The distribution and spread of naturally occurring Medea selfish genetic elements in the United States

Selfish genetic elements (SGEs) are DNA sequences that are transmitted to viable offspring in greater than Mendelian frequencies. Medea SGEs occur naturally in some populations of red flour beetle (Tribolium castaneum) and are expected to increase in frequency within populations and spread among populations. The large‐scale U.S. distributions of Medea‐4 (M⁴) had been mapped based on samples from 1993 to 1995. We sampled beetles in 2011–2014 and show that the distribution of M⁴ in the United States is dynamic and has shifted southward. By using a genetic marker of Medea‐1 (M¹), we found five unique geographic clusters with high and low M¹ frequencies in a pattern not predicted by microsatellite‐based analysis of population structure. Our results indicate the absence of rigid barriers to Medea spread in the United States, so assessment of what factors have limited its current distribution requires further investigation. There is great interest in using synthetic SGEs, including synthetic Medea, to alter or suppress pest populations, but there is concern about unpredicted spread of these SGEs and potential for populations to become resistant to them. The finding of patchy distributions of Medea elements suggests that released synthetic SGEs cannot always be expected to spread uniformly, especially in target species with limited dispersal.     

Modeling the Dynamics of a Non-Limited and a Self-Limited Gene Drive System in Structured Aedes aegypti Populations

Recently there have been significant advances in research on genetic strategies to control populations of disease-vectoring insects. Some of these strategies use the gene drive properties of selfish genetic elements to spread physically linked anti-pathogen genes into local vector populations. Because of the potential of these selfish elements to spread through populations, control approaches based on these strategies must be carefully evaluated to ensure a balance between the desirable spread of the refractoriness-conferring genetic cargo and the avoidance of potentially unwanted outcomes such as spread to non-target populations. There is also a need to develop better estimates of the economics of such releases. We present here an evaluation of two such strategies using a biologically realistic mathematical model that simulates the resident Aedes aegypti mosquito population of Iquitos, Peru. One strategy uses the selfish element Medea, a non-limited element that could permanently spread over a large geographic area; the other strategy relies on Killer-Rescue genetic constructs, and has been predicted to have limited spatial and temporal spread. We simulate various operational approaches for deploying these genetic strategies, and quantify the optimal number of released transgenic mosquitoes needed to achieve definitive spread of Medea-linked genes and/or high frequencies of Killer-Rescue-associated elements. We show that for both strategies the most efficient approach for achieving spread of anti-pathogen genes within three years is generally to release adults of both sexes in multiple releases over time. Even though females in these releases should not transmit disease, there could be public concern over such releases, making the less efficient male-only release more practical. This study provides guidelines for operational approaches to population replacement genetic strategies, as well as illustrates the use of detailed spatial models to assist in safe and efficient implementation of such novel genetic strategies.     

Genetics at the verge of extinction: insights from the Iberian lynx

Population viability might become compromised by the loss of genetic diversity and the accumulation of inbreeding resulting from population decline and fragmentation. The Iberian lynx (Lynx pardinus) provides a paradigmatic example of a species at the verge of extinction, and because of the well‐documented and different demographic histories of the two remaining populations (Doñana and Andújar), it provides the opportunity to evaluate the performance of analytical methods commonly applied to recently declined populations. We used mitochondrial sequences and 36 microsatellite markers to evaluate the current genetic status of the species and to assess the genetic signatures of its past history. Mitochondrial diversity was extremely low with only two haplotypes, alternatively fixed in each population. Both remnant populations have low levels of genetic diversity at microsatellite markers, particularly the population from Doñana, and genetic differentiation between the two populations is high. Bayesian coalescent‐based methods suggest an earlier decline starting hundreds of years ago, while heterozygosity excess and M‐ratio tests did not provide conclusive and consistent evidence for recent bottlenecks. Also, a model of gene flow received overwhelming support over a model of pure drift. Results that are in conflict with the known recent demography of the species call for caution in the use of these methods, especially when no information on previous demographic history is available. Overall, our results suggest that current genetic patterns in the Iberian lynx are mainly the result of its recent decline and fragmentation and alerts on possible genetic risks for its persistence. Conservation strategies should explicitly consider this threat and incorporate an integrated genetic management of wild, captive and re‐introduced populations, including genetic restoration through translocations.     

Using soil seed banks to assess temporal patterns of genetic variation in invasive plant populations

Most research on the genetics of invasive plant species has focused on analyzing spatial differences among existing populations. Using a long‐established Gunnera tinctoria population from Ireland, we evaluated the potential of using plants derived from seeds associated with different soil layers to track genetic variation through time. This species and site were chosen because (1) G. tinctoria produces a large and persistent seed bank; (2) it has been present in this locality, Sraheens, for ~90 years; (3) the soil is largely undisturbed; and (4) the soil's age can be reliably determined radiometrically at different depths. Amplified fragment length polymorphic markers (AFLPs) were used to assess differences in the genetic structure of 75 individuals sampled from both the standing population and from four soil layers, which spanned 18 cm (estimated at ~90 years based on ²¹⁰Pb and ¹³⁷Cs dating). While there are difficulties in interpreting such data, including accounting for the effects of selection, seed loss, and seed migration, a clear pattern of lower total allele counts, percentage polymorphic loci, and genetic diversity was observed in deeper soils. The greatest percentage increase in the measured genetic variables occurred prior to the shift from the lag to the exponential range expansion phases and may be of adaptive significance. These findings highlight that seed banks in areas with long‐established invasive populations can contain valuable genetic information relating to invasion processes and as such, should not be overlooked.     

Contrasting genetic responses to population fragmentation in a coevolving fig and fig wasp across a mainland–island archipelago

Interacting species of pollinator–host systems, especially the obligate ones, are sensitive to habitat fragmentation, due to the nature of mutual dependence. Comparative studies of genetic structure can provide insights into how habitat fragmentation contributes to patterns of genetic divergence among populations of the interacting species. In this study, we used microsatellites to analyse genetic variation in Chinese populations of a typical mutualistic system – Ficus pumila and its obligate pollinator Wiebesia sp. 1 – in a naturally fragmented landscape. The plants and wasps showed discordant patterns of genetic variation and geographical divergence. There was no significant positive relationship in genetic diversity between the two species. Significant isolation‐by‐distance (IBD) patterns occurred across the populations of F. pumila and Wiebesia sp. 1 as whole, and IBD also occurred among island populations of the wasps, but not the plants. However, there was no significant positive relationship in genetic differentiation between them. The pollinator populations had significantly lower genetic variation in small habitat patches than in larger patches, and three island pollinator populations showed evidence of a recent bottleneck event. No effects of patch size or genetic bottlenecks were evident in the plant populations. Collectively, the results indicate that, in more fragmented habitats, the pollinators, but not the plants, have experienced reduced genetic variation. The contrasting patterns have multiple potential causes, including differences in longevity and hence number of generations experiencing fragmentation; different dispersal patterns, with the host's genes dispersed as seeds as well as a result of pollen dispersal via the pollinator; asymmetrical responses to fluctuations in partner populations; and co‐existence of a rare second pollinating wasp on some islands. These results indicate that strongly interdependent species may respond in markedly different ways to habitat fragmentation.     

Effects of iron deficiency on the growth and photosynthesis of three bloom‐forming cyanobacterial species isolated from Lake Taihu

Cyanobacterial blooms are found in many freshwater ecosystems around the world, but the effect of environmental factors on their growth and the proportion of species still require more investigation. In this study, the physiological responses of bloom‐forming cyanobacteria M icrocystis aeruginosa FACHB912, M icrocystis flos‐aquae FACHB1028 and P seudanabaena sp. FACHB1282 to iron deficiency were investigated. Their specific growth rates were found to decrease as the available iron concentration decreased. At low available iron concentrations of 1 × 10^?7 M (pFe 21.3) and 5 × 10^?8 M (pFe 21.6), M . aeruginosa had the lowest specific growth rate among three studied species. The cell sizes of M . flos‐aquae and Pseudanabaena sp. were significantly smaller under the lowest iron concentration. The chlorophyll a content of the three species decreased at the lowest iron concentration. The maximal relative electron transport rate, photosynthetic efficiency, and light‐saturation parameter of M . aeruginosa were lower than the other two cyanobacteria at pFe 21.3. Therefore, M . aeruginosa was the least able to adapt to iron deficiency. Under iron deficiency, the functional absorption cross‐section of PSII and electron transport rate on the acceptor side of PSII decreased in M . aeruginosa, while the connectivity factor between individual photosynthetic units increased in M . flos‐aquae, and the electron transport rate on the acceptor side of PSII and between PSII and PSI decreased in P seudanabaena sp. The ability to store iron was highest in M . flos‐aquae, followed by P seudanabaena sp. and M . aeruginosa. Thus, these results provide necessary information for detecting the role of iron in the succession of cyanobacterial species in Lake Taihu, the third largest freshwater lake in China, because all three species were isolated from this lake.     

Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species

Prioritizing and making efficient conservation plans for threatened populations requires information at both evolutionary and ecological timescales. Nevertheless, few studies integrate multidisciplinary approaches, mainly because of the difficulty for conservationists to assess simultaneously the evolutionary and ecological status of populations. Here, we sought to demonstrate how combining genetic and demographic analyses allows prioritizing and initiating conservation plans. To do so, we combined snapshot microsatellite data and a 30‐year‐long demographic survey on a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale. Our results revealed low levels of genetic diversity and weak effective population sizes (<63 individuals) in all populations. We further detected severe bottlenecks dating back to the last centuries (200–800 years ago), which may explain the differentiation of certain populations. The demographic survey revealed a general decrease in the spatial distribution and abundance of P. toxostoma over the last three decades. We conclude that demo‐genetic approaches are essential for (1) identifying populations for which both evolutionary and ecological extinction risks are high; and (2) proposing conservation plans targeted toward these at risk populations, and accounting for the evolutionary history of populations. We suggest that demo‐genetic approaches should be the norm in conservation practices.     

Signatures of de‐domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi‐extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo‐Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south‐eastern European countries (Kr?kopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro‐geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild‐type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of ‘de‐domestication’ process, and wild resources are challenged by a ‘domestication’ drift. Both need to be further investigated and managed.     

Mapping the genomic architecture of ecological speciation in the wild: does linkage disequilibrium hold the key?

The hunt for the genes underlying ecological speciation has now closed in on a number of candidates, but making the link from genotype to phenotype continues to pose a significant challenge. This is partly because genetic studies in many systems remain impeded by long generation times or an inability to perform controlled crosses. Now, in this issue of Molecular Ecology, Malek et al. (2012) demonstrate the utility of a novel admixture mapping approach that can be used to identify genomic regions contributing to adaptive trait divergence between natural populations. Remarkably, they validate their approach by mapping traits associated with mate choice in a wild limnetic and benthic threespine stickleback (Gasterosteus aculeatus) species pair, finding several loci associated with male nuptial coloration and shape. While this study benefited from tried‐and‐true microsatellites in a well‐characterized species with a detailed genetic map (and genome sequence), the field is quickly moving towards the use of next‐generation sequencing, especially for nonmodel systems. The ability to characterize molecular polymorphisms for any system suggests that molecular ecologists working on virtually any species may benefit from applying Malek et al.'s approach, if naturally admixed populations are available.     

Environmental gradients shape the genetic structure of two medicinal Salvia species in Jordan

• Environmental gradients, and particularly climatic variables, exert a strong influence on plant distribution and, potentially, population genetic diversity and differentiation. Differences in water availability can cause among‐population variation in ecological processes and can thus interrupt populations’ connectivity and isolate them environmentally. The present study examines the effect of environmental heterogeneity on plant populations due to environmental isolation unrelated to geographic distance.
• Using AFLP markers, we analyzed genetic diversity and differentiation among 12 Salvia spinosa populations and 13 Salvia syriaca populations from three phytogeographical regions (Mediterranean, Irano‐Turanian and Saharo‐Arabian) representing the extent of the species’ geographic range in Jordan. Differences in geographic location and climate were considered in the analyses.
• For both species, flowering phenology varied among populations and regions. Irano‐Turanian and Saharo‐Arabian populations had higher genetic diversity than Mediterranean populations, and genetic diversity increased significantly with increasing temperature. Genetic diversity in Salvia syriaca was affected by population size, while genetic diversity responded to drought in S. spinosa. For both species, high levels of genetic differentiation were found as well as two well‐supported phytogeographical groups of populations, with Mediterranean populations clustering in one group and the Irano‐Turanian and Saharo‐Arabian populations in another. Genetic distance was significantly correlated to environmental distance, but not to geographic distance.
• Our data indicate that populations from moist vs. arid environments are environmentally isolated, where environmental gradients affect their flowering phenology, limit gene flow and shape their genetic structure. We conclude that environmental heterogeneity may act as driver for the observed variation in genetic diversity.

     

Transposable elements in individual genotypes of Drosophila simulans

Transposable elements are abundant, dynamic components of the genome that affect organismal phenotypes and fitness. In Drosophila melanogaster, they have increased in abundance as the species spread out of Africa, and different populations differ in their transposable element content. However, very little is currently known about how transposable elements differ between individual genotypes, and how that relates to the population dynamics of transposable elements overall. The sister species of D. melanogaster, D. simulans, has also recently become cosmopolitan, and panels of inbred genotypes exist from cosmopolitan and African flies. Therefore, we can determine whether the differences in colonizing populations are repeated in D. simulans, what the dynamics of transposable elements are in individual genotypes, and how that compares to wild flies. After estimating copy number in cosmopolitan and African D. simulans, I find that transposable element load is higher in flies from cosmopolitan populations. In addition, transposable element load varies considerably between populations, between genotypes, but not overall between wild and inbred lines. Certain genotypes either contain active transposable elements or are more permissive of transposition and accumulate copies of particular transposable elements. Overall, it is important to quantify genotype‐specific transposable element dynamics as well as population averages to understand the dynamics of transposable element accumulation over time.     

Genome‐scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human‐induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide‐rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.     

Postglacial colonization history reflects in the genetic structure of natural populations of Festuca rubra in Europe

We conducted a large‐scale population genetic survey of genetic diversity of the host grass Festuca rubra s.l., which fitness can be highly dependent on its symbiotic fungus Epichloë festucae, to evaluate genetic variation and population structure across the European range. The 27 studied populations have previously been found to differ in frequencies of occurrence of the symbiotic fungus E. festucae and ploidy levels. As predicted, we found decreased genetic diversity in previously glaciated areas in comparison with nonglaciated regions and discovered three major maternal genetic groups: southern, northeastern, and northwestern Europe. Interestingly, host populations from Greenland were genetically similar to those from the Faroe Islands and Iceland, suggesting gene flow also between those areas. The level of variation among populations within regions is evidently highly dependent on the postglacial colonization history, in particular on the number of independent long‐distance seed colonization events. Yet, also anthropogenic effects may have affected the population structure in F. rubra. We did not observe higher fungal infection rates in grass populations with lower levels of genetic variability. In fact, the fungal infection rates of E. festucae in relation to genetic variability of the host populations varied widely among geographical areas, which indicate differences in population histories due to colonization events and possible costs of systemic fungi in harsh environmental conditions. We found that the plants of different ploidy levels are genetically closely related within geographic areas indicating independent formation of polyploids in different maternal lineages.     

Trends in bowhead whales in West Greenland: Aerial surveys vs. genetic capture‐recapture analyses

We contrast two methods for estimating the trends of bowhead whales (Balaena mysticetus) in West Greenland: (1) double platform visual aerial survey, corrected for missed sightings and the time the whales are available at the surface; and (2) a genetic capture‐recapture approach based on a 14‐yr‐long biopsy sampling program in Disko Bay. The aerial survey covered 39,000 km² and resulted in 58 sightings, yielding an abundance estimate of 744 whales (CV = 0.34, 95% CI: 357–1,461). The genetic method relied on determining sex, mitochondrial haplotypes and genotypes of nine microsatellite markers. Based on samples from a total of 427 individuals, with 11 recaptures from previous years in 2013, this resulted in an estimate of 1,538 whales (CV = 0.24, 95% CI: 827–2,249). While the aerial survey is considered a snapshot of the local spring aggregation in Disko Bay, the genetic approach estimates the abundance of the source of this aggregation. As the whales in Disko Bay primarily are adult females that do not visit the bay annually, the genetic method would presumably yield higher estimates. The studies indicate that an increase in abundance observed between 1998 and 2006 has leveled off.     

Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape

Habitat fragmentation and climate change are both prominent manifestations of global change, but there is little knowledge on the specific mechanisms of how climate change may modify the effects of habitat fragmentation, for example, by altering dynamics of spatially structured populations. The long‐term viability of metapopulations is dependent on independent dynamics of local populations, because it mitigates fluctuations in the size of the metapopulation as a whole. Metapopulation viability will be compromised if climate change increases spatial synchrony in weather conditions associated with population growth rates. We studied a recently reported increase in metapopulation synchrony of the Glanville fritillary butterfly (Melitaea cinxia) in the Finnish archipelago, to see if it could be explained by an increase in synchrony of weather conditions. For this, we used 23 years of butterfly survey data together with monthly weather records for the same period. We first examined the associations between population growth rates within different regions of the metapopulation and weather conditions during different life‐history stages of the butterfly. We then examined the association between the trends in the synchrony of the weather conditions and the synchrony of the butterfly metapopulation dynamics. We found that precipitation from spring to late summer are associated with the M. cinxia per capita growth rate, with early summer conditions being most important. We further found that the increase in metapopulation synchrony is paralleled by an increase in the synchrony of weather conditions. Alternative explanations for spatial synchrony, such as increased dispersal or trophic interactions with a specialist parasitoid, did not show paralleled trends and are not supported. The climate driven increase in M. cinxia metapopulation synchrony suggests that climate change can increase extinction risk of spatially structured populations living in fragmented landscapes by altering their dynamics.     

Skeletal adaptations and phylogeny of the oldest mole Eotalpa (Talpidae,Lipotyphla, Mammalia) from the UK Eocene: the beginning of fossoriality in moles

The oldest talpid, Eotalpa, was previously known only from isolated cheek teeth from the European late Middle Eocene to earliest Oligocene. Screenwashing of Late Eocene sediments of the Hampshire Basin, UK, has yielded cranial and postcranial elements: maxilla, dentary, ulna, metacarpals, distal tibia, astragalus, calcaneum, metatarsals and phalanges. In addition to M_1–2 myotodonty, typical talpid features are as follows: ulna with long medially curved olecranon and deep abductor fossa and astragalar body with lateral process. However, Eotalpa retains certain soricid‐like primitive states (M¹ preparacrista, P⁴ with prominent mesiolingual protocone lobe, strongly angled astragalar neck and calcaneum with no space for a cuboid medial process) not found in modern talpids. Eotalpa is more derived than the most primitive living talpid Uropsilus in having lost the M^1–2 talon shelf, developed a convex radial facet on the ulna, an incipient proximal olecranon crest, relatively shorter metapodials and depressed manual unguals. Its astragalus with medial trochlear ridge taller than the lateral one and massive medial plantar process is typical of the Lipotyphla. Eotalpa lacks synostosis of tibia and fibula, found in other Talpidae, Soricidae and Erinaceidae, suggesting that synostosis in these groups has been independently acquired. Cladistic analysis places Eotalpa as stem member of the Talpidae and shows that much homoplasy arose during the early evolution of the family. Ground dwelling in Eotalpa is indicated by the following: astragalus with a medially dipping head, curved in a single plane; calcaneum with distal peroneal process and strongly overlapping ectal and sustentacular facets; and matching sized ectal and sustentacular facets on calcaneum and astragalus. These features would have restricted ankle mobility. Ungual and metatarsal shape and ulnar structure suggest a primitive stage in fossorial evolution and argue against a semiaquatic precursor stage in talpid fossoriality. Shrew‐moles may represent a reversal to surface foraging rather than an intermediate stage in fossoriality.     

Evidence of Cotton leaf curl Burewala virus Variant and its Associate Betasatellite Causing Yellow Mosaic of Eggplant (Solanum melongena) in Pakistan

Eggplant (Solanum melongena L.) plants with severe leaf mosaic and mottling were found in a kitchen garden near cotton fields in Pakistan. Rolling Circle Amplification products from six of the naturally infected eggplant plants, subjected to PCR, successfully amplified expected products of 2.8 and 1.4 kb using begomovirus and betasatellite‐specific primers, respectively. Based on 99% nucleotide sequence identity, the virus was identified as a variant of Cotton leaf curl Burewala virus (CLCuBuV) (GenBank Accession No. HG428709). Likewise, the sequenced betasatellite with a maximum of 97% nucleotide sequence identity was recognized as a new variant of Cotton leaf curl Multan betasatellite (CLCuMuB^Mul) (GenBank Accession No. HG428708). The symptomatic induction of Cotton leaf curl disease in CLCuBuV susceptible cotton genotype CIM‐496 by back‐indexing further confirmed the presence of CLCuBuV in eggplant. This is the first report of CLCuBuV and its associate betasatellite in naturally infected plants of eggplant.     

Population genetics of Wolbachia‐infected,parthenogenetic and uninfected,sexual populations of Tetrastichus coeruleus (Hymenoptera: Eulophidae)

Wolbachia are endosymbiotic bacteria known to manipulate the reproduction of their hosts. These manipulations are expected to have consequences on the population genetics of the host, such as heterozygosity levels, genetic diversity and gene flow. The parasitoid wasp Tetrastichus coeruleus has populations that are infected with parthenogenesis‐inducing Wolbachia and populations that are not infected. We studied the population genetics of T. coeruleus between and within Wolbachia‐infected and uninfected populations, using nuclear microsatellites and mitochondrial DNA. We expected reduced genetic diversity in both DNA types in infected populations. However, migration and gene flow could introduce new DNA variants into populations. We therefore paid special attention to individuals with unexpected (genetic) characteristics. Based on nuclear and mitochondrial DNA, two genetic clusters were evident: a thelytokous cluster containing all Wolbachia‐infected, parthenogenetic populations and an arrhenotokous cluster containing all uninfected, sexual populations. Nuclear and mitochondrial DNA did not exhibit concordant patterns of variation, although there was reduced genetic diversity in infected populations for both DNA types. Within the thelytokous cluster, there was nuclear DNA variation, but no mitochondrial DNA variation. This nuclear DNA variation may be explained by occasional sex between infected females and males, by horizontal transmission of Wolbachia, and/or by novel mutations. Several females from thelytokous populations were uninfected and/or heterozygous for microsatellite loci. These unexpected characteristics may be explained by migration, by inefficient transmission of Wolbachia, by horizontal transmission of Wolbachia, and/or by novel mutations. However, migration has not prevented the build‐up of considerable genetic differentiation between thelytokous and arrhenotokous populations.     

SIGN EPISTASIS LIMITS EVOLUTIONARY TRADE‐OFFS AT THE CONFLUENCE OF SINGLE‐ AND MULTI‐CARBON METABOLISM IN METHYLOBACTERIUM EXTORQUENS AM1

Adaptation of one set of traits is often accompanied by attenuation of traits important in other selective environments, leading to fitness trade‐offs. The mechanisms that either promote or prevent the emergence of trade‐offs remain largely unknown, and are difficult to discern in most systems. Here, we investigate the basis of trade‐offs that emerged during experimental evolution of Methylobacterium extorquens AM1 to distinct growth substrates. After 1500 generations of adaptation to a multi‐carbon substrate, succinate (S), many lineages had lost the ability to use one‐carbon compounds such as methanol (M), generating a mixture of M⁺ and M^? evolved phenotypes. We show that trade‐offs in M^? strains consistently arise via antagonistic pleiotropy through recurrent selection for loss‐of‐function mutations to ftfL (formate‐tetrahydrofolate ligase), which improved growth on S while simultaneously eliminating growth on M. But if loss of FtfL was beneficial, why were M trade‐offs not found in all populations? We discovered that eliminating FtfL was not universally beneficial on S, as it was neutral or even deleterious in certain evolved lineages that remained M⁺. This suggests that sign epistasis with earlier arising mutations prevented the emergence of mutations that drove trade‐offs through antagonistic pleiotropy, limiting the evolution of metabolic specialists in some populations.