Serena E. Dool

Phylogeography and postglacial recolonization of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers

The demographic history of Rhinolophus hipposideros (lesser horseshoe bat) was reconstructed across its European, North African and Middle‐Eastern distribution prior to, during and following the most recent glaciations by generating and analysing a multimarker data set. This data set consisted of an X‐linked nuclear intron (Bgn; 543 bp), mitochondrial DNA (cytb‐tRNA‐control region; 1630 bp) and eight variable microsatellite loci for up to 373 individuals from 86 localities. Using this data set of diverse markers, it was possible to determine the species’ demography at three temporal stages. Nuclear intron data revealed early colonization into Europe from the east, which pre‐dates the Quaternary glaciations. The mtDNA data supported multiple glacial refugia across the Mediterranean, the largest of which were found in the Ibero‐Maghreb region and an eastern location (Anatolia/Middle East)–that were used by R. hipposideros during the most recent glacial cycles. Finally, microsatellites provided the most recent information on these species’ movements since the Last Glacial Maximum and suggested that lineages that had diverged into glacial refugia, such as in the Ibero‐Maghreb region, have remained isolated. These findings should be used to inform future conservation management strategies for R. hipposideros and show the power of using a multimarker data set for phylogeographic studies.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera).

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

The Effects of Human-Mediated Habitat Fragmentation on a Sedentary Woodland-Associated Species (Rhinolophus hipposideros) at Its Range Margin

Among the many anthropogenic modifications to earths ecosystems, habitat loss and degradation pose the most immediate threat to many biota. The predicted consequences of fragmented habitats include lower species diversity, smaller population sizes, disrupted gene flow, increased drift and inbreeding and increased differentiation between neighbouring populations; all of which are thought to be further enhanced in species with low dispersal abilities. These factors, especially when occurring in tandem, can lead to an increased risk of extinction. To examine the genetic consequences of habitat fragmentation we selected an isolated population of a sedentary woodland specialist species (Rhinolophus hipposideros) to act as an indicator of disruptions to landscape level connectivity. Based on 491 individuals from 37 colonies our results revealed the presence of a broad North-Range/South-Range differentiation within this species in Ireland; a finding supported across datasets (mtDNA and nuclear microsatellites) and analyses. Analyses of echolocation data and microsatellites suggested further differentiation of the northern-most colonies. A landscape genetics framework to assess the impact of habitat versus geographic distance on population differentiation showed that habitat features (at a five km resolution) were equally likely to be correlated with differentiation as geographic distance considered alone. Further differentiation of the geographically disjunct groups is likely to occur in the future. The viability of either group alone is uncertain given their restricted distribution, small population sizes (based on census data and Ne estimates) and isolation. Roost provision and habitat restoration in the geographic region separating the differentiated groups will be fundamental to the recolonization of this area and the reestablishment of connectivity between the regional groups.

Mycobiomes of sympatric Amorphophallus albispathus (Araceae) and Camellia sinensis (Theaceae) – a case study reveals clear tissue preferences and differences in diversity and composition

Multiple biotic and abiotic parameters influence the dynamics of individual fungal species and entire communities. Major drivers for tropical plant endophytes are undoubtedly seasonality, local habitat conditions and biogeography. However, host specialization and tissue preferences also contribute to the structuring of endophytic mycobiomes. To elucidate such specializations and preferences, we sampled two commercially important, unrelated plant species, Amorphophallus albispathus and Camellia sinensis (tea plant) simultaneously at close proximity. The mycobiomes of different tissue types were assessed with high-throughput amplicon sequencing of the internal transcribed spacer DNA region. Both plants hosted different fungal communities and varied in α- and β-diversity, despite their neighboring occurrence. However, the fungal assemblages of Amorphophallus leaflets shared taxa with the mycobiomes of tea leaves, thereby suggesting common driving forces for leaf-inhabiting fungi irrespective of host plant identity. The mycobiome composition and diversity of tea leaves was clearly driven by leaf age. We suggest that the very youngest tea leaves are colonized by stochastic processes, while mycobiomes of old leaves are rather similar as the result of progressive succession. The biodiversity of fungi associated with A. albispathus was characterized by a large number of unclassified OTUs (at genus and species level) and by tissue-specific composition.This study is the first cultivation-independent high-throughput assessment of fungal biodiversity of an Amorphophallus species, and additionally expands the knowledge base on fungi associated with tea plants.

Variable Molecular Markers for the Order Mantophasmatodea (Insecta)

The recently discovered insect order Mantophasmatodea currently comprises 19 Southern African species. These mainly occur in allopatry, have high levels of color polymorphism and communicate via species- and gender-specific vibratory signals. High levels of interspecific morphological conservatism mean that cryptic species are likely to be uncovered. These aspects of Mantophasmatodean biology make them an ideal group in which to investigate population divergence due to habitat-specific adaptation, sexual selection, and potentially sensory speciation. Lack of appropriate genetic markers has thus far rendered such studies unfeasible. To address this need, the first microsatellite loci for this order were developed. Fifty polymorphic loci were designed specifically for Karoophasma biedouwense (Austrophasmatidae), out of which 23 were labeled and tested for amplification across the order using 2-3 individuals from 10 species, representing all 4 currently known families. A Bayesian mitochondrially encoded cytochrome c oxidase I (COI) topology was reconstructed and divergence dates within the order were estimated for the first time. Amplification success and levels of polymorphism were compared with genetic divergence and time since divergence. In agreement with studies on vertebrate taxa, both amplification and variability were negatively correlated with distance (temporal and genetic). The high number of informative loci will offer sufficient resolution for both broad level population genetic analysis and individual based pedigree or parentage analyses for most species in Austrophasmatidae, with at least some loci available for the other families. This resource will facilitate research into the evolutionary biology of this understudied but fascinating group.

Identification of the same polyomavirus species in different African horseshoe bat species is indicative of short-range host-switching events

Polyomaviruses (PyVs) are considered to be highly host-specific in different mammalian species, with no well-supported evidence for host-switching events. We examined the species diversity and host specificity of PyVs in horseshoe bats (Rhinolophus spp.), a broadly distributed and highly speciose mammalian genus. We annotated six PyV genomes, comprising four new PyV species, based on pairwise identity within the large T antigen (LTAg) coding region. Phylogenetic comparisons revealed two instances of highly related PyV species, one in each of the Alphapolyomavirus and Betapolyomavirus genera, present in different horseshoe bat host species (Rhinolophus blasii and R. simulator), suggestive of short-range host-switching events. The two pairs of Rhinolophus PyVs in different horseshoe bat host species were 99.9 and 88.8 % identical with each other over their respective LTAg coding sequences and thus constitute the same virus species. To corroborate the species identification of the bat hosts, we analysed mitochondrial cytb and a large nuclear intron dataset derived from six independent and neutrally evolving loci for bat taxa of interest. Bayesian estimates of the ages of the most recent common ancestors suggested that the near-identical and more distantly related PyV species diverged approximately 9.1E4 (5E3-2.8E5) and 9.9E6 (4E6-18E6) years before the present, respectively, in contrast to the divergence times of the bat host species: 12.4E6 (10.4E6-15.4E6). Our findings provide evidence that short-range host-switching of PyVs is possible in horseshoe bats, suggesting that PyV transmission between closely related mammalian species can occur.