Origin and demographic history of the endemic Taiwan spruce (Picea morrisonicola)

Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level

The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.

Assisted migration and the rare endemic plant species: the case of two endangered Mexican spruces

Background

In the projected climate change scenarios, assisted migration might play an important role in the ex situ conservation of the threatened plant species, by translocate them to similar suitable habitats outside their native distributions. However, it is unclear if such habitats will be available for the Rare Endemic Plant Species (REPS), because of their very restricted habitats. The aims of this study were to perform a population size assessment for the REPS Picea martinezii Patterson and Picea mexicana Martínez, and to evaluate the potential species distributions and their possibilities for assisted migration inside México and worldwide.

Methods

We performed demographic censuses, field surveys in search for new stands, and developed distribution models for Last Glacial Maximum (22,000 years ago), Middle Holocene (6,000 years ago), current (1961-1990) and future (2050 and 2070) periods, for the whole Mexican territory (considering climatic, soil, geologic and topographic variables) and for all global land areas (based only on climate).

Results

Our censuses showed populations of 89,266 and 39,059 individuals for P. martinezii and P. mexicana, respectively, including known populations and new stands. Projections for México indicated somewhat larger suitable areas in the past, now restricted to the known populations and new stands, where they will disappear by 2050 in a pessimistic climatic scenario, and scarce marginal areas (p = 0.5-0.79) remaining only for P. martinezii by 2070. Worldwide projections (based only on climate variables) revealed few marginal areas in 2050 only in México for P. martinezii, and several large areas (p ≥ 0.5) for P. mexicana around the world (all outside México), especially on the Himalayas in India and the Chungyang mountains in Taiwan with highly suitable (p ≥ 0.8) climate habitats in current and future (2050) conditions. However, those suitable areas are currently inhabited by other endemic spruces: Picea smithiana (Wall.) Boiss and Picea morrisonicola Hayata, respectively.

Conclusions

Assisted migration would only be an option for P. martinezii on scarce marginal sites in México, and the possibilities for P. mexicana would be continental and transcontinental translocations. This rises two possible issues for future ex situ conservation programs: the first is related to whether or not consider assisted migration to marginal sites which do not cover the main habitat requirements for the species; the second is related to which species (the local or the foreign) should be prioritized for conservation when suitable habitat is found elsewhere but is inhabited by other endemic species. This highlights the necessity to discuss new policies, guidelines and mechanisms of international cooperation to deal with the expected high species extinction rates, linked to projected climate change.

Complete plastome sequences of Picea asperata and P. crassifolia and comparative analyses with P. abies and P. morrisonicola

Picea asperata and P. crassifolia have sympatric ranges and are closely related, but the differences between these species at the plastome level are unknown. To better understand the patterns of variation among Picea plastomes, the complete plastomes of P. asperata and P. crassifolia were sequenced. Then, the plastomes were compared with the complete plastomes of P. abies and P. morrisonicola, which are closely and distantly related to the focal species, respectively. We also used these sequences to construct phylogenetic trees to determine the relationships among and between the four species as well as additional taxa from Pinaceae and other gymnosperms. Analysis of our sequencing data allowed us to identify 438 single nucleotide polymorphism (SNPs) point mutation events, 95 indel events, four inversion events, and seven highly variable regions, including six gene spacer regions (psbJ-petA, trnT-psaM, trnS-trnD, trnL-rps4, psaC-ccsA, and rps7-trnL) and one gene (ycf1). The highly variable regions are appropriate targets for future use in the phylogenetic reconstructions of closely related, sympatric species of Picea as well as Pinaceae in general.

Genetic structure of Sakhalin spruce (Picea glehnii) in northern Japan and adjacent regions revealed by nuclear microsatellites and mitochondrial gene sequences

The genetic structure of Sakhalin spruce (Picea glehnii) was studied across the natural range of the species, including two small isolated populations in south Sakhalin and Hayachine, by using six microsatellite loci and maternally inherited mitochondrial gene sequences. We also analyzed P. jezoensis, a sympatric spruce in the range. Genetic diversity of P. glehnii was higher in central Hokkaido and the lowest in the Hayachine. Bayesian clustering and principal coordinate analysis by using the microsatellites indicated that the Hayachine was clearly distinct from other populations, implying that it had undergone strong genetic drift since the last glacial period. P. glehnii harbored four mitochondrial haplotypes, two of which were shared with P. jezoensis. One of the two was observed without geographical concentration, suggesting its derivation from ancestral polymorphism. Another was observed in south Sakhalin and in P. jezoensis across Sakhalin. The Bayesian clustering--by using four microsatellite loci, including P. jezoensis populations--indicated unambiguous species delimitation, but with possible admixture of P. jezoensis genes into P. glehnii in south Sakhalin, where P. glehnii is abundantly overwhelmed by P. jezoensis; this might explain the occurrence of introgression of the haplotype of P. jezoensis into P. glehnii.

Amount of information needed for model choice in Approximate Bayesian Computation

Approximate Bayesian Computation (ABC) has become a popular technique in evolutionary genetics for elucidating population structure and history due to its flexibility. The statistical inference framework has benefited from significant progress in recent years. In population genetics, however, its outcome depends heavily on the amount of information in the dataset, whether that be the level of genetic variation or the number of samples and loci. Here we look at the power to reject a simple constant population size coalescent model in favor of a bottleneck model in datasets of varying quality. Not only is this power dependent on the number of samples and loci, but it also depends strongly on the level of nucleotide diversity in the observed dataset. Whilst overall model choice in an ABC setting is fairly powerful and quite conservative with regard to false positives, detecting weaker bottlenecks is problematic in smaller or less genetically diverse datasets and limits the inferences possible in non-model organism where the amount of information regarding the two models is often limited. Our results show it is important to consider these limitations when performing an ABC analysis and that studies should perform simulations based on the size and nature of the dataset in order to fully assess the power of the study.