Development and Transferability of EST-SSR Markers for Pinus koraiensis from Cold-Stressed Transcriptome through Illumina Sequencing

Widely Targeted Metabolomic Profiling Combined with Transcriptome Analysis Provides New Insights into Lipid Biosynthesis in Seed Kernels of Pinus koraiensis

Lipid-rich Pinus koraiensis seed kernels are highly regarded for their nutritional and health benefits. To ascertain the molecular mechanism of lipid synthesis, we conducted widely targeted metabolomic profiling together with a transcriptome analysis of the kernels in P. koraiensis cones at various developmental stages. The findings reveal that 148 different types of lipid metabolites, or 29.6% of total metabolites, are present in kernels. Among those metabolites, the concentrations of linoleic acid, palmitic acid, and α-linolenic acid were higher, and they steadily rose as the kernels developed. An additional 10 hub genes implicated in kernel lipid synthesis were discovered using weighted gene co-expression network analysis (WGCNA), gene interaction network analysis, oil body biosynthesis, and transcriptome analysis. This study used lipid metabolome and transcriptome analyses to investigate the mechanisms of key regulatory genes and lipid synthesis molecules during kernel development, which served as a solid foundation for future research on lipid metabolism and the creation of P. koraiensis kernel food.

De novo full-length transcriptome analysis of two ecotypes of Phragmites australis (swamp reed and dune reed) provides new insights into the transcriptomic complexity of dune reed and its long-term adaptation to desert environments

BACKGROUND

The extremely harsh environment of the desert is changing dramatically every moment, and the rapid adaptive stress response in the short term requires enormous energy expenditure to mobilize widespread regulatory networks, which is all the more detrimental to the survival of the desert plants themselves. The dune reed, which has adapted to desert environments with complex and variable ecological factors, is an ideal type of plant for studying the molecular mechanisms by which Gramineae plants respond to combinatorial stress of the desert in their natural state. But so far, the data on the genetic resources of reeds is still scarce, therefore most of their research has focused on ecological and physiological studies.

RESULTS

In this study, we obtained the first De novo non-redundant Full-Length Non-Chimeric (FLNC) transcriptome databases for swamp reeds (SR), dune reeds (DR) and the All of Phragmites australis (merged of iso-seq data from SR and DR), using PacBio Iso-Seq technology and combining tools such as Iso-Seq3 and Cogent. We then identified and described long non-coding RNAs (LncRNA), transcription factor (TF) and alternative splicing (AS) events in reeds based on a transcriptome database. Meanwhile, we have identified and developed for the first time a large number of candidates expressed sequence tag-SSR (EST-SSRs) markers in reeds based on UniTransModels. In addition, through differential gene expression analysis of wild-type and homogenous cultures, we found a large number of transcription factors that may be associated with desert stress tolerance in the dune reed, and revealed that members of the Lhc family have an important role in the long-term adaptation of dune reeds to desert environments.

CONCLUSIONS

Our results provide a positive and usable genetic resource for Phragmites australis with a widespread adaptability and resistance, and provide a genetic database for subsequent reeds genome annotation and functional genomic studies.

Development of unigene-derived SSR markers from RNA-seq data of Uraria lagopodioides (Fabaceae) and their application in the genus Uraria Desv. (Fabaceae)

BACKGROUND

Uraria Desv. belongs to the tribe Desmodieae (Fabaceae), a group of legume plants, some of which have medicinal properties. However, due to a lack of genomic information, the interspecific relationships, genetic diversity, population genetics, and identification of functional genes within Uraria species are still unclear.

RESULTS

Using RNA-Seq, a total of 66,026 Uraria lagopodioides unigenes with a total sequence content of 52,171,904 bp were obtained via de novo assembly and annotated using GO, KEGG, and KOG databases. 17,740 SSRs were identified from a set of 66,026 unigenes. Cross-species amplification showed that 54 out of 150 potential unigene-derived SSRs were transferable in Uraria, of which 19 polymorphic SSRs were developed. Cluster analysis based on polymorphisms successfully distinguished seven Uraria species and revealed their interspecific relationships. Seventeen samples of seven Uraria species were clustered into two monophyletic clades, and phylogenetic relationships of Uraria species based on unigene-derived SSRs were consistent with classifications based on morphological characteristics.

CONCLUSIONS

Unigenes annotated in the present study will provide new insights into the functional genomics of Uraria species. Meanwhile, the unigene-derived SSR markers developed here will be invaluable for assessing the genetic diversity and evolutionary history of Uraria and relatives.

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

Genetic diversity analysis and fingerprint construction of Korean pine (Pinus koraiensis) clonal seed orchard

Korean pine is a native tree species in Northeast China. In order to meet the needs of germplasm resource evaluation and molecular marker-assisted breeding of Korean pine, we collected Korean pine clones from 7 populations in Northeast China, analyzed the genetic diversity and genetic structure by SSR molecular marker technology and clustered them to revealed the inter- and intrapopulation differentiation characteristics of each clone. The fingerprint profiles of 161 Korean pine clones were also constructed. 77 alleles were detected for 11 markers, and 18 genotypes were identified on average for each marker. The PIC of the different markers ranged from 0.155-0.855, and the combination of PI and PIsibs for the 11 markers was 3.1 × 10^-8 and 1.14 × 10^-3, respectively. MANOVA showed that genetic variation existed mainly within populations, accounting for 98% of the total variation. The level of genetic differentiation among populations was low, with an average Nm between populations of 11.036. Genetic diversity is lower in the Lushuihe population and higher in the Tieli population. The 161 Korean pine clones were divided into 4 or 7 populations, and the 7 populations were not clearly distinguished from each other, with only the Lushuihe population showing partial differentiation. There is no significant correlation between the genetic distance of Korean pine populations and the geographical distance of their superior tree sources. This result can provide recommendations for future Korean pine breeding programs. The combination of 11 markers could completely distinguish 161 clones and establish the fingerprint. Genetic diversity of Korean pine clones from the 7 populations was abundant, and the genetic distances of individuals and populations were evenly dispersed. The fingerprint map can be used for the identification of Korean pine clones.

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.

Antioxidant processes involving epicatechin decreased symptoms of pine wilt disease

Since the pine wood nematode (PWN, Bursaphelenchus xylophilus) invasion of Northeast China, both symptomatic and asymptomatic PWN carriers have been found. Asymptomatic PWN carriers, which are more dangerous than symptomatic carriers, constitute a source of infection in the following spring. The simultaneous presence of symptomatic and asymptomatic PWN carriers indicates that Pinus koraiensis has different tolerance levels to PWN. In this study, validity of susceptibility testing discovered differential types of P. koraiensis including Latent Reservoirs, Low Susceptibles, High Susceptibles and Bell Ringers. Among those types, the Low Susceptibles and Latent Reservoirs were asymptomatic PWN carriers, and Latent Reservoirs were the most dangerous. Transcriptome and metabolomic data showed that 5 genes (3 ans and 2 anr gene) involved in the epicatechin (EC) synthesis pathway were significantly upregulated, which increased the content of EC antioxidants in Latent Reservoirs. Hydrogen peroxide (H₂O₂) staining and content determination showed that the hypersensitive response (HR) and H₂O₂, which functions as a signaling molecule in systemic acquired resistance, decreased in Latent Reservoirs. However, low contents of EC and high contents of H₂O₂ were found in the High Susceptibles of P. koraiensis. RT-PCR results showed that the expression of ans and anr was upregulated together only in Latent Reservoirs. These results show that the susceptibility of P. koraiensis to PWN differed among different individuals, although no resistant individuals were found. Latent Reservoirs, in which more PWNs resided without visible symptoms via prolonged incubation period, inhibited the symptoms caused by H₂O₂ because of increased contents of the EC antioxidants.

Analysis of the genetic diversity and population structure of Monochasma savatieri Franch. ex Maxim using novel EST-SSR markers

Background

Monochasma savatieri Franch. ex Maxim is a medicinally valuable herb. However, the collection and protection of the wild germplasm resources of M. savatieri are still insufficient, and their genetic diversity and population structure have been poorly studied.

Results

We collected and examined 46 M. savatieri individuals from Fujian, Hunan, Jiangxi, and Zhejiang provinces for genetic diversity and population structure, using 33 newly developed expressed sequence tag-simple sequence repeat (EST-SSR) markers. Applying these markers, we detected a total of 208 alleles, with an average of 6.303 alleles per locus. The polymorphic information content varied from 0.138 to 0.884 (average: 0.668), indicating a high level of polymorphism. At the population level, there was a low degree of genetic diversity among populations (I = 0.535, He = 0.342), with Zhejiang individuals showing the highest genetic diversity among the four populations (Fst = 0.497), which indicated little gene flow within the M. savatieri populations (Nm = 0.253). Mantel test analysis revealed a significant positive correlation between geographical and genetic distance among populations (R² = 0.3304, p < 0.05), and structure and principal coordinate analyses supported classification of populations into three clusters, which was consistent with the findings of cluster analysis.

Conclusions

As a rare medicinal plants, the protection of M. savatieri does not look optimistic, and accordingly, protective efforts should be beefed up on the natural wild populations. This study provided novel tools and insights for designing effective collection and conservation strategies for M. savatieri.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08832-x.

Root Breeding in the Post-Genomics Era: From Concept to Practice in Apple

The development of rootstocks with a high-quality dwarf-type root system is a popular research topic in the apple industry. However, the precise breeding of rootstocks is still challenging, mainly because the root system is buried deep underground, roots have a complex life cycle, and research on root architecture has progressed slowly. This paper describes ideas for the precise breeding and domestication of wild apple resources and the application of key genes. The primary goal of this research is to combine the existing rootstock resources with molecular breeding and summarize the methods of precision breeding. Here, we reviewed the existing rootstock germplasm, high-quality genome, and genetic resources available to explain how wild resources might be used in modern breeding. In particular, we proposed the 'from genotype to phenotype' theory and summarized the difficulties in future breeding processes. Lastly, the genetics governing root diversity and associated regulatory mechanisms were elaborated on to optimize the precise breeding of rootstocks.

Development and Application of EST-SSR Markers in Cephalotaxus oliveri From Transcriptome Sequences

Cephalotaxus oliveri is an endemic conifer of China, which has medicinal and ornamental value. However, the limited molecular markers and genetic information are insufficient for further genetic studies of this species. In this study, we characterized and developed the EST-SSRs from transcriptome sequences for the first time. The results showed that a total of 5089 SSRs were identified from 36446 unigenes with a density of one SSR per 11.1 kb. The most common type was trinucleotide repeats, excluding mononucleotide repeats, followed by dinucleotide repeats. AAG/CTT and AT/AT exhibited the highest frequency in the trinucleotide and dinucleotide repeats, respectively. Of the identified SSRs, 671, 1125, and 1958 SSRs were located in CDS, 3′UTR, and 5′UTR, respectively. Functional annotation showed that the SSR-containing unigenes were involved in growth and development with various biological functions. Among successfully designed primer pairs, 238 primer pairs were randomly selected for amplification and validation of EST-SSR markers and 47 primer pairs were identified as polymorphic. Finally, 28 high-polymorphic primers were used for genetic analysis and revealed a moderate level of genetic diversity. Seven natural C. oliveri sampling sites were divided into two genetic groups. Furthermore, the 28 EST-SSRs had 96.43, 71.43, and 78.57% of transferability rate in Cephalotaxus fortune, Ametotaxus argotaenia, and Pseudotaxus chienii, respectively. These markers developed in this study lay the foundation for further genetic and adaptive evolution studies in C. oliveri and related species.

An Overview of the Practices and Management Methods for Enhancing Seed Production in Conifer Plantations for Commercial Use

Flowering, the beginning of the reproductive growth, is a significant stage in the growth and development of plants. Conifers are economically and ecologically important, characterized by straight trunks and a good wood quality and, thus, conifer plantations are widely distributed around the world. In addition, conifer species have a good tolerance to biotic and abiotic stress, and a stronger survival ability. Seeds of some conifer species, such as Pinus koraiensis, are rich in vitamins, amino acids, mineral elements and other nutrients, which are used for food and medicine. Although conifers are the largest (giant sequoia) and oldest living plants (bristlecone pine), their growth cycle is relatively long, and the seed yield is unstable. In the present work, we reviewed selected literature and provide a comprehensive overview on the most influential factors and on the methods and techniques that can be adopted in order to improve flowering and seed production in conifers species. The review revealed that flowering and seed yields in conifers are affected by a variety of factors, such as pollen, temperature, light, water availability, nutrients, etc., and a number of management techniques, including topping off, pruning, fertilization, hormone treatment, supplementary pollination, etc. has been developed for improving cone yields. Furthermore, several flowering-related genes (FT, Flowering locus T and MADS-box, MCMI, AGAMOUS, DEFICIENCES and SRF) that play a crucial role in flowering in coniferous trees were identified. The results of this study can be useful for forest managers and for enhancing seed yields in conifer plantations for commercial use.

De Novo Transcriptome Assembly, Functional Annotation and SSR Marker Discovery of Qinling Takin (Budorcas taxicolor bedfordi)

The takin (Budorcas taxicolor) is an endemic ruminant species belonging to the bovine family. The International Union for Conservation of Nature (IUCN) has listed it as an endangered and vulnerable species. However, little is known about its molecular characterization since it lacks a reference genome. This study used RNA sequencing followed by de novo assembly, annotation and simple sequence repeats (SSRs) prediction to assess the transcriptome of Qinling takin (Budorcas taxicolor bedfordi) muscles. In total, 21,648 unigenes with an N50 and mean length of 1388 bp and 817 bp, respectively, were successfully detected and annotated against the public databases (NR, GO, KEGG, and EggNOG). Furthermore, 6222 SSRs were identified using the MIcroSAtellite (MISA) identification tool software. Taken together, these findings will provide valuable information for genetic, genomic, and evolutionary studies on takin.

Genetic Diversity and Population Differentiation of Pinus koraiensis in China

Pinus koraiensis is a well-known precious tree species in East Asia with high economic, ornamental and ecological value. More than fifty percent of the P. koraiensis forests in the world are distributed in northeast China, a region with abundant germplasm resources. However, these natural P. koraiensis sources are in danger of genetic erosion caused by continuous climate changes, natural disturbances such as wildfire and frequent human activity. Little work has been conducted on the population genetic structure and genetic differentiation of P. koraiensis in China because of the lack of genetic information. In this study, 480 P. koraiensis individuals from 16 natural populations were sampled and genotyped. Fifteen polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers were used to evaluate genetic diversity, population structure and differentiation in P. koraiensis. Analysis of molecular variance (AMOVA) of the EST-SSR marker data showed that 33% of the total genetic variation was among populations and 67% was within populations. A high level of genetic diversity was found across the P. koraiensis populations, and the highest levels of genetic diversity were found in HH, ZH, LS and TL populations. Moreover, pairwise Fst values revealed significant genetic differentiation among populations (mean Fst = 0.177). According to the results of the STRUCTURE and Neighbor-joining (NJ) tree analyses and principal component analysis (PCA), the studied geographical populations cluster into two genetic clusters: cluster 1 from Xiaoxinganling Mountains and cluster 2 from Changbaishan Mountains. These results are consistent with the geographical distributions of the populations. The results provide new genetic information for future genome-wide association studies (GWAS), marker-assisted selection (MAS) and genomic selection (GS) in natural P. koraiensis breeding programs and can aid the development of conservation and management strategies for this valuable conifer species.

Characterization and Development of Microsatellite Markers in Pseudotaxus chienii (Taxaceae) Based on Transcriptome Sequencing

Pseudotaxus chienii (Taxaceae) is an endangered conifer species endemic to China. However, a lack of suitable molecular markers hinders the genomic and genetic studies on this species. Here, we characterized and developed the microsatellite markers from a newly sequenced P. chienii transcriptome. A total of 21,835 microsatellite loci were detected from 161,131 non-redundant unigene sequences, and the frequency of SSRs was 13.55%, with an average of one SSR loci per 9.18 kb. Mono-nucleotide, di-nucleotide, and tri-nucleotide were the dominant repeat types, accounting for 50.06, 13.49, and 29.39% of the total SSRs, respectively. In terms of distribution location, the coding regions (CDS) with few microsatellites and mainly consisted of tri-nucleotides. There were significant differences in the length of microsatellite among genic regions and motif types. Functional annotation showed that the unigenes containing microsatellites had a wide range of biological functions, most of which were related to basic metabolism, and a few might be involved in expression regulation of gene and response to environmental stress. In addition, 375 primer pairs were randomly selected and synthesized for the amplification and validation of microsatellite markers. Seventy-seven primer pairs were successfully amplified and 40 primer pairs were found to be polymorphic. Finally, 20 pairs of primers with high polymorphism were selected to assess the genetic diversity in four P. chienii populations. In addition, the newly developed microsatellite markers exhibited high transferability (70%) in Amentotaxus argotaenia. Our study could enable further genetic diversity analysis and functional gene mining on Taxaceae.