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

Chromosome-scale genome assembly of Zoysia japonica uncovers cold tolerance candidate genes

Zoysiagrass stands out as a crucial native turfgrass due to its exceptional abiotic stress tolerance, extensive adaptability, and high ornamental value. In this study, we generated a high-quality chromosome-level genome assembly of Compadre (COM) zoysiagrass, leveraging PacBio SMRT sequencing and Hi-C scaffolding technologies. The resulting genome assembly (312.42 Mb) is anchored on 20 chromosomes, with a Scaffold N50 of 18.72 Mb. In total, 49,074 genes and 306,768 repeat sequences were annotated in the assembled genome. The first chromosome-scale genome of Zoysia japonica 'Compadre' provides a critical genetic resource for cold-tolerant turfgrass breeding through identifying stress-responsive candidate genes. Additionally, we have successfully established a cell nucleus extraction and library construction protocol tailored for zoysiagrass ATAC-seq technology, and a total of 80 low temperature tolerance candidate genes were preliminarily identified via ATAC-seq and RNA-seq profiling, thereby initiating the exploration of turfgrass epigenomics.

Integration of ATAC-seq and RNA-seq reveals the dynamics of chromatin accessibility and gene expression in zoysiagrass response to drought

KEY MESSAGE

The 'X4' accession of zoysiagrass demonstrated superior drought tolerance compared to other accessions. Integration analysis of transcriptomics and epigenomics revealed a positive correlation between ATAC-seq peak intensity and gene expression levels. Six motifs involved in regulating drought responses were identified, which are similar to the domains of the ERF and C2H2 transcription factor families. Heterologous expression of Zja11G000860 in yeast enhanced tolerance to drought stress, allowing robust growth even at high PEG6000 concentrations. Zoysiagrass is renowned for its drought tolerance and serves as an exceptional domestic turfgrass in China. However, the changes in chromatin accessibility during drought in zoysiagrass are not well understood. We conducted a preliminary evaluation of the phenotypic changes in drought tolerance for six zoysiagrass cultivars, taking into account their growth characteristics and physiological traits under drought conditions. Additionally, we utilized an integrated multi-omics strategy, encompassing RNA sequencing (RNA-seq), Assay for Transposase Accessible Chromatin using high-throughput sequencing (ATAC-seq), and reverse transcription quantitative PCR (RT-qPCR) verification experiments, to gain deeper understanding of the chromatin accessibility patterns linked to gene expression in response to drought stress in zoysiagrass. Preliminary analysis of the trends in relative water content and proline content suggested that the variety 'X4' exhibited superior drought tolerance compared to the other five accessions. The KEGG pathway enrichment analysis revealed that zoysiagrass responded to environmental stress by regulating stress response and antioxidant defense pathways. Notably, the expression levels of genes Zja03G031540 and Zja11G000860 were significantly increased in the 'X4' zoysiagrass genotype, which exhibited improved drought tolerance, compared to the 'X1' zoysiagrass genotype with reduced drought tolerance. This study suggested that 63 high-confidence genes are related to drought stress, including Zja03G031540 and Zja11G000860. Additionally, six motifs regulating drought responses were unearthed. Furthermore, the heterologous expression of Zja11G000860 in yeast enhanced tolerance to drought stress. The study discovered a positive correlation between ATAC-seq peak intensity and gene expression levels. The expression of high-confidence genes was linked to zoysiagrass resistance evaluation and phenotypic traits, implying that these genes are involved in responding to external drought stress. This study combined ATAC-seq and RNA-seq technologies for the first time to identify drought-related gene expression in zoysiagrass, elucidating the grass adaptation to environmental stress and the regulatory mechanisms underlying stress responses, and laying the groundwork for zoysiagrass improvement and breeding.

On the way to diploidization and unexpected ploidy in the grass Sporobolus section Spartina mesopolyploids

Plant history is characterized by cyclical whole genome duplication and diploidization with important biological and ecological consequences. Here, we explore the genome history of two related iconic polyploid grasses (Sporobolus alterniflorus and S. maritimus), involved in a well-known example of neopolyploid speciation. We report particular genome dynamics where an ancestral Sporobolus genome (n = 2x = 20) duplicated 9.6-24.4 million years ago (MYA), which was followed by descending dysploidy resulting in a genome with an unexpected base chromosome number (n = 15). This diploidized genome duplicated again 2.1-6.2 MYA to form a tetraploid lineage (2n = 4x = 60), thus reshuffling the ploidy of these species previously thought hexaploids. We also elucidate the mechanism accompanying the speciation between S. maritimus (2n = 60) and S. alterniflorus (2n = 62), resulting from chromosome restructuring, and identify key adaptive genes in the corresponding regions. This represents critical findings to decipher molecular mechanisms underlying species expansion, adaptation to environmental challenge and invasiveness.

Genome editing for grass improvement and future agriculture

Grasses, including turf and forage, cover most of the earth's surface; predominantly important for land, water, livestock feed, soil, and water conservation, as well as carbon sequestration. Improved production and quality of grasses by modern molecular breeding is gaining more research attention. Recent advances in genome-editing technologies are helping to revolutionize plant breeding and also offering smart and efficient acceleration on grass improvement. Here, we reviewed all recent researches using (CRISPR)/CRISPR-associated protein (Cas)-mediated genome editing tools to enhance the growth and quality of forage and turf grasses. Furthermore, we highlighted emerging approaches aimed at advancing grass breeding program. We assessed the CRISPR-Cas effectiveness, discussed the challenges associated with its application, and explored future perspectives primarily focusing on turf and forage grasses. Despite the promising potential of genome editing in grasses, its current efficiency remains limited due to several bottlenecks, such as the absence of comprehensive reference genomes, the lack of efficient gene delivery tools, unavailability of suitable vector and delivery for grass species, high polyploidization, and multiple homoeoalleles, etc. Despite these challenges, the CRISPR-Cas system holds great potential to fully harness its benefits in grass breeding and genetics, aiming to improve and sustain the quantity and quality of turf and forage grasses.

A comprehensive assessment of photosynthetic acclimation to shade in C4 grass (Cynodon dactylon (L.) Pers.)

BACKGROUND

Light deficit in shaded environment critically impacts the growth and development of turf plants. Despite this fact, past research has predominantly concentrated on shade avoidance rather than shade tolerance. To address this, our study examined the photosynthetic adjustments of Bermudagrass when exposed to varying intensities of shade to gain an integrative understanding of the shade response of C4 turfgrass.

RESULTS

We observed alterations in photosynthetic pigment-proteins, electron transport and its associated carbon and nitrogen assimilation, along with ROS-scavenging enzyme activity in shaded conditions. Mild shade enriched Chl b and LHC transcripts, while severe shade promoted Chl a, carotenoids and photosynthetic electron transfer beyond QA- (ET0/RC, φE0, Ψ0). The study also highlighted differential effects of shade on leaf and root components. For example, Soluble sugar content varied between leaves and roots as shade diminished SPS, SUT1 but upregulated BAM. Furthermore, we observed that shading decreased the transcriptional level of genes involving in nitrogen assimilation (e.g. NR) and SOD, POD, CAT enzyme activities in leaves, even though it increased in roots.

CONCLUSIONS

As shade intensity increased, considerable changes were noted in light energy conversion and photosynthetic metabolism processes along the electron transport chain axis. Our study thus provides valuable theoretical groundwork for understanding how C4 grass acclimates to shade tolerance.

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.

Breeding for reduced seed dormancy to domesticate new grass species

Introducing new grass species into cultivation has long been proposed as beneficial to increase the sustainability and diversity of productive systems. However, wild species with potential tend to show high seed dormancy, causing slow, poor, and unsynchronized seedling emergence. Meanwhile, domesticated species, such as cereals, show lower seed dormancy, facilitating their successful establishment. In this work, we conduct a review of phenotypic variation on seed dormancy and its genetic and molecular basis. This quantitative and highly heritable trait shows phenotype plasticity which is modulated by environmental factors. The level of dormancy depends on the expression of genes associated with the metabolism and sensitivity to the hormones abscisic acid (ABA) and gibberellins (GA), along with other dormancy-specific genes. The genetic regulation of these traits is highly conserved across species. The low seed dormancy observed in cereals and some temperate forages was mostly unconsciously selected during various domestication processes. Emphasis is placed on selecting materials with low seed dormancy for warm-season forage grasses to improve their establishment and adoption. Finally, we review advances in the domestication of dallisgrass, where seed dormancy was considered a focus trait throughout the process.

Development of genomic resources for Rhodes grass (Chloris gayana), draft genome and annotated variant discovery

Genomic resources for grasses, especially warm-season grasses are limited despite their commercial and environmental importance. Here, we report the first annotated draft whole genome sequence for diploid Rhodes grass (Chloris gayana), a tropical C4 species. Generated using long read nanopore sequencing and assembled using the Flye software package, the assembled genome is 603 Mbp in size and comprises 5,233 fragments that were annotated using the GenSas pipeline. The annotated genome has 46,087 predicted genes corresponding to 92.0% of the expected genomic content present via BUSCO analysis. Gene ontology terms and repetitive elements are identified and discussed. An additional 94 individual plant genotypes originating from three diploid and two tetraploid Rhodes grass cultivars were short-read whole genome resequenced (WGR) to generate a single nucleotide polymorphism (SNP) resource for the species that can be used to elucidate inter- and intra-cultivar relationships across both ploidy levels. A total of 75,777 high quality SNPs were used to generate a phylogenetic tree, highlighting the diversity present within the cultivars which agreed with the known breeding history. Differentiation was observed between diploid and tetraploid cultivars. The WGR data were also used to provide insights into the nature and evolution of the tetraploid status of the species, with results largely agreeing with the published literature that the tetraploids are autotetraploid.

Genetic diversity and population structure of bermudagrass (Cynodon spp.) revealed by genotyping-by-sequencing

Bermudagrass (Cynodon spp.) breeding and cultivar development is hampered by limited information regarding its genetic and phenotypic diversity. To explore diversity in bermudagrass, a total of 206 Cynodon accessions consisting of 193 common bermudagrass (C. dactylon var. dactylon) and 13 African bermudagrass (C. transvaalensis) accessions of worldwide origin were assembled for genetic characterization. Genotyping-by-sequencing (GBS) was employed for genetic marker development. With a minor allele frequency of 0.05 and a minimum call rate of 0.5, a total of 37,496 raw single nucleotide polymorphisms (SNPs) were called de novo and were used in the genetic diversity characterization. Population structure analysis using ADMIXTURE revealed four subpopulations in this germplasm panel, which was consistent with principal component analysis (PCA) and phylogenetic analysis results. The first three principal components explained 15.6%, 10.1%, and 3.8% of the variance in the germplasm panel, respectively. The first subpopulation consisted of C. dactylon accessions from various continents; the second subpopulation was comprised mainly of C. transvaalensis accessions; the third subpopulation contained C. dactylon accessions primarily of African origin; and the fourth subpopulation represented C. dactylon accessions obtained from the Oklahoma State University bermudagrass breeding program. Genetic diversity parameters including Nei's genetic distance, inbreeding coefficient, and Fst statistic revealed substantial genetic variation in the Cynodon accessions, demonstrating the potential of this germplasm panel for further genetic studies and cultivar development in breeding programs.

Chromosome-level genome assembly of a high-altitude-adapted frog (Rana kukunoris) from the Tibetan plateau provides insight into amphibian genome evolution and adaptation

BACKGROUND

The high-altitude-adapted frog Rana kukunoris, occurring on the Tibetan plateau, is an excellent model to study life history evolution and adaptation to harsh high-altitude environments. However, genomic resources for this species are still underdeveloped constraining attempts to investigate the underpinnings of adaptation.

RESULTS

The R. kukunoris genome was assembled to a size of 4.83 Gb and the contig N50 was 1.80 Mb. The 6555 contigs were clustered and ordered into 12 pseudo-chromosomes covering ~ 93.07% of the assembled genome. In total, 32,304 genes were functionally annotated. Synteny analysis between the genomes of R. kukunoris and a low latitude species Rana temporaria showed a high degree of chromosome level synteny with one fusion event between chr11 and chr13 forming pseudo-chromosome 11 in R. kukunoris. Characterization of features of the R. kukunoris genome identified that 61.5% consisted of transposable elements and expansions of gene families related to cell nucleus structure and taste sense were identified. Ninety-five single-copy orthologous genes were identified as being under positive selection and had functions associated with the positive regulation of proteins in the catabolic process and negative regulation of developmental growth. These gene family expansions and positively selected genes indicate regions for further interrogation to understand adaptation to high altitude.

CONCLUSIONS

Here, we reported a high-quality chromosome-level genome assembly of a high-altitude amphibian species using a combination of Illumina, PacBio and Hi-C sequencing technologies. This genome assembly provides a valuable resource for subsequent research on R. kukunoris genomics and amphibian genome evolution in general.

A Novel Beta-Glucosidase Gene for Plant Type Was Identified by Genome-Wide Association Study and Gene Co-Expression Analysis in Widespread Bermudagrass

Bermudagrass (Cynodon spp.) is one of the most widely distributed warm-season grasses globally. The growth habits and plant type of bermudagrass are strongly associated with the applied purpose of the landscape, livestock, and eco-remediation. Therefore, persistent efforts are made to investigate the genetic basis of plant type and growth habits of bermudagrass. Here, we dissect the genetic diversity of 91 wild bermudagrass resources by genome-wide association studies (GWAS) combined with weighted gene co-expression analysis (WGCNA). This work is based on the RNA-seq data and the genome of African bermudagrass (Cynodon transvaalensis Burtt Davy). Sixteen reliable single-nucleotide polymorphisms (SNPs) in transcribed regions were identified to be associated with the plant height and IAA content in diverse bermudagrass by GWAS. The integration of the results from WGCNA indicates that beta-glucosidase 31 (CdBGLU31) is a candidate gene underlying a G/A SNP signal. Furthermore, both qRT-PCR and correlation coefficient analyses indicate that CdBGLU31 might play a comprehensive role in plant height and IAA biosynthesis and signal. In addition, we observe lower plant height in Arabidopsis bglu11 mutants (homologs of CdBGLU31). It uncovers the breeding selection history of different plant types from diverse bermudagrass and provides new insights into the molecular function of CdBGLU31 both in plant types and in IAA biosynthetic pathways.

A de novo assembled high-quality chromosome-scale Trifolium pratense genome and fine-scale phylogenetic analysis

BACKGROUND

Red clover (Trifolium pratense L.) is a diploid perennial temperate legume with 14 chromosomes (2n = 14) native to Europe and West Asia, with high nutritional and economic value. It is a very important forage grass and is widely grown in marine climates, such as the United States and Sweden. Genetic research and molecular breeding are limited by the lack of high-quality reference genomes. In this study, we used Illumina, PacBio HiFi, and Hi-C to obtain a high-quality chromosome-scale red clover genome and used genome annotation results to analyze evolutionary relationships among related species.

RESULTS

The red clover genome obtained by PacBio HiFi assembly sequencing was 423 M. The assembly quality was the highest among legume genome assemblies published to date. The contig N50 was 13 Mb, scaffold N50 was 55 Mb, and BUSCO completeness was 97.9%, accounting for 92.8% of the predicted genome. Genome annotation revealed 44,588 gene models with high confidence and 52.81% repetitive elements in red clover genome. Based on a comparison of genome annotation results, red clover was closely related to Trifolium medium and distantly related to Glycine max, Vigna radiata, Medicago truncatula, and Cicer arietinum among legumes. Analyses of gene family expansions and contractions and forward gene selection revealed gene families and genes related to environmental stress resistance and energy metabolism.

CONCLUSIONS

We report a high-quality de novo genome assembly for the red clover at the chromosome level, with a substantial improvement in assembly quality over those of previously published red clover genomes. These annotated gene models can provide an important resource for molecular genetic breeding and legume evolution studies. Furthermore, we analyzed the evolutionary relationships among red clover and closely related species, providing a basis for evolutionary studies of clover leaf and legumes, genomics analyses of forage grass, the improvement of agronomic traits.

Progress and Challenges in China Turfgrass Abiotic Stress Resistance Research

Turfgrasses are an important vehicle for urban ecology and one of the most important indicators of economy and civilization. The biological characteristics of different turfgrass species affect the productivity and quality of the turf and its potential use in landscapes, slopes, and sports fields. Cultivation and management techniques can assist turfgrasses to meet the challenges of climate change, while the development of molecular breeding will provide a broader platform for the application of turfgrasses. The turfgrass industry of China has developed considerably in the last three decades; however, there is still an objective gap with developed countries. This manuscript reviewed the research progress of turfgrass resistance breeding, analyzed the bottlenecks in the development of turfgrass resistance breeding, and put forward the strategies to cope with the bottlenecks. Our review aims to promote research and utilization of turfgrasses.

High-quality chromosome-level genome assembly of Litsea coreana L. provides insights into Magnoliids evolution and flavonoid biosynthesis

The magnoliid Litsea coreana has been the subject of a substantial amount of research owing to its production of many flavonoid metabolites, high food processing value, and a controversial phylogenetic position. For this study, we assembled a high-grade genome at the chromosome scale and annotation of L. coreana that was anchored to 12 chromosomes. The total genome was 1139.45 Mb, while the N50 scaffold was 97.18 Mb long. The analysis of phylogenetic trees constructed by different methods show that the phylogeny of Magnoliids is inconsistent, indicating that the differentiation process of monocots, eudicots, and Magnoliids still remains in dispute. An ancient whole-genome duplication (WGD) event was shown to have occurred before the Magnoliales and Laurels had differentiated. Subsequently, an independent WGD appeared in the Lauralean lineage. A total of 27 types of flavonoids were detected in all five tissues of L. coreana. Chalcone synthases (CHSs) that are responsible for production of flavonoids have been validated at the bioinformatics level. The retention of comparative genomic analyses of the CHS gene family showed that this family had contracted significantly in L. coreana. Our research further elaborated the evolution of Lauraceae and perfected the genetic basis of flavonoid biosynthesis in L. coreana. SIGNIFICANCE STATEMENT: Provides evidence that determines the evolutionary status of Magnoliids. The chalcone synthase gene family was significantly contracted in Litsea coreana.

Expression of a Cytochrome P450 Gene from Bermuda Grass Cynodon dactylon in Soybean Confers Tolerance to Multiple Herbicides

Bermuda grass (Cynodon dactylon) is notoriously difficult to control with some commonly used herbicides. We cloned a cytochrome P450 gene from Bermuda grass, named P450-N-Z1, which was found to confer tolerance to multiple herbicides in transgenic Arabidopsis. These herbicides include: (1) acetolactate synthase (ALS) inhibitor herbicides nicosulfuron and penoxsulam; (2) p-hydroxyphenylpyruvate dioxygenase (HPPD)-inhibiting herbicide mesotrione; (3) synthetic auxin herbicide dicamba; (4) photosynthesis inhibitor bentazon. We further generated transgenic soybean plants expressing P450-N-Z1, and found that these transgenic soybean plants gained robust tolerance to nicosulfuron, flazasulfuron, and 2,4-dichlorophenoxyacetic acid (2,4-D) in greenhouse assays. A field trial demonstrated that transgenic soybean is tolerant to flazasulfuron and 2,4-D at 4-fold and 2-fold the recommended rates, respectively. Furthermore, we also demonstrated that flazasulfuron and dicamba are much more rapidly degraded in vivo in the transgenic soybean than in non-transgenic soybean. Therefore, P450-N-Z1 may be utilized for engineering transgenic crops for herbicide tolerance.

A High-Quality Haplotype-Resolved Genome of Common Bermudagrass (Cynodon dactylon L.) Provides Insights Into Polyploid Genome Stability and Prostrate Growth

Common bermudagrass (Cynodon dactylon L.) is an important perennial warm-season turfgrass species with great economic value. However, the reference genome is still deficient in C. dactylon, which severely impedes basic studies and breeding studies. In this study, a high-quality haplotype-resolved genome of C. dactylon cultivar Yangjiang was successfully assembled using a combination of multiple sequencing strategies. The assembled genome is approximately 1.01 Gb in size and is comprised of 36 pseudo chromosomes belonging to four haplotypes. In total, 76,879 protein-coding genes and 529,092 repeat sequences were annotated in the assembled genome. Evolution analysis indicated that C. dactylon underwent two rounds of whole-genome duplication events, whereas syntenic and transcriptome analysis revealed that global subgenome dominance was absent among the four haplotypes. Genome-wide gene family analyses further indicated that homologous recombination-regulating genes and tiller-angle-regulating genes all showed an adaptive evolution in C. dactylon, providing insights into genome-scale regulation of polyploid genome stability and prostrate growth. These results not only facilitate a better understanding of the complex genome composition and unique plant architectural characteristics of common bermudagrass, but also offer a valuable resource for comparative genome analyses of turfgrasses and other plant species.

Integration of the metabolome and transcriptome reveals the mechanism of resistance to low nitrogen supply in wild bermudagrass (Cynodon dactylon (L.) Pers.) roots

BACKGROUND

Nitrogen (N) is an essential macronutrient that significantly affects turf quality. Commercial cultivars of bermudagrass (Cynodon dactylon (L.) Pers.) require large amounts of nitrogenous fertilizer. Wild bermudagrass germplasm from natural habitats with poor nutrition and diverse N distributions is an important source for low-N-tolerant cultivated bermudagrass breeding. However, the mechanisms underlying the differences in N utilization among wild germplasm resources of bermudagrass are not clear.

RESULTS

To clarify the low N tolerance mechanism in wild bermudagrass germplasm, the growth, physiology, metabolome and transcriptome of two wild accessions, C291 (low-N-tolerant) and C716 (low-N-sensitive), were investigated. The results showed that root growth was less inhibited in low-N-tolerant C291 than in low-N-sensitive C716 under low N conditions; the root dry weight, soluble protein content and free amino acid content of C291 did not differ from those of the control, while those of C716 were significantly decreased. Down-regulation of N acquisition, primary N assimilation and amino acid biosynthesis was less pronounced in C291 than in C716 under low N conditions; glycolysis and the tricarboxylic acid (TCA) cycle pathway were also down-regulated, accompanied by a decrease in the biosynthesis of amino acids; strikingly, processes such as translation, biosynthesis of the structural constituent of ribosome, and the expression of individual aminoacyl-tRNA synthetase genes, most of genes associated with ribosomes related to protein synthesis were all up-regulated in C291, but down-regulated in C716.

CONCLUSIONS

Overall, low-N-tolerant wild bermudagrass tolerated low N nutrition by reducing N primary assimilation and amino acid biosynthesis, while promoting the root protein synthesis process and thereby maintaining root N status and normal growth.

Genome-Wide Analysis of the HSP20 Gene Family and Expression Patterns of HSP20 Genes in Response to Abiotic Stresses in Cynodon transvaalensis

African bermudagrass (Cynodon transvaalensis Burtt-Davy) is an important warm-season turfgrass and forage grass species. Heat shock protein 20 (HSP20) is a diverse, ancient, and important protein family. To date, HSP20 genes have not been characterized genome-widely in African bermudagrass. Here, we confirmed 41 HSP20 genes in African bermudagrass genome. On the basis of the phylogenetic tree and cellular locations, the HSP20 proteins were classified into 12 subfamilies. Motif composition was consistent with the phylogeny. Moreover, we identified 15 pairs of paralogs containing nine pairs of tandem duplicates and six pairs of WGD/segmental duplicates of HSP20 genes. Unsurprisingly, the syntenic genes revealed that African bermudagrass had a closer evolutionary relationship with monocots (maize and rice) than dicots (Arabidopsis and soybean). The expression patterns of HSP20 genes were identified with the transcriptome data under abiotic stresses. According to the expression profiles, HSP20 genes could be clustered into three groups (Groups I, II, and III). Group I was the largest, and these genes were up-regulated in response to heat stress as expected. In Group II, one monocot-specific HSP20, CtHSP20-14 maintained higher expression levels under optimum temperature and low temperature, but not high temperature. Moreover, a pair of WGD/segmental duplicates CtHSP20-9 and CtHSP20-10 were among the most conserved HSP20s across different plant species, and they seemed to be positively selected in response to extreme temperatures during evolution. A total of 938 cis-elements were captured in the putative promoters of HSP20 genes. Almost half of the cis-elements were stress responsive, indicating that the expression pattern of HSP20 genes under abiotic stresses might be largely regulated by the cis-elements. Additionally, three-dimensional structure simulations and protein-protein interaction networks were incorporated to resolve the function mechanism of HSP20 proteins. In summary, the findings fulfilled the HSP20 family analysis and could provide useful information for further functional investigations of the specific HSP20s (e.g., CtHSP20-9, CtHSP20-10, and CtHSP20-14) in African bermudagrass.