Nucleosome reorganisation in breast cancer tissues.

BACKGROUND: Nucleosome repositioning in cancer is believed to cause many changes in genome organisation and gene expression. Understanding these changes is important to elucidate fundamental aspects of cancer. It is also important for medical diagnostics based on cell-free DNA (cfDNA), which originates from genomic DNA regions protected from digestion by nucleosomes. RESULTS: We have generated high-resolution nucleosome maps in paired tumour and normal tissues from the same breast cancer patients using MNase-assisted histone H3 ChIP-seq and compared them with the corresponding cfDNA from blood plasma. This analysis has detected single-nucleosome repositioning at key regulatory regions in a patient-specific manner and common cancer-specific patterns across patients. The nucleosomes gained in tumour versus normal tissue were particularly informative of cancer pathways, with ~ 20-fold enrichment at CpG islands, a large fraction of which marked promoters of genes encoding DNA-binding proteins. The tumour tissues were characterised by a 5-10 bp decrease in the average distance between nucleosomes (nucleosome repeat length, NRL), which is qualitatively similar to the differences between pluripotent and differentiated cells. This effect was correlated with gene activity, differential DNA methylation and changes in local occupancy of linker histone variants H1.4 and H1X. CONCLUSIONS: Our study offers a novel resource of high-resolution nucleosome maps in breast cancer patients and reports for the first time the effect of systematic decrease of NRL in paired tumour versus normal breast tissues from the same patient. Our findings provide a new mechanistic understanding of nucleosome repositioning in tumour tissues that can be valuable for patient diagnostics, stratification and monitoring.

Nucleosome repositioning in chronic lymphocytic leukemia.

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression.

Hydrophobic nanoporous carbon scaffolds reveal the origin of polarity-dependent electrocapillary imbibition.

An engineered nanoporous carbon scaffold (NCS) consisting of a 3-D interconnected 85 nm nanopore network was used here as a model material to investigate the nanoscale transport of liquids as a function of the polarity and magnitude of an applied potential ('electro-imbibition'), all in 1 M KCl solution. A camera was used to track both meniscus formation and meniscus jump, front motion dynamics, and droplet expulsion, while also quantifying the electrocapillary imbibition height (H) as a function of the applied potential of the NCS material. Although no imbibition was seen over a wide range of potentials, at positive potentials (+1.2 V vs. the potential of zero charge (pzc)), imbibition was correlated with carbon surface electro-oxidation, as confirmed by both electrochemistry and post-imbibition surface analysis, with gas evolution (O2, CO2) seen visually only after imbibition was well underway. At negative potentials, vigorous hydrogen evolution reaction was observed at the NCS/KCl solution interface, well before imbibition began at -0.5 Vpzc, proposed to be nucleated by an electrical double layer charging-driven meniscus jump, followed by processes such as Marangoni flow, adsorption induced deformation, and hydrogen pressure driven flow. This study improves the understanding of electrocapillary imbibition at the nanoscale, being highly relevant in a wide range of multidisciplinary practical applications, including in energy storage and conversion devices, energy-efficient desalination, and electrical-integrated nanofluidics design.

DNA sequence-dependent formation of heterochromatin nanodomains.

The mammalian epigenome contains thousands of heterochromatin nanodomains (HNDs) marked by di- and trimethylation of histone H3 at lysine 9 (H3K9me2/3), which have a typical size of 3-10 nucleosomes. However, what governs HND location and extension is only partly understood. Here, we address this issue by introducing the chromatin hierarchical lattice framework (ChromHL) that predicts chromatin state patterns with single-nucleotide resolution. ChromHL is applied to analyse four HND types in mouse embryonic stem cells that are defined by histone methylases SUV39H1/2 or GLP, transcription factor ADNP or chromatin remodeller ATRX. We find that HND patterns can be computed from PAX3/9, ADNP and LINE1 sequence motifs as nucleation sites and boundaries that are determined by DNA sequence (e.g. CTCF binding sites), cooperative interactions between nucleosomes as well as nucleosome-HP1 interactions. Thus, ChromHL rationalizes how patterns of H3K9me2/3 are established and changed via the activity of protein factors in processes like cell differentiation.

Fas-threshold signalling in MSCs promotes pancreatic cancer progression and metastasis.

Mesenchymal stem cells (MSCs) belong to the tumour microenvironment and have been implicated in tumour progression. We found that the number of MSCs significantly increased in tumour-burdened mice driven by Fas-threshold signalling. Consequently, MSCs lacking Fas lost their ability to induce metastasis development in a pancreatic cancer model. Mixing of MSCs with pancreatic cancer cells led to sustained production of the pro-metastatic cytokines CCL2 and IL6 by the stem cells. The levels of these cytokines were dependent on the number of MSCs, linking Fas-mediated MSC-proliferation to their capacity to promote tumour progression. Furthermore, we discovered that CCL2 and IL6 were induced by pancreatic cancer cell-derived IL1. Importantly, analysis of patient transcriptomic data revealed that high FasL expression correlates with high levels of MSC markers as well as increased IL6 and CCL2 levels in pancreatic tumours. Moreover, both FasL and CCL2 are linked to elevated levels of markers specific for monocytes known to possess further pro-metastatic activities. These results confirm our experimental findings of a FasL-MSC-IL1-CCL2/IL6 axis in pancreatic cancer and highlights the role of MSCs in tumour progression.

Experimental and computational evaluation of cyclic solvent injection in fractured tight hydrocarbon reservoirs.

Multi-fractured horizontal wells have enabled commercial production from low-permeability ('tight') hydrocarbon reservoirs but recoveries remain exceedingly small (< 5-10%). As a result, operators have investigated the use of solvent (gas) injection schemes, such as huff-n-puff (HNP), to improve oil recovery. Previous HNP laboratory approaches, classified primary as 'flow-through-matrix' and 'flow-around-matrix' typically (1) are not fully representative of field conditions at near-fracture regions and (2) require long test times, even when performed on fractured cores. The objectives of this proof-of-concept study are to (1) design and implement a new experimental procedure that better reproduces HNP schemes in near-fracture regions and (2) use the results, simulated with a compositional lab-calibrated model, to explore the controls on enhanced hydrocarbon recovery in depleted tight oil plays. Performing multiple CO2 and (simplified) lean gas HNP cycles, the integrated experimental and simulation approach proposed herein achieves the ultimate recovery factors in a significantly shorter time frame (25-50%) compared to previous studies. The integrated experimental and computational approach proposed herein is valuable for core-based evaluation of cyclic solvent (CO2, CH4) injection in tight hydrocarbon reservoirs for (1) hydrocarbon recovery and (2) subsurface greenhouse (CO2, CH4) gas disposal/storage applications.

Genetic variation at the Cyp6m2 putative insecticide resistance locus in Anopheles gambiae and Anopheles coluzzii.

BACKGROUND: The emergence of insecticide resistance is a major threat to malaria control programmes in Africa, with many different factors contributing to insecticide resistance in its vectors, Anopheles mosquitoes. CYP6M2 has previously been recognized as an important candidate in cytochrome P450-mediated detoxification in Anopheles. As it has been implicated in resistance against pyrethroids, organochlorines and carbamates, its broad metabolic activity makes it a potential agent in insecticide cross-resistance. Currently, allelic variation within the Cyp6m2 gene remains unknown. METHODS: Here, Illumina whole-genome sequence data from Phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) was used to examine genetic variation in the Cyp6m2 gene across 16 populations in 13 countries comprising Anopheles gambiae and Anopheles coluzzii mosquitoes. To identify whether these alleles show evidence of selection either through potentially modified enzymatic function or by being linked to variants that change the transcriptional profile of the gene, hierarchical clustering of haplotypes, linkage disequilibrium, median joining networks and extended haplotype homozygosity analyses were performed. RESULTS: Fifteen missense biallelic substitutions at high frequency (defined as > 5% frequency in one or more populations) are found, which fall into five distinct haplotype groups that carry the main high frequency variants: A13T, D65A, E328Q, Y347F, I359V and A468S. Despite consistent reports of Cyp6m2 upregulation and metabolic activity in insecticide resistant Anophelines, no evidence of directional selection is found occurring on these variants or on the haplotype clusters in which they are found. CONCLUSION: These results imply that emerging resistance associated with Cyp6m2 is potentially driven by distant regulatory loci such as transcriptional factors rather than by its missense variants, or that other genes are playing a more significant role in conferring metabolic resistance.

Wetting dynamics of nanoliter water droplets in nanoporous media.

HYPOTHESIS: The Lucas-Washburn (L-W) equation is the classical theory to describe the dynamics of spontaneous imbibition in single micro-channels and micro-scale porous media. However, for nanoliter droplets imbibition in nanoporous media, the L-W equation may not be suitable, due to the nanoscale liquid-solid interactions, e.g., contact line pinning and capillary condensation. In addition, for an intrinsically hydrophobic nanoporous substrate, spontaneous imbibition of a nanoliter droplet is hypothesized to occur if capillary condensation had occurred internally already. EXPERIMENTS: A nanoporous carbon scaffold was synthesized and used as a model nanoporous medium. A recently-developed micro-injection technique was used to generate a series of nanoliter water droplets (2.8-34 nL); the entire wetting dynamics (i.e., apparent contact angle and droplet volume as a function of time) were observed inside an environmental scanning electron microscope. FINDINGS: The L-W equation does not describe the wetting dynamics of nanoliter water droplets in nanoporous media. A new theoretical model is developed to characterize the corresponding dynamics. It is demonstrated that, even for an intrinsically hydrophobic nanoporous substrate, spontaneous imbibition of a nanoliter droplet can occur if capillary condensation had occurred internally already.

Linker histone epitopes are hidden by in situ higher-order chromatin structure.

BACKGROUND: Histone H1 is the most mobile histone in the cell nucleus. Defining the positions of H1 on chromatin in situ, therefore, represents a challenge. Immunoprecipitation of formaldehyde-fixed and sonicated chromatin, followed by DNA sequencing (xChIP-seq), is traditionally the method for mapping histones onto DNA elements. But since sonication fragmentation precedes ChIP, there is a consequent loss of information about chromatin higher-order structure. Here, we present a new method, xxChIP-seq, employing antibody binding to fixed intact in situ chromatin, followed by extensive washing, a second fixation, sonication and immunoprecipitation. The second fixation is intended to prevent the loss of specifically bound antibody during washing and subsequent sonication and to prevent antibody shifting to epitopes revealed by the sonication process. In many respects, xxChIP-seq is comparable to immunostaining microscopy, which also involves interaction of the primary antibody with fixed and permeabilized intact cells. The only epitopes displayed after immunostaining are the "exposed" epitopes, not "hidden" by the fixation of chromatin higher-order structure. Comparison of immunoprecipitated fragments between xChIP-seq versus xxChIP-seq should indicate which epitopes become inaccessible with fixation and identify their associated DNA elements. RESULTS: We determined the genomic distribution of histone variants H1.2 and H1.5 in human myeloid leukemia cells HL-60/S4 and compared their epitope exposure by both xChIP-seq and xxChIP-seq, as well as high-resolution microscopy, illustrating the influences of preserved chromatin higher-order structure in situ. We found that xChIP and xxChIP H1 signals are in general negatively correlated, with differences being more pronounced near active regulatory regions. Among the intriguing observations, we find that transcription-related regions and histone PTMs (i.e., enhancers, promoters, CpG islands, H3K4me1, H3K4me3, H3K9ac, H3K27ac and H3K36me3) exhibit significant deficiencies (depletions) in H1.2 and H1.5 xxChIP-seq reads, compared to xChIP-seq. These observations suggest the existence of in situ transcription-related chromatin higher-order structures stabilized by formaldehyde. CONCLUSION: Comparison of H1 xxChIP-seq to H1 xChIP-seq allows the development of hypotheses on the chromosomal localization of (stabilized) higher-order structure, indicated by the generation of "hidden" H1 epitopes following formaldehyde crosslinking. Changes in H1 epitope exposure surrounding averaged chromosomal binding sites or epigenetic modifications can also indicate whether these sites have chromatin higher-order structure. For example, comparison between averaged active or inactive promoter regions suggests that both regions can acquire stabilized higher-order structure with hidden H1 epitopes. However, the H1 xChIP-seq comparison cannot define their differences. Application of the xxChIP-seq versus H1 xChIP-seq method is particularly relevant to chromatin-associated proteins, such as linker histones, that play dynamic roles in establishing chromatin higher-order structure.

CTCF-dependent chromatin boundaries formed by asymmetric nucleosome arrays with decreased linker length.

The CCCTC-binding factor (CTCF) organises the genome in 3D through DNA loops and in 1D by setting boundaries isolating different chromatin states, but these processes are not well understood. Here we investigate chromatin boundaries in mouse embryonic stem cells, defined by the regions with decreased Nucleosome Repeat Length (NRL) for ∼20 nucleosomes near CTCF sites, affecting up to 10% of the genome. We found that the nucleosome-depleted region (NDR) near CTCF is asymmetrically located >40 nucleotides 5'-upstream from the centre of CTCF motif. The strength of CTCF binding to DNA and the presence of cohesin is correlated with the decrease of NRL near CTCF, and anti-correlated with the level of asymmetry of the nucleosome array. Individual chromatin remodellers have different contributions, with Snf2h having the strongest effect on the NRL decrease near CTCF and Chd4 playing a major role in the symmetry breaking. Upon differentiation, a subset of preserved, common CTCF sites maintains asymmetric nucleosome pattern and small NRL. The sites which lost CTCF upon differentiation are characterized by nucleosome rearrangement 3'-downstream, with unchanged NDR 5'-upstream of CTCF motifs. Boundaries of topologically associated chromatin domains frequently contain several inward-oriented CTCF motifs whose effects, described above, add up synergistically.

A working method for estimating dynamic shear velocity in the montney formation.

In this paper, we present a customized method for estimating sonic shear velocity (Vs) from compressional velocity (Vp) logs in the Montney Formation, in wells lacking dipole sonic data. Following a multi-scenario analysis that comprised of assessing empirical Vs estimation relations [including lithology, porosity (Ø), and volume of clay (Vsh)-based Vs estimation techniques], bivariate statistics, and machine learning, we found that the Greenberg & Castagna (1992) shale lithology constants yield Vs log estimates that best match the measured Montney Formation Vs in our study area, with a regional correlation coefficient of 0.8. We have therefore customized the Vs estimation method in our study to use the Greenberg & Castagna (1992) shale lithology constants. Our working method: •Improves the efficacy of Vs log estimation from Vp logs in the study area•Demonstrates the importance of calibrating empirical relations for Vs estimation to a specific formation, and•Provides a more accurate complementary Vs log dataset for subsequent regional reservoir characterization studies.

DNA (de)methylation in embryonic stem cells controls CTCF-dependent chromatin boundaries.

Coordinated changes of DNA (de)methylation, nucleosome positioning, and chromatin binding of the architectural protein CTCF play an important role for establishing cell-type-specific chromatin states during differentiation. To elucidate molecular mechanisms that link these processes, we studied the perturbed DNA modification landscape in mouse embryonic stem cells (ESCs) carrying a double knockout (DKO) of the Tet1 and Tet2 dioxygenases. These enzymes are responsible for the conversion of 5-methylcytosine (5mC) into its hydroxymethylated (5hmC), formylated (5fC), or carboxylated (5caC) forms. We determined changes in nucleosome positioning, CTCF binding, DNA methylation, and gene expression in DKO ESCs and developed biophysical models to predict differential CTCF binding. Methylation-sensitive nucleosome repositioning accounted for a significant portion of CTCF binding loss in DKO ESCs, whereas unmethylated and nucleosome-depleted CpG islands were enriched for CTCF sites that remained occupied. A number of CTCF sites also displayed direct correlations with the CpG modification state: CTCF was preferentially lost from sites that were marked with 5hmC in wild-type (WT) cells but not from 5fC-enriched sites. In addition, we found that some CTCF sites can act as bifurcation points defining the differential methylation landscape. CTCF loss from such sites, for example, at promoters, boundaries of chromatin loops, and topologically associated domains (TADs), was correlated with DNA methylation/demethylation spreading and can be linked to down-regulation of neighboring genes. Our results reveal a hierarchical interplay between cytosine modifications, nucleosome positions, and DNA sequence that determines differential CTCF binding and regulates gene expression.

Slow polymer diffusion on brush-patterned surfaces in aqueous solution.

A model system for the investigation of diffusional transport in compartmentalized nanosystems is described. Arrays of "corrals" enclosed within poly[oligo(ethylene glycol)methyl ether methacrylate] (POEGMA) "walls" were fabricated using double-exposure interferometric lithography to deprotect aminosilane films protected by a nitrophenyl group. In exposed regions, removal of the nitrophenyl group enabled attachment of an initiator for the atom-transfer radical polymerization of end-grafted POEGMA (brushes). Diffusion coefficients for poly(ethylene glycol) in these corrals were obtained by fluorescence correlation spectroscopy. Two modes of surface diffusion were observed: one which is similar to diffusion on the unpatterned surface and a very slow mode of surface diffusion that becomes increasingly important as confinement increases. Diffusion within the POEGMA brushes does not significantly contribute to the results.

Poly(ADP-ribosyl)ation associated changes in CTCF-chromatin binding and gene expression in breast cells.

CTCF is an evolutionarily conserved and ubiquitously expressed architectural protein regulating a plethora of cellular functions via different molecular mechanisms. CTCF can undergo a number of post-translational modifications which change its properties and functions. One such modifications linked to cancer is poly(ADP-ribosyl)ation (PARylation). The highly PARylated CTCF form has an apparent molecular mass of 180 kDa (referred to as CTCF180), which can be distinguished from hypo- and non-PARylated CTCF with the apparent molecular mass of 130 kDa (referred to as CTCF130). The existing data accumulated so far have been mainly related to CTCF130. However, the properties of CTCF180 are not well understood despite its abundance in a number of primary tissues. In this study we performed ChIP-seq and RNA-seq analyses in human breast cells 226LDM, which display predominantly CTCF130 when proliferating, but CTCF180 upon cell cycle arrest. We observed that in the arrested cells the majority of sites lost CTCF, whereas fewer sites gained CTCF or remain bound (i.e. common sites). The classical CTCF binding motif was found in the lost and common, but not in the gained sites. The changes in CTCF occupancies in the lost and common sites were associated with increased chromatin densities and altered expression from the neighboring genes. Based on these results we propose a model integrating the CTCF130/180 transition with CTCF-DNA binding and gene expression changes. This study also issues an important cautionary note concerning the design and interpretation of any experiments using cells and tissues where CTCF180 may be present.

Live Imaging of Micro-Wettability Experiments Performed for Low-Permeability Oil Reservoirs.

Low-permeability (unconventional) hydrocarbon reservoirs exhibit a complex nanopore structure and micro (µm) -scale variability in composition which control fluid distribution, displacement and transport processes. Conventional methods for characterizing fluid-rock interaction are however typically performed at a macro (mm) -scale on rock sample surfaces. In this work, innovative methods for the quantification of micro-scale variations in wettability and fluid distribution in a low-permeability oil reservoir was enabled by using an environmental scanning electron microscope. Live imaging of controlled water condensation/evaporation experiments allowed micro-droplet contact angles to be evaluated, while imaging combined with x-ray mapping of cryogenically frozen samples facilitated the evaluation of oil and water micro-droplet contact angles after successive fluid injection. For the first time, live imaging of fluids injected through a micro-injection system has enabled quantification of sessile and dynamic micro-droplet contact angles. Application of these combined methods has revealed dramatic spatial changes in fluid contact angles at the micro-scale, calling into question the applicability of macro-scale observations of fluid-rock interaction.

Massive introgression drives species radiation at the range limit of Anopheles gambiae.

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

Cognitive analysis of schizophrenia risk genes that function as epigenetic regulators of gene expression.

Epigenetic mechanisms are an important heritable and dynamic means of regulating various genomic functions, including gene expression, to orchestrate brain development, adult neurogenesis, and synaptic plasticity. These processes when perturbed are thought to contribute to schizophrenia pathophysiology. A core feature of schizophrenia is cognitive dysfunction. For genetic disorders where cognitive impairment is more severe such as intellectual disability, there are a disproportionally high number of genes involved in the epigenetic regulation of gene transcription. Evidence now supports some shared genetic aetiology between schizophrenia and intellectual disability. GWAS have identified 108 chromosomal regions associated with schizophrenia risk that span 350 genes. This study identified genes mapping to those loci that have epigenetic functions, and tested the risk alleles defining those loci for association with cognitive deficits. We developed a list of 350 genes with epigenetic functions and cross-referenced this with the GWAS loci. This identified eight candidate genes: BCL11B, CHD7, EP300, EPC2, GATAD2A, KDM3B, RERE, SATB2. Using a dataset of Irish psychosis cases and controls (n = 1235), the schizophrenia risk SNPs at these loci were tested for effects on IQ, working memory, episodic memory, and attention. Strongest associations were for rs6984242 with both measures of IQ (P = 0.001) and episodic memory (P = 0.007). We link rs6984242 to CHD7 via a long range eQTL. These associations were not replicated in independent samples. Our study highlights that a number of genes mapping to risk loci for schizophrenia may function as epigenetic regulators of gene expression but further studies are required to establish a role for these genes in cognition. © 2016 Wiley Periodicals, Inc.

Characterization of Diblock Copolymer Order-Order Transitions in Semidilute Aqueous Solution Using Fluorescence Correlation Spectroscopy.

The temperature and pH-dependent diffusion of poly(glycerol monomethacrylate)-block-poly(2-hydroxypropyl methacrylate) nanoparticles prepared via polymerization-induced self-assembly in water is characterized using fluorescence correlation spectroscopy (FCS). Lowering the solution temperature or raising the solution pH induces a worm-to-sphere transition and hence an increase in diffusion coefficient by a factor of between four and eight. FCS enables morphological transitions to be monitored at relatively high copolymer concentrations (10% w/w) compared to those required for dynamic light scattering (0.1% w/w). This is important because such transitions are reversible at the former concentration, whereas they are irreversible at the latter. Furthermore, the FCS data suggest that the thermal transition takes place over a very narrow temperature range (less than 2 °C). These results demonstrate the application of FCS to characterize order-order transitions, as opposed to order-disorder transitions.

Contemporary evolution of resistance at the major insecticide target site gene Ace-1 by mutation and copy number variation in the malaria mosquito Anopheles gambiae.

Functionally constrained genes are ideal insecticide targets because disruption is often fatal, and resistance mutations are typically costly. Synaptic acetylcholinesterase (AChE) is an essential neurotransmission enzyme targeted by insecticides used increasingly in malaria control. In Anopheles and Culex mosquitoes, a glycine-serine substitution at codon 119 of the Ace-1 gene confers both resistance and fitness costs, especially for 119S/S homozygotes. G119S in Anopheles gambiae from Accra (Ghana) is strongly associated with resistance, and, despite expectations of cost, resistant 119S alleles are increasing significantly in frequency. Sequencing of Accra females detected only a single Ace-1 119S haplotype, whereas 119G diversity was high overall but very low at non-synonymous sites, evidence of strong purifying selection driven by functional constraint. Flanking microsatellites showed reduced diversity, elevated linkage disequilibrium and high differentiation of 119S, relative to 119G homozygotes across up to two megabases of the genome. Yet these signals of selection were inconsistent and sometimes weak tens of kilobases from Ace-1. This unexpected finding is attributable to apparently ubiquitous amplification of 119S alleles as part of a large copy number variant (CNV) far exceeding the size of the Ace-1 gene, whereas 119G alleles were unduplicated. Ace-1 CNV was detectable in archived samples collected when the 119S allele was rare in Ghana. Multicopy amplification of resistant alleles has not been observed previously and is likely to underpin the recent increase in 119S frequency. The large CNV compromised localization of the strong selective sweep around Ace-1, emphasizing the need to integrate CNV analysis into genome scans for selection.

Evolving the world's most dangerous animal.

The devastating consequences of malaria are well known but many mysteries remain about its key protagonists, a handful of Anopheles species. New work provides a framework for solving such puzzles, by generation and analysis of whole genome assemblies for 16 Anopheles species, with genomic flexibility a key emergent theme.