Molecular markers and human history: A tale of two haplogroup systems

To chronicle demographic movement across African Asian corridors, a variety of molecular (sequence analysis, restriction mapping and denaturing high performance liquid chromatography etc.) and statistical (correspondence analysis, AMOVA, calculation of diversity indices and phylogenetic inference, etc.) techniques were employed to assess the phylogeographic patterns of mtDNA control region and Y chromosomal variation among 14 sub-Saharan, North African and Middle Eastern populations. The patterns of genetic diversity revealed evidence of multiple migrations across several African Asian passageways as well within the African continent itself. The two-part analysis uncovered several interesting results which include the following: (1) a north (Egypt and Middle East Asia) to south (sub-Saharan Africa) partitioning of both mtDNA and Y chromosomal haplogroup diversity, (2) a genetic diversity gradient in sub-Saharan Africa from east to west, (3) evidence in favor of the Levantine Corridor over the Horn of Africa as the major genetic conduit since the Last Glacial Maximum, (4) a substantially higher mtDNA versus Y chromosomal sub-Saharan component in the Middle East collections, (5) a higher representation of East versus West African mtDNA haplotypes in the Arabian Peninsula populations versus no such bias in the Levant groups and lastly, (6) genetic remnants of the Bantu demographic expansion in sub-Saharan Africa. ^

Algorithmes de construction et correction d'arbres de gènes par la réconciliation

Les gènes, qui servent à encoder les fonctions biologiques des êtres vivants, forment l'unité moléculaire de base de l'hérédité. Afin d'expliquer la diversité des espèces que l'on peut observer aujourd'hui, il est essentiel de comprendre comment les gènes évoluent. Pour ce faire, on doit recréer le passé en inférant leur phylogénie, c'est-à-dire un arbre de gènes qui représente les liens de parenté des régions codantes des vivants. Les méthodes classiques d'inférence phylogénétique ont été élaborées principalement pour construire des arbres d'espèces et ne se basent que sur les séquences d'ADN. Les gènes sont toutefois riches en information, et on commence à peine à voir apparaître des méthodes de reconstruction qui utilisent leurs propriétés spécifiques. Notamment, l'histoire d'une famille de gènes en terme de duplications et de pertes, obtenue par la réconciliation d'un arbre de gènes avec un arbre d'espèces, peut nous permettre de détecter des faiblesses au sein d'un arbre et de l'améliorer. Dans cette thèse, la réconciliation est appliquée à la construction et la correction d'arbres de gènes sous trois angles différents: 1) Nous abordons la problématique de résoudre un arbre de gènes non-binaire. En particulier, nous présentons un algorithme en temps linéaire qui résout une polytomie en se basant sur la réconciliation. 2) Nous proposons une nouvelle approche de correction d'arbres de gènes par les relations d'orthologie et paralogie. Des algorithmes en temps polynomial sont présentés pour les problèmes suivants: corriger un arbre de gènes afin qu'il contienne un ensemble d'orthologues donné, et valider un ensemble de relations partielles d'orthologie et paralogie. 3) Nous montrons comment la réconciliation peut servir à "combiner'' plusieurs arbres de gènes. Plus précisément, nous étudions le problème de choisir un superarbre de gènes selon son coût de réconciliation.

Crassaminicella profunda gen. nov., sp nov., an anaerobic marine bacterium isolated from deep-sea sediments

A novel, anaerobic, chemo-organotrophic bacterium, designated strain Ra1766HT, was isolated from sediments of the Guaymas basin (Gulf of California, Mexico) taken from a depth of 2002 m. Cells were thin, motile, Gram-stain-positive, flexible rods forming terminal endospores. Strain Ra1766H(T) grew at temperatures of 25-45 degrees C (optimum 30 degrees C), pH 6.7-8.1 (optimum 7.5) and in a salinity of 5-60 g l(-1) NaCl (optimum 30 g l(-1)). It was an obligate heterotrophic bacterium fermenting carbohydrates (glucose and mannose) and organic acids (pyruvate and succinate). Casamino acids and amino acids (glutamate, aspartate and glycine) were also fermented. The main end products from glucose fermentation were acetate, butyrate, ethanol, H-2 and CO2. Sulfate, sulfite, thiosulfate, elemental sulfur, fumarate, nitrate, nitrite and Fe(III) were not used as terminal electron acceptors. The predominant cellular fatty acids were C-14 : 0, C-16:1 omega 7, C-16:1 omega 7 DMA and C-16:0. The main polar lipids consisted of phosphatidylglycerol, diphosphatidylglycerol, phosphatidylethanolamine and phospholipids. The G +C content of the genomic DNA was 33.7 molo/o. Phylogenetic analysis of the 16S rRNA gene sequence indicated that strain Ra1766H(T) was affiliated to cluster XI of the order Clostridia les, phylum Firmicutes. The closest phylogenetic relative of Ra1766H(T) was Geosporobacter subterraneus (94.2% 16S rRNA gene sequence similarity). On the basis of phylogenetic inference and phenotypic properties, strain Ra1766H(T) (=DSM 27501(T)=JCM 19377(T)) is proposed to be the type strain of a novel species of a novel genus, named Crassaminicella pro funda.

Algorithmes de construction et correction d'arbres de gènes par la réconciliation

Les gènes, qui servent à encoder les fonctions biologiques des êtres vivants, forment l'unité moléculaire de base de l'hérédité. Afin d'expliquer la diversité des espèces que l'on peut observer aujourd'hui, il est essentiel de comprendre comment les gènes évoluent. Pour ce faire, on doit recréer le passé en inférant leur phylogénie, c'est-à-dire un arbre de gènes qui représente les liens de parenté des régions codantes des vivants. Les méthodes classiques d'inférence phylogénétique ont été élaborées principalement pour construire des arbres d'espèces et ne se basent que sur les séquences d'ADN. Les gènes sont toutefois riches en information, et on commence à peine à voir apparaître des méthodes de reconstruction qui utilisent leurs propriétés spécifiques. Notamment, l'histoire d'une famille de gènes en terme de duplications et de pertes, obtenue par la réconciliation d'un arbre de gènes avec un arbre d'espèces, peut nous permettre de détecter des faiblesses au sein d'un arbre et de l'améliorer. Dans cette thèse, la réconciliation est appliquée à la construction et la correction d'arbres de gènes sous trois angles différents: 1) Nous abordons la problématique de résoudre un arbre de gènes non-binaire. En particulier, nous présentons un algorithme en temps linéaire qui résout une polytomie en se basant sur la réconciliation. 2) Nous proposons une nouvelle approche de correction d'arbres de gènes par les relations d'orthologie et paralogie. Des algorithmes en temps polynomial sont présentés pour les problèmes suivants: corriger un arbre de gènes afin qu'il contienne un ensemble d'orthologues donné, et valider un ensemble de relations partielles d'orthologie et paralogie. 3) Nous montrons comment la réconciliation peut servir à "combiner'' plusieurs arbres de gènes. Plus précisément, nous étudions le problème de choisir un superarbre de gènes selon son coût de réconciliation.

Protein markers of B. xylophilus Steiner & Buhrer, 1934 (Nickle, 1970) populations using quantitative proteomics and character compatibility.

The Pine Wood Nematode (PWN) Bursaphelenchus xylophilus is a severe forest pathogen in countries where it has been introduced and is considered a worldwide quarantine organism. In this study, protein markers for differentiating populations of this nematode were identified by studying differences among four selected Iberian and one American population. These populations were compared by quantitative proteomics (iTRAQ). From a total of 2860 proteins identified using the public database from the B. xylophilus genome project, 216 were unambiguous and significantly differentially regulated in the studied populations. Comparisons of their pairwise ratio were statistically treated and supported in order to convert them into discrete character states, suggesting that 141 proteins were not informative as population specific markers. Application of the Character Compatibility methodology on the remaining 75 proteins (belonging to families with different biological functions) excludes 27 which are incompatible among them. Considering only the compatible proteins, the method selects a subset of 30 specific unique protein markers which allowed the compared classification of the Iberian isolates. This approach makes it easier search for diagnostic tools and phylogenetic inference within species and populations of a pathogen exhibiting a high level of genetic diversity.

Phylogenetic relationships of Danio within the order cypriniformes: A framework for comparative and evolutionary studies of a model species

The evolutionary relationships of species of Danio and the monophyly and phylogenetic placement of the genus within the family Cyprinidae and subfamily Rasborinae provide fundamentally important phyloinformatics necessary for direct evaluations of an array of pertinent questions in modern comparative biology. Although the genus Danio is not one of the most diverse within the family, Danio rerio is one of the most important model species in biology. Many investigations have used this species or presumed close relatives to address specific questions that have lasting impact on the hypothesis and theory of development in vertebrates. Largely lacking from this approach has been a holistic picture of the exact phylogenetic or evolutionary relationships of this species and its close relatives. One thing that has been learned over the previous century is that many organismal attributes (e.g., developmental pathways, ecologies, behaviors, speciation) are historically constrained and their origins and functions are best explained via a phylogenetic approach. Herein, we provide a molecular evaluation of the phylogenetic placement of the model species Danio rerio within the genus Danio and among hypothesized closely related species and genera. Our analysis is derived from data using two nuclear genes (RAG1, rhodopsin) and five mitochondrial genes (ND4, ND4L, ND5, COI, cyt b) evaluated using parsimony, maximum likelihood, and Bayesian analyses. The family Cyprinidae is resolved as monophyletic but the subfamily Rasborinae (priority over Danioinae) is an unnatural assemblage. Danio is identified as a monophyletic group sister to a clade inclusive of the genera Chela, Microrasbora, Devario, and Inlecypris, not Devario nor Esomus as hypothesized in previous studies. Danio rerio is sister to D. kyathit among the species of Danio evaluated in this analysis. Microrasbora and Rasbora are non-monophyletic assemblages; however, Boraras is monophyletic.

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

Introduction: Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

Methods: We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

Results: We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

Conclusions: This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Phylogenetic Analyses of Eriotheca and Related Genera (Bombacoideae, Malvaceae)

Molecular and morphological data have shown that Bombacoideae and Malvoideae together form a well-supported Malvatheca clade. Phylogenetic relationships in Bombacoideae have been studied, but some genera in Bombax s. I. have not been adequately sampled for sufficiently variable molecular markers. The relationships of Eriotheca, for example, have yet to be resolved. Here, nuclear (ITS) and chloroplast (trnL-Fand matK) sequence data from 50 exemplars of Bombacoideae and seven additional taxa from other genera of Malvatheca were used to test monophyly of Eriotheca and its relationships with related genera of Bombax s. I. Parsimony and Bayesian analyses of individual and combined sequence data suggest that Eriotheca is not monophyletic as currently circumscribed but forms a paraphyletic grade containing Pachira s. 1. The newly discovered Eriotheca + Pachira clade has a probable synapomorphy of striate seeds. In addition, two other moderately supported clades emerged within the core Bombacoideae: Pseudobombax + Ceiba s. 1. and Bombax + Spirotheca + Pachira quinata. These three clades, and the African Rhodognaphalon together constitute the major clade of core Bombacoideae, whereas Adansonia appears to be more closely related to Catostemma, Scleronema, and Cavanillesia. The phylogenetic results imply three independent migrations from the New to Old World and homoplasy in staminal morphology.

Phylogenetic analysis of the order Pleuronectiformes (Teleostei) based on sequences of 12S and 16S mitochondrial genes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phylogenetic relationships within anuran clade Terrarana, with emphasis on the placement of Brazilian Atlantic rainforest frogs genus Ischnocnema (Anura: Brachycephalidae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phylogenetic Relationships in Eleocharis (Cyperaceae): C-4 Photosynthesis Origins and Patterns of Diversification in the Spikerushes

Eleocharis (Cyperaceae) is a morphologically and physiologically diverse lineage of 250 + species with a cosmopolitan distribution. We here explore phylogenetic relationships in this lineage using maximum parsimony and Bayesian inference analyses of nrDNA ITS and cpDNA trnC-ycf6 and ycf6-psbM sequence data with the goals of comparing our phylogenetic hypotheses to previous classifications, morphological variation, and photosynthetic pathway variation. Our results suggest that in Eleocharis C, photosynthesis has been derived at least three times, with several cases of possible reversion to C-3-like or intermediate pathways and several additional origins of C-3-C-4 intermediate photosynthetic pathways, as inferred by carbon isotope ratio measurements. Many classification units currently recognized in Eleocharis are not monophyletic, however, E. subgenus Limnochloa and E. subgenus Scirpidium are monophyletic. Other classification units largely corresponding to clades include E. subgenus Zinserlingia, E. subseries Chaetariae, and E. series Maculosae. Problems with species circumscription and morphological variation in several groups are discussed in light of the phylogeny, particularly in the context of species membership of seven focal clades found in the analyses.

Reinforcing and expanding the predictions of the disturbance vicariance hypothesis in Amazonian harlequin frogs: A molecular phylogenetic and climate envelope modelling approach

The disturbance vicariance hypothesis (DV) has been proposed to explain speciation in Amazonia, especially its edge regions, e. g. in eastern Guiana Shield harlequin frogs (Atelopus) which are suggested to have derived from a cool-adapted Andean ancestor. In concordance with DV predictions we studied that (i) these amphibians display a natural distribution gap in central Amazonia; (ii) east of this gap they constitute a monophyletic lineage which is nested within a pre-Andean/western clade; (iii) climate envelopes of Atelopus west and east of the distribution gap show some macroclimatic divergence due to a regional climate envelope shift; (iv) geographic distributions of climate envelopes of western and eastern Atelopus range into central Amazonia but with limited spatial overlap. We tested if presence and apparent absence data points of Atelopus were homogenously distributed with Ripley's K function. A molecular phylogeny (mitochondrial 16S rRNA gene) was reconstructed using Maximum Likelihood and Bayesian Inference to study if Guianan Atelopus constitute a clade nested within a larger genus phylogeny. We focused on climate envelope divergence and geographic distribution by computing climatic envelope models with MaxEnt based on macroscale bioclimatic parameters and testing them by using Schoener's index and modified Hellinger distance. We corroborated existing DV predictions and, for the first time, formulated new DV predictions aiming on species' climate envelope change. Our results suggest that cool-adapted Andean Atelopus ancestors had dispersed into the Amazon basin and further onto the eastern Guiana Shield where, under warm conditions, they were forced to change climate envelopes. © 2010 The Author(s).

Phylogenetic study of Class Armophorea (Alveolata, Ciliophora) based on 18S-rDNA data

The 18S rDNA phylogeny of Class Armophorea, a group of anaerobic ciliates, is proposed based on an analysis of 44 sequences (out of 195) retrieved from the NCBI/GenBank database. Emphasis was placed on the use of two nucleotide alignment criteria that involved variation in the gap-opening and gap-extension parameters and the use of rRNA secondary structure to orientate multiple-alignment. A sensitivity analysis of 76 data sets was run to assess the effect of variations in indel parameters on tree topologies. Bayesian inference, maximum likelihood and maximum parsimony phylogenetic analyses were used to explore how different analytic frameworks influenced the resulting hypotheses. A sensitivity analysis revealed that the relationships among higher taxa of the Intramacronucleata were dependent upon how indels were determined during multiple-alignment of nucleotides. The phylogenetic analyses rejected the monophyly of the Armophorea most of the time and consistently indicated that the Metopidae and Nyctotheridae were related to the Litostomatea. There was no consensus on the placement of the Caenomorphidae, which could be a sister group of the Metopidae + Nyctorheridae, or could have diverged at the base of the Spirotrichea branch or the Intramacronucleata tree.