A Mixture of Chemicals Found in Human Amniotic Fluid Disrupts Brain Gene Expression and Behavior in Xenopus laevis

Thyroid hormones (TH) are essential for normal brain development, influencing neural cell differentiation, migration, and synaptogenesis. Multiple endocrine-disrupting chemicals (EDCs) are found in the environment, raising concern for their potential effects on TH signaling and the consequences on neurodevelopment and behavior. While most research on EDCs investigates the effects of individual chemicals, human health may be adversely affected by a mixture of chemicals. The potential consequences of EDC exposure on human health are far-reaching and include problems with immune function, reproductive health, and neurological development. We hypothesized that embryonic exposure to a mixture of chemicals (containing phenols, phthalates, pesticides, heavy metals, and perfluorinated, polychlorinated, and polybrominated compounds) identified as commonly found in the human amniotic fluid could lead to altered brain development. We assessed its effect on TH signaling and neurodevelopment in an amphibian model (Xenopus laevis) highly sensitive to thyroid disruption. Fertilized eggs were exposed for eight days to either TH (thyroxine, T₄ 10 nM) or the amniotic mixture (at the actual concentration) until reaching stage NF47, where we analyzed gene expression in the brains of exposed tadpoles using both RT-qPCR and RNA sequencing. The results indicate that whilst some overlap on TH-dependent genes exists, T₄ and the mixture have different gene signatures. Immunohistochemistry showed increased proliferation in the brains of T₄-treated animals, whereas no difference was observed for the amniotic mixture. Further, we demonstrated diminished tadpoles' motility in response to T₄ and mixture exposure. As the individual chemicals composing the mixture are considered safe, these results highlight the importance of examining the effects of mixtures to improve risk assessment.

Crosstalk between Thyroid Hormone and Corticosteroid Signaling Targets Cell Proliferation in Xenopus tropicalis Tadpole Liver

Thyroid hormones (TH) and glucocorticoids (GC) are involved in numerous developmental and physiological processes. The effects of individual hormones are well documented, but little is known about the joint actions of the two hormones. To decipher the crosstalk between these two hormonal pathways, we conducted a transcriptional analysis of genes regulated by TH, GC, or both hormones together in liver of Xenopus tropicalis tadpoles using RNA-Seq. Among the differentially expressed genes (DE), 70.5% were regulated by TH only, 0.87% by GC only, and 15% by crosstalk between the two hormones. Gene ontology analysis of the crosstalk-regulated genes identified terms referring to DNA replication, DNA repair, and cell-cycle regulation. Biological network analysis identified groups of genes targeted by the hormonal crosstalk and corroborated the gene ontology analysis. Specifically, we found two groups of functionally linked genes (chains) mainly composed of crosstalk-regulated hubs (highly interactive genes), and a large subnetwork centred around the crosstalk-regulated genes psmb6 and cdc7. Most of the genes in the chains are involved in cell-cycle regulation, as are psmb6 and cdc7, which regulate the G2/M transition. Thus, the biological action of these two hormonal pathways acting together in the liver targets cell-cycle regulation.

Marine gregarine genomes reveal the breadth of apicomplexan diversity with a partially conserved glideosome machinery

Our current view of the evolutionary history, coding and adaptive capacities of Apicomplexa, protozoan parasites of a wide range of metazoan, is currently strongly biased toward species infecting humans, as data on early diverging apicomplexan lineages infecting invertebrates is extremely limited. Here, we characterized the genome of the marine eugregarine Porospora gigantea, intestinal parasite of Lobsters, remarkable for the macroscopic size of its vegetative feeding forms (trophozoites) and its gliding speed, the fastest so far recorded for Apicomplexa. Two highly syntenic genomes named A and B were assembled. Similar in size (~ 9 Mb) and coding capacity (~ 5300 genes), A and B genomes are 10.8% divergent at the nucleotide level, corresponding to 16-38 My in divergent time. Orthogroup analysis across 25 (proto)Apicomplexa species, including Gregarina niphandrodes, showed that A and B are highly divergent from all other known apicomplexan species, revealing an unexpected breadth of diversity. Phylogenetically these two species branch sisters to Cephaloidophoroidea, and thus expand the known crustacean gregarine superfamily. The genomes were mined for genes encoding proteins necessary for gliding, a key feature of apicomplexans parasites, currently studied through the molecular model called glideosome. Sequence analysis shows that actin-related proteins and regulatory factors are strongly conserved within apicomplexans. In contrast, the predicted protein sequences of core glideosome proteins and adhesion proteins are highly variable among apicomplexan lineages, especially in gregarines. These results confirm the importance of studying gregarines to widen our biological and evolutionary view of apicomplexan species diversity, and to deepen our understanding of the molecular bases of key functions such as gliding, well known to allow access to the intracellular parasitic lifestyle in Apicomplexa.

Adhesive gland transcriptomics uncovers a diversity of genes involved in glue formation in marine tube-building polychaetes

Tube-building sabellariid polychaetes are hermatypic organisms capable of forming vast reefs in highly turbulent marine habitats. Sabellariid worms assemble their tube by gluing together siliceous and calcareous clastic particles using a polyelectrolytic biocement. Here, we performed transcriptomic analyses to investigate the genes that are differentially expressed in the parathorax region, which contains the adhesive gland and tissues, from the rest of the body. We found a large number of candidate genes to be involved in the composition and formation of biocement in two species: Sabellaria alveolata and Phragmatopoma caudata. Our results indicate that the glue is likely to be composed by a large diversity of cement-related proteins, including Poly(S), GY-rich, H-repeat and miscellaneous categories. However, sequences divergence and differences in expression profiles between S. alveolata and P. caudata of cement-related proteins may reflect adaptation to the type of substratum used to build their tube, and/or to their habitat (temperate vs tropical, amplitude of pH, salinity …). Related to the L-DOPA metabolic pathways and linked with the genes that were differentially expressed in the parathorax region, we found that tyrosinase and peroxidase gene families may have undergone independent expansion in the two Sabellariidae species investigated. Our data also reinforce the importance of protein modifications in cement formation. Altogether these new genomic resources help to identify novel transcripts encoding for cement-related proteins, but also important enzymes putatively involved in the chemistry of the adhesion process, such as kinases, and may correspond to new targets to develop biomimetic approaches.

STATEMENTS OF SIGNIFICANCE

The diversity of bioadhesives elaborated by marine invertebrates is a tremendous source of inspiration to develop biomimetic approaches for biomedical and technical applications. Recent studies on the adhesion system of mussel, barnacle and sea star had highlighted the usefulness of high-throughput RNA sequencing in accelerating the development of biomimetic adhesives. Adhesion in sandcastle worms, which involves catechol and phosphate chemistries, polyelectrolyte complexes, supramolecular architectures, and a coacervation process, is a useful model to develop multipurpose wet adhesives. Using transcriptomic tools, we have explored the diversity of genes encoding for structural and catalytic proteins involved in cement formation of two sandcastle worm species, Sabellaria alveolata and Phragmatopoma caudata. The important genomic resource generated should help to design novel "blue" adhesives.

Comparative analysis of molecular signatures suggests the use of gabapentin for the management of endometriosis-associated pain

Background

It has been repetitively shown that the transcription factors DLX5 and DLX6 are drastically downregulated in endometriotic lesions when compared with eutopic endometrium. These findings suggest that regulatory cascades involving DLX5/6 might be at the origin of endometriosis symptoms such as chronic pelvic pain (CPP). We have shown that inactivation of Dlx5 and Dlx5/6 in the mouse uterus results in an endometrial phenotype reminiscent of endometriosis.

Methods

We focused on genes that present a similar deregulation in endometriosis and in Dlx5/6-null mice in search of new endometriosis targets.

Results

We confirmed a strong reduction of DLX5 expression in endometriosis implants. We identified a signature of 30 genes similarly deregulated in human endometriosis implants and in Dlx5/6-null mouse uteri, reinforcing the notion that the downregulation of Dlx5/6 is an early event in the progress of endometriosis. CACNA2D3, a component of the α2δ family of voltage-dependent calcium channel complex, was strongly overexpressed both in mutant mouse uteri and in endometriosis implants, were also CACNA2D1 and CACNA2D2, other members of the α2δ family involved in nociception, are upregulated.

Conclusion

Comparative analysis of gene expression signatures from endometriosis and mouse models showed that calcium channel subunits α2δ involved in nociception can be targets for the treatment of endometriosis-associated pain. CACNA2D3 has been associated with pain sensitization and heat nociception in animal models. In patients, CACNA2D3 variants were associated with reduced sensitivity to acute noxious stimuli. As α2δs were targets of gabapentinoid analgesics, the results suggested the use of these drugs for the treatment of endometriosis-associated pain. Indeed, recent small-scale clinical studies have shown that gabapentin could be effective in women with CPP. The findings of this study reinforce the need for a large definitive trial.

Insights from the Shell Proteome: Biomineralization to Adaptation

Bivalves have evolved a range of complex shell forming mechanisms that are reflected by their incredible diversity in shell mineralogy and microstructures. A suite of proteins exported to the shell matrix space plays a significant role in controlling these features, in addition to underpinning some of the physical properties of the shell itself. Although, there is a general consensus that a minimum basic protein tool kit is required for shell construction, to date, this remains undefined. In this study, the shell matrix proteins (SMPs) of four highly divergent bivalves (The Pacific oyster, Crassostrea gigas; the blue mussel, Mytilus edulis; the clam, Mya truncata, and the king scallop, Pecten maximus) were analyzed in an identical fashion using proteomics pipeline. This enabled us to identify the critical elements of a "basic tool kit" for calcification processes, which were conserved across the taxa irrespective of the shell morphology and arrangement of the crystal surfaces. In addition, protein domains controlling the crystal layers specific to aragonite and calcite were also identified. Intriguingly, a significant number of the identified SMPs contained domains related to immune functions. These were often are unique to each species implying their involvement not only in immunity, but also environmental adaptation. This suggests that the SMPs are selectively exported in a complex mix to endow the shell with both mechanical protection and biochemical defense.

An integrated omic analysis of hepatic alteration in medaka fish chronically exposed to cyanotoxins with possible mechanisms of reproductive toxicity

Cyanobacterial blooms threaten human health as well as the population of other living organisms in the aquatic environment, particularly due to the production of natural toxic components, the cyanotoxin. So far, the most studied cyanotoxins are microcystins (MCs). In this study, the hepatic alterations at histological, proteome and transcriptome levels were evaluated in female and male medaka fish chronically exposed to 1 and 5 μg L^-1 microcystin-LR (MC-LR) and to the extract of MC-producing Microcystis aeruginosa PCC 7820 (5 μg L^-1 of equivalent MC-LR) by balneation for 28 days, aiming at enhancing our understanding of the potential reproductive toxicity of cyanotoxins in aquatic vertebrate models. Indeed, both MC and Microcystis extract adversely affect reproductive parameters including fecundity and egg hatchability. The liver of toxin treated female fish present glycogen storage loss and cellular damages. The quantitative proteomics analysis revealed that the quantities of 225 hepatic proteins are dysregulated. In particular, a notable decrease in protein quantities of vitellogenin and choriogenin was observed, which could explain the decrease in reproductive output. Liver transcriptome analysis through Illumina RNA-seq reveals that over 100-400 genes are differentially expressed under 5 μg L^-1 MC-LR and Microcystis extract treatments, respectively. Ingenuity pathway analysis of the omic data attests that various metabolic pathways, such as energy production, protein biosynthesis and lipid metabolism, are disturbed by both MC-LR and the Microcystis extract, which could provoke the observed reproductive impairment. The transcriptomics analysis also constitutes the first report of the impairment of circadian rhythm-related gene induced by MCs. This study contributes to a better understanding of the potential consequences of chronic exposure of fish to environmental concentrations of cyanotoxins, suggesting that Microcystis extract could impact a wider range of biological pathways, compared with pure MC-LR, and even 1 μg L^-1 MC-LR potentially induces a health risk for aquatic organisms.

Deep sexual dimorphism in adult medaka fish liver highlighted by multi-omic approach

Sexual dimorphism describes the features that discriminate between the two sexes at various biological levels. Especially, during the reproductive phase, the liver is one of the most sexually dimorphic organs, because of different metabolic demands between the two sexes. The liver is a key organ that plays fundamental roles in various physiological processes, including digestion, energetic metabolism, xenobiotic detoxification, biosynthesis of serum proteins, and also in endocrine or immune response. The sex-dimorphism of the liver is particularly obvious in oviparous animals, as the female liver is the main organ for the synthesis of oocyte constituents. In this work, we are interested in identifying molecular sexual dimorphism in the liver of adult medaka fish and their sex-variation in response to hepatotoxic exposures. By developing an integrative approach combining histology and different high-throughput omic investigations (metabolomics, proteomics and transcriptomics), we were able to globally depict the strong sexual dimorphism that concerns various cellular and molecular processes of hepatocytes comprising protein synthesis, amino acid, lipid and polysaccharide metabolism, along with steroidogenesis and detoxification. The results of this work imply noticeable repercussions on the biology of oviparous organisms environmentally exposed to chemical or toxin issues.

Improved Genome Editing Efficiency and Flexibility Using Modified Oligonucleotides with TALEN and CRISPR-Cas9 Nucleases

Genome editing has now been reported in many systems using TALEN and CRISPR-Cas9 nucleases. Precise mutations can be introduced during homology-directed repair with donor DNA carrying the wanted sequence edit, but efficiency is usually lower than for gene knockout and optimal strategies have not been extensively investigated. Here, we show that using phosphorothioate-modified oligonucleotides strongly enhances genome editing efficiency of single-stranded oligonucleotide donors in cultured cells. In addition, it provides better design flexibility, allowing insertions more than 100 bp long. Despite previous reports of phosphorothioate-modified oligonucleotide toxicity, clones of edited cells are readily isolated and targeted sequence insertions are achieved in rats and mice with very high frequency, allowing for homozygous loxP site insertion at the mouse ROSA locus in particular. Finally, when detected, imprecise knockin events exhibit indels that are asymmetrically positioned, consistent with genome editing taking place by two steps of single-strand annealing.

Dlx5 and Dlx6 control uterine adenogenesis during post-natal maturation: possible consequences for endometriosis

Dlx5 and Dlx6 are two closely associated homeobox genes which code for transcription factors involved in the control of steroidogenesis and reproduction. Inactivation of Dlx5/6 in the mouse results in a Leydig cell defect in the male and in ovarian insufficiency in the female. DLX5/6 are also strongly expressed by the human endometrium but their function in the uterus is unknown. The involvement of DLX5/6 in human uterine pathology is suggested by their strong downregulation in endometriotic lesions and upregulation in endometrioïd adenocarcinomas. We first show that Dlx5/6 expression begins in Müllerian ducts epithelia and persists then in the uterine luminal and glandular epithelia throughout post-natal maturation and in the adult. We then use a new mouse model in which Dlx5 and Dlx6 can be simultaneously inactivated in the endometrium using a Pgr(cre/+) allele. Post-natal inactivation of Dlx5/6 in the uterus results in sterility without any obvious ovarian involvement. The uteri of Pgr(cre/+); Dlx5/6(flox/flox) mice present very few uterine glands and numerous abnormally large and branched invaginations of the uterine lumen. In Dlx5/6 mutant uteri, the expression of genes involved in gland formation (Foxa2) and in epithelial remodelling during implantation (Msx1) is significantly reduced. Furthermore, we show that DLX5 is highly expressed in human endometrial glandular epithelium and that its expression is affected in endometriosis. We conclude that Dlx5 and Dlx6 expression determines uterine architecture and adenogenesis and is needed for implantation. Given their importance for female reproduction, DLX5 and DLX6 must be regarded as interesting targets for future clinical research.

Genome-wide negative feedback drives transgenerational DNA methylation dynamics in Arabidopsis

Epigenetic variations of phenotypes, especially those associated with DNA methylation, are often inherited over multiple generations in plants. The active and inactive chromatin states are heritable and can be maintained or even be amplified by positive feedback in a transgenerational manner. However, mechanisms controlling the transgenerational DNA methylation dynamics are largely unknown. As an approach to understand the transgenerational dynamics, we examined long-term effect of impaired DNA methylation in Arabidopsis mutants of the chromatin remodeler gene DDM1 (Decrease in DNA Methylation 1) through whole genome DNA methylation sequencing. The ddm1 mutation induces a drastic decrease in DNA methylation of transposable elements (TEs) and repeats in the initial generation, while also inducing ectopic DNA methylation at hundreds of loci. Unexpectedly, this ectopic methylation can only be seen after repeated self-pollination. The ectopic cytosine methylation is found primarily in the non-CG context and starts from 3’ regions within transcription units and spreads upstream. Remarkably, when chromosomes with reduced DNA methylation were introduced from a ddm1 mutant into a DDM1 wild-type background, the ddm1-derived chromosomes also induced analogous de novo accumulation of DNA methylation in trans. These results lead us to propose a model to explain the transgenerational DNA methylation redistribution by genome-wide negative feedback. The global negative feedback, together with local positive feedback, would ensure robust and balanced differentiation of chromatin states within the genome.

DNA methylation is important for controlling activity of transposable elements and genes. An intriguing feature of DNA methylation in plants is that its pattern can be inherited over multiple generations at high fidelity in a Mendelian manner. However, mechanisms controlling the trans-generational DNA methylation dynamics are largely unknown. Arabidopsis mutants of a chromatin remodeler gene DDM1 (Decrease in DNA Methylation 1) show drastic reduction of DNA methylation in transposons and repeats, and also show progressive changes in developmental phenotypes during propagation through self-pollination. We now show using whole genome DNA methylation sequencing that upon repeated selfing, the ddm1 mutation induces an ectopic accumulation of DNA methylation at hundreds of loci. Remarkably, even in the wild type background, the analogous de novo increase of DNA methylation can be induced in trans by chromosomes with reduced DNA methylation. Collectively, our findings support a model to explain the transgenerational DNA methylation redistribution by genome-wide negative feedback, which should be important for balanced differentiation of DNA methylation states within the genome.

ncPRO-seq: a tool for annotation and profiling of ncRNAs in sRNA-seq data

SUMMARY

Non-coding RNA (ncRNA) PROfiling in small RNA (sRNA)-seq (ncPRO-seq) is a stand-alone, comprehensive and flexible ncRNA analysis pipeline. It can interrogate and perform detailed profiling analysis on sRNAs derived from annotated non-coding regions in miRBase, Rfam and RepeatMasker, as well as specific regions defined by users. The ncPRO-seq pipeline performs both gene-based and family-based analyses of sRNAs. It also has a module to identify regions significantly enriched with short reads, which cannot be classified under known ncRNA families, thus enabling the discovery of previously unknown ncRNA- or small interfering RNA (siRNA)-producing regions. The ncPRO-seq pipeline supports input read sequences in fastq, fasta and color space format, as well as alignment results in BAM format, meaning that sRNA raw data from the three current major platforms (Roche-454, Illumina-Solexa and Life technologies-SOLiD) can be analyzed with this pipeline. The ncPRO-seq pipeline can be used to analyze read and alignment data, based on any sequenced genome, including mammals and plants.

AVAILABILITY

Source code, annotation files, manual and online version are available at http://ncpro.curie.fr/.

CONTACT

bioinfo.ncproseq@curie.fr or cciaudo@ethz.ch

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.