Phylogeny and circumscription of Sapindaceae revisited: molecular sequence data, morphology and biogeography support recognition of a new family, Xanthoceraceae

Background and aims Recent studies have adopted a broad definition of Sapindaceae that includes taxa traditionally placed in Aceraceae and Hippocastanaceae, achieving monophyly but yielding a family difficult to characterize and for which no obvious morphological synapomorphy exists. This expanded circumscription was necessitated by the finding that the monotypic, temperate Asian genus Xanthoceras, historically placed in Sapindaceae tribe Harpullieae, is basal within the group. Here we seek to clarify the relationships of Xanthoceras based on phylogenetic analyses using a dataset encompassing nearly 3/4 of sapindaceous genera, comparing the results with information from morphology and biogeography, in particular with respect to the other taxa placed in Harpullieae. We then re-examine the appropriateness of maintaining the current broad, morphologically heterogeneous definition of Sapindaceae and explore the advantages of an alternative family circumscription. Methods Using 243 samples representing 104 of the 142 currently recognized genera of Sapindaceae s. lat. (including all in Harpullieae), sequence data were analyzed for nuclear (ITS) and plastid (matK, rpoB, trnD-trnT, trnK-matK, trnL-trnF and trnS-trnG) markers, adopting the methodology of a recent family-wide study, performing single-gene and total evidence analyses based on maximum likelihood (ML) and maximum parsimony (MP) criteria, and applying heuristic searches developed for large datasets, viz, a new strategy implemented in RAxML (for ML) and the parsimony ratchet (for MP). Bootstrap analyses were performed for each method to test for congruence between markers. Key results Our findings support earlier suggestions that Harpullieae are polyphyletic: Xanthoceras is confirmed as sister to all other sampled taxa of Sapindaceae s. lat.; the remaining members belong to three other clades within Sapindaceae s. lat., two of which correspond respectively to the groups traditionally treated as Aceraceae and Hippocastanaceae, together forming a clade sister to the largely tropical Sapindaceae s. str., which is monophyletic and morphologically coherent provided Xanthoceras is excluded. Conclusion To overcome the difficulties of a broadly circumscribed Sapindaceae, we resurrect the historically recognized temperate families Aceraceae and Hippocastanaceae, and describe a new family, Xanthoceraceae, thus adopting a monophyletic and easily characterized circumscription of Sapindaceae nearly identical to that used for over a century.

Conflict resolution and evolution of social structures in insect societies

In colonies of social Hymenoptera (which include all ants, as well as some wasp and bee species), only queens reproduce whereas workers generally perform other tasks. The evolution of worker's reproductive altruism can be explained by kin selection, which states that workers can indirectly transmit copies of their genes by helping the reproduction of relatives. The relatedness between queens and workers may however be low, particularly when there are multiple queens per colony, which limits the transmission of copies of workers genes and increases potential conflicts between colony members. In this thesis, we investigated the link between social structure variations and conflicts, and explored the mechanisms involved in variation of colony queen number in ants. According to kin selection, workers should rear the brood they are most related to. In social Hymenoptera, males are haploid whereas females (workers and queens) are diploid. As a result, workers can be up to three times more related to females than males in some colonies, where they should consequently favour the production of females. In contrast, queens are equally related to daughters and sons in all types of colonies and therefore should favour a balanced sex ratio. In a meta-analysis across all studies of social Hymenoptera, we showed that colony sex ratio is generally largely influenced by workers. Hence, the evolution of social structures where queens and workers are equally related to males and females may contribute to decrease the conflict between the two castes over colony sex ratio. Another conflict between queens and workers can occur over male production. Many species contain workers that still have the ability to lay haploid eggs. In some social structures, workers are on average more related to sons of queens than to sons of other workers. As a result, workers should eliminate worker-laid eggs to favour queen-laid eggs. We showed that in the ant Formica selysi, workers eliminate more worker-laid than queen-laid eggs, independently of colony social structure. These results therefore suggest that worker policing can evolve independently from relatedness, potentially because of costs of worker reproduction at the colony-level. Colony queen number is a key parameter that influences relatedness between group members. Queen body size is generally linked to the success of independent colony foundation by single queens and may influence the number of queens in the new colony. In the ant F. selysi, single-queen colonies produce larger queens than multiple-queen colonies. We showed that this association results from genes or maternal effects transmitted to the eggs. However, we also found that queens produced in colonies of the two social forms did not differ in their general ability to found new colonies independently. Queen body size may also influence queen dispersal ability and constrain small queens to be re-adopted in their original nest after mating at proximity. We tested the acceptance of new queens in another ant species, Formica paralugubris, which has numerous queens per colony. Our results show that workers do not discriminate between nestmate and foreign queens, and more generally accept new queens at a limited rate. To conclude, this thesis shows that mechanisms influencing variation in colony queen number and the influence of these changes on conflict resolution are complex. Data gathered in this thesis therefore constitute a solid background for further research on the evolution and the maintenance of complex organisations in insect societies.

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.

Little effect of seasonal constraints on population genetic structure in eusocial paper wasps.

Climate has long been suggested to affect population genetic structures of eusocial insect societies. For instance, Hamilton [Journal of Theoretical Biology7 (1964) 17] discusses whether temperate and tropical eusocial insects may show differences in population-level genetic structure and viscosity, and how this might relate to differences in the degree of synchrony in their life cycles or modes of nest founding. Despite the importance of Hamilton's 1964 papers, this specific idea has not been tested in actual populations of wasps, probably due to the paucity of studies on tropical species. Here, we compare colony and population genetic structures in two species of primitively eusocial paper wasps with contrasting ecologies: the tropical species Polistes canadensis and the temperate species P. dominulus. Our results provide important clarifications of Hamilton's discussion. Specifically, we show that the genetic structures of the temperate and tropical species were very similar, indicating that seasonality does not greatly affect population viscosity or inbreeding. For both species, the high genetic differentiation between nests suggests strong selection at the nest level to live with relatives, whereas low population viscosity and low genetic differentiation between nest aggregations might reflect balancing selection to disperse, avoiding competition with relatives. Overall, our study suggests no prevalence of seasonal constraints of the life cycle in affecting the population genetic structure of eusocial paper wasps. These conclusions are likely to apply also to other primitively eusocial insects, such as halictine bees. They also highlight how selection for a kin structure that promotes altruism can override potential effects of ecology in eusocial insects.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.

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.

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.

Nature versus nurture in social insect caste determination

Recent evidence for genetic effects on royal and worker caste differentiation from diverse social insect taxa has put an end to the view that these phenotypes stem solely from a developmental switch controlled by environmental factors. Instead, the relative influences of genotypic and environmental effects on caste vary among species, ranging from largely environmentally controlled phenotypes to almost purely genetic systems. Disentangling the selective forces that generate variation for caste predisposition will require characterizing the genetic mechanisms underlying this variation, and identifying particular life-history strategies and kin structures associated with strong genetic effects on caste.

Alternative reproductive strategies in Formica lugubris Zett (Hymenoptera, Formicidae)

A long term study on the reproductive strategies of Formica lugubris (Formica rufa group) has revealed many new important facts. Within a supercolony the sex-ratio is strongly female biased. Some of the sexuals mate on the nest surface, whereas others fly from the nest. Those sexuals participating in flights aggregate in particular meadow sites (mating places). On these places females attract males by releasing a sex pheromone.

Unifying niche shift studies: insights from biological invasions.

Assessing whether the climatic niche of a species may change between different geographic areas or time periods has become increasingly important in the context of ongoing global change. However, approaches and findings have remained largely controversial so far, calling for a unification of methods. Here, we build on a review of empirical studies of invasion to formalize a unifying framework that decomposes niche change into unfilling, stability, and expansion situations, taking both a pooled range and range-specific perspective on the niche, while accounting for climatic availability and climatic analogy. This framework provides new insights into the nature of climate niche shifts and our ability to anticipate invasions, and may help in guiding the design of experiments for assessing causes of niche changes.

Genesis and evolution of the Evx and Mox genes and the extended Hox and ParaHox gene clusters

Background: Hox and ParaHox gene clusters are thought to have resulted from the duplication of a ProtoHox gene cluster early in metazoan evolution. However, the origin and evolution of the other genes belonging to the extended Hox group of homeobox-containing genes, that is, Mox and Evx, remains obscure. We constructed phylogenetic trees with mouse, amphioxus and Drosophila extended Hox and other related Antennapedia-type homeobox gene sequences and analyzed the linkage data available for such genes.Results: We claim that neither Mox nor Evx is a Hox or ParaHox gene. We propose a scenariothat reconciles phylogeny with linkage data, in which an Evx/Mox ancestor gene linked to aProtoHox cluster was involved in a segmental tandem duplication event that generated an arrayof all Hox-like genes, referred to as the `coupled¿ cluster. A chromosomal breakage within thiscluster explains the current composition of the extended Hox cluster (with Evx, Hox and Moxgenes) and the ParaHox cluster.Conclusions: Most studies dealing with the origin and evolution of Hox and ParaHox clustershave not included the Hox-related genes Mox and Evx. Our phylogenetic analyses and theavailable linkage data in mammalian genomes support an evolutionary scenario in which anancestor of Evx and Mox was linked to the ProtoHox cluster, and that a tandem duplication of alarge genomic region early in metazoan evolution generated the Hox and ParaHox clusters, plusthe cluster-neighbors Evx and Mox. The large `coupled¿ Hox-like cluster EvxHox/MoxParaHox wassubsequently broken, thus grouping the Mox and Evx genes to the Hox clusters, and isolating theParaHox cluster.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.