Unexplored microbial diversity from 2,500 food metagenomes and links with the human microbiome.

Complex microbiomes are part of the food we eat and influence our own microbiome, but their diversity remains largely unexplored. Here, we generated the open access curatedFoodMetagenomicData (cFMD) resource by integrating 1,950 newly sequenced and 583 public food metagenomes. We produced 10,899 metagenome-assembled genomes spanning 1,036 prokaryotic and 108 eukaryotic species-level genome bins (SGBs), including 320 previously undescribed taxa. Food SGBs displayed significant microbial diversity within and between food categories. Extension to >20,000 human metagenomes revealed that food SGBs accounted on average for 3% of the adult gut microbiome. Strain-level analysis highlighted potential instances of food-to-gut transmission and intestinal colonization (e.g., Lacticaseibacillus paracasei) as well as SGBs with divergent genomic structures in food and humans (e.g., Streptococcus gallolyticus and Limosilactobabillus mucosae). The cFMD expands our knowledge on food microbiomes, their role in shaping the human microbiome, and supports future uses of metagenomics for food quality, safety, and authentication.

Discovering and exploring the hidden diversity of human gut viruses using highly enriched virome samples.

Viruses are an abundant and crucial component of the human microbiome, but accurately discovering them via metagenomics is still challenging. Currently, the available viral reference genomes poorly represent the diversity in microbiome samples, and expanding such a set of viral references is difficult. As a result, many viruses are still undetectable through metagenomics even when considering the power of de novo metagenomic assembly and binning, as viruses lack universal markers. Here, we describe a novel approach to catalog new viral members of the human gut microbiome and show how the resulting resource improves metagenomic analyses. We retrieved >3,000 viral-like particles (VLP) enriched metagenomic samples (viromes), evaluated the efficiency of the enrichment in each sample to leverage the viromes of highest purity, and applied multiple analysis steps involving assembly and comparison with hundreds of thousands of metagenome-assembled genomes to discover new viral genomes. We reported over 162,000 viral sequences passing quality control from thousands of gut metagenomes and viromes. The great majority of the retrieved viral sequences (~94.4%) were of unknown origin, most had a CRISPR spacer matching host bacteria, and four of them could be detected in >50% of a set of 18,756 gut metagenomes we surveyed. We included the obtained collection of sequences in a new MetaPhlAn 4.1 release, which can quantify reads within a metagenome matching the known and newly uncovered viral diversity. Additionally, we released the viral database for further virome and metagenomic studies of the human microbiome.

A multidisciplinary approach to tackling invasive species: barcoding, morphology, and metataxonomy of the leafhopper Arboridia adanae.

The leafhopper genus Arboridia includes several species that feed on Vitis vinifera and cause leaf chlorosis. We report the first alien Arboridia infestation in Italy in 2021 in an Apulian vineyard. To confirm the taxonomic status of the species responsible for crop damage, and reconstruct its demographic history, we barcoded individuals from Apulia together with Arboridia spp. from Crete (Greece), A. adanae from Central Turkey and other specimens of the presumed sister species, A. dalmatina from Dalmatia (Croatia). Molecular phylogenies and barcoding gap analysis identified clades not associated with sampling locations. This result is incongruent with classical specimen assignment and is further supported by morphological analyses, which did not reveal significant differences among the populations. Therefore, we propose A. dalmatina as a junior synonym of A. adanae, which would become the only grapevine-related Arboridia species in the eastern Mediterranean. To further characterise A. adanae evolution, we performed a molecular clock analysis that suggested a radiation during the Pleistocene glaciations. Finally, to assess whether the Apulian individuals carried microorganisms of agricultural relevance, we sequenced their bacterial microbiota using 16S rRNA amplicon sequencing identifying three phytopathogens not generally associated with Arboridia activities as well as Wolbachia in one Apulian haplogroup. We discuss the agricultural implications of this infestation.

Shallow Whole-Genome Sequencing of Aedes japonicus and Aedes koreicus from Italy and an Updated Picture of Their Evolution Based on Mitogenomics and Barcoding.

Aedes japonicus and Aedes koreicus are two invasive mosquitoes native to East Asia that are quickly establishing in temperate regions of Europe. Both species are vectors of arboviruses, but we currently lack a clear understanding of their evolution. Here, we present new short-read, shallow genome sequencing of A. japonicus and A. koreicus individuals from northern Italy, which we used for downstream phylogenetic and barcode analyses. We explored associated microbial DNA and found high occurrences of Delftia bacteria in both samples, but neither Asaia nor Wolbachia. We then assembled complete mitogenomes and used these data to infer divergence times estimating the split of A. japonicus from A. koreicus in the Oligocene, which was more recent than that previously reported using mitochondrial markers. We recover a younger age for most other nodes within Aedini and other Culicidae. COI barcoding and phylogenetic analyses indicate that A. japonicus yaeyamensis, A. japonicus amamiensis, and the two A. koreicus sampled from Europe should be considered as separate species within a monophyletic species complex. Our studies further clarify the evolution of A. japonicus and A. koreicus, and indicate the need to obtain whole-genome data from putative species in order to disentangle their complex patterns of evolution.

Spatial and temporal dynamics of West Nile virus between Africa and Europe.

It is unclear whether West Nile virus (WNV) circulates between Africa and Europe, despite numerous studies supporting an African origin and high transmission in Europe. We integrated genomic data with geographic observations and phylogenetic and phylogeographic inferences to uncover the spatial and temporal viral dynamics of WNV between these two continents. We focused our analysis towards WNV lineages 1 (L1) and 2 (L2), the most spatially widespread and pathogenic WNV lineages. Our study shows a Northern-Western African origin of L1, with back-and-forth exchanges between West Africa and Southern-Western Europe; and a Southern African origin of L2, with one main introduction from South Africa to Europe, and no back introductions observed. We also noticed a potential overlap between L1 and L2 Eastern and Western phylogeography and two Afro-Palearctic bird migratory flyways. Future studies linking avian and mosquito species susceptibility, migratory connectivity patterns, and phylogeographic inference are suggested to elucidate the dynamics of emerging viruses.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

First genome-wide data from Italian European beech (Fagus sylvatica L.): Strong and ancient differentiation between Alps and Apennines.

The European beech (Fagus sylvatica L.) is one of the most widespread forest trees in Europe whose distribution and intraspecific diversity has been largely shaped by repeated glacial cycles. Previous studies, mainly based on palaeobotanical evidence and a limited set of chloroplast and nuclear genetic markers, highlighted a complex phylogeographic scenario, with southern and western Europe characterized by a rather heterogeneous genetic structure, as a result of recolonization from different glacial refugia. Despite its ecological and economic importance, the genome of this broad-leaved tree has only recently been assembled, and its intra-species genomic diversity is still largely unexplored. Here, we performed whole-genome resequencing of nine Italian beech individuals sampled from two stands located in the Alpine and Apennine mountain ranges. We investigated patterns of genetic diversity at chloroplast, mitochondrial and nuclear genomes and we used chloroplast genomes to reconstruct a temporally-resolved phylogeny. Results allowed us to test European beech differentiation on a whole-genome level and to accurately date their divergence time. Our results showed comparable, relatively high levels of genomic diversity in the two populations and highlighted a clear differentiation at chloroplast, mitochondrial and nuclear genomes. The molecular clock analysis indicated an ancient split between the Alpine and Apennine populations, occurred between the Günz and the Riss glaciations (approximately 660 kyrs ago), suggesting a long history of separation for the two gene pools. This information has important conservation implications in the context of adaptation to ongoing climate changes.

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae).

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

Investigating the microbial community of Cacopsylla spp. as potential factor in vector competence of phytoplasma.

Phytoplasmas are obligatory intracellular bacteria that colonize the phloem of many plant species and cause hundreds of plant diseases worldwide. In nature, phytoplasmas are primarily transmitted by hemipteran vectors. While all phloem-feeding insects could in principle transmit phytoplasmas, only a limited number of species have been confirmed as vectors. Knowledge about factors that might determine the vector capacity is currently scarce. Here, we characterized the microbiomes of vector and non-vector species of apple proliferation (AP) phytoplasma 'Candidatus Phytoplasma mali' to investigate their potential role in the vector capacity of the host. We performed high-throughput 16S rRNA metabarcoding of the two principal AP-vectors Cacopsylla picta and Cacopsylla melanoneura and eight Cacopsylla species, which are not AP-vectors but co-occur in apple orchards. The microbiomes of all species are dominated by Carsonella, the primary endosymbiont of psyllids and a second uncharacterized Enterobacteriaceae endosymbiont. Each Cacopsylla species harboured a species-specific phylotype of both symbionts. Moreover, we investigated differences between the microbiomes of AP-vector versus non-vector species and identified the predominant endosymbionts but also Wolbachia and several minor taxa as potential indicator species. Our study highlights the importance of considering the microbiome in future investigations of potential factors influencing host vector competence. We investigated the potential role of symbiotic bacteria in the acquisition and transmission of phytoplasma. By comparing the two main psyillid vector species of Apple proliferation (AP) phytoplasma and eight co-occurring species, which are not able to vector AP-phytoplasma, we found differences in the microbial communities of AP-vector and non-vector species, which appear to be driven by the predominant symbionts in both vector species and Wolbachia and several minor taxa in the non-vector species. In contrast, infection with AP-phytoplasma did not affect microbiome composition in both vector species. Our study provides new insights into the endosymbiont diversity of Cacopsylla spp. and highlights the importance of considering the microbiome when investigating potential factors influencing host vector competence.

The Impact of Fast Radiation on the Phylogeny of Bactrocera Fruit Flies as Revealed by Multiple Evolutionary Models and Mutation Rate-Calibrated Clock.

Several true fruit flies (Tephritidae) cause major damage to agriculture worldwide. Among them, species of the genus Bactrocera are extensively studied to understand the traits associated with their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but several nodes of the Bactrocera phylogeny are still controversial, especially concerning the reciprocal affinities of the two major pests B. dorsalis and B. tryoni. Here, we analyzed a newly assembled genomic-scaled dataset using different models of evolution to infer a phylogenomic backbone of ten representative Bactrocera species and two outgroups. We further provide the first genome-scaled inference of their divergence by calibrating the clock using fossil records and the spontaneous mutation rate. The results reveal a closer relationship of B. dorsalis with B. latifrons than to B. tryoni, contrary to what was previously supported by mitochondrial-based phylogenies. By employing coalescent-aware and heterogeneous evolutionary models, we show that this incongruence likely derives from a hitherto undetected systematic error, exacerbated by incomplete lineage sorting and possibly hybridization. This agrees with our clock analysis, which supports a rapid and recent radiation of the clade to which B. dorsalis, B. latifrons and B. tryoni belong. These results provide a new picture of Bactrocera phylogeny that can serve as the basis for future comparative analyses.

Large-scale mitogenome sequencing reveals consecutive expansions of domestic taurine cattle and supports sporadic aurochs introgression.

The contribution of domestic cattle in human societies is enormous, making cattle, along with other essential benefits, the economically most important domestic animal in the world today. To expand existing knowledge on cattle domestication and mitogenome diversity, we performed a comprehensive complete mitogenome analysis of the species (802 sequences, 114 breeds). A large sample was collected in South-east Europe, an important agricultural gateway to Europe during Neolithization and a region rich in cattle biodiversity. We found 1725 polymorphic sites (810 singletons, 853 parsimony-informative sites and 57 indels), 701 unique haplotypes, a haplotype diversity of 0.9995 and a nucleotide diversity of 0.0015. In addition to the dominant T3 and several rare haplogroups (Q, T5, T4, T2 and T1), we have identified maternal line in Austrian Murbodner cattle that possess surviving aurochs' mitochondria haplotype P1 that diverged prior to the Neolithization process. This is convincing evidence for rare female-mediated adaptive introgression of wild aurochs into domesticated cattle in Europe. We revalidated the existing haplogroup classification and provided Bayesian phylogenetic inference with a more precise estimated divergence time than previously available. Occasionally, classification based on partial mitogenomes was not reliable; for example, some individuals with haplogroups P and T5 were not recognized based on D-loop information. Bayesian skyline plot estimates (median) show that the earliest population growth began before domestication in cattle with haplogroup T2, followed by Q (~10.0-9.5 kyBP), whereas cattle with T3 (~7.5 kyBP) and T1 (~3.0-2.5 kyBP) expanded later. Overall, our results support the existence of interactions between aurochs and cattle during domestication and dispersal of cattle in the past, contribute to the conservation of maternal cattle diversity and enable functional analyses of the surviving aurochs P1 mitogenome.

On the Origin and Propagation of the COVID-19 Outbreak in the Italian Province of Trento, a Tourist Region of Northern Italy.

BACKGROUND: Trentino is an Italian province with a tourism-based economy, bordering the regions of Lombardy and Veneto, where the two earliest and largest outbreaks of COVID-19 occurred in Italy. The earliest cases in Trentino were reported in the first week of March 2020, with most of the cases occurring in the winter sport areas in the Dolomites mountain range. The number of reported cases decreased over the summer months and was followed by a second wave in the autumn and winter of 2020. METHODS: we performed high-coverage Oxford Nanopore sequencing of 253 positive SARS-CoV-2 swabs collected in Trentino between March and December 2020. RESULTS: in this work, we analyzed genome sequences to trace the routes through which the virus entered the area, and assessed whether the autumnal resurgence could be attributed to lineages persisting undetected during summer, or as a consequence of new introductions. CONCLUSIONS: Comparing the draft genomes analyzed with a large selection of European sequences retrieved from GISAID we found that multiple introductions of the virus occurred at the early stage of the epidemics; the two epidemic waves were unrelated; the second wave was due to reintroductions of the virus in summer when traveling restrictions were uplifted.

Epidemiological and Evolutionary Analysis of West Nile Virus Lineage 2 in Italy.

West Nile virus (WNV) is a mosquito-borne virus potentially causing serious illness in humans and other animals. Since 2004, several studies have highlighted the progressive spread of WNV Lineage 2 (L2) in Europe, with Italy being one of the countries with the highest number of cases of West Nile disease reported. In this paper, we give an overview of the epidemiological and genetic features characterising the spread and evolution of WNV L2 in Italy, leveraging data obtained from national surveillance activities between 2011 and 2021, including 46 newly assembled genomes that were analysed under both phylogeographic and phylodynamic frameworks. In addition, to better understand the seasonal patterns of the virus, we used a machine learning model predicting areas at high-risk of WNV spread. Our results show a progressive increase in WNV L2 in Italy, clarifying the dynamics of interregional circulation, with no significant introductions from other countries in recent years. Moreover, the predicting model identified the presence of suitable conditions for the 2022 earlier and wider spread of WNV in Italy, underlining the importance of using quantitative models for early warning detection of WNV outbreaks. Taken together, these findings can be used as a reference to develop new strategies to mitigate the impact of the pathogen on human and other animal health in endemic areas and new regions.

Hallstatt miners consumed blue cheese and beer during the Iron Age and retained a non-Westernized gut microbiome until the Baroque period.

We subjected human paleofeces dating from the Bronze Age to the Baroque period (18th century AD) to in-depth microscopic, metagenomic, and proteomic analyses. The paleofeces were preserved in the underground salt mines of the UNESCO World Heritage site of Hallstatt in Austria. This allowed us to reconstruct the diet of the former population and gain insights into their ancient gut microbiome composition. Our dietary survey identified bran and glumes of different cereals as some of the most prevalent plant fragments. This highly fibrous, carbohydrate-rich diet was supplemented with proteins from broad beans and occasionally with fruits, nuts, or animal food products. Due to these traditional dietary habits, all ancient miners up to the Baroque period have gut microbiome structures akin to modern non-Westernized individuals whose diets are also mainly composed of unprocessed foods and fresh fruits and vegetables. This may indicate a shift in the gut community composition of modern Westernized populations due to quite recent dietary and lifestyle changes. When we extended our microbial survey to fungi present in the paleofeces, in one of the Iron Age samples, we observed a high abundance of Penicillium roqueforti and Saccharomyces cerevisiae DNA. Genome-wide analysis indicates that both fungi were involved in food fermentation and provides the first molecular evidence for blue cheese and beer consumption in Iron Age Europe.

Metagenomic analysis of ancient dental calculus reveals unexplored diversity of oral archaeal Methanobrevibacter.

BACKGROUND: Dental calculus (mineralised dental plaque) preserves many types of microfossils and biomolecules, including microbial and host DNA, and ancient calculus are thus an important source of information regarding our ancestral human oral microbiome. In this study, we taxonomically characterised the dental calculus microbiome from 20 ancient human skeletal remains originating from Trentino-South Tyrol, Italy, dating from the Neolithic (6000-3500 BCE) to the Early Middle Ages (400-1000 CE). RESULTS: We found a high abundance of the archaeal genus Methanobrevibacter in the calculus. However, only a fraction of the sequences showed high similarity to Methanobrevibacter oralis, the only described Methanobrevibacter species in the human oral microbiome so far. To further investigate the diversity of this genus, we used de novo metagenome assembly to reconstruct 11 Methanobrevibacter genomes from the ancient calculus samples. Besides the presence of M. oralis in one of the samples, our phylogenetic analysis revealed two hitherto uncharacterised and unnamed oral Methanobrevibacter species that are prevalent in ancient calculus samples sampled from a broad range of geographical locations and time periods. CONCLUSIONS: We have shown the potential of using de novo metagenomic assembly on ancient samples to explore microbial diversity and evolution. Our study suggests that there has been a possible shift in the human oral microbiome member Methanobrevibacter over the last millennia. Video abstract.

The evolution of insect biodiversity.

Insects comprise over half of all described animal species. Together with the Protura (coneheads), Collembola (springtails) and Diplura (two-pronged bristletails), insects form the Hexapoda, a terrestrial arthropod lineage characterised by possessing six legs. Exponential growth of genome-scale data for the hexapods has substantially altered our understanding of the origin and evolution of insect biodiversity. Phylogenomics has provided a new framework for reconstructing insect evolutionary history, resolving their position among the arthropods and some long-standing internal controversies such as the placement of the termites, twisted-winged insects, lice and fleas. However, despite the greatly increased size of phylogenomic datasets, contentious relationships among key insect clades remain unresolved. Further advances in insect phylogeny cannot rely on increased depth and breadth of genome and taxon sequencing. Improved modelling of the substitution process is fundamental to countering tree-reconstruction artefacts, while gene content, modelling of duplications and deletions, and comparative morphology all provide complementary lines of evidence to test hypotheses emerging from the analysis of sequence data. Finally, the integration of molecular and morphological data is key to the incorporation of fossil species within insect phylogeny. The emerging integrated framework of insect evolution will help explain the origins of insect megadiversity in terms of the evolution of their body plan, species diversity and ecology. Future studies of insect phylogeny should build upon an experimental, hypothesis-driven approach where the robustness of hypotheses generated is tested against increasingly realistic evolutionary models as well as complementary sources of phylogenetic evidence.

Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Viviparity and habitat restrictions may influence the evolution of male reproductive genes in tsetse fly (Glossina) species.

BACKGROUND: Glossina species (tsetse flies), the sole vectors of African trypanosomes, maintained along their long evolutionary history a unique reproductive strategy, adenotrophic viviparity. Viviparity reduces their reproductive rate and, as such, imposes strong selective pressures on males for reproductive success. These species live in sub-Saharan Africa, where the distributions of the main sub-genera Fusca, Morsitans, and Palpalis are restricted to forest, savannah, and riverine habitats, respectively. Here we aim at identifying the evolutionary patterns of the male reproductive genes of six species belonging to these three main sub-genera. We then interpreted the different patterns we found across the species in the light of viviparity and the specific habitat restrictions, which are known to shape reproductive behavior. RESULTS: We used a comparative genomic approach to build consensus evolutionary trees that portray the selective pressure acting on the male reproductive genes in these lineages. Such trees reflect the long and divergent demographic history that led to an allopatric distribution of the Fusca, Morsitans, and Palpalis species groups. A dataset of over 1700 male reproductive genes remained conserved over the long evolutionary time scale (estimated at 26.7 million years) across the genomes of the six species. We suggest that this conservation may result from strong functional selective pressure on the male imposed by viviparity. It is noteworthy that more than half of these conserved genes are novel sequences that are unique to the Glossina genus and are candidates for selection in the different lineages. CONCLUSIONS: Tsetse flies represent a model to interpret the evolution and differentiation of male reproductive biology under different, but complementary, perspectives. In the light of viviparity, we must take into account that these genes are constrained by a post-fertilization arena for genomic conflicts created by viviparity and absent in ovipositing species. This constraint implies a continuous antagonistic co-evolution between the parental genomes, thus accelerating inter-population post-zygotic isolation and, ultimately, favoring speciation. Ecological restrictions that affect reproductive behavior may further shape such antagonistic co-evolution.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species.

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe.

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.