Evidence Of Recent Positive Selection Research Articles

Genome-wide association study of leprosy in Malawi and Mali.

Leprosy is a chronic infection of the skin and peripheral nerves caused by Mycobacterium leprae. Despite recent improvements in disease control, leprosy remains an important cause of infectious disability globally. Large-scale genetic association studies in Chinese, Vietnamese and Indian populations have identified over 30 susceptibility loci for leprosy. There is a significant burden of leprosy in Africa, however it is uncertain whether the findings of published genetic association studies are generalizable to African populations. To address this, we conducted a genome-wide association study (GWAS) of leprosy in Malawian (327 cases, 436 controls) and Malian (247 cases, 368 controls) individuals. In that analysis, we replicated four risk loci previously reported in China, Vietnam and India; MHC Class I and II, LACC1 and SLC29A3. We further identified a novel leprosy susceptibility locus at 10q24 (rs2015583; combined p = 8.81 × 10-9; OR = 0.51 [95% CI 0.40 - 0.64]). Using publicly-available data we characterise regulatory activity at this locus, identifying ACTR1A as a candidate mediator of leprosy risk. This locus shows evidence of recent positive selection and demonstrates pleiotropy with established risk loci for inflammatory bowel disease and childhood-onset asthma. A shared genetic architecture for leprosy and inflammatory bowel disease has been previously described. We expand on this, strengthening the hypothesis that selection pressure driven by leprosy has shaped the evolution of autoimmune and atopic disease in modern populations. More broadly, our data highlights the importance of defining the genetic architecture of disease across genetically diverse populations, and that disease insights derived from GWAS in one population may not translate to all affected populations.

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities.

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Ancestral genomic contributions to complex traits in contemporary Europeans

The contemporary European genetic makeup formed in the last 8,000 years when local Western Hunter-Gatherers (WHGs) mixed with incoming Anatolian Neolithic farmers and Pontic Steppe pastoralists.1-3 This encounter combined genetic variants with distinct evolutionary histories and, together with new environmental challenges faced by the post-Neolithic Europeans, unlocked novel adaptations.4 Previous studies inferred phenotypes in these source populations, using either a few single loci5-7 or polygenic scores based on genome-wide association studies,8-10 and investigated the strength and timing of natural selection on lactase persistence or height, among others.6,11,12 However, how ancient populations contributed to present-day phenotypic variation is poorly understood. Here, we investigate how the unique tiling of genetic variants inherited from different ancestral components drives the complex traits landscape of contemporary Europeans and quantify selection patterns associated with these components. Using matching individual-level genotype and phenotype data for 27 traits in the Estonian biobank13 and genotype data directly from the ancient source populations, we quantify the contributions from each ancestry to present-day phenotypic variation in each complex trait. We find substantial differences in ancestry for eye and hair color, body mass index, waist/hip circumferences, and their ratio, height, cholesterol levels, caffeine intake, heart rate, and age at menarche. Furthermore, we find evidence for recent positive selection linked to four of these traits and, in addition, sleep patterns and blood pressure. Our results show that these ancient components were differentiated enough to contribute ancestry-specific signatures to the complex trait variability displayed by contemporary Europeans.

Identity-by-descent detection across 487,409 British samples reveals fine scale population structure and ultra-rare variant associations

Detection of Identical-By-Descent (IBD) segments provides a fundamental measure of genetic relatedness and plays a key role in a wide range of analyses. We develop FastSMC, an IBD detection algorithm that combines a fast heuristic search with accurate coalescent-based likelihood calculations. FastSMC enables biobank-scale detection and dating of IBD segments within several thousands of years in the past. We apply FastSMC to 487,409 UK Biobank samples and detect ~214 billion IBD segments transmitted by shared ancestors within the past 1500 years, obtaining a fine-grained picture of genetic relatedness in the UK. Sharing of common ancestors strongly correlates with geographic distance, enabling the use of genomic data to localize a sample’s birth coordinates with a median error of 45 km. We seek evidence of recent positive selection by identifying loci with unusually strong shared ancestry and detect 12 genome-wide significant signals. We devise an IBD-based test for association between phenotype and ultra-rare loss-of-function variation, identifying 29 association signals in 7 blood-related traits.

The impact of genetic adaptation on chimpanzee subspecies differentiation.

Chimpanzees, humans' closest relatives, are in danger of extinction. Aside from direct human impacts such as hunting and habitat destruction, a key threat is transmissible disease. As humans continue to encroach upon their habitats, which shrink in size and grow in density, the risk of inter-population and cross-species viral transmission increases, a point dramatically made in the reverse with the global HIV/AIDS pandemic. Inhabiting central Africa, the four subspecies of chimpanzees differ in demographic history and geographical range, and are likely differentially adapted to their particular local environments. To quantitatively explore genetic adaptation, we investigated the genic enrichment for SNPs highly differentiated between chimpanzee subspecies. Previous analyses of such patterns in human populations exhibited limited evidence of adaptation. In contrast, chimpanzees show evidence of recent positive selection, with differences among subspecies. Specifically, we observe strong evidence of recent selection in eastern chimpanzees, with highly differentiated SNPs being uniquely enriched in genic sites in a way that is expected under recent adaptation but not under neutral evolution or background selection. These sites are enriched for genes involved in immune responses to pathogens, and for genes inferred to differentiate the immune response to infection by simian immunodeficiency virus (SIV) in natural vs. non-natural host species. Conversely, central chimpanzees exhibit an enrichment of signatures of positive selection only at cytokine receptors, due to selective sweeps in CCR3, CCR9 and CXCR6 -paralogs of CCR5 and CXCR4, the two major receptors utilized by HIV to enter human cells. Thus, our results suggest that positive selection has contributed to the genetic and phenotypic differentiation of chimpanzee subspecies, and that viruses likely play a predominate role in this differentiation, with SIV being a likely selective agent. Interestingly, our results suggest that SIV has elicited distinctive adaptive responses in these two chimpanzee subspecies.

An Aptitude for Altitude: Are Epigenomic Processes Involved?

In recent years, high-throughput genomic technologies and computational advancements have invigorated efforts to identify the molecular mechanisms regulating human adaptation to high altitude. Although exceptional progress regarding the identification of genomic regions showing evidence of recent positive selection has been made, many of the key “hypoxia tolerant” phenotypes of highland populations have not yet been linked to putative adaptive genetic variants. As a result, fundamental questions regarding the biological processes by which such adaptations are acquired remain unanswered. This Mini Review discusses the hypothesis that the epigenome works in coordination with underlying genomic sequence to govern adaptation to the chronic hypoxia of high altitude by influencing adaptive capacity and phenotypic variation under conditions of environmental hypoxia. Efforts to unravel the complex interactions between the genome, epigenome, and environmental exposures are essential to more fully appreciate the mechanisms underlying human adaptation to hypoxia.

Signatures of Recent Positive Selection in Enhancers Across 41 Human Tissues.

Evolutionary changes in enhancers are widely associated with variation in human traits and diseases. However, studies comprehensively quantifying levels of selection on enhancers at multiple evolutionary periods during recent human evolution and how enhancer evolution varies across human tissues are lacking. To address these questions, we integrated a dataset of 41,561 transcribed enhancers active in 41 different human tissues (FANTOM Consortium) with whole genome sequences of 1,668 individuals from the African, Asian, and European populations (1000 Genomes Project). Our analyses based on four different metrics (Tajima’s D, FST, H12, nSL) showed that ∼5.90% of enhancers showed evidence of recent positive selection and that genes associated with enhancers under very recent positive selection are enriched for diverse immune-related functions. The distributions of these metrics for brain and testis enhancers were often statistically significantly different and in the direction suggestive of less positive selection compared to those of other tissues; the same was true for brain and testis enhancers that are tissue-specific compared to those that are tissue-broad and for testis enhancers associated with tissue-enriched and non-tissue-enriched genes. These differences varied considerably across metrics and tissues and were generally in the form of changes in distributions’ shapes rather than shifts in their values. Collectively, these results suggest that many human enhancers experienced recent positive selection throughout multiple time periods in human evolutionary history, that this selection occurred in a tissue-dependent and immune-related functional context, and that much like the evolution of their protein-coding gene counterparts, the evolution of brain and testis enhancers has been markedly different from that of enhancers in other tissues.

Adaptation to milking agropastoralism in Chilean goat herders and nutritional benefit of lactase persistence.

The genetic trait of lactase persistence (LP) evolved as an adaptation to milking pastoralism in the Old World and is a well‐known example of positive natural selection in humans. However, the specific mechanisms conferring this selective advantage are unknown. To understand the relationship between milk drinking, LP, growth, reproduction, and survival, communities of the Coquimbo Region in Chile, with recent adoption of milking agropastoralism, were used as a model population.DNA samples and data on stature, reproduction, and diet were collected from 451 participants. Lactose tolerance tests were done on 41 of them. The European −13,910*T (rs4988235) was the only LP causative variant found, showing strong association (99.6%) with LP phenotype.Models of associations of inferred LP status and milk consumption, with fertility, mortality, height, and weight were adjusted with measures of ancestry and relatedness to control for population structure. Although we found no statistically significant effect of LP on fertility, a significant effect (P = 0.002) was observed of LP on body mass index (BMI) in males and of BMI on fertility (P = 0.003). These results fail to support a causal relationship between LP and fertility yet suggest the idea of a nutritional advantage of LP. Furthermore, the proportion of European ancestry around the genetic region of −13,910*T is significantly higher (P = 0.008) than the proportion of European ancestry genome‐wide, providing evidence of recent positive selection since European–Amerindian admixture. This signature was absent in nonpastoralist Latin American populations, supporting the hypothesis of specific adaptation to milking agropastoralism in the Coquimbo communities.

Systematic detection of positive selection in the human-pathogen interactome and lasting effects on infectious disease susceptibility.

Infectious disease has shaped the natural genetic diversity of humans throughout the world. A new approach to capture positive selection driven by pathogens would provide information regarding pathogen exposure in distinct human populations and the constantly evolving arms race between host and disease-causing agents. We created a human pathogen interaction database and used the integrated haplotype score (iHS) to detect recent positive selection in genes that interact with proteins from 26 different pathogens. We used the Human Genome Diversity Panel to identify specific populations harboring pathogen-interacting genes that have undergone positive selection. We found that human genes that interact with 9 pathogen species show evidence of recent positive selection. These pathogens are Yersenia pestis, human immunodeficiency virus (HIV) 1, Zaire ebolavirus, Francisella tularensis, dengue virus, human respiratory syncytial virus, measles virus, Rubella virus, and Bacillus anthracis. For HIV-1, GWAS demonstrate that some naturally selected variants in the host-pathogen protein interaction networks continue to have functional consequences for susceptibility to these pathogens. We show that selected human genes were enriched for HIV susceptibility variants (identified through GWAS), providing further support for the hypothesis that ancient humans were exposed to lentivirus pandemics. Human genes in the Italian, Miao, and Biaka Pygmy populations that interact with Y. pestis show significant signs of selection. These results reveal some of the genetic footprints created by pathogens in the human genome that may have left lasting marks on susceptibility to infectious disease.

Fast and robust detection of ancestral selective sweeps.

There are few methods tailored for detecting signals of positive selection in populations directly ancestral to multiple descendent populations. We introduce the ancestral branch statistic (ABS), a four-population summary statistic for identifying selective sweeps occurring in the direct ancestor of a pair of populations. Simulations show that ABS performs at least as well as, and often better under model violations, than the complementary likelihood approach of 3P-CLR across diverse selection scenarios and parameter values. We first applied ABS to contemporary human genomic data to identify genes that may have been adaptive in ancestral East Asian populations, uncovering the well-established candidate EDAR, as well as a novel candidate SLC35F3, which encodes a putative thiamine transporter that may have been involved in adaptation to eating polished grains. Next, we performed scans with ancient European genomic data to reexamine evidence of recent positive selection in ancestral Europeans. The MCM6/LCT cluster and the SLC45A2 and HERC2 genes are strong outliers, agreeing with previous studies. Novel candidates, such as SLC30A9 and CYP1A2, may have been involved in adaptation to local nutrient sufficiency and lifestyle changes. Finally, we provide open-source software, CalcABS, which can perform genomic scans of ancestral sweeps with ABS from population allele frequency data.

Discovery of Genetic Variation on Chromosome 5q22 Associated with Mortality in Heart Failure.

Failure of the human heart to maintain sufficient output of blood for the demands of the body, heart failure, is a common condition with high mortality even with modern therapeutic alternatives. To identify molecular determinants of mortality in patients with new-onset heart failure, we performed a meta-analysis of genome-wide association studies and follow-up genotyping in independent populations. We identified and replicated an association for a genetic variant on chromosome 5q22 with 36% increased risk of death in subjects with heart failure (rs9885413, P = 2.7x10-9). We provide evidence from reporter gene assays, computational predictions and epigenomic marks that this polymorphism increases activity of an enhancer region active in multiple human tissues. The polymorphism was further reproducibly associated with a DNA methylation signature in whole blood (P = 4.5x10-40) that also associated with allergic sensitization and expression in blood of the cytokine TSLP (P = 1.1x10-4). Knockdown of the transcription factor predicted to bind the enhancer region (NHLH1) in a human cell line (HEK293) expressing NHLH1 resulted in lower TSLP expression. In addition, we observed evidence of recent positive selection acting on the risk allele in populations of African descent. Our findings provide novel genetic leads to factors that influence mortality in patients with heart failure.

Novel genomic signals of recent selection in an Ethiopian population.

The recent feasibility of genome-wide studies of adaptation in human populations has provided novel insights into biological pathways that have been affected by adaptive pressures. However, only a few African populations have been investigated using these genome-wide approaches. Here, we performed a genome-wide analysis for evidence of recent positive selection in a sample of 120 individuals of Wolaita ethnicity belonging to Omotic-speaking people who have inhabited the mid- and high-land areas of southern Ethiopia for millennia. Using the 11 HapMap populations as the comparison group, we found Wolaita-specific signals of recent positive selection in several human leukocyte antigen (HLA) loci. Notably, the selected loci overlapped with HLA regions that we previously reported to be associated with podoconiosis-a geochemical lymphedema of the lower legs common in the Wolaita area. We found selection signals in PPARA, a gene involved in energy metabolism during prolonged food deficiency. This finding is consistent with the dietary use of enset, a crop with high-carbohydrate and low-fat and -protein contents domesticated in Ethiopia subsequent to food deprivation 10 000 years ago, and with metabolic adaptation to high-altitude hypoxia. We observed novel selection signals in CDKAL1 and NEGR1, well-known diabetes and obesity susceptibility genes. Finally, the SLC24A5 gene locus known to be associated with skin pigmentation was in the top selection signals in the Wolaita, and the alleles of single-nucleotide polymorphisms rs1426654 and rs1834640 (SLC24A5) associated with light skin pigmentation in Eurasian populations were of high frequency (47.9%) in this Omotic-speaking indigenous Ethiopian population.

Insights into the Development and Evolution of Exaggerated Traits Using De Novo Transcriptomes of Two Species of Horned Scarab Beetles

Scarab beetles exhibit an astonishing variety of rigid exo-skeletal outgrowths, known as “horns”. These traits are often sexually dimorphic and vary dramatically across species in size, shape, location, and allometry with body size. In many species, the horn exhibits disproportionate growth resulting in an exaggerated allometric relationship with body size, as compared to other traits, such as wings, that grow proportionately with body size. Depending on the species, the smallest males either do not produce a horn at all, or they produce a disproportionately small horn for their body size. While the diversity of horn shapes and their behavioural ecology have been reasonably well studied, we know far less about the proximate mechanisms that regulate horn growth. Thus, using 454 pyrosequencing, we generated transcriptome profiles, during horn growth and development, in two different scarab beetle species: the Asian rhinoceros beetle, Trypoxylus dichotomus, and the dung beetle, Onthophagus nigriventris. We obtained over half a million reads for each species that were assembled into over 6,000 and 16,000 contigs respectively. We combined these data with previously published studies to look for signatures of molecular evolution. We found a small subset of genes with horn-biased expression showing evidence for recent positive selection, as is expected with sexual selection on horn size. We also found evidence of relaxed selection present in genes that demonstrated biased expression between horned and horn-less morphs, consistent with the theory of developmental decoupling of phenotypically plastic traits.

Genes associated with SLE are targets of recent positive selection.

The reasons for the ethnic disparities in the prevalence of systemic lupus erythematosus (SLE) and the relative high frequency of SLE risk alleles in the population are not fully understood. Population genetic factors such as natural selection alter allele frequencies over generations and may help explain the persistence of such common risk variants in the population and the differential risk of SLE. In order to better understand the genetic basis of SLE that might be due to natural selection, a total of 74 genomic regions with compelling evidence for association with SLE were tested for evidence of recent positive selection in the HapMap and HGDP populations, using population differentiation, allele frequency, and haplotype-based tests. Consistent signs of positive selection across different studies and statistical methods were observed at several SLE-associated loci, including PTPN22, TNFSF4, TET3-DGUOK, TNIP1, UHRF1BP1, BLK, and ITGAM genes. This study is the first to evaluate and report that several SLE-associated regions show signs of positive natural selection. These results provide corroborating evidence in support of recent positive selection as one mechanism underlying the elevated population frequency of SLE risk loci and supports future research that integrates signals of natural selection to help identify functional SLE risk alleles.

Editors’ Pick: milk sugar, migration and pastoralism in Africa

Lactose persistence (LP), the ability to digest milk sugar (lactose), is one of the best examples of selection-based evolutionary change in humans from milk-drinking cultures. LP has been documented in populations with northern European ancestry and those originating in Central Asia, Middle East, the Arabian Peninsula, and Africa. Based on molecular genetic studies, it has been estimated that the selective change in the ability to digest lactose beyond childhood occurred within approximately the past 5,000-10,000 years [1]. The estimation is consistent with an advantage to LP in dairy farming-based subsistence, which came to prominence in certain geographical areas around the same time. The important enzyme in the breakdown of milk sugar, lactase (or more specifically lactase-phlorizin hydrolase (LPH)), is coded by the LCT gene. The molecular genetic data reported in recent years show several interesting aspects of LP. The first genetic variation described for European LP was not in the LCT gene, but in the in the intron 13 of the MCM6 gene, residing ~14 kb upstream of MCM6 [2]. However, it soon was realized that the MCM6 intron 13 T-13910 allele, which explained most LP in Europe, was absent in Africa [3,4]. Instead, three other SNPs in the same intron, C-14010, G-13907 and G-13915, were associated with lactase persistence in Africa. Also, these SNPs originated on different haplotype backgrounds from the European T-13910 allele and from each other providing a major model of convergent evolution in LP due to strong selective pressure [3,5]. To date, several SNPs associated with LP are known. In Europe, in addition to MCM6 intron 13 T-13910, intron 9 G-22018 is significantly associated with the LP trait. In Africa and the Arabian Peninsula, four SNPs (those above and G-14009), which are all located within 100 bp of T-13910 have been shown to be associated with the LP predominantly in pastoralist populations. Those variants do not, however, account for all the phenotypic variance of LP in Africa, implying that additional genetic variation may play a role. A new study by Ranciaro et al. in American Journal of Human Genetics[6] sought out new genetic variants associated with LP in Africa. They studied 819 individuals from 63 African populations, and 154 individuals from non-African populations from Europe, the Middle East and Asia. The authors sequenced MCM6 gene introns 9 and 13 and ~2 kb of the LCT gene promoter region in these samples. Additionally, in an effort to reconstruct the origin and spread of LP-associated alleles in Africa, four microsatellites were genotyped in a ~198 kb region in a subset of 252 individuals. The participants’ LP status was tested using a standardized blood test that tracks lactose digestion. With the new larger data set, and sampling from outside Africa, the authors confirmed the earlier observed association between the LP trait in Africa and three variants in intron 13, and they also found two additional LP-associated SNPs, one in intron 13 (G-12962) and one in the LCT promoter region (T-956). Further functional studies are naturally needed to prove whether these new SNPs are causative or just merely in linkage disequilibrium with the earlier identified variants. Ranciaro et al. [6] then used allele frequency and long-range linkage disequilibrium based neutrality tests to detect any signatures of selection. They found evidence for recent positive selection in eastern African populations and the Fulani from central Africa. While looking for the new variants for LP within Africa, the study by Ranciaro et al. [6] also considered that if LP patterns coincide with the rise in pastoralism, the LP-related alleles found in present-day populations should reveal past migrations of those populations. In order to do that, they analyzed a subset of 252 study participants using four microsatellites and combined these with the SNPs to create haplotypes. The haplotype analysis has some interesting insights about the migration of pastoralism groups in Africa. For example, the haplotype analysis supported an eastern African origin (Kenya and Tanzania) of the C-14010 LP-associated mutation now found in southern Africa in Bantu-speaking Xhosa and the hunter-gatherer San population. Also, the presence of Middle Eastern origin G-13915 allele in North African populations suggests migration between these regions. Interestingly, the data are supported by known historical interactions between these populations. Some of the data in Ranciaro et al. [6] is in agreement with earlier findings of distinct genetic LP variants that have arisen in various parts of Africa through convergent evolution. For example, the G-13907 was almost exclusively restricted to northern Kenya, Ethiopia and northern Sudan, as reported earlier. Last but not least, the authors also detected the traditionally European LP haplotype background in some West-and North-Central African populations. But, even with these new data from Ranciaro et al. [6], the genetic variants found and analyzed in various populations to date, are far from fully explaining LP in Africa.

Runs of homozygosity associated with speech delay in autism in a taiwanese han population: evidence for the recessive model.

Runs of homozygosity (ROH) may play a role in complex diseases. In the current study, we aimed to test if ROHs are linked to the risk of autism and related language impairment. We analyzed 546,080 SNPs in 315 Han Chinese affected with autism and 1,115 controls. ROH was defined as an extended homozygous haplotype spanning at least 500 kb. Relative extended haplotype homozygosity (REHH) for the trait-associated ROH region was calculated to search for the signature of selection sweeps. Totally, we identified 676 ROH regions. An ROH region on 11q22.3 was significantly associated with speech delay (corrected p = 1.73×10−8). This region contains the NPAT and ATM genes associated with ataxia telangiectasia characterized by language impairment; the CUL5 (culin 5) gene in the same region may modulate the neuronal migration process related to language functions. These three genes are highly expressed in the cerebellum. No evidence for recent positive selection was detected on the core haplotypes in this region. The same ROH region was also nominally significantly associated with speech delay in another independent sample (p = 0.037; combinatorial analysis Stouffer’s z trend = 0.0005). Taken together, our findings suggest that extended recessive loci on 11q22.3 may play a role in language impairment in autism. More research is warranted to investigate if these genes influence speech pathology by perturbing cerebellar functions.

Selection of Phototransduction Genes inHomo sapiens

We investigated the evidence of recent positive selection in the human phototransduction system at single nucleotide polymorphism (SNP) and gene level. SNP genotyping data from the International HapMap Project for European, Eastern Asian, and African populations was used to discover differences in haplotype length and allele frequency between these populations. Numeric selection metrics were computed for each SNP and aggregated into gene-level metrics to measure evidence of recent positive selection. The level of recent positive selection in phototransduction genes was evaluated and compared to a set of genes shown previously to be under recent selection, and a set of highly conserved genes as positive and negative controls, respectively. Six of 20 phototransduction genes evaluated had gene-level selection metrics above the 90th percentile: RGS9, GNB1, RHO, PDE6G, GNAT1, and SLC24A1. The selection signal across these genes was found to be of similar magnitude to the positive control genes and much greater than the negative control genes. There is evidence for selective pressure in the genes involved in retinal phototransduction, and traces of this selective pressure can be demonstrated using SNP-level and gene-level metrics of allelic variation. We hypothesize that the selective pressure on these genes was related to their role in low light vision and retinal adaptation to ambient light changes. Uncovering the underlying genetics of evolutionary adaptations in phototransduction not only allows greater understanding of vision and visual diseases, but also the development of patient-specific diagnostic and intervention strategies.

Detecting Loci under Recent Positive Selection in Dairy and Beef Cattle by Combining Different Genome-Wide Scan Methods

As the methodologies available for the detection of positive selection from genomic data vary in terms of assumptions and execution, weak correlations are expected among them. However, if there is any given signal that is consistently supported across different methodologies, it is strong evidence that the locus has been under past selection. In this paper, a straightforward frequentist approach based on the Stouffer Method to combine P-values across different tests for evidence of recent positive selection in common variations, as well as strategies for extracting biological information from the detected signals, were described and applied to high density single nucleotide polymorphism (SNP) data generated from dairy and beef cattle (taurine and indicine). The ancestral Bovinae allele state of over 440,000 SNP is also reported. Using this combination of methods, highly significant (P<3.17×10−7) population-specific sweeps pointing out to candidate genes and pathways that may be involved in beef and dairy production were identified. The most significant signal was found in the Cornichon homolog 3 gene (CNIH3) in Brown Swiss (P = 3.82×10−12), and may be involved in the regulation of pre-ovulatory luteinizing hormone surge. Other putative pathways under selection are the glucolysis/gluconeogenesis, transcription machinery and chemokine/cytokine activity in Angus; calpain-calpastatin system and ribosome biogenesis in Brown Swiss; and gangliosides deposition in milk fat globules in Gyr. The composite method, combined with the strategies applied to retrieve functional information, may be a useful tool for surveying genome-wide selective sweeps and providing insights in to the source of selection.

Extended homozygous haplotypes at genes involved in brain development associated with autism

Background Recent evidence has suggested that extended homozygous haplotypes (EHH) in several genomic regions may be associated with risk of psychiatric disorders (for example, schizophrenia). The phenomenon of EHHs may arise from recent positive selection, inbreeding, as well as recessive models. Materials and methods We used Affymetrix 500K SNP arrays to search for EHH in 1,385 affected individuals and 1,498 unaffected individuals. The EHH was defined as at least 100 contiguous homozygous SNPs. To interrogate the associations between EHHs and autism, logistic regression analysis with generalized equation estimation model to adjust for intrafamily correlation was performed. We also examined if any of the associated EHHs were related to deletions by examining the data of copy number variants. Finally, we checked if any associated EHHs contained genes with signatures of recent positive selection in the Hapmap sample. We used the Sidak method to correct for multiple tests. Results The best finding was obtained at the HMLGCL1 gene (P =3×1 0 –5 , odds ratio 0.15). The HMLGCL1 gene has been found to be highly expressed in some brain regions. Other genes harbored in regions enriched with EHHs associated with autism include ZFP91, CNTF, NAPL1 and TLE4. We also used the webtool Panther to assess if this set of genes is over-represented in any pathways. The results suggest that these genes collectively may be involved in the Wnt signaling pathway and neurodevelopment process. No remarkable evidence for recent positive selection was obtained for most of these loci, except the TLE4 gene (P = 0.008). Conclusions Taken together, these genes involved in brain development may harbor variants that exert a recessive effect on the risk of autism. Further work is warranted to replicate these findings in other populations.

Crohn's Disease and Genetic Hitchhiking at IBD5

Inflammatory bowel disease 5 (IBD5) is a 250 kb haplotype on chromosome 5 that is associated with an increased risk of Crohn's disease in Europeans. The OCTN1 gene is centrally located on IBD5 and encodes a transporter of the antioxidant ergothioneine (ET). The 503F variant of OCTN1 is strongly associated with IBD5 and is a gain-of-function mutation that increases absorption of ET. Although 503F has been implicated as the variant potentially responsible for Crohn's disease susceptibility at IBD5, there is little evidence beyond statistical association to support its role in disease causation. We hypothesize that 503F is a recent adaptation in Europeans that swept to relatively high frequency and that disease association at IBD5 results not from 503F itself, but from one or more nearby hitchhiking variants, in the genes IRF1 or IL5. To test for evidence of recent positive selection on the 503F allele, we employed the iHS statistic, which was significant in the European CEU HapMap population (P=0.0007) and European Human Genome Diversity Panel populations (P≤0.01). To evaluate the hypothesis of disease-variant hitchhiking, we performed haplotype association tests on high-density microarray data in a sample of 1,868 Crohn's disease cases and 5,550 controls. We found that 503F haplotypes with recombination breakpoints between OCTN1 and IRF1 or IL5 were not associated with disease (odds ratio [OR]: 1.05, P=0.21). In contrast, we observed strong disease association for 503F haplotypes with no recombination between these three genes (OR: 1.24, P=2.6×10(-8)), as expected if the sweeping haplotype harbored one or more disease-causing mutations in IRF1 or IL5. To further evaluate these disease-gene candidates, we obtained expression data from lower gastrointestinal biopsies of healthy individuals and Crohn's disease patients. We observed a 72% increase in gene expression of IRF1 among Crohn's disease patients (P=0.0006) and no significant difference in expression of OCTN1. Collectively, these data indicate that the 503F variant has increased in frequency due to recent positive selection and that disease-causing variants in linkage disequilibrium with 503F have hitchhiked to relatively high frequency, thus forming the IBD5 risk haplotype. Finally, our association results and expression data support IRF1 as a strong candidate for Crohn's disease causation.