Identifying branch-specific positive selection throughout the regulatory genome using an appropriate proxy neutral.

BACKGROUND: Adaptive changes in cis-regulatory elements are an essential component of evolution by natural selection. Identifying adaptive and functional noncoding DNA elements throughout the genome is therefore crucial for understanding the relationship between phenotype and genotype. RESULTS: We used ENCODE annotations to identify appropriate proxy neutral sequences and demonstrate that the conservativeness of the test can be modulated during the filtration of reference alignments. We applied the method to noncoding Human Accelerated Elements as well as open chromatin elements previously identified in 125 human tissues and cell lines to demonstrate its utility. Then, we evaluated the impact of query region length, proxy neutral sequence length, and branch count on test sensitivity and specificity. We found that the length of the query alignment can vary between 150 bp and 1 kb without affecting the estimation of selection, while for the reference alignment, we found that a length of 3 kb is adequate for proper testing. We also simulated sequence alignments under different classes of evolution and validated our ability to distinguish positive selection from relaxation of constraint and neutral evolution. Finally, we re-confirmed that a quarter of all non-coding Human Accelerated Elements are evolving by positive selection. CONCLUSION: Here, we introduce a method we called adaptiPhy, which adds significant improvements to our earlier method that tests for branch-specific directional selection in noncoding sequences. The motivation for these improvements is to provide a more sensitive and better targeted characterization of directional selection and neutral evolution across the genome.

Gene expression and adaptive noncoding changes during human evolution.

BACKGROUND: Despite evidence for adaptive changes in both gene expression and non-protein-coding, putatively regulatory regions of the genome during human evolution, the relationship between gene expression and adaptive changes in cis-regulatory regions remains unclear. RESULTS: Here we present new measurements of gene expression in five tissues of humans and chimpanzees, and use them to assess this relationship. We then compare our results with previous studies of adaptive noncoding changes, analyzing correlations at the level of gene ontology groups, in order to gain statistical power to detect correlations. CONCLUSIONS: Consistent with previous studies, we find little correlation between gene expression and adaptive noncoding changes at the level of individual genes; however, we do find significant correlations at the level of biological function ontology groups. The types of function include processes regulated by specific transcription factors, responses to genetic or chemical perturbations, and differentiation of cell types within the immune system. Among functional categories co-enriched with both differential expression and noncoding adaptation, prominent themes include cancer, particularly epithelial cancers, and neural development and function.

Spatio-Temporal Variation in Landscape Composition May Speed Resistance Evolution of Pests to Bt Crops.

Transgenic crops that express insecticide genes from Bacillus thuringiensis (Bt) are used worldwide against moth and beetle pests. Because these engineered plants can kill over 95% of susceptible larvae, they can rapidly select for resistance. Here, we use a model for a pyramid two-toxin Bt crop to explore the consequences of spatio-temporal variation in the area of Bt crop and non-Bt refuge habitat. We show that variability over time in the proportion of suitable non-Bt breeding habitat, Q, or in the total area of Bt and suitable non-Bt habitat, K, can increase the overall rate of resistance evolution by causing short-term surges of intense selection. These surges can be exacerbated when temporal variation in Q and/or K cause high larval densities in refuges that increase density-dependent mortality; this will give resistant larvae in Bt fields a relative advantage over susceptible larvae that largely depend on refuges. We address the effects of spatio-temporal variation in a management setting for two bollworm pests of cotton, Helicoverpa armigera and H. punctigera, and field data on landscape crop distributions from Australia. Even a small proportion of Bt fields available to egg-laying females when refuges are sparse may result in high exposure to Bt for just a single generation per year and cause a surge in selection. Therefore, rapid resistance evolution can occur when Bt crops are rare rather than common in the landscape. These results highlight the need to understand spatio-temporal fluctuations in the landscape composition of Bt crops and non-Bt habitats in order to design effective resistance management strategies.

The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network.

Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.

Population genetics of cis-regulatory sequences that operate during embryonic development in the sea urchin Strongylocentrotus purpuratus.

Despite the fact that noncoding sequences comprise a substantial fraction of functional sites within all genomes, the evolutionary mechanisms that operate on genetic variation within regulatory elements remain poorly understood. In this study, we examine the population genetics of the core, upstream cis-regulatory regions of eight genes (AN, CyIIa, CyIIIa, Endo16, FoxB, HE, SM30 a, and SM50) that function during the early development of the purple sea urchin, Strongylocentrotus purpuratus. Quantitative and qualitative measures of segregating variation are not conspicuously different between cis-regulatory and closely linked "proxy neutral" noncoding regions containing no known functional sites. Length and compound mutations are common in noncoding sequences; conventional descriptive statistics ignore such mutations, under-representing true genetic variation by approximately 28% for these loci in this population. Patterns of variation in the cis-regulatory regions of six of the genes examined (CyIIa, CyIIIa, Endo16, FoxB, AN, and HE) are consistent with directional selection. Genetic variation within annotated transcription factor binding sites is comparable to, and frequently greater than, that of surrounding sequences. Comparisons of two paralog pairs (CyIIa/CyIIIa and AN/HE) suggest that distinct evolutionary processes have operated on their cis-regulatory regions following gene duplication. Together, these analyses provide a detailed view of the evolutionary mechanisms operating on noncoding sequences within a natural population, and underscore how little is known about how these processes operate on cis-regulatory sequences.

A potential role for glucose transporters in the evolution of human brain size.

Differences in cognitive abilities and the relatively large brain are among the most striking differences between humans and their closest primate relatives. The energy trade-off hypothesis predicts that a major shift in energy allocation among tissues occurred during human origins in order to support the remarkable expansion of a metabolically expensive brain. However, the molecular basis of this adaptive scenario is unknown. Two glucose transporters (SLC2A1 and SLC2A4) are promising candidates and present intriguing mutations in humans, resulting, respectively, in microcephaly and disruptions in whole-body glucose homeostasis. We compared SLC2A1 and SLC2A4 expression between humans, chimpanzees and macaques, and found compensatory and biologically significant expression changes on the human lineage within cerebral cortex and skeletal muscle, consistent with mediating an energy trade-off. We also show that these two genes are likely to have undergone adaptation and participated in the development and maintenance of a larger brain in the human lineage by modulating brain and skeletal muscle energy allocation. We found that these two genes show human-specific signatures of positive selection on known regulatory elements within their 5'-untranslated region, suggesting an adaptation of their regulation during human origins. This study represents the first case where adaptive, functional and genetic lines of evidence implicate specific genes in the evolution of human brain size.

Whole-genome positive selection and habitat-driven evolution in a shallow and a deep-sea urchin.

Comparisons of genomic sequence between divergent species can provide insight into the action of natural selection across many distinct classes of proteins. Here, we examine the extent of positive selection as a function of tissue-specific and stage-specific gene expression in two closely-related sea urchins, the shallow-water Strongylocentrotus purpuratus and the deep-sea Allocentrotus fragilis, which have diverged greatly in their adult but not larval habitats. Genes that are expressed specifically in adult somatic tissue have significantly higher dN/dS ratios than the genome-wide average, whereas those in larvae are indistinguishable from the genome-wide average. Testis-specific genes have the highest dN/dS values, whereas ovary-specific have the lowest. Branch-site models involving the outgroup S. franciscanus indicate greater selection (ω(FG)) along the A. fragilis branch than along the S. purpuratus branch. The A. fragilis branch also shows a higher proportion of genes under positive selection, including those involved in skeletal development, endocytosis, and sulfur metabolism. Both lineages are approximately equal in enrichment for positive selection of genes involved in immunity, development, and cell-cell communication. The branch-site models further suggest that adult-specific genes have experienced greater positive selection than those expressed in larvae and that ovary-specific genes are more conserved (i.e., experienced greater negative selection) than those expressed specifically in adult somatic tissues and testis. Our results chart the patterns of protein change that have occurred after habitat divergence in these two species and show that the developmental or functional context in which a gene acts can play an important role in how divergent species adapt to new environments.

Contrasts between adaptive coding and noncoding changes during human evolution.

Changes in non-protein-coding regulatory DNA sequences have been proposed to play distinctive roles in adaptive evolution. We analyzed correlations between gene functions and evidence for positive selection in a common statistical framework across several large surveys of coding and noncoding sequences throughout the human genome. Strong correlations with both classifications in gene ontologies and measurements of gene expression indicate that neural development and function have adapted mainly through noncoding changes. In contrast, adaptation via coding changes is dominated by immunity, olfaction, and male reproduction. Genes with highly tissue-specific expression have undergone more adaptive coding changes, suggesting that pleiotropic constraints inhibit such changes in broadly expressed genes. In contrast, adaptive noncoding changes do not exhibit this pattern. Our findings underscore the probable importance of noncoding changes in the evolution of human traits, particularly cognitive traits.

Multiple Functional Variants in cis Modulate PDYN Expression.

Understanding genetic variation and its functional consequences within cis-regulatory regions remains an important challenge in human genetics and evolution. Here, we present a fine-scale functional analysis of segregating variation within the cis-regulatory region of prodynorphin, a gene that encodes an endogenous opioid precursor with roles in cognition and disease. In order to characterize the functional consequences of segregating variation in cis in a region under balancing selection in different human populations, we examined associations between specific polymorphisms and gene expression in vivo and in vitro. We identified five polymorphisms within the 5' flanking region that affect transcript abundance: a 68-bp repeat recognized in prior studies, as well as two microsatellites and two single nucleotide polymorphisms not previously implicated as functional variants. The impact of these variants on transcription differs by brain region, sex, and cell type, implying interactions between cis genotype and the differentiated state of cells. The effects of individual variants on expression level are not additive in some combinations, implying epistatic interactions between nearby variants. These data reveal an unexpectedly complex relationship between segregating genetic variation and its expression-trait consequences and highlights the importance of close functional scrutiny of natural genetic variation within even relatively well-studied cis-regulatory regions.

Functional consequences of genetic variation in primates on tyrosine hydroxylase (TH) expression in vitro.

Tyrosine hydroxylase, the rate-limiting enzyme in catecholamine synthesis, is known to contain naturally occurring genetic variation in it's promoter region that associates with a number of neuropsychological disorders. As such, examining non-coding regions is important for understanding tyrosine hydroxylase function in human health and disease. We examined approximately 2 kb upstream of the translation start site within humans and non-human primates to obtain a fine resolution map of evolutionarily and functionally relevant cis-regulatory differences. Our study investigated Macaca mulatta, Pan troglodytes, Gorilla gorilla, and Homo sapiens haplotypes using transient dual-luciferase transfection in three neuroblastoma cell lines to assay the impact of naturally occurring sequence variation on expression level. In addition to trans effects between cell lines, there are several significant expression differences between primate species, but the most striking difference was seen between human haplotypes in one cell line. Underlying this variation are numerous sequence polymorphisms, two of which influence expression within humans in a non-additive and cell line-specific manner. This study highlights functional consequences of tyrosine hydroxylase genetic variation in primates. Additionally, the results emphasize the importance of examining more than one cell line, the existence of multiple functional variants in a given promoter region and the presence of non-additive cis-interactions.

Genomic features that predict allelic imbalance in humans suggest patterns of constraint on gene expression variation.

Variation in gene expression is an important contributor to phenotypic diversity within and between species. Although this variation often has a genetic component, identification of the genetic variants driving this relationship remains challenging. In particular, measurements of gene expression usually do not reveal whether the genetic basis for any observed variation lies in cis or in trans to the gene, a distinction that has direct relevance to the physical location of the underlying genetic variant, and which may also impact its evolutionary trajectory. Allelic imbalance measurements identify cis-acting genetic effects by assaying the relative contribution of the two alleles of a cis-regulatory region to gene expression within individuals. Identification of patterns that predict commonly imbalanced genes could therefore serve as a useful tool and also shed light on the evolution of cis-regulatory variation itself. Here, we show that sequence motifs, polymorphism levels, and divergence levels around a gene can be used to predict commonly imbalanced genes in a human data set. Reduction of this feature set to four factors revealed that only one factor significantly differentiated between commonly imbalanced and nonimbalanced genes. We demonstrate that these results are consistent between the original data set and a second published data set in humans obtained using different technical and statistical methods. Finally, we show that variation in the single allelic imbalance-associated factor is partially explained by the density of genes in the region of a target gene (allelic imbalance is less probable for genes in gene-dense regions), and, to a lesser extent, the evenness of expression of the gene across tissues and the magnitude of negative selection on putative regulatory regions of the gene. These results suggest that the genomic distribution of functional cis-regulatory variants in the human genome is nonrandom, perhaps due to local differences in evolutionary constraint.

Evolution of incompatibility-inducing microbes in subdivided host populations.

Many insects, other arthropods, and nematodes harbor maternally inherited bacteria inducing "cytoplasmic incompatibility" (CI), reduced egg hatch when infected males mate with uninfected females. Although CI drives the spread of these microbes, selection on alternative, mutually compatible strains in panmictic host populations does not act directly on CI intensity but favors higher "effective fecundity," the number of infected progeny an infected female produces. We analyze the consequences of host population subdivision using deterministic and stochastic models. In subdivided populations, effective fecundity remains the primary target of selection. For strains of equal effective fecundity, if population density is regulated locally (i.e., "soft selection"), variation among patches in infection frequencies may induce change in the relative frequencies of the strains. However, whether this change favors stronger incompatibility depends on initial frequencies. Demographic fluctuations maintain frequency variation that tends to favor stronger incompatibility. However, this effect is weak; even with small patches, minute increases in effective fecundity can offset substantial decreases in CI intensity. These results are insensitive to many details of host life cycle and migration and to systematic outbreeding or inbreeding within patches. Selection acting through transfer between host species may be required to explain the prevalence of CI.

Evolutionary analysis of the cis-regulatory region of the spicule matrix gene SM50 in strongylocentrotid sea urchins.

An evolutionary analysis of transcriptional regulation is essential to understanding the molecular basis of phenotypic diversity. The sea urchin is an ideal system in which to explore the functional consequence of variation in cis-regulatory sequences. We are particularly interested in the evolution of genes involved in the patterning and synthesis of its larval skeleton. This study focuses on the cis-regulatory region of SM50, which has already been characterized to a considerable extent in the purple sea urchin, Strongylocentrotus purpuratus. We have isolated the cis-regulatory region from 15 individuals of S. purpuratus as well as seven closely related species in the family Strongylocentrotidae. We have performed a variety of statistical tests and present evidence that the cis-regulatory elements upstream of the SM50 gene have been subject to positive selection along the lineage leading to S. purpuratus. In addition, we have performed electrophoretic mobility shift assays (EMSAs) and demonstrate that nucleotide substitutions within Element C affect the ability of nuclear proteins to bind to this cis-regulatory element among members of the family Strongylocentrotidae. We speculate that such changes in SM50 and other genes could accumulate to produce altered patterns of gene expression with functional consequences during skeleton formation.

When two is better than one.

Gene duplication and divergence has long been considered an important route to adaptation and phenotypic evolution. Reporting in Nature, Hittinger and Carroll (2007) provide the first clear example of adaptations in both regulatory regions and protein-coding regions after gene duplication. This combination of evolutionary changes appears to have resolved an adaptive conflict, leading to increased organismal fitness.

Migration load and the coexistence of ecologically similar sexuals and asexuals.

The frozen niche variation hypothesis suggests that sexuals can coexist with closely related, ecologically similar asexuals because sexuals and narrowly adapted asexual clones use different resources. However, because a collection of clones can potentially dominate the entire resource axis, such coexistence is not stable. We show that if the sexual population inhabits multiple selection regimes and asexuals are intrinsically slightly less fit than sexuals, migration load in the sexual population allows sexuals and asexuals to coexist stably at the regional level. By decreasing sexuals' fitness, migration load allows asexuals to invade the sexual population. However, as the sexuals' range contracts, migration load decreases, preventing asexuals from driving sexuals to extinction. This "buffering" effect of migration load is even more relevant in models that include more realistic conditions, such as demographic asymmetries or explicit spatial structure.

Promoter regions of many neural- and nutrition-related genes have experienced positive selection during human evolution.

Surveys of protein-coding sequences for evidence of positive selection in humans or chimpanzees have flagged only a few genes known to function in neural or nutritional processes, despite pronounced differences between humans and chimpanzees in behavior, cognition and diet. It may be that most such differences are due to changes in gene regulation rather than protein structure. Here, we present the first survey of promoter (5'-flanking) regions, which are rich in cis-regulatory sequences, for evidence of positive selection in humans. Our results indicate that positive selection has targeted the regulation of many genes known to be involved in neural development and function, both in the brain and elsewhere in the nervous system, and in nutrition, particularly in glucose metabolism.

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Mutation rate and the cost of complexity.

Two recent theoretical studies of adaptation suggest that more complex organisms tend to adapt more slowly. Specifically, in Fisher's "geometric" model of a finite population where multiple traits are under optimizing selection, the average progress ensuing from a single mutation decreases as the number of traits increases--the "cost of complexity." Here, I draw on molecular and histological data to assess the extent to which on a large phylogenetic scale, this predicted decrease in the rate of adaptation per mutation is mitigated by an increase in the number of mutations per generation as complexity increases. As an index of complexity for multicellular organisms, I use the number of visibly distinct types of cell in the body. Mutation rate is the product of mutational target size and population mutation rate per unit target. Despite much scatter, genome size appears to be positively correlated with complexity (as indexed by cell-type number), which along with other considerations suggests that mutational target size tends to increase with complexity. In contrast, effective population mutation rate per unit target appears to be negatively correlated with complexity. The net result is that mutation rate probably does tend to increase with complexity, although probably not fast enough to eliminate the cost of complexity.

Sexual conflict and protein polymorphism.

Sexual conflict, where male and female reproductive interests differ, is probably widespread and often mediated by male or sperm proteins and female or egg proteins that bind to each other during mating or fertilization. One potential consequence is maintenance of polymorphism in these proteins, which might result in reproductive isolation between sympatric subpopulations. I investigate the conditions for polymorphism maintenance in a series of mathematical models of sexual conflict over mating or fertilization frequency. The models represent a male or sperm ligand and a female or egg receptor, and they differ in whether expression of either protein is haploid or diploid. For diploid expression, the conditions imply that patterns of dominance, which involve neither overdominance nor underdominance, can determine whether polymorphism is maintained. For example, suppose ligand expression is diploid, and consider ligand alleles L1 and L2 in interactions with a given receptor genotype; if L1/L1 males are fitter than L2/L2 males in these interactions, then polymorphism is more likely to be maintained when L1/L2 males more closely resemble L1/L1 males in these interactions. Such fitter-allele dominance might be typical of a ligand or its receptor due to their biochemistry, in which case polymorphism might be typical of the pair.

Consequences of recurrent gene flow from crops to wild relatives.

Concern about gene flow from crops to wild relatives has become widespread with the increasing cultivation of transgenic crops. Possible consequences of such gene flow include genetic assimilation, wherein crop genes replace wild ones, and demographic swamping, wherein hybrids are less fertile than their wild parents, and wild populations shrink. Using mathematical models of a wild population recurrently receiving pollen from a genetically fixed crop, we find that the conditions for genetic assimilation are not stringent, and progress towards replacement can be fast, even for disfavoured crop genes. Demographic swamping and genetic drift relax the conditions for genetic assimilation and speed progress towards replacement. Genetic assimilation can involve thresholds and hysteresis, such that a small increase in immigration can lead to fixation of a disfavoured crop gene that had been maintained at a moderate frequency, even if the increase in immigration is cancelled before the gene fixes. Demographic swamping can give rise to 'migrational meltdown', such that a small increase in immigration can lead to not only fixation of a disfavoured crop gene but also drastic shrinkage of the wild population. These findings suggest that the spread of crop genes in wild populations should be monitored more closely.

Coexistence in MacArthur-style consumer-resource models.

The nature of and conditions for permanent coexistence of consumers and resources are characterized in a family of models that generalize MacArthur's consumer-resource model. The generalization is of the resource dynamics, which need not be of Lotka-Volterra form but are subject only to certain restrictions loose enough to admit many resource dynamics of biological interest. For any such model, (1) if there is an interior equilibrium, then it is globally attracting, else some boundary equilibrium is globally attracting-thus permanent coexistence is coexistence at a globally attracting equilibrium; (2) there is an interior equilibrium if and only if for any species, the equilibrium approached in the absence of that species and the presence of the others is invasible by that species--thus permanent coexistence is equivalent to mutual invasibility; (3) for resources without direct interactions, the conditions for permanent coexistence of the consumers admit an instructive formulation in terms of regression statistics. The significance and limitations of the models and results are discussed.