Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L

Development of SLAF-Sequence and Multiplex SNaPshot Panels for Population Genetic Diversity Analysis and Construction of DNA Fingerprints for Sugarcane

A genetic diversity analysis and identification of plant germplasms and varieties are important and necessary for plant breeding. Deoxyribonucleotide (DNA) fingerprints based on genomic molecular markers play an important role in accurate germplasm identification. In this study, Specific-Locus Amplified Fragment Sequencing (SLAF-seq) was conducted for a sugarcane population with 103 cultivated and wild accessions. In total, 105,325 genomic single nucleotide polymorphisms (SNPs) were called successfully to analyze population components and genetic diversity. The genetic diversity of the population was complex and clustered into two major subpopulations. A principal component analysis (PCA) showed that these accessions could not be completely classified based on geographical origin. After filtration, screening, and comparison, 192 uniformly-distributed SNP loci were selected for the 32 chromosomes of sugarcane. An SNP complex genotyping detection system was established using the SNaPshot typing method and used for the precise genotyping and identification of 180 sugarcane germplasm samples. According to the stability and polymorphism of the SNPs, 32 high-quality SNP markers were obtained and successfully used to construct the first SNP fingerprinting and quick response codes (QR codes) for sugarcane. The results provide new insights for genotyping, classifying, and identifying germplasm and resources for sugarcane breeding

A High-Quality Genome Assembly of Sorghum dochna

Sweet sorghum (Sorghum dochna) is a high-quality bio-energy crop that also serves as food for humans and animals. However, there is little information on the genomic characteristics of S. dochna. In this study, we presented a high-quality assembly of S. dochna with PacBio long reads, Illumina short reads, high-throughput chromosome capture technology (Hi-C) sequencing data, gene annotation, and a comparative genome analysis. The results showed that the genome of S. dochna was assembled to 777 Mb with a contig N50 of 553.47 kb and a scaffold N50 of 727.11 kb. In addition, the gene annotation predicted 37,971 genes and 39,937 transcripts in the genome of S. dochna. A Venn analysis revealed a set of 7,988 common gene annotations by integrating five databases. A Cafe software analysis showed that 191 gene families were significantly expanded, while 3,794 were significantly contracted in S. dochna. A GO enrichment analysis showed that the expanded gene families were primarily clustered in the metabolic process, DNA reconstruction, and DNA binding among others. The high-quality genome map constructed in this study provides a biological basis for the future analysis of the biological characteristics of S. dochna, which is crucial for its breeding.

Expression Profiling and MicroRNA Regulatory Networks of Homeobox Family Genes in Sugarcane Saccharum spontaneum L

Homeobox (HB) genes play important roles in plant growth and development processes, particularly in the formation of lateral organs. Thus, they could influence leaf morphogenesis and biomass formation in plants. However, little is known about HBs in sugarcane, a crucial sugar crop, due to its complex genetic background. Here, 302 allelic sequences for 104 HBs were identified and divided into 13 subfamilies in sugarcane Saccharum spontaneum. Comparative genomics revealed that whole-genome duplication (WGD)/segmental duplication significantly promoted the expansion of the HB family in S. spontaneum, with SsHB26, SsHB63, SsHB64, SsHB65, SsHB67, SsHB95, and SsHB96 being retained from the evolutionary event before the divergence of dicots and monocots. Based on the analysis of transcriptome and degradome data, we speculated that SsHB15 and SsHB97 might play important roles in regulating sugarcane leaf morphogenesis, with miR166 and SsAGO10 being involved in the regulation of SsHB15 expression. Moreover, subcellular localization and transcriptional activity detection assays demonstrated that these two genes, SsHB15 and SsHB97, were functional transcription factors. This study demonstrated the evolutionary relationship and potential functions of SsHB genes and will enable the further investigation of the functional characterization and the regulatory mechanisms of SsHBs.

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.

Nogueira E, Baldani JI, Hemerly AS. Sugarcane Genotypes with Contrasting Biological Nitrogen Fixation Efficiencies Differentially Modulate Nitrogen Metabolism, Auxin Signaling, and Microorganism Perception Pathways

Sugarcane is an economically important crop that is used for the production of fuel ethanol. Diazotrophic bacteria have been isolated from sugarcane tissues, without causing visible plant anatomical changes or disease symptoms. These bacteria can be beneficial to the plant by promoting root growth and an increase in plant yield. Different rates of Biological Nitrogen Fixation (BNF) were observed in different genotypes. The aim of this work was to conduct a comprehensive molecular and physiological analysis of two model genotypes for contrasting BNF efficiency in order to unravel plant genes that are differentially regulated during a natural association with diazotrophic bacteria. A next-generation sequencing of RNA samples from the genotypes SP70-1143 (high-BNF) and Chunee (low-BNF) was performed. A differential transcriptome analysis showed that several pathways were differentially regulated among the two BNF-contrasting genotypes, including nitrogen metabolism, hormone regulation and bacteria recognition. Physiological analyses, such as nitrogenase and GS activity quantification, bacterial colonization, auxin response and root architecture evaluation, supported the transcriptome expression analyses. The differences observed between the genotypes may explain, at least in part, the differences in BNF contributions. Some of the identified genes might be involved in key regulatory processes for a beneficial association and could be further used as tools for obtaining more efficient BNF genotypes.

Root nitrate uptake in sugarcane (Saccharum spp.) is modulated by transcriptional and presumably posttranscriptional regulation of the NRT2.1/NRT3.1 transport system

KEY MESSAGE

Nitrate uptake in sugarcane roots is regulated at the transcriptional and posttranscriptional levels based on the physiological status of the plant and is likely a determinant mechanism for discrimination against nitrate. Sugarcane (Saccharum spp.) is one of the most suitable energy crops for biofuel feedstock, but the reduced recovery of nitrogen (N) fertilizer by sugarcane roots increases the crop carbon footprint. The low nitrogen use efficiency (NUE) of sugarcane has been associated with the significantly low nitrate uptake, which limits the utilization of the large amount of nitrate available in agricultural soils. To understand the regulation of nitrate uptake in sugarcane roots, we identified the major canonical nitrate transporter genes (NRTs-NITRATE TRANSPORTERS) and then determined their expression profiles in roots under contrasting N conditions. Correlation of gene expression with ¹⁵N-nitrate uptake revealed that under N deprivation or inorganic N (ammonium or nitrate) supply in N-sufficient roots, the regulation of ScNRT2.1 and ScNRT3.1 expression is the predominant mechanism for the modulation of the activity of the nitrate high-affinity transport system. Conversely, in N-deficient roots, the induction of ScNRT2.1 and ScNRT3.1 transcription is not correlated with the marked repression of nitrate uptake in response to nitrate resupply or high N provision, which suggested the existence of a posttranscriptional regulatory mechanism. Our findings suggested that high-affinity nitrate uptake is regulated at the transcriptional and presumably at the posttranscriptional levels based on the physiological N status and that the regulation of NRT2.1 and NRT3.1 activity is likely a determinant mechanism for the discrimination against nitrate uptake observed in sugarcane roots, which contributes to the low NUE in this crop species.

Characterization and Functional Implications of the Nonexpressor of Pathogenesis-Related Genes 1 (NPR1) in Saccharum

Sugarcane (Saccharum spp.) is an important sugar and energy crop worldwide. As a core regulator of the salicylic acid (SA) signaling pathway, nonexpressor of pathogenesis-related genes 1 (NPR1) plays a significant role in the response of the plant to biotic and abiotic stresses. However, there is currently no report on the NPR1-like gene family in sugarcane. In this study, a total of 18 NPR1-like genes were identified in Saccharum spontaneum and classified into three clades (clade I, II, and III). The cis-elements predicted in the promotors revealed that the sugarcane NPR1-like genes may be involved in various phytohormones and stress responses. RNA sequencing and quantitative real-time PCR analysis demonstrated that NPR1-like genes were differentially expressed in sugarcane tissues and under Sporisorium scitamineum stress. In addition, a novel ShNPR1 gene from Saccharum spp. hybrid ROC22 was isolated by homologous cloning and validated to be a nuclear-localized clade II member. The ShNPR1 gene was constitutively expressed in all the sugarcane tissues, with the highest expression level in the leaf and the lowest in the bud. The expression level of ShNPR1 was decreased by the plant hormones salicylic acid (SA) and abscisic acid (ABA). Additionally, the transient expression showed that the ShNPR1 gene plays a positive role in Nicotiana benthamiana plants’ defense response to Ralstonia solanacearum and Fusarium solani var. coeruleum. This study provided comprehensive information for the NPR1-like family in sugarcane, which should be helpful for functional characterization of sugarcane NPR1-like genes in the future.

Equilibrated evolution of the mixed auto-/allopolyploid haplotype-resolved genome of the invasive hexaploid Prussian carp

Understanding genome evolution of polyploids requires dissection of their often highly similar subgenomes and haplotypes. Polyploid animal genome assemblies so far restricted homologous chromosomes to a 'collapsed' representation. Here, we sequenced the genome of the asexual Prussian carp, which is a close relative of the goldfish, and present a haplotype-resolved chromosome-scale assembly of a hexaploid animal. Genome-wide comparisons of the 150 chromosomes with those of two ancestral diploid cyprinids and the allotetraploid goldfish and common carp revealed the genomic structure, phylogeny and genome duplication history of its genome. It consists of 25 syntenic, homeologous chromosome groups and evolved by a recent autoploid addition to an allotetraploid ancestor. We show that de-polyploidization of the alloploid subgenomes on the individual gene level occurred in an equilibrated fashion. Analysis of the highly conserved actinopterygian gene set uncovered a subgenome dominance in duplicate gene loss of one ancestral chromosome set.

Impact of Agroclimatic Variables on Proteogenomics in Sugar Cane (Saccharum spp.) Plant Productivity

Sugar cane (Saccharum spp. hybrids) is a major crop for sugar and renewable bioenergy worldwide, grown in arid and semiarid regions. China, the world's fourth-largest sugar producer after Brazil, India, and the European Union, all share ∼80% of the global production, and the remaining ∼20% of sugar comes from sugar beets, mostly grown in the temperate regions of the Northern Hemisphere, also used as a raw material in production of bioethanol for renewable energy. In view of carboxylation strategies, sugar cane qualifies as one of the best C4 crop. It has dual CO2 concentrating mechanisms located in its unique Krantz anatomy, having dimorphic chloroplasts located in mesophylls and bundle sheath cells for integrated operation of C4 and C3 carbon fixation cycles, regulated by enzymes to upgrade/sustain an ability for improved carbon assimilation to acquire an optimum carbon economy by producing enhanced plant biomass along with sugar yield under elevated temperature and strong irradiance with improved water-use efficiency. These superior intrinsic physiological carbon metabolisms encouraged us to reveal and recollect the facts for moving ahead with the molecular approaches to reveal the expression of proteogenomics linked with plant productivity under abiotic stress during its cultivation in specific agrizones globally.

Genome architecture and tetrasomic inheritance of autotetraploid potato

Potato (Solanum tuberosum) is the most consumed non-cereal food crop. Most commercial potato cultivars are autotetraploids with highly heterozygous genomes, severely hampering genetic analyses and improvement. By leveraging the state-of-the-art sequencing technologies and polyploid graph binning, we achieved a chromosome-scale, haplotype-resolved genome assembly of a cultivated potato, Cooperation-88 (C88). Intra-haplotype comparative analyses revealed extensive sequence and expression differences in this tetraploid genome. We identified haplotype-specific pericentromeres on chromosomes, suggesting a distinct evolutionary trajectory of potato homologous centromeres. Furthermore, we detected double reduction events that are unevenly distributed on haplotypes in 1021 of 1034 selfing progeny, a feature of autopolyploid inheritance. By distinguishing maternal and paternal haplotype sets in C88, we simulated the origin of heterosis in cultivated tetraploid with a survey of 3110 tetra-allelic loci with deleterious mutations, which were masked in the heterozygous condition by two parents. This study provides insights into the genomic architecture of autopolyploids and will guide their breeding.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

Taxonomically Restricted Genes Are Associated With Responses to Biotic and Abiotic Stresses in Sugarcane (Saccharum spp.)

Orphan genes (OGs) are protein-coding genes that are restricted to particular clades or species and lack homology with genes from other organisms, making their biological functions difficult to predict. OGs can rapidly originate and become functional; consequently, they may support rapid adaptation to environmental changes. Extensive spread of mobile elements and whole-genome duplication occurred in the Saccharum group, which may have contributed to the origin and diversification of OGs in the sugarcane genome. Here, we identified and characterized OGs in sugarcane, examined their expression profiles across tissues and genotypes, and investigated their regulation under varying conditions. We identified 319 OGs in the Saccharum spontaneum genome without detected homology to protein-coding genes in green plants, except those belonging to Saccharinae. Transcriptomic analysis revealed 288 sugarcane OGs with detectable expression levels in at least one tissue or genotype. We observed similar expression patterns of OGs in sugarcane genotypes originating from the closest geographical locations. We also observed tissue-specific expression of some OGs, possibly indicating a complex regulatory process for maintaining diverse functional activity of these genes across sugarcane tissues and genotypes. Sixty-six OGs were differentially expressed under stress conditions, especially cold and osmotic stresses. Gene co-expression network and functional enrichment analyses suggested that sugarcane OGs are involved in several biological mechanisms, including stimulus response and defence mechanisms. These findings provide a valuable genomic resource for sugarcane researchers, especially those interested in selecting stress-responsive genes.

Plant pan-genomics: recent advances, new challenges, and roads ahead

Pan-genomics can encompass most of the genetic diversity of a species or population and has proved to be a powerful tool for studying genomic evolution and the origin and domestication of species, and for providing information for plant improvement. Plant genomics has greatly progressed because of improvements in sequencing technologies and the rapid reduction of sequencing costs. Nevertheless, pan-genomics still presents many challenges, including computationally intensive assembly methods, high costs with large numbers of samples, ineffective integration of big data, and difficulty in applying it to downstream multi-omics analysis and breeding research. In this review, we summarize the definition and recent achievements of plant pan-genomics, computational technologies used for pan-genome construction, and the applications of pan-genomes in plant genomics and molecular breeding. We also discuss challenges and perspectives for future pan-genomics studies and provide a detailed pipeline for sample selection, genome assembly and annotation, structural variation identification, and construction and application of graph-based pan-genomes. The aim is to provide important guidance for plant pan-genome research and a better understanding of the genetic basis of genome evolution, crop domestication, and phenotypic diversity for future studies.

Chromosome-level and haplotype-resolved genome provides insight into the tetraploid hybrid origin of patchouli

Patchouli (Pogostemon cablin (Blanco) Benth.), a member of the Lamiaceae family, is an important aromatic plant that has been widely used in medicine and perfumery. Here, we report a 1.94 Gb chromosome-scale assembly of the patchouli genome (contig N50 = 7.97 Mb). The gene annotation reveals that tandem duplication of sesquiterpene biosynthetic genes may be a major contributor to the biosynthesis of patchouli bioactivity components. We further phase the genome into two distinct subgenomes (A and B), and identify a chromosome substitution event that have occurred between them. Further investigations show that a burst of universal LTR-RTs in the A subgenome lead to the divergence between two subgenomes. However, no significant subgenome dominance is detected. Finally, we track the evolutionary scenario of patchouli including whole genome tetraploidization, subgenome divergency, hybridization, and chromosome substitution, which are the key forces to determine the complexity of patchouli genome. Our work sheds light on the evolutionary history of patchouli and offers unprecedented genomic resources for fundamental patchouli research and elite germplasm development.

The ploidy level of patchouli, an aromatic plant in the Lamiaceae family, remain unclear. Here, the authors assemble a chromosome-level and haplotype-resolved genome for patchouli and reveal that it is tetraploid hybrid as well as compensated aneuploidy.

Repetitive Sequence Barcode Probe for Karyotype Analysis in Tripidium arundinaceum

The barcode probe is a convenient and efficient tool for molecular cytogenetics. Tripidium arundinaceum, as a polyploid wild allied genus of Saccharum, is a useful genetic resource that confers biotic and abiotic stress resistance for sugarcane breeding. Unfortunately, the basic cytogenetic information is still unclear due to the complex genome. We constructed the Cot-20 library for screening moderately and highly repetitive sequences from T. arundinaceum, and the chromosomal distribution of these repetitive sequences was explored. We used the barcode of repetitive sequence probes to distinguish the ten chromosome types of T. arundinaceum by fluorescence in situ hybridization (FISH) with Ea-0907, Ea-0098, and 45S rDNA. Furthermore, the distinction among homology chromosomes based on repetitive sequences was constructed in T. arundinaceum by the repeated FISH using the barcode probes including Ea-0663, Ea-0267, EaCent, 5S rDNA, Ea-0265, Ea-0070, and 45S rDNA. We combined these probes to distinguish 37 different chromosome types, suggesting that the repetitive sequences may have different distributions on homologous chromosomes of T. arundinaceum. In summary, this method provide a basis for the development of similar applications for cytogenetic analysis in other species.

A Comprehensive Identification and Expression Analysis of VQ Motif-Containing Proteins in Sugarcane (Saccharum spontaneum L.) under Phytohormone Treatment and Cold Stress

The VQ motif-containing proteins play a vital role in various processes such as growth, resistance to biotic and abiotic stresses and development. However, there is currently no report on the VQ genes in sugarcane (Saccharum spp.). Herein, 78 VQ genes in Saccharum spontaneum were identified and classified into nine subgroups (I-IX) by comparative genomic analyses. Each subgroup had a similar structural and conservative motif. These VQ genes expanded mainly through whole-genome segmental duplication. The cis-regulatory elements (CREs) of the VQ genes were widely involved in stress responses, phytohormone responses and physiological regulation. The RNA-seq data showed that SsVQ gene expression patterns in 10 different samples, including different developmental stages, revealed distinct temporal and spatial patterns. A total of 23 SsVQ genes were expressed in all tissues, whereas 13 SsVQ genes were not expressed in any tissues. Sequence Read Archive (SRA) data showed that the majority of SsVQs responded to cold and drought stress. In addition, quantitative real-time PCR analysis showed that the SsVQs were variously expressed under salicylic acid (SA), jasmonic acid (JA), abscisic acid (ABA) and cold treatment. This study conducted a full-scale analysis of the VQ gene family in sugarcane, which could be beneficial for the functional characterization of sugarcane VQ genes and provide candidate genes for molecular resistance breeding in cultivated sugarcane in the future.

ScGAIL, a sugarcane N-terminal truncated DELLA-like protein, participates in gibberellin signaling in Arabidopsis

The hormone gibberellin (GA) is crucial for internode elongation in sugarcane. DELLA proteins are critical negative regulators of the GA signaling pathway. ScGAI encodes a DELLA protein that was previously implicated in the regulation of sugarcane culm development. Here, we characterized ScGAI-like (ScGAIL) in sugarcane, which lacked the N-terminal region but was otherwise homologous to ScGAI. ScGAIL differed from ScGAI in its chromosomal location, expression patterns, and cellular localization. Although transgenic Arabidopsis overexpressing ScGAIL were insensitive to GAs, GA synthesis was affected in these plants, suggesting that ScGAIL disrupted the GA signaling pathway. After GA treatment, the expression patterns of GA-associated genes differed between ScGAIL-overexpressing and wild-type Arabidopsis, and the degradation of AtDELLA proteins in transgenic lines was significantly inhibited compared with wild-type lines. A sugarcane GID1 gene (ScGID1) encoding a putative GA receptor was isolated and interacted with ScGAIL in a GA-independent manner. Five ScGAIL-interacting proteins were verified by yeast two-hybrid assays, and only one interacted with ScGAI. Therefore, ScGAIL may inhibit plant growth by modulating the GA signaling pathway.

Genomic insights into the recent chromosome reduction of autopolyploid sugarcane Saccharum spontaneum

Saccharum spontaneum is a founding Saccharum species and exhibits wide variation in ploidy levels. We have assembled a high-quality autopolyploid genome of S. spontaneum Np-X (2n = 4x = 40) into 40 pseudochromosomes across 10 homologous groups, that better elucidates recent chromosome reduction and polyploidization that occurred circa 1.5 million years ago (Mya). One paleo-duplicated chromosomal pair in Saccharum, NpChr5 and NpChr8, underwent fission followed by fusion accompanied by centromeric split around 0.80 Mya. We inferred that Np-X, with x = 10, most likely represents the ancestral karyotype, from which x = 9 and x = 8 evolved. Resequencing of 102 S. spontaneum accessions revealed that S. spontaneum originated in northern India from an x = 10 ancestor, which then radiated into four major groups across the Indian subcontinent, China, and Southeast Asia. Our study suggests new directions for accelerating sugarcane improvement and expands our knowledge of the evolution of autopolyploids.

A chromosome-level genome assembly of Paracymoriza distinctalis (Lepidoptera: Crambidae: Acentropinae)

Paracymoriza distinctalis is a semiaquatic lepidopteran insect, which is of great value for studying the differentiation of the Pyraloidea superfamily. However, the understanding of heredity, evolution, and functional genomics of P. distinctalis are limited by few genome-wide resources. Here, we applied PacBio sequencing and the chromosome capture technique to assemble the first P. distinctalis genome from a single female individual. The genome size is 1.2 Gb with 32 chromosomes and the N50 is 38.91 Mb. Approximately 576.37 Mb, accounting for 48.93% of the genome, was identified as repeats. The genome comprises 39,003 protein-coding genes, 66.56% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of P. distinctalis and Chilo suppressalis lived ~83.5 million years ago. This chromosome-level genome assembly work is not only conducive to the understanding of P. distinctalis, but also may promote the study of the genomes of other lepidopteran species.

De novo phasing resolves haplotype sequences in complex plant genomes

Genome phasing is a recently developed assembly method that separates heterozygous eukaryotic genomic regions and builds haplotype-resolved assemblies. Because differences between haplotypes are ignored in most published de novo genomes, assemblies are available as consensus genomes consisting of haplotype mixtures, thus increasing the need for genome phasing. Here, we review the operating principles and characteristics of several freely available and widely used phasing tools (TrioCanu, FALCON-Phase, and ALLHiC). An examination of downstream analyses using haplotype-resolved genome assemblies in plants indicated significant differences among haplotypes regarding chromosomal rearrangements, sequence insertions, and expression of specific alleles that contribute to the acquisition of the biological characteristics of plant species. Finally, we suggest directions to solve addressing limitations of current genome-phasing methods. This review provides insights into the current progress, limitations, and future directions of de novo genome phasing, which will enable researchers to easily access and utilize genome-phasing in studies involving highly heterozygous complex plant genomes.

Identification of Genetic Loci for Sugarcane Leaf Angle at Different Developmental Stages by Genome-Wide Association Study

Sugarcane (Saccharum spp.) is an efficient crop mainly used for sugar and bioethanol production. High yield and high sucrose of sugarcane are always the fundamental demands in sugarcane growth worldwide. Leaf angle and size of sugarcane can be attributed to planting density, which was associated with yield. In this study, we performed genome-wide association studies (GWAS) with a panel of 216 sugarcane core parents and their derived lines (natural population) to determine the genetic basis of leaf angle and key candidate genes with +2, +3, and +4 leaf at the seedling, elongation, and mature stages. A total of 288 significantly associated loci of sugarcane leaf angle at different developmental stages (eight phenotypes) were identified by GWAS with 4,027,298 high-quality SNP markers. Among them, one key locus and 11 loci were identified in all three stages and two stages, respectively. An InDel marker (SNP Ss6A_102766953) linked to narrow leaf angle was obtained. Overall, 4,089 genes were located in the confidence interval of significant loci, among which 3,892 genes were functionally annotated. Finally, 13 core parents and their derivatives tagged with SNPs were selected for marker-assisted selection (MAS). These candidate genes are mainly related to MYB transcription factors, auxin response factors, serine/threonine protein kinases, etc. They are directly or indirectly associated with leaf angle in sugarcane. This research provided a large number of novel genetic resources for the improvement of leaf angles and simultaneously to high yield and high bioethanol production.

The integrated genomics of crop domestication and breeding

As a major event in human civilization, wild plants were successfully domesticated to be crops, largely owing to continuing artificial selection. Here, we summarize new discoveries made during the past decade in crop domestication and breeding. The construction of crop genome maps and the functional characterization of numerous trait genes provide foundational information. Approaches to read, interpret, and write complex genetic information are being leveraged in many plants for highly efficient de novo or re-domestication. Understanding the underlying mechanisms of crop microevolution and applying the knowledge to agricultural productions will give possible solutions for future challenges in food security.

Allele expression biases in mixed-ploid sugarcane accessions

Allele-specific expression (ASE) represents differences in the magnitude of expression between alleles of the same gene. This is not straightforward for polyploids, especially autopolyploids, as knowledge about the dose of each allele is required for accurate estimation of ASE. This is the case for the genomically complex Saccharum species, characterized by high levels of ploidy and aneuploidy. We used a Beta-Binomial model to test for allelic imbalance in Saccharum, with adaptations for mixed-ploid organisms. The hierarchical Beta-Binomial model was used to test if allele expression followed the expectation based on genomic allele dosage. The highest frequencies of ASE occurred in sugarcane hybrids, suggesting a possible influence of interspecific hybridization in these genotypes. For all accessions, genes showing ASE (ASEGs) were less frequent than those with balanced allelic expression. These genes were related to a broad range of processes, mostly associated with general metabolism, organelles, responses to stress and responses to stimuli. In addition, the frequency of ASEGs in high-level functional terms was similar among the genotypes, with a few genes associated with more specific biological processes. We hypothesize that ASE in Saccharum is largely a genotype-specific phenomenon, as a large number of ASEGs were exclusive to individual accessions.

Identification of sex determination locus in sea cucumber Apostichopus japonicus using genome-wide association study

Background

Sex determination mechanisms are complicated and diverse across taxonomic categories. Sea cucumber Apostichopus japonicus is a benthic echinoderm, which is the closest group of invertebrates to chordate, and important economic and ecologically aquaculture species in China. A. japonicus is dioecious, and no phenotypic differences between males and females can be detected before sexual maturation. Identification of sex determination locus will broaden knowledge about sex-determination mechanism in echinoderms, which allows for the identification of sex-linked markers and increases the efficiency of sea cucumber breeding industry.

Results

Here, we integrated assembly of a novel chromosome-level genome and resequencing of female and male populations to investigate the sex determination mechanisms of A. japonicus. We built a chromosome-level genome assembly AJH1.0 using Hi-C technology. The assembly AJH1.0 consists of 23 chromosomes ranging from 22.4 to 60.4 Mb. To identify the sex-determination locus of A. japonicus, we conducted genome-wide association study (GWAS) and analyses of distribution characteristics of sex-specific SNPs and fixation index F_ST. The GWAS analysis showed that multiple sex-associated loci were located on several chromosomes, including chromosome 4 (24.8%), followed by chromosome 9 (10.7%), chromosome 17 (10.4%), and chromosome 18 (14.1%). Furthermore, analyzing the homozygous and heterozygous genotypes of plenty of sex-specific SNPs in females and males confirmed that A. japonicus might have a XX/XY sex determination system. As a physical region of 10 Mb on chromosome 4 included the highest number of sex-specific SNPs and higher F_ST values, this region was considered as the candidate sex determination region (SDR) in A. japonicus.

Conclusions

In the present study, we integrated genome-wide association study and analyses of sex-specific variations to investigate sex determination mechanisms. This will bring novel insights into gene regulation during primitive gonadogenesis and differentiation and identification of master sex determination gene in sea cucumber. In the sea cucumber industry, investigation of molecular mechanisms of sex determination will be helpful for artificial fertilization and precise breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08632-3.

Characterization of silicon transporter gene family in Saccharum and functional analysis of the ShLsi6 gene in biotic stress

Silicon, one of the most prevalent elements in the soil, is beneficial for plant growth and defense against different stresses. The silicon transporter gene (Lsi) plays an important role in the uptake and transport of silicon in higher plants. In this study, a total of 32 Lsi genes, including 20 SsLsi in sugarcane wild species Saccharum spontaneum, 5 ShLsi in Saccharum hybrid cultivar R570 and 7 SbLsi in sugarcane related species Sorghum bicolor, were identified and classified into three groups. Bioinformatics analysis showed that instability, hydrophobicity, localization of cell membranes and vacuoles were the main features of the Lsi proteins. Whole genome and segmental duplication contributed to the main expansion of Lsi gene family. Collinearity analysis of the Lsi genes showed that S. spontanum and R570 had a collinear relationship with monocotyledonous plants S. bicolor and Oryza sativa, but not with dicotyledonous plants Arabidopsis thaliana and Vitis vinifera. The replicated Lsi genes were mainly subjected to strong selection pressure for purification. The diverse cis-regulatory elements in the promoter of SsLsi, ShLsi and SbLsi genes suggested that they were widely involved in the response of plants to various stresses and the regulation of the growth and development. Transcriptome data and real time quantitative PCR analysis showed that the Lsi genes exhibited different expression profiles in sugarcane tissues and under Sporisorium scitamineum, drought and cold stresses. In addition, the cDNA and genomic DNA sequences of ShLsi6 that was homologous to SsLsi1b gene was cloned from Saccharum hybrid cultivar ROC22. Transient expression analysis showed that, compared with the control, Nicotiana benthamiana leaves which overexpressed the ShLsi6 gene showed a high sensitivity after inoculation with tobacco pathogens Ralstonia solanacearum and Fusarium solani var. coeruleum. This study provides important information for further functional analysis of Lsi genes and resistant breeding in sugarcane.

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.

The flying spider-monkey tree fern genome provides insights into fern evolution and arborescence

To date, little is known about the evolution of fern genomes, with only two small genomes published from the heterosporous Salviniales. Here we assembled the genome of Alsophila spinulosa, known as the flying spider-monkey tree fern, onto 69 pseudochromosomes. The remarkable preservation of synteny, despite resulting from an ancient whole-genome duplication over 100 million years ago, is unprecedented in plants and probably speaks to the uniqueness of tree ferns. Our detailed investigations into stem anatomy and lignin biosynthesis shed new light on the evolution of stem formation in tree ferns. We identified a phenolic compound, alsophilin, that is abundant in xylem, and we provided the molecular basis for its biosynthesis. Finally, analysis of demographic history revealed two genetic bottlenecks, resulting in rapid demographic declines of A. spinulosa. The A. spinulosa genome fills a crucial gap in the plant genomic landscape and helps elucidate many unique aspects of tree fern biology.

Characterization and Phylogenetic Analyses of the Complete Mitochondrial Genome of Sugarcane (Saccharum spp. Hybrids) Line A1

Modern sugarcane cultivars are highly polyploid with complex nuclear genomic genetic background, while their mitochondrion (mt) genomes are much simpler, smaller and more manageable and could provide useful phylogenetic information. In this study, the mt genome of a modern commercial cultivar A1 was sequenced via Illumina Hiseq XTen and PacBio Sequel platform. The assembled and annotated mitochondrial genomes of A1 were composed of two circular DNA molecules, one large and one small, which were named Chromosome 1 and Chromosome 2. The two distinct circular chromosomes of mitogenome construct is consisted with other sugarcane cultivars i.e., Saccharum officinarum Khon Kaen 3 and Saccharum spp. hybrids ROC22 and FN15. The Chromosome 1 of A1 mitogenome is 300,822 bp in length with the GC content of 43.94%, and 7.14% of Chromosome 1 sequences (21,468 nucleotides) are protein coding genes (PCGs) while 92.86% (279,354 nucleotides) are intergenic region. The length of Chromosome 2 is 144,744 bp with the GC content of 43.57%, and 8.20% of Chromosome 2 sequences (11,865 nucleotides) are PCGs while 91.80% (132,879 nucleotides) are intergenic region. A total of 43 genes are located on Chromosome 1, which contains 22 PCGs (six nad genes, four rps genes, four atp genes, three ccm genes, three cox genes, one mat gene and one mtt gene) and 21 non-coding genes including 15 tRNAs and 6 rRNAs. Chromosome 2 includes 18 genes in total, which contains 13 PCGs (four nad genes, three rps genes, two atp genes, one ccm gene, one cob gene, one cox gene and one rpl gene) and five non-coding genes (tRNA genes). Analysis of codon usage of 35 PCGs showed that codon ending in A/U was preferred. Investigation of gene composition indicated that the types and copy numbers of CDS genes, tRNAs and rRNAs of A1 and FN15 were identical. The cox1 gene has two copies and the trnP gene has one copy in A1, FN15 and ROC22 three lines, while there is only one copy of cox1 and two copies of trnP in S. officinarum Khon Kaen 3. In addition, S. officinarum Khon Kaen 3 have no nad1 gene and rps7 gene. 100 sequence repeats, 38 SSRs and 444 RNA editing sites in A1 mt genome were detected. Moreover, the maximum likelihood phylogenetic analysis found that A1 were more closely related to S. spp. hybrid (ROC22 and FN15) and S. officinarum (Khon Kaen 3). Herein, the complete mt genome of A1 will provide essential DNA molecular information for further phylogenetic and evolutionary analysis for Saccharum and Poaceae.

A Reappraisal of Polyploidy Events in Grasses (Poaceae) in a Rapidly Changing World

Around 80% of megaflora species became extinct at the Cretaceous–Paleogene (K–Pg) boundary. Subsequent polyploidy events drove the survival of thousands of plant species and played a significant historical role in the development of the most successful modern cereal crops. However, current and rapid global temperature change poses an urgent threat to food crops worldwide, including the world’s big three cereals: rice, wheat, and maize, which are members of the grass family, Poaceae. Some minor cereals from the same family (such as teff) have grown in popularity in recent years, but there are important knowledge gaps regarding the similarities and differences between major and minor crops, including how polyploidy affects their biological processes under natural and (a)biotic stress conditions and thus the potential to harness polyploidization attributes for improving crop climate resilience. This review focuses on the impact of polyploidy events on the Poaceae family, which includes the world’s most important food sources, and discusses the past, present, and future of polyploidy research for major and minor crops. The increasing accessibility to genomes of grasses and their wild progenitors together with new tools and interdisciplinary research on polyploidy can support crop improvement for global food security in the face of climate change.

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.

Twenty years of plant genome sequencing: achievements and challenges

Publication of the complete genome sequence of Arabidopsis thaliana, the first plant reference genome, in December 2000 heralded the beginning of the plant genome era. Over the past 20 years reference genomes have been generated for hundreds of plant species, spanning non-vascular to flowering plants. Releasing these plant genomes has dramatically advanced studies in all disciplines of plant biology. Importantly, multiple reference-level genomes have been generated for the major crops and their progenitors, enabling the creation of pan-genomes and exploration of domestication history and natural variations that can be adopted by modern crop breeding. We summarize the progress of plant genome sequencing and the challenges of sequencing more complex plant genomes and generating pan-genomes.

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.

Elbaiomy R, Chen L, Wang Z, Jiao J, Zhu J, Zhou W, Chen B, Soaud SA, Abbas M, Lin N, El-Sappah AH. Evolution analysis of FRIZZY PANICLE (FZP) orthologs explored the mutations in DNA coding sequences in the grass family (Poaceae)

FRIZZY PANICLE (FZP), an essential gene that controls spikelet differentiation and development in the grass family (Poaceae), prevents the formation of axillary bud meristems and is closely associated with crop yields. It is unclear whether the FZP gene or its orthologs were selected during the evolutionary process of grass species, which possess diverse spike morphologies. In the present study, we adopted bioinformatics methods for the evolutionary analysis of FZP orthologs in species of the grass family. Thirty-five orthologs with protein sequences identical to that of the FZP gene were identified from 29 grass species. Analysis of conserved domains revealed that the AP2/ERF domains were highly conserved with almost no amino acid mutations. However, species of the tribe Triticeae, genus Oryza, and C4 plants exhibited more significant amino acid mutations in the acidic C-terminus region. Results of the phylogenetic analysis showed that the 29 grass species could be classified into three groups, namely, Triticeae, Oryza, and C4 plants. Within the Triticeae group, the FZP genes originating from the same genome were classified into the same sub-group. When selection pressure analysis was performed, significant positive selection sites were detected in species of the Triticeae and Oryza groups. Our results show that the FZP gene was selected during the grass family's evolutionary process, and functional divergence may have already occurred among the various species. Therefore, researchers investigating the FZP gene's functions should take note of the possible presence of various roles in other grass species.

Identification of potential MTAs and candidate genes for juice quality- and yield-related traits in Saccharum clones: a genome-wide association and comparative genomic study

Sugarcane is an economically important commercial crop which provides raw material for the production of sugar, jaggery, bioethanol, biomass and other by-products. Sugarcane breeding till today heavily relies on conventional breeding approaches which is time consuming, laborious and costly. Integration of marker-assisted selection (MAS) in sugarcane genetic improvement programs for difficult to select traits like sucrose content, resistance to pests and diseases and tolerance to abiotic stresses will accelerate varietal development. In the present study, association mapping approach was used to identify QTLs and genes associated with sucrose and other important yield-contributing traits. A mapping panel of 110 diverse sugarcane genotypes and 148 microsatellite primers were used for structured association mapping study. An optimal subpopulation number (ΔK) of 5 was identified by structure analysis. GWAS analysis using TASSEL identified a total of 110 MTAs which were localized into 27 QTLs by GLM and MLM (Q + K, PC + K) approaches. Among the 24 QTLs sequenced, 12 were able to identify potential candidate genes, viz., starch branching enzyme, starch synthase 4, sugar transporters and G3P-DH related to carbohydrate metabolism and hormone pathway-related genes ethylene insensitive 3-like 1, reversion to ethylene sensitive1-like, and auxin response factor associated to juice quality- and yield-related traits. Six markers, NKS 5_185, SCB 270_144, SCB 370_256, NKS 46_176 and UGSM 648_245, associated with juice quality traits and marker SMC31CUQ_304 associated with NMC were validated and identified as significantly associated to the traits by one-way ANOVA analysis. In conclusion, 24 potential QTLs identified in the present study could be used in sugarcane breeding programs after further validation in larger population. The candidate genes from carbohydrate and hormone response pathway presented in this study could be manipulated with genome editing approaches to further improve sugarcane crop.

Chromosome-scale and haplotype-resolved genome assembly of a tetraploid potato cultivar

Potato is the most widely produced tuber crop worldwide. However, reconstructing the four haplotypes of its autotetraploid genome remained an unsolved challenge. Here, we report the 3.1 Gb haplotype-resolved (at 99.6% precision), chromosome-scale assembly of the potato cultivar ‘Otava’ based on high-quality long reads, single-cell sequencing of 717 pollen genomes and Hi-C data. Unexpectedly, ~50% of the genome was identical-by-descent due to recent inbreeding, which was contrasted by highly abundant structural rearrangements involving ~20% of the genome. Among 38,214 genes, only 54% were present in all four haplotypes with an average of 3.2 copies per gene. Taking the leaf transcriptome as an example, 11% of the genes were differently expressed in at least one haplotype, where 25% of them were likely regulated through allele-specific DNA methylation. Our work sheds light on the recent breeding history of potato, the functional organization of its tetraploid genome and has the potential to strengthen the future of genomics-assisted breeding.

Haplotype-resolved genome assembly of the tetraploid potato cultivar ‘Otava’ sheds light on functional organization of the tetraploid genome and provides the potential for genomics-assisted breeding.

Allele-aware chromosome-scale assembly of the allopolyploid genome of hexaploid Ma bamboo (Dendrocalamus latiflorus Munro)

Dendrocalamus latiflorus Munro is a woody clumping bamboo with rapid shoot growth. Both genetic transformation and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing techniques are available for D. latiflorus, enabling reverse genetic approaches. Thus, D. latiflorus has the potential to be a model bamboo species. However, the genome sequence of D. latiflorus has remained unreported due to its polyploidy and large genome size. Here, we sequenced the D. latiflorus genome and assembled it into three allele-aware subgenomes (AABBCC), representing the largest genome of a major bamboo species. We assembled 70 allelic chromosomes (2, 737 Mb) for hexaploid D. latiflorus using both single-molecule sequencing from the Pacific Biosciences (PacBio) Sequel platform and chromosome conformation capture sequencing (Hi-C). Repetitive sequences comprised 52.65% of the D. latiflorus genome. We annotated 135 231 protein-coding genes in the genome based on transcriptomes from eight different tissues. Transcriptome sequencing using RNA-Seq and PacBio single-molecule real-time long-read isoform sequencing revealed highly differential alternative splicing (AS) between non-abortive and abortive shoots, suggesting that AS regulates the abortion rate of bamboo shoots. This high-quality hexaploid genome and comprehensive strand-specific transcriptome datasets for this Poaceae family member will pave the way for bamboo research using D. latiflorus as a model species.

Divergence and conservation of defensins and lipid transfer proteins (LTPs) from sugarcane wild species and modern cultivar genomes

Plant defensins and lipid transfer proteins (LTPs) constitute a large and evolutionarily diverse family of antimicrobial peptides. Defensins and LTPs are two pathogenesis-related proteins (PR proteins) whose characterization may help to uncover aspects about the sugarcane response to pathogens attack. LTPs have also been investigated for their participation in the response to different types of stress. Despite the important roles of defensins and LTPs in biotic and abiotic stresses, scarce knowledge is found about these proteins in sugarcane. By using bioinformatics approaches, we characterized defensins and LTPs in the sugarcane wild species and modern cultivar genomes. The identification of defensins and LTPs showed that all five defensins groups and eight of the nine LTPs have their respective genes loci, although some was only identified in the cultivar genome. Phylogenetic analysis showed that defensins appear to be more conserved among groups of plants than LTPs. Some defensins and LTPs showed opposite expression during pathogenic and benefic bacterial interactions. Interestingly, the expression of defensins and LTPs in shoots and roots was completely different in plants submitted to benefic bacteria or water depletion. Finally, the modeling and comparison of isoforms of LTPs and defensins in wild species and cultivars revealed a high conservation of tertiary structures, with variation of amino acids in different regions of proteins, which could impact their antimicrobial activity. Our data contributed to the characterization of defensins and LTPs in sugarcane and provided new elements for understanding the involvement of these proteins in sugarcane response to different types of stress.

Characterization of Repetitive DNA in Saccharum officinarum and Saccharum spontaneum by Genome Sequencing and Cytological Assays

In most plant species, DNA repeated elements such as satellites and retrotransposons are composing the majority of their genomes. Saccharum officinarum (2n = 8x = 80) and S. spontaneum (2n = 40-128) are the two fundamental donors of modern sugarcane cultivars. These two species are polyploids with large genome sizes and are enriched in repetitive elements. In this work, we adopted a de novo strategy to isolate highly repetitive and abundant sequences in S. officinarum LA Purple and S. spontaneum SES208. The findings obtained from alignment to the genome assemblies revealed that the vast majority of the repeats (97.9% in LA Purple and 96.5% in SES208) were dispersed in the respective genomes. Fluorescence in situ hybridization assays were performed on 27 representative repeats to investigate their distributions and abundances. The results showed that the copies of some highly repeated sequences, including rDNA and centromeric or telomeric repeats, were underestimated in current genome assemblies. The analysis of the raw read mapping strategy showed more copy numbers for all studied repeats, suggesting that copy number underestimation is common for highly repeated sequences in current genome assemblies of LA Purple and SES208. In addition, the data showed that the centromeric retrotransposons in all SES208 centromeres were absent in certain S. spontaneum clones with different ploidies. This rapid turnover of centromeric DNA in sugarcane provides new clues regarding the pattern of centromeric retrotransposon formation and accumulation.

The highly continuous reference genome of a leaf-chimeric red pineapple (Ananas comosus var. bracteatus f. tricolor) provides insights into elaboration of leaf color

Ananas comosus var. bracteatus f. tricolor (GL1) is a red pineapple accession whose mostly green leaves with chimeric white leaf margins turn red in spring and autumn and during flowering. It is an important ornamental plant and ideal plant research model for anthocyanin metabolism, chimeric leaf development, and photosynthesis. Here, we generated a highly contiguous chromosome-scale genome assembly for GL1 and compared it with other 3 published pineapple assemblies (var. comosus accessions MD2 and F153, and var. bracteatus accession CB5). The GL1 assembly has a total size of ∼461 Mb, with a contig N50 of ∼2.97 Mb and Benchmarking Universal Single-Copy Ortholog score of 97.3%. More than 99% of the contigs are anchored to 25 pseudochromosomes. Compared with the other 3 published pineapple assemblies, the GL1 assembly was confirmed to be more continuous. Our evolutionary analysis showed that the Bromeliaceae and Poaceae diverged from their nearest common ancestor ∼82.36 million years ago (MYA). Population structure analysis showed that while GL1 has not undergone admixture, bracteatus accession CB5 has resulted from admixture of 3 species of Ananas. Through classification of orthogroups, analysis of genes under positive selection, and analysis of presence/absence variants, we identified a series of genes related to anthocyanin metabolism and development of chimeric leaves. The structure and evolution of these genes were compared among the published pineapple assemblies with reveal candidate genes for these traits. The GL1 genome assembly and its comparisons with other 3 pineapple genome assemblies provide a valuable resource for the genetic improvement of pineapple and serve as a model for understanding the genomic basis of important traits in different pineapple varieties and other pan-cereal crops.

Chromosome-specific painting unveils chromosomal fusions and distinct allopolyploid species in the Saccharum complex

Karyotypes provide key cytogenetic information on the phylogenetic relationships and evolutionary origins in related eukaryotic species. Despite our knowledge of the chromosome numbers of sugarcane and its wild relatives, the chromosome composition and evolution among the species in the Saccharum complex have been elusive owing to the complex polyploidy and the large numbers of chromosomes of these species. Oligonucleotide-based chromosome painting has become a powerful tool of cytogenetic studies especially for plant species with large numbers of chromosomes. We developed oligo-based chromosome painting probes for all 10 chromosomes in Saccharum officinarum (2n = 8x = 80). The 10 painting probes generated robust fluorescence in situ hybridization signals in all plant species within the Saccharum complex, including species in the genera Saccharum, Miscanthus, Narenga and Erianthus. We conducted comparative chromosome analysis using the same set of probes among species from four different genera within the Saccharum complex. Excitingly, we discovered several novel cytotypes and chromosome rearrangements in these species. We discovered that fusion from two different chromosomes is a common type of chromosome rearrangement associated with the species in the Saccharum complex. Such fusion events changed the basic chromosome number and resulted in distinct allopolyploids in the Saccharum complex.

Limited allele-specific gene expression in highly polyploid sugarcane

Polyploidy is widespread in plants, allowing the different copies of genes to be expressed differently in a tissue-specific or developmentally specific way. This allele-specific expression (ASE) has been widely reported, but the proportion and nature of genes showing this characteristic have not been well defined. We now report an analysis of the frequency and patterns of ASE at the whole-genome level in the highly polyploid sugarcane genome. Very high depth whole-genome sequencing and RNA sequencing revealed strong correlations between allelic proportions in the genome and in expressed sequences. This level of sequencing allowed discrimination of each of the possible allele doses in this 12-ploid genome. Most genes were expressed in direct proportion to the frequency of the allele in the genome with examples of polymorphisms being found with every possible discrete level of dose from 1:11 for single-copy alleles to 12:0 for monomorphic sites. The rarer cases of ASE were more frequent in the expression of defense-response genes, as well as in some processes related to the biosynthesis of cell walls. ASE was more common in genes with variants that resulted in significant disruption of function. The low level of ASE may reflect the recent origin of polyploid hybrid sugarcane. Much of the ASE present can be attributed to strong selection for resistance to diseases in both nature and domestication.

Interspecific complementation-restoration of phenotype in Arabidopsis cuc2cuc3 mutant by sugarcane CUC2 gene

BACKGROUND

In plants, a critical balance between differentiation and proliferation of stem cells at the shoot apical meristem zone is essential for proper growth. The spatiotemporal regulation of some crucial genes dictates the formation of a boundary within and around budding organs. The boundary plays a pivotal role in distinguishing one tissue type from another and provides a defined shape to the organs at their developed stage. NAM/CUC subfamily of the NAC transcription factors control the boundary formation during meristematic development.

RESULTS

Here, we have identified the CUP-SHAPED COTYLEDON (CUC) genes in sugarcane and named SsCUC2 (for the orthologous gene of CUC1 and CUC2) and SsCUC3. The phylogenetic reconstruction showed that SsCUCs occupy the CUC2 and CUC3 clade together with monocots, whereas eudicot CUC2 and CUC3 settled separately in the different clade. The structural analysis of CUC genes showed that most of the CUC3 genes were accompanied by an intron gain during eudicot divergence. Besides, the study of SsCUCs expression in the RNA-seq obtained during different stages of ovule development revealed that SsCUCs express in developing young tissues, and the expression of SsCUC2 is regulated by miR164. We also demonstrate that SsCUC2 (a monocot) could complement the cuc2cuc3 mutant phenotype of Arabidopsis (eudicot).

CONCLUSIONS

This study further supports that CUC2 has diverged in CUC1 and CUC2 during the evolution of monocots and eudicots from ancestral plants. The functional analysis of CUC expression patterns during sugarcane ovule development and ectopic expression of SsCUC2 in Arabidopsis showed that SsCUC2 has a conserved role in boundary formation. Overall, these findings improve our understanding of the functions of sugarcane CUC genes. Our results reveal the crucial functional role of CUC genes in sugarcane.

Genome assembly of alfalfa cultivar zhongmu-4 and identification of SNPs associated with agronomic traits

Alfalfa (Medicago sativa L.) is the most important legume forage crop worldwide with high nutritional value and yield. For a long time, the breeding of alfalfa was hampered by lacking reliable information on the autotetraploid genome and molecular markers linked to important agronomic traits. We herein reported the de novo assembly of the allele-aware chromosome-level genome of Zhongmu-4, a cultivar widely cultivated in China, and a comprehensive database of genomic variations based on resequencing of 220 germplasms. Approximate 2.74 Gb contigs (N50 of 2.06 Mb), accounting for 88.39% of the estimated genome, were assembled, and 2.56 Gb contigs were anchored to 32 pseudo-chromosomes. A total of 34,922 allelic genes were identified from the allele-aware genome. We observed the expansion of gene families, especially those related to the nitrogen metabolism, and the increase of repetitive elements including transposable elements, which probably resulted in the increase of Zhongmu-4 genome compared with Medicago truncatula. Population structure analysis revealed that the accessions from Asia and South America had relatively lower genetic diversity than those from Europe, suggesting that geography may influence alfalfa genetic divergence during local adaption. Genome-wide association studies identified 101 single nucleotide polymorphisms (SNPs) associated with 27 agronomic traits. Two candidate genes were predicted to be correlated with fall dormancy and salt response. We believe that the allele-aware chromosome-level genome sequence of Zhongmu-4 combined with the resequencing data of the diverse alfalfa germplasms will facilitate genetic research and genomics-assisted breeding in variety improvement of alfalfa.

Signatures of selection in recently domesticated macadamia

Macadamia is a high value nut crop that is recently domesticated, ideal for testing the effect of artificial selection. Here, we sequence the genome of Hawaiian cultivar ‘Kau’ and assemble into 794 Mb in 14 pseudo-chromosomes with 37,728 genes. Genome analysis reveals a whole-genome duplication event, occurred 46.8 million years ago. Gene expansions occurred in gene families involves in fatty acid biosynthesis. Gene duplication of MADS-Box transcription factors in proanthocyanidin biosynthesis are relevant for seed coat development. Genome re-sequencing of 112 accessions reveals the origin of Hawaiian cultivars from Mount Bauple in southeast Queensland in Australia. Selective sweeps are detected in macadamia cultivars, including genes involved in fatty acid biosynthesis, seed coat development, and heat stress response. Such strong effects of artificial selection in few generations reveals the genomic basis for ‘one-step operation’ for clonal crop domestication. The knowledge gained could accelerate domestication of new crops from wild species.

Macadamia is a recently domesticated nut crop. Here, the authors report the genome assembly of Hawaiian cultivar ‘Kau’ and conduct population genomic analyses to reveal the origin of Hawaiian cultivars and the genomic basis for one-step operation for the clonal crop domestication.

Genome-wide identification and expression analysis of the coronatine-insensitive 1 (COI1) gene family in response to biotic and abiotic stresses in Saccharum

Background

The coronatine insensitive 1 (COI1) gene is the core member of jasmonate signaling pathway, which is closely related to plant biotic and abiotic resistance. However, there have been no reports on COI1 in sugarcane (Sacharum spp.). Hence, systematically investigating the characteristics of the COI1 multigene family in sugarcane can provide a means to study and manipulate the jasmonic acid signaling pathway.

Results

A total of 156 COI1 proteins were obtained from the genomes of 19 land plants, while none were obtained from five algae species. A phylogenetic tree demonstrated that these COI1 proteins were classified into four groups, while 31 proteins of SsCOI1 from Saccharum spontaneum, SbCOI1 from Sorghum bicolor, and ShCOI1 from Saccharum spp. hybrid cultivar R570 clustered into three groups. Synteny analysis and duplication patterns revealed that COI1 genes expanded through various genome replication events and could have experienced strong purifying selective pressure during evolution in S. spontaneum, S. bicolor, and R570. An investigation of cis-acting elements suggests that COI1 genes may be involved in plant growth and development and response to various stresses. Expression analysis implied that 21 SsCOI1 genes were constitutively expressed, and had positive responses to drought, cold, and Sporisorium scitamineum stresses with different expression patterns. Among them, seven SsCOI1 haplotype genes may play different roles in response to methyl jasmonate. Furthermore, the ShCOI1–4, ShCOI1–5, and ShCOI1–6 genes were cloned from Saccharum spp. hybrid cultivar ROC22. Real-time quantitative PCR (RT-qPCR) analysis demonstrated that these three ShCOI1 genes had divergent expression profiles in response to salicylic acid, abscisic acid, polyethylene glycol, cold, and S. scitamineum.

Conclusions

These results suggest that COI1 genes may act in sugarcane growth, development, and response to various stresses via different regulatory mechanisms, which laying a foundation for the functional identification of the sugarcane COI1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08255-0.

Peripheral structures in unlabelled trees and the accumulation of subgenomes in the evolution of polyploids

Many angiosperms have undergone some series of polyploidization events over the course of their evolutionary history. In these genomes, especially those resulting from multiple autopolyploidization, it may be relatively easy to recognize all the ξ sets of n homeologous chromosomes, but it is much harder, if not impossible, to partition these chromosomes into n subgenomes, each representing one distinct genomic component of ξ chromosomes making up the original polyploid. Thus, if we wish to infer the polyploidization history of the genome, we could make use of all the gene trees inferred from the genes in one set of homeologous chromosomes to construct a consensus tree, but there is no evident way of combining the trees from the ξ different sets, because we have no labelling of the chromosomes that is known to be consistent across these sets. We suggest here that lacking a consistent leaf-labelling, the topological structure of the trees may display sufficient resemblance so that a higher level consensus could be revealing of evolutionary history. This would be especially true of the peripheral structures of the tree, likely representing events that occurred more recently and have thus been less obscured by subsequent evolutionary processes. Here, we present a statistical test to assess whether the subgenomes in a polyploid genome could have been added one at a time. The null hypothesis is that the accumulation of chromosomes follows a stochastic process in which transition from one generation to the next is through randomly choosing an edge, and then subdividing this edge in order to link the new internal vertex to a new external vertex. We analyze the probability distributions of a number of peripheral tree substructures, namely leaf- or terminal-pairs, triples and quadruples, arising from this stochastic process, in terms of some exact recurrences. We propose some conjectures regarding the asymptotic behaviours of these distributions. Applying our analysis to a sugarcane genome, we demonstrate that it is unlikely that the accumulation of subgenomes has occurred one at a time.

GreeNC 2.0: a comprehensive database of plant long non-coding RNAs

The Green Non-Coding Database (GreeNC) is one of the reference databases for the study of plant long non-coding RNAs (lncRNAs). Here we present our most recent update where 16 species have been updated, while 78 species have been added, resulting in the annotation of more than 495 000 lncRNAs. Moreover, sequence clustering was applied providing information about sequence conservation and gene families. The current version of the database is available at: http://greenc.sequentiabiotech.com/wiki2/Main_Page.

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.

Quantitative Trait Loci Mapping and Development of KASP Marker Smut Screening Assay Using High-Density Genetic Map and Bulked Segregant RNA Sequencing in Sugarcane (Saccharum spp.)

Sugarcane is one of the most important industrial crops globally. It is the second largest source of bioethanol, and a major crop for biomass-derived electricity and sugar worldwide. Smut, caused by Sporisorium scitamineum, is a major sugarcane disease in many countries, and is managed by smut-resistant varieties. In China, smut remains the single largest constraint for sugarcane production, and consequently it impacts the value of sugarcane as an energy feedstock. Quantitative trait loci (QTLs) associated with smut resistance and linked diagnostic markers are valuable tools for smut resistance breeding. Here, we developed an F1 population (192 progeny) by crossing two sugarcane varieties with contrasting smut resistance and used for genome-wide single nucleotide polymorphism (SNP) discovery and mapping, using a high-throughput genotyping method called "specific locus amplified fragment sequencing (SLAF-seq) and bulked-segregant RNA sequencing (BSR-seq). SLAF-seq generated 148,500 polymorphic SNP markers. Using SNP and previously identified SSR markers, an integrated genetic map with an average 1.96 cM marker interval was produced. With this genetic map and smut resistance scores of the F1 individuals from four crop years, 21 major QTLs were mapped, with a phenotypic variance explanation (PVE) > 8.0%. Among them, 10 QTLs were stable (repeatable) with PVEs ranging from 8.0 to 81.7%. Further, four QTLs were detected based on BSR-seq analysis. aligning major QTLs with the genome of a sugarcane progenitor Saccharum spontaneum, six markers were found co-localized. Markers located in QTLs and functional annotation of BSR-seq-derived unigenes helped identify four disease resistance candidate genes located in major QTLs. 77 SNPs from major QTLs were then converted to Kompetitive Allele-Specific PCR (KASP) markers, of which five were highly significantly linked to smut resistance. The co-localized QTLs, candidate resistance genes, and KASP markers identified in this study provide practically useful tools for marker-assisted sugarcane smut resistance breeding.