DOI: 10.5061/dryad.3f355
title: Data from: Mitochondrial DNA variants help monitor the dynamics of Wolbachia invasion into host populations

description: Wolbachia is the most widespread endosymbiotic bacterium of insects and
 other arthropods that can rapidly invade host populations. Deliberate
 releases of Wolbachia into natural populations of the dengue fever
 mosquito, Aedes aegypti, are used as a novel biocontrol strategy for
 dengue suppression. Invasion of Wolbachia through the host population
 relies on factors such as high fidelity of the endosymbiont transmission
 and limited immigration of uninfected individuals, but these factors can
 be difficult to measure. One way of acquiring relevant information is to
 consider mitochondrial DNA (mtDNA) variation alongside Wolbachia in
 field-caught mosquitoes. Here we used diagnostic mtDNA markers to
 differentiate infection-associated mtDNA haplotypes from those of the
 uninfected mosquitoes at release sites. Unique haplotypes associated with
 Wolbachia were found at locations outside Australia. We also performed
 mathematical and qualitative analyses including modelling the expected
 dynamics of the Wolbachia and mtDNA variants during and after a release.
 Our analyses identified key features in haplotype frequency patterns to
 infer the presence of imperfect maternal transmission of Wolbachia,
 presence of immigration and possibly incomplete cytoplasmic
 incompatibility. We demonstrate that ongoing screening of the mtDNA
 variants should provide information on maternal leakage and immigration,
 particularly in releases outside Australia. As we demonstrate in a case
 study, our models to track the Wolbachia dynamics can be successfully
 applied to temporal studies in natural populations or Wolbachia release
 programs, as long as there is co-occurring mtDNA variation that
 differentiates infected and uninfected populations.
description: RADseq_Loci_in_all_individualsSequences and hamming distance from assembled mitochondrial sequences in RAD sequencingNucleotide alignment_CO1Fasta file of partial sequences. Description of samples can be obtained from Excel spreadsheet provided. See Column G of "Sample_details_updated.xlsx".Nucleotide alignment_ND4Partial ND4 sequences. Labels match column B of "Sample_details_updated.xlsx". See comments in cell H1 for explanation of sequences labelled "...Assembly..."Nucleotide alignment_ND5Partial ND5 sequences. See column F of "Sample_details_updated.xlsx" to trace back actual sample details.Primer1Sequence from region amplified by primer pair 1 (see Behura et al. 2011). Labels match column B of "Sample_details_updated.xlsx". See comments in cell H1 for explanation of sequences labelled "...Assembly..."Primer4Sequence from region amplified by primer pair 4 (see Behura et al. 2011). Labels match column B of "Sample_details_updated.xlsx". See comments in cell H1 for explanation of sequences labelled "...Assembly..."Primer6Sequence from region amplified by primer pair 6 (see Behura et al. 2011). Labels match column B of "Sample_details_updated.xlsx". See comments in cell H1 for explanation of sequences labelled "...Assembly..."R-scriptR codes for modelsSample_details(relevant columns for ID). Column B lists all ID of samples that were considered, but not all were sequenced. Column F lists the labels in the ND5 fasta file. Column G lists the labels in the COI fasta file.