Search

The question of how interactions between the gut microbiome and vertebrate hosts contribute to host adaptation and speciation is one of the major problems in current evolutionary research. Using bacteriome and mycobiome metabarcoding, we examined how these two components of the gut microbiota vary with the degree of host admixture in secondary contact between two house mouse subspecies (Mus musculus musculusandM. m. domesticus). We used a large dataset collected at two replicates of the hybrid zone and model-based statistical analyses to ensure the robustness of our results. Assuming that the microbiota of wild hosts suffers from spatial autocorrelation, we directly compared the results of statistical models that were spatially naive with those that accounted for spatial autocorrelation. We showed that neglecting spatial autocorrelation can drastically affect the results and lead to misleading conclusions. The spatial analyses showed little difference between subspecies, both in microbiome composition and in individual bacterial lineages. Similarly, the degree of admixture had minimal effects on the gut bacteriome and mycobiome and was caused by changes in a few microbial lineages that correspond to the common symbionts of free-living house mice. In contrast to previous studies, these data do not support the hypothesis that the microbiota plays an important role in host reproductive isolation in this particular model system.

Data on the gut microbiota (GM) of wild animals are key to studies on evolutionary biology (host–GM interactions under natural selection), ecology and conservation biology (GM as a fitness component closely connected to the environment). Wildlife GM sampling often requires non-invasive techniques or sampling from dead animals. In a controlled experiment profiling microbial 16S rRNA in 52 house mice (Mus musculus) from eight families and four genetic backgrounds, we studied the effects of live- and snap-trapping on small mammal GM and evaluated the suitability of microbiota from non-fresh faeces as a proxy for caecal GM. We compared CM from individuals sampled 16–18 h after death with those in live traps and caged controls, and caecal and faecal GM collected from mice in live-traps. Sampling delay did not affect GM composition, validating data from fresh cadavers or snap-trapped animals. Animals trapped overnight displayed a slight but significant difference in GM composition to the caged controls, though the change only had negligible effect on GM diversity, composition and inter-individual divergence. Hence, the trapping process appears not to bias GM profiling. Despite their significant difference, caecal and faecal microbiota were correlated in composition and, to a lesser extent, diversity. Both showed congruent patterns of inter-individual divergence following the natural structure of the dataset. Thus, the faecal microbiome represents a good non-invasive proxy of the caecal microbiome, making it suitable for detecting biologically relevant patterns. However, care should be taken when analysing mixed datasets containing both faecal and caecal samples.

Crocidura (Eulipotyphla, Soricidae) is the most species-rich genus among mammals, with high cryptic diversity and complicated taxonomy. The hirta-flavescens group of Crocidura represents the most abundant and widespread shrews in savannahs of eastern and southern Africa, making them a suitable phylogeographical model for assessing the role of paleoclimatic changes on current biodiversity in open African habitats. We present the first comprehensive study on the phylogeography, evolutionary history, geographical distribution, systematics, and taxonomy of the group, using the integration of mitochondrial, genome-wide (ddRAD sequencing), morpho-logical and morphometrical data collected from specimens over most of the known geographic distribution. Our genomic data confirmed the monophyly of this group and its sister relationship with the olivieri group of Cro-cidura. There is a substantial genetic variation within the hirta-flavescens group, with three highly supported clades showing parapatric distribution and which can be distinguished morphologically: C. hirta, distributed in both the Zambezian and Somali-Masai bioregions, C. flavescens, known from South Africa and south-western Zambia, and C. cf. flavescens, which is known to occur only in central and western Tanzania. Morphometric data revealed relatively minor differences between C. hirta and C. cf. flavescens, but they differ in the colouration of the pelage. Diversification of the hirta-flavescens group has most likely happened during phases of grassland expansion and contraction during Plio-Pleistocene climatic cycles. Eastern African Rift system, rivers, and the distinctiveness of Zambezian and Somali-Masai bioregions seem to have also shaped the pattern of their di-versity, which is very similar to sympatric rodent species living in open habitats. Finally, we review the group's taxonomy and propose to revalidate C. bloyeti, currently a synonym of C. hirta, including the specimens treated as C. cf. flavescens.

Mammarenaviruses are bi-segmented RNA viruses. They encompass viruses responsible for several severe diseases in humans. While performing a de novo assembly of a new virus found in a wild single-striped grass mouse in Tanzania, we found a single S but two divergent L segments. Natural co-infections, common within reptarenaviruses in captivity, were never reported for mammarenaviruses and never in a wild sample. This finding can have implications for virus evolution as co-infection could trigger viral recombination/reassortment in natural reservoirs.

The gastrointestinal microbiota (GM) is considered an important component of the vertebrate holobiont. GM–host interactions influence the fitness of holobionts and are, therefore, an integral part of evolution. The house mouse is a prominent model for GM–host interactions, and evidence suggests a role for GM in mouse speciation. However, previous studies based on short 16S rRNA GM profiles of wild house mouse subspecies failed to detect GM divergence, which is a prerequisite for the inclusion of GM in Dobzhansky–Muller incompatibilities. Here, we used standard 16S rRNA GM profiling in two mouse subspecies, Mus musculus musculus and M. m. domesticus, including the intestinal mucosa and content of three gut sections (ileum, caecum, and colon). We reduced environmental variability by sampling GM in the offspring of wild mice bred under seminatural conditions. Although the breeding conditions allowed a contact between the subspecies, we found a clear differentiation of GM between them, in all three gut sections. Differentiation was mainly driven by several Helicobacters and two H. ganmani variants showed a signal of codivergence with their hosts. Helicobacters represent promising candidates for studying GM–host coadaptations and the fitness effects of their interactions.

Polymorphism of the major histocompatibility complex (MHC), DAB1 gene was characterized for the first time in the European bitterling (Rhodeus amarus), a freshwater fish employed in studies of host-parasite coevolution and mate choice, taking advantage of newly designed primers coupled with high-throughput amplicon sequencing. Across 221 genotyped individuals, we detected 1–4 variants per fish, with 28% individuals possessing 3–4 variants. We identified 36 DAB1 variants, and they showed high sequence diversity mostly located within predicted antigen-binding sites, and both global and codon-specific excess of non-synonymous mutations. Despite deep divergence between two major allelic lineages, functional diversity was surprisingly low (3 supertypes). Overall, these findings suggest the role of positive and balancing selection in promotion and long-time maintenance of DAB1 polymorphism. Further investigations will clarify the role of pathogen-mediated selection to drive the evolution of DAB1 variation.

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

The close phylogenetic relationship between humans and nonhuman primates (NHPs) can result in a high potential for pathogen exchange. In recent decades, NHP and human interactions have become more frequent due to increasing habitat encroachment and ecotourism. Strongylid communities, which include members of several genera, are typically found in NHPs. Using optimized high-throughput sequencing for strain-level identification of primate strongylids, we studied the structure of strongylid communities in NHPs and humans co-habiting a tropical forest ecosystem in the Central African Republic. General taxonomic assignment of 85 ITS-2 haplotypes indicated that the studied primates harbour at least nine genera of strongylid nematodes, with Oesophagostomum and Necator being the most prevalent. We detected both host-specific and shared strongylid haplotypes. Skin-penetrating Necator gorillaehaplotypes were shared between humans and gorillas but Necator americanus were much more restricted to humans. Strongylid communities of local hunter-gatherers employed as trackers were more similar to those of gorillas compared to their relatives, who spent more time in villages. This was due to lower abundance of human-origin N. americanus in both gorillas and trackers. Habituated gorillas or those under habituation did not show larger overlap of strongylids with humans compared to unhabituated. We concluded that the occurrence of the human-specific strongylids in gorillas does not increase with direct contact between gorillas and humans due to the habituation. Overall, our results indicate that the degree of habitat sharing between hosts, together with mode of parasite transmission, are important factors for parasite spillover among primates.

The Ethiopian highlands represent a remarkable biodiversity ‘hot spot’ with a very high number of endemic species, even among vertebrates. Ethiopian representatives of a species complex of speckled brush‐furred rats (Lophuromys flavopunctatus sensu lato) inhabit highland habitats ranging from low‐elevation forests to Afroalpine grasslands. These may serve as a suitable model for understanding evolutionary processes leading to high genetic and ecological diversity in montane biodiversity hot spots. Here, we analyse the most comprehensive genetic data set of this group, comprising 315 specimens (all nine putative Ethiopian Lophuromys taxa sampled across most of their distribution ranges) genotyped at one mitochondrial and four nuclear markers, and thousands of SNPs from ddRAD sequencing. We performed phylogenetic analyses, delimited species and mapped their distribution and estimated divergence time between species (under the species‐tree framework) and mitochondrial lineages. We found significant incongruence between mitochondrial and nuclear phylogenies, most probably caused by multiple interspecific introgression events. We discuss alternative scenarios of Ethiopian Lophuromys evolution, from retention of ancestral polymorphism to hybridization upon secondary contact of partially reproductively isolated lineages leading to reticulate evolution. Finally, we use the diversity of the speckled brush‐furred rats for the description of the main biogeographic patterns in the fauna of the Ethiopian highlands.

Vertebrate gut microbiota (GM) is comprised of a taxonomically diverse consortium of symbiotic and commensal microorganisms that have a pronounced effect on host physiology, immune system function and health status. Despite much research on interactions between hosts and their GM, the factors affecting inter- and intraspecific GM variation in wild populations are still poorly known. We analysed data on faecal microbiota composition in 51 passerine species (319 individuals) using Illumina MiSeq sequencing of bacterial 16S rRNA (V3-V4 variable region). Despite pronounced interindividual variation, GM composition exhibited significant differences at the interspecific level, accounting for approximately 20%-30% of total GM variation. We also observed a significant correlation between GM composition divergence and host's phylogenetic divergence, with strength of correlation higher than that of GM vs. ecological or life history traits and geographic variation. The effect of host's phylogeny on GM composition was significant, even after statistical control for these confounding factors. Hence, our data do not support codiversification of GM and passerine phylogeny solely as a by-product of their ecological divergence. Furthermore, our findings do not support that GM vs. host's phylogeny codiversification is driven primarily through trans-generational GM transfer as the GM vs. phylogeny correlation does not increase with higher sequence similarity used when delimiting operational taxonomic units. Instead, we hypothesize that the GM vs. phylogeny correlation may arise as a consequence of interspecific divergence of genes that directly or indirectly modulate composition of GM.

Hantaviruses, well-known human pathogens, have only recently been identified on the African continent. Tigray virus (TIGV) was found in Ethiopia in 2012 in a Murinae species, Stenocephalemys albipes, but the genetic data obtained at that time were too limited to correctly assess its phylogenetic position within the hantavirus tree. We used high throughput sequencing to determine the complete genome of TIGV, which showed a typical hantavirus organisation. The large (L), medium (M), and small (S) genome segments were found to be 6532, 3594 and 1908 nucleotides long, respectively, and the 5' and 3' termini for all three segments were predicted to form the panhandle-like structure typical for bunyaviruses. Nucleotide-based phylogenetic analyses revealed that all three coding segments cluster in the phylogroup III sensu Guo et al. (2013). However, while TIGV S segment is basal to the Murinae-associated hantaviruses, the M and L segments are basal to the Soricomorpha-associated hantaviruses. TIGV is the first Murinae-borne hantavirus showing this inconsistent segmental clustering in the hantavirus phylogenetic tree. We finally propose non-exclusive scenarios that could explain the original phylogenetic position of TIGV. (C) 2016 Elsevier B.V. All rights reserved.

Strongylid nematodes in large terrestrial herbivores such as great apes, equids, elephants, and humans tend to occur in complex communities. However, identification of all species within strongylid communities using traditional methods based on coproscopy or single nematode amplification and sequencing is virtually impossible. High-throughput sequencing (HTS) technologies provide opportunities to generate large amounts of sequence data and enable analyses of samples containing a mixture of DNA from multiple species/genotypes. We designed and tested an HTS approach for strain-level identification of gastrointestinal strongylids using ITS-2 metabarcoding at the MiSeq Illumina platform in samples from two free-ranging non-human primate species inhabiting the same environment, but differing significantly in their host traits and ecology. Although we observed overlapping of particular haplotypes, overall the studied primate species differed in their strongylid nematode community composition. Using HTS, we revealed hidden diversity in the strongylid nematode communities in non-human primates, more than one haplotype was found in more than 90% of samples and coinfections of more than one putative species occurred in 80% of samples. In conclusion, the HTS approach on strongylid nematodes, preferably using fecal samples, represents a time and cost-efficient way of studying strongylid communities and provides a resolution superior to traditional approaches.

Murine cytomegalovirus (MCMV) has been reported from house mice (Mus musculus) worldwide, but only recently from Eastern house mice (M. m. musculus), of particular interest because they form a semi-permeable species barrier in Europe with Western house mice, M. m. domesticus. Here we report genome sequences of EastMCMV (from Eastern mice), and set these in the context of MCMV genomes from genus Mus hosts. We show EastMCMV and WestMCMV are genetically distinct. Phylogeny splitting analyses show a genome wide (94%) pattern consistent with no West-East introgression, the major exception (3.8%) being a genome-terminal region of duplicated genes involved in host immune system evasion. As expected from its function, this is a region of maintenance of ancestral polymorphism: The lack of clear splitting signal cannot be interpreted as evidence of introgression. The EastMCMV genome sequences reported here can therefore serve as a well-described resource for exploration of murid MCMV diversity.

Understanding the complex Entamoeba communities in the mammalian intestine has been, to date, complicated by the lack of a suitable approach for molecular detection of multiple variants co-occurring in mixed infections. Here, we report on the application of a high throughput sequencing approach based on partial 18S rDNA using the Illumina MiSeq platform. We describe, to our knowledge, for the first time, the Entamoeba communities in humans, free-ranging western lowland gorillas and central chimpanzees living in the Dja Faunal Reserve in Cameroon. We detected 36 Entamoeba haplotypes belonging to six haplotype clusters, containing haplotypes possessing high and low host specificity. Most of the detected haplotypes belonged to commensal Entamoeba, however, the pathogenic species (Entamoeba histolytica and Entamoeba nuttalli) were also detected. We observed that some Entamoeba haplotypes are shared between humans and other hosts, indicating their zoonotic potential. The findings are important not only for understanding the epidemiology of amoebiasis in humans in rural African localities, but also in the context of wild great ape conservation.

Captive-bred mallards Anas platyrhynchos have been released for hunting purposes at a very large scale in Europe since the mid-1970s. In spite of a potential genetic impact, the actual contribution of restocked mallards to the genome of the target population has received little attention. The genetic structure of modern wild mallards in the Camargue, Southern France, was assessed from two samples: one originating from shot birds in hunting bags and one from presumed wild ducks captured alive in a hunting-free reserve. Reference samples originated from five mallard farms, as well as from museum samples collected before the mid-1970s (i.e. before massive mallard releases started). Our results revealed that the genetic signature of wild wintering mallards has not changed significantly because museum and presumed wild samples from the Camargue hunting-free nature reserve were genetically similar, and clearly differentiated from the farm mallards. This suggests that mallard releases in the Camargue or elsewhere in France, although massive, have not actually translated into complete admixture of wild and captive genomes, most likely due to low survival of released birds once in the wild. Nevertheless, although genetic introgression of the wild population by captive-bred was contained, we found significant rates of hybridization between wild and captive-bred mallards in modern samples. This result suggests that long-term releases of captive-bred mallards, if carried on at such large scale, could compromise irreversibly the genetic structure and composition of European mallards. This work contributes to fill in the gap on the monitoring of the genetic consequences of large-scale game releases for exploitation.

The genetic integrity of natural populations can be threatened through large-scale introduction of farmed stocks with different genetic or geographic origin. Huge numbers of farm-reared mallard (Anas platyrhynchos, Anatidae) have been introduced into the wild in many European countries since 1970. Czech breeding facilities currently produce around 200–300,000 ducks annually, exceeding wild numbers by around 10 times. Such facilities, however, were founded with hybrid ducks from outside the Czech Republic. Three types of DNA markers, two neutral (14 microsatellite DNA loci and the mitochondrial DNA control region) and one under selection (MHC class I locus), were used to genotype mallards from six Czech breeding facilities (n = 131) and seven wild nesting localities (n = 139). We found marked genetic divergence between wild and captive-bred populations, the latter having significantly lower genetic diversity. Released captive-bred mallards were integrated into breeding wild population through hybridization mediated by high frequency nesting. Overall, our data suggest that release of captive-bred individuals threatens the genetic integrity of wild population. Massive restocking may also be undesirable as regards public health. Waterfowl are known reservoirs of transmittable pathogens and large-scale restocking could alter immune defence gene frequencies in wild population. We propose the establishment of a national genetic monitoring programme for breeding facilities.

PCR diagnostics detected 100% prevalence of Helicobacter in 425 wild house mice ( Mus musculus ) from across central Europe. Of seven species identified, the five most frequent were Helicobacter rodentium (78%), H. typhlonius (53%), H. hepaticus (41%), H. bilis (30%), and H. muridarum (1%). Double infections were more common (42%) than single (30%) and triple (21%) infections. Wild house mice could be considered potential reservoirs of Helicobacter strains for both humans and other vertebrates.

Syphacia obvelata is a common gastro-intestinal parasitic nematode of the house mouse (Mus musculus), a prime model rodent species. Investigations of the genetic structure, variability of parasite populations and other biological aspects of this host-parasite system are limited due to the lack of genetic resources for S. obvelata. To fill this gap, we developed a set of microsatellite markers for S. obvelata, using a 454 pyrosequencing approach. We designed three multiplex panels allowing genotyping of 10 polymorphic loci and scrutinized them on 42 samples from two different regions inhabited by two different house mouse subspecies (Mus musculus musculus and M. m. domesticus). The numbers of alleles ranged from 2 to 6 with mean observed heterozygosities 0.1476 and 0.2095 for domesticus and musculus worms, respectively. The described markers will facilitate further studies on population biology and co-evolution of this host-parasite system. (C) 2012 Elsevier B.V. All rights reserved.

In this study, we present 10 new microsatellite markers for a carnivorous bat, Megaderma lyra (family Megadermatidae). Out of 10 loci, only one locus showed a statistically significant value (P < 0.05) for FIS coefficient, and one was found to be X-linked. These loci were tested on 189 individuals representing four populations from southern India, and all loci were highly polymorphic. The mean number of observed alleles per locus was 10.1 (range from 6 to 16). Expected heterozygosity values ranged from 0.54 to 0.91. These polymorphic markers will be used to examine population structure, mating and dispersal behaviour of M. lyra in the field.

The effects of gastrointestinal tract microbiota (GTM) on host physiology and health have been the subject of considerable interest in recent years. While a variety of captive bred species have been used in experiments, the extent to which GTM of captive and/or inbred individuals resembles natural composition and variation in wild populations is poorly understood. Using 454 pyrosequencing, we performed 16S rDNA GTM barcoding for 30 wild house mice (Mus musculus) and wild-derived inbred strain mice belonging to two subspecies (M. m. musculus and M. m. domesticus). Sequenced individuals were selected according to a 2 × 2 experimental design: wild (14) vs. inbred origin (16) and M. m. musculus (15) vs. M. m. domesticus (15). We compared alpha diversity (i.e. number of operational taxonomic units - OTUs), beta diversity (i.e. interindividual variability) and microbiota composition across the four groups. We found no difference between M. m. musculus and M. m. domesticus subspecies, suggesting low effect of genetic differentiation between these two subspecies on GTM structure. Both inbred and wild populations showed the same level of microbial alpha and beta diversity; however, we found strong differentiation in microbiota composition between wild and inbred populations. Relative abundance of ~ 16% of OTUs differed significantly between wild and inbred individuals. As laboratory mice represent the most abundant model for studying the effects of gut microbiota on host metabolism, immunity and neurology, we suggest that the distinctness of laboratory-kept mouse microbiota, which differs from wild mouse microbiota, needs to be considered in future biomedical research.

P. Abelló ...et al.-- 3 pages, 1 table, This article documents the addition of 123 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Brenthis ino, Cichla orinocensis, Cichla temensis, Epinephelus striatus, Gobio gobio, Liocarcinus depurator, Macrolophus pygmaeus, Monilinia vaccinii-corymbosi, Pelochelys cantorii, Philotrypesis josephi, Romanogobio vladykovi, Takydromus luyeanus and Takydromus viridipunctatus. These loci were cross-tested on the following species: Cichla intermedia, Cichla ocellaris, Cichla pinima, Epinephelus acanthistius, Gobio carpathicus, Gobio obtusirostris, Gobio sp. 1, Gobio volgensis, Macrolophus costalis, Macrolophus melanotoma, Macrolophus pygmaeus, Romanogobio albipinnatus, Romanogobio banaticus, Romanogobio belingi, Romanogobio kesslerii, Romanogobio parvus, Romanogobio pentatrichus, Romanogobio uranoscopus, Takydromus formosanus, Takydromus hsuehshanesis and Takydromus stejnegeri

Ten polymorphic microsatellites for the common pochard were isolated from microsatellite enriched libraries. Seven microsatellites were obtained based on cross-species amplification. These 17 microsatellites exhibited polymorphism within a population sample, with numbers of alleles per locus ranging from 2 to 17 and expected heterozygosities from 0.053 to 0.916. These markers will be helpful for studying reproductive strategies in the pochard, as well as in addressing population genetic questions. Four of the microsatellite loci showed significant departures from Hardy–Weinberg equilibrium and/or high estimates of null allele frequencies in our study population. Hence, they should be used with caution.

The mammalian major histocompatibility complex (MHC) is a tightly linked cluster of immune genes, and is often thought of as inherited as a unit. This has led to the hope that studying a single MHC gene will reveal patterns of evolution representative of the MHC as a whole. In this study we analyse a 1000-km transect of MHC variation traversing the European house mouse hybrid zone to compare signals of selection and patterns of diversification at two closely linked MHC class II genes, H-2Aa and H-2Eb. We show that although they are 0.01 cM apart (that is, recombination is expected only once in 10 000 meioses), disparate evolutionary patterns were detected. H-2Aa shows higher allelic polymorphism, faster allelic turnover due to higher mutation rates, stronger positive selection at antigen-binding sites and higher population structuring than H-2Eb. H-2Eb alleles are maintained in the gene pool for longer, including over separation of the subspecies, some H-2Eb alleles are positively and others negatively selected and some of the alleles are not expressed. We conclude that studies on MHC genes in wild-living vertebrates can give substantially different results depending on the MHC gene examined and that the level of polymorphism in a related species is a poor criterion for gene choice.

Effects of vertebrate-associated microbiota on physiology and health are of significant interest in current biological research. Most previous studies have focused on host-microbiota interactions in captive-bred mammalian models. These interactions and their outcomes are still relatively understudied, however, in wild populations and non-mammalian taxa. Using deep pyrosequencing, we described the cloacal microbiome (CM) composition in free living barn swallows Hirundo rustica, a long-distance migratory passerine bird. Barn swallow CM was dominated by bacteria of the Actinobacteria, Proteobacteria and Firmicutes phyla. Bacteroidetes, which represent an important proportion of the digestive tract microbiome in many vertebrate species, was relatively rare in barn swallow CM (< 5%). CM composition did not differ between males and females. A significant correlation of CM within breeding pair members is consistent with the hypothesis that cloacal contact during within-pair copulation may promote transfer of bacterial assemblages. This effect on CM composition had a relatively low effect size, however, possibly due to the species' high level of sexual promiscuity.