Uncovering the genetic basis of adaptive change: on the intersection of landscape genomics and theoretical population genetics

A workshop recently held at the Ecole Polytechnique Federale de Lausanne (EPFL, Switzerland) was dedicated to understanding the genetic basis of adaptive change, taking stock of the different approaches developed in theoretical population genetics and landscape genomics and bringing together knowledge accumulated in both research fields. Indeed, an important challenge in theoretical population genetics is to incorporate effects of demographic history and population structure. But important design problems (e.g. focus on populations as units, focus on hard selective sweeps, no hypothesis-based framework in the design of the statistical tests) reduce their capability of detecting adaptive genetic variation. In parallel, landscape genomics offers a solution to several of these problems and provides a number of advantages (e.g. fast computation, landscape heterogeneity integration). But the approach makes several implicit assumptions that should be carefully considered (e.g. selection has had enough time to create a functional relationship between the allele distribution and the environmental variable, or this functional relationship is assumed to be constant). To address the respective strengths and weaknesses mentioned above, the workshop brought together a panel of experts from both disciplines to present their work and discuss the relevance of combining these approaches, possibly resulting in a joint software solution in the future.

The role of parasite-driven selection in shaping landscape genomic structure in red grouse (Lagopus lagopus scotica)

This article is protected by copyright. All rights reserved. Acknowledgements This study was funded by a BBSRC studentship (MAW) and NERC grants NE/H00775X/1 and NE/D000602/1 (SBP). The authors are grateful to Mario Röder and Keliya Bai for fieldwork assistance, and all estate owners, factors and keepers for access to field sites, most particularly MJ Taylor and Mike Nisbet (Airlie), Neil Brown (Allargue), RR Gledson and David Scrimgeour (Delnadamph), Andrew Salvesen and John Hay (Dinnet), Stuart Young and Derek Calder (Edinglassie), Kirsty Donald and David Busfield (Glen Dye), Neil Hogbin and Ab Taylor (Glen Muick), Alistair Mitchell (Glenlivet), Simon Blackett, Jim Davidson and Liam Donald (Invercauld), Richard Cooke and Fred Taylor† (Invermark), Shaila Rao and Christopher Murphy (Mar Lodge), and Ralph Peters and Philip Astor (Tillypronie). Data accessibility • Genotype data (DataDryad: doi:10.5061/dryad.4t7jk) • Metadata (information on sampling sites, phenotypes and medication regimen) (DataDryad: doi:10.5061/dryad.4t7jk)

Phylogeography and population genetics of north European Atlantic salmon (Salmo salar L.)

Although abundant in the number of individuals, the Atlantic salmon may be considered as a threatened species in many areas of its native distribution range. Human activities such as building of power plant dams, offshore overfishing, pollution, clearing of riverbeds for timber floating and badly designed stocking regimes have diminished the distribution of Atlantic salmon. As a result of this, many of the historical populations both in Europe and northern America have gone extinct or are severely depressed. In fact, only 1% of Atlantic salmon existing today are of natural origin, the rest being farmed salmon. All of this has lead to a vast amount of research and many restoration programmes aiming to bring Atlantic salmon back to rivers from where it has vanished. However, many of the restoration programmes conducted thus far have been unsuccessful due to inadequate scientific research or lack of its implementation, highlighting the fact that more research is needed to fully understand the biology of this complex species. The White and Barents Seas in northwest Russia are among the last regions in Europe where Atlantic salmon populations are still stable, thus forming an important source of biodiversity for the entire European region. Salmon stocks from this area are also of immense economic and social importance for the local people in the form of fishing tourism. The main aim of this thesis was to elucidate the post-glacial history and population genetic structure of north European and particularly northwest Russian Atlantic salmon, both of which are aspects of great importance for the management and conservation of the species. Throughout the whole thesis, these populations were studied by utilizing microsatellites as the main molecular tool. One of the most important discoveries of the thesis was the division of Atlantic salmon from the White and Barents Seas into four separate clusters, which has not been observed in previous studies employing nuclear markers although is supported by mtDNA studies. Populations from the western Barents Sea clustered together with the northeast Atlantic populations into a clearly distinguishable group while populations from the White Sea and eastern Barents Sea were separated into three additional groups. This has important conservation implications as this thesis clearly indicates that conservation of populations from all of the observed clusters is warranted in order to conserve as much of the genetic diversity as possible in this area. The thesis also demonstrates how differences in population life histories within a species, migratory behaviour in this case, and in their phylogeographic origin affect the genetic characteristics of populations, namely diversity and divergence levels. The anadromous populations from the Atlantic Ocean, White Sea and Barents Sea possessed higher levels of genetic diversity than the anadromous populations form the Baltic Sea basin. Among the non-anadromous populations the result was the opposite: the Baltic freshwater populations were more variable. This emphasises the importance of taking the life history of a population into consideration when developing conservation strategies: due to the limited possibilities for new genetic diversity to be generated via gene flow, it is expected that freshwater Atlantic salmon populations would be more vulnerable to extinction following a population crash and thus deserve a high conservation status. In the last chapter of this thesis immune relevant marker loci were developed and screened for signatures of natural selection along with loci linked to genes with other functions or no function at all. Also, a novel landscape genomics method, which combines environmental information with molecular data, was employed to investigate whether immune relevant markers displayed significant correlations to various environmental variables more frequently than other loci. Indications of stronger selection pressure among immune-relevant loci compared to non-immune relevant EST-linked loci was found but further studies are needed to evaluate whether it is a common phenomenon in Atlantic salmon.

Update of the FANTOM web resource: from mammalian transcriptional landscape to its dynamic regulation.

The international Functional Annotation Of the Mammalian Genomes 4 (FANTOM4) research collaboration set out to better understand the transcriptional network that regulates macrophage differentiation and to uncover novel components of the transcriptome employing a series of high-throughput experiments. The primary and unique technique is cap analysis of gene expression (CAGE), sequencing mRNA 5'-ends with a second-generation sequencer to quantify promoter activities even in the absence of gene annotation. Additional genome-wide experiments complement the setup including short RNA sequencing, microarray gene expression profiling on large-scale perturbation experiments and ChIP-chip for epigenetic marks and transcription factors. All the experiments are performed in a differentiation time course of the THP-1 human leukemic cell line. Furthermore, we performed a large-scale mammalian two-hybrid (M2H) assay between transcription factors and monitored their expression profile across human and mouse tissues with qRT-PCR to address combinatorial effects of regulation by transcription factors. These interdependent data have been analyzed individually and in combination with each other and are published in related but distinct papers. We provide all data together with systematic annotation in an integrated view as resource for the scientific community (http://fantom.gsc.riken.jp/4/). Additionally, we assembled a rich set of derived analysis results including published predicted and validated regulatory interactions. Here we introduce the resource and its update after the initial release.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

Combining niche modelling and landscape genetics to study local adaptation: A novel approach illustrated using alpine plants

Understanding the factors that shape adaptive genetic variation across species niches has become of paramount importance in evolutionary ecology, especially to understand how adaptation to changing climate affects the geographic range of species. The distribution of adaptive alleles in the ecological niche is determined by the emergence of novel mutations, their fitness consequences and gene flow that connects populations across species niches. Striking demographical differences and source sink dynamics of populations between the centre and the margin of the niche can play a major role in the emergence and spread of adaptive alleles. Although some theoretical predictions have long been proposed, the origin and distribution of adaptive alleles within species niches remain untested. In this paper, we propose and discuss a novel empirical approach that combines landscape genetics with species niche modelling, to test whether alleles that confer local adaptation are more likely to occur in either marginal or central populations of species niches. We illustrate this new approach by using a published data set of 21 alpine plant species genotyped with a total of 2483 amplified fragment length polymorphisms (AFLP), distributed over more than 1733 sampling sites across the Alps. Based on the assumption that alleles that were statistically associated with environmental variables were adaptive, we found that adaptive alleles in the margin of a species niche were also present in the niche centre, which suggests that adaptation originates in the niche centre. These findings corroborate models of species range evolution, in which the centre of the niche contributes to the emergence of novel adaptive alleles, which diffuse towards niche margins and facilitate niche and range expansion through subsequent local adaptation. Although these results need to be confirmed via fitness measurements in natural populations and functionally characterised genetic sequences, this study provides a first step towards understanding how adaptive genetic variation emerges and shapes species niches and geographic ranges along environmental gradients.

Using gene-network landscape to dissect genotype effects of TCF7L2 genetic variant on diabetes and cardiovascular risk

Vaquero AR, Ferreira NE, Omae SV, Rodrigues MV, Teixeira SK, Krieger JE, Pereira AC. Using gene-network landscape to dissect genotype effects of TCF7L2 genetic variant on diabetes and cardiovascular risk. Physiol Genomics 44: 903-914, 2012. First published August 7, 2012; doi:10.1152/physiolgenomics.00030.2012.-The single nucleotide polymorphism (SNP) within the TCF7L2 gene, rs7903146, is, to date, the most significant genetic marker associated with Type 2 diabetes mellitus (T2DM) risk. Nonetheless, its functional role in disease pathology is poorly understood. The aim of the present study was to investigate, in vascular smooth muscle cells from 92 patients undergoing aortocoronary bypass surgery, the contribution of this SNP in T2DM using expression levels and expression correlation comparison approaches, which were visually represented as gene interaction networks. Initially, the expression levels of 41 genes (seven TCF7L2 splice forms and 40 other T2DM relevant genes) were compared between rs7903146 wild-type (CC) and T2DM-risk (CT + TT) genotype groups. Next, we compared the expression correlation patterns of these 41 genes between groups to observe if the relationships between genes were different. Five TCF7L2 splice forms and nine genes showed significant expression differences between groups. RXR alpha gene was pinpointed as showing the most different expression correlation pattern with other genes. Therefore, T2DM risk alleles appear to be influencing TCF7L2 splice form's expression in vascular smooth muscle cells, and RXR alpha gene is pointed out as a treatment target candidate for risk reduction in individuals with high risk of developing T2DM, especially individuals harboring TCF7L2 risk genotypes.

Analysis of opo cis-regulatory landscape uncovers Vsx2 requirement in early eye morphogenesis

The self-organized morphogenesis of the vertebrate optic cup entails coupling the activation of the retinal gene regulatory network to the constriction-driven infolding of the retinal epithelium. Yet the genetic mechanisms underlying this coordination remain largely unexplored. Through phylogenetic footprinting and transgenesis in zebrafish, here we examine the cis-regulatory landscape of opo, an endocytosis regulator essential for eye morphogenesis. Among the different conserved enhancers identified, we isolate a single retina-specific element (H6_10137) and show that its activity depends on binding sites for the retinal determinant Vsx2. Gain- and loss-of-function experiments and ChIP analyses reveal that Vsx2 regulates opo expression through direct binding to this retinal enhancer. Furthermore, we show that vsx2 knockdown impairs the primary optic cup folding. These data support a model by which vsx2, operating through the effector gene opo, acts as a central transcriptional node that coordinates neural retina patterning and optic cup invagination in zebrafish.

Genomics and proteomics approaches to the study of cancer-stroma interactions

Background: The development and progression of cancer depend on its genetic characteristics as well as on the interactions with its microenvironment. Understanding these interactions may contribute to diagnostic and prognostic evaluations and to the development of new cancer therapies. Aiming to investigate potential mechanisms by which the tumor microenvironment might contribute to a cancer phenotype, we evaluated soluble paracrine factors produced by stromal and neoplastic cells which may influence proliferation and gene and protein expression. Methods: The study was carried out on the epithelial cancer cell line (Hep-2) and fibroblasts isolated from a primary oral cancer. We combined a conditioned-medium technique with subtraction hybridization approach, quantitative PCR and proteomics, in order to evaluate gene and protein expression influenced by soluble paracrine factors produced by stromal and neoplastic cells. Results: We observed that conditioned medium from fibroblast cultures (FCM) inhibited proliferation and induced apoptosis in Hep-2 cells. In neoplastic cells, 41 genes and 5 proteins exhibited changes in expression levels in response to FCM and, in fibroblasts, 17 genes and 2 proteins showed down-regulation in response to conditioned medium from Hep-2 cells (HCM). Nine genes were selected and the expression results of 6 down-regulated genes (ARID4A, CALR, GNB2L1, RNF10, SQSTM1, USP9X) were validated by real time PCR. Conclusions: A significant and common denominator in the results was the potential induction of signaling changes associated with immune or inflammatory response in the absence of a specific protein.

Deep weathering and landscape evolution in a tropical plateau

The thick weathering profiles of humid tropical areas are an important, yet generally neglected, source of information on landscape evolution. Six complete profiles of the weathering mantle were sampled by drilling on the three stepped levels of the Campos do Jordao Plateau, on the NW flank of the Continental Rift of Southeastern Brazil. Mineralogical and micromorphological analyses of drill core samples, complemented by geochemical interpretations and by previous data on the upper saprolite, indicate continuity of a general lateritic trend during the entire process of mantle formation. Lateritization phases of different intensity were defined and considered to reflect adjustment to changes in environmental conditions created by the gradual uplift of the plateau to its present position. Older and more superficial materials related to intense lateritic weathering are characterized by allitization with direct formation of gibbsite from silicates, probably related to tropical climates existing immediately before the formation of the continental rift, during the Paleogene, and also before any significant increase in altitude. Monosialitization phase with general kaolinitization and restricted indirect formation of gibbsite from silicates could be associated to less aggressive climates that followed the Neogene (Miocene?) accentuation of uplift rates along the continental rift. The changes produced by uplift in the tropical climate eventually favored the development of a podzolization trend in soils above 1800m. (C) 2011 Elsevier BM. All rights reserved.

Landscape unit discrimination for pedological surveys by orbital spectral response

Landscape unit discrimination for pedological surveys by orbital spectral response. The objective of tins study was compare two soil survey methods. The first was performed by methods traditionally used to distinguish landscape units and soil class discrimination. The second was based on soil class distinction through orbital spectral response. In order to establish soil characteristics and their classification, soil samples were collected at two depths in a grid system, with a distance of 500 meters between points. With these samples, physical and chemical analyses were carried out. In the sampling points, the apparent reflectance of the soil, front the orbital image, was determined and, through cluster analysis landscape units were established. In order to evaluate the resemblance reliability between the landscape units established in each method, the Kappa index was used, the value set for the confusion matrix was 0.43, indicating high quality in the comparison, showing that the non-conventional method was as close as the one carried out by photointerpretation.

Genomics of pyrrolnitrin biosynthetic loci: evidence for conservation and whole-operon mobility within Gram-negative bacteria

Pyrrolnitrin (PRN) is a tryptophan-derived secondary metabolite produced by a narrow range of Gram-negative bacteria. The PRN biosynthesis by rhizobacteria presumably has a key role in their life strategies and in the biocontrol of plant diseases. The biosynthetic operon that encodes the pathway that converts tryptophan to PRN is composed of four genes, prnA through D, whose diversity, genomic context and spread over bacterial genomes are poorly understood. Therefore, we launched an endeavour aimed at retrieving, by in vitro and in silico means, diverse bacteria carrying the prnABCD biosynthetic loci in their genomes. Analysis of polymorphisms of the prnD gene sequences revealed a high level of conservation between Burkholderia, Pseudomonas and Serratia spp. derived sequences. Whole-operon- and prnD-based phylogeny resulted in tree topologies that are incongruent with the taxonomic status of the evaluated strains as predicted by 16S rRNA gene phylogeny. The genomic composition of c. 20 kb DNA fragments containg the PRN operon varied in different strains. Highly conserved and distinct transposase-encoding genes surrounding the PRN biosynthetic operons of Burkholderia pseudomallei strains were found. A prnABCD-deprived genomic region in B. pseudomallei strain K96243 contained the same gene composition as, and shared high homology with, the flanking regions of the PRN operon in B. pseudomallei strains 668, 1106a and 1710b. Our results strongly suggest that the PRN biosynthetic operon is mobile. The extent, frequency and promiscuity of this mobility remain to be understood.

Geochronology of weathering and landscape evolution, Dugald River valley, NW Queensland, Australia

40Ar/39Ar laser incremental heating analyses of individual grains of supergene jarosite, alunite, and cryptomelane from weathering profiles in the Dugald River area, Queensland, Australia, show a strong positive correlation between a sample’s age and its elevation. We analyzed 125 grains extracted from 35 hand specimens collected from weathering profiles at 11 sites located at 3 distinct elevations. The highest elevation profile hosts the oldest supergene minerals, whereas progressively younger samples occur at lower positions in the landscape. The highest elevation sampling sites (three sites), located on top of an elongated mesa (255 to 275 m elevation), yield ages in the 16 to 12 Ma range. Samples from an intermediate elevation site (225 to 230 m elevation) yield ages in the 6 to 4 Ma range. Samples collected at the lowest elevation sites (200 to 220 m elevation) yield ages in the 2.2 to 0.8 Ma interval. Grains of supergene alunite, jarosite, and cryptomelane analyzed from individual single hand specimens yield reproducible results, confirming the suitability of these minerals to 40Ar/39Ar geochronology. Multiple samples collected from the same site also yield reproducible results, indicating that the ages measured are true precipitation ages for the samples analyzed. Different sites, up to 3 km apart, sampled from weathering profiles at the same elevation again yield reproducible results. The consistency of results confirms that 40Ar/39Ar geochronology of supergene jarosite, alunite, and cryptomelane yields ages of formation of weathering profiles, providing a reliable numerical basis for differentiating and correlating these profiles. The age versus elevation relationship obtained suggest that the stepped landscapes in the Dugald River area record a progressive downward migration of a relatively flat weathering front. The steps in the landscape result from differential erosion of previously weathered bedrock displaying different susceptibility to weathering and contrasting resistance to erosion. Combined, the age versus elevation relationships measured yield a weathering rate of 3.8 m. Myr−1 (for the past 15 Ma) if a descending subhorizontal weathering front is assumed. The results also permit the calculation of the erosion rate of the more easily weathered and eroded lithologies, assuming an initially flat landscape as proposed in models of episodic landscape development. The average erosion rate for the past 15 Ma is 3.3 m. Myr−1, consistent with erosion rates obtained by cosmogenic isotope studies in the region.