Search

Pyrethroid-impregnated bed nets have driven considerable reductions in malaria-associated morbidity and mortality in Africa since the beginning of the century1. The intense selection pressure exerted by bed nets has precipitated widespread and escalating resistance to pyrethroids in African Anopheles populations, threatening to reverse the gains that been made by malaria control2. Here we show that expression of a sensory appendage protein (SAP2), which is enriched in the legs, confers pyrethroid resistance to Anopheles gambiae. Expression of SAP2 is increased in insecticide-resistant populations and is further induced after the mosquito comes into contact with pyrethroids. SAP2 silencing fully restores mortality of the mosquitoes, whereas SAP2 overexpression results in increased resistance, probably owing to high-affinity binding of SAP2 to pyrethroid insecticides. Mining of genome sequence data reveals a selective sweep near the SAP2 locus in the mosquito populations of three West African countries (Cameroon, Guinea and Burkina Faso) with the observed increase in haplotype-associated single-nucleotide polymorphisms mirroring the increasing resistance of mosquitoes to pyrethroids reported in Burkina Faso. Our study identifies a previously undescribed mechanism of insecticide resistance that is likely to be highly relevant to malaria control efforts. The leg-enriched sensory appendage protein, SAP2, confers pyrethroid resistance to Anopheles gambiae, through high-affinity binding of pyrethroid insecticides; an observed selective sweep in field mosquitoes mirrors the increasing resistance reported in Africa.

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae, caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so.

Summary Background Clinical prognostic groupings for localised prostate cancers are imprecise, with 30–50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. Methods We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. Findings Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1–9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65–0·76]) and radical prostatectomy (4·0 [1·6–9·7]; p=0·0024; AUC 0·57 [0·52–0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2–12]; p=0·019; AUC 0·67 [0·61–0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0–19]; p=0·0015; AUC 0·74 [95% CI 0·65–0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4–6·0]; p=0·0039; AUC 0·68 [95% CI 0·63–0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. Interpretation This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. Funding Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

Background 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent congener of the dioxin class of environmental contaminants. Exposure to TCDD causes a wide range of toxic outcomes, ranging from chloracne to acute lethality. The severity of toxicity is highly dependent on the aryl hydrocarbon receptor (AHR). Binding of TCDD to the AHR leads to changes in transcription of numerous genes. Studies evaluating the transcriptional changes brought on by TCDD may provide valuable insight into the role of the AHR in human health and disease. We therefore compiled a collection of transcriptomic datasets that can be used to aid the scientific community in better understanding the transcriptional effects of ligand-activated AHR. Results Specifically, we have created a datasets package – TCDD.Transcriptomics – for the R statistical environment, consisting of 63 unique experiments comprising 377 samples, including various combinations of 3 species (human derived cell lines, mouse and rat), 4 tissue types (liver, kidney, white adipose tissue and hypothalamus) and a wide range of TCDD exposure times and doses. These datasets have been fully standardized using consistent preprocessing and annotation packages (available as of September 14, 2015). To demonstrate the utility of this R package, a subset of “AHR-core” genes were evaluated across the included datasets. Ahrr, Nqo1 and members of the Cyp family were significantly induced following exposure to TCDD across the studies as expected while Aldh3a1 was induced specifically in rat liver. Inmt was altered only in liver tissue and primarily by rat-AHR. Conclusions Analysis of the “AHR-core” genes demonstrates a continued need for studies surrounding the impact of AHR-activity on the transcriptome; genes believed to be consistently regulated by ligand-activated AHR show surprisingly little overlap across species and tissues. Until now, a comprehensive assessment of the transcriptome across these studies was challenging due to differences in array platforms, processing methods and annotation versions. We believe that this package, which is freely available for download (http://labs.oicr.on.ca/boutros-lab/tcdd-transcriptomics) will prove to be a highly beneficial resource to the scientific community evaluating the effects of TCDD exposure as well as the variety of functions of the AHR. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3446-z) contains supplementary material, which is available to authorized users.

Prostate tumours are highly variable in their response to therapies, but clinically available prognostic factors can explain only a fraction of this heterogeneity. Here we analysed 200 whole-genome sequences and 277 additional whole-exome sequences from localized, non-indolent prostate tumours with similar clinical risk profiles, and carried out RNA and methylation analyses in a subset. These tumours had a paucity of clinically actionable single nucleotide variants, unlike those seen in metastatic disease. Rather, a significant proportion of tumours harboured recurrent non-coding aberrations, large-scale genomic rearrangements, and alterations in which an inversion repressed transcription within its boundaries. Local hypermutation events were frequent, and correlated with specific genomic profiles. Numerous molecular aberrations were prognostic for disease recurrence, including several DNA methylation events, and a signature comprised of these aberrations outperformed well-described prognostic biomarkers. We suggest that intensified treatment of genomically aggressive localized prostate cancer may improve cure rates.

© 2014 Nature America, Inc. All rights reserved., Polymorphous low-grade adenocarcinoma (PLGA) is the second most frequent type of malignant tumor of the minor salivary glands. We identified PRKD1 hotspot mutations encoding p.Glu710Asp in 72.9% of PLGAs but not in other salivary gland tumors. Functional studies demonstrated that this kinase-activating alteration likely constitutes a driver of PLGA., This work was supported in part by an IDEAS grant from Princess Margaret Hospital, the Head and Neck Translational Research Program (I.W., B.A.C., P.C.B. and J.D.M.), the Ontario Institute for Cancer Research and the government of Ontario (P.C.B. and J.D.M.) and by a Terry Fox Research Institute New Investigator Award (P.C.B.). C.H. and F.-F.L. acknowledge support from the Wharton family, Joe's Team, Gordon Tozer, the Campbell Family Institute for Cancer Research and the Ministry of Health and Long-Term Planning, Canada.

The sustainability of malaria control in Africa is threatened by rising levels of insecticide resistance, and new tools to prevent malaria transmission are urgently needed. To gain a better understanding of the mosquito populations that transmit malaria, we sequenced the genomes of 765 wild specimens ofAnopheles gambiaeandAnopheles coluzziisampled from 15 locations across Africa. The data reveal high levels of genetic diversity, with over 50 million single nucleotide polymorphisms across the 230 Mbp genome. We observe complex patterns of population structure and marked variations in local population size, some of which may be due at least in part to malaria control interventions. Insecticide resistance genes show strong signatures of recent selection associated with multiple independent mutations spreading over large geographical distances and between species. The genetic variability of natural populations substantially reduces the target space for novel gene-drive strategies for mosquito control. This large dataset provides a foundation for tracking the emergence and spread of insecticide resistance and developing new vector control tools.

Risk factors for central nervous system (CNS) tumours in children remain largely unknown. Evidence of an inverse relationship between atopy and tumour development exists in adults but little is known about childhood tumours. This study aims to examine the risk of childhood CNS tumours given a history of eczema and asthma. Cases of children diagnosed with CNS tumours (n=575) and controls (n=6292) from the UK Childhood Cancer Study (UKCCS) were analysed using conditional logistic regression comparing reported histories of allergic disease. Asthma was statistically significantly and negatively associated with all CNS tumours (odds ratios, OR 0.75, confidence of interval, CI(95%): 0.58-0.97), though this was not observed for eczema (OR 0.94, CI(95%): 0.74-1.18). Individuals who had suffered both asthma and eczema showed the most significant reduction in risk (OR 0.48, CI(95%): 0.28-0.81). Analysis by tumour subtype showed the strongest effect for the medulloblastoma/PNET group. These results may have a biological explanation with raised immunosurveillance in atopic individuals protecting against the development of brain tumours. Alternative explanations might include bias, reverse causality or confounding.

The emergence of next generation DNA sequencing technology is enabling high-resolution cancer genome analysis. Large-scale projects like the International Cancer Genome Consortium (ICGC) are systematically scanning cancer genomes to identify recurrent somatic mutations. Second generation DNA sequencing, however, is still an evolving technology and procedures, both experimental and analytical, are constantly changing. Thus the research community is still defining a set of best practices for cancer genome data analysis, with no single protocol emerging to fulfil this role. Here we describe an extensive benchmark exercise to identify and resolve issues of somatic mutation calling. Whole genome sequence datasets comprising tumor-normal pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, were shared within the ICGC and submissions of somatic mutation calls were compared to verified mutations and to each other. Varying strategies to call mutations, incomplete awareness of sources of artefacts, and even lack of agreement on what constitutes an artefact or real mutation manifested in widely varying mutation call rates and somewhat low concordance among submissions. We conclude that somatic mutation calling remains an unsolved problem. However, we have identified many issues that are easy to remedy that are presented here. Our study highlights critical issues that need to be addressed before this valuable technology can be routinely used to inform clinical decision-making.

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

We describe the presentation, management, and clinical outcome of a massive acinic cell carcinoma of the parotid gland. The primary tumor and blood underwent exome sequencing which revealed deletions in CDKN2A as well as PPP1R13B, which induces p53. A damaging nonsynonymous mutation was noted in EP300, a histone acetylase which plays a role in cellular proliferation. This study provides the first insights into the genetic underpinnings of this cancer. Future large-scale efforts will be necessary to define the mutational landscape of salivary gland malignancies to identify therapeutic targets and biomarkers of treatment failure.

Little is known regarding the aetiology of central nervous system tumors in children. Recent studies have speculated on a potential infectious aetiology, but no clear associations have been found. This article uses parent reported questionnaire data from the UK Childhood Cancer Study (UKCCS), a population-based case-control study, to examine the relationship between the infectious exposure in the first year of life and the likelihood of developing a CNS tumor. The variables representing infectious exposure were social contact (including social contact with other infants and attendance at informal and formal day care), sharing a bedroom with another child, birth order, and exposure to a school-age child within the home. Children reported to have had no social contact with other infants in the first year of life displayed an increased risk of developing a CNS tumor when compared to those who had (OR 1.37, 95% CI 1.08-1.75). This effect was most prominent in the primitive neuroectodermal tumor/medulloblastoma subgroup (OR 1.78, 95% CI 1.12-2.83). Those who had attended informal (OR 0.86, 95% CI 0.68-1.09) or formal day care (OR 0.93, 95% CI 0.68-1.26) showed slightly non-statistically significant reduced risks when compared to those reporting social contact only. No association with any of the other variables was observed. Overall, the inconsistent findings by variable and tumor subtype suggest that an early exposure to infections is not strongly implicated in the aetiology of CNS tumors. However, the effect for social contact outside the home, particularly for PNET/medulloblastomas warrants further investigation.

We investigated infant feeding habits in relation to risk of childhood central nervous system tumours among 633 cases in the UK Childhood Cancer Study (UKCCS). No significant effect of breastfeeding was detected overall (odds ratio 1.01, confidence interval: 0.85–1.21) nor in any morphological subgroup. Similarly, no effect for the duration of breastfeeding or any other feeding practices was observed.

Prostate cancer (CaP) remains the most common male malignancy worldwide, leading to over 300,000 deaths per year. In Western countries, most prostate tumours are diagnosed while they are confined to the prostate and have relatively indolent histology, as assessed by the Gleason Score (GS). CaP is a C-class tumour, characterized by large number of driver copy-number aberrations and genomic rearrangements. Therefore, while previous sequencing studies have focused largely on the coding regions of late-stage disease, herein we comprehensively characterized the copy-number profiles of 250 localized prostate cancers and analyzed the whole genomes of 124 matched tumour/normal pairs derived from patients with GS6 and GS7 prostate cancer. Using this – the largest whole-genome sequencing dataset of prostate cancer to date – we confirm the C-class character of the disease and identify strong genomic subtypes that stretch across multiple types of somatic alteration, including SNVs, CNAs and genomic rearrangements. We provide the first assessments of localized hyper-mutation phenomena (chromothripsis and kataegis) in prostate cancer, and identify specific genes driving higher levels of these hyper-mutations. We identify unexpected biases in the location and role of both non-coding SNVs and genomic rearrangements, including clear association with epigenetic processes, and with genome-wide profiling of methylation in 92 samples. Finally, we demonstrate a stark paucity of clinically-actionable mutations in localized GS6 and GS7 disease, even lacking those common in high-risk localized disease, indicating that novel therapeutic development against the recurrent targets identified here will be key to allowing less-aggressive, targeted treatment of early-stage disease. Citation Format: Michael E. Fraser, Veronica Y. Sabelnykova, Takafumi N. Yamaguchi, Alice Meng, Lawrence E. Heisler, Junyan Zhang, Julie Livingstone, Vincent Huang, Andre P. Masella, Fouad Yousif, Michael Xie, Nicholas J. Harding, Xihui Lin, Haiying Kong, Stephenie D. Prokopec, Alejandro Berlin, Dominique Trudel, Xuemei Luo, Timothy E. Beck, Richard de Borja, Alister D'Costa, Robert E. Denroche, Natalie S. Fox, Emilie Lalonde, Ada Wong, Taryne Chong, Michelle Sam, Jeremy Johns, Lee Timms, Nicholas Buchner, Michele Orain, Valerie Picard, Helene Hovington, Kenneth C. Chu, Christine P'ng, Bryan Lo, Francis Nguyen, Kathleen E. Houlahan, Christopher Cooper, Shaylan K. Govind, Clement Fung, Louis Lacombe, Colin C. Collins, Yves Fradet, Bernard Tetu, Theodorus van der Kwast, John McPherson, Thomas J. Hudson, Rob G. Bristow, Paul Boutros. The mutational landscape of localized gleason 6 and 7 prostate cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2966. doi:10.1158/1538-7445.AM2015-2966

Men with localized prostate cancer vary widely in clinical outcome, with a 30-50% failure rate after primary treatment. There is thus significant interest in developing genomically refined prognostic groups. We sought to evaluate the extent of genetic heterogeneity, both between patients (inter-prostate) and within individual prostate glands (intra-prostate) to assess the impact of genetic heterogeneity on risk stratification within a tight clinical cohort. Copy number aberrations (CNAs) from 75 Gleason 7 patients were determined by OncoScan SNP microarrays. We measure the percentage of genome involved in a CNA, termed percent genome aberration (PGA), a proxy for genomic instability. Additionally, whole genome sequencing was applied to 10 intermediate-risk prostate tumours and matched blood, including multiple manually macro-dissected regions from 5 of the prostates (range 2 to 9). Somatic single nucleotide variants (SNVs) and genomic rearrangements (GR) were extracted from each patient. We find a high degree of inter-prostatic heterogeneity between the 75 Gleason 7 patients, with the number of CNAs per patient ranging from 0 to 929, corresponding to PGA 0 to 16.7%. Known prognostic markers can differentiate between patients at higher risk for biochemical recurrence, but only account for a fraction of the cohort. Notably, when these prognostic genes are examined within multiple regions of five independent tumours, they differ in copy number between cancerous regions of the same prostate. For example, TP53 is deleted in 1/2, 1/3, 4/9, 0/4, and 4/5 prostate regions. Indeed, phylogenetic analysis of geographically distinct regions revealed multi-clonal disease in two of the five patients; separate analyses based on SNVs, CNAs, and GRs all concluded that these patient have two genetically distinct cancers within their prostate. We demonstrate dramatic levels of inter- and intra- prostate genetic heterogeneity within pathologically identical or similar cancers. The observed intra-prostatic genomic heterogeneity, both in terms of multi-focal and multi-clonal disease, has critical implications for clinical management. Prognostic information obtained by biopsy may be inconsistent depending on the site of biopsy, and applying personalized medicine to prostate cancer will be challenging. This study highlights the need for further evaluation of how intra-prostatic heterogeneity is related to patient prognosis. Citation Information: Mol Cancer Ther 2013;12(11 Suppl):B129. Citation Format: Emilie Lalonde, Paul C. Boutros, Michael Fraser, Richard de Borja, Nicholas J. Harding, Dominique Trudel, Alice Meng, Pablo H. Hennings-Yeomans, Andrew McPherson, Amin Zia, Jianxin Wang, Timothy Beck, Natalie S. Fox, Taryne Chong, Michelle Sam, Jeremy Johns, Lee Timms, Nicholas Buchner, Sohrab Shah, Cenk Sahinalp, Thomas J. Hudson, John D. McPherson, Theodorus van der Kwast, Robert G. Bristow. Clinical implications of inter- and intra- prostatic heterogeneity. [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2013 Oct 19-23; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2013;12(11 Suppl):Abstract nr B129.

Intermediate risk prostate cancer (CaP) with Gleason score (GS) of 7 show up to 100x variability in genetic instability. As CaP is multifocal and likely multiclonal, there is a need to characterize heterogeneity for patient stratification, which would increase the ability to act on genomic information by adding adjuvant therapies to offset systemic occult metastases that currently limit cure in ∼30% of patients. Individual genetic portraits could be used to improve cure on combined clinical-molecular staging criteria. We undertook a pilot study to assess the genetic heterogeneity of potentially curable GS=7 CaP. We selected 10 men with GS=7 CaP; 5 treated with external beam radiotherapy (frozen pre-treatment biopsies) and 5 treated with radical prostatectomy (RadP, frozen tumour). Additionally, DNA from 18 distinct formalin-fixed, paraffin-embedded (FFPE) foci from the 5 RadP were analysed. Each of these 28 foci were subjected to whole-genome sequencing (WGS) and OncoScan SNP arrays to yield comprehensive genetic profiles. mRNA expression was evaluated on frozen RadP by microarray. Germline DNA from whole-blood was also analysed. Following independent pathology reviews and manual macro-dissection of tumour areas of ≥70% cellularity, WGS (≥50x tumour, ≥30x germline) was performed on as little as 50 ng genomic DNA, and OncoScan arrays were performed using as little as 30ng DNA using either amplified or innate genomic DNA. Regions of CaP in FFPE RadP were recorded using a tissue map to identify independent malignant foci, and ERG immunostaining was performed to assist in the identification. In cases where ERG-positive and -negative foci were adjacent, ERG staining was repeated on an un-stained slide to confirm separate foci based on 3D multi-section analyses. ERG fusion status was also assessed in frozen samples by aCGH or IHC. Validation of SNVs via SNP array and deep-resequencing showed ∼99% accuracy. Tumour cellularity was estimated using Qpure and was >60% for all samples. Phylogenetic techniques were used to demonstrate clear multi-clonality in two tumours. Across all tumours, ∼50% of SNVs were specific to an individual tumour-region. Phylogenies were confirmed with both SNVs and CNAs, but CNAs generally exhibited greater concordance amongst different regions of the same tumour. Some previously observed recurrent mutations were previously identified as recurrent in CaP (e.g. SPOP), and the overall mutation rate for intermediate-risk CaP was only somewhat below that reported for castrate-resistant disease (11,230 somatic SNVs per tumour). Our studies support the concept that a complete characterization of inter- and intra-CaP heterogeneity is possible in fresh and archival tissues; the latter is important for correlations to clinical outcome. These approaches can then be streamlined for high-throughput analyses within personalized medicine laboratories leading to “point of care” molecular tests and individualization of therapy. Citation Format: Michael E. Fraser, Richard de Borja, Dominique Trudel, Nicholas J. Harding, Pablo H. Hennings-Yeomans, Alice Meng, Emilie R. Lalonde, Andrew Brown, Natalie S. Fox, Taryne Chong, Amin Zia, Michelle Sam, Jianxin Wang, Michelle A. Chan-Seng-Yue, Jeremy Johns, Lee Timms, Nicholas Buchner, Ada Wong, Fouad Yousif, Rob Denroche, Gaetano Zafarana, Maud HW Starmans, Hanbert Chen, Shaylan Govind, Francis Nguyen, Melania Pintilie, Neil Fleshner, Stanislav Volik, Lakshmi Muthuswamy, Colin C. Collins, Thomas J. Hudson, Lincoln D. Stein, Timothy Beck, John D. McPherson, Theodorus van der Kwast, Paul C. Boutros, Rob G. Bristow. A molecular portrait of potentially curable prostate cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2003. doi:10.1158/1538-7445.AM2013-2003

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii. We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures.