A chromosome-scale assembly of the model desiccation tolerant grass Oropetium thomaeum

Development of a Reference Transcriptome and Identification of Differentially Expressed Genes Linked to Salt Stress in Salt Marsh Grass (Sporobolus alterniflorus) along Delaware Coastal Regions

Salt marsh grass (Sporobolus alterniflorus) plays a crucial role in Delaware coastal regions by serving as a physical barrier between land and water along the inland bays and beaches. This vegetation helps to stabilize the shoreline and prevent erosion, protecting the land from the powerful forces of the waves and tides. In addition to providing a physical barrier, salt marsh grass is responsible for filtering nutrients in the water, offering an environment for aquatic species and presenting a focal point of study for high salt tolerance in plants. As seawater concentrations vary along the Delaware coast from low to medium to high salinity, our study seeks to identify the impact of salt tolerance in marsh grass and to identify genes associated with salt tolerance levels. We developed more than 211,000 next-generation-sequencing (Illumina) transcriptomic reads to create a reference transcriptome from low-, medium-, and high-salinity marsh grass leaf samples collected from the Delaware coastline. Contiguous sequences were annotated based on a homology search using BLASTX against rice (Oryza sativa), foxtail millet (Setaria italica), and non-redundant species within the Viridiplantae database. Additionally, we identified differentially expressed genes related to salinity stress as candidates for salt stress qPCR analysis. The data generated from this study may help to elucidate the genetic signatures and physiological responses of plants to salinity stress, thereby offering valuable insight into the use of innovative approaches for gene expression studies in crops that are less salt tolerant.

Convergent evolution of desiccation tolerance in grasses

Desiccation tolerance has evolved repeatedly in plants as an adaptation to survive extreme environments. Plants use similar biophysical and cellular mechanisms to survive life without water, but convergence at the molecular, gene and regulatory levels remains to be tested. Here we explore the evolutionary mechanisms underlying the recurrent evolution of desiccation tolerance across grasses. We observed substantial convergence in gene duplication and expression patterns associated with desiccation. Syntenic genes of shared origin are activated across species, indicative of parallel evolution. In other cases, similar metabolic pathways are induced but using different gene sets, pointing towards phenotypic convergence. Species-specific mechanisms supplement these shared core mechanisms, underlining the complexity and diversity of evolutionary adaptations to drought. Our findings provide insight into the evolutionary processes driving desiccation tolerance and highlight the roles of parallel and convergent evolution in response to environmental challenges.

Higher order polyploids exhibit enhanced desiccation tolerance in the grass Microchloa caffra

Desiccation tolerance evolved recurrently across diverse plant lineages to enable survival in water-limited conditions. Many resurrection plants are polyploid, and several groups have hypothesized that polyploidy contributed to the evolution of desiccation tolerance. However, due to the vast phylogenetic distance between resurrection plant lineages, the rarity of desiccation tolerance, and the prevalence of polyploidy in plants, this hypothesis has been difficult to test. Here, we surveyed natural variation in morphological, reproductive, and desiccation tolerance traits across several cytotypes of a single species to test for links between polyploidy and increased resilience. We sampled multiple natural populations of the resurrection grass Microchloa caffra across an environmental gradient ranging from mesic to xeric in South Africa. We describe two distinct ecotypes of M. caffra that occupy different extremes of the environmental gradient and exhibit consistent differences in ploidy, morphological, reproductive, and desiccation tolerance traits in both field and common growth conditions. Interestingly, plants with more polyploid genomes exhibited consistently higher recovery from desiccation, were less reproductive, and were larger than plants with smaller genomes and lower ploidy. These data indicate that selective pressures in increasingly xeric sites may play a role in maintaining and increasing desiccation tolerance and are mediated by changes in ploidy.

Whole-genome sequencing of allotetraploid bermudagrass reveals the origin of Cynodon and candidate genes for salt tolerance

Bermudagrass (Cynodon dactylon) is a globally distributed, extensively used warm-season turf and forage grass with high tolerance to salinity and drought stress in alkaline environments. However, the origin of the species and genetic mechanisms for salinity tolerance in the species are basically unknown. Accordingly, we set out to study evolution divergence events in the Cynodon genome and to identify genes for salinity tolerance. We developed a 604.0 Mb chromosome-level polyploid genome sequence for bermudagrass 'A12359' (n = 18). The C. dactylon genome comprises 2 complete sets of homoeologous chromosomes, each with approximately 30 000 genes, and most genes are conserved as syntenic pairs. Phylogenetic study showed that the initial Cynodon species diverged from Oropetium thomaeum approximately 19.7-25.4 million years ago (Mya), the A and B subgenomes of C. dactylon diverged approximately 6.3-9.1 Mya, and the bermudagrass polyploidization event occurred 1.5 Mya on the African continent. Moreover, we identified 82 candidate genes associated with seven agronomic traits using a genome-wide association study, and three single-nucleotide polymorphisms were strongly associated with three salt resistance genes: RAP2-2, CNG channels, and F14D7.1. These genes may be associated with enhanced bermudagrass salt tolerance. These bermudagrass genomic resources, when integrated, may provide fundamental insights into evolution of diploid and tetraploid genomes and enhance the efficacy of comparative genomics in studying salt tolerance in Cynodon.

Identification and characterization of ACR gene family in maize for salt stress tolerance

Background

Members of the ACR gene family are commonly involved in various physiological processes, including amino acid metabolism and stress responses. In recent decades, significant progress has been made in the study of ACR genes in plants. However, little is known about their characteristics and function in maize.

Methods

In this study, ACR genes were identified from the maize genome, and their molecular characteristics, gene structure, gene evolution, gene collinearity analysis, cis-acting elements were analyzed. qRT-PCR technology was used to verify the expression patterns of the ZmACR gene family in different tissues under salt stress. In addition, Ectopic expression technique of ZmACR5 in Arabidopsis thaliana was utilized to identify its role in response to salt stress.

Results

A total of 28 ZmACR genes were identified, and their molecular characteristics were extensively described. Two gene pairs arising from segmented replication events were detected in maize, and 18 collinear gene pairs were detected between maize and 3 other species. Through phylogenetic analysis, three subgroups were revealed, demonstrating distinct divergence between monocotyledonous and dicotyledonous plants. Analysis of ZmACR cis-acting elements revealed the optional involvement of ZmACR genes in light response, hormone response and stress resistance. Expression analysis of 8 ZmACR genes under salt treatment clearly revealed their role in the response to salt stress. Ectopic overexpression of ZmACR5 in Arabidopsis notably reduced salt tolerance compared to that of the wild type under salt treatment, suggesting that ZmACR5 has a negative role in the response to salt stress.

Conclusion

Taken together, these findings confirmed the involvement of ZmACR genes in regulating salt stress and contributed significantly to our understanding of the molecular function of ACR genes in maize, facilitating further research in this field.

Genome analyses reveal population structure and a purple stigma color gene candidate in finger millet

Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a 'poor man's crop', finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21^st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

Genome-wide comparative analysis of the valine glutamine motif containing genes in four Ipomoea species

BACKGROUND

Genes with valine glutamine (VQ) motifs play an essential role in plant growth, development, and resistance to biotic and abiotic stresses. However, little information on the VQ genes in sweetpotato and other Ipomoea species is available.

RESULTS

This study identified 55, 58, 50 and 47 VQ genes from sweetpotato (I. batatas), I.triflida, I. triloba and I. nil, respectively. The phylogenetic analysis revealed that the VQ genes formed eight clades (I-VII), and the members in the same group exhibited similar exon-intron structure and conserved motifs distribution. The distribution of the VQ genes among the chromosomes of Ipomoea species was disproportional, with no VQ genes mapped on a few of each species' chromosomes. Duplication analysis suggested that segmental duplication significantly contributes to their expansion in sweetpotato, I.trifida, and I.triloba, while the segmental and tandem duplication contributions were comparable in I.nil. Cis-regulatory elements involved in stress responses, such as W-box, TGACG-motif, CGTCA-motif, ABRE, ARE, MBS, TCA-elements, LTR, and WUN-motif, were detected in the promoter regions of the VQ genes. A total of 30 orthologous groups were detected by syntenic analysis of the VQ genes. Based on the analysis of RNA-seq datasets, it was found that the VQ genes are expressed distinctly among different tissues and hormone or stress treatments. A total of 40 sweetpotato differentially expressed genes (DEGs) refer to biotic (sweetpotato stem nematodes and Ceratocystis fimbriata pathogen infection) or abiotic (cold, salt and drought) stress treatments were detected. Moreover, IbVQ8, IbVQ25 and IbVQ44 responded to the five stress treatments and were selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis.

CONCLUSIONS

Our study may provide new insights into the evolution of VQ genes in the four Ipomoea genomes and contribute to the future molecular breeding of sweetpotatoes.

Core cellular and tissue-specific mechanisms enable desiccation tolerance in Craterostigma

Resurrection plants can survive prolonged life without water (anhydrobiosis) in regions with seasonal drying. This desiccation tolerance requires the coordination of numerous cellular processes across space and time, and individual plant tissues face unique constraints related to their function. Here, we analyzed the complex, octoploid genome of the model resurrection plant Craterostigma (C. plantagineum), and surveyed spatial and temporal expression dynamics to identify genetic elements underlying desiccation tolerance. Homeologous genes within the Craterostigma genome have divergent expression profiles, suggesting the subgenomes contribute differently to desiccation tolerance traits. The Craterostigma genome contains almost 200 tandemly duplicated early light-induced proteins, a hallmark trait of desiccation tolerance, with massive upregulation under water deficit. We identified a core network of desiccation-responsive genes across all tissues, but observed almost entirely unique expression dynamics in each tissue during recovery. Roots and leaves have differential responses related to light and photoprotection, autophagy and nutrient transport, reflecting their divergent functions. Our findings highlight a universal set of likely ancestral desiccation tolerance mechanisms to protect cellular macromolecules under anhydrobiosis, with secondary adaptations related to tissue function.

Comparative analysis of the MYB gene family in seven Ipomoea species

The MYB transcription factors regulate plant growth, development, and defense responses. However, information about the MYB gene family in Ipomoea species is rare. Herein, we performed a comprehensive genome-wide comparative analysis of this gene family among seven Ipomoea species, sweet potato (I. batatas), I. trifida, I. triloba, I. nil, I. purpurea, I. cairica, and I. aquatic, and identified 296, 430, 411, 291, 226, 281, and 277 MYB genes, respectively. The identified MYB genes were classified into five types: 1R-MYB (MYB-related), 2R-MYB (R2R3-MYB), 3R-MYB (R1R2R3-MYB), 4R-MYB, and 5R-MYB, and the MYB-related or R2R3-MYB type was the most abundant MYB genes in the seven species. The Ipomoea MYB genes were classed into distinct subgroups based on the phylogenetic topology and the classification of the MYB superfamily in Arabidopsis. Analysis of gene structure and protein motifs revealed that members within the same phylogenetic group presented similar exon/intron and motif organization. The identified MYB genes were unevenly mapped on the chromosomes of each Ipomoea species. Duplication analysis indicated that segmental and tandem duplications contribute to expanding the Ipomoea MYB genes. Non-synonymous substitution (Ka) to synonymous substitution (Ks) [Ka/Ks] analysis showed that the duplicated Ipomoea MYB genes are mainly under purifying selection. Numerous cis-regulatory elements related to stress responses were detected in the MYB promoters. Six sweet potato transcriptome datasets referring to abiotic and biotic stresses were analyzed, and MYB different expression genes' (DEGs') responses to stress treatments were detected. Moreover, 10 sweet potato MYB DEGs were selected for qRT-PCR analysis. The results revealed that four responded to biotic stress (stem nematodes and Ceratocystis fimbriata pathogen infection) and six responded to the biotic stress (cold, drought, and salt). The results may provide new insights into the evolution of MYB genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

From tetraploid to diploid, a pangenomic approach to identify genes lost during synthetic diploidization of Eragrostis curvula

INTRODUCTION

In Eragrostis curvula, commonly known as weeping lovegrass, a synthetic diploidization event of the facultative apomictic tetraploid Tanganyika INTA cv. originated from the sexual diploid Victoria cv. Apomixis is an asexual reproduction by seeds in which the progeny is genetically identical to the maternal plant.

METHODS

To assess the genomic changes related to ploidy and to the reproductive mode occurring during diploidization, a mapping approach was followed to obtain the first E. curvula pangenome assembly. In this way, gDNA of Tanganyika INTA was extracted and sequenced in 2x250 Illumina pair-end reads and mapped against the Victoria genome assembly. The unmapped reads were used for variant calling, while the mapped reads were assembled using Masurca software.

RESULTS

The length of the assembly was 28,982,419 bp distributed in 18,032 contigs, and the variable genes annotated in these contigs rendered 3,952 gene models. Functional annotation of the genes showed that the reproductive pathway was differentially enriched. PCR amplification in gDNA and cDNA of Tanganyika INTA and Victoria was conducted to validate the presence/absence variation in five genes related to reproduction and ploidy. The polyploid nature of the Tanganyika INTA genome was also evaluated through the variant calling analysis showing the single nucleotide polymorphism (SNP) coverage and allele frequency distribution with a segmental allotetraploid pairing behavior.

DISCUSSION

The results presented here suggest that the genes were lost in Tanganyika INTA during the diploidization process that was conducted to suppress the apomictic pathway, affecting severely the fertility of Victoria cv.

Regulatory dynamics distinguishing desiccation tolerance strategies within resurrection grasses

Desiccation tolerance has evolved recurrently in grasses using two unique strategies of either protecting or dismantling the photosynthetic apparatus to minimize photooxidative damage under life without water (anhydrobiosis). Here, we surveyed chromatin architecture and gene expression during desiccation in two closely related grasses with distinguishing desiccation tolerance strategies to identify regulatory dynamics underlying these unique adaptations. In both grasses, we observed a strong association between nearby chromatin accessibility and gene expression in desiccated tissues compared to well-watered, reflecting an unusual chromatin stability under anhydrobiosis. Integration of chromatin accessibility (ATACseq) and expression data (RNAseq) revealed a core desiccation response across these two grasses. This includes many genes with binding sites for the core seed development transcription factor ABI5, supporting the long-standing hypothesis that vegetative desiccation tolerance evolved from rewiring seed pathways. Oropetium thomaeum has a unique set of desiccation induced genes and regulatory elements associated with photoprotection, pigment biosynthesis, and response to high light, reflecting its adaptation of protecting the photosynthetic apparatus under desiccation (homoiochlorophyly). By contrast, Eragrostis nindensis has unique accessible and expressed genes related to chlorophyll catabolism, scavenging of amino acids, and hypoxia, highlighting its poikilochlorophyllous adaptations of dismantling the photosynthetic apparatus and degrading chlorophyll under desiccation. Together, our results highlight the complex regulatory and expression dynamics underlying desiccation tolerance in grasses.

Chromosome-scale genome assemblies and annotations for Poales species Carex cristatella, Carex scoparia, Juncus effusus, and Juncus inflexus

The majority of sequenced genomes in the monocots are from species belonging to Poaceae, which include many commercially important crops. Here, we expand the number of sequenced genomes from the monocots to include the genomes of 4 related cyperids: Carex cristatella and Carex scoparia from Cyperaceae and Juncus effusus and Juncus inflexus from Juncaceae. The high-quality, chromosome-scale genome sequences from these 4 cyperids were assembled by combining whole-genome shotgun sequencing of Nanopore long reads, Illumina short reads, and Hi-C sequencing data. Some members of the Cyperaceae and Juncaceae are known to possess holocentric chromosomes. We examined the repeat landscapes in our sequenced genomes to search for potential repeats associated with centromeres. Several large satellite repeat families, comprising 3.2–9.5% of our sequenced genomes, showed dispersed distribution of large satellite repeat clusters across all Carex chromosomes, with few instances of these repeats clustering in the same chromosomal regions. In contrast, most large Juncus satellite repeats were clustered in a single location on each chromosome, with sporadic instances of large satellite repeats throughout the Juncus genomes. Recognizable transposable elements account for about 20% of each of the 4 genome assemblies, with the Carex genomes containing more DNA transposons than retrotransposons while the converse is true for the Juncus genomes. These genome sequences and annotations will facilitate better comparative analysis within monocots.

Comparative analysis reveals chromosome number reductions in the evolution of African bermudagrass (Cynodon transvaalensis Burtt-Davy)

African bermudagrass (Cynodon transvaalensis Burtt-Davy) (2n = 2x = 18) belongs to the genus Cynodon, tribe Cynodonteae, subfamily Chloridoideae in the grass family Poaceae. The species is frequently crossed with common bermudagrass (Cynodon dactylon Pers.) in developing high-quality hybrid turf cultivars. Molecular resources for C. transvaalensis are scarce; thus, its genomic evolution is unknown. Recently, a linkage map consisting of 1278 markers provided a powerful tool for African bermudagrass genomic research. The objective of this study was to investigate chromosome number reduction events that resulted in the nine haploid chromosomes in this species. Tag sequences of mapped single nucleotide polymorphism markers in C. transvaalensis were compared against genome sequences of Oropetium thomaeum (L.f.) Trin. (2n = 2x = 20), a genomic model in the Cynodonteae tribe. The comparative genomic analyses revealed broad collinearity between the genomes of these two species. The analyses further revealed that two major interchromosomal rearrangements of the paleochromosome ρ12 (ρ1-ρ12-ρ1 and ρ6-ρ12-ρ6) resulted in nine chromosomes in the genome of C. transvaalensis. The findings provide novel information regarding the formation of the initial diploid species in the Cynodon genus.

Genomic insights into the origin, adaptive evolution, and herbicide resistance of Leptochloa chinensis, a devastating tetraploid weedy grass in rice fields

Chinese sprangletop (Leptochloa chinensis), belonging to the grass subfamily Chloridoideae, is one of the most notorious weeds in rice ecosystems. Here, we report a chromosome-scale reference genome assembly and a genomic variation map of the tetraploid L. chinensis. The L. chinensis genome is derived from two diploid progenitors that diverged ∼10.9 million years ago, and its two subgenomes display neither fractionation bias nor overall gene expression dominance. Comparative genomic analyses reveal substantial genome rearrangements in L. chinensis after its divergence from the common ancestor of Chloridoideae and, together with transcriptome profiling, demonstrate the important contribution of tetraploidization to the gene sources for the herbicide resistance of L. chinensis. Population genomic analyses of 89 accessions from China reveal that L. chinensis accessions collected from southern/southwestern provinces have substantially higher nucleotide diversity than those from the middle and lower reaches of the Yangtze River, suggesting that L. chinensis spread in China from the southern/southwestern provinces to the middle and lower reaches of the Yangtze River. During this spread, L. chinensis developed significantly increased herbicide resistance, accompanied by the selection of numerous genes involved in herbicide resistance. Taken together, our study generated valuable genomic resources for future fundamental research and agricultural management of L. chinensis, and provides significant new insights into the herbicide resistance as well as the origin and adaptive evolution of L. chinensis.

High-quality chromosome-scale de novo assembly of the Paspalum notatum 'Flugge' genome

BACKGROUND

Paspalum notatum 'Flugge' is a diploid with 20 chromosomes (2n = 20) multi-purpose subtropical herb native to South America and has a high ecological significance. It is currently widely planted in tropical and subtropical regions. Despite the gene pool of P. notatum 'Flugge' being unearthed to a large extent in the past decade, no details about the genomic information of relevant species in Paspalum have been reported. In this study, the complete genome information of P. notatum was established and annotated through sequencing and de novo assembly of its genome.

RESULTS

The latest PacBio third-generation HiFi assembly and sequencing revealed that the genome size of P. notatum 'Flugge' is 541 M. The assembly result is the higher index among the genomes of the gramineous family published so far, with a contig N50 = 52Mbp, scaffold N50 = 49Mbp, and BUSCOs = 98.1%, accounting for 98.5% of the estimated genome. Genome annotation revealed 36,511 high-confidence gene models, thus providing an important resource for future molecular breeding and evolutionary research. A comparison of the genome annotation results of P. notatum 'Flugge' with other closely related species revealed that it had a close relationship with Zea mays but not close compared to Brachypodium distachyon, Setaria viridis, Oryza sativa, Puccinellia tenuiflora, Echinochloa crusgalli. An analysis of the expansion and contraction of gene families suggested that P. notatum 'Flugge' contains gene families associated with environmental resistance, increased reproductive ability, and molecular evolution, which explained its excellent agronomic traits.

CONCLUSION

This study is the first to report the high-quality chromosome-scale-based genome of P. notatum 'Flugge' assembled using the latest PacBio third-generation HiFi sequencing reads. The study provides an excellent genetic resource bank for gramineous crops and invaluable perspectives regarding the evolution of gramineous plants.

A comparative genomics examination of desiccation tolerance and sensitivity in two sister grass species

This is a significant sister group contrast comparative study of the underpinning genomics and evolution of desiccation tolerance (DT), a critical trait in the evolution of land plants. Our results revealed that the DT grass Sporobolus stapfianus is transcriptionally primed to tolerate a dehydration/desiccation event and that the desiccation response in the DT S. stapfianus is distinct from the water stress response of the desiccation-sensitive Sporobolus pyramidalis. Our results also show that the desiccation response is largely unique, indicating a recent evolution of this trait within the angiosperms, and that inhibition of senescence during dehydration is likely critical in rendering a plant desiccation tolerant.

Desiccation tolerance is an ancient and complex trait that spans all major lineages of life on earth. Although important in the evolution of land plants, the mechanisms that underlay this complex trait are poorly understood, especially for vegetative desiccation tolerance (VDT). The lack of suitable closely related plant models that offer a direct contrast between desiccation tolerance and sensitivity has hampered progress. We have assembled high-quality genomes for two closely related grasses, the desiccation-tolerant Sporobolus stapfianus and the desiccation-sensitive Sporobolus pyramidalis. Both species are complex polyploids; S. stapfianus is primarily tetraploid, and S. pyramidalis is primarily hexaploid. S. pyramidalis undergoes a major transcriptome remodeling event during initial exposure to dehydration, while S. stapfianus has a muted early response, with peak remodeling during the transition between 1.5 and 1.0 grams of water (gH₂O) g⁻¹ dry weight (dw). Functionally, the dehydration transcriptome of S. stapfianus is unrelated to that for S. pyramidalis. A comparative analysis of the transcriptomes of the hydrated controls for each species indicated that S. stapfianus is transcriptionally primed for desiccation. Cross-species comparative analyses indicated that VDT likely evolved from reprogramming of desiccation tolerance mechanisms that evolved in seeds and that the tolerance mechanism of S. stapfianus represents a recent evolution for VDT within the Chloridoideae. Orthogroup analyses of the significantly differentially abundant transcripts reconfirmed our present understanding of the response to dehydration, including the lack of an induction of senescence in resurrection angiosperms. The data also suggest that failure to maintain protein structure during dehydration is likely critical in rendering a plant desiccation sensitive.

Genome-wide comparative analysis of the nucleotide-binding site-encoding genes in four Ipomoea species

The nucleotide-binding site (NBS)-encoding gene is a major type of resistance (R) gene, and its diverse evolutionary patterns were analyzed in different angiosperm lineages. Until now, no comparative studies have been done on the NBS encoding genes in Ipomoea species. In this study, various numbers of NBS-encoding genes were identified across the whole genome of sweet potato (Ipomoea batatas) (#889), Ipomoea trifida (#554), Ipomoea triloba (#571), and Ipomoea nil (#757). Gene analysis showed that the CN-type and N-type were more common than the other types of NBS-encoding genes. The phylogenetic analysis revealed that the NBS-encoding genes formed three monophyletic clades: CNL, TNL, and RNL, which were distinguished by amino acid motifs. The distribution of the NBS-encoding genes among the chromosomes was non-random and uneven; 83.13, 76.71, 90.37, and 86.39% of the genes occurred in clusters in sweet potato, I. trifida, I. triloba, and I. nil, respectively. The duplication pattern analysis reveals the presence of higher segmentally duplicated genes in sweet potatoes than tandemly duplicated ones. The opposite trend was found for the other three species. A total of 201 NBS-encoding orthologous genes were found to form synteny gene pairs between any two of the four Ipomea species, suggesting that each of the synteny gene pairs was derived from a common ancestor. The gene expression patterns were acquired by analyzing using the published datasets. To explore the candidate resistant genes in sweet potato, transcriptome analysis has been carried out using two resistant (JK20 and JK274) and susceptible cultivars (Tengfei and Santiandao) of sweet potato for stem nematodes and Ceratocystis fimbriata pathogen, respectively. A total of 11 differentially expressed genes (DEGs) were found in Tengfei and JK20 for stem nematodes and 19 DEGs in Santiandao and JK274 for C. fimbriata. Moreover, six DEGs were further selected for quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis, and the results were consistent with the transcriptome analysis. The results may provide new insights into the evolution of NBS-encoding genes in the Ipomoea genome and contribute to the future molecular breeding of sweet potatoes.

Systems biology of resurrection plants

Plant species that exhibit vegetative desiccation tolerance can survive extreme desiccation for months and resume normal physiological activities upon re-watering. Here we survey the recent knowledge gathered from the sequenced genomes of angiosperm and non-angiosperm desiccation-tolerant plants (resurrection plants) and highlight some distinct genes and gene families that are central to the desiccation response. Furthermore, we review the vast amount of data accumulated from analyses of transcriptomes and metabolomes of resurrection species exposed to desiccation and subsequent rehydration, which allows us to build a systems biology view on the molecular and genetic mechanisms of desiccation tolerance in plants.

The first draft genome of feather grasses using SMRT sequencing and its implications in molecular studies of Stipa

The Eurasian plant Stipa capillata is the most widespread species within feather grasses. Many taxa of the genus are dominants in steppe plant communities and can be used for their classification and in studies related to climate change. Moreover, some species are of economic importance mainly as fodder plants and can be used for soil remediation processes. Although large-scale molecular data has begun to appear, there is still no complete or draft genome for any Stipa species. Thus, here we present a single-molecule long-read sequencing dataset generated using the Pacific Biosciences Sequel System. A draft genome of about 1004 Mb was obtained with a contig N50 length of 351 kb. Importantly, here we report 81,224 annotated protein-coding genes, present 77,614 perfect and 58 unique imperfect SSRs, reveal the putative allopolyploid nature of S. capillata, investigate the evolutionary history of the genus, demonstrate structural heteroplasmy of the chloroplast genome and announce for the first time the mitochondrial genome in Stipa. The assembled nuclear, mitochondrial and chloroplast genomes provide a significant source of genetic data for further works on phylogeny, hybridisation and population studies within Stipa and the grass family Poaceae.

The genome of the warm-season turfgrass African bermudagrass (Cynodon transvaalensis)

Cynodon species can be used for multiple purposes and have high economic and ecological significance. However, the genetic basis of the favorable agronomic traits of Cynodon species is poorly understood, partially due to the limited availability of genomic resources. In this study, we report a chromosome-scale genome assembly of a diploid Cynodon species, C. transvaalensis, obtained by combining Illumina and Nanopore sequencing, BioNano, and Hi-C. The assembly contains 282 scaffolds (~423.42 Mb, N50 = 5.37 Mb), which cover ~93.2% of the estimated genome of C. transvaalensis (~454.4 Mb). Furthermore, 90.48% of the scaffolds (~383.08 Mb) were anchored to nine pseudomolecules, of which the largest was 60.78 Mb in length. Evolutionary analysis along with transcriptome comparison provided a preliminary genomic basis for the adaptation of this species to tropical and/or subtropical climates, typically with dry summers. The genomic resources generated in this study will not only facilitate evolutionary studies of the Chloridoideae subfamily, in particular, the Cynodonteae tribe, but also facilitate functional genomic research and genetic breeding in Cynodon species for new leading turfgrass cultivars in the future.

Genetic Diversity of Tef [Eragrostis tef (Zucc.)Trotter] as Revealed by Microsatellite Markers

Genetic variability is the fundamental prerequisite of any crop-breeding program to develop superior cultivars. There are about 350 Eragrostis1 species, of which, tef is the only species cultivated for human consumption. Currently, the Ethiopian Biodiversity Institute (EBI) collected over five thousand tef accessions from different geographical regions, diverse in terms of climate and elevation, which are uncharacterized yet. The objective of this study was to evaluate the genetic diversity among 64 tef accessions using 10 selected polymorphic simple sequence repeats (SSRs) markers. A total of 314 alleles were detected with an average of 14.5 alleles per locus and amplicon size ranged from 90 bp-320 bp. The mean value of polymorphic information content (PIC) was 0.87, appearing polymorphic for all loci. The lowest Fst value (0.05) was recorded among the studied tef populations. The mean value of major allele frequency and the number of effective alleles were 0.33 and 3.32, respectively. The mean value of gene flow (Nm) and Shannon's information index (I) was 4.74 and 1.65, respectively. The observed (Ho) and expected (He) heterozygosities varied from 0.34 to 0.56 and from 0.58 to 0.76, respectively. The cluster analysis has grouped the 64 tef accessions into three distinct clusters based on their similarity. The PCoA analysis showed that clustering is basing on the geographical origin of accessions. Analysis of molecular variance revealed 56%, 39% and 5% of the total variation due to variation within populations, among individuals and among populations, respectively. Structure bar-plot also inferred three gene pools, but with high level of admixtures. Thus, the present study shows that the identified tef accessions could be of great interest for the initiation of a planned breeding and conservation programs.

To dry perchance to live: Insights from the genome of the desiccation-tolerant biocrust moss Syntrichia caninervis

With global climate change, water scarcity threatens whole agro/ecosystems. The desert moss Syntrichia caninervis, an extremophile, offers novel insights into surviving desiccation and heat. The sequenced S. caninervis genome consists of 13 chromosomes containing 16 545 protein-coding genes and 2666 unplaced scaffolds. Syntenic relationships within the S. caninervis and Physcomitrella patens genomes indicate the S. caninervis genome has undergone a single whole genome duplication event (compared to two for P. patens) and evidence suggests chromosomal or segmental losses in the evolutionary history of S. caninervis. The genome contains a large sex chromosome composed primarily of repetitive sequences with a large number of Copia and Gypsy elements. Orthogroup analyses revealed an expansion of ELIP genes encoding proteins important in photoprotection. The transcriptomic response to desiccation identified four structural clusters of novel genes. The genomic resources established for this extremophile offer new perspectives for understanding the evolution of desiccation tolerance in plants.

Plant Co-expression Annotation Resource: a web server for identifying targets for genetically modified crop breeding pipelines

The development of genetically modified crops (GM) includes the discovery of candidate genes through bioinformatics analysis using genomics data, gene expression, and others. Proteins of unknown function (PUFs) are interesting targets for GM crops breeding pipelines for the novelty associated with such targets and also to avoid copyright protection. One method of inferring the putative function of PUFs is by relating them to factors of interest such as abiotic stresses using orthology and co-expression networks, in a guilt-by-association manner. In this regard, we have downloaded, analyzed, and processed genomics data of 53 angiosperms, totaling 1,862,010 genes and 2,332,974 RNA. Diamond and InterproScan were used to discover 72,266 PUFs for all organisms. RNA-seq datasets related to abiotic stresses were downloaded from NCBI/GEO. The RNA-seq data was used as input to the LSTrAP software to construct co-expression networks. LSTrAP also created clusters of transcripts with correlated expression, whose members are more probably related to the molecular mechanisms associated with abiotic stresses in the plants. Orthologous groups were created (OrhtoMCL) using all 2,332,974 proteins in order to associate PUFs to abiotic stress-related clusters of co-expression and therefore infer their function in a guilt-by-association manner. A freely available web resource named "Plant Co-expression Annotation Resource" ( https://www.machado.cnptia.embrapa.br/plantannot ), Plantannot, was created to provide indexed queries to search for PUF putatively associated with abiotic stresses. The web interface also allows browsing, querying, and retrieving of public genomics data from 53 plants. We hope Plantannot to be useful for researchers trying to obtain novel GM crops resistant to climate change hazards.

Sequence-based genetic mapping of Cynodon dactylon Pers. reveals new insights into genome evolution in Poaceae

Bermudagrass (Cynodon dactylon Pers.) is an important warm-season perennial used extensively for turf, forage, soil conservation and remediation worldwide. However, limited genomic information has hindered the application of molecular tools towards understanding genome evolution and in breeding new cultivars. We genotype a first-generation selfed population derived from the tetraploid (4x = 36) ‘A12359’ using genotyping-by-sequencing. A high-density genetic map of 18 linkage groups (LGs) is constructed with 3,544 markers. Comparative genomic analyses reveal that each of nine homeologous LG pairs of C. dactylon corresponds to one of the first nine chromosomes of Oropetium thomaeum. Two nested paleo-ancestor chromosome fusions (ρ6-ρ9-ρ6, ρ2-ρ10-ρ2) may have resulted in a 12-to-10 chromosome reduction. A segmental dissemination of the paleo-chromosome ρ12 (ρ1-ρ12-ρ1, ρ6-ρ12-ρ6) leads to the 10-to-9 chromosome reduction in C. dactylon genome. The genetic map will assist in an ongoing whole genome sequence assembly and facilitate marker-assisted selection (MAS) in developing new cultivars.

Tilin Fang et al. study the genome of Bermudagrass (Cynodon dactylon Pers.). They use genotyping-by-sequencing and provide a genetic map with a 10-fold increase in genetic marker density. Comparative genomics analyses reveal chromosome rearrangements. Their work will contribute to whole genome assembly efforts which will be beneficial for developing new cultivars.

Convergent Loss of an EDS1/PAD4 Signaling Pathway in Several Plant Lineages Reveals Coevolved Components of Plant Immunity and Drought Response

Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses.

Intertwined signatures of desiccation and drought tolerance in grasses

Grasses are among the most resilient plants, and some can survive prolonged desiccation in semiarid regions with seasonal rainfall. However, the genetic elements that distinguish grasses that are sensitive versus tolerant to extreme drying are largely unknown. Here, we leveraged comparative genomic approaches with the desiccation-tolerant grass Eragrostis nindensis and the related desiccation-sensitive cereal Eragrostis tef to identify changes underlying desiccation tolerance. These analyses were extended across C4 grasses and cereals to identify broader evolutionary conservation and divergence. Across diverse genomic datasets, we identified changes in chromatin architecture, methylation, gene duplications, and expression dynamics related to desiccation in E. nindensis It was previously hypothesized that transcriptional rewiring of seed desiccation pathways confers vegetative desiccation tolerance. Here, we demonstrate that the majority of seed-dehydration-related genes showed similar expression patterns in leaves of both desiccation-tolerant and -sensitive species. However, we identified a small set of seed-related orthologs with expression specific to desiccation-tolerant species. This supports a broad role for seed-related genes, where many are involved in typical drought responses, with only a small subset of crucial genes specifically induced in desiccation-tolerant plants.

Exceptional subgenome stability and functional divergence in the allotetraploid Ethiopian cereal teff

Teff (Eragrostis tef) is a cornerstone of food security in the Horn of Africa, where it is prized for stress resilience, grain nutrition, and market value. Here, we report a chromosome-scale assembly of allotetraploid teff (variety Dabbi) and patterns of subgenome dynamics. The teff genome contains two complete sets of homoeologous chromosomes, with most genes maintaining as syntenic gene pairs. TE analysis allows us to estimate that the teff polyploidy event occurred ~1.1 million years ago (mya) and that the two subgenomes diverged ~5.0 mya. Despite this divergence, we detect no large-scale structural rearrangements, homoeologous exchanges, or biased gene loss, in contrast to many other allopolyploids. The two teff subgenomes have partitioned their ancestral functions based on divergent expression across a diverse expression atlas. Together, these genomic resources will be useful for accelerating breeding of this underutilized grain crop and for fundamental insights into polyploid genome evolution.

IsoSeq transcriptome assembly of C3 panicoid grasses provides tools to study evolutionary change in the Panicoideae

The number of plant species with genomic and transcriptomic data has been increasing rapidly. The grasses-Poaceae-have been well represented among species with published reference genomes. However, as a result the genomes of wild grasses are less frequently targeted by sequencing efforts. Sequence data from wild relatives of crop species in the grasses can aid the study of domestication, gene discovery for breeding and crop improvement, and improve our understanding of the evolution of C4 photosynthesis. Here, we used long-read sequencing technology to characterize the transcriptomes of three C3 panicoid grass species: Dichanthelium oligosanthes, Chasmanthium laxum, and Hymenachne amplexicaulis. Based on alignments to the sorghum genome, we estimate that assembled consensus transcripts from each species capture between 54.2% and 65.7% of the conserved syntenic gene space in grasses. Genes co-opted into C4 were also well represented in this dataset, despite concerns that because these genes might play roles unrelated to photosynthesis in the target species, they would be expressed at low levels and missed by transcript-based sequencing. A combined analysis using syntenic orthologous genes from grasses with published reference genomes and consensus long-read sequences from these wild species was consistent with previously published phylogenies. It is hoped that these data, targeting underrepresented classes of species within the PACMAD grasses-wild species and species utilizing C3 photosynthesis-will aid in future studies of domestication and C4 evolution by decreasing the evolutionary distance between C4 and C3 species within this clade, enabling more accurate comparisons associated with evolution of the C4 pathway.

Massive Tandem Proliferation of ELIPs Supports Convergent Evolution of Desiccation Tolerance across Land Plants

Desiccation tolerance was a critical adaptation for the colonization of land by early nonvascular plants. Resurrection plants have maintained or rewired these ancestral protective mechanisms, and desiccation-tolerant species are dispersed across the land plant phylogeny. Although common physiological, biochemical, and molecular signatures are observed across resurrection plant lineages, features underlying the recurrent evolution of desiccation tolerance are unknown. Here we used a comparative approach to identify patterns of genome evolution and gene duplication associated with desiccation tolerance. We identified a single gene family with dramatic expansion in all sequenced resurrection plant genomes and no expansion in desiccation-sensitive species. This gene family of early light-induced proteins (ELIPs) expanded in resurrection plants convergent through repeated tandem gene duplication. ELIPs are universally highly expressed during desiccation in all surveyed resurrection plants and may play a role in protecting against photooxidative damage of the photosynthetic apparatus during prolonged dehydration. Photosynthesis is particularly sensitive to dehydration, and the increased abundance of ELIPs may help facilitate the rapid recovery observed for most resurrection plants. Together, these observations support convergent evolution of desiccation tolerance in land plants through tandem gene duplication.