Background: Iron (Fe) deficiency is a common problem in citrus production. As the second largest superfamily of transcription factors (TFs), the basic/helix-loop-helix (bHLH) proteins have been shown to participate in the regulation of Fe homeostasis and a series of other biological and developmental processes in plants. However, this family of members in citrus and their functions in citrus Fe deficiency are still largely unknown. Results: In this study, we identified a total of 128 CgbHLHs from pummelo (Citrus grandis) genome that were classified into 18 subfamilies by phylogenetic comparison with Arabidopsis thaliana bHLH proteins. All of these CgbHLHs were randomly distributed on nine known (125 genes) and one unknown (3 genes) chromosomes, and 12 and 47 of them were identified to be tandem and segmental duplicated genes, respectively. Sequence analysis showed detailed characteristics of their intron-exon structures, bHLH domain and conserved motifs. Gene ontology (GO) analysis suggested that most of CgbHLHs were annotated to the nucleus, DNA-binding transcription factor activity, response to abiotic stimulus, reproduction, post-embryonic development, flower development and photosynthesis. In addition, 27 CgbHLH proteins were predicted to have direct or indirect protein-protein interactions. Based on GO annotation, RNA sequencing data in public database and qRT-PCR results, several of CgbHLHs were identified as the key candidates that respond to iron deficiency. Conclusions: In total, 128 CgbHLH proteins were identified from pummelo, and their detailed sequence and structure characteristics and putative functions were analyzed. This study provides comprehensive information for further functional elucidation of CgbHLH genes in citrus.

BackgroundCopper (Cu) toxicity has become a potential threat for citrus production, but little is known about related mechanisms. This study aims to uncover the global landscape of mRNAs, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in response to Cu toxicity so as to construct a regulatory network of competing endogenous RNAs (ceRNAs) and to provide valuable knowledge pertinent to Cu response in citrus.ResultsTolerance of four commonly used rootstocks to Cu toxicity was evaluated, and ‘Ziyang Xiangcheng’ (Citrus junos) was found to be the most tolerant genotype. Then the roots and leaves sampled from ‘Ziyang Xiangcheng’ with or without Cu treatment were used for whole-transcriptome sequencing. In total, 5734 and 222 mRNAs, 164 and 5 lncRNAs, 45 and 17 circRNAs, and 147 and 130 miRNAs were identified to be differentially expressed (DE) in Cu-treated roots and leaves, respectively, in comparison with the control. Gene ontology enrichment analysis showed that most of the DEmRNAs and targets of DElncRNAs and DEmiRNAs were annotated to the categories of ‘oxidation-reduction’, ‘phosphorylation’, ‘membrane’, and ‘ion binding’. The ceRNA network was then constructed with the predicted pairs of DEmRNAs-DEmiRNAs and DElncRNAs-DEmiRNAs, which further revealed regulatory roles of these DERNAs in Cu toxicity.ConclusionsA large number of mRNAs, lncRNAs, circRNAs, and miRNAs in ‘Ziyang Xiangcheng’ were altered in response to Cu toxicity, which may play crucial roles in mitigation of Cu toxicity through the ceRNA regulatory network in this Cu-tolerant rootstock.