An RNA-Seq-based reference transcriptome for Citrus

Citrus genomes: past, present and future

Over the past decade, genome sequencing and assembly approaches have been greatly improved, resulting in the assembly of many genomes for citrus, including wild, domesticated, and citrus-related genomes. Improvements in technologies have led to assembled genomes with higher completeness, contiguity, quality, and accuracy that have greatly facilitated annotation and analysis. This review summarizes the evolution of the sequencing, assembly, and annotation technologies leading to citrus genomes over the past 11 years, a comprehensive evaluation of their quality, contiguity, and completeness, and the major findings and applications. Of the 50 genomes now available, 35 have been assembled to chromosome level and 15 to draft level, and 14 were haplotype-resolved assemblies. To date there have been four pangenome-wide studies for citrus. The very recent genomes assembled with long-read sequencing have achieved >99% and >98% assembly and annotation completeness (BUSCO), respectively. However, some early genomes are not of the same high quality as more recently sequenced genomes and would benefit from re-sequencing. A more comprehensive pangenome based upon a larger set of species and genotypes assembled at the haplotype level would allow genomics to deliver the maximum benefits for citrus improvement and research.

OrangeExpDB: an integrative gene expression database for Citrus spp

BACKGROUND

Citrus is a major fruit crop, and RNA-sequencing (RNA-seq) data can be utilized to investigate its gene functions, heredity, evolution, development, and the detection of genes linked to essential traits or resistance to pathogens. However, it is challenging to use the public RNA-seq datasets for researchers without bioinformatics training, and expertise.

RESULTS

OrangeExpDB is a web-based database that integrates transcriptome data of various Citrus spp., including C. limon (L.) Burm., C. maxima (Burm.) Merr., C. reticulata Blanco, C. sinensis (L.) Osbeck, and Poncirus trifoliata (L.) Raf., downloaded from the NCBI SRA database. It features a blast tool for browsing and searching, enabling quick download of expression matrices for different transcriptome samples. Expression of genes of interest can be easily generated by searching gene IDs or sequence similarity. Expression data in text format can be downloaded and presented as a heatmap, with additional sample information provided at the bottom of the webpage.

CONCLUSIONS

Researchers can utilize OrangeExpDB to facilitate functional genomic analysis and identify key candidate genes, leveraging publicly available citrus RNA-seq datasets. OrangeExpDB can be accessed at http://www.orangeexpdb.com/ .

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Current insights into posttranscriptional regulation of fleshy fruit ripening

Fruit ripening is a complicated process that is accompanied by the formation of fruit quality. It is not only regulated at the transcriptional level via transcription factors or DNA methylation but also fine-tuned after transcription occurs. Here, we review recent advances in our understanding of key regulatory mechanisms of fleshy fruit ripening after transcription. We mainly highlight the typical mechanisms by which fruit ripening is controlled, namely, alternative splicing, mRNA N6-methyladenosine RNA modification methylation, and noncoding RNAs at the posttranscriptional level; regulation of translation efficiency and upstream open reading frame-mediated translational repression at the translational level; and histone modifications, protein phosphorylation, and protein ubiquitination at the posttranslational level. Taken together, these posttranscriptional regulatory mechanisms, along with transcriptional regulation, constitute the molecular framework of fruit ripening. We also critically discuss the potential usage of some mechanisms to improve fruit traits.

De novo transcriptome analysis and identification of defensive genes in garlic (Allium sativum L.) using high-throughput sequencing

BACKGROUND

Garlic (Allium sativum L.) is the second most widely cultivated Allium which is mainly grown in temperate regions and used as a flavoring agent in a wide variety of foods. Garlic contains various bioactive compounds whose metabolic pathways, plant-pathogen interactions, defensive genes, identify interaction networks, and functional genomics were not previously predicted in the garlic at the genomic level. To address this issue, we constructed two garlic Illumina 2000 libraries from tissues of garlic clove and leaf.

RESULTS

Approximately 43 million 125 bp paired-end reads were obtained in the two libraries. A total of 239,973 contigs were generated by de novo assembly of both samples and were compared with the sequences in the NCBI non-redundant protein database (Nr). In total, 42% of contigs were matched to known proteins in public databases including Nr, Gene Ontology (GO), and Cluster Orthologous Gene Database (COG), and then, contigs were mapped to 138 via functional annotation against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, a number of regulatory genes including the CCHC (Zn) family, followed by WD40, bromodomain, bZIP, AP2-EREBP, BED-type (Zn) proteins, and defense response proteins related to different conserved domains, such as RGA3, NBS-LRR, TIR-NBS-LRR, LRR, NBS-ARC, and CC-NBS-LRR were discovered based on the transcriptome dataset. We compared the ortholog gene family of the A. sativum transcriptome to A. thaliana, O. sativa, and Z. mays and found that 12,077 orthologous gene families are specific to A. sativum L. Furthermore, we identified genes involved in plant defense mechanisms, their protein-protein interaction network, and plant-pathogen interaction pathways.

CONCLUSIONS

Our study contains an extensive sequencing and functional gene-annotation analysis of A. sativum L. The findings provide insights into the molecular basis of TFs, defensive genes, and a reference for future studies on the genetics and breeding of A. sativum L.

Conserving Citrus Diversity: From Vavilov's Early Explorations to Genebanks around the World

Citrus is among the most economically important fruit crops. Its vast species diversity and global production was observed by N.I. Vavilov during his international plant explorations from the early to mid-1900s. Currently, ex situ citrus collections located around the world conserve and protect citrus genetic resources, as revealed in a survey conducted in 2021. Responses were received from 43 collections in 27 countries, of which 35 provided data regarding collection composition, management practices, and security, as well as other information. The six largest citrus collections have between 1000 and 1735 accessions. The largest accession holdings are mandarins and sweet oranges, although all citrus fruit types are maintained: mandarin, sweet orange, lemon, pummelo, grapefruit, hybrids, lime, sour orange, citron, kumquat, papeda, finger lime, and crop wild relatives. Diseases pose significant threats to collections, though some collections are maintained in a clean-plant state as a result of intensive sanitation efforts. National and regional quarantine regulations often limit the export and import of citrus plants or propagative materials, thus limiting the availability of materials at an international level. Resources, both financial and human, are necessary to ensure the long-term safety and security of citrus collections on a global scale. Future efforts to develop citrus genebanking communities will provide opportunities for improved conservation, as well as collaborations and training.

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.

CitrusKB: a comprehensive knowledge base for transcriptome and interactome of Citrus spp. infected by Xanthomonas citri subsp. citri at different infection stages

Citrus canker type A is a serious disease caused by Xanthomonas citri subsp. citri (X. citri), which is responsible for severe losses to growers and to the citrus industry worldwide. To date, no canker-resistant citrus genotypes are available, and there is limited information regarding the molecular and genetic mechanisms involved in the early stages of the citrus canker development. Here, we present the CitrusKB knowledge base. This is the first in vivo interactome database for different citrus cultivars, and it was produced to provide a valuable resource of information on citrus and their interaction with the citrus canker bacterium X. citri. CitrusKB provides tools for a user-friendly web interface to let users search and analyse a large amount of information regarding eight citrus cultivars with distinct levels of susceptibility to the disease, with controls and infected plants at different stages of infection by the citrus canker bacterium X. citri. Currently, CitrusKB comprises a reference citrus genome and its transcriptome, expressed transcripts, pseudogenes and predicted genomic variations (SNPs and SSRs). The updating process will continue over time by the incorporation of novel annotations and analysis tools. We expect that CitrusKB may substantially contribute to the field of citrus genomics. CitrusKB is accessible at http://bioinfo.deinfo.uepg.br/citrus. Users can download all the generated raw sequences and generated datasets by this study from the CitrusKB website.

Genome sequencing and CRISPR/Cas9 gene editing of an early flowering Mini-Citrus (Fortunella hindsii)

Hongkong kumquat (Fortunella hindsii) is a wild citrus species characterized by dwarf plant height and early flowering. Here, we identified the monoembryonic F. hindsii (designated as 'Mini-Citrus') for the first time and constructed its selfing lines. This germplasm constitutes an ideal model for the genetic and functional genomics studies of citrus, which have been severely hindered by the long juvenility and inherent apomixes of citrus. F. hindsii showed a very short juvenile period (~8 months) and stable monoembryonic phenotype under cultivation. We report the first de novo assembled 373.6 Mb genome sequences (Contig-N50 2.2 Mb and Scaffold-N50 5.2 Mb) for F. hindsii. In total, 32 257 protein-coding genes were annotated, 96.9% of which had homologues in other eight Citrinae species. The phylogenomic analysis revealed a close relationship of F. hindsii with cultivated citrus varieties, especially with mandarin. Furthermore, the CRISPR/Cas9 system was demonstrated to be an efficient strategy to generate target mutagenesis on F. hindsii. The modifications of target genes in the CRISPR-modified F. hindsii were predominantly 1-bp insertions or small deletions. This genetic transformation system based on F. hindsii could shorten the whole process from explant to T1 mutant to about 15 months. Overall, due to its short juvenility, monoembryony, close genetic background to cultivated citrus and applicability of CRISPR, F. hindsii shows unprecedented potentials to be used as a model species for citrus research.

Facing Climate Change: Biotechnology of Iconic Mediterranean Woody Crops

The Mediterranean basin is especially sensitive to the adverse outcomes of climate change and especially to variations in rainfall patterns and the incidence of extremely high temperatures. These two concurring adverse environmental conditions will surely have a detrimental effect on crop performance and productivity that will be particularly severe on woody crops such as citrus, olive and grapevine that define the backbone of traditional Mediterranean agriculture. These woody species have been traditionally selected for traits such as improved fruit yield and quality or alteration in harvesting periods, leaving out traits related to plant field performance. This is currently a crucial aspect due to the progressive and imminent effects of global climate change. Although complete genome sequence exists for sweet orange (Citrus sinensis) and clementine (Citrus clementina), olive tree (Olea europaea) and grapevine (Vitis vinifera), the development of biotechnological tools to improve stress tolerance still relies on the study of the available genetic resources including interspecific hybrids, naturally occurring (or induced) polyploids and wild relatives under field conditions. To this respect, post-genomic era studies including transcriptomics, metabolomics and proteomics provide a wide and unbiased view of plant physiology and biochemistry under adverse environmental conditions that, along with high-throughput phenotyping, could contribute to the characterization of plant genotypes exhibiting physiological and/or genetic traits that are correlated to abiotic stress tolerance. The ultimate goal of precision agriculture is to improve crop productivity, in terms of yield and quality, making a sustainable use of land and water resources under adverse environmental conditions using all available biotechnological tools and high-throughput phenotyping. This review focuses on the current state-of-the-art of biotechnological tools such as high throughput -omics and phenotyping on grapevine, citrus and olive and their contribution to plant breeding programs.

Transcriptome analysis of the endangered Notopterygium incisum: Cold-tolerance gene discovery and identification of EST-SSR and SNP markers

Notopterygium incisum C. C. Ting ex H. T. Chang (Apiaceae) is an endangered perennial herb in China. The lack of transcriptomic and genomic resources for N. incisum greatly hinders studies of its population genetics and conservation. In this study, we employed RNA-seq technology to characterize transcriptomes for the flowers, leaves, and stems of this endangered herb. A total of 56 million clean reads were assembled into 120,716 unigenes with an N50 length of 850 bp. Among these unigenes, 70,245 (58.19%) were successfully annotated and 65,965 (54.64%) were identified as coding sequences based on their similarities with sequences in public databases. We identified 21 unigenes that had significant relationships with cold tolerance in N. incisum according to gene ontology (GO) annotation analysis. In addition, 13,149 simple sequence repeats (SSRs) and 85,681 single nucleotide polymorphisms were detected as potential molecular genetic markers. Ninety-six primer pairs of SSRs were randomly selected to validate their amplification efficiency and polymorphism. Nineteen SSR loci exhibited polymorphism in three natural populations of N. incisum. These results provide valuable resources to facilitate future functional genomics and conservation genetics studies of N. incisum.

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.

Comparative transcriptome analysis of Poncirus trifoliata identifies a core set of genes involved in arbuscular mycorrhizal symbiosis

The perennial woody plants of citrus are one of the most important fruit crops in the world and largely depends on arbuscular mycorrhizal symbiosis (AMS) to obtain essential nutrients from soil. However, the molecular aspects of AMS in citrus and perennial woody plants in general have largely been understudied. We used RNA-sequencing to identify differentially expressed genes in roots of Poncirus trifoliata upon mycorrhization by the AM fungus Glomus versiforme and evaluated their conservation by comparative transcriptome analyses with four herbaceous model plants. We identified 282 differentially expressed genes in P. trifoliata, including orthologs of 21 genes with characterized roles in AMS and 83 genes that are considered to be conserved in AM-host plants. Comparative transcriptome analysis revealed a 'core set' of 156 genes from P. trifoliata whose orthologous genes from at least three of the five species also exhibited similar transcriptional changes during AMS. Functional analysis of one of these conserved AM-induced genes, a 3-keto-acyl-ACP reductase (FatG) involved in fatty acid biosynthesis, confirmed its involvement in AMS in Medicago truncatula. Our results identify a core transcriptional program for AMS that is largely conserved between P. trifoliata and other plants. The comparative transcriptomics approach adds to previous phylogenomics studies to identify conserved genes required for AMS.

Transcriptome Analysis of Bael (Aegle marmelos (L.) Corr.) a Member of Family Rutaceae

Aegle marmelos (L.) Corr. is a medicinally and horticulturally important tree member of the family Rutaceae. It is native to India, where it is also known as Bael. Despite its importance, the genomic resources of this plant are scarce. This study presented the first-ever report of expressed transcripts in the leaves of Aegle marmelos. A total of 133,616 contigs were assembled to 46,335 unigenes with minimum and maximum lengths of 201 bp and 14,853 bp, respectively. There were 7002 transcription factors and 94,479 simple sequence repeat (SSR) markers. The A. marmelos transcripts were also annotated based on information from other members of Rutaceae; namely Citrus clementina and Citrus sinensis. A total of 482 transcripts were annotated as cytochrome p450s (CYPs), and 314 transcripts were annotated as glucosyltransferases (GTs). In the A. marmelos leaves, the monoterpenoid biosynthesis pathway was predominant. This study provides an important genomic resource along with useful information about A. marmelos.

The HTPmod Shiny application enables modeling and visualization of large-scale biological data

The wave of high-throughput technologies in genomics and phenomics are enabling data to be generated on an unprecedented scale and at a reasonable cost. Exploring the large-scale data sets generated by these technologies to derive biological insights requires efficient bioinformatic tools. Here we introduce an interactive, open-source web application (HTPmod) for high-throughput biological data modeling and visualization. HTPmod is implemented with the Shiny framework by integrating the computational power and professional visualization of R and including various machine-learning approaches. We demonstrate that HTPmod can be used for modeling and visualizing large-scale, high-dimensional data sets (such as multiple omics data) under a broad context. By reinvestigating example data sets from recent studies, we find not only that HTPmod can reproduce results from the original studies in a straightforward fashion and within a reasonable time, but also that novel insights may be gained from fast reinvestigation of existing data by HTPmod.

Dijun Chen et al. present HTPmod, a Shiny web application for modeling and visualization of large-scale genomic and phenomic datasets. The authors show that HTPmod can quickly reproduce analyses of high-throughput biological datasets and produce publication-quality figures.