Genomics of the origin and evolution of Citrus.

The genus Citrus, comprising some of the most widely cultivated fruit crops worldwide, includes an uncertain number of species. Here we describe ten natural citrus species, using genomic, phylogenetic and biogeographic analyses of 60 accessions representing diverse citrus germ plasms, and propose that citrus diversified during the late Miocene epoch through a rapid southeast Asian radiation that correlates with a marked weakening of the monsoons. A second radiation enabled by migration across the Wallace line gave rise to the Australian limes in the early Pliocene epoch. Further identification and analyses of hybrids and admixed genomes provides insights into the genealogy of major commercial cultivars of citrus. Among mandarins and sweet orange, we find an extensive network of relatedness that illuminates the domestication of these groups. Widespread pummelo admixture among these mandarins and its correlation with fruit size and acidity suggests a plausible role of pummelo introgression in the selection of palatable mandarins. This work provides a new evolutionary framework for the genus Citrus.

Comparative transcriptomic analyses of citrus cold-resistant vs. sensitive rootstocks might suggest a relevant role of ABA signaling in triggering cold scion adaption.

BACKGROUND: The citrus genus comprises a number of sensitive tropical and subtropical species to cold stress, which limits global citrus distribution to certain latitudes and causes major economic loss. We used RNA-Seq technology to analyze changes in the transcriptome of Valencia delta seedless orange in response to long-term cold stress grafted on two frequently used citrus rootstocks: Carrizo citrange (CAR), considered one of the most cold-tolerant accessions; C. macrophylla (MAC), a very sensitive one. Our objectives were to identify the genetic mechanism that produce the tolerant or sensitive phenotypes in citrus, as well as to gain insights of the rootstock-scion interactions that induce the cold tolerance or sensitivity in the scion. RESULTS: Plants were kept at 1 ºC for 30 days. Samples were taken at 0, 15 and 30 days. The metabolomic analysis showed a significant increase in the concentration of free sugars and proline, which was higher for the CAR plants. Hormone quantification in roots showed a substantially increased ABA concentration during cold exposure in the CAR roots, which was not observed in MAC. Different approaches were followed to analyze gene expression. During the stress treatment, the 0-15-day comparison yielded the most DEGs. The functional characterization of DEGs showed enrichment in GO terms and KEGG pathways related to abiotic stress responses previously described in plant cold adaption. The DEGs analysis revealed that several key genes promoting cold adaption were up-regulated in the CAR plants, and those repressing it had higher expression levels in the MAC samples. CONCLUSIONS: The metabolomic and transcriptomic study herein performed indicates that the mechanisms activated in plants shortly after cold exposure remain active in the long term. Both the hormone quantification and differential expression analysis suggest that ABA signaling might play a relevant role in promoting the cold hardiness or sensitiveness of Valencia sweet orange grafted onto Carrizo citrange or Macrophylla rootstocks, respectively. Our work provides new insights into the mechanisms by which rootstocks modulate resistance to abiotic stress in the production variety grafted onto them.

Comparative transcriptomics of wild and commercial Citrus during early ripening reveals how domestication shaped fruit gene expression.

BACKGROUND: Interspecific hybridizations and admixtures were key in Citrus domestication, but very little is known about their impact at the transcriptomic level. To determine the effects of genome introgressions on gene expression, the transcriptomes of the pulp and flavedo of three pure species (citron, pure mandarin and pummelo) and four derived domesticated genetic admixtures (sour orange, sweet orange, lemon and domesticated mandarin) have been analyzed at color break. RESULTS: Many genes involved in relevant physiological processes for domestication, such sugar/acid metabolism and carotenoid/flavonoid synthesis, were differentially expressed among samples. In the low-sugar, highly acidic species lemon and citron, many genes involved in sugar metabolism, the TCA cycle and GABA shunt displayed a reduced expression, while the P-type ATPase CitPH5 and most subunits of the vacuolar ATPase were overexpressed. The red-colored species and admixtures were generally characterized by the overexpression in the flavedo of specific pivotal genes involved in the carotenoid biosynthesis, including phytoene synthase, ζ-carotene desaturase, ß-lycopene cyclase and CCD4b, a carotenoid cleavage dioxygenase. The expression patterns of many genes involved in flavonoid modifications, especially the flavonoid and phenylpropanoid O-methyltransferases showed extreme diversity. However, the most noticeable differential expression was shown by a chalcone synthase gene, which catalyzes a key step in the biosynthesis of flavonoids. This chalcone synthase was exclusively expressed in mandarins and their admixed species, which only expressed the mandarin allele. In addition, comparisons between wild and domesticated mandarins revealed that the major differences between their transcriptomes concentrate in the admixed regions. CONCLUSION: In this work we present a first study providing broad evidence that the genome introgressions that took place during citrus domestication largely shaped gene expression in their fruits.

IDA (INFLORESCENCE DEFICIENT IN ABSCISSION)-like peptides and HAE (HAESA)-like receptors regulate corolla abscission in Nicotiana benthamiana flowers.

BACKGROUND: Abscission is an active, organized, and highly coordinated cell separation process enabling the detachment of aerial organs through the modification of cell-to-cell adhesion and breakdown of cell walls at specific sites on the plant body known as abscission zones. In Arabidopsis thaliana, abscission of floral organs and cauline leaves is regulated by the interaction of the hormonal peptide INFLORESCENCE DEFICIENT IN ABSCISSION (IDA), a pair of redundant receptor-like protein kinases, HAESA (HAE) and HAESA-LIKE2 (HSL2), and SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) co-receptors. However, the functionality of this abscission signaling module has not yet been demonstrated in other plant species. RESULTS: The expression of the pair of NbenIDA1 homeologs and the receptor NbenHAE.1 was supressed at the base of the corolla tube by the inoculation of two virus-induced gene silencing (VIGS) constructs in Nicotiana benthamiana. These gene suppression events arrested corolla abscission but did not produce any obvious effect on plant growth. VIGS plants retained a higher number of corollas attached to the flowers than control plants, an observation related to a greater corolla breakstrength. The arrest of corolla abscission was associated with the preservation of the parenchyma tissue at the base of the corolla tube that, in contrast, was virtually collapsed in normal corollas. In contrast, the inoculation of a viral vector construct that increased the expression of NbenIDA1A at the base of the corolla tube negatively affected the growth of the inoculated plants accelerating the timing of both corolla senescence and abscission. However, the heterologous ectopic overexpression of citrus CitIDA3 and Arabidopsis AtIDA in N. benthamiana did not alter the standard plant phenotype suggesting that the proteolytic processing machinery was unable to yield active peptides. CONCLUSION: Here, we demonstrate that the pair of NbenIDA1 homeologs encoding small peptides of the IDA-like family and the receptor NbenHAE.1 control cellular breakdown at the base of the corolla tube awhere an adventitious AZ should be formed and, therefore, corolla abscission in N. benthamiana flowers. Altogether, our results provide the first evidence supporting the notion that the IDA-HAE/HSL2 signaling module is conserved in angiosperms.

Integration of mandarin (Citrus reticulata) cytogenetic map with its genome sequence.

Citrus is an extremely important genus in terms of world fruit production. Despite its economic importance and the small genome sizes of its species (2n = 18, 1C = 430 ± 68 Mbp), entire genomic assemblies have only recently become available for some of its representatives. Together with the previous CMA/DAPI banding and fluorescence in situ hybridization (FISH) in the group, these data are important for understanding the complex relationships between its species and for assisting breeding programs. To anchor genomic data with the cytogenetic map of mandarin (Citrus reticulata), the parental species of several economically important hybrids such as sweet orange and clementine, 18 BAC (bacterial artificial chromosome) clones were used. Eleven clementine BACs were positioned by BAC-FISH, doubling the number of chromosome markers so far available for BAC-FISH in citrus. Additionally, six previously mapped BACs were end-sequenced, allowing, together with one BAC previously sequenced, their assignment to scaffolds and the subsequent integration of chromosomes and the genome assembly. This study therefore established correlations between mandarin scaffolds and chromosomes, allowing further structural genomic and comparative study with the sweet orange genome, as well as insights into the chromosomal evolution of the group.

Transcriptomic analysis of Citrus clementina mandarin fruits maturation reveals a MADS-box transcription factor that might be involved in the regulation of earliness.

BACKGROUND: Harvest time is a relevant economic trait in citrus, and selection of cultivars with different fruit maturity periods has a remarkable impact in the market share. Generation of early- and late-maturing cultivars is an important target for citrus breeders, therefore, generation of knowledge regarding the genetic mechanisms controlling the ripening process and causing the early and late phenotypes is crucial. In this work we analyze the evolution of the transcriptome during fruit ripening in 3 sport mutations derived from the Fina clementine (Citrus clementina) mandarin: Clemenules (CLE), Arrufatina (ARR) and Hernandina (HER) that differ in their harvesting periods. CLE is considered a mid-season cultivar while ARR and HER are early- and late-ripening mutants, respectively. RESULTS: We used RNA-Seq technology to carry out a time course analysis of the transcriptome of the 3 mutations along the ripening period. The results indicated that in these mutants, earliness and lateness during fruit ripening correlated with the advancement or delay in the expression of a set of genes that may be implicated in the maturation process. A detailed analysis of the transcription factors known to be involved in the regulation of fruit ripening identified a member of the MADS box family whose expression was lower in ARR, the early-ripening mutant, and higher in HER, the late-ripening mutant. The pattern of expression of this gene during the maturation period was basically contrary to those of the ethylene biosynthetic genes, SAM and ACC synthases and ACC oxidase. The gene was present in hemizygous dose in the early-ripening mutant. CONCLUSIONS: Our analysis provides new clues about the genetic control of fruit ripening in citrus and allowed the identification of a transcription factor that could be involved in the early phenotype.

A Phylogenetic Analysis of 34 Chloroplast Genomes Elucidates the Relationships between Wild and Domestic Species within the Genus Citrus.

Citrus genus includes some of the most important cultivated fruit trees worldwide. Despite being extensively studied because of its commercial relevance, the origin of cultivated citrus species and the history of its domestication still remain an open question. Here, we present a phylogenetic analysis of the chloroplast genomes of 34 citrus genotypes which constitutes the most comprehensive and detailed study to date on the evolution and variability of the genus Citrus. A statistical model was used to estimate divergence times between the major citrus groups. Additionally, a complete map of the variability across the genome of different citrus species was produced, including single nucleotide variants, heteroplasmic positions, indels (insertions and deletions), and large structural variants. The distribution of all these variants provided further independent support to the phylogeny obtained. An unexpected finding was the high level of heteroplasmy found in several of the analyzed genomes. The use of the complete chloroplast DNA not only paves the way for a better understanding of the phylogenetic relationships within the Citrus genus but also provides original insights into other elusive evolutionary processes, such as chloroplast inheritance, heteroplasmy, and gene selection.

An RNA-Seq-based reference transcriptome for Citrus.

Previous RNA-Seq studies in citrus have been focused on physiological processes relevant to fruit quality and productivity of the major species, especially sweet orange. Less attention has been paid to vegetative or reproductive tissues, while most Citrus species have never been analysed. In this work, we characterized the transcriptome of vegetative and reproductive tissues from 12 Citrus species from all main phylogenetic groups. Our aims were to acquire a complete view of the citrus transcriptome landscape, to improve previous functional annotations and to obtain genetic markers associated with genes of agronomic interest. 28 samples were used for RNA-Seq analysis, obtained from 12 Citrus species: C. medica, C. aurantifolia, C. limon, C. bergamia, C. clementina, C. deliciosa, C. reshni, C. maxima, C. paradisi, C. aurantium, C. sinensis and Poncirus trifoliata. Four different organs were analysed: root, phloem, leaf and flower. A total of 3421 million Illumina reads were produced and mapped against the reference C. clementina genome sequence. Transcript discovery pipeline revealed 3326 new genes, the number of genes with alternative splicing was increased to 19,739, and a total of 73,797 transcripts were identified. Differential expression studies between the four tissues showed that gene expression is overall related to the physiological function of the specific organs above any other variable. Variants discovery analysis revealed the presence of indels and SNPs in genes associated with fruit quality and productivity. Pivotal pathways in citrus such as those of flavonoids, flavonols, ethylene and auxin were also analysed in detail.

Involvement of a citrus meiotic recombination TTC-repeat motif in the formation of gross deletions generated by ionizing radiation and MULE activation.

BACKGROUND: Transposable-element mediated chromosomal rearrangements require the involvement of two transposons and two double-strand breaks (DSB) located in close proximity. In radiobiology, DSB proximity is also a major factor contributing to rearrangements. However, the whole issue of DSB proximity remains virtually unexplored. RESULTS: Based on DNA sequencing analysis we show that the genomes of 2 derived mutations, Arrufatina (sport) and Nero (irradiation), share a similar 2 Mb deletion of chromosome 3. A 7 kb Mutator-like element found in Clemenules was present in Arrufatina in inverted orientation flanking the 5' end of the deletion. The Arrufatina Mule displayed "dissimilar" 9-bp target site duplications separated by 2 Mb. Fine-scale single nucleotide variant analyses of the deleted fragments identified a TTC-repeat sequence motif located in the center of the deletion responsible of a meiotic crossover detected in the citrus reference genome. CONCLUSIONS: Taken together, this information is compatible with the proposal that in both mutants, the TTC-repeat motif formed a triplex DNA structure generating a loop that brought in close proximity the originally distinct reactive ends. In Arrufatina, the loop brought the Mule ends nearby the 2 distinct insertion target sites and the inverted insertion of the transposable element between these target sites provoked the release of the in-between fragment. This proposal requires the involvement of a unique transposon and sheds light on the unresolved question of how two distinct sites become located in close proximity. These observations confer a crucial role to the TTC-repeats in fundamental plant processes as meiotic recombination and chromosomal rearrangements.

SNP mining in C. clementina BAC end sequences; transferability in the Citrus genus (Rutaceae), phylogenetic inferences and perspectives for genetic mapping.

BACKGROUND: With the increasing availability of EST databases and whole genome sequences, SNPs have become the most abundant and powerful polymorphic markers. However, SNP chip data generally suffers from ascertainment biases caused by the SNP discovery and selection process in which a small number of individuals are used as discovery panels. The ongoing International Citrus Genome Consortium sequencing project of the highly heterozygous Clementine and sweet orange genomes will soon result in the release of several hundred thousand SNPs. The primary goals of this study were: (i) to estimate the transferability within the genus Citrus of SNPs discovered from Clementine BACend sequencing (BES), (ii) to estimate bias associated with the very narrow discovery panel, and (iii) to evaluate the usefulness of the Clementine-derived SNP markers for diversity analysis and comparative mapping studies between the different cultivated Citrus species. RESULTS: Fifty-four accessions covering the main Citrus species and 52 interspecific hybrids between pummelo and Clementine were genotyped on a GoldenGate array platform using 1,457 SNPs mined from Clementine BES and 37 SNPs identified between and within C. maxima, C. medica, C. reticulata and C. micrantha. Consistent results were obtained from 622 SNP loci. Of these markers, 116 displayed incomplete transferability primarily in C. medica, C. maxima and wild Citrus species. The two primary biases associated with the SNP mining in Clementine were an overestimation of the C. reticulata diversity and an underestimation of the interspecific differentiation. However, the genetic stratification of the gene pool was high, with very frequent significant linkage disequilibrium. Furthermore, the shared intraspecific polymorphism and accession heterozygosity were generally enough to perform interspecific comparative genetic mapping. CONCLUSIONS: A set of 622 SNP markers providing consistent results was selected. Of the markers mined from Clementine, 80.5% were successfully transferred to the whole Citrus gene pool. Despite the ascertainment biases in relation to the Clementine origin, the SNP data confirm the important stratification of the gene pools around C. maxima, C. medica and C. reticulata as well as previous hypothesis on the origin of secondary species. The implemented SNP marker set will be very useful for comparative genetic mapping in Citrus and genetic association in C. reticulata.

A reference genetic map of C. clementina hort. ex Tan.; citrus evolution inferences from comparative mapping.

BACKGROUND: Most modern citrus cultivars have an interspecific origin. As a foundational step towards deciphering the interspecific genome structures, a reference whole genome sequence was produced by the International Citrus Genome Consortium from a haploid derived from Clementine mandarin. The availability of a saturated genetic map of Clementine was identified as an essential prerequisite to assist the whole genome sequence assembly. Clementine is believed to be a 'Mediterranean' mandarin × sweet orange hybrid, and sweet orange likely arose from interspecific hybridizations between mandarin and pummelo gene pools. The primary goals of the present study were to establish a Clementine reference map using codominant markers, and to perform comparative mapping of pummelo, sweet orange, and Clementine. RESULTS: Five parental genetic maps were established from three segregating populations, which were genotyped with Single Nucleotide Polymorphism (SNP), Simple Sequence Repeats (SSR) and Insertion-Deletion (Indel) markers. An initial medium density reference map (961 markers for 1084.1 cM) of the Clementine was established by combining male and female Clementine segregation data. This Clementine map was compared with two pummelo maps and a sweet orange map. The linear order of markers was highly conserved in the different species. However, significant differences in map size were observed, which suggests a variation in the recombination rates. Skewed segregations were much higher in the male than female Clementine mapping data. The mapping data confirmed that Clementine arose from hybridization between 'Mediterranean' mandarin and sweet orange. The results identified nine recombination break points for the sweet orange gamete that contributed to the Clementine genome. CONCLUSIONS: A reference genetic map of citrus, used to facilitate the chromosome assembly of the first citrus reference genome sequence, was established. The high conservation of marker order observed at the interspecific level should allow reasonable inferences of most citrus genome sequences by mapping next-generation sequencing (NGS) data in the reference genome sequence. The genome of the haploid Clementine used to establish the citrus reference genome sequence appears to have been inherited primarily from the 'Mediterranean' mandarin. The high frequency of skewed allelic segregations in the male Clementine data underline the probable extent of deviation from Mendelian segregation for characters controlled by heterozygous loci in male parents.

Development of indel markers from Citrus clementina (Rutaceae) BAC-end sequences and interspecific transferability in Citrus.

PREMISE OF THE STUDY: Indel markers were developed from BAC-end sequences of Citrus clementina cv. Nules. Transferability and polymorphism were tested in the Citrus genus to estimate the potential of indel markers mined from a single genotype for use in genetic studies. METHODS AND RESULTS: Using polyacrylamide gel electrophoresis and DNA silver staining, 89 indel markers were tested for their transferability and polymorphism. Thirty-eight markers were selected. Heterozygosity in C. clementina cv. Nules was confirmed for 33 of these indel pairs. A preliminary diversity study using a capillary electrophoresis fragment analyzer was conducted with 21 indels using 45 accessions representing Citrus genus diversity. Intraspecific and interspecific polymorphisms were observed. CONCLUSIONS: These results indicate the utility of indel markers developed from sequence data of a single genotype of interspecific origin. In Citrus, these markers will be useful for genetic mapping, germplasm characterization, and phylogenetic assignment of DNA fragments.

Shaping the biology of citrus: I. Genomic determinants of evolution.

We performed genomic analyses on wild species of the genus Citrus to identify major determinants of evolution. The most notable effect occurred on the pathogen-defense genes, as observed in many other plant genera. The gene space was also characterized by changes in gene families intimately related to relevant biochemical properties of citrus fruit, such as pectin modifying enzymes, HDR (4-hydroxy-3-methylbut-2-enyl diphosphate reductase) genes, and O-methyltransferases. Citrus fruits are highly abundant on pectins and secondary metabolites such as terpenoids and flavonoids, the targets of these families. Other gene types under positive selection, expanded through tandem duplications and retained as triplets from whole genome duplications, codified for purple acid phosphatases and MATE-efflux proteins. Although speciation has not been especially rapid in the genus, analyses of selective pressure at the codon level revealed that the extant species evolved from the ancestral citrus radiation show signatures of pervasive adaptive evolution and is therefore potentially responsible for the vast phenotypic differences observed among current species.

Shaping the biology of citrus: II. Genomic determinants of domestication.

We performed genomic analyses on species and varieties of the genus Citrus to identify several determinants of domestication, based on the pattern of pummelo [Citrus maxima (Burr. f) Merr] and mandarin (Citrus reticulata Blanco) admixture into the ancestral genome, as well as population genetic tests at smaller scales. Domestication impacted gene families regulating pivotal components of citrus flavor (such as acidity) because in edible mandarin varieties, chromosome areas with negative Tajimas values were enriched with genes associated with the regulation of citric acid. Detection of sweeps in edible mandarins that diverged from wild relatives indicated that domestication reduced chemical defenses involving cyanogenesis and alkaloid synthesis, thus increasing palatability. Also, a cluster of SAUR genes in domesticated mandarins derived from the pummelo genome appears to contain candidate genes controlling fruit size. Similarly, conserved stretches of pure mandarin areas were likely important as well for domestication, as, for example, a fragment in chromosome 1 that is involved in the apomictic reproduction of most edible mandarins. Interestingly, our results also support the hypothesis that various genes subject to selective pressure during evolution or derived from whole genome duplication events later became potential targets of domestication.

The aconitate hydratase family from Citrus.

BACKGROUND: Research on citrus fruit ripening has received considerable attention because of the importance of citrus fruits for the human diet. Organic acids are among the main determinants of taste and organoleptic quality of fruits and hence the control of fruit acidity loss has a strong economical relevance. In citrus, organic acids accumulate in the juice sac cells of developing fruits and are catabolized thereafter during ripening. Aconitase, that transforms citrate to isocitrate, is the first step of citric acid catabolism and a major component of the citrate utilization machinery. In this work, the citrus aconitase gene family was first characterized and a phylogenetic analysis was then carried out in order to understand the evolutionary history of this family in plants. Gene expression analyses of the citrus aconitase family were subsequently performed in several acidic and acidless genotypes to elucidate their involvement in acid homeostasis. RESULTS: Analysis of 460,000 citrus ESTs, followed by sequencing of complete cDNA clones, identified in citrus 3 transcription units coding for putatively active aconitate hydratase proteins, named as CcAco1, CcAco2 and CcAco3. A phylogenetic study carried on the Aco family in 14 plant species, shows the presence of 5 Aco subfamilies, and that the ancestor of monocot and dicot species shared at least one Aco gene. Real-time RT-PCR expression analyses of the three aconitase citrus genes were performed in pulp tissues along fruit development in acidic and acidless citrus varieties such as mandarins, oranges and lemons. While CcAco3 expression was always low, CcAco1 and CcAco2 genes were generally induced during the rapid phase of fruit growth along with the maximum in acidity and the beginning of the acid reduction. Two exceptions to this general pattern were found: 1) Clemenules mandarin failed inducing CcAco2 although acid levels were rapidly reduced; and 2) the acidless "Sucreña" orange showed unusually high levels of expression of both aconitases, an observation correlating with the acidless phenotype. However, in the acidless "Dulce" lemon aconitase expression was normal suggesting that the acidless trait in this variety is not dependent upon aconitases. CONCLUSIONS: Phylogenetic studies showed the occurrence of five different subfamilies of aconitate hydratase in plants and sequence analyses identified three active genes in citrus. The pattern of expression of two of these genes, CcAco1 and CcAco2, was normally associated with the timing of acid content reduction in most genotypes. Two exceptions to this general observation suggest the occurrence of additional regulatory steps of citrate homeostasis in citrus.

Identification of genes associated with bud dormancy release in Prunus persica by suppression subtractive hybridization.

To better understand the molecular and physiological mechanisms underlying maintenance and release of seasonal bud dormancy in perennial trees, we identified differentially expressed genes during dormancy progression in reproductive buds from peach (Prunus persica [L.] Batsch) by suppression subtractive hybridization (SSH) and microarray hybridization. Four SSH libraries were constructed, which were respectively enriched in cDNA highly expressed in dormant buds (named DR), in dormancy-released buds (RD) and in the cultivars with different chilling requirement, 'Zincal 5' (ZS) and 'Springlady' (SZ), sampled after dormancy release. About 2500 clones picked from the four libraries were loaded on a glass microarray. Hybridization of microarrays with the final products of SSH procedure was performed in order to validate the selected clones that were effectively enriched in their respective sample. Nearly 400 positive clones were sequenced, which corresponded to 101 different unigenes with diverse functional annotation. We obtained DAM4, 5 and 6 genes coding for MADS-box transcription factors previously related to growth cessation and terminal bud formation in the evergrowing mutant of peach. Several other cDNAs are similar to dormancy factors described in other species, and others have been related to bud dormancy for the first time in this study. Quantitative reverse transcription polymerase chain reaction analysis confirmed differential expression of cDNAs coding for a Zn-finger transcription factor, a GRAS-like regulator, a DNA-binding protein and proteins similar to forisome subunits involved in the reversible occlusion of sieve elements in Fabaceae, among others.

High-throughput functional annotation and data mining with the Blast2GO suite.

Functional genomics technologies have been widely adopted in the biological research of both model and non-model species. An efficient functional annotation of DNA or protein sequences is a major requirement for the successful application of these approaches as functional information on gene products is often the key to the interpretation of experimental results. Therefore, there is an increasing need for bioinformatics resources which are able to cope with large amount of sequence data, produce valuable annotation results and are easily accessible to laboratories where functional genomics projects are being undertaken. We present the Blast2GO suite as an integrated and biologist-oriented solution for the high-throughput and automatic functional annotation of DNA or protein sequences based on the Gene Ontology vocabulary. The most outstanding Blast2GO features are: (i) the combination of various annotation strategies and tools controlling type and intensity of annotation, (ii) the numerous graphical features such as the interactive GO-graph visualization for gene-set function profiling or descriptive charts, (iii) the general sequence management features and (iv) high-throughput capabilities. We used the Blast2GO framework to carry out a detailed analysis of annotation behaviour through homology transfer and its impact in functional genomics research. Our aim is to offer biologists useful information to take into account when addressing the task of functionally characterizing their sequence data.

Reprogramming of Retrotransposon Activity during Speciation of the Genus Citrus.

Speciation of the genus Citrus from a common ancestor has recently been established to begin ∼8 Ma during the late Miocene, a period of major climatic alterations. Here, we report the changes in activity of Citrus LTR retrotransposons during the process of diversification that gave rise to the current Citrus species. To reach this goal, we analyzed four pure species that diverged early during Citrus speciation, three recent admixtures derived from those species and an outgroup of the Citrus clade. More than 30,000 retrotransposons were grouped in ten linages. Estimations of LTR insertion times revealed that retrotransposon activity followed a species-specific pattern of change that could be ascribed to one of three different models. In some genomes, the expected pattern of gradual transposon accumulation was suddenly arrested during the radiation of the ancestor that gave birth to the current Citrus species. The individualized analyses of retrotransposon lineages showed that in each and every species studied, not all lineages follow the general pattern of the species itself. For instance, in most of the genomes, the retrotransposon activity of elements from the SIRE lineage reached its highest level just before Citrus speciation, while for Retrofit elements, it has been steadily growing. Based on these observations, we propose that Citrus retrotransposons may respond to stressful conditions driving speciation as a part of the genetic response involved in adaptation. This proposal implies that the evolving conditions of each species interact with the internal regulatory mechanisms of the genome controlling the proliferation of mobile elements.

Development of genomic resources for Citrus clementina: characterization of three deep-coverage BAC libraries and analysis of 46,000 BAC end sequences.

BACKGROUND: Citrus species constitute one of the major tree fruit crops of the subtropical regions with great economic importance. However, their peculiar reproductive characteristics, low genetic diversity and the long-term nature of tree breeding mostly impair citrus variety improvement. In woody plants, genomic science holds promise of improvements and in the Citrus genera the development of genomic tools may be crucial for further crop improvements. In this work we report the characterization of three BAC libraries from Clementine (Citrus clementina), one of the most relevant citrus fresh fruit market cultivars, and the analyses of 46.000 BAC end sequences. Clementine is a diploid plant with an estimated haploid genome size of 367 Mb and 2n = 18 chromosomes, which makes feasible the use of genomics tools to boost genetic improvement. RESULTS: Three genomic BAC libraries of Citrus clementina were constructed through EcoRI, MboI and HindIII digestions and 56,000 clones, representing an estimated genomic coverage of 19.5 haploid genome-equivalents, were picked. BAC end sequencing (BES) of 28,000 clones produced 28.1 Mb of genomic sequence that allowed the identification of the repetitive fraction (12.5% of the genome) and estimation of gene content (31,000 genes) of this species. BES analyses identified 3,800 SSRs and 6,617 putative SNPs. Comparative genomic studies showed that citrus gene homology and microsyntheny with Populus trichocarpa was rather higher than with Arabidopsis thaliana, a species phylogenetically closer to citrus. CONCLUSION: In this work, we report the characterization of three BAC libraries from C. clementina, and a new set of genomic resources that may be useful for isolation of genes underlying economically important traits, physical mapping and eventually crop improvement in Citrus species. In addition, BAC end sequencing has provided a first insight on the basic structure and organization of the citrus genome and has yielded valuable molecular markers for genetic mapping and cloning of genes of agricultural interest. Paired end sequences also may be very helpful for whole-genome sequencing programs.

Transferability of the EST-SSRs developed on Nules clementine (Citrus clementina Hort ex Tan) to other Citrus species and their effectiveness for genetic mapping.

BACKGROUND: During the last decade, numerous microsatellite markers were developed for genotyping and to identify closely related plant genotypes. In citrus, previously developed microsatellite markers were arisen from genomic libraries and more often located in non coding DNA sequences. To optimize the use of these EST-SSRs as genetic markers in genome mapping programs and citrus systematic analysis, we have investigated their polymorphism related to the type (di or trinucleotide) or their position in the coding sequences. RESULTS: Among 11000 unigenes from a Clementine EST library, we have found at least one microsatellite sequence (repeated units size ranged from 2 to 6 nucleotides) in 1500 unigenes (13.6%). More than 95% of these SSRs were di or trinucleotides. If trinucleotide microsatellites were encountered trough all part of EST sequences, dinucleotide microsatellites were preferentially (50%) concentrated in the 5' 100th nucleotides. We assessed the polymorphism of 41 EST-SSR, by PCR amplification droved with flanking primers among ten Citrus species plus 3 from other genera. More than 90% of EST-SSR markers were polymorphic. Furthermore, dinucleotide microsatellite markers were more polymorphic than trinucleotide ones, probably related to their distribution that was more often located in the 5' UnTranslated Region (UTR). We obtained a good agreement of diversity relationships between the citrus species and relatives assessed with EST-SSR markers with the established taxonomy and phylogeny. To end, the heterozygosity of each genotype and all dual combinations were studied to evaluate the percentage of mappable markers. Higher values (> 45%) were observed for putative Citrus inter-specific hybrids (lime lemon, or sour orange) than for Citrus basic true species (mandarin, pummelo and citron) (<30%). Most favorable combinations for genome mapping were observed in those involving interspecific hybrid genotypes. Those gave higher levels of mappable markers (>70%) with a significant proportion suitable for synteny analysis. CONCLUSION: Fourty one new EST-SSR markers were produced and were available for citrus genetic studies. Whatever the position of the SSR in the ESTs the EST-SSR markers we developed are powerful to investigate genetic diversity and genome mapping in citrus.