Genotype Imputation Reference Panel Selection Using Maximal Phylogenetic Diversity

The recent dramatic cost reduction of next-generation sequencing technology enables investigators to assess most variants in the human genome to identify risk variants for complex diseases. However, sequencing large samples remains very expensive. For a study sample with existing genotype data, such as array data from genome-wide association studies, a cost-effective approach is to sequence a subset of the study sample and then to impute the rest of the study sample, using the sequenced subset as a reference panel. The use of such an internal reference panel identifies population-specific variants and avoids the problem of a substantial mismatch in ancestry background between the study population and the reference population. To efficiently select an internal panel, we introduce an idea of phylogenetic diversity from mathematical phylogenetics and comparative genomics. We propose the “most diverse reference panel”, defined as the subset with the maximal “phylogenetic diversity”, thereby incorporating individuals that span a diverse range of genotypes within the sample. Using data both from simulations and from the 1000 Genomes Project, we show that the most diverse reference panel can substantially improve the imputation accuracy compared to randomly selected reference panels, especially for the imputation of rare variants. The improvement in imputation accuracy holds across different marker densities, reference panel sizes, and lengths for the imputed segments. We thus propose a novel strategy for planning sequencing studies on samples with existing genotype data.     

Host-Parasite Incongruences in Rodent Eimeria Suggest Significant Role of Adaptation Rather than Cophylogeny in Maintenance of Host Specificity

The degree of host specificity, its phylogenetic conservativeness and origin are virtually unknown in Eimeria. This situation is largely due to the inadequate sample of eimerian molecular data available for reliable phylogenetic analyses. In this study, we extend the data set by adding 71 new sequences of coccidia infecting 16 small-mammal genera, mostly rodents. According to the respective feasibility of PCR gene amplification, the new samples are represented by one or more of the following genes: nuclear 18S rRNA, plastid ORF 470, and mitochondrial COI. Phylogenetic analyses of these sequences confirm the previous hypothesis that Eimeria, in its current morphology-based delimitation, is not a monophyletic group. Several samples of coccidia corresponding morphologically to other genera are scattered among the Eimeria lineages. More importantly, the distribution of eimerians from different hosts indicates that the clustering of eimerian species is influenced by their host specificity, but does not arise from a cophylogenetic/cospeciation process; while several clusters are specific to a particular host group, inner topologies within these clusters do not reflect host phylogeny. This observation suggests that the host specificity of Eimeria is caused by adaptive rather than cophylogenetic processes.     

Evolutionary Diversification of New Caledonian Araucaria

New Caledonia is a global biodiversity hotspot. Hypotheses for its biotic richness suggest either that the island is a ‘museum’ for an old Gondwana biota or alternatively it has developed following relatively recent long distance dispersal and in situ radiation. The conifer genus Araucaria (Araucariaceae) comprises 19 species globally with 13 endemic to this island. With a typically Gondwanan distribution, Araucaria is particularly well suited to testing alternative biogeographic hypotheses concerning the origins of New Caledonian biota. We derived phylogenetic estimates using 11 plastid and rDNA ITS2 sequence data for a complete sampling of Araucaria (including multiple accessions of each of the 13 New Caledonian Araucaria species). In addition, we developed a dataset comprising 4 plastid regions for a wider taxon sample to facilitate fossil based molecular dating. Following statistical analyses to identify a credible and internally consistent set of fossil constraints, divergence times estimated using a Bayesian relaxed clock approach were contrasted with geological scenarios to explore the biogeographic history of Araucaria. The phylogenetic data resolve relationships within Araucariaceae and among the main lineages in Araucaria, but provide limited resolution within the monophyletic New Caledonian species group. Divergence time estimates suggest a Late Cretaceous-Cenozoic radiation of extant Araucaria and a Neogene radiation of the New Caledonian lineage. A molecular timescale for the evolution of Araucariaceae supports a relatively recent radiation, and suggests that earlier (pre-Cenozoic) fossil types assigned to Araucaria may have affinities elsewhere in Araucariaceae. While additional data will be required to adequately resolve relationships among the New Caledonian species, their recent origin is consistent with overwater dispersal following Eocene emersion of New Caledonia but is too old to support a single dispersal from Australia to Norfolk Island for the radiation of the Pacific Araucaria sect. Eutacta clade.     

Degrees of freedom in interspecific allometry: An adjustment for the effects of phylogenetic constraint

The data used in studies of bivariate interspecific allometry usually violate the assumption of statistical independence. Although the traits of each species are commonly treated as independent, the expression of a trait among species within a genus may covary because of shared common ancestry. The same effect exists for genera within a family and so on up the phylogenetic hierarchy. Determining sample size by counting data points overestimates the effective sample size, which then leads to overestimating the degrees of freedom that should be used in calculating probabilities and confidence intervals. This results in an inflated Type 1 error rate. Although some workers (e.g., Felsenstein [1985] Am. Nat. 125:1–15) have suggested that this issue may invalidate interspecific allometry as a comparative method, a correction for the problem can be approximated with variance components from a nested analysis of variance. Variance components partition the total variation in the data set among the levels of the nested hierarchy. If the variance component for each nested level is weighted by the number of groups at that level, the sum of these values is an estimate of an effective sample size for the data set which reflects the effects of phylogenetic constraint. Analysis of two data sets, using taxonomy to define levels of the nested hierarchy, suggests that it has been common for published studies of interspecific allometry to severely overestimate the number of degrees of freedom. Interspecific allometry remains an important comparative method for evaluating questions concerning individual species that are not similarly addressed by the format of most of the newer comparative methods. With the correction proposed here for estimating degrees of freedom, the major statistical weakness of the procedure is substantially reduced. © 1994 Wiley-Liss, Inc.     

Edge Principal Components and Squash Clustering: Using the Special Structure of Phylogenetic Placement Data for Sample Comparison

Principal components analysis (PCA) and hierarchical clustering are two of the most heavily used techniques for analyzing the differences between nucleic acid sequence samples taken from a given environment. They have led to many insights regarding the structure of microbial communities. We have developed two new complementary methods that leverage how this microbial community data sits on a phylogenetic tree. Edge principal components analysis enables the detection of important differences between samples that contain closely related taxa. Each principal component axis is a collection of signed weights on the edges of the phylogenetic tree, and these weights are easily visualized by a suitable thickening and coloring of the edges. Squash clustering outputs a (rooted) clustering tree in which each internal node corresponds to an appropriate “average” of the original samples at the leaves below the node. Moreover, the length of an edge is a suitably defined distance between the averaged samples associated with the two incident nodes, rather than the less interpretable average of distances produced by UPGMA, the most widely used hierarchical clustering method in this context. We present these methods and illustrate their use with data from the human microbiome.     

Metagenomic Sequencing of an In Vitro-Simulated Microbial Community

Background

Microbial life dominates the earth, but many species are difficult or even impossible to study under laboratory conditions. Sequencing DNA directly from the environment, a technique commonly referred to as metagenomics, is an important tool for cataloging microbial life. This culture-independent approach involves collecting samples that include microbes in them, extracting DNA from the samples, and sequencing the DNA. A sample may contain many different microorganisms, macroorganisms, and even free-floating environmental DNA. A fundamental challenge in metagenomics has been estimating the abundance of organisms in a sample based on the frequency with which the organism''s DNA was observed in reads generated via DNA sequencing.

Methodology/Principal Findings

We created mixtures of ten microbial species for which genome sequences are known. Each mixture contained an equal number of cells of each species. We then extracted DNA from the mixtures, sequenced the DNA, and measured the frequency with which genomic regions from each organism was observed in the sequenced DNA. We found that the observed frequency of reads mapping to each organism did not reflect the equal numbers of cells that were known to be included in each mixture. The relative organism abundances varied significantly depending on the DNA extraction and sequencing protocol utilized.

Conclusions/Significance

We describe a new data resource for measuring the accuracy of metagenomic binning methods, created by in vitro-simulation of a metagenomic community. Our in vitro simulation can be used to complement previous in silico benchmark studies. In constructing a synthetic community and sequencing its metagenome, we encountered several sources of observation bias that likely affect most metagenomic experiments to date and present challenges for comparative metagenomic studies. DNA preparation methods have a particularly profound effect in our study, implying that samples prepared with different protocols are not suitable for comparative metagenomics.     

Application of the phylogenetic informativeness method to chloroplast markers: A test case of closely related species in tribe Hydrangeeae (Hydrangeaceae)

In evolutionary biology appropriate marker selection for the reconstruction of solid phylogenetic hypotheses is fundamental. One of the most challenging tasks addresses the appropriate choice of genomic regions in studies of closely related species. Robust phylogenetic frameworks are central to studies dealing with questions ranging from evolutionary and conservation biology, biogeography to plant breeding. Phylogenetic informativeness profiles provide a quantitative measure of the phylogenetic signal in markers and therefore a method for locus prioritization. The present work profiles phylogenetic informativeness of mostly non-coding chloroplast regions in an angiosperm lineage of closely related species: the popular ornamental tribe Hydrangeeae (Hydrangeaceae, Cornales, Asterids). A recent phylogenetic study denoted a case of resolution contrast between the two strongly supported clades within tribe Hydrangeeae. We evaluate the phylogenetic signal of 13 highly variable plastid markers for estimating relationships within and among the currently recognized monophyletic groups of this tribe. A selection of combined loci based on their phylogenetic informativeness retrieved more robust phylogenetic hypotheses than simply combining individual markers performing best with respect to resolution, nodal support and accuracy or those presenting the highest number of parsimony informative characters. We propose the rpl32–ndhF intergenic spacer (IGS), trnV–ndhC IGS, trnL–rpl32 IGS, psbT–petB region and ndhA intron as the best candidates for future phylogenetic studies in Hydrangeeae and potentially in other Asterids. We also contrasted the phylogenetic informativeness of coded indels against substitutions concluding that, despite their low phylogenetic informativeness, coded indels provide additional phylogenetic signal that is nearly free of noise. Phylogenetic relationships obtained from our total combined analyses showed improved resolution and nodal support with respect to recently published results.     

Double‐digest RADseq loci using standard Illumina indexes improve deep and shallow phylogenetic resolution of Lophodermium,a widespread fungal endophyte of pine needles

The phylogenetic and population genetic structure of symbiotic microorganisms may correlate with important ecological traits that can be difficult to directly measure, such as host preferences or dispersal rates. This study develops and tests a low‐cost double‐digest restriction site‐associated DNA sequencing (ddRADseq) protocol to reveal among‐ and within‐species genetic structure for Lophodermium, a genus of fungal endophytes whose evolutionary analyses have been limited by the scarcity of informative markers. The protocol avoids expensive barcoded adapters and incorporates universal indexes for multiplexing. We tested for reproducibility and functionality by comparing shared loci from sample replicates and assessed the effects of numbers of ambiguous sites and clustering thresholds on coverage depths, number of shared loci among samples, and phylogenetic reconstruction. Errors between technical replicates were minimal. Relaxing the quality‐filtering criteria increased the mean coverage depth per locus and the number of loci recovered within a sample, but had little effect on the number of shared loci across samples. Increasing clustering threshold decreased the mean coverage depth per cluster and increased the number of loci recovered within a sample but also decreased the number of shared loci across samples, especially among distantly related species. The combination of low similarity clustering (70%) and relaxed quality‐filtering (allowing up to 30 ambiguous sites per read) performed the best in phylogenetic analyses at both recent and deep genetic divergences. Hence, this method generated sufficient number of shared homologous loci to investigate the evolutionary relationships among divergent fungal lineages with small haploid genomes. The greater genetic resolution also revealed new structure within species that correlated with ecological traits, providing valuable insights into their cryptic life histories.     

Partitioned Gene-Tree Analyses and Gene-Based Topology Testing Help Resolve Incongruence in a Phylogenomic Study of Host-Specialist Bees (Apidae: Eucerinae)

Incongruence among phylogenetic results has become a common occurrence in analyses of genome-scale data sets. Incongruence originates from uncertainty in underlying evolutionary processes (e.g., incomplete lineage sorting) and from difficulties in determining the best analytical approaches for each situation. To overcome these difficulties, more studies are needed that identify incongruences and demonstrate practical ways to confidently resolve them. Here, we present results of a phylogenomic study based on the analysis 197 taxa and 2,526 ultraconserved element (UCE) loci. We investigate evolutionary relationships of Eucerinae, a diverse subfamily of apid bees (relatives of honey bees and bumble bees) with >1,200 species. We sampled representatives of all tribes within the group and >80% of genera, including two mysterious South American genera, Chilimalopsis and Teratognatha. Initial analysis of the UCE data revealed two conflicting hypotheses for relationships among tribes. To resolve the incongruence, we tested concatenation and species tree approaches and used a variety of additional strategies including locus filtering, partitioned gene-trees searches, and gene-based topological tests. We show that within-locus partitioning improves gene tree and subsequent species-tree estimation, and that this approach, confidently resolves the incongruence observed in our data set. After exploring our proposed analytical strategy on eucerine bees, we validated its efficacy to resolve hard phylogenetic problems by implementing it on a published UCE data set of Adephaga (Insecta: Coleoptera). Our results provide a robust phylogenetic hypothesis for Eucerinae and demonstrate a practical strategy for resolving incongruence in other phylogenomic data sets.     

Phylobetadiversity among Forest Types in the Brazilian Atlantic Forest Complex

Phylobetadiversity is defined as the phylogenetic resemblance between communities or biomes. Analyzing phylobetadiversity patterns among different vegetation physiognomies within a single biome is crucial to understand the historical affinities between them. Based on the widely accepted idea that different forest physiognomies within the Southern Brazilian Atlantic Forest constitute different facies of a single biome, we hypothesize that more recent phylogenetic nodes should drive phylobetadiversity gradients between the different forest types within the Atlantic Forest, as the phylogenetic divergence among those forest types is biogeographically recent. We compiled information from 206 checklists describing the occurrence of shrub/tree species across three different forest physiognomies within the Southern Brazilian Atlantic Forest (Dense, Mixed and Seasonal forests). We analyzed intra-site phylogenetic structure (phylogenetic diversity, net relatedness index and nearest taxon index) and phylobetadiversity between plots located at different forest types, using five different methods differing in sensitivity to either basal or terminal nodes (phylogenetic fuzzy weighting, COMDIST, COMDISTNT, UniFrac and Rao’s H). Mixed forests showed higher phylogenetic diversity and overdispersion than the other forest types. Furthermore, all forest types differed from each other in relation phylobetadiversity patterns, particularly when phylobetadiversity methods more sensitive to terminal nodes were employed. Mixed forests tended to show higher phylogenetic differentiation to Dense and Seasonal forests than these latter from each other. The higher phylogenetic diversity and phylobetadiversity levels found in Mixed forests when compared to the others likely result from the biogeographical origin of several taxa occurring in these forests. On one hand, Mixed forests shelter several temperate taxa, like the conifers Araucaria and Podocarpus. On the other hand, tropical groups, like Myrtaceae, are also very representative of this forest type. We point out to the need of more attention to Mixed forests as a conservation target within the Brazilian Atlantic Forest given their high phylogenetic uniqueness.     

Genetic analyses reveal further cryptic lineages within the Myotis nattereri species complex

In recent years, new cryptic mammalian species have been discovered in Europe, many of which belong to the order Chiroptera. Within this order, some species such as Myotis nattereri contain several cryptic lineages/species, especially in the Mediterranean region. Here we present genetic, phylogenetic and morphological analysis on the Myotis nattereri species complex, focusing on France which is thought to be a contact zone for the different lineages/species. We sequenced the full Cytochrome b gene from individuals from 23 localities and investigated diagnostic morphological characteristics. Our results reveal new phylogenetic relationships within the Myotis nattereri species complex and among closely related species. We discuss morphological characters which are synapomorphic to each lineage and can be used to differentiate them. Our results also demonstrate the presence of a new lineage within the Myotis nattereri species complex. This lineage, endemic to Corsica, possibly represents a new cryptic species for which we present preliminary ecological data. We further identify the presence of three lineages/species in France and detail their distribution with potential contact zones.     

A review of the correlation of tergites, sternites, and leg pairs in diplopods

Background

Poorly preserved biological tissues have become an important source of DNA for a wide range of zoological studies. Measuring the quality of DNA obtained from these samples is often desired; however, there are no widely used techniques available for quantifying damage in highly degraded DNA samples. We present a general method that can be used to determine the frequency of polymerase blocking DNA damage in specific gene-regions in such samples. The approach uses quantitative PCR to measure the amount of DNA present at several fragment sizes within a sample. According to a model of random degradation the amount of available template will decline exponentially with increasing fragment size in damaged samples, and the frequency of DNA damage (λ) can be estimated by determining the rate of decline.

Results

The method is illustrated through the analysis of DNA extracted from sea lion faecal samples. Faeces contain a complex mixture of DNA from several sources and different components are expected to be differentially degraded. We estimated the frequency of DNA damage in both predator and prey DNA within individual faecal samples. The distribution of fragment lengths for each target fit well with the assumption of a random degradation process and, in keeping with our expectations, the estimated frequency of damage was always less in predator DNA than in prey DNA within the same sample (mean λ_predator = 0.0106 per nucleotide; mean λ_prey = 0.0176 per nucleotide). This study is the first to explicitly define the amount of template damage in any DNA extracted from faeces and the first to quantify the amount of predator and prey DNA present within individual faecal samples.

Conclusion

We present an approach for characterizing mixed, highly degraded PCR templates such as those often encountered in ecological studies using non-invasive samples as a source of DNA, wildlife forensics investigations and ancient DNA research. This method will allow researchers to measure template quality in order to evaluate alternate sources of DNA, different methods of sample preservation and different DNA extraction protocols. The technique could also be applied to study the process of DNA decay.     

Confidence Intervals for Population Allele Frequencies: The General Case of Sampling from a Finite Diploid Population of Any Size

The estimation of population allele frequencies using sample data forms a central component of studies in population genetics. These estimates can be used to test hypotheses on the evolutionary processes governing changes in genetic variation among populations. However, existing studies frequently do not account for sampling uncertainty in these estimates, thus compromising their utility. Incorporation of this uncertainty has been hindered by the lack of a method for constructing confidence intervals containing the population allele frequencies, for the general case of sampling from a finite diploid population of any size. In this study, we address this important knowledge gap by presenting a rigorous mathematical method to construct such confidence intervals. For a range of scenarios, the method is used to demonstrate that for a particular allele, in order to obtain accurate estimates within 0.05 of the population allele frequency with high probability (%), a sample size of is often required. This analysis is augmented by an application of the method to empirical sample allele frequency data for two populations of the checkerspot butterfly (Melitaea cinxia L.), occupying meadows in Finland. For each population, the method is used to derive % confidence intervals for the population frequencies of three alleles. These intervals are then used to construct two joint % confidence regions, one for the set of three frequencies for each population. These regions are then used to derive a % confidence interval for Jost''s D, a measure of genetic differentiation between the two populations. Overall, the results demonstrate the practical utility of the method with respect to informing sampling design and accounting for sampling uncertainty in studies of population genetics, important for scientific hypothesis-testing and also for risk-based natural resource management.     

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom‐feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre‐ and post‐mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi‐locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.     

The First Mitogenome of the Cyprus Mouflon (Ovis gmelini ophion): New Insights into the Phylogeny of the Genus Ovis

Sheep are thought to have been one of the first livestock to be domesticated in the Near East, thus playing an important role in human history. The current whole mitochondrial genome phylogeny for the genus Ovis is based on: the five main domestic haplogroups occurring among sheep (O. aries), along with molecular data from two wild European mouflons, three urials, and one argali. With the aim to shed some further light on the phylogenetic relationship within this genus, the first complete mitochondrial genome sequence of a Cypriot mouflon (O. gmelini ophion) is here reported. Phylogenetic analyses were performed using a dataset of whole Ovis mitogenomes as well as D-loop sequences. The concatenated sequence of 28 mitochondrial genes of one Cypriot mouflon, and the D-loop sequence of three Cypriot mouflons were compared to sequences obtained from samples representatives of the five domestic sheep haplogroups along with samples of the extant wild and feral sheep. The sample included also individuals from the Mediterranean islands of Sardinia and Corsica hosting remnants of the first wave of domestication that likely went then back to feral life. The divergence time between branches in the phylogenetic tree has been calculated using seven different calibration points by means of Bayesian and Maximum Likelihood inferences. Results suggest that urial (O. vignei) and argali (O. ammon) diverged from domestic sheep about 0.89 and 1.11 million years ago (MYA), respectively; and dates the earliest radiation of domestic sheep common ancestor at around 0.3 MYA. Additionally, our data suggest that the rise of the modern sheep haplogroups happened in the span of time between six and 32 thousand years ago (KYA). A close phylogenetic relationship between the Cypriot and the Anatolian mouflon carrying the X haplotype was detected. The genetic distance between this group and the other ovine haplogroups supports the hypothesis that it may be a new haplogroup never described before. Furthermore, the updated phylogenetic tree presented in this study determines a finer classification of ovine species and may help to classify more accurately new mitogenomes within the established haplogroups so far identified.     

Phylogenetic estimation of Mytilidae in the East China Sea inferred from mitochondrial genes and nuclear DNA sequence variation

We performed a phylogenetic estimation of the family Mytilidae in the East China Sea based on nuclear internal transcribed spacer (ITS) genes and two mitochondrial genes (COI and 16S RNA). Analysis of five mytilid species based on each of the three genes resulted in mostly congruent trees, although there were some discrepancies in the classification of these species. We combine the results obtained from the three separate analyses to provide a phylogenetic estimation of Mytilidae. We found that the Mytilidae was divided into two major lineages: in one clade, Mytilus galloprovincialis was grouped with Mytilus coruscus; in the second clade, Septifer bilocularis was placed at the basal position in an individual clade, and Perna viridis and Musculista senhousia were recovered as a monophyletic group. Although these finding provide important insights into the taxonomic relationships among the Mytilidae, many aspects of Mytilidae phylogeny remain unresolved. Further analysis based on more molecular information and extensive taxon sampling is necessary to elucidate the phylogenetic relationships among the major lineages within the Mytilidae.     

An ITS-based phylogenetic framework for the genus Vorticella: finding the molecular and morphological gaps in a taxonomically difficult group

Vorticella includes more than 100 currently recognized species and represents one of the most taxonomically challenging genera of ciliates. Molecular phylogenetic analysis of Vorticella has been performed so far with only sequences coding for small subunit ribosomal RNA (SSU rRNA); only a few of its species have been investigated using other genetic markers owing to a lack of similar sequences for comparison. Consequently, phylogenetic relationships within the genus remain unclear, and molecular discrimination between morphospecies is often difficult because most regions of the SSU rRNA gene are too highly conserved to be helpful. In this paper, we move molecular systematics for this group of ciliates to the infrageneric level by sequencing additional molecular markers—fast-evolving internal transcribed spacer (ITS) regions—in a broad sample of 66 individual samples of 28 morphospecies of Vorticella collected from Asia, North America and Europe. Our phylogenies all featured two strongly supported, highly divergent, paraphyletic clades (I, II) comprising the morphologically defined genus Vorticella. Three major lineages made up clade I, with a relatively well-resolved branching order in each one. The marked divergence of clade II from clade I confirms that the former should be recognized as a separate taxonomic unit as indicated by SSU rRNA phylogenies. We made the first attempt to elucidate relationships between species in clade II using both morphological and multi-gene approaches, and our data supported a close relationship between some morphospecies of Vorticella and Opisthonecta, indicating that relationships between species in the clade are far more complex than would be expected from their morphology. Different patterns of helix III of ITS2 secondary structure were clearly specific to clades and subclades of Vorticella and, therefore, may prove useful for resolving phylogenetic relationships in other groups of ciliates.     

Using multilocus approach to uncover cryptic diversity within Pseudotrapelus lizards from Saudi Arabia

The agamid Pseudotrapelus lizards inhabit the mountainous areas of the Arabian Peninsula and eastern North Africa. Currently six Pseudotrapelus species are recognised, though diagnostic morphological characters are still lacking, creating great difficulty in describing new species. Recently, two specimens of Pseudotrapelus were collected from the vicinity of Riyadh in central Saudi Arabia, an area that was not sampled in previous phylogenetic studies. In here we used both mitochondrial and nuclear data to investigate the phylogenetic position of the new samples, and assess their phylogenetic relationships with the other recognised species of Pseudotrapelus from across the distribution range of the genus. We used a multilocus approach of haplotype networks, concatenated datasets and species trees, performed mitochondrial and nuclear species delimitation analyses, and estimated divergence times. In general, our results support previous molecular studies and uncover the presence of cryptic diversity within Pseudotrapelus. The phylogenetic structure of the genus is of two major clades and within them seven distinct, delimited phylogenetic groups belonging to the six recognised species and the seventh to the individuals from Riyadh. The Riyadh specimens were distinct in all analyses performed. We suggest that the new specimens from the Riyadh area are a distinct lineage, forming a clade with their phylogenetic relatives, P. sinaitus and P. chlodnickii. The clade formed by these three species diverged during the Late Miocene around 6.4 Ma, with cladogenesis possibly facilitated by vicariance and isolation caused due to climatic fluctuations and the progression of sandy areas. Our results suggest further morphological research is necessary to revise the taxonomic status of this lineage and of the entire genus.     

Diversity of Diazotrophic Unicellular Cyanobacteria in the Tropical North Atlantic Ocean

We present data on the genetic diversity and phylogenetic affinities of N₂-fixing unicellular cyanobacteria in the plankton of the tropical North Atlantic Ocean. Our dinitrogenase gene (nifH) sequences grouped together with a group of cyanobacteria from the subtropical North Pacific; another subtropical North Pacific group was only distantly related. Most of the 16S ribosomal DNA sequences from our tropical North Atlantic samples were closely allied with sequences from a symbiont of the diatom Climacodium frauenfeldianum. These findings suggest a complex pattern of evolutionary and ecological divergence among unicellular cyanobacteria within and between ocean basins.