Haplotype analysis of sucrose synthase gene family in three Saccharum species

Transcriptomic Analysis for Diurnal Temperature Differences Reveals Gene-Regulation-Network Response to Accumulation of Bioactive Ingredients of Protocorm-like Bodies in Dendrobium officinale

Dendrobium officinale Kimura et Migo (D. officinale) is one of the most important traditional Chinese medicinal herbs, celebrated for its abundant bioactive ingredients. This study demonstrated that the diurnal temperature difference (DIF) (T1: 13/13 °C, T2: 25/13 °C, and T3: 25/25 °C) was more favorable for high chlorophyll, increased polysaccharide, and total flavonoid contents compared to constant temperature treatments in D. officinale PLBs. The transcriptome analysis revealed 4251, 4404, and 4536 differentially expressed genes (DEGs) in three different comparisons (A: 25/13 °C vs. 13/13 °C, B: 13/13 °C vs. 25/25 °C, and C: 25/13 °C vs. 25/25 °C, respectively). The corresponding up-/down-regulated DEGs were 1562/2689, 2825/1579, and 2310/2226, respectively. GO and KEGG enrichment analyses of DEGs showed that the pathways of biosynthesis of secondary metabolites, carotenoid biosynthesis, and flavonoid biosynthesis were enriched in the top 20; further analysis of the sugar- and flavonol-metabolism pathways in D. officinale PLBs revealed that the DIF led to a differential gene expression in the enzymes linked to sugar metabolism, as well as to flavonol metabolism. Certain key metabolic genes related to ingredient accumulation were identified, including those involved in polysaccharide metabolism (SUS, SUT, HKL1, HGL, AMY1, and SS3) and flavonol (UGT73C and UGT73D) metabolism. Therefore, these findings indicated that these genes may play an important role in the regulatory network of the DIF in the functional metabolites of D. officinale PLBs. In a MapMan annotation of abiotic stress pathways, the DEGs with significant changes in their expression levels were mainly concentrated in the heat-stress pathways, including heat-shock proteins (HSPs) and heat-shock transcription factors (HSFs). In particular, the expression levels of HSP18.2, HSP70, and HSF1 were significantly increased under DIF treatment, which suggested that HSF1, HSP70 and HSP18.2 may respond to the DIF. In addition, they can be used as candidate genes to study the effect of the DIF on the PLBs of D. officinale. The results of our qPCR analysis are consistent with those of the transcriptome-expression analysis, indicating the reliability of the sequencing. The results of this study revealed the transcriptome mechanism of the DIF on the accumulation of the functional metabolic components of D. officinale. Furthermore, they also provide an important theoretical basis for improving the quality of D. officinale via the DIF in production.

Current perspectives on the regulatory mechanisms of sucrose accumulation in sugarcane

Sugars transported from leaves (source) to stems (sink) energize cell growth, elongation, and maintenance. which are regulated by a variety of genes. This review reflects progress and prospects in the regulatory mechanism for maximum sucrose accumulation, including the role of sucrose metabolizing enzymes, sugar transporters and the elucidation of post-transcriptional control of sucrose-induced regulation of translation (SIRT) in the accumulation of sucrose. The current review suggests that SIRT is emerging as a significant mechanism controlling Scbzip44 activities in response to endogenous sugar signals (via the negative feedback mechanism). Sucrose-controlled upstream open reading frame (SC-uORF) exists at the 5' leader region of Scbzip44's main ORF, which inhibits sucrose accumulation through post-transcriptional regulatory mechanisms. Sucrose transporters (SWEET1a/4a/4b/13c, TST, SUT1, SUT4 and SUT5) are crucial for sucrose translocation from source to sink. Particularly, SWEET13c was found to be a major contributor to the efflux in the transportation of stems. Tonoplast sugar transporters (TSTs), which import sucrose into the vacuole, suggest their tissue-specific role from source to sink. Sucrose cleavage has generally been linked with invertase isozymes, whereas sucrose synthase (SuSy)-catalyzed metabolism has been associated with biosynthetic processes such as UDP-Glc, cellulose, hemicellulose and other polymers. However, other two key sucrose-metabolizing enzymes, such as sucrose-6-phosphate phosphohydrolase (S6PP) and sucrose phosphate synthase (SPS) isoforms, have been linked with sucrose biosynthesis. These findings suggest that manipulation of genes, such as overexpression of SPS genes and sucrose transporter genes, silencing of the SC-uORF of Scbzip44 (removing the 5' leader region of the main ORF that is called SIRT-Insensitive) and downregulation of the invertase genes, may lead to maximum sucrose accumulation. This review provides an overview of sugarcane sucrose-regulating systems and baseline information for the development of cultivars with higher sucrose accumulation.

Sucrose synthase gene SUS3 could enhance cold tolerance in tomato

Tomatoes are susceptible to damage from cold temperatures in all stages of growth. Therefore, it is important to identify genetic resources and genes that can enhance tomato's ability to tolerate cold. In this study, a population of 223 tomato accessions was used to identify the sensitivity or tolerance of plants to cold stress. Transcriptome analysis of these accessions revealed that SUS3, a member of the sucrose synthase gene family, was induced by cold stress. We further investigated the role of SUS3 in cold stress by overexpression (OE) and RNA interference (RNAi). Compared with the wild type, SUS3-OE lines accumulated less MDA and electrolyte leakage and more proline and soluble sugar, maintained higher activities of SOD and CAT, reduced superoxide radicals, and suffered less membrane damage under cold. Thus, our findings indicate that SUS3 plays a crucial role in the response to cold stress. This study indicates that SUS3 may serve as a direct target for genetic engineering and improvement projects, which aim to augment the cold tolerance of tomato crops.

Relationship between Sugarcane eIF4E Gene and Resistance against Sugarcane Streak Mosaic Virus

Sugarcane mosaic disease, mainly caused by Sugarcane streak mosaic virus (SCSMV), has serious adverse effects on the yield and quality of sugarcane. Eukaryotic translation initiation factor 4E (eIF4E) is a natural resistance gene in plants. The eIF4E-mediated natural recessive resistance results from non-synonymous mutations of the eIF4E protein. In this study, two sugarcane varieties, CP94-1100 and ROC22, were selected for analysis of their differences in resistance to SCSMV. Four-base missense mutations in the ORF region of eIF4E resulted in different conserved domains. Therefore, the differences in resistance to SCSMV are due to the inherent differences in eIF4E of the sugarcane varieties. The coding regions of eIF4E included 28 SNP loci and no InDel loci, which were affected by negative selection and were relatively conserved. A total of 11 haploids encoded 11 protein sequences. Prediction of the protein spatial structure revealed three non-synonymous mutation sites for amino acids located in the cap pocket of eIF4E; one of these sites existed only in a resistant material (Yuetang 55), whereas the other site existed only in a susceptible material (ROC22), suggesting that these two sites might be related to the resistance to SCSMV. The results provide a strong basis for further analysis of the functional role of eIF4E in regulating mosaic resistance in sugarcane.

Haplotype variations of sucrose phosphate synthase B gene among sugarcane accessions with different sucrose content

BACKGROUND

Sucrose phosphate synthase B (SPSB) gene encoding an important rate-limiting enzyme for sucrose synthesis in sugarcane is mainly expressed on leaves, where its alleles control sucrose synthesis. In this study, genetic variation of SPSB gene represented by different haplotypes in sugarcane was investigated in hybrid clones with high and low sugar content and various accessory species.

RESULTS

A total of 39 haplotypes of SPSB gene with 2, 824 bp in size were identified from 18 sugarcane accessions. These haplotypes mainly distributed on Chr3B, Chr3C, and Chr3D according to the AP85-441 reference genome. Single nucleotide polymorphisms (SNPs) and insertion/deletion (InDels) were very dense (42 bp/sequence variation) including 44 transitional and 23 transversional SNPs among the 39 haplotypes. The sequence diversity related Hd, Eta, and Pi values were all lower in clones of high sucrose content (HS) than those in clones of low sucrose content (LS). The evolutionary network analysis showed that about half SPSB haplotypes (19 out of 39) were clustered into one group, including 6 (HAP4, HAP6, HAP7, HAP9, HAP17 and HAP20) haplotypes under positive selection in comparison to HAP26 identified in Badila (S. officinarum), an ancestry noble cane species and under purification selection (except HAP19 under neutral selection) in comparison to HAP18 identified from India1 (S. spontaneum), an ancestry species with low sugar content but high stress tolerance. The average number of haplotypes under positive selection in HS clones was twice as that in LS. Most of the SNPs and InDels sequence variation sites were positively correlated with sucrose and fiber content and negatively correlated with reducing sugar.

CONCLUSIONS

A total of 39 haplotypes of SPSB gene were identified in this study. Haplotypes potentially associated with high sucrose synthesis efficiency were identified. The mutations of SPSB haplotypes in HS were favorable and tended to be selected and fixed. The results of this study are informative and beneficial to the molecular assisted breeding of sucrose synthesis in sugarcane in the future.

A short review on sugarcane: its domestication, molecular manipulations and future perspectives

Sugarcane (Saccharum spp.) is a special crop plant that underwent anthropogenic evolution from a wild grass species to an important food, fodder, and energy crop. Unlike any other grass species which were selected for their kernels, sugarcane was selected for its high stem sucrose accumulation. Flowering in sugarcane is not favored since flowering diverts the stored sugar resources for the reproductive and developmental energy needs. Cultivars are vegetatively propagated and sugarcane breeding is still essentially focused on conventional methods, since the knowledge of sugarcane genetics has lagged that of other major crops. Cultivar improvement has been extremely challenging due to its polyploidy and aneuploidy nature derived from a few interspecific hybridizations between Saccharum officinarum and Saccharum spontaneum, revealing the coexistence of two distinct genome organization modes in the modern variety. Alongside implementation of modern agricultural techniques, generation of hybrid clones, transgenics and genome edited events will help to meet the ever-growing bioenergy needs. Additionally, there are two common biotechnological approaches to improve plant stress tolerance, which includes marker-assisted selection (MAS) and genetic transformation. During the past two decades, the use of molecular approaches has contributed greatly to a better understanding of the genetic and biochemical basis of plant stress-tolerance and in some cases, it led to the development of plants with enhanced tolerance to abiotic stress. Hence, this review mainly intends on the events that shaped the sugarcane as what it is now and what challenges ahead and measures taken to further improve its yield, production and maximize utilization to beat the growing demands.

Transcriptome Dynamics Underlying Magnesium Deficiency Stress in Three Founding Saccharum Species

Modern sugarcane cultivars were generated through interspecific crossing of the stress resistance Saccharum spontaneum and the high sugar content Saccharum officinarum which was domesticated from Saccharum robustum. Magnesium deficiency (MGD) is particularly prominent in tropical and subtropical regions where sugarcane is grown, but the response mechanism to MGD in sugarcane remains unknown. Physiological and transcriptomic analysis of the three founding Saccharum species under different magnesium (Mg) levels was performed. Our result showed that MGD decreased chlorophyll content and photosynthetic efficiency of three Saccharum species but led to increased starch in leaves and lignin content in roots of Saccharum robustum and Saccharum spontaneum. We identified 12,129, 11,306 and 12,178 differentially expressed genes (DEGs) of Saccharum officinarum, Saccharum robustum and Saccharum spontaneum, respectively. In Saccharum officinarum, MGD affected signal transduction by up-regulating the expression of xylan biosynthesis process-related genes. Saccharum robustum, responded to the MGD by regulating the expression of transcription and detoxification process-related genes. Saccharum spontaneum, avoids damage from MGD by regulating the expression of the signing transduction process and the transformation from growth and development to reproductive development. This novel repertoire of candidate genes related to MGD response in sugarcane will be helpful for engineering MGD tolerant varieties.

Functional characterization and analysis of transcriptional regulation of sugar transporter SWEET13c in sugarcane Saccharum spontaneum

BACKGROUND

Sugarcane is an important crop for sugar production worldwide. The Sugars Will Eventually be Exported Transporters (SWEETs) are a group of sugar transporters recently identified in sugarcane. In Saccharum spontaneum, SsSWEET13c played a role in the sucrose transportation from the source to the sink tissues, which was found to be mainly active in the mature leaf. However, the function and regulation of SWEETs in sugarcane remain elusive despite extensive studies performed on sugar metabolism.

RESULTS

In this study, we showed that SsSWEET13c is a member of SWEET gene family in S. spontaneum, constituting highest circadian rhythm-dependent expression. It is a functional gene that facilitates plant root elongation and increase fresh weight of Arabidopsis thaliana, when overexpressed. Furthermore, yeast one-hybrid assays indicate that 20 potential transcription factors (TFs) could bind to the SsSWEET13c promoter in S. spontaneum. We combined transcriptome data from developmental gradient leaf with distinct times during circadian cycles and stems/leaves at different growth stages. We have uncovered that 14 out of 20 TFs exhibited positive/negative gene expression patterns relative to SsSWEET13c. In the source tissues, SsSWEET13c was mainly positively regulated by SsbHLH34, SsTFIIIA-a, SsMYR2, SsRAP2.4 and SsbHLH035, while negatively regulated by SsABS5, SsTFIIIA-b and SsERF4. During the circadian rhythm, it was noticed that SsSWEET13c was more active in the morning than in the afternoon. It was likely due to the high level of sugar accumulation at night, which was negatively regulated by SsbZIP44, and positively regulated by SsbHLH34. Furthermore, in the sink tissues, SsSWEET13c was also active for sugar accumulation, which was positively regulated by SsbZIP44, SsTFIIIA-b, SsbHLH34 and SsTFIIIA-a, and negatively regulated by SsERF4, SsHB36, SsDEL1 and SsABS5. Our results were further supported by one-to-one yeast hybridization assay which verified that 12 potential TFs could bind to the promoter of SsSWEET13c.

CONCLUSIONS

A module of the regulatory network was proposed for the SsSWEET13c in the developmental gradient of leaf and circadian rhythm in S. spontaneum. These results provide a novel understanding of the function and regulation of SWEET13c during the sugar transport and biomass production in S. spontaneum.

Expression characterization and cross-species complementation uncover the functional conservation of YABBY genes for leaf abaxial polarity and carpel polarity establishment in Saccharum spontaneum

BACKGROUND

Cell polarity establishment and maintenance is indispensable for plant growth and development. In plants, the YABBY transcription factor family has a distinct role in leaf asymmetric polarity establishment and lateral organ initiation. However, for the important sugar crop Saccharum, little information on YABBY genes is available.

RESULTS

In this study, a total of 20 sequences for 7 SsYABBY genes were identified in the sugarcane genome, designated as SsYABBY1-7 based on their chromosome locations, and characterized by phylogenetic analysis. We provided a high-resolution map of SsYABBYs' global expression dynamics during vegetative and reproductive organ morphogenesis and revealed that SsYABBY3/4/5 are predominately expressed at the seedling stage of stem and leaf basal zone; SsYABBY2/5/7 are highly expressed in ovules. Besides, cross-species overexpression and/or complementation verified the conserved function of SsYABBY2 in establishing leaf adaxial-abaxial polarity and ovules development. We found that the SsYABBY2 could successfully rescue the leaves curling, carpel dehiscence, and ovule abortion defects in Arabidopsis crc mutant.

CONCLUSIONS

Collectively, our study demonstrates that SsYABBY genes retained a conserved function in establishing and preserving leaf adaxial-abaxial polarity and lateral organ development during evolution.

Genomic analysis of SBP gene family in Saccharum spontaneum reveals their association with vegetative and reproductive development

Background

SQUAMOSA promoter binding proteins (SBPs) genes encode a family of plant-specific transcription factors involved in various growth and development processes, including flower and fruit development, leaf initiation, phase transition, and embryonic development. The SBP gene family has been identified and characterized in many species, but no systematic analysis of the SBP gene family has been carried out in sugarcane.

Results

In the present study, a total of 50 sequences for 30 SBP genes were identified by the genome-wide analysis and designated SsSBP1 to SsSBP30 based on their chromosomal distribution. According to the phylogenetic tree, gene structure and motif features, the SsSBP genes were classified into eight groups (I to VIII). By synteny analysis, 27 homologous gene pairs existed in SsSBP genes, and 37 orthologous gene pairs between sugarcane and sorghum were found. Expression analysis in different tissues, including vegetative and reproductive organs, showed differential expression patterns of SsSBP genes, indicating their functional diversity in the various developmental processes. Additionally, 22 SsSBP genes were predicted as the potential targets of miR156. The differential expression pattern of miR156 exhibited a negative correlation of transcription levels between miR156 and the SsSBP gene in different tissues.

Conclusions

The sugarcane genome possesses 30 SsSBP genes, and they shared similar gene structures and motif features in their subfamily. Based on the transcriptional and qRT-PCR analysis, most SsSBP genes were found to regulate the leaf initial and female reproductive development. The present study comprehensively and systematically analyzed SBP genes in sugarcane and provided a foundation for further studies on the functional characteristics of SsSBP genes during different development processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08090-3.

Genome-wide identification and expression profiling of DREB genes in Saccharum spontaneum

Background

The dehydration-responsive element-binding proteins (DREBs) are important transcription factors that interact with a DRE/CRT (C-repeat) sequence and involve in response to multiple abiotic stresses in plants. Modern sugarcane are hybrids from the cross between Saccharum spontaneum and Saccharum officinarum, and the high sugar content is considered to the attribution of S. officinaurm, while the stress tolerance is attributed to S. spontaneum. To understand the molecular and evolutionary characterization and gene functions of the DREBs in sugarcane, based on the recent availability of the whole genome information, the present study performed a genome-wide in silico analysis of DREB genes and transcriptome analysis in the polyploidy S. spontaneum.

Results

Twelve DREB1 genes and six DREB2 genes were identified in S. spontaneum genome and all proteins contained a conserved AP2/ERF domain. Eleven SsDREB1 allele genes were assumed to be originated from tandem duplications, and two of them may be derived after the split of S. spontaneum and the proximal diploid species sorghum, suggesting tandem duplication contributed to the expansion of DREB1-type genes in sugarcane. Phylogenetic analysis revealed that one DREB2 gene was lost during the evolution of sugarcane. Expression profiling showed different SsDREB genes with variable expression levels in the different tissues, indicating seven SsDREB genes were likely involved in the development and photosynthesis of S. spontaneum. Furthermore, SsDREB1F, SsDREB1L, SsDREB2D, and SsDREB2F were up-regulated under drought and cold condition, suggesting that these four genes may be involved in both dehydration and cold response in sugarcane.

Conclusions

These findings demonstrated the important role of DREBs not only in the stress response, but also in the development and photosynthesis of S. spontaneum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07799-5.

Genome-wide identification and analysis of the sucrose synthase gene family in cassava (Manihot esculenta Crantz)

Sucrose synthase (SUS), a key enzyme of the sucrose metabolism pathway, is encoded by a multi-gene family in plants. To date, dozens of SUS gene families have been characterized in various plant genomes. However, only a few studies have performed comprehensive analyses in tropical crops like cassava (Manihot esculenta Crantz). In the present study, seven non-redundant members of the SUS gene family (MeSUS1-7) were identified and characterized from the cassava genome. The MeSUS genes were distributed on five chromosomes (Chr1, Chr2, Chr3, Chr14, and Chr16) and the encoded proteins could be classified into three major groups with other SUS proteins from both dicot and monocot species (SUS I, SUS II, and SUS III). The spatio-temporal expression profiles of MeSUS genes showed a developmental stage-dependent, partially overlapping pattern, mainly expressed in the source and sink tissues. Cold and drought treatments significantly induced the expressions of MeSUS2, MeSUS4, MeSUS6, and MeSUS7 and the activities of the encoded enzymes, indicating that these genes may play crucial roles in resistance against abiotic stresses. These results provide new insights into the multifaceted role of the SUS gene family members in various physiological processes, especially sucrose transport and starch accumulation in cassava roots.

Comparative analysis of sucrose phosphate synthase (SPS) gene family between Saccharum officinarum and Saccharum spontaneum

BACKGROUND

Sucrose phosphate synthase (SPS) genes play vital roles in sucrose production across various plant species. Modern sugarcane cultivar is derived from the hybridization between the high sugar content species Saccharum officinarum and the high stress tolerance species Saccharum spontaneum, generating one of the most complex genomes among all crops. The genomics of sugarcane SPS remains under-studied despite its profound impact on sugar yield.

RESULTS

In the present study, 8 and 6 gene sequences for SPS were identified from the BAC libraries of S. officinarum and S. spontaneum, respectively. Phylogenetic analysis showed that SPSD was newly evolved in the lineage of Poaceae species with recently duplicated genes emerging from the SPSA clade. Molecular evolution analysis based on Ka/Ks ratios suggested that polyploidy reduced the selection pressure of SPS genes in Saccharum species. To explore the potential gene functions, the SPS expression patterns were analyzed based on RNA-seq and proteome dataset, and the sugar content was detected using metabolomics analysis. All the SPS members presented the trend of increasing expression in the sink-source transition along the developmental gradient of leaves, suggesting that the SPSs are involved in the photosynthesis in both Saccharum species as their function in dicots. Moreover, SPSs showed the higher expression in S. spontaneum and presented expressional preference between stem (SPSA) and leaf (SPSB) tissue, speculating they might be involved in the differentia of carbohydrate metabolism in these two Saccharum species, which required further verification from experiments.

CONCLUSIONS

SPSA and SPSB genes presented relatively high expression and differential expression patterns between the two Saccharum species, indicating these two SPSs are important in the formation of regulatory networks and sucrose traits in the two Saccharum species. SPSB was suggested to be a major contributor to the sugar accumulation because it presented the highest expressional level and its expression positively correlated with sugar content. The recently duplicated SPSD2 presented divergent expression levels between the two Saccharum species and the relative protein content levels were highest in stem, supporting the neofunctionalization of the SPSD subfamily in Saccharum.

Genome-Wide Identification and Expression Profile Analysis of WRKY Family Genes in the Autopolyploid Saccharum spontaneum

WRKY is one of the largest transcription factor families in plants and plays important roles in the regulation of developmental and physiological processes. To date, the WRKY gene family has not been identified in Saccharum species because of its complex polyploid genome. In this study, a total of 294 sequences for 154 SsWRKY genes were identified in the polyploid Saccharum spontaneum genome and then named on the basis of their chromosome locations, including 13 (8.4%) genes with four alleles, 29 (18.8%) genes with three alleles and 41 (26.6%) genes with two alleles. Among them, 73.8% and 16.0% of the SsWRKY genes originated from segmental duplications and tandem duplications, respectively. The WRKY members exhibited conserved gene structures and amino acid sequences among the allelic haplotypes, which were accompanied by variations in intron sizes. Phylogenetic and collinearity analyses revealed that 27 SsWRKYs originated after the split of sorghum and Saccharum, resulting in a significantly higher number of WRKYs in sugarcane than in the proximal diploid species sorghum. The analysis of RNA-seq data revealed that SsWRKYs' expression profiles in 46 different samples including different developmental stages revealed distinct temporal and spatial patterns with 52 genes expressed in all tissues, four genes not expressed in any tissues and 21 SsWRKY genes likely to be involved in photosynthesis. The comprehensive analysis of SsWRKYs' expression will provide an important and valuable foundation for further investigation of the regulatory mechanisms of WRKYs in physiological roles in sugarcane S. spontaneum.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop

BACKGROUND

Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing.

RESULTS

Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences.

CONCLUSIONS

This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

Analysis of the diversity and tissue specificity of sucrose synthase genes in the long read transcriptome of sugarcane

BACKGROUND

Sugarcane accumulates very high levels of sucrose in the culm. Elucidation of the molecular mechanisms that allows such high sucrose synthesis and accumulation (up to 650 mM) is made difficult by the complexity of the highly polyploid genome. Here we report the use of RNA Seq data to characterize the sucrose synthase (SuSy) genes expressed in the transcriptome of the mature sugarcane plant.

RESULTS

Four SuSy gene families were identified in the sugarcane Iso-Seq long read transcriptome (SUGIT) through gene annotation of transcripts that mapped to reference SuSy genes from sorghum and maize. In total, 38, 19, 14, and 2 transcripts were identified for the four corresponding SuSy genes 1, 2, 4 and 7, respectively. Comparative studies using available SuSy genes from sorghum (1, 2, 4, 6, 7) and maize (1-7) revealed that the sugarcane SuSy genes were interrupted by multiple introns and that they share a highly conserved gene structure. Spatial expression of the four SuSy genes in sugarcane genotypes and in the progenitor species, Saccharum spontaneum and Saccharum officinarum, was studied in the leaf and root tissues and also in three regions of the culm tissue; top, middle and bottom internodes. Expression profiles indicated that all SuSy transcripts were differentially expressed between the top and bottom tissues, with high expression in the top tissues, lower expression in the bottom and moderate expression in the middle, indicating a gradient of SuSy activity in the sugarcane culm. Further, the root tissue had similar expression levels to that of the top internodes while leaf tissues showed lower expression. In the progenitors, SuSy7 was found to be highly expressed in S. officinarum while the other three SuSy genes had moderate expression in both the progenitors.

CONCLUSIONS

The high expression of the SuSy genes in sink tissues, the top internodes and the roots suggests functional roles in sucrose utilization to support growth. The SuSy7 gene has not been previously reported in sugarcane. As sugarcane is unique in storing such high amounts of sucrose, it is possible that there are more SuSy genes/isoforms with specific expression patterns to be discovered in this complex system.

Genome-wide identification and expression analysis of sucrose synthase genes in allotetraploid Brassica juncea

Sucrose plays pivotal role in energy metabolism and regulating gene expression of several physiological processes in higher plants. Here, fourteen sucrose synthase (SUS) genes have been identified in the allotetraploid genome of Indian mustard, Brassica juncea. The identified SUS genes in B. juncea (BjSUS) were derived from the two-progenitor species, B. rapa and B. nigra. Intron-exon analysis indicated loss or gain of 1-3 introns in diversification of SUS gene family. Phylogenetic analysis revealed discrete evolutionary paths for the BjSUS genes, originating from three ancestor groups, SUS I, SUS II and SUS III. Gene expression study revealed significant variability in expression of the BjSUS paralogs across the different tissues. BjSUS genes showed transcriptional activation in response to defense hormones and a late response to wounding. Tissue and temporal specificity of expression revealed importance of specific SUS paralogs at different developmental stages and under different stress conditions. The study highlighted differential involvement of SUS paralogs in sucrose metabolism across the tissues and stress-responses, in a major oilseed crop B. juncea.

An Overview of Sucrose Synthases in Plants

Sucrose is the end product of photosynthesis and the primary sugar transported in the phloem of most plants. Sucrose synthase (SuSy) is a glycosyl transferase enzyme that plays a key role in sugar metabolism, primarily in sink tissues. SuSy catalyzes the reversible cleavage of sucrose into fructose and either uridine diphosphate glucose (UDP-G) or adenosine diphosphate glucose (ADP-G). The products of sucrose cleavage by SuSy are available for many metabolic pathways, such as energy production, primary-metabolite production, and the synthesis of complex carbohydrates. SuSy proteins are usually homotetramers with an average monomeric molecular weight of about 90 kD (about 800 amino acids long). Plant SuSy isozymes are mainly located in the cytosol or adjacent to plasma membrane, but some SuSy proteins are found in the cell wall, vacuoles, and mitochondria. Plant SUS gene families are usually small, containing between four to seven genes, with distinct exon-intron structures. Plant SUS genes are divided into three separate clades, which are present in both monocots and dicots. A comprehensive phylogenetic analysis indicates that a first SUS duplication event may have occurred before the divergence of the gymnosperms and angiosperms and a second duplication event probably occurred in a common angiosperm ancestor, leading to the existence of all three clades in both monocots and dicots. Plants with reduced SuSy activity have been shown to have reduced growth, reduced starch, cellulose or callose synthesis, reduced tolerance to anaerobic-stress conditions and altered shoot apical meristem function and leaf morphology. Plants overexpressing SUS have shown increased growth, increased xylem area and xylem cell-wall width, and increased cellulose and starch contents, making SUS high-potential candidate genes for the improvement of agricultural traits in crop plants. This review summarizes the current knowledge regarding plant SuSy, including newly discovered possible developmental roles for SuSy in meristem functioning that involve sugar and hormonal signaling.

Comparative genomics revealed the gene evolution and functional divergence of magnesium transporter families in Saccharum

Background

Sugarcane served as the model plant for discovery of the C₄ photosynthetic pathway. Magnesium is the central atom of chlorophyll, and thus is considered as a critical nutrient for plant development and photosynthesis. In plants, the magnesium transporter (MGT) family is composed of a number of membrane proteins, which play crucial roles in maintaining Mg homeostasis. However, to date there is no information available on the genomics of MGTs in sugarcane due to the complexity of the Saccharum genome.

Results

Here, we identified 10 MGTs from the Saccharum spontaneum genome. Phylogenetic analysis of MGTs suggested that the MGTs contained at least 5 last common ancestors before the origin of angiosperms. Gene structure analysis suggested that MGTs family of dicotyledon may be accompanied by intron loss and pseudoexon phenomena during evolution. The pairwise synonymous substitution rates corresponding to a divergence time ranged from 142.3 to 236.6 Mya, demonstrating that the MGTs are an ancient gene family in plants. Both the phylogeny and Ks analyses indicated that SsMGT1/SsMGT2 originated from the recent ρWGD, and SsMGT7/SsMGT8 originated from the recent σ WGD. These 4 recently duplicated genes were shown low expression levels and assumed to be functionally redundant. MGT6, MGT9 and MGT10 weredominant genes in the MGT family and werepredicted to be located inthe chloroplast. Of the 3 dominant MGTs, SsMGT6 expression levels were found to be induced in the light period, while SsMGT9 and SsMTG10 displayed high expression levels in the dark period. These results suggested that SsMGT6 may have a function complementary to SsMGT9 and SsMTG10 that follows thecircadian clock for MGT in the leaf tissues of S. spontaneum. MGT3, MGT7 and MGT10 had higher expression levels Insaccharum officinarum than in S. spontaneum, suggesting their functional divergence after the split of S. spontaneum and S. officinarum.

Conclusions

This study of gene evolution and expression of MGTs in S. spontaneum provided basis for the comprehensive genomic study of the entire MGT genes family in Saccharum. The results are valuable for further functional analyses of MGT genes and utilization of the MGTs for Saccharum genetic improvement.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5437-3) contains supplementary material, which is available to authorized users.

Recent polyploidization events in three Saccharum founding species

The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F₁ individuals derived from S. officinarumLA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole-genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole-genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long-standing hypothesis and clarified the timing of speciation and whole-genome duplication events in Saccharum.

New insights into the evolution and functional divergence of the SWEET family in Saccharum based on comparative genomics

BACKGROUND

The SWEET (Sugars Will Eventually be Exported Transporters) gene family is a recently identified group of sugar transporters that play an indispensable role in sugar efflux, phloem loading, plant-pathogen interaction, nectar secretion, and reproductive tissue development. However, little information on Saccharum SWEET is available for this crop with a complex genetic background.

RESULTS

In this study, 22 SWEET genes were identified from Saccharum spontaneum Bacterial Artificial Chromosome libraries sequences. Phylogenetic analyses of SWEETs from 11 representative plant species showed that gene expansions of the SWEET family were mainly caused by the recent gene duplication in dicot plants, while these gene expansions were attributed to the ancient whole genome duplication (WGD) in monocot plant species. Gene expression profiles were obtained from RNA-seq analysis. SWEET1a and SWEET2s had higher expression levels in the transitional zone and maturing zone than in the other analyzed zones. SWEET1b was mainly expressed in the leaf tissues and the mature zone of the leaf of both S. spontaneum and S. officinarum, and displayed a peak in the morning and was undetectable in both sclerenchyma and parenchyma cells from the mature stalks of S. officinarum. SsSWEET4a\4b had higher expression levels than SWEET4c and were mainly expressed in the stems of seedlings and mature plants. SWEET13s are recently duplicated genes, and the expression of SWEET13s dramatically increased from the maturing to mature zones. SWEET16b's expression was not detected in S. officinarum, but displayed a rhythmic diurnal expression pattern.

CONCLUSIONS

Our study revealed the gene evolutionary history of SWEETs in Saccharum and SWEET1b was found to be a sucrose starvation-induced gene involved in the sugar transportation in the high photosynthetic zones. SWEET13c was identified as the key player in the efflux of sugar transportation in mature photosynthetic tissues. SWEET4a\4b were found to be mainly involved in sugar transportation in the stalk. SWEET1a\2a\4a\4b\13a\16b were suggested to be the genes contributing to the differences in sugar contents between S. spontaneum and S. officinarum. Our results are valuable for further functional analysis of SWEET genes and utilization of the SWEET genes for genetic improvement of Saccharum for biofuel production.

Luxurious Nitrogen Fertilization of Two Sugar Cane Genotypes Contrasting for Lignin Composition Causes Changes in the Stem Proteome Related to Carbon, Nitrogen, and Oxidant Metabolism but Does Not Alter Lignin Content

Sugar cane is an important crop for sugar and biofuel production. Its lignocellulosic biomass represents a promising option as feedstock for second-generation ethanol production. Nitrogen fertilization can affect differently tissues and its biopolymers, including the cell-wall polysaccharides and lignin. Lignin content and composition are the most important factors associated with biomass recalcitrance to convert cell-wall polysaccharides into fermentable sugars. Thus it is important to understand the metabolic relationship between nitrogen fertilization and lignin in this feedstock. In this study, a large-scale proteomics approach based on GeLC-MS/MS was employed to identify and relatively quantify proteins differently accumulated in two contrasting genotypes for lignin composition after excessive nitrogen fertilization. From the ∼1000 nonredundant proteins identified, 28 and 177 were differentially accumulated in response to nitrogen from IACSP04-065 and IACSP04-627 lines, respectively. These proteins were associated with several functional categories, including carbon metabolism, amino acid metabolism, protein turnover, and oxidative stress. Although nitrogen fertilization has not changed lignin content, phenolic acids and lignin composition were changed in both species but not in the same way. Sucrose and reducing sugars increased in plants of the genotype IACSP04-065 receiving nitrogen.

Molecular cloning, structure, phylogeny and expression analysis of the invertase gene family in sugarcane

BACKGROUND

Invertases (INVs) are key enzymes regulating sucrose metabolism and are here revealed to be involved in responses to environmental stress in plants. To date, individual members of the invertase gene family and their expression patterns are unknown in sugarcane due to its complex genome despite their significance in sucrose metabolism.

RESULTS

In this study, based on comparative genomics, eleven cDNA and twelve DNA sequences belonging to 14 non-redundant members of the invertase gene family were successfully cloned from sugarcane. A comprehensive analysis of the invertase gene family was carried out, including gene structures, phylogenetic relationships, functional domains, conserved motifs of proteins. The results revealed that the 14 invertase members from sugarcane could be clustered into three subfamilies, including 6 neutral/alkaline invertases (ShN/AINVs), and 8 acid invertases (ShAINVs). Faster divergence occurred in acid INVs than in neutral/alkaline INVs after the split of sugarcane and sorghum. At least a one-time gene duplication event was observed to have occurred in the four groups of acid INVs, whereas ShN/AINV1 and ShN/AINV2 in the β8 lineage were revealed to be the most recently duplicated genes among their paralogous genes in the β group of N/AINVs. Furthermore, comprehensive expression analysis of these genes was performed in sugarcane seedlings subjected to five abiotic stresses (drought, low temperature, glucose, fructose, and sucrose) using Quantitative Real-time PCR. The results suggested a functional divergence of INVs and their potential role in response to the five different treatments. Enzymatic activity in sugarcane seedlings was detected under five abiotic stresses treatments, and showed that the activities of all INVs were significantly inhibited in response to five different abiotic stresses, and that the neutral/alkaline INVs played a more prominent role in abiotic stresses than the acid INVs.

CONCLUSIONS

In this study, we determined the INV gene family members of sugarcane by PCR cloning using sorghum as a reference, providing the first study of the INV gene family in sugarcane. Combining existing INV gene data from 7 plants with a comparative approach including a series of comprehensive analyses to isolate and identify INV gene family members proved to be highly successful. Moreover, the expression levels of INV genes and the variation of enzymatic activities associated with drought, low temperature, glucose, fructose, and sucrose are reported in sugarcane for the first time. The results offered useful foundation and framework for future research for understanding the physiological roles of INVs for sucrose accumulation in sugarcane.

A survey of the complex transcriptome from the highly polyploid sugarcane genome using full-length isoform sequencing and de novo assembly from short read sequencing

BACKGROUND

Despite the economic importance of sugarcane in sugar and bioenergy production, there is not yet a reference genome available. Most of the sugarcane transcriptomic studies have been based on Saccharum officinarum gene indices (SoGI), expressed sequence tags (ESTs) and de novo assembled transcript contigs from short-reads; hence knowledge of the sugarcane transcriptome is limited in relation to transcript length and number of transcript isoforms.

RESULTS

The sugarcane transcriptome was sequenced using PacBio isoform sequencing (Iso-Seq) of a pooled RNA sample derived from leaf, internode and root tissues, of different developmental stages, from 22 varieties, to explore the potential for capturing full-length transcript isoforms. A total of 107,598 unique transcript isoforms were obtained, representing about 71% of the total number of predicted sugarcane genes. The majority of this dataset (92%) matched the plant protein database, while just over 2% was novel transcripts, and over 2% was putative long non-coding RNAs. About 56% and 23% of total sequences were annotated against the gene ontology and KEGG pathway databases, respectively. Comparison with de novo contigs from Illumina RNA-Sequencing (RNA-Seq) of the internode samples from the same experiment and public databases showed that the Iso-Seq method recovered more full-length transcript isoforms, had a higher N50 and average length of largest 1,000 proteins; whereas a greater representation of the gene content and RNA diversity was captured in RNA-Seq. Only 62% of PacBio transcript isoforms matched 67% of de novo contigs, while the non-matched proportions were attributed to the inclusion of leaf/root tissues and the normalization in PacBio, and the representation of more gene content and RNA classes in the de novo assembly, respectively. About 69% of PacBio transcript isoforms and 41% of de novo contigs aligned with the sorghum genome, indicating the high conservation of orthologs in the genic regions of the two genomes.

CONCLUSIONS

The transcriptome dataset should contribute to improved sugarcane gene models and sugarcane protein predictions; and will serve as a reference database for analysis of transcript expression in sugarcane.

Co-expression network analysis reveals transcription factors associated to cell wall biosynthesis in sugarcane

Sugarcane is a hybrid of Saccharum officinarum and Saccharum spontaneum, with minor contributions from other species in Saccharum and other genera. Understanding the molecular basis of cell wall metabolism in sugarcane may allow for rational changes in fiber quality and content when designing new energy crops. This work describes a comparative expression profiling of sugarcane ancestral genotypes: S. officinarum, S. spontaneum and S. robustum and a commercial hybrid: RB867515, linking gene expression to phenotypes to identify genes for sugarcane improvement. Oligoarray experiments of leaves, immature and intermediate internodes, detected 12,621 sense and 995 antisense transcripts. Amino acid metabolism was particularly evident among pathways showing natural antisense transcripts expression. For all tissues sampled, expression analysis revealed 831, 674 and 648 differentially expressed genes in S. officinarum, S. robustum and S. spontaneum, respectively, using RB867515 as reference. Expression of sugar transporters might explain sucrose differences among genotypes, but an unexpected differential expression of histones were also identified between high and low Brix° genotypes. Lignin biosynthetic genes and bioenergetics-related genes were up-regulated in the high lignin genotype, suggesting that these genes are important for S. spontaneum to allocate carbon to lignin, while S. officinarum allocates it to sucrose storage. Co-expression network analysis identified 18 transcription factors possibly related to cell wall biosynthesis while in silico analysis detected cis-elements involved in cell wall biosynthesis in their promoters. Our results provide information to elucidate regulatory networks underlying traits of interest that will allow the improvement of sugarcane for biofuel and chemicals production.

Transcript Analysis and Regulative Events during Flower Development in Olive (Olea europaea L.)

The identification and characterization of transcripts involved in flower organ development, plant reproduction and metabolism represent key steps in plant phenotypic and physiological pathways, and may generate high-quality transcript variants useful for the development of functional markers. This study was aimed at obtaining an extensive characterization of the olive flower transcripts, by providing sound information on the candidate MADS-box genes related to the ABC model of flower development and on the putative genetic and molecular determinants of ovary abortion and pollen-pistil interaction. The overall sequence data, obtained by pyrosequencing of four cDNA libraries from flowers at different developmental stages of three olive varieties with distinct reproductive features (Leccino, Frantoio and Dolce Agogia), included approximately 465,000 ESTs, which gave rise to more than 14,600 contigs and approximately 92,000 singletons. As many as 56,700 unigenes were successfully annotated and provided gene ontology insights into the structural organization and putative molecular function of sequenced transcripts and deduced proteins in the context of their corresponding biological processes. Differentially expressed genes with potential regulatory roles in biosynthetic pathways and metabolic networks during flower development were identified. The gene expression studies allowed us to select the candidate genes that play well-known molecular functions in a number of biosynthetic pathways and specific biological processes that affect olive reproduction. A sound understanding of gene functions and regulatory networks that characterize the olive flower is provided.

Structure, phylogeny, allelic haplotypes and expression of sucrose transporter gene families in Saccharum

BACKGROUND

Sugarcane is an economically important crop contributing to about 80% of the world sugar production. Increasing efforts in molecular biological studies have been performed for improving the sugar yield and other relevant important agronomic traits. However, due to sugarcane's complicated genomes, it is still challenging to study the genetic basis of traits, such as sucrose accumulation. Sucrose transporters (SUTs) are critical for both phloem loading in source tissue and sucrose uptaking in sink tissue, and are considered to be the control points for regulating sucrose storage. However, no genomic study for sugarcane sucrose transporter (SsSUT) families has been reported up to date.

RESULTS

By using comparative genomics and bacterial artificial chromosomes (BACs), six SUT genes were identified and characterized in S. spontaenum. Phylogenetic analyses revealed that the two pairs SsSUTs (SsSUT1/SsSUT3 and SsSUT5/SsSUT6) could be clustered together into two separate monocot specific SUT groups, while SsSUT2 and SsSUT4 were separated into the other two groups, with members from both dicot and monocot species. Gene structure comparison demonstrated that the number and position of exons/introns in SUTs were highly conserved among the close orthologs; in contrast, there were variations among the paralogous SUTs in Sacchuarm. Though with the high polyploidy level, gene allelic haplotype comparative analysis showed that the examined four SsSUT members exhibited conservations of gene structures and amino acid sequences among the allelic haplotypes accompanied by variations of intron sizes. Gene expression analyses were performed for tissues from seedlings under drought stress and mature plants of three Saccharum species (S.officinarnum, S.spotaneum and S.robustum). Both SUT1 and SUT4 expressed abundantly at different conditions. SUT2 had similar expression level in all of the examined tissues, but SUT3 was undetectable. Both of SUT5 and SUT6 had lower expression level than other gene member, and expressed stronger in source leaves and are likely to play roles in phloem loading. In the seeding plant leave under water stress, four genes SUT1, SUT2, SUT4 and SUT5 were detectable. In these detectable genes, SUT1 and SUT4 were down regulated, while, SUT2 and SUT5 were up regulated.

CONCLUSIONS

In this study, we presented the first comprehensive genomic study for a whole gene family, the SUT family, in Saccharum. We speculated that there were six SUT members in the S. spotaneum genome. Out of the six members, SsSUTs, SsSUT5 and SsSUT6 were recent duplication genes accompanied by rapid evolution, while, SsSUT2 and SsSUT4 were the ancient members in the families. Despite the high polypoidy genome, functional redundancy may not exist among the SUTs allelic haplotypes supported by the evidence of strong purifying selection of the gene allele. SUT3 could be a low active member in the family because it is undetectable in our study, but it might not be a pseudogene because it harbored integrated gene structure. SUT1 and SUT4 were the main members for the sucrose transporter, while, these SUTs had sub-functional divergence in response to sucrose accumulation and plant development in Saccharum.

Natural Allelic Variations in Highly Polyploidy Saccharum Complex

Sugarcane (Saccharum spp.) is an important sugar and biofuel crop with high polyploid and complex genomes. The Saccharum complex, comprised of Saccharum genus and a few related genera, are important genetic resources for sugarcane breeding. A large amount of natural variation exists within the Saccharum complex. Though understanding their allelic variation has been challenging, it is critical to dissect allelic structure and to identify the alleles controlling important traits in sugarcane. To characterize natural variations in Saccharum complex, a target enrichment sequencing approach was used to assay 12 representative germplasm accessions. In total, 55,946 highly efficient probes were designed based on the sorghum genome and sugarcane unigene set targeting a total of 6 Mb of the sugarcane genome. A pipeline specifically tailored for polyploid sequence variants and genotype calling was established. BWA-mem and sorghum genome approved to be an acceptable aligner and reference for sugarcane target enrichment sequence analysis, respectively. Genetic variations including 1,166,066 non-redundant SNPs, 150,421 InDels, 919 gene copy number variations, and 1,257 gene presence/absence variations were detected. SNPs from three different callers (Samtools, Freebayes, and GATK) were compared and the validation rates were nearly 90%. Based on the SNP loci of each accession and their ploidy levels, 999,258 single dosage SNPs were identified and most loci were estimated as largely homozygotes. An average of 34,397 haplotype blocks for each accession was inferred. The highest divergence time among the Saccharum spp. was estimated as 1.2 million years ago (MYA). Saccharum spp. diverged from Erianthus and Sorghum approximately 5 and 6 MYA, respectively. The target enrichment sequencing approach provided an effective way to discover and catalog natural allelic variation in highly polyploid or heterozygous genomes.

Full-length enriched cDNA libraries and ORFeome analysis of sugarcane hybrid and ancestor genotypes

Sugarcane is a major crop used for food and bioenergy production. Modern cultivars are hybrids derived from crosses between Saccharum officinarum and Saccharum spontaneum. Hybrid cultivars combine favorable characteristics from ancestral species and contain a genome that is highly polyploid and aneuploid, containing 100–130 chromosomes. These complex genomes represent a huge challenge for molecular studies and for the development of biotechnological tools that can facilitate sugarcane improvement. Here, we describe full-length enriched cDNA libraries for Saccharum officinarum, Saccharum spontaneum, and one hybrid genotype (SP803280) and analyze the set of open reading frames (ORFs) in their genomes (i.e., their ORFeomes). We found 38,195 (19%) sugarcane-specific transcripts that did not match transcripts from other databases. Less than 1.6% of all transcripts were ancestor-specific (i.e., not expressed in SP803280). We also found 78,008 putative new sugarcane transcripts that were absent in the largest sugarcane expressed sequence tag database (SUCEST). Functional annotation showed a high frequency of protein kinases and stress-related proteins. We also detected natural antisense transcript expression, which mapped to 94% of all plant KEGG pathways; however, each genotype showed different pathways enriched in antisense transcripts. Our data appeared to cover 53.2% (17,563 genes) and 46.8% (937 transcription factors) of all sugarcane full-length genes and transcription factors, respectively. This work represents a significant advancement in defining the sugarcane ORFeome and will be useful for protein characterization, single nucleotide polymorphism and splicing variant identification, evolutionary and comparative studies, and sugarcane genome assembly and annotation.

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene-rich regions. Gene-enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene-enrichment strategy, we have compared assemblies using methyl-filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA-seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single-nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.