Phylogeny Inference Research Articles

DNA Multiple Sequence Alignment Guided by Protein Domains: The MSA-PAD 2.0 Method.

Multiple sequence alignment (MSA) is a fundamental component in many DNA sequence analyses including metagenomics studies and phylogeny inference. When guided by protein profiles, DNA multiple alignments assume a higher precision and robustness. Here we present details of the use of the upgraded version of MSA-PAD (2.0), which is a DNA multiple sequence alignment framework able to align DNA sequences coding for single/multiple protein domains guided by PFAM or user-defined annotations. MSA-PAD has two alignment strategies, called "Gene" and "Genome," accounting for coding domains order and genomic rearrangements, respectively. Novel options were added to the present version, where the MSA can be guided by protein profiles provided by the user. This allows MSA-PAD 2.0 to run faster and to add custom protein profiles sometimes not present in PFAM database according to the user's interest. MSA-PAD 2.0 is currently freely available as a Web application at https://recasgateway.cloud.ba.infn.it/ .

Phylogenomics of Bartheletia paradoxa reveals its basal position in Agaricomycotina and that the early evolutionary history of basidiomycetes was rapid and probably not strictly bifurcating

The higher level phylogeny of fungi has been addressed in previous studies, but for those analyses, either taxon sampling or gene sampling was low, or some basal lineages important for the inference of basidiomycete phylogeny were lacking. Here, a phylogenomic analysis based on highly conserved genes and including the enigmatic species Bartheletia paradoxa from Ginkgo biloba is presented. While phylogenetic analyses including also less conserved parts of core eukaryotic genes yielded a basal position for the extremophile genus Wallemia with low support, an exclusion of highly variable parts of these genes suggested Bartheletia paradoxa as the most basal member of the Agaricomycotina, but again with low support. Network analyses suggest a network-like evolution at the base of the Basidiomycota, supported by phylogenies based on single genes and gene clusters with shared topology. When further removing noise by removing poorly resolving genes, strong but not maximum support was obtained for Bartheletia paradoxa being the sister lineage to all other Agaricomycotina. We speculate that the lack of support for the early splits in Agaricomycotina and Basidiomycota can probably be explained by rapid radiation, linked to major evolutionary developments, such as, in the case of Basidiomycota, the advent of basidia in the last common ancestor.

Non-Mendelian Single-Nucleotide Polymorphism Inheritance and Atypical Meiotic Configurations are Prevalent in Hop.

Hop ( L.) breeding programs seek to exploit genetic resources for bitter flavor, aroma, and disease resistance. However, these efforts have been thwarted by segregation distortion including female-biased sex ratios. To better understand the transmission genetics of hop, we genotyped 4512 worldwide accessions of hop, including cultivars, landraces, and over 100 wild accessions using a genotyping-by-sequencing (GBS) approach. From the resulting ∼1.2 million single-nucleotide polymorphisms (SNPs), prequalified GBS markers were validated by inferences in population structures and phylogeny. Analysis of pseudo-testcross (Pt) mapping data from F families revealed mixed patterns of Mendelian and non-Mendelian segregation. Three-dimensional (3D) cytogenetic analysis of late meiotic prophase nuclei from two wild and two cultivated hop revealed conspicuous and prevalent occurrences of multiple, atypical, nondisomic chromosome complexes including autosomes. We used genome-wide association studies (GWAS) and fixation index (F) analysis to demonstrate selection mapping of genetic loci for key traits including sex, bitter acids, and drought tolerance. Among the possible mechanisms underlying the observed segregation distortion from the genomic data analysis, the cytogenetic analysis points to meiotic chromosome behavior as one of the contributing factors. The findings shed light on long-standing questions on the unusual transmission genetics and phenotypic variation in hop, with major implications for breeding, cultivation, and the natural history of .

Genomics of adaptive divergence with chromosome-scale heterogeneity in crossover rate.

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome-scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low-crossover centres relative to the high-crossover peripheries of chromosomes ("Chromosome Centre-Biased Differentiation", CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome-scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome-scale heterogeneity in crossover rate in evolutionary genomics.

Phylogeny analysis from gene-order data with massive duplications

BackgroundGene order changes, under rearrangements, insertions, deletions and duplications, have been used as a new type of data source for phylogenetic reconstruction. Because these changes are rare compared to sequence mutations, they allow the inference of phylogeny further back in evolutionary time. There exist many computational methods for the reconstruction of gene-order phylogenies, including widely used maximum parsimonious methods and maximum likelihood methods. However, both methods face challenges in handling large genomes with many duplicated genes, especially in the presence of whole genome duplication.MethodsIn this paper, we present three simple yet powerful methods based on maximum-likelihood (ML) approaches that encode multiplicities of both gene adjacency and gene content information for phylogenetic reconstruction.ResultsExtensive experiments on simulated data sets show that our new method achieves the most accurate phylogenies compared to existing approaches. We also evaluate our method on real whole-genome data from eleven mammals. The package is publicly accessible at http://www.geneorder.org.ConclusionsOur new encoding schemes successfully incorporate the multiplicity information of gene adjacencies and gene content into an ML framework, and show promising results in reconstruct phylogenies for whole-genome data in the presence of massive duplications.

SiFit: inferring tumor trees from single-cell sequencing data under finite-sites models

Single-cell sequencing enables the inference of tumor phylogenies that provide insights on intra-tumor heterogeneity and evolutionary trajectories. Recently introduced methods perform this task under the infinite-sites assumption, violations of which, due to chromosomal deletions and loss of heterozygosity, necessitate the development of inference methods that utilize finite-sites models. We propose a statistical inference method for tumor phylogenies from noisy single-cell sequencing data under a finite-sites model. The performance of our method on synthetic and experimental data sets from two colorectal cancer patients to trace evolutionary lineages in primary and metastatic tumors suggests that employing a finite-sites model leads to improved inference of tumor phylogenies.

Genome-wide RAD sequencing data provide unprecedented resolution of the phylogeny of temperate bamboos (Poaceae: Bambusoideae)

The temperate bamboos (tribe Arundinarieae, Poaceae) are strongly supported as monophyly in recent molecular studies, but taxonomic delineation and phylogenetic relationships within the tribe lack resolution. Here, we sampled 39 species (36 temperate bamboos and 3 outgroups) for restriction-site associated DNA sequencing (RAD-seq) with an emphasis on Phyllostachys clade and related clades. Using the largest data matrix for the bamboos to date, we were able to infer phylogenetic relationships with unparalleled resolution. The Phyllostachys, Shibataea, and Arundinaria clades defined from plastid phylogeny, were not supported as monophyletic group. However, the RAD-seq phylogeny largely agreed with the morphology-based taxonomy, with two clades having leptomorph rhizomes strongly supported as monophyletic group. We also explored two approaches, BWA-GATK (a mapping system) and Stacks (a grouping system), for differences in SNP calling and phylogeny inference. For the same level of missing data, the BWA-GATK pipeline produced much more SNPs in comparison with Stacks. Phylogenetic analyses of the largest data matrices from both pipelines, using concatenation and coalescent methods provided similar tree topologies, despite the presence of missing data. Our study demonstrates the utility of RAD-seq data for elucidating phylogenetic relationships between genera and higher taxonomic levels in this important but phylogenetically challenging group.

Evolutionary analysis of Chironius snakes unveils cryptic diversity and provides clues to diversification in the Neotropics

Recent hypotheses to explain tropical diversity involves the Neogene and Quaternary geoclimatic dynamics, but the absence of unambiguous data permitting the choice between alternative hypotheses makes a general theory for the origin of tropical biodiversity far to be achieved. The occurrence of Chironius snakes in well-defined biogeographical regions led us to adopt Chironius as a model to unveil patterns of vertebrate diversification in the Neotropics. Here, we used molecular markers and records on geographic distribution to investigate Chironius evolution and, subsequently, providing hints on diversification in the Neotropics. To avoid analyzing nominal species that do not constitute exclusive evolutionary lineages, we firstly conducted a species delimitation study prior to carrying out the species distribution modeling analysis. We generated 161 sequences of 12S, 16S, c-mos and rag2 for 15 species and 50 specimens, and included additional data from GenBank yielding a matrix of 137 terminals, and performed the following evolutionary analyses: inference of a concatenated gene tree, estimation of gene divergence times, inference of the coalescent-based phylogeny of Chironius, estimation of effective population sizes and modeling potential distribution of species across the last millennia. We tested for species boundaries within Chironius by implementing a coalescent-based Bayesian species delimitation approach. Our analyses supported the monophyly of Chironius, although our findings underscored cryptic candidate species in C. flavolineatus and C. exoletus. The inferred timetree suggested that Chironius snakes have evolved in the early Miocene (ca. 20.2Mya) and began to diversify from the late Miocene to the early Pliocene, values that are much older than previously reported. Following genetic divergence of virtually all extant Chironius species investigated, the effective sizes of the populations have expanded when compared to their MRCAs. The evolutionary and SDM data from C. brazili and C. diamantina provided additional evidence of the origin of species in the Neotropics. We argue that temperature, instead of altitude, has been the major driving factor in the evolution of both species, and thus we present a case for the consequences of global warming.

Phase transition in the sample complexity of likelihood-based phylogeny inference

Reconstructing evolutionary trees from molecular sequence data is a fundamental problem in computational biology. Stochastic models of sequence evolution are closely related to spin systems that have been extensively studied in statistical physics and that connection has led to important insights on the theoretical properties of phylogenetic reconstruction algorithms as well as the development of new inference methods. Here, we study maximum likelihood, a classical statistical technique which is perhaps the most widely used in phylogenetic practice because of its superior empirical accuracy. At the theoretical level, except for its consistency, that is, the guarantee of eventual correct reconstruction as the size of the input data grows, much remains to be understood about the statistical properties of maximum likelihood in this context. In particular, the best bounds on the sample complexity or sequence-length requirement of maximum likelihood, that is, the amount of data required for correct reconstruction, are exponential in the number, n, of tips—far from known lower bounds based on information-theoretic arguments. Here we close the gap by proving a new upper bound on the sequence-length requirement of maximum likelihood that matches up to constants the known lower bound for some standard models of evolution. More specifically, for the r-state symmetric model of sequence evolution on a binary phylogeny with bounded edge lengths, we show that the sequence-length requirement behaves logarithmically in n when the expected amount of mutation per edge is below what is known as the Kesten-Stigum threshold. In general, the sequence-length requirement is polynomial in n. Our results imply moreover that the maximum likelihood estimator can be computed efficiently on randomly generated data provided sequences are as above. Our main technical contribution, which may be of independent interest, relates the total variation distance between the leaf state distributions of two trees with a notion of combinatorial distance between the trees. In words we show in a precise quantitative manner that the more different two evolutionary trees are, the easier it is to distinguish their output.

Striking pseudogenization in avian phylogenetics: Numts are large and common in falcons

Nuclear copies of mitochondrial genes (numts) are a well-known feature of eukaryotic genomes and a concern in systematics, as they can mislead phylogenetic inferences when inadvertently used. Studies on avian numts initially based on the chicken genome suggest that numts may be uncommon and relatively short among birds. Here we ask how common numts are in falcons, based on recently sequenced genomes of the Saker falcon (Falco cherrug) and Peregrine falcon (F. peregrinus). We identified numts by BLASTN searches and then extracted CYTB, ND2 and COI sequences from them, which were then used for phylogeny inference along with several sequences from other species in Falconiformes. Our results indicate that avian numts may be much more frequent and longer than previously thought. Phylogenetic inferences revealed multiple independent nuclear insertions throughout the history of the Falconiformes, including cases of sequences available in public databases and wrongly identified as authentic mtDNA. New sequencing technologies and ongoing efforts for whole genome sequencing will provide exciting opportunities for avian numt research in the near future.

A new distributed alignment-free approach to compare whole proteomes

Phylogeny inference has moved in recent years from the analysis of a single or few proteins to that of whole proteomes. However, the reconstruction of evolutionary trees for big number of species poses a significant computational challenge when using complete proteomes, even when relatively fast pairwise sequence comparison algorithms are used. We present a distributed approach that relies on the computation of distance measures based on maximal shared substrings within a bounded Hamming distance. The distributed system we built to implement this approach is flexible in that it supports a variety of design choices. It is based on the Spark framework and covers all the steps required by our approach, starting from the initial indexing of a set of FASTA sequences up to producing a report detailing the distances among these sequences, ranked according to a user-defined measure. Here we apply it to compare all proteins of selected organisms, divide them into groups and perform the comparisons within each group separately. The groups include: the functionally characterized proteins, the ribosomal proteins, and the unannotated proteins. We compute the average distances within the groups and evaluate their relationship and ability to capture the evolutionary closeness of organisms. We run experiments on selected species using a Hadoop computing cluster running Spark. The results show that the system implementing our approach is scalable and accurate.

P4P: a peptidome-based strain-level genome comparison web tool.

Peptidome similarity analysis enables researchers to gain insights into differential peptide profiles, providing a robust tool to discriminate strain-specific peptides, true intra-species differences among biological replicates or even microorganism-phenotype variations. However, no in silico peptide fingerprinting software existed to facilitate such phylogeny inference. Hence, we developed the Peptidomes for Phylogenies (P4P) web tool, which enables the survey of similarities between microbial proteomes and simplifies the process of obtaining new biological insights into their phylogeny. P4P can be used to analyze different peptide datasets, i.e. bacteria, viruses, eukaryotic species or even metaproteomes. Also, it is able to work with whole proteome datasets and experimental mass-to-charge lists originated from mass spectrometers. The ultimate aim is to generate a valid and manageable list of peptides that have phylogenetic signal and are potentially sample-specific. Sample-to-sample comparison is based on a consensus peak set matrix, which can be further submitted to phylogenetic analysis. P4P holds great potential for improving phylogenetic analyses in challenging taxonomic groups, biomarker identification or epidemiologic studies. Notably, P4P can be of interest for applications handling large proteomic datasets, which it is able to reduce to small matrices while maintaining high phylogenetic resolution. The web server is available at http://sing-group.org/p4p.

Species delimitation in the crypsis‐defended and polymorphic stick insects of the genus Libethra (Phasmatodea, Diapheromeridae)

Crypsis is a widespread defensive mechanism for prey ambushed by visually oriented predators. One well‐documented consequence of this interaction is the formation of search images by the predator that leads to the coexistence of multiple morphs in the prey population. Libethra is a genus of stick insects whose ornamentation and coloration are decidedly protective and exhibit exuberant intraspecific variation. Despite this fact, these traits have been extensively used to support the boundaries between the species it includes. This is the case in the north Andean species Libethra rabdota and Libethra rabdotula, delimited on the basis of a specific arrangement of ornamental characters. To state whether such taxonomical practice is meaningful and convenient, we evaluated the grouping of specimens of both species by means of the independent inference of phylogenies obtained with mitochondrial DNA sequences and ootaxonomy. This article synonymizes L. rabdota and L. rabdotula, redescribes the adults and describes the egg of the species. Additionally, we demonstrate that the characters related to crypsis used for the diagnosis of L. rabdota are not exclusive to this species by describing Libethra rioblanco sp. n. and Libethra ucumariensis sp. n. Our findings suggest that this genus requires extensive revision, but also provide support for the usage of different tools that might produce more accurate and stable classification. Evolutionary scenarios for both the extreme morphological variation and possible homoplastic state of camouflage‐related traits are discussed herein.

Live phylogeny with polytomies: Finding the most compact parsimonious trees

Construction of phylogenetic trees has traditionally focused on binary trees where all species appear on leaves, a problem for which numerous efficient solutions have been developed. Certain application domains though, such as viral evolution and transmission, paleontology, linguistics, and phylogenetic stemmatics, often require phylogeny inference that involves placing input species on ancestral tree nodes (live phylogeny), and polytomies. These requirements, despite their prevalence, lead to computationally harder algorithmic solutions and have been sparsely examined in the literature to date. In this article we prove some unique properties of most parsimonious live phylogenetic trees with polytomies, and their mapping to traditional binary phylogenetic trees. We show that our problem reduces to finding the most compact parsimonious tree for n species, and describe a novel efficient algorithm to find such trees without resorting to exhaustive enumeration of all possible tree topologies.

Abstract P1-06-06: Evidence for tumor heterogeneity and clonal evolution during invasive progression of breast cancer

Abstract Purpose: Intratumoral heterogeneity is well recognized to be an important driver of treatment resistance and metastasis. We undertook this N of one study to measure the degree of heterogeneity in a single large preinvasive lesion with an invasive component to determine the relationship between tumor heterogeneity, spatial distribution, clonal evolution, and invasive progression. Methods: We identified a patient with extensive DCIS measuring 7.5 cm, associated with 1.5 cm of an invasive component. We segregated the tumor sample into 32 unique blocks with precise geospatial localization; invasive cancer was identified in 3 of 32 blocks. NGS libraries were made from FFPE derived DNA (20-40ng) for full exome sequencing and hybridization to a 4.8 million element SNP array. All data were analyzed and a phylogentic tree was constructed. Results: The sequence data was analyzed with Platypus[1] and 3922 somatic mutation sites were found in total. These sites were concatenated into one sequence for each sample. Then a Neighbor-Joining tree was built with a Jukes-Cantor model and 1000 bootstrap replicates using FastME 2.0[2], to assess the reliability of the tree (Figure 1). Phylogenetic analysis revealed that invasive cancer evolved twice, independently. Dense sampling allowed reconstruction of the temporal order of mutations that accumulated in the cell lineage of the invasive cancers. Furthermore, the phylogeny revealed that distant regions may be closely genetically related, the oldest parts of the tumor were in the interior of the tumor, and that the invasive tumors evolved near the oldest parts of the tumor, rather than at the expanding front. Conclusions: Extensive sampling and sequencing of a single tumor yields important insights about tumor heterogeneity and tumor progression for DCIS to invasive cancer. Foci of invasion were geospatially associated with preinvasive regions of progressively higher mutational load. 1. Rimmer A, Phan H, Mathieson I, Iqbal Z, Twigg SR, Wilkie AO, McVean G, Lunter G, Consortium W: Integrating mapping-, assembly-and haplotype-based approaches for calling variants in clinical sequencing applications. Nature genetics 2014, 46(8):912-918. 2. Lefort V, Desper R, Gascuel O: FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. Molecular biology and evolution 2015:msv150. Citation Format: Ding Y, Marks JR, King LM, Hall AH, Mardis ER, Rodrigo AG, Maley CC, Hwang E-S. Evidence for tumor heterogeneity and clonal evolution during invasive progression of breast cancer [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P1-06-06.

VMCMC: a graphical and statistical analysis tool for Markov chain Monte Carlo traces

BackgroundMCMC-based methods are important for Bayesian inference of phylogeny and related parameters. Although being computationally expensive, MCMC yields estimates of posterior distributions that are useful for estimating parameter values and are easy to use in subsequent analysis. There are, however, sometimes practical difficulties with MCMC, relating to convergence assessment and determining burn-in, especially in large-scale analyses. Currently, multiple software are required to perform, e.g., convergence, mixing and interactive exploration of both continuous and tree parameters.ResultsWe have written a software called VMCMC to simplify post-processing of MCMC traces with, for example, automatic burn-in estimation. VMCMC can also be used both as a GUI-based application, supporting interactive exploration, and as a command-line tool suitable for automated pipelines.ConclusionsVMCMC is a free software available under the New BSD License. Executable jar files, tutorial manual and source code can be downloaded from https://bitbucket.org/rhali/visualmcmc/.

Towards a higher-level Ensifera phylogeny inferred from mitogenome sequences

Although mitogenomes are useful tools for inferring evolutionary history, only a few representative ones can be used for most Ensifera lineages. Thirty-two ensiferan mitogenomes were determined using ABI Sanger sequencing and standard primer walking of 2–3 overlapping Long-PCR fragments, or Illumina® HiSeq2000 for “shotgun” sequenced long-PCR-amplified mitochondrial or total genomic DNA. Six patterns of gene arrangements, including the novel trnR-trnSAGN-trnA-trnN-trnG-nad3 in Lipotactes tripyrga (Lipotactinae), were identified from 59 ensiferan mitogenomes. The results suggest that trnM-trnI-trnQ and trnA-trnR-trnE-trnSAGN-trnN-trnF rearrangements might be a shared derived character in Pseudophyllinae and Gryllidae, respectively. We found base composition biases in our dataset, which potentially complicate the inference of higher-level ensiferan phylogeny. Site-heterogeneous Bayesian inference (BI) and site-homogeneous maximum likelihood (ML) analyses recovered all ensiferan superfamilies as monophyletic. The site-homogeneous BI analysis failed to recover the monophyly of Stenopelmatoidea. As Schizodactyloidea was only represented by Comicus campestris, its monophyly could not be tested. In the Triassic/Jurassic boundary, Ensifera diverged into grylloid and non-grylloid clades. All analyses confirmed Grylloidea and Gryllotalpoidea as sister groups. Site-heterogeneous BI analysis found Schizodactyloidea as the most basal lineage and sister to the clade formed by Grylloidea and Gryllotalpoidea, but the site-homogeneous analyses placed it basally to the non-grylloid clade and recovered a sister relationship between Tettigonioidea and (Hagloidea, Rhaphidophoroidea, Stenopelmatoidea), although this clade had a low support. The site-heterogeneous BI analysis found Tettigonioidea and Hagloidea were sister groups (posterior probability (PP)=0.99), Stenopelmatoidea was sister to (Tettigonioidea, Hagloidea) (PP>0.91), and Rhaphidophoroidea was basal to the non-grylloid clade. At a lower level, all analyses divided Tettigonioidea into Phaneropteridae and Tettigoniidae.

Identification and Characterization of Rhizoctonia Species Associated with Soybean Seedling Disease.

In an effort to identify the Rhizoctonia spp. associated with seedling diseases of soybean, Rhizoctonia isolates were recovered from soybean seedlings with damping off and root and hypocotyl rot symptoms from Arkansas, Illinois, Kansas, Michigan, Minnesota, and the Canadian province of Ontario between 2012 and 2014. Based on cultural morphology, polymerase chain reaction restriction fragment length polymorphism, and phylogenetic analysis of the internal transcribed spacer (ITS) region of the ribosomal RNA genes, 80 isolates were confirmed to be Rhizoctonia solani, 24 were binucleate Rhizoctonia spp., and 10 were R. zeae. Of the 80 R. solani isolates, one belonged to anastomosis group (AG) 2-1, 52 belonged to AG-2-2IIIB, five belonged to AG-3 PT, three belonged to AG-4 HGI, two belonged to AG-4 HGIII, nine belonged to AG-7, and eight belonged to AG-11. Bayesian inference of phylogeny using the ITS region revealed two clades of R. solani AG-7 that possibly correspond to different AG-7 subgroups. Phylogenetic analysis also provided evidence for genetic relatedness between certain binucleate Rhizoctonia and some R. solani isolates. On 'Williams 82' soybean, isolates of AG-2-2IIIB were the most aggressive, followed by isolates of AG-7, AG-4, and AG-11. On 'Jubilee', a sweet corn cultivar, AG-2-2IIIB and AG-4 isolates caused significant stunting and root damage, whereas the damage caused by the AG-11 isolates was mostly restricted to the mesocotyl. Isolates of R. zeae and the binucleate Rhizoctonia spp. were not pathogenic on soybean or corn. Our results indicate that soybean and corn are hosts to the predominant and aggressive AG of R. solani, implying that rotation between these two crops may not be an effective management practice.

A New Species of Pit Mite (Trombidiformes: Harpirhynchidae) from the South American Rattlesnake (Viperidae): Morphological and Molecular Analysis

Background: Mites of the genus Ophioptes, parasitize a wide range of snakes’ species worldwide. Pit mites develop in capsules inside the connective tissue or scales of their hosts and all stages have a genital-anal opening with no connection to the midgut. To this date, there are 15 known species, of which five occur in the Neotropical region. In South America four species have been described from Colubrid snakes. Methods: Mites were collected from the chin shields and infralabial area of the head, and the anterior third portion of the snake. Comparisons of South American species of pit mites are provided for identification purposes. SEM imaging and illustration were made to provide morphological details of the new species. DNA extraction, sequencing, and phylogeny inference were performed of the new mite species and other species of Trombidiformes mites found on reptiles and amphibians. Results: Ophioptes ekans n. sp. is described from the pits made by the mite on the scales and skin of a South American rattlesnake (Crotalus durissus terrificus) in Campo Limpo Paulista, Sao Paulo state, Brazil, captured on January 2014. The Genbank accession numbers of the new species are KU891263, KU891264 and KU891265. DNA sequences were used for molecular phylogenetic inference. Three nymhpal stages were observed for this species. Conclusion: This is the first record of a viper snake from the sub-family Crotaline parasitized by Ophioptes mites. Molecular analyses showed that molecular systematic of Trombidiformes mites is still unclear and more sequences and other genes are needed do better elucidate the relationships within the group. These are the first DNA sequences (18rRNA V4 region) of mites from the Ophioptinae subfamily.

Identification ofArmillariaspecies on declined oak in Britain: implications for oak health

The identity of 51 isolates of Armillaria from 15 Quercus robur trees in poor health, and a single healthy tree, at nine sites in England, was determined using multilocus sequence analysis (MLSA) of three gene regions. Sequences of the ITS-1, IGS-1 and EF-1α gene regions were obtained by polymerase chain reaction (PCR) amplification and sequencing, and phylogenetic trees were generated based on maximum likelihood and Bayesian inference of phylogenies. Four Armillaria species were isolated: Armillaria gallica, A. mellea, A. ostoyae and A. tabescens. Armillaria gallica was most frequently isolated (40/51 isolates), but only from woodland trees. Armillaria mellea was isolated infrequently (3/51), from garden trees; A. tabescens was isolated infrequently (4/51), from trees either in a garden or a parkland location. Armillaria ostoyae (4/51 isolates) was coisolated with A. gallica, raising interesting questions about the synecology of these species, suggesting that more thorough investigations are required to detect all species present on a single host. The distribution of these Armillaria species in Britain and historical information about them on oak are described. It is concluded that further studies are necessary to determine the role of Armillaria in oak declines; A. gallica should be a key focus, but investigations should include polymicrobial interactions with other microorganisms, including other Armillaria species.