Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Evolution of the correlation between expression divergence and protein divergence in mammals.

Divergence of protein sequences and gene expression patterns are two fundamental mechanisms that generate organismal diversity. Here, we have used genome and transcriptome data from eight mammals and one bird to study the positive correlation of these two processes throughout mammalian evolution. We demonstrate that the correlation is stable over time and most pronounced in neural tissues, which indicates that it is the result of strong negative selection. The correlation is not driven by genes with specific functions and may instead best be viewed as an evolutionary default state, which can nevertheless be evaded by certain gene types. In particular, genes with developmental and neural functions are skewed toward changes in gene expression, consistent with selection against pleiotropic effects associated with changes in protein sequences. Surprisingly, we find that the correlation between expression divergence and protein divergence is not explained by between-gene variation in expression level, tissue specificity, protein connectivity, or other investigated gene characteristics, suggesting that it arises independently of these gene traits. The selective constraints on protein sequences and gene expression patterns also fluctuate in a coordinate manner across phylogenetic branches: We find that gene-specific changes in the rate of protein evolution in a specific mammalian lineage tend to be accompanied by similar changes in the rate of expression evolution. Taken together, our findings highlight many new aspects of the correlation between protein divergence and expression divergence, and attest to its role as a fundamental property of mammalian genome evolution.

Gene duplication, population genomics, and species-level differentiation within a tropical mountain shrub.

Gene duplication leads to paralogy, which complicates the de novo assembly of genotyping-by-sequencing (GBS) data. The issue of paralogous genes is exacerbated in plants, because they are particularly prone to gene duplication events. Paralogs are normally filtered from GBS data before undertaking population genomics or phylogenetic analyses. However, gene duplication plays an important role in the functional diversification of genes and it can also lead to the formation of postzygotic barriers. Using populations and closely related species of a tropical mountain shrub, we examine 1) the genomic differentiation produced by putative orthologs, and 2) the distribution of recent gene duplication among lineages and geography. We find high differentiation among populations from isolated mountain peaks and species-level differentiation within what is morphologically described as a single species. The inferred distribution of paralogs among populations is congruent with taxonomy and shows that GBS could be used to examine recent gene duplication as a source of genomic differentiation of nonmodel species.

Aalenian to Cenomanian Radiolaria of the Bermeja Complex (Puerto Rico) and Pacific origin of radiolarites on the Caribbean Plate

The study of the radiolarian ribbon chert is a key in determining the origins of associated Mesozoic oceanic terranes and may help to achieve a general agreement regarding the basic principles on the evolution of the Caribbean Plate. The Bermeja Complex of Puerto Rico, which contains serpentinized peridotite, altered basalt, amphibolite, and chert (Mariquita Chert Formation), is one of these crucial oceanic terranes. The radiolarian biochronology presented in this work is mainly based by correlation on the biozonations of Baumgartner et al. (1995) and O'Dogherty (1994) and indicates an early Middle Jurassic to early Late Cretaceous (late Bajocian-early Callovian to late early Albian-early middle Cenomanian) age. The illustrated assemblages contain about 120 species, of which one is new (Pantanellium karinae), and belonging to about 50 genera. A review of the previous radiolarian published works on the Mariquita Chert Formation and the results of this study suggest that this formation ranges in age from Middle Jurassic to early Late Cretaceous (late Aalenian to early-middle Cenomanian) and also reveal a possible feature of the Bermeja Complex, which is the younging of radiolarian cherts from north to south, evoking a polarity of accretion. On the basis of a currently exhaustive inventory of the radiolarite facies s.s. on the Caribbean Plate, a re-examination of the regional distribution of Middle Jurassic sediments associated with oceanic crust, and a paleoceanographic argumentation on the water currents, we come to the conclusion that the radiolarite and associated Mesozoic oceanic terranes of the Caribbean Plate are of Pacific origin. Eventually, a discussion on the origin of the cherts of the Mariquita Formation illustrated by Middle Jurassic to middle Cretaceous geodynamic models of the Pacific and Caribbean realms bring up the possibility that the rocks of the Bermeja Complex are remnants of two different oceans.

Evolutionary patterns of MHC class II B in owls and their implications for the understanding of avian MHC evolution

Owing to its special mode of evolution and central role in the adaptive immune system, the major histocompatibility complex (MHC) has become the focus of diverse disciplines such as immunology, evolutionary ecology, and molecular evolution. MHC evolution has been studied extensively in diverse vertebrate lineages over the last few decades, and it has been suggested that birds differ from the established mammalian norm. Mammalian MHC genes evolve independently, and duplication history (i.e., orthology) can usually be traced back within lineages. In birds, this has been observed in only 3 pairs of closely related species. Here we report strong evidence for the persistence of orthology of MHC genes throughout an entire avian order. Phylogenetic reconstructions of MHC class II B genes in 14 species of owls trace back orthology over tens of thousands of years in exon 3. Moreover, exon 2 sequences from several species show closer relationships than sequences within species, resembling transspecies evolution typically observed in mammals. Thus, although previous studies suggested that long-term evolutionary dynamics of the avian MHC was characterized by high rates of concerted evolution, resulting in rapid masking of orthology, our results question the generality of this conclusion. The owl MHC thus opens new perspectives for a more comprehensive understanding of avian MHC evolution.

Genetic causes of transitions from sexual reproduction to asexuality in plants and animals.

The persistence of sexual reproduction in the face of competition from asexual invaders is more likely if asexual lineages are produced infrequently or have low fitness. The generation rate and success of new asexual lineages will be influenced by the proximate mechanisms underlying transitions to asexuality. As such, characterization of these mechanisms can help explain the distribution of reproductive modes among natural populations. Here, we synthesize the literature addressing proximate causes of transitions from sexual to asexual reproduction in plants and animals. In cyclical and facultatively asexual taxa, individual mutations can cause obligate asexuality. The evolution of asexuality in obligately sexual groups is more complex, requiring the simultaneous acquisition of two traits generally controlled by different genetic factors: unreduced gamete formation and spontaneous development of unfertilized gametes. At least three 'pre-adaptations' could favour transitions to obligate asexuality in obligate sexuals. First, linkage among loci affecting separate key components of asexuality facilitates its spread, with evidence for these linkage blocks in plants. Second, asexuality should evolve more readily in haplodiploids; support for this hypothesis comes from two examples where a single locus causes transitions to asexuality. Third, standing genetic variation for the production of unreduced gametes could facilitate transitions to asexuality, but whether the ability to produce unreduced gametes contributes to the evolution of obligate asexuality remains unclear. We close by reviewing the associations between asexuality, hybridization and polyploidy, and argue that current data suggest that hybridization is more likely to play a causal role in transitions to asexuality than polyploidy.

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.

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Analyses of six homologous proteins of Protochlamydia amoebophila UWE25 encoded by large GC-rich genes (lgr): a model of evolution and concatenation of leucine-rich repeats.

BACKGROUND: Along the chromosome of the obligate intracellular bacteria Protochlamydia amoebophila UWE25, we recently described a genomic island Pam100G. It contains a tra unit likely involved in conjugative DNA transfer and lgrE, a 5.6-kb gene similar to five others of P. amoebophila: lgrA to lgrD, lgrF. We describe here the structure, regulation and evolution of these proteins termed LGRs since encoded by "Large G+C-Rich" genes. RESULTS: No homologs to the whole protein sequence of LGRs were found in other organisms. Phylogenetic analyses suggest that serial duplications producing the six LGRs occurred relatively recently and nucleotide usage analyses show that lgrB, lgrE and lgrF were relocated on the chromosome. The C-terminal part of LGRs is homologous to Leucine-Rich Repeats domains (LRRs). Defined by a cumulative alignment score, the 5 to 18 concatenated octacosapeptidic (28-meric) LRRs of LGRs present all a predicted alpha-helix conformation. Their closest homologs are the 28-residue RI-like LRRs of mammalian NODs and the 24-meres of some Ralstonia and Legionella proteins. Interestingly, lgrE, which is present on Pam100G like the tra operon, exhibits Pfam domains related to DNA metabolism. CONCLUSION: Comparison of the LRRs, enable us to propose a parsimonious evolutionary scenario of these domains driven by adjacent concatenations of LRRs. Our model established on bacterial LRRs can be challenged in eucaryotic proteins carrying less conserved LRRs, such as NOD proteins and Toll-like receptors.

Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

Total lipid and fatty acid accumulation during basidiospore formation in the ectomycorrhizal fungus Pisolithus sp.

The basidiospores of Pisolithus sp. contain large amounts of lipids, indicating provision for future germination in the host rhizosphere. However, the accumulation, composition, and mobilization of lipids during formation and germination of these spores are largely unknown. In this study, lipid storage and fatty acid composition during basidiosporogenesis were analyzed in fresh basidiocarps using bright-field microscopy and gas chromatography. Abundant lipid bodies are found in the hyphae, basidia, and basidiospores of fungal basidiocarps. This evidences a considerable C transport in the basidiocarp to meet the C demand during basidiospore formation. Fatty acid composition analysis revealed the presence of 24 compounds with chains of 9 to 18 C atoms, either saturated or insaturated, with one or two insaturations. The fatty acid composition and content varied according to the developmental stage of the peridioles. In free basidiospores, the predominant compounds were 16:0, 16:1w5c, 18:1w9c, and 18:2w6,9c/18:0ante, at concentrations of 76, 46, 192, and 51 µg g-1 dry matter, respectively. Our results indicate that oleic acid is the major constituent of lipid reserves in Pisolithus sp. basidiospores. Further studies are being conducted to determine the factors that induce lipid mobilization during spore germination.

Biofilm formation on bone grafts and bone graft substitutes: comparison of different materials by a standard in vitro test and microcalorimetry.

We analyzed the initial adhesion and biofilm formation of Staphylococcus aureus (ATCC 29213) and S. epidermidis RP62A (ATCC 35984) on various bone grafts and bone graft substitutes under standardized in vitro conditions. In parallel, microcalorimetry was evaluated as a real-time microbiological assay in the investigation of biofilm formation and material science research. The materials beta-tricalcium phosphate (beta-TCP), processed human spongiosa (Tutoplast) and poly(methyl methacrylate) (PMMA) were investigated and compared with polyethylene (PE). Bacterial counts (log(10) cfu per sample) were highest on beta-TCP (S. aureus 7.67 +/- 0.17; S. epidermidis 8.14 +/- 0.05) while bacterial density (log(10) cfu per surface) was highest on PMMA (S. aureus 6.12 +/- 0.2, S. epidermidis 7.65 +/- 0.13). Detection time for S. aureus biofilms was shorter for the porous materials (beta-TCP and processed human spongiosa, p < 0.001) compared to the smooth materials (PMMA and PE), with no differences between beta-TCP and processed human spongiosa (p > 0.05) or PMMA and PE (p > 0.05). In contrast, for S. epidermidis biofilms the detection time was different (p < 0.001) between all materials except between processed human spongiosa and PE (p > 0.05). The quantitative analysis by quantitative culture after washing and sonication of the material demonstrated the importance of monitoring factors like specific surface or porosity of the test materials. Isothermal microcalorimetry proved to be a suitable tool for an accurate, non-invasive and real-time microbiological assay, allowing the detection of bacterial biomass without removing the biofilm from the surface.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Iron oxidation on the surface of adventitious roots and its relation to aerenchyma formation in rice genotypes

Establishment of the water layer in an irrigated rice crop leads to consumption of free oxygen in the soil which enters in a chemical reduction process mediated by anaerobic microorganisms, changing the crop environment. To maintain optimal growth in an environment without O2, rice plants develop pore spaces (aerenchyma) that allow O2 transport from air to the roots. Carrying capacity is determined by the rice genome and it may vary among cultivars. Plants that have higher capacity for formation of aerenchyma should theoretically carry more O2 to the roots. However, part of the O2 that reaches the roots is lost due to permeability of the roots and the O2 gradient created between the soil and roots. The O2 that is lost to the outside medium can react with chemically reduced elements present in the soil; one of them is iron, which reacts with oxygen and forms an iron plaque on the outer root surface. Therefore, evaluation of the iron plaque and of the formation of pore spaces on the root can serve as a parameter to differentiate rice cultivars in regard to the volume of O2 transported via aerenchyma. An experiment was thus carried out in a greenhouse with the aim of comparing aerenchyma and iron plaque formation in 13 rice cultivars grown in flooded soils to their formation under growing conditions similar to a normal field, without free oxygen. The results indicated significant differences in the volume of pore spaces in the roots among cultivars and along the root segment in each cultivar, indicating that under flooded conditions the genetic potential of the plant is crucial in induction of cell death and formation of aerenchyma in response to lack of O2. In addition, the amount of Fe accumulated on the root surface was different among genotypes and along the roots. Thus, we concluded that the rice genotypes exhibit different responses for aerenchyma formation, oxygen release by the roots and iron plaque formation, and that there is a direct relationship between porosity and the amount of iron oxidized on the root surface.

Large N phases, gravitational instantons and the nuts and bolts of AdS holography

Recent results in the literature concerning holography indicate that the thermodynamics of quantum gravity (at least with a negative cosmological constant) can be modeled by the large N thermodynamics of quantum field theory. We emphasize that this suggests a completely unitary evolution of processes in quantum gravity, including black hole formation and decay, and even more extreme examples involving topology change. As concrete examples which show that this correspondence holds even when the space-time is only locally asymptotically AdS, we compute the thermodynamical phase structure of the AdS-Taub-NUT and AdS-Taub-bolt spacetimes, and compare them to a (2+1)-dimensional conformal field theory (at large N) compactified on a squashed three-sphere and on the twisted plane.