Relative strength of fine-scale spatial genetic structure in paternally vs biparentally inherited DNA in a dioecious plant depends on both sex proportions and pollen-to-seed dispersal ratio

Paternal Inheritance of Mitochondrial DNA May Lead to Dioecy in Conifers

In angiosperms cytoplasmic DNA is typically passed on maternally through ovules. Genes in the mtDNA may cause male sterility. When male-sterile (female) cytotypes produce more seeds than cosexuals, they pass on more copies of their mtDNA and will co-occur with cosexuals with a neutral cytotype. Cytoplasmic gynodioecy is a well-known phenomenon in angiosperms, both in wild and crop plants. In some conifer families (e.g. Pinaceae) mitochondria are also maternally inherited. However in some other families (e.g. Taxaceae and Cupressaceae) mtDNA is paternally inherited through the pollen. With paternal mtDNA inheritance, male cytotypes that produce more pollen than cosexuals are expected to co-occur with cosexuals. This is uncharted territory. An ESS model shows that the presence of male cytotypes selects for more female allocation in the cosexual, i.e. for sexual specialisation. An allele that switches sex from male to female can then invade. This leads to rapid loss of the neutral cytotype of the cosexual, fixation of the male cytotype and dioecy with 50% males and 50% females. The models suggest that paternal inheritance of mtDNA facilitates the evolution dioecy. Consistent with this hypothesis the Pinaceae are 100% monoecious, while dioecy is common in the Taxaceae family and in the genus Juniperus (Cupressaceae). However, no reliable data are yet available on both mode of inheritance of mtDNA and gender variation of the same species. When cosexuals benefit from reproductive assurance (high selfing rate, low inbreeding depression, low fertilisation) they maintain themselves next to males and females. This predicted pattern with three sex types present in the same population is observed in conifers in nature.

ORTHOSKIM: in silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in non-model organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: 1) global sequence assembly, 2) mapping against reference sequences, and 3) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and non-coding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g. ultra-conserved elements) or collections of single-copy ortholog genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and 4 RNAseq libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 license and is available at: https://github.com/cpouchon/ORTHOSKIM.

Multiscale spatial genetic structure within and between populations of wild cherry trees in nuclear genotypes and chloroplast haplotypes

Spatial genetic structure (SGS) of plants mainly depends on the effective population size and gene dispersal. Maternally inherited loci are expected to have higher genetic differentiation between populations and more intensive SGS within populations than biparentally inherited loci because of smaller effective population sizes and fewer opportunities of gene dispersal in the maternally inherited loci. We investigated biparentally inherited nuclear genotypes and maternally inherited chloroplast haplotypes of microsatellites in 17 tree populations of three wild cherry species under different conditions of tree distribution and seed dispersal. As expected, interpopulation genetic differentiation was 6-9 times higher in chloroplast haplotypes than in nuclear genotypes. This difference indicated that pollen flow 4-7 times exceeded seed flow between populations. However, no difference between nuclear and chloroplast loci was detected in within-population SGS intensity due to their substantial variation among the populations. The SGS intensity tended to increase as trees became more aggregated, suggesting that tree aggregation biased pollen and seed dispersal distances toward shorter. The loss of effective seed dispersers, Asian black bears, did not affect the SGS intensity probably because of mitigation of the bear loss by other vertebrate dispersers and too few tree generations after the bear loss to alter SGS. The findings suggest that SGS is more variable in smaller spatial scales due to various ecological factors in local populations.

Recent Fragmentation May Not Alter Genetic Patterns in Endangered Long-Lived Species: Evidence From Taxus cuspidata

Forestland fragmentation caused by overexploitation of forest resources can in principle reduce genetic diversity, limit gene flow and eventually lead to species developing strong genetic structure. However, the genetic consequences of recent anthropogenic fragmentation of tree species remain unclear. Taxus cuspidata, which has extremely small populations distributed mainly in Changbai Mt. in Northeast (NE) China, has recently endured severe habitat fragmentation. Here, we investigate the pattern of genetic diversity and structure, identify risk factors, predict the future distribution and finally provide guidelines for the conservation and management of this species. We used three chloroplast and two mitochondrial DNA fragments, which are both paternally inherited in yews but differ in mutation rates, to genotype a total of 265 individuals from 26 populations covering the distribution of the species in China. Both chloroplast and mitochondrial data showed high degrees of genetic diversity, extensive gene flow over the entire geographical range and historical stability of both effective population size and distribution of the species. However, ecological niche modeling suggests a decrease in suitable areas for this species by the years 2050 and 2070. The maintenance of high genetic diversity and the existence of sufficient gene flow suggest that recent fragmentation has not affected the genetic composition of the long-lived tree T. cuspidata. However, severe impacts of anthropogenic activities are already threatening the species. Conservation and management strategies should be implemented in order to protect the remnant populations.

Seed and pollen gene dispersal in Taxus baccata, a dioecious conifer in the face of strong population fragmentation

Background and Aims

Dispersal is crucial due to its direct impact on dynamics of a species' distribution as well as having a role in shaping adaptive potential through gene flow. In plants forming scarce and small populations, knowledge about the dispersal process is required to assess the potential for colonizing new habitats and connectivity of present and future populations. This study aimed to assess dispersal potential in Taxus baccata, a dioecious gymnosperm tree with a wide but highly fragmented distribution.

Methods

Seed and pollen dispersal kernels were estimated directly in the framework of the spatially explicit mating model, where genealogies of naturally established seedlings were reconstructed with the help of microsatellite markers. In this way, six differently shaped dispersal functions were compared.

Key Results

Seed dispersal followed a leptokurtic distribution, with the Exponential-Power, the Power-law and Weibull being almost equally best-fitting models. The pollen dispersal kernel appeared to be more fat-tailed than the seed dispersal kernel, and the Lognormal and the Exponential-Power function showed the best fit. The rate of seed immigration from the background sources was not significantly different from the rate of pollen immigration (13.1 % vs. 19.7 %) and immigration rates were in agreement with or below maximum predictions based on the estimated dispersal kernels. Based on the multimodel approach, 95 % of seeds travel <109 m, while 95 % of pollen travels <704 m from the source.

Conclusions

The results showed that, at a local spatial scale, yew seeds travel shorter distances than pollen, facilitating a rapid development of a kinship structure. At the landscape level, however, although yew exhibits some potential to colonize new habitats through seed dispersal, genetic connectivity between different yew remnants is strongly limited. Taking into account strong population fragmentation, the study suggests that gene dispersal may be a limiting factor of the adaptability of the species.