A Positively Selected MAGEE2 LoF Allele Is Associated with Sexual Dimorphism in Human Brain Size and Shows Similar Phenotypes in Magee2 Null Mice.

A nonsense allele at rs1343879 in human MAGEE2 on chromosome X has previously been reported as a strong candidate for positive selection in East Asia. This premature stop codon causing ∼80% protein truncation is characterized by a striking geographical pattern of high population differentiation: common in Asia and the Americas (up to 84% in the 1000 Genomes Project East Asians) but rare elsewhere. Here, we generated a Magee2 mouse knockout mimicking the human loss-of-function mutation to study its functional consequences. The Magee2 null mice did not exhibit gross abnormalities apart from enlarged brain structures (13% increased total brain area, P = 0.0022) in hemizygous males. The area of the granular retrosplenial cortex responsible for memory, navigation, and spatial information processing was the most severely affected, exhibiting an enlargement of 34% (P = 3.4×10-6). The brain size in homozygous females showed the opposite trend of reduced brain size, although this did not reach statistical significance. With these insights, we performed human association analyses between brain size measurements and rs1343879 genotypes in 141 Chinese volunteers with brain MRI scans, replicating the sexual dimorphism seen in the knockout mouse model. The derived stop gain allele was significantly associated with a larger volume of gray matter in males (P = 0.00094), and smaller volumes of gray (P = 0.00021) and white (P = 0.0015) matter in females. It is unclear whether or not the observed neuroanatomical phenotypes affect behavior or cognition, but it might have been the driving force underlying the positive selection in humans.

Genomic analyses inform on migration events during the peopling of Eurasia.

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago.

Prioritising positively selected variants in whole-genome sequencing data using FineMAV.

BACKGROUND: In population genomics, polymorphisms that are highly differentiated between geographically separated populations are often suggestive of Darwinian positive selection. Genomic scans have highlighted several such regions in African and non-African populations, but only a handful of these have functional data that clearly associates candidate variations driving the selection process. Fine-Mapping of Adaptive Variation (FineMAV) was developed to address this in a high-throughput manner using population based whole-genome sequences generated by the 1000 Genomes Project. It pinpoints positively selected genetic variants in sequencing data by prioritizing high frequency, population-specific and functional derived alleles. RESULTS: We developed a stand-alone software that implements the FineMAV statistic. To graphically visualise the FineMAV scores, it outputs the statistics as bigWig files, which is a common file format supported by many genome browsers. It is available as a command-line and graphical user interface. The software was tested by replicating the FineMAV scores obtained using 1000 Genomes Project African, European, East and South Asian populations and subsequently applied to whole-genome sequencing datasets from Singapore and China to highlight population specific variants that can be subsequently modelled. The software tool is publicly available at https://github.com/fadilla-wahyudi/finemav . CONCLUSIONS: The software tool described here determines genome-wide FineMAV scores, using low or high-coverage whole-genome sequencing datasets, that can be used to prioritize a list of population specific, highly differentiated candidate variants for in vitro or in vivo functional screens. The tool displays these scores on the human genome browsers for easy visualisation, annotation and comparison between different genomic regions in worldwide human populations.

Evolutionary and functional analysis of RBMY1 gene copy number variation on the human Y chromosome.

Human RBMY1 genes are located in four variable-sized clusters on the Y chromosome, expressed in male germ cells and possibly associated with sperm motility. We have re-investigated the mutational background and evolutionary history of the RBMY1 copy number distribution in worldwide samples and its relevance to sperm parameters in an Estonian cohort of idiopathic male factor infertility subjects. We estimated approximate RBMY1 copy numbers in 1218 1000 Genomes Project phase 3 males from sequencing read-depth, then chose 14 for valid ation by multicolour fibre-FISH. These fibre-FISH samples provided accurate calibration standards for the entire panel and led to detailed insights into population variation and mutational mechanisms. RBMY1 copy number worldwide ranged from 3 to 13 with a mode of 8. The two larger proximal clusters were the most variable, and additional duplications, deletions and inversions were detected. Placing the copy number estimates onto the published Y-SNP-based phylogeny of the same samples suggested a minimum of 562 mutational changes, translating to a mutation rate of 2.20 × 10-3 (95% CI 1.94 × 10-3 to 2.48 × 10-3) per father-to-son Y-transmission, higher than many short tandem repeat (Y-STRs), and showed no evidence for selection for increased or decreased copy number, but possible copy number stabilizing selection. An analysis of RBMY1 copy numbers among 376 infertility subjects failed to replicate a previously reported association with sperm motility and showed no significant effect on sperm count and concentration, serum follicle stimulating hormone (FSH), luteinizing hormone (LH) and testosterone levels or testicular and semen volume. These results provide the first in-depth insights into the structural rearrangements underlying RBMY1 copy number variation across diverse human lineages.

Transcriptome analysis of Escherichia coli K1 after therapy with hesperidin conjugated with silver nanoparticles.

BACKGROUNDS: Escherichia coli K1 causes neonatal meningitis. Transcriptome studies are indispensable to comprehend the pathology and biology of these bacteria. Recently, we showed that nanoparticles loaded with Hesperidin are potential novel antibacterial agents against E. coli K1. Here, bacteria were treated with and without Hesperidin conjugated with silver nanoparticles, and silver alone, and 50% minimum inhibitory concentration was determined. Differential gene expression analysis using RNA-seq, was performed using Degust software and a set of genes involved in cell stress response and metabolism were selected for the study. RESULTS: 50% minimum inhibitory concentration with silver-conjugated Hesperidin was achieved with 0.5 µg/ml of Hesperidin conjugated with silver nanoparticles at 1 h. Differential genetic analysis revealed the expression of 122 genes (≥ 2-log FC, P< 0.01) in both E. coli K1 treated with Hesperidin conjugated silver nanoparticles and E. coli K1 treated with silver alone, compared to untreated E. coli K1. Of note, the expression levels of cation efflux genes (cusA and copA) and translocation of ions, across the membrane genes (rsxB) were found to increase 2.6, 3.1, and 3.3- log FC, respectively. Significant regulation was observed for metabolic genes and several genes involved in the coordination of flagella. CONCLUSIONS: The antibacterial mechanism of nanoparticles maybe due to disruption of the cell membrane, oxidative stress, and metabolism in E. coli K1. Further studies will lead to a better understanding of the genetic mechanisms underlying treatment with nanoparticles and identification of much needed novel antimicrobial drug candidates.

Mitochondrial DNA Profiling Reveals Two Lineages of Sun Bears in East and West Malaysia.

Sun bear populations are fragmented and at risk from habitat loss and exploitation for body parts. These threats are made worse by significant gaps in knowledge of sun bear population genetic diversity, population connectivity, and taxonomically significant management units. Using a complete sun bear mitochondrial genome, we developed a set of mitochondrial markers to assess haplotype variation and the evolutionary history of sun bears from Peninsular (West) Malaysia and Sabah (East Malaysia). Genetic samples from 28 sun bears from Peninsular Malaysia, 36 from Sabah, and 18 from Thailand were amplified with primers targeting a 1800 bp region of the mitochondrial genome including the complete mitochondrial control region and adjacent genes. Sequences were analyzed using phylogenetic methods. We identified 51 mitochondrial haplotypes among 82 sun bears. Phylogenetic and network analyses provided strong support for a deep split between Malaysian sun bears and sun bears in East Thailand and Yunnan province in China. The Malaysian lineage was further subdivided into two clades: Peninsular Malaysian and Malaysian Borneo (Sabah). Sun bears from Thailand occurred in both Sabah and Peninsular Malaysian clades. Our study supports recent findings that sun bears from Sundaland form a distinct clade from those in China and Indochina with Thailand possessing lineages from the three clades. Importantly, we demonstrate a more recent and clear genetic delineation between sun bears from the Malay Peninsula and Sabah indicating historical barriers to gene flow within the Sundaic region.

Investigation of culturable human gut mycobiota from the segamat community in Johor, Malaysia.

Although several studies have already been carried out in investigating the general profile of the gut mycobiome across several countries, there has yet to be an officially established baseline of a healthy human gut mycobiome, to the best of our knowledge. Microbial composition within the gastrointestinal tract differ across individuals worldwide, and most human gut fungi studies concentrate specifically on individuals from developed countries or diseased cohorts. The present study is the first culture-dependent community study assessing the prevalence and diversity of gut fungi among different ethnic groups from South East Asia. Samples were obtained from a multi-ethnic semi-rural community from Segamat in southern Malaysia. Faecal samples were screened for culturable fungi and questionnaire data analysis was performed. Culturable fungi were present in 45% of the participants' stool samples. Ethnicity had an impact on fungal prevalence and density in stool samples. The prevalence of resistance to fluconazole, itraconazole, voriconazole and 5-fluorocytosine, from the Segamat community, were 14%, 14%, 11% and 7% respectively. It was found that Jakun individuals had lower levels of antifungal resistance irrespective of the drug tested, and male participants had more fluconazole resistant yeast in their stool samples. Two novel point mutations were identified in the ERG11 gene from one azole resistant Candida glabrata, suggesting a possible cause of the occurrence of antifungal resistant isolates in the participant's faecal sample.

Monitoring of heavy metal pollution in urban and rural environments across Pakistan using House crows (Corvus splendens) as bioindicator.

A widely distributed urban bird, the house crow (Corvus splendens), was used to assess bioavailable heavy metals in urban and rural environments across Pakistan. Bioaccumulation of arsenic (As), zinc (Zn), lead (Pb), cadmium (Cd), nickel (Ni), iron (Fe), manganese (Mn), chromium (Cr), and copper (Cu) was investigated in wing feathers of 96 crows collected from eight locations and categorized into four groups pertaining to their geographical and environmental similarities. Results revealed that the concentrations of Pb, Ni, Mn, Cu, and Cr were positively correlated and varied significantly among the four groups. Zn, Fe, Cr, and Cu regarded as industrial outputs, were observed in birds both in industrialized cities and in adjoining rural agricultural areas irrigated through the Indus Basin Irrigation System. Birds in both urban regions accrued Pb more than the metal toxicity thresholds for birds. The house crow was ranked in the middle on the metal accumulation levels in feathers between highly accumulating raptor and piscivore and less contaminated insectivore and granivore species in the studied areas,. This study suggests that the house crow is an efficient bioindicator and supports the feasibility of using feathers to discriminate the local pollution differences among terrestrial environments having different levels and kinds of anthropogenic activities.

The Genetic Legacy of the Indian Ocean Slave Trade: Recent Admixture and Post-admixture Selection in the Makranis of Pakistan.

From the eighth century onward, the Indian Ocean was the scene of extensive trade of sub-Saharan African slaves via sea routes controlled by Muslim Arab and Swahili traders. Several populations in present-day Pakistan and India are thought to be the descendants of such slaves, yet their history of admixture and natural selection remains largely undefined. Here, we studied the genome-wide diversity of the African-descent Makranis, who reside on the Arabian Sea coast of Pakistan, as well that of four neighboring Pakistani populations, to investigate the genetic legacy, population dynamics, and tempo of the Indian Ocean slave trade. We show that the Makranis are the result of an admixture event between local Baluch tribes and Bantu-speaking populations from eastern or southeastern Africa; we dated this event to ∼300 years ago during the Omani Empire domination. Levels of parental relatedness, measured through runs of homozygosity, were found to be similar across Pakistani populations, suggesting that the Makranis rapidly adopted the traditional practice of endogamous marriages. Finally, we searched for signatures of post-admixture selection at traits evolving under positive selection, including skin color, lactase persistence, and resistance to malaria. We demonstrate that the African-specific Duffy-null blood group-believed to confer resistance against Plasmodium vivax infection-was recently introduced to Pakistan through the slave trade and evolved adaptively in this P. vivax malaria-endemic region. Our study reconstructs the genetic and adaptive history of a neglected episode of the African Diaspora and illustrates the impact of recent admixture on the diffusion of adaptive traits across human populations.

Cytosine methylation of rice mitochondrial DNA from grain and leaf tissues.

MAIN CONCLUSION: The rice leaf mitochondrial DNA is  more methylated compared to the rice grain mitochondrial DNA. The old rice leaf mitochondrial DNA has also a higher methylation level than the young rice leaf mitochondrial DNA. The presence of DNA methylation in rice organelles has not been well characterized. We have previously shown that cytosine methylation of chloroplast DNA is different between leaf and grain, and varies between young and old leaves in rice. However, the variation in cytosine methylation of mitochondrial DNA is still poorly characterized. In this study, we have investigated cytosine methylation of mitochondrial DNA in the rice grain and leaf. Based on CpG, CHG, and CHH methylation analyses, the leaf mitochondrial DNA was found to be  more methylated compared to the grain mitochondrial DNA. The methylation of the leaf mitochondrial DNA was also higher in old compared to young leaves. Differences in methylation were observed at different cytosine positions of the mitochondrial DNA between grain and leaf, although there were also positions with a similar level of high methylation in all the tissues examined. The differentially methylated cytosine positions in rice mitochondrial DNA were observed mostly in the intergenic region and in some mitochondrial-specific genes involved in ATP production, transcription, and translation. The functional importance of cytosine methylation in the life cycle of rice mitochondria is still to be determined.

Selected neuropeptide genes show genetic differentiation between Africans and non-Africans.

BACKGROUND: Publicly available genome data provides valuable information on the genetic variation patterns across different modern human populations. Neuropeptide genes are crucial to the nervous, immune, endocrine system, and physiological homeostasis as they play an essential role in communicating information in neuronal functions. It remains unclear how evolutionary forces, such as natural selection and random genetic drift, have affected neuropeptide genes among human populations. To date, there are over 100 known human neuropeptides from the over 1000 predicted peptides encoded in the genome. The purpose of this study is to analyze and explore the genetic variation in continental human populations across all known neuropeptide genes by examining highly differentiated SNPs between African and non-African populations. RESULTS: We identified a total of 644,225 SNPs in 131 neuropeptide genes in 6 worldwide population groups from a public database. Of these, 5163 SNPs that had ΔDAF |(African - non-African)| ≥ 0.20 were identified and fully annotated. A total of 20 outlier SNPs that included 19 missense SNPs with a moderate impact and one stop lost SNP with high impact, were identified in 16 neuropeptide genes. Our results indicate that an overall strong population differentiation was observed in the non-African populations that had a higher derived allele frequency for 15/20 of those SNPs. Highly differentiated SNPs in four genes were particularly striking: NPPA (rs5065) with high impact stop lost variant; CHGB (rs6085324, rs236150, rs236152, rs742710 and rs742711) with multiple moderate impact missense variants; IGF2 (rs10770125) and INS (rs3842753) with moderate impact missense variants that are in linkage disequilibrium. Phenotype and disease associations of these differentiated SNPs indicated their association with hypertension and diabetes and highlighted the pleiotropic effects of these neuropeptides and their role in maintaining physiological homeostasis in humans. CONCLUSIONS: We compiled a list of 131 human neuropeptide genes from multiple databases and literature survey. We detect significant population differentiation in the derived allele frequencies of variants in several neuropeptide genes in African and non-African populations. The results highlights SNPs in these genes that may also contribute to population disparities in prevalence of diseases such as hypertension and diabetes.

Naegleria fowleri: differential genetic expression following treatment with Hesperidin conjugated with silver nanoparticles using RNA-Seq.

Naegleria fowleri causes a deadly infection known as primary amoebic meningoencephalitis (PAM). To our knowledge, there are very few transcriptome studies conducted on these brain-eating amoebae, despite rise in the number of cases. Although the Naegleria genome has been sequenced, currently, it is not well annotated. Transcriptome level studies are needed to help understand the pathology and biology of this fatal parasitic infection. Recently, we showed that nanoparticles loaded with the flavonoid Hesperidin (HDN) are potential novel antimicrobial agents. N. fowleri trophozoites were treated with and without HDN-conjugated with silver nanoparticles (AgNPs) and silver only, and then, 50% minimum inhibitory concentration (MIC) was determined. The results revealed that the MIC of HDN-conjugated AgNPs was 12.5 microg/mL when treated for 3 h. As no reference genome exists for N. fowleri, de novo RNA transcriptome analysis using RNA-Seq and differential gene expression analysis was performed using the Trinity software. Analysis revealed that more than 2000 genes were differentially expressed in response to N. fowleri treatment with HDN-conjugated AgNPs. Some of the genes were linked to oxidative stress response, DNA repair, cell division, cell signalling and protein synthesis. The downregulated genes were linked with processes such as protein modification, synthesis of aromatic amino acids, when compared with untreated N. fowleri. Further transcriptome studies will lead to understanding of genetic mechanisms of the biology and pathogenesis and/or the identification of much needed drug candidates.

Comparative sequence and methylation analysis of chloroplast and amyloplast genomes from rice.

KEY MESSAGE: Grain amyloplast and leaf chloroplast DNA sequences are identical in rice plants but are differentially methylated. The leaf chloroplast DNA becomes more methylated as the rice plant ages. Rice is an important crop worldwide. Chloroplasts and amyloplasts are critical organelles but the amyloplast genome is poorly studied. We have characterised the sequence and methylation of grain amyloplast DNA and leaf chloroplast DNA in rice. We have also analysed the changes in methylation patterns in the chloroplast DNA as the rice plant ages. Total genomic DNA from grain, old leaf and young leaf tissues were extracted from the Oryza sativa ssp. indica cv. MR219 and sequenced using Illumina Miseq. Sequence variant analysis revealed that the amyloplast and chloroplast DNA of MR219 were identical to each other. However, comparison of CpG and CHG methylation between the identical amyloplast and chloroplast DNA sequences indicated that the chloroplast DNA from rice leaves collected at early ripening stage was more methylated than the amyloplast DNA from the grains of the same plant. The chloroplast DNA became more methylated as the plant ages so that chloroplast DNA from young leaves was less methylated overall than amyloplast DNA. These differential methylation patterns were primarily observed in organelle-encoded genes related to photosynthesis followed by those involved in transcription and translation.

Demographic History and Genetic Adaptation in the Himalayan Region Inferred from Genome-Wide SNP Genotypes of 49 Populations.

We genotyped 738 individuals belonging to 49 populations from Nepal, Bhutan, North India, or Tibet at over 500,000 SNPs, and analyzed the genotypes in the context of available worldwide population data in order to investigate the demographic history of the region and the genetic adaptations to the harsh environment. The Himalayan populations resembled other South and East Asians, but in addition displayed their own specific ancestral component and showed strong population structure and genetic drift. We also found evidence for multiple admixture events involving Himalayan populations and South/East Asians between 200 and 2,000 years ago. In comparisons with available ancient genomes, the Himalayans, like other East and South Asian populations, showed similar genetic affinity to Eurasian hunter-gatherers (a 24,000-year-old Upper Palaeolithic Siberian), and the related Bronze Age Yamnaya. The high-altitude Himalayan populations all shared a specific ancestral component, suggesting that genetic adaptation to life at high altitude originated only once in this region and subsequently spread. Combining four approaches to identifying specific positively selected loci, we confirmed that the strongest signals of high-altitude adaptation were located near the Endothelial PAS domain-containing protein 1 and Egl-9 Family Hypoxia Inducible Factor 1 loci, and discovered eight additional robust signals of high-altitude adaptation, five of which have strong biological functional links to such adaptation. In conclusion, the demographic history of Himalayan populations is complex, with strong local differentiation, reflecting both genetic and cultural factors; these populations also display evidence of multiple genetic adaptations to high-altitude environments.

Copy number variation arising from gene conversion on the human Y chromosome.

We describe the variation in copy number of a ~ 10 kb region overlapping the long intergenic noncoding RNA (lincRNA) gene, TTTY22, within the IR3 inverted repeat on the short arm of the human Y chromosome, leading to individuals with 0-3 copies of this region in the general population. Variation of this CNV is common, with 266 individuals having 0 copies, 943 (including the reference sequence) having 1, 23 having 2 copies, and two having 3 copies, and was validated by breakpoint PCR, fibre-FISH, and 10× Genomics Chromium linked-read sequencing in subsets of 1234 individuals from the 1000 Genomes Project. Mapping the changes in copy number to the phylogeny of these Y chromosomes previously established by the Project identified at least 20 mutational events, and investigation of flanking paralogous sequence variants showed that the mutations involved flanking sequences in 18 of these, and could extend over > 30 kb of DNA. While either gene conversion or double crossover between misaligned sister chromatids could formally explain the 0-2 copy events, gene conversion is the more likely mechanism, and these events include the longest non-allelic gene conversion reported thus far. Chromosomes with three copies of this CNV have arisen just once in our data set via another mechanism: duplication of 420 kb that places the third copy 230 kb proximal to the existing proximal copy. Our results establish gene conversion as a previously under-appreciated mechanism of generating copy number changes in humans and reveal the exceptionally large size of the conversion events that can occur.

The Kalash genetic isolate: ancient divergence, drift, and selection.

The Kalash represent an enigmatic isolated population of Indo-European speakers who have been living for centuries in the Hindu Kush mountain ranges of present-day Pakistan. Previous Y chromosome and mitochondrial DNA markers provided no support for their claimed Greek descent following Alexander III of Macedon's invasion of this region, and analysis of autosomal loci provided evidence of a strong genetic bottleneck. To understand their origins and demography further, we genotyped 23 unrelated Kalash samples on the Illumina HumanOmni2.5M-8 BeadChip and sequenced one male individual at high coverage on an Illumina HiSeq 2000. Comparison with published data from ancient hunter-gatherers and European farmers showed that the Kalash share genetic drift with the Paleolithic Siberian hunter-gatherers and might represent an extremely drifted ancient northern Eurasian population that also contributed to European and Near Eastern ancestry. Since the split from other South Asian populations, the Kalash have maintained a low long-term effective population size (2,319-2,603) and experienced no detectable gene flow from their geographic neighbors in Pakistan or from other extant Eurasian populations. The mean time of divergence between the Kalash and other populations currently residing in this region was estimated to be 11,800 (95% confidence interval = 10,600-12,600) years ago, and thus they represent present-day descendants of some of the earliest migrants into the Indian sub-continent from West Asia.

A recent bottleneck of Y chromosome diversity coincides with a global change in culture.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Insights into hominid evolution from the gorilla genome sequence.

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

Revisiting the thrifty gene hypothesis via 65 loci associated with susceptibility to type 2 diabetes.

We have investigated the evidence for positive selection in samples of African, European, and East Asian ancestry at 65 loci associated with susceptibility to type 2 diabetes (T2D) previously identified through genome-wide association studies. Selection early in human evolutionary history is predicted to lead to ancestral risk alleles shared between populations, whereas late selection would result in population-specific signals at derived risk alleles. By using a wide variety of tests based on the site frequency spectrum, haplotype structure, and population differentiation, we found no global signal of enrichment for positive selection when we considered all T2D risk loci collectively. However, in a locus-by-locus analysis, we found nominal evidence for positive selection at 14 of the loci. Selection favored the protective and risk alleles in similar proportions, rather than the risk alleles specifically as predicted by the thrifty gene hypothesis, and may not be related to influence on diabetes. Overall, we conclude that past positive selection has not been a powerful influence driving the prevalence of T2D risk alleles.

A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.