Challenges of phylogenetic analyses of aDNA sequences

One of the crucial steps of authentication of aDNA sequences is phylogenetic consistency. Amplified sequences should fit into the phylogenetic framework of their supposed origin. An inherent property of aDNA sequences however, is their short sequence length. Additionally, genes for aDNA studies are often chosen by their preservation potential rather than by phylogenetically informative content. This poses potential challenges regarding their analyses, and might result in an inaccurate reflection of the supposed phylogenetic history of the sequence or organism under study. In this paper some fundamental problems of phylogenetic analysis and interpretation of aDNA datasets are discussed. Suggestions for character sampling and treatment of missing data are made. The publication is the result of a talk from the 1st PAMINSA Meeting in Rio de Janeiro, July 2005.

Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism.

The tendency for more closely related species to share similar traits and ecological strategies can be explained by their longer shared evolutionary histories and represents phylogenetic conservatism. How strongly species traits co-vary with phylogeny can significantly impact how we analyze cross-species data and can influence our interpretation of assembly rules in the rapidly expanding field of community phylogenetics. Phylogenetic conservatism is typically quantified by analyzing the distribution of species values on the phylogenetic tree that connects them. Many phylogenetic approaches, however, assume a completely sampled phylogeny: while we have good estimates of deeper phylogenetic relationships for many species-rich groups, such as birds and flowering plants, we often lack information on more recent interspecific relationships (i.e., within a genus). A common solution has been to represent these relationships as polytomies on trees using taxonomy as a guide. Here we show that such trees can dramatically inflate estimates of phylogenetic conservatism quantified using S. P. Blomberg et al.'s K statistic. Using simulations, we show that even randomly generated traits can appear to be phylogenetically conserved on poorly resolved trees. We provide a simple rarefaction-based solution that can reliably retrieve unbiased estimates of K, and we illustrate our method using data on first flowering times from Thoreau's woods (Concord, Massachusetts, USA).

Inteins in pathogenic fungi: a phylogenetic tool and perspectives for therapeutic applications

Inteins or "internal proteins" are coding sequences that are transcribed and translated with flanking sequences (exteins). After translation, the inteins are excised by an autocatalytic process and the host protein assumes its normal conformation and develops its expected function. These parasitic genetic elements have been found in important, conserved proteins in all three domains of life. Most of the eukaryotic inteins are present in the fungi kingdom and the PRP8 intein is one of the most widespread inteins, occurring in important pathogens such as Cryptococcus neoformans (varieties grubii and neoformans), Cryptococcus gattii, Histoplasma capsulatum and Paracoccidioides brasiliensis. The knowledge of conserved and non-conserved domains in inteins have opened up new opportunities for the study of population variability in pathogenic fungi, including their phylogenetic relationships and recognition or diagnoses of species. Furthermore, inteins in pathogenic fungi should also be considered a promising therapeutic drug target, since once the autocatalytic splicing is inhibited, the host protein, which is typically vital, will not be able to perform its normal function and the fungal cell will not survive or reproduce.

Dengue virus type 3 in Brazil: a phylogenetic perspective

Circulation of a new dengue virus (DENV)-3 genotype was recently described in Brazil and Colombia, but the precise classification of this genotype has been controversial. Here we perform phylogenetic and nucleotide-distance analyses of the envelope gene, which support the subdivision of DENV-3 strains into five distinct genotypes (GI to GV) and confirm the classification of the new South American genotype as GV. The extremely low genetic distances between Brazilian GV strains and the prototype Philippines/L11423 GV strain isolated in 1956 raise important questions regarding the origin of GV in South America.

Phylogenetic plant community structure along elevation is lineage specific

The trend of closely related taxa to retain similar environmental preferences mediated by inherited traits suggests that several patterns observed at the community scale originate from longer evolutionary processes. While the effects of phylogenetic relatedness have been previously studied within a single genus or family, lineage-specific effects on the ecological processes governing community assembly have rarely been studied for entire communities or flora. Here, we measured how community phylogenetic structure varies across a wide elevation gradient for plant lineages represented by thirty-five families, using a co-occurrence index and net relatedness index (NRI). We propose a framework that analyses each lineage separately and reveals the trend of ecological assembly at tree nodes. We found prevailing phylogenetic clustering for more ancient nodes and overdispersion in more recent tree nodes. Closely related species may thus rapidly evolve new environmental tolerances to radiate into distinct communities, while older lineages likely retain inherent environmental tolerances to occupy communities in similar environments, either through efficient dispersal mechanisms or the exclusion of older lineages with more divergent environmental tolerances. Our study illustrates the importance of disentangling the patterns of community assembly among lineages to better interpret the ecological role of traits. It also sheds light on studies reporting absence of phylogenetic signal, and opens new perspectives on the analysis of niche and trait conservatism across lineages.

Phylogenetic reconstruction of transmission events from individuals with acute HIV infection: toward more-rigorous epidemiological definitions.

Phylogenetic reconstructions of transmission events from individuals with acute human immunodeficiency virus (HIV) infection are conducted to illustrate this group's heightened infectivity. Varied definitions of acute infection and assumptions about observed phylogenetic clusters may produce misleading results. We conducted a phylogenetic analysis of HIV pol sequences from 165 European patients with estimated infection dates and calculated the difference between dates within clusters. Nine phylogenetic clusters were observed. Comparison of dates within clusters revealed that only 2 could have been generated during acute infection. Previous analyses may have incorrectly assigned transmission events to the acutely HIV infected when they were more likely to have occurred during chronic infection.

Phylogenetic analysis of Biomphalaria tenagophila (Orbigny, 1835) (Mollusca: Gastropoda)

Mitochondrial DNA of Biomphalaria tenagophila, a mollusc intermediate host of Schistosoma mansoni in Brazil, was sequenced and characterised. The genome size found for B. tenagophila was 13,722 bp and contained 13 messenger RNAs, 22 transfer RNAs (tRNA) and two ribosomal RNAs (rRNA). In addition to sequencing, the mitochondrial DNA (mtDNA) genome organization of B. tenagophila was analysed based on its content and localization of both coding and non-coding regions, regions of gene overlap and tRNA nucleotide sequences. Sequences of protein, rRNA 12S and rRNA 16S nucleotides as well as gene organization were compared between B. tenagophila and Biomphalaria glabrata, as the latter is the most important S. mansoni intermediate host in Brazil. Differences between such species were observed regarding rRNA composition. The complete sequence of the B. tenagophila mitochondrial genome was deposited in GenBank (accession EF433576). Furthermore, phylogenetic relationships were estimated among 28 mollusc species, which had their complete mitochondrial genome deposited in GenBank, using the neighbour-joining method, maximum parsimony and maximum likelihood bootstrap. B. tenagophila was positioned at a branch close to B. glabrata and Pulmonata molluscs, collectively comprising a paraphyletic group, contrary to Opistobranchia, which was positioned at a single branch and constituted a monophyletic group.

Phylogenetic alpha and beta diversities of butterfly communities correlate with climate in the western Swiss Alps

The observation of non-random phylogenetic distribution of traits in communities provides evidence for niche-based community assembly. Environment may influence the phylogenetic structure of communities because traits determining how species respond to prevailing conditions can be phylogenetically conserved. In this study, we investigate the variation of butterfly species richness and of phylogenetic - and -diversities along temperature and plant species richness gradients. Our study indicates that butterfly richness is independently positively correlated to temperature and plant species richness in the study area. However, the variation of phylogenetic - and -diversities is only correlated to temperature. The significant phylogenetic clustering at high elevation suggests that cold temperature filters butterfly lineages, leading to communities mostly composed of closely related species adapted to those climatic conditions. These results suggest that in colder and more severe conditions at high elevations deterministic processes and not purely stochastic events drive the assemblage of butterfly communities.

Inference of distributional parameters from compositional samples containing nondetects

Low concentrations of elements in geochemical analyses have the peculiarity of beingcompositional data and, for a given level of significance, are likely to be beyond thecapabilities of laboratories to distinguish between minute concentrations and completeabsence, thus preventing laboratories from reporting extremely low concentrations of theanalyte. Instead, what is reported is the detection limit, which is the minimumconcentration that conclusively differentiates between presence and absence of theelement. A spatially distributed exhaustive sample is employed in this study to generateunbiased sub-samples, which are further censored to observe the effect that differentdetection limits and sample sizes have on the inference of population distributionsstarting from geochemical analyses having specimens below detection limit (nondetects).The isometric logratio transformation is used to convert the compositional data in thesimplex to samples in real space, thus allowing the practitioner to properly borrow fromthe large source of statistical techniques valid only in real space. The bootstrap method isused to numerically investigate the reliability of inferring several distributionalparameters employing different forms of imputation for the censored data. The casestudy illustrates that, in general, best results are obtained when imputations are madeusing the distribution best fitting the readings above detection limit and exposes theproblems of other more widely used practices. When the sample is spatially correlated, itis necessary to combine the bootstrap with stochastic simulation

DNA multigene sequencing of topotypic specimens of the fascioliasis vector Lymnaea diaphana and phylogenetic analysis of the genus Pectinidens (Gastropoda)

Freshwater lymnaeid snails are crucial in defining transmission and epidemiology of fascioliasis. In South America, human endemic areas are related to high altitudes in Andean regions. The species Lymnaea diaphana has, however, been involved in low altitude areas of Chile, Argentina and Peru where human infection also occurs. Complete nuclear ribosomal DNA 18S, internal transcribed spacer (ITS)-2 and ITS-1 and fragments of mitochondrial DNA 16S and cytochrome c oxidase (cox)1 genes of L. diaphana specimens from its type locality offered 1,848, 495, 520, 424 and 672 bp long sequences. Comparisons with New and Old World Galba/Fossaria, Palaearctic stagnicolines, Nearctic stagnicolines, Old World Radix and Pseudosuccinea allowed to conclude that (i) L. diaphana shows sequences very different from all other lymnaeids, (ii) each marker allows its differentiation, except cox1 amino acid sequence, and (iii) L. diaphana is not a fossarine lymnaeid, but rather an archaic relict form derived from the oldest North American stagnicoline ancestors. Phylogeny and large genetic distances support the genus Pectinidens as the first stagnicoline representative in the southern hemisphere, including colonization of extreme world regions, as most southern Patagonia, long time ago. The phylogenetic link of L. diaphana with the stagnicoline group may give light to the aforementioned peculiar low altitude epidemiological scenario of fascioliasis.

Cardenolides in nectar may be more than a consequence of allocation to other plant parts: a phylogenetic study of Asclepias

The primary function of secondary plant metabolites is thought to be defence against herbivores. The frequent occurrence of these same noxious compounds in floral nectar, which functions primarily to attract pollinators, has been seen as paradoxical. Although these compounds may have an adaptive purpose in nectar, they may also occur as a nonadaptive consequence of chemical defence in other plant parts. If nectar chemistry reflects physiological constraints or passive leakage from other tissues, we expect that the identity and relative concentration of nectar cardenolides to be correlated with those of other plant parts; in contrast, discordant distributions of compounds in nectar and other tissues may suggest adaptive roles in nectar. We compared the concentrations and identities of cardenolides in the nectar, leaves and flowers of 12 species from a monophyletic clade of Asclepias. To measure putative toxicity of nectar cardenolides, we then examined the effects of a standard cardenolide (digoxin) on the behaviour of bumblebees, a common generalist pollinator of Asclepias. We found that the average cardenolide concentrations in nectar, leaves and flowers of the 12 Asclepias species were positively correlated as predicted by nonadaptive hypotheses. However, significant differences in the identities and concentrations of individual cardenolides between nectar and leaves suggest that the production or allocation of cardenolides may be independently regulated at each plant part. In addition, cardenolide concentrations in leaves and nectar exhibited no phylogenetic signal. Surprisingly, bumblebees did not demonstrate an aversion to digoxin-rich nectar, which may indicate that nectar cardenolides have little effect on pollination. Although the idea that discordant patterns of secondary metabolites across tissue types may signal adaptive functions is attractive, there is evidence to suggest constraint contributes to nectar secondary chemistry. Further work testing the ecological impacts of such patterns will be critical in determining the functional significance of nectar cardenolides

Applications of stochastic analysis in finance and statistical inference

First: A continuous-time version of Kyle's model (Kyle 1985), known as the Back's model (Back 1992), of asset pricing with asymmetric information, is studied. A larger class of price processes and of noise traders' processes are studied. The price process, as in Kyle's model, is allowed to depend on the path of the market order. The process of the noise traders' is an inhomogeneous Lévy process. Solutions are found by the Hamilton-Jacobi-Bellman equations. With the insider being risk-neutral, the price pressure is constant, and there is no equilibirium in the presence of jumps. If the insider is risk-averse, there is no equilibirium in the presence of either jumps or drifts. Also, it is analised when the release time is unknown. A general relation is established between the problem of finding an equilibrium and of enlargement of filtrations. Random announcement time is random is also considered. In such a case the market is not fully efficient and there exists equilibrium if the sensitivity of prices with respect to the global demand is time decreasing according with the distribution of the random time. Second: Power variations. it is considered, the asymptotic behavior of the power variation of processes of the form _integral_0^t u(s-)dS(s), where S_ is an alpha-stable process with index of stability 0&alpha&2 and the integral is an Itô integral. Stable convergence of corresponding fluctuations is established. These results provide statistical tools to infer the process u from discrete observations. Third: A bond market is studied where short rates r(t) evolve as an integral of g(t-s)sigma(s) with respect to W(ds), where g and sigma are deterministic and W is the stochastic Wiener measure. Processes of this type are particular cases of ambit processes. These processes are in general not of the semimartingale kind.

Inference about the number of contributors to a DNA mixture: Comparative analyses of a Bayesian network approach and the maximum allele count method.

In the forensic examination of DNA mixtures, the question of how to set the total number of contributors (N) presents a topic of ongoing interest. Part of the discussion gravitates around issues of bias, in particular when assessments of the number of contributors are not made prior to considering the genotypic configuration of potential donors. Further complication may stem from the observation that, in some cases, there may be numbers of contributors that are incompatible with the set of alleles seen in the profile of a mixed crime stain, given the genotype of a potential contributor. In such situations, procedures that take a single and fixed number contributors as their output can lead to inferential impasses. Assessing the number of contributors within a probabilistic framework can help avoiding such complication. Using elements of decision theory, this paper analyses two strategies for inference on the number of contributors. One procedure is deterministic and focuses on the minimum number of contributors required to 'explain' an observed set of alleles. The other procedure is probabilistic using Bayes' theorem and provides a probability distribution for a set of numbers of contributors, based on the set of observed alleles as well as their respective rates of occurrence. The discussion concentrates on mixed stains of varying quality (i.e., different numbers of loci for which genotyping information is available). A so-called qualitative interpretation is pursued since quantitative information such as peak area and height data are not taken into account. The competing procedures are compared using a standard scoring rule that penalizes the degree of divergence between a given agreed value for N, that is the number of contributors, and the actual value taken by N. Using only modest assumptions and a discussion with reference to a casework example, this paper reports on analyses using simulation techniques and graphical models (i.e., Bayesian networks) to point out that setting the number of contributors to a mixed crime stain in probabilistic terms is, for the conditions assumed in this study, preferable to a decision policy that uses categoric assumptions about N.