Constructing Phylogenies in the Presence Of Intra-Individual Site Polymorphisms (2ISPs) with a Focus on the Nuclear Ribosomal Cistron

Abstract

Nuclear DNA is widely used to estimate phylogenetic and phylogeographic relationships. Nuclear gene variants may be present in an individual's genome, and these result in Intra-Individual Site Polymorphisms (2ISP; pronounced "twisp") in direct-PCR or individual-consensus sequences based on a sample of clones or fragment sequences from next generation sequencing (NGS). 2ISPs can occur fairly often, especially within, but not restricted to, high-copy-number regions such as the widely used internal transcribed spacers of the nuclear ribosomal cistron. Dealing with 2ISPs has been problematic as phylogeny reconstruction optimality criteria generally do not take account of this variation. Here we test whether an approach that treats 2ISPs as additional (termed "informative"), rather than ambiguous, characters offers improved support in three common criteria used for phylogenetic inference: Minimum Evolution (via Neighbour Joining), Maximum Parsimony, and Maximum Likelihood. We demonstrate significant improvements using the 2ISP-informative treatment with simulated, real-world, and case-study data sets. We envisage that this 2ISP-informative approach will greatly aid phylogenetic inference using any nuclear DNA regions that contain polymorphisms within individuals (including consensus sequences generated from NGS), especially at the intrageneric or intraspecific level.

Species relationships and divergence times in beeches: new insights from the inclusion of 53 young and old fossils in a birth–death clock model

“…We used a bipartition network as an alternative means to visualize the output (electronic supplementary material, figure S3). Bipartition networks allow efficient inspection of topological alternatives [11,62,65] and help in assessing whether the fossil placements were as intended and which fossils 'moved' between placements, based on their ages (morphology obviously could not have played a role).…”

Contrasting genetic structure between mitochondrial and nuclear markers in the dengue fever mosquito from Rio de Janeiro: implications for vector control

“…W for adenine or cytosine, Y for cytosine or thymine, etc. ), and the RAXML program considers these ambiguous sites in the phylogenetic reconstruction (Potts et al 2014). We tested 44 substitution models with the Gamma-distributed rate variation among sites (+G) and calculated the AIC weight as a measure of model support in JMODELTEST2.…”

The Colchic region as refuge for relict tree lineages: cryptic speciation in field maples

“…Bootstrap (BS) support was established using the fast bootstrap implementation in RAxML (Stamatakis et al, 2008); number of necessary bootstrap replicates were determined by the extended majority rule bootstrap criterion (Pattengale et al, 2009). The matrix of strict individual-consensus sequences was (a) analysed in the standard setting (GTR-A) using a 4 × 4 general-time reversible substitution model (GTR; Rodriguez et al, 1990) allowing for site-rate variation modelled via a Gamma distribution (+a), followed by (b) treating all nucleotides and polymorphisms as distinct characters (Potts et al, 2014) using the multistate model implemented in RAxML (-m MULTIGAMMA). Two models were applied for (b): a GTR model allowing for unequal substitution rates between all possible states (GTR-I) and Lewis' (2001) Markov-model, which optimises a single substitution probability (MK-I).…”