Evolution of a new function by degenerative mutation in cephalochordate steroid receptors

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers-molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here, we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), an α2β2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition were simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous, binding via the same surface patch on interacting subunits, but rotated 180° relative to each other. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric and homomeric interactions. Our findings suggest that elaborate and specific symmetrical molecular complexes may often evolve via simple genetic and physical mechanisms.

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers - molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), aα 2 β 2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition were simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous - involving the same surface patch on interacting subunits, rotated 180° relative to each other. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric vs homomeric interactions. Our findings suggest that elaborate and specific symmetrical molecular complexes may often evolve via simple genetic and physical mechanisms.

The interplay of oxytocin and sex hormones

The neuropeptide oxytocin has historically been associated with reproduction and maternal behavior. However, more recent research has uncovered that oxytocin has a much wider range of roles in physiology and behavior. Despite the excitement surrounding potential therapeutical applications of intranasally administered oxytocin, the results of these intervention studies have been inconsistent. Various reasons for these mixed results have been proposed, which tend to focus on methodological issues, such as study design. While methodological issues are certainly important, emerging evidence suggests that the interaction between oxytocin and sex hormones may also account for these varied findings. To better understand the purpose and function of the interaction of oxytocin with sex hormones, with a focus on estrogens, progesterone, and testosterone, we conducted a comprehensive thematic review via four perspectives: evolutionary, developmental, mechanistic, and survival. Altogether, this synergistic approach highlights the critical function of sex hormone activity for accomplishing the diverse roles of oxytocin via the modulation of oxytocin release and oxytocin receptor activity, which is also likely to contribute to the heterogeneity of outcomes after oxytocin administration.

Implication of Molecular Constraints Facilitating the Functional Evolution of Pseudomonas aeruginosa KPR2 into a Versatile α-Keto-Acid Reductase

Theoretical concepts linking the structure, function, and evolution of a protein, while often intuitive, necessitate validation through investigations in real-world systems. Our study empirically explores the evolutionary implications of multiple gene copies in an organism by shedding light on the structure-function modulations observed in Pseudomonas aeruginosa's second copy of ketopantoate reductase (PaKPR2). We demonstrated with two apo structures that the typical active site cleft of the protein transforms into a two-sided pocket where a molecular gate made up of two residues controls the substrate entry site, resulting in its inactivity toward the natural substrate ketopantoate. Strikingly, this structural modification made the protein active against several important α-keto-acid substrates with varied efficiency. Structural constraints at the binding site for this altered functional trait were analyzed with two binary complexes that show the conserved residue microenvironment faces restricted movements due to domain closure. Finally, its mechanistic highlights gathered from a ternary complex structure help in delineating the molecular perspectives behind its kinetic cooperativity toward these broad range of substrates. Detailed structural characteristics of the protein presented here also identified four key amino acid residues responsible for its versatile α-keto-acid reductase activity, which can be further modified to improve its functional properties through protein engineering.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Functional dissection and assembly of a small, newly evolved, W chromosome-specific genomic region of the African clawed frog Xenopus laevis

Genetic triggers for sex determination are frequently co-inherited with other linked genes that may also influence one or more sex-specific phenotypes. To better understand how sex-limited regions evolve and function, we studied a small W chromosome-specific region of the frog Xenopus laevis that contains only three genes (dm-w, scan-w, ccdc69-w) and that drives female differentiation. Using gene editing, we found that the sex-determining function of this region requires dm-w but that scan-w and ccdc69-w are not essential for viability, female development, or fertility. Analysis of mesonephros+gonad transcriptomes during sexual differentiation illustrates masculinization of the dm-w knockout transcriptome, and identifies mostly non-overlapping sets of differentially expressed genes in separate knockout lines for each of these three W-specific gene compared to wildtype sisters. Capture sequencing of almost all Xenopus species and PCR surveys indicate that the female-determining function of dm-w is present in only a subset of species that carry this gene. These findings map out a dynamic evolutionary history of a newly evolved W chromosome-specific genomic region, whose components have distinctive functions that frequently degraded during Xenopus diversification, and evidence the evolutionary consequences of recombination suppression.

Comparison and quantification of estrogen receptor-mediated responsiveness to endocrine disruptors in bivalves by using complementary model and a novel yeast assay approach

Endocrine disrupting chemicals (EDCs) in estuaries and coastal habitats have been widely detected over the world and caused global concern. Bivalves have been shown to be vulnerable to endocrine disruption. However, estrogen receptors (ERs) sensitivity to steroids and EDCs has long been considered to be restricted to vertebrates. In the present study, a computational simulation docking model was applied to qualitatively predict the binding behavior of two bivalve ERs to estradiol and compared the docking activity with zebra fish ERa. A novel reconstituted yeast system was constructed by using transcriptional activator GAL-4 consists of ER-expressing plasmid and ERE (estrogen responsive element)-containing plasmid. The assays showed that bivalve ER specifically activate transcription in response to tested steroids and EDCs, but the activation ability is weaker compared to zebra fish ERa. The results corroborate the presence of an active ER in bivalve molluscs and provide a promising tool for screening of marine environmental pollutants active in disturbing ERs of bivalves, as well as understanding the underlying mechanism across taxonomic groups and phyla.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

Divergent Evolution of Progesterone and Mineralocorticoid Receptors in Terrestrial Vertebrates and Fish Influences Endocrine Disruption

There is much concern about disruption of endocrine physiology regulated by steroid hormones in humans, other terrestrial vertebrates and fish by industrial chemicals, such as bisphenol A, and pesticides, such as DDT. These endocrine-disrupting chemicals influence steroid-mediated physiology in humans and other vertebrates by competing with steroids for receptor binding sites, disrupting diverse responses involved in reproduction, development and differentiation. Here I discuss that due to evolution of the progesterone receptor (PR) and mineralocorticoid receptor (MR) after ray-finned fish and terrestrial vertebrates diverged from a common ancestor, each receptor evolved to respond to different steroids in ray-finned fish and terrestrial vertebrates. In elephant shark, a cartilaginous fish that diverged before the separation between ray-finned fish and terrestrial vertebrates, both progesterone and 17,20β-dihydroxy-progesterone activate the PR. During the evolution of ray-finned fish and terrestrial vertebrates, the PR in terrestrial vertebrates continued responding to progesterone and evolved to weakly respond to 17,20β-dihydroxy-progesterone. In contrast, the physiological progestin for the PR in zebrafish and other ray-finned fish is 17,20β-dihydroxy-progesterone, and ray-finned fish PR responds weakly to progesterone. The MR in fish and terrestrial vertebrates also diverged to have different responses to progesterone. Progesterone is a potent agonist for elephant shark MR, zebrafish MR and other fish MRs, in contrast to progesterone's opposite activity as an antagonist for aldosterone, the physiological mineralocorticoid for human MR. These different physiological ligands for fish and terrestrial vertebrate PR and MR need to be considered in applying data for their disruption by chemicals in fish and terrestrial vertebrates to each other.

Yap5 Competes With Hap4 for the Regulation of Iron Homeostasis Genes in the Human Pathogen Candida glabrata

The CCAAT-binding complex (CBC) is a conserved heterotrimeric transcription factor which, in fungi, requires additional regulatory subunits to act on transcription. In the pathogenic yeast Candida glabrata, CBC has a dual role. Together with the Hap4 regulatory subunit, it activates the expression of genes involved in respiration upon growth with non-fermentable carbon sources, while its association with the Yap5 regulatory subunit is required for the activation of iron tolerance genes in response to iron excess. In the present work, we investigated further the interplay between CBC, Hap4 and Yap5. We showed that Yap5 regulation requires a specific Yap Response Element in the promoter of its target gene GRX4 and that the presence of Yap5 considerably strengthens the binding of CBC to the promoters of iron tolerance genes. Chromatin immunoprecipitation (ChIP) and transcriptome experiments showed that Hap4 can also bind these promoters but has no impact on the expression of those genes when Yap5 is present. In the absence of Yap5 however, GRX4 is constitutively regulated by Hap4, similarly to the genes involved in respiration. Our results suggest that the distinction between the two types of CBC targets in C. glabrata is mainly due to the dependency of Yap5 for very specific DNA sequences and to the competition between Hap4 and Yap5 at the promoter of the iron tolerance genes.

Identification of discriminative gene-level and protein-level features associated with pathogenic gain-of-function and loss-of-function variants

Identifying whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity is an important step toward understanding disease mechanisms because they may result in markedly different clinical phenotypes. Here, we generated an extensive database of documented germline GOF and LOF pathogenic variants by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. We then investigated various gene- and protein-level features of GOF and LOF variants and applied machine learning and statistical analyses to identify discriminative features. We found that GOF variants were enriched in essential genes, for autosomal-dominant inheritance, and in protein binding and interaction domains, whereas LOF variants were enriched in singleton genes, for protein-truncating variants, and in protein core regions. We developed a user-friendly web-based interface that enables the extraction of selected subsets from the GOF/LOF database by a broad set of annotated features and downloading of up-to-date versions. These results improve our understanding of how variants affect gene/protein function and may ultimately guide future treatment options.

A Crab Is Not a Fish: Unique Aspects of the Crustacean Endocrine System and Considerations for Endocrine Toxicology

Crustaceans-and arthropods in general-exhibit many unique aspects to their physiology. These include the requirement to moult (ecdysis) in order to grow and reproduce, the ability to change color, and multiple strategies for sexual differentiation. Accordingly, the endocrine regulation of these processes involves hormones, receptors, and enzymes that differ from those utilized by vertebrates and other non-arthropod invertebrates. As a result, environmental chemicals known to disrupt endocrine processes in vertebrates are often not endocrine disruptors in crustaceans; while, chemicals that disrupt endocrine processes in crustaceans are often not endocrine disruptors in vertebrates. In this review, we present an overview of the evolution of the endocrine system of crustaceans, highlight endocrine endpoints known to be a target of disruption by chemicals, and identify other components of endocrine signaling that may prove to be targets of disruption. This review highlights that crustaceans need to be evaluated for endocrine disruption with consideration of their unique endocrine system and not with consideration of the endocrine system of vertebrates.

Development of an activity assay for characterizing deoxyhypusine synthase and its diverse reaction products

Deoxyhypusine synthase transfers an aminobutyl moiety from spermidine to the eukaryotic translation initiation factor 5A (eIF5A) in the first step of eIF5A activation. This exclusive post-translational modification is conserved in all eukaryotes. Activated eIF5A has been shown to be essential for cell proliferation and viability. Recent reports have linked the activation of eIF5A to several human diseases. Deoxyhypusine synthase, which is encoded by a single gene copy in most eukaryotes, was duplicated in several plant lineages during evolution, the copies being repeatedly recruited to pyrrolizidine alkaloid biosynthesis. However, the function of many of these duplicates is unknown. Notably, deoxyhypusine synthase is highly promiscuous and can catalyze various reactions, often of unknown biological relevance. To facilitate in-depth biochemical studies of this enzyme, we report here the development of a simple and robust in vitro enzyme assay. It involves precolumn derivatization of the polyamines taking part in the reaction and avoids the need for the previously used radioactively labeled tracers. The derivatized polyamines are quantified after high-performance liquid chromatography coupled to diode array and fluorescence detectors. By performing kinetic analyses of deoxyhypusine synthase and its paralog from the pyrrolizidine alkaloid-producing plant Senecio vernalis, we demonstrate that the assay unequivocally differentiates the paralogous enzymes. Furthermore, it detects and quantifies, in a single assay, the side reactions that occur in parallel to the main reaction. The presented assay thus provides a detailed biochemical characterization of deoxyhypusine synthase and its paralogs.

Similarities in biological processes can be used to bridge ecology and molecular biology

Much of the research in biology aims to understand the origin of diversity. Naturally, ecological diversity was the first object of study, but we now have the necessary tools to probe diversity at molecular scales. The inherent differences in how we study diversity at different scales caused the disciplines of biology to be organized around these levels, from molecular biology to ecology. Here, we illustrate that there are key properties of each scale that emerge from the interactions of simpler components and that these properties are often shared across different levels of organization. This means that ideas from one level of organization can be an inspiration for novel hypotheses to study phenomena at another level. We illustrate this concept with examples of events at the molecular level that have analogs at the organismal or ecological level and vice versa. Through these examples, we illustrate that biological processes at different organization levels are governed by general rules. The study of the same phenomena at different scales could enrich our work through a multidisciplinary approach, which should be a staple in the training of future scientists.

A Review on Sex Steroid Hormone Estrogen Receptors in Mammals and Fish

Steroid hormones play essential roles in the reproductive biology of vertebrates. Estrogen exercises its effect through estrogen receptors and is not only a female reproductive hormone but acts virtually in all vertebrates, including fish, and is involved in the physiological and pathological states in all males and females. Estrogen has been implicated in mandible conservation and circulatory and central nervous systems as well as the reproductive system. This review intended to understand the structure, function, binding affinities, and activations of estrogens and estrogen receptors and to discuss the understanding of the role of sex steroid hormone estrogen receptors in mammals and fish.

Steroid receptors and vertebrate evolution

Considering that life on earth evolved about 3.7 billion years ago, vertebrates are young, appearing in the fossil record during the Cambrian explosion about 542 to 515 million years ago. Results from sequence analyses of genomes from bacteria, yeast, plants, invertebrates and vertebrates indicate that receptors for adrenal steroids (aldosterone, cortisol), and sex steroids (estrogen, progesterone, testosterone) also are young, with an estrogen receptor and a 3-ketosteroid receptor first appearing in basal chordates (cephalochordates: amphioxus), which are close ancestors of vertebrates. Duplication and divergence of the 3-ketosteroid receptor yielded an ancestral progesterone receptor and an ancestral corticoid receptor, the common ancestor of the glucocorticoid and mineralocorticoid receptors, in jawless vertebrates (cyclostomes: lampreys, hagfish). This was followed by evolution of an androgen receptor, distinct glucocorticoid and mineralocorticoid receptors and estrogen receptor-α and -β in cartilaginous fishes (Chondrichthyes: sharks). Further evolution of mineralocorticoid signaling occurred with the evolution of aldosterone synthase in lungfish, a forerunner of terrestrial vertebrates. Adrenal and sex steroid receptors are not found in echinoderms and hemichordates, which are ancestors in the lineage of cephalochordates and vertebrates. The evolution of steroid receptors at key nodes in the evolution of vertebrates, in which steroid receptors act as master switches to regulate differentiation, development, reproduction, immune responses, electrolyte homeostasis and stress responses, suggests an important role for steroid receptors in the evolutionary success of vertebrates, considering that the human genome contains about 22,000 genes, which is not much larger than genomes of invertebrates, such as Caenorhabditis elegans (~18,000 genes) and Drosophila (~14,000 genes).

Paralog dependency indirectly affects the robustness of human cells

The protective redundancy of paralogous genes partly relies on the fact that they carry their functions independently. However, a significant fraction of paralogous proteins may form functionally dependent pairs, for instance, through heteromerization. As a consequence, one could expect these heteromeric paralogs to be less protective against deleterious mutations. To test this hypothesis, we examined the robustness landscape of gene loss-of-function by CRISPR-Cas9 in more than 450 human cell lines. This landscape shows regions of greater deleteriousness to gene inactivation as a function of key paralog properties. Heteromeric paralogs are more likely to occupy such regions owing to their high expression and large number of protein-protein interaction partners. Further investigation revealed that heteromers may also be under stricter dosage balance, which may also contribute to the higher deleteriousness upon gene inactivation. Finally, we suggest that physical dependency may contribute to the deleteriousness upon loss-of-function as revealed by the correlation between the strength of interactions between paralogs and their higher deleteriousness upon loss of function.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

Genomic sequence analyses of classical and non-classical lamprey progesterone receptor genes and the inference of homologous gene evolution in metazoans

BACKGROUND

Nuclear progesterone receptor (nPR) is an evolutionary innovation in vertebrates that mediates genomic responses to progesterone. Vertebrates also respond to progesterone via membrane progesterone receptors (mPRs) or membrane associated progesterone receptors (MAPRs) through rapid nongenomic mechanisms. Lampreys are extant agnathan vertebrates, residing at the evolutionary juncture where vertebrates diverged from invertebrates. A survey of the progesterone receptor (PR) gene sequences in lamprey genomes would inform PR gene evolutionary events during the transition from invertebrates to vertebrates.

RESULTS

In this study, we annotated sequences of one nPR, four mPR (β, γ, δ and ε) and four MAPR genes from genomes of two lamprey species (Petromyzon marinus and Lethenteron japonicum). To infer the origin and evolutionary history of PR genes, we constructed phylogenetic trees of PR homologous sequences across representative species of metazoans. Phylogenetic analyses revealed that the mPRγ gene first appeared in non-bilaterians, and the mPRβ gene likely arose from a duplication of mPRγ. On the other hand, the mPRγ gene gave rise to the mPRδ and ε genes much later in the vertebrate lineage. In addition, the mPRα gene first appeared in cartilaginous fishes, likely derived from duplication of mPRβ after the agnathan-gnathostome divergence. All known MAPR genes were present in the lamprey genomes. Progesterone receptor membrane component 1 (PGRMC1), neudesin and neuferricin genes probably evolved in parallel in non-bilaterians, whereas two copies of PGRMC genes probably derived from duplication of ancestral PGRMC1 sequence and appeared before the speciation of lampreys.

CONCLUSIONS

Non-classical mPR and MAPR genes first evolved in non-bilaterians and classical nPR genes evolved later in basal vertebrates. Sequence repertoires for membrane progesterone receptor genes in vertebrates likely originated from an ancestral metazoan sequence and expanded via several duplication events.

Fatty-acylation target sequence in the ligand-binding domain of vertebrate steroid receptors demarcates evolution from estrogen-related receptors

Present-day nuclear receptors (NRs) responding to adrenal and sex steroids are key regulators of reproduction and growth in mammals, and are thought to have evolved from an ancestral NR most closely related to extant estrogen-related receptors (ERRs). The molecular events (and ligands) that distinguish steroid-activated NRs (SRs) from their inferred ancestor, that gave rise to both the ERRs and SRs, remain unknown. We report that target sequences for fatty-acylation (palmitoylation) at a key cysteine residue (corresponding to Cys447 in human estrogen receptor ERα) in helix 8 of the ligand-binding domain accurately demarcate SRs from ERRs. Docking studies are consistent with the hypothesis that palmitate embeds into a key groove in the receptor surface. The implications of lipidation, and of potential alternative ligands for the key cysteine residue, for receptor function and the evolution of SRs are discussed.

The promiscuous estrogen receptor: Evolution of physiological estrogens and response to phytochemicals and endocrine disruptors

Many actions of estradiol (E2), the principal physiological estrogen in vertebrates, are mediated by estrogen receptor-α (ERα) and ERβ. An important physiological feature of vertebrate ERs is their promiscuous response to several physiological steroids, including estradiol (E2), Δ⁵-androstenediol, 5α-androstanediol, and 27-hydroxycholesterol. A novel structural characteristic of Δ⁵-androstenediol, 5α-androstanediol, and 27-hydroxycholesterol is the presence of a C19 methyl group, which precludes the presence of an aromatic A ring with a C3 phenolic group that is a defining property of E2. The structural diversity of these estrogens can explain the response of the ER to synthetic chemicals such as bisphenol A and DDT, which disrupt estrogen physiology in vertebrates, and the estrogenic activity of a variety of plant-derived chemicals such as genistein, coumestrol, and resveratrol. Diversity in the A ring of physiological estrogens also expands potential structures of industrial chemicals that can act as endocrine disruptors. Compared to E2, synthesis of 27-hydroxycholesterol and Δ⁵-androstenediol is simpler, leading us, based on parsimony, to propose that one or both of these steroids or a related metabolite was a physiological estrogen early in the evolution of the ER, with E2 assuming this role later as the canonical estrogen. In addition to the well-studied role of the ER in reproductive physiology, the ER also is an important transcription factor in non-reproductive tissues such as the cardiovascular system, kidney, bone, and brain. Some of these ER actions in non-reproductive tissues appeared early in vertebrate evolution, long before the emergence of mammals.

Increased Affinity for RNA Targets Evolved Early in Animal and Plant Dicer Lineages through Different Structural Mechanisms

Understanding the structural basis for evolutionary changes in protein function is central to molecular evolutionary biology and can help determine the extent to which functional convergence occurs through similar or different structural mechanisms. Here, we combine ancestral sequence reconstruction with functional characterization and structural modeling to directly examine the evolution of sequence-structure-function across the early differentiation of animal and plant Dicer/DCL proteins, which perform the first molecular step in RNA interference by identifying target RNAs and processing them into short interfering products. We found that ancestral Dicer/DCL proteins evolved similar increases in RNA target affinities as they diverged independently in animal and plant lineages. In both cases, increases in RNA target affinities were associated with sequence changes that anchored the RNA’s 5′phosphate, but the structural bases for 5′phosphate recognition were different in animal versus plant lineages. These results highlight how molecular-functional evolutionary convergence can derive from the evolution of unique protein structures implementing similar biochemical mechanisms.

The effects of repeated whole genome duplication events on the evolution of cytokinin signaling pathway

BACKGROUND

It is thought that after whole-genome duplications (WGDs), a large fraction of the duplicated gene copies is lost over time while few duplicates are retained. Which factors promote survival or death of a duplicate remains unclear and the underlying mechanisms are poorly understood. According to the model of gene dosage balance, genes encoding interacting proteins are predicted to be preferentially co-retained after WGDs. Among these are genes encoding proteins involved in complexes or in signal transduction.

RESULTS

We have investigated the way that repeated WGDs during land plant evolution have affected cytokinin signaling to study patterns of gene duplicability and co-retention in this important signal transduction pathway. Through the integration of phylogenetic analyses with comparisons of genome collinearity, we have found that signal input mediated by cytokinin receptors proved to be highly conserved over long evolutionary time-scales, with receptors showing predominantly gene loss after repeated WGDs. However, the downstream elements, e,g. response regulators, were mainly retained after WGDs and thereby formed gene families in most plant lineages.

CONCLUSIONS

Gene dosage balance between the interacting components indicated by co-retention after WGDs seems to play a minor role in the evolution of cytokinin signaling pathway. Overall, core genes of cytokinin signaling show a highly heterogeneous pattern of gene retention after WGD, reflecting complex relationships between the various factors that shape the long-term fate of a duplicated gene.

LXRα and LXRβ Nuclear Receptors Evolved in the Common Ancestor of Gnathostomes

Nuclear receptors (NRs) regulate numerous aspects of the endocrine system. They mediate endogenous and exogenous cues, ensuring a homeostatic control of development and metabolism. Gene duplication, loss and mutation have shaped the repertoire and function of NRs in metazoans. Here, we examine the evolution of a pivotal orchestrator of cholesterol metabolism in vertebrates, the liver X receptors (LXRs). Previous studies suggested that LXRα and LXRβ genes emerged in the mammalian ancestor. However, we show through genome analysis and functional assay that bona fide LXRα and LXRβ orthologues are present in reptiles, coelacanth and chondrichthyans but not in cyclostomes. These findings show that LXR duplicated before gnathostome radiation, followed by asymmetric paralogue loss in some lineages. We suggest that a tighter control of cholesterol levels in vertebrates was achieved through the exploitation of a wider range of oxysterols, an ability contingent on ligand-binding pocket remodeling.

Genome-wide identification of nuclear receptor (NR) genes and the evolutionary significance of the NR1O subfamily in the monogonont rotifer Brachionus spp

Nuclear receptors (NRs) are a large family of transcription factors that are involved in many fundamental biological processes. NRs are considered to have originated from a common ancestor, and are highly conserved throughout the whole animal taxa. Therefore, the genome-wide identification of NR genes in an animal taxon can provide insight into the evolutionary tendencies of NRs. Here, we identified all the NR genes in the monogonont rotifer Brachionus spp., which are considered an ecologically key species due to their abundance and world-wide distribution. The NR family was composed of 40, 32, 29, and 32 genes in the genomes of the rotifers B. calyciflorus, B. koreanus, B. plicatilis, and B. rotundiformis, respectively, which were classified into seven distinct subfamilies. The composition of each subfamily was highly conserved between species, except for NR1O genes, suggesting that they have undergone sporadic evolutionary processes for adaptation to their different environmental pressures. In addition, despite the dynamics of NR evolution, the significance of the conserved endocrine system, particularly for estrogen receptor (ER)-signaling, in rotifers was discussed on the basis of phylogenetic analyses. The results of this study may help provide a better understanding the evolution of NRs, and expand our knowledge of rotifer endocrine systems.

Gene duplication can impart fragility, not robustness, in the yeast protein interaction network

The maintenance of duplicated genes is thought to protect cells from genetic perturbations, but the molecular basis of this robustness is largely unknown. By measuring the interaction of yeast proteins with their partners in wild-type cells and in cells lacking a paralog, we found that 22 out of 56 paralog pairs compensate for the lost interactions. An equivalent number of pairs exhibit the opposite behavior and require each other's presence for maintaining their interactions. These dependent paralogs generally interact physically, regulate each other's abundance, and derive from ancestral self-interacting proteins. This reveals that gene duplication may actually increase mutational fragility instead of robustness in a large number of cases.

The evolution, structure and function of the ray finned fish (Actinopterygii) glucocorticoid receptors

Basal ray-finned fish (Actinopterygii) possess a single glucocorticoid receptor (GR) and when compared to the lobe-finned vertebrate (Sarcopterygii) GR possess nine additional amino acids between the zinc-finger of the DNA binding domain. A whole genome duplication event which occurred between 320 and 350MYA in the teleost lineage following the split from the basal ray-finned fish resulted in 2 GRs: one GR group, GR1, has retained the 9 amino acids insert whereas the other group, GR2, has not. The exception to this is the zebrafish, that have lost one of the GRs, but they do possess 2 GRs with a splice variant that lacks the C-terminal portion of the GR to form GRβ which acts as a dominant-repressor of the wildtype GR. Another splice variant sees the basal ray-finned GR and teleost GR1 without the 9 amino acids insert. The molecular basis for GRs retention is beginning to be unravelled. In Pantadon buchholzi, rainbow trout, carp, marine and Japanese medaka GR2 is more sensitive to glucocorticoids (GC), thus potentially playing a more significant role in regulating gene expression at basal circulatory GC concentrations. However, this division in GC sensitivity is not seen in other species. The few studies to evaluate the significance of the 9 amino acid insert have shown that it affect maximal transactivational activity the extent to which is dependent on the number of glucocorticoid response elements (GREs) present in the reporter plasmid. The retention of these GRs would suggest there was an evolutionary advantage, which saw the development of a complex regulatory process to mediate the actions of the glucocorticoids.

30 YEARS OF THE MINERALOCORTICOID RECEPTOR: Evolution of the mineralocorticoid receptor: sequence, structure and function

The mineralocorticoid receptor (MR) is descended from a corticoid receptor (CR), which has descendants in lamprey and hagfish, cyclostomes (jawless fish), a taxon that evolved at the base of the vertebrate line. A distinct MR and GR first appear in cartilaginous fishes (Chondrichthyes), such as sharks, skates, rays and chimeras. Skate MR has a strong response to corticosteroids that are mineralocorticoids and glucocorticoids in humans. The half-maximal responses (EC50s) for skate MR for the mineralocorticoids aldosterone and 11-deoxycorticosterone are 0.07 nM and 0.03 nM, respectively. EC50s for the glucocorticoids cortisol and corticosterone are 1 nM and 0.09 nM, respectively. The physiological mineralocorticoid in ray-finned fish, which do not synthesize aldosterone, is not fully understood because several 3-ketosteroids, including cortisol, 11-deoxycortisol, corticosterone, 11-deoxycorticosterone and progesterone are transcriptional activators of fish MR. Further divergence of the MR and GR in terrestrial vertebrates, which synthesize aldosterone, led to emergence of aldosterone as a selective ligand for the MR. Here, we combine sequence analysis of the CR and vertebrate MRs and GRs, analysis of crystal structures of human MR and GR and data on transcriptional activation by 3-ketosteroids of wild-type and mutant MRs and GRs to investigate the evolution of selectivity for 3-ketosteroids by the MR in terrestrial vertebrates and ray-finned fish, as well as the basis for binding of some glucocorticoids by human MR and other vertebrate MRs.

Convergence of Domain Architecture, Structure, and Ligand Affinity in Animal and Plant RNA-Binding Proteins

Reconstruction of ancestral protein sequences using phylogenetic methods is a powerful technique for directly examining the evolution of molecular function. Although ancestral sequence reconstruction (ASR) is itself very efficient, downstream functional, and structural studies necessary to characterize when and how changes in molecular function occurred are often costly and time-consuming, currently limiting ASR studies to examining a relatively small number of discrete functional shifts. As a result, we have very little direct information about how molecular function evolves across large protein families. Here we develop an approach combining ASR with structure and function prediction to efficiently examine the evolution of ligand affinity across a large family of double-stranded RNA binding proteins (DRBs) spanning animals and plants. We find that the characteristic domain architecture of DRBs-consisting of 2-3 tandem double-stranded RNA binding motifs (dsrms)-arose independently in early animal and plant lineages. The affinity with which individual dsrms bind double-stranded RNA appears to have increased and decreased often across both animal and plant phylogenies, primarily through convergent structural mechanisms involving RNA-contact residues within the β1-β2 loop and a small region of α2. These studies provide some of the first direct information about how protein function evolves across large gene families and suggest that changes in molecular function may occur often and unassociated with major phylogenetic events, such as gene or domain duplications.

Origin of an ancient hormone/receptor couple revealed by resurrection of an ancestral estrogen

The origin of ancient ligand/receptor couples is often analyzed via reconstruction of ancient receptors and, when ligands are products of metabolic pathways, they are not supposed to evolve. However, because metabolic pathways are inherited by descent with modification, their structure can be compared using cladistic analysis. Using this approach, we studied the evolution of steroid hormones. We show that side-chain cleavage is common to most vertebrate steroids, whereas aromatization was co-opted for estrogen synthesis from a more ancient pathway. The ancestral products of aromatic activity were aromatized steroids with a side chain, which we named "paraestrols." We synthesized paraestrol A and show that it effectively binds and activates the ancestral steroid receptor. Our study opens the way to comparative studies of biologically active small molecules.

Molecular cloning and characterization of hagfish estrogen receptors

One or more distinct forms of the nuclear estrogen receptor (ER) have been isolated from many vertebrates to date. To better understand the molecular evolution of ERs, we cloned and characterized er cDNAs from the inshore hagfish, Eptatretus burgeri, a modern representative of the most primitive vertebrates, the agnathans. Two er cDNAs, er1 and er2, were isolated from the liver of a reproductive female hagfish. A phylogenetic analysis placed hagfish ER1 into a position prior to the divergence of vertebrate ERs. Conversely, hagfish ER2 was placed at the base of the vertebrate ERβ clade. The tissue distribution patterns of both ER subtype mRNAs appeared to be different, suggesting that each subtype has different physiological roles associated with estrogen actions. An estrogen responsive-luciferase reporter assay using mammalian HEK293 cells was used to functionally characterize these hagfish ERs. Both ER proteins displayed estrogen-dependent activation of transcription. These results clearly demonstrate that the hagfish has two functional ER subtypes.

Resurrecting ancestral structural dynamics of an antiviral immune receptor: adaptive binding pocket reorganization repeatedly shifts RNA preference

Background

Although resurrecting ancestral proteins is a powerful tool for understanding the molecular-functional evolution of gene families, nearly all studies have examined proteins functioning in relatively stable biological processes. The extent to which more dynamic systems obey the same ‘rules’ governing stable processes is unclear. Here we present the first detailed investigation of the functional evolution of the RIG-like receptors (RLRs), a family of innate immune receptors that detect viral RNA in the cytoplasm.

Results

Using kinetic binding assays and molecular dynamics simulations of ancestral proteins, we demonstrate how a small number of adaptive protein-coding changes repeatedly shifted the RNA preference of RLRs throughout animal evolution by reorganizing the shape and electrostatic distribution across the RNA binding pocket, altering the hydrogen bond network between the RLR and its RNA target. In contrast to observations of proteins involved in metabolism and development, we find that RLR-RNA preference ‘flip flopped’ between two functional states, and shifts in RNA preference were not always coupled to gene duplications or speciation events. We demonstrate at least one reversion of RLR-RNA preference from a derived to an ancestral function through a novel structural mechanism, indicating multiple structural implementations of similar functions.

Conclusions

Our results suggest a model in which frequent shifts in selection pressures imposed by an evolutionary arms race preclude the long-term functional optimization observed in stable biological systems. As a result, the evolutionary dynamics of immune receptors may be less constrained by structural epistasis and historical contingency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0818-6) contains supplementary material, which is available to authorized users.

Evaluation of the scientific underpinnings for identifying estrogenic chemicals in nonmammalian taxa using mammalian test systems

The US Environmental Protection Agency has responsibility for assessing endocrine activity of more than 10 000 chemicals, a task that cannot reasonably be achieved solely through use of available mammalian and nonmammalian in vivo screening assays. Hence, it has been proposed that chemicals be prioritized for in vivo testing using data from in vitro high-throughput assays for specific endocrine system targets. Recent efforts focused on potential estrogenic chemicals-specifically those that activate estrogen receptor-alpha (ERα)-have broadly demonstrated feasibility of the approach. However, a major uncertainty is whether prioritization based on mammalian (primarily human) high-throughput assays accurately reflects potential chemical-ERα interactions in nonmammalian species. The authors conducted a comprehensive analysis of cross-species comparability of chemical-ERα interactions based on information concerning structural attributes of estrogen receptors, in vitro binding and transactivation data for ERα, and the effects of a range of chemicals on estrogen-signaling pathways in vivo. Overall, this integrated analysis suggests that chemicals with moderate to high estrogenic potency in mammalian systems also should be priority chemicals in nonmammalian vertebrates. However, the degree to which the prioritization approach might be applicable to invertebrates is uncertain because of a lack of knowledge of the biological role(s) of possible ERα orthologs found in phyla such as annelids. Further, comparative analysis of in vitro data for fish and reptiles suggests that mammalian-based assays may not effectively capture ERα interactions for low-affinity chemicals in all vertebrate classes. Environ Toxicol Chem 2016;35:2806-2816. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US Government work and, as such, is in the public domain in the United States of America.

Potential mechanisms underlying estrogen-induced expression of the molluscan estrogen receptor (ER) gene

In vertebrates, estrogens and estrogen mimicking chemicals modulate gene expression mainly through a genomic pathway mediated by the estrogen receptors (ERs). Although the existence of an ER orthologue in the mollusc genome has been known for some time, its role in estrogen signalling has yet to be deciphered. This is largely due to its constitutive (ligand-independent) activation and a limited mechanistic understanding of its regulation. To fill this knowledge gap, we cloned and characterised an ER cDNA (sgER) and the 5'-flanking region of the gene from the Sydney rock oyster Saccostrea glomerata. The sgER cDNA is predicted to encode a 477-amino acid protein that contains a DNA-binding domain (DBD) and a ligand-binding domain (LBD) typically conserved among both vertebrate and invertebrate ERs. A comparison of the sgER LBD sequence with those of other ligand-dependent ERs revealed that the sgER LBD is variable at several conserved residues known to be critical for ligand binding and receptor activation. Ligand binding assays using fluorescent-labelled E2 and purified sgER protein confirmed that sgER is devoid of estrogen binding. In silico analysis of the sgER 5'-flanking sequence indicated the presence of three putative estrogen responsive element (ERE) half-sites and several putative sites for ER-interacting transcription factors, suggesting that the sgER promoter may be autoregulated by its own gene product. sgER mRNA is ubiquitously expressed in adult oyster tissues, with the highest expression found in the ovary. Ovarian expression of sgER mRNA was significantly upregulated following in vitro and in vivo exposure to 17β-estradiol (E2). Notably, the activation of sgER expression by E2 in vitro was abolished by the specific ER antagonist ICI 182, 780. To determine whether sgER expression is epigenetically regulated, the in vivo DNA methylation status of the putative proximal promoter in ovarian tissues was assessed using bisulfite genomic sequencing. The results showed that the promoter is predominantly hypomethylated (with 0-3.3% methylcytosines) regardless of sgER mRNA levels. Overall, our investigations suggest that the estrogen responsiveness of sgER is regulated by a novel ligand-dependent receptor, presumably via a non-genomic pathway(s) of estrogen signalling.

A second estrogen receptor from Japanese lamprey (Lethenteron japonicum) does not have activities for estrogen binding and transcription

Estrogens regulate many physiological responses in vertebrates by binding to the estrogen receptor (ER), a ligand-activated transcription factor. To understand the evolution of vertebrate ERs and to investigate how estrogen acts in a jawless vertebrate, we used degenerate primer sets and PCR to isolate DNA fragments encoding two distinct ER subtypes, Esr1a and Esr1b from the Japanese lamprey, Lethenteron japonicum. Phylogenetic analysis indicates that these two ERs are the result of lineage-specific gene duplication within the jawless fishes, different from the previous duplication event of Esr1 (ERα) and Esr2 (ERβ) within the jawed vertebrates. Reporter gene assays show that lamprey Esr1a displays both constitutive and estrogen-dependent activation of gene transcription. Domain swapping experiments indicate that constitutive activity resides in the A/B domain of lamprey Esr1a. Unexpectedly, lamprey Esr1b does not bind estradiol and is not stimulated by other estrogens, androgens or corticosteroids. A 3D model of lamprey Esr1b suggests that although estradiol fits into the steroid binding site, some stabilizing contacts between the ligand and side chains that are found in human Esr1 and Esr2 are missing in lamprey Esr1b.

Evolution of Gene Duplication in Plants

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

PhyloBot: A Web Portal for Automated Phylogenetics, Ancestral Sequence Reconstruction, and Exploration of Mutational Trajectories

The method of phylogenetic ancestral sequence reconstruction is a powerful approach for studying evolutionary relationships among protein sequence, structure, and function. In particular, this approach allows investigators to (1) reconstruct and “resurrect” (that is, synthesize in vivo or in vitro) extinct proteins to study how they differ from modern proteins, (2) identify key amino acid changes that, over evolutionary timescales, have altered the function of the protein, and (3) order historical events in the evolution of protein function. Widespread use of this approach has been slow among molecular biologists, in part because the methods require significant computational expertise. Here we present PhyloBot, a web-based software tool that makes ancestral sequence reconstruction easy. Designed for non-experts, it integrates all the necessary software into a single user interface. Additionally, PhyloBot provides interactive tools to explore evolutionary trajectories between ancestors, enabling the rapid generation of hypotheses that can be tested using genetic or biochemical approaches. Early versions of this software were used in previous studies to discover genetic mechanisms underlying the functions of diverse protein families, including V-ATPase ion pumps, DNA-binding transcription regulators, and serine/threonine protein kinases. PhyloBot runs in a web browser, and is available at the following URL: http://www.phylobot.com. The software is implemented in Python using the Django web framework, and runs on elastic cloud computing resources from Amazon Web Services. Users can create and submit jobs on our free server (at the URL listed above), or use our open-source code to launch their own PhyloBot server.

Mechanistic insights into induction of vitellogenin gene expression by estrogens in Sydney rock oysters, Saccostrea glomerata

Marine molluscs, such as oysters, respond to estrogenic compounds with the induction of the egg yolk protein precursor, vitellogenin (Vtg), availing a biomarker for estrogenic pollution. Despite this application, the precise molecular mechanism through which estrogens exert their action to induce molluscan vitellogenesis is unknown. As a first step to address this question, we cloned a gene encoding Vtg from the Sydney rock oyster Saccostrea glomerata (sgVtg). Using primers designed from a partial sgVtg cDNA sequence available in Genbank, a full-length sgVtg cDNA of 8498bp was obtained by 5'- and 3'-RACE. The open reading frame (ORF) of sgVtg was determined to be 7980bp, which is substantially longer than the orthologs of other oyster species. Its deduced protein sequence shares the highest homology at the N- and C-terminal regions with other molluscan Vtgs. The full-length genomic DNA sequence of sgVtg was obtained by genomic PCR and genome walking targeting the gene body and flanking regions, respectively. The genomic sequence spans 20kb and consists of 30 exons and 29 introns. Computer analysis identified three closely spaced half-estrogen responsive elements (EREs) in the promoter region and a 210-bp CpG island 62bp downstream of the transcription start site. Upregulation of sgVtg mRNA expression was observed in the ovaries following in vitro (explants) and in vivo (tank) exposure to 17β-estradiol (E2). Notably, treatment with an estrogen receptor (ER) antagonist in vitro abolished the upregulation, suggesting a requirement for an estrogen-dependent receptor for transcriptional activation. DNA methylation of the 5' CpG island was analysed using bisulfite genomic sequencing of the in vivo exposed ovaries. The CpG island was found to be hypomethylated (with 0-3% methylcytosines) in both control and E2-exposed oysters. However, no significant differential methylation or any correlation between methylation and sgVtg expression levels was observed. Overall, the results support the possible involvement of an ERE-containing promoter and an estrogen-activated receptor in estrogen signalling in marine molluscs.

Evolution of estrogen receptors in ray-finned fish and their comparative responses to estrogenic substances

In vertebrates, estrogens play fundamental roles in regulating reproductive activities through estrogen receptors (ESRs), and disruption of estrogen signaling is now of global concern for both wildlife and human health. To date, ESRs of only a limited number of species have been characterized. We investigated the functional diversity and molecular basis or ligand sensitivity of ESRs among ray-finned fish species (Actinopterygii), the most variable group within vertebrates. We cloned and characterized ESRs from several key species in the evolution of ray-finned fish including bichir (Polypteriformes, ESR1 and ESR2) at the basal lineage of ray-finned fish, and arowana (Osteoglossiformes, ESR1 and ESR2b) and eel (Anguilliformes, ESR1, ESR2a and ESR2b) both belonging to ancient early-branching lineages of teleosts, and suggest that ESR2a and ESR2b emerged through teleost-specific whole genome duplication, but an ESR1 paralogue has been lost in the early lineage of euteleost fish species. All cloned ESR isoforms showed similar responses to endogenous and synthetic steroidal estrogens, but they responded differently to non-steroidal estrogenic endocrine disrupting chemicals (EDCs) (e.g., ESR2a exhibits a weaker reporter activity compared with ESR2b). We show that variation in ligand sensitivity of ESRs can be attributed to phylogeny among species of different taxonomic groups in ray-finned fish. The molecular information provided contributes both to understanding of the comparative role of ESRs in the reproductive biology of fish and their comparative responses to EDCs.

Gene Duplicability of Core Genes Is Highly Consistent across All Angiosperms

Gene duplication is an important mechanism for adding to genomic novelty. Hence, which genes undergo duplication and are preserved following duplication is an important question. It has been observed that gene duplicability, or the ability of genes to be retained following duplication, is a nonrandom process, with certain genes being more amenable to survive duplication events than others. Primarily, gene essentiality and the type of duplication (small-scale versus large-scale) have been shown in different species to influence the (long-term) survival of novel genes. However, an overarching view of "gene duplicability" is lacking, mainly due to the fact that previous studies usually focused on individual species and did not account for the influence of genomic context and the time of duplication. Here, we present a large-scale study in which we investigated duplicate retention for 9178 gene families shared between 37 flowering plant species, referred to as angiosperm core gene families. For most gene families, we observe a strikingly consistent pattern of gene duplicability across species, with gene families being either primarily single-copy or multicopy in all species. An intermediate class contains gene families that are often retained in duplicate for periods extending to tens of millions of years after whole-genome duplication, but ultimately appear to be largely restored to singleton status, suggesting that these genes may be dosage balance sensitive. The distinction between single-copy and multicopy gene families is reflected in their functional annotation, with single-copy genes being mainly involved in the maintenance of genome stability and organelle function and multicopy genes in signaling, transport, and metabolism. The intermediate class was overrepresented in regulatory genes, further suggesting that these represent putative dosage-balance-sensitive genes.

Paralogue Interference Affects the Dynamics after Gene Duplication

Proteins tend to form homomeric complexes of identical subunits, which act as functional units. By definition, the subunits are encoded from a single genetic locus. When such a gene is duplicated, the gene products are suggested initially to cross-interact when coexpressed, thus resulting in the phenomenon of paralogue interference. In this opinion article, we explore how paralogue interference can shape the fate of a duplicated gene. One important outcome is a prolonged time window in which both copies remain under selection increasing the chance to accumulate mutations and to develop new properties. Thereby, paralogue interference can mediate the coevolution of duplicates and here we illustrate the potential of this phenomenon in light of recent new studies.

Evolutionary history exposes radical diversification among classes of interaction partners of the MLLE domain of plant poly(A)-binding proteins

Background

Poly(A)-binding proteins (PABPs) are evolutionarily conserved proteins that have important functions in the regulation of translation and the control of mRNA stability in eukaryotes. Most PABPs encode a C-terminal domain known as the MLLE domain (previously PABC or CTC), which can mediate protein interactions. In earlier work we identified and predicted that four classes of MLLE-interacting proteins were present in Arabidopsis thaliana, which we named CID A, B, C, and D. These proteins encode transcription-activating domains (CID A), the Lsm and LsmAD domains of ataxin-2 (CID B), the CUE and small MutS-related domains (CID C), and two RNA recognition domains (CID D). We recently found that a novel class that lacks the LsmAD domain is present in CID B proteins.

Results

We extended our analysis to other classes of CIDs present in the viridiplantae. We found that novel variants also evolved in classes CID A and CID C. A specific transcription factor domain is present in a distinct lineage in class A, and a variant that lacks at least two distinct domains was also identified in a divergent lineage in class C. We did not detect any variants in Class D CIDs. This class often consists of four to six highly conserved RNA-binding proteins, which suggests that major redundancy is present in this class.

Conclusions

CIDs are likely to operate as components of posttranscriptional regulatory assemblies. The evident diversification of CIDs may be neutral or may be important for plant adaptation to the environment and for acquisition of specific traits during evolution. The fact that CIDs subclasses are maintained in early lineages suggest that a presumed interference between duplicates was resolved, and a defined function for each subclass was achieved.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0475-1) contains supplementary material, which is available to authorized users.

Interaction between the trout mineralocorticoid and glucocorticoid receptors in vitro

The salmonid corticosteroid receptors (CRs), glucocorticoid receptors 1 and 2 (GR1 and GR2) and the mineralocorticoid receptor (MR) share a high degree of homology with regard to structure, ligand- and DNA response element-binding, and cellular co-localization. Typically, these nuclear hormone receptors homodimerize to confer transcriptional activation of target genes, but a few studies using mammalian receptors suggest some degree of heterodimerization. We observed that the trout MR confers a several fold lower transcriptional activity compared to the trout GRs. This made us question the functional relevance of the MR when this receptor is located in the same cells as the GRs and activated by cortisol. A series of co-transfection experiments using different glucocorticoid response elements (GREs) containing promoter-reporter constructs were carried out to investigate any possible interaction between the piscine CRs. Co-transfection of the GRs with the MR significantly reduced the total transcriptional activity even at low MR levels, suggesting interaction between these receptors. Co-transfection of GR1 or GR2 with the MR did not affect the subcellular localization of the GRs, and the MR-mediated inhibition seemed to be independent of specific activation or inhibition of the MR. Site-directed mutagenesis of the DNA-binding domain and dimerization interface of the MR showed that the inhibition was dependent on DNA binding but not necessarily on dimerization ability. Thus, we suggest that the interaction between MR and the GRs may regulate the cortisol response in cell types where the receptors co-localize and propose a dominant-negative role for the MR in cortisol-mediated transcriptional activity.

Origin of the response to adrenal and sex steroids: Roles of promiscuity and co-evolution of enzymes and steroid receptors

Many responses to adrenal and sex steroids are mediated by receptors that belong to the nuclear receptor family of transcription factors. We investigated the co-evolution of these vertebrate steroid receptors and the enzymes that synthesize adrenal and sex steroids through data mining of genomes from cephalochordates [amphioxus], cyclostomes [lampreys, hagfish], chondrichthyes [sharks, rays, skates], actinopterygii [ray-finned fish], sarcopterygii [coelacanths, lungfishes and terrestrial vertebrates]. An ancestor of the estrogen receptor and 3-ketosteroid receptors evolved in amphioxus. A corticoid receptor and a progesterone receptor evolved in cyclostomes, and an androgen receptor evolved in gnathostomes. Amphioxus contains CYP11, CYP17, CYP19, 3β/Δ5-4-HSD and 17β-HSD14, which suffice for the synthesis of estradiol and Δ5-androstenediol. Amphioxus also contains CYP27, which catalyzes the synthesis of 27-hydroxy-cholesterol, another estrogen. Lamprey contains, in addition, CYP21, which catalyzes the synthesis of 11-deoxycortisol. Chondrichthyes contain, in addition, CYP11A, CYP11C, CYP17A1, CYP17A2. Coelacanth also contains CYP11C1, the current descendent from a common ancestor with modern land vertebrate CYP11B genes, which catalyze the synthesis of cortisol, corticosterone and aldosterone. Interestingly, CYP11B2, aldosterone synthase, evolved from separate gene duplications in at least old world monkeys and two suborders of rodents. Sciurognathi (including mice and rats) and Hystricomorpha (including guinea pigs). Thus, steroid receptors and steroidogenic enzymes co-evolved at key transitions in the evolution of vertebrates. Together, this suite of receptors and enzymes through their roles in transcriptional regulation of reproduction, development, homeostasis and the response to stress contributed to the evolutionary diversification of vertebrates. This article is part of a Special Issue entitled 'Steroid/Sterol signaling'.

Intermolecular epistasis shaped the function and evolution of an ancient transcription factor and its DNA binding sites

Complexes of specifically interacting molecules, such as transcription factor proteins (TFs) and the DNA response elements (REs) they recognize, control most biological processes, but little is known concerning the functional and evolutionary effects of epistatic interactions across molecular interfaces. We experimentally characterized all combinations of genotypes in the joint protein-DNA sequence space defined by an historical transition in TF-RE specificity that occurred some 500 million years ago in the DNA-binding domain of an ancient steroid hormone receptor. We found that rampant epistasis within and between the two molecules was essential to specific TF-RE recognition and to the evolution of a novel TF-RE complex with unique derived specificity. Permissive and restrictive epistatic mutations across the TF-RE interface opened and closed potential evolutionary paths accessible by the other, making the evolution of each molecule contingent on its partner's history and allowing a molecular complex with novel specificity to evolve.

The nuclear receptors of Biomphalaria glabrata and Lottia gigantea: implications for developing new model organisms

Nuclear receptors (NRs) are transcription regulators involved in an array of diverse physiological functions including key roles in endocrine and metabolic function. The aim of this study was to identify nuclear receptors in the fully sequenced genome of the gastropod snail, Biomphalaria glabrata, intermediate host for Schistosoma mansoni and compare these to known vertebrate NRs, with a view to assessing the snail's potential as a invertebrate model organism for endocrine function, both as a prospective new test organism and to elucidate the fundamental genetic and mechanistic causes of disease. For comparative purposes, the genome of a second gastropod, the owl limpet, Lottia gigantea was also investigated for nuclear receptors. Thirty-nine and thirty-three putative NRs were identified from the B. glabrata and L. gigantea genomes respectively, based on the presence of a conserved DNA-binding domain and/or ligand-binding domain. Nuclear receptor transcript expression was confirmed and sequences were subjected to a comparative phylogenetic analysis, which demonstrated that these molluscs have representatives of all the major NR subfamilies (1-6). Many of the identified NRs are conserved between vertebrates and invertebrates, however differences exist, most notably, the absence of receptors of Group 3C, which includes some of the vertebrate endocrine hormone targets. The mollusc genomes also contain NR homologues that are present in insects and nematodes but not in vertebrates, such as Group 1J (HR48/DAF12/HR96). The identification of many shared receptors between humans and molluscs indicates the potential for molluscs as model organisms; however the absence of several steroid hormone receptors indicates snail endocrine systems are fundamentally different.