Optimal Intron Analyses in the Trimeresurus Radiation of Asian Pitvipers

The Artefactual Branch Effect and Phylogenetic Conflict: Species Delimitation with Gene Flow in Mangrove Pit Vipers (Trimeresurus purpureomaculatus-erythrurus Complex)

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to naturally inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate. Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent nonmonophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict.

A new species of pit-viper from the Ayeyarwady and Yangon regions in Myanmar (Viperidae, Trimeresurus)

In a genomic study by Chan and colleagues, pit-vipers of the Trimeresuruserythrurus-purpureomaculatus complex from the Ayeyarwady and Yangon regions in Myanmar were demonstrated to be a distinct species based on robust population genetic and species delimitation analyses. Here, we provide morphological characterizations and a formal description of those populations as a new species. The new species, Trimeresurusayeyarwadyensissp. nov., is most closely related to T.erythrurus and T.purpureomaculatus and shares morphological characteristics with both of those species. Some specimens of T.ayeyarwadyensissp. nov. have green dorsal coloration and no distinct dorsal blotches (a trait shared with T.erythrurus but not T.purpureomaculatus), while others have dark dorsal blotches (a trait shared with T.purpureomaculatus but not T.erythrurus). The distinct evolutionary trajectory of the new species, coupled with the lack of obvious morphological differentiation, represents a classic example of the cryptic nature of species commonly found in the Trimeresurus group of Asian pit-vipers and underscores the need for data-rich analyses to verify species' boundaries more broadly within this genus.

De Novo Assembly of Venom Gland Transcriptome of Tropidolaemus wagleri (Temple Pit Viper, Malaysia) and Insights into the Origin of Its Major Toxin, Waglerin

The venom proteome of Temple Pit Viper (Tropidolaemus wagleri) is unique among pit vipers, characterized by a high abundance of a neurotoxic peptide, waglerin. To further explore the genetic diversity of its toxins, the present study de novo assembled the venom gland transcriptome of T. wagleri from west Malaysia. Among the 15 toxin gene families discovered, gene annotation and expression analysis reveal the dominating trend of bradykinin-potentiating peptide/angiotensin-converting enzyme inhibitor-C-type natriuretic peptide (BPP/ACEI-CNP, 76.19% of all-toxin transcription) in the transcriptome, followed by P-III snake venom metalloproteases (13.91%) and other toxins. The transcript TwBNP01 of BPP/ACEI-CNP represents a large precursor gene (209 amino acid residues) containing the coding region for waglerin (24 residues). TwBNP01 shows substantial sequence variations from the corresponding genes of its sister species, Tropidolaemus subannulatus of northern Philippines, and other viperid species which diversely code for proline-rich small peptides such as bradykinin-potentiating peptides (BPPs). The waglerin/waglerin-like peptides, BPPs and azemiopsin are proline-rich, evolving de novo from multiple highly diverged propeptide regions within the orthologous BPP/ACEI-CNP genes. Neofunctionalization of the peptides results in phylogenetic constraints consistent with a phenotypic dichotomy, where Tropidolaemus spp. and Azemiops feae convergently evolve a neurotoxic trait while vasoactive BPPs evolve only in other species.

T. Lalremsanga EO, Bhosale H, Gowande G, Patel H, Idiatullina SS, Poyarkov NA. Systematics of Trimeresurus popeiorum Smith, 1937 with a revised molecular phylogeny of Asian pitvipers of the genus Trimeresurus Lacépède, 1804 sensu lato

The pit viper snake genus Trimeresurus Lacépède, 1804 sensu lato, is a diverse group of nocturnal serpents comprising over 61 species. The genus is morphologically heterogeneous and has been divided into several subgenera. We present an updated phylogeny of Asian pitvipers and propose a revised classification. Additionally, we revise the taxonomy of T. popeiorum Smith, 1937 and propose taxonomic changes with support from molecular and morphological data. We restrict T. popeiorumsensu stricto to northeastern India, Bangladesh, southern China, and northern Myanmar; populations beyond these areas require further assessment. We also synonymize T. yingjiangensisChen et al., 2019 with T. popeiorum based on overlapping morphological characters, molecular data, and distribution. The findings shed new light on the taxonomy of T. popeiorum, warranting the need for assessing the population of T. popeiorum from southeast Asia.

A large-scale systematic framework of Chinese snakes based on a unified multilocus marker system

Snakes are one of the most diverse groups of terrestrial vertebrates, with approximately 3500 extant species. A robust phylogeny and taxonomy of snakes is crucial for us to know, study and protect them. For a large group such as snakes, broad-scale phylogenetic reconstructions largely rely on data integration. Increasing the compatibility of the data from different researches is thus important, which can be facilitated by standardization of the loci used in systematic analyses. In this study, we proposed a unified multilocus marker system for snake systematics by conflating 5 mitochondrial markers, 19 vertebrate-universal nuclear protein coding (NPC) markers and 72 snake-specific noncoding intron markers. This marker system is an addition to the large squamate conserved locus set (SqCL) for studies preferring a medium-scale data set. We applied this marker system to over 440 snake samples and constructed the currently most comprehensive systematic framework of the snakes in China. Robust snake phylogenetic relationships were recovered at both deep and shallow evolutionary depths, demonstrating the usefulness of this multilocus marker system. Discordance was revealed by a parallel comparison between the snake tree based on the multilocus marker system and that based on only the mitochondrial loci, highlighting the necessity of using multiple types of markers to better understand the snake evolutionary histories. The divergence times of different snake groups were estimated with the nuclear data set. Our comprehensive snake tree not only confirms many important nodes inferred in previous studies but also contributes new insights into many snake phylogenetic relationships. Suggestions are made for the current Chinese snake taxonomy.

Multilocus approach reveals discordant molecular markers and corridors for gene flow between North African populations of Fasciola hepatica

Fasciolosis is a foodborne trematodosis characterised by a worldwide distribution. Various approaches have been developed for the study of the causative agents of this parasitic infection: Fasciola hepatica, Fasciola gigantica and the aspermic intermediated forms (hybrid and introgressed). In the present study, novel and common molecular markers (pepck and pold, ITS, CO1, ND1 and CO1-trnT-rrnL) were used to characterise Fasciola flukes from the Tunisian-Algerian border, to estimate the gene flow between these populations and to evaluate the reliability of different molecular markers. All nuclear and mitochondrial markers, apart from pepck, supported the monophyly of the studied flukes identified as F. hepatica. Multiplex PCR for pepck revealed three different genotypes corresponding to F. hepatica (pepck-Fh), F. gigantica (pepck-Fg) and the aspermic Fasciola flukes (pepck-Fh/Fg). Sequence analysis of pepck revealed high polymorphism, length variation, within this intronic marker. The observed inconsistencies were due to the position of the forward primer within the intronic region. Pepck sequences showed different level of heterozygosity and homozygosity with length polymorphisms in the introns. Pepck multiplex PCR patterns could not differentiate between Fasciola species. All studies based on only pepck multiplex PCR with mitochondrial markers should be revised. Nuclear and mitochondrial markers revealed an important gene flow between Tunisian and Algerian populations of F. hepatica. The combination of nuclear and mitochondrial sequence analysis is still the best method to distinguish these taxa. Effective measures are needed in order to better control cross-country illegal trade of vector.

A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein-coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single-copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next-generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p-distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.

Choosing and Using Introns in Molecular Phylogenetics

Introns are now commonly used in molecular phylogenetics in an attempt to recover gene trees that are concordant with species trees, but there are a range of genomic, logistical and analytical considerations that are infrequently discussed in empirical studies that utilize intron data. This review outlines expedient approaches for locus selection, overcoming paralogy problems, recombination detection methods and the identification and incorporation of LVHs in molecular systematics. A range of parsimony and Bayesian analytical approaches are also described in order to highlight the methods that can currently be employed to align sequences and treat indels in subsequent analyses. By covering the main points associated with the generation and analysis of intron data, this review aims to provide a comprehensive introduction to using introns (or any non-coding nuclear data partition) in contemporary phylogenetics.

Complete mitochondrial genome of Camponotus atrox (Hymenoptera: Formicidae): a new tRNA arrangement in Hymenoptera

We sequenced the complete mitochondrial (mt) genome of Camponotus atrox (Hymenoptera: Formicidae), which is only distributed in Korea. The genome was 16 540 bp in size and contained typical sets of genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs). The C. atrox A+T-rich region, at 1402 bp, was the longest of all sequenced ant genomes and was composed of an identical tandem repeat consisting of six 100-bp copies and one 96-bp copy. A total of 315 bp of intergenic spacer sequence was spread over 23 regions. An alignment of the spacer sequences in ants was largely feasible among congeneric species, and there was substantial sequence divergence, indicating their potential use as molecular markers for congeneric species. The A/T contents at the first and second codon positions of protein-coding genes (PCGs) were similar for ant species, including C. atrox (73.9% vs. 72.3%, on average). With increased taxon sampling among hymenopteran superfamilies, differences in the divergence rates (i.e., the non-synonymous substitution rates) between the suborders Symphyta and Apocrita were detected, consistent with previous results. The C. atrox mt genome had a unique gene arrangement, trnI-trnM-trnQ, at the A+T-rich region and ND2 junction (underline indicates inverted gene). This may have originated from a tandem duplication of trnM-trnI, resulting in trnM-trnI-trnM-trnI-trnQ, and the subsequent loss of the first trnM and second trnI, resulting in trnI-trnM-trnQ.

Incorporating indels as phylogenetic characters: impact for interfamilial relationships within Arctoidea (Mammalia: Carnivora)

Insertion and deletion events (indels) provide a suite of markers with enormous potential for molecular phylogenetics. Using many more indel characters than those in previous studies, we here for the first time address the impact of indel inclusion on the phylogenetic inferences of Arctoidea (Mammalia: Carnivora). Based on 6843 indel characters from 22 nuclear intron loci of 16 species of Arctoidea, our analyses demonstrate that when the indels were not taken into consideration, the monophyly of Ursidae and Pinnipedia tree and the monophyly of Pinnipedia and Musteloidea tree were both recovered, whereas inclusion of indels by using three different indel coding schemes give identical phylogenetic tree topologies supporting the monophyly of Ursidae and Pinnipedia. Our work brings new perspectives on the previously controversial placements among Arctoidea families, and provides another example demonstrating the importance of identifying and incorporating indels in the phylogenetic analyses of introns. In addition, comparison of indel incorporation methods revealed that the three indel coding methods are all advantageous over treating indels as missing data, given that incorporating indels produces consistent results across methods. This is the first report of the impact of different indel coding schemes on phylogenetic reconstruction at the family level in Carnivora, which indicates that indels should be taken into account in the future phylogenetic analyses.

Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia

Were molecular data available for extinct taxa, questions regarding the origins of many groups could be settled in short order. As this is not the case, various strategies have been proposed to combine paleontological and neontological data sets. The use of fossil dates as node age calibrations for divergence time estimation from molecular phylogenies is commonplace. In addition, simulations suggest that the addition of morphological data from extinct taxa may improve phylogenetic estimation when combined with molecular data for extant species, and some studies have merged morphological and molecular data to estimate combined evidence phylogenies containing both extinct and extant taxa. However, few, if any, studies have attempted to estimate divergence times using phylogenies containing both fossil and living taxa sampled for both molecular and morphological data. Here, I infer both the phylogeny and the time of origin for Lissamphibia and a number of stem tetrapods using Bayesian methods based on a data set containing morphological data for extinct taxa, molecular data for extant taxa, and molecular and morphological data for a subset of extant taxa. The results suggest that Lissamphibia is monophyletic, nested within Lepospondyli, and originated in the late Carboniferous at the earliest. This research illustrates potential pitfalls for the use of fossils as post hoc age constraints on internal nodes and highlights the importance of explicit phylogenetic analysis of extinct taxa. These results suggest that the application of fossils as minima or maxima on molecular phylogenies should be supplemented or supplanted by combined evidence analyses whenever possible.

On the phylogeny of Mustelidae subfamilies: analysis of seventeen nuclear non-coding loci and mitochondrial complete genomes

BACKGROUND

Mustelidae, as the largest and most-diverse family of order Carnivora, comprises eight subfamilies. Phylogenetic relationships among these Mustelidae subfamilies remain argumentative subjects in recent years. One of the main reasons is that the mustelids represent a typical example of rapid evolutionary radiation and recent speciation event. Prior investigation has been concentrated on the application of different mitochondrial (mt) sequence and nuclear protein-coding data, herein we employ 17 nuclear non-coding loci (>15 kb), in conjunction with mt complete genome data (>16 kb), to clarify these enigmatic problems.

RESULTS

The combined nuclear intron and mt genome analyses both robustly support that Taxidiinae diverged first, followed by Melinae. Lutrinae and Mustelinae are grouped together in all analyses with strong supports. The position of Helictidinae, however, is enigmatic because the mt genome analysis places it to the clade uniting Lutrinae and Mustelinae, whereas the nuclear intron analysis favors a novel view supporting a closer relationship of Helictidinae to Martinae. This finding emphasizes a need to add more data and include more taxa to resolve this problem. In addition, the molecular dating provides insights into the time scale of the origin and diversification of the Mustelidae subfamilies. Finally, the phylogenetic performances and limits of nuclear introns and mt genes are discussed in the context of Mustelidae phylogeny.

CONCLUSION

Our study not only brings new perspectives on the previously obscured phylogenetic relationships among Mustelidae subfamilies, but also provides another example demonstrating the effectiveness of nuclear non-coding loci for reconstructing evolutionary histories in a group that has undergone rapid bursts of speciation.

Second-generation environmental sequencing unmasks marine metazoan biodiversity

Biodiversity is of crucial importance for ecosystem functioning, sustainability and resilience, but the magnitude and organization of marine diversity at a range of spatial and taxonomic scales are undefined. In this paper, we use second-generation sequencing to unmask putatively diverse marine metazoan biodiversity in a Scottish temperate benthic ecosystem. We show that remarkable differences in diversity occurred at microgeographical scales and refute currently accepted ecological and taxonomic paradigms of meiofaunal identity, rank abundance and concomitant understanding of trophic dynamics. Richness estimates from the current benchmarked Operational Clustering of Taxonomic Units from Parallel UltraSequencing analyses are broadly aligned with those derived from morphological assessments. However, the slope of taxon rarefaction curves for many phyla remains incomplete, suggesting that the true alpha diversity is likely to exceed current perceptions. The approaches provide a rapid, objective and cost-effective taxonomic framework for exploring links between ecosystem structure and function of all hitherto intractable, but ecologically important, communities.

Recent developments in sequencing technologies have provided the opportunity to investigate the biodiversity of ecosystems. Such a metagenomic approach, combined with taxon clustering, is used here to demonstrate that the species richness of a marine community in Scotland is much greater than anticipated.

Genomic scans support repetitive continental colonization events during the rapid radiation of voles (Rodentia: Microtus): the utility of AFLPs versus mitochondrial and nuclear sequence markers

Single locus studies might not resolve phylogenetic relationships and the evolutionary history of taxa. The analysis of multiple markers promises higher resolution, and congruence among loci may indicate that the phylogenies represent the underlying species history. Here, we examine the utility of a genome-wide approach based on amplified fragment length polymorphisms (AFLP) and several DNA sequence markers in resolving phylogenetic signals in the rapidly radiating rodent genus Microtus which produced about 70 vole species within the last 1.2-2 myr. The current Holarctic distribution of Microtus is assumed to have resulted from three independent colonization events out of Asia to North America, Europe, and northern Asia without subsequent colonization, which would have led to deep splits between species from different continents. We investigated this hypothesis of three single colonization events by reconstructing the phylogenetic relationships among species from all three continents based on data from the first exon of the nuclear arginine vasopressin receptor 1a gene (EXON1), an adjacent noncoding region and the mitochondrial cytochrome b gene. The phylogenetic patterns obtained from these sequence markers are contrasted to genome-wide data on more than 1800 amplified fragment length polymorphisms (AFLP) analyzed for the same samples. Our results show that the single sequence markers partially resolve the phylogenetic relationships within Microtus, but with some incongruence mostly between EXON1 and the other loci. However, deeper nodes of the radiation are only weakly supported and neither the combination of the markers nor additional nuclear sequences improved the resolution significantly. AFLPs provided much stronger support for major continent-specific clades, and show also that reciprocal monophyly of American and European voles is incomplete. Our results demonstrate that Microtus voles colonized the American and European continents each repeatedly in several independent events on similar colonization routes during their radiation. More generally, this study supports the suitability of AFLPs as an alternative to sequence markers to resolve the evolutionary history of rapidly radiating taxa.

[Intra-individual allele heterozygotes in phylogenetic analysis]

The study of evolutionary relationships among organisms is vital in evolutionary biology. To reconstruct a reliable species phylogeny, one of the most important issues is to choose proper molecular markers and take full advantage of phylogenetic information contained in these markers. Intra-individual allele heterozygotes (IIAHs) have been commonly detected in intron phylogenetic studies. How to incorporate IIAHs into phylogenetic framework has been a focus in current studies. In this review, the conception, isolation, and analytic methods of IIAHs in phylogeny were summarized.

Toward understanding the distribution of Laurasian frogs: a test of Savage's biogeographical hypothesis using the genus Bombina

Several anuran groups of Laurasian origin are each co-distributed in four isolated regions of the Northern Hemisphere: central/southern Europe and adjacent areas, Korean Peninsula and adjacent areas, Indo-Malaya, and southern North America. Similar distribution patterns have been observed in diverse animal and plant groups. Savage [Savage, J.M., 1973. The geographic distribution of frogs: patterns and predictions. In: Vial, J.L. (Ed.), Evolutionary Biology of the Anurans. University of Missouri Press, Columbia, pp. 351-445] hypothesized that the Miocene global cooling and increasing aridities in interiors of Eurasia and North America caused a southward displacement and range contraction of Laurasian frogs (and other groups). We use the frog genus Bombina to test Savage's biogeographical hypothesis. A phylogeny of Bombina is reconstructed based on three mitochondrial and two nuclear gene fragments. The genus is divided into three major clades: an Indo-Malaya clade includes B. fortinuptialis, B. lichuanensis, B. maxima, and B. microdeladigitora; a European clade includes B. bombina, B. pachypus, and B. variegata; and a Korean clade contains B. orientalis. The European and Korean clades form sister-group relationship. Molecular dating of the phylogenetic tree using the penalized likelihood and Bayesian analyses suggests that the divergence between the Indo-Malaya clade and other Bombina species occurred 5.9-28.6 million years ago. The split time between the European clade and the Korean clade is estimated at 5.1-20.9 million years ago. The divergence times of these clades are not significantly later than the timing of Miocene cooling and drying, and therefore can not reject Savage's hypothesis. Some other aspects of biogeography of Bombina also are discussed. The Korean Peninsula and the Shandong Peninsula might have supplied distinct southern refugia for B. orientalis during the Pleistocene glacial maxima. In the Indo-Malaya clade, the uplift of the Tibetan Plateau might have promoted the split between B. maxima and the other species.

Spliceosomal introns as tools for genomic and evolutionary analysis

Over the past 5 years, the availability of dozens of whole genomic sequences from a wide variety of eukaryotic lineages has revealed a very large amount of information about the dynamics of intron loss and gain through eukaryotic history, as well as the evolution of intron sequences. Implicit in these advances is a great deal of information about the structure and evolution of surrounding sequences. Here, we review the wealth of ways in which structures of spliceosomal introns as well as their conservation and change through evolution may be harnessed for evolutionary and genomic analysis. First, we discuss uses of intron length distributions and positions in sequence assembly and annotation, and for improving alignment of homologous regions. Second, we review uses of introns in evolutionary studies, including the utility of introns as indicators of rates of sequence evolution, for inferences about molecular evolution, as signatures of orthology and paralogy, and for estimating rates of nucleotide substitution. We conclude with a discussion of phylogenetic methods utilizing intron sequences and positions.

Deep genetic divergences among Indo-Pacific populations of the coral reef sponge Leucetta chagosensis (Leucettidae): founder effects, vicariance, or both?

Background

An increasing number of studies demonstrate that genetic differentiation and speciation in the sea occur over much smaller spatial scales than previously appreciated given the wide distribution range of many morphologically defined coral reef invertebrate species and the presumed dispersal-enhancing qualities of ocean currents. However, knowledge about the processes that lead to population divergence and speciation is often lacking despite being essential for the understanding, conservation, and management of marine biodiversity. Sponges, a highly diverse, ecologically and economically important reef-invertebrate taxon, exhibit spatial trends in the Indo-West Pacific that are not universally reflected in other marine phyla. So far, however, processes generating those unexpected patterns are not understood.

Results

We unraveled the phylogeographic structure of the widespread Indo-Pacific coral reef sponge Leucetta chagosensis across its known geographic range using two nuclear markers: the rDNA internal transcribed spacers (ITS 1&2) and a fragment of the 28S gene, as well as the second intron of the ATP synthetase beta subunit-gene (ATPSb-iII). This enabled the detection of several deeply divergent clades congruent over both loci, one containing specimens from the Indian Ocean (Red Sea and Maldives), another one from the Philippines, and two other large and substructured NW Pacific and SW Pacific clades with an area of overlap in the Great Barrier Reef/Coral Sea. Reciprocally monophyletic populations were observed from the Philippines, Red Sea, Maldives, Japan, Samoa, and Polynesia, demonstrating long-standing isolation. Populations along the South Equatorial Current in the south-western Pacific showed isolation-by-distance effects. Overall, the results pointed towards stepping-stone dispersal with some putative long-distance exchange, consistent with expectations from low dispersal capabilities.

Conclusion

We argue that both founder and vicariance events during the late Pliocene and Pleistocene were responsible to varying degrees for generating the deep phylogeographic structure. This structure was perpetuated largely as a result of the life history of L. chagosensis, resulting in high levels of regional isolation. Reciprocally monophyletic populations constitute putative sibling (cryptic) species, while population para- and polyphyly may indicate incipient speciation processes. The genetic diversity and biodiversity of tropical Indo-Pacific sponges appears to be substantially underestimated since the high level of genetic divergence is not necessarily manifested at the morphological level.

Molecular phylogenetics of the lizard genus Microlophus (squamata:tropiduridae): aligning and retrieving indel signal from nuclear introns

We use a multigene data set (the mitochondrial locus and nine nuclear gene regions) to test phylogenetic relationships in the South American "lava lizards" (genus Microlophus) and describe a strategy for aligning noncoding sequences that accounts for differences in tempo and class of mutational events. We focus on seven nuclear introns that vary in size and frequency of multibase length mutations (i.e., indels) and present a manual alignment strategy that incorporates insertions and deletions (indels) for each intron. Our method is based on mechanistic explanations of intron evolution that does not require a guide tree. We also use a progressive alignment algorithm (Probabilistic Alignment Kit; PRANK) and distinguishes insertions from deletions and avoids the "gapcost" conundrum. We describe an approach to selecting a guide tree purged of ambiguously aligned regions and use this to refine PRANK performance. We show that although manual alignment is successful in finding repeat motifs and the most obvious indels, some regions can only be subjectively aligned, and there are limits to the size and complexity of a data matrix for which this approach can be taken. PRANK alignments identified more parsimony-informative indels while simultaneously increasing nucleotide identity in conserved sequence blocks flanking the indel regions. When comparing manual and PRANK with two widely used methods (CLUSTAL, MUSCLE) for the alignment of the most length-variable intron, only PRANK recovered a tree congruent at deeper nodes with the combined data tree inferred from all nuclear gene regions. We take this concordance as an objective function of alignment quality and present a strongly supported phylogenetic hypothesis for Microlophus relationships. From this hypothesis we show that (1) a coded indel data partition derived from the PRANK alignment contributed significantly to nodal support and (2) the indel data set permitted detection of significant conflict between mitochondrial and nuclear data partitions, which we hypothesize arose from secondary contact of distantly related taxa, followed by hybridization and mtDNA introgression.

Secondary structure of ITS2 rRNA provides taxonomic characters for systematic studies--a case in Lycoperdaceae (Basidiomycota)

The secondary structure of the ITS2 rDNA transcript (pre-rRNA) could provide information for identifying homologous nucleotide characters useful for cladistic inference of relationships. Such structure data could become taxonomic characters. This work compares the effect of several modern nucleotide alignment strategies, including those making use of structure data, on phylogenetic inference. From both the phylogenetic analyses and comparative secondary structure, implications for taxonomy and evolution of puffball fungi are discussed. Lycoperdaceae remain insufficiently resolved with present taxon and data sampling. Neither alignment allows statistically robust phylogenetic hypotheses under any current optimality criterion. The secondary structure data at this time are best used as accessory taxonomic characters as their phylogenetic resolving power and confidence in validity is limited compared with underlying nucleotide characters. We introduce a preliminary nomenclature convention to describe secondary structure for defining consensus features. These consensus structures are illustrated for the clades /Calvatia, /Handkea-Echinatum, /Vascellum, /Morganella, and /Plumbea-Paludosa (Bovista).