Assessment of genetic diversity among four orchids based on ddRAD sequencing data for conservation purposes

Conservation genetics and potential geographic distribution modeling of Corybas taliensis, a small 'sky Island' orchid species in China

BACKGROUND

Corybas taliensis is an endemic species of sky islands in China. Its habitat is fragile and unstable, and it is likely that the species is threatened. However, it is difficult to determine the conservation priority or unit without knowing the genetic background and the overall distribution of this species. In this study, we used double digest restriction-site associated DNA-sequencing (ddRAD-seq) to investigate the conservation genomics of C. taliensis. At the same time, we modeled the extent of suitable habitat for C. taliensis in present and future (2030 and 2090) habitat using the maximum-entropy (MaxEnt) model.

RESULTS

The results suggested that the related C. fanjingshanensis belongs to C. taliensis and should not be considered a separate species. All the sampling locations were divided into three genetic groups: the Sichuan & Guizhou population (SG population), the Hengduan Mountains population (HD population) and Himalayan population (HM population), and we found that there was complex gene flow between the sampling locations of HD population. MT was distinct genetically from the other sampling locations due to the unique environment in Motuo. The genetic diversity (π, He) of C. taliensis was relatively high, but its contemporary effective population size (Ne) was small. C. taliensis might be currently affected by inbreeding depression, although its large population density may be able to reduce the effect of this. The predicted areas of suitable habitat currently found in higher mountains will not change significantly in the future, and these suitable habitats are predicted to spread to other higher mountains under future climate change. However, suitable habitat in relatively low altitude areas may disappear in the future. This suggests that C. taliensis will be caught in a 'summit trap' in low altitude areas, however, in contrast, the high altitude of the Himalaya and the Hengduan Mountains are predicted to act as 'biological refuges' for C. taliensis in the future.

CONCLUSIONS

These results not only provide a new understanding of the genetic background and potential resource distribution of C. taliensis, but also lay the foundation for its conservation and management.

The Genomics of Isolated Populations of Gampsocleis glabra (Orthoptera: Tettigoniidae) in Central and Western Europe

Habitat destruction and fragmentation are among the major current threats to global biodiversity. Fragmentation may also affect species with good dispersal abilities. We study the heath bushcricket Gampsocleis glabra, a specialist of steppe-like habitats across Europe that are highly fragmented, investigating if these isolated populations can be distinguished using population genomics and if there are any traces of admixture or dispersal among them. We try to answer these questions using genome-wide SNP data generated with ddRAD sequencing. We calculated F-statistics and visualized differentiation using STRUCTURE plots. While limited by the difficulty of sampling this threatened species, our results show that all populations except one that was represented by a singleton were clearly distinct, with pairwise FST values between 0.010 and 0.181. STRUCTURE indicated limited but visible admixture across most populations and probably also an exchange of individuals between populations of Germany and The Netherlands. We conclude that in G. glabra, a certain amount of gene flow has persisted, at least in the past, also among populations that are isolated today. We also detect a possibly more recent dispersal event between a population in The Netherlands and one in Germany, which may be human aided. We suggest that the conservation of larger populations should be maintained, that efforts should be taken to restore abandoned habitat, that the preservation even of small habitat fragments may be beneficial for the conservation of this species, and that these habitats should be regularly monitored for possible (re-)colonization.

Genomic evidence reveals high genetic diversity in a narrowly distributed species and natural hybridization risk with a widespread species in the genus Geodorum

BACKGROUND

Understanding genetic diversity is a core issue in conservation genetics. However, previous genetic diversity evaluations of narrowly distributed species have rarely used closely related widespread species as a reference. Furthermore, identifying natural hybridization signals between narrowly and widely distributed sympatric species is of great importance for the development of species conservation programs.

METHODS

In this study, population genotyping by sequencing (GBS) was performed for a narrowly distributed species, Geodorum eulophioides (endemic and endangered in Southwest China), and a widespread species, G. densiflorum. A total of 18,490 high-quality single nucleotide polymorphisms (SNPs) were identified at the whole-genome level.

RESULTS

The results showed that the nucleotide diversity and heterozygosity of G. eulophioides were significantly higher than those of G. densiflorum, confirming that narrowly distributed species can still preserve high genetic diversity. Consistent with taxonomic boundaries, all sampled individuals from the two species were divided into two genetic clusters and showed high genetic differentiation between species. However, in a sympatric population, a few G. eulophioides individuals were detected with genetic components from G. densiflorum, suggesting potential interspecific natural hybridization. This hypothesis was supported by Treemix analysis and hand-hybridization trials. Invasion of the habitat of G. eulophioides invasion by G. densiflorum under anthropogenic disturbance may be the main factor causing interspecific hybridization.

CONCLUSIONS

Therefore, reducing or avoiding habitat disturbance is a key measure to protect the G. eulophioides populations. This study provides valuable information for future conservation programs for narrowly distributed species.

Construction of a high-density genetic linkage map and QTL mapping for bioenergy-related traits in sweet sorghum [Sorghum bicolor (L.) Moench]

Sorghum is an important but arguably undervalued cereal crop, grown in large areas in Asia and Africa due to its natural resilience to drought and heat. There is growing demand for sweet sorghum as a source of bioethanol as well as food and feed. The improvement of bioenergy-related traits directly affects bioethanol production from sweet sorghum; therefore, understanding the genetic basis of these traits would enable new cultivars to be developed for bioenergy production. In order to reveal the genetic architecture behind bioenergy-related traits, we generated an F₂ population from a cross between sweet sorghum cv. 'Erdurmus' and grain sorghum cv. 'Ogretmenoglu'. This was used to construct a genetic map from SNPs discovered by double-digest restriction-site associated DNA sequencing (ddRAD-seq). F₃ lines derived from each F₂ individual were phenotyped for bioenergy-related traits in two different locations and their genotypes were analyzed with the SNPs to identify QTL regions. On chromosomes 1, 7, and 9, three major plant height (PH) QTLs (qPH1.1, qPH7.1, and qPH9.1) were identified, with phenotypic variation explained (PVE) ranging from 10.8 to 34.8%. One major QTL (qPJ6.1) on chromosome 6 was associated with the plant juice trait (PJ) and explained 35.2% of its phenotypic variation. For fresh biomass weight (FBW), four major QTLs (qFBW1.1, qFBW6.1, qFBW7.1, and qFBW9.1) were determined on chromosomes 1, 6, 7, and 9, which explained 12.3, 14.5, 10.6, and 11.9% of the phenotypic variation, respectively. Moreover, two minor QTLs (qBX3.1 and qBX7.1) of Brix (BX) were mapped on chromosomes 3 and 7, explaining 8.6 and 9.7% of the phenotypic variation, respectively. The QTLs in two clusters (qPH7.1/qBX7.1 and qPH7.1/qFBW7.1) overlapped for PH, FBW and BX. The QTL, qFBW6.1, has not been previously reported. In addition, eight SNPs were converted into cleaved amplified polymorphic sequences (CAPS) markers, which can be easily detected by agarose gel electrophoresis. These QTLs and molecular markers can be used for pyramiding and marker-assisted selection studies in sorghum, to develop advanced lines that include desirable bioenergy-related traits.

Genetic diversity and structure of Rhododendron meddianum, a plant species with extremely small populations

Rhododendron meddianum is a critically endangered species with important ornamental value and is also a plant species with extremely small populations. In this study, we used double digest restriction-site-associated DNA sequencing (ddRAD) technology to assess the genetic diversity, genetic structure and demographic history of the three extant populations of R. meddianum. Analysis of SNPs indicated that R. meddianum populations have a high genetic diversity (π = 0.0772 ± 0.0024, H _E  = 0.0742 ± 0.002). Both F _ST values (0.1582-0.2388) and AMOVA showed a moderate genetic differentiation among the R. meddianum populations. Meanwhile, STRUCTURE, PCoA and NJ trees indicated that the R. meddianum samples were clustered into three distinct genetic groups. Using the stairway plot, we found that R. meddianum underwent a population bottleneck about 70,000 years ago. Furthermore, demographic models of R. meddianum and its relative, Rhododendron cyanocarpum, revealed that these species diverged about 3.05 (2.21-5.03) million years ago. This divergence may have been caused by environmental changes that occurred after the late Pliocene, e.g., the Asian winter monsoon intensified, leading to a drier climate. Based on these findings, we recommend that R. meddianum be conserved through in situ, ex situ approaches and that its seeds be collected for germplasm.

Molecular characterization of endangered Iranian terrestrial orchids using ISSR markers and association with floral and tuber-related phenotypic traits

We aimed to study the genetic diversity, population structure, and phylogeny of Iranian orchids using inter-simple sequence repeat (ISSR) markers to find markers associated with phenotypic traits. Based on the phenotypic analysis, the inflorescence length and the flower number of studied accessions ranged from 3.92 to 27.13 cm and 5 to 50, respectively. On the other hand, the tuber length ranged from 1.80 to 9.35 cm. A total of 310 reproducible ISSR fragments with a size range of 150 to 3000 bp were amplified. ISSR primers provided an average polymorphism information content of 0.391, varied from 0.488 for UBC-876 to 0.351 for UBC-842. Os.J population showed the lowest genetic diversity (H = 0.057 and I = 0.075), while Oyst.JA population showed the highest genetic diversity (H = 0.114 and I = 0.158). At species level, the average coefficient of genetic differentiation (G ST) ranged from 0.265 for Orchis simia to 0.587 for Himantoglossum affine. Gene flow (Nm) varied from 1.38 (O. simia) to 0.756 (Anacamptis collina). The UPGMA genetic similarity dendrogram using Jaccard coefficients (r = 0.973) revealed six main clusters. Based on the Bayesian clustering method, the highest probability of the data was achieved when accessions were divided into eight groups. Floral and tuber-related phenotypic traits represented high correlations together, and they were associated with some ISSR bands based on the multiple association analysis. Altogether, ISSR markers proved to be useful for discrimination and clarification of the relationships among species and populations collected from geographically different locations. Furthermore, it could identify the polymorphism among accessions within each population and species.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s12298-020-00920-0) contains supplementary material, which is available to authorized users.

A High-Density SNP Genetic Map Construction Using ddRAD-Seq and Mapping of Capsule Shattering Trait in Sesame

The seed-bearing capsule of sesame shatters at harvest. This wildish trait makes the crop unsuitable for mechanized harvesting and also restricts its commercial potential by limiting the cultivation for countries that have no access to low-cost labor. Therefore, the underlying genetic basis of the capsule shattering trait is highly important in order to develop mechanization-ready varieties for sustainable sesame farming. In the present study, we generated a sesame F₂ population derived from a cross between a capsule shattering cultivar (Muganli-57) and a non-shattering mutant (PI 599446), which was used to construct a genetic map based on double-digest restriction-site-associated DNA sequencing. The resulting high-density genetic map contained 782 single-nucleotide polymorphisms (SNPs) and spanned a length of 697.3 cM, with an average marker interval of 0.89 cM. Based on the reference genome, the capsule shattering trait was mapped onto SNP marker S8_5062843 (78.9 cM) near the distal end of LG8 (chromosome 8). In order to reveal genes potentially controlling the shattering trait, the marker region (S8_5062843) was examined, and a candidate gene including six CDSs was identified. Annotation showed that the gene encodes a protein with 440 amino acids, sharing ∼99% homology with transcription repressor KAN1. Compared with the capsule shattering allele, the SNP change and altered splicing in the flanking region of S8_5062843 caused a frameshift mutation in the mRNA, resulting in the loss of function of this gene in the mutant parent and thus in non-shattering capsules and leaf curling. With the use of genomic data, InDel and CAPS markers were developed to differentiate shattering and non-shattering capsule genotypes in marker-assisted selection studies. The obtained results in the study can be beneficial in breeding programs to improve the shattering trait and enhance sesame productivity.

Genome-wide SNP discovery from Darjeeling tea cultivars - their functional impacts and application toward population structure and trait associations

Genotyping by sequencing and identification of functionally relevant nucleotide variations in crop accessions are the key steps to unravel genetic control of desirable traits. Elite cultivars of Darjeeling tea were undergone SNP genotyping by double-digest restriction-site associated DNA sequencing method. This study reports a set of 54,206 high-quality SNP markers discovered from ~10.4 GB sequence data, encompassing 15 chromosomes of the reference tea genome. Genetic relatedness among the accessions conforms to the analyses of Bayesian clustering, UPGMA, and PCoA methods. Genomic positions of the discovered SNPs and their putative effect on annotated genes designated a thoughtful understanding of their functional aspects in tea system biology. A group of 95 genes was identified to be affected by high impact variants. Genome-wide association analyses of 21 agronomic and biochemical phenotypes resulted in trait-linked polymorphic loci with strong confidence (p < 0.05 and 0.001).

Using genomic information for management planning of an endangered perennial, Viola uliginosa

Species occupying habitats subjected to frequent natural and/or anthropogenic changes are a challenge for conservation management. We studied one such species, Viola uliginosa, an endangered perennial wetland species typically inhabiting sporadically flooded meadows alongside rivers/lakes. In order to estimate genomic diversity, population structure, and history, we sampled five sites in Finland, three in Estonia, and one each in Slovenia, Belarus, and Poland using genomic SNP data with double-digest restriction site-associated DNA sequencing (ddRAD-seq). We found monophyletic populations, high levels of inbreeding (mean population F SNP = 0.407-0.945), low effective population sizes (N e = 0.8-50.9), indications of past demographic expansion, and rare long-distance dispersal. Our results are important in implementing conservation strategies for V. uliginosa, which should include founding of seed banks, ex situ cultivations, and reintroductions with individuals of proper origin, combined with continuous population monitoring and habitat management.

Genetic diversity and population structure of the Mediterranean sesame core collection with use of genome-wide SNPs developed by double digest RAD-Seq

The Mediterranean sesame core collection contains agro-morphologically superior sesame accessions from geographically diverse regions in four continents. In the present investigation, the genetic diversity and population structure of this collection was analyzed with 5292 high-quality SNPs discovered by double-digest restriction site associated DNA (ddRAD) sequencing, a cost-effective and flexible next-generation sequencing method. The genetic distance between pairs of accessions varied from 0.023 to 0.524. The gene diversity was higher in accessions from Asia than from America, Africa, and Europe. The highest genetic differentiation was observed between accessions collected from America and Europe. Structure analysis showed the presence of three subpopulations among the sesame accessions, and only six accessions were placed in an admixture group. Phylogenetic tree and principal coordinate analysis clustered the accessions based on their countries of origin. However, no clear division was evident among the sesame accessions with regard to their continental locations. This result was supported by an AMOVA analysis, which revealed a genetic variation among continental groups of 5.53% of the total variation. The large number of SNPs clearly indicated that the Mediterranean sesame core collection is a highly diverse genetic resource. The collection can be exploited by breeders to select appropriate accessions that will provide high genetic gain in sesame improvement programs. The high-quality SNP data generated here should also be used in genome-wide association studies to explore qualitative trait loci and SNPs related to economically and agronomically important traits in sesame.

Optimizing ddRADseq in Non-Model Species: A Case Study in Eucalyptus dunnii Maiden

Restriction site-associated DNA sequencing (RADseq) and its derived protocols, such as double digest RADseq (ddRADseq), offer a flexible and highly cost-effective strategy for efficient plant genome sampling. This has become one of the most popular genotyping approaches for breeding, conservation, and evolution studies in model and non-model plant species. However, universal protocols do not always adapt well to non-model species. Herein, this study reports the development of an optimized and detailed ddRADseq protocol in Eucalyptus dunnii, a non-model species, which combines different aspects of published methodologies. The initial protocol was established using only two samples by selecting the best combination of enzymes and through optimal size selection and simplifying lab procedures. Both single nucleotide polymorphisms (SNPs) and simple sequence repeats (SSRs) were determined with high accuracy after applying stringent bioinformatics settings and quality filters, with and without a reference genome. To scale it up to 24 samples, we added barcoded adapters. We also applied automatic size selection, and therefore obtained an optimal number of loci, the expected SNP locus density, and genome-wide distribution. Reliability and cross-sequencing platform compatibility were verified through dissimilarity coefficients of 0.05 between replicates. To our knowledge, this optimized ddRADseq protocol will allow users to go from the DNA sample to genotyping data in a highly accessible and reproducible way.

Genome-wide researches and applications on Dendrobium

This review summarizes current knowledge of chromosome characterization, genetic mapping, genomic sequencing, quality formation, floral transition, propagation, and identification in Dendrobium. The widely distributed Dendrobium has been studied for a long history, due to its important economic values in both medicine and ornamental. In recent years, some species of Dendrobium and other orchids had been reported on genomic sequences, using the next-generation sequencing technology. And the chloroplast genomes of many Dendrobium species were also revealed. The chromosomes of most Dendrobium species belong to mini-chromosomes, and showed 2n = 38. Only a few of genetic studies were reported in Dendrobium. After revealing of genomic sequences, the techniques of transcriptomics, proteomics and metabolomics could be employed on Dendrobium easily. Some other molecular biological techniques, such as gene cloning, gene editing, genetic transformation and molecular marker developing, had also been applied on the basic research of Dendrobium, successively. As medicinal plants, insights into the biosynthesis of some medicinal components were the most important. As ornamental plants, regulation of flower related characteristics was the most important. More, knowledge of growth and development, environmental interaction, evolutionary analysis, breeding of new cultivars, propagation, and identification of species and herbs were also required for commercial usage. All of these studies were improved using genomic sequences and related technologies. To answer some key scientific issues in Dendrobium, quality formation, flowering, self-incompatibility and seed germination would be the focus of future research. And genome related technologies and studies would be helpful.