Strong reproductive isolation among sympatric lineages of an international mosquito pest

Using population genomics and a chromosome-level assembly, researchers identified three cryptic, strongly reproductively isolated lineages of the mosquito *Aedes notoscriptus*, revealing that a single lineage (VIC1) drove international invasions while coexisting with a distinct lineage (VIC2) in Australia that exhibits hybridization barriers and a massive structural variant.

The Gist

The Big Picture: A Mosquito Mystery

Imagine you walk into a room and see a group of people who all look exactly the same. You assume they are all from the same family. But then, you discover that hidden inside their DNA, they are actually three completely different families who have been living in the same house for a long time without mixing much.

That is exactly what scientists found with the Australian backyard mosquito (Aedes notoscriptus). This mosquito is a pest that carries diseases like Ross River virus and even the bacteria that causes Buruli ulcer. For years, scientists thought there was just one type of this mosquito. This study used high-tech "genomic detective work" to reveal that there are actually three distinct, hidden lineages (let's call them Family A, Family B, and Family C).

The Three Families

1. Family A (VIC1): This is the "Global Traveler." It lives in Victoria (Australia) but is also the one that escaped to invade New Zealand and California. It's the one causing trouble internationally.
2. Family B (VIC2): This is the "Local Resident." It also lives in Victoria, right next door to Family A. It has been there a long time and has a very diverse family tree.
3. Family C (NT): This family lives way up north in the Northern Territory and is quite different from the others.

The big surprise? Family A and Family B are living in the same neighborhoods in Melbourne. They are "sympatric," meaning they share the same backyard. You'd expect them to mix, have babies together, and become one big blended family. But they aren't.

The "No-Go" Zone: Why They Don't Mix

Even though Family A and Family B live side-by-side, they act like two different species that refuse to date.

• The Analogy: Imagine two different high schools in the same town. The students walk past each other every day, but they rarely hang out. If they do, the "mixed" kids (hybrids) are rare.
• The Science: The researchers found that when they looked at the mosquitoes, there were far fewer "mixed" babies than they expected if the two families were mating randomly. It's like if you flipped a coin 100 times and expected 50 heads and 50 tails, but you only got 10 heads. Something is stopping them from mixing.
• The Reason: It's likely a mix of pre-mating barriers (they might dance to different rhythms or look for food at different times) and post-mating barriers (if they do mate, the babies might be weak or sterile).

The "Genetic Bottleneck": A Family That Lost Its Way

Family A (the one that invaded California and New Zealand) has a very interesting history.

• The Analogy: Imagine a family that used to be huge and diverse, but then a disaster happened, and only a few members survived to start a new town. That new town has very little genetic variety compared to the original home.
• The Science: Family A has much less genetic diversity than Family B. This suggests that Family A went through a "bottleneck"—a population crash—before expanding into Melbourne and then jumping to the US and New Zealand. They are essentially the "survivors" of a genetic crash.

The "Super-Block" on the Chromosome

Here is the coolest discovery. Inside Family B's DNA, there is a massive chunk of genetic code (about 25 million letters long, which is huge for a mosquito!) that is flipped upside down compared to Family A.

• The Analogy: Imagine two instruction manuals for building a house. Family A has the manual printed normally. Family B has a 50-page section of their manual printed upside down.
• The Effect: Because this section is flipped, the genes inside it get "locked" together. They can't swap parts with Family A easily. This "super-block" contains genes that control how the mosquito finds hosts, how it reacts to stress, and how it behaves. This locked-up section might be the secret weapon that keeps the two families separate and helps them survive in different ways.

Why Should You Care?

This isn't just about mosquito gossip; it matters for public health.

1. Disease Control: If you try to stop mosquitoes using a new method (like releasing sterile males or bacteria), you need to know which "family" you are targeting. If you target Family A but miss Family B, the disease could still spread.
2. Invasion Tracking: We now know that the mosquitoes invading California and New Zealand are specifically from Family A. This helps us trace where they came from and how they travel.
3. The "One Size Fits All" Myth: This study proves that just because two bugs look the same, they aren't the same. Treating them as one group could be a mistake in fighting disease.

The Bottom Line

Scientists built a brand-new, high-definition map of the mosquito's genome (like upgrading from a blurry sketch to a 4K photo). Using this map, they discovered that the "Australian backyard mosquito" is actually a complex of hidden families. Two of these families live together in Melbourne but refuse to mix, likely due to a massive genetic "flip" in their DNA. Understanding these differences is the key to stopping the diseases they carry.

Technical Summary

1. Problem Statement

The Australian mosquito Aedes notoscriptus is a significant public health concern, acting as a vector for arboviruses (Ross River, Barmah Forest), dog heartworm, and Mycobacterium ulcerans (Buruli ulcer). While previous mitochondrial DNA (mtDNA) studies suggested the existence of cryptic lineages within this species, the extent of reproductive isolation, nuclear genomic structure, and the mechanisms driving lineage divergence remained unclear. Specifically, it was unknown whether these lineages represented distinct species, how they interact in sympatry (specifically in Greater Melbourne), and what role they play in international invasions (New Zealand and California). The reliance on mtDNA markers was suspected to be insufficient due to potential mito-nuclear discordance.

2. Methodology

The study employed a comprehensive population genomics approach coupled with high-quality reference genome assembly:

• Reference Genome Assembly: The authors generated a chromosome-level reference genome for Ae. notoscriptus using PacBio HiFi long-read sequencing (52.1 Gb) and Hi-C chromatin conformation capture data. The assembly spans ~900 Mb with three primary chromosome scaffolds.
• Sampling: Genome-wide sequence data (ddRADseq) was generated for 193 individuals collected from the native range (Victoria, NSW, Queensland, Northern Territory, Australia) and invasive populations (New Zealand, California, USA).
• Population Structure Analysis:
  • PCA and DAPC: Used to identify genetic clusters and visualize differentiation.
  • Phylogenetics: Constructed using genome-wide absolute divergence ($d_{XY}$) to minimize bias from recent drift.
  • Hybridization Detection: Utilized ADMIXTURE, Treemix (for migration edges), and fineRADstructure (for microhaplotype similarity) to detect admixture.
  • Demographic History: Analyzed autosomal heterozygosity and Tajima's D to infer population bottlenecks and expansion.
  • Structural Variant Detection: Employed local PCA (using the lostruct package) to identify genomic regions with abnormal population structure indicative of inversions.
  • Functional Annotation: Candidate genes within structural variants were annotated via BLAST against Ae. aegypti and analyzed for Gene Ontology (GO) enrichment.
• Additional Screening: All individuals were screened for Wolbachia infection (via qPCR and alignment) to rule out cytoplasmic incompatibility as a reproductive barrier. mtDNA COI gene sequencing was performed to compare nuclear and mitochondrial phylogenies.

3. Key Contributions

• First Chromosome-Level Assembly: Provided the first high-quality, chromosome-scale reference genome for Ae. notoscriptus, a critical resource for future evolutionary and vector control research.
• Discovery of Cryptic Lineages: Identified three deeply divergent nuclear lineages (VIC1, VIC2, and NT) that are discordant with previously defined mtDNA clades, highlighting the limitations of mitochondrial markers for this species.
• Quantification of Reproductive Isolation: Demonstrated that despite extensive sympatry in Greater Melbourne, two major lineages (VIC1 and VIC2) exhibit strong but incomplete reproductive isolation, with significantly fewer hybrids than expected under random mating.
• Structural Variant Discovery: Uncovered a massive structural variant (~25 Mbp, ~10% of a chromosome) segregating in the VIC2 lineage, likely an inversion, which is fixed in VIC1.

4. Key Results

• Lineage Identification:
  • VIC1: The lineage responsible for all international invasions (NZ and USA). It is widespread across southern Australia but shows signs of a recent demographic bottleneck (reduced heterozygosity, elevated Tajima's D) suggesting a recent range expansion into Greater Melbourne.
  • VIC2: A distinct lineage co-occurring with VIC1 in Greater Melbourne. It retains higher genetic diversity, suggesting it represents the ancestral population in the region.
  • NT: A third lineage restricted to the Northern Territory.
• Mito-Nuclear Discordance: There is a lack of concordance between nuclear and mitochondrial phylogenies. mtDNA clades are distributed across nuclear lineages, likely due to historical hybridization and backcrossing followed by selection against nuclear introgression.
• Reproductive Isolation:
  • In the 11 sites where VIC1 and VIC2 co-occur, hybrids were observed at a frequency of 22%, significantly lower than the 41% expected under random mating ($p = 0.005$).
  • Only 11 individuals were identified as recent hybrids (F1 or early generation) based on high heterozygosity and admixture signals, suggesting strong selection against hybrids or pre-mating barriers.
  • No Wolbachia infections were found, ruling out cytoplasmic incompatibility as the driver of isolation.
• Demographic History: VIC1 shows a ~21.5% reduction in heterozygosity compared to VIC2 and elevated Tajima's D, consistent with a recent bottleneck and range expansion into Melbourne. VIC2 shows no such signature, indicating long-term stability in the region.
• Structural Variant (Inversion):
  • A ~25 Mbp region on Chromosome 2 was identified as a putative inversion segregating in VIC2 but fixed in VIC1.
  • This region contains ~112 genes enriched for neuronal signaling (ligand-gated ion channels, neurotransmitter receptors), oxidoreductase activity (stress response), and chromatin regulation (Polycomb complex).
  • The variant is not geographically clinal within Melbourne but is fixed in VIC1 and segregating in VIC2, suggesting it may act as a "supergene" maintaining adaptive differences between the lineages.
• Dispersal: Neither lineage showed significant isolation-by-distance within Greater Melbourne, indicating high effective dispersal and near-panmixia within each lineage, despite the reproductive barriers between them.

5. Significance

• Vector Control Implications: The discovery of cryptic lineages with different demographic histories and potential behavioral/physiological differences (suggested by the inversion's gene content) implies that vector competence, insecticide resistance, and dispersal capabilities may vary between lineages. Control strategies (e.g., Wolbachia releases or gene drives) must account for these distinct populations.
• Evolutionary Insight: The study provides a compelling model for studying reproductive isolation in sympatry. It demonstrates that strong reproductive barriers can evolve and persist even in the absence of geographic separation, potentially driven by large-scale structural variants (inversions) that suppress recombination and maintain adaptive allele combinations.
• Methodological Advancement: The work underscores the necessity of genome-wide nuclear data over mtDNA for accurate species delineation in mosquitoes and establishes a foundational genomic resource for Ae. notoscriptus to support future comparative genomics and epidemiological surveillance.