Structural organization,classification and phylogenetic relationship of cytochrome P450 genes in Citrus clementina and Citrus sinensis

The genus Citrus is an important fruit crop and nutritional source for the good health of humans. Cytochrome P450s represent about 1 % of the proteome and mediate diverse biochemical reactions pertaining to both primary and secondary metabolism. Analysis of Citrus genomic resources identified 296 plant cytochrome P450s (CYP) coding genes in Citrus clementina, 272 in double haploid (dh) Citrus sinensis, and 202 in C. sinensis. In C. clementina and dh C. sinensis, CYP genes are distributed into nine clans. In the three genomes, single intron containing CYP genes are predominant in the A-type families. Among non-A-type CYP families, multiple intron containing genes are predominant. More number of genes in CYP A-type families over non-A-type families is attributed to rapid evolution of A-type genes facilitated by their gene organization. Further, complex gene organization of non-A-type genes with the presence of multiple introns might have contributed to the slower evolvement of paralogs. Majority of introns (1,660) from three genomes showed canonical GT-AG splice sites. However, 33 introns showed non-conventional GC… PyAG splice sites and functionality of these splice sites is confirmed by the ESTs lacking this intron. Across the families, gene organization is conserved between the three genomes. In dh C. sinensis, 22 genes were identified to have alternate splicing. Examination of scaffolds in C. clementina revealed that majority of the Citrus CYP genes are solitary and a few of them are in clusters of 3–8 genes. PCR amplification of C. sinensis genomic DNA with gene-specific primers failed to amplify out-grouped genes Ccl-CYP706A16 and Ccl-CYP706B1, confirming that they are specific to C. clementina. Differential number of CYP genes observed between C. clementina and C. sinensis is attributed to the extent of variability between their parents representing ancestral taxa.     

Complete Chloroplast Genome Sequence of Omani Lime (Citrus aurantiifolia) and Comparative Analysis within the Rosids

The genus Citrus contains many economically important fruits that are grown worldwide for their high nutritional and medicinal value. Due to frequent hybridizations among species and cultivars, the exact number of natural species and the taxonomic relationships within this genus are unclear. To compare the differences between the Citrus chloroplast genomes and to develop useful genetic markers, we used a reference-assisted approach to assemble the complete chloroplast genome of Omani lime (C. aurantiifolia). The complete C. aurantiifolia chloroplast genome is 159,893 bp in length; the organization and gene content are similar to most of the rosids lineages characterized to date. Through comparison with the sweet orange (C. sinensis) chloroplast genome, we identified three intergenic regions and 94 simple sequence repeats (SSRs) that are potentially informative markers with resolution for interspecific relationships. These markers can be utilized to better understand the origin of cultivated Citrus. A comparison among 72 species belonging to 10 families of representative rosids lineages also provides new insights into their chloroplast genome evolution.     

Use of an expressed sequence tag-based method for single nucleotide polymorphism identification and discrimination of Citrus species and cultivars

Sweet orange [Citrus sinensis (L.) Osbeck] represents the most important Citrus species, followed by clementine (C. clementina Hort. ex Tan.). Citrus species and genotypes are difficult to recognize as they have a moderate level of diversity due to nucellar selection, vegetative propagation and origin by single spontaneous mutation. Despite the large number of available sequences and the existence of a draft assembly of sweet orange and clementine, there are currently no single nucleotide polymorphism (SNP) databases for Citrus species. For this purpose, the QualitySNP software was used to discover SNPs in 19 Citrus species starting from 540,000 expressed sequence tags (ESTs) assembled in 52,000 contigs. The vast majority of ESTs, contigs and SNPs were found in C. clementina and C. sinensis: 4,400 out of 16,000 contigs (27 %) of C. clementina and 4,100 out of 17,000 contigs (24 %) of C. sinensis contained putative SNPs. A total of 3,634 sequences were associated with enzymes belonging to 121 metabolic KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, among which the secondary metabolite pathway was the most represented. A total of 163 SNPs from 52 contigs and genes of specific functional categories were validated and 81 polymorphic sites were found. Thirty-seven selected SNPs, validated by Sanger sequencing, confirmed that polymorphisms were mainly between species, while poor within-species variability was discovered. This work provides a collection of 15,879 putative SNP markers that could be exploited by the Citrus community. Furthermore, the validated SNPs associated with specific genes could be used for functional genetic studies in germplasm diversity analysis, mapping and breeding.     

Phylogenetic Relationships of Citrus and Its Relatives Based on matK Gene Sequences

The genus Citrus includes mandarin, orange, lemon, grapefruit and lime, which have high economic and nutritional value. The family Rutaceae can be divided into 7 subfamilies, including Aurantioideae. The genus Citrus belongs to the subfamily Aurantioideae. In this study, we sequenced the chloroplast matK genes of 135 accessions from 22 genera of Aurantioideae and analyzed them phylogenetically. Our study includes many accessions that have not been examined in other studies. The subfamily Aurantioideae has been classified into 2 tribes, Clauseneae and Citreae, and our current molecular analysis clearly discriminate Citreae from Clauseneae by using only 1 chloroplast DNA sequence. Our study confirms previous observations on the molecular phylogeny of Aurantioideae in many aspects. However, we have provided novel information on these genetic relationships. For example, inconsistent with the previous observation, and consistent with our preliminary study using the chloroplast rbcL genes, our analysis showed that Feroniella oblata is not nested in Citrus species and is closely related with Feronia limonia. Furthermore, we have shown that Murraya paniculata is similar to Merrillia caloxylon and is dissimilar to Murraya koenigii. We found that “true citrus fruit trees” could be divided into 2 subclusters. One subcluster included Citrus, Fortunella, and Poncirus, while the other cluster included Microcitrus and Eremocitrus. Compared to previous studies, our current study is the most extensive phylogenetic study of Citrus species since it includes 93 accessions. The results indicate that Citrus species can be classified into 3 clusters: a citron cluster, a pummelo cluster, and a mandarin cluster. Although most mandarin accessions belonged to the mandarin cluster, we found some exceptions. We also obtained the information on the genetic background of various species of acid citrus grown in Japan. Because the genus Citrus contains many important accessions, we have comprehensively discussed the classification of this genus.     

Genome-wide identification of NBS genes in<Emphasis Type="Italic"> japonica</Emphasis> rice reveals significant expansion of divergent non-TIR NBS-LRR genes

A complete set of candidate disease resistance ( R) genes encoding nucleotide-binding sites (NBSs) was identified in the genome sequence of japonica rice ( Oryza sativa L. var. Nipponbare). These putative R genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution, and compared with those in Arabidopsis thaliana. We found 535 NBS-coding sequences, including 480 non-TIR (Toll/IL-1 receptor) NBS-LRR (Leucine Rich Repeat) genes. TIR NBS-LRR genes, which are common in A. thaliana, have not been identified in the rice genome. The number of non-TIR NBS-LRR genes in rice is 8.7 times higher than that in A. thaliana, and they account for about 1% of all of predicted ORFs in the rice genome. Some 76% of the NBS genes were located in 44 gene clusters or in 57 tandem arrays, and 16 apparent gene duplications were detected in these regions. Phylogenetic analyses based both NBS and N-terminal regions classified the genes into about 200 groups, but no deep clades were detected, in contrast to the two distinct clusters found in A. thaliana. The structural and genetic diversity that exists among NBS-LRR proteins in rice is remarkable, and suggests that diversifying selection has played an important role in the evolution of R genes in this agronomically important species. (Supplemental material is available online at .)Communicated by R. HagemannThe first three authors contributed equally to this work     

Nuclear Species-Diagnostic SNP Markers Mined from 454 Amplicon Sequencing Reveal Admixture Genomic Structure of Modern Citrus Varieties

Most cultivated Citrus species originated from interspecific hybridisation between four ancestral taxa (C. reticulata, C. maxima, C. medica, and C. micrantha) with limited further interspecific recombination due to vegetative propagation. This evolution resulted in admixture genomes with frequent interspecific heterozygosity. Moreover, a major part of the phenotypic diversity of edible citrus results from the initial differentiation between these taxa. Deciphering the phylogenomic structure of citrus germplasm is therefore essential for an efficient utilization of citrus biodiversity in breeding schemes. The objective of this work was to develop a set of species-diagnostic single nucleotide polymorphism (SNP) markers for the four Citrus ancestral taxa covering the nine chromosomes, and to use these markers to infer the phylogenomic structure of secondary species and modern cultivars. Species-diagnostic SNPs were mined from 454 amplicon sequencing of 57 gene fragments from 26 genotypes of the four basic taxa. Of the 1,053 SNPs mined from 28,507 kb sequence, 273 were found to be highly diagnostic for a single basic taxon. Species-diagnostic SNP markers (105) were used to analyse the admixture structure of varieties and rootstocks. This revealed C. maxima introgressions in most of the old and in all recent selections of mandarins, and suggested that C. reticulata × C. maxima reticulation and introgression processes were important in edible mandarin domestication. The large range of phylogenomic constitutions between C. reticulata and C. maxima revealed in mandarins, tangelos, tangors, sweet oranges, sour oranges, grapefruits, and orangelos is favourable for genetic association studies based on phylogenomic structures of the germplasm. Inferred admixture structures were in agreement with previous hypotheses regarding the origin of several secondary species and also revealed the probable origin of several acid citrus varieties. The developed species-diagnostic SNP marker set will be useful for systematic estimation of admixture structure of citrus germplasm and for diverse genetic studies.     

LTR retrotransposons from the <Emphasis Type="Italic">Citrus x clementina</Emphasis> genome: characterization and application

Long terminal repeat retrotransposons (LTR-RTs) are a large portion of most plant genomes, and can be used as a powerful molecular marker system. The first citrus reference genome (Citrus x clementina) has been publicly available since 2011; however, previous studies in citrus have not utilized the whole genome for LTR-RT marker development. In this study, 3959 full-length LTR-RTs were identified in the C. x clementina genome using structure-based (LTR_FINDER) and homology-based (RepeatMasker) methods. LTR-RTs were first classified by protein domain into Gypsy and Copia superfamilies, and then clustered into 1074 families based on LTR sequence similarity. Three hundred fifty Copia families were grouped into four lineages: Retrofit, Tork, Sire, and Oryco. One hundred seventy-eight Gypsy families were sorted into six lineages: Athila, Tat, Renia, CRM, Galadriel, and Del. Most LTR-RTs (3218 or 81.3%) were anchored to the nine Clementine mandarin linkage groups, accounting for 9.74% of chromosomes currently assembled. Accessions of 25 Rutaceae species were genotyped using 17 inter-retrotransposon amplified polymorphism (IRAP) markers developed from conserved LTR regions. Sequence-specific amplified polymorphism (SSAP) makers were used to distinguish ‘Valencia’ and ‘Pineapple’ sweet oranges (C. x sinensis), and 24 sweet orange clones. LTR-RT markers developed from the Clementine genome can be transferred within the Rutaceae family demonstrating that they are an excellent tool for citrus and Rutaceae genetic analysis.     

Phylogeny and genomic organization of the TIR and non-tIR NBS-LRR resistance gene family in Medicago truncatula

Sequences homologous to the nucleotide binding site (NBS) domain of NBS-leucine-rich repeat (LRR) resistance genes were retrieved from the model legume M. truncatula through several methods. Phylogenetic analysis classified these sequences into TIR (toll and interleukin-1 receptor) and non-TIR NBS subfamilies and further subclassified them into several well-defined clades within each subfamily. Comparison of M. truncatula NBS sequences with those from several closely related legumes, including members of the tribes Trifoleae, Viceae, and Phaseoleae, reveals that most clades contain sequences from multiple legume species. Moreover, sequences from species within the closely related Trifoleae and Viceae tribes (e.g., Medicago and Pisum spp.) tended to be cophyletic and distinct from sequences of Phaseoleae species (e.g., soybean and bean). These results suggest that the origin of major clades within the NBS-LRR family predate radiation of these Papilionoid legumes, while continued diversification of these sequences mirrors speciation within this legume subfamily. Detailed genetic and physical mapping of both TIR and non-TIR NBS sequences in M. truncatula reveals that most NBS sequences are organized into clusters, and few, if any, clusters contain both TIR and non-TIR sequences. Examples were found, however, of physical clusters that contain sequences from distinct phylogenetic clades within the TIR or non-TIR subfamilies. Comparative mapping reveals several blocks of resistance gene loci that are syntenic between M. truncatula and soybean and between M. truncatula and pea.     

Cytoplasmic diversity,phylogenetic relationships and molecular evolution of Tunisian Citrus species as inferred from mutational events and pseudogene of chloroplast trnL-trnF spacer

The genus Citrus L. is among the most important fruit trees in the world. In this report, cytoplasmic polymorphism of twenty seven Tunisian Citrus cultivars was explored using the chloroplast trnL-trnF intergenic spacer. Chloroplast sequences showed variation in length and nucleotide content. Haplotype and nucleotide diversity showed low variations. Molecular phylogenetic tree identifies Citrus maternal origins and demonstrates two major groups distinguishing between mandarin and pummelo groups. The trnL-trnF intergenic spacer showed one copy of pseudogene of the original trnF gene in 27 Citrus species at position 275 bp with a size varying from 49 to 63 bp. The anticodon domain was identified as the most conserved element, but one transversion (T−>C) was found in the D-domain. Meanwhile, one transversion (T−>A) and one transition (T−>G) were found in the T-domain. Neutrality tests (Tajima, Fu & Li and Fu) which revealed positive and non-significant values and Pi and θ_W assume a neutral model of evolution and advocated a constant population size. The study demonstrates the resolving power of trnL-trnF sequence data to prove both pummelo and mandarin gene pool’s contribution in the development of Tunisian secondary species and inferring their genetic and phylogenetic relationships.     

Genome-wide identification and expression analysis of the molecular chaperone binding protein <Emphasis Type="Italic">BiP</Emphasis> genes in <Emphasis Type="Italic">Citrus</Emphasis>

Here, we report for the first time the genome-wide identification and expression analysis of the molecular chaperone BiP genes in Citrus. Six genes encoding the conserved protein domain family GPR78/BiP/KAR2 were identified in the genome of Citrus sinensis and C. clementina. Two of them, named here as CsBiP1 and CsBiP2, were classified as true BiPs based on their deduced amino acid sequences. Alignment of the deduced amino acid sequences of CsBiP1 and CsBiP2 with BiP homologs from soybean and Arabidopsis showed that they contain all the conserved functional motifs of BiPs. Analysis of the promoter region of CsBiPs revealed the existence of cis-acting regulatory sequences involved in abiotic, heat-shock, and endoplasmic reticulum (ER) stress responses. Publicly available RNA-seq data indicated that CsBiP1 is abundantly expressed in leaf, flower, fruit, and callus, whereas CsBiP2 expression is rarely detected in any tissues under normal conditions. Comparative quantitative real-time PCR (qPCR) analysis of expression of these genes between C. sinensis grafted on the drought-tolerant “Rangpur” lime (C. limonia) and -sensitive “Flying Dragon” trifoliate orange (Poncirus trifoliata) rootstocks showed that CsBiP1 was upregulated by drought stress on the former but downregulated on the latter, whereas the CsBiP2 mRNA levels were downregulated on drought-stressed “Flying Dragon,” but remained constant on “Rangpur.” CsBiP2 upregulation was only observed in C. sinensis seedlings subjected to osmotic and cold treatments. Taken together, these results indicate the existence of two highly conserved BiP genes in Citrus that are differentially regulated in the different tissues and in response to abiotic stresses.     

Characterization of genes associated with polyembryony and in vitro somatic embryogenesis in Citrus

Apomixis is a common reproduction system in the genus Citrus. To investigate the molecular mechanism of apomictic embryogenesis in Citrus, genes expressed specifically in an apomictic genotype were isolated by PCR-selected suppression subtractive hybridization with total RNAs obtained from the ovule at anthesis. Several genes showing conspicuously different expressions between polyembryonic (apomictic) and monoembryonic (nonapomictic) genotypes were selected, and their expression profiles during ovule development were analyzed in detail. This analysis identified two apomictic and three nonapomictic genotype-specific genes. Among the latter, msg-2 was highly expressed in the late stage of somatic embryogenesis. Specific expression during ovule development in monoembryonic cultivars and in the late stage of somatic embryogenesis indicated that msg-2 is not expressed in the initiation stage of polyembryogenesis and somatic embryogenesis, suggesting its role in suppressing initial cell formation of somatic embryos. The full-length complementary DNA of msg-2 contained small open reading frames in its sequence but showed no homology to functionally known genes in the public databases. As sequences similar to msg-2 were frequently found among Citrus expressed sequence tags, msg-2 may be associated with polyembryogenesis and somatic embryogenesis in a Citrus-specific manner.     

The plastome of Citrus. Physical map,variation among Citrus cultivars and species and comparison with related genera

Summary A physical plastome map was constructed for Citrus aurantium, and the plastomes of species and cultivars of Citrus and of two Citrus relatives were analysed by Southern blot-hybridisation of labelled total tobacco cpDNA to digests of total Citrus DNA. A resemblance was found between the plastomes of cultivars of C. limon (lemon), C. sinensis (orange), C. aurantium (sour orange), C. paradisii (grapefruit) and C. grandis (pomello). The plastomes of other Citrus types such as mandarin (C. reticulata) and citron (C. medico) differed from each other as well as from the plastomes of the aforementioned group. The plastomes of Poncirus trifoliata and Microcitrus sp. are distinct from each other as well as from the Citrus types.     

Large-scale analysis of NBS domain-encoding resistance gene analogs in Triticeae

Proteins containing nucleotide binding sites (NBS) encoded by plant resistance genes play an important role in the response of plants to a wide array of pathogens. In this paper, an in silico search was conducted in order to identify and characterize members of NBS-encoding gene family in the tribe of Triticeae. A final dataset of 199 sequences was obtained by four search methods. Motif analysis confirmed the general structural organization of the NBS domain in cereals, characterized by the presence of the six commonly conserved motifs: P-loop, RNBS-A, Kinase-2, Kinase-3a, RNBS-C and GLPL. We revealed the existence of 11 distinct distribution patterns of these motifs along the NBS domain. Four additional conserved motifs were shown to be significantly present in all 199 sequences. Phylogenetic analyses, based on genetic distance and parsimony, revealed a significant overlap between Triticeae sequences and Coiled coil-Nucleotide binding site-Leucine rich repeat (CNL)-type functional genes from monocotyledons. Furthermore, several Triticeae sequences belonged to clades containing functional homologs from non Triticeae species, which has allowed for these sequences to be functionally assigned. The findings reported, in this study, will provide a strong groundwork for the isolation of candidate R-genes in Triticeae crops and the understanding of their evolution.     

Diversity and evolution of a non-TIR-NBS sequence family that clusters to a chromosomal ’’hotspot” for disease resistance genes in soybean

In soybean, genes controlling resistance to numerous diseases have been shown to cluster to regions on several chromosomes. One such vital chromosomal region is on the soybean molecular linkage group (MLG) F flanked by the RFLP markers K644 and B212. Here, genes controlling resistance to bacterial blight, Phytophthora root rot, and several viral diseases, as well as QTLs conditioning resistance to corn earworm, root knot nematode, and white mold have been mapped. We have previously identified two classes (b and j) of disease resistance-related nucleotide binding site (NBS) sequences that localize to this cluster of genes. Using both cDNA and genomic analyses, we have studied one multi-gene family of sequences representing the previously reported class j NBS of soybean. This class of NBS resembles the RPS2-like NBS sequences. RPS2 and similar resistance genes are referred to as non-TIR because they do not encode motifs homologous to the Toll-Interleukin-1 region (TIR). By designing PCR primers that specifically target these non-TIR-NBS encoding sequences, we have amplified at least six class j sequence members from soybean. In addition, we have conducted genomic and cDNA library screenings to identify additional class j members. In all, we have characterized 12 class j NBS sequence members. These members have been mapped within a 2-cM region of the soybean F linkage group. We have also identified homoeologous chromosomal regions on linkage groups A2 and E that contain class j NBS sequences. A BLAST search of the GenBank database has shown that non-TIR NBS sequences are present across the legume family. We have compared these non-TIR sequences from other legumes with the soybean clones to assess the level of diversity within this class of disease resistance-related sequences.     

In vitro anther culture of sweet orange (Citrus sinensis L. Osbeck) genotypes and of a C. clementina × C. sinensis ‘Hamlin’ hybrid

Citrus, and particularly sweet oranges, are very recalcitrant to anther culture. In this paper it was evaluated for the first time the response of 27 genotypes of Citrus sinensis and of one hybrid C. clementina × C. sinensis, to in vitro anther culture. Ten genotypes of sweet oranges showed embryogenic callus induction, mostly blood sweet oranges genotypes, such as Tarocco, Moro and Sanguinelli. In vitro microspore developmental switches from the gamethophytic to the sporophytic pathway were shown by DAPI staining in microspores of these responsive genotypes, after 10 months in culture. However, microsatellite marker analyses showed that these calli were heterozygous. The flow-cytometric analysis of these embryogenic calli showed the presence of two peaks, corresponding to haploid (n) and diploid (2n) genotypes. Differently, anther cultures of the hybrid C. clementina × C. sinensis produced tri-haploid (3n) embryogenic calli and the embryos obtained were homozygous when analyzed by molecular markers (sample sequence repeats), confirming the more responsive characteristic of clementine to microspore embryogenesis through anther culture.     

Comparison of Five Populations of Tylenchulus semipenetrans to Citrus,Poncirus, and their Hybrids

The infectivity of five populations of Tylenchulus semipenetrans were compared and differentiated on 10 hosts (5 Citrus spp., 1 Poncirus trifoliata, and 4 hybrids of Citrus spp. X P. trifoliata). Differences in levels of infection and development (P = 0.01) occurred between Citrus spp. and P. trifoliata cv. ''Pomeroy'' and their three hybrids, C. paradisi X P. trifoliata cv. ''Swingle'' citruntelo and C. sinensis, cv. ''Ruby'' orange X P. trifoliata cv. ''Webber Fawcett 14-7'', and ''15-7''. Poncirus trifoliata cv. Pomeroy was susceptible to a California biotype 3 and highly resistant to the other citrus nematode populations. Low infection levels with California biotype 1, Arizona, and Florida populations on Swingle citrumelo, and the two Ruby orange hybrids indicated inherited resistance. Reproduction of the nematode population from Texas was greatest on the three hybrids, Swingle citrumelo, Ruby orange 14-7, and 15-7, from the California 1, Arizona, and Florida populations, but its comparable densities on P. trifoliata and Citrus spp. were not sufficiently different from these populations to consider it a separate biotype. California biotype 3 was sufficiently different from all other populations to be considered a different biotype, and it was named the "Poncirus biotype."