Assessing the phylogenetic utility of four mitochondrial genes and a nuclear intron in the asian pit viper genus, Trimeresurus: separate, simultaneous, and conditional data combination analyses

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.

Venomics of Trimeresurus (Popeia) nebularis, the Cameron Highlands Pit Viper from Malaysia: Insights into Venom Proteome, Toxicity and Neutralization of Antivenom

Trimeresurus nebularis is a montane pit viper that causes bites and envenomation to various communities in the central highland region of Malaysia, in particular Cameron’s Highlands. To unravel the venom composition of this species, the venom proteins were digested by trypsin and subjected to nano-liquid chromatography-tandem mass spectrometry (LC-MS/MS) for proteomic profiling. Snake venom metalloproteinases (SVMP) dominated the venom proteome by 48.42% of total venom proteins, with a characteristic distribution of P-III: P-II classes in a ratio of 2:1, while P-I class was undetected. Snaclecs constituted the second most venomous protein family (19.43%), followed by snake venom serine proteases (SVSP, 14.27%), phospholipases A₂ (5.40%), disintegrins (5.26%) and minor proteins including cysteine-rich secretory proteins, L-amino acid oxidases, phosphodiesterases, 5′-nucleotidases. The venomic profile correlates with local (painful progressive edema) and systemic (hemorrhage, coagulopathy, thrombocytopenia) manifestation of T. nebularis envenoming. As specific antivenom is unavailable for T. nebularis, the hetero-specific Thai Green Pit viper Monovalent Antivenom (GPVAV) was examined for immunological cross-reactivity. GPVAV exhibited good immunoreactivity to T. nebularis venom and the antivenom effectively cross-neutralized the hemotoxic and lethal effects of T. nebularis (lethality neutralizing potency = 1.6 mg venom per mL antivenom). The findings supported GPVAV use in treating T. nebularis envenoming.

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.

Venomics of Tropidolaemus wagleri, the sexually dimorphic temple pit viper: Unveiling a deeply conserved atypical toxin arsenal

Tropidolaemus wagleri (temple pit viper) is a medically important snake in Southeast Asia. It displays distinct sexual dimorphism and prey specificity, however its venomics and inter-sex venom variation have not been thoroughly investigated. Applying reverse-phase HPLC, we demonstrated that the venom profiles were not significantly affected by sex and geographical locality (Peninsular Malaya, insular Penang, insular Sumatra) of the snakes. Essentially, venoms of both sexes share comparable intravenous median lethal dose (LD₅₀) (0.56-0.63 μg/g) and cause neurotoxic envenomation in mice. LCMS/MS identified six waglerin forms as the predominant lethal principles, comprising 38.2% of total venom proteins. Fourteen other toxin-protein families identified include phospholipase A₂, serine proteinase, snaclec and metalloproteinase. In mice, HPLC fractions containing these proteins showed insignificant contribution to the overall venom lethality. Besides, the unique elution pattern of approximately 34.5% of non-lethal, low molecular mass proteins (3-5 kDa) on HPLC could be potential biomarker for this primitive crotalid species. Together, the study unveiled the venom proteome of T. wagleri that is atypical among many pit vipers as it comprises abundant neurotoxic peptides (waglerins) but little hemotoxic proteinases. The findings also revealed that the venom is relatively well conserved intraspecifically despite the drastic morphological differences between sexes.

Isolation and characterisation of 20 microsatellite loci from the Chinese green tree viper Trimeresurus stejnegeri stejnegeri

The Chinese green tree viper (Trimeresurus stejnegeri stejnegeri), one of the most common snakes in Southeast Asia, can be a good model species for evolutionary and behavioral research. However, there is no high polymorphic co-dominant marker that can be used for individual-based genetic analyses available for this species. Therefore, we developed 20 polymorphic microsatellite loci for T. s. stejnegeri in Taiwan by screening a microsatellite-enriched DNA library. The allele numbers of these loci ranged from 3 to 22, and the observed heterozygosity were 0.042-1.000. The probability of false parent non-exclusion was 2 × 10−5 for those loci that conformed to the Hardy–Weinberg equilibrium. These highly polymorphic microsatellite markers are useful for both individual-based and population-level analysis. Furthermore, this is the first set of microsatellite markers developed for Trimeresurus, and could be applied to closely related species to address various research questions on evolution and animal behavior.

[Summary of phylogeny in family Felidae of Carnivora]

Felidae (cats) is one of the strict carnivorous groups in the order Carnivora, many of which are most familiar and spectacular to us. They are the top predators in the world. Thirty-six of 37 living cat species are considered as either "endangered" or "threatened". The relationships among species of the family Felidae, which evolved recently and rapidly, are difficult to resolve, and have been the subject of debate. Construction of a reliable Felidae phylogeny will be of evolutionarily significance and conservation value. In this paper, we summarized phylogeny of Felidae, including cytological, morphological and molecular evidence, and pointed out the existing phylogenetic problems. This review is expected to guide future researches of Felidae phylogeny, and to lay a theoretic foundation for the protection of this animal group.

New insights into the evolution of intronic sequences of the beta-fibrinogen gene and their application in reconstructing mustelid phylogeny

Mustelidae is the largest and most diverse family in the order Carnivora. The phylogenetic relationships among the subfamilies have especially long been a focus of study. Herein we are among the first to employ two new introns (4 and 7) of the nuclear beta-fibrinogen gene to clarify these enigmatic problems. In addition, two previously available nuclear (IRBP exon 1 and TTR intron 1) and one mt (ND2) data sets were also combined and analyzed simultaneously with the newly obtained sequence data in this study. Detailed characterizations of the two intronic regions not only reveal the remarkable occurrences of short interspersed element (SINE) insertion events, providing a new example supporting the attractive hypothesis that attrition of an earlier retroposition may offer a proper environment for successive retropositions by forming a "dimer-like" structure, but also demonstrate their utility in the resolution of mustelid phylogeny. All of our analyses confirm the assemblage of Mustelinae, Lutrinae, and Melinae with confidence; moreover, two clades within Mustelinae were clearly recognized, i.e., genera Mustela and Martes. Notably, genus Martes of Mustelinae was found to branch off first, followed by Melinae and then a clade containing Lutrinae and genus Mustela of Mustelinae, indicating paraphyly of Mustelinae. In addition, Mephitinae diverges before the other mustelids and the monophyletic Procyonidae in all cases, supporting its elevation to a separate family. Additional independent genetic markers are still in need to resolve the trichotomy among Mephitinae and the other two carnivoran clades, Ailuridae and Procyonidae/non-mephitine Mustelidae.

Phylogenetics of the woodrat genus Neotoma (Rodentia: Muridae): implications for the evolution of phenotypic variation in male external genitalia

Interspecific morphological variation in animal genitalia has long attracted the attention of evolutionary biologists because of the role genital form may play in the generation and/or maintenance of species boundaries. Here we examine the origin and evolution of genital variation in rodents of the muroid genus Neotoma. We test the hypothesis that a relatively rare genital form has evolved only once in Neotoma. We use four mitochondrial and four nuclear markers to evaluate this hypothesis by establishing a phylogenetic framework in which to examine genital evolution. We find intron seven of the beta-fibrinogen gene to be a highly informative nuclear marker for the levels of differentiation that characterize Neotoma with this locus evolving at a rate slower than cytochrome b but faster than 12S. We estimate phylogenetic relationships within Neotoma using both maximum parsimony and maximum likelihood-based Bayesian methods. Our Bayesian and parsimony reconstructions differ in significant ways, but we show that our parsimony analysis may be influenced by long-branch attraction. Furthermore, our estimate of Neotoma phylogeny remains consistent across various data partitioning strategies in the Bayesian analyses. Using ancestral state reconstruction, we find support for the monophyly of taxa that possess the relatively rare genital form. However, we also find support for the independent evolution of the common genital form and discuss possible underlying developmental shifts that may have contributed to our observed patterns of morphological evolution.

Parallel evolution of larval morphology and habitat in the snail-killing fly genus Tetanocera

In this study, we sequenced one nuclear and three mitochondrial DNA loci to construct a robust estimate of phylogeny for all available species of Tetanocera. Character optimizations suggested that aquatic habitat was the ancestral condition for Tetanocera larvae, and that there were at least three parallel transitions to terrestrial habitat, with one reversal. Maximum likelihood analyses of character state transformations showed significant correlations between habitat transitions and changes in four larval morphological characteristics (cuticular pigmentation and three characters associated with the posterior spiracular disc). We provide evidence that phylogenetic niche conservatism has been responsible for the maintenance of aquatic-associated larval morphological character states, and that concerted convergence and/or gene linkage was responsible for parallel morphological changes that were derived in conjunction with habitat transitions. These habitat-morphology associations were consistent with the action of natural selection in facilitating the morphological changes that occurred during parallel aquatic to terrestrial habitat transitions in Tetanocera.

Attributes and congruence of three molecular data sets: Inferring phylogenies among Septoria-related species from woody perennial plants

To improve our understanding of phylogenetic relationships within the anamorphic genus Septoria, three molecular data sets representing 2,417 bp of nuclear and mitochondrial genes were evaluated. Separate gene analyses and combined analyses were performed using first, the maximum parsimony criterion and second, a Bayesian framework. The homogeneity of data partitions was evaluated via a combination of homogeneity partition tests and tree topology incongruence tests before conducting combined analyses. A last incongruence re-evaluation using partitioned Bremer support was performed on the combined tree, which corroborated the previous estimates. After each separate data set attributes were examined, simple explanations were advocated as the causes of the significant incongruences detected. The analysis of multiple gene partitions showed unprecedented phylogenetic resolution within the genus Septoria that supported the results from previously published single gene phylogenies. Specifically, we have delimited distinct but closely related species representing monophyletic groups that frequently correlated with their respective host families. Conversely, the occurrence of well-supported groups including closely related but distinct molecular taxa sampled on unrelated host-plants allowed us to reject, in these particular cases, the co-evolutionary concept expected between a parasite and its host and to discuss alternative evolutionary models recently proposed for these pathogens.

Phylogeny of the caniform carnivora: evidence from multiple genes

The monophyletic group Caniformia in the order Carnivora currently comprises seven families whose relationships remain contentious. The phylogenetic positions of the two panda species within the Caniformia have also been evolutionary puzzles over the past decades, especially for Ailurus fulgens (the red panda). Here, new nuclear sequences from two introns of the beta-fibrinogen gene (beta-fibrinogen introns 4 and 7) and a complete mitochondrial (mt) gene (ND2) from 17 caniform representatives were explored for their utilities in resolving higher-level relationships in the Caniformia. In addition, two previously available nuclear (IRBP exon 1 and TTR intron 1) data sets were also combined and analyzed simultaneously with the newly obtained sequence data in this study. Combined analyses of four nuclear and one mt genes (4417 bp) recover a branching order in which almost all nodes were strongly supported. The present analyses provide evidence in favor of Ailurus fulgens as the closest taxon to the procyonid-mustelid (i.e., Musteloidea sensu stricto) clade, followed by pinnipeds (i.e., Otariidae and Phocidae), Ursidae (including Ailuropoda melanoleuca), and Canidae, the most basal lineage in the Caniformia. The potential utilities of different genes in the context of caniform phylogeny were also evaluated, with special attention to the previously unexplored beta-fibrinogen intron 4 and 7 genes.

Bayesian mixed models and the phylogeny of pitvipers (Viperidae: Serpentes)

The subfamily Crotalinae (pitvipers) contains over 190 species of venomous snakes distributed in both the Old and New World. We incorporated an extensive sampling of taxa (including 28 of 29 genera), and sequences of four mitochondrial gene fragments (2.3kb) per individual, to estimate the phylogeny of pitvipers based on maximum parsimony and Bayesian phylogenetic methods. Our Bayesian analyses incorporated complex mixed models of nucleotide evolution that allocated independent models to various partitions of the dataset within combined analyses. We compared results of unpartitioned versus partitioned Bayesian analyses to investigate how much unpartitioned (versus partitioned) models were forced to compromise estimates of model parameters, and whether complex models substantially alter phylogenetic conclusions to the extent that they appear to extract more phylogenetic signal than simple models. Our results indicate that complex models do extract more phylogenetic signal from the data. We also address how differences in phylogenetic results (e.g., bipartition posterior probabilities) obtained from simple versus complex models may be interpreted in terms of relative credibility. Our estimates of pitviper phylogeny suggest that nearly all recently proposed generic reallocations appear valid, although certain Old and New World genera (Ovophis, Trimeresurus, and Bothrops) remain poly- or paraphyletic and require further taxonomic revision. While a majority of nodes were resolved, we could not confidently estimate the basal relationships among New World genera and which lineage of Old World species is most closely related to this New World group.

Optimal Intron Analyses in the Trimeresurus Radiation of Asian Pitvipers

Nuclear introns are commonly used as phylogenetic markers, but a number of issues related to alignment strategies, indel treatments, and the incorporation of length-variant heterozygotes (LVHs) are not routinely addressed when generating phylogenetic hypotheses. Topological congruence in relation to an extensive mitochondrial DNA multigene phylogeny (derived from 2,423 bp of 12S, 16S, ND4, and CYTB genes) of the Asian pitviper Trimeresurus radiation was used to compare combinations of "by eye" and edited and unedited ClustalX 1.8 alignments of two nuclear introns. Indels were treated as missing data, fifth character states, and assigned simple and multistate codes. Upon recovery of the optimal alignment and indel treatment strategy, a total evidence approach was used to investigate the phylogenetic utility of the indels and test new generic arrangements within Trimeresurus. Approximately one third of the intron data partitions exhibited LVHs, suggesting that they are common in introns. Furthermore, a simple concatenation approach can facilitate the incorporation of LVHs into phylogenetic analyses to make use of all available data and investigate mechanisms of molecular evolution. Analyses of ClustalX 1.8-assisted alignments were generally more congruent than the "by eye" alignment and the analysis of a simple coded, edited ClustalX 1.8 (gap opening cost 5, gap extension cost 1) alignment revealed the most congruent tree. The total evidence approach supported the new arrangements within Trimeresurus, suggesting that the phylogeny should be considered as a working benchmark in Asian pitviper systematics. Finally, a critical appraisal of the diverse array of indels (56 to 57 per intron, ranging from 1 to 151 bp in length) suggested that they are a combination of Hennigian and homoplasious events unrelated to indel size or location within the intron. [Alignment; indels; intron analysis; length-variant heterozygotes; Trimeresurus.].

Phylogenetic studies of pantherine cats (Felidae) based on multiple genes, with novel application of nuclear beta-fibrinogen intron 7 to carnivores

The pantherine lineage of the cat family Felidae (order: Carnivora) includes five big cats of genus Panthera and a great many midsized cats known worldwide. Presumably because of their recent and rapid radiation, the evolutionary relationship among pantherines remains ambiguous. We provide an independent assessment of the evolutionary history of pantherine lineage using two complete mitochondrial (mt) genes (ND2 and ND4) and the nuclear beta-fibrinogen intron 7 gene, whose utility in carnivoran phylogeny was first explored. The available four mt (ND5, cytb, 12S, and 16SrRNA) and two nuclear (IRBP and TTR) sequence loci were also combined to reconstruct phylogeny of 14 closely related cat species. Our analyses of combined mt data (six genes; approximately 3750 bp) and combined mt and nuclear data (nine genes; approximately 6500 bp) obtained identical tree topologies, which were well-resolved and strongly supported for almost all nodes. Monophyly of Panthera genus in pantherine lineage was confirmed and interspecific affinities within this genus revealed a novel branching pattern, with P. tigris diverging first in Panthera genus, followed by P. onca, P. leo, and last two sister species P. pardus and P. uncia. In addition, close association of Neofelis nebulosa to Panthera, the phylogenetic redefinition of Otocolobus manul within the domestic cat group, and the relatedness of Acinonyx jubatus and Puma concolor were all important findings in the resulting phylogenies. The potential utilities of nine different genes for phylogenetic resolution of closely related pantherine species were also evaluated, with special interest in that of the novel nuclear beta-fibrinogen intron 7.

Unusual venom phospholipases A2 of two primitive tree vipers Trimeresurus puniceus and Trimeresurus borneensis

To explore the venom diversity of Asian pit vipers, we investigated the structure and function of venom phospholipase A2 (PLA2) derived from two primitive tree vipers Trimeresurus puniceus and Trimeresurus borneensis. We purified six novel PLA2s from T. puniceus venom and another three from T. borneensis venom. All cDNAs encoding these PLA2s except one were cloned, and the molecular masses and N-terminal sequences of the purified enzymes closely matched those predicted from the cDNA. Three contain K49 and lack a disulfide bond at C61-C91, in contrast with the D49-containing PLA2s in both venom species. They are less thermally stable than other K49-PLA2s which contain seven disulfide bonds, as indicated by a decrease of 8.8 degrees C in the melting temperature measured by CD spectroscopy. The M110D mutation in one of the K49-PLA2s apparently reduced its edematous potency. A phylogenetic tree based on the amino-acid sequences of 17 K49-PLA2s from Asian pit viper venoms illustrates close relationships among the Trimeresurus species and intergeneric segregations. Basic D49-PLA2s with a unique Gly6 substitution were also purified from both venoms. They showed edema-inducing and anticoagulating activities. It is notable that acidic PLA2s from both venoms inhibited blood coagulation rather than platelet aggregation, and this inhibition was only partially dependent on enzyme activity. These results contribute to our understanding of the evolution of Trimeresurus pit vipers and the structure-function relationships between various subtypes of crotalid venom PLA2.

Comparative proteomics and subtyping of venom phospholipases A2 and disintegrins of Protobothrops pit vipers

To explore the venom diversity and systematics of pit vipers under the genus Protobothrops, the venom phospholipases A2 (PLA2s) of P. mangshanensis, P. elegans and P. tokarensis were purified and characterized for the first time. The results were compared with the corresponding venom data of other co-generic species including P. mucrosquamatus, P. flavoviridis and P. jerdonii. Based on sequence features at the N-terminal regions, we identified five PLA2 subtypes, i.e., the Asp49-PLA2s with N6, E6 or R6 substitution and the Lys49-PLA2. However, not all subtypes were expressed in each of the species. Venom N6-PLA2s from P. mangshanensis and P. tokarensis venom were weakly neurotoxic toward chick biventer cervicis tissue preparations. The venoms of P. tokarensis and P. flavoviridis contained identical PLA2 isoforms. In most Protobothrop disintegrins, sequences flanking the RGD-motif are conserved. Phylogenetic analyses based on amino acid sequences of both families of the acidic PLA2s and the disintegrins clarify that these species could belong to a monophyletic group.

Molecular phylogeny of the hexagrammid fishes using a multi-locus approach

Ideally, organisms are grouped into monophyletic assemblages reflecting their evolutionary histories. Single (molecular) markers can reflect the evolutionary history of the marker, rather than the species in question, therefore, phylogenetic relationships should be inferred from adequate sampling of characters. Because the use of multiple loci greatly improves the resolving power of the molecular assay, we constructed a molecular phylogeny of the family Hexagrammidae based on six loci, including two mitochondrial and four nuclear loci. The resulting molecular phylogeny, from the combined data, was significantly different from the morphological topology suggested by Shinohara [Memoirs of the Faculty of Fisheries, Hokkaido University 41 (1994) 1]. Our data support a monophyletic assemblage for the genera Hexagrammos and Pleurogrammus. However, other taxa traditionally included in the family Hexagrammidae did not form a monophyletic assemblage. The monotypic genus Ophiodon was more closely associated with cottids than with other hexagrammids. Our data concur with the morphological topology in that the genera Zaniolepis and Oxylebius formed a monophyletic clade, which was distinct and basal to the remaining hexagrammids, seven cottids and one agonid.

A phylogeny of four mitochondrial gene regions suggests a revised taxonomy for Asian pitvipers (Trimeresurus and Ovophis)

We present a phylogeny of the Asian pitvipers, based on 2403 bp of four mitochondrial gene regions. All but six known species of Trimeresurus sensu stricto (s.s.) as currently defined, as well as multiple populations of widespread species, which may yet be described as full species, and representatives of all other Asian pitviper genera, are included. Both the greater sampling and larger dataset provide improved resolution over previous studies and support the existence of distinct species groups within Trimeresurus s.s. Although all but two species currently referred to this genus form a monophyletic group, morphological and molecular analyses identify four subgroups that warrant recognition at the generic level. We propose a new generic arrangement to reflect these findings. We also highlight the non-monophyly of Ovophis, and propose a new genus to accommodate a species formerly assigned to Ovophis.