Citrus Genome Research Articles

Comparative Analysis of Transposable Elements in the Genomes of Citrus and Citrus-Related Genera.

Transposable elements (TEs) significantly contribute to the evolution and diversity of plant genomes. In this study, we explored the roles of TEs in the genomes of Citrus and Citrus-related genera by constructing a pan-genome TE library from 20 published genomes of Citrus and Citrus-related accessions. Our results revealed an increase in TE content and the number of TE types compared to the original annotations, as well as a decrease in the content of unclassified TEs. The average length of TEs per assembly was approximately 194.23 Mb, representing 41.76% (Murraya paniculata) to 64.76% (Citrus gilletiana) of the genomes, with a mean value of 56.95%. A significant positive correlation was found between genome size and both the number of TE types and TE content. Consistent with the difference in mean whole-genome size (39.83 Mb) between Citrus and Citrus-related genera, Citrus genomes contained an average of 34.36 Mb more TE sequences than Citrus-related genomes. Analysis of the estimated insertion time and half-life of long terminal repeat retrotransposons (LTR-RTs) suggested that TE removal was not the primary factor contributing to the differences among genomes. These findings collectively indicate that TEs are the primary determinants of genome size and play a major role in shaping genome structures. Principal coordinate analysis (PCoA) of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) identifiers revealed that the fragmented TEs were predominantly derived from ancestral genomes, while intact TEs were crucial in the recent evolutionary diversification of Citrus. Moreover, the presence or absence of intact TEs near the AdhE superfamily was closely associated with the bitterness trait in the Citrus species. Overall, this study enhances TE annotation in Citrus and Citrus-related genomes and provides valuable data for future genetic breeding and agronomic trait research in Citrus.

Chromosome-Scale, De Novo, Phased Genome Assemblies of Three Australian Limes: Citrus australasica, C. inodora, and C. glauca.

Huanglongbing (HLB) is a severe citrus disease worldwide. Wild Australian limes like Citrus australasica, C. inodora, and C. glauca possess beneficial HLB resistance traits. Individual trees of the three taxa were extensively used in a breeding program for over a decade to introgress resistance traits into commercial-quality citrus germplasm. We generated high-quality, phased, de novo genome assemblies of the three Australian limes using PacBio long-read sequencing. The genome assembly sizes of the primary and alternate haplotypes were determined for C. australasica (337 Mb/335 Mb), C. inodora (304 Mb/299 Mb), and C. glauca (376 Mb/379 Mb). The nine chromosome-scale scaffolds included 86-91% of the genome sequences generated. The integrity and completeness of the assembled genomes were estimated to be at 97.2-98.8%. Gene annotation studies identified 25,461 genes in C. australasica, 27,665 in C. inodora, and 30,067 in C. glauca. Genes belonging to 118 orthogroups were specific to Australian lime genomes compared to other citrus genomes analyzed. Significantly fewer canonical resistance (R) genes were found in C. inodora and C. glauca (319 and 449, respectively) compared to C. australasica (576), C. clementina (579), and C. sinensis (651). Similar patterns were observed for other gene families associated with potential HLB resistance, including Phloem protein 2 (PP2) and Callose synthase (CalS) genes predicted in the Australian lime genomes. The genomic information on Australian limes developed in the present study will help understand the genetic basis of HLB resistance.

Multiplexed gene editing in citrus by using a multi-intron containing Cas9 gene.

Several expression systems have been developed in clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (CRISPR/Cas9) framework allowing for gene editing of disease-associated genes across diverse citrus varieties. In this study, we present a new approach employing a multi-intron containing Cas9 gene plus multiple gRNAs separated with tRNA sequences to target the phytoene desaturase gene in both 'Carrizo' citrange and 'Duncan' grapefruit. Notably, using this unified vector significantly boosted editing efficiency in both citrus varieties, showcasing mutations in all three designated targets. The implementation of this multiplex gene editing system with a multi-intron-containing Cas9 plus a gRNA-tRNA array demonstrates a promising avenue for efficient citrus genome editing, equipping us with potent tools in the ongoing battle against several diseases such as canker and huanglongbing.

Genome-Wide Analysis of the GLK Gene Family and Its Expression at Different Leaf Ages in the Citrus Cultivar Kanpei.

The GLK gene family plays a crucial role in the regulation of chloroplast development and participates in chlorophyll synthesis. However, the precise mechanism by which GLK contributes to citrus's chlorophyll synthesis remains elusive. The GLK gene family causes variations in the photosynthetic capacity and chlorophyll synthesis of different citrus varieties. In this study, we identified tissue-specific members and the key CcGLKs involved in chlorophyll synthesis. A total of thirty CcGLK transcription factors (TFs) were discovered in the citrus genome, distributed across all nine chromosomes. The low occurrence of gene tandem duplication events and intronic variability suggests that intronic variation may be the primary mode of evolution for CcGLK TFs. Tissue-specific expression patterns were observed for various GLK family members; for instance, CcGLK12 and CcGLK15 were specifically expressed in the skin, while CcGLK30 was specific to the ovary, and CcGLK10, CcGLK6, CcGLK21, CcGLK2, CcGLK18, CcGLK9, CcGLK28, and CcGLK8 were specifically expressed in the leaves. CcGLK4, CcGLK5, CcGLK11, CcGLK23, CcGLKl7, CcGLK26, and CcGLK20 may participate in the regulation of the ALA, prochlorophylate, protoporphyrin IX, Mg-protoporphyrin IX, Chl b, T-Chl, MG-ProtoIX ME, and POR contents in citrus.

AGRAMP: machine learning models for predicting antimicrobial peptides against phytopathogenic bacteria.

Antimicrobial peptides (AMPs) are promising alternatives to traditional antibiotics for combating plant pathogenic bacteria in agriculture and the environment. However, identifying potent AMPs through laborious experimental assays is resource-intensive and time-consuming. To address these limitations, this study presents a bioinformatics approach utilizing machine learning models for predicting and selecting AMPs active against plant pathogenic bacteria. N-gram representations of peptide sequences with 3-letter and 9-letter reduced amino acid alphabets were used to capture the sequence patterns and motifs that contribute to the antimicrobial activity of AMPs. A 5-fold cross-validation technique was used to train the machine learning models and to evaluate their predictive accuracy and robustness. The models were applied to predict putative AMPs encoded by intergenic regions and small open reading frames (ORFs) of the citrus genome. Approximately 7% of the 10,000-peptide dataset from the intergenic region and 7% of the 685,924-peptide dataset from the whole genome were predicted as probable AMPs. The prediction accuracy of the reported models range from 0.72 to 0.91. A subset of the predicted AMPs was selected for experimental test against Spiroplasma citri, the causative agent of citrus stubborn disease. The experimental results confirm the antimicrobial activity of the selected AMPs against the target bacterium, demonstrating the predictive capability of the machine learning models. Hydrophobic amino acid residues and positively charged amino acid residues are among the key features in predicting AMPs by the Random Forest Algorithm. Aggregation propensity appears to be correlated with the effectiveness of the AMPs. The described models would contribute to the development of effective AMP-based strategies for plant disease management in agricultural and environmental settings. To facilitate broader accessibility, our model is publicly available on the AGRAMP (Agricultural Ngrams Antimicrobial Peptides) server.

Biochemical characterization of a regiospecific flavonoid 3'-O-methyltransferase from orange

Citrus plants have diverse methoxyflavonoids including, chrysoeriol, isosakuranetin, and nobiletin. In plants, O-methyltransferases (OMTs) participate in the methylation of a vast array of secondary metabolites, including flavonoids, phenylpropanoids, and alkaloids. To identify functional OMTs involved in the formation of methoxyflavonoids, orange (Citrus sinensis) OMT (CsOMT) genes were retrieved from the Citrus Genome Database. The phylogenetic relationships with functional OMTs suggested that three CsOMTs, CsOMT15, CsOMT16, and CsOMT30, are possible candidates for flavonoid OMTs (FOMTs). These CsOMTs were heterologously expressed in Escherichia coli, and their OMT activity was examined with flavonoid substrates. Of the examined CsOMTs, CsOMT16 catalyzed the regiospecific 3'-O-methylation of flavonoids to the respective 3'-methoxyflavonoids. A kinetic study demonstrated that CsOMT16 accepts diverse flavonoids as a substrate with a comparable preference. The flavonoids eriodictyol, luteolin, and quercetin were efficiently converted to homoeriodictyol, chrysoeriol, and isorhamnetin by CsOMT16-transformed E. coli cells, respectively. These findings suggest that CsOMT16 contributes to the methoxyflavonoid formation in orange and is applicable to the biotechnological production of 3'-methoxyflavonoids.

Comprehensive analyses of the citrus WRKY gene family involved in the metabolism of fruit sugars and organic acids.

Sugars and organic acids are the main factors determining the flavor of citrus fruit. The WRKY transcription factor family plays a vital role in plant growth and development. However, there are still few studies about the regulation of citrus WRKY transcription factors (CsWRKYs) on sugars and organic acids in citrus fruit. In this work, a genome-wide analysis of CsWRKYs was carried out in the citrus genome, and a total of 81 CsWRKYs were identified, which contained conserved WRKY motifs. Cis-regulatory element analysis revealed that most of the CsWRKY promoters contained several kinds of hormone-responsive and abiotic-responsive cis-elements. Furthermore, gene expression analysis and fruit quality determination showed that multiple CsWRKYs were closely linked to fruit sugars and organic acids with the development of citrus fruit. Notably, transcriptome co-expression network analysis further indicated that three CsWRKYs, namely, CsWRKY3, CsWRKY47, and CsWRKY46, co-expressed with multiple genes involved in various pathways, such as Pyruvate metabolism and Citrate cycle. These CsWRKYs may participate in the metabolism of fruit sugars and organic acids by regulating carbohydrate metabolism genes in citrus fruit. These findings provide comprehensive knowledge of the CsWRKY family on the regulation of fruit quality.

Proteomic landscape presents cues for vegetative to reproductive transition in mango

Proteome-based vegetative and flower bud formation characterization was utilized to identify and differentiate protein species with significant variable abundance during floral transition in mango cv. Dashehari using 2DE and corroborated the identified protein spots using gene expression analysis. Total soluble proteins were phenol-extracted from mango cv's vegetative and floral flush. Dashehari and separated on 2D gels at pH 4-7. The protein spots with variable intensity were identified through SameSpots software. The protein sequences of differentially accumulated spots were identified based on PI and MW using Citrus sinensis proteome isoelectric focusing database. Furthermore, these protein sequences were used to conduct (tBLASTn) against Mangifera indica to predict the protein. Real time gene expression was done to corroborate identified proteins. Total 301 spots were detected, out of which 16 were identified as differentially expressed (P?0.05) and a 2-fold change. These 16 protein spots were identified on the basis of in silico comparative mapping protein against genome of mango and citrus, a close relative. They were classified into eight categories: transcriptional regulation, phenylpropanoid pathway and cell wall /cytoskeleton metabolism-related proteins, hormone signalling, flowering time, signal-transduction, transport and protein synthesis to flowering. Five genes coding for shortlisted proteins were used for validation of results using gene expression analysis. SAM (S adenosyl methonine synthase) was found up-regulated in floral flush, involved in the biosynthesis of polyamines has association with flowering, and stress responses. Furthermore, ARF (Auxin Response Factor), serine/threonine kinase gene members were also found to play critical role in determining floral development process in mango, consistent with results obtained through 2DE. Protein species that are putatively involved in phenylpropanoid pathway were also identified, showing the process of mango flowering from a new perspective beyond the conventional view. This flowering related proteomics study provides an overview of the biological pathways and regulatory mechanisms associated with flowering developmental physiology.

Mining key genes related to root morphogenesis through genome-wide identification and expression analysis of RR gene family in citrus.

Morphogenesis of root is a vital factor to determine the root system architecture. Cytokinin response regulators (RRs) are the key transcription factors in cytokinin signaling, which play important roles in regulating the root morphogenesis. In this study, 29 RR proteins, including 21 RRs and 8 pseudo RRs, were identified from the genome of citrus, and termed as CcRR1-21 and CcPRR1-8, respectively. Phylogenetic analysis revealed that the 29 CcRRs could be classified into four types according to their representative domains. Analysis of cis-elements of CcRRs indicated that they were possibly involved in the regulation of growth and abiotic stress resistance in citrus. Within the type A and type B CcRRs, CcRR4, CcRR5, CcRR6 and CcRR16 highly expressed in roots and leaves, and dramatically responded to the treatments of hormones and abiotic stresses. CcRR2, CcRR10, CcRR14 and CcRR19 also highly expressed in roots under different treatments. Characteristic analysis revealed that the above 8 CcRRs significantly and differentially expressed in the three zones of root, suggesting their functional differences in regulating root growth and development. Further investigation of the 3 highly and differentially expressed CcRRs, CcRR5, CcRR10 and CcRR14, in 9 citrus rootstocks showed that the expression of CcRR5, CcRR10 and CcRR14 was significantly correlated to the length of primary root, the number of lateral roots, and both primary root and the number of lateral roots, respectively. Results of this study indicated that CcRRs were involved in regulating the growth and development of the root in citrus with different functions among the members.

Genome-Wide Analysis of the GLK Gene Family and the Expression under Different Growth Stages and Dark Stress in Sweet Orange (Citrus sinensis)

The GLK gene family is of great significance in regulating chloroplast development and participating in chlorophyll synthesis. However, the mechanism of GLK involvement in Citrus chlorophyll synthesis remains unclear. In this study, bioinformatics methods were used to analyze the gene structure, protein evolution, chromosome distribution, promoter elements and expression profile of GLK gene family in Citrus. Overall, 27 CsGLK TFs were identified from Citrus genome and divided into three subgroups according to the conserved domains. All members were distributed on nine chromosomes. The tandem replication events (ka/ks < 1) indicated that CsGLK TFs underwent a purification selection evolutionary process. The intron variation might be a vital configuration for the evolution of CsGLK genes. The expression pattern of CsGLKs showed that family members had higher expression levels in different tissues and at different growth stages and could actively respond to dark stress. CsGLK TFs of the same group had similar structures, but their expression patterns were quite different, indicating that they may have different functions and not be redundant. Correlation analysis showed that CsGLK2, CsGLK9, CsGLK10, CsGLK11, CsGLK20 and CsGLK24 were significantly positive correlations with Chl a and Chl b contents. In addition, CsGLK2, CsGLK5, CsGLK10, CsGLK11, CsGLK12, CsGLK15, CsGLK20 and CsGLK24 were significantly positive related to Mg-Proto IX, Proto IX and Pchl.

Genome-wide identification of the CsPAL gene family and functional analysis for strengthening green mold resistance in citrus fruit

Phenylalanine aminolyase (PAL) is a key rate-limiting enzyme. It catalyzes the initial step of the phenylpropane metabolic pathway. PAL is involved in the regulation of fruit diseases, but the specific function of the PAL gene in the regulation of postharvest fruit diseases is unknown. In this study, we identified four CsPAL members in citrus, including CsPAL1-CsPAL4. Based on physicochemical properties, gene structure, protein structure, protein sequence, and phylogenetic relationship analysis, all four CsPALs are predicted to have catalytic activity and take part in the disease resistance of citrus fruit to green mold. In addition, using Agrobacterium-mediated transient overexpression, all four CsPALs were shown to improve the tolerance of citrus fruit to green mold. Measurements of secondary metabolites in peel samples from the four transient overexpression groups showed a significant increase in both flavonoids and total phenols and a significant decrease in lignin content. For the first time, four CsPAL members were identified from the citrus genome, demonstrating that they are all implicated in the gene function of citrus fruit resistance to green mold through the elevation of flavonoids and total phenols but not lignin.

Genome-wide identification and expression analysis of MATE gene family in citrus fruit (Citrus clementina)

Multidrug and toxic compound extrusion (MATE) proteins are a class of secondary active multidrug transporters. In plants, this family has significantly expanded and is involved in numerous plant physiological processes. Although MATE proteins have been identified in an increasing number of species, the understanding about this family in citrus remains unclear. In this study, a total of 69 MATE transporters were identified in the citrus genome (Citrus clementina) and classified into four groups by phylogenetic analysis. Tandem and segmental duplication events were the main causes of the citrus MATE family expansion. RNA-seq and qRT-PCR analyses were performed during citrus fruit development. The results indicated that CitMATE genes showed specific expression profiles in citrus peels and flesh at different developmental stages. Combined with the variations of flavonoids and citrate levels in citrus fruit, we suggested that CitMATE43 and CitMATE66 may be involved in the transport process of flavonoids and citrate in citrus fruit, respectively. In addition, two flavonoids positive regulators, CitERF32 and CitERF33, both directly bind to and activated the CitMATE43 promoter. Our results provide comprehensive information on citrus MATE genes and valuable understanding for the flavonoids and citrate metabolism in citrus fruit.

Genome-wide identification of ovate family in Citrus and functional characterization of CitOFP19

Fruit shape is an important trait for fruit appearance and commercial value. Diversity of fruit-shape has been utilized in the breeding of pummelo (Citrus maxima), a basic species in Citrus. However, little is known about genetic basis of fruit shape in citrus. In this study, we identified 16 OVATE family protein (OFP) genes in the pummelo genome. Phylogenetically, they were classified into three subfamilies, which was consistent with the classification of their Arabidopsis orthologs. Synteny analysis suggested that segment and tandem duplications were responsible for their expansion in pummelo. Expression pattern analysis of Citrus OFPs (CitOFPs) showed that CitOFP19 had significantly higher expression level in the ovaries of round pummelo than in those of pear-shaped pummelo. Heterologous overexpression of CitOFP19 in tomato resulted in pear-shaped ovary and fruit shape. Taken together, this study characterized OVATE gene family in Citrus genome and assessed the function of CitOFP19.

Not Just a Cycle: Three gab Genes Enable the Non-Cyclic Flux Toward Succinate via GABA Shunt in 'Candidatus Liberibacter asiaticus'-Infected Citrus.

Although the mitochondria retain all required enzymes for an intact tricarboxylic acid (TCA) cycle, plants might shift the cyclic flux from the TCA cycle to an alternative noncyclic pathway via γ-aminobutyric acid (GABA) shunt under specific physiological conditions. We hypothesize that several genes may ease this noncyclic flux and contribute to the citrus response to the phytopathogenic bacterium 'Candidatus Liberibacter asiaticus', the causal agent of Huanglongbing in citrus. To test this hypothesis, we used multiomics techniques (metabolomics, fluxomics, and transcriptomics) to investigate the potential roles of putative gab homologies from Valencia sweet orange (Citrus sinensis). Our findings showed that 'Ca. L. asiaticus' significantly increased the endogenous GABA and succinate content but decreased ketoglutarate in infected citrus plants. Citrus genome harbors three putative gab genes, including amino-acid permease (also known as GABA permease; CsgabP), GABA transaminase (CsgabT), and succinate-semialdehyde dehydrogenase (also known as GABA dehydrogenase; CsgabD). The transcript levels of CsgabP, CsgabT, and CsgabD were upregulated in citrus leaves upon the infection with 'Ca. L. asiaticus' and after the exogenous application of GABA or deuterium-labeled GABA isotope (GABA-D6). Moreover, our finding showed that exogenously applied GABA is quickly converted to succinate and fed into the TCA cycle. Likewise, the fluxomics study showed that GABA-D6 is rapidly metabolized to succinate-D4. Our work proved that GABA shunt and three predicated gab genes from citrus, support the upstream noncyclic flux toward succinate rather than an intact TCA cycle and contribute to citrus defense responses to 'Ca. L. asiaticus'.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Highly Efficient Generation of Canker-Resistant Sweet Orange Enabled by an Improved CRISPR/Cas9 System.

Sweet orange (Citrus sinensis) is the most economically important species for the citrus industry. However, it is susceptible to many diseases including citrus bacterial canker caused by Xanthomonas citri subsp. citri (Xcc) that triggers devastating effects on citrus production. Conventional breeding has not met the challenge to improve disease resistance of sweet orange due to the long juvenility and other limitations. CRISPR-mediated genome editing has shown promising potentials for genetic improvements of plants. Generation of biallelic/homozygous mutants remains difficult for sweet orange due to low transformation rate, existence of heterozygous alleles for target genes, and low biallelic editing efficacy using the CRISPR technology. Here, we report improvements in the CRISPR/Cas9 system for citrus gene editing. Based on the improvements we made previously [dicot codon optimized Cas9, tRNA for multiplexing, a modified sgRNA scaffold with high efficiency, citrus U6 (CsU6) to drive sgRNA expression], we further improved our CRISPR/Cas9 system by choosing superior promoters [Cestrum yellow leaf curling virus (CmYLCV) or Citrus sinensis ubiquitin (CsUbi) promoter] to drive Cas9 and optimizing culture temperature. This system was able to generate a biallelic mutation rate of up to 89% for Carrizo citrange and 79% for Hamlin sweet orange. Consequently, this system was used to generate canker-resistant Hamlin sweet orange by mutating the effector binding element (EBE) of canker susceptibility gene CsLOB1, which is required for causing canker symptoms by Xcc. Six biallelic Hamlin sweet orange mutant lines in the EBE were generated. The biallelic mutants are resistant to Xcc. Biallelic mutation of the EBE region abolishes the induction of CsLOB1 by Xcc. This study represents a significant improvement in sweet orange gene editing efficacy and generating disease-resistant varieties via CRISPR-mediated genome editing. This improvement in citrus genome editing makes genetic studies and manipulations of sweet orange more feasible.

On how to generalize species-specific conceptual schemes to generate a species-independent Conceptual Schema of the Genome

BackgroundUnderstanding the genome, with all of its components and intrinsic relationships, is a great challenge. Conceptual modeling techniques have been used as a means to face this challenge. The heterogeneity and idiosyncrasy of genomic use cases mean that conceptual modeling techniques are used to generate conceptual schemes that focus on too specific scenarios (i.e., they are species-specific conceptual schemes). Our research group developed two different conceptual schemes. The first one is the Conceptual Schema of the Human Genome, which is intended to improve Precision Medicine and genetic diagnosis. The second one is the Conceptual Schema of the Citrus Genome, which is intended to identify the genetic cause of relevant phenotypes in the agri-food field.MethodsOur two conceptual schemes have been ontologically compared to identify their similarities and differences. Based on this comparison, several changes have been performed in the Conceptual Schema of the Human Genome in order to obtain the first version of a species-independent Conceptual Schema of the Genome. Identifying the different genome information items used in each genomic case study has been essential in achieving our goal. The changes needed to provide an expanded, more generic version of the Conceptual Schema of the Human Genome are analyzed and discussed.ResultsThis work presents a new CS called the Conceptual Schema of the Genome that is ready to be adapted to any specific working genome-based context (i.e., species-independent).ConclusionThe generated Conceptual Schema of the Genome works as a global, generic element from which conceptual views can be created in order to work with any specific species. This first working version can be used in the human use case, in the citrus use case, and, potentially, in more use cases of other species.

The Breeding Information Management System (BIMS): an online resource for crop breeding.

In this era of big data, breeding programs are producing ever larger amounts of data. This necessitates access to efficient management systems to keep track of cross, performance, pedigree, geographical and image-based data, as well as genotyping data. In this article, we report the progress on the Breeding Information Management System (BIMS), a free, secure and online breeding management system that allows breeders to store, manage, archive and analyze their private breeding data. BIMS is the first publicly available database system that enables individual breeders to integrate their private phenotypic and genotypic data with public data and, at the same time, have complete control of their own breeding data along with access to tools such as data import/export, data analysis and data archiving. The integration of breeding data with publicly available genomic and genetic data enhances genetic understanding of important traits and maximizes the marker-assisted breeding utility for breeders and allied scientists. BIMS incorporates the use of the Android App Field Book, open-source phenotype data collection software for phones and tablets that allows breeders to replace hard copy field books, thus alleviating the possibility of transcription errors while providing faster access to the collected data. BIMS comes with training materials and support for individual or small group training and is currently implemented in the Genome Database for Rosaceae, CottonGEN, the Citrus Genome Database, the Pulse Crop Database, and the Genome Database for Vaccinium.Database URLs: (https://www.rosaceae.org/), (https://www.cottongen.org/), (https://www.citrusgenomedb.org/), (https://www.pulsedb.org/) and (https://www.vaccinium.org/)

High-Density Genetic Map Construction and Identification of QTLs Controlling Leaf Abscission Trait in Poncirus trifoliata.

A high-density genetic linkage map is essential for genetic and genomic studies including QTL mapping, genome assembly, and comparative genomic analysis. Here, we constructed a citrus high-density linkage map using SSR and SNP markers, which are evenly distributed across the citrus genome. The integrated linkage map contains 4163 markers with an average distance of 1.12 cM. The female and male linkage maps contain 1478 and 2976 markers with genetic lengths of 1093.90 cM and 1227.03 cM, respectively. Meanwhile, a genetic map comparison demonstrates that the linear order of common markers is highly conserved between the clementine mandarin and Poncirus trifoliata. Based on this high-density integrated citrus genetic map and two years of deciduous phenotypic data, two loci conferring leaf abscission phenotypic variation were detected on scaffold 1 (including 36 genes) and scaffold 8 (including 107 genes) using association analysis. Moreover, the expression patterns of 30 candidate genes were investigated under cold stress conditions because cold temperature is closely linked with the deciduous trait. The developed high-density genetic map will facilitate QTL mapping and genomic studies, and the localization of the leaf abscission deciduous trait will be valuable for understanding the mechanism of this deciduous trait and citrus breeding.

Association of T2/S-RNase With Self-Incompatibility of Japanese Citrus Accessions Examined by Transcriptomic, Phylogenetic, and Genetic Approaches.

Several citrus varieties show gametophytic self-incompatibility (GSI), which can contribute to seedless fruit production in several cultivars. This study investigated the genes regulating this trait through RNA-seq performed using styles collected from the flowers of Japanese citrus cultivars ‘Hyuganatsu,' ‘Tosabuntan,' ‘Hassaku,' ‘Banpeiyu,' and ‘Sweet Spring'. We screened the transcripts of putative T2 RNases, i.e., the protein family including all S-RNases from S-RNase-based GSI plants, and constructed a phylogenetic tree using the screened T2 RNases and S-RNases retrieved from citrus genome databases and a public database. Three major clusters (class I–III) were formed, among which, the class III cluster contained family specific subclusters formed by S-RNase and a citrus-specific cluster monophyletic to the S-RNase clusters. From the citrus class III cluster, six transcripts were consistent with the S haplotypes previously determined in Japanese citrus accessions, sharing characteristics such as isoelectric point, extracellular localization, molecular weight, intron number and position, and tissue-specific expression with S-RNases. One T2 RNase gene in self-incompatible Hyuganatsu was significantly down-regulated in the styles of a self-compatible mutant of Hyuganatsu in RNA-seq and qPCR analyses. In addition, the inheritance pattern of some T2 RNase genes was consistent with the pattern of the S haplotype in the progeny population of Hyuganatsu and Tosabuntan. As all results supported citrus self-incompatibility being based on S-RNase, we believe that six T2 RNase genes were S-RNases. The homology comparison between the six T2 RNases and S-RNases recently reported in Chinese citrus revealed that three out of six T2 RNases were identical to S-RNases from Chinese citrus. Thus, the other three T2 RNases were finally concluded to be novel citrus S-RNases involved in self-incompatibility.

Genome-wide identification of citrus histone acetyltransferase and deacetylase families and their expression in response to arbuscular mycorrhizal fungi and drought

ABSTRACT Histone acetyltransferase (HAT) and histone deacetylase (HDAC) family members control histone acetylation levels to regulate gene expression and play significant roles in microbe interactions and abioticstress responses. However, the HAT and HDAC family members in citrus have not yet been identified, and their responses to drought (abiotic stress) and arbuscular mycorrhizal fungi (AMF) colonization, remain unknown. In the present study, HAT and HDAC family genes were characterized in the citrus genome. A total of 14 CsHATs and 11 CsHDACs were identified in the citrus genome, of which 8 CsHATs and 7 CsHDACs showed syntenic relationships with Arabidopsis HATs and HDACs, respectively. Furthermore, the 14 CsHATs were classified into GNAT-MYST, TAFII250, MYST, CBP, and one citrus-unique family. The 11 CsHDACs were classified into SIR2, HD2, and RPD3/HDA1. Gene expression analysis showed AMF colonizationor drought induced the expression of CsHAT1, 2, 6, 7, 9, 11, 12, and 14, and repressed CsHDAC1, 2, 5, 6, 8, and 10. The regulation of HAT and HDAC expression was not always consistent with the cis-elements identified in their promoters. The present study offers useful insights into the citrus HAT and HDAC gene families and their expression in response to AMF colonization and drought stress.