A machine-learning approach to map landscape connectivity in Aedes aegypti with genetic and environmental data.

Mapping landscape connectivity is important for controlling invasive species and disease vectors. Current landscape genetics methods are often constrained by the subjectivity of creating resistance surfaces and the difficulty of working with interacting and correlated environmental variables. To overcome these constraints, we combine the advantages of a machine-learning framework and an iterative optimization process to develop a method for integrating genetic and environmental (e.g., climate, land cover, human infrastructure) data. We validate and demonstrate this method for the Aedes aegypti mosquito, an invasive species and the primary vector of dengue, yellow fever, chikungunya, and Zika. We test two contrasting metrics to approximate genetic distance and find Cavalli-Sforza-Edwards distance (CSE) performs better than linearized FST The correlation (R) between the model's predicted genetic distance and actual distance is 0.83. We produce a map of genetic connectivity for Ae. aegypti's range in North America and discuss which environmental and anthropogenic variables are most important for predicting gene flow, especially in the context of vector control.

Infection with endosymbiotic Spiroplasma disrupts tsetse (Glossina fuscipes fuscipes) metabolic and reproductive homeostasis.

Tsetse flies (Glossina spp.) house a population-dependent assortment of microorganisms that can include pathogenic African trypanosomes and maternally transmitted endosymbiotic bacteria, the latter of which mediate numerous aspects of their host's metabolic, reproductive, and immune physiologies. One of these endosymbionts, Spiroplasma, was recently discovered to reside within multiple tissues of field captured and laboratory colonized tsetse flies grouped in the Palpalis subgenera. In various arthropods, Spiroplasma induces reproductive abnormalities and pathogen protective phenotypes. In tsetse, Spiroplasma infections also induce a protective phenotype by enhancing the fly's resistance to infection with trypanosomes. However, the potential impact of Spiroplasma on tsetse's viviparous reproductive physiology remains unknown. Herein we employed high-throughput RNA sequencing and laboratory-based functional assays to better characterize the association between Spiroplasma and the metabolic and reproductive physiologies of G. fuscipes fuscipes (Gff), a prominent vector of human disease. Using field-captured Gff, we discovered that Spiroplasma infection induces changes of sex-biased gene expression in reproductive tissues that may be critical for tsetse's reproductive fitness. Using a Gff lab line composed of individuals heterogeneously infected with Spiroplasma, we observed that the bacterium and tsetse host compete for finite nutrients, which negatively impact female fecundity by increasing the length of intrauterine larval development. Additionally, we found that when males are infected with Spiroplasma, the motility of their sperm is compromised following transfer to the female spermatheca. As such, Spiroplasma infections appear to adversely impact male reproductive fitness by decreasing the competitiveness of their sperm. Finally, we determined that the bacterium is maternally transmitted to intrauterine larva at a high frequency, while paternal transmission was also noted in a small number of matings. Taken together, our findings indicate that Spiroplasma exerts a negative impact on tsetse fecundity, an outcome that could be exploited for reducing tsetse population size and thus disease transmission.

Phylogenomic analysis of Syngnathidae reveals novel relationships, origins of endemic diversity and variable diversification rates.

BACKGROUND: Seahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved. RESULTS: Here, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea. CONCLUSIONS: Our results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events.

Big Data in Conservation Genomics: Boosting Skills, Hedging Bets, and Staying Current in the Field.

A current challenge in the fields of evolutionary, ecological, and conservation genomics is balancing production of large-scale datasets with additional training often required to handle such datasets. Thus, there is an increasing need for conservation geneticists to continually learn and train to stay up-to-date through avenues such as symposia, meetings, and workshops. The ConGen meeting is a near-annual workshop that strives to guide participants in understanding population genetics principles, study design, data processing, analysis, interpretation, and applications to real-world conservation issues. Each year of ConGen gathers a diverse set of instructors, students, and resulting lectures, hands-on sessions, and discussions. Here, we summarize key lessons learned from the 2019 meeting and more recent updates to the field with a focus on big data in conservation genomics. First, we highlight classical and contemporary issues in study design that are especially relevant to working with big datasets, including the intricacies of data filtering. We next emphasize the importance of building analytical skills and simulating data, and how these skills have applications within and outside of conservation genetics careers. We also highlight recent technological advances and novel applications to conservation of wild populations. Finally, we provide data and recommendations to support ongoing efforts by ConGen organizers and instructors-and beyond-to increase participation of underrepresented minorities in conservation and eco-evolutionary sciences. The future success of conservation genetics requires both continual training in handling big data and a diverse group of people and approaches to tackle key issues, including the global biodiversity-loss crisis.

The population genomics of multiple tsetse fly (Glossina fuscipes fuscipes) admixture zones in Uganda.

Understanding the mechanisms that enforce, maintain or reverse the process of speciation is an important challenge in evolutionary biology. This study investigates the patterns of divergence and discusses the processes that form and maintain divergent lineages of the tsetse fly Glossina fuscipes fuscipes in Uganda. We sampled 251 flies from 18 sites spanning known genetic lineages and the four admixture zones between them. We apply population genomics, hybrid zone and approximate Bayesian computation to the analysis of three types of genetic markers: 55,267 double-digest restriction site-associated DNA (ddRAD) SNPs to assess genome-wide admixture, 16 microsatellites to provide continuity with published data and accurate biogeographic modelling, and a 491-bp fragment of mitochondrial cytochrome oxidase I and II to infer maternal inheritance patterns. Admixture zones correspond with regions impacted by the reorganization of Uganda's river networks that occurred during the formation of the West African Rift system over the last several hundred thousand years. Because tsetse fly population distributions are defined by rivers, admixture zones likely represent both old and new regions of secondary contact. Our results indicate that older hybrid zones contain mostly parental types, while younger zones contain variable hybrid types resulting from multiple generations of interbreeding. These findings suggest that reproductive barriers are nearly complete in the older admixture zones, while nearly absent in the younger admixture zones. Findings are consistent with predictions of hybrid zone theory: Populations in zones of secondary contact transition rapidly from early to late stages of speciation or collapse all together.

Ectomycorrhizas and tree seedling establishment are strongly influenced by forest edge proximity but not soil inoculum.

Reforestation is challenging when timber harvested areas have been degraded, invaded by nonnative species, or are of marginal suitability to begin with. Conifers form mutualistic partnerships with ectomycorrhizal fungi (EMF) to obtain greater access to soil resources, and these partnerships may be especially important in degraded areas. However, timber harvest can impact mycorrhizal fungi by removing or compacting topsoil, removing host plants, and warming and drying the soil. We used a field experiment to evaluate the role of EMF in Douglas-fir reforestation in clearcuts invaded by Cytisus scoparius (Scotch broom) where traditional reforestation approaches have repeatedly failed. We tested how planting distance from intact Douglas-fir forest edges influenced reforestation success and whether inoculation with forest soils can be used to restore EMF relationships. We used an Illumina DNA sequencing approach to measure the abundance, richness and composition of ectomycorrhizal fungi on Douglas-fir roots, and assessed differences in Douglas-fir seedling survival and growth near to and far from forest edges with and without forest soil inoculum. Planting Douglas-fir seedlings near forest edges increased seedling survival, growth, and EMF root colonization. Edge proximity had no effect on EMF richness but did change fungal community composition. Inoculations with forest soil did not increase EMF abundance or richness or change community composition, nor did it improve seedling establishment. With Illumina sequencing, we identified two to three times greater species richness than described in previous edge effects studies. Of the 95 EMF species we identified, 40% of the species occurred on less than 5% of the seedlings. The ability to detect fungi at low abundance may explain why we did not detect differences in EMF richness with distance to hosts as previous studies. Our findings suggest that forest edges are suitable for reforestation, even when the interiors of deforested areas are not. We advocate for timber harvest designs that maximize edge habitat where ectomycorrhizal fungi contribute to tree establishment. However, this study does not support the use of inoculation with forest soil as a simple method to enhance EMF and seedling survival.

Molecular phylogeny and patterns of diversification in syngnathid fishes.

The family Syngnathidae is a large and diverse clade of morphologically unique bony fishes, with 57 genera and 300 described species of seahorses, pipefishes, pipehorses, and seadragons. They primarily inhabit shallow coastal waters in temperate and tropical oceans, and are characterized by a fused jaw, male brooding, and extraordinary crypsis. Phylogenetic relationships within the Syngnathidae remain poorly resolved due to lack of generic taxon sampling, few diagnostic morphological characters, and limited molecular data. The phylogenetic placement of the threatened, commercially exploited seahorses remains a topic of intense interest, with conflicting topologies based on morphology and predominantly mitochondrial genetic data. In this study, we integrate eight nuclear and mitochondrial markers and 17 morphological characters to investigate the phylogenetic structure of the family Syngnathidae at the generic level. We include 91 syngnathid species representing 48 of the 57 recognized genera, all major ocean basins, and a broad array of temperate and tropical habitats including rocky and coral reefs, sand and silt, mangroves, seagrass beds, estuaries, and rivers. Maximum likelihood and Bayesian analyses of 5160bp from eight loci produced high congruence among alternate topologies, defining well-supported and sometimes novel clades. We present a hypothesis that confirms a deep phylogenetic split between lineages with trunk- or tail-brood pouch placement, and provides significant new insights into the morphological evolution and biogeography of this highly derived fish clade. Based on the fundamental division between lineages - the tail brooding "Urophori" and the trunk brooding "Gastrophori" - we propose a revision of Syngnathidae classification into only two subfamilies: the Nerophinae and the Syngnathinae. We find support for distinct principal clades within the trunk-brooders and tail-brooders, the latter of which include seahorses, seadragons, independent lineages of pipehorses, and clades that originated in southern Australia and the Western Atlantic. We suggest the seahorse genus Hippocampus is of Indo-Pacific origin and its sister clade is an unexpected grouping of several morphologically disparate Indo-Pacific genera, including the Pacific pygmy pipehorses. Taxonomic revision is required for multiple genera, particularly to reflect deep evolutionary splits in nominal lineages from the Atlantic versus the Indo-Pacific.

Introgression between invasive and native blue mussels (genus Mytilus) in the central California hybrid zone.

The ecological and genetic factors determining the extent of introgression between species in secondary contact zones remain poorly understood. Here, we investigate the relative importance of isolating barriers and the demographic expansion of invasive Mytilus galloprovincialis on the magnitude and the direction of introgression with the native Mytilus trossulus in a hybrid zone in central California. We use double-digest restriction-site-associated DNA sequencing (ddRADseq) to genotype 1337 randomly selected single nucleotide polymorphisms and accurately distinguish early and advanced generation hybrids for the first time in the central California Mytilus spp. hybrid zone. Weak levels of introgression were observed in both directions but were slightly more prevalent from the native M. trossulus into the invasive M. galloprovincialis. Few early and advanced backcrossed individuals were observed across the hybrid zone confirming the presence of strong barriers to interbreeding. Heterogeneous patterns of admixture across the zone of contact were consistent with the colonization history of M. galloprovincialis with more extensive introgression in northern localities furthest away from the putative site of introduction in southern California. These observations reinforce the importance of dynamic spatial and demographic expansions in determining patterns of introgression between close congeners, even in those with high dispersal potential and well-developed reproductive barriers. Our results suggest that the threat posed by invasive M. galloprovincialis is more ecological than genetic as it has displaced, and continues to displace the native M. trossulus from much of central and southern California.

Genomics and conservation: Guidance from training to analyses and applications.

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Genomic evidence of sex chromosome aneuploidy and infection-associated genotypes in the tsetse fly Glossina fuscipes, the major vector of African trypanosomiasis in Uganda.

The primary vector of the trypanosome parasite causing human and animal African trypanosomiasis in Uganda is the riverine tsetse fly Glossina fuscipes fuscipes (Gff). Our study improved the Gff genome assembly with whole genome 10× Chromium sequencing of a lab reared pupae, identified autosomal versus sex-chromosomal regions of the genome with ddRAD-seq data from 627 field caught Gff, and identified SNPs associated with trypanosome infection with genome-wide association (GWA) analysis in a subset of 351 flies. Results from 10× Chromium sequencing greatly improved Gff genome assembly metrics and assigned a full third of the genome to the sex chromosome. Results from ddRAD-seq suggested possible sex-chromosome aneuploidy in Gff and identified a single autosomal SNP to be highly associated with trypanosome infection. The top associated SNP was ∼1100 bp upstream of the gene lecithin cholesterol acyltransferase (LCAT), an important component of the molecular pathway that initiates trypanosome lysis and protection in mammals. Results suggest that there may be naturally occurring genetic variation in Gff in genomic regions in linkage disequilibrium with LCAT that can protect against trypanosome infection, thereby paving the way for targeted research into novel vector control strategies that can promote parasite resistance in natural populations.

Multiple introductions and overwintering shape the progressive invasion of Aedes albopictus beyond the Alps.

Aedes albopictus originates from Southeast Asia and is considered one of the most invasive species globally. This mosquito is a nuisance and a disease vector of significant public health relevance. In Europe, Ae. albopictus is firmly established and widespread south of the Alps, a mountain range that forms a formidable biogeographic barrier to many organisms. Recent reports of Ae. albopictus north of the Alps raise questions of (1) the origins of its recent invasion, and (2) if this mosquito has established overwintering populations north of the Alps. To answer these questions, we analyzed population genomic data from >4000 genome-wide SNPs obtained through double-digest restriction site-associated DNA sequencing. We collected SNP data from specimens from six sites in Switzerland, north and south of the Alps, and analyzed them together with specimens from other 33 European sites, five from the Americas, and five from its Asian native range. At a global level, we detected four genetic clusters with specimens from Indonesia, Brazil, and Japan as the most differentiated, whereas specimens from Europe, Hong Kong, and USA largely overlapped. Across the Alps, we detected a weak genetic structure and high levels of genetic admixture, supporting a scenario of rapid and human-aided dispersal along transportation routes. While the genetic pattern suggests frequent re-introductions into Switzerland from Italian sources, the recovery of a pair of full siblings in two consecutive years in Strasbourg, France, suggests the presence of an overwintering population north of the Alps. The suggestion of overwintering populations of Ae. albopictus north of the Alps and the expansion patterns identified points to an increased risk of further northward expansion and the need for increased surveillance of mosquito populations in Northern Europe.

A machine learning approach to integrating genetic and ecological data in tsetse flies (Glossina pallidipes) for spatially explicit vector control planning.

Vector control is an effective strategy for reducing vector-borne disease transmission, but requires knowledge of vector habitat use and dispersal patterns. Our goal was to improve this knowledge for the tsetse species Glossina pallidipes, a vector of human and animal African trypanosomiasis, which are diseases that pose serious health and socioeconomic burdens across sub-Saharan Africa. We used random forest regression to (i) build and integrate models of G. pallidipes habitat suitability and genetic connectivity across Kenya and northern Tanzania and (ii) provide novel vector control recommendations. Inputs for the models included field survey records from 349 trap locations, genetic data from 11 microsatellite loci from 659 flies and 29 sampling sites, and remotely sensed environmental data. The suitability and connectivity models explained approximately 80% and 67% of the variance in the occurrence and genetic data and exhibited high accuracy based on cross-validation. The bivariate map showed that suitability and connectivity vary independently across the landscape and was used to inform our vector control recommendations. Post hoc analyses show spatial variation in the correlations between the most important environmental predictors from our models and each response variable (e.g., suitability and connectivity) as well as heterogeneity in expected future climatic change of these predictors. The bivariate map suggests that vector control is most likely to be successful in the Lake Victoria Basin and supports the previous recommendation that G. pallidipes from most of eastern Kenya should be managed as a single unit. We further recommend that future monitoring efforts should focus on tracking potential changes in vector presence and dispersal around the Serengeti and the Lake Victoria Basin based on projected local climatic shifts. The strong performance of the spatial models suggests potential for our integrative methodology to be used to understand future impacts of climate change in this and other vector systems.

Phylogeography and population structure of the tsetse fly Glossina pallidipes in Kenya and the Serengeti ecosystem.

Glossina pallidipes is the main vector of animal African trypanosomiasis and a potential vector of human African trypanosomiasis in eastern Africa where it poses a large economic burden and public health threat. Vector control efforts have succeeded in reducing infection rates, but recent resurgence in tsetse fly population density raises concerns that vector control programs require improved strategic planning over larger geographic and temporal scales. Detailed knowledge of population structure and dispersal patterns can provide the required information to improve planning. To this end, we investigated the phylogeography and population structure of G. pallidipes over a large spatial scale in Kenya and northern Tanzania using 11 microsatellite loci genotyped in 600 individuals. Our results indicate distinct genetic clusters east and west of the Great Rift Valley, and less distinct clustering of the northwest separate from the southwest (Serengeti ecosystem). Estimates of genetic differentiation and first-generation migration indicated high genetic connectivity within genetic clusters even across large geographic distances of more than 300 km in the east, but only occasional migration among clusters. Patterns of connectivity suggest isolation by distance across genetic breaks but not within genetic clusters, and imply a major role for river basins in facilitating gene flow in G. pallidipes. Effective population size (Ne) estimates and results from Approximate Bayesian Computation further support that there has been recent G. pallidipes population size fluctuations in the Serengeti ecosystem and the northwest during the last century, but also suggest that the full extent of differences in genetic diversity and population dynamics between the east and the west was established over evolutionary time periods (tentatively on the order of millions of years). Findings provide further support that the Serengeti ecosystem and northwestern Kenya represent independent tsetse populations. Additionally, we present evidence that three previously recognized populations (the Mbeere-Meru, Central Kenya and Coastal "fly belts") act as a single population and should be considered as a single unit in vector control.

Spatio-temporal distribution of Spiroplasma infections in the tsetse fly (Glossina fuscipes fuscipes) in northern Uganda.

Tsetse flies (Glossina spp.) are vectors of parasitic trypanosomes, which cause human (HAT) and animal African trypanosomiasis (AAT) in sub-Saharan Africa. In Uganda, Glossina fuscipes fuscipes (Gff) is the main vector of HAT, where it transmits Gambiense disease in the northwest and Rhodesiense disease in central, southeast and western regions. Endosymbionts can influence transmission efficiency of parasites through their insect vectors via conferring a protective effect against the parasite. It is known that the bacterium Spiroplasma is capable of protecting its Drosophila host from infection with a parasitic nematode. This endosymbiont can also impact its host's population structure via altering host reproductive traits. Here, we used field collections across 26 different Gff sampling sites in northern and western Uganda to investigate the association of Spiroplasma with geographic origin, seasonal conditions, Gff genetic background and sex, and trypanosome infection status. We also investigated the influence of Spiroplasma on Gff vector competence to trypanosome infections under laboratory conditions. Generalized linear models (GLM) showed that Spiroplasma probability was correlated with the geographic origin of Gff host and with the season of collection, with higher prevalence found in flies within the Albert Nile (0.42 vs 0.16) and Achwa River (0.36 vs 0.08) watersheds and with higher prevalence detected in flies collected in the intermediate than wet season. In contrast, there was no significant correlation of Spiroplasma prevalence with Gff host genetic background or sex once geographic origin was accounted for in generalized linear models. Additionally, we found a potential negative correlation of Spiroplasma with trypanosome infection, with only 2% of Spiroplasma infected flies harboring trypanosome co-infections. We also found that in a laboratory line of Gff, parasitic trypanosomes are less likely to colonize the midgut in individuals that harbor Spiroplasma infection. These results indicate that Spiroplasma infections in tsetse may be maintained by not only maternal but also via horizontal transmission routes, and Spiroplasma infections may also have important effects on trypanosome transmission efficiency of the host tsetse. Potential functional effects of Spiroplasma infection in Gff could have impacts on vector control approaches to reduce trypanosome infections.

Genetic Differentiation of Glossina pallidipes Tsetse Flies in Southern Kenya.

The tsetse fly Glossina pallidipes, the major vector of the parasite that causes animal African trypanosomiasis in Kenya, has been subject to intense control measures with only limited success. The G. pallidipes population dynamics and dispersal patterns that underlie limited success in vector control campaigns remain unresolved, and knowledge on genetic connectivity can provide insights, and thereby improve control and monitoring efforts. We therefore investigated the population structure and estimated migration and demographic parameters in G. pallidipes using genotypic data from 11 microsatellite loci scored in 250 tsetse flies collected from eight localities in Kenya. Clustering analysis identified two genetically distinct eastern and western clusters (mean between-cluster F ST = 0.202) separated by the Great Rift Valley. We also found evidence of admixture and migration between the eastern and western clusters, isolation by distance, and a widespread signal of inbreeding. We detected differences in population dynamics and dispersal patterns between the western and eastern clusters. These included lower genetic diversity (allelic richness; 7.48 versus 10.99), higher relatedness (percent related individuals; 21.4% versus 9.1%), and greater genetic differentiation (mean within-cluster F ST; 0.183 versus 0.018) in the western than the eastern cluster. Findings are consistent with the presence of smaller, less well-connected populations in Western relative to eastern Kenya. These data suggest that recent anthropogenic influences such as land use changes and vector control programs have influenced population dynamics in G. pallidipes in Kenya, and that vector control efforts should include some region-specific strategies to effectively control this disease vector.

A spatial genetics approach to inform vector control of tsetse flies (Glossina fuscipes fuscipes) in Northern Uganda.

Tsetse flies (genus Glossina) are the only vector for the parasitic trypanosomes responsible for sleeping sickness and nagana across sub-Saharan Africa. In Uganda, the tsetse fly Glossina fuscipes fuscipes is responsible for transmission of the parasite in 90% of sleeping sickness cases, and co-occurrence of both forms of human-infective trypanosomes makes vector control a priority. We use population genetic data from 38 samples from northern Uganda in a novel methodological pipeline that integrates genetic data, remotely sensed environmental data, and hundreds of field-survey observations. This methodological pipeline identifies isolated habitat by first identifying environmental parameters correlated with genetic differentiation, second, predicting spatial connectivity using field-survey observations and the most predictive environmental parameter(s), and third, overlaying the connectivity surface onto a habitat suitability map. Results from this pipeline indicated that net photosynthesis was the strongest predictor of genetic differentiation in G. f. fuscipes in northern Uganda. The resulting connectivity surface identified a large area of well-connected habitat in northwestern Uganda, and twenty-four isolated patches on the northeastern margin of the G. f. fuscipes distribution. We tested this novel methodological pipeline by completing an ad hoc sample and genetic screen of G. f. fuscipes samples from a model-predicted isolated patch, and evaluated whether the ad hoc sample was in fact as genetically isolated as predicted. Results indicated that genetic isolation of the ad hoc sample was as genetically isolated as predicted, with differentiation well above estimates made in samples from within well-connected habitat separated by similar geographic distances. This work has important practical implications for the control of tsetse and other disease vectors, because it provides a way to identify isolated populations where it will be safer and easier to implement vector control and that should be prioritized as study sites during the development and improvement of vector control methods.

Sequence-Based Analysis of Thermal Adaptation and Protein Energy Landscapes in an Invasive Blue Mussel (Mytilus galloprovincialis).

Adaptive responses to thermal stress in poikilotherms plays an important role in determining competitive ability and species distributions. Amino acid substitutions that affect protein stability and modify the thermal optima of orthologous proteins may be particularly important in this context. Here, we examine a set of 2,770 protein-coding genes to determine if proteins in a highly invasive heat tolerant blue mussel (Mytilus galloprovincialis) contain signals of adaptive increases in protein stability relative to orthologs in a more cold tolerant M. trossulus. Such thermal adaptations might help to explain, mechanistically, the success with which the invasive marine mussel M. galloprovincialis has displaced native species in contact zones in the eastern (California) and western (Japan) Pacific. We tested for stabilizing amino acid substitutions in warm tolerant M. galloprovincialis relative to cold tolerant M. trossulus with a generalized linear model that compares in silico estimates of recent changes in protein stability among closely related congeners. Fixed substitutions in M. galloprovincialis were 3,180.0 calories per mol per substitution more stabilizing at genes with both elevated dN/dS ratios and transcriptional responses to heat stress, and 705.8 calories per mol per substitution more stabilizing across all 2,770 loci investigated. Amino acid substitutions concentrated in a small number of genes were more stabilizing in M. galloprovincialis compared with cold tolerant M. trossulus. We also tested for, but did not find, enrichment of a priori GO terms in genes with elevated dN/dS ratios in M. galloprovincialis. This might indicate that selection for thermodynamic stability is generic across all lineages, and suggests that the high change in estimated protein stability that we observed in M. galloprovincialis is driven by selection for extra stabilizing substitutions, rather than by higher incidence of selection in a greater number of genes in this lineage. Nonetheless, our finding of more stabilizing amino acid changes in the warm adapted lineage is important because it suggests that adaption for thermal stability has contributed to M. galloprovincialis' superior tolerance to heat stress, and that pairing tests for positive selection and tests for transcriptional response to heat stress can identify candidates of protein stability adaptation.

Multiple evolutionary origins of Trypanosoma evansi in Kenya.

Trypanosoma evansi is the parasite causing surra, a form of trypanosomiasis in camels and other livestock, and a serious economic burden in Kenya and many other parts of the world. Trypanosoma evansi transmission can be sustained mechanically by tabanid and Stomoxys biting flies, whereas the closely related African trypanosomes T. brucei brucei and T. b. rhodesiense require cyclical development in tsetse flies (genus Glossina) for transmission. In this study, we investigated the evolutionary origins of T. evansi. We used 15 polymorphic microsatellites to quantify levels and patterns of genetic diversity among 41 T. evansi isolates and 66 isolates of T. b. brucei (n = 51) and T. b. rhodesiense (n = 15), including many from Kenya, a region where T. evansi may have evolved from T. brucei. We found that T. evansi strains belong to at least two distinct T. brucei genetic units and contain genetic diversity that is similar to that in T. brucei strains. Results indicated that the 41 T. evansi isolates originated from multiple T. brucei strains from different genetic backgrounds, implying independent origins of T. evansi from T. brucei strains. This surprising finding further suggested that the acquisition of the ability of T. evansi to be transmitted mechanically, and thus the ability to escape the obligate link with the African tsetse fly vector, has occurred repeatedly. These findings, if confirmed, have epidemiological implications, as T. brucei strains from different genetic backgrounds can become either causative agents of a dangerous, cosmopolitan livestock disease or of a lethal human disease, like for T. b. rhodesiense.

Genetic diversity and population structure of the tsetse fly Glossina fuscipes fuscipes (Diptera: Glossinidae) in Northern Uganda: Implications for vector control.

Uganda is the only country where the chronic and acute forms of human African Trypanosomiasis (HAT) or sleeping sickness both occur and are separated by < 100 km in areas north of Lake Kyoga. In Uganda, Glossina fuscipes fuscipes is the main vector of the Trypanosoma parasites responsible for these diseases as well for the animal African Trypanosomiasis (AAT), or Nagana. We used highly polymorphic microsatellite loci and a mitochondrial DNA (mtDNA) marker to provide fine scale spatial resolution of genetic structure of G. f. fuscipes from 42 sampling sites from the northern region of Uganda where a merger of the two disease belts is feared. Based on microsatellite analyses, we found that G. f. fuscipes in northern Uganda are structured into three distinct genetic clusters with varying degrees of interconnectivity among them. Based on genetic assignment and spatial location, we grouped the sampling sites into four genetic units corresponding to northwestern Uganda in the Albert Nile drainage, northeastern Uganda in the Lake Kyoga drainage, western Uganda in the Victoria Nile drainage, and a transition zone between the two northern genetic clusters characterized by high level of genetic admixture. An analysis using HYBRIDLAB supported a hybrid swarm model as most consistent with tsetse genotypes in these admixed samples. Results of mtDNA analyses revealed the presence of 30 haplotypes representing three main haplogroups, whose location broadly overlaps with the microsatellite defined clusters. Migration analyses based on microsatellites point to moderate migration among the northern units located in the Albert Nile, Achwa River, Okole River, and Lake Kyoga drainages, but not between the northern units and the Victoria Nile drainage in the west. Effective population size estimates were variable with low to moderate sizes in most populations and with evidence of recent population bottlenecks, especially in the northeast unit of the Lake Kyoga drainage. Our microsatellite and mtDNA based analyses indicate that G. f. fuscipes movement along the Achwa and Okole rivers may facilitate northwest expansion of the Rhodesiense disease belt in Uganda. We identified tsetse migration corridors and recommend a rolling carpet approach from south of Lake Kyoga northward to minimize disease dispersal and prevent vector re-colonization. Additionally, our findings highlight the need for continuing tsetse monitoring efforts during and after control.

Effective population sizes of a major vector of human diseases, Aedes aegypti.

The effective population size (Ne ) is a fundamental parameter in population genetics that determines the relative strength of selection and random genetic drift, the effect of migration, levels of inbreeding, and linkage disequilibrium. In many cases where it has been estimated in animals, Ne is on the order of 10%-20% of the census size. In this study, we use 12 microsatellite markers and 14,888 single nucleotide polymorphisms (SNPs) to empirically estimate Ne in Aedes aegypti, the major vector of yellow fever, dengue, chikungunya, and Zika viruses. We used the method of temporal sampling to estimate Ne on a global dataset made up of 46 samples of Ae. aegypti that included multiple time points from 17 widely distributed geographic localities. Our Ne estimates for Ae. aegypti fell within a broad range (~25-3,000) and averaged between 400 and 600 across all localities and time points sampled. Adult census size (Nc) estimates for this species range between one and five thousand, so the Ne /Nc ratio is about the same as for most animals. These Ne values are lower than estimates available for other insects and have important implications for the design of genetic control strategies to reduce the impact of this species of mosquito on human health.