Building the sugarcane genome for biotechnology and identifying evolutionary trends

The extent of multiallelic, co-editing of LIGULELESS1 in highly polyploid sugarcane tunes leaf inclination angle and enables selection of the ideotype for biomass yield

Sugarcane (Saccharum spp. hybrid) is a prime feedstock for commercial production of biofuel and table sugar. Optimizing canopy architecture for improved light capture has great potential for elevating biomass yield. LIGULELESS1 (LG1) is involved in leaf ligule and auricle development in grasses. Here, we report CRISPR/Cas9-mediated co-mutagenesis of up to 40 copies/alleles of the putative LG1 in highly polyploid sugarcane (2n = 100-120, x = 10-12). Next generation sequencing revealed co-editing frequencies of 7.4%-100% of the LG1 reads in 16 of the 78 transgenic lines. LG1 mutations resulted in a tuneable leaf angle phenotype that became more upright as co-editing frequency increased. Three lines with loss of function frequencies of ~12%, ~53% and ~95% of lg1 were selected following a randomized greenhouse trial and grown in replicated, multi-row field plots. The co-edited LG1 mutations were stably maintained in vegetative progenies and the extent of co-editing remained constant in field tested lines L26 and L35. Next generation sequencing confirmed the absence of potential off targets. The leaf inclination angle corresponded to light transmission into the canopy and tiller number. Line L35 displaying loss of function in ~12% of the lg1 NGS reads exhibited an 18% increase in dry biomass yield supported by a 56% decrease in leaf inclination angle, a 31% increase in tiller number, and a 25% increase in internode number. The scalable co-editing of LG1 in highly polyploid sugarcane allows fine-tuning of leaf inclination angle, enabling the selection of the ideotype for biomass yield.

Leveraging the sugarcane CRISPR/Cas9 technique for genetic improvement of non-cultivated grasses

Under changing climatic scenarios, grassland conservation and development have become imperative to impart functional sustainability to their ecosystem services. These goals could be effectively and efficiently achieved with targeted genetic improvement of native grass species. To the best of our literature search, very scant research findings are available pertaining to gene editing of non-cultivated grass species (switch grass, wild sugarcane, Prairie cordgrass, Bermuda grass, Chinese silver grass, etc.) prevalent in natural and semi-natural grasslands. Thus, to explore this novel research aspect, this study purposes that gene editing techniques employed for improvement of cultivated grasses especially sugarcane might be used for non-cultivated grasses as well. Our hypothesis behind suggesting sugarcane as a model crop for genetic improvement of non-cultivated grasses is the intricacy of gene editing owing to polyploidy and aneuploidy compared to other cultivated grasses (rice, wheat, barley, maize, etc.). Another reason is that genome editing protocols in sugarcane (x = 10-13) have been developed and optimized, taking into consideration the high level of genetic redundancy. Thus, as per our knowledge, this review is the first study that objectively evaluates the concept and functioning of the CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 technique in sugarcane regarding high versatility, target specificity, efficiency, design simplicity, and multiplexing capacity in order to explore novel research perspectives for gene editing of non-cultivated grasses against biotic and abiotic stresses. Additionally, pronounced challenges confronting sugarcane gene editing have resulted in the development of different variants (Cas9, Cas12a, Cas12b, and SpRY) of the CRISPR tool, whose technicalities have also been critically assessed. Moreover, different limitations of this technique that could emerge during gene editing of non-cultivated grass species have also been highlighted.

Identification and functional characterization of the sugarcane (Saccharum spp.) AMT2-type ammonium transporter ScAMT3;3 revealed a presumed role in shoot ammonium remobilization

Sugarcane (Saccharum spp.) is an important crop for sugar and bioethanol production worldwide. To maintain and increase sugarcane yields in marginal areas, the use of nitrogen (N) fertilizers is essential, but N overuse may result in the leaching of reactive N to the natural environment. Despite the importance of N in sugarcane production, little is known about the molecular mechanisms involved in N homeostasis in this crop, particularly regarding ammonium (NH4 +), the sugarcane's preferred source of N. Here, using a sugarcane bacterial artificial chromosome (BAC) library and a series of in silico analyses, we identified an AMMONIUM TRANSPORTER (AMT) from the AMT2 subfamily, sugarcane AMMONIUM TRANSPORTER 3;3 (ScAMT3;3), which is constitutively and highly expressed in young and mature leaves. To characterize its biochemical function, we ectopically expressed ScAMT3;3 in heterologous systems (Saccharomyces cerevisiae and Arabidopsis thaliana). The complementation of triple mep mutant yeast demonstrated that ScAMT3;3 is functional for NH3/H+ cotransport at high availability of NH4 + and under physiological pH conditions. The ectopic expression of ScAMT3;3 in the Arabidopsis quadruple AMT knockout mutant restored the transport capacity of 15N-NH4 + in roots and plant growth under specific N availability conditions, confirming the role of ScAMT3;3 in NH4 + transport in planta. Our results indicate that ScAMT3;3 belongs to the low-affinity transport system (Km 270.9 µM; Vmax 209.3 µmol g-1 root DW h-1). We were able to infer that ScAMT3;3 plays a presumed role in NH4 + source-sink remobilization in the shoots via phloem loading. These findings help to shed light on the functionality of a novel AMT2-type protein and provide bases for future research focusing on the improvement of sugarcane yield and N use efficiency.

Functional and comparative analysis of THI1 gene in grasses with a focus on sugarcane

De novo synthesis of thiamine (vitamin B1) in plants depends on the action of thiamine thiazole synthase, which synthesizes the thiazole ring, and is encoded by the THI1 gene. Here, we investigated the evolution and diversity of THI1 in Poaceae, where C4 and C3 photosynthetic plants co-evolved. An ancestral duplication of THI1 is observed in Panicoideae that remains in many modern monocots, including sugarcane. In addition to the two sugarcane copies (ScTHI1-1 and ScTHI1-2), we identified ScTHI1-2 alleles showing differences in their sequence, indicating divergence between ScTHI1-2a and ScTHI1-2b. Such variations are observed only in the Saccharum complex, corroborating the phylogeny. At least five THI1 genomic environments were found in Poaceae, two in sugarcane, M. sinensis, and S. bicolor. The THI1 promoter in Poaceae is highly conserved at 300 bp upstream of the start codon ATG and has cis-regulatory elements that putatively bind to transcription factors associated with development, growth, development and biological rhythms. An experiment set to compare gene expression levels in different tissues across the sugarcane R570 life cycle showed that ScTHI1-1 was expressed mainly in leaves regardless of age. Furthermore, ScTHI1 displayed relatively high expression levels in meristem and culm, which varied with the plant age. Finally, yeast complementation studies with THI4-defective strain demonstrate that only ScTHI1-1 and ScTHI1-2b isoforms can partially restore thiamine auxotrophy, albeit at a low frequency. Taken together, the present work supports the existence of multiple origins of THI1 harboring genomic regions in Poaceae with predicted functional redundancy. In addition, it questions the contribution of the levels of the thiazole ring in C4 photosynthetic plant tissues or potentially the relevance of the THI1 protein activity.

Redox homeostasis at SAM: a new role of HINT protein

ScHINT1 was identified at sugarcane SAM using subtractive libraries. Here, by bioinformatic tools, two-hybrid approach, and biochemical assays, we proposed that its role might be associated to control redox homeostasis. Such control is important for plant development and flowering transition, and this is ensured with some protein partners such as PAL and SBT that interact with ScHINT1. The shoot apical meristem transition from vegetative to reproductive is a crucial step for plants. In sugarcane (Saccharum spp.), this process is not well known, and it has an important impact on production due to field reduction. In view of this, ScHINT1 (Sugarcane HISTIDINE TRIAD NUCLEOTIDE-BINDING PROTEIN) was identified previously by subtractive cDNA libraries using Shoot Apical Meristem (SAM) by our group. This protein is a member of the HIT superfamily that was composed of hydrolase with an AMP site ligation. To better understand the role of ScHINT1 in sugarcane flowering, here its function in SAM was characterized using different approaches such as bioinformatics, two-hybrid assays, transgenic plants, and biochemical assays. ScHINT1 was conserved in plants, and it was grouped into four clades (HINT1, HINT2, HINT3, and HINT4). The 3D model proposed that ScHINT1 might be active as it was able to ligate to AMP subtract. Moreover, the two-hybrid approach identified two protein interactions: subtilase and phenylalanine ammonia-lyase. The evolutionary tree highlighted the relationships that each sequence has with specific subfamilies and different proteins. The 3D models constructed reveal structure conservation when compared with other PDB-related crystals, which indicates probable functional activity for the sugarcane models assessed. The interactome analysis showed a connection to different proteins that have antioxidative functions in apical meristems. Lastly, the transgenic plants with 35S::ScHINT1_AS (anti-sense orientation) produced more flowers than wild-type or 35S::ScHINT1_S (sense). Alpha-tocopherol and antioxidant enzymes measurement showed that their levels were higher in 35S::ScHINT_S plants than in 35S::ScHINT1_AS or wild-type plants. These results proposed that ScHINT1 might have an important role with other proteins in orchestrating this complex network for plant development and flowering.

Functional characterization of the sugarcane (Saccharum spp.) ammonium transporter AMT2;1 suggests a role in ammonium root-to-shoot translocation

AMMONIUM TRANSPORTER/METHYLAMMONIUM PERMEASE/RHESUS (AMT) family members transport ammonium across membranes in all life domains. Plant AMTs can be categorized into AMT1 and AMT2 subfamilies. Functional studies of AMTs, particularly AMT1-type, have been conducted using model plants but little is known about the function of AMTs from crops. Sugarcane (Saccharum spp.) is a major bioenergy crop that requires heavy nitrogen fertilization but depends on a low carbon-footprint for competitive sustainability. Here, we identified and functionally characterized sugarcane ScAMT2;1 by complementing ammonium uptake-defective mutants of Saccharomyces cerevisiae and Arabidopsis thaliana. Reporter gene driven by the ScAMT2;1 promoter in A. thaliana revealed preferential expression in the shoot vasculature and root endodermis/pericycle according to nitrogen availability and source. Arabidopsis quadruple mutant plants expressing ScAMT2;1 driven by the CaMV35S promoter or by a sugarcane endogenous promoter produced significantly more biomass than mutant plants when grown in NH₄ ⁺ and showed more ¹⁵N-ammonium uptake by roots and nitrogen translocation to shoots. In A. thaliana, ScAMT2;1 displayed a K_m of 90.17 µM and V_max of 338.99 µmoles h^-1 g^-1 root DW. Altogether, our results suggest that ScAMT2;1 is a functional high-affinity ammonium transporter that might contribute to ammonium uptake and presumably to root-to-shoot translocation under high NH₄ ⁺ conditions.

Recent Advances in Sugarcane Genomics, Physiology, and Phenomics for Superior Agronomic Traits

Advances in sugarcane breeding have contributed significantly to improvements in agronomic traits and crop yield. However, the growing global demand for sugar and biofuel in the context of climate change requires further improvements in cane and sugar yields. Attempts to achieve the desired rates of genetic gain in sugarcane by conventional breeding means are difficult as many agronomic traits are genetically complex and polygenic, with each gene exerting small effects. Unlike those of many other crops, the sugarcane genome is highly heterozygous due to its autopolyploid nature, which further hinders the development of a comprehensive genetic map. Despite these limitations, many superior agronomic traits/genes for higher cane yield, sugar production, and disease/pest resistance have been identified through the mapping of quantitative trait loci, genome-wide association studies, and transcriptome approaches. Improvements in traits controlled by one or two loci are relatively easy to achieve; however, this is not the case for traits governed by many genes. Many desirable phenotypic traits are controlled by quantitative trait nucleotides (QTNs) with small and variable effects. Assembling these desired QTNs by conventional breeding methods is time consuming and inefficient due to genetic drift. However, recent developments in genomics selection (GS) have allowed sugarcane researchers to select and accumulate desirable alleles imparting superior traits as GS is based on genomic estimated breeding values, which substantially increases the selection efficiency and genetic gain in sugarcane breeding programs. Next-generation sequencing techniques coupled with genome-editing technologies have provided new vistas in harnessing the sugarcane genome to look for desirable agronomic traits such as erect canopy, leaf angle, prolonged greening, high biomass, deep root system, and the non-flowering nature of the crop. Many desirable cane-yielding traits, such as single cane weight, numbers of tillers, numbers of millable canes, as well as cane quality traits, such as sucrose and sugar yield, have been explored using these recent biotechnological tools. This review will focus on the recent advances in sugarcane genomics related to genetic gain and the identification of favorable alleles for superior agronomic traits for further utilization in sugarcane breeding programs.

Bioinformatic analyses to uncover genes involved in trehalose metabolism in the polyploid sugarcane

Trehalose-6-phosphate (T6P) is an intermediate of trehalose biosynthesis that plays an essential role in plant metabolism and development. Here, we comprehensively analyzed sequences from enzymes of trehalose metabolism in sugarcane, one of the main crops used for bioenergy production. We identified protein domains, phylogeny, and in silico expression levels for all classes of enzymes. However, post-translational modifications and residues involved in catalysis and substrate binding were analyzed only in trehalose-6-phosphate synthase (TPS) sequences. We retrieved 71 putative full-length TPS, 93 trehalose-6-phosphate phosphatase (TPP), and 3 trehalase (TRE) of sugarcane, showing all their conserved domains, respectively. Putative TPS (Classes I and II) and TPP sugarcane sequences were categorized into well-known groups reported in the literature. We measured the expression levels of the sequences from one sugarcane leaf transcriptomic dataset. Furthermore, TPS Class I has specific N-glycosylation sites inserted in conserved motifs and carries catalytic and binding residues in its TPS domain. Some of these residues are mutated in TPS Class II members, which implies loss of enzyme activity. Our approach retrieved many homo(eo)logous sequences for genes involved in trehalose metabolism, paving the way to discover the role of T6P signaling in sugarcane.

OsAT1, an anion transporter, negatively regulates grain size and yield in rice

Improving the grain yield of rice is a central goal of basic and applied scientific research. Here, we identified an anion transporter, OsAT1, localized in the endoplasmic reticulum and Golgi. OsAT1 is highly expressed in flag, stem, and sheath as monitored using qRT-PCR and pOsAT1::GUS. Thousand-grain weight, grain weight per plant, and content of starch were significantly increased in OsAT1 knock-down mutants (OsAT1-Ri) but significantly decreased in OsAT1 overexpressed lines (OsAT1-OE). In addition, the grain weight per plant increased by 6.17% to 6.78% in OsAT1-RNAi lines, whereas it decreased by 45.93% to 46.76% in OsAT1-OE lines, compared to wild-type. Moreover, the copper content was noticeably reduced in flag leaf of OsAT1-Ri lines and increased in OsAT1-OE lines. RNA-sequencing analysis of OsAT1-OE lines revealed that the genes related to starch biosynthesis and metabolism pathway were enriched in the down-regulated category. Thus, our results suggest that knock-down of OsAT1 in rice possibly reduces copper accumulation and improves the accumulation of storage starch, hence, increasing the grain size and weight. OsAT1 may be a useful gene to consider for cereal breeding programs.

Contributions of sugar transporters to crop yield and fruit quality

The flux, distribution, and storage of soluble sugars regulate crop yield in terms of starch, oil, protein, and total carbohydrates, and affect the quality of many horticultural products. Sugar transporters contribute to phloem loading and unloading. The mechanisms of phloem loading have been studied in detail, but the complex and diverse mechanisms of phloem unloading and sugar storage in sink organs are less explored. Unloading and subsequent transport mechanisms for carbohydrates vary in different sink organs. Analyzing the transport and storage mechanisms of carbohydrates in important storage organs, such as cereal seeds, fruits, or stems of sugarcane, will provide information for genetic improvements to increase crop yield and fruit quality. This review discusses current research progress on sugar transporters involved in carbohydrate unloading and storage in sink organs. The roles of sugar transporters in crop yield and the accumulation of sugars are also discussed to highlight their contribution to efficient breeding.

A comprehensive molecular cytogenetic analysis of the genome architecture in modern sugarcane cultivars

Modern sugarcane cultivars are derived from the hybridization of Saccharum officinarum (2n = 80) and S. spontaneum (2n = 40-128), leading to a variety of complex genomes with highly polyploid and varied chromosome structures. These complex genomes have hindered deciphering the genome structure and marker-assisted selection in sugarcane breeding. Ten cultivars were analyzed by fluorescence in situ hybridization adopting chromosome painting and S. spontaneum-specific probes. The results showed six types of chromosomes in the studied cultivars, including S. spontaneum or S. officinarum chromosomes, interspecific recombinations from homoeologous or nonhomoeologous chromosomes, and translocations of S. spontaneum or S. officinarum chromosomes. The results showed unexpectedly high proportions of interspecific recombinations in these cultivars (11.9-40.9%), which renew our knowledge that less than 13% of chromosomes result from interspecific exchanges. Also, the results showed a high frequency of translocations (an average of 2.15 translocations per chromosome) between S. officinarum chromosomes. The diverse types of chromosomes in cultivars imply that hybrid gametes of S. spontaneum and S. officinarum may form unusual chromosome pairs, including homoeologous or nonhomoeologous chromosomes either between or within S. spontaneum and S. officinarum. Moreover, we consistently observed 11 or 12 copies for the four studied chromosomes, i.e., chromosomes 1, 2, 7, and 8, suggesting steady transmission during the breeding program. By comparison, we found a relatively fewer copies of S. spontaneum chromosome 1 than those of S. spontaneum chromosomes 2, 7, and 8. These results provide deep insights into the structure of cultivars and may facilitate chromosome-assisted selection in sugarcane breeding.

A Sugarcane G-Protein-Coupled Receptor, ShGPCR1, Confers Tolerance to Multiple Abiotic Stresses

Sugarcane (Saccharum spp.) is a prominent source of sugar and serves as bioenergy/biomass feedstock globally. Multiple biotic and abiotic stresses, including drought, salinity, and cold, adversely affect sugarcane yield. G-protein-coupled receptors (GPCRs) are components of G-protein-mediated signaling affecting plant growth, development, and stress responses. Here, we identified a GPCR-like protein (ShGPCR1) from sugarcane and energy cane (Saccharum spp. hybrids) and characterized its function in conferring tolerance to multiple abiotic stresses. ShGPCR1 protein sequence contained nine predicted transmembrane (TM) domains connected by four extracellular and four intracellular loops, which could interact with various ligands and heterotrimeric G proteins in the cells. ShGPCR1 sequence displayed other signature features of a GPCR, such as a putative guanidine triphosphate (GTP)-binding domain, as well as multiple myristoylation and protein phosphorylation sites, presumably important for its biochemical function. Expression of ShGPCR1 was upregulated by drought, salinity, and cold stresses. Subcellular imaging and calcium (Ca²⁺) measurements revealed that ShGPCR1 predominantly localized to the plasma membrane and enhanced intracellular Ca²⁺ levels in response to GTP, respectively. Furthermore, constitutive overexpression of ShGPCR1 in sugarcane conferred tolerance to the three stressors. The stress-tolerance phenotype of the transgenic lines corresponded with activation of multiple drought-, salinity-, and cold-stress marker genes, such as Saccharum spp. LATE EMBRYOGENESIS ABUNDANT, DEHYDRIN, DROUGHT RESPONSIVE 4, GALACTINOL SYNTHASE, ETHYLENE RESPONSIVE FACTOR 3, SALT OVERLY SENSITIVE 1, VACUOLAR Na+/H+ ANTIPORTER 1, NAM/ATAF1/2/CUC2, COLD RESPONSIVE FACTOR 2, and ALCOHOL DEHYDROGENASE 3. We suggest that ShGPCR1 plays a key role in conferring tolerance to multiple abiotic stresses, and the engineered lines may be useful to enhance sugarcane production in marginal environments with fewer resources.

Characterization of full-length transcriptome in Saccharum officinarum and molecular insights into tiller development

BACKGROUND

Although extensive breeding efforts are ongoing in sugarcane (Saccharum officinarum L.), the average yield is far below the theoretical potential. Tillering is an important component of sugarcane yield, however, the molecular mechanism underlying tiller development is still elusive. The limited genomic data in sugarcane, particularly due to its complex and large genome, has hindered in-depth molecular studies.

RESULTS

Herein, we generated full-length (FL) transcriptome from developing leaf and tiller bud samples based on PacBio Iso-Seq. In addition, we performed RNA-seq from tiller bud samples at three developmental stages (T0, T1 and T2) to uncover key genes and biological pathways involved in sugarcane tiller development. In total, 30,360 and 20,088 high-quality non-redundant isoforms were identified in leaf and tiller bud samples, respectively, representing 41,109 unique isoforms in sugarcane. Likewise, we identified 1063 and 1037 alternative splicing events identified in leaf and tiller bud samples, respectively. We predicted the presence of coding sequence for 40,343 isoforms, 98% of which was successfully annotated. Comparison with previous FL transcriptomes in sugarcane revealed 2963 unreported isoforms. In addition, we characterized 14,946 SSRs from 11,700 transcripts and 310 lncRNAs. By integrating RNA-seq with the FL transcriptome, 468 and 57 differentially expressed genes (DEG) were identified in T1vsT0 and T2vsT0, respectively. Strong up-regulation of several pyruvate phosphate dikinase and phosphoenolpyruvate carboxylase genes suggests enhanced carbon fixation and protein synthesis to facilitate tiller growth. Similarly, up-regulation of linoleate 9S-lipoxygenase and lipoxygenase genes in the linoleic acid metabolism pathway suggests high synthesis of key oxylipins involved in tiller growth and development.

CONCLUSIONS

Collectively, we have enriched the genomic data available in sugarcane and provided candidate genes for manipulating tiller formation and development, towards productivity enhancement in sugarcane.

Three founding ancestral genomes involved in the origin of sugarcane

BACKGROUND AND AIMS

Modern sugarcane cultivars (Saccharum spp.) are high polyploids, aneuploids (2n = ~12x = ~120) derived from interspecific hybridizations between the domesticated sweet species Saccharum officinarum and the wild species S. spontaneum.

METHODS

To analyse the architecture and origin of such a complex genome, we analysed the sequences of all 12 hom(oe)ologous haplotypes (BAC clones) from two distinct genomic regions of a typical modern cultivar, as well as the corresponding sequence in Miscanthus sinense and Sorghum bicolor, and monitored their distribution among representatives of the Saccharum genus.

KEY RESULTS

The diversity observed among haplotypes suggested the existence of three founding genomes (A, B, C) in modern cultivars, which diverged between 0.8 and 1.3 Mya. Two genomes (A, B) were contributed by S. officinarum; these were also found in its wild presumed ancestor S. robustum, and one genome (C) was contributed by S. spontaneum. These results suggest that S. officinarum and S. robustum are derived from interspecific hybridization between two unknown ancestors (A and B genomes). The A genome contributed most haplotypes (nine or ten) while the B and C genomes contributed one or two haplotypes in the regions analysed of this typical modern cultivar. Interspecific hybridizations likely involved accessions or gametes with distinct ploidy levels and/or were followed by a series of backcrosses with the A genome. The three founding genomes were found in all S. barberi, S. sinense and modern cultivars analysed. None of the analysed accessions contained only the A genome or the B genome, suggesting that representatives of these founding genomes remain to be discovered.

CONCLUSIONS

This evolutionary model, which combines interspecificity and high polyploidy, can explain the variable chromosome pairing affinity observed in Saccharum. It represents a major revision of the understanding of Saccharum diversity.

Multiallelic, Targeted Mutagenesis of Magnesium Chelatase With CRISPR/Cas9 Provides a Rapidly Scorable Phenotype in Highly Polyploid Sugarcane

Genome editing with sequence-specific nucleases, such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), is revolutionizing crop improvement. Developing efficient genome-editing protocols for highly polyploid crops, including sugarcane (x = 10-13), remains challenging due to the high level of genetic redundancy in these plants. Here, we report the efficient multiallelic editing of magnesium chelatase subunit I (MgCh) in sugarcane. Magnesium chelatase is a key enzyme for chlorophyll biosynthesis. CRISPR/Cas9-mediated targeted co-mutagenesis of 49 copies/alleles of magnesium chelatase was confirmed via Sanger sequencing of cloned PCR amplicons. This resulted in severely reduced chlorophyll contents, which was scorable at the time of plant regeneration in the tissue culture. Heat treatment following the delivery of genome editing reagents elevated the editing frequency 2-fold and drastically promoted co-editing of multiple alleles, which proved necessary to create a phenotype that was visibly distinguishable from the wild type. Despite their yellow leaf color, the edited plants were established well in the soil and did not show noticeable growth retardation. This approach will facilitate the establishment of genome editing protocols for recalcitrant crops and support further optimization, including the evaluation of alternative RNA-guided nucleases to overcome the limitations of the protospacer adjacent motif (PAM) site or to develop novel delivery strategies for genome editing reagents.

Role of miRNAs in the host-pathogen interaction between sugarcane and Colletotrichum falcatum, the red rot pathogen

KEY MESSAGE

Sugarcane microRNAs specifically involved during compatible and incompatible interactions with red rot pathogen Colletotrichum falcatum were identified. We have identified how the miRNAs regulate their gene targets and elaborated evidently on the underlying molecular mechanism of sugarcane defense response to C. falcatum for the first time. Resistance against the fungal pathogen Colletotrichum falcatum causing red rot is one of the most desirable traits for sustainable crop cultivation in sugarcane. To gain new insight into the host defense mechanism against C. falcatum, we studied the role of sugarcane microRNAs during compatible and incompatible interactions by adopting the NGS platform. We have sequenced a total of 80 miRNA families that comprised 980 miRNAs, and the putative targets of the miRNAs include transcription factors, membrane-bound proteins, glutamate receptor proteins, lignin biosynthesis proteins, signaling cascade proteins, transporter proteins, mitochondrial proteins, ER proteins, defense-related, stress response proteins, translational regulation proteins, cell proliferation, and ubiquitination proteins. Further, qRT-PCR analyses of 8 differentially regulated miRNAs and 26 gene transcript targets expression indicated that these miRNAs have a regulatory effect on the expression of respective target genes in most of the cases. Also, the results suggest that certain miRNA regulates many target genes that are involved in inciting early responses to the pathogen infection, signaling pathways, endoplasmic reticulum stress, and resistance gene activation through feedback response from various cellular processes during the compatible and incompatible interaction with the red rot pathogen C. falcatum. The present study revealed the role of sugarcane miRNAs and their target genes during sugarcane-C. falcatum interaction and provided new insight into the miRNA-mediated defense mechanism in sugarcane for the first time.

Error-free recombination in sugarcane mediated by only 30 nucleotides of homology and CRISPR/Cas9 induced DNA breaks or Cre-recombinase

Precision genome editing by homology directed repair has tremendous potential for crop improvement. This study describes in planta homologous recombination mediated by CRISPR/Cas9 induced DNA double strand break in proximity to a single short (∼30 nt) homology arm. The efficiency of CRISPR/Cas9-mediated recombination between two loxP sites was compared with Cre (Cyclization recombination enzyme) and codon-optimized Cre-mediated site-specific recombination in sugarcane. A transgenic locus was generated with a selectable nptII coding sequence with terminator between two loxP sites located downstream of a constitutive promoter and acting as transcription block for the downstream promoter-less gusA coding sequence with terminator. Recombination between the two loxP sites resulted in deletion of the transcription block and restored gus activity. This transgenic locus provided an efficient screen for identification of recombination events in sugarcane callus following biolistic delivery of Cre, codon-optimized Cre, or the combination of sgRNA and Cas9 targeting the 5' loxP site. The Cre codon optimized for sugarcane displayed the highest efficiency in mediating the recombination that restored gus activity followed by cre and CRISPR/Cas9. Remarkably the short region of homology of the loxP site cleaved by Cas9 (30 nt)-mediated error-free recombination in all 21 events from three different experiments that were analyzed by Sanger sequencing consistent with homology directed repair. These findings will inform rational design of strategies for precision genome editing in plants.

Sucrose-phosphate phosphatase from sugarcane reveals an ancestral tandem duplication

BACKGROUND

Sugarcane is capable to store large amounts of sucrose in the culm at maturity hence it became a major source of sucrose for the food and the renewable energy industries. Sucrose, the main disaccharide produced by photosynthesis, is mainly stored in the vacuole of the cells of non-photosynthetic tissues. Two pathways are known to release free sucrose in plant cells, one is de novo synthesis dependent on sucrose phosphate synthase (SPS) and sucrose phosphate phosphatase (S6PP) while the other is regulatory and dependent on sucrose synthase (SuSy) activity. The molecular understanding of genes that give rise to the expression of the enzyme sucrose phosphate phosphatase, responsible for the release of sucrose in the last synthetic step lag behind the regulatory SuSy gene.

RESULTS

Sugarcane genome sequencing effort disclosed the existence of a tandem duplication and the present work further support that both S6PP.1 and S6PP_2D isoforms are actively transcribed in young sugarcane plants but significantly less at maturity. Two commercial hybrids (SP80-3280 and R570) and both Saccharum spontaneum (IN84-58) and S.officinarum (BADILLA) exhibit transcriptional activity at three-month-old plants of the tandem S6PP_2D in leaves, culm, meristem and root system with a cultivar-specific distribution. Moreover, this tandem duplication is shared with other grasses and is ancestral in the group.

CONCLUSION

Detection of a new isoform of S6PP resulting from the translation of 14 exon-containing transcript (S6PP_2D) will contribute to the knowledge of sucrose metabolism in plants. In addition, expression varies along plant development and between sugarcane cultivars and parental species.

Multiallelic, Targeted Mutagenesis of Magnesium Chelatase With CRISPR/Cas9 Provides a Rapidly Scorable Phenotype in Highly Polyploid Sugarcane

Genome editing with sequence-specific nucleases, such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), is revolutionizing crop improvement. Developing efficient genome-editing protocols for highly polyploid crops, including sugarcane (x = 10-13), remains challenging due to the high level of genetic redundancy in these plants. Here, we report the efficient multiallelic editing of magnesium chelatase subunit I (MgCh) in sugarcane. Magnesium chelatase is a key enzyme for chlorophyll biosynthesis. CRISPR/Cas9-mediated targeted co-mutagenesis of 49 copies/alleles of magnesium chelatase was confirmed via Sanger sequencing of cloned PCR amplicons. This resulted in severely reduced chlorophyll contents, which was scorable at the time of plant regeneration in the tissue culture. Heat treatment following the delivery of genome editing reagents elevated the editing frequency 2-fold and drastically promoted co-editing of multiple alleles, which proved necessary to create a phenotype that was visibly distinguishable from the wild type. Despite their yellow leaf color, the edited plants were established well in the soil and did not show noticeable growth retardation. This approach will facilitate the establishment of genome editing protocols for recalcitrant crops and support further optimization, including the evaluation of alternative RNA-guided nucleases to overcome the limitations of the protospacer adjacent motif (PAM) site or to develop novel delivery strategies for genome editing reagents.

Applying Molecular Phenotyping Tools to Explore Sugarcane Carbon Potential

Sugarcane (Saccharum spp.), a C₄ grass, has a peculiar feature: it accumulates, gradient-wise, large amounts of carbon (C) as sucrose in its culms through a complex pathway. Apart from being a sustainable crop concerning C efficiency and bioenergetic yield per hectare, sugarcane is used as feedstock for producing ethanol, sugar, high-value compounds, and products (e.g., polymers and succinate), and bioelectricity, earning the title of the world's leading biomass crop. Commercial cultivars, hybrids bearing high levels of polyploidy, and aneuploidy, are selected from a large number of crosses among suitable parental genotypes followed by the cloning of superior individuals among the progeny. Traditionally, these classical breeding strategies have been favoring the selection of cultivars with high sucrose content and resistance to environmental stresses. A current paradigm change in sugarcane breeding programs aims to alter the balance of C partitioning as a means to provide more plasticity in the sustainable use of this biomass for metabolic engineering and green chemistry. The recently available sugarcane genetic assemblies powered by data science provide exciting perspectives to increase biomass, as the current sugarcane yield is roughly 20% of its predicted potential. Nowadays, several molecular phenotyping tools can be applied to meet the predicted sugarcane C potential, mainly targeting two competing pathways: sucrose production/storage and biomass accumulation. Here we discuss how molecular phenotyping can be a powerful tool to assist breeding programs and which strategies could be adopted depending on the desired final products. We also tackle the advances in genetic markers and mapping as well as how functional genomics and genetic transformation might be able to improve yield and saccharification rates. Finally, we review how "omics" advances are promising to speed up plant breeding and reach the unexplored potential of sugarcane in terms of sucrose and biomass production.

Sugarcane mosaic virus mediated changes in cytosine methylation pattern and differentially transcribed fragments in resistance-contrasting sugarcane genotypes

Sugarcane mosaic virus (SCMV) is the causal agent of sugarcane mosaic disease (SMD) in Brazil; it is mainly controlled by using resistant cultivars. Studies on the changes in sugarcane transcriptome provided the first insights about the molecular basis underlying the genetic resistance to SMD; nonetheless, epigenetic modifications such as cytosine methylation is also informative, considering its roles in gene expression regulation. In our previous study, differentially transcribed fragments (DTFs) were obtained using cDNA-amplified fragment length polymorphism by comparing mock- and SCMV-inoculated plants from two sugarcane cultivars with contrasting responses to SMD. In this study, the identification of unexplored DTFs was continued while the same leaf samples were used to evaluate SCMV-mediated changes in the cytosine methylation pattern by using methylation-sensitive amplification polymorphism. This analysis revealed minor changes in cytosine methylation in response to SCMV infection, but distinct changes between the cultivars with contrasting responses to SMD, with higher hypomethylation events 24 and 72 h post-inoculation in the resistant cultivar. The differentially methylated fragments (DMFs) aligned with transcripts, putative promoters, and genomic regions, with a preponderant distribution within CpG islands. The transcripts found were associated with plant immunity and other stress responses, epigenetic changes, and transposable elements. The DTFs aligned with transcripts assigned to stress responses, epigenetic changes, photosynthesis, lipid transport, and oxidoreductases, in which the transcriptional start site is located in proximity with CpG islands and tandem repeats. Real-time quantitative polymerase chain reaction results revealed significant upregulation in the resistant cultivar of aspartyl protease and VQ protein, respectively, selected from DMF and DTF alignments, suggesting their roles in genetic resistance to SMD and supporting the influence of cytosine methylation in gene expression. Thus, we identified new candidate genes for further validation and showed that the changes in cytosine methylation may regulate important mechanisms underlying the genetic resistance to SMD.

The sugarcane ShMYB78 transcription factor activates suberin biosynthesis in Nicotiana benthamiana

KEY MESSAGE

A sugarcane MYB present in the culm induces suberin biosynthesis and is involved both with fatty acid and phenolics metabolism. Few transcription factors have been described as regulators of cell wall polymers deposition in C4 grasses. Particularly, regulation of suberin biosynthesis in this group of plants remains poorly understood. Here, we showed that the sugarcane MYB transcription factor ShMYB78 is an activator of suberin biosynthesis and deposition. ShMYB78 was identified upon screening genes whose expression was upregulated in sugarcane internodes undergoing suberization during culm development or triggered by wounding. Agrobacterium-mediated transient expression of ShMYB78 in Nicotiana benthamiana leaves induced the ectopic deposition of suberin and its aliphatic and aromatic monomers. Further, the expression of suberin-related genes was induced by ShMYB78 heterologous expression in Nicotiana benthamiana leaves. ShMYB78 was shown to be a nuclear protein based on its presence in sugarcane internode nuclear protein extracts, and protoplast transactivation assays demonstrated that ShMYB78 activates the promoters of the sugarcane suberin biosynthetic genes β-ketoacyl-CoA synthase (ShKCS20) and caffeic acid-O-methyltransferase (ShCOMT). Our results suggest that ShMYB78 may be involved in the transcriptional regulation of suberin deposition, from fatty acid metabolism to phenylpropanoid biosynthesis, in sugarcane internodes.

Differential expression in leaves of Saccharum genotypes contrasting in biomass production provides evidence of genes involved in carbon partitioning

BACKGROUND

The development of biomass crops aims to meet industrial yield demands, in order to optimize profitability and sustainability. Achieving these goals in an energy crop like sugarcane relies on breeding for sucrose accumulation, fiber content and stalk number. To expand the understanding of the biological pathways related to these traits, we evaluated gene expression of two groups of genotypes contrasting in biomass composition.

RESULTS

First visible dewlap leaves were collected from 12 genotypes, six per group, to perform RNA-Seq. We found a high number of differentially expressed genes, showing how hybridization in a complex polyploid system caused extensive modifications in genome functioning. We found evidence that differences in transposition and defense related genes may arise due to the complex nature of the polyploid Saccharum genomes. Genotypes within both biomass groups showed substantial variability in genes involved in photosynthesis. However, most genes coding for photosystem components or those coding for phosphoenolpyruvate carboxylases (PEPCs) were upregulated in the high biomass group. Sucrose synthase (SuSy) coding genes were upregulated in the low biomass group, showing that this enzyme class can be involved with sucrose synthesis in leaves, similarly to sucrose phosphate synthase (SPS) and sucrose phosphate phosphatase (SPP). Genes in pathways related to biosynthesis of cell wall components and expansins coding genes showed low average expression levels and were mostly upregulated in the high biomass group.

CONCLUSIONS

Together, these results show differences in carbohydrate synthesis and carbon partitioning in the source tissue of distinct phenotypic groups. Our data from sugarcane leaves revealed how hybridization in a complex polyploid system resulted in noticeably different transcriptomic profiles between contrasting genotypes.

Elucidating the Genetic Architecture of Fiber Quality in Hemp (Cannabis sativa L.) Using a Genome-Wide Association Study

Hemp (Cannabis sativa L.) is a bast-fiber crop with a great potential in the emerging bio-based economy. Yet, hemp breeding for fiber quality is restricted and that is mainly due to the limited knowledge of the genetic architecture of its fiber quality. A panel of 123 hemp accessions, with large phenotypic variability, was used to study the genetic basis of seven cell wall and bast fiber traits relevant to fiber quality. These traits showed large genetic variance components and high values of broad sense heritability in this hemp panel, as concluded from the phenotypic evaluation across three test locations with contrasting environments. The hemp panel was genotyped using restriction site associated DNA sequencing (RAD-seq). Subsequently, a large set (> 600,000) of selected genome-wide single nucleotide polymorphism (SNP) markers was used for a genome-wide association study (GWAS) approach to get insights into quantitative trait loci (QTLs) controlling fiber quality traits. In absence of a complete hemp genome sequence, identification of QTLs was based on the following characteristics: (i) association level to traits, (ii) fraction of explained trait variance, (iii) collinearity between QTLs, and (iv) detection across different environments. Using this approach, 16 QTLs were identified across locations for different fiber quality traits, including contents of glucose, glucuronic acid, mannose, xylose, lignin, and bast fiber content. Among them, six were found across the three environments. The genetic markers composing the QTLs that are common across locations are valuable tools to develop novel genotypes of hemp with improved fiber quality. Underneath the QTLs, 12 candidate genes were identified which are likely to be involved in the biosynthesis and modification of monosaccharides, polysaccharides, and lignin. These candidate genes were suggested to play an important role in determining fiber quality in hemp. This study provides new insights into the genetic architecture of fiber traits, identifies QTLs and candidate genes that form the basis for molecular breeding for high fiber quality hemp cultivars.

Time-series expression profiling of sugarcane leaves infected with Puccinia kuehnii reveals an ineffective defense system leading to susceptibility

Successful orange rust development on sugarcane can potentially be explained as suppression of the plant immune system by the pathogen or delayed plant signaling to trigger defense responses. Puccinia kuehnii is an obligate biotrophic fungus that infects sugarcane leaves causing a disease called orange rust. It spread out to other countries resulting in reduction of crop yield since its first outbreak. One of the knowledge gaps of that pathosystem is to understand the molecular mechanisms altered in susceptible plants by this biotic stress. Here, we investigated the changes in temporal expression of transcripts in pathways associated with the immune system. To achieve this purpose, we used RNA-Seq to analyze infected leaf samples collected at five time points after inoculation. Differential expression analyses of adjacent time points revealed substantial changes at 12, 48 h after inoculation and 12 days after inoculation, coinciding with the events of spore germination, haustoria post-penetration and post-sporulation, respectively. During the first 24 h, a lack of transcripts involved with resistance mechanisms was revealed by underrepresentation of hypersensitive and defense response related genes. However, two days after inoculation, upregulation of genes involved with immune response regulation provided evidence of some potential defense response. Events related to biotic stress responses were predominantly downregulated in the initial time points, but expression was later restored to basal levels. Genes involved in carbohydrate metabolism showed evidence of repression followed by upregulation, possibly to ensure the pathogen nutritional requirements were met. Our results support the hypothesis that P. kuehnii initially suppressed sugarcane genes involved in plant defense systems. Late overexpression of specific regulatory pathways also suggests the possibility of an inefficient recognition system by a susceptible sugarcane genotype.

Role of NGS and SNP genotyping methods in sugarcane improvement programs

Sugarcane (Saccharum spp.) is one of the most economically significant crops because of its high sucrose content and it is a promising biomass feedstock for biofuel production. Sugarcane genome sequencing and analysis is a difficult task due to its heterozygosity and polyploidy. Long sequence read technologies, PacBio Single-Molecule Real-Time (SMRT) sequencing, the Illumina TruSeq, and the Oxford Nanopore sequencing could solve the problem of genome assembly. On the applications side, next generation sequencing (NGS) technologies played a major role in the discovery of single nucleotide polymorphism (SNP) and the development of low to high throughput genotyping platforms. The two mainstream high throughput genotyping platforms are the SNP microarray and genotyping by sequencing (GBS). This paper reviews the NGS in sugarcane genomics, genotyping methodologies, and the choice of these methods. Array-based SNP genotyping is robust, provides consistent SNPs, and relatively easier downstream data analysis. The GBS method identifies large scale SNPs across the germplasm. A combination of targeted GBS and array-based genotyping methods should be used to increase the accuracy of genomic selection and marker-assisted breeding.

Alternative Splicing of Circadian Clock Genes Correlates With Temperature in Field-Grown Sugarcane

Alternative Splicing (AS) is a mechanism that generates different mature transcripts from precursor mRNAs (pre-mRNAs) of the same gene. In plants, a wide range of physiological and metabolic events are related to AS, as well as fast responses to changes in temperature. AS is present in around 60% of intron-containing genes in Arabidopsis, 46% in rice, and 38% in maize and it is widespread among the circadian clock genes. Little is known about how AS influences the circadian clock of C4 plants, like commercial sugarcane, a C4 crop with a complex hybrid genome. This work aims to test if the daily dynamics of AS forms of circadian clock genes are regulated by environmental factors, such as temperature, in the field. A systematic search for AS in five sugarcane clock genes, ScLHY, ScPRR37, ScPRR73, ScPRR95, and ScTOC1 using different organs of sugarcane sampled during winter, with 4 months old plants, and during summer, with 9 months old plants, revealed temperature- and organ-dependent expression of at least one alternatively spliced isoform in all genes. Expression of AS isoforms varied according to the season. Our results suggest that AS events in circadian clock genes are correlated with temperature.

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.

Genomic resources for energy cane breeding in the post genomics era

Sugarcane is one of the most sustainable energy crops among cultivated crops presenting the highest tonnage of cultivated plants. Its high productivity of sugar, bioethanol and bioelectricity make it a promising green alternative to petroleum. Furthermore, the myriad of products that can be derived from sugarcane biomass has been driving breeding programs towards varieties with a higher yield of fiber and a more vigorous and sustainable performance: the energy cane. Here we provide an overview of the energy cane including plant description, breeding efforts, types, and end-uses. In addition, we describe recently published genomic resources for the development of this crop, discuss current knowledge of cell wall metabolism, bioinformatic tools and databases available for the community.

Sugarcane Omics: An Update on the Current Status of Research and Crop Improvement

Sugarcane is an important crop from Poaceae family, contributing about 80% of the total world's sucrose with an annual value of around US$150 billion. In addition, sugarcane is utilized as a raw material for the production of bioethanol, which is an alternate source of renewable energy. Moving towards sugarcane omics, a remarkable success has been achieved in gene transfer from a wide variety of plant and non-plant sources to sugarcane, with the accessibility of efficient transformation systems, selectable marker genes, and genetic engineering gears. Genetic engineering techniques make possible to clone and characterize useful genes and also to improve commercially important traits in elite sugarcane clones that subsequently lead to the development of an ideal cultivar. Sugarcane is a complex polyploidy crop, and hence no single technique has been found to be the best for the confirmation of polygenic and phenotypic characteristics. To better understand the application of basic omics in sugarcane regarding agronomic characters and industrial quality traits as well as responses to diverse biotic and abiotic stresses, it is important to explore the physiology, genome structure, functional integrity, and collinearity of sugarcane with other more or less similar crops/plants. Genetic improvements in this crop are hampered by its complex genome, low fertility ratio, longer production cycle, and susceptibility to several biotic and abiotic stresses. Biotechnology interventions are expected to pave the way for addressing these obstacles and improving sugarcane crop. Thus, this review article highlights up to date information with respect to how advanced data of omics (genomics, transcriptomic, proteomics and metabolomics) can be employed to improve sugarcane crops.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

A genome-wide association study identified loci for yield component traits in sugarcane (Saccharum spp.)

Sugarcane (Saccharum spp.) has a complex genome with variable ploidy and frequent aneuploidy, which hampers the understanding of phenotype and genotype relations. Despite this complexity, genome-wide association studies (GWAS) may be used to identify favorable alleles for target traits in core collections and then assist breeders in better managing crosses and selecting superior genotypes in breeding populations. Therefore, in the present study, we used a diversity panel of sugarcane, called the Brazilian Panel of Sugarcane Genotypes (BPSG), with the following objectives: (i) estimate, through a mixed model, the adjusted means and genetic parameters of the five yield traits evaluated over two harvest years; (ii) detect population structure, linkage disequilibrium (LD) and genetic diversity using simple sequence repeat (SSR) markers; (iii) perform GWAS analysis to identify marker-trait associations (MTAs); and iv) annotate the sequences giving rise to SSR markers that had fragments associated with target traits to search for putative candidate genes. The phenotypic data analysis showed that the broad-sense heritability values were above 0.48 and 0.49 for the first and second harvests, respectively. The set of 100 SSR markers produced 1,483 fragments, of which 99.5% were polymorphic. These SSR fragments were useful to estimate the most likely number of subpopulations, found to be four, and the LD in BPSG, which was stronger in the first 15 cM and present to a large extension (65 cM). Genetic diversity analysis showed that, in general, the clustering of accessions within the subpopulations was in accordance with the pedigree information. GWAS performed through a multilocus mixed model revealed 23 MTAs, six, three, seven, four and three for soluble solid content, stalk height, stalk number, stalk weight and cane yield traits, respectively. These MTAs may be validated in other populations to support sugarcane breeding programs with introgression of favorable alleles and marker-assisted selection.

Genome-wide alternative splicing landscapes modulated by biotrophic sugarcane smut pathogen

Alternative splicing (AS) promotes transcriptome and proteome diversity during growth, development, and stress responses in eukaryotes. Genome-wide studies of AS in sugarcane (Saccharum spp.) are lacking, mainly due to the absence of a high-quality sequenced reference genome, sugarcane's large, complex genome, and the variable chromosome numbers and polyploidy of sugarcane cultivars. Here, we analyzed changes in the sugarcane isoform-level transcriptome and AS landscape during infection with the smut fungus (Sporisorium scitamineum) using a hybrid approach involving Sorghum bicolor reference-based and Trinity de novo mapping tools. In total, this analysis detected 16,039 and 15,379 transcripts (≥2 FPKM) at 5 and 200 days after infection, respectively. A conservative estimate of isoform-level expression suggested that approximately 5,000 (14%) sugarcane genes undergo AS. Differential expression analysis of the alternatively spliced genes in healthy and smut-infected sugarcane revealed 896 AS events modulated at different stages of infection. Gene family and gene ontology functional enrichment analysis of the differentially spliced genes revealed overrepresentation of functional categories related to the cell wall, defense, and redox homeostasis pathways. Our study provides novel insight into the AS landscape of sugarcane during smut disease interactions.

The population genetic structure approach adds new insights into the evolution of plant LTR retrotransposon lineages

Long terminal repeat retrotransposons (LTR-RTs) in plant genomes differ in abundance, structure and genomic distribution, reflecting the large number of evolutionary lineages. Elements within lineages can be considered populations, in which each element is an individual in its genomic environment. In this way, it would be reasonable to apply microevolutionary analyses to understand transposable element (TE) evolution, such as those used to study the genetic structure of natural populations. Here, we applied a Bayesian method to infer genetic structure of populations together with classical phylogenetic and dating tools to analyze LTR-RT evolution using the monocot Setaria italica as a model species. In contrast to a phylogeny, the Bayesian clusterization method identifies populations by assigning individuals to one or more clusters according to the most probabilistic scenario of admixture, based on genetic diversity patterns. In this work, each LTR-RT insertion was considered to be one individual and each LTR-RT lineage was considered to be a single species. Nine evolutionary lineages of LTR-RTs were identified in the S. italica genome that had different genetic structures with variable numbers of clusters and levels of admixture. Comprehensive analysis of the phylogenetic, clusterization and time of insertion data allowed us to hypothesize that admixed elements represent sequences that harbor ancestral polymorphic sequence signatures. In conclusion, application of microevolutionary concepts in genome evolution studies is suitable as a complementary approach to phylogenetic analyses to address the evolutionary history and functional features of TEs.

Gene Duplication in the Sugarcane Genome: A Case Study of Allele Interactions and Evolutionary Patterns in Two Genic Regions

Sugarcane (Saccharum spp.) is highly polyploid and aneuploid. Modern cultivars are derived from hybridization between S. officinarum and S. spontaneum. This combination results in a genome exhibiting variable ploidy among different loci, a huge genome size (~10 Gb) and a high content of repetitive regions. An approach using genomic, transcriptomic, and genetic mapping can improve our knowledge of the behavior of genetics in sugarcane. The hypothetical HP600 and Centromere Protein C (CENP-C) genes from sugarcane were used to elucidate the allelic expression and genomic and genetic behaviors of this complex polyploid. The physically linked side-by-side genes HP600 and CENP-C were found in two different homeologous chromosome groups with ploidies of eight and ten. The first region (Region01) was a Sorghum bicolor ortholog region with all haplotypes of HP600 and CENP-C expressed, but HP600 exhibited an unbalanced haplotype expression. The second region (Region02) was a scrambled sugarcane sequence formed from different noncollinear genes containing partial duplications of HP600 and CENP-C (paralogs). This duplication resulted in a non-expressed HP600 pseudogene and a recombined fusion version of CENP-C and the orthologous gene Sobic.003G299500 with at least two chimeric gene haplotypes expressed. It was also determined that it occurred before Saccharum genus formation and after the separation of sorghum and sugarcane. A linkage map was constructed using markers from nonduplicated Region01 and for the duplication (Region01 and Region02). We compare the physical and linkage maps, demonstrating the possibility of mapping markers located in duplicated regions with markers in nonduplicated region. Our results contribute directly to the improvement of linkage mapping in complex polyploids and improve the integration of physical and genetic data for sugarcane breeding programs. Thus, we describe the complexity involved in sugarcane genetics and genomics and allelic dynamics, which can be useful for understanding complex polyploid genomes.

A miniature inverted-repeat transposable element, AddIn-MITE, located inside a WD40 gene is conserved in Andropogoneae grasses

Miniature inverted-repeat transposable elements (MITEs) have been associated with genic regions in plant genomes and may play important roles in the regulation of nearby genes via recruitment of small RNAs (sRNA) to the MITEs loci. We identified eight families of MITEs in the sugarcane genome assembly with MITE-Hunter pipeline. These sequences were found to be upstream, downstream or inserted into 67 genic regions in the genome. The position of the most abundant MITE (Stowaway-like) in genic regions, which we call AddIn-MITE, was confirmed in a WD40 gene. The analysis of four monocot species showed conservation of the AddIn-MITE sequence, with a large number of copies in their genomes. We also investigated the conservation of the AddIn-MITE’ position in the WD40 genes from sorghum, maize and, in sugarcane cultivars and wild Saccharum species. In all analyzed plants, AddIn-MITE has located in WD40 intronic region. Furthermore, the role of AddIn-MITE-related sRNA in WD40 genic region was investigated. We found sRNAs preferentially mapped to the AddIn-MITE than to other regions in the WD40 gene in sugarcane. In addition, the analysis of the small RNA distribution patterns in the WD40 gene and the structure of AddIn-MITE, suggests that the MITE region is a proto-miRNA locus in sugarcane. Together, these data provide insights into the AddIn-MITE role in Andropogoneae grasses.

The control of endopolygalacturonase expression by the sugarcane RAV transcription factor during aerenchyma formation

The development of lysigenous aerenchyma starts with cell expansion and degradation of pectin from the middle lamella, leading to cell wall modification, and culminating with cell separation. Here we report that nutritional starvation of sugarcane induced gene expression along sections of the first 5 cm of the root and between treatments. We selected two candidate genes: a RAV transcription factor, from the ethylene response factors superfamily, and an endopolygalacturonase (EPG), a glycosyl hydrolase related to homogalacturonan hydrolysis from the middle lamella. epg1 and rav1 transcriptional patterns suggest they are essential genes at the initial steps of pectin degradation during aerenchyma development in sugarcane. Due to the high complexity of the sugarcane genome, rav1 and epg1 were sequenced from 17 bacterial artificial chromosome clones containing hom(e)ologous genomic regions, and the sequences were compared with those of Sorghum bicolor. We used one hom(e)olog sequence from each gene for transactivation assays in tobacco. rav1 was shown to bind to the epg1 promoter, repressing β-glucuronidase activity. RAV repression upon epg1 transcription is the first reported link between ethylene regulation and pectin hydrolysis during aerenchyma formation. Our findings may help to elucidate cell wall degradation in sugarcane and therefore contribute to second-generation bioethanol production.

A gene-rich fraction analysis of the Passiflora edulis genome reveals highly conserved microsyntenic regions with two related Malpighiales species

Passiflora edulis is the most widely cultivated species of passionflowers, cropped mainly for industrialized juice production and fresh fruit consumption. Despite its commercial importance, little is known about the genome structure of P. edulis. To fill in this gap in our knowledge, a genomic library was built, and now completely sequenced over 100 large-inserts. Sequencing data were assembled from long sequence reads, and structural sequence annotation resulted in the prediction of about 1,900 genes, providing data for subsequent functional analysis. The richness of repetitive elements was also evaluated. Microsyntenic regions of P. edulis common to Populus trichocarpa and Manihot esculenta, two related Malpighiales species with available fully sequenced genomes were examined. Overall, gene order was well conserved, with some disruptions of collinearity identified as rearrangements, such as inversion and translocation events. The microsynteny level observed between the P. edulis sequences and the compared genomes is surprising, given the long divergence time that separates them from the common ancestor. P. edulis gene-rich segments are more compact than those of the other two species, even though its genome is much larger. This study provides a first accurate gene set for P. edulis, opening the way for new studies on the evolutionary issues in Malpighiales genomes.

A mosaic monoploid reference sequence for the highly complex genome of sugarcane

Sugarcane (Saccharum spp.) is a major crop for sugar and bioenergy production. Its highly polyploid, aneuploid, heterozygous, and interspecific genome poses major challenges for producing a reference sequence. We exploited colinearity with sorghum to produce a BAC-based monoploid genome sequence of sugarcane. A minimum tiling path of 4660 sugarcane BAC that best covers the gene-rich part of the sorghum genome was selected based on whole-genome profiling, sequenced, and assembled in a 382-Mb single tiling path of a high-quality sequence. A total of 25,316 protein-coding gene models are predicted, 17% of which display no colinearity with their sorghum orthologs. We show that the two species, S. officinarum and S. spontaneum, involved in modern cultivars differ by their transposable elements and by a few large chromosomal rearrangements, explaining their distinct genome size and distinct basic chromosome numbers while also suggesting that polyploidization arose in both lineages after their divergence.

Metabolite Profiles of Sugarcane Culm Reveal the Relationship Among Metabolism and Axillary Bud Outgrowth in Genetically Related Sugarcane Commercial Cultivars

Metabolic composition is known to exert influence on several important agronomic traits, and metabolomics, which represents the chemical composition in a cell, has long been recognized as a powerful tool for bridging phenotype-genotype interactions. In this work, sixteen truly representative sugarcane Brazilian varieties were selected to explore the metabolic networks in buds and culms, the tissues involved in the vegetative propagation of this species. Due to the fact that bud sprouting is a key trait determining crop establishment in the field, the sprouting potential among the genotypes was evaluated. The use of partial least square discriminant analysis indicated only mild differences on bud outgrowth potential under controlled environmental conditions. However, primary metabolite profiling provided information on the variability of metabolic features even under a narrow genetic background, typical for modern sugarcane cultivars. Metabolite-metabolite correlations within and between tissues revealed more complex patterns for culms in relation to buds, and enabled the recognition of key metabolites (e.g., sucrose, putrescine, glutamate, serine, and myo-inositol) affecting sprouting ability. Finally, those results were associated with the genetic background of each cultivar, showing that metabolites can be potentially used as indicators for the genetic background.

Revisiting Meiosis in Sugarcane: Chromosomal Irregularities and the Prevalence of Bivalent Configurations

Traditional sugarcane cultivars (Saccharum officinarum) proved highly susceptible to diseases, and this led breeders to progress to interspecific crosses resulting in disease resistance. A backcrossing program to S. officinarum was then required to boost sucrose content. Clonal selection across generations and incorporation of other germplasm into cultivated backgrounds established the (narrow) genetic base of modern cultivars (Saccharum spp.), which have a man-made genome. The genome complexity has inspired several molecular studies that have elucidated aspects of sugarcane genome constitution, architecture, and cytogenetics. However, there is a critical shortage of information on chromosome behavior throughout meiosis in modern cultivars. In this study, we examined the microsporogenesis of a contemporary variety, providing a detailed analysis of the meiotic process and chromosome association at diakinesis, using FISH with centromeric probes. Chromosomal abnormalities were documented by examining high quality preparations of pollen mother cells (700 in total). Approximately 70% of the cells showed abnormalities, such as metaphase chromosomes not lined up at the plate, lagging chromosomes and chromosomal bridges, and tetrad cells with micronuclei. Some dyads with asynchronous behavior were also observed. Due to the hybrid composition of the sugarcane genome, we suggest that bivalent incomplete pairing may occur in the first prophase leading to univalency. The presence of rod bivalents showing the lagging tendency is consistent with a reduction in chiasma frequency. Finally, the presence of chromatin bridges indicates the indirect occurrence of chromosomal inversions, although chromosome fragments were not clearly recognized. Possible reasons for such meiotic abnormalities and the large prevalence of bivalent formation are discussed.

TALEN-mediated targeted mutagenesis of more than 100 COMT copies/alleles in highly polyploid sugarcane improves saccharification efficiency without compromising biomass yield

Sugarcane is the world's most efficient feedstock for commercial production of bioethanol due to its superior biomass production and accumulation of sucrose in stems. Integrating first- and second-generation ethanol conversion processes will enhance the biofuel yield per unit area by utilizing both sucrose and cell wall-bound sugars for fermentation. RNAi suppression of the lignin biosynthetic gene caffeic acid O-methyltransferase (COMT) has been demonstrated to improve bioethanol production from lignocellulosic biomass. Genome editing has been used in a number of crops for creation of loss of function phenotypes but is very challenging in sugarcane due to its highly polyploid genome. In this study, a conserved region of COMT was targeted with a single-transcription activator-like effector nuclease (TALEN) pair for multi-allelic mutagenesis to modify lignin biosynthesis in sugarcane. Field-grown TALEN-mediated COMT mutants showed up to 19.7% lignin reduction and significantly decreased syringyl to guaiacyl (S/G) ratio resulting in an up to 43.8% improved saccharification efficiency. Biomass production of COMT mutant lines with superior saccharification efficiency did not differ significantly from the original cultivar under replicated field conditions. Sanger sequencing of cloned COMT amplicons (1351-1657 bp) revealed co-editing of 107 of the 109 unique COMT copies/alleles in vegetative progeny of line CB6 using a single TALEN pair. Line CB6 combined altered cell wall composition and drastically improved saccharification efficiency with good agronomic performance. These findings confirm the feasibility of co-mutagenesis of a very large number of target alleles/copies for improvement in crops with complex genomes.

"Targeted Sequencing by Gene Synteny," a New Strategy for Polyploid Species: Sequencing and Physical Structure of a Complex Sugarcane Region

Sugarcane exhibits a complex genome mainly due to its aneuploid nature and high ploidy level, and sequencing of its genome poses a great challenge. Closely related species with well-assembled and annotated genomes can be used to help assemble complex genomes. Here, a stable quantitative trait locus (QTL) related to sugar accumulation in sorghum was successfully transferred to the sugarcane genome. Gene sequences related to this QTL were identified in silico from sugarcane transcriptome data, and molecular markers based on these sequences were developed to select bacterial artificial chromosome (BAC) clones from the sugarcane variety SP80-3280. Sixty-eight BAC clones containing at least two gene sequences associated with the sorghum QTL were sequenced using Pacific Biosciences (PacBio) technology. Twenty BAC sequences were found to be related to the syntenic region, of which nine were sufficient to represent this region. The strategy we propose is called "targeted sequencing by gene synteny," which is a simpler approach to understanding the genome structure of complex genomic regions associated with traits of interest.

The Challenge of Analyzing the Sugarcane Genome

Reference genome sequences have become key platforms for genetics and breeding of the major crop species. Sugarcane is probably the largest crop produced in the world (in weight of crop harvested) but lacks a reference genome sequence. Sugarcane has one of the most complex genomes in crop plants due to the extreme level of polyploidy. The genome of modern sugarcane hybrids includes sub-genomes from two progenitors Saccharum officinarum and S. spontaneum with some chromosomes resulting from recombination between these sub-genomes. Advancing DNA sequencing technologies and strategies for genome assembly are making the sugarcane genome more tractable. Advances in long read sequencing have allowed the generation of a more complete set of sugarcane gene transcripts. This is supporting transcript profiling in genetic research. The progenitor genomes are being sequenced. A monoploid coverage of the hybrid genome has been obtained by sequencing BAC clones that cover the gene space of the closely related sorghum genome. The complete polyploid genome is now being sequenced and assembled. The emerging genome will allow comparison of related genomes and increase understanding of the functioning of this polyploidy system. Sugarcane breeding for traditional sugar and new energy and biomaterial uses will be enhanced by the availability of these genomic resources.

Development and Applications of Chromosome-Specific Cytogenetic BAC-FISH Probes in S. spontaneum

Saccharum spontaneum is a major Saccharum species that contributed to the origin of modern sugarcane cultivars, and due to a high degree of polyploidy is considered to be a plant species with one of the most complex genetics. Fluorescence in situ hybridization (FISH) is a powerful and widely used tool in genome studies. Here, we demonstrated that FISH based on bacterial artificial chromosome (BAC) clones can be used as a specific cytological marker to identify S. spontaneum individual chromosomes and study the relationship between S. spontaneum and other related species. We screened low-copy BACs as probes from the sequences of a high coverage of S. spontaneum BAC library based on BLAST search of the sorghum genome. In total, we isolated 49 positive BAC clones, and identified 27 BAC clones that can give specific signals on the S. spontaneum chromosomes. Of the 27 BAC probes, 18 were confirmed to be able to discriminate the eight basic chromosomes of S. spontaneum. Moreover, BAC-24, BAC-66, BAC-78, BAC-69, BAC-71, BAC-73, and BAC-77 probes were used to construct physical maps of chromosome 1 and chromosome 2 of S. spontaneum, which indicated synteny in Sb01 between S. spontaneum and sorghum. Furthermore, we found that BAC-14 and BAC-19 probes, corresponding to the sorghum chromosomes 2 and 8, respectively, localized to different arms of the same S. spontaneum chromosome, suggesting that there was an inter-chromosomal rearrangement event between S. spontaneum and sorghum. Our study provides the first set of chromosome-specific cytogenetic markers in Saccharum and is critical for future advances in cytogenetics and genome sequencing studies in Saccharum.

Analysis of Three Sugarcane Homo/Homeologous Regions Suggests Independent Polyploidization Events of Saccharum officinarum and Saccharum spontaneum

Whole genome duplication has played an important role in plant evolution and diversification. Sugarcane is an important crop with a complex hybrid polyploid genome, for which the process of adaptation to polyploidy is still poorly understood. In order to improve our knowledge about sugarcane genome evolution and the homo/homeologous gene expression balance, we sequenced and analyzed 27 BACs (Bacterial Artificial Chromosome) of sugarcane R570 cultivar, containing the putative single-copy genes LFY (seven haplotypes), PHYC (four haplotypes), and TOR (seven haplotypes). Comparative genomic approaches showed that these sugarcane loci presented a high degree of conservation of gene content and collinearity (synteny) with sorghum and rice orthologous regions, but were invaded by transposable elements (TE). All the homo/homeologous haplotypes of LFY, PHYC, and TOR are likely to be functional, because they are all under purifying selection (dN/dS ≪ 1). However, they were found to participate in a nonequivalently manner to the overall expression of the corresponding gene. SNPs, indels, and amino acid substitutions allowed inferring the S. officinarum or S. spontaneum origin of the TOR haplotypes, which further led to the estimation that these two sugarcane ancestral species diverged between 2.5 and 3.5 Ma. In addition, analysis of shared TE insertions in TOR haplotypes suggested that two autopolyploidization may have occurred in the lineage that gave rise to S. officinarum, after its divergence from S. spontaneum.