Population genetic simulation: Benchmarking frameworks for non-standard models of natural selection.

Population genetic simulation has emerged as a common tool for investigating increasingly complex evolutionary and demographic models. Software capable of handling high-level model complexity has recently been developed, and the advancement of tree sequence recording now allows simulations to merge the efficiency and genealogical insight of coalescent simulations with the flexibility of forward simulations. However, frameworks utilizing these features have not yet been compared and benchmarked. Here, we evaluate various simulation workflows using the coalescent simulator msprime and the forward simulator SLiM, to assess resource efficiency and determine an optimal simulation framework. Three aspects were evaluated: (1) the burn-in, to establish an equilibrium level of neutral diversity in the population; (2) the forward simulation, in which temporally fluctuating selection is acting; and (3) the final computation of summary statistics. We provide typical memory and computation time requirements for each step. We find that the fastest framework, a combination of coalescent and forward simulation with tree sequence recording, increases simulation speed by over twenty times compared to classical forward simulations without tree sequence recording, although it does require six times more memory. Overall, using efficient simulation workflows can lead to a substantial improvement when modelling complex evolutionary scenarios-although the optimal framework ultimately depends on the available computational resources.

Genetic Risk Assessment of Degenerative Eye Disease (GRADE): study protocol of a prospective assessment of polygenic risk scores to predict diagnosis of glaucoma and age-related macular degeneration.

BACKGROUND: Glaucoma and age-related macular degeneration (AMD) account for a substantial portion of global blindness. Both conditions are highly heritable, with recognised monogenic and polygenic inheritance patterns. Current screening guidelines lack decisive recommendations. Polygenic risk scores (PRS) allow for cost-effective broad population risk stratification for these conditions. The predictive potential of PRS could facilitate earlier diagnosis and treatment, and prevent unnecessary vision loss. METHODS: The Genetic Risk Assessment of Degenerative Eye disease (GRADE) study is a prospective study designed to generate high-quality evidence about the feasibility of PRS to stratify individuals from the general population, enabling identification of those at highest risk of developing glaucoma or AMD. The targeted recruitment is 1000 individuals aged over 50 years, from which blood or saliva samples will be used for genotyping and an individual PRS for glaucoma and AMD will be derived. Individuals with PRS values in the bottom decile (n = 100), top decile (n = 100) and middle 80% (n = 100) for both glaucoma and AMD will undergo a detailed eye examination for glaucoma and/or AMD. DISCUSSION: The primary objective will be to compare the prevalence of glaucoma and AMD cases between low, intermediate, and high PRS risk groups. We expect to find a higher prevalence of both diseases in the high PRS risk group, as compared to the middle and low risk groups. This prospective study will assess the clinical validity of a PRS for glaucoma and AMD in the general Australian population. Positive findings will support the implementation of PRS into clinical practice.

No Strong Association between the Apolipoprotein E E4 Allele and Glaucoma: A Multicohort Study.

Purpose: To elucidate a potential association between the apolipoprotein E (APOE) E4 allele and glaucoma prevalence in large cohorts. Design: A cross-sectional analysis of baseline and prospectively collected cohort data. Participants: UK Biobank (UKBB) participants of genetically determined European ancestry (n = 438 711). Replication analyses were performed using clinical and genotyping data collected from European participants recruited to the Canadian Longitudinal Study of Aging (CLSA; n = 18 199), the Australian and New Zealand Registry of Advanced Glaucoma (ANZRAG; n = 1970), and the Blue Mountains Eye Study (BMES; n = 2440). Methods: Apolipoprotein E alleles and genotypes were determined, and their distributions were compared on the basis of glaucoma status. Similar analyses were performed using positive control outcomes associated with the APOE E4 allele (death, dementia, age-related macular degeneration) and negative control outcomes not associated with the APOE E4 allele (cataract, diabetic eye disease). Outcome phenotypes were also correlated with Alzheimer's dementia (AD), a clinical outcome highly associated with the APOE E4 allele. Main Outcome Measures: Results of APOE E4 genotype-phenotype comparisons were reported as odds ratios (ORs) with 95% confidence intervals (CIs). Replication analyses investigated APOE E4 associations in 2 replication cohorts (CLSA and ANZRAG/BMES). Results: The APOE E4 allele was inversely associated with glaucoma (OR, 0.96; 95% CI, 0.93-0.99; P = 0.016) and both negative controls (cataract: OR, 0.98; 95% CI, 0.96-0.99; P = 0.015; diabetic eye disease: OR, 0.92; 95% CI, 0.87-0.97; P = 0.003) in the UKBB cohort. A paradoxical positive association was observed between AD and both glaucoma (OR, 1.30; 95% CI, 1.08-1.54; P < 0.01) and cataract (OR, 1.15; 1.04-1.28; P = 0.018). No association between the APOE E4 allele and glaucoma was observed in either replication cohort (CLSA: OR, 1.03; 95% CI, 0.89-1.19; P = 0.66; ANZRAG/BMES: OR, 0.97; 95% CI, 0.84-1.12; P = 0.65). Conclusions: A small negative association observed between APOE E4 and glaucoma within the UKBB was not evident in either replication cohort and may represent an artifact of glaucoma underdiagnosis in APOE E4 carriers. Financial Disclosures: The author(s) have no proprietary or commercial interest in any materials discussed in this article.

Fluctuating selection and the determinants of genetic variation.

Recent studies of cosmopolitan Drosophila populations have found hundreds to thousands of genetic loci with seasonally fluctuating allele frequencies, bringing temporally fluctuating selection to the forefront of the historical debate surrounding the maintenance of genetic variation in natural populations. Numerous mechanisms have been explored in this longstanding area of research, but these exciting empirical findings have prompted several recent theoretical and experimental studies that seek to better understand the drivers, dynamics, and genome-wide influence of fluctuating selection. In this review, we evaluate the latest evidence for multilocus fluctuating selection in Drosophila and other taxa, highlighting the role of potential genetic and ecological mechanisms in maintaining these loci and their impacts on neutral genetic variation.

Mitochondrial TXNRD2 and ME3 Genetic Risk Scores Are Associated with Specific Primary Open-Angle Glaucoma Phenotypes.

PURPOSE: Genetic variants in regions that include the mitochondrial genes thioredoxin reductase 2 (TXNRD2) and malic enzyme 3 (ME3) are associated with primary open-angle glaucoma (POAG) in genome-wide association studies (GWASs). To assess their clinical impact, we investigated whether TXNRD2 and ME3 genetic risk scores (GRSs) are associated with specific glaucoma phenotypes. DESIGN: Cross-sectional study. PARTICIPANTS: A total of 2617 patients with POAG and 2634 control participants from the National Eye Institute Glaucoma Human Genetics Collaboration Hereditable Overall Operational Database (NEIGHBORHOOD) consortium. METHODS: All POAG-associated single nucleotide polymorphisms (SNPs) in the TXNRD2 and ME3 loci were identified using GWAS data (P < 0.05). Of these, 20 TXNRD2 and 24 ME3 SNPs were selected after adjusting for linkage disequilibrium. The correlation between SNP effect size and gene expression levels was investigated using the Gene-Tissue Expression database. Genetic risk scores were constructed for each individual using the unweighted sum of TXNRD2, ME3, and TXNRD2 + ME3 combined risk alleles. Age- and sex-adjusted odds ratios (ORs) for POAG diagnosis were calculated per decile for each GRS. Additionally, the clinical features of patients with POAG in the top 1%, 5%, and 10% of each GRS were compared with those in the bottom 1%, 5%, and 10%, respectively. MAIN OUTCOME MEASURES: Primary open-angle glaucoma OR per GRS decile, maximum treated intraocular pressure (IOP), and prevalence of paracentral visual field loss among patients with POAG with high versus low GRSs. RESULTS: A larger SNP effect size strongly correlated with higher TXNRD2 and lower ME3 expression levels (r = 0.95 and r = -0.97, respectively; P < 0.05 for both). Individuals in decile 10 of the TXNRD2 + ME3 GRS had the highest odds of POAG diagnosis (OR, 1.79 compared with decile 1; 95% confidence interval, 1.39-2.30; P < 0.001). Patients with POAG in the top 1% of the TXNRD2 GRS showed higher mean maximum treated IOP compared with the bottom 1% (19.9 mmHg vs. 15.6 mmHg; adjusted P = 0.03). Patients with POAG in the top 1% of the ME3 and TXNRD2 + ME3 GRS showed a higher prevalence of paracentral field loss than the bottom 1% (72.7% vs. 14.3% for ME3 GRS and 88.9% vs. 33.3% for TXNRD2+ME3 GRS; adjusted P = 0.03 for both). CONCLUSIONS: Patients with POAG with higher TXNRD2 and ME3 GRSs showed higher treated IOP and a greater prevalence of paracentral field loss. Functional studies exploring how these variants impact mitochondrial function in patients with glaucoma are warranted. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

The APOE E4 Allele Is Associated with Faster Rates of Neuroretinal Thinning in a Prospective Cohort Study of Suspect and Early Glaucoma.

Purpose: To investigate the association between the apolipoprotein E (APOE) E4 dementia-risk allele and prospective longitudinal retinal thinning in a cohort study of suspect and early manifest glaucoma. Design: Retrospective analysis of prospective cohort data. Participants: This study included all available eyes from participants recruited to the Progression Risk of Glaucoma: Relevant SNPs [single nucleotide polymorphisms] with Significant Association (PROGRESSA) study with genotyping data from which APOE genotypes could be determined. Methods: Apolipoprotein E alleles and genotypes were determined in PROGRESSA, and their distributions were compared with an age-matched and ancestrally matched normative cohort, the Blue Mountains Eye Study. Structural parameters of neuroretinal atrophy measured using spectral-domain OCT were compared within the PROGRESSA cohort on the basis of APOE E4 allele status. Main Outcome Measures: Longitudinal rates of thinning in the macular ganglion cell-inner plexiform layer (mGCIPL) complex and the peripapillary retinal nerve fiber layer (pRNFL). Results: Rates of mGCIPL complex thinning were faster in participants harboring ≥1 copies of the APOE E4 allele (ß = -0.13 µm/year; P ≤0.001). This finding was strongest in eyes affected by normal-tension glaucoma (NTG; ß = -0.20 µm/year; P = 0.003). Apolipoprotein E E4 allele carriers were also more likely to be lost to follow-up (P = 0.01) and to demonstrate a thinner average mGCIPL complex (70.9 µm vs. 71.9 µm; P = 0.011) and pRNFL (77.6 µm vs. 79.2 µm; P = 0.045) after a minimum of 3 years of monitoring. Conclusions: The APOE E4 allele was associated with faster rates of mCGIPL complex thinning, particularly in eyes with NTG. These results suggest that the APOE E4 allele may be a risk factor for retinal ganglion cell degeneration in glaucoma.

Genetic adaptations to SIV across chimpanzee populations.

Central and eastern chimpanzees are infected with Simian Immunodeficiency Virus (SIV) in the wild, typically without developing acute immunodeficiency. Yet the recent zoonotic transmission of chimpanzee SIV to humans, which were naïve to the virus, gave rise to the Human Immunodeficiency Virus (HIV), which causes AIDS and is responsible for one of the deadliest pandemics in human history. Chimpanzees have likely been infected with SIV for tens of thousands of years and have likely evolved to reduce its pathogenicity, becoming semi-natural hosts that largely tolerate the virus. In support of this view, central and eastern chimpanzees show evidence of positive selection in genes involved in SIV/HIV cell entry and immune response to SIV, respectively. We hypothesise that the population first infected by SIV would have experienced the strongest selective pressure to control the lethal potential of zoonotic SIV, and that population genetics will reveal those first critical adaptations. With that aim we used population genetics to investigate signatures of positive selection in the common ancestor of central-eastern chimpanzees. The genes with signatures of positive selection in the ancestral population are significantly enriched in SIV-related genes, especially those involved in the immune response to SIV and those encoding for host genes that physically interact with SIV/HIV (VIPs). This supports a scenario where SIV first infected the central-eastern ancestor and where this population was under strong pressure to adapt to zoonotic SIV. Interestingly, integrating these genes with candidates of positive selection in the two infected subspecies reveals novel patterns of adaptation to SIV. Specifically, we observe evidence of positive selection in numerous steps of the biological pathway responsible for T-helper cell differentiation, including CD4 and multiple genes that SIV/HIV use to infect and control host cells. This pathway is active only in CD4+ cells which SIV/HIV infects, and it plays a crucial role in shaping the immune response so it can efficiently control the virus. Our results confirm the importance of SIV as a selective factor, identify specific genetic changes that may have allowed our closest living relatives to reduce SIV's pathogenicity, and demonstrate the potential of population genomics to reveal the evolutionary mechanisms used by naïve hosts to reduce the pathogenicity of zoonotic pathogens.

Diagnostic yield of candidate genes in an Australian corneal dystrophy cohort.

Corneal dystrophies describe a clinically and genetically heterogeneous group of inherited disorders. The International Classification of Corneal Dystrophies (IC3D) lists 22 types of corneal dystrophy, 17 of which have been demonstrated to result from pathogenic variants in 19 identified genes. In this study, we investigated the diagnostic yield of genetic testing in a well-characterised cohort of 58 individuals from 44 families with different types of corneal dystrophy. Individuals diagnosed solely with Fuchs endothelial corneal dystrophy were excluded. Clinical details were obtained from the treating ophthalmologist. Participants and their family members were tested using a gene candidate and exome sequencing approach. We identified a likely molecular diagnosis in 70.5% families (31/44). The detection rate was significantly higher among probands with a family history of corneal dystrophy (15/16, 93.8%) than those without (16/28, 57.1%, p = .015), and among those who had undergone corneal graft surgery (9/9, 100.0%) compared to those who had not (22/35, 62.9%, p = .041). We identified eight novel variants in five genes and identified five families with syndromes associated with corneal dystrophies. Our findings highlight the genetic heterogeneity of corneal dystrophies and the clinical utility of genetic testing in reaching an accurate clinical diagnosis.

Population dynamics and genetic connectivity in recent chimpanzee history.

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

RNA Sequencing of Lens Capsular Epithelium Implicates Novel Pathways in Pseudoexfoliation Syndrome.

Purpose: Pseudoexfoliation syndrome (PEX) is a common systemic disease that results in severe and often irreversible vision loss. Despite considerable research effort, PEX remains incompletely understood. This study sought to perform the first RNAseq study in elucidate the pathophysiology of PEX, and contribute a publicly available transcriptomic data resource for future research. Methods: Human ocular lens capsular epithelium samples were collected from 25 patients with PEX and 39 non-PEX controls undergoing cataract surgery. RNA extracted from these specimens was subjected to polyadenylated (mRNA) selection and deep bulk RNA sequencing. Differential expression analysis investigated protein-coding gene transcripts. Exploratory analyses used pathway analysis tools, and curated class- and disease-specific gene sets. Results: Differential expression analysis demonstrated that 2882 genes were differentially expressed according to PEX status. Genes associated with viral gene expression pathways were among the most upregulated, alongside genes encoding ribosomal and mitochondrial respiratory transport chain proteins. Cell adhesion protein transcripts including type 4 collagen subunits were downregulated. Conclusions: This comparative transcriptomic dataset highlights novel and previously recognized pathogenic pathways in PEX and provides the first comprehensive transcriptomic resource, adding an additional layer to build further understanding of PEX pathophysiology.

Identification of Structural Variation in Chimpanzees Using Optical Mapping and Nanopore Sequencing.

Recent efforts to comprehensively characterize great ape genetic diversity using short-read sequencing and single-nucleotide variants have led to important discoveries related to selection within species, demographic history, and lineage-specific traits. Structural variants (SVs), including deletions and inversions, comprise a larger proportion of genetic differences between and within species, making them an important yet understudied source of trait divergence. Here, we used a combination of long-read and -range sequencing approaches to characterize the structural variant landscape of two additional Pan troglodytes verus individuals, one of whom carries 13% admixture from Pan troglodytes troglodytes. We performed optical mapping of both individuals followed by nanopore sequencing of one individual. Filtering for larger variants (>10 kbp) and combined with genotyping of SVs using short-read data from the Great Ape Genome Project, we identified 425 deletions and 59 inversions, of which 88 and 36, respectively, were novel. Compared with gene expression in humans, we found a significant enrichment of chimpanzee genes with differential expression in lymphoblastoid cell lines and induced pluripotent stem cells, both within deletions and near inversion breakpoints. We examined chromatin-conformation maps from human and chimpanzee using these same cell types and observed alterations in genomic interactions at SV breakpoints. Finally, we focused on 56 genes impacted by SVs in >90% of chimpanzees and absent in humans and gorillas, which may contribute to chimpanzee-specific features. Sequencing a greater set of individuals from diverse subspecies will be critical to establish the complete landscape of genetic variation in chimpanzees.

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.

Structural Variants and Selective Sweep Foci Contribute to Insecticide Resistance in the Drosophila Genetic Reference Panel.

Patterns of nucleotide polymorphism within populations of Drosophila melanogaster suggest that insecticides have been the selective agents driving the strongest recent bouts of positive selection. However, there is a need to explicitly link selective sweeps to the particular insecticide phenotypes that could plausibly account for the drastic selective responses that are observed in these non-target insects. Here, we screen the Drosophila Genetic Reference Panel with two common insecticides; malathion (an organophosphate) and permethrin (a pyrethroid). Genome-wide association studies map survival on malathion to two of the largest sweeps in the D. melanogaster genome; Ace and Cyp6g1 Malathion survivorship also correlates with lines which have high levels of Cyp12d1, Jheh1 and Jheh2 transcript abundance. Permethrin phenotypes map to the largest cluster of P450 genes in the Drosophila genome, however in contrast to a selective sweep driven by insecticide use, the derived allele seems to be associated with susceptibility. These results underscore previous findings that highlight the importance of structural variation to insecticide phenotypes: Cyp6g1 exhibits copy number variation and transposable element insertions, Cyp12d1 is tandemly duplicated, the Jheh loci are associated with a Bari1 transposable element insertion, and a Cyp6a17 deletion is associated with susceptibility.

Human local adaptation of the TRPM8 cold receptor along a latitudinal cline.

Ambient temperature is a critical environmental factor for all living organisms. It was likely an important selective force as modern humans recently colonized temperate and cold Eurasian environments. Nevertheless, as of yet we have limited evidence of local adaptation to ambient temperature in populations from those environments. To shed light on this question, we exploit the fact that humans are a cosmopolitan species that inhabit territories under a wide range of temperatures. Focusing on cold perception-which is central to thermoregulation and survival in cold environments-we show evidence of recent local adaptation on TRPM8. This gene encodes for a cation channel that is, to date, the only temperature receptor known to mediate an endogenous response to moderate cold. The upstream variant rs10166942 shows extreme population differentiation, with frequencies that range from 5% in Nigeria to 88% in Finland (placing this SNP in the 0.02% tail of the FST empirical distribution). When all populations are jointly analyzed, allele frequencies correlate with latitude and temperature beyond what can be explained by shared ancestry and population substructure. Using a Bayesian approach, we infer that the allele originated and evolved neutrally in Africa, while positive selection raised its frequency to different degrees in Eurasian populations, resulting in allele frequencies that follow a latitudinal cline. We infer strong positive selection, in agreement with ancient DNA showing high frequency of the allele in Europe 3,000 to 8,000 years ago. rs10166942 is important phenotypically because its ancestral allele is protective of migraine. This debilitating disorder varies in prevalence across human populations, with highest prevalence in individuals of European descent-precisely the population with the highest frequency of rs10166942 derived allele. We thus hypothesize that local adaptation on previously neutral standing variation may have contributed to the genetic differences that exist in the prevalence of migraine among human populations today.

Signatures of Long-Term Balancing Selection in Human Genomes.

Balancing selection maintains advantageous diversity in populations through various mechanisms. Although extensively explored from a theoretical perspective, an empirical understanding of its prevalence and targets lags behind our knowledge of positive selection. Here, we describe the Non-central Deviation (NCD), a simple yet powerful statistic to detect long-term balancing selection (LTBS) that quantifies how close frequencies are to expectations under LTBS, and provides the basis for a neutrality test. NCD can be applied to a single locus or genomic data, and can be implemented considering only polymorphisms (NCD1) or also considering fixed differences with respect to an outgroup (NCD2) species. Incorporating fixed differences improves power, and NCD2 has higher power to detect LTBS in humans under different frequencies of the balanced allele(s) than other available methods. Applied to genome-wide data from African and European human populations, in both cases using chimpanzee as an outgroup, NCD2 shows that, albeit not prevalent, LTBS affects a sizable portion of the genome: ∼0.6% of analyzed genomic windows and 0.8% of analyzed positions. Significant windows (P < 0.0001) contain 1.6% of SNPs in the genome, which disproportionally fall within exons and change protein sequence, but are not enriched in putatively regulatory sites. These windows overlap ∼8% of the protein-coding genes, and these have larger number of transcripts than expected by chance even after controlling for gene length. Our catalog includes known targets of LTBS but a majority of them (90%) are novel. As expected, immune-related genes are among those with the strongest signatures, although most candidates are involved in other biological functions, suggesting that LTBS potentially influences diverse human phenotypes.

Insights into DDT Resistance from the Drosophila melanogaster Genetic Reference Panel.

Insecticide resistance is considered a classic model of microevolution, where a strong selective agent is applied to a large natural population, resulting in a change in frequency of alleles that confer resistance. While many insecticide resistance variants have been characterized at the gene level, they are typically single genes of large effect identified in highly resistant pest species. In contrast, multiple variants have been implicated in DDT resistance in Drosophila melanogaster; however, only the Cyp6g1 locus has previously been shown to be relevant to field populations. Here we use genome-wide association studies (GWAS) to identify DDT-associated polygenes and use selective sweep analyses to assess their adaptive significance. We identify and verify two candidate DDT resistance loci. A largely uncharacterized gene, CG10737, has a function in muscles that ameliorates the effects of DDT, while a putative detoxifying P450, Cyp6w1, shows compelling evidence of positive selection.

Chimpanzee genomic diversity reveals ancient admixture with bonobos.

Our closest living relatives, chimpanzees and bonobos, have a complex demographic history. We analyzed the high-coverage whole genomes of 75 wild-born chimpanzees and bonobos from 10 countries in Africa. We found that chimpanzee population substructure makes genetic information a good predictor of geographic origin at country and regional scales. Multiple lines of evidence suggest that gene flow occurred from bonobos into the ancestors of central and eastern chimpanzees between 200,000 and 550,000 years ago, probably with subsequent spread into Nigeria-Cameroon chimpanzees. Together with another, possibly more recent contact (after 200,000 years ago), bonobos contributed less than 1% to the central chimpanzee genomes. Admixture thus appears to have been widespread during hominid evolution.

Genomic and Transcriptomic Associations Identify a New Insecticide Resistance Phenotype for the Selective Sweep at the Cyp6g1 Locus of Drosophila melanogaster.

Scans of the Drosophila melanogaster genome have identified organophosphate resistance loci among those with the most pronounced signature of positive selection. In this study, the molecular basis of resistance to the organophosphate insecticide azinphos-methyl was investigated using the Drosophila Genetic Reference Panel, and genome-wide association. Recently released full transcriptome data were used to extend the utility of the Drosophila Genetic Reference Panel resource beyond traditional genome-wide association studies to allow systems genetics analyses of phenotypes. We found that both genomic and transcriptomic associations independently identified Cyp6g1, a gene involved in resistance to DDT and neonicotinoid insecticides, as the top candidate for azinphos-methyl resistance. This was verified by transgenically overexpressing Cyp6g1 using natural regulatory elements from a resistant allele, resulting in a 6.5-fold increase in resistance. We also identified four novel candidate genes associated with azinphos-methyl resistance, all of which are involved in either regulation of fat storage, or nervous system development. In Cyp6g1, we find a demonstrable resistance locus, a verification that transcriptome data can be used to identify variants associated with insecticide resistance, and an overlap between peaks of a genome-wide association study, and a genome-wide selective sweep analysis.

Gene duplication in the major insecticide target site, Rdl, in Drosophila melanogaster.

The Resistance to Dieldrin gene, Rdl, encodes a GABA-gated chloride channel subunit that is targeted by cyclodiene and phenylpyrazole insecticides. The gene was first characterized in Drosophila melanogaster by genetic mapping of resistance to the cyclodiene dieldrin. The 4,000-fold resistance observed was due to a single amino acid replacement, Ala(301) to Ser. The equivalent change was subsequently identified in Rdl orthologs of a large range of resistant insect species. Here, we report identification of a duplication at the Rdl locus in D. melanogaster. The 113-kb duplication contains one WT copy of Rdl and a second copy with two point mutations: an Ala(301) to Ser resistance mutation and Met(360) to Ile replacement. Individuals with this duplication exhibit intermediate dieldrin resistance compared with single copy Ser(301) homozygotes, reduced temperature sensitivity, and altered RNA editing associated with the resistant allele. Ectopic recombination between Roo transposable elements is involved in generating this genomic rearrangement. The duplication phenotypes were confirmed by construction of a transgenic, artificial duplication integrating the 55.7-kb Rdl locus with a Ser(301) change into an Ala(301) background. Gene duplications can contribute significantly to the evolution of insecticide resistance, most commonly by increasing the amount of gene product produced. Here however, duplication of the Rdl target site creates permanent heterozygosity, providing unique potential for adaptive mutations to accrue in one copy, without abolishing the endogenous role of an essential gene.