Vestigialization of an allosteric switch: genetic and structural mechanisms for the evolution of constitutive activity in a steroid hormone receptor

Vertebrate endocrine disruptors induce sex-reversal in blue mussels

Mollusks are the second most diverse animal phylum, yet little is known about their endocrinology or how they respond to endocrine disrupting compound (EDC) pollution. Characteristic effects of endocrine disruption are reproductive impairment, skewed sex ratios, development of opposite sex characteristics, and population decline. However, whether classical vertebrate EDCs, such as steroid hormone-like chemicals and inhibitors of steroidogenesis, exert effects on mollusks is controversial. In the blue mussel, Mytilus edulis, EDC exposure is correlated with feminized sex ratios in wild and laboratory mussels, but sex reversal has not been confirmed. Here, we describe a non-destructive qPCR assay to identify the sex of M. edulis allowing identification of males and females prior to experimentation. We exposed male mussels to 17α-ethinylestradiol and female mussels to ketoconazole, EDCs that mimic vertebrate steroid hormones or inhibit their biosynthesis. Both chemicals changed the sex of individual mussels, interfered with gonadal development, and disrupted gene expression of the sex differentiation pathway. Impacts from ketoconazole treatment, including changes in steroid levels, confirmed a role for steroidogenesis and steroid-like hormones in mollusk endocrinology. The present study expands the possibilities for laboratory and field monitoring of mollusk species and provides key insights into endocrine disruption and sexual differentiation in bivalves.

Identification of an ionic mechanism for ERα-mediated rapid excitation in neurons

The major female ovarian hormone, 17β-estradiol (E2), can alter neuronal excitability within milliseconds to regulate a variety of physiological processes. Estrogen receptor-α (ERα), classically known as a nuclear receptor, exists as a membrane-bound receptor to mediate this rapid action of E2, but the ionic mechanisms remain unclear. Here, we show that a membrane channel protein, chloride intracellular channel protein-1 (Clic1), can physically interact with ERα with a preference to the membrane-bound ERα. Clic1-mediated currents can be enhanced by E2 and reduced by its depletion. In addition, Clic1 currents are required to mediate the E2-induced rapid excitations in multiple brain ERα populations. Further, genetic disruption of Clic1 in hypothalamic ERα neurons blunts the regulations of E2 on female body weight balance. In conclusion, we identified the Clic1 chloride channel as a key mediator for E2-induced rapid neuronal excitation, which may have a broad impact on multiple neurobiological processes regulated by E2.

Estrogen receptor knockdown suggests its role in gonadal development regulation in Manila clam Ruditapes philippinarum

The estrogen receptor (ER), a ligand-dependent transcription factor, is critical for vertebrate reproduction. However, its role in bivalves is not well understood, with ongoing debates regarding its function in regulating reproduction similarly to vertebrates. To investigate ER's function, we conducted a 21-day RNA interference experiment focusing on its role in gonadal development in bivalves. Histological analyses revealed that ER inhibition significantly suppressed ovarian development in females and, conversely, promoted gonadal development in males. Additionally, levels of 17β-estrogen (E2) were markedly reduced in the gonads of both sexes following ER suppression. Transcriptomic analysis from RNA-seq of testes and ovaries after ER interference showed changes in the expression of key genes such as Vtg, CYP17, 3β-HSD, and 17β-HSD. These genes are involved in the estrogen signaling pathway and steroid hormone biosynthesis. Furthermore, ER suppression significantly affected the expression of genes linked to gametogenesis and the reproductive cycle. Our findings highlight ER's crucial, yet complex and sex-specific roles in gonadal development in bivalves, emphasizing the need for further detailed studies.

Is there potential for estradiol receptor signaling in lophotrochozoans?

Estrogen receptors (ERs) are thought to be the ancestor of all steroid receptors and are present in most lophotrochozoans studied to date, including molluscs, annelids, and rotifers. A number of studies have investigated the functional role of estrogen receptors in invertebrate species, although most are in molluscs, where the receptor is constitutively active. In vitro experiments provided evidence for ligand-activated estrogen receptors in annelids, raising important questions about the role of estrogen signalling in lophotrochozoan lineages. Here, we review the concordant and discordant evidence of estradiol receptor signalling in lophotrochozoans, with a focus on annelids and rotifers. We explore the de novo synthesis of estrogens, the evolution and expression of estrogen receptors, and physiological responses to activation of estrogen receptors in the lophotrochozoan phyla Annelida and Rotifera. Key data are missing to determine if de novo biosynthesis of estradiol in non-molluscan lophotrochozoans is likely. For example, an ortholog for the CYP11 gene is present, but confirmation of substrate conversion and measured tissue products is lacking. Orthologs CYP17 and CYP19 are lacking, yet intermediates or products (e.g. estradiol) in tissues have been measured. Estrogen receptors are present in multiple species, and for a limited number, in vitro data show agonist binding of estradiol and/or transcriptional activation. The expression patterns of the lophotrochozoan ERs suggest developmental, reproductive, and digestive roles but are highly species dependent. E2 exposures suggest that lophotrochozoan ERs may play a role in reproduction, but no strong dose-response relationship has been established. Therefore, we expect most lophotrochozoan species, outside of perhaps platyhelminths, to have an ER but their physiological role remains elusive. Mining genomes for orthologs gene families responsible for steroidogenesis, coupled with in vitro and in vivo studies of the steroid pathway are needed to better assess whether lophotrochozoans are capable of estradiol biosynthesis. One major challenge is that much of the data are divided across a diversity of species. We propose that the polychaetes Capitella teleta or Platyneris dumerilii, and rotifer Brachionus manjavacas may be strong species choices for studies of estrogen receptor signalling, because of available genomic data, established laboratory culture techniques, and gene knockout potential.

Ioxynil and diethylstilbestrol impair cardiac performance and shell growth in the mussel Mytilus coruscus

The herbicide ioxynil (IOX) and the synthetic estrogen diethylstilbestrol (DES) are environmentally relevant contaminants that act as endocrine disruptors (EDCs) and have recently been shown to be cardiovascular disruptors in vertebrates. Mussels, Mytilus coruscus, were exposed to low doses of IOX (0.37, 0.037 and 0.0037 mg/L) and DES (0.27, 0.027 and 0.0027 mg/L) via the water and the effect monitored by generating whole animal transcriptomes and measuring cardiac performance and shell growth. One day after IOX (0.37 and 0.037 mg/L) and DES (0.27 and 0.027 mg/L) exposure heart rate frequency was decreased in both groups and 0.27 mg/L DES significantly reduced heart rate frequency with increasing time of exposure (P < 0.05) and no acclimatization occurred. The functional effects were coupled to significant differential expression of genes of the serotonergic synapse pathway and cardiac-related genes at 0.027 mg/L DES, which suggests that impaired heart function may be due to interference with neuroendocrine regulation and direct cardiac effect genes. Multiple genes related to detoxifying xenobiotic substances were up regulated and genes related to immune function were down regulated in the DES group (vs. control), indicating that detoxification processes were enhanced, and the immune response was depressed. In contrast, IOX had a minor disrupting effect at a molecular level. Of note was a significant suppression (P < 0.05) by DES of shell growth in juveniles and lower doses (< 0.0027 mg/L) had a more severe effect. The shell growth depression in 0.0027 mg/L DES-treated juveniles was not accompanied by abundant differential gene expression, suggesting that the effect of 0.0027 mg/L DES on shell growth may be direct. The results obtained in the present study reveal for the first time that IOX and DES may act as neuroendocrine disrupters with a broad spectrum of effects on cardiac performance and shell growth, and that DES exposure had a much more pronounced effect than IOX in a marine bivalve.

Nuclear receptor superfamily structural diversity in pacific oyster: In silico identification of estradiol binding candidates

The increasing presence of anthropogenic contaminants in aquatic environments poses challenges for species inhabiting contaminated sites. Due to their structural binding characteristics to ligands that inhibit or activate gene transcription, these xenobiotic compounds frequently target the nuclear receptor superfamily. The present work aims to understand the potential interaction between the hormone 17-β-estradiol, an environmental contaminant, and the nuclear receptors of Crassostrea gigas, the Pacific oyster. This filter-feeding, sessile oyster species is subject to environmental changes and exposure to contaminants. In the Pacific oyster, the estrogen-binding nuclear receptor is not able to bind this hormone as it does in vertebrates. However, another receptor may exhibit responsiveness to estrogen-like molecules and derivatives. We employed high-performance in silico methodologies, including three-dimensional modeling, molecular docking and atomistic molecular dynamics to identify likely binding candidates with the target moecule. Our approach revealed that among the C. gigas nuclear receptor superfamily, candidates with the most favorable interaction with the molecule of interest belonged to the NR1D, NR1H, NR1P, NR2E, NHR42, and NR0B groups. Interestingly, NR1H and NR0B were associated with planktonic/larval life cycle stages, while NR1P, NR2E, and NR0B were associated with sessile/adult life stages. The application of this computational methodological strategy demonstrated high performance in the virtual screening of candidates for binding with the target xenobiotic molecule and can be employed in other studies in the field of ecotoxicology in non-model organisms.

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.

Topiary: Pruning the manual labor from ancestral sequence reconstruction

Ancestral sequence reconstruction (ASR) is a powerful tool to study the evolution of proteins and thus gain deep insight into the relationships among protein sequence, structure, and function. A major barrier to its broad use is the complexity of the task: it requires multiple software packages, complex file manipulations, and expert phylogenetic knowledge. Here we introduce topiary, a software pipeline that aims to overcome this barrier. To use topiary, users prepare a spreadsheet with a handful of sequences. Topiary then: (1) Infers the taxonomic scope for the ASR study and finds relevant sequences by BLAST; (2) Does taxonomically informed sequence quality control and redundancy reduction; (3) Constructs a multiple sequence alignment; (4) Generates a maximum-likelihood gene tree; (5) Reconciles the gene tree to the species tree; (6) Reconstructs ancestral amino acid sequences; and (7) Determines branch supports. The pipeline returns annotated evolutionary trees, spreadsheets with sequences, and graphical summaries of ancestor quality. This is achieved by integrating modern phylogenetics software (Muscle5, RAxML-NG, GeneRax, and PastML) with online databases (NCBI and the Open Tree of Life). In this paper, we introduce non-expert readers to the steps required for ASR, describe the specific design choices made in topiary, provide a detailed protocol for users, and then validate the pipeline using datasets from a broad collection of protein families. Topiary is freely available for download: https://github.com/harmslab/topiary.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Resurrecting Enzymes by Ancestral Sequence Reconstruction

Ancestral Sequence Reconstruction (ASR) allows one to infer the sequences of extinct proteins using the phylogeny of extant proteins. It consists of disclosing the evolutionary history-i.e., the phylogeny-of a protein family of interest and then inferring the sequences of its ancestors-i.e., the nodes in the phylogeny. Assisted by gene synthesis, the selected ancestors can be resurrected in the lab and experimentally characterized. The crucial step to succeed with ASR is starting from a reliable phylogeny. At the same time, it is of the utmost importance to have a clear idea on the evolutionary history of the family under study and the events that influenced it. This allows us to implement ASR with well-defined hypotheses and to apply the appropriate experimental methods. In the last years, ASR has become popular to test hypotheses about the origin of functionalities, changes in activities, understanding physicochemical properties of proteins, among others. In this context, the aim of this chapter is to present the ASR approach applied to the reconstruction of enzymes-i.e., proteins with catalytic roles. The spirit of this contribution is to provide a basic, hands-to-work guide for biochemists and biologists who are unfamiliar with molecular phylogenetics.

Epistasis and intramolecular networks in protein evolution

Proteins are molecular machines composed of complex, highly connected amino acid networks. Their functional optimization requires the reorganization of these intramolecular networks by evolution. In this review, we discuss the mechanisms by which epistasis, that is, the dependence of the effect of a mutation on the genetic background, rewires intramolecular interactions to alter protein function. Deciphering the biophysical basis of epistasis is crucial to our understanding of evolutionary dynamics and the elucidation of sequence-structure-function relationships. We featured recent studies that provide insights into the molecular mechanisms giving rise to epistasis, particularly at the structural level. These studies illustrate the convoluted and fascinating nature of the intramolecular networks co-opted by epistasis during the evolution of protein function.

Nuclear Receptors and Development of Marine Invertebrates

Nuclear Receptors (NRs) are a superfamily of transcription factors specific to metazoans that have the unique ability to directly translate the message of a signaling molecule into a transcriptional response. In vertebrates, NRs are pivotal players in countless processes of both embryonic and adult physiology, with embryonic development being one of the most dynamic periods of NR activity. Accumulating evidence suggests that NR signaling is also a major regulator of development in marine invertebrates, although ligands and transactivation dynamics are not necessarily conserved with respect to vertebrates. The explosion of genome sequencing projects and the interpretation of the resulting data in a phylogenetic context allowed significant progress toward an understanding of NR superfamily evolution, both in terms of molecular activities and developmental functions. In this context, marine invertebrates have been crucial for characterizing the ancestral states of NR-ligand interactions, further strengthening the importance of these organisms in the field of evolutionary developmental biology.

Estrogen receptor regulates immune defense by suppressing NF-κB signaling in the Crassostrea hongkongensis

The crosstalk between the estrogen receptor (ER) and NF-κB signalling pathways has merged in vertebrates and plays a key role in the control of genes involved in inflammation, cell proliferation and apoptosis. However, such crosstalk between the endocrine and immune systems needs to be explored in lower invertebrates. In this study, we identified a 2856-bp homologue of the estrogen receptor from Hong Kong oyster (ChER), containing a 5' untranslated region (UTR) of 234 bp, a 3' UTR of 387 bp, and an open reading frame (ORF) of 2235 bp. We observed that overexpression of ChER suppressed ChRel-dependent NF-kappaB (NF-κB) activation in the HEK293T (human embryonic kidney 293T) cell line, and depletion of ChER in vivo resulted in upregulation of two NF-κB-responsive marker genes, namely, TNF-α and IL-17, which confirmed its potential role in controlling NF-κB signalling. Furthermore, an EMSA (electrophoretic mobility shift assay) showed that ChER could negatively regulate the binding of ChRel to NF-κB probe-responsive elements. Serial domain requirement analysis showed that both region C (DNA-binding domain) and region E (ligand-binding domain) of ChER were essential for mediating the crosstalk underlying ChER-dependent NF-κB suppression. In conclusion, we demonstrate for the first time the negative regulatory role of the ER in NF-κB signalling in oysters, strongly indicating the presence of complex crosstalk between the endocrine and immune systems in lower marine molluscs.

Effects of Whole-Lake Additions of Ethynylestradiol on Leech Populations

Leeches are widespread, found in many freshwater habitats, and have diverse dietary habits. Despite their close phylogenetic relationships to Mollusca, a phylum with species affected by exogenous estrogens, it is unclear whether Hirudinea may also be impacted. A whole-lake experiment was done at the Experimental Lakes Area in Ontario, Canada, to assess whether 17α-ethynylestradiol (EE2) affected fishes and other species. Herein, we examined whether EE2 impacted leech community composition, species abundance, growth rates, gonad size, and cocoon production, when compared with reference lakes using a before-after-control-impact design. Each month baited leech traps were set overnight in the littoral zone at 10 sites around experimental Lake 260 and 2 reference lakes, and individuals were identified, weighed, and measured. Male and female organs of common species Haemopis marmorata were measured. Across all lakes, 9 species representing 3 families were collected. There were no apparent effects of EE2 on numbers, species richness, or community composition; however, condition in 1 of 5 species increased significantly after EE2 exposures. Total gonadosomatic index (GSI) and the GSI for all male or all female organs combined for H. marmorata were not affected by EE2 additions. However, some individual reproductive organs including relative sperm sac length (+), relative epididymis weight (-), relative vaginal bulb length (+), and relative ovisac + albumen gland length (+) changed after EE2 additions. Finally, whereas overall cocoon production was similar, it occurred earlier in the EE2-amended lake. In summary, few individual through community measures of leeches responded to low ng/L concentrations of EE2, suggesting that Hirudinea are less sensitive to this endocrine disruptor than other invertebrates and vertebrates. Environ Toxicol Chem 2020;39:1608-1619. © 2020 SETAC.

A neural mechanism of nuclear receptor expression and regionalization

Spatially restricted expression of genes by global circulating inducers (hormones, secreted proteins, growth factors, neuromodulators, etc.) was a prerequisite for the evolution of animals. Far from a random occurrence, it is a systematically occurring, certain event, implying that specific information is invested for it to happen. In this minireview, we show for the first time that the expression and regionalization takes place at the level of receptors via a neural mechanism and make an attempt to reconstruct the causal chain from neural signaling to expression of nuclear receptors.

Invertebrates facing environmental contamination by endocrine disruptors: Novel evidences and recent insights

The crisis of biodiversity we currently experience raises the question of the impact of anthropogenic chemicals on wild life health. Endocrine disruptors are notably incriminated because of their possible effects on development and reproduction, including at very low doses. As commonly recorded in the field, the burden they impose on wild species also concerns invertebrates, with possible specificities linked with the specific physiology of these animals. A better understanding of chemically-mediated endocrine disruption in these species has clearly gained from knowledge accumulated on vertebrate models. But the molecular pathways specific to invertebrates also need to be reckoned, which implies dedicated research efforts to decipher their basic functioning in order to be able to assess its possible disruption. The recent rising of omics technologies opens the way to an intensification of these efforts on both aspects, even in species almost uninvestigated so far.

How to resurrect ancestral proteins as proxies for ancient biogeochemistry

Throughout the history of life, enzymes have served as the primary molecular mediators of biogeochemical cycles by catalyzing the metabolic pathways that interact with geochemical substrates. The byproducts of enzymatic activities have been preserved as chemical and isotopic signatures in the geologic record. However, interpretations of these signatures are limited by the assumption that such enzymes have remained functionally conserved over billions of years of molecular evolution. By reconstructing ancient genetic sequences in conjunction with laboratory enzyme resurrection, preserved biogeochemical signatures can instead be related to experimentally constrained, ancestral enzymatic properties. We may thereby investigate instances within molecular evolutionary trajectories potentially tied to significant biogeochemical transitions evidenced in the geologic record. Here, we survey recent enzyme resurrection studies to provide a reasoned assessment of areas of success and common pitfalls relevant to ancient biogeochemical applications. We conclude by considering the Great Oxidation Event, which provides a constructive example of a significant biogeochemical transition that warrants investigation with ancestral enzyme resurrection. This event also serves to highlight the pitfalls of facile interpretation of paleophenotype models and data, as applied to two examples of enzymes that likely both influenced and were influenced by the rise of atmospheric oxygen - RuBisCO and nitrogenase.

Estrogenic compounds as exogenous modulators of physiological functions in molluscs: Signaling pathways and biological responses

Molluscs have been widely utilized to evaluate the effects of estrogenic compounds, one of the most widespread classes of Endocrine Disrupting Chemicals-EDCs. However, knowledge on steroid signaling and metabolism in molluscs has considerably increased in the last decade: from these studies, a considerable debate emerged on the role of 'natural' steroids in physiology, in particular in reproduction, of this invertebrate group. In this work, available information on the effects and mechanisms of action of estrogens in molluscs will be reviewed, with particular emphasis on bivalves that, widespread in aquatic ecosystems, are most likely affected by exposure to estrogenic EDCs. Recent advances in steroid uptake and metabolism, and estrogen receptors-ERs in molluscs, as well as in estrogen signaling in vertebrates, will be considered. The results so far obtained with 17β-estradiol and different estrogenic compounds in the model bivalve Mytilus spp., demonstrate specific effects on immune function, development and metabolism. Transcriptomic data reveal non genomic estrogen signaling pathways in mussel tissues that are supported by new observations at the cellular level. In vitro and in vivo data show, through independent lines of evidence, that estrogens act through non-genomic signaling pathways in bivalves. In this light, regardless of whether molluscs synthesize estrogens de novo or not, and despite their ERs are not directly activated by ligand binding, estrogens can interact with multiple signaling components, leading to modulation of different physiological functions. Increasing knowledge in endocrine physiology of molluscs will provide a framework for a better evaluation and interpretation of data on the impact of estrogenic EDCs in this invertebrate group.

Aromatisation of steroids in the bivalve Mytilus trossulus

In this study, we demonstrated the presence of the enzymatic complex able to perform aromatization (estrogen synthesis) in both, the microsomal and mitochondrial fractions of gills and gonads from Mytilus trossulus. Based on in vitro experiments, we highlighted the importance of temperature as the limiting factor of aromatisation efficiency (AE) in mussels. After testing range of temperatures (4–23 °C), the highest AE was found during incubation at 8 °C and pH 7.6 (41.66 pmol/h/mg protein in gills and 58.37 pmol/h/mg protein in gonads). The results were confirmed during field studies where the most efficient aromatisation occurred in bivalves collected in spring while the least effective in those collected in winter. During in vitro studies, AE turned out to be more intensive in female gonads than in male gonads. The process was also more intensive in mitochondrial fraction than in microsomal one (62.97 pmol/h/mg protein in male gills and 73.94 pmol/h/mg protein in female gonads). Enzymatic complex (aromatase-like enzyme) catalysing aromatisation in mussels was found to be insensitive to inhibitory effect of selective inhibitors of mammalian aromatase such as letrozole and anastrazole, suggesting its different structure from vertebrate aromatase. Further in vivo studies using ¹³C-labeled steroids at 8 °C temperature window confirmed that bivalves are able to uptake testosterone and androstenedione from the ambient environment and metabolise them to estrone and 17β-estradiol thus confirming endogenous estrogen’ synthesis.

Freshwater snail Parafossarulus striatulus estrogen receptor: Characteristics and expression profiles under lab and field exposure

The modes of action by which putative endocrine disrupting chemicals (EDCs) elicit toxicity in mollusks remains unclear due to our limited understanding of the molluscan endocrine system. We identified and partially characterised the estrogen receptor (ER) of the mollusk Parafossarulus striatulus. The full-length cDNA of the ER of P. striatulus (psER) was isolated and found to have an ORF of 1386 bp which corresponded to 461 amino acids. Phylogenetic analysis revealed that psER is an orthologue of ER of other mollusks. Moreover, the DNA-binding domain, ligand-binding domain, P-box, D-box, and AF2 domain were also identified in psER. Exposure of females and males to 17β-estradiol (E2, 100 ng/L) for 24 h and 72 h did not alter psER transcription, but exposure to 17α-methyltestosterone (MT, 100 μg/L) for 72 h significantly decreased ER transcription in females only (p < 0.05). psER transcription was surveyed in males and females seeded in different regions in Taihu Lake, China. psER transcription were elevated among females and males maintained at site ML. This elevation was statistically significant (p < 0.05) among male snails as compared to snails held at the more pristine site of SZ. This was different to the results from lab, implying that some unknown chemicals or other environmental factors in field could affect psER transcription level in snails. Furthermore, females and males held at site ML also exhibited a significant elevation in vitellogenin transcription as compared to snails held at site SZ, suggesting that vitellogenin production may be directly regulated by psER or co-regulated with psER in this species.

NR3E receptors in cnidarians: A new family of steroid receptor relatives extends the possible mechanisms for ligand binding

Steroid hormone receptors are important regulators of development and physiology in bilaterian animals, but the role of steroid signaling in cnidarians has been contentious. Cnidarians produce steroids, including A-ring aromatic steroids with a side-chain, but these are probably made through pathways different than the one used by vertebrates to make their A-ring aromatic steroids. Here we present comparative genomic analyses indicating the presence of a previously undescribed nuclear receptor family within medusozoan cnidarians, that we propose to call NR3E. This family predates the diversification of ERR/ER/SR in bilaterians, indicating that the first NR3 evolved in the common ancestor of the placozoan and cnidarian-bilaterian with lineage-specific loss in the anthozoans, even though multiple species in this lineage have been shown to produce aromatic steroids, whose function remain unclear. We discovered serendipitously that a cytoplasmic factor within epidermal cells of transgenic Hydra vulgaris can trigger the nuclear translocation of heterologously expressed human ERα. This led us to hypothesize that aromatic steroids may also be present in the medusozoan cnidarian lineage, which includes Hydra, and may explain the translocation of human ERα. Docking experiments with paraestrol A, a cnidarian A-ring aromatic steroid, into the ligand-binding pocket of Hydra NR3E indicates that, if an aromatic steroid is indeed the true ligand, which remains to be demonstrated, it would bind to the pocket through a partially distinct mechanism from the manner in which estradiol binds to vertebrate ER.

Transcriptomic and Network Analyses Reveal Mechanistic-Based Biomarkers of Endocrine Disruption in the Marine Mussel, Mytilus edulis

Transcriptomics, high-throughput assays, and adverse outcome pathways (AOP) are promising approaches applied to toxicity monitoring in the 21st century, but development of these methods is challenging for nonmodel organisms and emerging contaminants. For example, Endocrine Disrupting Compounds (EDCs) may cause reproductive impairments and feminization of male bivalves; however, the mechanism linked to this adverse outcome is unknown. To develop mechanism-based biomarkers that may be linked through an AOP, we exposed Mytilus edulis to 17-alpha-ethinylestradiol (5 and 50 ng/L) and 4-nonylphenol (1 and 100 μg/L) for 32 and 39 days. When mussels were exposed to these EDCs, we found elevated female specific transcripts and significant female-skewed sex ratios using a RT-qPCR assay. We performed gene expression analysis on digestive gland tissue using an M. edulis microarray and through network and targeted analyses identified the nongenomic estrogen signaling pathway and steroidogenesis pathway as the likely mechanisms of action for a putative AOP. We also identified several homologues to genes within the vertebrate steroidogenesis pathway including the cholesterol side chain cleavage complex. From this AOP, we designed the Coastal Biosensor for Endocrine Disruption (C-BED) assay which was confirmed in the laboratory and tested in the field.

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.

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.

Conservation of estrogen receptor function in invertebrate reproduction

Background

Rotifers are microscopic aquatic invertebrates that reproduce both sexually and asexually. Though rotifers are phylogenetically distant from humans, and have specialized reproductive physiology, this work identifies a surprising conservation in the control of reproduction between humans and rotifers through the estrogen receptor. Until recently, steroid signaling has been observed in only a few invertebrate taxa and its role in regulating invertebrate reproduction has not been clearly demonstrated. Insights into the evolution of sex signaling pathways can be gained by clarifying how receptors function in invertebrate reproduction.

Results

In this paper, we show that a ligand-activated estrogen-like receptor in rotifers binds human estradiol and regulates reproductive output in females. In other invertebrates characterized thus far, ER ligand binding domains have occluded ligand-binding sites and the ERs are not ligand activated. We have used a suite of computational, biochemical and biological techniques to determine that the rotifer ER binding site is not occluded and can bind human estradiol.

Conclusions

Our results demonstrate that this mammalian hormone receptor plays a key role in reproduction of the ancient microinvertebrate Brachinous manjavacas. The presence and activity of the ER within the phylum Rotifera indicates that the ER structure and function is highly conserved throughout animal evolution.

Electronic supplementary material

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

Robustness of Reconstructed Ancestral Protein Functions to Statistical Uncertainty

Hypotheses about the functions of ancient proteins and the effects of historical mutations on them are often tested using ancestral protein reconstruction (APR)-phylogenetic inference of ancestral sequences followed by synthesis and experimental characterization. Usually, some sequence sites are ambiguously reconstructed, with two or more statistically plausible states. The extent to which the inferred functions and mutational effects are robust to uncertainty about the ancestral sequence has not been studied systematically. To address this issue, we reconstructed ancestral proteins in three domain families that have different functions, architectures, and degrees of uncertainty; we then experimentally characterized the functional robustness of these proteins when uncertainty was incorporated using several approaches, including sampling amino acid states from the posterior distribution at each site and incorporating the alternative amino acid state at every ambiguous site in the sequence into a single "worst plausible case" protein. In every case, qualitative conclusions about the ancestral proteins' functions and the effects of key historical mutations were robust to sequence uncertainty, with similar functions observed even when scores of alternate amino acids were incorporated. There was some variation in quantitative descriptors of function among plausible sequences, suggesting that experimentally characterizing robustness is particularly important when quantitative estimates of ancient biochemical parameters are desired. The worst plausible case method appears to provide an efficient strategy for characterizing the functional robustness of ancestral proteins to large amounts of sequence uncertainty. Sampling from the posterior distribution sometimes produced artifactually nonfunctional proteins for sequences reconstructed with substantial ambiguity.

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.

Impact of bisphenol A (BPA) on early embryo development in the marine mussel Mytilus galloprovincialis: Effects on gene transcription

Bisphenol A (BPA), a monomer used in plastic manufacturing, is weakly estrogenic and a potential endocrine disruptor in mammals. Although it degrades quickly, it is pseudo-persistent in the environment because of continual inputs, with reported concentrations in aquatic environments between 0.0005 and 12 μg/L. BPA represents a potential concern for aquatic ecosystems, as shown by its reproductive and developmental effects in aquatic vertebrates. In invertebrates, endocrine-related effects of BPA were observed in different species and experimental conditions, with often conflicting results, indicating that the sensitivity to this compound can vary considerably among related taxa. In the marine mussel Mytilus galloprovincialis BPA was recently shown to affect early development at environmental concentrations. In this work, the possible effects of BPA on mussel embryos were investigated at the molecular level by evaluating transcription of 13 genes, selected on the basis of their biological functions in adult mussels. Gene expression was first evaluated in trocophorae and D-veligers (24 and 48 h post fertilization) grown in physiological conditions, in comparison with unfertilized eggs. Basal expressions showed a general up-regulation during development, with distinct transcript levels in trocophorae and D-veligers. Exposure of fertilized eggs to BPA (10 μg/L) induced a general upregulation at 24 h pf, followed by down regulation at 48 h pf. Mytilus Estrogen Receptors, serotonin receptor and genes involved in biomineralization (Carbonic Anydrase and Extrapallial Protein) were the most affected by BPA exposure. At 48 h pf, changes in gene expression were associated with irregularities in shell formation, as shown by scanning electron microscopy (SEM), indicating that the formation of the first shelled embryo, a key step in mussel development, represents a sensitive target for BPA. Similar results were obtained with the natural estrogen 17β-estradiol. The results demonstrate that BPA and E₂ can affect Mytilus early development through dysregulation of gene transcription.

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.

Epistasis in protein evolution

The structure, function, and evolution of proteins depend on physical and genetic interactions among amino acids. Recent studies have used new strategies to explore the prevalence, biochemical mechanisms, and evolutionary implications of these interactions-called epistasis-within proteins. Here we describe an emerging picture of pervasive epistasis in which the physical and biological effects of mutations change over the course of evolution in a lineage-specific fashion. Epistasis can restrict the trajectories available to an evolving protein or open new paths to sequences and functions that would otherwise have been inaccessible. We describe two broad classes of epistatic interactions, which arise from different physical mechanisms and have different effects on evolutionary processes. Specific epistasis-in which one mutation influences the phenotypic effect of few other mutations-is caused by direct and indirect physical interactions between mutations, which nonadditively change the protein's physical properties, such as conformation, stability, or affinity for ligands. In contrast, nonspecific epistasis describes mutations that modify the effect of many others; these typically behave additively with respect to the physical properties of a protein but exhibit epistasis because of a nonlinear relationship between the physical properties and their biological effects, such as function or fitness. Both types of interaction are rampant, but specific epistasis has stronger effects on the rate and outcomes of evolution, because it imposes stricter constraints and modulates evolutionary potential more dramatically; it therefore makes evolution more contingent on low-probability historical events and leaves stronger marks on the sequences, structures, and functions of protein families.

Molecular evolution of anthocyanin pigmentation genes following losses of flower color

Background

Phenotypic transitions, such as trait gain or loss, are predicted to carry evolutionary consequences for the genes that control their development. For example, trait losses can result in molecular decay of the pathways underlying the trait. Focusing on the Iochrominae clade (Solanaceae), we examine how repeated losses of floral anthocyanin pigmentation associated with flower color transitions have affected the molecular evolution of three anthocyanin pathway genes (Chi, F3h, and Dfr).

Results

We recovered intact coding regions for the three genes in all of the lineages that have lost floral pigmentation, suggesting that molecular decay is not associated with these flower color transitions. However, two of the three genes (Chi, F3h) show significantly elevated dN/dS ratios in lineages without floral pigmentation. Maximum likelihood analyses suggest that this increase is due to relaxed constraint on anthocyanin genes in the unpigmented lineages as opposed to positive selection. Despite the increase, the values for dN/dS in both pigmented and unpigmented lineages were consistent overall with purifying selection acting on these loci.

Conclusions

The broad conservation of anthocyanin pathway genes across lineages with and without floral anthocyanins is consistent with the growing consensus that losses of pigmentation are largely achieved by changes in gene expression as opposed to structural mutations. Moreover, this conservation maintains the potential for regain of flower color, and indicates that evolutionary losses of floral pigmentation may be readily reversible.

Electronic supplementary material

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

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.

Genetic and structural analyses of cytochrome P450 hydroxylases in sex hormone biosynthesis: Sequential origin and subsequent coevolution

Biosynthesis of steroid hormones in vertebrates involves three cytochrome P450 hydroxylases, CYP11A1, CYP17A1 and CYP19A1, which catalyze sequential steps in steroidogenesis. These enzymes are conserved in the vertebrates, but their origin and existence in other chordate subphyla (Tunicata and Cephalochordata) have not been clearly established. In this study, selected protein sequences of CYP11A1, CYP17A1 and CYP19A1 were compiled and analyzed using multiple sequence alignment and phylogenetic analysis. Our analyses show that cephalochordates have sequences orthologous to vertebrate CYP11A1, CYP17A1 or CYP19A1, and that echinoderms and hemichordates possess CYP11-like but not CYP19 genes. While the cephalochordate sequences have low identity with the vertebrate sequences, reflecting evolutionary distance, the data show apparent origin of CYP11 prior to the evolution of CYP19 and possibly CYP17, thus indicating a sequential origin of these functionally related steroidogenic CYPs. Co-occurrence of the three CYPs in early chordates suggests that the three genes may have coevolved thereafter, and that functional conservation should be reflected in functionally important residues in the proteins. CYP19A1 has the largest number of conserved residues while CYP11A1 sequences are less conserved. Structural analyses of human CYP11A1, CYP17A1 and CYP19A1 show that critical substrate binding site residues are highly conserved in each enzyme family. The results emphasize that the steroidogenic pathways producing glucocorticoids and reproductive steroids are several hundred million years old and that the catalytic structural elements of the enzymes have been conserved over the same period of time. Analysis of these elements may help to identify when precursor functions linked to these enzymes first arose.

Environmental Effects of BPA: Focus on Aquatic Species

Research on bisphenol A (BPA) as an environmental contaminant has now major regulatory implications toward the ecosystem health, and hence it is incumbent on scientists to do their research to the highest standards possible, in order that the most appropriate decisions are made to mitigate the impacts to aquatic wildlife. However, the contribution given so far appears rather fragmented. The present overview aims to collect available information on the effects of BPA on aquatic vertebrates and invertebrates to provide a general scenario and to suggest future developments toward more comprehensive approaches useful for aquatic species protection.

Molecular characterization of an estrogen receptor and estrogen-related receptor and their autoregulatory capabilities in two Mytilus species

Vertebrate-like sex steroid hormones have been widely detected in mollusks, and numerous experiments have shown the importance of steroids in gonad development. Nevertheless, their signaling pathways in invertebrates have not been uncovered yet. Steroid receptors are an ancient class of transcription factors with multiple roles in not only vertebrates but also invertebrates. Estrogen signaling is thought to have major roles in mollusk physiology, but the full repertoire of estrogen receptors is unknown. We presented the successful cloning of two novel forms of estrogen receptor-like genes. These receptors are present in two closely related species of Mytilus: Mytilus edulis and Mytilus galloprovincialis, commonly known and widely distributed sentinel species. Our phylogenetic analysis revealed that one of these receptors is an estrogen receptor (ER) and the other one is an estrogen-related receptor (ERR). Studies of expression analysis showed that both receptor mRNAs were localized in the oocytes and follicle cells in contact with developing oocytes in the ovary and Sertoli cells in the testis, and in the ciliated cells of the gill. In addition, we have evidence that one (ER) of these may have a capacity to autoregulate its own expression in the gonadal cells by estrogen (E2) and that this gene is responsive to estrogenic compounds.

Marine invertebrate xenobiotic-activated nuclear receptors: their application as sensor elements in high-throughput bioassays for marine bioactive compounds

Developing high-throughput assays to screen marine extracts for bioactive compounds presents both conceptual and technical challenges. One major challenge is to develop assays that have well-grounded ecological and evolutionary rationales. In this review we propose that a specific group of ligand-activated transcription factors are particularly well-suited to act as sensors in such bioassays. More specifically, xenobiotic-activated nuclear receptors (XANRs) regulate transcription of genes involved in xenobiotic detoxification. XANR ligand-binding domains (LBDs) may adaptively evolve to bind those bioactive, and potentially toxic, compounds to which organisms are normally exposed to through their specific diets. A brief overview of the function and taxonomic distribution of both vertebrate and invertebrate XANRs is first provided. Proof-of-concept experiments are then described which confirm that a filter-feeding marine invertebrate XANR LBD is activated by marine bioactive compounds. We speculate that increasing access to marine invertebrate genome sequence data, in combination with the expression of functional recombinant marine invertebrate XANR LBDs, will facilitate the generation of high-throughput bioassays/biosensors of widely differing specificities, but all based on activation of XANR LBDs. Such assays may find application in screening marine extracts for bioactive compounds that could act as drug lead compounds.

A mollusk retinoic acid receptor (RAR) ortholog sheds light on the evolution of ligand binding

Nuclear receptors are transcription factors that regulate networks of target genes in response to small molecules. There is a strong bias in our knowledge of these receptors because they were mainly characterized in classical model organisms, mostly vertebrates. Therefore, the evolutionary origins of specific ligand-receptor couples still remain elusive. Here we present the identification and characterization of a retinoic acid receptor (RAR) from the mollusk Nucella lapillus (NlRAR). We show that this receptor specifically binds to DNA response elements organized in direct repeats as a heterodimer with retinoid X receptor. Surprisingly, we also find that NlRAR does not bind all-trans retinoic acid or any other retinoid we tested. Furthermore, NlRAR is unable to activate the transcription of reporter genes in response to stimulation by retinoids and to recruit coactivators in the presence of these compounds. Three-dimensional modeling of the ligand-binding domain of NlRAR reveals an overall structure that is similar to vertebrate RARs. However, in the ligand-binding pocket (LBP) of the mollusk receptor, the alteration of several residues interacting with the ligand has apparently led to an overall decrease in the strength of the interaction with the ligand. Accordingly, mutations of NlRAR at key positions within the LBP generate receptors that are responsive to retinoids. Altogether our data suggest that, in mollusks, RAR has lost its affinity for all-trans retinoic acid, highlighting the evolutionary plasticity of its LBP. When put in an evolutionary context, our results reveal new structural and functional features of nuclear receptors validated by millions of years of evolution that were impossible to reveal in model organisms.