Evolution based on domain combinations: the case of glutaredoxins

Background: Protein domains represent the basic units in the evolution of proteins. Domain duplication and shuffling by recombination and fusion, followed by divergence are the most common mechanisms in this process. Such domain fusion and recombination events are predicted to occur only once for a given multidomain architecture. However, other scenarios may be relevant in the evolution of specific proteins, such as convergent evolution of multidomain architectures. With this in mind, we study glutaredoxin (GRX) domains, because these domains of approximately one hundred amino acids are widespread in archaea, bacteria and eukaryotes and participate in fusion proteins. GRXs are responsible for the reduction of protein disulfides or glutathione-protein mixed disulfides and are involved in cellular redox regulation, although their specific roles and targets are often unclear. Results: In this work we analyze the distribution and evolution of GRX proteins in archaea,bacteria and eukaryotes. We study over one thousand GRX proteins, each containing at least one GRX domain, from hundreds of different organisms and trace the origin and evolution of the GRX domain within the tree of life. Conclusion: Our results suggest that single domain GRX proteins of the CGFS and CPYC classes have, each, evolved through duplication and divergence from one initial gene that was present in the last common ancestor of all organisms. Remarkably, we identify a case of convergent evolution in domain architecture that involves the GRX domain. Two independent recombination events of a TRX domain to a GRX domain are likely to have occurred, which is an exception to the dominant mechanism of domain architecture evolution.

Genomic evidence for role of inversion 3RP of Drosophila melanogaster in facilitating climate change adaptation.

Chromosomal inversion polymorphisms are common in animals and plants, and recent models suggest that alternative arrangements spread by capturing different combinations of alleles acting additively or epistatically to favour local adaptation. It is also thought that inversions typically maintain favoured combinations for a long time by suppressing recombination between alternative chromosomal arrangements. Here, we consider patterns of linkage disequilibrium and genetic divergence in an old inversion polymorphism in Drosophila melanogaster (In(3R)Payne) known to be associated with climate change adaptation and a recent invasion event into Australia. We extracted, karyotyped and sequenced whole chromosomes from two Australian populations, so that changes in the arrangement of the alleles between geographically separated tropical and temperate areas could be compared. Chromosome-wide linkage disequilibrium (LD) analysis revealed strong LD within the region spanned by In(3R)Payne. This genomic region also showed strong differentiation between the tropical and the temperate populations, but no differentiation between different karyotypes from the same population, after controlling for chromosomal arrangement. Patterns of differentiation across the chromosome arm and in gene ontologies were enhanced by the presence of the inversion. These data support the notion that inversions are strongly selected by bringing together combinations of genes, but it is still not clear if such combinations act additively or epistatically. Our data suggest that climatic adaptation through inversions can be dynamic, reflecting changes in the relative abundance of different forms of an inversion and ongoing evolution of allelic content within an inversion.

Diversity and genetic struture in cajá tree (spondias mombin l.) populations in northeastern Brazil

Spondias mombin L. is a fruit tree from the American continent from the Anacardiaceae family. In Brazil it is common in different vegetation types but is more frequent in the Atlantic and Amazonian rainforests. It is economically important because of its fruits, which are widely consumed raw or processed as fruit jellies, juices and ice creams. The leaves have great importance in the pharmaceutical industry because of their antibacterial properties. In the state of Pernambuco, cajá tree is widely distributed in the Zona da Mata region and less frequently in the Agreste and Sertão areas. In this work diversity and genetic structure were studied in four populations of cajá tree from Pernambuco's Zona da Mata, Northeast Brazil, using isozymes polymorphism analyses from electrophoreses. The result showed 100% of polymorphism (P) for nine alleles and the average of alleles per locus s was 2.4. The expected heterozygosity ranged from 0.530 to 0.574 and the observed heterozygosity , from 0.572 to 0.735. It was not observed inbreeding and the average F IT was -0.175, whereas within population inbreeding (f) varied from -0.08 to- 0.37. The genetic divergence among the populations (F ST) ranged from 0.006 to 0.028 and the average was 0.026. The average of estimated gene flow (Nm) was high (5.27). The CG-IPA population, corresponding to the germplasm collection of IPA, showed more than 96% of genetic similarity with other populations; therefore, it is a good representative of the existent genetic diversity in the Zona da Mata region.

The evolution of alternative splicing patterns and regulatory mechanisms

Alternative splicing produces multiple isoforms from the same gene, thus increasing the number of transcripts of the species. Alternative splicing is a virtually ubiquitous mechanism in eukaryotes, for example more than 90% of protein-coding genes in human are alternatively spliced. Recent evolutionary studies showed that alternative splicing is a fast evolving and highly species- specific mechanism. The rapid evolution of alternative splicing was considered as a contribution to the phenotypic diversity between species. However, the function of many isoforms produced by alternative splicing remains unclear and they might be the result of noisy splicing. Thus, the functional relevance of alternative splicing and the evolutionary mechanisms of its rapid divergence among species are still poorly understood. During my thesis, I performed a large-scale analysis of the regulatory mechanisms that drive the rapid evolution of alternative splicing. To study the evolution of alternative splicing regulatory mechanisms, I used an extensive RNA-sequencing dataset comprising 12 tetrapod species (human, chimpanzee and bonobo, gorilla, orangutan, macaque, marmoset, mouse, opossum, platypus, chicken and frog) and 8 tissues (cerebellum, brain, heart, kidney, liver, testis, placenta and ovary). To identify the catalogue of alternative splicing eis-acting regulatory elements in the different tetrapod species, I used a previously defined computational approach. This approach is a statistical analysis of exons/introns and splice sites composition and relies on a principle of compensation between splice sites strength and the presence of additional regulators. With an evolutionary comparative analysis of the exonic eis-acting regulators, I showed that these regulatory elements are generally shared among primates and more conserved than non-regulatory elements. In addition, I showed that the usage of these regulatory elements is also more conserved than expected by chance. In addition to the identification of species- specific eis-acting regulators, these results may explain the rapid evolution of alternative splicing. I also developed a new approach based on evolutionary sequence changes and corresponding alternative splicing changes to identify potential splicing eis-acting regulators in primates. The identification of lineage-specific substitutions and corresponding lineage-specific alternative splicing changes, allowed me to annotate the genomic sequences that might have played a role in the alternative splicing pattern differences among primates. Finally, I showed that the identified splicing eis-acting regulator datasets are enriched in human disease-causing mutations, thus confirming their biological relevance.

Vanishing chromosomal inversion clines in Drosophila subobscura from Chile: Is behavioral thermoregulation to blame?

Chromosomal inversion clines paralleling the long-standing ones in native Palearctic populations of Drosophila subobscura evolved swiftly after this species invaded the Americas in the late 1970s and early 1980s. However, the new clines did not consistently continue to converge on the Old World baseline. Our recent survey of Chilean populations of D. subobscura shows that inversion clines have faded or even changed sign with latitude. Here, we investigate the hypothesis that this fading of inversion clines might be due to the Bogert effect, namely, that flies' thermoregulatory behavior has eventually compensated for environmental variation in temperature, thus buffering selection on thermal-related traits. We show that latitudinal divergence in thermal preference (T-p) has evolved in Chile for females, with higher-latitude flies having a lower mean T-p. Plastic responses in T-p also lessen latitudinal thermal variation because flies developed at colder temperatures prefer warmer microclimates. Our results are consistent with the idea that active behavioral thermoregulation might buffer environmental variation and reduce the potential effect of thermal selection on other traits as chromosomal arrangements.

The UlaG protein familly defines novel structural and functional motifs grafted on an ancient RNase fold

Background: Bacterial populations are highly successful at colonizing new habitats and adapting to changing environmental conditions, partly due to their capacity to evolve novel virulence and metabolic pathways in response to stress conditions and to shuffle them by horizontal gene transfer (HGT). A common theme in the evolution of new functions consists of gene duplication followed by functional divergence. UlaG, a unique manganese-dependent metallo-b-lactamase (MBL) enzyme involved in L-ascorbate metabolism by commensal and symbiotic enterobacteria, provides a model for the study of the emergence of new catalytic activities from the modification of an ancient fold. Furthermore, UlaG is the founding member of the so-called UlaG-like (UlaGL) protein family, a recently established and poorly characterized family comprising divalent (and perhaps trivalent)metal-binding MBLs that catalyze transformations on phosphorylated sugars and nucleotides. Results: Here we combined protein structure-guided and sequence-only molecular phylogenetic analyses to dissect the molecular evolution of UlaG and to study its phylogenomic distribution, its relatedness with present-day UlaGL protein sequences and functional conservation. Phylogenetic analyses indicate that UlaGL sequences are present in Bacteria and Archaea, with bona fide orthologs found mainly in mammalian and plant-associated Gramnegative and Gram-positive bacteria. The incongruence between the UlaGL tree and known species trees indicates exchange by HGT and suggests that the UlaGL-encoding genes provided a growth advantage under changing conditions. Our search for more distantly related protein sequences aided by structural homology has uncovered that UlaGL sequences have a common evolutionary origin with present-day RNA processing and metabolizing MBL enzymes widespread in Bacteria, Archaea, and Eukarya. This observation suggests an ancient origin for the UlaGL family within the broader trunk of the MBL superfamily by duplication, neofunctionalization and fixation. Conclusions: Our results suggest that the forerunner of UlaG was present as an RNA metabolizing enzyme in the last common ancestor, and that the modern descendants of that ancestral gene have a wide phylogenetic distribution and functional roles. We propose that the UlaGL family evolved new metabolic roles among bacterial and possibly archeal phyla in the setting of a close association with metazoans, such as in the mammalian gastrointestinal tract or in animal and plant pathogens, as well as in environmental settings. Accordingly, the major evolutionary forces shaping the UlaGL family include vertical inheritance and lineage-specific duplication and acquisition of novel metabolic functions, followed by HGT and numerous lineage-specific gene loss events.

Patterns of nucleotide diversity at the regions encompassing the Drosophila insulin-like peptide (dilp) genes: demography vs positive selection in Drosophila melanogaster.

In Drosophila, the insulin-signaling pathway controls some life history traits, such as fertility and lifespan, and it is considered to be the main metabolic pathway involved in establishing adult body size. Several observations concerning variation in body size in the Drosophila genus are suggestive of its adaptive character. Genes encoding proteins in this pathway are, therefore, good candidates to have experienced adaptive changes and to reveal the footprint of positive selection. The Drosophila insulin-like peptides (DILPs) are the ligands that trigger the insulin-signaling cascade. In Drosophila melanogaster, there are several peptides that are structurally similar to the single mammalian insulin peptide. The footprint of recent adaptive changes on nucleotide variation can be unveiled through the analysis of polymorphism and divergence. With this aim, we have surveyed nucleotide sequence variation at the dilp1-7 genes in a natural population of D. melanogaster. The comparison of polymorphism in D. melanogaster and divergence from D. simulans at different functional classes of the dilp genes provided no evidence of adaptive protein evolution after the split of the D. melanogaster and D. simulans lineages. However, our survey of polymorphism at the dilp gene regions of D. melanogaster has provided some evidence for the action of positive selection at or near these genes. The regions encompassing the dilp1-4 genes and the dilp6 gene stand out as likely affected by recent adaptive events.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group).

BACKGROUND: Hybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group. RESULTS: We found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts. CONCLUSIONS: Our findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

Molecular phylogenetics and temporal diversification in the genus Aeromonas based on the sequences of five housekeeping genes

Several approaches have been developed to estimate both the relative and absolute rates of speciation and extinction within clades based on molecular phylogenetic reconstructions of evolutionary relationships, according to an underlying model of diversification. However, the macroevolutionary models established for eukaryotes have scarcely been used with prokaryotes. We have investigated the rate and pattern of cladogenesis in the genus Aeromonas (γ-Proteobacteria, Proteobacteria, Bacteria) using the sequences of five housekeeping genes and an uncorrelated relaxed-clock approach. To our knowledge, until now this analysis has never been applied to all the species described in a bacterial genus and thus opens up the possibility of establishing models of speciation from sequence data commonly used in phylogenetic studies of prokaryotes. Our results suggest that the genus Aeromonas began to diverge between 248 and 266 million years ago, exhibiting a constant divergence rate through the Phanerozoic, which could be described as a pure birth process.

Annonaceae substitution rates: a codon model perspective

The Annonaceae includes cultivated species of economic interest and represents an important source of information for better understanding the evolution of tropical rainforests. In phylogenetic analyses of DNA sequence data that are used to address evolutionary questions, it is imperative to use appropriate statistical models. Annonaceae are cases in point: Two sister clades, the subfamilies Annonoideae and Malmeoideae, contain the majority of Annonaceae species diversity. The Annonoideae generally show a greater degree of sequence divergence compared to the Malmeoideae, resulting in stark differences in branch lengths in phylogenetic trees. Uncertainty in how to interpret and analyse these differences has led to inconsistent results when estimating the ages of clades in Annonaceae using molecular dating techniques. We ask whether these differences may be attributed to inappropriate modelling assumptions in the phylogenetic analyses. Specifically, we test for (clade-specific) differences in rates of non-synonymous and synonymous substitutions. A high ratio of nonsynonymous to synonymous substitutions may lead to similarity of DNA sequences due to convergence instead of common ancestry, and as a result confound phylogenetic analyses. We use a dataset of three chloroplast genes (rbcL, matK, ndhF) for 129 species representative of the family. We find that differences in branch lengths between major clades are not attributable to different rates of non-synonymous and synonymous substitutions. The differences in evolutionary rate between the major clades of Annonaceae pose a challenge for current molecular dating techniques that should be seen as a warning for the interpretation of such results in other organisms.

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions

Phylogeny and biogeography of Primula sect. Armerina: implications for plant evolution under climate change and the uplift of the Qinghai-Tibet Plateau.

BACKGROUND: The historical orogenesis and associated climatic changes of mountain areas have been suggested to partly account for the occurrence of high levels of biodiversity and endemism. However, their effects on dispersal, differentiation and evolution of many groups of plants are still unknown. In this study, we examined the detailed diversification history of Primula sect. Armerina, and used biogeographic analysis and macro-evolutionary modeling to investigate a series of different questions concerning the evolution of the geographical and ecological distribution of the species in this section. RESULTS: We sequenced five chloroplast and one nuclear genes for species of Primula sect. Armerina. Neither chloroplast nor nuclear trees support the monophyly of the section. The major incongruences between the two trees occur among closely related species and may be explained by hybridization. Our dating analyses based on the chloroplast dataset suggest that this section began to diverge from its relatives around 3.55 million years ago, largely coinciding with the last major uplift of the Qinghai-Tibet Plateau (QTP). Biogeographic analysis supports the origin of the section in the Himalayan Mountains and dispersal from the Himalayas to Northeastern QTP, Western QTP and Hengduan Mountains. Furthermore, evolutionary models of ecological niches show that the two P. fasciculata clades have significantly different climatic niche optima and rates of niche evolution, indicating niche evolution under climatic changes and further providing evidence for explaining their biogeographic patterns. CONCLUSION: Our results support the hypothesis that geologic and climatic events play important roles in driving biological diversification of organisms in the QTP area. The Pliocene uplift of the QTP and following climatic changes most likely promoted both the inter- and intraspecific divergence of Primula sect. Armerina. This study also illustrates how niche evolution under climatic changes influences biogeographic patterns.

Climate niche conservatism does not explain restricted distribution patterns in Tynanthus (Bignonieae, Bignoniaceae)

Studies on niche evolution allow us to establish how species niches have changed over time as well as to identify how long-term evolutionary processes have led to present-day species distributions. Here, we investigate the patterns of climatic niche evolution in Tynanthus (Bignonieae, Bignoniaceae), a genus comprising narrowly distributed species. We test the hypothesis that niche conservatism has played an important role in the diversification history of this group of Neotropical lianas. For that, we perform univariate and multivariate comparisons between species' climatic niches and associated environmental data with information on species' phylogenetic relationships. We encountered considerable divergence in niches among species, indicating that niche conservatism in climatic variables has does not seem to havenot played a key role in the diversification of the genus. Our results are used as a basis to discuss patterns of ecological niche evolution in the group and to suggest novel approaches for future analyses.

HETEROTIC GROUP FORMATION IN PSIDIUM GUAJAVA L. BY ARTIFICIAL NEURAL NETWORK AND DISCRIMINANT ANALYSIS

ABSTRACT The present study aimed at evaluating the heterotic group formation in guava based on quantitative descriptors and using artificial neural network (ANN). For such, we evaluated eight quantitative descriptors. Large genetic variability was found for the eight quantitative traits in the 138 genotypes of guava. The artificial neural network technique determined that the optimal number of groups was three. The grouping consistency was determined by linear discriminant analysis, which obtained classification percentage of the groups, with a value of 86 %. It was concluded that the artificial neural network method is effective to detect genetic divergence and heterotic group formation.

The dynamics of coexistence: habitat sharing versus segregation patterns among three sympatric montane vipers

Contact zones of closely related and ecologically similar species constitute rare opportunities to study the evolutionary consequences of past speciation processes. They represent natural laboratories in which strong competition could lead to the exclusion of one species, or the various species may switch into distinct ecological niches. Alternatively, if reproductive isolation has not yet been achieved, they may hybridize. We elucidate the degree of taxon integrity by comparing genetics and habitat use of three similar-sized congeneric viper species, Vipera ammodytes, Viperaaspis, and Viperaberus, of Nadiza Valley in western Slovenia. No hybridization was detected for either mitochondrial or nuclear genomes. Similarly, external intermediacy by a single prestudy viper (probably V.ammodytesxV. aspis) indicates that hybridization occasionally occurs, but should be very rare. Populations of the three related viperids are partially allopatric in Nadiza Valley, but they also coexist in a narrow contact zone in the montane grassland along the south-exposed slope of Mount Stol (1673m a.s.l.). Here, the three species that occupy areas in or near patches of rocky microhabitats (e.g. stone piles, slides, and walls) live in syntopy. However, fine-scale measurements of structural components show partial habitat segregation, in which V.berus becomes more dominant at elevations above 1400m and occupies mostly the mountain ridge and north-exposed slopes of Mount Stol, V.aspis occurs below 1300m and is the only species to inhabit stoneless patches of grass and bushes around 1000m and lower, and V.ammodytes occurs at all elevations up to 1500m, but is restricted to a rocky microhabitat. We suggest that a high degree of microstructure divergence, slightly different environmental niches, and a generally favourable habitat for all three viper species, keep the pressure for mis-mating and hybridization low, although mechanisms such as reduced hybrid inferiority and temporal mating segregation cannot yet be excluded.