The molecular population genetics of regulatory genes

Regulatory loci, which may encode both trans acting proteins as well as cis acting promoter regions, are crucial components of an organism's genetic architecture. Although evolution of these regulatory loci is believed to underlie the evolution of numerous adaptive traits, there is little information on natural variation of these genes. Recent molecular population genetic studies, however, have provided insights into the extent of natural variation at regulatory genes, the evolutionary forces that shape them and the phenotypic effects of molecular regulatory variants. These recent analyses suggest that it may be possible to study the molecular evolutionary ecology of regulatory diversification by examining both the extent and patterning of regulatory gene diversity, the phenotypic effects of molecular variation at these loci and their ecological consequences.     

Footprints of divergent selection in natural populations of Castanopsis fargesii (Fagaceae)

Given predicted rapid climate change, an understanding of how environmental factors affect genetic diversity in natural populations is important. Future selection pressures are inherently unpredictable, so forest management policies should maintain both overall diversity and identify genetic markers associated with the environmental factors expected to change most rapidly, like temperature and rainfall. In this study, we genotyped 648 individuals in 28 populations of Castanopsis fargesii (Fagaceae) using 32 expressed sequence tag (EST)-derived microsatellite markers. After removing six loci that departed from Hardy–Weinberg equilibrium, we measured genetic variation, population structure and identified candidate loci putatively under selection by temperature and precipitation. We found that C. fargesii populations possessed high genetic diversity and moderate differentiation among them, indicating predominant outcrossing and few restrictions to gene flow. These patterns reduce the possible impact of stochastic effects or the influence of genetic isolation. Clear footprints of divergent selection at four loci were discovered. Frequencies of five alleles at these loci were strongly correlated with environmental factors, particularly extremes in precipitation. These alleles varied from being near fixation at one end of the gradient to being completely absent at the other. Our study species is an important forest tree in the subtropical regions of China and could have a major role in future management and reforestation plans. Our results demonstrate that the gene flow is widespread and abundant in natural populations, maintaining high diversity, while diversifying selection is acting on specific genomic regions.     

Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence–absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.     

Genetic diversity and gene flow in collapsed and healthy abalone fisheries

Overexploitation of marine species invariably results in population decline but can also have indirect effects on ecological processes such as larval dispersal and recruitment that ultimately affect genetic diversity and population resilience. We compared microsatellite DNA variation among depleted and healthy populations of the black-lip abalone Haliotis rubra from Tasmania, Australia, to determine if over-fishing had affected genetic diversity. We also used genetic data to assess whether variation in the scale and frequency of larval dispersal was linked to greater population decline in some regions than in others, and if larval dispersal was sufficient to facilitate natural recovery of depleted populations. Surprisingly, allelic diversity was higher in depleted populations than in healthy populations ( P <  0.05). Significant subdivision across hundreds of metres among our sampling sites ( F _ST = 0.026, P  < 0.01), coupled with assignment tests, indicated that larval dispersal is restricted in all regions studied, and that abalone populations across Tasmania are largely self-recruiting. Low levels of larval exchange appear to occur at the meso-scale (7–20 km), but age estimates based on shell size indicated that successful migration of larvae between any two sites may happen only once every few years. We suggest that genetic diversity may be higher in depleted populations due to the higher relative ratio of migrant to self-recruiting larvae. In addition, we expect that recovery of depleted abalone populations will be reliant on sources of larvae at the meso-scale (tens of km), but that natural recovery is only likely to occur on a timescale unacceptable to fishers and resource managers.     

Complex selection on a regulator of social cognition: Evidence of balancing selection,regulatory interactions and population differentiation in the prairie vole Avpr1a locus

Adaptive variation in social behaviour depends upon standing genetic variation, but we know little about how evolutionary forces shape genetic diversity relevant to brain and behaviour. In prairie voles (Microtus ochrogaster), variants at the Avpr1a locus predict expression of the vasopressin 1a receptor in the retrosplenial cortex (RSC), a brain region that mediates spatial and contextual memory; cortical V1aR abundance in turn predicts diversity in space use and sexual fidelity in the field. To examine the potential contributions of adaptive and neutral forces to variation at the Avpr1a locus, we explore sequence diversity at the Avpr1a locus and throughout the genome in two populations of wild prairie voles. First, we refine results demonstrating balancing selection at the locus by comparing the frequency spectrum of variants at the locus to a random sample of the genome. Next, we find that the four single nucleotide polymorphisms that predict high V1aR expression in the RSC are in stronger linkage disequilibrium than expected by chance despite high recombination among intervening variants, suggesting that epistatic selection maintains their association. Analysis of population structure and a haplotype network for two populations revealed that this excessive LD was unlikely to be due to admixture alone. Furthermore, the two populations differed considerably in the region shown to be a regulator of V1aR expression despite the extremely low levels of genomewide genetic differentiation. Together, our data suggest that complex selection on Avpr1a locus favours specific combinations of regulatory polymorphisms, maintains the resulting alleles at population‐specific frequencies, and may contribute to unique patterns of spatial cognition and sexual fidelity among populations.     

Minimal increase in genetic diversity enhances predation resistance

The importance of species diversity to emergent, ecological properties of communities is increasingly appreciated, but the importance of within‐species genetic diversity for analogous emergent properties of populations is only just becoming apparent. Here, the properties and effects of genetic variation on predation resistance in populations were assessed and the molecular mechanism underlying these emergent effects was investigated. Using biofilms of the ubiquitous bacterium Serratia marcescens, we tested the importance of genetic diversity in defending biofilms against protozoan grazing, a main source of mortality for bacteria in all natural ecosystems. S. marcescens biofilms established from wild‐type cells produce heritable, stable variants, which when experimentally combined, persist as a diverse assemblage and are significantly more resistant to grazing than either wild type or variant biofilms grown in monoculture. This diversity effect is biofilm‐specific, a result of either facilitation or resource partitioning among variants, with equivalent experiments using planktonic cultures and grazers resulting in dominance by a single resistant strain. The variants studied are all the result of single nucleotide polymorphisms in one regulatory gene suggesting that the benefits of genetic diversity in clonal biofilms can occur through remarkably minimal genetic change. The findings presented here provide a new insight on the integration of genetics and population ecology, in which diversity arising through minimal changes in genotype can have major ecological implications for natural populations.     

GM haplotype diversity of 82 populations over the world suggests a centrifugal model of human migrations

This study investigates the GM genetic relationships of 82 human populations, among which 10 represent original data, within and among the main broad geographic areas of the world. Different approaches are used: multidimensional scaling analysis and test for isolation by distance, to assess the correlation between genetic variation and spatial distributions; analysis of variance, to investigate the genetic structure at different hierarchical levels of population subdivision; genetic similarity map (geographic map distorted by available genetic information), to identify regions of high and low genetic variation; and minimal spanning network, to point out possible migration routes across continental areas. The results show that the GM polymorphism is characterized by one of the highest amounts of genetic variation observed so far among populations of different continents (Fct=0.3915, P < 0.0001). GM diversity can be explained by a model of isolation by distance (IBD) at most continental levels, with a particularly significant fit to IBD for the Middle East and Europe. Five peripheral regions of the world (Europe, west and south sub-Saharan Africa, Southeast Asia, and America) exhibit a low level of genetic diversity both within and among populations. By contrast, East and North African, Southwest Asian, and Northeast Asian populations are highly diverse and interconnected genetically by large genetic distances. Therefore, the observed GM variation can be explained by a "centrifugal model" of modern humans peopling history, involving ancient dispersals across a large intercontinental area spanning from East Africa to Northeast Asia, followed by recent migrations in peripheral geographic regions.     

The epigenetic footprint of poleward range‐expanding plants in apomictic dandelions

Epigenetic modifications, such as DNA methylation variation, can generate heritable phenotypic variation independent of the underlying genetic code. However, epigenetic variation in natural plant populations is poorly documented and little understood. Here, we test whether northward range expansion of obligate apomicts of the common dandelion (Taraxacum officinale) is associated with DNA methylation variation. We characterized and compared patterns of genetic and DNA methylation variation in greenhouse‐reared offspring of T. officinale that were collected along a latitudinal transect of northward range expansion in Europe. Genetic AFLP and epigenetic MS‐AFLP markers revealed high levels of local diversity and modest but significant heritable differentiation between sampling locations and between the southern, central and northern regions of the transect. Patterns of genetic and epigenetic variation were significantly correlated, reflecting the genetic control over epigenetic variation and/or the accumulation of lineage‐specific spontaneous epimutations, which may be selectively neutral. In addition, we identified a small component of DNA methylation differentiation along the transect that is independent of genetic variation. This epigenetic differentiation might reflect environment‐specific induction or, in case the DNA methylation variation affects relevant traits and fitness, selection of heritable DNA methylation variants. Such generated epigenetic variants might contribute to the adaptive capacity of individual asexual lineages under changing environments. Our results highlight the potential of heritable DNA methylation variation to contribute to population differentiation along ecological gradients. Further studies are needed using higher resolution methods to understand the functional significance of such natural occurring epigenetic differentiation.     

Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping

The extraordinary phenotypic diversity of dog breeds has been sculpted by a unique population history accompanied by selection for novel and desirable traits. Here we perform a comprehensive analysis using multiple test statistics to identify regions under selection in 509 dogs from 46 diverse breeds using a newly developed high-density genotyping array consisting of >170,000 evenly spaced SNPs. We first identify 44 genomic regions exhibiting extreme differentiation across multiple breeds. Genetic variation in these regions correlates with variation in several phenotypic traits that vary between breeds, and we identify novel associations with both morphological and behavioral traits. We next scan the genome for signatures of selective sweeps in single breeds, characterized by long regions of reduced heterozygosity and fixation of extended haplotypes. These scans identify hundreds of regions, including 22 blocks of homozygosity longer than one megabase in certain breeds. Candidate selection loci are strongly enriched for developmental genes. We chose one highly differentiated region, associated with body size and ear morphology, and characterized it using high-throughput sequencing to provide a list of variants that may directly affect these traits. This study provides a catalogue of genomic regions showing extreme reduction in genetic variation or population differentiation in dogs, including many linked to phenotypic variation. The many blocks of reduced haplotype diversity observed across the genome in dog breeds are the result of both selection and genetic drift, but extended blocks of homozygosity on a megabase scale appear to be best explained by selection. Further elucidation of the variants under selection will help to uncover the genetic basis of complex traits and disease.     

Variation in the human TAS1R taste receptor genes

We have performed a comprehensive evaluation of single-nucleotide polymorphisms (SNPs) and haplotypes in the human TAS1R gene family, which encodes receptors for sweet and umami tastes. Complete DNA sequences of TAS1R1-, TAS1R2-, and TAS1R3-coding regions, obtained from 88 individuals of African, Asian, European, and Native American origin, revealed substantial coding and noncoding diversity: polymorphisms are common in these genes, and polymorphic sites and SNP frequencies vary widely in human populations. The genes TAS1R1 and TAS1R3, which encode proteins that act as a dimer to form the umami (glutamate) taste receptor, showed less variation than the TAS1R2 gene, which acts as a dimer with TAS1R3 to form the sweet taste receptor. The TAS1R3 gene, which encodes a subunit common to both the sweet and umami receptors, was the most conserved. Evolutionary genetic analysis indicates that these variants have come to their current frequencies under natural selection during population growth and support the view that the coding sequence variants affect receptor function. We propose that human populations likely vary little with respect to umami perception, which is controlled by one major form of the receptor that is optimized for detecting glutamate but may vary much more with respect to sweet perception.     

Low Frequency Variants,Collapsed Based on Biological Knowledge,Uncover Complexity of Population Stratification in 1000 Genomes Project Data

Analyses investigating low frequency variants have the potential for explaining additional genetic heritability of many complex human traits. However, the natural frequencies of rare variation between human populations strongly confound genetic analyses. We have applied a novel collapsing method to identify biological features with low frequency variant burden differences in thirteen populations sequenced by the 1000 Genomes Project. Our flexible collapsing tool utilizes expert biological knowledge from multiple publicly available database sources to direct feature selection. Variants were collapsed according to genetically driven features, such as evolutionary conserved regions, regulatory regions genes, and pathways. We have conducted an extensive comparison of low frequency variant burden differences (MAF<0.03) between populations from 1000 Genomes Project Phase I data. We found that on average 26.87% of gene bins, 35.47% of intergenic bins, 42.85% of pathway bins, 14.86% of ORegAnno regulatory bins, and 5.97% of evolutionary conserved regions show statistically significant differences in low frequency variant burden across populations from the 1000 Genomes Project. The proportion of bins with significant differences in low frequency burden depends on the ancestral similarity of the two populations compared and types of features tested. Even closely related populations had notable differences in low frequency burden, but fewer differences than populations from different continents. Furthermore, conserved or functionally relevant regions had fewer significant differences in low frequency burden than regions under less evolutionary constraint. This degree of low frequency variant differentiation across diverse populations and feature elements highlights the critical importance of considering population stratification in the new era of DNA sequencing and low frequency variant genomic analyses.     

Genetic Evaluation of the Efficacy of In Situ and Ex Situ Conservation of <Emphasis Type="Italic">Parashorea chinensis</Emphasis> (Dipterocarpaceae) in Southwestern China

The majority of research in genetic diversity yields recommendations rather than actual conservation achievements. We assessed the efficacy of actual in situ and ex situ efforts to conserve Parashorea chinensis (Dipterocarpaceae) against the background of the geographic pattern of genetic variation of this species. Samples from seven natural populations, including three in a nature reserve, and one ex situ conservation population were studied. Across the natural populations, 47.8% of RAPD loci were polymorphic; only 20.8% on average varied at the population level. Mean population genetic diversity was 0.787 within natural populations and 1.410 for the whole species. Significant genetic differentiation among regions and isolation by distance were present on larger scales (among regions). AMOVA revealed that the majority of the among-population variation occurred among regions rather than among populations within regions. Regression analysis, Mantel test, principal coordinates analysis, and cluster analysis consistently demonstrated increasing genetic isolation with increasing geographic distance. Genetic differentiation within the region was quite low compared to that among regions. Multilocus spatial autocorrelation analysis of these three populations revealed random distribution of genetic variation in two populations, but genetic clustering was detected in the third population. The ex situ conserved population contained a medium level of genetic variation compared with the seven natural populations; it contained 77.1% of the total genetic variation of this species and 91% of the moderate to high frequency RAPD fragments (f > 0.05). Exclusive bands were detected in natural populations, but none were found in the ex situ conserved population. The populations protected in the nature reserve contained most of the genetic variation of the whole species, with 81.4% of the total genetic variation and 95.7% of the fragments with moderate to high frequency (f > 0.05) of this species conserved. The results show that the ex situ conserved population does not contain enough genetic variation to meet the need of release in the future, and that more extensive ex situ sampling in natural populations TY, NP, HK, and MG is needed. The in situ conserved population contains representative genetic variation to maintain long-term survival and evolutionary processes of P. chinensis.