Metagenomic analysis of 16S ribosomal RNA has been used to profile microbial communities at high resolution, and to examine their association with host diet or diseases. We examined the oral and gut microbiome composition of two captive koalas to determine whether bacterial communities are unusual in this species, given that their diet consists almost exclusively of Eucalyptus leaves. Despite a highly specialized diet, koala oral and gut microbiomes were similar in composition to the microbiomes from the same body regions of other mammals. Rectal swabs contained all of the diversity present in faecal samples, along with additional taxa, suggesting that faecal bacterial communities may merely subsample the gut bacterial diversity. Furthermore, the faecal microbiomes of the captive koalas were similar to those reported for wild koalas, suggesting that captivity may not compromise koala microbial health. Since koalas frequently suffer from ocular diseases caused by Chlamydia infection, we also examined the eye microbiome composition of two captive koalas, establishing the healthy baseline for this body part. The eye microbial community was very diverse, similar to other mammalian ocular microbiomes but with an unusually high representation of bacteria from the family Phyllobacteriaceae.

Tick-borne encephalitis (TBE) is a severe zoonotic neurological disease endemic in northeast Italy since 1992. In the Province of Trento, a sharp increase in TBE incidence has been recorded since 2012, despite the vaccination efforts. To assess current TBE infection hazard in this area, we applied an integrated approach combining the distribution of human cases, the seroprevalence of tick-borne encephalitis virus (TBEV) in sentinel hosts and the screening of questing ticks for TBEV. A total of 706 goat sera from 69 farms were screened for TBEV-specific antibodies resulting in 5 positive farms, while the location of human cases was provided by the local Public Health Agency. Tick sampling was concentrated in areas where TBEV circulation was suggested by either seroprevalence in goats or human cases, resulting in 2,410 Ixodes ricinus collected and analyzed by real-time RT-PCR. Four tick samples from 2 areas with record of human cases were positive to TBEV corresponding to a 0.17% prevalence in the region, while risk areas suggested by serology on goats were not confirmed by tick screening. Our results revealed an increase in TBEV prevalence in ticks and the emergence of new active TBE foci, compared to previous surveys, and demonstrated the importance of an integrated approach for TBE risk assessment. A phylogenetic analysis of the partial E gene confirmed that the European TBEV subtype is circulating in northeast Italy and suggested that the different Italian TBEV strains originated independently as a result of different introductions from neighbouring countries, presumably through migratory birds.

Gibbon ape leukemia virus (GALV) and koala retrovirus (KoRV) most likely originated from a cross-species transmission of an ancestral retrovirus into koalas and gibbons via one or more intermediate as-yet-unknown hosts. A virus highly similar to GALV has been identified in an Australian native rodent (Melomys burtoni) after extensive screening of Australian wildlife. GALV-like viruses have also been discovered in several Southeast Asian species, although screening has not been extensive and viruses discovered to date are only distantly related to GALV. We therefore screened 26 Southeast Asian rodent species for KoRV-and GALV-like sequences, using hybridization capture and high-throughput sequencing, in the attempt to identify potential GALV and KoRV hosts. Only the individuals belonging to a newly discovered subspecies of Melomys burtoni from Indonesia were positive, yielding an endogenous provirus very closely related to a strain of GALV. The sequence of the critical receptor domain for GALV infection in the Indonesian M. burtoni subsp. was consistent with the susceptibility of the species to GALV infection. The second record of a GALV in M. burtoni provides further evidence that M. burtoni, and potentially other lineages within the widespread subfamily Murinae, may play a role in the spread of GALV-like viruses. The discovery of a GALV in the most western part of the Australo-Papuan distribution of M. burtoni, specifically in a transitional zone between Asia and Australia (Wallacea), may be relevant to the cross-species transmission to gibbons in Southeast Asia and broadens the known distribution of GALVs in wild rodents. IMPORTANCEGibbon ape leukemia virus (GALV) and the koala retrovirus (KoRV) are very closely related, yet their hosts neither are closely related nor overlap geographically. Direct cross-species infection between koalas and gibbons is unlikely. Therefore, GALV and KoRV may have arisen via a cross-species transfer from an intermediate host whose range overlaps those of both gibbons and koalas. Using hybridization capture and high-throughput sequencing, we have screened a wide range of rodent candidate hosts from Southeast Asia for KoRV-and GALV-like sequences. Only a Melomys burtoni subspecies from Wallacea (Indonesia) was positive for GALV. We report the genome sequence of this newly identified GALV, the critical domain for infection of its potential cellular receptor, and its phylogenetic relationships with the other previously characterized GALVs. We hypothesize that Melomys burtoni, and potentially related lineages with an Australo-Papuan distribution, may have played a key role in crossspecies transmission to other taxa. The evolutionary mechanisms involved in cross-species transmissions (CSTs) of viruses are complex and generally poorly understood. Viral evolution, host contact rates, biological similarity in host defense systems (receptors, viral restriction factors), and host evolutionary relationships have been proposed as key factors in CST rates and outcomes (1). However, the...

Since it has been understood that gut microbiota of vector mosquitoes can influence their vector competence, efforts have been undertaken to develop new control strategies based on host microbiota manipulation, and aimed at suppressing the vector population or replacing it with a less competent one. For the proper design of such control strategies it is necessary to know the microbiota composition of the target vector species, how it is acquired, and how it changes throughout the host's life cycle. In this study, 16S rRNA amplicon high-throughput sequencing was used to characterize the changes in microbiota from the aquatic environment (larval breeding sites) to the different developmental stages of field-collected Aedes koreicus in Italy, an emerging invasive mosquito species in Europe and a potential vector of several pathogens. The bacterial communities of the aquatic breeding sites, larvae, pupae and adults showed distinctive structures to one another. Indeed, 84% of community members were unique to a given sample type. Nevertheless, almost 40% of the sequences generated were assigned to bacteria detected in all sample types, suggesting the importance of bacteria transstadially transmitted from water to the adult stage in constituting mosquito microbiota. Among these, genus C39 largely constituted water microbiota, family Burkholderiaceae was the most abundant in larvae and pupae, and genus Asaia dominated adult communities. In addition, Asaia constituted a core microbiota across all sample types. Our results suggest that the microbiota of Ae. koreicus mosquitoes is composed by a community which derives from the aquatic bacteria of the larval breeding sites, is then filtered by the larval gut, where only certain members are able to persist, rearranged by metamorphosis and finally modified by the change in diet at the adult stage. Understanding how the microbiota of Ae. koreicus changes through the mosquito life cycle represents a first step in selecting bacterial candidates for use in microbiotabased intervention measures for this species. The properties which Asaia exhibits in this species, such as dominance, high prevalence and transstadial transmission, prevent the use of Wolbachia but make Asaia an ideal candidate for paratransgenesis.

Gibbon ape leukemia viruses (GALVs) are part of a larger group of pathogenic gammaretroviruses present across phylogenetically diverse host species of Australasian mammals. Despite the biomedical utility of GALVs as viral vectors and in cancer gene therapy, full genome sequences have not been determined for all of the five identified GALV strains, nor has a comprehensive evolutionary analysis been performed. We therefore generated complete genomic sequences for each GALV strain using hybridization capture and high-throughput sequencing. The four strains of GALV isolated from gibbons formed a monophyletic clade that was closely related to the woolly monkey virus (WMV), which is a GALV strain that likely originated in a gibbon host. The GALV-WMV clade in turn formed a sister group to the koala retroviruses (KoRVs). Genomic signatures of episodic diversifying selection were detected among the gammaretroviruses with concentration in the env gene across the GALV strains that were particularly oncogenic and KoRV strains that were potentially exogenous, likely reflecting their adaptation to the host immune system. In vitro studies involving vectors chimeric between GALV and KoRV-B established that variable regions A and B of the surface unit of the envelope determine which receptor is used by a viral strain to enter host cells. IMPORTANCEThe gibbon ape leukemia viruses (GALVs) are among the most medically relevant retroviruses due to their use as viral vectors for gene transfer and in cancer gene therapy. Despite their importance, full genome sequences have not been determined for the majority of primate isolates, nor has comprehensive evolutionary analysis been performed, despite evidence that the viruses are facing complex selective pressures associated with cross-species transmission. Using hybridization capture and high-throughput sequencing, we report here the full genome sequences of all the GALV strains and demonstrate that diversifying selection is acting on them, particularly in the envelope gene in functionally important domains, suggesting that host immune pressure is shaping GALV evolution. G ibbon ape leukemia virus (GALV) is an exogenous gammaretrovirus associated with hematopoietic neoplasms in captive colonies of white-handed gibbon (Hylobates lar). Five strains of GALV have been isolated from gibbons. The first was isolated from an animal with lymphocytic leukemia in a colony at the San Francisco Medical Center (strain SF) (1, 2). GALV was later isolated from gibbons displaying malignant tumors, notably an individual gibbon with granulocytic leukemia, at the Southeast Asia Treaty Organization Medical Research Laboratory in Bangkok, Thailand (strain SEATO) (3, 4), and another gibbon with lymphocytic leukemia from a colony on Hall's Island, near Bermuda (strain GALV-H) (5, 6). The Brain strain was isolated from two healthy gibbons injected with brain extracts from human patients with kuru and from an uninoculated cage mate (7). The SEATO strain has been shown to cause chronic myelogenous leukemia...

SummaryEnvironmental DNA (eDNA) and its subdiscipline, invertebrate-derived DNA (iDNA) have been used to survey biodiversity non-invasively [1,2]. Water is ubiquitous in most ecosystems, and, among invertebrates, terrestrial haematophagous leeches are abundant and can be easily collected in many tropical rainforests [3,4]. Such non-invasive nucleic acid sources can mitigate difficulties of obtaining wildlife samples, particularly in remote areas or for rare species. Recently, eDNA/iDNA sources have been applied to monitoring specific wildlife pathogens [5,6]. However, previous studies have focused on known pathogens, whereas most wildlife pathogens are uncharacterized and unknown. Non-invasive approaches to monitoring known and novel pathogens may be of particular benefit in ecosystems prone to viral emergence, many of which occur in areas where invasive sampling is challenging, such as tropical rainforests. Here, we show that both eDNA from natural waterholes, and iDNA from terrestrial haematophagous leeches, can be used to detect unknown viruses circulating in mammalian hosts (Figure 1). Using a curated set of RNA oligonucleotides based on the ViroChip microarray assay [7] as baits in a hybridization capture system, multiple mammalian RNA and DNA viruses were detected from both eDNA and iDNA samples. Congruence was found between host DNA assignment and viruses identified in leeches, and between animals observed visiting the waterholes and the viruses detected. Our results demonstrate that eDNA/iDNA samples may represent an effective non-invasive resource for studying wildlife viral diversity. Several of the detected viruses were novel, highlighting the potential of eDNA/iDNA for epidemiological analysis of emerging viruses prior to their emergence.HighlightsEnvironmental DNA (water and blood-sucking leeches) provided a non-invasive method of screening wildlife for virusesA comprehensive viral RNA oligonucleotide bait set was developed to capture known and unknown mammalian virus diversityLeech blood meal host determination and viruses identified were congruentViruses determined from water correlated with known and observed species visiting the water sourcesIn briefAlfano, Dayaram, et al. demonstrate that environmental DNA from southeast Asian leech bloodmeals and waterholes from Africa and Mongolia can be used as to detect viruses circulating in wildlife. These nucleic acid sources may represent an effective non-invasive resource for studying wildlife viral diversity and emerging viruses pre-emergence.

1. Environmental DNA (eDNA) and invertebrate-derived DNA (iDNA) have been used to survey biodiversity non-invasively to mitigate difficulties of obtaining wildlife samples, particularly in remote areas or for rare species. Recently, eDNA/iDNA have been applied to monitor known wildlife pathogens, however, most wildlife pathogens are unknown and often evolutionarily divergent.2. To detect and identify known and novel mammalian viruses from eDNA/iDNA sources, we used a curated set of RNA oligonucleotides as viral baits in a hybridization capture system coupled with high throughput sequencing.3. We detected multiple known and novel mammalian RNA and DNA viruses from multiple viral families from both waterhole eDNA and leech derived iDNA. Congruence was found between detected hosts and viruses identified in leeches and waterholes. 4. Our results demonstrate that eDNA/iDNA samples represent an effective non-invasive resource for studying wildlife viral diversity and for detecting novel potentially zoonotic viruses prior to their emergence.

Background: Domestication is a complex, multi-stage and species-specific process that results in organisms living close to humans. In the arboviral vectorAedes aegyptiadaptation to living in proximity with anthropogenic environments has been recognized as a major evolutionary shift, separating a generalist form,Aedes aegypti formosus(Aaf), from the domestic formAedes aegypti aegypti(Aaa), which tends to deposit eggs artificial containers and bite humans for a blood meal. These behaviors enhance the mosquito vectorial capacity. The extent to which domestication has impacted theAe. aegyptigenome has not been thoroughly investigated yet. Results: Taking advantage of two forms' distinct and historically documented geographic distributions, we analyzed the genomes of 634 worldwideAe. aegyptimosquitoes. Using more than 300 million high-confidence SNPs, we found a unique origin for all out-of-AfricaAe. aegyptimosquitoes, with no evidence of admixture events in Africa, apart from Kenya. A group of 1033 genes were under positive selection only in out-of-Africa mosquitoes and 236 genes had nonsynonymous mutations, occurring at statistically different frequencies in Aaa and Aaf mosquitoes. Conclusion: We identified a clear signal of genetic differentiation between Aaa and Aaf, circumscribed to a catalogue of candidate genes. TheseAaa molecular signaturegenes extend beyond chemosensory genes to genes linked to neuronal and hormonal functions. This suggests that the behavioral shift to domestication may rely on the fine regulation of metabolic and neuronal functions, more than the role of a few significant genes. Our results also provide the foundation to investigate new targets for the control ofAe. aegyptipopulations.