Geographic differentiation in the house fly estimated by microsatellite and mitochondrial variation

Integrated genome-wide investigations of the housefly, a global vector of diseases reveal unique dispersal patterns and bacterial communities across farms

Background

Houseflies (Musca domestica L.) live in intimate association with numerous microorganisms and is a vector of human pathogens. In temperate areas, houseflies will overwinter in environments constructed by humans and recolonize surrounding areas in early summer. However, the dispersal patterns and associated bacteria across season and location are unclear. We used genotyping-by-sequencing (GBS) for the simultaneous identification and genotyping of thousands of Single Nucleotide Polymorphisms (SNPs) to establish dispersal patterns of houseflies across farms. Secondly, we used 16S rRNA gene amplicon sequencing to establish the variation and association between bacterial communities and the housefly across farms.

Results

Using GBS we identified 18,000 SNPs across 400 individuals sampled within and between 11 dairy farms in Denmark. There was evidence for sub-structuring of Danish housefly populations and with genetic structure that differed across season and sex. Further, there was a strong isolation by distance (IBD) effect, but with large variation suggesting that other hidden geographic barriers are important. Large individual variations were observed in the community structure of the microbiome and it was found to be dependent on location, sex, and collection time. Furthermore, the relative prevalence of putative pathogens was highly dependent on location and collection time.

Conclusion

We were able to identify SNPs for the determination of the spatiotemporal housefly genetic structure, and to establish the variation and association between bacterial communities and the housefly across farms using novel next-generation sequencing (NGS) techniques. These results are important for disease prevention given the fine-scale population structure and IBD for the housefly, and that individual houseflies carry location specific bacteria including putative pathogens.

Is Multifactorial Sex Determination in the House Fly, Musca domestica (L.), Stable Over Time?

Sex determination pathways evolve rapidly, usually because of turnover of master regulatory genes at the top of the developmental pathway. Polygenic sex determination is expected to be a transient state between ancestral and derived conditions. However, polygenic sex determination has been observed in numerous animal species, including the house fly, Musca domestica House fly males carry a male-determining factor (M) that can be located on any chromosome, and an individual male may have multiple M factors. Females lack M and/or have a dominant allele of the Md-tra gene (Md-tra ^D ) that acts as a female-determining locus even in the presence of multiple copies of M. We found the frequency and linkage of M in house flies collected in Chino, CA (USA) was relatively unchanged between 1982 and 2014. The frequency of females with Md-tra ^D in the 2014 collection was 33.6% (n = 140). Analysis of these results, plus previously published data, revealed a strong correlation between the frequencies of Md-tra ^D and multiple M males, and we find that these populations are expected to have balanced sex ratios. We also find that fitness values that allow for the invasion and maintenance of multiple sex determining loci suggest that sexually antagonistic selection could be responsible for maintaining polygenic sex determination in house fly populations. The stability over time and equilibrium frequencies within populations suggest the house fly polygenic sex determination system is not in transition, and provide guidance for future investigations on the factors responsible for the polymorphism.

Mitochondrial genetic variations in natural house fly (Musca domestica L.) populations from the western and southern parts of Turkey

The house fly Musca domestica Linnaeus (Diptera) is one of the most studied species that is globally distributed and well known to everyone. In order to ensure baseline knowledge for the genetic resources of the species, genetic variation in M. domestica populations from western and southern parts of Turkey was investigated using nucleotide sequence analysis of 348 base pairs (bp) in the mitochondrial cytochrome oxidase subunit I gene (COI). Samples of 192 individuals were collected from 16 localities of Turkey. There were 10 variable sites defining two haplotypes of COI in this species. There was no difference in geographical distribution frequency between the two regions of Turkey. Overall, haplotype diversity (h) was low, ranging from 0 to 0.5606 with the average overall value of 0.178 ± 0.04 and nucleotide diversity (π), ranged from 0 to 0.0056 with the overall mean of 0.0016. Analysis of molecular variance (AMOVA) indicated that genetic differentiation within individuals and populations was low and significant (p < 0.05). Except Afyon population, conventional population statistic FST showed no significant genetic structure along the range of M. domestica populations. Sixteen populations clustered under six haplotypes and two of them are unique to Turkey. Haplotype networks suggested that house fly populations in Turkey are grouped with the Palearctic region, which is the most probable place for the origin of this species.

Phenotypic variation of the housefly, Musca domestica: amounts and patterns of wing shape asymmetry in wild populations and laboratory colonies

Musca domestica L. (Diptera: Muscidae) is a vector of a range variety of pathogens infecting humans and animals. During a year, housefly experiences serial population bottlenecks resulted in reduction of genetic diversity. Population structure has also been subjected to different selection regimes created by insect control programs and pest management. Both environmental and genetic disturbances can affect developmental stability, which is often reflected in morphological traits as asymmetry. Since developmental stability is of great adaptive importance, the aim of this study was to examine fluctuating asymmetry (FA), as a measure of developmental instability, in both wild populations and laboratory colonies of M. domestica. The amount and pattern of wing shape FA was compared among samples within each of two groups (laboratory and wild) and between groups. Firstly, the amount of FA does not differ significantly among samples within the group and neither does it differ between groups. Regarding the mean shape of FA, contrary to non-significant difference within the wild population group and among some colonies, the significant difference between groups was found. These results suggest that the laboratory colonies and wild samples differ in buffering mechanisms to perturbations during development. Hence, inbreeding and stochastic processes, mechanisms dominating in the laboratory-bred samples contributed to significant changes in FA of wing shape. Secondly, general patterns of left-right displacements of landmarks across both studied sample groups are consistent. Observed consistent direction of FA implies high degrees of wing integration. Thus, our findings shed light on developmental buffering processes important for population persistence in the environmental change and genetic stress influence on M. domestica.

Multiple origins of kdr-type resistance in the house fly, Musca domestica

Insecticide resistance is a model phenotype that can be used to investigate evolutionary processes underlying the spread of alleles across a global landscape, while offering valuable insights into solving the problems that resistant pests present to human health and agriculture. Pyrethroids are one of the most widely used classes of insecticides world-wide and they exert their toxic effects through interactions with the voltage-sensitive sodium channel (Vssc). Specific mutations in Vssc (kdr, kdr-his and super-kdr) are known to cause resistance to pyrethroid insecticides in house flies. In order to determine the number of evolutionary origins of kdr, kdr-his and super-kdr, we sequenced a region of Vssc from house flies collected in the USA, Turkey and China. Our phylogenetic analysis of Vssc unequivocally supports the hypothesis of multiple independent origins of kdr, super-kdr and kdr-his on an unprecedented geographic scale. The implications of these evolutionary processes on pest management are discussed.

Molecular evidence to suggest the origin of a colonization: Drosophila subobscura in America

The recent colonization of America by Drosophila subobscura represents a great opportunity for evolutionary biology studies. Knowledge of the populations from which the colonization started would provide an understanding of how genetic composition changed during adaptation to the new environment. Thus, a 793 nucleotide fragment of the Odh (Octanol dehydrogenase) gene was sequenced in 66 chromosomal lines from Barcelona (western Mediterranean) and in 66 from Mt. Parnes (Greece, eastern Mediterranean). No sequence of Odh fragment in Barcelona or Mt. Parnes was identical to any of those previously detected in America. However, an Odh sequence from Barcelona differed in only one nucleotide from another found in American populations. In both cases, the chromosomal lines presented the same inversion: O(7), and the Odh gene was located within this inversion. This evidence suggests a possible western Mediterranean origin for the colonization. Finally, the molecular and inversion data indicate that the colonization was not characterized by multiple reintroductions.

Phylogenetic characterization of bacteria in the gut of house flies (Musca domestica L.)

House flies (Musca domestica L.) are cosmopolitan, ubiquitous, synanthropic insects that serve as mechanical or biological vectors for various microorganisms. To fully assess the role of house flies in the epidemiology of human diseases, it is essential to understand the diversity of microbiota harbored by natural fly populations. This study aimed to identify the diversity of house fly gut bacteria by both culture-dependent and culture-independent approaches. A total of 102 bacterial strains were isolated from the gut of 65 house flies collected from various public places including a garden, public park, garbage/dump area, public toilet, hospital, restaurant/canteen, mutton shop/market, and house/human habitation. Molecular phylogenetic analyses placed these isolates into 22 different genera. The majority of bacteria identified were known potential pathogens of the genera Klebsiella, Aeromonas, Shigella, Morganella, Providencia, and Staphylococcus. Culture-independent methods involved the construction of a 16S rRNA gene clone library, and sequence analyses supported culture recovery results. However, additional bacterial taxa not determined via culture recovery were revealed using this methodology and included members of the classes Alphaproteobacteria, Deltaproteobacteria, and the phylum Bacteroidetes. Here, we show that the house fly gut is an environmental reservoir for a vast number of bacterial species, which may have impacts on vector potential and pathogen transmission.

A microsatellite marker linkage map of the housefly, Musca domestica: evidence for male recombination

We present the first molecular marker linkage map for Musca domestica containing 35 microsatellite plus six visible markers. We report the development of 33 new microsatellite markers of which 19 are included in the linkage map. Two hundred and thirty-six F2 individuals were genotyped from three crosses yielding a linkage map consisting of five linkage groups that represent the five autosomes of the housefly. The map covers a total of 229.6 cM with an average marker spacing of 4.4 cM spanning approximately 80.2% of the genome. We found up to 29% recombination in male houseflies in contrast to most previous studies. The linkage map will add to genetic studies of the housefly.

Tsetse flies: genetics, evolution, and role as vectors

Tsetse flies (Diptera: Glossinidae) are an ancient taxon of one genus, Glossina, and limited species diversity. All are exclusively haematophagous and confined to sub-Saharan Africa. The Glossina are the principal vectors of African trypanosomes Trypanosoma sp. (Kinetoplastida: Trypanosomatidae) and as such, are of great medical and economic importance. Clearly tsetse flies and trypanosomes are coadapted and evolutionary interactions between them are manifest. Numerous clonally reproducing strains of Trypanosoma sp. exist and their genetic diversities and spatial distributions are inadequately known. Here I review the breeding structures of the principle trypanosome vectors, G. morsitans s.l., G. pallidipes, G. palpalis s.l. and G. fuscipes fuscipes. All show highly structured populations among which there is surprisingly little detectable gene flow. Rather less is known of the breeding structure of T. brucei sensu lato vis à vis their vector tsetse flies but many genetically differentiated strains exist in nature. Genetic recombination in Trypanosoma via meiosis has recently been demonstrated in the laboratory thereby furnishing a mechanism of strain differentiation in addition to that of simple mutation. Spatially and genetically representative sampling of both trypanosome species and strains and their Glossina vectors is a major barrier to a comprehensive understanding of their mutual relationships.

Complex population history of two Anopheles dirus mosquito species in Southeast Asia suggests the influence of Pleistocene climate change rather than human-mediated effects

Anopheles dirus and Anopheles baimaii are closely related species which feed on primates, particularly humans, and transmit malaria in the tropical forests of mainland Southeast Asia. Here, we report an in-depth phylogeographic picture based on 269 individuals from 21 populations from mainland Southeast Asia. Analysis of 1537 bp of mtDNA sequence revealed that the population history of A. baimaii is far more complex than previously thought. An old expansion (pre-300 kyr BP) was inferred in northern India/Bangladesh with a wave of south-eastwards expansion arriving at the Thai border (ca 135-173 kyr BP) followed by leptokurtic dispersal very recently (ca 16 kyr BP) into peninsular Thailand. The long and complex population history of these anthropophilic species suggests their expansions are not in response to the relatively recent (ca 40 kyr BP) human expansions in mainland Southeast Asia but, rather, fit well with our understanding of Pleistocene climatic change there.

Microsatellite evolution: Mutations, sequence variation, and homoplasy in the hypervariable avian microsatellite locus HrU10

Background

Microsatellites are frequently used genetic markers in a wide range of applications, primarily due to their high length polymorphism levels that can easily be genotyped by fragment length analysis. However, the mode of microsatellite evolution is yet not fully understood, and the role of interrupting motifs for the stability of microsatellites remains to be explored in more detail. Here we present a sequence analysis of mutation events and a description of the structure of repeated regions in the hypervariable, pentanucleotide microsatellite locus HrU10 in barn swallows (Hirundo rustica) and tree swallows (Tachycineta bicolor).

Results

In a large-scale parentage analysis in barn swallows and tree swallows, broods were screened for mutations at the HrU10 locus. In 41 cases in the barn swallows and 15 cases in the tree swallows, mutations corresponding to the loss or gain of one or two repeat units were detected. The parent and mutant offspring alleles were sequenced for 33 of these instances (26 in barn swallows and 7 in tree swallows). Replication slippage was considered the most likely mutational process. We tested the hypothesis that HrU10, a microsatellite with a wide allele size range, has an increased probability of introductions of interruptive motifs (IMs) with increasing length of the repeated region. Indeed, the number and length of the IMs was strongly positively correlated with the total length of the microsatellite. However, there was no significant correlation with the length of the longest stretch of perfectly repeated units, indicating a threshold level for the maximum length of perfectly repeated pentanucleotide motifs in stable HrU10 alleles. The combination of sequence and pedigree data revealed that 15 barn swallow mutations (58%) produced alleles that were size homoplasic to other alleles in the data set.

Conclusion

Our results give further insights into the mode of microsatellite evolution, and support the assumption of increased slippage rate with increased microsatellite length and a stabilizing effect of interrupting motifs for microsatellite regions consisting of perfect repeats. In addition, the observed extent of size homoplasy may impose a general caution against using hypervariable microsatellites in genetic diversity measures when alleles are identified by fragment length analysis only.

Phylogeography of stable fly (Diptera: Muscidae) estimated by diversity at ribosomal 16S and cytochrome oxidase I mitochondrial genes

The blood-feeding cosmopolitan stable fly, Stomoxys calcitrans L. (Diptera: Muscidae), is thought to disperse rapidly and widely, and earlier studies of allozyme variation were consistent with high vagility in this species. The geographic origins of New World populations are unknown. Diversity at mitochondrial loci r16S and cytochrome oxidase I was examined in 277 stable flies from 11 countries, including five zoogeographical regions. Of 809 nucleotides, 174 were polymorphic and 133 were parsimony informative. Seventy-six haplotypes were found in frequencies consistent with the Wright-Fisher infinite allele model. None were shared among four or more zoogeographical regions. The null hypothesis of mutation neutrality was not rejected, thereby validating the observed distribution. Fifty-nine haplotypes were singular, eight were private and confined to the Old World, and three of 76 haplotypes were shared between the Old and New World. Only 19 haplotypes were found in the New World, 14 of which were singletons. Haplotype and nucleotide diversities were heterogeneous among countries and regions. The most diversity was observed in sub-Saharan Africa. Regional differentiation indices were C(RT) = 0.26 and N(RT) = 0.31, indicating populations were highly structured macrogeographically. Palearctic and New World flies were the least differentiated from each other. There were strong genetic similarities among populations in the Nearctic, Neotropical, and Palearctic regions, and it is most likely that New World populations were derived from the Palearctic after 1492 CE, in the colonial era.

The influence of altitude and topography on genetic structure in the long-toed salamander (Ambystoma macrodactulym)

A primary goal of molecular ecology is to understand the influence of abiotic factors on the spatial distribution of genetic variation. Features including altitudinal clines, topography and landscape characteristics affect the proportion of suitable habitat, influence dispersal patterns, and ultimately structure genetic differentiation among populations. We studied the effects of altitude and topography on genetic variation of long-toed salamanders (Ambystoma macrodactylum), a geographically widespread amphibian species throughout northwestern North America. We focused on 10 low altitude sites (< 1200 m) and 11 high-altitude sites in northwestern Montana and determined multilocus genotypes for 549 individuals using seven microsatellite loci. We tested four hypotheses: (1) gene flow is limited between high- and low-altitude sites; and, (2) gene flow is limited among high-altitude sites due to harsh habitat and extreme topographical relief between sites; (3) low-altitude sites exhibit higher among-site gene flow due to frequent flooding events and low altitudinal relief; and (4) there is a negative correlation between altitude and genetic variation. Overall F(ST) values were moderate (0.08611; P < 0.001). Pairwise F(ST) estimates between high and low populations and a population graphing method supported the hypothesis that low-altitude and high-altitude sites, taken together, are genetically differentiated from each other. Also as predicted, gene flow is more prominent among low-altitude sites than high-altitude sites; low-altitude sites had a significantly lower F(ST) (0.03995; P < 0.001) than high altitude sites (F(ST) = 0.10271; P < 0.001). Use of Bayesian analysis of population structure (BAPS) resulted in delineation of 10 genetic groups, two among low-altitude populations and eight among high-altitude populations. In addition, within high altitude populations, basin-level genetic structuring was apparent. A nonequilibrium algorithm for detecting current migration rates supported these population distinctions. Finally, we also found a significant negative correlation between genetic diversity and altitude. These results are consistent with the hypothesis that topography and altitudinal gradients shape the spatial distribution of genetic variation in a species with a broad geographical range and diverse life history. Our study sheds light on which key factors limit dispersal and ultimately species' distributions.