Parallel evolution in mosquito vectors - a duplicated esterase locus is associated with resistance to pirimiphos-methyl in An. gambiae.

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in Anopheles gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to organophosphates in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in An. gambiae, An. coluzzii and An. arabiensis. As in Cx. pipiens, copy number variants have arisen at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in Anopheles gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programmes.

Discovery of knock-down resistance in the major African malaria vector Anopheles funestus.

A major mechanism of insecticide resistance in insect pests is knock-down resistance (kdr) caused by mutations in the voltage-gated sodium channel (Vgsc) gene. Despite being common in most malaria Anopheles vector species, kdr mutations have never been observed in Anopheles funestus, the principal malaria vector in Eastern and Southern Africa. While monitoring 10 populations of An. funestus in Tanzania, we unexpectedly found resistance to DDT, a banned insecticide, in one location. Through whole-genome sequencing of 333 An. funestus samples from these populations, we found 8 novel amino acid substitutions in the Vgsc gene, including the kdr variant, L976F (L1014F in An. gambiae), in tight linkage disequilibrium with another (P1842S). The mutants were found only at high frequency in one region, with a significant decline between 2017 and 2023. Notably, kdr L976F was strongly associated with survivorship to the exposure to DDT insecticide, while no clear association was noted with a pyrethroid insecticide (deltamethrin). Further study is necessary to identify the origin and spread of kdr in An. funestus, and the potential threat to current insecticide-based vector control in Africa.

Copy number variants underlie the major selective sweeps in insecticide resistance genes in Anopheles arabiensis from Tanzania.

To keep ahead of the evolution of resistance to insecticides in mosquitoes, national malaria control programmes must make use of a range of insecticides, both old and new, while monitoring resistance mechanisms. Knowledge of the mechanisms of resistance remains limited in Anopheles arabiensis, which in many parts of Africa is of increasing importance because it is apparently less susceptible to many indoor control interventions. Furthermore, comparatively little is known in general about resistance to non-pyrethroid insecticides such as pirimiphos-methyl (PM), which are crucial for effective control in the context of resistance to pyrethroids. We performed a genome-wide association study to determine the molecular mechanisms of resistance to deltamethrin (commonly used in bednets) and PM, in An. arabiensis from two regions in Tanzania. Genomic regions of positive selection in these populations were largely driven by copy number variants (CNVs) in gene families involved in resistance to these two insecticides. We found evidence of a new gene cluster involved in resistance to PM, identifying a strong selective sweep tied to a CNV in the Coeae2g-Coeae6g cluster of carboxylesterase genes. Using complementary data from An. coluzzii in Ghana, we show that copy number at this locus is significantly associated with PM resistance. Similarly, for deltamethrin, resistance was strongly associated with a novel CNV allele in the Cyp6aa / Cyp6p cluster. Against this background of metabolic resistance, target site resistance was very rare or absent for both insecticides. Mutations in the pyrethroid target site Vgsc were at very low frequency in Tanzania, yet combining these samples with three An. arabiensis individuals from West Africa revealed a startling diversity of evolutionary origins of target site resistance, with up to 5 independent origins of Vgsc-995 mutations found within just 8 haplotypes. Thus, despite having been first recorded over 10 years ago, Vgsc resistance mutations in Tanzanian An. arabiensis have remained at stable low frequencies. Overall, our results provide a new copy number marker for monitoring resistance to PM in malaria mosquitoes, and reveal the complex picture of resistance patterns in An. arabiensis.

Parallel evolution in mosquito vectors - a duplicated esterase locus is associated with resistance to pirimiphos-methyl in An. gambiae.

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in Anopheles gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to organophosphates in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in An. gambiae, An. coluzzii and An. arabiensis. As in Cx. pipiens, copy number variation seems to play a role in the evolution of insecticide resistance at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in Anopheles gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programmes.

Adaptive Long-Read Sequencing Reveals GGC Repeat Expansion in ZFHX3 Associated with Spinocerebellar Ataxia Type 4.

BACKGROUND: Spinocerebellar ataxia type 4 (SCA4) is an autosomal dominant ataxia with invariable sensory neuropathy originally described in a family with Swedish ancestry residing in Utah more than 25 years ago. Despite tight linkage to the 16q22 region, the molecular diagnosis has since remained elusive. OBJECTIVES: Inspired by pathogenic structural variation implicated in other 16q-ataxias with linkage to the same locus, we revisited the index SCA4 cases from the Utah family using novel technologies to investigate structural variation within the candidate region. METHODS: We adopted a targeted long-read sequencing approach with adaptive sampling on the Oxford Nanopore Technologies (ONT) platform that enables the detection of segregating structural variants within a genomic region without a priori assumptions about any variant features. RESULTS: Using this approach, we found a heterozygous (GGC)n repeat expansion in the last coding exon of the zinc finger homeobox 3 (ZFHX3) gene that segregates with disease, ranging between 48 and 57 GGC repeats in affected probands. This finding was replicated in a separate family with SCA4. Furthermore, the estimation of this GGC repeat size in short-read whole genome sequencing (WGS) data of 21,836 individuals recruited to the 100,000 Genomes Project in the UK and our in-house dataset of 11,258 exomes did not reveal any pathogenic repeats, indicating that the variant is ultrarare. CONCLUSIONS: These findings support the utility of adaptive long-read sequencing as a powerful tool to decipher causative structural variation in unsolved cases of inherited neurological disease. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Early human occupation of Australia's eastern seaboard.

Secure archaeological evidence for human occupation on the eastern seaboard of Australia before ~ 25,000 years ago has proven elusive. This has prompted some researchers to argue that the coastal margins remained uninhabited prior to 25 ka. Here we show evidence for human occupation beginning between 30 ± 6 and 49 ± 8 ka at Wallen Wallen Creek (WWC), and at Middle Canalpin Creek (MCA20) between 38 ± 8 and 41 ± 8 ka. Both sites are located on the western side of Minjerribah (North Stradbroke Island), the second largest sand island in the world, isolated by rising sea levels in the early Holocene. The earliest occupation phase at both sites consists of charcoal and heavily retouched stone artefacts made from exotic raw materials. Heat-treatment of imported silcrete artefacts first appeared in sediment dated to ~ 30,000 years ago, making these amongst Australia's oldest dated heat-treated artefacts. An early human presence on Minjerribah is further suggested by palaeoenvironmental records of anthropogenic burning beginning by 45,000 years ago. These new chronologies from sites on a remnant portion of the continental margin confirm early human occupation along Sahul's now-drowned eastern continental shelf.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression.

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis-regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae and Anopheles coluzzii.

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of malaria control, but the genetics of resistance are only partially understood. We performed a large scale multi-country genome-wide association study of resistance to two widely used insecticides: deltamethrin and pirimiphos-methyl, using sequencing data from An. gambiae and An. coluzzii from ten locations in West Africa. Resistance was highly multi-genic, multi-allelic and variable between populations. While the strongest and most consistent association with deltamethrin resistance came from Cyp6aa1, this was based on several independent copy number variants (CNVs) in An. coluzzii, and on a non-CNV haplotype in An. gambiae. For pirimiphos-methyl, signals included Ace1, cytochrome P450s, glutathione S-transferases and the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes showed evidence of cross-resistance to both insecticides. These locally-varying, multi-allelic patterns highlight the challenges involved in genomic monitoring of resistance, and may form the basis for improved surveillance methods.

Druggable proteins influencing cardiac structure and function: Implications for heart failure therapies and cancer cardiotoxicity.

Dysfunction of either the right or left ventricle can lead to heart failure (HF) and subsequent morbidity and mortality. We performed a genome-wide association study (GWAS) of 16 cardiac magnetic resonance (CMR) imaging measurements of biventricular function and structure. Cis-Mendelian randomization (MR) was used to identify plasma proteins associating with CMR traits as well as with any of the following cardiac outcomes: HF, non-ischemic cardiomyopathy, dilated cardiomyopathy (DCM), atrial fibrillation, or coronary heart disease. In total, 33 plasma proteins were prioritized, including repurposing candidates for DCM and/or HF: IL18R (providing indirect evidence for IL18), I17RA, GPC5, LAMC2, PA2GA, CD33, and SLAF7. In addition, 13 of the 25 druggable proteins (52%; 95% confidence interval, 0.31 to 0.72) could be mapped to compounds with known oncological indications or side effects. These findings provide leads to facilitate drug development for cardiac disease and suggest that cardiotoxicities of several cancer treatments might represent mechanism-based adverse effects.

Genome-wide association studies reveal novel loci associated with pyrethroid and organophosphate resistance in Anopheles gambiae s.l.

Resistance to insecticides in Anopheles mosquitoes threatens the effectiveness of the most widespread tools currently used to control malaria. The genetic underpinnings of resistance are still only partially understood, with much of the variance in resistance phenotype left unexplained. We performed a multi-country large scale genome-wide association study of resistance to two insecticides widely used in malaria control: deltamethrin and pirimiphos-methyl. Using a bioassay methodology designed to maximise the phenotypic difference between resistant and susceptible samples, we sequenced 969 phenotyped female An. gambiae and An. coluzzii from ten locations across four countries in West Africa (Benin, Côte d'Ivoire, Ghana and Togo), identifying single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) segregating in the populations. The patterns of resistance association were highly multiallelic and variable between populations, with different genomic regions contributing to resistance, as well as different mutations within a given region. While the strongest and most consistent association with deltamethrin resistance came from the region around Cyp6aa1 , this resistance was based on a combination of several independent CNVs in An. coluzzii , and on a non-CNV bearing haplotype in An. gambiae . Further signals involved a range of cytochrome P450, mitochondrial, and immunity genes. Similarly, for pirimiphos-methyl, while the strongest signal came from the region of Ace1 , more widespread signals included cytochrome P450s, glutathione S-transferases, and a subunit of the nAChR target site of neonicotinoid insecticides. The regions around Cyp9k1 and the Tep family of immune genes were associated with resistance to both insecticide classes, suggesting possible cross-resistance mechanisms. These locally-varying, multigenic and multiallelic patterns highlight the challenges involved in genomic monitoring and surveillance of resistance, and form the basis for improvement of methods used to detect and predict resistance. Based on simulations of resistance variants, we recommend that yet larger scale studies, exceeding 500 phenotyped samples per population, are required to better identify associated genomic regions.

65,000-years of continuous grinding stone use at Madjedbebe, Northern Australia.

Grinding stones and ground stone implements are important technological innovations in later human evolution, allowing the exploitation and use of new plant foods, novel tools (e.g., bone points and edge ground axes) and ground pigments. Excavations at the site of Madjedbebe recovered Australia's (if not one of the world's) largest and longest records of Pleistocene grinding stones, which span the past 65 thousand years (ka). Microscopic and chemical analyses show that the Madjedbebe grinding stone assemblage displays the earliest known evidence for seed grinding and intensive plant use, the earliest known production and use of edge-ground stone hatchets (aka axes), and the earliest intensive use of ground ochre pigments in Sahul (the Pleistocene landmass of Australia and New Guinea). The Madjedbebe grinding stone assemblage reveals economic, technological and symbolic innovations exemplary of the phenotypic plasticity of Homo sapiens dispersing out of Africa and into Sahul.

Stone toolmaking difficulty and the evolution of hominin technological skills.

Stone tools are a manifestation of the complex cognitive and dexterous skills of our hominin ancestors. As such, much research has been devoted to understanding the skill requirements of individual lithic technologies. Yet, comparing skill across different technologies, and thus across the vast timespan of the Palaeolithic, is an elusive goal. We seek to quantify a series of commensurable metrics of knapping skill across four different lithic technologies (discoids, handaxes, Levallois, and prismatic blades). To compare the requisite dexterity, coordination, and care involved in each technology, we analysed video footage and lithic material from a series of replicative knapping experiments to quantify deliberation (strike time), precision (platform area), intricacy (flake size relative to core size), and success (relative blank length). According to these four metrics, discoidal knapping appears to be easiest among the sample. Levallois knapping involved an intricate reduction sequence, but did not require as much motor control as handaxes and especially prismatic blades. Compared with the other Palaeolithic technologies, we conclude that prismatic blade knapping is set apart by being a skill intensive means of producing numerous standardised elongate end-products.

Stochastic models support rapid peopling of Late Pleistocene Sahul.

The peopling of Sahul (the combined continent of Australia and New Guinea) represents the earliest continental migration and settlement event of solely anatomically modern humans, but its patterns and ecological drivers remain largely conceptual in the current literature. We present an advanced stochastic-ecological model to test the relative support for scenarios describing where and when the first humans entered Sahul, and their most probable routes of early settlement. The model supports a dominant entry via the northwest Sahul Shelf first, potentially followed by a second entry through New Guinea, with initial entry most consistent with 50,000 or 75,000 years ago based on comparison with bias-corrected archaeological map layers. The model's emergent properties predict that peopling of the entire continent occurred rapidly across all ecological environments within 156-208 human generations (4368-5599 years) and at a plausible rate of 0.71-0.92 km year-1. More broadly, our methods and approaches can readily inform other global migration debates, with results supporting an exit of anatomically modern humans from Africa 63,000-90,000 years ago, and the peopling of Eurasia in as little as 12,000-15,000 years via inland routes.

The genetic architecture of target-site resistance to pyrethroid insecticides in the African malaria vectors Anopheles gambiae and Anopheles coluzzii.

Resistance to pyrethroid insecticides is a major concern for malaria vector control. Pyrethroids target the voltage-gated sodium channel (VGSC), an essential component of the mosquito nervous system. Substitutions in the amino acid sequence can induce a resistance phenotype. We use whole-genome sequence data from phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) to provide a comprehensive account of genetic variation in the Vgsc gene across 13 African countries. In addition to known resistance alleles, we describe 20 other non-synonymous nucleotide substitutions at appreciable population frequency and map these variants onto a protein model to investigate the likelihood of pyrethroid resistance phenotypes. Thirteen of these novel alleles were found to occur almost exclusively on haplotypes carrying the known L995F kdr (knock-down resistance) allele and may enhance or compensate for the L995F resistance genotype. A novel mutation I1527T, adjacent to a predicted pyrethroid-binding site, was found in tight linkage with V402L substitutions, similar to allele combinations associated with resistance in other insect species. We also analysed genetic backgrounds carrying resistance alleles, to determine which alleles have experienced recent positive selection, and describe ten distinct haplotype groups carrying known kdr alleles. Five of these groups are observed in more than one country, in one case separated by over 3000 km, providing new information about the potential for the geographical spread of resistance. Our results demonstrate that the molecular basis of target-site pyrethroid resistance in malaria vectors is more complex than previously appreciated, and provide a foundation for the development of new genetic tools for insecticide resistance management.

Pandanus nutshell generates a palaeoprecipitation record for human occupation at Madjedbebe, northern Australia.

Little is known about the Pleistocene climatic context of northern Australia at the time of early human settlement. Here we generate a palaeoprecipitation proxy using stable carbon isotope analysis of modern and archaeological pandanus nutshell from Madjedbebe, Australia's oldest known archaeological site. We document fluctuations in precipitation over the last 65,000 years and identify periods of lower precipitation during the penultimate and last glacial stages, Marine Isotope Stages 4 and 2. However, the lowest effective annual precipitation is recorded at the present time. Periods of lower precipitation, including the earliest phase of occupation, correspond with peaks in exotic stone raw materials and artefact discard at the site. This pattern is interpreted as suggesting increased group mobility and intensified use of the region during drier periods.

Papuan mitochondrial genomes and the settlement of Sahul.

New Guineans represent one of the oldest locally continuous populations outside Africa, harboring among the greatest linguistic and genetic diversity on the planet. Archeological and genetic evidence suggest that their ancestors reached Sahul (present day New Guinea and Australia) by at least 55,000 years ago (kya). However, little is known about this early settlement phase or subsequent dispersal and population structuring over the subsequent period of time. Here we report 379 complete Papuan mitochondrial genomes from across Papua New Guinea, which allow us to reconstruct the phylogenetic and phylogeographic history of northern Sahul. Our results support the arrival of two groups of settlers in Sahul within the same broad time window (50-65 kya), each carrying a different set of maternal lineages and settling Northern and Southern Sahul separately. Strong geographic structure in northern Sahul remains visible today, indicating limited dispersal over time despite major climatic, cultural, and historical changes. However, following a period of isolation lasting nearly 20 ky after initial settlement, environmental changes postdating the Last Glacial Maximum stimulated diversification of mtDNA lineages and greater interactions within and beyond Northern Sahul, to Southern Sahul, Wallacea and beyond. Later, in the Holocene, populations from New Guinea, in contrast to those of Australia, participated in early interactions with incoming Asian populations from Island Southeast Asia and continuing into Oceania.

Human occupation of northern India spans the Toba super-eruption ~74,000 years ago.

India is located at a critical geographic crossroads for understanding the dispersal of Homo sapiens out of Africa and into Asia and Oceania. Here we report evidence for long-term human occupation, spanning the last ~80 thousand years, at the site of Dhaba in the Middle Son River Valley of Central India. An unchanging stone tool industry is found at Dhaba spanning the Toba eruption of ~74 ka (i.e., the Youngest Toba Tuff, YTT) bracketed between ages of 79.6 ± 3.2 and 65.2 ± 3.1 ka, with the introduction of microlithic technology ~48 ka. The lithic industry from Dhaba strongly resembles stone tool assemblages from the African Middle Stone Age (MSA) and Arabia, and the earliest artefacts from Australia, suggesting that it is likely the product of Homo sapiens as they dispersed eastward out of Africa.

The first Australian plant foods at Madjedbebe, 65,000-53,000 years ago.

There is little evidence for the role of plant foods in the dispersal of early modern humans into new habitats globally. Researchers have hypothesised that early movements of human populations through Island Southeast Asia and into Sahul were driven by the lure of high-calorie, low-handling-cost foods, and that the use of plant foods requiring processing was not common in Sahul until the Holocene. Here we present the analysis of charred plant food remains from Madjedbebe rockshelter in northern Australia, dated to between 65 kya and 53 kya. We demonstrate that Australia's earliest known human population exploited a range of plant foods, including those requiring processing. Our finds predate existing evidence for such subsistence practices in Sahul by at least 23ky. These results suggest that dietary breadth underpinned the success of early modern human populations in this region, with the expenditure of labour on the processing of plants guaranteeing reliable access to nutrients in new environments.

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes.

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.