Cytonuclear disequilibrium and genetic drift in a natural population of ponderosa pine

Integrative Pre-Breeding for Biotic Resistance in Forest Trees

Climate change is unleashing novel biotic antagonistic interactions for forest trees that may jeopardize populations' persistence. Therefore, this review article envisions highlighting major opportunities from ecological evolutionary genomics to assist the identification, conservation, and breeding of biotic resistance in forest tree species. Specifically, we first discuss how assessing the genomic architecture of biotic stress resistance enables us to recognize a more polygenic nature for a trait typically regarded Mendelian, an expectation from the Fisherian runaway pathogen-host concerted arms-race evolutionary model. Secondly, we outline innovative pipelines to capture and harness natural tree pre-adaptations to biotic stresses by merging tools from the ecology, phylo-geography, and omnigenetics fields within a predictive breeding platform. Promoting integrative ecological genomic studies promises a better understanding of antagonistic co-evolutionary interactions, as well as more efficient breeding utilization of resistant phenotypes.

Ancient events and climate adaptive capacity shaped distinct chloroplast genetic structure in the oak lineages

BACKGROUND

Understanding the origin of genetic variation is the key to predict how species will respond to future climate change. The genus Quercus is a species-rich and ecologically diverse woody genus that dominates a wide range of forests and woodland communities of the Northern Hemisphere. Quercus thus offers a unique opportunity to investigate how adaptation to environmental changes has shaped the spatial genetic structure of closely related lineages. Furthermore, Quercus provides a deep insight into how tree species will respond to future climate change. This study investigated whether closely related Quercus lineages have similar spatial genetic structures and moreover, what roles have their geographic distribution, ecological tolerance, and historical environmental changes played in the similar or distinct genetic structures.

RESULTS

Despite their close relationships, the three main oak lineages (Quercus sections Cyclobalanopsis, Ilex, and Quercus) have different spatial genetic patterns and occupy different climatic niches. The lowest level and most homogeneous pattern of genetic diversity was found in section Cyclobalanopsis, which is restricted to warm and humid climates. The highest genetic diversity and strongest geographic genetic structure were found in section Ilex, which is due to their long-term isolation and strong local adaptation. The widespread section Quercus is distributed across the most heterogeneous range of environments; however, it exhibited moderate haplotype diversity. This is likely due to regional extinction during Quaternary climatic fluctuation in Europe and North America.

CONCLUSIONS

Genetic variations of sections Ilex and Quercus were significantly predicted by geographic and climate variations, while those of section Cyclobalanopsis were poorly predictable by geographic or climatic diversity. Apart from the different historical environmental changes experienced by different sections, variation of their ecological or climatic tolerances and physiological traits induced varying responses to similar environment changes, resulting in distinct spatial genetic patterns.

Heteroplasmy and Patterns of Cytonuclear Linkage Disequilibrium in Wild Carrot

Organellar genomes are considered to be strictly uniparentally-inherited. Uniparental inheritance allows for cytonuclear coevolution and the development of highly coordinated cytonuclear interactions. Yet, instances of biparental inheritance have been documented across eukaryotes. Biparental inheritance in otherwise uniparentally-inherited organelles is termed leakage (maternal or paternal) and allows for the presence of multiple variants of the same organellar genome within an individual, called heteroplasmy. It is unclear what, if any, evolutionary consequences are placed on nuclear and/or organellar genomes due to heteroplasmy. One way of accessing cytonuclear interactions and potential coevolution is through calculating cytonuclear linkage disequilibrium (cnLD), or the non-random association of alleles between nuclear and organellar genomes. Patterns of cnLD can indicate positive or negative cytonuclear selection, coevolution between the nuclear and organellar genomes, non-traditional organellar inheritance, or instances of ancestral heteroplasmy. In plants, cytonuclear interactions have been shown to play a role in cytoplasmic male sterility which occurs in gynodioecious species and is associated with leakage. We used the gynodioecious species, Daucus carota L. spp. carota, or wild carrot, to investigate cnLD. We genotyped a total of 265 individuals from two regions of the USA at 15 nuclear microsatellites, the mitochondrial genes cox1 and atp9, and an intergenic region between trnS and trnG (StoG) in the plastid genome to calculate nuclear–nuclear LD (nucLD), cnLD, and organellar LD (i.e., within the mtDNA and between mtDNA and ptDNA) within the two regions. We were further able to identify cox1 and StoG heteroplasmy and calculate some of the same LD measures within heteroplasmic and homoplasmic (non-heteroplasmic) datasets. We used a Z-transformation test to demonstrate that heteroplasmic individuals display significantly higher levels of cnLD within both regions. In spite of this, within and between organellar LD is low to moderate. Given these patterns of LD in two regions of the USA in which gene flow has been shown to occur between crop and wild carrot, we suggest that heteroplasmy is an evolutionary mechanism which permits the maintenance of cnLD while also acting to disrupt organellar LD.

Genetic Diversity among Age Classes of a Pinus sylvestris (L.) Population from the Białowieża Primeval Forest, Poland

The Białowieża Primeval Forest (BPF) is Europe’s last primeval forest and an irreplaceable area for biodiversity conservation due to its size, protection status, and substantially undisturbed nature. There is no other forest in Europe with such a large surface representing highly-advanced natural succession. This article reports on the first analysis of the genetic variability and demographic structure of a self-renewed Pinus sylvestris population located in BPF, using both chloroplast and mitochondrial DNA markers. The analysis of molecular variance (AMOVA) for chloroplast simple sequence repeats (cpSSRs) revealed a significant genetic differentiation among age classes that accounted for about 2% of the total variance, comparable to those reported among different populations of Scots pine. None of the 117 detected chloroplast haplotypes were common to all age classes. Haplotype diversity ranged from 0.370 to 0.415 for cpSSRs and from 0.320 to 0.455 for mitochondrial markers. The genetic variation of the studied age classes—represented by mitochondrial markers—strongly depicts the maternal genetic structure, indicating limited seed dispersal. Temporal genetic substructuring is maintained within a self-renewed population of Scots pine from the BPF.

Spatial genetic structure in wild cardoon, the ancestor of cultivated globe artichoke: Limited gene flow, fragmentation and population history

Nuclear and chloroplast markers and phenotypic characters were integrated to analyse the population genetic structure of wild cardoon, Cynara cardunculus var. sylvestris, the ancestor of cultivated globe artichoke, Cynara cardunculus var. scolymus on the island of Sardinia, Italy. The spatial scale ranged from a few metres to ∼200km. Wild cardoon appears to be genetically fragmented, with significant genetic divergence at various scales, indicating that gene flow is insufficient to counterbalance the effects of genetic drift or founder effects. Divergence between populations was higher for chloroplast (40%) than for nuclear markers (15%), suggesting that gene flow via seed was lower than via pollen. Two main genetic groups were detected; these correlated with differences in flowering time, capitula size, glossiness, and anthocyanin pigmentation. A complex population structure of wild cardoon emerged over small spatial scales, likely resulting from the interplay between gene dispersal, colonisation history and selective forces. Indeed, Sardinia appears to be a 'hybrid zone' of different gene pools. The island has unique diverse germplasm that has originated from hybridisation among different gene pools. The sampling of seeds from a few plants but from many sites is suggested as the best strategy to harvest the genetic diversity of wild cardoon.

Genetic evidence for multiple sources of the non-native fish Cichlasoma urophthalmus (Günther; Mayan Cichlids) in southern Florida

The number and diversity of source populations may influence the genetic diversity of newly introduced populations and affect the likelihood of their establishment and spread. We used the cytochrome b mitochondrial gene and nuclear microsatellite loci to identify the sources of a successful invader in southern Florida, USA, Cichlasoma urophthalmus (Mayan cichlid). Our cytochrome b data supported an introduction from Guatemala, while our microsatellite data suggested movement of Mayan Cichlids from the upper Yucatán Peninsula to Guatemala and introductions from Guatemala and Belize to Florida. The mismatch between mitochondrial and nuclear genomes suggests admixture of a female lineage from Guatemala, where all individuals were fixed for the mitochondrial haplotype found in the introduced population, and a more diverse but also relatively small number of individuals from Belize. The Florida cytochrome b haplotype appears to be absent from Belize (0 out of 136 fish screened from Belize had this haplotype). Genetic structure within the Florida population was minimal, indicating a panmictic population, while Mexican and Central American samples displayed more genetic subdivision. Individuals from the Upper Yucatán Peninsula and the Petén region of Guatemala were more genetically similar to each other than to fish from nearby sites and movement of Mayan Cichlids between these regions occurred thousands of generations ago, suggestive of pre-Columbian human transportation of Mayan Cichlids through this region. Mayan Cichlids present a rare example of cytonuclear disequilibrium and reduced genetic diversity in the introduced population that persists more than 30 years (at least 7-8 generations) after introduction. We suggest that hybridization occurred in ornamental fish farms in Florida and may contribute their establishment in the novel habitat. Hybridization prior to release may contribute to other successful invasions.

Patterns of cyto-nuclear linkage disequilibrium in Silene latifolia: genomic heterogeneity and temporal stability

Non-random association of alleles in the nucleus and cytoplasmic organelles, or cyto-nuclear linkage disequilibrium (LD), is both an important component of a number of evolutionary processes and a statistical indicator of others. The evolutionary significance of cyto-nuclear LD will depend on both its magnitude and how stable those associations are through time. Here, we use a longitudinal population genetic data set to explore the magnitude and temporal dynamics of cyto-nuclear disequilibria through time. We genotyped 135 and 170 individuals from 16 and 17 patches of the plant species Silene latifolia in Southwestern VA, sampled in 1993 and 2008, respectively. Individuals were genotyped at 14 highly polymorphic microsatellite markers and a single-nucleotide polymorphism (SNP) in the mitochondrial gene, atp1. Normalized LD (D') between nuclear and cytoplasmic loci varied considerably depending on which nuclear locus was considered (ranging from 0.005-0.632). Four of the 14 cyto-nuclear associations showed a statistically significant shift over approximately seven generations. However, the overall magnitude of this disequilibrium was largely stable over time. The observed origin and stability of cyto-nuclear LD is most likely caused by the slow admixture between anciently diverged lineages within the species' newly invaded range, and the local spatial structure and metapopulation dynamics that are known to structure genetic variation in this system.

SNP detection from de novo transcriptome sequencing in the bivalve Macoma balthica: marker development for evolutionary studies

Hybrid zones are noteworthy systems for the study of environmental adaptation to fast-changing environments, as they constitute reservoirs of polymorphism and are key to the maintenance of biodiversity. They can move in relation to climate fluctuations, as temperature can affect both selection and migration, or remain trapped by environmental and physical barriers. There is therefore a very strong incentive to study the dynamics of hybrid zones subjected to climate variations. The infaunal bivalve Macoma balthica emerges as a noteworthy model species, as divergent lineages hybridize, and its native NE Atlantic range is currently contracting to the North. To investigate the dynamics and functioning of hybrid zones in M. balthica, we developed new molecular markers by sequencing the collective transcriptome of 30 individuals. Ten individuals were pooled for each of the three populations sampled at the margins of two hybrid zones. A single 454 run generated 277 Mb from which 17K SNPs were detected. SNP density averaged 1 polymorphic site every 14 to 19 bases, for mitochondrial and nuclear loci, respectively. An scan detected high genetic divergence among several hundred SNPs, some of them involved in energetic metabolism, cellular respiration and physiological stress. The high population differentiation, recorded for nuclear-encoded ATP synthase and NADH dehydrogenase as well as most mitochondrial loci, suggests cytonuclear genetic incompatibilities. Results from this study will help pave the way to a high-resolution study of hybrid zone dynamics in M. balthica, and the relative importance of endogenous and exogenous barriers to gene flow in this system.

Sympatric bromeliad species (Pitcairnia spp.) facilitate tests of mechanisms involved in species cohesion and reproductive isolation in Neotropical inselbergs

The roles of intra- and interspecific gene flow in speciation and species evolution are topics of great current interest in molecular ecology and evolutionary biology. Recent modelling studies call for new empirical data to test hypotheses arising from the recent shift from a 'whole-genome reproductive isolation' view to a 'genic' view of species and speciation. Particularly scarce (and thus of particular interest) are molecular genetic data on recently radiated, naturally hybridizing species in strongly structured and species-rich environments. Here, we studied four sympatric plant species (Pitcairnia spp.; Bromeliaceae) adapted to Neotropical inselbergs (isolated outcrops resembling habitat 'islands' in tropical rainforests) using nuclear and plastid DNA. Patterns of plastid DNA haplotype sharing and nuclear genomic admixture suggest the presence of both, incomplete lineage sorting and interspecific gene flow over extended periods of time. Integrity and cohesion of inselberg species of Pitcairnia are maintained despite introgression and in the face of extremely low within-species migration rates (N(e)m < 1 migrant per generation). Cross-evaluation of our genetic data against published pollination experiments indicate that species integrity is maintained by the simultaneous action of multiple prezygotic barriers, including flowering phenology, pollinator isolation and divergent mating systems. Postzygotic Bateson-Dobzhansky-Muller incompatibilities appear to contribute to isolation, as suggested by asymmetric introgression rates of single loci. Our results suggest that incomplete lineage sorting, hybridization and introgression form integral aspects of adaptive radiation in Neotropical inselberg 'archipelagos'. Inselbergs with multiple closely related co-occurring species should be of special interest to students of speciation in mountain systems, and to ongoing conservation programmes in the Atlantic Rainforest biodiversity hotspot.

F(ST) in the cytonuclear system

Selection on nuclear (or organelle) sites inevitably affects the spatial distribution of a neutral organelle (or nuclear) allele via transient cytonuclear disequilibrium. Here I examine this effect in terms of F(st) for a neutral allele by bringing together cytonuclear genomes with contrasting modes of inheritance. The relationships between cytonuclear disequilibrium and increment in F(st) are explored and confirmed through Monte Carlo simulations. Results show that the transient increment in F(st) for a neutral allele is not only related to the vectors of seed and pollen dispersal but also to the mode of its inheritance. Such increments can be substantial under certain conditions. Seed dispersal is more effective than pollen dispersal in changing the transient increment. The cumulative effects from multiple selective nuclear sites can amplify the transient increment in F(st) for a neutral paternal or maternal organelle allele. Selection on selective organelle sites facilitates the transient increment in F(st) for a neutral nuclear allele. Partial selfing can significantly reinforce the transient increment in F(st). These theoretical insights highlight the roles of transient cytonuclear disequilibrium as a biological factor in evolving population differentiation and refine our practical interpretations of F(st) with cytonuclear markers.

Joint match probabilities for Y chromosomal and autosomal markers

Empirical tests of association between Y chromosome and autosomal markers are presented and a theoretical framework for determining a joint match probability is recommended. Statistical analyses of association were performed in 16 US populations between the autosomal genotypes from loci CSF1PO, FGA, THO1, TPOX, vWA, D3S1358, D5S818, D7S820, D8S1179, D13S317, D16S539, D18S512, D21S11 and Y chromosome haplotypes from loci DYS19, DYS385ab, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS438, and DYS439. The sample populations include individuals of European-, African-, Hispanic-, Native-, and Asian-American ancestry. The results are consistent with independence of Y and autosomal markers, although small amounts of dependence would likely have escaped our tests. Given the data in hand, we suggest it is appropriate to compute joint match probabilities by multiplying the Y haplotype frequency with the appropriately corrected autosomal frequency. In addition to correcting for autosomal frequency differences between groups, a further correction may be required. Since two individuals sharing the same Y haplotype are likely to be more recently related than two randomly chosen individuals, the autosomal frequencies have to be adjusted to account for this, akin to the theta correction used to account for population substructure. The structure imposed on the autosomal frequencies conditioned in a Y match is a function of the number of markers scored and their mutation rate. However, in most settings theta<0.01. When population structure is already present in the autosomes, the additional effect due to conditioning on the Y is small. For example, if the amount of structure in the population is theta=0.01 or 0.03 (the NRCII range), then the effect of conditioning on the Y results in only a trivial increase in theta to 0.02-0.04, respectively.

Importance of genetic drift during Pleistocene divergence as revealed by analyses of genomic variation

Determining what factors affect the structuring of genetic variation is key to deciphering the relative roles of different evolutionary processes in species differentiation. Such information is especially critical to understanding how the frequent shifts and fragmentation of species distributions during the Pleistocene translates into species differences, and why the effect of such rapid climate change on patterns of species diversity varies among taxa. Studies of mitochondrial DNA (mtDNA) have detected significant population structure in many species, including those directly impacted by the glacial cycles. Yet, understanding the ultimate consequence of such structure, as it relates to how species divergence occurs, requires demonstration that such patterns are also shared with genomic patterns of differentiation. Here we present analyses of amplified fragment length polymorphisms (AFLPs) in the montane grasshopper Melanoplus oregonensis to assess the evolutionary significance of past demographic events and associated drift-induced divergence as inferred from mtDNA. As an inhabitant of the sky islands of the northern Rocky Mountains, this species was subject to repeated and frequent shifts in species distribution in response to the many glacial cycles. Nevertheless, significant genetic structuring of M. oregonensis is evident at two different geographic and temporal scales: recent divergence associated with the recolonization of the montane meadows in individual sky islands, as well as older divergence associated with displacements into regional glacial refugia. The genomic analyses indicate that drift-induced divergence, despite the lack of long-standing geographic barriers, has significantly contributed to species divergence during the Pleistocene. Moreover, the finding that divergence associated with past demographic events involves the repartitioning of ancestral variation without significant reductions of genomic diversity has intriguing implications - namely, the further amplification of drift-induced divergence by selection.

Distortion of symmetrical introgression in a hybrid zone: evidence for locus-specific selection and uni-directional range expansion

The fate of species integrity upon natural hybridization depends on the interaction between selection and dispersal. The relative significance of these processes may be studied in the initial phase of contact before selection and gene flow reach equilibrium. Here we study a hybrid zone of two salamander species, Lyciasalamandra antalyana and Lyciasalamandra billae, at the initial phase of hybridization. We quantify the degree and mode of introgression using nuclear and mtDNA markers. The hybrid zone can be characterized as an abrupt transition zone, the central hybrid zone being only c. 400 m, but introgressed genes were traced up to 3 km. Introgression was traced in both sexes but gene flow may be slightly male-biased. Indirect evidence suggests that hybrid males are less viable than females. Introgression occurred at two levels: (1) locus-specific selection led to different allelic introgression patterns independent of species, while (2) asymmetrical species-level introgression occurred predominately from L. antalyana to L. billae due to range expansion of the former. This indicates that foreign genes can be incorporated into novel genomic environments, which in turn may contribute to the great diversity of morphological variants in Lyciasalamandra.

Species delimitation in native South American fire ants

The taxonomy of fire ants has been plagued by difficulties in recognizing species on the basis of morphological characters. We surveyed allozyme markers and sequences of the mtDNA COI gene in several closely related nominal species from two areas of sympatry in the native ranges to learn whether the morphology-based delimitation of these species is supported by genetic data. We found that Solenopsis invicta and Solenopsis richteri, pest species whose distinctiveness has been debated, appear to be fully reproductively isolated at both study sites. This isolation contrasts with the extensive hybridization occurring between them in the USA, where both have been introduced. We also found strong genetic differentiation consistent with barriers to gene flow between Solenopsis quinquecuspis and the other two species. However, several lines of evidence suggest that nuclear and mitochondrial genes of S. invicta and S. richteri are introgressing into S. quinquecuspis. The latter apparently is a recently derived member of the clade that includes all three species, suggesting that there has been insufficient time for its full development of intrinsic isolating mechanisms. Finally, our discovery of genetically distinct populations within both S. invicta and S. richteri suggests the presence of previously unrecognized (cryptic) species. Their existence, together with the difficulties in developing diagnostic morphological characters for described species, imply that the group is actively radiating species and that morphological divergence generally does not keep pace with the development of reproductive isolation and neutral genetic divergence in this process.

The evolutionary processes of mitochondrial and chloroplast genomes differ from those of nuclear genomes

This paper first introduces our present knowledge of the origin of mitochondria and chloroplasts, and the organization and inheritance patterns of their genomes, and then carries on to review the evolutionary processes influencing mitochondrial and chloroplast genomes. The differences in evolutionary phenomena between the nuclear and cytoplasmic genomes are highlighted. It is emphasized that varying inheritance patterns and copy numbers among different types of genomes, and the potential advantage achieved through the transfer of many cytoplasmic genes to the nucleus, have important implications for the evolution of nuclear, mitochondrial and chloroplast genomes. Cytoplasmic genes transferred to the nucleus have joined the more strictly controlled genetic system of the nuclear genome, including also sexual recombination, while genes retained within the cytoplasmic organelles can be involved in selection and drift processes both within and among individuals. Within-individual processes can be either intra- or intercellular. In the case of heteroplasmy, which is attributed to mutations or biparental inheritance, within-individual selection on cytoplasmic DNA may provide a mechanism by which the organism can adapt rapidly. The inheritance of cytoplasmic genomes is not universally maternal. The presence of a range of inheritance patterns indicates that different strategies have been adopted by different organisms. On the other hand, the variability occasionally observed in the inheritance mechanisms of cytoplasmic genomes reduces heritability and increases environmental components in phenotypic features and, consequently, decreases the potential for adaptive evolution.

Estimating seed vs. pollen dispersal from spatial genetic structure in the common ash

Spatial genetic structure was analysed with five highly polymorphic microsatellite loci in a Romanian population of common ash (Fraxinus excelsior L.), a wind-pollinated and wind-dispersed tree species occurring in mixed deciduous forests over almost all of Europe. Contributions of seed and pollen dispersal to total gene flow were investigated by analysing the pattern of decrease in kinship coefficients among pairs of individuals with geographical distance and comparing it with simulation results. Plots of kinship against the logarithm of distance were decomposed into a slope and a shape component. Simulations showed that the slope is informative about the global level of gene flow, in agreement with theoretical expectations, whereas the shape component was correlated with the relative importance of seed vs. pollen dispersal. Hence, our results indicate that insights into the relative contributions of seed and pollen dispersal to overall gene flow can be gained from details of the pattern of spatial genetic structure at biparentally inherited loci. In common ash, the slope provided an estimate of total gene dispersal in terms of Wright's neighbourhood size of Nb = 519 individuals. No precise estimate of seed vs. pollen flow could be obtained from the shape because of the stochasticity inherent to the data, but the parameter combinations that best fitted the data indicated restricted seed flow, sigmas pound 14 m, and moderate pollen flow, 70 m pound sigmap pound 140 m.

Mitochondrial haplotype distribution, seed dispersal and patterns of postglacial expansion of ponderosa pine

Maternally inherited mtDNA in a secondary contact zone of ponderosa pine revealed a cline less than 10 km wide - much narrower than previously described. A survey of 76 populations gave no evidence either of intermixing or of a mosaic contact zone. Such sharp contact zones are consistent with diffusive range expansion, rather than long distance colonization. However, evidence for long distance seed dispersal events was found in two populations where haplotypes were observed far from their main area of occurrence. The results suggest a small number of long distance colonists with diffusive dispersal from these centres.

Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure

We describe a method for comparing nuclear and organelle population differentiation (F(ST)) in seed plants to test the hypothesis that pollen and seed gene flow rates are equal. Wright's infinite island model is used, with arbitrary levels of self-fertilization and biparental organelle inheritance. The comparison can also be applied to gene flow in animals. Since effective population sizes are smaller for organelle genomes than for nuclear genomes and organelles are often uniparentally inherited, organelle F(ST) is expected to be higher at equilibrium than nuclear F(ST) even if pollen and seed gene flow rates are equal. To reject the null hypothesis of equal seed and pollen gene flow rates, nuclear and organelle F(ST)'s must differ significantly from their expected values under this hypothesis. Finite island model simulations indicate that infinite island model expectations are not greatly biased by finite numbers of populations (>/=100 subpopulations). The power to distinguish dissimilar rates of pollen and seed gene flow depends on confidence intervals for fixation index estimates, which shrink as more subpopulations and loci are sampled. Using data from the tropical tree Corythophora alta, we rejected the null hypothesis that seed and pollen gene flow rates are equal but cannot reject the alternative hypothesis that pollen gene flow is 200 times greater than seed gene flow.