Mosquito genomics. Extensive introgression in a malaria vector species complex revealed by phylogenomics

Ultra-conserved sequences in the genomes of highly diverse Anopheles mosquitoes, with implications for malaria vector control

DNA sequences that are exactly conserved over long evolutionary time scales have been observed in a variety of taxa. Such sequences are likely under strong functional constraint and they have been useful in the field of comparative genomics for identifying genome regions with regulatory function. A potential new application for these ultra-conserved elements (UCEs) has emerged in the development of gene drives to control mosquito populations. Many gene drives work by recognizing and inserting at a specific target sequence in the genome, often imposing a reproductive load as a consequence. They can therefore select for target sequence variants that provide resistance to the drive. Focusing on highly conserved, highly constrained sequences lowers the probability that variant, gene drive-resistant alleles can be tolerated. Here, we search for conserved sequences of 18 bp and over in an alignment of 21 Anopheles genomes, spanning an evolutionary timescale of 100 million years, and characterize the resulting sequences according to their location and function. Over 8000 UCEs were found across the alignment, with a maximum length of 164 bp. Length-corrected gene ontology analysis revealed that genes containing Anopheles UCEs were over-represented in categories with structural or nucleotide-binding functions. Known insect transcription factor binding sites were found in 48% of intergenic Anopheles UCEs. When we looked at the genome sequences of 1142 wild-caught mosquitoes, we found that 15% of the Anopheles UCEs contained no polymorphisms. Our list of Anopheles UCEs should provide a valuable starting point for the selection and testing of new targets for gene-drive modification in the mosquitoes that transmit malaria.

Chromosome-level genome assemblies of the malaria vectors Anopheles coluzzii and Anopheles arabiensis

BACKGROUND

Anopheles coluzzii and Anopheles arabiensis belong to the Anopheles gambiae complex and are among the major malaria vectors in sub-Saharan Africa. However, chromosome-level reference genome assemblies are still lacking for these medically important mosquito species.

FINDINGS

In this study, we produced de novo chromosome-level genome assemblies for A. coluzzii and A. arabiensis using the long-read Oxford Nanopore sequencing technology and the Hi-C scaffolding approach. We obtained 273.4 and 256.8 Mb of the total assemblies for A. coluzzii and A. arabiensis, respectively. Each assembly consists of 3 chromosome-scale scaffolds (X, 2, 3), complete mitochondrion, and unordered contigs identified as autosomal pericentromeric DNA, X pericentromeric DNA, and Y sequences. Comparison of these assemblies with the existing assemblies for these species demonstrated that we obtained improved reference-quality genomes. The new assemblies allowed us to identify genomic coordinates for the breakpoint regions of fixed and polymorphic chromosomal inversions in A. coluzzii and A. arabiensis.

CONCLUSION

The new chromosome-level assemblies will facilitate functional and population genomic studies in A. coluzzii and A. arabiensis. The presented assembly pipeline will accelerate progress toward creating high-quality genome references for other disease vectors.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Total Ortholog Median Matrix as an alternative unsupervised approach for phylogenomics based on evolutionary distance between protein coding genes

The increasing number of available genomic data allowed the development of phylogenomic analytical tools. Current methods compile information from single gene phylogenies, whether based on topologies or multiple sequence alignments. Generally, phylogenomic analyses elect gene families or genomic regions to construct phylogenomic trees. Here, we presented an alternative approach for Phylogenomics, named TOMM (Total Ortholog Median Matrix), to construct a representative phylogram composed by amino acid distance measures of all pairwise ortholog protein sequence pairs from desired species inside a group of organisms. The procedure is divided two main steps, (1) ortholog detection and (2) creation of a matrix with the median amino acid distance measures of all pairwise orthologous sequences. We tested this approach within three different group of organisms: Kinetoplastida protozoa, hematophagous Diptera vectors and Primates. Our approach was robust and efficacious to reconstruct the phylogenetic relationships for the three groups. Moreover, novel branch topologies could be achieved, providing insights about some phylogenetic relationships between some taxa.

Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe)

Inferences about past processes of adaptation and speciation require a gene-scale and genome-wide understanding of the evolutionary history of diverging taxa. In this study, we use genome-wide capture of nuclear gene sequences, plus skimming of organellar sequences, to investigate the phylogenomics of monkeyflowers in Mimulus section Erythranthe (27 accessions from seven species). Taxa within Erythranthe, particularly the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated), have been a model system for investigating the ecological genetics of speciation and adaptation for over five decades. Across >8000 nuclear loci, multiple methods resolve a predominant species tree in which M. cardinalis groups with other hummingbird-pollinated taxa (37% of gene trees), rather than being sister to M. lewisii (32% of gene trees). We independently corroborate a single evolution of hummingbird pollination syndrome in Erythranthe by demonstrating functional redundancy in genetic complementation tests of floral traits in hybrids; together, these analyses overturn a textbook case of pollination-syndrome convergence. Strong asymmetries in allele sharing (Patterson's D-statistic and related tests) indicate that gene tree discordance reflects ancient and recent introgression rather than incomplete lineage sorting. Consistent with abundant introgression blurring the history of divergence, low-recombination and adaptation-associated regions support the new species tree, while high-recombination regions generate phylogenetic evidence for sister status for M. lewisii and M. cardinalis. Population-level sampling of core taxa also revealed two instances of chloroplast capture, with Sierran M. lewisii and Southern Californian M. parishii each carrying organelle genomes nested within respective sympatric M. cardinalis clades. A recent organellar transfer from M. cardinalis, an outcrosser where selfish cytonuclear dynamics are more likely, may account for the unexpected cytoplasmic male sterility effects of selfer M. parishii organelles in hybrids with M. lewisii. Overall, our phylogenomic results reveal extensive reticulation throughout the evolutionary history of a classic monkeyflower radiation, suggesting that natural selection (re-)assembles and maintains species-diagnostic traits and barriers in the face of gene flow. Our findings further underline the challenges, even in reproductively isolated species, in distinguishing re-use of adaptive alleles from true convergence and emphasize the value of a phylogenomic framework for reconstructing the evolutionary genetics of adaptation and speciation.

Dsuite - Fast D-statistics and related admixture evidence from VCF files

Patterson's D, also known as the ABBA-BABA statistic, and related statistics such as the f4 -ratio, are commonly used to assess evidence of gene flow between populations or closely related species. Currently available implementations often require custom file formats, implement only small subsets of the available statistics, and are impractical to evaluate all gene flow hypotheses across data sets with many populations or species due to computational inefficiencies. Here, we present a new software package Dsuite, an efficient implementation allowing genome scale calculations of the D and f4 -ratio statistics across all combinations of tens or hundreds of populations or species directly from a variant call format (VCF) file. Our program also implements statistics suited for application to genomic windows, providing evidence of whether introgression is confined to specific loci, and it can also aid in interpretation of a system of f4 -ratio results with the use of the "f-branch" method. Dsuite is available at https://github.com/millanek/Dsuite, is straightforward to use, substantially more computationally efficient than comparable programs, and provides a convenient suite of tools and statistics, including some not previously available in any software package. Thus, Dsuite facilitates the assessment of evidence for gene flow, especially across larger genomic data sets.

The Gene and Gene Expression (GAGE) Species Concept: An Universal Approach for All Eukaryotic Organisms

The Gene and Gene Expression (GAGE) species concept, a new version of the Pragmatic Species Concept of Seifert (2014), is proposed as a concept applicable to any described recent or fossil eukaryotic organism independent from its mode of reproduction or evolutionary history. In addition to presenting the concept as such, the article also provides practical recommendations for taxonomists when delimiting species and describing taxa. The wording of the new concept contains a heading core sentence plus five attached sentences addressing essential conditions for its translation into a sound taxonomic practice: "Species are separable clusters that have passed a threshold of evolutionary divergence and are exclusively defined by nuclear DNA sequences and/or their expression products. Nuclear DNA sequences and their expression products are different character systems but have a highly correlated indicative function. Character systems with the least risk of epigenetic or ontogenetic modification have superior indicative value when conflicts between character systems of integrative studies arise. All character systems have to be described by an adequate numerics allowing cluster formation and determination of thresholds. Thresholds for each character system should be fixed by consensus among the experts under the principle of avoiding oversplitting or lumping. Clusters must not be the expression of intraspecific polymorphism." Recognizing the distortions and conflicts caused to taxonomy through barcoding or through assessment on the basis of association with other organisms, the GAGE species concept strongly downgrades the use of cytoplasmic DNA of endosymbiotic origin (mtDNA, cpDNA) or DNA of closely associated microbes (e.g., Wolbachia bacteria) for final taxonomic decision-making. Recognizing the distortion of phylogenies by the high frequency of reticulate evolution, it is argued that delimiting and naming species has to be separated from constructing bifurcating phylogenetic trees. [Cytoplasmic DNA; lumping; nuclear DNA; numeric taxonomy; oversplitting; reticulate evolution.].

Resistance to pirimiphos-methyl in West African Anopheles is spreading via duplication and introgression of the Ace1 locus

Vector population control using insecticides is a key element of current strategies to prevent malaria transmission in Africa. The introduction of effective insecticides, such as the organophosphate pirimiphos-methyl, is essential to overcome the recurrent emergence of resistance driven by the highly diverse Anopheles genomes. Here, we use a population genomic approach to investigate the basis of pirimiphos-methyl resistance in the major malaria vectors Anopheles gambiae and A. coluzzii. A combination of copy number variation and a single non-synonymous substitution in the acetylcholinesterase gene, Ace1, provides the key resistance diagnostic in an A. coluzzii population from Côte d'Ivoire that we used for sequence-based association mapping, with replication in other West African populations. The Ace1 substitution and duplications occur on a unique resistance haplotype that evolved in A. gambiae and introgressed into A. coluzzii, and is now common in West Africa primarily due to selection imposed by other organophosphate or carbamate insecticides. Our findings highlight the predictive value of this complex resistance haplotype for phenotypic resistance and clarify its evolutionary history, providing tools to for molecular surveillance of the current and future effectiveness of pirimiphos-methyl based interventions.

Determining the probability of hemiplasy in the presence of incomplete lineage sorting and introgression

The incongruence of character states with phylogenetic relationships is often interpreted as evidence of convergent evolution. However, trait evolution along discordant gene trees can also generate these incongruences - a phenomenon known as hemiplasy. Classic comparative methods do not account for discordance, resulting in incorrect inferences about the number, timing, and direction of trait transitions. Biological sources of discordance include incomplete lineage sorting (ILS) and introgression, but only ILS has received theoretical consideration in the context of hemiplasy. Here, we present a model that shows introgression makes hemiplasy more likely, such that methods that account for ILS alone will be conservative. We also present a method and software (HeIST) for making statistical inferences about the probability of hemiplasy and homoplasy in large datasets that contain both ILS and introgression. We apply our methods to two empirical datasets, finding that hemiplasy is likely to contribute to the observed trait incongruences in both.

Maximum Likelihood Estimation of Species Trees from Gene Trees in the Presence of Ancestral Population Structure

Though large multilocus genomic data sets have led to overall improvements in phylogenetic inference, they have posed the new challenge of addressing conflicting signals across the genome. In particular, ancestral population structure, which has been uncovered in a number of diverse species, can skew gene tree frequencies, thereby hindering the performance of species tree estimators. Here we develop a novel maximum likelihood method, termed TASTI (Taxa with Ancestral structure Species Tree Inference), that can infer phylogenies under such scenarios, and find that it has increasing accuracy with increasing numbers of input gene trees, contrasting with the relatively poor performances of methods not tailored for ancestral structure. Moreover, we propose a supertree approach that allows TASTI to scale computationally with increasing numbers of input taxa. We use genetic simulations to assess TASTI's performance in the three- and four-taxon settings and demonstrate the application of TASTI on a six-species Afrotropical mosquito data set. Finally, we have implemented TASTI in an open-source software package for ease of use by the scientific community.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.

Radiation with reticulation marks the origin of a major malaria vector

Advances in genomics have led to an appreciation that introgression is common, but its evolutionary consequences are poorly understood. In recent species radiations the sharing of genetic variation across porous species boundaries can facilitate adaptation to new environments and generate novel phenotypes, which may contribute to further diversification. Most Anopheles mosquito species that are of major importance as human malaria vectors have evolved within recent and rapid radiations of largely nonvector species. Here, we focus on one of the most medically important yet understudied anopheline radiations, the Afrotropical Anopheles funestus complex (AFC), to investigate the role of introgression in its diversification and the possible link between introgression and vector potential. The AFC comprises at least seven morphologically similar species, yet only An. funestus sensu stricto is a highly efficient malaria vector with a pan-African distribution. Based on de novo genome assemblies and additional whole-genome resequencing, we use phylogenomic and population genomic analyses to establish species relationships. We show that extensive interspecific gene flow involving multiple species pairs has shaped the evolutionary history of the AFC since its diversification. The most recent introgression event involved a massive and asymmetrical movement of genes from a distantly related AFC lineage into An. funestus, an event that predated and plausibly facilitated its subsequent dramatic geographic range expansion across most of tropical Africa. We propose that introgression may be a common mechanism facilitating adaptation to new environments and enhancing vectorial capacity in Anopheles mosquitoes.

Unveiling mosquito cryptic species and their reproductive isolation

Mosquitoes are major vectors of many infectious pathogens or parasites. Understanding cryptic species and the speciation of disease vectors has important implications for vector management, evolution and host-pathogen and/or host-parasite interactions. Currently, mosquito cryptic species have been reported in many studies, most of which focus on the reproductive isolation of cryptic species and mainly on Anopheles gambiae sensu lato complex. Emerging species within the primary malaria vector Anopheles gambiae show different ecological preferences and significant prezygotic reproductive isolation, while Aedes mariae and Aedes zammitii show postmating reproductive isolation. However, data reporting the reproductive isolation in Culex and Aedes albopictus mosquito cryptic species is absent. The lack of systematic studies leaves many questions open, such as whether cryptic species are more common in particular habitats, latitudes or taxonomic groups; what mosquito cryptic species evolutionary processes bring about reproductive isolation in the absence of morphological differentiation? How does Wolbachia infection affect in mosquitoes' reproductive isolation? In this review, we provide a summary of recent advances in the discovery and identification of sibling or cryptic species within mosquito genera.

Introgressive hybridization between two non-native apple snails in China: widespread hybridization and homogenization in egg morphology

BACKGROUND

Apple snails from the genus Pomacea have spread widely in paddy fields and other wetlands of southern China since their introduction in the 1980s. Pomacea spp. are commonly identified using mitochondrial COI sequences. However, sequencing the nuclear elongation factor 1-alpha (EF1α) gene revealed genetic introgression between field populations of P. canaliculata and P. maculata, which produce surviving hybrids in laboratory crossbreeding experiments.

RESULTS

In this study, we sequenced 1054 EF1α clones to design specific primers and established a fast and accurate multiplex polymerase chain reaction (PCR) method for genotyping EF1α. Combined with genotyping P. canaliculata and P. maculata based on mitochondrial COI and nuclear EF1α, we revealed the genetic introgression patterns of 30 apple snail populations in China. Purebred and hybrid individuals of P. canaliculata were widely distributed, while pure maculata-EF1α type was detected only in a few individuals identified as P. canaliculata based on COI sequences. Each egg clutch had one to three genetic patterns, indicating multiple paternity or segregation in the progeny of hybrids. The higher percentages of hybrids in both wild populations and progeny than the homozygotes indicated a potential heterosis in the apple snail populations. Additionally, egg size and clutch size of the apple snails became homogeneous among the non-native populations exhibiting introgression hybridization.

CONCLUSION

Our findings emphasize the value of apple snails as a model to study the mechanisms and impacts of introgressive hybridization on fitness traits. © 2020 Society of Chemical Industry.

Impact of homologous recombination on core genome phylogenies

BACKGROUND

Core genome phylogenies are widely used to build the evolutionary history of individual prokaryote species. By using hundreds or thousands of shared genes, these approaches are the gold standard to reconstruct the relationships of large sets of strains. However, there is growing evidence that bacterial strains exchange DNA through homologous recombination at rates that vary widely across prokaryote species, indicating that core genome phylogenies might not be able to reconstruct true phylogenies when recombination rate is high. Few attempts have been made to evaluate the robustness of core genome phylogenies to recombination, but some analyses suggest that reconstructed trees are not always accurate.

RESULTS

In this study, we tested the robustness of core genome phylogenies to various levels of recombination rates. By analyzing simulated and empirical data, we observed that core genome phylogenies are relatively robust to recombination rates; nevertheless, our results suggest that many reconstructed trees are not completely accurate even when bootstrap supports are high. We found that some core genome phylogenies are highly robust to recombination whereas others are strongly impacted by it, and we identified that the robustness of core genome phylogenies to recombination is highly linked to the levels of selective pressures acting on a species. Stronger selective pressures lead to less accurate tree reconstructions, presumably because selective pressures more strongly bias the routes of DNA transfers, thereby causing phylogenetic artifacts.

CONCLUSIONS

Overall, these results have important implications for the application of core genome phylogenies in prokaryotes.

Phylogenomics and the Genetic Architecture of the Placental Mammal Radiation

The genomes of placental mammals are being sequenced at an unprecedented rate. Alignments of hundreds, and one day thousands, of genomes spanning the rich living and extinct diversity of species offer unparalleled power to resolve phylogenetic controversies, identify genomic innovations of adaptation, and dissect the genetic architecture of reproductive isolation. We highlight outstanding questions about the earliest phases of placental mammal diversification and the promise of newer methods, as well as remaining challenges, toward using whole genome data to resolve placental mammal phylogeny. The next phase of mammalian comparative genomics will see the completion and application of finished-quality, gapless genome assemblies from many ordinal lineages and closely related species. Interspecific comparisons between the most hypervariable genomic loci will likely reveal large, but heretofore mostly underappreciated, effects on population divergence, morphological innovation, and the origin of new species.

A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly

BACKGROUND

The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics.

RESULTS

Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS).

CONCLUSIONS

Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

Culicidae evolutionary history focusing on the Culicinae subfamily based on mitochondrial phylogenomics

Mosquitoes are insects of medical importance due their role as vectors of different pathogens to humans. There is a lack of information about the evolutionary history and phylogenetic positioning of the majority of mosquito species. Here we characterized the mitogenomes of mosquito species through low-coverage whole genome sequencing and data mining. A total of 37 draft mitogenomes of different species were assembled from which 16 are newly-sequenced species. We datamined additional 49 mosquito mitogenomes, and together with our 37 mitogenomes, we reconstructed the evolutionary history of 86 species including representatives from 15 genera and 7 tribes. Our results showed that most of the species clustered in clades with other members of their own genus with exception of Aedes genus which was paraphyletic. We confirmed the monophyletic status of the Mansoniini tribe including both Coquillettidia and Mansonia genus. The Aedeomyiini and Uranotaeniini were consistently recovered as basal to other tribes in the subfamily Culicinae, although the exact relationships among these tribes differed between analyses. These results demonstrate that low-coverage sequencing is effective to recover mitogenomes, establish phylogenetic knowledge and hence generate basic fundamental information that will help in the understanding of the role of these species as pathogen vectors.

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.

EFSA Panel on Genetically Modified Organisms (GMO), Naegeli H, Bresson J, Dalmay T, Dewhurst IC, Epstein MM, Guerche P, Hejatko J, Moreno FJ, Mullins E, Nogué F, Rostoks N, Sánchez Serrano JJ, Savoini G, Veromann E, Veronesi F, Bonsall MB, Mumford J, Wimmer EA, Devos Y, Paraskevopoulos K, Firbank LG. Adequacy and sufficiency evaluation of existing EFSA guidelines for the molecular characterisation, environmental risk assessment and post-market environmental monitoring of genetically modified insects containing engineered gene drives

Advances in molecular and synthetic biology are enabling the engineering of gene drives in insects for disease vector/pest control. Engineered gene drives (that bias their own inheritance) can be designed either to suppress interbreeding target populations or modify them with a new genotype. Depending on the engineered gene drive system, theoretically, a genetic modification of interest could spread through target populations and persist indefinitely, or be restricted in its spread or persistence. While research on engineered gene drives and their applications in insects is advancing at a fast pace, it will take several years for technological developments to move to practical applications for deliberate release into the environment. Some gene drive modified insects (GDMIs) have been tested experimentally in the laboratory, but none has been assessed in small-scale confined field trials or in open release trials as yet. There is concern that the deliberate release of GDMIs in the environment may have possible irreversible and unintended consequences. As a proactive measure, the European Food Safety Authority (EFSA) has been requested by the European Commission to review whether its previously published guidelines for the risk assessment of genetically modified animals (EFSA, 2012 and 2013), including insects (GMIs), are adequate and sufficient for GDMIs, primarily disease vectors, agricultural pests and invasive species, for deliberate release into the environment. Under this mandate, EFSA was not requested to develop risk assessment guidelines for GDMIs. In this Scientific Opinion, the Panel on Genetically Modified Organisms (GMO) concludes that EFSA's guidelines are adequate, but insufficient for the molecular characterisation (MC), environmental risk assessment (ERA) and post-market environmental monitoring (PMEM) of GDMIs. While the MC,ERA and PMEM of GDMIs can build on the existing risk assessment framework for GMIs that do not contain engineered gene drives, there are specific areas where further guidance is needed for GDMIs.

Assessment of the ecologically dependent post-zygotic isolation between Anopheles coluzzii and Anopheles gambiae

Within the Anopheles gambiae complex, the sibling species An. coluzzii and An. gambiae are undergoing sympatric speciation. These species are characterized by rare hybrids in most of their geographical distribution. A strong assortative mating mediated by spatial swarm segregation has been shown whereas no intrinsic post-zygotic barriers have been found in laboratory conditions. To test the role of the hybridisation in reproductive isolation in natural populations transplant experiment are therefore needed to establish the significance of post-zygotic barriers. Previous studies indicated that predation is one of the major forces driving ecological divergence between An. gambiae and An. coluzzii. Here we extended these studies to their hybrids. Parental species and their F1 hybrids from reciprocal crosses were generated by the forced-mating technique as follows: female An. coluzzii × male An. coluzzii; female An. coluzzii × male An. gambiae; female An. gambiae × male An. coluzzii and female An. gambiae × Male An. gambiae. First instar larvae of each group from the crossing (here after An. coluzzii, Hybrid COL/GAM, Hybrid GAM/COL and An. gambiae, respectively) were transplanted in a field experiment with predation effect. Emergence success, development time of larvae and body size of the newly emerging adults were estimated as fitness components and then compared between parental species and F1 hybrids in absence and in presence of predators. Our findings confirm that An. coluzzii had higher fitness than An. gambiae in presence of predators versus in absence of predators. Moreover, the fitness of the F1 hybrid COL/GAM whose female parent was An. coluzzii matched that of An. coluzzii while that of the F1 reciprocal hybrid GAM/COL was similar to An. gambiae.

Detecting inversions with PCA in the presence of population structure

Chromosomal inversions can lead to reproductive isolation and adaptation in insects such as Drosophila melanogaster and the non-model malaria vector Anopheles gambiae. Inversions can be detected and characterized using principal component analysis (PCA) of single nucleotide polymorphisms (SNPs). To aid in developing such methods, we formed a new benchmark derived from three publicly-available insect data. We then used this benchmark to perform an extended validation of our software for inversion analysis (Asaph). Through that process, we identified and characterized several problematic test cases liable to misinterpretation that can help guide PCA-based inversion detection. Lastly, we re-analyzed the 2R chromosome arm of 150 An. gambiae and coluzzii samples and observed two inversions (2Rc and 2Rd) that were previously known but not annotated in these particular individuals. The resulting benchmark data set and methods will be useful for future inversion detection based solely on SNP data.

Evolution of the Insecticide Target Rdl in African Anopheles Is Driven by Interspecific and Interkaryotypic Introgression

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across nonconcordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors.

Speciation in North American Junonia from a genomic perspective

Delineating species boundaries in phylogenetic groups undergoing recent radiation is a daunting challenge akin to discretizing continuity. Here, we propose a general approach exemplified by American butterflies from the genus Junonia Hübner notorious for the variety of similar phenotypes, ease of hybridization, and the lack of consensus about their classification. We obtain whole-genome shotgun sequences of about 200 specimens. We reason that discreteness emerges from continuity by means of a small number of key players, and search for the proteins that diverged markedly between sympatric populations of different species, while keeping low polymorphism within these species. Being 0.25% of the total number, these three dozen 'speciation' proteins indeed partition pairs of Junonia populations into two clusters with a prominent break in between, while all proteins taken together fail to reveal this discontinuity. Populations with larger divergence from each other, comparable to that between two sympatric species, form the first cluster and correspond to different species. The other cluster is characterized by smaller divergence, similar to that between allopatric populations of the same species and comprise conspecific pairs. Using this method, we conclude that J. genoveva (Cramer), J. litoralis Brévignon, J. evarete (Cramer), and J. divaricata C. & R. Felder are restricted to South America. We find that six species of Junonia are present in the United States, one of which is new: Junonia stemosa Grishin, sp.n. (i), found in south Texas and phenotypically closest to J. nigrosuffusa W. Barnes & McDunnough (ii) in its dark appearance. In the pale nudum of the antennal club, these two species resemble J. zonalis C. & R. Felder (iii) from Florida and the Caribbean Islands. The pair of sister species, J. grisea Austin & J. Emmel (iv) and J. coenia Hübner (v), represent the classic west/east U.S.A. split. The mangrove feeder (as caterpillar), dark nudum J. neildi Brévignon (vi) enters south Texas as a new subspecies Junonia neildi varia Grishin ssp.n. characterized by more extensive hybridization with and introgression from J. coenia, and, as a consequence, more variable wing patterns compared with the nominal J. n. neildi in Florida. Furthermore, a new mangrove-feeding species from the Pacific Coast of Mexico is described as Junonia pacoma Grishin sp.n. Finally, genomic analysis suggests that J. nigrosuffusa may be a hybrid species formed by the ancestors of J. grisea and J. stemosa sp.n.

Climate and Urbanization Drive Mosquito Preference for Humans

The majority of mosquito-borne illness is spread by a few mosquito species that have evolved to specialize in biting humans, yet the precise causes of this behavioral shift are poorly understood. We address this gap in the arboviral vector Aedes aegypti. We first collect and characterize the behavior of mosquitoes from 27 sites scattered across the species' ancestral range in sub-Saharan Africa, revealing previously unrecognized variation in preference for human versus animal odor. We then use modeling to show that over 80% of this variation can be predicted by two ecological factors-dry season intensity and human population density. Finally, we integrate this information with whole-genome sequence data from 375 individual mosquitoes to identify a single underlying ancestry component linked to human preference. Genetic changes associated with human specialist ancestry were concentrated in a few chromosomal regions. Our findings suggest that human-biting in this important disease vector originally evolved as a by-product of breeding in human-stored water in areas where doing so provided the only means to survive the long, hot dry season. Our model also predicts that the rapid urbanization currently taking place in Africa will drive further mosquito evolution, causing a shift toward human-biting in many large cities by 2050.

Mitochondrial genomics reveals the evolutionary history of the porpoises (Phocoenidae) across the speciation continuum

Historical variation in food resources is expected to be a major driver of cetacean evolution, especially for the smallest species like porpoises. Despite major conservation issues among porpoise species (e.g., vaquita and finless), their evolutionary history remains understudied. Here, we reconstructed their evolutionary history across the speciation continuum. Phylogenetic analyses of 63 mitochondrial genomes suggest that porpoises radiated during the deep environmental changes of the Pliocene. However, all intra-specific subdivisions were shaped during the Quaternary glaciations. We observed analogous evolutionary patterns in both hemispheres associated with convergent evolution to coastal versus oceanic environments. This suggests that similar mechanisms are driving species diversification in northern (harbor and Dall's) and southern species (spectacled and Burmeister's). In contrast to previous studies, spectacled and Burmeister's porpoises shared a more recent common ancestor than with the vaquita that diverged from southern species during the Pliocene. The low genetic diversity observed in the vaquita carried signatures of a very low population size since the last 5,000 years. Cryptic lineages within Dall's, spectacled and Pacific harbor porpoises suggest a richer evolutionary history than previously suspected. These results provide a new perspective on the mechanisms driving diversification in porpoises and an evolutionary framework for their conservation.

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.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Biased Gene Retention in the Face of Introgression Obscures Species Relationships

Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.

The Genetic Architecture of Post-Zygotic Reproductive Isolation Between Anopheles coluzzii and An. quadriannulatus

The Anopheles gambiae complex is comprised of eight morphologically indistinguishable species and has emerged as a model system for the study of speciation genetics due to the rapid radiation of its member species over the past two million years. Male hybrids between most An. gambiae complex species pairs are sterile, and some genotype combinations in hybrid males cause inviability. We investigated the genetic basis of hybrid male inviability and sterility between An. coluzzii and An. quadriannulatus by measuring segregation distortion and performing a QTL analysis of sterility in a backcross population. Hybrid males were inviable if they inherited the An. coluzzii X chromosome and were homozygous at one or more loci in 18.9 Mb region of chromosome 3. The An. coluzzii X chromosome has a disproportionately large effect on hybrid sterility when introgressed into an An. quadriannulatus genetic background. Additionally, an epistatic interaction between the An. coluzzii X and a 1.12 Mb, pericentric region of the An. quadriannulatus 3L chromosome arm has a statistically significant contribution to the hybrid sterility phenotype. This same epistatic interaction occurs when the An. coluzzii X is introgressed into the genetic background of An. arabiensis, the sister species of An. quadriannulatus, suggesting that this may represent one of the first Dobzhansky-Muller incompatibilities to evolve early in the radiation of the Anopheles gambiae species complex. We describe the additive effects of each sterility QTL, epistatic interactions between them, and genes within QTL with protein functions related to mating behavior, reproduction, spermatogenesis, and microtubule morphogenesis, whose divergence may contribute to post-zygotic reproductive isolation between An. coluzzii and An. quadriannulatus.

Genes and speciation: is it time to abandon the biological species concept?

The biological species concept (BSC) is the cornerstone of neo-Darwinian thinking. In BSC, species do not exchange genes either during or after speciation. However, as gene flow during speciation is increasingly being reported in a substantial literature, it seems time to reassess the revered, but often doubted, BSC. Contrary to the common perception, BSC should expect substantial gene flow at the onset of speciation, not least because geographical isolation develops gradually. Although BSC does not stipulate how speciation begins, it does require a sustained period of isolation for speciation to complete its course. Evidence against BSC must demonstrate that the observed gene flow does not merely occur at the onset of speciation but continues until its completion. Importantly, recent genomic analyses cannot reject this more realistic version of BSC, although future analyses may still prove it wrong. The ultimate acceptance or rejection of BSC is not merely about a historical debate; rather, it is about the fundamental nature of species - are species (and, hence, divergent adaptations) driven by a relatively small number of genes, or by thousands of them? Many levels of biology, ranging from taxonomy to biodiversity, depend on this resolution.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Defining Species When There is Gene Flow

Whatever one's definition of species, it is generally expected that individuals of the same species should be genetically more similar to each other than they are to individuals of another species. Here, we show that in the presence of cross-species gene flow, this expectation may be incorrect. We use the multispecies coalescent model with continuous-time migration or episodic introgression to study the impact of gene flow on genetic differences within and between species and highlight a surprising but plausible scenario in which different population sizes and asymmetrical migration rates cause a genetic sequence to be on average more closely related to a sequence from another species than to a sequence from the same species. Our results highlight the extraordinary impact that even a small amount of gene flow may have on the genetic history of the species. We suggest that contrasting long-term migration rate and short-term hybridization rate, both of which can be estimated using genetic data, may be a powerful approach to detecting the presence of reproductive barriers and to define species boundaries.[Gene flow; introgression; migration; multispecies coalescent; species concept; species delimitation.].

Reticulate evolution as a management challenge: Patterns of admixture with phylogenetic distance in endemic fishes of western North America

Admixture in natural populations is a long-standing management challenge, with population genomic approaches offering means for adjudication. We now more clearly understand the permeability of species boundaries and the potential of admixture for promoting adaptive evolution. These issues particularly resonate in western North America, where tectonism and aridity have fragmented and reshuffled rivers over millennia, in turn promoting reticulation among endemic fishes, a situation compounded by anthropogenic habitat modifications and non-native introductions. The melding of historic and contemporary admixture has both confused and stymied management. We underscore this situation with a case study that quantifies basin-wide admixture among a group of native and introduced fishes by employing double-digest restriction site-associated DNA (ddRAD) sequencing. Our approach: (a) quantifies the admixed history of 343 suckers (10 species of Catostomidae) across the Colorado River Basin; (b) gauges admixture within the context of phylogenetic distance and "ecological specialization"; and (c) extrapolates potential drivers of introgression across hybrid crosses that involve endemic as well as invasive species. Our study extends across an entire freshwater basin and expands previous studies more limited in scope both geographically and taxonomically. Our results detected admixture involving all 10 species, with habitat alterations not only accelerating the breakdown of reproductive isolation, but also promoting introgression. Hybridization occurred across the genus despite phylogenetic distance, whereas introgression was only detected within subgenera, implicating phylogenetic distance and/or ecological specialization as drivers of reproductive isolation. Understanding the extent of admixture and reproductive isolation across multiple species serves to disentangle their reticulate evolutionary histories and provides a broadscale perspective for basin-wide conservation and management.

High-density genetic variation maps reveal the correlation between asymmetric interspecific introgressions and improvement of agronomic traits in Upland and Pima cotton varieties developed in Xinjiang, China

The two new world tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, are cultivated worldwide and are characterised by a high yield and superior fibre quality, respectively. Historical genetic introgression has been reported between them; however, the existence of introgression and its genetic effects on agronomic traits remain unclear with regard to independent breeding of G. hirsutum (Upland cotton) and G. barbadense (Pima cotton) elite cultivars. We collected 159 G. hirsutum and 70 G. barbadense cultivars developed in Xinjiang, China, along with 30 semi-wild accessions of G. hirsutum, to perform interspecific introgression tests, intraspecific selection analyses and genome-wide association studies (GWAS) with fibre quality and yield component traits in multiple environments. In total, we identified seven interspecific introgression events and 52 selective sweep loci in G. hirsutum, as well as 17 interspecific introgression events and 19 selective sweep loci in G. barbadense. Correlation tests between agronomic traits and introgressions showed that introgression loci were mutually beneficial for the improvement of fibre quality and yield traits in both species. In addition, the phenotypic effects of four interspecific introgression events could be detected by intraspecific GWAS, with Gb_INT13 significantly improving fibre yield in G. barbadense. The present study describes the landscape of genetic introgression and selection between the two species, and highlights the genetic effects of introgression among populations, which can be used for future improvement of fibre yield and quality in G. barbadense and G. hirsutum, respectively.

Gene flow and species delimitation in fishes of Western North America: Flannelmouth (Catostomus latipinnis) and Bluehead sucker (C. Pantosteus discobolus)

The delimitation of species boundaries, particularly those obscured by reticulation, is a critical step in contemporary biodiversity assessment. It is especially relevant for conservation and management of indigenous fishes in western North America, represented herein by two species with dissimilar life histories codistributed in the highly modified Colorado River (i.e., flannelmouth sucker, Catostomus latipinnis; bluehead sucker, C. (Pantosteus) discobolus). To quantify phylogenomic patterns and examine proposed taxonomic revisions, we first employed double-digest restriction site-associated DNA sequencing (ddRAD), yielding 39,755 unlinked SNPs across 139 samples. These were subsequently evaluated with multiple analytical approaches and by contrasting life history data. Three phylogenetic methods and a Bayesian assignment test highlighted similar phylogenomic patterns in each, but with considerable difference in presumed times of divergence. Three lineages were detected in bluehead sucker, supporting elevation of C. (P.) virescens to species status and recognizing C. (P.) discobolus yarrowi (Zuni bluehead sucker) as a discrete entity. Admixture in the latter necessitated a reevaluation of its contemporary and historic distributions, underscoring how biodiversity identification can be confounded by complex evolutionary histories. In addition, we defined three separate flannelmouth sucker lineages as ESUs (evolutionarily significant units), given limited phenotypic and genetic differentiation, contemporary isolation, and lack of concordance (per the genealogical concordance component of the phylogenetic species concept). Introgression was diagnosed in both species, with the Little Colorado and Virgin rivers in particular. Our diagnostic methods, and the agreement of our SNPs with previous morphological, enzymatic, and mitochondrial work, allowed us to partition complex evolutionary histories into requisite components, such as isolation versus secondary contact.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.]

A mitochondrial genetic divergence proxy predicts the reproductive compatibility of mammalian hybrids

Numerous pairs of evolutionarily divergent mammalian species have been shown to produce hybrid offspring. In some cases, F1 hybrids are able to produce F2s through matings with F1s. In other instances, the hybrids are only able to produce offspring themselves through backcrosses with a parent species owing to unisexual sterility (Haldane's Rule). Here, we explicitly tested whether genetic distance, computed from mitochondrial and nuclear genes, can be used as a proxy to predict the relative fertility of the hybrid offspring resulting from matings between species of terrestrial mammals. We assessed the proxy's predictive power using a well-characterized felid hybrid system, and applied it to modern and ancient hominins. Our results revealed a small overlap in mitochondrial genetic distance values that distinguish species pairs whose calculated distances fall within two categories: those whose hybrid offspring follow Haldane's Rule, and those whose hybrid F1 offspring can produce F2s. The strong correlation between genetic distance and hybrid fertility demonstrated here suggests that this proxy can be employed to predict whether the hybrid offspring of two mammalian species will follow Haldane's Rule.

Nuclear phylogenomic analyses of asterids conflict with plastome trees and support novel relationships among major lineages

PREMISE

Discordance between nuclear and organellar phylogenies (cytonuclear discordance) is a well-documented phenomenon at shallow evolutionary levels but has been poorly investigated at deep levels of plant phylogeny. Determining the extent of cytonuclear discordance across major plant lineages is essential not only for elucidating evolutionary processes, but also for evaluating the currently used framework of plant phylogeny, which is largely based on the plastid genome.

METHODS

We present a phylogenomic examination of a major angiosperm clade (Asteridae) based on sequence data from the nuclear, plastid, and mitochondrial genomes as a means of evaluating currently accepted relationships inferred from the plastome and exploring potential sources of genomic conflict in this group.

RESULTS

We recovered at least five instances of well-supported cytonuclear discordance concerning the placements of major asterid lineages (i.e., Ericales, Oncothecaceae, Aquifoliales, Cassinopsis, and Icacinaceae). We attribute this conflict to a combination of incomplete lineage sorting and hybridization, the latter supported in part by previously inferred whole-genome duplications.

CONCLUSIONS

Our results challenge several long-standing hypotheses of asterid relationships and have implications for morphological character evolution and for the importance of ancient whole-genome duplications in early asterid evolution. These findings also highlight the value of reevaluating broad-scale angiosperm and green-plant phylogeny with nuclear genomic data.

Polyploidy breaks speciation barriers in Australian burrowing frogs Neobatrachus

Polyploidy has played an important role in evolution across the tree of life but it is still unclear how polyploid lineages may persist after their initial formation. While both common and well-studied in plants, polyploidy is rare in animals and generally less understood. The Australian burrowing frog genus Neobatrachus is comprised of six diploid and three polyploid species and offers a powerful animal polyploid model system. We generated exome-capture sequence data from 87 individuals representing all nine species of Neobatrachus to investigate species-level relationships, the origin and inheritance mode of polyploid species, and the population genomic effects of polyploidy on genus-wide demography. We describe rapid speciation of diploid Neobatrachus species and show that the three independently originated polyploid species have tetrasomic or mixed inheritance. We document higher genetic diversity in tetraploids, resulting from widespread gene flow between the tetraploids, asymmetric inter-ploidy gene flow directed from sympatric diploids to tetraploids, and isolation of diploid species from each other. We also constructed models of ecologically suitable areas for each species to investigate the impact of climate on differing ploidy levels. These models suggest substantial change in suitable areas compared to past climate, which correspond to population genomic estimates of demographic histories. We propose that Neobatrachus diploids may be suffering the early genomic impacts of climate-induced habitat loss, while tetraploids appear to be avoiding this fate, possibly due to widespread gene flow. Finally, we demonstrate that Neobatrachus is an attractive model to study the effects of ploidy on the evolution of adaptation in animals.

The Earth BioGenome project: opportunities and challenges for plant genomics and conservation

Sequencing them all. That is the ambitious goal of the recently launched Earth BioGenome project (Proceedings of the National Academy of Sciences of the United States of America, 115, 4325-4333), which aims to produce reference genomes for all eukaryotic species within the next decade. In this perspective, we discuss the opportunities of this project with a plant focus, but highlight also potential limitations. This includes the question of how to best capture all plant diversity, as the green taxon is one of the most complex clades in the tree of life, with over 300 000 species. For this, we highlight four key points: (i) the unique biological insights that could be gained from studying plants, (ii) their apparent underrepresentation in sequencing efforts given the number of threatened species, (iii) the necessity of phylogenomic methods that are aware of differences in genome complexity and quality, and (iv) the accounting for within-species genetic diversity and the historical aspect of conservation genetics.

Integrated likelihood for phylogenomics under a no-common-mechanism model

Background

Multi-locus species phylogeny inference is based on models of sequence evolution on gene trees as well as models of gene tree evolution within the branches of species phylogenies. Almost all statistical methods for this inference task assume a common mechanism across all loci as captured by a single value of each branch length of the species phylogeny.

Results

In this paper, we pursue a “no common mechanism" (NCM) model, where every gene tree evolves according to its own parameters of the species phylogeny. Based on this model, we derive an analytically integrated likelihood of both species trees and networks given the gene trees of multiple loci under an NCM model. We demonstrate the performance of inference under this integrated likelihood on both simulated and biological data.

Conclusions

The model presented here will afford opportunities for exploring connections among various criteria for estimating species phylogenies from multiple, independent loci. Furthermore, further development of this model could potentially result in more efficient methods for searching the space of species phylogenies by focusing solely on the topology of the phylogeny.

Divergence-Based Introgression Polarization

Introgressive hybridization results in the transfer of genetic material between species, often with fitness implications for the recipient species. The development of statistical methods for detecting the signatures of historical introgression in whole-genome data has been a major area of focus. Although existing techniques are able to identify the taxa that exchanged genes during introgression using a four-taxon system, most methods do not explicitly distinguish which taxon served as donor and which as recipient during introgression (i.e., polarization of introgression directionality). Existing methods that do polarize introgression are often only able to do so when there is a fifth taxon available and that taxon is sister to one of the taxa involved in introgression. Here, we present divergence-based introgression polarization (DIP), a method for polarizing introgression using patterns of sequence divergence across whole genomes, which operates in a four-taxon context. Thus, DIP can be applied to infer the directionality of introgression when additional taxa are not available. We use simulations to show that DIP can polarize introgression and identify potential sources of bias in the assignment of directionality, and we apply DIP to a well-described hominin introgression event.

Triad hybridization via a conduit species

Introgressive hybridization can affect the evolution of populations in several important ways. It may retard or reverse divergence of species, enable the development of novel traits, enhance the potential for future evolution by elevating levels of standing variation, create new species, and alleviate inbreeding depression in small populations. Most of what is known of contemporary hybridization in nature comes from the study of pairs of species, either coexisting in the same habitat or distributed parapatrically and separated by a hybrid zone. More rarely, three species form an interbreeding complex (triad), reported in vertebrates, insects, and plants. Often, one species acts as a genetic link or conduit for the passage of genes (alleles) between two others that rarely, if ever, hybridize. Demographic and genetic consequences are unknown. Here we report results of a long-term study of interbreeding Darwin's finches on Daphne Major island, Galápagos. Geospiza fortis acted as a conduit for the passage of genes between two others that have never been observed to interbreed on Daphne: Geospiza fuliginosa, a rare immigrant, and Geospiza scandens, a resident. Microsatellite gene flow from G. fortis into G. scandens increased in frequency during 30 y of favorable ecological conditions, resulting in genetic and morphological convergence. G. fortis, G. scandens, and the derived dihybrids and trihybrids experienced approximately equal fitness. Especially relevant to young adaptive radiations, where species differ principally in ecology and behavior, these findings illustrate how new combinations of genes created by hybridization among three species can enhance the potential for evolutionary change.

Linked-read sequencing identifies abundant microinversions and introgression in the arboviral vector Aedes aegypti

Background

Aedes aegypti is the principal mosquito vector of Zika, dengue, and yellow fever viruses. Two subspecies of Ae. aegypti exhibit phenotypic divergence with regard to habitat, host preference, and vectorial capacity. Chromosomal inversions have been shown to play a major role in adaptation and speciation in dipteran insects and would be of great utility for studies of Ae. aegypti. However, the large and highly repetitive genome of Ae. aegypti makes it difficult to detect inversions with paired-end short-read sequencing data, and polytene chromosome analysis does not provide sufficient resolution to detect chromosome banding patterns indicative of inversions.

Results

To characterize chromosomal diversity in this species, we have carried out deep Illumina sequencing of linked-read (10X Genomics) libraries in order to discover inversion loci as well as SNPs. We analyzed individuals from colonies representing the geographic limits of each subspecies, one contact zone between subspecies, and a closely related sister species. Despite genome-wide SNP divergence and abundant microinversions, we do not find any inversions occurring as fixed differences between subspecies. Many microinversions are found in regions that have introgressed and have captured genes that could impact behavior, such as a cluster of odorant-binding proteins that may play a role in host feeding preference.

Conclusions

Our study shows that inversions are abundant and widely shared among subspecies of Aedes aegypti and that introgression has occurred in regions of secondary contact. This library of 32 novel chromosomal inversions demonstrates the capacity for linked-read sequencing to identify previously intractable genomic rearrangements and provides a foundation for future population genetics studies in this species.

Toward the Definition of Efficacy and Safety Criteria for Advancing Gene Drive-Modified Mosquitoes to Field Testing

Mosquitoes containing gene drive systems are being developed as complementary tools to prevent transmission of malaria and other mosquito-borne diseases. As with any new tool, decision makers and other stakeholders will need to balance risks (safety) and benefits (efficacy) when considering the rationale for testing and deploying gene drive-modified mosquito products. Developers will benefit from standards for judging whether an investigational gene drive product meets acceptability criteria for advancing to field trials. Such standards may be formalized as preferred product characteristics and target product profiles, which describe the desired attributes of the product category and of a particular product, respectively. This report summarizes discussions from two scientific workshops aimed at identifying efficacy and safety characteristics that must be minimally met for an investigational gene drive-modified mosquito product to be deemed viable to move from contained testing to field release and the data that will be needed to support an application for first field release.

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome-wide patterns of within- and between- species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole-genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long-term balancing selection have also been crucial components in shaping patterns of genome-wide variation during the speciation process.