Phylogenetic analyses of C4 photosynthesis evolution in grasses: genetics of convergence and convergence of genes

Summary : During the evolutionary diversification of organisms, similar ecological constraints led to the recurrent appearances of the same traits (phenotypes) in distant lineages, a phenomenon called convergence. In most cases, the genetic origins of the convergent traits remain unknown, but recent studies traced the convergent phenotypes to recurrent alterations of the same gene or, in a few cases, to identical genetic changes. However, these cases remain anecdotal and there is a need for a study system that evolved several times independently and whose genetic determinism is well resolved and straightforward, such as C4 photosynthesis. This adaptation to warm environments, possibly driven by past atmospheric CO2 decreases, consists in a CO2-concentrating pump, created by numerous morphological and biochemical novelties. All genes encoding C4 enzymes already existed in C3 ancestors, and are supposed to have been recruited through gene duplication followed by neo-functionalization, to acquire the cell specific expression pattern and altered kinetic properties that characterize Ca-specific enzymes. These predictions have so far been tested only in species-poor and ecologically marginal C4 dicots. The monocots, and especially the grass family (Poaceae), the most important C4 family in terms of species number, ecological dominance and economical importance, have been largely under-considered as suitable study systems. This thesis aimed at understanding the evolution of the C4 trait in grasses at a molecular level and to use the genetics of C4 photosynthesis to infer the evolutionary history of the C4 phenotype and its driving selective pressures. A molecular phylogeny of grasses and affiliated monocots identified 17 to 18 independent acquisitions of the C4 pathway in the grass family. A relaxed molecular clock was used to date these events and the first C4 evolution was estimated in the Chloridoideae subfamily, between 32-25 million years ago, at a period when atmospheric CO2 abruptly declined. Likelihood models showed that after the COZ decline the probability of evolving the C4 pathway strongly increased, confirming low CO2 as a likely driver of C4 photosynthesis evolution. In order to depict the genetic changes linked to the numerous C4 origins, genes encoding phopshoenolpyruvate carboxylase (PEPC), the key-enzyme responsible for the initial fixation of atmospheric CO2 in the C4 pathway, were isolated from a large sample of C3 and C4 grasses. Phylogenetic analyses were used to reconstruct the evolutionary history of the PEPC multigene family and showed that the evolution of C4-specific PEPC had been driven by positive selection on 21 codons simultaneously in up to eight C4 lineages. These selective pressures led to numerous convergent genetic changes in many different C4 clades, highlighting the repeatability of some evolutionary processes, even at the molecular level. PEPC C4-adaptive changes were traced and used to show multiple appearances of the C, pathway in clades where species tree inferences were unable to differentiate multiple C4 appearances and a single appearance followed by C4 to C3 reversion. Further investigations of genes involved in some of the C4 subtypes only (genes encoding decarboxylating enzymes NADP-malic enzyme and phosphoenolpyruvate carboxykinase) showed that these C4-enzymes also evolved through strong positive selection and underwent parallel genetic changes during the different Ca origins. The adaptive changes on these subtype-specific C4 genes were used to retrace the history of the C4-subtypes phenotypes, which revealed that the evolution of C4-PEPC and C4-decarboxylating enzymes was in several cases disconnected, emphasizing the multiplicity of the C4 trait and the gradual acquisition of the features that create the CO2-pump. Finally, phylogenetic analyses of a gene encoding the Rubisco (the enzyme responsible for the fixation of CO2 into organic compounds in all photosynthetic organisms) showed that C4 evolution switched the selective pressures on this gene. Five codons were recurrently mutated to adapt the enzyme kinetics to the high CO2 concentrations of C4 photosynthetic cells. This knowledge could be used to introgress C4-like Rubisco in C3 crops, which could lead to an increased yield under predicted future high CO2 atmosphere. Globally, the phylogenetic framework adopted during this thesis demonstrated the widespread occurrence of genetic convergence on C4-related enzymes. The genetic traces of C4 photosynthesis evolution allowed reconstructing events that happened during the last 30 million years and proved the usefulness of studying genes directly responsible for phenotype variations when inferring evolutionary history of a given trait. Résumé Durant la diversification évolutive des organismes, des pressions écologiques similaires ont amené à l'apparition récurrente de certains traits (phénotypes) dans des lignées distantes, un phénomène appelé évolution convergente. Dans la plupart des cas, l'origine génétique des traits convergents reste inconnue mais des études récentes ont montré qu'ils étaient dus dans certains cas à des changements répétés du même gène ou, dans de rares cas, à des changements génétiques identiques. Malgré tout, ces cas restent anecdotiques et il y a un réel besoin d'un système d'étude qui ait évolué indépendamment de nombreuses fois et dont le déterminisme génétique soit clairement identifié. La photosynthèse dite en Ça répond à ces critères. Cette adaptation aux environnements chauds, dont l'évolution a pu être encouragé par des baisses passées de la concentration atmosphérique en CO2, est constituée de nombreuses nouveautés morphologiques et biochimiques qui créent une pompe à CO2. La totalité des gènes codant les enzymes Ç4 étaient déjà présents dans les ancêtres C3. Leur recrutement pour la photosynthèse Ç4 est supposé s'être fait par le biais de duplications géniques suivies par une néo-fonctionnalisation pour leur conférer l'expression cellule-spécifique et les propriétés cinétiques qui caractérisent les enzymes C4. Ces prédictions n'ont jusqu'à présent été testées que dans des familles C4 contenant peu d'espèces et ayant un rôle écologique marginal. Les graminées (Poaceae), qui sont la famille C4 la plus importante, tant en termes de nombre d'espèces que de dominance écologique et d'importance économique, ont toujours été considérés comme un système d'étude peu adapté et ont fait le sujet de peu d'investigations évolutives. Le but de cette thèse était de comprendre l'évolution de la photosynthèse en C4 chez les graminées au niveau génétique et d'utiliser les gènes pour inférer l'évolution du phénotype C4 ainsi que les pressions de sélection responsables de son évolution. Une phylogénie moléculaire de la famille des graminées et des monocotylédones apparentés a identifié 17 à 18 acquisitions indépendantes de la photosynthèse chez les graminées. Grâce à une méthode d'horloge moléculaire relâchée, ces évènements ont été datés et la première apparition C4 a été estimée dans la sous-famille des Chloridoideae, il y a 32 à 25 millions d'années, à une période où les concentrations atmosphériques de CO2 ont décliné abruptement. Des modèles de maximum de vraisemblance ont montré qu'à la suite du déclin de CO2, la probabilité d'évoluer la photosynthèse C4 a fortement augmenté, confirmant ainsi qu'une faible concentration de CO2 est une cause potentielle de l'évolution de la photosynthèse C4. Afin d'identifier les mécanismes génétiques responsables des évolutions répétées de la photosynthèse C4, un segment des gènes codant pour la phosphoénolpyruvate carboxylase (PEPC), l'enzyme responsable de la fixation initiale du CO2 atmosphérique chez les plantes C4, ont été séquencés dans une centaine de graminées C3 et C4. Des analyses phylogénétiques ont permis de reconstituer l'histoire évolutive de la famille multigénique des PEPC et ont montré que l'évolution de PEPC spécifiques à la photosynthèse Ça a été causée par de la sélection positive agissant sur 21 codons, et ce simultanément dans huit lignées C4 différentes. Cette sélection positive a conduit à un grand nombre de changements génétiques convergents dans de nombreux clades différents, ce qui illustre la répétabilité de certains phénomènes évolutifs, et ce même au niveau génétique. Les changements sur la PEPC liés au C4 ont été utilisés pour confirmer des évolutions indépendantes du phénotype C4 dans des clades où l'arbre des espèces était incapable de différencier des apparitions indépendantes d'une seule apparition suivie par une réversion de C4 en C3. En considérant des gènes codant des protéines impliquées uniquement dans certains sous-types C4 (deux décarboxylases, l'enzyme malique à NADP et la phosphoénolpyruvate carboxykinase), des études ultérieures ont montré que ces enzymes C4 avaient elles-aussi évolué sous forte sélection positive et subi des changements génétiques parallèles lors des différentes origines de la photosynthèse C4. Les changements adaptatifs sur ces gènes liés seulement à certains sous-types C4 ont été utilisés pour retracer l'histoire des phénotypes de sous-types C4, ce qui a révélé que les caractères formant le trait C4 ont, dans certains cas, évolué de manière déconnectée. Ceci souligne la multiplicité du trait C4 et l'acquisition graduelle de composants participant à la pompe à CO2 qu'est la photosynthèse C4. Finalement, des analyses phylogénétiques des gènes codant pour la Rubisco (l'enzyme responsable de la fixation du CO2 en carbones organiques dans tous les organismes photosynthétiques) ont montré que l'évolution de la photosynthèse Ça a changé les pressions de sélection sur ce gène. Cinq codons ont été mutés de façon répétée afin d'adapter les propriétés cinétiques de la Rubisco aux fortes concentrations de CO2 présentes dans les cellules photosynthétiques des plantes C4. Globalement, l'approche phylogénétique adoptée durant cette thèse de doctorat a permis de démontré des phénomène fréquents de convergence génétique sur les enzymes liées à la photosynthèse C4. Les traces génétiques de l'évolution de la photosynthèse C4 ont permis de reconstituer des évènements qui se sont produits durant les derniers 30 millions d'années et ont prouvé l'utilité d'étudier des gènes directement responsables des variations phénotypiques pour inférer l'histoire évolutive d'un trait donné.

Evidence for adaptive radiation from a phylogenetic study of plant defenses.

One signature of adaptive radiation is a high level of trait change early during the diversification process and a plateau toward the end of the radiation. Although the study of the tempo of evolution has historically been the domain of paleontologists, recently developed phylogenetic tools allow for the rigorous examination of trait evolution in a tremendous diversity of organisms. Enemy-driven adaptive radiation was a key prediction of Ehrlich and Raven's coevolutionary hypothesis [Ehrlich PR, Raven PH (1964) Evolution 18:586-608], yet has remained largely untested. Here we examine patterns of trait evolution in 51 North American milkweed species (Asclepias), using maximum likelihood methods. We study 7 traits of the milkweeds, ranging from seed size and foliar physiological traits to defense traits (cardenolides, latex, and trichomes) previously shown to impact herbivores, including the monarch butterfly. We compare the fit of simple random-walk models of trait evolution to models that incorporate stabilizing selection (Ornstein-Ulenbeck process), as well as time-varying rates of trait evolution. Early bursts of trait evolution were implicated for 2 traits, while stabilizing selection was implicated for several others. We further modeled the relationship between trait change and species diversification while allowing rates of trait evolution to vary during the radiation. Species-rich lineages underwent a proportionately greater decline in latex and cardenolides relative to species-poor lineages, and the rate of trait change was most rapid early in the radiation. An interpretation of this result is that reduced investment in defensive traits accelerated diversification, and disproportionately so, early in the adaptive radiation of milkweeds.

Phylogenetic analysis of Tomato mosaic virus from Hemerocallis sp. and Impatiens hawkeri

The culture and commercialization of ornamental plants have considerably increased in the last years. To supply the commercial demand, several Hemerocallis and Impatiens varieties have been bred for appreciated qualities such as flowers with a diversity of shapes and colors. With the aim of characterizing the tobamovirus isolated from Hemerocallis sp. (tobamo-H) and Impatiens hawkeri (tobamo-I) from the USA and São Paulo, respectively, as well as to establish phylogenetic relationships between them and other Tobamovirus species, the viruses were submitted to RNA extraction, RT-PCR amplification, coat-protein gene sequencing and phylogenetic analyses. Comparison of tobamovirus homologous sequences yielded values superior to 98.5% of identity with Tomato mosaic virus (ToMV) isolates at the nucleotide level. In relation to tobamo-H, 100% of identity with ToMV from tomatoes from Australia and Peru was found. Based on maximum likelihood (ML) analysis it was suggested that tobamo-H and tobamo-I share a common ancestor with ToMV, Tobacco mosaic virus, Odontoglossum ringspot virus and Pepper mild mottle virus. The tree topology reconstructed under ML methodology shows a monophyletic group, supported by 100% of bootstrap, consisting of various ToMV isolates from different hosts, including some ornamentals, from different geographical locations. The results indicate that Hemerocallis sp. and I. hawkeri are infected by ToMV. This is the first report of the occurrence of this virus in ornamental species in Brazil.

Coevolution of Heliconius spp. and Passiflora spp. : a phylogenetic comparison

Although a substantial amount of research has been done on all aspects ofHeliconius biology and their ecological interactions with Passiflora, there has not hitherto been a phylogenetic examination of this association for coevolution. To test the HeliconiuslPassilfora association for coevolutionary congruence, phylogenies for each group were established and compared. The phylogeny for 14 species ofHeliconiinae from Costa Rica was based on combined sequence data from rRNA ITS 2 and partial EF-1a gene regions. For the Passifloraceae, 17 host plant species were utilized to establish a phylogeny based on tRNALeucine and ITS 1/5.8S1 ITS 2 sequence data. The phylogenies for both groups were largely in agreement with current classification (for Passifloraceae) and previously established phylogenies. Associations with the large subgenera Passiflora and Decaloba correspond with the two major Advanced Radiation groups in Heliconius. Although strict congruence above subgenus level was not observed, broad scale congruence was evident. One main host shift as well as other possible explanations for lack of strict congruence are suggested.

Phylogenetic variation in the shoot mineral concentration of angiosperms

The calcium (Ca) concentration of plant shoot tissues varies systematically between angiosperm orders. The phylogenetic variation in the shoot concentration of other mineral nutrients has not yet been described at an ordinal level. The aims of this study were (1) to quantify the shoot mineral concentration of different angiosperm orders, (2) to partition the phylogenetic variation in shoot mineral concentration between and within orders, (3) to determine if the shoot concentration of different minerals are correlated across angiosperm species, and (4) to compare experimental data with published ecological survey data on 81 species sampled from their natural habitats. Species, selected pro rata from different angiosperm orders, were grown in a hydroponic system under a constant external nutrient regime. Shoots of 117 species were sampled during vegetative growth. Significant variation in shoot carbon (C), calcium (Ca), potassium (K), and magnesium (Mg) concentration occurred between angiosperm orders. There was no evidence for systematic differences in shoot phosphorus (P) or organic-nitrogen (N) concentration between orders. At a species level, there were strong positive correlations between shoot Ca and Mg concentration, between shoot P and organic-N concentration, and between shoot K concentration and shoot fresh weight:dry weight ratio. Shoot C and cation concentration correlated negatively at a species level. Species within the Poales and the Caryophyllales had distinct shoot mineralogies in both the designed experiment and in the ecological survey.

Phylogenetic relationships in Kumanoa (Batrachospermales, Rhodophyta) species in Brazil with the proposal of Kumanoa amazonensis sp nov.

This Study evaluated the species-level taxonomy and phylogenetic relationship among Kumanoa species from Brazil with other regions of the world based on the plastid-encoded RUBISCO large Subunit gene (rhcL). Partial rbcL sequences were obtained for 11 Kulnanoa specimens. Eight species are recognised from Brazil on the basis of molecular and morphological data: seven previously described (K abilii, K ambignia, K. breviarticulata, K. cipoensis, K. equisetoidea, K. globospora and K procarpa) and a new species here proposed (K. amazonensis sp. nov. Necchi & Vis). The new species has reduced and dense whorls but differs from the two closest related species in lacking secondary fascicles. Previously proposed infrageneric categories were not supported by the molecular data. Species described and endemic (K. breviarticulata, K. cipoensis, K equiseloidea and K. procarpa) to Brazil are not grouped together but are variously related to other species from North America, Europe and Australasia. With the species recognised in this study using molecular and morphological data and those previously distinguished by morphology, 13 species of Kumanoa are Currently documented from Brazil.

Biogeography of Glandulocaudinae (Teleostei : Characiformes : Characidae) revisited: phylogenetic patterns, historical geology and genetic connectivity

The biogeography of the Glandulocaudinae ( former Glandulocaudini) is reviewed. The major pattern of diversification presented by this group of freshwater fishes can be clearly associated to the main aspects of the tectonic evolution of the southern portion of the Cis-Andean South American Platform. The phylogenetic relationships within the group suggest that the clade represented by Lophiobrycon is the sister-group of the more derived clade represented by the genus Glandulocauda and Mimagoniates. Lophiobrycon and Glandulocauda occur in areas of the ancient crystalline shield of southeastern Brazil and their present allopatric distribution is probably due to relict survival and tectonic vicariant events. Populations of Glandulocauda melanogenys are found in contiguous drainages in presently isolated upper parts of the Tiete, Guaratuba, Itatinga, and Ribeira de Iguape basins and this pattern of distribution is probably the result of river capture caused by tectonic processes that affected a large area in eastern and southeastern Brazil. The species of Mimagoniates are predominantly distributed along the eastern and southeastern coastal areas, but M. microlepis is additionally found in the rio Iguacu and Tibagi basins. Mimagoniates barberi occurs in both SW margin of the upper rio Parana basin and the lower Paraguay and Mimagoniates sp. occurs in the upper Paraguay river basin. Tectonic activations of the Continental Rift of Southeastern Brazil along the eastern margin of the Upper Parana basin promoted population fragmentation responsible of the present day distribution presented by Glandulocauda melanogenys. We hypothesize that occurrence of Mimagoniates along the lowland area around the Parana basin was due to a single or a multiple fragmentation of populations along the W-SW border of the upper Parana Basin, probably due to the major tectonic origin of the Chaco-Pantanal wetland foreland basins since the Miocene as well as Cenozoic tectonic activity along the borders of the upper Parana basin, such as in the eastern Paraguay, in the Asuncion Rift. Distributional pattern of Mimagoniates suggests that its initial diversification may be related to the tectonic evolution of the Chaco-Pantanal foreland basin system and a minimum age of 2.5 M.Y are proposed for this monophyletic group. Previous hypotheses on sea level fluctuations of the late Quaternary as being the main causal mechanism promoting cladogenesis and speciation of the group are critically reviewed. Phylogeographic studies based on molecular data indicate significant differences among the isolated populations of M. microlepis. These findings suggest that a much longer period of time and a paleogeographic landscape configuration of the Brazilian southeastern coastal region explain the present observed phylogenetic and biogeographic patterns.

Phylogenetic analysis of PRP8 intein in Paracoccidioides brasiliensis species

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

Humidity levels drive reproductive modes and phylogenetic diversity of amphibians in the Brazilian Atlantic Forest

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

A phylogenetic analysis of Pleurodema (Anura: Leptodactylidae: Leiuperinae) based on mitochondrial and nuclear gene sequences, with comments on the evolution of anuran foam nests

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

Phylogenetic relationships in Kumanoa (Batrachospermales, Rhodophyta) species in Brazil with the proposal of Kumanoa amazonensis sp nov.

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Phylogenetic relationship of Sirodotia species (Batrachospermales, Rhodophyta) in North and South America

Sequence data from the RUBISCO large subunit (rbcL) plastid gene and nuclear small-subunit ribosomal DNA (SSU rDNA) were examined for five samples of Sirodotia delicatula from southeastern Brazil. Data from six North American samples previously identified as S. huillensis and S. suecica were also included in the analysis. Molecular data supported the continued recognition of these three species as separate entities, although one of the North American collections was misidentified. These results were shown to be congruent with morphology, chromosome number and geographic distribution. S. delicatula is more closely related to S. huillensis, both occurring in tropical-subtropical regions, than either to S. suecica with a temperate-boreal distribution. There was little rbcL variation within S. delicatula from Brazil and Costa Rica (the latter a collection previously identified as S. huillensis), with the six samples sequenced diverging from each other by 0-8 bp (0-0.67%). SSU rDNA data set did not provide sufficient resolution to infer phylogenetic relationships among the species of this group due to the low rates of variation (5 bp). Sirodotia was a well-supported clade (100% bootstrap or 1.00 a posteriori probability) based on rbcL sequences. Thus, the results confirm that Sirodotia is a monophyletic group within the Batrachospermales and we continue to recognize it at the generic level. The species S. delicatula, S. huillensis and S. suecica are morphologically and genetically distinct.

Phylogenetic systematics of dart-poison frogs and their relatives (Amphibia: Athesphatanura: Dendrobatidae)

The known diversity of dart-poison frog species has grown from 70 in the 1960s to 247 at present, with no sign that the discovery of new species will wane in the foreseeable future. Although this growth in knowledge of the diversity of this group has been accompanied by detailed investigations of many aspects of the biology of dendrobatids, their phylogenetic relationships remain poorly understood. This study was designed to test hypotheses of dendrobatid diversification by combining new and prior genotypic and phenotypic evidence in a total evidence analysis. DNA sequences were sampled for five mitochondrial and six nuclear loci (approximately 6,100 base pairs [bp]; x=3,740 bp per terminal; total dataset composed of approximately 1.55 million bp), and 174 phenotypic characters were scored from adult and larval morphology, alkaloid profiles, and behavior. These data were combined with relevant published DNA sequences. Ingroup sampling targeted several previously unsampled species, including Aromobates nocturnus, which was hypothesized previously to be the sister of all other dendrobatids. Undescribed and problematic species were sampled from multiple localities when possible. The final dataset consisted of 414 terminals: 367 ingroup terminals of 156 species and 47 outgroup terminals of 46 species. Direct optimization parsimony analysis of the equally weighted evidence resulted in 25,872 optimal trees. Forty nodes collapse in the strict consensus, with all conflict restricted to conspecific terminals. Dendrobatids were recovered as monophyletic, and their sister group consisted of Crossodactylus, Hylodes, and Megaelosia, recognized herein as Hylodidae. Among outgroup taxa, Centrolenidae was found to be the sister group of all athesphatanurans except Hylidae, Leptodactyidae was polyphyletic, Thoropa was nested within Cycloramphidae, and Ceratophryinae was paraphyletic with respect to Telmatobiinae. Among dendrobatids, the monophyly and content of Mannophryne and Phyllobates were corroborated. Aromobates nocturnus and Colostethus saltuensis were found to be nested within Nephelobates, and Minyobates was paraphyletic and nested within Dendrobates. Colostethus was shown to be rampantly nonmonophyletic, with most species falling into two unrelated cis- and trans-Andean clades. A morphologically and behaviorally diverse clade of median lingual process-possessing species was discovered. In light of these findings and the growth in knowledge of the diversity of this large clade over the past 40 years, we propose a new, monophyletic taxonomy for dendrobatids, recognizing the inclusive clade as a superfamily (Dendrobatoidea) composed of two families (one of which is new), six subfamilies (three new), and 16 genera (four new). Although poisonous frogs did not form a monophyletic group, the three poisonous lineages are all confined to the revised family Dendrobatidae, in keeping with the traditional application of this name. We also propose changes to achieve a monophyletic higher-level taxonomy for the athesphatanuran outgroup taxa. Analysis of character evolution revealed multiple origins of phytotelm-breeding, parental provisioning of nutritive oocytes for larval consumption (larval oophagy), and endotrophy. Available evidence indicates that transport of tadpoles on the dorsum of parent nurse frogs-a dendrobatid synapomorphy-is carried out primitively by male nurse frogs, with three independent origins of female transport and five independent origins of biparental transport. Reproductive amplexus is optimally explained as having been lost in the most recent common ancestor of Dendrobatoidea, with cephalic amplexus arising independently three times. © American Museum of Natural History 2006.

Phylogenetic relationships of the genus Sibynophis (Serpentes: Colubroidea)

We present the results of the first molecular analysis of the phylogenetic affinities of the Asian colubroid genus Sibynophis. We recovered a sister-group relationship between Sibynophis and the New World Scaphiodontophis. Although Liophidium sometimes is associated with these genera, the relationship is distant. Morphological characters that Liophidium shares with Sibynophis and Scaphiodontophis are resolved as homoplasies that probably reflect the similarities of their specialized feeding habits. The traditional subfamily Sibynophiinae is elevated to the family-level, and Scaphiodontophiinae is placed in its synonymy.