Investigating the Genetic Diversity, Population Differentiation and Population Dynamics of Cycas segmentifida (Cycadaceae) Endemic to Southwest China by Multiple Molecular Markers

Identification and genetic diversity analysis of Rickettsia in Dermacentor nuttalli within inner Mongolia, China

Background

The genus Rickettsia contains the lineages spotted fever group (SFG), typhus group (TG), and transitional group (TRG). The spotted fever group Rickettsia (SFGR) is transmitted by ticks. The tick species Dermacentor nuttalli is considered the main vector carrying SFGR in Inner Mongolia. Studying the genetic diversity and population structure of Rickettsia is essential for developing effective control strategies and predicting evolutionary trends of Rickettsia.

Methods

In 2019 we collected 408 D. nuttalli in the Inner Mongolia Autonomous Region, detected the percentage of Rickettsia-positive specimens, and characterized the haplotypes. From the Rickettsia-positive ticks, the gltA and ompA genes were extracted, amplified, and sequenced.

Results

Ten haplotypes of the gltA gene and 22 haplotypes of the ompA gene were obtained. The phylogenetic analysis showed that the haplotypes G1–G7 and G9 of the gltA gene cluster with Rickettsia raoultii, while G8 and G10 cluster with Rickettsia sibirica. Haplotypes O1–O15, O18 and O20–O22 of the ompA gene cluster with R. raoultii, while O16 and O19 cluster with R. sibirica. The average haplotype diversity was 0.3 for gltA and 0.7 for ompA. The average nucleotide diversity was greater than 0.05. Neutrality tests were nonsignificant for Tajima’s D results and Fu’s Fs results. The fixation index values (FST) showed that the degree of genetic differentiation between most sampled populations was small (FST < 0.05), whereas some populations showed a medium (FST > 0.05) or large (FST > 0.15) degree of differentiation. Analysis of molecular variance (AMOVA) revealed that the variation within populations was greater than that between populations. The mismatch analysis of Rickettsia showed double peaks.

Conclusions

We found two Rickettsia spp. (R. raoultii and R. sibirica). The high genetic disparity of Rickettsia allows for easy adaption to different environments. Genetic differentiation between populations is small, and Rickettsia populations do not show a geographically differentiated structure. The high rates of retention and infection of Rickettsia in D. nuttalli together with the animal husbandry exchange in Inner Mongolia gradually led to the harmonization of genetic characteristics of Rickettsia across various regions. Overall, the significant genetic diversity and geographical structure of Rickettsia in D. nuttalli are critical for SFGR control.

Graphical Abstract

Genome-Wide Evidence for Complex Hybridization and Demographic History in a Group of Cycas From China

Cycads represent one of the most ancestral living seed plants as well as one of the most threatened plant groups in the world. South China is a major center and potential origin of Cycas, the most rapidly diversified lineage of cycads. However, genomic-wide diversity of Cycas remains poorly understood due to the challenge of generating genomic markers associated with their inherent large genomes. Here, we perform a comprehensive conservation genomic study based on restriction-site associated DNA sequencing (RADseq) data in six representative species of Cycas in South China. Consistently low genetic diversity and strong genetic differentiation were detected across species. Both phylogenetic inference and genetic structure analysis via several methods revealed generally congruent groups among the six Cycas species. The analysis with ADMIXTURE showed low mixing of genetic composition among species, while individuals of C. dolichophylla exhibited substantial genetic admixture with C. bifida, C. changjiangensis, and C. balansae. Furthermore, the results from Treemix, f₄-statistic, and ABBA-BABA test were generally consistent and revealed the complex patterns of interspecific gene flow. Relatively strong signals of hybridization were detected between C. dolichophylla and C. szechuanensis, and the ancestor of C. taiwaniana and C. changjiangensis. Distinct patterns of demographic history were inferred for these species by Stairway Plot, and our results suggested that both climate fluctuation and frequent geological activities during the late Pleistocene exerted deep impacts on the population dynamics of these species in South China. Finally, we explore the practical implications of our findings for the development of conservation strategies in Cycas. The present study demonstrates the efficiency of RADseq for conservation genomic studies on non-model species with large and complex genomes. Given the great significance of cycads as a radical transition in the evolution of plant biodiversity, our study provides important insights into the mechanisms of diversification in such recently radiated living fossil taxa.

Genetic characterization of the entire range of Cycas panzhihuaensis (Cycadaceae)

Cycas panzhihuaensis L. Zhou & S. Y. Yang (Cycadaceae) is an endangered gymnosperm species endemic to the dry-hot valley of the Jinsha River basin in southwest China. Although the wild C. panzhihuaensis population from Panzhihua Cycad Natural Reserve is well protected and its genetic diversity has been well assessed, the genetic characteristics of populations outside the nature reserve, which face larger risks of extinction, remain unknown. Furthermore, the population genetics and historical dynamics of this endemic and endangered species have not been examined across its entire range. In this study, to analyze the genetic diversity, phylogeographical structure and demographic history of C. panzhihuaensis from all its seven known locations, we sequenced and compared molecular data from chloroplastic DNA (psbA-trnH, psbM-trnD, and trnS-trnG), single-copy nuclear genes (PHYP, AC5, HSP70, and AAT) from 61 individuals, as well as 11 nuclear microsatellite loci (SSR) from 102 individuals. We found relatively high genetic diversity within populations and high genetic differentiation among populations of C. panzhihuaensis, which is consistent with the patterns of other Asian inland cycads. Although no significant phylogeographical structure was detected, we found that small and unprotected populations possess higher genetic diversity and more unique haplotypes, which revises our understanding of diversity within this species and deserves due attention. Analysis of demographic dynamics suggest that human activity might be the key threat to C. panzhihuaensis. Based on the genetic characterization of C. panzhihuaensis, we propose several practical guidelines for the conservation of this species, especially for the populations with small sizes.

Population Differentiation and Demographic History of the Cycas taiwaniana Complex (Cycadaceae) Endemic to South China as Indicated by DNA Sequences and Microsatellite Markers

Historical geology, climatic oscillations, and seed dispersal capabilities are thought to influence the population dynamics and genetics of plants, especially for distribution-restricted and threatened species. Investigating the genetic resources within and among taxa is a prerequisite for conservation management. The Cycas taiwaniana complex consists of six endangered species that are endemic to South China. In this study, we investigated the relationship between phylogeographic history and the genetic structure of the C. taiwaniana complex. To estimate the phylogeographic history of the complex, we assessed the genetic structure and divergence time, and performed phylogenetic and demographic historical analyses. Two chloroplast DNA intergenic regions (cpDNA), two single-copy nuclear genes (SCNGs), and six microsatellite loci (SSR) were sequenced for 18 populations. The SCNG data indicated a high genetic diversity within populations, a low genetic diversity among populations, and significant genetic differentiation among populations. Significant phylogeographical structure was detected. Structure and phylogenetic analyses both revealed that the 18 populations of the C. taiwaniana complex have two main lineages, which were estimated to diverge in the Middle Pleistocene. We propose that Cycas fairylakea was incorporated into Cycas szechuanensis and that the other populations, which are mainly located on Hainan Island, merged into one lineage. Bayesian skyline plot analyses revealed that the C. taiwaniana complex experienced a recent decline, suggesting that the complex probably experienced a bottleneck event. We infer that the genetic structure of the C. taiwaniana complex has been affected by Pleistocene climate shifts, sea-level oscillations, and human activities. In addition to providing new insights into the evolutionary legacy of the genus, the genetic characterizations will be useful for the conservation of Cycas species.

Recent Trends in Research on the Genetic Diversity of Plants: Implications for Conservation

Genetic diversity and its distribution, both within and between populations, may be determined by micro-evolutionary processes, such as the demographic history of populations, natural selection, and gene flow. In plants, indices of genetic diversity (e.g., k, h and π) and structure (e.g., FST) are typically inferred from sequences of chloroplast markers. Given the recent advances and popularization of molecular techniques for research in population genetics, phylogenetics, phylogeography, and ecology, we adopted a scientometric approach to compile evidence on the recent trends in the use of cpDNA sequences as markers for the analysis of genetic diversity in botanical studies, over the years. We also used phylogenetic modeling to assess the relative contribution of relatedness or ecological and reproductive characters to the genetic diversity of plants. We postulated that genetic diversity could be defined not only by microevolutionary factors and life history traits, but also by relatedness, so that species more closely related phylogenetically would have similar genetic diversities. We found a clear tendency for an increase in the number of studies over time, confirming the hypothesis that the advances in the area of molecular genetics have supported the accumulation of data on the genetic diversity of plants. However, we found that the vast majority of these data have been produced by Chinese authors, and refer specifically to populations of Chinese plants. Most of the data on genetic diversity have been obtained for species in the International Union for Conservation of Nature (IUCN) category NE (Not Evaluated), which indicates a relative lack of attention on threatened species. In general, we observed very high FST values in the groups analyzed and, as we focused primarily on species that have not been evaluated by the IUCN, the number of plant species that are threatened with extinction may be much greater than that indicated by the listing of this organization. We also found that the number of haplotypes (k) was influenced by the type of geographic distribution of the plant, while haplotype diversity (h) was affected by the type of flower, and the fixation index (FST), by the type of habitat. The plant species most closely-related phylogenetically have similar levels of genetic diversity. Overall, then, it will important to consider phylogenetic dependence in future studies that evaluate the effects of life-history traits on plant genetic diversity.