Search

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.

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.

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.

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.