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.

Description and Phylogenetic Relationships of a New Genus and Species of Lizard (Squamata, Gymnophthalmidae) from the Amazonian Rainforest of Northern Brazil

Marinussaurus curupira, a new genus and species of Gymnophthalmidae lizard is described from Iranduba, state of Amazonas, Brazil. The genus is characterized by an elongate body; short and stout pentadactyl limbs; all digits clawed; single frontonasal; two prefrontals; absence of frontoparietals; interparietal and parietals forming a straight posterior margin, with interparietal shorter than parietals; distinctive ear opening and eyelid; few temporals; three pairs of chin shields; nasal divided; a distinct collar; smooth, mainly hexagonal, dorsal scales; smooth quadrangular ventral scales; two precloacal and three femoral pores on each side in males; pores between three or four scales. Parsimony (PAR) and partitioned Bayesian (BA) phylogenetic analyses with morphological and molecular data recovered the new genus as a member of the Ecpleopodini radiation of the Cercosaurinae. A close relationship of the new genus with Arthrosaura is postulated.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

Background: Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei. Results: The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infraordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes. Conclusions: Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Inference of phylogeny and taxonomy within the Didymozoidae (Digenea) from the second internal transcribed spacer (ITS2) of ribosomal DNA

DNA sequences of the second internal transcribed spacer (ITS2) of ribosomal DNA (rDNA) were determined for 11 species from four genera of Didymozoinae (Indodidymozoon, Helicodidymozoon, Rhopalotrema and Neometadidymozoon) and a species of the Lecithasteridae, Lecithaster stellatus. Sequences were used to test the validity of species recognised on morphological criteria and to infer phylogenetic relationships. Sequences of the 11 didymozoids differed by 0.5% to 19%. Our phylogenetic analyses: (i) indicate that species in the genera Helicodidymozoon and Rhopalotrema are a monophyletic group; (ii) support separation of the genus Helicodidymozoon from the genera Indodidymozoon and Neometadidymozoon; and (iii) support recognition of Rhopalotrema as a genus distinct from Neometadidymozoon. We found the gonochoristic species, I. pearsoni and I. suttiei, to be genetically similar to the hermaphroditic species in the genus Indodidymozoon and found no evidence to indicate that they belong in a separate genus.

Molecular characterization of the Hepatitis B virus genotypes in Colombia: A Bayesian inference on the genotype F

Hepatitis B is a worldwide health problem affecting about 2 billion people and more than 350 million are chronic carriers of the virus. Nine HBV genotypes (A to I) have been described. The geographical distribution of HBV genotypes is not completely understood due to the limited number of samples from some parts of the world. One such example is Colombia, in which few studies have described the HBV genotypes. In this study, we characterized HBV genotypes in 143 HBsAg-positive volunteer blood donors from Colombia. A fragment of 1306 bp partially comprising HBsAg and the DNA polymerase coding regions (S/POL) was amplified and sequenced. Bayesian phylogenetic analyses were conducted using the Markov Chain Monte Carlo (MCMC) approach to obtain the maximum clade credibility (MCC) tree using BEAST v.1.5.3. Of all samples, 68 were positive and 52 were successfully sequenced. Genotype F was the most prevalent in this population (77%) - subgenotypes F3 (75%) and Fib (2%). Genotype G (7.7%) and subgenotype A2 (15.3%) were also found. Genotype G sequence analysis suggests distinct introductions of this genotype in the country. Furthermore, we estimated the time of the most recent common ancestor (TMRCA) for each HBV/F subgenotype and also for Colombian F3 sequences using two different datasets: (i) 77 sequences comprising 1306 bp of S/POL region and (ii) 283 sequences comprising 681 bp of S/POL region. We also used two other previously estimated evolutionary rates: (i) 2.60 x 10(-4) s/s/y and (ii) 1.5 x 10(-5) s/s/y. Here we report the HBV genotypes circulating in Colombia and estimated the TMRCA for the four different subgenotypes of genotype F. (C) 2010 Elsevier B.V. All rights reserved.

Human synthetic lethal inference as potential anti-cancer target gene detection

Background: Two genes are called synthetic lethal (SL) if mutation of either alone is not lethal, but mutation of both leads to death or a significant decrease in organism's fitness. The detection of SL gene pairs constitutes a promising alternative for anti-cancer therapy. As cancer cells exhibit a large number of mutations, the identification of these mutated genes' SL partners may provide specific anti-cancer drug candidates, with minor perturbations to the healthy cells. Since existent SL data is mainly restricted to yeast screenings, the road towards human SL candidates is limited to inference methods. Results: In the present work, we use phylogenetic analysis and database manipulation (BioGRID for interactions, Ensembl and NCBI for homology, Gene Ontology for GO attributes) in order to reconstruct the phylogenetically-inferred SL gene network for human. In addition, available data on cancer mutated genes (COSMIC and Cancer Gene Census databases) as well as on existent approved drugs (DrugBank database) supports our selection of cancer-therapy candidates.Conclusions: Our work provides a complementary alternative to the current methods for drug discovering and gene target identification in anti-cancer research. Novel SL screening analysis and the use of highly curated databases would contribute to improve the results of this methodology.

An evaluation of new parsimony-based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae

Aim  Recently developed parametric methods in historical biogeography allow researchers to integrate temporal and palaeogeographical information into the reconstruction of biogeographical scenarios, thus overcoming a known bias of parsimony-based approaches. Here, we compare a parametric method, dispersal-extinction-cladogenesis (DEC), against a parsimony-based method, dispersal-vicariance analysis (DIVA), which does not incorporate branch lengths but accounts for phylogenetic uncertainty through a Bayesian empirical approach (Bayes-DIVA). We analyse the benefits and limitations of each method using the cosmopolitan plant family Sapindaceae as a case study.Location  World-wide.Methods  Phylogenetic relationships were estimated by Bayesian inference on a large dataset representing generic diversity within Sapindaceae. Lineage divergence times were estimated by penalized likelihood over a sample of trees from the posterior distribution of the phylogeny to account for dating uncertainty in biogeographical reconstructions. We compared biogeographical scenarios between Bayes-DIVA and two different DEC models: one with no geological constraints and another that employed a stratified palaeogeographical model in which dispersal rates were scaled according to area connectivity across four time slices, reflecting the changing continental configuration over the last 110 million years.Results  Despite differences in the underlying biogeographical model, Bayes-DIVA and DEC inferred similar biogeographical scenarios. The main differences were: (1) in the timing of dispersal events - which in Bayes-DIVA sometimes conflicts with palaeogeographical information, and (2) in the lower frequency of terminal dispersal events inferred by DEC. Uncertainty in divergence time estimations influenced both the inference of ancestral ranges and the decisiveness with which an area can be assigned to a node.Main conclusions  By considering lineage divergence times, the DEC method gives more accurate reconstructions that are in agreement with palaeogeographical evidence. In contrast, Bayes-DIVA showed the highest decisiveness in unequivocally reconstructing ancestral ranges, probably reflecting its ability to integrate phylogenetic uncertainty. Care should be taken in defining the palaeogeographical model in DEC because of the possibility of overestimating the frequency of extinction events, or of inferring ancestral ranges that are outside the extant species ranges, owing to dispersal constraints enforced by the model. The wide-spanning spatial and temporal model proposed here could prove useful for testing large-scale biogeographical patterns in plants.

Phylogenetic relationships and divergence times among dormice (Rodentia, Gliridae) based on three nuclear genes

We examined phylogenetic relationships among six species representing three subfamilies, Glirinae, Graphiurinae and Leithiinae with sequences from three nuclear protein-coding genes (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). Phylogenetic trees reconstructed from maximum-parsimony (MP), maximum-likelihood (ML) and Bayesian-inference (BI) analyses showed the monophyly of Glirinae (Glis and Glirulus) and Leithiinae (Dryomys, Eliomys and Muscardinus) with strong support, although the branch length maintaining this relationship was very short, implying rapid diversification among the three subfamilies. Divergence time estimates were calculated from ML (local clock model) and Bayesian-dating method using a calibration point of 25 Myr (million years) ago for the divergence between Glis and Glirulus, and 55 Myr ago for the split between lineages of Gliridae and Sciuridae on the basis of fossil records. The results showed that each lineage of Graphiuros, Glis, Glirulus and Muscardinus dates from the Late Oligocene to the Early Miocene period, which is mostly in agreement with fossil records. Taking into account that warm climate harbouring a glirid-favoured forest dominated from Europe to Asia during this period, it is considered that this warm environment triggered the prosperity of the glirid species through the rapid diversification. Glirulus japonicas is suggested to be a relict of this ancient diversification during the warm period.

Phylogenetic analysis of plastid and nuclear DNA sequences reveals a rapid late Miocene radiation of the temperate bamboo tribe Arundinarieae (Poaceae, Bambusoideae)

Background: Arundinarieae are a large tribe of temperate woody bamboos for which phylogenetics are poorly understood because of limited taxon sampling and lack of informative characters. Aims: This study assessed phylogenetic relationships, origins and classification of Arundinarieae. Methods: DNA sequences (plastid trnL-F; nuclear ITS) were used for parsimony and Bayesian inference including 41 woody bamboo taxa. Divergence dates were estimated using a relaxed Bayesian clock. Results: Arundinarieae were monophyletic but their molecular divergence was low compared to the tropical Bambuseae. Ancestors of the Arundinarieae lineage were estimated to have diverged from the other bamboos 23 (15-30) million years ago (Mya). However, the Arundinarieae radiation occurred 10 (6-16) Mya compared to 18 (11-25) Mya for the tropical Bambuseae. Some groups could be defined within Arundinarieae, but they do not correspond to recognised subtribes such as Arundinariinae or Shibataeinae. Conclusions: Arundinarieae are a relatively ancient bambusoid lineage that underwent a rapid radiation in the late Miocene. The radiation coincides with the continental collision of the Indo-Australian and Eurasian Plates. Arundinarieae are distributed primarily in East Asia and the Himalayas to northern Southeast Asia. It is unknown whether they were present in Asia long before their radiation, but we believe recent dispersal is a more likely scenario. Keywords: Arundinarieae; Bambuseae; internal transcribed spacer (ITS); molecular clock; phylogenetics; radiation; temperate bamboos; Thamnocalaminae; trnL-F

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)