Comprehensive analysis of SSRs and database construction using all complete gene-coding sequences in major horticultural and representative plants

TriticeaeSSRdb: a comprehensive database of simple sequence repeats in Triticeae

Microsatellites, known as simple sequence repeats (SSRs), are short tandem repeats of 1 to 6 nucleotide motifs found in all genomes, particularly eukaryotes. They are widely used as co-dominant markers in genetic analyses and molecular breeding. Triticeae, a tribe of grasses, includes major cereal crops such as bread wheat, barley, and rye, as well as abundant forage and lawn grasses, playing a crucial role in global food production and agriculture. To enhance genetic work and expedite the improvement of Triticeae crops, we have developed TriticeaeSSRdb, an integrated and user-friendly database. It contains 3,891,705 SSRs from 21 species and offers browsing options based on genomic regions, chromosomes, motif types, and repeat motif sequences. Advanced search functions allow personalized searches based on chromosome location and length of SSR. Users can also explore the genes associated with SSRs, design customized primer pairs for PCR validation, and utilize practical tools for whole-genome browsing, sequence alignment, and in silico SSR prediction from local sequences. We continually update TriticeaeSSRdb with additional species and practical utilities. We anticipate that this database will greatly facilitate trait genetic analyses and enhance molecular breeding strategies for Triticeae crops. Researchers can freely access the database at http://triticeaessrdb.com/.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Flowering genes identification, network analysis, and database construction for 837 plants

Flowering is one of the most important biological phenomena in the plant kingdom, which not only has important ecological significance, but also has substantial horticultural ornamental value. In this study, we undertook an exhaustive review of the advancements in our understanding of plant flowering genes. We delved into the identification and conducted comparative analyses of flowering genes across virtually all sequenced angiosperm plant genomes. Furthermore, we established an extensive angiosperm flowering atlas, encompassing a staggering 183 720 genes across eight pathways, along with 10 155 ABCDE mode genes, which play a pivotal role in plant flowering regulation. Through the examination of expression patterns, we unveiled the specificities of these flowering genes. An interaction network between flowering genes of the ABCDE model and their corresponding upstream genes offered a blueprint for comprehending their regulatory mechanisms. Moreover, we predicted the miRNA and target genes linked to the flowering processes of each species. To culminate our efforts, we have built a user-friendly web interface, named the Plant Flowering-time Gene Database (PFGD), accessible at http://pfgd.bio2db.com/. We firmly believe that this database will serve as a cornerstone in the global research community, facilitating the in-depth exploration of flowering genes in the plant kingdom. In summation, this pioneering endeavor represents the first comprehensive collection and comparative analysis of flowering genes in plants, offering valuable resources for the study of plant flowering genetics.

The genome of Stephania japonica provides insights into the biosynthesis of cepharanthine

Stephania japonica is an early-diverging eudicotyledon plant with high levels of cepharanthine, proven to be effective in curing coronavirus infections. Here, we report a high-quality S. japonica genome. The genome size is 688.52 Mb, and 97.37% sequences anchor to 11 chromosomes. The genome comprises 67.46% repetitive sequences and 21,036 genes. It is closely related to two Ranunculaceae species, which diverged from their common ancestor 55.90-71.02 million years ago (Mya) with a whole-genome duplication 85.59-96.75 Mya. We further reconstruct ancestral karyotype of Ranunculales. Several cepharanthine biosynthesis genes are identified and verified by western blot. Two genes (Sja03G0243 and Sja03G0241) exhibit catalytic activity as shown by liquid chromatography-mass spectrometry. Then, cepharanthine biosynthesis genes, transcription factors, and CYP450 family genes are used to construct a comprehensive network. Finally, we construct an early-diverging eudicotyledonous genome resources (EEGR) database. As the first genome of the Menispermaceae family to be released, this study provides rich resources for genomic studies.

Microsatellite Content in 397 Nuclear Exons and Their Flanking Regions in the Fern Family Ophioglossaceae

Microsatellites or SSRs are small tandem repeats that are 1-6 bp long. They are usually highly polymorphic and form important portions of genomes. They have been extensively analyzed in humans, animals and model plants; however, information from non-flowering plants is generally lacking. Here, we examined 29 samples of Ophioglossaceae ferns, mainly from the genera Botrychium and Sceptridium. We analyzed the SSR distribution, density and composition in almost 400 nuclear exons and their flanking regions. We detected 45 SSRs in exons and 1475 SSRs in the flanking regions. In the exons, only di-, tri- and tetranucleotides were found, and all of them were 12 bp long. The annotation of the exons containing SSRs showed that they were related to various processes, such as metabolism, catalysis, transportation or plant growth. The flanking regions contained SSRs from all categories, with the most numerous being dinucleotides, followed by tetranucleotides. More than one-third of all the SSRs in the flanking regions were 12 bp long. The SSR densities in the exons were very low, ranging from 0 to 0.07 SSRs/kb, while those in the flanking regions ranged from 0.24 to 0.81 SSRs/kb; and those in the combined dataset ranged from 0.2 to 0.81 SSRs/kb. The majority of the detected SSRs in the flanking regions were polymorphic and present at the same loci across two or more samples but differing in the number of repeats. The SSRs detected here may serve as a basis for further population genetic, phylogenetic or evolutionary genetic studies, as well as for further studies focusing on SSRs in the genomes and their roles in adaptation, evolution and diseases.

PHGD: An integrative and user-friendly database for plant hormone-related genes

Plant Hormone Gene Database (PHGD) database platform construction pipeline. First, we collected all reported hormone-related genes in the model plant Arabidopsis thaliana, and combined with the existing experimental background, mapped the hormone-gene interaction network to provide a blueprint. Next, we collected 469 high-quality plant genomes. Then, bioinformatics was used to identify hormone-related genes in these plants. Finally, these genetic data were programmed to be stored in a database and a platform website PHGD was built. PHGD was divided into eight modules, namely Home, Browse, Search, Resources, Download, Tools, Help, and Contact. We provided data resources and platform services to facilitate the study of plant hormones.

Construction of a high-density genetic map for yardlong bean and identification of ANT1 as a regulator of anthocyanin biosynthesis

Because its long, tender pods supply essential proteins, vitamins, and fibers to humans, yardlong bean (Vigna unguiculata ssp. sesquipedalis) is a commonly consumed vegetable, especially in Southeast Asia. To provide insights into the genetic bases of key agricultural traits in yardlong bean, we here created a high-density bin-map with 2084 bin markers using 514 227 SNPs from a recombinant-inbred line (RIL) population. Quantitative trait loci (QTL) mapping was carried out to identify loci associated with anthocyanin content (ANT), vitamin E content (VE), total soluble protein content (TSP), pod length (PL), hundred-seed weight (HSW), seed length and width (SL and SW, respectively), and seed coat color (SCC). In total, 20 related QTLs were isolated, explaining 7.58-56.03% of the phenotypic variation. Of these, five major QTLs (qANT5, qTSP11, qVE7, qPL3, and qSCC9) were detected in 2020, 2021, and the combined environment, explaining 11.96-56.03% of the phenotypic variation. VuANT1 was identified as a causal gene for the QTL qANT5, which regulated anthocyanin content; VuANT1 was highly expressed in immature purple pods but barely detectable in white pods. VuANT1 overexpression in tobacco leaves and yardlong bean hairy roots resulted in purple coloration as a result of anthocyanin accumulation. These findings suggested that VuANT1 was a key regulator of anthocyanin accumulation in yardlong bean. Our results lay a firm foundation for target agricultural trait improvement and clarification of the genetic mechanisms underlying agricultural traits in yardlong bean.

The high-quality sequencing of the Brassica rapa 'XiangQingCai' genome and exploration of genome evolution and genes related to volatile aroma

'Vanilla' (XQC, brassica variety chinensis) is an important vegetable crop in the Brassica family, named for its strong volatile fragrance. In this study, we report the high-quality chromosome-level genome sequence of XQC. The assembled genome length was determined as 466.11 Mb, with an N50 scaffold of 46.20 Mb. A total of 59.50% repetitive sequences were detected in the XQC genome, including 47 570 genes. Among all examined Brassicaceae species, XQC had the closest relationship with B. rapa QGC ('QingGengCai') and B. rapa Pakchoi. Two whole-genome duplication (WGD) events and one recent whole-genome triplication (WGT) event occurred in the XQC genome in addition to an ancient WGT event. The recent WGT was observed to occur during 21.59-24.40 Mya (after evolution rate corrections). Our findings indicate that XQC experienced gene losses and chromosome rearrangements during the genome evolution of XQC. The results of the integrated genomic and transcriptomic analyses revealed critical genes involved in the terpenoid biosynthesis pathway and terpene synthase (TPS) family genes. In summary, we determined a chromosome-level genome of B. rapa XQC and identified the key candidate genes involved in volatile fragrance synthesis. This work can act as a basis for the comparative and functional genomic analysis and molecular breeding of B. rapa in the future.

Genome-Wide Identification of the WD40 Gene Family in Tomato (Solanum lycopersicum L.)

WD40 proteins are a superfamily of regulatory proteins widely found in eukaryotes that play an important role in regulating plant growth and development. However, the systematic identification and characterization of WD40 proteins in tomato (Solanum lycopersicum L.) have not been reported. In the present study, we identified 207 WD40 genes in the tomatoes genome and analyzed their chromosomal location, gene structure and evolutionary relationships. A total of 207 tomato WD40 genes were classified by structural domain and phylogenetic tree analyses into five clusters and 12 subfamilies and were found to be unevenly distributed across the 12 tomato chromosomes. We identified six tandem duplication gene pairs and 24 segmental duplication pairs in the WD40 gene family, with segmental duplication being the major mode of expansion in tomatoes. Ka/Ks analysis revealed that paralogs and orthologs of WD40 family genes underwent mainly purifying selection during the evolutionary process. RNA-seq data from different tissues and developmental periods of tomato fruit development showed tissue-specific expression of WD40 genes. In addition, we constructed four coexpression networks according to the transcriptome and metabolome data for WD40 proteins involved in fruit development that may be related to total soluble solid formation. The results provide a comprehensive overview of the tomato WD40 gene family and will provide valuable information for the validation of the function of tomato WD40 genes in fruit development.

High-quality ice plant reference genome analysis provides insights into genome evolution and allows exploration of genes involved in the transition from C3 to CAM pathways

Ice plant (Mesembryanthemum crystallinum), a member of the Aizoaceae family, is a typical halophyte crop and a model plant for studying the mechanism of transition from C3 photosynthesis to crassulacean acid metabolism (CAM). Here, we report a high-quality chromosome-level ice plant genome sequence. This 98.05% genome sequence is anchored to nine chromosomes, with a total length of 377.97 Mb and an N50 scaffold of 40.45 Mb. Almost half of the genome (48.04%) is composed of repetitive sequences, and 24 234 genes have been annotated. Subsequent to the ancient whole-genome triplication (WGT) that occurred in eudicots, there has been no recent whole-genome duplication (WGD) or WGT in ice plants. However, we detected a novel WGT event that occurred in the same order in Simmondsia chinensis, which was previously overlooked. Our findings revealed that ice plants have undergone chromosome rearrangements and gene removal during evolution. Combined with transcriptome and comparative genomic data and expression verification, we identified several key genes involved in the CAM pathway and constructed a comprehensive network. As the first genome of the Aizoaceae family to be released, this report will provide a rich data resource for comparative and functional genomic studies of Aizoaceae, especially for studies on salt tolerance and C3-to-CAM transitions to improve crop yield and resistance.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

The Brassicaceae genome resource (TBGR): A comprehensive genome platform for Brassicaceae plants

The Brassicaceae is an important plant family. We built a user-friendly, web-based, comparative, and functional genomic database, The Brassicaceae Genome Resource (TBGR, http://www.tbgr.org.cn), based on 82 released genomes from 27 Brassicaceae species. The TBGR database contains a large number of important functional genes, including 4,096 glucosinolate genes, 6,625 auxin genes, 13,805 flowering genes, 36,632 resistance genes, 1,939 anthocyanin genes, and 1,231 m6A genes. A total of 1,174,049 specific guide sequences for clustered regularly interspaced short palindromic repeats and 5,856,479 transposable elements were detected in Brassicaceae. TBGR also provides information on synteny, duplication, and orthologs for 27 Brassicaceae species. The TBGR database contains 1,183,851 gene annotations obtained using the TrEMBL, Swiss-Prot, Nr, GO, and Pfam databases. The BLAST, Synteny, Primer Design, Seq_fetch, and JBrowse tools are provided to help users perform comparative genomic analyses. All the genome assemblies, gene models, annotations, and bioinformatics results can be easily downloaded from the TBGR database. We plan to improve and continuously update the database with newly assembled genomes and comparative genomic studies. We expect the TBGR database to become a key resource for the study of the Brassicaceae.

Comprehensive Analyses of Simple Sequence Repeat (SSR) in Bamboo Genomes and Development of SSR Markers with Peroxidase Genes

Simple sequence repeats (SSRs) are one of the most important molecular markers, which are widespread in plants. Bamboos are important forest resources worldwide. Here, the comprehensive identification and comparative analysis of SSRs were performed in three woody and two herbaceous bamboo species. Altogether 567,175 perfect SSRs and 71,141 compound SSRs were identified from 5737.8 Mb genome sequences of five bamboo species. Di-nucleotide SSRs were the most predominant type, with an average of ~50,152.2 per species. Most SSRs were located in intergenic regions, while those located in genic regions were relatively less. Moreover, the results of annotation distribution indicated that terms with P450, peroxidase and ATP-binding cassette transporter related to lignin biosynthesis might play important roles in woody and herbaceous bamboos under the mediation of SSRs. Furthermore, the peroxidase gene family consisted of a large number of genes containing SSRs was selected for the evolutionary relationship analysis and SSR markers development. Fifteen SSR markers derived from peroxidase family genes of Phyllostachys edulis were identified as polymorphic in 34 accessions belonging to seven genera in Bambusoideae. These results provided a comprehensive insight of SSR markers into bamboo genomes, which would facilitate bamboo research related to comparative genomics, evolution and marker-assisted selection.

Large-scale analyses of heat shock transcription factors and database construction based on whole-genome genes in horticultural and representative plants

Heat shock transcription factor (Hsf) plays a critical role in regulating heat resistance. Here, 2950 Hsf family genes were identified from 111 horticultural and representative plants. More Hsf genes were detected in higher plants than lower plants. Based on all Hsf genes, we constructed a phylogenetic tree, which indicated that Hsf genes of each branch evolved independently after species differentiation. Furthermore, we uncovered the evolutionary trajectories of Hsf genes by motif analysis. There were only 6 motifs (M1 to M6) in lower plants, and then 4 novel motifs (M7-M10) appeared in higher plants. However, the motifs of some Hsf genes were lost in higher plant, indicating that Hsf genes have undergone sequence variation during the evolution. The number of Hsf gene loss was more than duplication after whole-genome duplication in higher plants. The heat response network was constructed using 24 Hsf genes, 2421 downstream, and 222 upstream genes of Arabidopsis. Further enrichment analysis revealed that Hsf genes and other transcription factors interacted with each other to response heat resistance. The global expression maps were illustrated for Hsf genes under various abiotic, biotic stresses, and several developmental stages in Arabidopsis. The syntenic and phylogenetic analyses were conducted using Hsf genes of Arabidopsis and Pan-genome of 18 Brassica rapa accessions. We also performed the expression pattern analysis of Hsf and six Hsp family genes using expression values from different tissues and heat treatments in B. rapa. The interaction network between Hsf and Hsp gene families was constructed in B. rapa, and several core genes were detected in the network. Finally, we constructed a Hsf database (http://hsfdb.bio2db.com) for researchers to retrieve Hsf gene family information. Therefore, our study will provide rich resources for the evolution and functional study of Hsf genes.

TCMPG: an integrative database for traditional Chinese medicine plant genomes

Because of their great therapeutic and economic value, medicinal plants have attracted increasing scientific attention. With the rapid development of high-throughput sequencing technology, the genomes of many medicinal plants have been sequenced. Storing and analyzing the increasing volume of genomic data has become an urgent task. To solve this challenge, we have proposed the Traditional Chinese Medicine Plant Genome database (TCMPG, http://cbcb.cdutcm.edu.cn/TCMPG/), an integrative database for storing the scattered genomes of medicinal plants. TCMPG currently includes 160 medicinal plants, 195 corresponding genomes, and 255 herbal medicines. Detailed information on plant species, genomes, and herbal medicines is also integrated into TCMPG. Popular genomic analysis tools are embedded in TCMPG to facilitate the systematic analysis of medicinal plants. These include BLAST for identifying orthologs from different plants, SSR Finder for identifying simple sequence repeats, JBrowse for browsing genomes, Synteny Viewer for displaying syntenic blocks between two genomes, and HmmSearch for identifying protein domains. TCMPG will be continuously updated by integrating new data and tools for comparative and functional genomic analysis.

An empirical analysis of mtSSRs: could microsatellite distribution patterns explain the evolution of mitogenomes in plants?

Microsatellites (SSRs) are tandem repeat sequences in eukaryote genomes, including plant cytoplasmic genomes. The mitochondrial genome (mtDNA) has been shown to vary in size, number, and distribution of SSRs among different plant groups. Thus, SSRs contribute with genomic diversity in mtDNAs. However, the abundance, distribution, and evolutionary significance of SSRs in mtDNA from a wide range of algae and plants have not been explored. In this study, the mtDNAs of 204 plant and algal species were investigated related to the presence of SSRs. The number of SSRs was positively correlated with genome size. Its distribution is dependent on plant and algal groups analyzed, although the cluster analysis indicates the conservation of some common motifs in algal and terrestrial plants that reflect common ancestry of groups. Many SSRs in coding and non-coding regions can be useful for molecular markers. Moreover, mitochondrial SSRs are highly abundant, representing an important source for natural or induced genetic variation, i.e., for biotechnological approaches that can modulate mtDNA gene regulation. Thus, this comparative study increases the understanding of the plant and algal SSR evolution and brings perspectives for further studies.