Genetic diversity and population structure of black cottonwood (Populus deltoides) revealed using simple sequence repeat markers

Population genetic variation and geographic distribution of suitable areas of Coptis species in China

Introduction

The rhizomes of Coptis plants have been used in traditional Chinese medicine over 2000 years. Due to increasing market demand, the overexploitation of wild populations, habitat degradation and indiscriminate artificial cultivation of Coptis species have severely damaged the native germplasms of species in China.

Methods

Genome-wide simple-sequence repeat (SSR) markers were developed using the genomic data of C. chinensis. Population genetic diversity and structure of 345 Coptis accessions collected from 32 different populations were performed based on these SSRs. The distribution of suitable areas for three taxa in China was predicted and the effects of environmental variables on genetic diversity in relation to different population distributions were further analyzed.

Results

22 primer pairs were selected as clear, stable, and polymorphic SSR markers. These had an average of 16.41 alleles and an average polymorphism information content (PIC) value of 0.664. In the neighbor-joining (N-J) clustering analysis, the 345 individuals clustered into three groups, with C. chinensis, C. chinensis var. brevisepala and C. teeta being clearly separated. All C. chinensis accessions were further divided into four subgroups in the population structure analysis. The predicted distributions of suitable areas and the environmental variables shaping these distributions varied considerably among the three species.

Discussion

Overall, the amount of solar radiation, precipitation and altitude were the most important environmental variables influencing the distribution and genetic variation of three species. The findings will provide key information to guide the conservation of genetic resources and construction of a core reserve for species.

Interspecific variations in growth, physiology and Cd accumulation between Populus deltoides and P. × canadensis in response to Cd pollution under two soil types

Both acid and alkaline purple soils in China are increasingly affected by Cd contamination. The selection of fast-growing trees suitable for remediating different soil types is urgent, yet there is a severe lack of relevant knowledge. In this study, we conducted a controlled pot experiment to compare the growth, physiology, and Cd accumulation efficiency of two widely recognized poplar species, namely Populus deltoides and P. × canadensis, under Cd contamination (1 mg kg-1) in acid and alkaline purple soils. The objective was to determine which poplar species is best suited for remediating different soil types. Our findings are as follows: (1) the total biomass of both poplars remained largely unaffected by Cd pollution in both soil types. Notably, under Cd pollution, the total biomass of P. deltoides in acid purple soil was 1.53 times greater than that in alkaline purple soil. (2) Cd pollution did not significantly induce oxidative damage in the leaves of either poplar species in both soil types. However, in acid purple soil, Cd contamination led to a 21% increase in NO3- concentration and a 44% increase in NH4+ concentration in P. × canadensis leaves, whereas in alkaline purple soil, it led to a 59% increase in NH4+ concentration in P. deltoides leaves. (3) Cd concentrations in all root orders of P. × canadensis were significantly higher than those in P. deltoides, especially in the first three root orders, under alkaline purple soil. The total Cd accumulation by P. × canadensis in Cd-polluted alkaline purple soil was 2.18 times higher than that in Cd-polluted acid purple soil, a difference not observed in P. deltoides. (4) redundancy analysis indicated that the sequestration effect of higher soil organic matter on Cd availability in acid purple soil was more pronounced than the release effects caused by lower pH. In conclusion, P. × canadensis is better suited for remediating alkaline purple soil due to its higher capacity for Cd uptake, while P. deltoides is more suitable for remediating Cd-contaminated acid purple soil due to its better growth conditions and greater Cd enrichment capability.

De Novo Mining and Validating Novel Microsatellite Markers to Assess Genetic Diversity in Maruca vitrata (F.), a Legume Pod Borer

Maruca vitrata (Fabricius) is an invasive insect pest capable of causing enormous economic losses to a broad spectrum of leguminous crops. Microsatellites are valuable molecular markers for population genetic studies; however, an inadequate number of M. vitrata microsatellite loci are available to carry out population association studies. Thus, we utilized this insect's public domain databases for mining expressed sequence tags (EST)-derived microsatellite markers. In total, 234 microsatellite markers were identified from 10053 unigenes. We discovered that trinucleotide repeats were the most predominant microsatellite motifs (61.53%), followed by dinucleotide repeats (23.50%) and tetranucleotide repeats (14.95%). Based on the analysis, twenty-five markers were selected for validation in M. vitrata populations collected from various regions of India. The number of alleles (Na), observed heterozygosity (Ho), and expected heterozygosity (He) ranged from 2 to 5; 0.00 to 0.80; and 0.10 to 0.69, respectively. The polymorphic loci showed polymorphism information content (PIC), ranging from 0.09 to 0.72. Based on the genetic distance matrix, the unrooted neighbor-joining dendrogram differentiated the selected populations into two discrete groups. The SSR markers developed and validated in this study will be helpful in population-level investigations of M. vitrata to understand the gene flow, demography, dispersal patterns, biotype differentiation, and host dynamics.

eYGFPuv-Assisted Transgenic Selection in Populus deltoides WV94 and Multiplex Genome Editing in Protoplasts of P. trichocarpa × P. deltoides Clone '52-225'

Although CRISPR/Cas-based genome editing has been widely used for plant genetic engineering, its application in the genetic improvement of trees has been limited, partly because of challenges in Agrobacterium-mediated transformation. As an important model for poplar genomics and biotechnology research, eastern cottonwood (Populus deltoides) clone WV94 can be transformed by A. tumefaciens, but several challenges remain unresolved, including the relatively low transformation efficiency and the relatively high rate of false positives from antibiotic-based selection of transgenic events. Moreover, the efficacy of CRISPR-Cas system has not been explored in P. deltoides yet. Here, we first optimized the protocol for Agrobacterium-mediated stable transformation in P. deltoides WV94 and applied a UV-visible reporter called eYGFPuv in transformation. Our results showed that the transgenic events in the early stage of transformation could be easily recognized and counted in a non-invasive manner to narrow down the number of regenerated shoots for further molecular characterization (at the DNA or mRNA level) using PCR. We found that approximately 8.7% of explants regenerated transgenic shoots with green fluorescence within two months. Next, we examined the efficacy of multiplex CRISPR-based genome editing in the protoplasts derived from P. deltoides WV94 and hybrid poplar clone '52-225' (P. trichocarpa × P. deltoides clone '52-225'). The two constructs expressing the Trex2-Cas9 system resulted in mutation efficiency ranging from 31% to 57% in hybrid poplar clone 52-225, but no editing events were observed in P. deltoides WV94 transient assay. The eYGFPuv-assisted plant transformation and genome editing approach demonstrated in this study has great potential for accelerating the genome editing-based breeding process in poplar and other non-model plants species and point to the need for additional CRISPR work in P. deltoides.

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.

Genetic Diversity and Population Structure Analysis of Castanopsis hystrix and Construction of a Core Collection Using Phenotypic Traits and Molecular Markers

Castanopsis hystrix is a valuable native, broad-leaved, and fast-growing tree in South China. In this study, 15 phenotypic traits and 32 simple sequence repeat (SSR) markers were used to assess the genetic diversity and population structure of a natural population of C. hystrix and to construct a core germplasm collection by a set of 232 accessions. The results showed that the original population of C. hystrix had relatively high genetic diversity, with the number of alleles (Na), effective number of alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), Shannon's information index (I), and polymorphism information content (PIC) averaging at 26.188, 11.565, 0.863, 0.897, 2.660, and 0.889, respectively. Three sub-populations were identified based on a STRUCTURE analysis, indicating a strong genetic structure. The results from the phylogenetic and population structures showed a high level of agreement, with 232 germplasms being classified into three main groups. The analysis of molecular variance (AMOVA) test indicated that 96% of the total variance was derived from within populations, which revealed a low differentiation among populations. A core collection composed of 157 germplasms was firstly constructed thereafter, of which the diversity parameters non-significantly differed from the original population. These results revealed the genetic diversity and population structure of C. hystrix germplasms, which have implications for germplasm management and genome-wide association studies on C. hystrix, as well as for core collection establishment applications in other wood-producing hardwood species.

A Neural Network-Based Spectral Approach for the Assignment of Individual Trees to Genetically Differentiated Subpopulations

Studying population structure has made an essential contribution to understanding evolutionary processes and demographic history in forest ecology research. This inference process basically involves the identification of common genetic variants among individuals, then grouping the similar individuals into subpopulations. In this study, a spectral-based classification of genetically differentiated groups was carried out using a provenance–progeny trial of Eucalyptus cladocalyx. First, the genetic structure was inferred through a Bayesian analysis using single-nucleotide polymorphisms (SNPs). Then, different machine learning models were trained with foliar spectral information to assign individual trees to subpopulations. The results revealed that spectral-based classification using the multilayer perceptron method was very successful at classifying individuals into their respective subpopulations (with an average of 87% of correct individual assignments), whereas 85% and 81% of individuals were assigned to their respective classes correctly by convolutional neural network and partial least squares discriminant analysis, respectively. Notably, 93% of individual trees were assigned correctly to the class with the smallest size using the spectral data-based multi-layer perceptron classification method. In conclusion, spectral data, along with neural network models, are able to discriminate and assign individuals to a given subpopulation, which could facilitate the implementation and application of population structure studies on a large scale.

Contrasting Reproductive Strategies of Two Nymphaea Species Affect Existing Natural Genetic Diversity as Assessed by Microsatellite Markers: Implications for Conservation and Wetlands Restoration

Nymphaea, commonly known as water lily, is the largest and most widely distributed genus in the order Nymphaeales. The importance of Nymphaea in wetland ecosystems and their increased vulnerability make them a great choice for conservation and management. In this work, we studied genetic diversity in a collection of 90 N. micrantha and 92 N. nouchali individuals from six different states of India, i.e., Assam, Manipur, Meghalaya, Maharashtra, Goa, and Kerala, using simple sequence repeat (SSR) markers developed by low throughput Illumina sequencing (10X coverage of genome) of N. micrantha. Nymphaea nouchali is native to India, whereas N. micrantha is suggested to be introduced to the country for its aesthetic and cultural values. The study revealed extensive polymorphism in N. nouchali, while in N. micrantha, no apparent genetic divergence was detected prompting us to investigate the reason(s) by studying the reproductive biology of the two species. The study revealed that N. micrantha predominantly reproduces asexually which has impacted the genetic diversity of the species to a great extent. This observation is of immense importance for a successful re-establishment of Nymphaea species during restoration programs of wetlands. The information generated on reproductive behaviors and their association with genotypic richness can help in strategizing genetic resource conservation, especially for species with limited distribution. The study has also generated 22,268 non-redundant microsatellite loci, out of which, 143 microsatellites were tested for polymorphism and polymorphic markers were tested for transferability in five other Nymphaea species, providing genomic resources for further studies on this important genus.

EST-SSR Primer Development and Genetic Structure Analysis of Psathyrostachys juncea Nevski

Psathyrostachys juncea is a perennial forage grass which plays an important role in soil and water conservation and ecological maintenance in cold and dry areas of temperate regions. In P. juncea, a variety of biotic and abiotic stress related genes have been used in crop improvement, indicating its agronomic, economic, forage, and breeding value. To date, there have been few studies on the genetic structure of P. juncea. Here, the genetic diversity and population structure of P. juncea were analyzed by EST-SSR molecular markers to evaluate the genetic differentiation related to tillering traits in P. juncea germplasm resources. The results showed that 400 simple sequence repeat (SSR) loci were detected in 2,020 differentially expressed tillering related genes. A total of 344 scored bands were amplified using 103 primer pairs, out of which 308 (89.53%) were polymorphic. The Nei's gene diversity of 480 individuals was between 0.092 and 0.449, and the genetic similarity coefficient was between 0.5008 and 0.9111, with an average of 0.6618. Analysis of molecular variance analysis showed that 93% of the variance was due to differences within the population, and the remaining 7% was due to differences among populations. Psathyrostachys juncea materials were clustered into five groups based on population genetic structure, principal coordinate analysis and unweighted pair-group method with arithmetic means (UPGMA) analysis. The results were similar between clustering methods, but a few individual plants were distributed differently by the three models. The clustering results, gene diversity and genetic similarity coefficients showed that the overall genetic relationship of P. juncea individuals was relatively close. A Mantel test, UPGMA and structural analysis also showed a significant correlation between genetic relationship and geographical distribution. These results provide references for future breeding programs, genetic improvement and core germplasm collection of P. juncea.

Genetic diversity and population structure analysis in a large collection of Vicia amoena in China with newly developed SSR markers

Vicia amoena is a high-nutritional quality forage similar to alfalfa. However, studies on the genetic background of V. amoena are scarce. In the present study, the genetic variation of 24 V. amoena populations was assessed with newly developed simple sequence repeat (SSR) markers. A total of 8799 SSRs were identified in the V. amoena genomic-enriched sequences, and the most abundant repeat number was four. A total of 569 sampled individuals were assayed to evaluate the genetic diversity of the V. amoena populations based on 21 polymorphic SSR primers. The polymorphism information content (PIC) ranged from 0.896 to 0.968, with an average of 0.931, which indicated that the markers were highly informative. Based on analysis of molecular variance, 88% of the variance occurred within populations, and the remaining 12% of the variance occurred among populations. The high degree of gene flow (Nm= 4.958) also showed slight differentiation among the V. amoena populations. The V. amoena populations were mainly clustered by steppe and mountain habitats based on principal coordinate analysis (PCoA) and STRUCTURE analysis. This indicated that the elevation and special habitat of geographical origins may be important factors affecting the clustered pattern of V. amoena populations. Neighbour-joining (NJ) analysis did not separate the populations well by geographical origin, which indicated that the genetic structure of V. amoena was complex and needs further study. Overall, our results showed that the newly developed SSR markers could benefit the V. amoena research community by providing genetic background information to help establish a foundation for breeding improvement and germplasm resource conservation.

Morphological, physiological, and transcriptional responses to low nitrogen stress in Populus deltoides Marsh. clones with contrasting nitrogen use efficiency

Background

Nitrogen (N) is one of the main factors limiting the wood yield in poplar cultivation. Understanding the molecular mechanism of N utilization could play a guiding role in improving the nitrogen use efficiency (NUE) of poplar.

Results

In this study, three N-efficient genotypes (A1-A3) and three N-inefficient genotypes (C1-C3) of Populus deltoides were cultured under low N stress (5 μM NH₄NO₃) and normal N supply (750 μM NH₄NO₃). The dry matter mass, leaf morphology, and chlorophyll content of both genotypes decreased under N starvation. The low nitrogen adaptation coefficients of the leaves and stems biomass of group A were significantly higher than those of group C (p < 0.05). Interestingly, N starvation induced fine root growth in group A, but not in group C. Next, a detailed time-course analysis of enzyme activities and gene expression in leaves identified 2062 specifically differentially expressed genes (DEGs) in group A and 1118 in group C. Moreover, the sensitivity to N starvation of group A was weak, and DEGs related to hormone signal transduction and stimulus response played an important role in the low N response this group. Weighted gene co-expression network analysis identified genes related to membranes, catalytic activity, enzymatic activity, and response to stresses that might be critical for poplar’s adaption to N starvation and these genes participated in the negative regulation of various biological processes. Finally, ten influential hub genes and twelve transcription factors were identified in the response to N starvation. Among them, four hub genes were related to programmed cell death and the defense response, and PodelWRKY18, with high connectivity, was involved in plant signal transduction. The expression of hub genes increased gradually with the extension of low N stress time, and the expression changes in group A were more obvious than those in group C.

Conclusions

Under N starvation, group A showed stronger adaptability and better NUE than group C in terms of morphology and physiology. The discovery of hub genes and transcription factors might provide new information for the analysis of the molecular mechanism of NUE and its improvement in poplar.

Supplementary Information

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

Population Genetic Diversity and Structure of an Endangered Salicaceae Species in Northeast China: Chosenia arbutifolia (Pall.) A. Skv

Chosenia arbutifolia (Pall.) A. Skv. is a unique and endangered species belonging to the Salicaceae family. It has great potential for ornamental and industrial use. However, human interference has led to a decrease in and fragmentation of its natural populations in the past two decades. To effectively evaluate, utilize, and conserve available resources, the genetic diversity and population structure of C. arbutifolia were analyzed in this study. A total of 142 individuals from ten provenances were sampled and sequenced. Moderate diversity was detected among these, with a mean expected heterozygosity and Shannon’s Wiener index of 0.3505 and 0.5258, respectively. The inbreeding coefficient was negative, indicating a significant excess of heterozygotes. The fixation index varied from 0.0068 to 0.3063, showing a varied genetic differentiation between populations. Analysis of molecular variance demonstrated that differentiation accounted for 82.23% of the total variation among individuals, while the remaining 17.77% variation was between populations. Furthermore, the results of population structure analysis indicated that the 142 individuals originated from three primitive groups. To provide genetic information and help design conservation and management strategies, landscape genomics analysis was performed by investigating loci associated with environmental variables. Eighteen SNP markers were associated with altitude and annual average temperature, of which five were ascribed with specific functions. In conclusion, the current study furthers the understanding of C. arbutifolia genetic architecture and provides insights for germplasm protection.

Genetic Diversity and Population Structure of Capirona (Calycophyllum spruceanum Benth.) from the Peruvian Amazon Revealed by RAPD Markers

Capirona (Calycophyllum spruceanum Benth.) is a tree species of commercial importance widely distributed in South American forests that is traditionally used for its medicinal properties and wood quality. Studies on this tree species have been focused mainly on wood properties, propagation, and growth. However, genetic studies on capirona have been very limited to date. Currently, it is possible to explore genetic diversity and population structure in a fast and reliable manner by using molecular markers. We here used 10 random amplified polymorphic DNA (RAPD) markers to analyze the genetic diversity and population structure of 59 samples of capirona that were sampled from four provinces located in the eastern region of the Peruvian amazon. A total of 186 bands were manually scored, generating a 59 × 186 presence/absence matrix. A dendrogram was generated using the UPGMA clustering algorithm, and, similar to the principal coordinate analysis (PCoA), it showed four groups that correspond to the geographic origin of the capirona samples (LBS, Irazola, Masisea, Iñapari). Similarly, a discriminant analysis of principal components (DAPC) and STRUCTURE analysis confirmed that capirona is grouped into four clusters. However, we also noticed that a few samples were intermingled. Genetic diversity estimation was conducted considering the four groups (populations) identified by STRUCTURE software. AMOVA revealed the greatest variation within populations (71.56%) and indicated that variability among populations is 28.44%. Population divergence (Fst) between clusters 1 and 4 revealed the highest genetic difference (0.269), and the lowest Fst was observed between clusters 3 and 4 (0.123). RAPD markers were successful and effective. However, more studies are needed, employing other molecular tools. To the best of our knowledge, this is the first investigation employing molecular markers in capirona in Peru considering its natural distribution, and as such it is hoped that this helps to pave the way towards its genetic improvement and the urgent sustainable management of forests in Peru.

Genome skimming-based simple sequence repeat (SSR) marker discovery and characterization in Grevillea robusta

Proteaceae, a largely southern hemisphere family consisting of 80 genera distributed in Australia and southern Africa as its centres of greatest diversity, also extends well in northern and southern America. Under this family, Grevillea robusta is a fast-growing species got popularity in farm and avenue plantations. Despite the ecological and economic importance, the species has not yet been investigated for its genetic improvement and genome-based studies. Only a few molecular markers are available for the species or its close relatives, which hinders  genomic and population genetics studies. Genetic markers have been intensively applied for the main strategies in breeding programs, especially for the economically important traits. Hence, it is of utmost priority to develop genomic database resources and species-specific markers for studying quantitative genetics in G. robusta. Given this, the present study aimed to develop de novo genome sequencing, robust microsatellites markers, sequence annotation and their validation in different stands of G. robusta in northern India. Library preparation and sequencing were carried out using Illumina paired-end sequencing technology. Approximately, ten gigabases (Gb) sequence data with 70.87 million raw reads assembled into 425,923 contigs (read mapped to 76.48%) comprising 455 Mb genome size (23 × coverage) generated through genome skimming approach. In total, 9421 simple sequence repeat (SSR) primer pairs were successfully designed from 13,335 microsatellite repeats. Afterward, a subset of 161 primer pairs was randomly selected, synthesized and validated. All the tested primers showed successful amplification but only 13 showed polymorphisms. The polymorphic SSRs were further used to estimate the measures of genetic diversity in 12 genotypes each from the states of Punjab, Haryana, Himachal Pradesh and Uttarakhand. Importantly, the average number of alleles (Na), observed heterozygosity (Ho), expected heterozygosity (He), and the polymorphism information content (PIC) were recorded as 2.69, 0.356, 0.557 and 0.388, respectively. The availability of sequence information and newly developed SSR markers could potentially be used in various genetic analyses and improvements through molecular breeding strategies for G. robusta.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01035-w.

Genome Assembly of Salicaceae Populus deltoides (Eastern Cottonwood) I-69 Based on Nanopore Sequencing and Hi-C Technologies

Populus deltoides has important ecological and economic values, widely used in poplar breeding programs due to its superior characteristics such as rapid growth and resistance to disease. Although the genome sequence of P. deltoides WV94 is available, the assembly is fragmented. Here, we reported an improved chromosome-level assembly of the P. deltoides cultivar I-69 by combining Nanopore sequencing and chromosome conformation capture (Hi-C) technologies. The assembly was 429.3 Mb in size and contained 657 contigs with a contig N50 length of 2.62 Mb. Hi-C scaffolding of the contigs generated 19 chromosome-level sequences, which covered 97.4% (418 Mb) of the total assembly size. Moreover, repetitive sequences annotation showed that 39.28% of the P. deltoides genome was composed of interspersed elements, including retroelements (23.66%), DNA transposons (6.83%), and unclassified elements (8.79%). We also identified a total of 44 362 protein-coding genes in the current P. deltoides assembly. Compared with the previous genome assembly of P. deltoides WV94, the current assembly had some significantly improved qualities: the contig N50 increased 3.5-fold and the proportion of gaps decreased from 3.2% to 0.08%. This high-quality, well-annotated genome assembly provides a reliable genomic resource for identifying genome variants among individuals, mining candidate genes that control growth and wood quality traits, and facilitating further application of genomics-assisted breeding in populations related to P. deltoides.