Adaptive introgression of anticoagulant rodent poison resistance by hybridization between old world mice

Hybrids of two destructive subterranean termites established in the field, revealing a potential for gene flow between species

Hybridization between invasive pest species may lead to significant genetic and economic impacts that require close monitoring. The two most invasive and destructive termite species worldwide, Coptotermes formosanus Shiraki and Coptotermes gestroi (Wasmann), have the potential for hybridization in the field. A three-year field survey conducted during the dispersal flight season of Coptotermes in Taiwan identified alates with atypical morphology, which were confirmed as hybrids of the two Coptotermes species using microsatellite and mitochondrial analyses. Out of 27,601 alates collected over three years, 4.4% were confirmed as hybrid alates, and some advanced hybrids (>F1 generations) were identified. The hybrid alates had a dispersal flight season that overlapped with the two parental species 13 out of 15 times. Most of the hybrid alates were females, implying that mating opportunities beyond F1 may primarily be possible through female hybrids. However, the incipient colony growth results from all potential mating combinations suggest that only backcross colonies with hybrid males could sometimes lead to brood development. The observed asymmetrical viability and fertility of hybrid alates may critically reduce the probability of advanced-hybrid colonies being established in the field.

Inference of selective forces on house mouse genomes during secondary contact in East Asia

The house mouse (Mus musculus), which is commensal to humans, has spread globally via human activities, leading to secondary contact between genetically divergent subspecies. This pattern of genetic admixture can provide insights into the selective forces at play in this well-studied model organism. Our analysis of 163 house mouse genomes, with a particular focus on East Asia, revealed substantial admixture between the subspecies castaneus and musculus, particularly in Japan and southern China. We revealed, despite the different level of autosomal admixture among regions, that all Y Chromosomes in the East Asian samples belonged to the musculus-type haplogroup, potentially explained by genomic conflict under sex-ratio distortion owing to varying copy numbers of ampliconic genes on sex chromosomes, Slx and Sly Our computer simulations, designed to replicate the observed scenario, show that the preferential fixation of musculus-type Y Chromosomes can be achieved with a slight increase in the male-to-female birth ratio. We also investigated the influence of selection on the posthybridization of the subspecies castaneus and musculus in Japan. Even though the genetic background of most Japanese samples closely resembles the subspecies musculus, certain genomic regions overrepresented the castaneus-like genetic components, particularly in immune-related genes. Furthermore, a large genomic block (∼2 Mbp) containing a vomeronasal/olfactory receptor gene cluster predominantly harbored castaneus-type haplotypes in the Japanese samples, highlighting the crucial role of olfaction-based recognition in shaping hybrid genomes.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

Widespread distribution of rodenticide resistance-conferring mutations in the Vkorc1 gene among house mouse populations in Portuguese Macaronesian islands and Iberian Atlantic areas

Growing evidence of widespread resistance to anticoagulant rodenticides (ARs) in house mice pose significant challenges to pest control efforts. First-generation ARs were introduced in the early 1950s but resistance to these emerged later that decade. Second-generation rodenticides were then developed, with resistance being reported in the late 1970s. Research has linked resistance to ARs with mutations in the Vkorc1 gene, leading to the use of more toxic and environmentally harmful compounds. In this study, 243 tail tips of house mice from mainland Portugal and Southern Spain, the Azores and Madeira archipelagos were analysed for all 3 exons of the Vkorc1 gene. Mutations L128S, Y139C, along with the so-called spretus genotype Vkorc1spr are considered responsible for reduced susceptibility of house mice to ARs. All these sequence variants were broadly detected throughout the sampling regions. Vkorc1spr was the most often recorded among mainland populations, whereas Y139C was nearly ubiquitous among the insular populations. In contrast, L128S was only detected in mainland Portugal and four islands of the Azores archipelago. All first generation ARs such as warfarin and coumatetralyl are deemed ineffective against all Vkorc1 variants identified in this study. Second-generation bromadiolone and difenacoum should also be discarded to control populations carrying Vkorc1spr, Y139C or L128S mutations. Inadequate use of ARs in regions where resistant animals have been found in large proportions will result in the spreading of rodenticide resistance among rodent populations through the positive selection of non-susceptible individuals. Consequently, ineffectiveness of rodent control will increase and potentiate environmental contamination, hazarding non-target wildlife through secondary poisoning. We highlight the need for Vkorc1 screening as a crucial tool in rodent management, aiding in the selection of the most appropriate control/eradication method in order to prevent misuse of these toxic biocides and the spread of rodenticide resistance among house mouse populations.

Local Effects of Nest-Boxes for Avian Predators over Common Vole Abundance during a Mid-Density Outbreak

At the end of the 20th century, the common vole (Microtus arvalis) colonized the practical totality of agricultural ecosystems in the northern sub-plateau of the Iberian Peninsula. To prevent crop damage, chemical control campaigns using anticoagulant rodenticides have been employed. This approach has a high environmental impact, and it has been banned in most countries in the European Union, including Spain. It is therefore essential to analyze alternative methods with lower environmental impacts. Here we explored the efficacy of biological control by avian predators to reduce vole abundance by providing nest-boxes in croplands. We used an indirect index based on the presence/absence of vole activity signs to measure the effect of nest-boxes on common vole abundance. We found that vole abundance was significantly lower near occupied nest-boxes at distances less than 180 m, where vole abundance increases progressively with increasing distance to the nearest nest-box. We also observed that the predatory pressure negatively affects the vole abundance at the end of the breeding period, considering the total number of fledglings. However, the effect of nest-boxes was highly variable depending on the study area and more limited in alfalfa fields, the optimal habitat for voles in agrarian ecosystems. Thus, nest-box supplementation would be a feasible measure for the biological control of the common vole in Mediterranean ecosystems, but it needs improvements for vole control in alfalfa fields within an integrated pest control program. We provide several recommendations to improve the performance of biological control in alfalfa fields.

From chip to SNP: Rapid development and evaluation of a targeted capture genotyping-by-sequencing approach to support research and management of a plaguing rodent

The management of invasive species has been greatly enhanced by population genetic analyses of multilocus single-nucleotide polymorphism (SNP) datasets that provide critical information regarding pest population structure, invasion pathways, and reproductive biology. For many applications there is a need for protocols that offer rapid, robust and efficient genotyping on the order of hundreds to thousands of SNPs, that can be tailored to specific study populations and that are scalable for long-term monitoring schemes. Despite its status as a model laboratory species, there are few existing resources for studying wild populations of house mice (Mus musculus spp.) that strike this balance between data density and laboratory efficiency. Here we evaluate the utility of a custom targeted capture genotyping-by-sequencing approach to support research on plaguing house mouse populations in Australia. This approach utilizes 3,651 hybridization capture probes targeting genome-wide SNPs identified from a sample of mice collected in grain-producing regions of southeastern Australia genotyped using a commercially available microarray platform. To assess performance of the custom panel, we genotyped wild caught mice (N = 320) from two adjoining farms and demonstrate the ability to correctly assign individuals to source populations with high confidence (mean >95%), as well as robust kinship inference within sites. We discuss these results in the context of proposed applications for future genetic monitoring of house mice in Australia.

Large-scale identification of rodenticide resistance in Rattus norvegicus and Mus musculus in the Netherlands based on Vkorc1 codon 139 mutations

BACKGROUND

Resistance to rodenticides has been reported globally and poses a considerable problem for efficacy in pest control. The most-documented resistance to rodenticides in commensal rodents is associated with mutations in the Vkorc1 gene, in particular in codon 139. Resistance to anticoagulant rodenticides has been reported in the Netherlands since 1989. A study from 2013 showed that 25% of 169 Norway rats (Rattus norvegicus) had a mutation at codon 139 of the Vkorc1 gene. To gain insight in the current status of rodenticide resistance amongst R. norvegicus and house mice Mus musculus in the Netherlands, we tested these rodents for mutations in codon 139 of the Vkorc1 gene. In addition, we collected data from pest controllers on their use of rodenticides and experience with rodenticide resistance.

RESULTS

A total of 1801 rodent samples were collected throughout the country consisting of 1404 R. norvegicus and 397 M. musculus. In total, 15% of R. norvegicus [95% confidence interval (CI): 13-17%] and 38% of M. musculus (95% CI: 33-43%) carried a genetic mutation at codon 139 of the Vkorc1 gene.

CONCLUSION

This study demonstrates genetic mutations at codon 139 of the Vkorc1 gene in M. musculus in the Netherlands. Resistance to anticoagulant rodenticides is present in R. norvegicus and M. musculus in multiple regions in the Netherlands. The results of this comprehensive study provide a baseline and facilitate trend analyses of Vkorc1 codon 139 mutations and evaluation of integrated pest management (IPM) strategies as these are enrolled in the Netherlands. © 2022 The Dutch Pest and Wildlife. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

The demographic history of house mice (Mus musculus domesticus) in eastern North America

The Western European house mouse (Mus musculus domesticus) is a widespread human commensal that has recently been introduced to North America. Its introduction to the Americas is thought to have resulted from the transatlantic movements of Europeans that began in the early 16th century. To study the details of this colonization history, we examine population structure, explore relevant demographic models, and infer the timing of divergence among house mouse populations in the eastern United States using published exome sequences from five North American populations and two European populations. For North American populations of house mice, levels of nucleotide variation were lower, and low-frequency alleles were less common than for European populations. These patterns provide evidence of a mild bottleneck associated with the movement of house mice into North America. Several analyses revealed that one North American population is genetically admixed, which indicates at least two source populations from Europe were independently introduced to eastern North America. Estimated divergence times between North American and German populations ranged between ∼1,000 and 7,000 years ago and overlapped with the estimated divergence time between populations from Germany and France. Demographic models comparing different North American populations revealed that these populations diverged from each other mostly within the last 500 years, consistent with the timing of the arrival of Western European settlers to North America. Together, these results support a recent introduction of Western European house mice to eastern North America, highlighting the effects of human migration and colonization on the spread of an invasive human commensal.

Gene drive in species complexes: defining target organisms

Engineered gene drives, which bias their own inheritance to increase in frequency in target populations, are being developed to control mosquito malaria vectors. Such mosquitoes can belong to complexes of both vector and nonvector species that can produce fertile interspecific hybrids, making vertical gene drive transfer (VGDT) to sibling species biologically plausible. While VGDT to other vectors could positively impact human health protection goals, VGDT to nonvectors might challenge biodiversity ones. Therefore, environmental risk assessment of gene drive use in species complexes invites more nuanced considerations of target organisms and nontarget organisms than for transgenes not intended to increase in frequency in target populations. Incorporating the concept of target species complexes offers more flexibility when assessing potential impacts from VGDT.

A Phylogenomic Assessment of Processes Underpinning Convergent Evolution in Open-Habitat Chats

Insights into the processes underpinning convergent evolution advance our understanding of the contributions of ancestral, introgressed, and novel genetic variation to phenotypic evolution. Phylogenomic analyses characterizing genome-wide gene tree heterogeneity can provide first clues about the extent of ILS and of introgression and thereby into the potential of these processes or (in their absence) the need to invoke novel mutations to underpin convergent evolution. Here, we were interested in understanding the processes involved in convergent evolution in open-habitat chats (wheatears of the genus Oenanthe and their relatives). To this end, based on whole-genome resequencing data from 50 taxa of 44 species, we established the species tree, characterized gene tree heterogeneity, and investigated the footprints of ILS and introgression within the latter. The species tree corroborates the pattern of abundant convergent evolution, especially in wheatears. The high levels of gene tree heterogeneity in wheatears are explained by ILS alone only for 30% of internal branches. For multiple branches with high gene tree heterogeneity, D-statistics and phylogenetic networks identified footprints of introgression. Finally, long branches without extensive ILS between clades sporting similar phenotypes provide suggestive evidence for the role of novel mutations in the evolution of these phenotypes. Together, our results suggest that convergent evolution in open-habitat chats involved diverse processes and highlight that phenotypic diversification is often complex and best depicted as a network of interacting lineages.

MaLAdapt Reveals Novel Targets of Adaptive Introgression From Neanderthals and Denisovans in Worldwide Human Populations

Adaptive introgression (AI) facilitates local adaptation in a wide range of species. Many state-of-the-art methods detect AI with ad-hoc approaches that identify summary statistic outliers or intersect scans for positive selection with scans for introgressed genomic regions. Although widely used, approaches intersecting outliers are vulnerable to a high false-negative rate as the power of different methods varies, especially for complex introgression events. Moreover, population genetic processes unrelated to AI, such as background selection or heterosis, may create similar genomic signals to AI, compromising the reliability of methods that rely on neutral null distributions. In recent years, machine learning (ML) methods have been increasingly applied to population genetic questions. Here, we present a ML-based method called MaLAdapt for identifying AI loci from genome-wide sequencing data. Using an Extra-Trees Classifier algorithm, our method combines information from a large number of biologically meaningful summary statistics to capture a powerful composite signature of AI across the genome. In contrast to existing methods, MaLAdapt is especially well-powered to detect AI with mild beneficial effects, including selection on standing archaic variation, and is robust to non-AI selective sweeps, heterosis from deleterious mutations, and demographic misspecification. Furthermore, MaLAdapt outperforms existing methods for detecting AI based on the analysis of simulated data and the validation of empirical signals through visual inspection of haplotype patterns. We apply MaLAdapt to the 1000 Genomes Project human genomic data and discover novel AI candidate regions in non-African populations, including genes that are enriched in functionally important biological pathways regulating metabolism and immune responses.

Investigation of anticoagulant rodenticide resistance induced by Vkorc1 mutations in rodents in Lebanon

Anticoagulant rodenticides (AR) remain the most effective chemical substances used to control rodents in order to limit their agricultural and public health damage in both rural and urban environments. The emergence of genetically based resistance to AR worldwide has threatened effective rodent control. This study gives a first overview of the distribution and frequency of single nucleotide polymorphism in the vitamin K epoxide reductase subcomponent 1 (Vkorc1) gene in rodents in Lebanon. In the Mus genus, we detected two missense mutations Leu128Ser and Tyr139Cys, that confer resistance to anticoagulant rodenticides in house mice and a new missense mutation Ala72Val in the Mus macedonicus species, not previously described. In the Rattus genus, we found one missense mutation Leu90Ile in the roof rat and one missense mutation Ser149Ile in the Norway rat. This is the first study to demonstrate potential resistance to AR in Lebanese rodents and therefore it provides data to pest control practitioners to choose the most suitable AR to control rodents in order to keep their efficacy.

VKORC1 single nucleotide polymorphisms in rodents in Spain

Rodents are considered one of the animal pests with the greatest impact on agricultural production and public health. Anticoagulant rodenticides (ARs), used as one of the most effective ways to control rodent populations worldwide, inhibit the vitamin K 2,3-epoxide reductase (VKORC1) enzyme involved in blood coagulation. Resistances to ARs are mainly associated with mutations or single nucleotide polymorphisms (SNPs) in the vkorc1 gene. Since the information on this subject is scarce in Spain, we monitored and discovered rodent SNPs that could favour genetic resistance in its populations. For that, more than 200 samples of stools and tails from brown rat (Rattus norvegicus), black rat (Rattus rattus) and mouse (Mus musculus) were collected from 12 Spanish regions previously identified with low AR efficacy in coordination with the National Association of Environmental Sanitation Companies (ANECPLA) and the managing entities of four locations. We then sequenced their vkorc1 exon 3 corresponding genomic DNA. We identified genotypic vkorc1 variations corresponding to amino acid changes at the VKORC1 protein at the S149I - S149T and the E155K - E155Q mutations, depending on the rodent species. Computational analysis of binding predictions found out that the brown rat S149I mutation predicted a high reduction of the binding affinity of chlorophacinone and brodifacoum ARs while, the black rat S149T, E155K and E155Q mutations slightly reduced bromadiolone AR binding. These results suggest that these mutations may be one of the causes of the increased resistance to those ARs.

Revisiting the Baltimore-Bullock's Oriole hybrid zone reveals changing plumage colour in Bullock's Orioles

Hybrid zones are dynamic areas where populations of two or more interbreeding species may change through an influx of novel genetic material resulting from hybridization or selection on standing genetic variation. Documenting changes in populations through time, however, is challenging because repeated samples are often missing or because long-term storage can affect trait morphologies, especially colour traits that may fade through time. We document a change in carotenoid-based orange breast feathers of Bullock's Orioles (Icterus bullockii) from the Great Plains hybrid zone, USA. Contemporary Bullock's Orioles are more orange than historic individuals from the same location sampled approximately 60 years ago. Spectrophotometry revealed that contemporary Bullock's Orioles resemble orange colour profiles of Baltimore Orioles (I. galbula), the species with which they hybridize. Fading or changes in diet hypotheses do not appear to explain the shift in colour we report for Bullock's Orioles. We propose that these changes in colour are facilitated through introgression with Baltimore Orioles, and favoured by females that choose brighter, more orange males. Our study highlights the long memory of natural history collections and how they offer new insights to the dynamic roll of hybrid zones in trait evolution between interacting species.

Widespread resistance to anticoagulant rodenticides in Mus musculus domesticus in the city of Barcelona

Control of rodent populations is a big challenge because of the rapid evolution of resistance to commonly used rodenticides and the collateral negative impacts that these products may have on biodiversity. Second-generation anticoagulants are very efficient but different single nucleotide polymorphisms (SNPs) in the Vkorc1 gene may confer resistance in rodents. We sequenced exons 1, 2 and 3 of the Vkorc1 gene from 111 mice (Mus musculus domesticus) captured across the city of Barcelona and found SNPs associated with resistance to first- and second-generation anticoagulants in all of them. Although most of the SNPs were associated with resistance to bromadiolone, we also found SNPs associated with resistance to brodifacoum. Out of all the individuals analyzed, 94.59 % carried mutations associated to introgression events with Mus spretus, a sympatric rodent species. Currently most of the chemical products for rodent control commercialized in the area are based on bromadiolone, although recent public control campaigns have already shifted to other products. Thus, the widespread occurrence of resistant mice to bromadiolone represents a challenge for rodent control in Barcelona and may increase the risk of secondary poisoning of animals preying on this species. Public health managers, pest control companies and citizens should be aware that the use of bromadiolone based products is ineffective and represents a risk for the environment, including human and animal health.

Anticoagulant rodenticide blood-clotting dose-responses and resistance factors for Tyrosine139Cysteine (Y139C) heterozygous- and homozygous-resistant house mice (Mus musculus)

BACKGROUND

The house mouse (Mus musculus) is a globally distributed rodent pest species against which anticoagulant rodenticides are widely used for the protection of human and animal health and the conservation of threatened wildlife. Anticoagulant-resistant house mice have been known for more than half a century. A house mouse strain was developed in the laboratory that was homozygous resistant for the single nucleotide polymorphism (SNP) Tyrosine139Cysteine (Y139C) and, subsequently, heterozygous resistant animals were produced from this strain by crossing with the homozygous susceptible strain.

RESULTS

Using blood clotting response tests, resistance factors at the ED₅₀ level in the homozygous resistant strain for the first-generation anticoagulants warfarin, chlorophacinone, diphacinone and coumatetralyl were in the range 31.5 to 628.0 for males (M) and 21.6 to 628.0 for females (F), thus indicating that Y139C house mice are substantially resistant to all these substances. Resistance factors at the ED₅₀ level for the homozygous strain generated against the second-generation compounds were: brodifacoum (M, 1.7; F, 1.9), bromadiolone (M, 16.6; F, 21.0), difenacoum (M, 1.2; F, 2.7), difethialone (M, 1.5; F, 1.5), and flocoumafen (M, 0.9; F, 1.2). Equivalent values for the heterozygous strain were: brodifacoum (M, 1.6; F, 1.4), bromadiolone (M, 5.6; F, 6.5), difenacoum (M, 1.0; F, 1.3), difethialone (M, 1.1; F, 1.1), flocoumafen (M, 0.9; F, 1.1).

CONCLUSION

Y139C SNP homozygous resistant mice are more resistant to anticoagulants than heterozygous resistant animals. All first-generation anticoagulants are highly resisted and, among the second-generation compounds, Y139C mice are resistant to bromadiolone and sometimes to difenacoum. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Counteracting forces of introgressive hybridization and interspecific competition shape the morphological traits of cryptic Iberian Eptesicus bats

Cryptic species that coexist in sympatry are likely to simultaneously experience strong competition and hybridization. The first phenomenon would lead to character displacement, whereas the second can potentially promote morphological similarity through adaptive introgression. The main goal of this work was to investigate the effect of introgressive hybridization on the morphology of cryptic Iberian Eptesicus bats when facing counteracting evolutionary forces from interspecific competition. We found substantial overlap both in dentition and in wing morphology traits, though mainly in individuals in sympatry. The presence of hybrids contributes to a fifth of this overlap, with hybrids showing traits with intermediate morphometry. Thus, introgressive hybridization may contribute to species adaptation to trophic and ecological space responding directly to the macro-habitats characteristics of the sympatric zone and to local prey availability. On the other hand, fur shade tended to be browner and brighter in hybrids than parental species. Colour differences could result from partitioning of resources as an adaptation to environmental factors such as roost and microhabitats. We argue that a balance between adaptive introgression and niche partitioning shapes species interactions with the environment through affecting morphological traits under selection.

Accelerated evolution of Vkorc1 in desert rodent species reveals genetic preadaptation to anticoagulant rodenticides

BACKGROUND

Some rodent species living in arid areas show elevated physiological tolerance to anti-vitamin K rodenticides (AVKs), which seems to be due to some unknown selective pressures that rodents may experience in desert habitats. Genes involved in the ϒ-carboxylation of blood coagulation, including vitamin K epoxide reductase complex, subunit 1 (Vkorc1), ϒ-glutamyl-carboxylase (Ggcx) and NAD(P)H quinone one dehydrogenase (Nqo1) are associated with anticoagulant resistance, or some levels of elevated tolerance, in rodents. To detect whether the DNA sequences of the three genes are also under natural selection in the desert rodent species, we analyzed the Vkorc1, Ggcx and Nqo1 genes of the desert rodents and compared them with other rodent species.

RESULTS

We found an accelerated evolutionary rate in Vkorc1 of desert rodents, especially in Mus spretus, Nannospalax galili and Psammomys obesus. By contrast, signals of positive selection were absent for Ggcx and Nqo1 in all species. Mapping the amino acid variations on the VKORC1 protein three-dimensional model suggested most interspecific amino acid variations occur on the outer surface of the VKORC1 pocket, whereas most intraspecific amino acid changes and known AVK resistance mutations occurred on the inner surface and endoplasmic reticulum luminal loop regions. Some desert-species-specific amino acid variations were found on the positions where known resistance mutations occurred, indicating these variations might be related to the elevated physical tolerance to AVKs in desert rodents.

CONCLUSION

The evolution of Vkorc1 has been accelerated in some desert rodent species, indicating genetic preadaptation to anticoagulant rodenticides. Positive selection and relaxed selection have been detected in Psammomys obesus and Nannospalax galili, indicating the two rodent species might also show tolerance to AVKs, which needs further verification. © 2022 Society of Chemical Industry.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

Chromosome size affects sequence divergence between species through the interplay of recombination and selection

The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades-Peromyscus rodents, Mus mice, and great apes-from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.

Genomes of hybrid lineages are mosaics of those of their parent species and harbour variation that has the potential to facilitate adaptation when hybrids encounter diverse environments. However, genetic incompatibilities between parental species can also act to limit possible combinations of parental alleles, constraining hybrid genome formation. What is the relative importance of selection and constraints in form of admixture proportions and genomic architecture in this process? We investigated this in the Italian sparrow, a hybrid species resulting from past hybridization between the house and Spanish sparrow. Using three independent hybrid lineages, we addressed how their genomes, harbouring different parental combinations, have evolved. We examined the roles of selection due to divergent local adaptation, recombination and purging of genetic incompatibilities in predicting differentiation. We found that selection against incompatibilities may constrain hybrid genome composition. In addition, signals of local selection as well as estimates of differentiation were correlated across populations, and outliers were shared among the hybrid lineages more often than expected by chance. Overall, our results suggest that in the Italian sparrow selection interacts with constraints linked to genetic incompatibilities affecting which sections of the genome can readily diverge among hybrid lineages.

Bidirectional introgression between Mus musculus domesticus and Mus spretus

Introgressed variants from other species can be an important source of genetic variation because they may arise rapidly, can include multiple mutations on a single haplotype, and have often been pretested by selection in the species of origin. Although introgressed alleles are generally deleterious, several studies have reported introgression as the source of adaptive alleles—including the rodenticide-resistant variant of Vkorc1 that introgressed from Mus spretus into European populations of Mus musculus domesticus. Here, we conducted bidirectional genome scans to characterize introgressed regions into one wild population of M. spretus from Spain and three wild populations of M. m. domesticus from France, Germany, and Iran. Despite the fact that these species show considerable intrinsic postzygotic reproductive isolation, introgression was observed in all individuals, including in the M. musculus reference genome (GRCm38). Mus spretus individuals had a greater proportion of introgression compared with M. m. domesticus, and within M. m. domesticus, the proportion of introgression decreased with geographic distance from the area of sympatry. Introgression was observed on all autosomes for both species, but not on the X-chromosome in M. m. domesticus, consistent with known X-linked hybrid sterility and inviability genes that have been mapped to the M. spretus X-chromosome. Tract lengths were generally short with a few outliers of up to 2.7 Mb. Interestingly, the longest introgressed tracts were in olfactory receptor regions, and introgressed tracts were significantly enriched for olfactory receptor genes in both species, suggesting that introgression may be a source of functional novelty even between species with high barriers to gene flow.

Prevalence and Adaptive Impact of Introgression

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

Varying Intensities of Introgression Obscure Incipient Venom-Associated Speciation in the Timber Rattlesnake (Crotalus horridus)

Ecologically divergent selection can lead to the evolution of reproductive isolation through the process of ecological speciation, but the balance of responsible evolutionary forces is often obscured by an inadequate assessment of demographic history and the genetics of traits under selection. Snake venoms have emerged as a system for studying the genetic basis of adaptation because of their genetic tractability and contributions to fitness, and speciation in venomous snakes can be associated with ecological diversification such as dietary shifts and corresponding venom changes. Here, we explored the neurotoxic (type A)–hemotoxic (type B) venom dichotomy and the potential for ecological speciation among Timber Rattlesnake (Crotalus horridus) populations. Previous work identified the genetic basis of this phenotypic difference, enabling us to characterize the roles geography, history, ecology, selection, and chance play in determining when and why new species emerge or are absorbed. We identified significant genetic, proteomic, morphological, and ecological/environmental differences at smaller spatial scales, suggestive of incipient ecological speciation between type A and type B C. horridus. Range-wide analyses, however, rejected the reciprocal monophyly of venom type, indicative of varying intensities of introgression and a lack of reproductive isolation across the range. Given that we have now established the phenotypic distributions and ecological niche models of type A and B populations, genome-wide data are needed and capable of determining whether type A and type B C. horridus represent distinct, reproductively isolated lineages due to incipient ecological speciation or differentiated populations within a single species.

Finding unknown species in the genomes of extant species

Although many species have gone extinct, their genetic components might exist in extant species because of ancient hybridization. Via advances in genome sequencing and development of modern population genetics, one can find the legacy of unknown or extinct species in the context of available genomes from extant species. Such discovery can be used as a strategy to search for hidden species or fossils in conservation biology and archeology, gain novel insight into complex evolutionary history, and provide the new sources of genetic variation for breeding and trait improvement in agriculture.

The genic view of hybridization in the Anthropocene

Human impact is noticeable around the globe, indicating that a new era might have begun: the Anthropocene. Continuing human activities, including land-use changes, introduction of non-native species and rapid climate change, are altering the distributions of countless species, often giving rise to human-mediated hybridization events. While the interbreeding of different populations or species can have detrimental effects, such as genetic extinction, it can be beneficial in terms of adaptive introgression or an increase in genetic diversity. In this paper, I first review the different mechanisms and outcomes of anthropogenic hybridization based on literature from the last five years (2016-2020). The most common mechanisms leading to the interbreeding of previously isolated taxa include habitat change (51% of the studies) and introduction of non-native species (34% intentional and 19% unintentional). These human-induced hybridization events most often result in introgression (80%). The high incidence of genetic exchange between the hybridizing taxa indicates that the application of a genic view of speciation (and introgression) can provide crucial insights on how to address hybridization events in the Anthropocene. This perspective considers the genome as a dynamic collection of genetic loci with distinct evolutionary histories, giving rise to a heterogenous genomic landscape in terms of genetic differentiation and introgression. First, understanding this genomic landscape can lead to a better selection of diagnostic genetic markers to characterize hybrid populations. Second, describing how introgression patterns vary across the genome can help to predict the likelihood of negative processes, such as demographic and genetic swamping, as well as positive outcomes, such as adaptive introgression. It is especially important to not only quantify how much genetic material introgressed, but also what has been exchanged. Third, comparing introgression patterns in pre-Anthropocene hybridization events with current human-induced cases might provide novel insights into the likelihood of genetic swamping or species collapse during an anthropogenic hybridization event. However, this comparative approach remains to be tested before it can be applied in practice. Finally, the genic view of introgression can be combined with conservation genomic studies to determine the legal status of hybrids and take appropriate measures to manage anthropogenic hybridization events. The interplay between evolutionary and conservation genomics will result in the constant exchange of ideas between these fields which will not only improve our knowledge on the origin of species, but also how to conserve and protect them.

Human-modified canids in human-modified landscapes: The evolutionary consequences of hybridization for grey wolves and free-ranging domestic dogs

Introgressive hybridization between domestic animals and their wild relatives is an indirect form of human-induced evolution, altering gene pools and phenotypic traits of wild and domestic populations. Although this process is well documented in many taxa, its evolutionary consequences are poorly understood. In this study, we assess introgression patterns in admixed populations of Eurasian wolves and free-ranging domestic dogs (FRDs), identifying chromosomal regions with significantly overrepresented hybrid ancestry and assessing whether genes located within these regions show signatures of selection. Although the dog admixture proportion in West Eurasian wolves (2.7%) was greater than the wolf admixture proportion in FRDs (0.75%), the number and average length of chromosomal blocks showing significant overrepresentation of hybrid ancestry were smaller in wolves than FRDs. In wolves, 6% of genes located within these blocks showed signatures of positive selection compared to 23% in FRDs. We found that introgression from wolves may provide a considerable adaptive advantage to FRDs, counterbalancing some of the negative effects of domestication, which can include reduced genetic diversity and excessive tameness. In wolves, introgression from FRDs is mostly driven by drift, with a small number of positively selected genes associated with brain function and behaviour. The predominance of drift may be the consequence of small effective size of wolf populations, which reduces efficiency of selection for weakly advantageous or against weakly disadvantageous introgressed variants. Small wolf population sizes result largely from human-induced habitat loss and hunting, thus linking introgression rates to anthropogenic processes. Our results imply that maintenance of large population sizes should be an important element of wolf management strategies aimed at reducing introgression rates of dog-derived variants.

Demographic modeling informs functional connectivity and management interventions in Graham’s beardtongue

Functional connectivity (i.e., the movement of individuals across a landscape) is essential for the maintenance of genetic variation and persistence of rare species. However, illuminating the processes influencing functional connectivity and ultimately translating this knowledge into management practice remains a fundamental challenge. Here, we combine various population structure analyses with pairwise, population-specific demographic modeling to investigate historical functional connectivity in Graham’s beardtongue (Penstemon grahamii), a rare plant narrowly distributed across a dryland region of the western US. While principal component and population structure analyses indicated an isolation-by-distance pattern of differentiation across the species’ range, spatial inferences of effective migration exposed an abrupt shift in population ancestry near the range center. To understand these seemingly conflicting patterns, we tested various models of historical gene flow and found evidence for recent admixture (~ 3400 generations ago) between populations near the range center. This historical perspective reconciles population structure patterns and suggests management efforts should focus on maintaining connectivity between these previously isolated lineages to promote the ongoing transfer of genetic variation. Beyond providing species-specific knowledge to inform management options, our study highlights how understanding demographic history may be critical to guide conservation efforts when interpreting population genetic patterns and inferring functional connectivity.

Colonization and Authentication of the Pyrethroid-Resistant Anopheles gambiae s.s. Muleba-Kis Strain; an Important Test System for Laboratory Screening of New Insecticides

Simple Summary

Malaria control and prevention have traditionally relied on the use of insecticides in the form of treated bed nets or residual spraying in households. However, scaling up of these interventions—based on few available insecticide classes—resulted in the development and spread of insecticide resistance in malaria-transmitting mosquitoes. There is therefore an urgent need for introducing and applying new insecticides that are effective against these mosquitoes. Laboratories tasked with evaluating the efficacy of novel insecticides need to establish a large colony of resistant mosquitoes. In this study, we report the procedures used and challenges faced during the establishment and maintenance of a resistant mosquito strain in the laboratory which reflects the characteristics of the wild-resistant mosquito populations found in East Africa.

Abstract

Background: The emergence and spread of insecticide resistance in malaria vectors to major classes of insecticides call for urgent innovation and application of insecticides with novel modes of action. When evaluating new insecticides for public health, potential candidates need to be screened against both susceptible and resistant mosquitoes to determine efficacy and to identify potential cross-resistance to insecticides currently used for mosquito control. The challenges and lessons learned from establishing, maintaining, and authenticating the pyrethroid-resistant An. gambiae s.s. Muleba-Kis strain at the KCMUCo-PAMVERC Test Facility are described in this paper. Methods: Male mosquitoes from the F₁ generation of wild-pyrethroid resistant mosquitoes were cross-bred with susceptible female An. gambiae s.s. Kisumu laboratory strain followed by larval selection using a pyrethroid insecticide solution. Periodic screening for phenotypic and genotypic resistance was done. WHO susceptibility tests and bottle bioassays were used to assess the phenotypic resistance, while Taqman™ assays were used to screen for known target-site resistance alleles (kdr and ace-1). Additionally, the strains were periodically assessed for quality control by monitoring adult weight and wing length. Results: By out-crossing the wild mosquitoes with an established lab strain, a successful resistant insectary colony was established. Intermittent selection pressure using alphacypermethrin has maintained high kdr mutation (leucine-serine) frequencies in the selected colony. There was consistency in the wing length and weight measurements from the year 2016 to 2020, with the exception that one out of four years was significantly different. Mean annual wing length varied between 0.0142–0.0028 mm compared to values obtained in 2016, except in 2019 where it varied by 0.0901 mm. Weight only varied by approximately 0.001 g across four years, except in 2017 where it differed by 0.005 g. Routine phenotypic characterization on Muleba-Kis against pyrethroids using the WHO susceptibility test indicated high susceptibility when type I pyrethroids were used compared to type II pyrethroids. Dynamics on susceptibility status also depended on the lapse time when the selection was last done. Conclusions: This study described the procedure for introducing, colonizing, and maintaining a resistant An. gambiae s.s. strain in the laboratory with leucine to serine substitution kdr allele which reflects the features of the wild-resistant population in East Africa. Challenges in colonizing a wild-resistant mosquito strain were overcome by out-crossing between mosquito strains of desired traits followed by intermittent insecticide selection at the larval stage to select for the resistant phenotype.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O2 transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O2 affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.

Exchanging fluids The sociocultural implications of microbial, cultural, and ethnic admixture in Latin America

Knowledge of evolutionary influences on patterns of human mating, social interactions, and differential health is increasing, yet these insights have rarely been applied to historical analyses of human population dynamics. The genetic and evolutionary forces behind biases in interethnic mating and in the health of individuals of different ethnic groups in Latin America and the Caribbean since the European colonization of America are still largely ignored. We discuss how historical and contemporary sociocultural interactions and practices are strongly influenced by population-level evolutionary forces. Specifically, we discuss the historical implications of functional (de facto) polygyny, sex-biased admixture, and assortative mating in Latin America. We propose that these three evolutionary mechanisms influenced mating patterns, shaping the genetic and cultural landscape across Latin America and the Caribbean. Further, we discuss how genetic differences between the original populations that migrated at different times into Latin America contributed to their accommodation to and survival in the different local ecologies and interethnic interactions. Relevant medical and social implications follow from the genetic and cultural changes reviewed.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

Inferring Adaptive Introgression Using Hidden Markov Models

Adaptive introgression-the flow of adaptive genetic variation between species or populations-has attracted significant interest in recent years and it has been implicated in a number of cases of adaptation, from pesticide resistance and immunity, to local adaptation. Despite this, methods for identification of adaptive introgression from population genomic data are lacking. Here, we present Ancestry_HMM-S, a hidden Markov model-based method for identifying genes undergoing adaptive introgression and quantifying the strength of selection acting on them. Through extensive validation, we show that this method performs well on moderately sized data sets for realistic population and selection parameters. We apply Ancestry_HMM-S to a data set of an admixed Drosophila melanogaster population from South Africa and we identify 17 loci which show signatures of adaptive introgression, four of which have previously been shown to confer resistance to insecticides. Ancestry_HMM-S provides a powerful method for inferring adaptive introgression in data sets that are typically collected when studying admixed populations. This method will enable powerful insights into the genetic consequences of admixture across diverse populations. Ancestry_HMM-S can be downloaded from https://github.com/jesvedberg/Ancestry_HMM-S/.

Independent evolution toward larger body size in the distinctive Faroe Island mice

Most phenotypic traits in nature involve the collective action of many genes. Traits that evolve repeatedly are particularly useful for understanding how selection may act on changing trait values. In mice, large body size has evolved repeatedly on islands and under artificial selection in the laboratory. Identifying the loci and genes involved in this process may shed light on the evolution of complex, polygenic traits. Here, we have mapped the genetic basis of body size variation by making a genetic cross between mice from the Faroe Islands, which are among the largest and most distinctive natural populations of mice in the world, and a laboratory mouse strain selected for small body size, SM/J. Using this F2 intercross of 841 animals, we have identified 111 loci controlling various aspects of body size, weight and growth hormone levels. By comparing against other studies, including the use of a joint meta-analysis, we found that the loci involved in the evolution of large size in the Faroese mice were largely independent from those of a different island population or other laboratory strains. We hypothesize that colonization bottleneck, historical hybridization, or the redundancy between multiple loci have resulted in the Faroese mice achieving an outwardly similar phenotype through a distinct evolutionary path.

Rapid parallel adaptation despite gene flow in silent crickets

Gene flow is predicted to impede parallel adaptation via de novo mutation, because it can introduce pre-existing adaptive alleles from population to population. We test this using Hawaiian crickets (Teleogryllus oceanicus) in which 'flatwing' males that lack sound-producing wing structures recently arose and spread under selection from an acoustically-orienting parasitoid. Morphometric and genetic comparisons identify distinct flatwing phenotypes in populations on three islands, localized to different loci. Nevertheless, we detect strong, recent and ongoing gene flow among the populations. Using genome scans and gene expression analysis we find that parallel evolution of flatwing on different islands is associated with shared genomic hotspots of adaptation that contain the gene doublesex, but the form of selection differs among islands and corresponds to known flatwing demographics in the wild. We thus show how parallel adaptation can occur on contemporary timescales despite gene flow, indicating that it could be less constrained than previously appreciated.

Polyploidy and high environmental tolerance increase the invasive success of plants

Ploidy level and genome size (GS) could affect the invasive capacity of plants, although these parameters can be contradictory. While small GS seems to favor dispersion, polyploidy-which increases the GS-also seems to favor it. Using a phylogenetic path analysis, we evaluated the effects of both factors on the environmental tolerance and invasive success of plants. We selected 99 invasive plant species from public online databases and gathered information about invasive capacity (number of non-original countries in which each species occurs), tolerance to environmental factors, ploidy level, and GS. The invasive capacity varied depending on the ploidy level and tolerance to environmental factors. Polyploids and species with increased tolerance to elevated temperatures and rainfall values exhibited high invasive capacity. We found no evidence that GS affects the invasive capacity of plants. We suggest that the genetic variability provided by polyploidization has a positive impact on plant competitiveness, which may ultimately lead to an increased ability to colonize new environments. In a global warming scenario, integrative approaches using phenotypic, genetic, epigenetic, and ecological traits should be a productive route to unveil the aspects of invasive plants.

Gene Flow in Volant Vertebrates: Species Biology, Ecology and Climate Change

Gene flow, the exchange of genetic material between populations is an important biological process, which shapes and maintains biodiversity. The successful movement of individuals between populations depends on multiple factors determined by species biology and the environment. One of the most important factors regulating gene flow is the ability to move, and flight allows individuals to easily move across geographical barriers. Volant vertebrates are found on some of the remotest islands and contribute significantly to the biodiversity and ecosystem. The availability of next-generation sequencing data for non-model animals has substantially improved our understanding of gene flow and its consequences, allowing us to look at fine-scale patterns. However, most of our understanding regarding gene flow comes from the temperate regions and the Neotropics. The lack of studies from species-rich Asia is striking. In this review, we outline the importance of gene flow and the factors affecting gene flow, especially for volant vertebrates. We especially discuss research studies from tropical biomes of South and Southeast Asia, highlight the lacuna in literature and provide an outline for future studies in this species-rich region.

A VKORC1-based SNP survey of anticoagulant rodenticide resistance in the house mouse, Norway rat and roof rat in the USA

BACKGROUND

We conducted a vitamin K epoxide reductase subcomponent 1 (Vkorc1)-based nonsynonymous Single Nucleotide Polymorphism (nsSNP) screen with focus on the house mouse (Mus musculus domesticus), but that also considered the Norway rat (Rattus norvegicus) and roof rat (R. rattus) in the USA.

RESULTS

We detected six Vkorc1 nsSNPs underlying the amino-acid polymorphisms Ala21Thr, Trp59Leu, Ile104Val, Val118Leu, Leu128Ser and Tyr139Cys in house mice (average coverage/SNP; n = 182 individuals), two nsSNPs underlying Arg35Pro and Gly46Ser in the Norway rat (n = 93), with the notable absence of Tyr139Cys (n = 179), and one nsSNP underlying Tyr25Phe in the roof rat (n = 27). Inferred resistance frequency is 29.1% for mice (variability of states 0-98.8%), 6.5% (0-33.3%) for the Norway rat, and 39.3% (0-52.6%) for the roof rat based on Tyr25Phe frequencies.

CONCLUSIONS

Resistance detected in the USA in the 1980s likely was the consequence of Vkorc1 mutations in mice (Leu128Ser and Tyr139Cys), Norway rats (Arg35Pro) and roof rats (Tyr25Phe). Patterns of variant sharing between the USA and Europe indicate the importance of convergent evolution and gene flow in spreading resistance. The spread of nsSNPs in mice between continents appears to have been more effective than in Norway rats. We hypothesize that Arg35Pro may have originated in Norway rats in the USA, whereas Tyr139Cys variants originated in Europe. Tyr25Phe is the likely cause for resistance in roof rats. Further genetic testing in the USA is required to close sampling gaps, and population genomic data are needed to study the origin and spread of this adaptive trait.

The complex genomic basis of rapid convergent adaptation to pesticides across continents in a fungal plant pathogen

Convergent evolution leads to identical phenotypic traits in different species or populations. Convergence can be driven by standing variation allowing selection to favour identical alleles in parallel or the same mutations can arise independently. However, the molecular basis of such convergent adaptation remains often poorly resolved. Pesticide resistance in agricultural ecosystems is a hallmark of convergence in phenotypic traits. Here, we analyse the major fungal pathogen Zymoseptoria tritici causing serious losses on wheat and with fungicide resistance emergence across several continents. We sampled three population pairs each from a different continent spanning periods early and late in the application of fungicides. To identify causal loci for resistance, we combined knowledge from molecular genetics work and performed genome-wide association studies (GWAS) on a global set of isolates. We discovered yet unknown factors in azole resistance including a gene encoding membrane associated functions. We found strong support for the "hotspot" model of resistance evolution with convergent changes in a small set of loci but additional loci showed more population-specific allele frequency changes. Genome-wide scans of selection showed that half of all known resistance loci were overlapping a selective sweep region. Hence, the application of fungicides was one of the major selective agents acting on the pathogen over the past decades. Furthermore, loci identified through GWAS showed the highest overlap with selective sweep regions underlining the importance to map phenotypic trait variation in evolving populations. Our population genomic analyses highlighted that both de novo mutations and gene flow contributed to convergent pesticide adaptation.

Recessivity of pyrethroid resistance and limited interspecies hybridization across Hyalella clades supports rapid and independent origins of resistance

Several populations of the amphipod, Hyalella azteca, have developed resistance to pyrethroid insecticides due to non-target exposure, but the dominance of the resistance trait is unknown. The current study investigated the dominance level of point mutations in natural populations of insecticide-resistant H. azteca and determined whether H. azteca from different clades with and without resistant alleles can hybridize and produce viable offspring. A parent generation (P₀) of non-resistant homozygous wild type H. azteca was crossbred with pyrethroid-resistant homozygous mutant animals and the tolerance of the filial 1 (F₁) generation to the pyrethroid insecticide, permethrin, was measured. Then the genotypes of the F₁ generation was examined to assure heterozygosity. The resistant parents had permethrin LC₅₀ values that ranged from 52 to 82 times higher than the non-resistant animals and both crossbreeding experiments produced heterozygous hybrid offspring that had LC₅₀ values similar to the non-resistant H. azteca parent. Dominance levels calculated for each of the crosses showed values close to 0, confirming that the L925I and L925V mutations were completely recessive. The lack of reproduction by hybrids of the C x D breeding confirmed that these clades are reproductively isolated and therefore introgression of adaptive alleles across these clades is unlikely. Potential evolutionary consequences of this selection include development of population bottlenecks, which may arise leading to fitness costs and reduced genetic diversity of H. azteca.

Deep Convergence, Shared Ancestry, and Evolutionary Novelty in the Genetic Architecture of Heliconius Mimicry

Convergent evolution can occur through different genetic mechanisms in different species. It is now clear that convergence at the genetic level is also widespread, and can be caused by either (i) parallel genetic evolution, where independently evolved convergent mutations arise in different populations or species, or (ii) collateral evolution in which shared ancestry results from either ancestral polymorphism or introgression among taxa. The adaptive radiation of Heliconius butterflies shows color pattern variation within species, as well as mimetic convergence between species. Using comparisons from across multiple hybrid zones, we use signals of shared ancestry to identify and refine multiple putative regulatory elements in Heliconius melpomene and its comimics, Heliconius elevatus and Heliconius besckei, around three known major color patterning genes: optix, WntA, and cortex While we find that convergence between H. melpomene and H. elevatus is caused by a complex history of collateral evolution via introgression in the Amazon, convergence between these species in the Guianas appears to have evolved independently. Thus, we find adaptive convergent genetic evolution to be a key driver of regulatory changes that lead to rapid phenotypic changes. Furthermore, we uncover evidence of parallel genetic evolution at some loci around optix and WntA in H. melpomene and its distant comimic Heliconius erato Ultimately, we show that all three of convergence, conservation, and novelty underlie the modular architecture of Heliconius color pattern mimicry.

Mus musculus populations in Western Australia lack VKORC1 mutations conferring resistance to first generation anticoagulant rodenticides: Implications for conservation and biosecurity

Background

Humans routinely attempt to manage pest rodent populations with anticoagulant rodenticides (ARs). We require information on resistance to ARs within rodent populations to have effective eradication programs that minimise exposure in non-target species. Mutations to the VKORC1 gene have been shown to confer resistance in rodents with high proportions of resistance in mice found in all European populations tested. We screened mutations in Mus musculus within Western Australia, by sampling populations from the capital city (Perth) and a remote island (Browse Island). These are the first Australian mouse populations screened for resistance using this method. Additionally, the mitochondrial D-loop of house mice was sequenced to explore population genetic structure, identify the origin of Western Australian mice, and to elucidate whether resistance was linked to certain haplotypes.

Results

No resistance-related VKORC1 mutations were detected in either house mouse population. A genetic introgression in the intronic sequence of the VKORC1 gene of Browse Island house mouse was detected which is thought to have originated through hybridisation with the Algerian mouse (Mus spretus). Analysis of the mitochondrial D-loop reported two haplotypes in the house mouse population of Perth, and two haplotypes in the population of Browse Island.

Conclusions

Both house mouse populations exhibited no genetic resistance to ARs, in spite of free use of ARs in Western Australia. Therefore weaker anticoagulant rodenticides can be employed in pest control and eradication attempts, which will result in reduced negative impacts on non-target species. Biosecurity measures must be in place to avoid introduction of resistant house mice, and new house mouse subspecies to Western Australia.

Further resolution of the house mouse (Mus musculus) phylogeny by integration over isolation-with-migration histories

BACKGROUND

The three main subspecies of house mice, Mus musculus castaneus, Mus musculus domesticus, and Mus musculus musculus, are estimated to have diverged ~ 350-500KYA. Resolution of the details of their evolutionary history is complicated by their relatively recent divergence, ongoing gene flow among the subspecies, and complex demographic histories. Previous studies have been limited to some extent by the number of loci surveyed and/or by the scope of the method used. Here, we apply a method (IMa3) that provides an estimate of a population phylogeny while allowing for complex histories of gene exchange.

RESULTS

Results strongly support a topology with M. m. domesticus as sister to M. m. castaneus and M. m. musculus. In addition, we find evidence of gene flow between all pairs of subspecies, but that gene flow is most restricted from M. m. musculus into M. m. domesticus. Estimates of other key parameters are dependent on assumptions regarding generation time and mutation rate in house mice. Nevertheless, our results support previous findings that the effective population size, Ne, of M. m. castaneus is larger than that of the other two subspecies, that the three subspecies began diverging ~ 130 - 420KYA, and that the time between divergence events was short.

CONCLUSIONS

Joint demographic and phylogenetic analyses of genomic data provide a clearer picture of the history of divergence in house mice.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

Patterns and impacts of nonvertical evolution in eukaryotes: a paradigm shift

Evolution of eukaryotic species and their genomes has been traditionally understood as a vertical process in which genetic material is transmitted from parents to offspring along a lineage, and in which genetic exchange is restricted within species boundaries. However, mounting evidence from comparative genomics indicates that this paradigm is often violated. Horizontal gene transfer and mating between diverged lineages blur species boundaries and challenge the reconstruction of evolutionary histories of species and their genomes. Nonvertical evolution might be more restricted in eukaryotes than in prokaryotes, yet it is not negligible and can be common in certain groups. Recognition of such processes brings about the need to incorporate this complexity into our models, as well as to conceptually reframe eukaryotic diversity and evolution. Here, I review the recent work from genomics studies that supports the effects of nonvertical modes of evolution including introgression, hybridization, and horizontal gene transfer in different eukaryotic groups. I then discuss emerging patterns and effects, illustrated by specific examples, that support the conclusion that nonvertical processes are often at the root of important evolutionary transitions and adaptations. I will argue that a paradigm shift is needed to naturally accommodate nonvertical processes in eukaryotic evolution.