Bayesian Phylogenetic Inference Research Articles

Parallel power posterior analyses for fast computation of marginal likelihoods in phylogenetics.

In Bayesian phylogenetic inference, marginal likelihoods can be estimated using several different methods, including the path-sampling or stepping-stone-sampling algorithms. Both algorithms are computationally demanding because they require a series of power posterior Markov chain Monte Carlo (MCMC) simulations. Here we introduce a general parallelization strategy that distributes the power posterior MCMC simulations and the likelihood computations over available CPUs. Our parallelization strategy can easily be applied to any statistical model despite our primary focus on molecular substitution models in this study. Using two phylogenetic example datasets, we demonstrate that the runtime of the marginal likelihood estimation can be reduced significantly even if only two CPUs are available (an average performance increase of 1.96x). The performance increase is nearly linear with the number of available CPUs. We record a performance increase of 13.3x for cluster nodes with 16 CPUs, representing a substantial reduction to the runtime of marginal likelihood estimations. Hence, our parallelization strategy enables the estimation of marginal likelihoods to complete in a feasible amount of time which previously needed days, weeks or even months. The methods described here are implemented in our open-source software RevBayes which is available from http://www.RevBayes.com.

A new molecular phylogeny for the Tanypodinae (Diptera: Chironomidae) places the Australian diversity in a global context

The non-biting midge subfamily Tanypodinae (Diptera: Chironomidae) is species-rich, ecologically diverse, and near-globally distributed. Within the subfamily, aspects of generic and species-level taxonomy remain poorly understood, in particular the validity of assignment of Australian and New Zealand taxa to genera erected for northern hemisphere (Holarctic) fauna. Here, we place the austral diversity within this global context by extensive geographical and taxonomic sampling in concert with a multilocus phylogenetic approach. We incorporated sequence data for mitochondrial COI, and nuclear 28S and CAD, and conducted Bayesian and maximum likelihood phylogenetic inferences and Bayesian divergence time estimation. The resolved phylogeny supported many associations of Australian taxa with their proposed Holarctic congeners, with the exception of Apsectrotanypus Fittkau, and validates several taxa as endemic. Three of four New Zealand sampled taxa had their sister groups in Australia; New Zealand Monopelopia Fittkau was sister to a German congener. This included the first record of Procladius Kieffer from New Zealand. Most nodes connecting austral and Holarctic taxa clustered around the Cretaceous-Tertiary boundary (60–80 mya), whereas New Zealand-Australia nodes were generally slightly younger (53–57 mya). Together, these data contribute substantially to our understanding of the taxonomy, systematics and biogeography of the Australian Tanypodinae and more broadly to knowledge of Australia’s aquatic insect biodiversity.

Occurrence of Alternaria arborescens causing leaf spot on Pereskia aculeata in Brazil

Pereskia aculeata, commonly known as ora-pro-nobis, is a is a highly nutritive species used for human consumption with potential of expansion. However, leaf lesions that did not match the symptoms and signs of previously described pathogenic species to this culture were observed in an ora-pro-nobis cultivation in Curitiba, Paraná state, southern Brazil. Therefore, the objectives of this study were to isolate, test the pathogenicity and characterize morphologically and molecularly the causal agent of the symptoms of this species. For this purpose, isolation from leaf fragments was performed. Pathogenicity was tested in potted cuttings which were wounded to simulate the thorns natural injuries and inoculated with a conidia suspension. The colony characteristics and morphological features of conidia were determined on V8 medium, and the molecular characterization was performed with Bayesian phylogenetic inference of the partial concatenated alignment of the ITS, gapdh and tef1 genes of the isolate. Morphological characterization matched the described for Alternaria spp. The symptoms on inoculated cuttings appeared within 4 days in all inoculated leaves, while the control plants remained asymptomatic, and 10 days after inoculation all symptomatic leaves fell from the plant. The characteristics of the re-isolated fungus from the pathogenicity test lesions were morphologically similar to the isolate from naturally infected leaves. The sequences of the pathogen grouped in a reliable clade called Alternaria arborescens. To our knowledge, this is the first report of A.arborescens causing lesions on P.aculeata worldwide. Implications of these finding should be further studied.

Redefining the Distributional Boundaries and Phylogenetic Relationships for Ctenomids From Central Argentina.

With about 68 recognized living species, subterranean rodents of the genus Ctenomys are found in a multiplicity of habitats, from the dunes of the Atlantic coast to the Andes Mountains, including environments ranging from humid steppes of Pampas to the dry deserts of Chaco region. However, this genus needs an exhaustive reevaluation of its systematic and phylogenetic relationships regarding the different groups that compose it. This knowledge is essential to propose biodiversity conservation strategies both at species level and at higher hierarchical levels. In order to clarify the taxonomy and the recent evolutionary history from populations of Ctenomys in the Pampas region, Argentina, phylogenetic relationships among them were evaluated using mitochondrial DNA sequences: gene encoding cytochrome b protein (1,140 bp) and the non-coding D-loop region (434 bp). To infer the divergence times inside the Ctenomys clade, a Bayesian calibrate tree using fossil remains data from different families within Caviomorpha was performed at first. Secondly, that calibration data was used as priors in a new Bayesian phylogenetic inference within the genus Ctenomys. This phylogenetic tree emphasized on species currently distributed on the Pampas region, more precisely considering both the talarum and mendocinus groups. Bayesian inferences (BI) were integrated with the results of a Maximum Likelihood approach (ML). Based on these results, the distributional limits of the mendocinus and talarum groups appear to be related to the physiognomy of the Pampas region soils. On the other hand, the validity of C. pundti complex as a differentiated species of C. talarum is debated. According to previous evidence from morphological and chromosomal studies, these results show a very low divergence between those species that originally were classified within the talarum group. Mitochondrial DNA sequences from populations associated with these putative species have not recovered as reciprocal monophyletic groups in the phylogenetic analyses. In conclusion, C. talarum and C. pundti complex might be considered as the same biological species, or lineages going through a recent or incipient differentiation process. The results obtained in this study have important implications for conservation policies and practices, since both species are currently categorized as Vulnerable and Endangered, respectively.

New specimens from Mid-Cretaceous Myanmar amber illuminate the phylogenetic placement of Lagonomegopidae (Arachnida: Araneae)

Abstract New lagonomegopid spiders are described from Mid-Cretaceous Myanmar (Burmese) amber. Two new genera and species based on single specimens, Scopomegops fax gen. & sp. nov. and Hiatomegops spinalis gen. & sp. nov. are described. Two specimens belonging to Lineaburmops beigeli are further described. Additionally, after re-examining the holotype of Odontomegops titan, a detailed description of its basal ventral abdomen is added here. A phylogenetic analysis was performed to investigate the phylogenetic placement of Lagonomegopidae. A matrix of 79 morphological characters, scored for six lagonomegopid taxa and 26 non-lagonomegopid taxa, was analysed through parsimony and Bayesian phylogenetic inference. Our results recover extant Palpimanoidea as a monophyletic group and partly suggest that Lagonomegopidae is the sister-group to extant Palpimanoidea. The external sexual organs, retrolateral tibial apophysis on the male palp and tracheal spiracle in lagonomegopids are discussed.

Robustness of Phylogenetic Inference to Model Misspecification Caused by Pairwise Epistasis.

Likelihood-based phylogenetic inference posits a probabilistic model of character state change along branches of a phylogenetic tree. These models typically assume statistical independence of sites in the sequence alignment. This is a restrictive assumption that facilitates computational tractability, but ignores how epistasis, the effect of genetic background on mutational effects, influences the evolution of functional sequences. We consider the effect of using a misspecified site-independent model on the accuracy of Bayesian phylogenetic inference in the setting of pairwise-site epistasis. Previous work has shown that as alignment length increases, tree reconstruction accuracy also increases. Here, we present a simulation study demonstrating that accuracy increases with alignment size even if the additional sites are epistatically coupled. We introduce an alignment-based test statistic that is a diagnostic for pairwise epistasis and can be used in posterior predictive checks.

Historical biogeography, systematics, and integrative taxonomy of the non-Ethiopian speckled pelage brush-furred rats (Lophuromys flavopunctatus group)

BackgroundThe speckled-pelage brush-furred rats (Lophuromys flavopunctatus group) have been difficult to define given conflicting genetic, morphological, and distributional records that combine to obscure meaningful accounts of its taxonomic diversity and evolution. In this study, we inferred the systematics, phylogeography, and evolutionary history of the L. flavopunctatus group using maximum likelihood and Bayesian phylogenetic inference, divergence times, historical biogeographic reconstruction, and morphometric discriminant tests. We compiled comprehensive datasets of three loci (two mitochondrial [mtDNA] and one nuclear) and two morphometric datasets (linear and geometric) from across the known range of the genus Lophuromys.ResultsThe mtDNA phylogeny supported the division of the genus Lophuromys into three primary groups with nearly equidistant pairwise differentiation: one group corresponding to the subgenus Kivumys (Kivumys group) and two groups corresponding to the subgenus Lophuromys (L. sikapusi group and L. flavopunctatus group). The L. flavopunctatus group comprised the speckled-pelage brush-furred Lophuromys endemic to Ethiopia (Ethiopian L. flavopunctatus members [ETHFLAVO]) and the non-Ethiopian ones (non-Ethiopian L. flavopunctatus members [NONETHFLAVO]) in deeply nested relationships. There were distinctly geographically structured mtDNA clades among the NONETHFLAVO, which were incongruous with the nuclear tree where several clades were unresolved. The morphometric datasets did not systematically assign samples to meaningful taxonomic units or agree with the mtDNA clades. The divergence dating and ancestral range reconstructions showed the NONETHFLAVO colonized the current ranges over two independent dispersal events out of Ethiopia in the early Pleistocene.ConclusionThe phylogenetic associations and divergence times of the L. flavopunctatus group support the hypothesis that paleoclimatic impacts and ecosystem refugia during the Pleistocene impacted the evolutionary radiation of these rodents. The overlap in craniodental variation between distinct mtDNA clades among the NONETHFLAVO suggests unraveling underlying ecomorphological drivers is key to reconciling taxonomically informative morphological characters. The genus Lophuromys requires a taxonomic reassessment based on extensive genomic evidence to elucidate the patterns and impacts of genetic isolation at clade contact zones.

Discovery of a new mammal species (Soricidae: Eulipotyphla) from Narcondam volcanic island, India

We discovered a new Crocidura species of shrew (Soricidae: Eulipotyphla) from Narcondam Island, India by using both morphological and molecular approaches. The new species, Crocidura narcondamica sp. nov. is of medium size (head and body lengths) and has a distinct external morphology (darker grey dense fur with a thick, darker tail) and craniodental characters (braincase is rounded and elevated with weak lambdoidal ridges) in comparison to other close congeners. This is the first discovery of a shrew from this volcanic island and increases the total number of Crocidura species catalogued in the Indian checklist of mammals to 12. The newly discovered species shows substantial genetic distances (12.02% to 16.61%) to other Crocidura species known from the Indian mainland, the Andaman and Nicobar Archipelago, Myanmar, and from Sumatra. Both Maximum-Likelihood and Bayesian phylogenetic inferences, based on mitochondrial (cytochrome b) gene sequences showed distinct clustering of all included soricid species and exhibit congruence with the previous evolutionary hypothesis on this mammalian group. The present phylogenetic analyses also furnished the evolutionary placement of the newly discovered species within the genus Crocidura.

Population genetic structure of sharpbelly Hemiculter leucisculus (Basilesky, 1855) and morphological diversification along climate gradients in China.

Sharpbelly Hemiculter leucisculus (Basilewski, 1855) is a small, widespread, and native cyprinid fish with prominent habitat suitability and high invasive potential and is becoming the dominant species in freshwater ecosystems under intensified environmental disturbances. But how H. leucisculus acclimates to extremely heterogeneous environments remains unclear. In current study, the genetic structure of H. leucisculus was analyzed using Bayesian phylogenetic inference, haplotype network, and STRUCTURE base on cytb gene across 18 populations spanning 20 degrees of latitude and 18 degrees of longitude in China. The morphological diversification of body size and shape for H. leucisculus along the climate gradient was studied. The results showed that the 18 H. leucisculus populations were divided into 3 clusters: one cluster mainly from Huanghe River Basin, another cluster mainly from Yangzi River Basin, and H cluster containing Hainan and Beihai populations. The fish from southern populations were deeper bodied while individuals from northern populations were more slender. Inland individuals were more streamlined while coastal individuals were of deeper body. The partial Mantel test predicts that the potential mechanism underlining the intraspecies morphological diversification along climate gradients is primarily the divergent selection pressures among different environments, while genetic variation had less contribution to morphological differentiation. The formation of the Nanling Mountain Range could drive genetic differentiation between Beihai population and those from Yangzi River Basin. The present results highlight strong selective pressures of climate on widespread species and enrich morphological differentiation basis of acclimation for species with high habitat suitability and invasive potential.

Adaptive dating and fast proposals: Revisiting the phylogenetic relaxed clock model.

Relaxed clock models enable estimation of molecular substitution rates across lineages and are widely used in phylogenetics for dating evolutionary divergence times. Under the (uncorrelated) relaxed clock model, tree branches are associated with molecular substitution rates which are independently and identically distributed. In this article we delved into the internal complexities of the relaxed clock model in order to develop efficient MCMC operators for Bayesian phylogenetic inference. We compared three substitution rate parameterisations, introduced an adaptive operator which learns the weights of other operators during MCMC, and we explored how relaxed clock model estimation can benefit from two cutting-edge proposal kernels: the AVMVN and Bactrian kernels. This work has produced an operator scheme that is up to 65 times more efficient at exploring continuous relaxed clock parameters compared with previous setups, depending on the dataset. Finally, we explored variants of the standard narrow exchange operator which are specifically designed for the relaxed clock model. In the most extreme case, this new operator traversed tree space 40% more efficiently than narrow exchange. The methodologies introduced are adaptive and highly effective on short as well as long alignments. The results are available via the open source optimised relaxed clock (ORC) package for BEAST 2 under a GNU licence (https://github.com/jordandouglas/ORC).

The inner ear anatomy of glyptodonts and pampatheres (Xenarthra, Cingulata): Functional and phylogenetic implications

In recent years, the increase in studies on the inner ear anatomy of xenarthrans provided new insights regarding some locomotor and phylogenetic aspects. These works have begun to include fossil specimens, although so far, only two extinct representatives of the suborder Folivora (ground sloths) have been analyzed. In the present study, we present the first reconstruction and description of the bony labyrinth of five extinct specimens of the suborder Cingulata, the glyptodonts Glyptodon, Doedicurus, Panochthus, and Pseudoplohophorus, as well as the pampathere Holmesina, in order to analyze their morphology and evaluate functional and phylogenetic aspects. The anatomical and morphometric comparison shows that the inner ear of glyptodonts is largely similar to that of modern armadillos, particularly sharing many characteristics with Chlamyphorus. While the inner ear of Holmesina is similar to that of glyptodonts, it presents some characteristics that make it more similar to extant armadillos. Regarding some functional aspects inferred from the anatomy of the semicircular canals, glyptodonts would have had low levels of agility, although similar to most extant armadillos, and low head angular velocity. On the other hand, Holmesina would have had an agility level similar to that of the giant armadillo Priodontes, although it would also present low head angular velocity as seen in most extant armadillos. Finally, preliminary parsimony and Bayesian phylogenetic inference analyses based on a reduced set of inner ear characters showed interesting results, with glyptodonts and Holmesina recovered as closely related to Chlamyphorus. These results closely resemble the recent topologies obtained using ancient DNA and markedly contrast with the topologies obtained using more traditional morphological characters, highlighting the importance of inner ear data for the inference of phylogenetic relationships in cases where DNA preservation in fossil taxa is unaccessible.

Adaptive Tree Proposals for Bayesian Phylogenetic Inference.

Bayesian inference of phylogeny with Markov chain Monte Carlo plays a key role in the study of evolution. Yet, this method still suffers from a practical challenge identified more than two decades ago: designing tree topology proposals that efficiently sample tree spaces. In this article, I introduce the concept of adaptive tree proposals for unrooted topologies, that is, tree proposals adapting to the posterior distribution as it is estimated. I use this concept to elaborate two adaptive variants of existing proposals and an adaptive proposal based on a novel design philosophy in which the structure of the proposal is informed by the posterior distribution of trees. I investigate the performance of these proposals by first presenting a metric that captures the performance of each proposal within a mixture of proposals. Using this metric, I compare the performance of the adaptive proposals to the performance of standard and parsimony-guided proposals on 11 empirical data sets. Using adaptive proposals led to consistent performance gains and resulted in up to 18-fold increases in mixing efficiency and 6-fold increases in convergence rate without increasing the computational cost of these analyses. [Bayesian phylogenetic inference; Markov chain Monte Carlo; posterior probability distribution; tree proposals.]

Molecular species delimitation reveals high diversity in the mosquito Anopheles tessellatus Theobald, 1901 (Diptera, Culicidae) across its range

Anopheles tessellatus is a potentially important vector found across South, East and Southeast Asia. While it was formerly considered a formidable vector of human Plasmodium and filarial parasites in the Maldives, and of lesser importance as a vector of human Plasmodium in Sri Lanka and parts of Indonesia, it is currently of little or unknown health importance in many other parts of its range. This study describes the genetic diversity and evolutionary relationships among An. tessellatus populations in nine Asian countries at the COI gene using maximum-likelihood and Bayesian phylogenetic inference tree and cluster-based species delimitation approaches. These analyses reveal exceptional levels of genetic diversity in An. tessellatus populations across its known range, and identify up to six putative species in the newly determined Tessellatus Complex. The existence of such cryptic diversity has potentially important consequences for vector management and disease control. Differences in the ecologies and life histories among these species may have considerable impact on vectorial capacity and may go some way towards explaining why An. tessellatus s.l. has such varying degrees of public health importance across its range.

The Enigmatic Thelebolaceae (Thelebolales, Leotiomycetes): One New Genus Solomyces and Five New Species.

The family Thelebolaceae belongs to the order Thelebolales, class Leotiomycetes, and contains 22 genera. In this study, we introduce a new genus Solomyces gen. nov. in the family Thelebolaceae, which is supported by morphological observation and multilocus-based [internal transcribed spacers (ITS) + LSU and ITS + LSU+ MCM7+ EF1A+ RPB2] phylogenetic analysis. Maximum-likelihood and Bayesian phylogenetic inference analyses indicated that Solomyces is a distinct genus within this family. The new genus is compared against related Thelebolaceae genera, and its description and illustration are provided. This genus comprises one new species and one unnamed species (including two strains). We also report the addition of four new species – Pseudogymnoascus shaanxiensis, Pseudogymnoascus guizhouensis, Pseudogymnoascus sinensis, and Geomyces obovatus – in the family Thelebolaceae and present their morphological and phylogenetic characterizations.

Quantifying the impact of an inference model in Bayesian phylogenetics

Abstract Phylogenetic trees are currently routinely reconstructed from an alignment of character sequences (usually nucleotide sequences). Bayesian tools, such as MrBayes, RevBayes and BEAST2, have gained much popularity over the last decade, as they allow joint estimation of the posterior distribution of the phylogenetic trees and the parameters of the underlying inference model. An important ingredient of these Bayesian approaches is the species tree prior. In principle, the Bayesian framework allows for comparing different tree priors, which may elucidate the macroevolutionary processes underlying the species tree. In practice, however, only macroevolutionary models that allow for fast computation of the prior probability are used. The question is how accurate the tree estimation is when the real macroevolutionary processes are substantially different from those assumed in the tree prior. Here we present pirouette, a free and open‐source r package that assesses the inference error made by Bayesian phylogenetics for a given macroevolutionary diversification model. pirouette makes use of BEAST2, but its philosophy applies to any Bayesian phylogenetic inference tool. We describe pirouette’s usage providing full examples in which we interrogate a model for its power to describe another. Last, we discuss the results obtained by the examples and their interpretation.

Particle Gibbs sampling for Bayesian phylogenetic inference.

The combinatorial sequential Monte Carlo (CSMC) has been demonstrated to be an efficient complementary method to the standard Markov chain Monte Carlo (MCMC) for Bayesian phylogenetic tree inference using biological sequences. It is appealing to combine the CSMC and MCMC in the framework of the particle Gibbs (PG) sampler to jointly estimate the phylogenetic trees and evolutionary parameters. However, the Markov chain of the PG may mix poorly for high dimensional problems (e.g. phylogenetic trees). Some remedies, including the PG with ancestor sampling and the interacting particle MCMC, have been proposed to improve the PG. But they either cannot be applied to or remain inefficient for the combinatorial tree space. We introduce a novel CSMC method by proposing a more efficient proposal distribution. It also can be combined into the PG sampler framework to infer parameters in the evolutionary model. The new algorithm can be easily parallelized by allocating samples over different computing cores. We validate that the developed CSMC can sample trees more efficiently in various PG samplers via numerical experiments. The implementation of our method and the data underlying this articleare available at https://github.com/liangliangwangsfu/phyloPMCMC. Supplementary data are available at Bioinformatics online.

Phylogenetic inference enables reconstruction of a long-overlooked outbreak of almond leaf scorch disease (Xylella fastidiosa) in Europe

The recent introductions of the bacterium Xylella fastidiosa (Xf) into Europe are linked to the international plant trade. However, both how and when these entries occurred remains poorly understood. Here, we show how almond scorch leaf disease, which affects ~79% of almond trees in Majorca (Spain) and was previously attributed to fungal pathogens, was in fact triggered by the introduction of Xf around 1993 and subsequently spread to grapevines (Pierceʼs disease). We reconstructed the progression of almond leaf scorch disease by using broad phylogenetic evidence supported by epidemiological data. Bayesian phylogenetic inference predicted that both Xf subspecies found in Majorca, fastidiosa ST1 (95% highest posterior density, HPD: 1990–1997) and multiplex ST81 (95% HPD: 1991–1998), shared their most recent common ancestors with Californian Xf populations associated with almonds and grapevines. Consistent with this chronology, Xf-DNA infections were identified in tree rings dating to 1998. Our findings uncover a previously unknown scenario in Europe and reveal how Pierce’s disease reached the continent.

Regional drivers of diversification in the late Quaternary in a widely distributed generalist species, the common pheasant Phasianus colchicus

AbstractAimPleistocene climate and associated environmental changes have influenced phylogeographic patterns of many species. These not only depend on a species’ life history but also vary regionally. Consequently, populations of widespread species that occur in several biomes might display different evolutionary trajectories. We aimed to identify regional drivers of diversification in the common pheasant, a widely distributed ecological generalist.LocationAsia.TaxonCommon pheasant Phasianus colchicus.MethodsUsing a comprehensive geographical sampling of 204 individuals from the species’ entire range genotyped at seven nuclear and two mitochondrial loci, we reconstructed spatio‐temporal diversification and demographic history of the common pheasant. We applied Bayesian phylogenetic inference to describe phylogeographic structure, generated a species tree and inferred demographic history within and migration between lineages. Moreover, to establish a taxonomic framework, we conducted a species delimitation analysis.ResultsThe common pheasant diversified during the Late Pleistocene into eight distinct lineages. It originated at the edge of the Qinghai–Tibetan plateau and spread to East and Central Asia. Only the widely distributed lowland lineage of East Asia displayed recent range expansion. Greater phylogeographic structure was identified elsewhere, with lineages showing no sign of recent demographic changes. One lineage in south‐central China is the result of long‐term isolation within a climatically stable but topographically complex region. In lineages from arid Central Asia and China, range expansions were impeded by repeated population fragmentation during dry glacial periods and by recent aridification.Main conclusionsSpatio‐temporal phylogeographic frameworks of widespread taxa such as the common pheasant provide valuable opportunities to identify divergent drivers of regional diversification. Our results suggest that diversification and population histories in the eight distinct evolutionary lineages were shaped by regionally variable effects of past climate and associated environmental changes. The evolutionary history of the common pheasant is best reflected by its being split into three species.

A Bayesian phylogenetic hidden Markov model for B cell receptor sequence analysis.

The human body generates a diverse set of high affinity antibodies, the soluble form of B cell receptors (BCRs), that bind to and neutralize invading pathogens. The natural development of BCRs must be understood in order to design vaccines for highly mutable pathogens such as influenza and HIV. BCR diversity is induced by naturally occurring combinatorial “V(D)J” rearrangement, mutation, and selection processes. Most current methods for BCR sequence analysis focus on separately modeling the above processes. Statistical phylogenetic methods are often used to model the mutational dynamics of BCR sequence data, but these techniques do not consider all the complexities associated with B cell diversification such as the V(D)J rearrangement process. In particular, standard phylogenetic approaches assume the DNA bases of the progenitor (or “naive”) sequence arise independently and according to the same distribution, ignoring the complexities of V(D)J rearrangement. In this paper, we introduce a novel approach to Bayesian phylogenetic inference for BCR sequences that is based on a phylogenetic hidden Markov model (phylo-HMM). This technique not only integrates a naive rearrangement model with a phylogenetic model for BCR sequence evolution but also naturally accounts for uncertainty in all unobserved variables, including the phylogenetic tree, via posterior distribution sampling.

The complete mitochondrial genome of Quasipaa exilispinosa (Anura: Dicroglossidae)

The complete mitochondrial genome (mitogenome) of Quasipaa exilispinosa (Anura: Dicroglossidae) we sequenced can provide critical information for mitogenome evolution. This mitogenome with 17,046 bp long contains 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and one non-coding regions. The overall AT content of 58.6%. The data using Bayesian phylogenetic inference supports the close relationship with Q. spinosa from the genus Quasipaa. Our results will be helpful for detail study on phylogenetic relationships of the related taxa.