Richard P. Brown

The vaginal microbiome during pregnancy and the postpartum period in a European population

The composition and structure of the pregnancy vaginal microbiome may influence susceptibility to adverse pregnancy outcomes. Studies on the pregnant vaginal microbiome have largely been limited to Northern American populations. Using MiSeq sequencing of 16S rRNA gene amplicons, we characterised the vaginal microbiota of a mixed British cohort of women (n = 42) who experienced uncomplicated term delivery and who were sampled longitudinally throughout pregnancy (8–12, 20–22, 28–30 and 34–36 weeks gestation) and 6 weeks postpartum. We show that vaginal microbiome composition dramatically changes postpartum to become less Lactobacillus spp. dominant with increased alpha-diversity irrespective of the community structure during pregnancy and independent of ethnicity. While the pregnancy vaginal microbiome was characteristically dominated by Lactobacillus spp. and low alpha-diversity, unlike Northern American populations, a significant number of pregnant women this British population had a L. jensenii-dominated microbiome characterised by low alpha-diversity. L. jensenii was predominantly observed in women of Asian and Caucasian ethnicity whereas L. gasseri was absent in samples from Black women. This study reveals new insights into biogeographical and ethnic effects upon the pregnancy and postpartum vaginal microbiome and has important implications for future studies exploring relationships between the vaginal microbiome, host health and pregnancy outcomes.

Phylogeography of skinks (Chalcides) in the Canary Islands inferred from mitochondrial DNA sequences

Mitochondrial DNA (mtDNA) evolution was investigated in skinks of the genus Chalcides found in the Canary Islands (Ch. sexlineatus, Ch. viridanus and Ch. simonyi), together with some North African congenerics (Ch. polylepis and Ch. mionecton). Several sites were included within islands to cover areas of known within‐island geographical variation in morphology. Skinks from the islands of El Hierro and La Gomera appear to be sister taxa. The relationships between this clade and the Tenerife and Gran Canarian skinks were not fully resolved, although the best working hypothesis indicated monophyly with the former, with the latter forming a closely related outgroup. Ch. simonyi from Fuerteventura was more distantly related to the Western Canary Island skinks and did not show close relationships with the North African species Ch. mionecton and Ch. polylepis. Possible colonization sequences for the four most Western Canary Islands were considered. El Hierro appears to have been colonized relatively recently from La Gomera, commensurate with the recent origin of this island, while dispersal between La Gomera and Tenerife and between Gran Canaria and Tenerife or La Gomera appears to have taken place considerably earlier. Substantial within‐island haplotype divergence was found in Gran Canaria and Tenerife. This may be a result of recent periods of intense volcanic activity found within these two islands. Lower levels of within‐island differentiation are found in La Gomera and El Hierro and may be explained by lower levels of volcanic activity during recent geological history and a more recent colonization, respectively.

Divergence times and colonization of the Canary Islands by Gallotia lizards

The Canary Islands have become a model region for evolutionary studies. We obtained 1.8 Kbp of mtDNA sequence from all known island forms of the endemic lizard genus Gallotia and from its sister taxon Psammodromus in order to reanalyze phylogenetic relationships within the archipelago, estimate lineage divergence times, and reconstruct the colonization history of this group. Well-supported phylogenies were obtained using maximum parsimony and Bayesian inference. Previous studies have been unable to establish the branching pattern at the base of the tree. We found evidence that G. stehlini (Gran Canaria) originated from the most basal Gallotia node and G. atlantica from the subsequent node. Divergence times were estimated under a global clock using Bayesian Markov Chain Monte Carlo methods implemented by three different programs: BEAST, MCMCTREE, MULTIDIVTIME. Node constraints were derived from subaerial island appearance data and were incorporated into the analyses as soft or hard maximal bounds. Posterior node ages differed slightly between programs, possibly due to different priors on divergence times. The most eastern Canary Islands first emerged just over 20 mya and their colonization appears to have taken place relatively quickly, around 17-20 mya. The subsequent node is consistent with cladogenesis due to colonization of Gran Canaria from the eastern islands about 11-13 mya. The western islands appear to have been colonized by a dispersal event from Lanzarote/Fuerteventura in the east to either La Gomera or one of the ancient edifices that subsequently formed Tenerife in the west, about 9-10 mya. Within the western islands, the most recent node that is ancestral to both the G. intermedia/G. gomerana/G. simonyi and the G.galloti/G. caesaris clades is dated at about 5-6 mya. Subsequent dispersal events between ancient Tenerife islands and La Gomera are dated at around 3 mya in both clades, although the direction of dispersal cannot be determined. Finally, we show that G. galloti is likely to have colonized La Palma more than 0.5 Ma after emergence of the island 1.77 mya, while G. caesaris from the same clade may have colonized El Hierro very soon after it emerged 1.12 mya. There are tentative indications that the large-bodied endangered G. simonyi colonized El Hierro around the same time or even later than the smaller-bodied G. caesaris. This study demonstrates the effectiveness of Bayesian dating of a phylogeny in helping reconstruct the historical pattern of dispersal across an oceanic archipelago.

The Atlas mountains as a biogeographical divide in North-West Africa: evidence from mtDNA evolution in the agamid lizard Agama impalearis

Since the early Miocene there have been several physical events within NW Africa that are likely to have had a major impact on its faunal diversity. Phylogeographical studies will shed new light on the biogeography of the region. We analysed mitochondrial DNA diversity in the agamid lizard Agama impalearis (also called A. bibronii) based on sequences from mitochondrial genes with very different evolutionary rates (16S rRNA and ND2). Well-supported topologies of rooted maximum parsimony trees (with a Laudakia outgroup) and unrooted haplotype networks indicated two major clades with similar branch lengths. These clades have non-overlapping distributions representing respective areas to the North and West and South and East of the Atlas mountain chain and each could be given full species recognition. Nested clade analyses indicate that historical and possible present-day allopatry account for the primary phylogeographic pattern. Further evidence is provided by the estimated timing of cladogenesis, based on calibration of evolutionary rates in the ND2 gene of another continental Agamid. Sequence divergence between clades corresponds to 8.5-9.4mya, coinciding with the main period of orogenic uplift of the Atlas. Additional evidence of cladogenesis by allopatric fragmentation is also detected within the North/West Atlas clade, although contiguous range expansion is the most predominant explanation of more recent phylogeographic effects in this species. Miocene vicariance mediated by the Atlas may provide a general explanation of intra- and interspecific biogeographical patterns in NW African species.

Rate variation and estimation of divergence times using strict and relaxed clocks

BackgroundUnderstanding causes of biological diversity may be greatly enhanced by knowledge of divergence times. Strict and relaxed clock models are used in Bayesian estimation of divergence times. We examined whether: i) strict clock models are generally more appropriate in shallow phylogenies where rate variation is expected to be low, ii) the likelihood ratio test of the clock (LRT) reliably informs which model is appropriate for dating divergence times. Strict and relaxed models were used to analyse sequences simulated under different levels of rate variation. Published shallow phylogenies (Black bass, Primate-sucking lice, Podarcis lizards, Gallotiinae lizards, and Caprinae mammals) were also analysed to determine natural levels of rate variation relative to the performance of the different models.ResultsStrict clock analyses performed well on data simulated under the independent rates model when the standard deviation of log rate on branches, σ, was low (≤0.1), but were inappropriate when σ>0.1 (95% of rates fall within 0.0082-0.0121 subs/site/Ma when σ = 0.1, for a mean rate of 0.01). The independent rates relaxed clock model performed well at all levels of rate variation, although posterior intervals on times were significantly wider than for the strict clock. The strict clock is therefore superior when rate variation is low. The performance of a correlated rates relaxed clock model was similar to the strict clock. Increased numbers of independent loci led to slightly narrower posteriors under the relaxed clock while older root ages provided proportionately narrower posteriors. The LRT had low power for σ = 0.01-0.1, but high power for σ = 0.5-2.0. Posterior means of σ2 were useful for assessing rate variation in published datasets. Estimates of natural levels of rate variation ranged from 0.05-3.38 for different partitions. Differences in divergence times between relaxed and strict clock analyses were greater in two datasets with higher σ2 for one or more partitions, supporting the simulation results.ConclusionsThe strict clock can be superior for trees with shallow roots because of low levels of rate variation between branches. The LRT allows robust assessment of suitability of the clock model as does examination of posteriors on σ2.

Cladogenesis and phylogeography of the lizard Phrynocephalus vlangalii (Agamidae) on the Tibetan plateau

Phrynocephalus vlangalii is restricted to dry sand or Gobi desert highlands between major mountain ranges in the Qinghai (Tibetan) Plateau. Mitochondrial DNA (mtDNA) sequence (partial ND2, tRNATrp and partial tRNAAla) was obtained from 293 Phrynocephalus sampled from 34 sites across the plateau. Partitioned Bayesian and maximum parsimony phylogenetic analyses revealed that P. vlangalii and two other proposed species (P. erythrus and P. putjatia) together form a monophyletic mtDNA clade which, in contrast with previous studies, does not include P. theobaldi and P. zetangensis. The main P. vlangalli clade comprises seven well‐supported lineages that correspond to distinct geographical areas with little or no overlap, and share a most recent common ancestor at 5.06 ± 0.68 million years ago (mya). This is much older than intraspecific lineages in other Tibetan animal groups. Analyses of molecular variance indicated that most of the observed genetic variation occurred among populations/regions implying long‐term interruption of maternal gene flow. A combined approach based on tests of population expansion, estimation of node dates, and significance tests on clade areas indicated that phylogeographical structuring has been primarily shaped by three main periods of plateau uplift during the Pliocene and Pleistocene, specifically 3.4 mya, 2.5 mya and 1.7 mya. These periods corresponded to the appearance of several mountain ranges that formed physical barriers between lineages. Populations from the Qaidam Basin are shown to have undergone major demographic and range expansions in the early Pleistocene, consistent with colonization of areas previously covered by the huge Qaidam palaeolake, which desiccated at this time. The study represents one of the most detailed phylogeographical analyses of the Qinghai Plateau to date and shows how geological events have shaped current patterns of diversity.

Bayesian estimation of post-Messinian divergence times in Balearic island lizards

Phylogenetic relationships and timings of major cladogenesis events are investigated in the Balearic Island lizards Podarcislilfordi and P.pityusensis using 2675bp of mitochondrial and nuclear DNA sequences. Partitioned Bayesian and Maximum Parsimony analyses provided a well-resolved phylogeny with high node-support values. Bayesian MCMC estimation of node dates was investigated by comparing means of posterior distributions from different subsets of the sequence against the most robust analysis which used multiple partitions and allowed for rate heterogeneity among branches under a rate-drift model. Evolutionary rates were systematically underestimated and thus divergence times overestimated when sequences containing lower numbers of variable sites were used (based on ingroup node constraints). The following analyses allowed the best recovery of node times under the constant-rate (i.e., perfect clock) model: (i) all cytochrome b sequence (partitioned by codon position), (ii) cytochrome b (codon position 3 alone), (iii) NADH dehydrogenase (subunits 1 and 2; partitioned by codon position), (iv) cytochrome b and NADH dehydrogenase sequence together (six gene-codon partitions), (v) all unpartitioned sequence, (vi) a full multipartition analysis (nine partitions). Of these, only (iv) and (vi) performed well under the rate-drift model. These findings have significant implications for dating of recent divergence times in other taxa. The earliest P.lilfordi cladogenesis event (divergence of Menorcan populations), occurred before the end of the Pliocene, some 2.6Ma. Subsequent events led to a West Mallorcan lineage (2.0Ma ago), followed 1.2Ma ago by divergence of populations from the southern part of the Cabrera archipelago from a widely-distributed group from north Cabrera, northern and southern Mallorcan islets. Divergence within P.pityusensis is more recent with the main Ibiza and Formentera clades sharing a common ancestor at about 1.0Ma ago. Climatic and sea level changes are likely to have initiated cladogenesis, with lineages making secondary contact during periodic landbridge formation. This oscillating cross-archipelago pattern in which ancient divergence is followed by repeated contact resembles that seen between East-West refugia populations from mainland Europe.

Population density in Podarcis lilfordi (Squamata, Lacertidae), a lizard species endemic to small islets in the Balearic Islands (Spain)

The Balearic lizard, Podarcis lilfordi , is present in 43 insular populations in the Cabrera archipelago and around the coasts of Mallorca and Menorca islands (Spain). We studied lizard densities over the entire range of distribution, analyzing observed differences of density in relation to island area, habitat diversity, availability of resources, presence of predators, competitors and human disturbances. The density of the Balearic lizard varies from less than 35 to almost 8000 lizards ha –1 , with an average of around 1500 lizards ha –1 . In some very small islets we detected no more than 10 individuals. Using a subsample of nine coastal islets (Menorca) we did not find any significant correlation between ground arthropod biomass and lizard density. The combination of island area and its maximal altitude, its so-called biotic capacity, was also uncorrelated with lizard density. In addition, neither degree of island accessibility nor presence/absence of seagull breeding colonies, were able to explain lizard densities. Islands without ship rats ( Rattus rattus ) showed a significantly higher lizard density, but islands in which rat eradication programs were launched during the study period, showed lower densities than those with rats but no eradication actions. Genetic variability was significantly higher on bigger lizard populations, lacking a correlation with lizard densities. No single independent variable can explain density differences among populations under study. Our results are discussed in the light of available hypotheses on factors affecting population densities.

Geographic variation and population systematics: Distinguishing between ecogenetics and phylogenetics

Abstract A series of quantitative procedures (1 ‐ pattern of anagenesis in putative phylogentic trees; 2 ‐ ≪correlation≫ based tests comparing observed and hypothesized patterns: 3 ‐ congruence assessed by random re‐sampling) are tested against a series of ≪known≫ cases of ecogenetically or phylogenetically caused geographic variation to establish whether they can distinguish between these two causes. The former two procedures can distinguish between these causes in some circumstances, but the last procedure is not particularly useful for this purpose because independent morphological character systems can have congruent patterns of geographic variation even when they are ecogenetically caused.

Bayesian Dating of Shallow Phylogenies with a Relaxed Clock

Bayesian methods are increasingly being used to estimate divergence times without the restrictive assumption of a global clock. Little is known about their reliability for shallow phylogenies where DNA sequence divergence is low. We analyzed both simulated and real sequences to evaluate dating methods in phylogenies with mid-late Miocene roots. A large number of data sets (5000) with 10 taxa each were simulated under a rate-drift model for trees with 2 topologies (balanced or unbalanced) and with different sets of divergence times (characterized by long or short external branches). Data were analyzed using Bayesian Markov chain Monte Carlo methods in which the prior on divergence times was specified from a birth-death process with species sampling (BDS) or a Dirichlet distribution using the programs MCMCTREE and MULTIDIVTIME. The programs generally performed well on shallow phylogenies, but posterior mean node ages were biased and 95% posterior intervals included true ages in fewer than 95% of trees in some analyses. This typically occurred when the 95% prior interval did not include the true age and/or sequence lengths were </= 1 kbp. Widths of posterior intervals were also very dependent on the position of the calibrated node within the tree, particularly when sequences were short. Different divergence times priors within MCMCTREE, MULTIDIVTIME, and BEAST were used to analyze mitochondrial DNA data sets from a Bovid subfamily (the Caprinae) from Asian Laudakia and North African Chalcides lizards. Posterior divergence times were quite sensitive to different BDS priors but less sensitive to different Dirichlet priors. Our study demonstrates the impact of the prior on divergence times in shallow phylogenies and shows that 1) prior intervals on nodes should be assessed as a prerequisite to a dating analysis, 2) >or= 1 kbp of quite rapidly evolving sequence may be required to obtain accurate posterior means and usefully narrow posterior intervals.