An RNA Sequencing Transcriptome Analysis and Development of EST-SSR Markers in Chinese Hawthorn through Illumina Sequencing

Transcriptome analysis and development of EST-SSR markers in the mushroom Auricularia heimuer

Auricularia heimuer, the third most frequently cultivated edible mushroom species worldwide, has high medicinal value. However, a shortage of molecular marker hinders the efficiency and accuracy of genetic breeding efforts for A. heimuer. High-throughput transcriptome sequencing data are essential for gene discovery and molecular markers development. This study aimed to clarify the distribution of SSR loci across the A. heimuer transcriptome and to develop highly informative EST-SSR markers. These tools can be used for phylogenetic analysis, functional gene mining, and molecular marker-assisted breeding of A. heimuer. This study used Illumina high-throughput sequencing technology to obtain A. heimuer transcriptome data. The results revealed 37,538 unigenes in the A. heimuer transcriptome. Of these unigenes, 24,777 (66.01%) were annotated via comparison with the COG, Pfam, and NR databases. Overall, 2510 SSRs were identified from the unigenes, including 6 types of SSRs. The most abundant type of repeats were trinucleotides (1425, 56.77%), followed by mononucleotides (391, 15.58%) and dinucleotides (456, 18.17%). Primer pairs for 102 SSR loci were randomly designed for validity confirmation and polymorphism identification; this process yielded 53 polymorphic EST-SSR markers. Finally, 13 pairs of highly polymorphic EST-SSR primers were used to analyze the genetic diversity and population structure of 52 wild A. heimuer germplasms, revealing that the 52 germplasms could be divided into three categories. These results indicated that SSR loci were abundant in types, numbers, and frequencies, providing a potential basis for germplasm resource identification, genetic diversity analysis, and molecular marker-assisted breeding of A. heimuer.

Crataegus pinnatifida: A botanical, ethnopharmacological, phytochemical, and pharmacological overview

ETHNOPHARMACOLOGICAL RELEVANCE

Crataegus pinnatifida belongs to the Rosaceae family and extensively distribute in North China, Europe, and North America. Its usage was first described in "Xinxiu Ben Cao." The dried fruits of Crataegus pinnatifida Bunge or Crataegus pinnatifida var. major N. E. Br., also known as "Shanzha," is a famous medicine and food homology herb with a long history of medicinal usage in China. C. pinnatifida has the functions for digestive promotion, cardiovascular protection, and lipid reduction. It was traditionally used to treat indigestion, cardiodynia, thoracalgia, hernia, postpartum blood stagnation, and hemafecia. In recent years, C. pinnatifida has attracted worldwide attention as an important medicinal and economical crop due to its multiple and excellent health-promoting effects on cardiovascular, nervous, digestive, endocrine systems, and morbigenous microorganisms of the human body due to its medicinal and nutritional values.

AIM OF THE REVIEW

The current review aims to provide a comprehensive analysis of the geographical distribution, traditional usage, phytochemical components, pharmacological actions, clinical settings, and toxicities of C. pinnatifida. Moreover, the connection between the claimed biological activities and the traditional usage, along with the future perspectives for ongoing research on this plant, were also critically summarized.

MATERIALS AND METHODS

We collected the published literature on C. pinnatifida using a variety of scientific databases, including Web of Science, ScienceDirect, PubMed, Wiley, Springer, Taylor & Francis, ACS Publications, Google Scholar, Baidu Scholar, CNKI, The Plant List Database, and other literature sources (Ph.D. and MSc dissertations) from 2012 to 2022.

RESULTS

In the last decade, over 250 phytochemical compounds containing lignans, phenylpropanoids, flavonoids, triterpenoids, and their glycosides, as well as other compounds, have been isolated and characterized from different parts, including the fruit, leaves, and seeds of C. pinnatifida. Among these compounds, flavonoids and triterpenoids were major bioactive components of C. pinnatifida. They exhibited a broad spectrum of pharmacological actions with low toxicity in vitro and in vivo, such as cardiovascular protection, neuroprotection, anti-inflammatory, antioxidant, antibacterial, antiviral, anti-diabetes, anti-cancer, anti-mutagenic, anti-osteoporosis, anti-aging, anti-obesity, and hepatoprotection and other actions.

CONCLUSION

A long history of traditional uses and abundant pharmacochemical and pharmacological investigations have demonstrated that C. pinnatifida is an important medicine and food homology herb, which displays outstanding therapeutic potential, especially in the digestive system and cardiovascular disease. Nevertheless, the current studies on the active ingredients or crude extracts of C. pinnatifida and the possible mechanism of action are unclear. More evidence-based scientific studies are required to verify the traditional uses of C. pinnatifida. Furthermore, more efforts must be paid to selecting index components for quality control research and toxicity and safety studies of C. pinnatifida.

EST-SSR Markers' Development Based on RNA-Sequencing and Their Application in Population Genetic Structure and Diversity Analysis of Eleusine indica in China

Goosegrass (Eleusine indica) is one of the worst agricultural weeds in China. Molecular markers were developed for genetic diversity and population structure analyses. In this study, we identified 8391 expressed sequence tag-simple sequence repeat (EST-SSR) markers from the de novo assembled unigenes of E. indica. Mononucleotides were the most abundant type of repeats (3591, 42.79%), followed by trinucleotides (3162, 37.68%). The most dominant mononucleotide and trinucleotide repeat motifs were A/T (3406, 40.59%) and AAT/ATT (103, 1.5%), respectively. Fourteen pairs of EST-SSR primers were verified and used to analyze the genetic diversity and population structure of 59 goosegrass populations. A total of 49 alleles were amplified, with the number of alleles (N_a) ranging from two to eleven per locus, and the effective number of alleles (N_e) ranged from 1.07 to 4.53. The average polymorphic information content (PIC) was 0.36. Genetic structure analysis (K = 2) and principal coordinate analysis divided 59 E. indica populations into two groups in a manner similar to the unweighted pair-group method (Dice genetic similarity coefficient = 0.700). This study developed a set of EST-SSR markers in E. indica and successfully analyzed the diversity and population genetic structures of 59 E. indica populations in China.

New Insight into the Phylogeny and Taxonomy of Cultivated and Related Species of Crataegus in China, Based on Complete Chloroplast Genome Sequencing

Hawthorns (Crataegus L.) are one of the most important processing and table fruits in China, due to their medicinal properties and health benefits. However, the interspecific relationships and evolution history of cultivated Crataegus in China remain unclear. Our previously published data showed C. bretschneideri may be derived from the hybridization of C. pinnatifida with C. maximowiczii, and that introgression occurs between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. In the present study, chloroplast sequences were used to further elucidate the phylogenetic relationships of cultivated Crataegus native to China. The chloroplast genomes of three cultivated species and one related species of Crataegus were sequenced for comparative and phylogenetic analyses. The four chloroplast genomes of Crataegus exhibited typical quadripartite structures and ranged from 159,607 bp (C. bretschneideri) to 159,875 bp (C. maximowiczii) in length. The plastomes of the four species contained 113 genes consisting of 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Six hypervariable regions (ndhC-trnV(UAC)-trnM(CAU), ndhA, atpH-atpI, ndhF, trnR(UCU)-atpA, and ndhF-rpl32), 196 repeats, and a total of 386 simple sequence repeats were detected as potential variability makers for species identification and population genetic studies. In the phylogenomic analyses, we also compared the entire chloroplast genomes of three published Crataegus species: C. hupehensis (MW201730.1), C. pinnatifida (MN102356.1), and C. marshallii (MK920293.1). Our phylogenetic analyses grouped the seven Crataegus taxa into two main clusters. One cluster included C. bretschneideri, C. maximowiczii, and C. marshallii, whereas the other included C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. Taken together, our findings indicate that C. maximowiczii is the maternal origin of C. bretschneideri. This work provides further evidence of introgression between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major, and suggests that C. pinnatifida var. major might have been artificially selected and domesticated from hybrid populations, rather than evolved from C. pinnatifida.

Transcriptome sequencing and microsatellite marker discovery in Ailanthus altissima (Mill.) Swingle (Simaroubaceae)

Ailanthus altissima Swingle, is a tree species native to East Asia and has a great potential in decorative, bioenergy and industrial applications in many countries. To date, despite its commercial importance, the genomic and genetic resources available for this species are still insufficient. In this study, we characterized the transcriptome of A. altissima and developed thirteen EST-SSRs (expressed sequence tag-simple sequence repeats) based on Illumina paired-end RNA sequencing (RNA-seq). Besides, we developed ten polymorphic chloroplast microsatellite (cpSSR) markers using the available chloroplast genome of A. altissima. The transcriptome data produced 87,797 unigenes, of which 64,891 (73.91%) unigenes were successfully annotated in at least one protein database. For cpSSR markers the number of detected alleles (N) per marker varied from three at cpSSR12 to twelve at cpSSR8, the unbiased haploid diversity indices (uh) varied from 0.111 to 0.485, and haploid diversity indices (h) ranged from 0.101 to 0.444 with an average unbiased haploid diversity index (uh) of 0.274. Overall, a total of 65 different cpSSR alleles were identified at the ten loci among 165 individuals of A. altissima. The allele number per locus for EST-SSRs varied from 2.143 to 9.357, and the values of observed and expected heterozygosity ranged from 0.312 to 1.000 and 0.505 to 0.826, respectively. The molecular markers developed in this study will facilitate future genetic diversity, population structure, long distance-gene transfer and pollen-based gene flow analyses of A. altissima populations from its known distribution ranges in China focusing on planted and natural forest stands.

The complete chloroplast genome of Crataegus hupehensis Sarg. (Rosaceae), a medicinal and edible plant in China

Crataegus hupehensis Sarg. is well-known for its medicinal and nutritive value. In this study, the complete chloroplast genome sequence of C. hupehensis was determined by using Illumina high-throughput sequencing approach. The complete chloroplast genome is 159,766 bp with 36.6% GC content. It contained a pair of inverted repeat regions of 26,385 bp, a large single-copy region of 87,852 bp, and a small single-copy region of 19,144 bp. It contained 112 distinct genes, including 78 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. Phylogenetic analysis based on chloroplast genomes indicated that C. hupehensisis was closely related to C. kansuensis and C. marshallii in the subfamily Maloideae. This complete chloroplast genome will provide valuable insight into evolution, molecular breeding, and phylogenetic analysis of Crataegus species.

De novo transcriptomic analysis and identification of EST-SSR markers in Stephanandra incisa

Stephanandra incisa is a wild-type shrub with beautiful leaves and white flowers and is commonly used as a garden decoration accessory. However, the limited availability of genomic data of S. incisa has restricted its breeding process. Here, we identified EST-SSR markers using de novo transcriptome sequencing. In this study, a transcriptome database containing 35,251 unigenes, having an average length of 985 bp, was obtained from S. incisa. From these unigene sequences, we identified 5,555 EST-SSRs, with a distribution density of one SSR per 1.60 kb. Dinucleotides (52.96%) were the most detected SSRs, followed by trinucleotides (34.64%). From the EST-SSR loci, we randomly selected 100 sites for designing primer and used the DNA of 60 samples to verify the polymorphism. The average value of the effective number of alleles (Ne), Shannon's information index (I), and expective heterozygosity (He) was 1.969, 0.728, and 0.434, respectively. The polymorphism information content (PIC) value was in the range of 0.108 to 0.669, averaging 0.406, which represented a middle polymorphism level. Cluster analysis of S. incisa were also performed based on the obtained EST-SSR data in our work. As shown by structure analysis, 60 individuals could be classified into two groups. Thus, the identification of these novel EST-SSR markers provided valuable sequence information for analyzing the population structure, genetic diversity, and genetic resource assessment of S. incisa and other related species.

De novo transcriptome assembly and analysis of Phragmites karka, an invasive halophyte, to study the mechanism of salinity stress tolerance

With the rapidly deteriorating environmental conditions, the development of stress tolerant plants has become a priority for sustaining agricultural productivity. Therefore, studying the process of stress tolerance in naturally tolerant species hold significant promise. Phragmites karka is an invasive plant species found abundantly in tropical and sub tropical regions, fresh water regions and brackish marshy areas, such as river banks and lake shores. The plant possesses the ability to adapt and survive under conditions of high salinity. We subjected P. karka seedlings to salt stress and carried out whole transcriptome profiling of leaf and root tissues. Assessing the global transcriptome changes under salt stress resulted in the identification of several genes that are differentially regulated under stress conditions in root and leaf tissue. A total of 161,403 unigenes were assembled and used as a reference for digital gene expression analysis. A number of key metabolic pathways were found to be over-represented. Digital gene expression analysis was validated using qRT-PCR. In addition, a number of different transcription factor families including WRKY, MYB, CCCH, NAC etc. were differentially expressed under salinity stress. Our data will facilitate further characterisation of genes involved in salinity stress tolerance in P. karka. The DEGs from our results are potential candidates for understanding and engineering abiotic stress tolerance in plants.

Constructing a Core Collection of the Medicinal Plant Angelica biserrata Using Genetic and Metabolic Data

Angelica biserrata is an important medicinal plant in Chinese traditional medicine. Its roots, which are known as Duhuo in Chinese, are broadly applied to treat inflammation, arthritis, and headache. With increasing market demand, the wild resources of A. biserrata have been overexploited, and conservation, assessment of genetic resources and breeding for this species is needed. Here, we sequenced the transcriptome of A. biserrata and developed simple sequence repeat (SSR) markers from it to construct a core collection based on 208 samples collected from Changyang-related regions. A total of 132 alleles were obtained for 17 SSR loci used with the polymorphic information content (PIC) ranging from 0.44 to 0.83. Abundant genetic diversity was inferred by Shannon's information index (1.51), observed (0.57) and expected heterozygosity (0.72). The clustering analysis resulted into two sample groups and analysis of molecular variance (AMOVA) showed only 6% genetic variation existed among populations. A further metabolic analysis of these samples revealed the main coumarin contents, such as osthole and columbianadin. According to the genetic and metabolic data, we adopted the least distance stepwise sampling strategy to construct seven preliminary core collections, of which the 20CC collection, which possessed 42 A. biserrata individuals accounting for 90.20% of the genetic diversity of the original germplasm, represented the best core collection. This study will contribute to the conservation and management of A. biserrata wild germplasm resources and provide a material basis for future selection and breeding of this medicinal plant.

Plant virus interaction mechanism and associated pathways in mosaic disease of small cardamom (Elettaria cardamomum Maton) by RNA-Seq approach

Small cardamom (Elettaria cardamomum), grown in limited coastal tropical countries is one of the costliest and widely exported agri-produce having global turnover of >10 billion USD. Mosaic/marble disease is one of the major impediments that requires understanding of disease at molecular level. Neither whole genome sequence nor any genomic resources are available, thus RNA seq approach can be a rapid and economical alternative. De novo transcriptome assembly was done with Illumina Hiseq data. A total of 5317 DEGs, 2267 TFs, 114 pathways and 175,952 genic region putative markers were obtained. Gene regulatory network analysis deciphered molecular events involved in marble disease. This is the first transcriptomic report revealing disease mechanism mediated by perturbation in auxin homeostasis and ethylene signalling leading to senescence. The web-genomic resource (SCMVTDb) catalogues putative molecular markers, candidate genes and transcript information. SCMVTDb can be used in germplasm improvement against mosaic disease in endeavour of small cardamom productivity. Availability of genomic resource, SCMVTDb: http://webtom.cabgrid.res.in/scmvtdb/.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Development of SSR markers via de novo transcriptome assembly in Akebia trifoliata (Thunb.) Koidz

Owing to its high nutritive, economic, and medicinal values, Akebia trifoliata has received increased attention, making worthy of being used as a new fruit crop for further domestication and commercialization in China. However, molecular research of A. trifoliata has lagged as investigations of its genomic resources and molecular markers are rare. In this study, a cDNA library of A. trifoliata leaves was sequenced using the Illumina NovaSeq. 6000 sequencing system. In total, 101 417 transcripts, 63 757 unigenes, and 9494 simple sequence repeats were assembled and identified from the transcriptome datasets. The majority of the SSRs were di- and trinucleotide repeats. Length and number of SSR motifs ranged from 15 to 66, and 5 to 48 bp, respectively. Of which, the A/T mononucleotide motif and AG/TC and CT/GA dinucleotide motifs were the most abundant. Furthermore, 100 SSR primers were randomly selected to validate amplification and polymorphism, and 88 A. trifoliata accessions were definitively distinguished by 49 primers. With the Qinling mountains and Huaihe River line as the boundaries, the northern and southern accessions were clustered into different groups, but no clear geographical patterns (city or origin) were observed in the southern accessions. These newly identified molecular markers may provide a foundation for the genetic identification and diversity analysis and marker-assisted selection breeding in species of Akebia.

Transcriptome Profiling Provides Insight into the Genes in Carotenoid Biosynthesis during the Mesocarp and Seed Developmental Stages of Avocado (Persea americana)

Avocado (Persea americana Mill.) is an economically important crop because of its high nutritional value. However, the absence of a sequenced avocado reference genome has hindered investigations of secondary metabolism. For next-generation high-throughput transcriptome sequencing, we obtained 365,615,152 and 348,623,402 clean reads as well as 109.13 and 104.10 Gb of sequencing data for avocado mesocarp and seed, respectively, during five developmental stages. High-quality reads were assembled into 100,837 unigenes with an average length of 847.40 bp (N50 = 1725 bp). Additionally, 16,903 differentially expressed genes (DEGs) were detected, 17 of which were related to carotenoid biosynthesis. The expression levels of most of these 17 DEGs were higher in the mesocarp than in the seed during five developmental stages. In this study, the avocado mesocarp and seed transcriptome were also sequenced using single-molecule long-read sequencing to acquired 25.79 and 17.67 Gb clean data, respectively. We identified 233,014 and 238,219 consensus isoforms in avocado mesocarp and seed, respectively. Furthermore, 104 and 59 isoforms were found to correspond to the putative 11 carotenoid biosynthetic-related genes in the avocado mesocarp and seed, respectively. The isoform numbers of 10 out of the putative 11 genes involved in the carotenoid biosynthetic pathway were higher in the mesocarp than those in the seed. Besides, alpha- and beta-carotene contents in the avocado mesocarp and seed during five developmental stages were also measured, and they were higher in the mesocarp than in the seed, which validated the results of transcriptome profiling. Gene expression changes and the associated variations in gene dosage could influence carotenoid biosynthesis. These results will help to further elucidate carotenoid biosynthesis in avocado.

Transcriptome Sequencing of Different Avocado Ecotypes: de novo Transcriptome Assembly, Annotation, Identification and Validation of EST-SSR Markers

Avocado (Persea americana Mill.) could be considered as an important tropical and subtropical woody oil crop with high economic and nutritional value. Despite the importance of this species, genomic information is currently unavailable for avocado and closely related congeners. In this study, we generated more than 216 million clean reads from different avocado ecotypes using Illumina HiSeq high-throughput sequencing technology. The high-quality reads were assembled into 154,310 unigenes with an average length of 922 bp. A total of 55,558 simple sequence repeat (SSR) loci detected among the 43,270 SSR-containing unigene sequences were used to develop 74,580 expressed sequence tag (EST)-SSR markers. From these markers, a subset of 100 EST-SSR markers was randomly chosen to identify polymorphic EST-SSR markers in 28 avocado accessions. Sixteen EST-SSR markers with moderate to high polymorphism levels were detected, with polymorphism information contents ranging from 0.33 to 0.84 and averaging 0.63. These 16 polymorphic EST-SSRs could clearly and effectively distinguish the 28 avocado accessions. In summary, our study is the first presentation of transcriptome data of different avocado ecotypes and comprehensive study on the development and analysis of a set of EST-SSR markers in avocado. The application of next-generation sequencing techniques for SSR development is a potentially powerful tool for genetic studies.