De novo genome assembly of two tomato ancestors, Solanum pimpinellifolium and Solanum lycopersicum var. cerasiforme, by long-read sequencing

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Solanum pimpinellifolium exhibits complex genetic resistance to Pseudomonas syringae pv. tomato

Pseudomonas syringae pv. tomato (Pst) is the causal agent of bacterial speck disease in tomatoes. The Pto/Prf gene cluster from Solanum pimpinellifolium was introgressed into several modern tomato cultivars and provided protection against Pst race 0 strains for many decades. However, virulent Pst race 1 strains that evade Pto-mediated immunity now predominate in tomato-growing regions worldwide. Here we report the identification of resistance to a Pst race 1 strain (Pst19) in the wild tomato accession S. pimpinellifolium LA1589 (hereafter LA1589), using our rapid high-throughput seedling screen. LA1589 supports less bacterial growth than cultivars, and does not exhibit a hypersensitive response to Pst19. We tested an existing set of 87 Inbred Backcross Lines (IBLs) derived from a cross between susceptible Solanum lycopersicum E-6203 and Solanum pimpinellifolium LA1589 for resistance to Pst19. Using single-marker analysis, we identified three genomic regions associated with resistance. Bacterial growth assays on IBLs confirmed that these regions contribute to resistance in planta. We also mapped candidate genes associated with resistance in a cross between the Solanum lycopersicum var. lycopersicum cultivar Heinz BG-1706 and S. pimpinellifolium LA1589. By comparing candidates from the two mapping approaches, we were able to identify 3 QTL and 5 candidate genes in LA1589 for a role in resistance to Pst19. This work will assist in molecular marker-assisted breeding to protect tomato from bacterial speck disease.

Chromosome-level genome assembly of Solanum pimpinellifolium

Solanum pimpinellifolium, the closest wild relative of the domesticated tomato, has high potential for use in breeding programs aimed at developing multi-pathogen resistance and quality improvement. We generated a chromosome-level genome assembly of S. pimpinellifolium LA1589, with a size of 833 Mb and a contig N50 of 31 Mb. We anchored 98.80% of the contigs into 12 pseudo-chromosomes, and identified 74.47% of the sequences as repetitive sequences. The genome evaluation revealed BUSCO and LAI score of 98.3% and 14.49, respectively, indicating high quality of this assembly. A total of 41,449 protein-coding genes were predicted in the genome, of which 89.17% were functionally annotated. This high-quality genome assembly serves as a valuable resource for accelerating the biological discovery and molecular breeding of this important horticultural crop.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

A chromosome-scale and haplotype-resolved genome assembly of carnation (Dianthus caryophyllus) based on high-fidelity sequencing

Dianthus caryophyllus is an economic species often considered excellent cut flowers and is suitable for bouquets and gardens. Here, we assembled the haplotype-resolved genome of D. caryophyllus 'Aili' at the chromosome level for the first time. The total lengths of the two assembled haplotypes of carnation were 584.88 Mb for haplotype genome 1 (hap1) and 578.78 Mb for haplotype genome 2 (hap2), respectively. We predicted a total of 44,098 and 42,425 protein-coding genes, respectively. The remarkable structure variation was identified between two haplotypes. Moreover, we identified 403.80 Mb of transposable elements (TEs) in hap1, which accounted for 69.34% of the genome. In contrast, hap2 had 402.70 Mb of TEs, representing 69.61% of the genome. Long terminal repeats were the predominant transposable elements. Phylogenetic analysis showed that the species differentiation time between carnation and gypsophila was estimated to be ~54.43 MYA. The unique gene families of carnation genomes were identified in 'Aili' and previously published 'Francesco' and 'Scarlet Queen'. The assembled and annotated haplotype-resolved D. caryophyllus genome not only promises to facilitate molecular biology studies but also contributes to genome-level evolutionary studies.

Identification and Expression Analysis of the Alfin-like Gene Family in Tomato and the Role of SlAL3 in Salt and Drought Stresses

Alfin-like (AL) transcription factors are a family of plant-specific genes with a PHD-finger-like structural domain at the C-terminus and a DUF3594 structural domain at the N-terminus that play important roles in plant development and stress response. In the present study, genome-wide identification and analysis were performed of the AL protein family in cultivated tomato (Solanum lycopersicum) and three wild relatives (S. pennellii, S. pimpinellifolium, and S. lycopersicoides) to evaluate their response to different abiotic stresses. A total of 39 ALs were identified and classified into four groups and based on phylogenetic tree and evolutionary analysis were shown to have formed prior to the differentiation of monocotyledons and dicots. Moreover, cis-acting element analysis revealed that various phytohormone response and abiotic stress response elements were highly existed in tomato. In addition, further analysis of the SlAL3 gene revealed that its expression was induced by drought and salt stresses and localized to the nucleus. In conclusion, our findings concerning AL genes provide useful information for further studies on their functions and regulatory mechanisms and provide theoretical references for studying AL gene response to abiotic stresses in plants.

A large-scale assessment of the quality of plant genome assemblies using the LTR assembly index

Recent advances in genome sequencing have led to an increase in the number of sequenced genomes. However, the presence of repetitive sequences complicates the assembly of plant genomes. The LTR assembly index (LAI) has recently been widely used to assess the quality of genome assembly, as a higher LAI is associated with a higher quality of assembly. Here, we assessed the quality of assembled genomes of 1664 plant and algal genomes using LAI and reported the results as data repository called PlantLAI (https://bioinformatics.um6p.ma/PlantLAI). A number of 55 117 586 pseudomolecules/scaffolds with a total length of 988.11 gigabase-pairs were examined using the LAI workflow. A total of 46 583 551 accurate LTR-RTs were discovered, including 2 263 188 Copia, 2 933 052 Gypsy, and 1 387 311 unknown superfamilies. Consequently, only 1136 plant genomes are suitable for LAI calculation, with values ranging from 0 to 31.59. Based on the quality classification system, 476 diploid genomes were classified as draft, 472 as reference, and 135 as gold genomes. We also provide a free webtool to calculate the LAI of newly assembled genomes and the ability to save the result in the repository. The data repository is designed to fill in the gaps in the reported LAI of existing genomes, while the webtool is designed to help researchers calculate the LAI of their newly sequenced genomes.

Chromosome-level genome assembly and annotation of the prickly nightshade Solanum rostratum Dunal

The prickly nightshade Solanum rostratum, an annual malignant weed, is native to North America and has globally invaded 34 countries, causing serious threats to ecosystems, agriculture, animal husbandry, and human health. In this study, we constructed a chromosome-level genome assembly and annotation of S. rostratum. The contig-level genome was initially assembled in 898.42 Mb with a contig N50 of 62.00 Mb from PacBio high-fidelity reads. With Hi-C sequencing data scaffolding, 96.80% of the initially assembled sequences were anchored and orientated onto 12 pseudo-chromosomes, generating a genome of 869.69 Mb with a contig N50 of 72.15 Mb. We identified 649.92 Mb (72.26%) of repetitive sequences and 3,588 non-coding RNAs in the genome. A total of 29,694 protein-coding genes were predicted, with 28,154 (94.81%) functionally annotated genes. We found 99.5% and 91.3% complete embryophyta_odb10 genes in the pseudo-chromosomes genome and predicted gene datasets by BUSCO assessment. The present genomic resource provides essential information for subsequent research on the mechanisms of environmental adaptation of S. rostratum and host shift in Colorado potato beetles.

Super-pangenome analyses highlight genomic diversity and structural variation across wild and cultivated tomato species

Effective utilization of wild relatives is key to overcoming challenges in genetic improvement of cultivated tomato, which has a narrow genetic basis; however, current efforts to decipher high-quality genomes for tomato wild species are insufficient. Here, we report chromosome-scale tomato genomes from nine wild species and two cultivated accessions, representative of Solanum section Lycopersicon, the tomato clade. Together with two previously released genomes, we elucidate the phylogeny of Lycopersicon and construct a section-wide gene repertoire. We reveal the landscape of structural variants and provide entry to the genomic diversity among tomato wild relatives, enabling the discovery of a wild tomato gene with the potential to increase yields of modern cultivated tomatoes. Construction of a graph-based genome enables structural-variant-based genome-wide association studies, identifying numerous signals associated with tomato flavor-related traits and fruit metabolites. The tomato super-pangenome resources will expedite biological studies and breeding of this globally important crop.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Transcriptome analysis reveals the regulation of cyclic nucleotide-gated ion channels in response to exogenous abscisic acid and calcium treatment under drought stress in tomato

Background: Drought stress can limit the growth and development of tomato seedlings and cause considerable loss of tomato yield. Exogenous abscisic acid (ABA) and calcium (Ca²⁺) can effectively alleviate the damage of drought stress to plants in part because Ca²⁺ acts as a second messenger in the drought resistance pathway. Although cyclic nucleotide-gated ion channels (CNGCs) are common non-specific Ca²⁺ osmotic channels in cell membranes, a thorough understanding of the transcriptome characteristics of tomato treated with exogenous ABA and Ca²⁺ under drought stress is necessary to characterize the molecular mechanism of CNGC involved in tomato drought resistance. Results: There were 12,896 differentially expressed genes in tomato under drought stress, as well as 11,406 and 12,502 differentially expressed genes after exogenous ABA and Ca²⁺ application, respectively. According to functional annotations and reports, the 19 SlCNGC genes related to Ca²⁺ transport were initially screened, with 11 SlCNGC genes that were upregulated under drought stress and downregulated after exogenous ABA application. After exogenous Ca²⁺ application, the data showed that two of these genes were upregulated, while nine genes were downregulated. Based on these expression patterns, we predicted the role of SlCNGC genes in the drought resistance pathway and their regulation by exogenous ABA and Ca²⁺ in tomato. Conclusion: The results of this study provide foundational data for further study of the function of SlCNGC genes and a more comprehensive understanding of drought resistance mechanisms in tomato.

Next generation sequencing technologies to explore the diversity of germplasm resources: achievements and trends in tomato

Tomato is one of the major vegetable crops grown worldwide and a model species for genetic and biological research. Progress in genomic technologies made possible the development of forefront methods for high-scale sequencing, providing comprehensive insight into the genetic architecture of germplasm resources. This review revisits next-generation sequencing strategies and applications to investigate the diversity of tomato, describing the common platforms used for SNP genotyping of large collections, de novo sequencing, and whole genome resequencing. Significant findings in evolutionary history are outlined, thus discussing how genomics has provided new hints about the processes behind domestication. Finally, achievement and perspectives on pan-genome construction and graphical pan-genome development toward precise mining of the natural variation to be exploited for breeding purposes are presented.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43-49 and ∼40-46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

Genome-wide analysis of the WRKY gene family unveil evolutionary history and expression characteristics in tomato and its wild relatives

WRKY transcription factors (WRKYs) are one of the largest plant gene families in plants involved in various biotic and abiotic stress responses. Based on the conservation of WRKY proteins, we identified a total of 642 WRKYs in Amborella trichopoda (33), Vitis vinifera (64), Arabidopsis thaliana (48), Solanum lycopersicoides (88), S. pennellii (77), S. pimpinellifolium (80), S. lycopersicum var. cerasiforme (85), S. lycopersicum cv. Heinz1706 (85), and S. lycopersicum cv. M82 (82) genomes. Phylogenetic analysis clustered WRKYs from nine genomes above into two clusters (Cluster1 and Cluster2). Evolutionary analysis revealed that most of the WRKYs in tomato and its wild relatives were expanded after the whole genome triplication (WGT) event of Solanum ancestor. Effects of tandem duplication (TD) event for WRKYs revealed that several WRKYs have experienced TD event and drove the expansion of the WRKY gene family in tomato and its wild relatives. Comparative analysis of WRKYs derived from WGT and TD events indicated that the WGT event performed a stronger influence on the expansion of the WRKY gene family than the effects of the TD event. Transcriptome profiling of WRKYs in S. lycopersicum cv. Heinz1706 under the biotic stress condition relative to the control condition uncovered a number of up-regulated WRKYs in response to biotic stress. The diversified expression pattern among paralogs derived from TD and WGT implied the impact of gene duplication events on gene functional divergence and diversity in tomato. We hope that this project will supply novel knowledge for studying the evolutionary history and functional characteristics of WRKYs involved in biotic stress in tomato.

Genome-wide identification, phylogenetic analysis, and expression profiles of trihelix transcription factor family genes in quinoa (Chenopodium quinoa Willd.) under abiotic stress conditions

Background

The trihelix family of transcription factors plays essential roles in the growth, development, and abiotic stress response of plants. Although several studies have been performed on the trihelix gene family in several dicots and monocots, this gene family is yet to be studied in Chenopodium quinoa (quinoa).

Results

In this study, 47 C. quinoa trihelix (CqTH) genes were in the quinoa genome. Phylogenetic analysis of the CqTH and trihelix genes from Arabidopsis thaliana and Beta vulgaris revealed that the genes were clustered into five subfamilies: SIP1, GTγ, GT1, GT2, and SH4. Additionally, synteny analysis revealed that the CqTH genes were located on 17 chromosomes, with the exception of chromosomes 8 and 11, and 23 pairs of segmental duplication genes were detected. Furthermore, expression patterns of 10 CqTH genes in different plant tissues and at different developmental stages under abiotic stress and phytohormone treatment were examined. Among the 10 genes, CqTH02, CqTH25, CqTH18, CqTH19, CqTH25, CqTH31, and CqTH36, were highly expressed in unripe achenes 21 d after flowering and in mature achenes compared with other plant tissues. Notably, the 10 CqTH genes were upregulated in UV-treated leaves, whereas CqTH36 was consistently upregulated in the leaves under all abiotic stress conditions.

Conclusions

The findings of this study suggest that gene duplication could be a major driver of trihelix gene evolution in quinoa. These findings could serve as a basis for future studies on the roles of CqTH transcription factors and present potential genetic markers for breeding stress-resistant and high-yielding quinoa varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08726-y.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene β-cyclase whose function we demonstrate.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

The genome of hibiscus hamabo reveals its adaptation to saline and waterlogged habitat

Hibiscus hamabo is a semi-mangrove species with strong tolerance to salt and waterlogging stress. However, the molecular basis and mechanisms that underlie this strong adaptability to harsh environments remain poorly understood. Here, we assembled a high-quality, chromosome-level genome of this semi-mangrove plant and analyzed its transcriptome under different stress treatments to reveal regulatory responses and mechanisms. Our analyses suggested that H. hamabo has undergone two recent successive polyploidy events, a whole-genome duplication followed by a whole-genome triplication, resulting in an unusually large gene number (107 309 genes). Comparison of the H. hamabo genome with that of its close relative Hibiscus cannabinus, which has not experienced a recent WGT, indicated that genes associated with high stress resistance have been preferentially preserved in the H. hamabo genome, suggesting an underlying association between polyploidy and stronger stress resistance. Transcriptomic data indicated that genes in the roots and leaves responded differently to stress. In roots, genes that regulate ion channels involved in biosynthetic and metabolic processes responded quickly to adjust the ion concentration and provide metabolic products to protect root cells, whereas no such rapid response was observed from genes in leaves. Using co-expression networks, potential stress resistance genes were identified for use in future functional investigations. The genome sequence, along with several transcriptome datasets, provide insights into genome evolution and the mechanism of salt and waterlogging tolerance in H. hamabo, suggesting the importance of polyploidization for environmental adaptation.

Organelle genome assembly uncovers the dynamic genome reorganization and cytoplasmic male sterility associated genes in tomato

To identify cytoplasmic male sterility (CMS)-associated genes in tomato, we determined the genome sequences of mitochondria and chloroplasts in three CMS tomato lines derived from independent asymmetric cell fusions, their nuclear and cytoplasmic donors, and male fertile weedy cultivated tomato and wild relatives. The structures of the CMS mitochondrial genomes were highly divergent from those of the nuclear and cytoplasmic donors, and genes of the donors were mixed up in these genomes. On the other hand, the structures of CMS chloroplast genomes were moderately conserved across the donors, but CMS chloroplast genes were unexpectedly likely derived from the nuclear donors. Comparative analysis of the structures and contents of organelle genes and transcriptome analysis identified three genes that were uniquely present in the CMS lines, but not in the donor or fertile lines. RNA-sequencing analysis indicated that these three genes transcriptionally expressed in anther, and identified different RNA editing levels in one gene, orf265, that was partially similar to ATP synthase subunit 8, between fertile and sterile lines. The orf265 was a highly potential candidate for CMS-associated gene. This study suggests that organelle reorganization mechanisms after cell fusion events differ between mitochondria and chloroplasts, and provides insight into the development of new F1 hybrid breeding programs employing the CMS system in tomato.

Genome-wide identification and expression profile analysis of trihelix transcription factor family genes in response to abiotic stress in sorghum [Sorghum bicolor (L.) Moench]

Background

Transcription factors, including trihelix transcription factors, play vital roles in various growth and developmental processes and in abiotic stress responses in plants. The trihelix gene has been systematically studied in some dicots and monocots, including Arabidopsis, tomato, chrysanthemum, soybean, wheat, corn, rice, and buckwheat. However, there are no related studies on sorghum.

Results

In this study, a total of 40 sorghum trihelix (SbTH) genes were identified based on the sorghum genome, among which 34 were located in the nucleus, 5 in the chloroplast, 1 (SbTH38) in the cytoplasm, and 1 (SbTH23) in the extracellular membrane. Phylogenetic analysis of the SbTH genes and Arabidopsis and rice trihelix genes indicated that the genes were clustered into seven subfamilies: SIP1, GTγ, GT1, GT2, SH4, GTSb8, and orphan genes. The SbTH genes were located in nine chromosomes and none on chromosome 10. One pair of tandem duplication gene and seven pairs of segmental duplication genes were identified in the SbTH gene family. By qPCR, the expression of 14 SbTH members in different plant tissues and in plants exposed to six abiotic stresses at the seedling stage were quantified. Except for the leaves in which the genes were upregulated after only 2 h exposure to high temperature, the 12 SbTH genes were significantly upregulated in the stems of sorghum seedlings after 24 h under the other abiotic stress conditions. Among the selected genes, SbTH10/37/39 were significantly upregulated, whereas SbTH32 was significantly downregulated under different stress conditions.

Conclusions

In this study, we identified 40 trihelix genes in sorghum and found that gene duplication was the main force driving trihelix gene evolution in sorghum. The findings of our study serve as a basis for further investigation of the functions of SbTH genes and providing candidate genes for stress-resistant sorghum breeding programmes and increasing sorghum yield.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08000-7.