A novel conjugative transposon carrying an autonomously amplified plasmid.

Tetracyclines serve as broad-spectrum antibiotics to treat bacterial infections. The discovery of new tetracycline resistance genes has led to new questions about the underlying mechanisms of resistance, gene transfer, and their relevance to human health. We tracked changes in the abundance of a 55-kbp conjugative transposon (CTn214) carrying tetQ, a tetracycline resistance gene, within a Bacteroides fragilis metagenome-assembled genome derived from shotgun sequencing of microbial DNA extracted from the ileal pouch of a patient with ulcerative colitis. The mapping of metagenomic reads to CTn214 revealed the multi-copy nature of a 17,044-nt region containing tetQ in samples collected during inflammation and uninflamed visits. B. fragilis cultivars isolated from the same patient during periods of inflammation harbored CTn214 integrated into the chromosome or both a circular, multi-copy, extrachromosomal region of the CTn214 containing tetQ and the corresponding integrated form. The tetracycline-dependent mechanism for the transmission of CTn214 is nearly identical to a common conjugative transposon found in the genome of B. fragilis (CTnDOT), but the autonomously amplified nature of a circular 17,044-nt region of CTn214 that codes for tetQ and the integration of the same sequence in the linear chromosome within the same cell is a novel observation. Genome and transcriptome sequencing of B. fragilis cultivars grown under different concentrations of tetracycline and ciprofloxacin indicates that tetQ in strains containing the circular form remains actively expressed regardless of treatment, while the expression of tetQ in strains containing the linear form increases only in the presence of tetracycline.IMPORTANCEThe exchange of antibiotic production and resistance genes between microorganisms can lead to the emergence of new pathogens. In this study, short-read mapping of metagenomic samples taken over time from the illeal pouch of a patient with ulcerative colitis to a Bacteroides fragilis metagenome-assembled genome revealed two distinct genomic arrangements of a novel conjugative transposon, CTn214, that encodes tetracycline resistance. The autonomous amplification of a plasmid-like circular form from CTn214 that includes tetQ potentially provides consistent ribosome protection against tetracycline. This mode of antibiotic resistance offers a novel mechanism for understanding the emergence of pathobionts like B. fragilis and their persistence for extended periods of time in patients with inflammatory bowel disease.

Microbial chemolithoautotrophs are abundant in salt marsh sediment following long-term experimental nitrate enrichment.

Long-term anthropogenic nitrate (NO3-) enrichment is a serious threat to many coastal systems. Nitrate reduction coupled with the oxidation of reduced forms of sulfur is conducted by chemolithoautotrophic microbial populations in a process that decreases nitrogen (N) pollution. However, little is known about the diversity and distribution of microbes capable of carbon fixation within salt marsh sediment and how they respond to long-term NO3- loading. We used genome-resolved metagenomics to characterize the distribution, phylogenetic relationships, and adaptations important to microbial communities within NO3--enriched sediment. We found NO3- reducing sulfur oxidizers became dominant members of the microbial community throughout the top 25 cm of the sediment following long-term NO3- enrichment. We also found that most of the chemolithoautotrophic genomes recovered contained striking metabolic versatility, including the potential for complete denitrification and evidence of mixotrophy. Phylogenetic reconstruction indicated that similar carbon fixation strategies and metabolic versatility can be found in several phylogenetic groups, but the genomes recovered here represent novel organisms. Our results suggest that the role of chemolithoautotrophy within NO3--enriched salt marsh sediments may be quantitatively more important for retaining carbon and filtering NO3- than previously indicated and further inquiry is needed to explicitly measure their contribution to carbon turnover and removal of N pollution.

Sulphide addition favours respiratory ammonification (DNRA) over complete denitrification and alters the active microbial community in salt marsh sediments.

The balance between nitrate respiration pathways, denitrification and dissimilatory nitrate (NO3 - ) reduction to ammonium (DNRA), determines whether bioavailable nitrogen is removed as N2 gas or recycled as ammonium. Saltwater intrusion and organic matter enrichment may increase sulphate reduction leading to sulphide accumulation. We investigated the effects of sulphide on the partitioning of NO3 - between complete denitrification and DNRA and the microbial communities in salt marsh sediments. Complete denitrification significantly decreased with increasing sulphide, resulting in an increase in the contribution of DNRA to NO3 - respiration. Alternative fates of NO3 - became increasingly important at higher sulphide treatments, which could include N2 O production and/or transport into intracellular vacuoles. Higher 16S transcript diversity was observed in the high sulphide treatment, with clear shifts in composition. Generally, low and no sulphide, coupled with high NO3 - , favoured the activity of Campylobacterales, Oceanospirillales and Altermonadales, all of which include opportunistic denitrifiers. High ∑sulphide conditions promoted the activity of potential sulphide oxidizing nitrate reducers (Desulfobulbaceae, Acidiferrobacteraceae and Xanthomonadales) and sulphate reducers (Desulfomonadaceae, Desulfobacteraceae). Our study highlights the tight coupling between N and S cycling, and the implications of these dynamics on the fate of bioavailable N in coastal environments susceptible to intermittent saltwater inundation and organic matter enrichment.

The Toxicogenome of Hyalella azteca: A Model for Sediment Ecotoxicology and Evolutionary Toxicology.

Hyalella azteca is a cryptic species complex of epibenthic amphipods of interest to ecotoxicology and evolutionary biology. It is the primary crustacean used in North America for sediment toxicity testing and an emerging model for molecular ecotoxicology. To provide molecular resources for sediment quality assessments and evolutionary studies, we sequenced, assembled, and annotated the genome of the H. azteca U.S. Lab Strain. The genome quality and completeness is comparable with other ecotoxicological model species. Through targeted investigation and use of gene expression data sets of H. azteca exposed to pesticides, metals, and other emerging contaminants, we annotated and characterized the major gene families involved in sequestration, detoxification, oxidative stress, and toxicant response. Our results revealed gene loss related to light sensing, but a large expansion in chemoreceptors, likely underlying sensory shifts necessary in their low light habitats. Gene family expansions were also noted for cytochrome P450 genes, cuticle proteins, ion transporters, and include recent gene duplications in the metal sequestration protein, metallothionein. Mapping of differentially expressed transcripts to the genome significantly increased the ability to functionally annotate toxicant responsive genes. The H. azteca genome will greatly facilitate development of genomic tools for environmental assessments and promote an understanding of how evolution shapes toxicological pathways with implications for environmental and human health.

Patient-Specific Bacteroides Genome Variants in Pouchitis.

A 2-year longitudinal microbiome study of 22 patients who underwent colectomy with an ileal pouch anal anastomosis detected significant increases in distinct populations of Bacteroides during 9 of 11 patient visits that coincided with inflammation (pouchitis). Oligotyping and metagenomic short-read annotation identified Bacteroides populations that occurred in early samples, bloomed during inflammation, and reappeared after antibiotic treatment. Targeted cultivation of Bacteroides isolates from the same individual at multiple time points and from several patients detected subtle genomic changes, including the identification of rapidly evolving genomic elements that differentiate isogenic strains of Bacteroides fragilis from the mucosa versus lumen. Each patient harbored Bacteroides spp. that are closely related to commonly occurring clinical isolates, including Bacteroides ovatus, B. thetaiotaomicron, B. vulgatus, and B. fragilis, which contained unique loci in different patients for synthesis of capsular polysaccharides. The presence of unique Bacteroides capsular polysaccharide loci within different hosts and between the lumen and mucosa may represent adaptations to stimulate, suppress, and evade host-specific immune responses at different microsites of the ileal pouch. IMPORTANCE: This longitudinal study provides an opportunity to describe shifts in the microbiomes of individual patients who suffer from ulcerative colitis (UC) prior to and following inflammation. Pouchitis serves as a model for UC with a predictable incidence of disease onset and enables prospective longitudinal investigations of UC etiology prior to inflammation. Because of insufficient criteria for predicting which patients will develop UC or pouchitis, the interpretation of cross-sectional study designs suffers from lack of information about the microbiome structure and host gene expression patterns that directly correlate with the onset of disease. Our unique longitudinal study design allows each patient to serve as their own control, providing information about the state of the microbiome and host prior to and during the course of disease. Of significance to the broader community, this study identifies microbial strains that may have genetic elements that trigger the onset of disease in susceptible hosts.

Genome reconstructions indicate the partitioning of ecological functions inside a phytoplankton bloom in the Amundsen Sea, Antarctica.

Antarctica polynyas support intense phytoplankton blooms, impacting their environment by a substantial depletion of inorganic carbon and nutrients. These blooms are dominated by the colony-forming haptophyte Phaeocystis antarctica and they are accompanied by a distinct bacterial population. Yet, the ecological role these bacteria may play in P. antarctica blooms awaits elucidation of their functional gene pool and of the geochemical activities they support. Here, we report on a metagenome (~160 million reads) analysis of the microbial community associated with a P. antarctica bloom event in the Amundsen Sea polynya (West Antarctica). Genomes of the most abundant Bacteroidetes and Proteobacteria populations have been reconstructed and a network analysis indicates a strong functional partitioning of these bacterial taxa. Three of them (SAR92, and members of the Oceanospirillaceae and Cryomorphaceae) are found in close association with P. antarctica colonies. Distinct features of their carbohydrate, nitrogen, sulfur and iron metabolisms may serve to support mutualistic relationships with P. antarctica. The SAR92 genome indicates a specialization in the degradation of fatty acids and dimethylsulfoniopropionate (compounds released by P. antarctica) into dimethyl sulfide, an aerosol precursor. The Oceanospirillaceae genome carries genes that may enhance algal physiology (cobalamin synthesis). Finally, the Cryomorphaceae genome is enriched in genes that function in cell or colony invasion. A novel pico-eukaryote, Micromonas related genome (19.6 Mb, ~94% completion) was also recovered. It contains the gene for an anti-freeze protein, which is lacking in Micromonas at lower latitudes. These draft genomes are representative for abundant microbial taxa across the Southern Ocean surface.

Anvi'o: an advanced analysis and visualization platform for 'omics data.

Advances in high-throughput sequencing and 'omics technologies are revolutionizing studies of naturally occurring microbial communities. Comprehensive investigations of microbial lifestyles require the ability to interactively organize and visualize genetic information and to incorporate subtle differences that enable greater resolution of complex data. Here we introduce anvi'o, an advanced analysis and visualization platform that offers automated and human-guided characterization of microbial genomes in metagenomic assemblies, with interactive interfaces that can link 'omics data from multiple sources into a single, intuitive display. Its extensible visualization approach distills multiple dimensions of information about each contig, offering a dynamic and unified work environment for data exploration, manipulation, and reporting. Using anvi'o, we re-analyzed publicly available datasets and explored temporal genomic changes within naturally occurring microbial populations through de novo characterization of single nucleotide variations, and linked cultivar and single-cell genomes with metagenomic and metatranscriptomic data. Anvi'o is an open-source platform that empowers researchers without extensive bioinformatics skills to perform and communicate in-depth analyses on large 'omics datasets.

Comparison of Sewage and Animal Fecal Microbiomes by Using Oligotyping Reveals Potential Human Fecal Indicators in Multiple Taxonomic Groups.

Most DNA-based microbial source tracking (MST) approaches target host-associated organisms within the order Bacteroidales, but the gut microbiota of humans and other animals contain organisms from an array of other taxonomic groups that might provide indicators of fecal pollution sources. To discern between human and nonhuman fecal sources, we compared the V6 regions of the 16S rRNA genes detected in fecal samples from six animal hosts to those found in sewage (as a proxy for humans). We focused on 10 abundant genera and used oligotyping, which can detect subtle differences between rRNA gene sequences from ecologically distinct organisms. Our analysis showed clear patterns of differential oligotype distributions between sewage and animal samples. Over 100 oligotypes of human origin occurred preferentially in sewage samples, and 99 human oligotypes were sewage specific. Sequences represented by the sewage-specific oligotypes can be used individually for development of PCR-based assays or together with the oligotypes preferentially associated with sewage to implement a signature-based approach. Analysis of sewage from Spain and Brazil showed that the sewage-specific oligotypes identified in U.S. sewage have the potential to be used as global alternative indicators of human fecal pollution. Environmental samples with evidence of prior human fecal contamination had consistent ratios of sewage signature oligotypes that corresponded to the trends observed for sewage. Our methodology represents a promising approach to identifying new bacterial taxa for MST applications and further highlights the potential of the family Lachnospiraceae to provide human-specific markers. In addition to source tracking applications, the patterns of the fine-scale population structure within fecal taxa suggest a fundamental relationship between bacteria and their hosts.

Sewage reflects the microbiomes of human populations.

UNLABELLED: Molecular characterizations of the gut microbiome from individual human stool samples have identified community patterns that correlate with age, disease, diet, and other human characteristics, but resources for marker gene studies that consider microbiome trends among human populations scale with the number of individuals sampled from each population. As an alternative strategy for sampling populations, we examined whether sewage accurately reflects the microbial community of a mixture of stool samples. We used oligotyping of high-throughput 16S rRNA gene sequence data to compare the bacterial distribution in a stool data set to a sewage influent data set from 71 U.S. cities. On average, only 15% of sewage sample sequence reads were attributed to human fecal origin, but sewage recaptured most (97%) human fecal oligotypes. The most common oligotypes in stool matched the most common and abundant in sewage. After informatically separating sequences of human fecal origin, sewage samples exhibited ~3× greater diversity than stool samples. Comparisons among municipal sewage communities revealed the ubiquitous and abundant occurrence of 27 human fecal oligotypes, representing an apparent core set of organisms in U.S. populations. The fecal community variability among U.S. populations was significantly lower than among individuals. It clustered into three primary community structures distinguished by oligotypes from either: Bacteroidaceae, Prevotellaceae, or Lachnospiraceae/Ruminococcaceae. These distribution patterns reflected human population variation and predicted whether samples represented lean or obese populations with 81 to 89% accuracy. Our findings demonstrate that sewage represents the fecal microbial community of human populations and captures population-level traits of the human microbiome. IMPORTANCE: The gut microbiota serves important functions in healthy humans. Numerous projects aim to define a healthy gut microbiome and its association with health states. However, financial considerations and privacy concerns limit the number of individuals who can be screened. By analyzing sewage from 71 cities, we demonstrate that geographically distributed U.S. populations share a small set of bacteria whose members represent various common community states within U.S. adults. Cities were differentiated by their sewage bacterial communities, and the community structures were good predictors of a city's estimated level of obesity. Our approach demonstrates the use of sewage as a means to sample the fecal microbiota from millions of people and its potential to elucidate microbiome patterns associated with human demographics.

A single genus in the gut microbiome reflects host preference and specificity.

Delineating differences in gut microbiomes of human and animal hosts contributes towards understanding human health and enables new strategies for detecting reservoirs of waterborne human pathogens. We focused upon Blautia, a single microbial genus that is important for nutrient assimilation as preliminary work suggested host-related patterns within members of this genus. In our dataset of 57 M sequence reads of the V6 region of the 16S ribosomal RNA gene in samples collected from seven host species, we identified 200 high-resolution taxonomic units within Blautia using oligotyping. Our analysis revealed 13 host-specific oligotypes that occurred exclusively in fecal samples of humans (three oligotypes), swine (six oligotypes), cows (one oligotype), deer (one oligotype), or chickens (two oligotypes). We identified an additional 171 oligotypes that exhibited differential abundance patterns among all the host species. Blautia oligotypes in the human population obtained from sewage and fecal samples displayed remarkable continuity. Oligotypes from only 10 Brazilian human fecal samples collected from individuals in a rural village encompassed 97% of all Blautia oligotypes found in a Brazilian sewage sample from a city of three million people. Further, 75% of the oligotypes in Brazilian human fecal samples matched those in US sewage samples, implying that a universal set of Blautia strains may be shared among culturally and geographically distinct human populations. Such strains can serve as universal markers to assess human fecal contamination in environmental samples. Our results indicate that host-specificity and host-preference patterns of organisms within this genus are driven by host physiology more than dietary habits.

Minimum entropy decomposition: unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences.

Molecular microbial ecology investigations often employ large marker gene datasets, for example, ribosomal RNAs, to represent the occurrence of single-cell genomes in microbial communities. Massively parallel DNA sequencing technologies enable extensive surveys of marker gene libraries that sometimes include nearly identical sequences. Computational approaches that rely on pairwise sequence alignments for similarity assessment and de novo clustering with de facto similarity thresholds to partition high-throughput sequencing datasets constrain fine-scale resolution descriptions of microbial communities. Minimum Entropy Decomposition (MED) provides a computationally efficient means to partition marker gene datasets into 'MED nodes', which represent homogeneous operational taxonomic units. By employing Shannon entropy, MED uses only the information-rich nucleotide positions across reads and iteratively partitions large datasets while omitting stochastic variation. When applied to analyses of microbiomes from two deep-sea cryptic sponges Hexadella dedritifera and Hexadella cf. dedritifera, MED resolved a key Gammaproteobacteria cluster into multiple MED nodes that are specific to different sponges, and revealed that these closely related sympatric sponge species maintain distinct microbial communities. MED analysis of a previously published human oral microbiome dataset also revealed that taxa separated by less than 1% sequence variation distributed to distinct niches in the oral cavity. The information theory-guided decomposition process behind the MED algorithm enables sensitive discrimination of closely related organisms in marker gene amplicon datasets without relying on extensive computational heuristics and user supervision.

Comparison of brush and biopsy sampling methods of the ileal pouch for assessment of mucosa-associated microbiota of human subjects.

BACKGROUND: Mucosal biopsy is the most common sampling technique used to assess microbial communities associated with the intestinal mucosa. Biopsies disrupt the epithelium and can be associated with complications such as bleeding. Biopsies sample a limited area of the mucosa, which can lead to potential sampling bias. In contrast to the mucosal biopsy, the mucosal brush technique is less invasive and provides greater mucosal coverage, and if it can provide equivalent microbial community data, it would be preferable to mucosal biopsies. RESULTS: We compared microbial samples collected from the intestinal mucosa using either a cytology brush or mucosal biopsy forceps. We collected paired samples from patients with ulcerative colitis (UC) who had previously undergone colectomy and ileal pouch anal anastomosis (IPAA), and profiled the microbial communities of the samples by sequencing V4-V6 or V4-V5 16S rRNA-encoding gene amplicons. Comparisons of 177 taxa in 16 brush-biopsy sample pairs had a mean R2 of 0.94. We found no taxa that varied significantly between the brush and biopsy samples after adjusting for multiple comparisons (false discovery rate ≤0.05). We also tested the reproducibility of DNA amplification and sequencing in 25 replicate pairs and found negligible variation (mean R2 = 0.99). A qPCR analysis of the two methods showed that the relative yields of bacterial DNA to human DNA were several-fold higher in the brush samples than in the biopsies. CONCLUSIONS: Mucosal brushing is preferred to mucosal biopsy for sampling the epithelial-associated microbiota. Although both techniques provide similar assessments of the microbial community composition, the brush sampling method has relatively more bacterial to host DNA, covers a larger surface area, and is less traumatic to the epithelium than the mucosal biopsy.

A filtering method to generate high quality short reads using illumina paired-end technology.

Consensus between independent reads improves the accuracy of genome and transcriptome analyses, however lack of consensus between very similar sequences in metagenomic studies can and often does represent natural variation of biological significance. The common use of machine-assigned quality scores on next generation platforms does not necessarily correlate with accuracy. Here, we describe using the overlap of paired-end, short sequence reads to identify error-prone reads in marker gene analyses and their contribution to spurious OTUs following clustering analysis using QIIME. Our approach can also reduce error in shotgun sequencing data generated from libraries with small, tightly constrained insert sizes. The open-source implementation of this algorithm in Python programming language with user instructions can be obtained from https://github.com/meren/illumina-utils.

Multiphasic analysis of the temporal development of the distal gut microbiota in patients following ileal pouch anal anastomosis.

BACKGROUND: The indigenous gut microbiota are thought to play a crucial role in the development and maintenance of the abnormal inflammatory responses that are the hallmark of inflammatory bowel disease. Direct tests of the role of the gut microbiome in these disorders are typically limited by the fact that sampling of the microbiota generally occurs once disease has become manifest. This limitation could potentially be circumvented by studying patients who undergo total proctocolectomy with ileal pouch anal anastomosis (IPAA) for the definitive treatment of ulcerative colitis. A subset of patients who undergo IPAA develops an inflammatory condition known as pouchitis, which is thought to mirror the pathogenesis of ulcerative colitis. Following the development of the microbiome of the pouch would allow characterization of the microbial community that predates the development of overt disease. RESULTS: We monitored the development of the pouch microbiota in four patients who underwent IPAA. Mucosal and luminal samples were obtained prior to takedown of the diverting ileostomy and compared to samples obtained 2, 4 and 8 weeks after intestinal continuity had been restored. Through the combined analysis of 16S rRNA-encoding gene amplicons, targeted 16S amplification and microbial cultivation, we observed major changes in structure and function of the pouch microbiota following ileostomy. There is a relative increase in anaerobic microorganisms with the capacity for fermentation of complex carbohydrates, which corresponds to the physical stasis of intestinal contents in the ileal pouch. Compared to the microbiome structure encountered in the colonic mucosa of healthy individuals, the pouch microbial community in three of the four individuals was quite distinct. In the fourth patient, a community that was much like that seen in a healthy colon was established, and this patient also had the most benign clinical course of the four patients, without the development of pouchitis 2 years after IPAA. CONCLUSIONS: The microbiota that inhabit the ileal-anal pouch of patients who undergo IPAA for treatment of ulcerative colitis demonstrate significant structural and functional changes related to the restoration of fecal flow. Our preliminary results suggest once the pouch has assumed the physiologic role previously played by the intact colon, the precise structure and function of the pouch microbiome, relative to a normal colonic microbiota, will determine if there is establishment of a stable, healthy mucosal environment or the reinitiation of the pathogenic cascade that results in intestinal inflammation.

Detection of broadly distributed sodium channel alleles characteristic of insect pyrethroid resistance in West Nile virus vector Culex pipiens complex mosquitoes in the United States.

West Nile virus (WNV) has emerged as a health threat to the North American population since its initial outbreak in New York City in 1999. Culex (Culex) pipiens complex mosquitoes have been considered to play the primary role in the enzootic maintenance and transmission of WNV in North America. The voltage-gated sodium channel (NaCh) gene contains pyrethroid resistance-associated mutations in the coding region in many insect species. However, the knowledge of potential NaCh mutations was minimal in Culex. Seeking pyrethroid resistance alleles in Culex, we evaluated a transect along the east coast of the United States with an NaCh-based genotyping tool that amplified a portion of the transcribed sequence containing kdr mutations and the intron immediately downstream of the mutation site. Three genotypes that are typically associated with pyrethroid resistance in insects have been identified in Culex pipiens complex mosquitoes in this study: susceptible wild type kds, the classical knock-down resistance Leu --> Phe mutation (Phe/kdr), and a second resistance mechanism, a Leu --> Ser mutation (Ser/kdr). Moreover, we observed heterozygotic individual mosquitoes possessing both kdr alleles. Results of this study advance our knowledge of the potential for pyrethroid insecticide resistance among the populations of Cx. pipiens complex in the United States.

Microsatellite data suggest significant population structure and differentiation within the malaria vector Anopheles darlingi in Central and South America.

BACKGROUND: Anopheles darlingi is the most important malaria vector in the Neotropics. An understanding of A. darlingi's population structure and contemporary gene flow patterns is necessary if vector populations are to be successfully controlled. We assessed population genetic structure and levels of differentiation based on 1,376 samples from 31 localities throughout the Peruvian and Brazilian Amazon and Central America using 5-8 microsatellite loci. RESULTS: We found high levels of polymorphism for all of the Amazonian populations (mean RS = 7.62, mean HO = 0.742), and low levels for the Belize and Guatemalan populations (mean RS = 4.3, mean HO = 0.457). The Bayesian clustering analysis revealed five population clusters: northeastern Amazonian Brazil, southeastern and central Amazonian Brazil, western and central Amazonian Brazil, Peruvian Amazon, and the Central American populations. Within Central America there was low non-significant differentiation, except for between the populations separated by the Maya Mountains. Within Amazonia there was a moderate level of significant differentiation attributed to isolation by distance. Within Peru there was no significant population structure and low differentiation, and some evidence of a population expansion. The pairwise estimates of genetic differentiation between Central America and Amazonian populations were all very high and highly significant (FST = 0.1859 - 0.3901, P < 0.05). Both the DA and FST distance-based trees illustrated the main division to be between Central America and Amazonia. CONCLUSION: We detected a large amount of population structure in Amazonia, with three population clusters within Brazil and one including the Peru populations. The considerable differences in Ne among the populations may have contributed to the observed genetic differentiation. All of the data suggest that the primary division within A. darlingi corresponds to two white gene genotypes between Amazonia (genotype 1) and Central America, parts of Colombia and Venezuela (genotype 2), and are in agreement with previously published mitochondrial COI gene sequences interpreted as incipient species. Overall, it appears that two main factors have contributed to the genetic differentiation between the population clusters: physical distance between the populations and the differences in effective population sizes among the subpopulations.

Population structure of the malaria vector Anopheles darlingi in a malaria-endemic region of eastern Amazonian Brazil.

Anopheles darlingi is the primary malaria vector in Latin America, and is especially important in Amazonian Brazil. Historically, control efforts have been focused on indoor house spraying using a variety of insecticides, but since the mid-1990s there has been a shift to patient treatment and focal insecticide fogging. Anopheles darlingi was believed to have been significantly reduced in a gold-mining community, Peixoto de Azevedo (in Mato Grosso State), in the early 1990s by insecticide use during a severe malaria epidemic. In contrast, although An. darlingi was eradicated from some districts of the city of Belem (the capital of Para State) in 1968 to reduce malaria, populations around the water protection area in the eastern district were treated only briefly. To investigate the population structure of An. darlingi including evidence for a population bottleneck in Peixoto, we analyzed eight microsatellite loci of 256 individuals from seven locations in Brazil: three in Amapa State, three in Para State, and one in Mato Grosso State. Allelic diversity and mean expected heterozygosity were high for all populations (mean number alleles/locus and H(E) were 13.5 and 0.834, respectively) and did not differ significantly between locations. Significant heterozygote deficits were associated with linkage disequilibrium, most likely due to either the Wahlund effect or selection. We found no evidence for a population bottleneck in Peixoto, possibly because the reduction was not extreme enough to be detected. Overall estimates of long-term N(e) varied from 92.4 individuals under the linkage disequilibrium model to infinity under the heterozygote excess model. Fixation indices and analysis of molecular variance demonstrated significant differentiation between locations north and south of the Amazon River, suggesting a degree of genetic isolation between them, attributed to isolation by distance.