Assessing rapid evolution in a changing environment.

Climate change poses a serious threat to species persistence. Effective modelling of evolutionary responses to rapid climate change is therefore essential. In this review we examine recent advances in phylogenetic comparative methods, techniques normally used to study adaptation over long periods, which allow them to be applied to the study of adaptation over shorter time scales. This increased applicability is largely due to the emergence of more flexible models of character evolution and the parallel development of molecular technologies that can be used to assess adaptive variation at loci scattered across the genome. The merging of phylogenetic and population genetic approaches to the study of adaptation has significant potential to advance our understanding of rapid responses to environmental change.

Genome Biology and Evolution of the Animal Kingdom

Estudi realitzat a partir d’una estada a la Institut J.W. Jenkinson Laboratory for Evolution and Development of the University of Oxford, Regne Unit, entre 2010 i 2012. He estat membre del laboratori del Professor Peter W.H. Holland com a becari post-doctoral Beatriu de Pinós des de setembre de 2010 al setembre de 2012. El nostre projecte de recerca se centra en l'anàlisi genòmic comparatiu del Regne Animal, tot explorant el contingut dels genomes a través de totes les branques de l'arbre dels animals. Totes les referències a les meves publicacions durant aquest post-doc es poden trobar a http://about.me/jordi_paps. Crec que el nombre i la qualitat dels resultats del meu post-doc, un total de 8 publicacions incloent dos articles a la prestigiosa revista Nature, són prova de l'èxit d'aquest post-doc. Prof Peter W. H. Holland (Departament de Zoologia de la Universitat d'Oxford) i jo som coautors de tres articles de genòmica comparativa, resultats directes d'aquest projecte: 1) comparació de families gèniques entre vertebrats invertebrats (Briefings in Functional Genomics), 2) el genoma de l'ostra (publicat a la revista Nature), i 3) els genomes de 6 platihelmints paràsits (acceptat també a Nature). A més, tenim altres 2 treballs en preparació. Un d'ells analitza l'evolució, expressió i funció dels gens Hox al a la tènia Hymenolepis. El perfil fi d'aquests gens clau del desenvolupament esclareix els canvis d'estil de vida dels organismes. A més, durant aquest últim post-doc he participat en diverses col•laboracions, incloent anàlisi de gens d'envelliment a cucs plans, un estudi sobre la filogènia del grup Gastrotricha, una revisió de l'evolució phylum Platyhelminthes, així com un capítol d'un llibre sobre l'evolució dels animals bilaterals. Finalment, gràcies a la beca Beatriu de Pinós, el Prof. Peter W.H. Holland m'ha convidat a formar part del seu equip com un investigador post-doctoral en el seu projecte ERC Advance actual sobre duplicacions genòmiques.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are “genomic fossils” valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome’s structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction (∼80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

SESAM - a new framework integrating macroecological and species distribution models for predicting spatio-temporal patterns of species assemblages

Two different approaches currently prevail for predicting spatial patterns of species assemblages. The first approach (macroecological modelling, MEM) focuses directly on realised properties of species assemblages, whereas the second approach (stacked species distribution modelling, S-SDM) starts with constituent species to approximate assemblage properties. Here, we propose to unify the two approaches in a single 'spatially-explicit species assemblage modelling' (SESAM) framework. This framework uses relevant species source pool designations, macroecological factors, and ecological assembly rules to constrain predictions of the richness and composition of species assemblages obtained by stacking predictions of individual species distributions. We believe that such a framework could prove useful in many theoretical and applied disciplines of ecology and evolution, both for improving our basic understanding of species assembly across spatio-temporal scales and for anticipating expected consequences of local, regional or global environmental changes. In this paper, we propose such a framework and call for further developments and testing across a broad range of community types in a variety of environments.

Plant and arbuscular mycorrhizal fungal diversity - are we looking at the relevant levels of diversity and are we using the right techniques?

Fungal symbionts commonly occur in plants influencing host growth, physiology, and ecology (Carlile et al., 2001). However, while whole-plant growth responses to biotrophic fungi are readily demonstrated, it has been much more difficult to identify and detect the physiological mechanisms responsible. Previous work on the clonal grass Glyceria striata has revealed that the systemic fungal endophyte Epichloë glyceriae has a positive effect on clonal growth of its host (Pan & Clay, 2002; 2003). The latest study from these authors, in this issue (pp. 467- 475), now suggests that increased carbon movement in hosts infected by E. glyceriae may function as one mechanism by which endophytic fungi could increase plant growth. Given the widespread distribution of both clonal plants and symbiotic fungi, this research will have implications for our understanding of the ecology and evolution of fungus-plant associations in natural communities.

The role of gene duplication in the origin and evolution of new biological functions

Abstract : Gene duplication is an essential source of material for the origin of genetic novelty and the evolution of lineage- or species-specific phenotypic traits. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with a significant number of gene copies during the last ~63 million years (MYA) of primate evolution. We estimated that at least 1 new functional gene (retrogene) per MYA emerged by retroposition in the primate lineage leading to humans. Using a combination of comparative sequencing and evolutionary simulations, we obtained strong evidence of functionality for 7 primate specific retrogenes. Most of these genes are specifically expressed in testis suggesting that retroposition has contributed with genetic raw material necessary for the evolution ofmale-specific functions in primates. We characterized CDC14Bretro (identified in the previous survey) that originated from the retroposition of a cell cycle gene - CDC14B - in the common ancestor of humans and apes. We demonstrate that CDC14Bretro experienced a period of intense positive selection in the African ape ancestor. By virtue of the amino acid substitutions that occurred during this period CDC 14Bretro adapted to a new subcellular compartment in African apes. Further analyses indicate that this subcellular shift reflects the evolution of anew functional role of CDC 14Bretro. Prompted by this result, we used yeast (Saccharomyces cerevisiae) to investigate on a global scale the extent of functional diversification of duplicate genes through the subcellular adaptation of their encoded proteins. We found that duplicate proteins frequently evolved new cellular localization patterns, either by partitioning of ancestral localizations ("sublocalization"), or more frequently by relocalization to previously unoccupied compartments ("neolocalization"). Interestingly, proteins involved in processes with a wider subcellular distribution more frequently evolved new localization patterns suggesting that subcellular localization changes are dependent on progenitor gene functions. Relocated proteins adapted to their new subcellular environments and evolved new functional roles through changes of their physio-chemical properties, expression levels, and interaction partners. Our work suggests an important role of subcellular adaptation for the emergence of new gene functions.

Hybridization and Population Structure of the Culex Pipiens Complex in the Islands of Macaronesia

The Culex pipiens complex includes two widespread mosquito vector species, Cx. pipiens and Cx. quinquefasciatus. The distribution of these species varies in latitude, with the former being present in temperate regions and the latter in tropical and subtropical regions. However, their distribution range overlaps in certain areas and interspecific hybridization has been documented. Genetic introgression between these species may have epidemiological repercussions for West Nile virus (WNV) transmission. Bayesian clustering analysis based on multilocus genotypes of 12 microsatellites was used to determine levels of hybridization between these two species in Macaronesian islands, the only contact zone described in West Africa. The distribution of the two species reflects both the islands’ biogeography and historical aspects of human colonization. Madeira Island displayed a homogenous population of Cx. pipiens, whereas Cape Verde showed a more intriguing scenario with extensive hybridization. In the islands of Brava and Santiago, only Cx. quinquefasciatus was found, while in Fogo and Maio high hybrid rates (~40%) between the two species were detected. Within the admixed populations, second-generation hybrids (~50%) were identified suggesting a lack of isolation mechanisms. The observed levels of hybridization may locally potentiate the transmission to humans of zoonotic arboviruses such as WNV.

PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth-death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open-source command-line Python program available at http://sourceforge.net/projects/pyrate/.

Molecular epidemiology and evolution of HIV-1 in Portugal and portuguese speaking african countries

The aims of this thesis were to better characterize HIV-1 diversity in Portugal, Angola, Mozambique and Cape Verde and to investigate the origin and epidemiological history of HIV-1 in these countries. The impact of these issues in diagnosis, disease progression and susceptibility to ARV therapy was also investigated. Finally, the nature, dynamics and prevalence of transmitted drug resistance (TDR) was determined in untreated HIV-1 infected patients. In Angola, practically all HIV-1 genetic forms were found, including almost all subtypes, untypable (U) strains, CRFs and URFs. Recombinants (first and second generation) were present in 47.1% of the patients. HIV/AIDS epidemic in Angola probably started in 1961, the major cause being the independence war, subsequently spreading to Portugal. In Maputo, 81% of the patients were infected with subtype C viruses. Subtype G, U and recombinants such as CRF37_cpx, were also present. The results suggest that HIV-1 epidemic in Mozambique is evolving rapidly in genetic complexity. In Cape Verde, where HIV-1 and HIV-2 co-circulate, subtype G is the prevailed subtype. Subtypes B, C, F1, U, CRF02_AG and other recombinant strains were also found. HIV-2 isolates belonged to group A, some being closely related to the original ROD isolate. In all three countries numerous new polymorphisms were identified in the RT and PR of HIV-1 viruses. Mutations conferring resistance to the NRTIs or NNRTIs were found in isolates from 2 (2%) patients from Angola, 4 (6%) from Mozambique and 3 (12%) from Cape Verde. None of the isolates containing TDR mutations would be fully sensitive to the standard first-line therapeutic regimens used in these countries. Close surveillance in treated and untreated populations will be crucial to prevent further transmission of drug resistant strains and maximize the efficacy of ARV therapy. In Portugal, investigation of a seronegative case infection with rapid progression to AIDS and death revealed that the patient was infected with a CRF14_BG-like R5-tropic strain selectively transmitted by his seropositive sexual partner. The results suggest a massive infection with a highly aggressive CRF14_BG like strain and/or the presence of an unidentified immunological problem that prevented the formation of HIV-1-specific antibodies. Near full-length genomic sequences obtained from three unrelated patients enabled the first molecular and phylogenomic characterization of CRF14_BG from Portugal; all sequences were strongly related with CRF14_BG Spanish isolates. The mean date of origin of CRF14_BG was estimated to be 1992. We propose that CRF14_BG emerged in Portugal in the early 1990s, spread to Spain in late 1990s as a consequence of IDUs migration and then to the rest of Europe. Most CRF14_BG strains were predicted to use CXCR4 and were associated with rapid CD4 depletion and disease progression. Finally, we provide evidence suggesting that the X4 tropism of CRF14_BG may have resulted from convergent evolution of the V3 loop possibly driven by an effective escape from neutralizing antibody response.

The role of phylogeny in the spatial distributions and assembly of mountain plant communities.

A major challenge in community ecology is a thorough understanding of the processes that govern the assembly and composition of communities in time and space. The growing threat of climate change to the vascular plant biodiversity of fragile ecosystems such as mountains has made it equally imperative to develop comprehensive methodologies to provide insights into how communities are assembled. In this perspective, the primary objective of this PhD thesis is to contribute to the theoretical and methodological development of community ecology, by proposing new solutions to better detect the ecological and evolutionary processes that govern community assembly. As phylogenetic trees provide by far, the most advanced tools to integrate the spatial, ecological and evolutionary dynamics of plant communities, they represent the cornerstone on which this work was based. In this thesis, I proposed new solutions to: (i) reveal trends in community assembly on phylogenies, depicted by the transition of signals at the nodes of the different species and lineages responsible for community assembly, (ii) contribute to evidence the importance of evolutionarily labile traits in the distribution of mountain plant species. More precisely, I demonstrated that phylogenetic and functional compositional turnover in plant communities was driven by climate and human land use gradients mostly influenced by evolutionarily labile traits, (iii) predict and spatially project the phylogenetic structure of communities using species distribution models, to identify the potential distribution of phylogenetic diversity, as well as areas of high evolutionary potential along elevation. The altitudinal setting of the Diablerets mountains (Switzerland) provided an appropriate model for this study. The elevation gradient served as a compression of large latitudinal variations similar to a collection of islands within a single area, and allowed investigations on a large number of plant communities. Overall, this thesis highlights that stochastic and deterministic environmental filtering processes mainly influence the phylogenetic structure of plant communities in mountainous areas. Negative density-dependent processes implied through patterns of phylogenetic overdispersion were only detected at the local scale, whereas environmental filtering implied through phylogenetic clustering was observed at both the regional and local scale. Finally, the integration of indices of phylogenetic community ecology with species distribution models revealed the prospects of providing novel and insightful explanations on the potential distribution of phylogenetic biodiversity in high mountain areas. These results generally demonstrate the usefulness of phylogenies in inferring assembly processes, and are worth considering in the theoretical and methodological development of tools to better understand phylogenetic community structure.

Intraspecific genetic variation in arbuscular mycorrhizal fungi and its consequences for molecular biology, ecology, and development of inoculum

It has been known for some time that different arbuscular mycorrhizal fungal (AMF) taxa confer differences in plant growth. Although genetic variation within AMF species has been given less attention, it could potentially be an ecologically important source of variation. Ongoing studies on variability in AMF genes within Glomus intraradices indicate that at least for some genes, such as the BiP gene, sequence variability can be high, even in coding regions. This suggests that genetic variation within an AMF may not be selectively neutral. This clearly needs to be investigated in more detail for other coding regions of AMF genomes. Similarly, studies on AMF population genetics indicate high genetic variation in AMF populations, and a considerable amount of variation seen in phenotypes in the population can be attributed to genetic differences among the fungi. The existence of high within-species genetic variation could have important consequences for how investigations on AMF gene expression and function are conducted. Furthermore, studies of within-species genetic variability and how it affects variation in plant growth will help to identify at what level of precision ecological studies should be conducted to identify AMF in plant roots in the field. A population genetic approach to studying AMF genetic variability can also be useful for inoculum development. By knowing the amount of genetic variability in an AMF population, the maximum and minimum numbers of spores that will contain a given amount of genetic diversity can be estimated. This could be particularly useful for developing inoculum with high adaptability to different environments.

The fossil record and evolution of freshwater plants: A review

Palaeobotany applied to freshwater plants is an emerging field of palaeontology. Hydrophytic plants reveal evolutionary trends of their own, clearly distinct from those of the terrestrial and marine flora. During the Precambrian, two groups stand out in the fossil record of freshwater plants: the Cyanobacteria (stromatolites) in benthic environments and the prasinophytes (leiosphaeridian acritarchs) in transitional planktonic environments. During the Palaeozoic, green algae (Chlorococcales, Zygnematales, charophytes and some extinct groups) radiated and developed the widest range of morphostructural patterns known for these groups. Between the Permian and Early Cretaceous, charophytes dominated macrophytic associations, with the consequence that over tens of millions of years, freshwater flora bypassed the dominance of vascular plants on land. During the Early Cretaceous, global extension of the freshwater environments is associated with diversification of the flora, including new charophyte families and the appearance of aquatic angiosperms and ferns for the first time. Mesozoic planktonic assemblages retained their ancestral composition that was dominated by coenobial Chlorococcales, until the appearance of freshwater dinoflagellates in the Early Cretaceous. In the Late Cretaceous, freshwater angiosperms dominated almost all macrophytic communities worldwide. The Tertiary was characterised by the diversification of additional angiosperm and aquatic fern lineages, which resulted in the first differentiation of aquatic plant biogeoprovinces. Phytoplankton also diversified during the Eocene with the development of freshwater diatoms and chrysophytes. Diatoms, which were exclusively marine during tens of millions of years, were dominant over the Chlorococcales during Neogene and in later assemblages. During the Quaternary, aquatic plant communities suffered from the effects of eutrophication, paludification and acidification, which were the result of the combined impact of glaciation and anthropogenic disturbance.