Development of a CRISPR/Cas9-based gene drive system in Anopheles gambiae, the main vector for the malaria parasite, paves the way for control of this pest insect. Gene drive systems that enable super-Mendelian inheritance of a transgene have the potential to modify insect populations over a timeframe of a few years. We describe CRISPR-Cas9 endonuclease constructs that function as gene drive systems in Anopheles gambiae, the main vector for malaria. We identified three genes (AGAP005958, AGAP011377 and AGAP007280) that confer a recessive female-sterility phenotype upon disruption, and inserted into each locus CRISPR-Cas9 gene drive constructs designed to target and edit each gene. For each targeted locus we observed a strong gene drive at the molecular level, with transmission rates to progeny of 91.4 to 99.6%. Population modeling and cage experiments indicate that a CRISPR-Cas9 construct targeting one of these loci, AGAP007280, meets the minimum requirement for a gene drive targeting female reproduction in an insect population. These findings could expedite the development of gene drives to suppress mosquito populations to levels that do not support malaria transmission.

A continuous-time dynamic model of a network of $N$ nonlinear elements interacting via random asymmetric couplings is studied. A self-consistent mean-field theory, exact in the $N\to \infty$ limit, predicts a transition from a stationary phase to a chaotic phase occurring at a critical value of the gain parameter. The autocorrelations of the chaotic flow as well as the maximal Lyapunov exponent are calculated.

In the human malaria vector Anopheles gambiae, the gene doublesex (Agdsx) encodes two alternatively spliced transcripts, dsx-female (AgdsxF) and dsx-male (AgdsxM), that control differentiation of the two sexes. The female transcript, unlike the male, contains an exon (exon 5) whose sequence is highly conserved in all Anopheles mosquitoes so far analyzed. We found that CRISPR–Cas9-targeted disruption of the intron 4–exon 5 boundary aimed at blocking the formation of functional AgdsxF did not affect male development or fertility, whereas females homozygous for the disrupted allele showed an intersex phenotype and complete sterility. A CRISPR–Cas9 gene drive construct targeting this same sequence spread rapidly in caged mosquitoes, reaching 100% prevalence within 7–11 generations while progressively reducing egg production to the point of total population collapse. Owing to functional constraint of the target sequence, no selection of alleles resistant to the gene drive occurred in these laboratory experiments. Cas9-resistant variants arose in each generation at the target site but did not block the spread of the drive.

INTRODUCTION: Control of mosquito vectors has historically proven to be an effective means of eliminating malaria. Human malaria is transmitted only by mosquitoes in the genus Anopheles , but not all species within the genus, or even all members of each vector species, are efficient malaria vectors. Variation in vectorial capacity for human malaria among Anopheles mosquito species is determined by many factors, including behavior, immunity, and life history. RATIONALE: This variation in vectorial capacity suggests an underlying genetic/genomic plasticity that results in variation of key traits determining vectorial capacity within the genus. Sequencing the genome of Anopheles gambiae , the most important malaria vector in sub-Saharan Africa, has offered numerous insights into how that species became highly specialized to live among and feed upon humans and how susceptibility to mosquito control strategies is determined. Until very recently, similar genomic resources have not existed for other anophelines, limiting comparisons to individual genes or sets of genomic markers with no genome-wide data to investigate attributes associated with vectorial capacity across the genus. RESULTS: We sequenced and assembled the genomes and transcriptomes of 16 anophelines from Africa, Asia, Europe, and Latin America, spanning ~100 million years of evolution and chosen to represent a range of evolutionary distances from An. gambiae , a variety of geographic locations and ecological conditions, and varying degrees of vectorial capacity. Genome assembly quality reflected DNA template quality and homozygosity. Despite variation in contiguity, the assemblies were remarkably complete and searches for arthropod-wide single-copy orthologs generally revealed few missing genes. Genome annotation supported with RNA sequencing transcriptomes yielded between 10,738 and 16,149 protein-coding genes for each species. Relative to Drosophila, the closest dipteran genus for which equivalent genomic resources exist, Anopheles exhibits a dynamic genomic evolutionary profile. Comparative analyses show a fivefold faster rate of gene gain and loss, elevated gene shuffling on the X chromosome, and more intron losses in Anopheles . Some determinants of vectorial capacity, such as chemosensory genes, do not show elevated turnover but instead diversify through protein-sequence changes. We also document evidence of variation in important reproductive phenotypes, genes controlling immunity to Plasmodium malaria parasites and other microbes, genes encoding cuticular and salivary proteins, and genes conferring metabolic insecticide resistance. This dynamism of anopheline genes and genomes may contribute to their flexible capacity to take advantage of new ecological niches, including adapting to humans as primary hosts. CONCLUSIONS: Anopheline mosquitoes exhibit a molecular evolutionary profile very distinct from Drosophila , and their genomes harbor strong evidence of functional variation in traits that determine vectorial capacity. These 16 new reference genome assemblies provide a foundation for hypothesis generation and testing to further our understanding of the diverse biological traits that determine vectorial capacity. Geography, vector status, and molecular phylogeny of the 16 newly sequenced anopheline mosquitoes and selected other dipterans.: The maximum likelihood molecular phylogeny of all sequenced anophelines and two mosquito outgroups was constructed from the aligned protein sequences of 1085 single-copy orthologs. Shapes between branch termini and species names indicate vector status and are colored according to geographic ranges depicted on the map.

The average eigenvalue distribution $\rho \left(\lambda \right)$ of $N\times N$ real random asymmetric matrices $J_{\mathrm{ij}} \left(J_{\mathrm{ji}}\ne J_{\mathrm{ij}}\right)$ is calculated in the limit of $N\to \infty$. It is found that $\rho \left(\lambda \right)$ is uniform in an ellipse, in the complex plane, whose real and imaginary axes are $1+\tau$ and $1-\tau$, respectively. The parameter $\tau$ is given by $\tau =N{\left[J_{\mathrm{ij}}{J}_{\mathrm{ji}}\right]}_J$ and $N{\left[J_{\mathrm{ij}}^2\right]}_J$ is normalized to 1. In the $\tau =1\mathrm{}$ limit, Wigner's semicircle law is recovered. The results are extended to complex asymmetric matrices.