Detecting positive selection from genome scans of linkage disequilibrium

Background

During the lifetime of a fermenter culture, the soil bacterium S. coelicolor undergoes a major metabolic switch from exponential growth to antibiotic production. We have studied gene expression patterns during this switch, using a specifically designed Affymetrix genechip and a high-resolution time-series of fermenter-grown samples.

Results

Surprisingly, we find that the metabolic switch actually consists of multiple finely orchestrated switching events. Strongly coherent clusters of genes show drastic changes in gene expression already many hours before the classically defined transition phase where the switch from primary to secondary metabolism was expected. The main switch in gene expression takes only 2 hours, and changes in antibiotic biosynthesis genes are delayed relative to the metabolic rearrangements. Furthermore, global variation in morphogenesis genes indicates an involvement of cell differentiation pathways in the decision phase leading up to the commitment to antibiotic biosynthesis.

Conclusions

Our study provides the first detailed insights into the complex sequence of early regulatory events during and preceding the major metabolic switch in S. coelicolor, which will form the starting point for future attempts at engineering antibiotic production in a biotechnological setting.     

A new approach for using genome scans to detect recent positive selection in the human genome

Genome-wide scanning for signals of recent positive selection is essential for a comprehensive and systematic understanding of human adaptation. Here, we present a genomic survey of recent local selective sweeps, especially aimed at those nearly or recently completed. A novel approach was developed for such signals, based on contrasting the extended haplotype homozygosity (EHH) profiles between populations. We applied this method to the genome single nucleotide polymorphism (SNP) data of both the International HapMap Project and Perlegen Sciences, and detected widespread signals of recent local selection across the genome, consisting of both complete and partial sweeps. A challenging problem of genomic scans of recent positive selection is to clearly distinguish selection from neutral effects, given the high sensitivity of the test statistics to departures from neutral demographic assumptions and the lack of a single, accurate neutral model of human history. We therefore developed a new procedure that is robust across a wide range of demographic and ascertainment models, one that indicates that certain portions of the genome clearly depart from neutrality. Simulations of positive selection showed that our tests have high power towards strong selection sweeps that have undergone fixation. Gene ontology analysis of the candidate regions revealed several new functional groups that might help explain some important interpopulation differences in phenotypic traits.     

Scanning the human genome for signals of selection

The search for adaptive evolution in the human genome has reached a new era with the advent of genome-wide surveys of genetic variation. However, making sense, let alone use, of such experiments is far from straightforward. Key problems include the way in which the data have been collected, the need to control for factors such as population history and variable recombination rates, which influence the discovery rates for both true and false positives, and the inherent difficulty of falsification. Nevertheless, recent work has shown that genome scans can be used to identify both functional polymorphisms underlying selected traits and entire classes of genes enriched for signals of adaptation.     

Codon volatility as an indicator of positive selection: data from eukaryotic genome comparisons

It has been suggested that codon volatility (the proportion of the point-mutation neighbors of a codon that encode different amino acids) can be used as an index of past positive selection. We compared codon volatility with patterns of synonymous and nonsynonymous nucleotide substitution in genome-wide comparisons of orthologous genes between three pairs of related genomes: (1) the protists Plasmodium falciparum and P. yoelii, (2) the fungi Saccharomyces cerevisiae and S. paradoxus, and (3) the mammals mouse and rat. Codon volatility was not consistently associated with an elevated rate of nonsynonymous substitution, as would be expected under positive selection. Rather, the most consistent and powerful correlate of elevated codon volatility was nucleotide content at the second codon position, as expected, given the nature of the genetic code.     

A map of recent positive selection in the human genome

The identification of signals of very recent positive selection provides information about the adaptation of modern humans to local conditions. We report here on a genome-wide scan for signals of very recent positive selection in favor of variants that have not yet reached fixation. We describe a new analytical method for scanning single nucleotide polymorphism (SNP) data for signals of recent selection, and apply this to data from the International HapMap Project. In all three continental groups we find widespread signals of recent positive selection. Most signals are region-specific, though a significant excess are shared across groups. Contrary to some earlier low resolution studies that suggested a paucity of recent selection in sub-Saharan Africans, we find that by some measures our strongest signals of selection are from the Yoruba population. Finally, since these signals indicate the existence of genetic variants that have substantially different fitnesses, they must indicate loci that are the source of significant phenotypic variation. Though the relevant phenotypes are generally not known, such loci should be of particular interest in mapping studies of complex traits. For this purpose we have developed a set of SNPs that can be used to tag the strongest ∼250 signals of recent selection in each population.     

Searching for evidence of positive selection in the human genome using patterns of microsatellite variability

Both natural selection and nonequilibrium population-level processes can lead to a skew in the frequency distribution of polymorphisms. Population-level processes are expected to affect all loci in a roughly equal fashion, whereas selection will affect only some regions of the genome. We conducted a sliding-window analysis of the frequency distribution of microsatellite polymorphisms across the human genome to identify regions that may be under positive selection. The analysis was based on a published data set of 5,257 mapped microsatellites in individuals of European ancestry. Observed and expected numbers of alleles were compared under a stepwise mutation model (SMM) using analytical formulae. Observed and expected heterozygosities were compared under a SMM using coalescent simulations. The two sets of analyses gave similar results. Approximately one-fourth of all loci showed a significant deficit of heterozygosity, consistent with a recent population expansion. Forty-three windows were identified with extreme skews in the frequency distribution of polymorphisms (in the direction of a deficit of heterozygosity, given the number of alleles). If these extreme windows are tracking selection at linked sites, theory predicts that they should be more common in regions of the genome with less recombination. We tested this prediction by comparing recombination rates in these extreme windows and in other regions of the genome and found that extreme windows had a significantly lower recombination rate than the genomic average. The proportion of extreme windows was significantly higher on the X chromosome than on the autosomes. Moreover, all the windows with extreme skews on the X chromosome were found in two clusters near the centromere; both these clusters exhibit markedly reduced recombination rates. These analyses point to regions of the genome that may recently have been subject to positive selection. These results also suggest that the effects of positive selection may be more pronounced on the X chromosome than on the autosomes in humans.     

Empirical bayes method for incorporating data from multiple genome scans

Individual genome scans tend to have low power and can produce markedly biased estimates of QTL effects. Further, the confidence interval for their location is often prohibitively large for subsequent fine mapping and positional cloning. Given that a large number of genome scans have been conducted, not to mention the large number of variables and subsets tested, it is difficult to confidently rule out type 1 error as an explanation for significant effects even when there is apparent replication in a separate data set. We adapted Empirical Bayes (EB) methods [1] to analyze data from multiple genome scans simultaneously and alleviate each of these problems while still allowing for different QTL population effects across studies. We investigated the effects of using the EB method to include data from background studies to update the results of a single study of interest via simulation and demonstrated that it has a stable confidence level over a wide range of parameters defining the background studies and increased the power to detect linkage, even when some of the background studies were null or had QTL effect at other markers. This EB method for incorporating data from multiple studies into genome scan analyses seems promising.     

Uninformative polymorphisms bias genome scans for signatures of selection

Background

A female preference for intense sexual visual signals is widespread in animals. Although the preferences for a signal per se and for the intensity of the signal were often regarded to have the identical origin, no study has demonstrated if this is true. It was suggested that the female fiddler crabs prefer males with courtship structures because of direct benefit to escape predation. Here we tested if female preference for both components (i.e. presence and size) of the courtship structure in Uca lactea is from the sensory bias to escape predation. If both components have the identical origin, females should show the same response to different-sized courtship structures regardless of predation risk.

Results

First, we observed responses of mate-searching female U. lactea to courting males with full-sized, half-sized and no semidomes which were experimentally manipulated. Females had a directional preference for males with bigger semidomes within normal variation. Thereafter, we tested the effect of predation risk on the female bias in the non-courtship context. When threatened by an avian mock predator, females preferentially approached burrows with full-sized semidomes regardless of reproductive cycles (i.e. reproductive periods and non-reproductive periods). When the predator cue was absent, however, females preferred burrows with semidomes without discriminating structure size during reproductive periods but did not show any bias during non-reproductive periods.

Conclusions

Results indicate that selection for the size of courtship structures in U. lactea may have an origin in the function to reduce predation risk, but that the preference for males with structures may have evolved by female choice, independent of predation pressure.     

Transpositional landscape of the rice genome revealed by paired-end mapping of high-throughput re-sequencing data

Transposable elements (TEs) are mobile entities that densely populate most eukaryotic genomes and contribute to both their structural and functional dynamics. However, most TE-related sequences in both plant and animal genomes correspond to inactive, degenerated elements, due to the combined effect of silencing pathways and elimination through deletions. One of the major difficulties in fully characterizing the molecular basis of genetic diversity of a given species lies in establishing its genome-wide transpositional activity. Here, we provide an extensive survey of the transpositional landscape of a plant genome using a deep sequencing strategy. This was achieved through paired-end mapping of a fourfold coverage of the genome of rice mutant line derived from an in vitro callus culture using Illumina technology. Our study shows that at least 13 TE families are active in this genotype, causing 34 new insertions. This next-generation sequencing-based strategy provides new opportunities to quantify the impact of TEs on the genome dynamics of the species.     

Automating the identification of DNA variations using quality-based fluorescence re-sequencing: analysis of the human mitochondrial genome.

Diagnostic re-sequencing plays a central role in medical and evolutionary genetics. In this report we describe a process that applies fluorescence-based re-sequencing and an integrated set of analysis tools to automate and simplify the identification of DNA variations using the human mitochondrial genome as a model system. Two programs used in genome sequence analysis (Phred, a base-caller, and Phrap, a sequence assembler) are applied to assess the quality of each base call across the sequence. Potential DNA variants are automatically identified and 'tagged' by comparing the assembled sequence with a reference sequence. We also show that employing the Consed program to display a set of highly annotated reference sequences greatly simplifies data analysis by providing a visual database containing information on the location of the PCR primers, coding and regulatory sequences and previously known DNA variants. Among the 12 genomes sequenced 378 variants including 29 new variants were identified along with two heteroplasmic sites, automatically detected by the PolyPhred program. Overall we document the ease and speed of performing high quality and accurate fluorescence-based re-sequencing on long tracts of DNA as well as the application of new approaches to automatically find and view DNA variants among these sequences.     

Delineation of signals required for thymocyte positive selection

Peptide/MHC complexes capable of inducing positive selection in mouse fetal thymic organ cultures fail to do so in suspension culture. Furthermore, this type of culture does not promote initial stages of differentiation, such as coreceptor down-modulation, unless peptides used for stimulation have (at least) weak agonist activity. We show in this study that signals provided in suspension culture by nonagonist peptide/MHC complexes on the surface of macrophages, even though apparently silent, are sufficient to promote complete phenotypic differentiation when CD4+CD8+ thymocytes are subsequently placed in a proper anatomical setting. Furthermore, the synergistic actions of suboptimal concentrations of phorbol esters and nonagonist peptide/MHC complexes can make the initial stages of positive selection visible, without converting maturation into negative selection. Thus, the correlation between efficiency of positive selection and the degree of coreceptor down-modulation on CD4+CD8+ thymocytes is not linear. Furthermore, these results suggest that the unique role of thymic stromal cells in positive selection is related not to presentation of self-peptide/MHC complexes, but most likely to another ligand.     

Endothelial-like cells derived from human CD14 positive monocytes

In the present study, we show that endothelial-like cells (ELCs) can develop from human CD14-positive mononuclear cells (CD14 cells) in the presence of angiogenic growth factors. The CD14 cells became loosely adherent within 24 h of culture and subsequently underwent a distinct process of morphological transformation to caudated or oval cells with eccentric nuclei. After 1 week in culture the cells showed a clear expression of endothelial cell markers, including von Willebrand factor (vWF), CD144 (VE-cadherin), CD105 (endoglin), acetylated low-density lipoprotein (AC-LDL)-receptor, CD36 (thrombospondin receptor), FLT-1, which is vascular endothelial cell growth factor (VEGF) receptor-1, and, to a weaker extent, KDR (VEGF receptor-2). Furthermore, in these cells structures resembling Weibel-Palade bodies at different storage stages were identified by electron microscopy, and upon culturing on three-dimensional fibrin gels the cells build network-like structures. In addition, cell proliferation and vWF expression was stimulated by VEGF, and the endothelial cell adhesion molecules CD54 (ICAM-1), and CD106 (VCAM-1) became transiently inducible by tumor necrosis factor-alpha (TNF-alpha). In contrast, the dendritic markers CD1a, and CD83 were not expressed to any significant extent. The expression of CD68, CD80 (B7-1), CD86 (B7-2), HLA-DR and CD36 may also suggest that ELCs might be related to macrophages, sinus lining or microvascular endothelial cells. Taken together, our observations indicate that ELCs can differentiate from cells of the monocytic lineage, suggesting a closer relationship between the monocyte/macrophage- and the endothelial cell systems than previously supposed.     

Management, presentation and interpretation of genome scans using GSCANDB

MOTIVATION: Advances in high-throughput genotyping have made it possible to carry out genome-wide association studies using very high densities of genetic markers. This has led to the problem of the storage, management, quality control, presentation and interpretation of results. In order to achieve a successful outcome, it may be necessary to analyse the data in different ways and compare the results with genome annotations and other genome scans. RESULTS: We created GSCANDB, a database for genome scan data, using a MySQL backend and Perl-CGI web interface. It displays genome scans of multiple phenotypes analysed in different ways and projected onto genome annotations derived from EnsMart. The current version is optimized for analysis of mouse data, but is customizable to other species. AVAILABILITY: Source code and example data are available under the GPL, in versions tailored to either human or mouse association studies, from http://gscan.well.ox.ac.uk/software.     

Genome-wide scans provide evidence for positive selection of genes implicated in Lassa fever

Rapidly evolving viruses and other pathogens can have an immense impact on human evolution as natural selection acts to increase the prevalence of genetic variants providing resistance to disease. With the emergence of large datasets of human genetic variation, we can search for signatures of natural selection in the human genome driven by such disease-causing microorganisms. Based on this approach, we have previously hypothesized that Lassa virus (LASV) may have been a driver of natural selection in West African populations where Lassa haemorrhagic fever is endemic. In this study, we provide further evidence for this notion. By applying tests for selection to genome-wide data from the International Haplotype Map Consortium and the 1000 Genomes Consortium, we demonstrate evidence for positive selection in LARGE and interleukin 21 (IL21), two genes implicated in LASV infectivity and immunity. We further localized the signals of selection, using the recently developed composite of multiple signals method, to introns and putative regulatory regions of those genes. Our results suggest that natural selection may have targeted variants giving rise to alternative splicing or differential gene expression of LARGE and IL21. Overall, our study supports the hypothesis that selective pressures imposed by LASV may have led to the emergence of particular alleles conferring resistance to Lassa fever, and opens up new avenues of research pursuit.     

TESS3: fast inference of spatial population structure and genome scans for selection

Geography and landscape are important determinants of genetic variation in natural populations, and several ancestry estimation methods have been proposed to investigate population structure using genetic and geographic data simultaneously. Those approaches are often based on computer‐intensive stochastic simulations and do not scale with the dimensions of the data sets generated by high‐throughput sequencing technologies. There is a growing demand for faster algorithms able to analyse genomewide patterns of population genetic variation in their geographic context. In this study, we present TESS3 , a major update of the spatial ancestry estimation program TESS . By combining matrix factorization and spatial statistical methods, TESS3 provides estimates of ancestry coefficients with accuracy comparable to TESS and with run‐times much faster than the Bayesian version. In addition, the TESS3 program can be used to perform genome scans for selection, and separate adaptive from nonadaptive genetic variation using ancestral allele frequency differentiation tests. The main features of TESS3 are illustrated using simulated data and analysing genomic data from European lines of the plant species Arabidopsis thaliana.     

Shared and unique components of human population structure and genome-wide signals of positive selection in South Asia

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.     

Genetics of autism: from genome scans to candidate genes

The autistic disorder was firstly described by Leo Kanner sixty years ago. This complex developmental disability is characterized by social and communicative impairments and repetitive and stereotyped behaviours and interests. The prevalence of autism in the general population is about 1 in 1,000, with four males affected for one female. In approximately 15% of the cases, autism is associated with known genetic disorders, such as fragile X syndrome, tuberous sclerosis or Rett syndrome. Nevertheless, a recognised medical etiology can only be identified in a minority of cases. A higher recurrence risk in families with autistic subjects (45 times greater than the prevalence in the general population) and higher concordance for autism among monozygotic (60-90%) than dizygotic (0-10%) twins argue for a genetic predisposition to idiopathic autism. The past decade has been marked by an increased interest in the genetic basis of autism, with a series of multiple independent whole genome scans and chromosomal abnormalities studies. These analyses have pointed out several candidate regions on chromosomes 2q, 7q, 6q, 15q and sex chromosomes. These regions possess candidate genes that have been screened for mutations or association with autism. However, a clear involvement of a major susceptibility gene (or genes) in autism remains far from clear. The results from linkage studies and the clear drop in the concordance rates between monozygotic and dizygotic twins suggests that the genetic aetiology of autism is certainly heterogeneous (different genes in different families) and polygenic (more than one affected gene per individual). The almost finished sequence of the human genome and the generation of haplotype maps will shed light on the inter-individual genetic variability and will certainly increase the power and reliability of association studies for complex traits, such as autism.