Generational effects and abnormalities in craniofacial chondrogenesis in zebrafish (Danio rerio) embryos upon maternal exposure to estrogen endocrine disrupting chemicals

Bisphenol A (BPA) and diethyl phthalate (DEP) are estrogenic endocrine disrupting chemicals (EEDCs). The present study reconfirmed that the angle of the ceratohyal cartilage (CH) in embryos were larger from maternal BPA and E2, but smaller from DEP compared to the control. However, it is still unknown whether both the BPA and DEP chemicals disrupted the action of E2 and thereby influence the estrogen signaling pathways. Additionally, it remains unclear whether they also disrupted certain related genes in the migratory pathways of neural crest cells (NCCs) in their offspring. The present data showed that nuclear estrogen receptors and membrane estrogen receptors have different disrupted profiles among female zebrafish exposed to BPA (F-BPA), and DEP (F-DEP), and external E2 (F-E2). However, certain related genes in the migratory pathways of NCCs in embryos from F-BPA and F-E2 such as the sox10, chm1, and tgfbr1a mRNA expressions showed a positive relationship compared with CH angles; the gene expressions of sox9a, smad3, and col2a1a and the CH angles of embryos exhibited an opposite relationship upon F-DEP treatments. Thus, we suggested that the genes involved in NCCs migration were potentially induced by the residual maternal DEP contents. Two sets of genes, chm1/tgfb3 and chm1/gper1, exhibited an identical profile in the ovary and its offspring at 2 h of post fertilization upon F-E2 and F-BPA treatments, respectively. We suggested that the maternal mRNA from female to embryos were transferred before the maternal-to-zygotic transition stage.

The Role of Estrogen and Thyroid Hormones in Zebrafish Visual System Function

Visual system development is a highly complex process involving coordination of environmental cues, cell pathways, and integration of functional circuits. Consequently, a change to any step, due to a mutation or chemical exposure, can lead to deleterious consequences. One class of chemicals known to have both overt and subtle effects on the visual system is endocrine disrupting compounds (EDCs). EDCs are environmental contaminants which alter hormonal signaling by either preventing compound synthesis or binding to postsynaptic receptors. Interestingly, recent work has identified neuronal and sensory systems, particularly vision, as targets for EDCs. In particular, estrogenic and thyroidogenic signaling have been identified as critical modulators of proper visual system development and function. Here, we summarize and review this work, from our lab and others, focusing on behavioral, physiological, and molecular data collected in zebrafish. We also discuss different exposure regimes used, including long-lasting effects of developmental exposure. Overall, zebrafish are a model of choice to examine the impact of EDCs and other compounds targeting estrogen and thyroid signaling and the consequences of exposure in visual system development and function.

The Effects of 17α-Ethinylestradiol on the Heart Rate of Embryonic Japanese Medaka (Oryzias latipes)

Estrogen toxicity has been an area of priority in aquatic toxicology over the last 20 yr. Currently, estrogen toxicity is primarily linked to classical estrogen signaling, the interaction of estrogen receptors alpha and beta (ERα and ERβ). Recent evidence has indicated that a rapid, nongenomic, nonclassical estrogen signaling pathway exists via the G protein-coupled estrogen receptor (GPER), which is expressed in many biological systems, with roles in the cardiovascular system. The objective of the present study was to investigate the effect of 17α-ethinylestradiol (EE2) on the heart rate of embryonic Japanese medaka (Oryzias latipes). A significant decrease (bradycardia) in embryonic heart rate was observed at all treatment concentrations (0.1, 1, 10, 100, and 1000 ng/L EE2) at 144, 168, and 192 h postfertilization (hpf; p ≤ 0.05), whereas 120 and 216 hpf embryos experienced a significant decrease from the control at 10, 100, and 1000 ng/L EE2 and 0.1, 100, and 1000 ng/L EE2, respectively (p ≤ 0.05). In addition, using select estrogen receptor modulators, it was demonstrated that estrogen-induced bradycardia appears to be linked to GPER and not ERα and ERβ. The present study highlights GPER as a novel and alternative mode of action for EE2 toxicity at environmentally relevant concentrations. Environ Toxicol Chem 2020;39:904-912. © 2020 SETAC.

Transcriptome profiling of the fungus Aspergillus nidulans exposed to a commercial glyphosate-based herbicide under conditions of apparent herbicide tolerance

Glyphosate-based herbicides, such as Roundup®, are the most widely used non-selective, broad-spectrum herbicides. The release of these compounds in large amounts into the environment is susceptible to affect soil quality and health, especially because of the non-target effects on a large range of organisms including soil microorganisms. The soil filamentous fungus Aspergillus nidulans, a well-characterized experimental model organism that can be used as a bio-indicator for agricultural soil health, has been previously shown to be highly affected by Roundup GT Plus (R450: 450 g/L of glyphosate) at concentrations far below recommended agricultural application rate, including at a dose that does not cause any macroscopic effect. In this study, we determined alterations in the transcriptome of A. nidulans when exposed to R450 at a dose corresponding to the no-observed-adverse-effect level (NOAEL) for macroscopic parameters. A total of 1816 distinct genes had their expression altered. The most affected biological functions were protein synthesis, amino acids and secondary metabolisms, stress response, as well as detoxification pathways through cytochromes P450, glutathione-S-transferases, and ABC transporters. These results partly explain the molecular mechanisms underlying alterations in growth parameters detected at higher concentrations for this ascomycete fungus. In conclusion, our results highlight molecular disturbances in a soil fungus under conditions of apparent tolerance to the herbicide, and thus confirm the need to question the principle of "substantial equivalence" when applied to plants made tolerant to herbicides.

Estrogen Activation of G-Protein-Coupled Estrogen Receptor 1 Regulates Phosphoinositide 3-Kinase and mTOR Signaling to Promote Liver Growth in Zebrafish and Proliferation of Human Hepatocytes

BACKGROUND & AIMS

Patients with cirrhosis are at high risk for hepatocellular carcinoma (HCC) and often have increased serum levels of estrogen. It is not clear how estrogen promotes hepatic growth. We investigated the effects of estrogen on hepatocyte proliferation during zebrafish development, liver regeneration, and carcinogenesis. We also studied human hepatocytes and liver tissues.

METHODS

Zebrafish were exposed to selective modifiers of estrogen signaling at larval and adult stages. Liver growth was assessed by gene expression, fluorescent imaging, and histologic analyses. We monitored liver regeneration after hepatocyte ablation and HCC development after administration of chemical carcinogens (dimethylbenzanthrazene). Proliferation of human hepatocytes was measured in a coculture system. We measured levels of G-protein-coupled estrogen receptor (GPER1) in HCC and nontumor liver tissues from 68 patients by immunohistochemistry.

RESULTS

Exposure to 17β-estradiol (E2) increased proliferation of hepatocytes and liver volume and mass in larval and adult zebrafish. Chemical genetic and epistasis experiments showed that GPER1 mediates the effects of E2 via the phosphoinositide 3-kinase-protein kinase B-mechanistic target of rapamycin pathway: gper1-knockout and mtor-knockout zebrafish did not increase liver growth in response to E2. HCC samples from patients had increased levels of GPER1 compared with nontumor tissue samples; estrogen promoted proliferation of human primary hepatocytes. Estrogen accelerated hepatocarcinogenesis specifically in male zebrafish. Chemical inhibition or genetic loss of GPER1 significantly reduced tumor development in the zebrafish.

CONCLUSIONS

In an analysis of zebrafish and human liver cells and tissues, we found GPER1 to be a hepatic estrogen sensor that regulates liver growth during development, regeneration, and tumorigenesis. Inhibitors of GPER1 might be developed for liver cancer prevention or treatment.

TRANSCRIPT PROFILING

The accession number in the Gene Expression Omnibus is GSE92544.

Early phosphoproteomic changes for adverse outcome pathway development in the fathead minnow (Pimephales promelas) brain

Adverse outcome pathways (AOPs) are conceptual frameworks that organize and link contaminant-induced mechanistic molecular changes to adverse biological responses at the individual and population level. AOPs leverage molecular and high content mechanistic information for regulatory decision-making, but most current AOPs for hormonally active agents (HAAs) focus on nuclear receptor-mediated effects only despite the overwhelming evidence that HAAs also activate membrane receptors. Activation of membrane receptors triggers non-genomic signaling cascades often transduced by protein phosphorylation leading to phenotypic changes. We utilized label-free LC-MS/MS to identify proteins differentially phosphorylated in the brain of fathead minnows (Pimephales promelas) aqueously exposed for 30 minutes to two HAAs, 17α-ethinylestradiol (EE2), a strong estrogenic substance, and levonorgestrel (LNG), a progestin, both components of the birth control pill. EE2 promoted differential phosphorylation of proteins involved in neuronal processes such as nervous system development, synaptic transmission, and neuroprotection, while LNG induced differential phosphorylation of proteins involved in axon cargo transport and calcium ion homeostasis. EE2 and LNG caused similar enrichment of synaptic plasticity and neurogenesis. This study is the first to identify molecular changes in vivo in fish after short-term exposure and highlights transduction of rapid signaling mechanisms as targets of HAAs, in addition to nuclear receptor-mediated pathways.

The comprehensive expression analysis of the G protein-coupled receptor from Penaeus monodon indicating it participates in innate immunity and anti-ammonia nitrogen stress

The G protein-coupled receptors (GPCRs) composed a superfamily that played an important role in physiological processes of crustaceans, with multiple functions such as growth and development, acting as a defense against stimulations from external factors. In this paper, one kind of GPCRs were identified from Penaeus monodon, called PmGPCR, included an open reading frame (ORF) of 1113 bp. Bioinformatic analysis showed that PmGPCR protein had the typical structure of seven transmembrane domains (7TM), especially the special Asp-Arg-Try motif (DRY motif) between the third transmembrane structures (TM3) and the second intracellular loops (IL-2) which can prove that PmGPCR belongs to the rhodopsin-like family. The analyses of phylogenetic tree indicated that the amino acid sequence of PmGPCR should be merged into Procambarus clarkiic with high identity (98%). Quantitative real-time PCR (q RT-PCR) revealed that PmGPCR mRNA was highly expressed in hepatopancreas, abdominal ganglia and lymph, in which it was significantly higher than that of other tissues (P < 0.05). In addition, the expression of PmGPCR was analyzed during three days post-stimulation with the gram-positive/negative bacteria, the mRNA expression level increased after challenged with gram - positive bacteria in hepatopancreas, lymph and intestines. During the development stages, PmGPCR showed significantly higher expression in nauplius, zoea III, mysis III and post larvae stages than that in other development stages. Meanwhile, the highest transcripts expression of PmGPCR in abdominal ganglia, hepatopancreas, lymph and intestines respectively appeared at D0, D1, D2 and D3/D4 stages of molting. High or low concentration of ammonia nitrogen up-regulated the expression level of PmGPCR at the initial stage in hepatopancreas and gill, and then down-regulated at 48 h. These results indicated PmGPCR may mediate the pathways that involved in growth and development process, survival in the adversity, in addition, provided the useful data to research GPCR-mediated physiological and biological process and explain the mechanisms to defense pathogens and anti-stress in shrimp.

Steroid Transport, Local Synthesis, and Signaling within the Brain: Roles in Neurogenesis, Neuroprotection, and Sexual Behaviors

Sex steroid hormones are synthesized from cholesterol and exert pleiotropic effects notably in the central nervous system. Pioneering studies from Baulieu and colleagues have suggested that steroids are also locally-synthesized in the brain. Such steroids, called neurosteroids, can rapidly modulate neuronal excitability and functions, brain plasticity, and behavior. Accumulating data obtained on a wide variety of species demonstrate that neurosteroidogenesis is an evolutionary conserved feature across fish, birds, and mammals. In this review, we will first document neurosteroidogenesis and steroid signaling for estrogens, progestagens, and androgens in the brain of teleost fish, birds, and mammals. We will next consider the effects of sex steroids in homeostatic and regenerative neurogenesis, in neuroprotection, and in sexual behaviors. In a last part, we will discuss the transport of steroids and lipoproteins from the periphery within the brain (and vice-versa) and document their effects on the blood-brain barrier (BBB) permeability and on neuroprotection. We will emphasize the potential interaction between lipoproteins and sex steroids, addressing the beneficial effects of steroids and lipoproteins, particularly HDL-cholesterol, against the breakdown of the BBB reported to occur during brain ischemic stroke. We will consequently highlight the potential anti-inflammatory, anti-oxidant, and neuroprotective properties of sex steroid and lipoproteins, these latest improving cholesterol and steroid ester transport within the brain after insults.

Fish as a model to assess chemical toxicity in bone

Environmental toxicology has been expanding as growing concerns on the impact of produced and released chemical compounds over the environment and human health are being demonstrated. Among the toxic effects observed in organisms exposed to pollutants, those affecting skeletal tissues (osteotoxicity) have been somehow overlooked in comparison to hepato-, immune-, neuro- and/or reproductive toxicities. Nevertheless, sub-lethal effects of toxicants on skeletal development and/or bone maintenance may result in impaired growth, reduced survival rate, increased disease susceptibility and diminished welfare. Osteotoxicity may occur by acute or chronic exposure to different environmental insults. Because of biologically and technically advantagous features - easy to breed and inexpensive to maintain, external and rapid rate of development, translucent larvae and the availability of molecular and genetic tools - the zebrafish (Danio rerio) has emerged in the last decade as a vertebrate model system of choice to evaluate osteotoxicity. Different experimental approaches in fish species and analytical tools have been applied, from in vitro to in vivo systems, from specific to high throughput methodologies. Current knowledge on osteotoxicity and underlying mechanisms gained using fish, with a special emphasis on zebrafish systems, is reviewed here. Osteotoxicants have been classified into four categories according to the pathway involved in the transduction of the osteotoxic effects: activation/inhibition of membrane and/or nuclear receptors, alteration of redox condition, mimicking of bone constituents and unknown pathways. Knowledge on these pathways is also reported here as it may provide critical insights into the development, production and release of future chemical compounds with none or low osteotoxicity, thus promoting the green/environmental friendly chemistry.

G protein-coupled estrogen receptor regulates embryonic heart rate in zebrafish

Estrogens act by binding to estrogen receptors alpha and beta (ERα, ERβ), ligand-dependent transcription factors that play crucial roles in sex differentiation, tumor growth and cardiovascular physiology. Estrogens also activate the G protein-coupled estrogen receptor (GPER), however the function of GPER in vivo is less well understood. Here we find that GPER is required for normal heart rate in zebrafish embryos. Acute exposure to estrogens increased heart rate in wildtype and in ERα and ERβ mutant embryos but not in GPER mutants. GPER mutant embryos exhibited reduced basal heart rate, while heart rate was normal in ERα and ERβ mutants. We detected gper transcript in discrete regions of the brain and pituitary but not in the heart, suggesting that GPER acts centrally to regulate heart rate. In the pituitary, we observed gper expression in cells that regulate levels of thyroid hormone triiodothyronine (T3), a hormone known to increase heart rate. Compared to wild type, GPER mutants had reduced levels of T3 and estrogens, suggesting pituitary abnormalities. Exposure to exogenous T3, but not estradiol, rescued the reduced heart rate phenotype in gper mutant embryos, demonstrating that T3 acts downstream of GPER to regulate heart rate. Using genetic and mass spectrometry approaches, we find that GPER regulates maternal estrogen levels, which are required for normal embryonic heart rate. Our results demonstrate that estradiol plays a previously unappreciated role in the acute modulation of heart rate during zebrafish embryonic development and suggest that GPER regulates embryonic heart rate by altering maternal estrogen levels and embryonic T3 levels.

Estrogen hormones are important for the formation and function of the nervous, reproductive and cardiovascular systems. Here we report that acute exposure to estrogens increases heart rate, a previously unappreciated function of estrogens. Using zebrafish with mutations in genes that respond to estrogens, we found that heart rate is regulated not by the typical molecules that respond to estrogens–the nuclear estrogen receptors–but rather by a different molecule, the G protein-coupled estrogen receptor. We also show that estrogens increase heart rate by increasing levels of thyroid hormone. Our results reveal a new function for the G protein-coupled estrogen receptor and a new connection between estrogens and thyroid hormone. Environmental compounds that mimic estrogens can be harmful because they can influence gonad function. Our results suggest that endocrine disrupting compounds may also influence cardiac function.

Contribution of G protein-coupled estrogen receptor 1 (GPER) to 17β-estradiol-induced developmental toxicity in zebrafish

Exposure to 17β-estradiol (E2) influences the regulation of multiple signaling pathways, and E2-mediated disruption of signaling events during early development can lead to malformations such as cardiac defects. In this study, we investigated the potential role of the G-protein estrogen receptor 1 (GPER) in E2-induced developmental toxicity. Zebrafish embryos were exposed to E2 from 2h post fertilization (hpf) to 76 hpf with subsequent transcriptional measurements of heart and neural crest derivatives expressed 2 (hand2), leucine rich repeat containing 10 (lrrc10), and gper at 12, 28 and 76 hpf. Alteration in the expression of lrrc10, hand2 and gper was observed at 12 hpf and 76 hpf, but not at 28 hpf. Expression of these genes was also altered after exposure to G1 (a GPER agonist) at 76 hpf. Expression of lrrc10, hand2 and gper all coincided with the formation of cardiac edema at 76 hpf as well as other developmental abnormalities. While co-exposure of G1 with G36 (a GPER antagonist) rescued G1-induced abnormalities and altered gene expression, co-exposure of E2 with G36, or ICI 182,780 (an estrogen receptor antagonist) did not rescue E2-induced cardiac deformities or gene expression. In addition, no effects on the concentrations of downstream ER and GPER signaling molecules (cAMP or calcium) were observed in embryo homogenates after E2 treatment. These data suggest that the impacts of E2 on embryonic development at this stage are complex and may involve multiple receptor and/or signaling pathways.

Comparative Developmental Toxicity of Flavonoids Using an Integrative Zebrafish System

Flavonoids are a large, structurally diverse class of bioactive naturally occurring chemicals commonly detected in breast milk, soy based infant formulas, amniotic fluid, and fetal cord blood. The potential for pervasive early life stage exposures raises concerns for perturbation of embryogenesis, though developmental toxicity and bioactivity information is limited for many flavonoids. Therefore, we evaluated a suite of 24 flavonoid and flavonoid-like chemicals using a zebrafish embryo-larval toxicity bioassay-an alternative model for investigating developmental toxicity of environmentally relevant chemicals. Embryos were exposed to 1-50 µM of each chemical from 6 to 120 h postfertilization (hpf), and assessed for 26 adverse developmental endpoints at 24, 72, and 120 hpf. Behavioral changes were evaluated in morphologically normal animals at 24 and 72 hpf, at 120 hpf using a larval photomotor response (LPR) assay. Gene expression was comparatively evaluated for all compounds for effects on biomarker transcripts indicative of AHR (cyp1a) and ER (cyp19a1b, esr1, lhb, vtg) pathway bioactivity. Overall, 15 of 24 flavonoids elicited adverse effects on one or more of the developmental or behavioral endpoints. Hierarchical clustering and principle component analyses compared toxicity profiles and identified 3 distinct groups of bioactive flavonoids. Despite robust induction of multiple estrogen-responsive biomarkers, co-exposure with ER and GPER antagonists did not ameliorate toxicity, suggesting ER-independence and alternative modes of action. Taken together, these studies demonstrate that development is sensitive to perturbation by bioactive flavonoids in zebrafish that are not related to traditional estrogen receptor mode of action pathways. This integrative zebrafish platform provides a useful framework for evaluating flavonoid developmental toxicity and hazard prioritization.

Suberoylanilide Hydroxamic Acid Restores Estrogen Reduced-cTnI Expression in Neonatal Hearts of Mice

Diastolic cardiac dysfunction can be caused by abnormality in cTnI expression during cardiogenesis. In this study, we investigated the effects of estrogen on the abnormal expression of cTnI in the hearts of neonatal mice and its potential epigenetic mechanisms. We then evaluated suberoylanilide hydroxamic acid (SAHA), a HDAC inhibitor, as a new target treatment of diastolic cardiac dysfunction. Postnatal day 0.5 C57BL/6 mice were injected with estrogen for 1 week, then the hearts of 7-day-old neonatal mice were retrieved for examination. The activities of HDAC and HAT were assayed by colorimetry, and the interaction of cTnI with HDAC5 in mice hearts were examined using chromatin immunoprecipitation assays. The expression of cTnI was tested by quantitative real-time RT-PCR and Western blot. Estrogen treated groups displayed a significantly increased HDAC activity in the hearts of neonatal mice while HAT activity remained unchanged. Additionally, HDAC5 was higher at the cTnI promoter, as compared to the saline treated control groups. The acetylation of histone H3K9ac on cTnI promoter significantly decreased in the hearts of neonatal mice treated with estrogen, and the expression of cTnI at transcriptional and protein levels also decreased. SAHA was shown to increase the acetylation of histone H3K9ac and upregulate the expression of cTnI. The data demonstrated that SAHA can correct cTnI expression abnormality caused by estrogen through inhibiting the binding of HDAC5 to the promoter of cTnI. J. Cell. Biochem. 117: 2377-2384, 2016. © 2016 Wiley Periodicals, Inc.

Phenotype anchoring in zebrafish reveals a potential role for matrix metalloproteinases (MMPs) in tamoxifen's effects on skin epithelium

The zebrafish is a powerful alternative model used to link phenotypes with molecular effects to discover drug mode of action. Using a zebrafish embryo-larval toxicity bioassay, we evaluated the effects of tamoxifen--a widely used anti-estrogen chemotherapeutic. Zebrafish exposed to ≥ 10 μM tamoxifen exhibited a unique necrotic caudal fin phenotype that was rapidly induced regardless of developmental life-stage when treatment was applied. To define tamoxifen's bioactivity resulting in this phenotype, targeted gene expression was used to evaluate 100 transcripts involved in tissue remodeling, calcium signaling, cell cycle and cell death, growth factors, angiogenesis and hypoxia. The most robustly misregulated transcripts in the tail were matrix metalloproteinases mmp9 and mmp13a, induced 127 and 1145 fold, respectively. Expression of c-fos, c-jun, and ap1s1 were also moderately elevated (3-7 fold), consistent with AP-1 activity--a transcription factor that regulates MMP expression. Immunohistochemistry confirmed high levels of induction for MMP13a in affected caudal fin skin epithelial tissue. The necrotic caudal fin phenotype was significantly attenuated or prevented by three functionally unique MMP inhibitors: EDTA (metal chelator), GM 6001 (broad MMP inhibitor), and SR 11302 (AP-1 transcription factor inhibitor), suggesting MMP-dependence. SR 11302 also inhibited induction of mmp9, mmp13a, and a putative MMP target, igfbp1a. Overall, our studies suggest that tamoxifen's effect is the result of perturbation of the MMP system in the skin leading to ectopic expression, cytotoxicity, and the necrotic caudal fin phenotype. These studies help advance our understanding of tamoxifen's non-classical mode of action and implicate a possible role for MMPs in tissues such as skin.

Triphenyl phosphate-induced developmental toxicity in zebrafish: potential role of the retinoic acid receptor

Using zebrafish as a model, we previously reported that developmental exposure to triphenyl phosphate (TPP) - a high-production volume organophosphate-based flame retardant - results in dioxin-like cardiac looping impairments that are independent of the aryl hydrocarbon receptor. Using a pharmacologic approach, the objective of this study was to investigate the potential role of retinoic acid receptor (RAR) - a nuclear receptor that regulates vertebrate heart morphogenesis - in mediating TPP-induced developmental toxicity in zebrafish. We first revealed that static exposure of zebrafish from 5-72h post-fertilization (hpf) to TPP in the presence of non-toxic concentrations of an RAR antagonist (BMS493) significantly enhanced TPP-induced toxicity (relative to TPP alone), even though identical non-toxic BMS493 concentrations mitigated retinoic acid (RA)-induced toxicity. BMS493-mediated enhancement of TPP toxicity was not a result of differential TPP uptake or metabolism, as internal embryonic doses of TPP and diphenyl phosphate (DPP) - a primary TPP metabolite - were not different in the presence or absence of BMS493. Using real-time PCR, we then quantified the relative change in expression of cytochrome P450 26a1 (cyp26a1) - a major target gene for RA-induced RAR activation in zebrafish - and found that RA and TPP exposure resulted in a ∼5-fold increase and decrease in cyp26a1 expression, respectively, relative to vehicle-exposed embryos. To address whether TPP may interact with human RARs, we then exposed Chinese hamster ovary cells stably transfected with chimeric human RARα-, RARβ-, or RARγ to TPP in the presence of RA, and found that TPP significantly inhibited RA-induced luciferase activity in a concentration-dependent manner. Overall, our findings suggest that zebrafish RARs may be involved in mediating TPP-induced developmental toxicity, a mechanism of action that may have relevance to humans.

Aberrant gene expression in humans

Gene expression as an intermediate molecular phenotype has been a focus of research interest. In particular, studies of expression quantitative trait loci (eQTL) have offered promise for understanding gene regulation through the discovery of genetic variants that explain variation in gene expression levels. Existing eQTL methods are designed for assessing the effects of common variants, but not rare variants. Here, we address the problem by establishing a novel analytical framework for evaluating the effects of rare or private variants on gene expression. Our method starts from the identification of outlier individuals that show markedly different gene expression from the majority of a population, and then reveals the contributions of private SNPs to the aberrant gene expression in these outliers. Using population-scale mRNA sequencing data, we identify outlier individuals using a multivariate approach. We find that outlier individuals are more readily detected with respect to gene sets that include genes involved in cellular regulation and signal transduction, and less likely to be detected with respect to the gene sets with genes involved in metabolic pathways and other fundamental molecular functions. Analysis of polymorphic data suggests that private SNPs of outlier individuals are enriched in the enhancer and promoter regions of corresponding aberrantly-expressed genes, suggesting a specific regulatory role of private SNPs, while the commonly-occurring regulatory genetic variants (i.e., eQTL SNPs) show little evidence of involvement. Additional data suggest that non-genetic factors may also underlie aberrant gene expression. Taken together, our findings advance a novel viewpoint relevant to situations wherein common eQTLs fail to predict gene expression when heritable, rare inter-individual variation exists. The analytical framework we describe, taking into consideration the reality of differential phenotypic robustness, may be valuable for investigating complex traits and conditions.

The uniqueness of individuals is due to differences in the combination of genetic, epigenetic and environmental determinants. Understanding the genetic basis of phenotypic variation is a key objective in genetics. Gene expression has been considered as an intermediate phenotype, and the association between gene expression and commonly-occurring genetic variants in the general population has been convincingly established. However, there are few methods to assess the impact of rare genetic variants, such as private SNPs, on gene expression. Here we describe a systematic approach, based on the theory of multivariate outlier detection, to identify individuals that show unusual or aberrant gene expression, relative the rest of the study cohort. Through characterizing detected outliers and corresponding gene sets, we are able to identify which gene sets tend to be aberrantly expressed and which individuals show deviant gene expression within a population. One of our major findings is that private SNPs may contribute to aberrant expression in outlier individuals. These private SNPs are more frequently located in the enhancer and promoter regions of genes that are aberrantly expressed, suggesting a possible regulatory function of these SNPs. Overall, our results provide new insight into the determinants of inter-individual variation, which have not been evaluated by large population-level cohort studies.

Early expression of aromatase and the membrane estrogen receptor GPER in neuromasts reveals a role for estrogens in the development of the frog lateral line system

Estrogens and their receptors are present at very early stages of vertebrate embryogenesis before gonadal tissues are formed. However, the cellular source and the function of estrogens in embryogenesis remain major questions in developmental endocrinology. We demonstrate the presence of estrogen-synthesizing enzyme aromatase and G protein-coupled estrogen receptor (GPER) proteins throughout early embryogenesis in the model organism, Silurana tropicalis. We provide the first evidence of aromatase in the vertebrate lateral line. High levels of aromatase were detected in the mantle cells of neuromasts, the mechanosensory units of the lateral line, which persisted throughout the course of development (Nieuwkoop and Faber stages 34-47). We show that GPER is expressed in both the accessory and hair cells. Pharmacological activation of GPER with the agonist G-1 disrupted neuromast development and migration. Future study of this novel estrogen system in the amphibian lateral line may shed light on similar systems such as the mammalian inner ear.

Modulation of estrogen causes disruption of craniofacial chondrogenesis in Danio rerio

Estrogen is a steroid hormone that is ubiquitous in vertebrates, but its role in cartilage formation has not been extensively studied. Abnormalities of craniofacial cartilage and bone account for a large portion of birth defects in the United States. Zebrafish (Danio rerio) have been used as models of human disease, and their transparency in the embryonic period affords additional advantages in studying craniofacial development. In this study, zebrafish embryos were treated with 17-β estradiol (E2) or with an aromatase inhibitor and observed for defects in craniofacial cartilage. Concentrations of E2 greater than 2μM caused major disruptions in cartilage formation. Concentrations below 2μM caused subtle changed in cartilage morphology that were only revealed by measurement. The angles formed by cartilage elements in fish treated with 1.5 and 2μM E2 were increasingly wide, while the length of the primary anterior-posterior cartilage element in these fish decreased significantly from controls. These treatments resulted in fish with shorter, flatter faces as estrogen concentration increased. Inhibition of aromatase activity also resulted in similar craniofacial disruption indicating that careful control of estrogen signaling is required for appropriate development. Further investigation of the phenomena described in this study could lead to a better understanding of the etiology of craniofacial birth defects and endocrine disruption of cartilage formation.