Population Genetic Structure and Demographic History of Primula fasciculata in Southwest China

The Metasequoia genome and evolutionary relationships among redwoods

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

A comparative plastomic analysis of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow and implication of the origin of Chinese jujube

Comparative plastomics can be used to explicitly dissect various types of plastome variation. In the present study, the plastome variation pattern of Ziziphus jujuba var. spinosa (also called sour jujube) and its phylogenomic relationship with Chinese jujube were investigated. Plastomes of 21 sour jujube individuals were sequenced and assembled. The length of the sour jujube plastomes ranged between 159399 and 161279 bp. The plastomes exhibited collinearity of structure, gene order and content. The most divergent regions were located in the intergenic spacers, such as trnR-UCU-atpA and psbZ-trnG-UCC. Sliding window analysis demonstrated that the sequence variation among the sour jujube plastomes was relatively low. Sixty-two to 76 SSRs with 4 motif types were identified in the sour jujube plastomes with a predominant motif type of A/T. Three protein-coding genes exhibited higher nonsynonymous/synonymous substitution ratios, indicating that these genes may undergo positive selection. A total of 80 SNPs were detected and 1266 potential RNA editing sites of 23 protein-coding genes were predicted. In the phylogenomic tree constructed, sour jujube has a sister relationship to Chinese jujube, which indicates that Chinese jujube may have originated or been domesticated from sour jujube. The present study explicitly investigated the individual-level plastome variation of sour jujube and provides potential valuable molecular markers for future genetic-related study of this lineage.

Differentiated historical demography and ecological niche forming present distribution and genetic structure in coexisting two salamanders (Amphibia, Urodela, Hynobiidae) in a small island, Japan

Background

The climatic oscillations in the Quaternary period considerably shaped the distribution and population genetic structure of organisms. Studies on the historical dynamics of distribution and demography not only reflect the current geographic distribution but also allow us to understand the adaption and genetic differentiation of species. However, the process and factors affecting the present distribution and genetic structure of many taxa are still poorly understood, especially for endemic organisms to small islands.

Methods

Here, we integrated population genetic and ecological niche modelling approaches to investigate the historical distribution and demographic dynamics of two co-existing salamanders on Tsushima Island, Japan: the true H. tsuensis (Group A), and Hynobius sp. (Group B). We also examined the hypothesis on the equivalency and similarity of niches of these groups by identity and background tests for ecological niche space.

Results

Our result showed that Group A is considered to have undergone a recent population expansion after the Last Glacial Maximum while it is unlikely to have occurred in Group B. The highest suitability was predicted for Group A in southern Tsushima Island, whereas the northern part of Tsushima Island was the potential distribution of Group B. The results also suggested a restricted range of both salamanders during the Last Interglacial and Last Glacial Maximum, and recent expansion in Mid-Holocene. The genetic landscape-shape interpolation analysis and historical suitable area of ecological niche modelling were consistent, and suggested refugia used during glacial ages in southern part for Group A, and in northern part of Tsushima Island for Group B. Additionally, we found evidence of nonequivalence for the ecological niche of the two groups of the salamanders, although our test could not show either niche divergence or conservatism based on the background tests. The environmental predictors affecting the potential distribution of each group also showed distinctiveness, leading to differences in selecting suitable areas. Finally, the combination of population genetics and ecological modeling has revealed the differential demographic/historical response between coexisting two salamanders on a small island.

Chromosome-level de novo genome assembly and whole-genome resequencing of the threatened species Acanthochlamys bracteata (Velloziaceae) provide insights into alpine plant divergence in a biodiversity hotspot

The Hengduan Mountains region is an important hotspot of alpine plant diversity and endemism. Acanthochlamys bracteata is a species of a threatened monotypic genus endemic to the Hengduan Mountains. In this study, we present a high-quality, chromosome-level reference genome for A. bracteata, constructed using long reads, short reads and Hi-C technology. We characterized its genetic diversity, population structure, demographic history and gene flow by resequencing individuals collected across its distribution. Comparative genomics analyses based on sequence information from single-copy orthologous genes revealed that A. bracteata and Dioscorea rotundata diverged ~104.5 million years ago. Whole-genome resequencing based on population genetic analysis revealed that the division of the 14 populations into 10 distinct clusters reflected geographical divergence, and three separate high levels of gene flow occurred sequentially between isolated populations of the Hengduan Mountains, a finding which is consistent with the turnover between ice ages and interglacial periods. Our findings indicate that Quaternary climatic changes played an important role in shaping the genetic structure and demographic trajectories of A. bracteata, and provide critical insights into the genetic status and evolutionary history of this poorly understood species, and possibly other alpine plants with a similar distribution. This study demonstrates the usefulness of population genomics for evaluating the effects of past climatic changes and identifying conservation units for the conservation and management of threatened species. Our high-quality genome represents a valuable resource for future studies of the underlying molecular mechanisms of adaptive evolution and provides insight for further comparative genomic analysis with other Velloziaceae species.

Shennongjia-Wushan Mountains-One cryptic glacial refugium introduced by the phylogeographical study of the Geometridae moth Ourapteryx szechuana Wehrli

The origin and evolution of biodiversity in the Shennongjia and Wushan Mountains, located in central China, are little known. In this study, we used Ourapteryx szechuana, which is widely distributed in China and northern Nepal, to explore whether these mountains acted as glacial refugia during climate oscillations of the Quaternary. In total, 192 samples of O. szechuana were collected throughout much of the distribution range. Phylogenetic analysis, molecular dating, demographic history reconstructions, and MAXENT were used to investigate the evolutionary history and differentiation mechanisms and predict the potential species distributions during four different periods. The phylogenetic tree and the star-like median-joining network strongly supported two reciprocally monophyletic and allopatric lineages. Lineage I was restricted to the Shennongjia and Wushan Mountains. The divergence time of O. szechuana from its sister species O. thibetaria was approximately 1.94 Ma. The differentiation processes of the two intraspecific lineages occurred at approximately 0.47 Ma. The demographic history reconstruction and the ecological niche model suggested that Lineage II experienced an expansion after the LGM (Last Glacial Maximum), whereas Lineage I did not experience any expansion. Our results suggested the Naynayxungla glaciation promoted the divergence of the two lineages by restricting them to different refugia. The valleys of the Shennongjia-Wushan Mountains may have kept stable and warm (thus ice-free) environments during Quaternary glaciations, allowing this region to act as a glacial refugia. Our studies show that the Shennongjia and Wushan Mountains are likely to be important but little studied glacial refugia for the insect and thus worthy of more attention.

Genomic analysis of Medicago ruthenica provides insights into its tolerance to abiotic stress and demographic history

Medicago ruthenica has been recently cultivated as a new forage crop and has been recognized as a source of genes to improve abiotic stress tolerance in cultivated alfalfa because of its remarkable tolerance to drought, salinity-alkalinity, and cold and snowy winters. Here, we reveal a chromosome-scale genome sequence of M. ruthenica based on Illumina, PacBio, and Hi-C data. The assembled genome consists of 903.56 Mb with 50,268 annotated protein-coding genes, which is larger and contains relatively more genes than Medicago truncatula (420 Mb and 44,623 genes) and Medicago sativa spp. caerulea (793 Mb and 47,202 genes). All three species shared the ancestral Papilionoideae whole-genome duplication event before their divergence. The more recent expansion of repetitive elements compared to that in the other two species was determined to have contributed greatly to the larger genome size of M. ruthenica. We further found that multiple gene and transcription factor families (e.g., SOS homologous genes, NAC, C2H2, and CAMTA) have expanded in M. ruthenica, which might have led to its enhanced tolerance to abiotic stress. In addition, M. ruthenica harbors more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed two genetic lineages, reflecting the west and east of its geographical distribution, respectively. The two lineages probably diverged during the last glaciation and survived in multiple refugia at the last glacial maximum, followed by recent expansion. Our genomic data provide a genetic basis for further molecular breeding research on M. ruthenica and alfalfa.