From population genetics to population genomics of forest trees: integrated population genomics approach

Early works by Altukhov and his associates on pine and spruce laid the foundation for Russian population genetic studies on tree species with the use of molecular genetic markers. In recent years, these species have become especially popular as nontraditional eukaryotic models for population and evolutionary genomic research. Tree species with large, cross-pollinating native populations, high genetic and phenotypic variation, growing in diverse environments and affected by environmental changes during hundreds of years of their individual development, are an ideal model for studying the molecular genetic basis of adaptation. The great advance in this field is due to the rapid development of population genomics in the last few years. In the broad sense, population genomics is a novel, fast-developing discipline, combining traditional population genetic approaches with the genomic level of analysis. Thousands of genes with known function and sometimes known genomic localization can be simultaneously studied in many individuals. This opens new prospects for obtaining statistical estimates for a great number of genes and segregating elements. Mating system, gene exchange, reproductive population size, population disequilibrium, interaction among populations, and many other traditional problems of population genetics can be now studied using data on variation in many genes. Moreover, population genomic analysis allows one to distinguish factors that affect individual genes, alleles, or nucleotides (such as, for example, natural selection) from factors affecting the entire genome (e.g., demography). This paper presents a brief review of traditional methods of studying genetic variation in forest tree species and introduces a new, integrated population genomics approach. The main stages of the latter are : (1) selection of genes, which are tentatively involved in variation of adaptive traits, by means of a detailed examination of the regulation and the expression of individual genes and genotypes, with subsequent determination of their complete allelic composition by direct nucleotide sequencing; (2) examination of the phenotypic effects of individual alleles by, e.g., association mapping; and (3) determining the frequencies of the selected alleles in natural population for identification of the adaptive variation pattern in the heterogeneous environment. Through decoding the phenotypic effects of individual alleles and identification of adaptive variation patterns at the population level, population genomics in the future will serve as a very helpful, efficient, and economical tool, essential for developing a correct strategy for conserving and increasing forests and other commercially valuable plant and animal species.     

From population genetics to population genomics of forest trees: Integrated population genomics approach

Early works by Altukhov and his associates on pine and spruce laid the foundation for Russian population genetic studies on tree species with the use of molecular genetic markers. In recent years, these species have become especially popular as nontraditional eukaryotic models for population and evolutionary genome-wide research. Tree species with large, cross-pollinating native populations, high genetic and phenotypic variation, growing in diverse environments and affected by environmental changes during hundreds of years of their individual development, are an ideal model for studying the molecular genetic basis of adaptation. The great advance in this field is due to the rapid development of population genomics in the last few years. In the broad sense, population genomics is a novel, fast-developing discipline, combining traditional population genetic approaches with the genome-wide level of analysis. Thousands of genes with known function and sometimes known genome-wide localization can be simultaneously studied in many individuals. This opens new prospects for obtaining statistical estimates for a great number of genes and segregating elements. Mating system, gene exchange, reproductive population size, population disequilibrium, interaction among populations, and many other traditional problems of population genetics can be now studied using data on variation in many genes. Moreover, population genome-wide analysis allows one to distinguish factors that affect individual genes, allelles, or nucleotides (such as, for example, natural selection) from factors affecting the entire genome (e.g., demography). This paper presents a brief review of traditional methods of studying genetic variation in forest tree species and introduces a new, integrated population genomics approach. The main stages of the latter are: (1) selection of genes, which are tentatively involved in variation of adaptive traits, by means of a detailed examination of the regulation and the expression of individual genes and genotypes, with subsequent determination of their complete allelic composition by direct nucleotide sequencing; (2) examination of the phenotypic effects of individual alleles by, e.g., association mapping; and (3) determining the frequencies of the selected alleles in natural population for identification of the adaptive variation pattern in the heterogeneous environment. Through decoding the phenotypic effects of individual alleles and identification of adaptive variation patterns at the population level, population genomics in the future will serve as a very helpful, efficient, and economical tool, essential for developing a correct strategy for conserving and increasing forests and other commercially valuable plant and animal species.     

The population biology of coevolution

New populational approaches to the study of coevolution among species are confronting two major problems: the geographic scale at which coevolution proceeds, and the long-standing issue of how species may coevolve with more than one other species. By incorporating the ecological structure of life histories and populations into analyses of the coevolutionary process, these studies are indicating that coevolutionary change is much more ecologically dynamic than indicated by earlier work. Rather than simply a slow, stately process shaping species over long periods of time, parts of the coevolutionary process may proceed rapidly (sometimes observable in less than a decade), continually molding and remolding populations and communities locally and over broad geographic scales.     

Establishing the evolutionary compatibility of potential sources of colonizers for overfished stocks: a population genomics approach

Identifying fish stock structure is fundamental to pinpoint stocks that might contribute colonizers to overfished stocks. However, a stock's potential to contribute to rebuilding hinges on demographic connectivity, a challenging parameter to measure. With genomics as a new tool, fisheries managers can detect signatures of natural selection and thus identify fishing areas likely to contribute evolutionarily compatible colonizers to an overfished area (i.e. colonizers that are not at a fitness disadvantage in the overfished area and able to reproduce at optimal rates). Identifying evolutionarily compatible stocks would help narrow the focus on establishing demographic connectivity where it matters. Here, we genotype 4723 SNPs in 616 orange roughy (Hoplostethus atlanticus) across five fishing areas off the Tasmanian coast in Australia. We ask whether these areas form a single genetic unit, and test for signatures of local adaptation. Results from amova , structure , discriminant analysis of principal components, bayesass and isolation by distance suggest that sampled locations are subjected to geneflow amounts that are above what is needed to establish ‘drift connectivity’. However, it remains unclear whether there is a single panmictic population or several highly connected populations. Most importantly, we did not find any evidence of local adaptation, suggesting that the examined orange roughy stocks are evolutionarily compatible. The data have helped test an assumption of the orange roughy management programme and to formulate hypotheses regarding stock demographic connectivity. Overall, our results demonstrate the potential of genomics to inform fisheries management, even when evidence for stock structure is sparse.     

Genome divergence and diversification within a geographic mosaic of coevolution

Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome‐wide divergence was slight between ecotypes (F_ST = 0.011–0.035), we found evidence of relative genetic differentiation (as measured by F_ST) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome‐wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource‐based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000 years.     

Cytonuclear coevolution: the genomics of cooperation

Without mitochondria we would be in big trouble, and there would be a global biological energy crisis if it were not for chloroplasts. Fortunately, genomic evolution over the past two billion years has ensured that the functions of these key organelles are with us to stay. Whole-genome analyses have not only proven that mitochondria and chloroplasts are descended from formerly free-living bacteria, but have also shown that it is difficult to define eukaryotes without reference to the fusion and coevolution of host and endosymbiont genomes. Here, we review how the macro- and microevolutionary insights that follow from the genomics of cytonuclear interactions are uniting molecular evolution, structural proteomics, population genetics and problems in aging and disease. Our goals are to clarify the coevolutionary events that have governed nuclear and organelle evolution, and to encourage further critical analyses of these interactions as problems in the study of co-adapted gene complexes.     

The influence of a competitor on the geographic mosaic of coevolution between crossbills and lodgepole pine

The geographic mosaic theory of coevolution posits that the form of selection between interacting species varies across a landscape with coevolution important and active in some locations (i.e., coevolutionary hotspots) but not in others (i.e., coevolutionary coldspots). We tested the hypothesis that the presence of red squirrels (Tamiasciurus hudsonicus) affects the occurrence of coevolution between red crossbills (Loxia curvirostra complex) and Rocky Mountain lodgepole pine (Pinus contorta ssp. latifolia) and thereby provides a mechanism giving rise to a geographic mosaic of selection. Red squirrels are the predominant predispersal seed predator and selective agent on lodgepole pine cones. However, in four isolated mountain ranges east and west of the Rocky Mountains, red squirrels are absent and red crossbills are the main predispersal seed predator. These isolated populations of pine have apparently evolved without Tamiasciurus for about 10,000 to 12,000 years. Based on published morphological, genetic, and paleobotanical studies, we infer that cone traits in these isolated populations that show parallel differences from cones in the Rocky Mountains have changed in parallel. We used data on crossbill and conifer cone morphology and feeding preferences and efficiency to detect whether red crossbills and lodgepole pine exhibit reciprocal adaptations, which would imply coevolution. Cone traits that act to deter Tamiasciurus and result in high ratios of cone mass to seed mass were less developed in the isolated populations. Cone traits that act to deter crossbills include larger and thicker scales and perhaps increased overlap between successive scales and were enhanced in the isolated populations. In the larger, isolated mountain ranges crossbills have evolved deeper, shorter, and therefore more decurved bills to exploit these cones. This provides crossbills with higher feeding rates, and the change in bill shape has improved efficiency by reducing the concomitant increases in body mass and daily energy expenditures that would have resulted if only bill size had increased. These parallel adaptations and counter adaptations in red crossbills and lodgepole pine are interpreted as reciprocal adaptations and imply that these crossbills and pine are in coevolutionary arms races where red squirrels are absent (i.e., coevolutionary hotspots) but not where red squirrels are present (i.e., coevolutionary cold-spots).     

The rise of yeast population genomics

Genome sequences of multiple individuals are essential to determine the forces shaping sequence variation as well as to understand the relationship between genotype and phenotype. Because of their wide ecological, geographical and genetic diversity, yeast species represent an ideal model system for population genomics. Recently, there has been a renewed interest in characterizing the genetic diversity within yeast species such as Saccharomyces cerevisiae and Saccharomyces paradoxus. Here, we review recent progress in the exploration of the intraspecific diversity using large collections of yeast isolates. These recent large-scale polymorphism surveys have increased our understanding of the population structures as well as the evolutionary history of the species. In addition, these resources represent a powerful framework for dissecting the relationship between genotype and phenotype.     

The population genomics of within-host Mycobacterium tuberculosis

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.Subject terms: Genetic variation, Population genetics     

The population genomics of hepatitis B virus

Hepatitis B virus (HBV) infection is considered as the fifth leading cause of death due to infectious diseases and has a worldwide prevalence. The particular geographical distribution of the eight previously defined genotypes of HBV suggests that the viral population is highly structured. The presence of such population structure is likely to affect the geographical distribution of polymorphisms involved in disease progression. In this study, we determined the structure of the HBV population using a clustering approach based on the observed allele frequencies at the polymorphic loci. We used all full-genome sequences publicly available and obtained a significant clustering of the HBV population into four main clusters, strongly associated with the current classification into genotypes. One of these main clusters could itself be split into three well-supported subclusters, highlighting the hierarchical nature of the population differentiation between HBV strains. The extremely clear-cut subdivision of the HBV population further indicates that recombination in HBV is not as extensive as previously assumed.     

Dos and don'ts of testing the geographic mosaic theory of coevolution

The geographic mosaic theory of coevolution is stimulating much new research on interspecific interactions. We provide a guide to the fundamental components of the theory, its processes and main predictions. Our primary objectives are to clarify misconceptions regarding the geographic mosaic theory of coevolution and to describe how empiricists can test the theory rigorously. In particular, we explain why confirming the three main predicted empirical patterns (spatial variation in traits mediating interactions among species, trait mismatching among interacting species and few species-level coevolved traits) does not provide unequivocal support for the theory. We suggest that strong empirical tests of the geographic mosaic theory of coevolution should focus on its underlying processes: coevolutionary hot and cold spots, selection mosaics and trait remixing. We describe these processes and discuss potential ways each can be tested.     

Comparative genomics reveals different population structures associated with host and geographic origin in antimicrobial-resistant Salmonella enterica

Genetic variation in a pathogen, including the causative agent of salmonellosis, Salmonella enterica, can occur as a result of eco-evolutionary forces triggered by dissimilarities of ecological niches. Here, we applied comparative genomics to study 90 antimicrobial resistant (AMR) S. enterica isolates from bovine and human hosts in New York and Washington states to understand host- and geographic-associated population structure. Results revealed distinct presence/absence profiles of functional genes and pseudogenes (e.g., virulence genes) associated with bovine and human isolates. Notably, bovine isolates contained significantly more transposase genes but fewer transposase pseudogenes than human isolates, suggesting the occurrence of large-scale transposition in genomes of bovine and human isolates at different times. The high correlation between transposase genes and AMR genes, as well as plasmid replicons, highlights the potential role of horizontally transferred transposons in promoting adaptation to antibiotics. By contrast, a number of potentially geographic-associated single-nucleotide polymorphisms (SNPs), rather than geographic-associated genes, were identified. Interestingly, 38% of these SNPs were in genes annotated as cell surface protein-encoding genes, including some essential for antibiotic resistance and host colonization. Overall, different evolutionary forces and limited recent inter-population transmission appear to shape AMR S. enterica population structure in different hosts and geographic origins.     

The demography and population genomics of evolutionary transitions to self-fertilization in plants

The evolution of self-fertilization from outcrossing has occurred on numerous occasions in flowering plants. This shift in mating system profoundly influences the morphology, ecology, genetics and evolution of selfing lineages. As a result, there has been sustained interest in understanding the mechanisms driving the evolution of selfing and its environmental context. Recently, patterns of molecular variation have been used to make inferences about the selective mechanisms associated with mating system transitions. However, these inferences can be complicated by the action of linked selection following the transition. Here, using multilocus simulations and comparative molecular data from related selfers and outcrossers, we demonstrate that there is little evidence for strong bottlenecks associated with initial transitions to selfing, and our simulation results cast doubt on whether it is possible to infer the role of bottlenecks associated with reproductive assurance in the evolution of selfing. They indicate that the effects of background selection on the loss of diversity and efficacy of selection occur rapidly following the shift to high selfing. Future comparative studies that integrate explicit ecological and genomic details are necessary for quantifying the independent and joint effects of selection and demography on transitions to selfing and the loss of genetic diversity.