Transcriptome and DNA methylome analyses reveal underlying mechanisms for the racial disparity in uterine fibroids.

Uterine fibroids (leiomyomas) affect Black women disproportionately compared with women of other races and ethnicities in terms of prevalence, incidence, and severity of symptoms. The causes of this racial disparity are essentially unknown. We hypothesized that myometria of Black women are more susceptible to developing fibroids, and we examined the transcriptomic and DNA methylation profiles of myometria and fibroids from Black and White women for comparison. Myometrial samples cluster by race in both their transcriptome and DNA methylation profiles, whereas fibroid samples only cluster by race in the latter. More differentially expressed genes (DEGs) were detected in the Black and White myometrial sample comparison than in the fibroid comparison. Leiomyoma gene set expression analysis identified 4 clusters of DEGs, including a cluster of 24 genes with higher expression in myometrial samples from Black women. One of the DEGs in this group, von Willibrands factor (VWF), was significantly hypomethylated in both myometrial samples from Black women and in all fibroids at 2 CpG probes that are near a putative enhancer site and that are correlated with VWF expression levels. These results suggest that the molecular basis for the disparity in fibroid disease between Black and White women could be found in the myometria before fibroid development and not in the fibroids themselves.

Changes in nucleus accumbens gene expression accompany sex-specific suppression of spontaneous physical activity in aromatase knockout mice.

Aromatase catalyzes conversion of testosterone to estradiol and is expressed in a variety of tissues, including the brain. Suppression of aromatase adversely affects metabolism and physical activity behavior, but mechanisms remain uncertain. The hypothesis tested herein was that whole body aromatase deletion would cause gene expression changes in the nucleus accumbens (NAc), a brain regulating motivated behaviors such as physical activity, which is suppressed with loss of estradiol. Metabolic and behavioral assessments were performed in male and female wild-type (WT) and aromatase knockout (ArKO) mice. NAc-specific differentially expressed genes (DEGs) were identified with RNAseq, and associations between the measured phenotypic traits were determined. Female ArKO mice had greater percent body fat, reduced spontaneous physical activity (SPA), consumed less energy, and had lower relative resting energy expenditure (REE) than WT females. Such differences were not observed in ArKO males. However, in both sexes, a top DEG was Pts, a gene encoding an enzyme necessary for catecholamine (e.g., dopamine) biosynthesis. In comparing male and female WT mice, top DEGs were related to sexual development/fertility, immune regulation, obesity, dopamine signaling, and circadian regulation. SPA correlated strongly with Per3, a gene regulating circadian function, thermoregulation, and metabolism (r = -0.64, P = .002), which also correlated with adiposity (r = 0.54, P = .01). In conclusion, aromatase ablation leads to gene expression changes in NAc, which may in turn result in reduced SPA and related metabolic abnormalities. These findings may have significance to post-menopausal women and those treated with an aromatase inhibitor.

Sexual dimorphism in brain transcriptomes of Amami spiny rats (Tokudaia osimensis): a rodent species where males lack the Y chromosome.

BACKGROUND: Brain sexual differentiation is sculpted by precise coordination of steroid hormones during development. Programming of several brain regions in males depends upon aromatase conversion of testosterone to estrogen. However, it is not clear the direct contribution that Y chromosome associated genes, especially sex-determining region Y (Sry), might exert on brain sexual differentiation in therian mammals. Two species of spiny rats: Amami spiny rat (Tokudaia osimensis) and Tokunoshima spiny rat (T. tokunoshimensis) lack a Y chromosome/Sry, and these individuals possess an XO chromosome system in both sexes. Both Tokudaia species are highly endangered. To assess the neural transcriptome profile in male and female Amami spiny rats, RNA was isolated from brain samples of adult male and female spiny rats that had died accidentally and used for RNAseq analyses. RESULTS: RNAseq analyses confirmed that several genes and individual transcripts were differentially expressed between males and females. In males, seminal vesicle secretory protein 5 (Svs5) and cytochrome P450 1B1 (Cyp1b1) genes were significantly elevated compared to females, whereas serine (or cysteine) peptidase inhibitor, clade A, member 3 N (Serpina3n) was upregulated in females. Many individual transcripts elevated in males included those encoding for zinc finger proteins, e.g. zinc finger protein X-linked (Zfx). CONCLUSIONS: This method successfully identified several genes and transcripts that showed expression differences in the brain of adult male and female Amami spiny rat. The functional significance of these findings, especially differential expression of transcripts encoding zinc finger proteins, in this unusual rodent species remains to be determined.

The novel cyst nematode effector protein 30D08 targets host nuclear functions to alter gene expression in feeding sites.

Cyst nematodes deliver effector proteins into host cells to manipulate cellular processes and establish a metabolically hyperactive feeding site. The novel 30D08 effector protein is produced in the dorsal gland of parasitic juveniles, but its function has remained unknown. We demonstrate that expression of 30D08 contributes to nematode parasitism, the protein is packaged into secretory granules and it is targeted to the plant nucleus where it interacts with SMU2 (homolog of suppressor of mec-8 and unc-52 2), an auxiliary spliceosomal protein. We show that SMU2 is expressed in feeding sites and an smu2 mutant is less susceptible to nematode infection. In Arabidopsis expressing 30D08 under the SMU2 promoter, several genes were found to be alternatively spliced and the most abundant functional classes represented among differentially expressed genes were involved in RNA processing, transcription and binding, as well as in development, and hormone and secondary metabolism, representing key cellular processes known to be important for feeding site formation. In conclusion, we demonstrated that the 30D08 effector is secreted from the nematode and targeted to the plant nucleus where its interaction with a host auxiliary spliceosomal protein may alter the pre-mRNA splicing and expression of a subset of genes important for feeding site formation.

Hypothalamic transcriptomic alterations in male and female California mice (Peromyscus californicus) developmentally exposed to bisphenol A or ethinyl estradiol.

Bisphenol A (BPA) is an endocrine-disrupting chemical (EDC) prevalent in many household items. Rodent models and human epidemiological studies have linked this chemical to neurobehavior impairments. In California mice, developmental exposure to BPA results in sociosexual disorders at adulthood, including communication and biparental care deficits, behaviors that are primarily regulated by the hypothalamus. Thus, we sought to examine the transcriptomic profile in this brain region of juvenile male and female California mice offspring exposed from periconception through lactation to BPA or ethinyl estradiol (EE, estrogen present in birth control pills and considered a positive estrogen control for BPA studies). Two weeks prior to breeding, P0 females were fed a control diet, or this diet supplemented with 50 mg BPA/kg feed weight or 0.1 ppb EE, and continued on the diets through lactation. At weaning, brains from male and female offspring were collected, hypothalamic RNA isolated, and RNA-seq analysis performed. Results indicate that BPA and EE groups clustered separately from controls with BPA and EE exposure leading to unique set of signature gene profiles. Kcnd3 was downregulated in the hypothalamus of BPA- and EE-exposed females, whereas Tbl2, Topors, Kif3a, and Phactr2 were upregulated in these groups. Comparison of transcripts differentially expressed in BPA and EE groups revealed significant enrichment of gene ontology terms associated with microtubule-based processes. Current results show that perinatal exposure to BPA or EE can result in several transcriptomic alterations, including those associated with microtubule functions, in the hypothalamus of California mice. It remains to be determined whether these genes mediate BPA-induced behavioral disruptions.

Effects of exposure to bisphenol A and ethinyl estradiol on the gut microbiota of parents and their offspring in a rodent model.

Gut dysbiosis may result in various diseases, such as metabolic and neurobehavioral disorders. Exposure to endocrine disrupting chemicals (EDCs), including bisphenol A (BPA) and ethinyl estradiol (EE), especially during development, may also increase the risk for such disorders. An unexplored possibility is that EDC-exposure might alter the gut microbial composition. Gut flora and their products may thus be mediating factors for the disease-causing effects of these chemicals. To examine the effects of EDCs on the gut microbiome, female and male monogamous and biparental California mice (Peromyscus californicus) were exposed to BPA (50 mg/kg feed weight) or EE (0.1 ppb) or control diet from periconception through weaning. 16s rRNA sequencing was performed on bacterial DNA isolated from fecal samples, and analyses performed for P0 and F1 males and females. Both BPA and EE induced generational and sex-dependent gut microbiome changes. Many of the bacteria, e.g. Bacteroides, Mollicutes, Prevotellaceae, Erysipelotrichaceae, Akkermansia, Methanobrevibacter, Sutterella, whose proportions increase with exposure to BPA or EE in the P0 or F1 generation are associated with different disorders, such as inflammatory bowel disease (IBD), metabolic disorders, and colorectal cancer. However, the proportion of the beneficial bacterium, Bifidobacterium, was also elevated in fecal samples of BPA- and EE-exposed F1 females. Intestinal flora alterations were also linked to changes in various metabolic and other pathways. Thus, BPA and EE exposure may disrupt the normal gut flora, which may in turn result in systemic effects. Probiotic supplementation might be an effective means to mitigate disease-promoting effects of these chemicals.

Intermittent severe hypoxia induces plasticity within serotonergic and catecholaminergic neurons in the neonatal rat ventrolateral medulla.

5-HT neurons contribute to autoresuscitation and survival during intermittent severe hypoxia (IsH). In adults, catecholaminergic neurons in the ventrolateral medulla (VLM) contribute to the autonomic response to hypoxia. We hypothesized that 1) catecholaminergic neurons in the neonatal VLM are activated following IsH, 2) this activation is compromised following an acute loss of brain stem 5-HT, and 3) IsH induces cellular and/or transcriptomic plasticity within catecholaminergic and serotonergic neurons that are within or project to the VLM, respectively. To test these hypotheses, we treated rat pups with 6-fluorotryptophan, a tryptophan hydroxylase (TPH) inhibitor, and then exposed treated and vehicle controls to IsH or air. Along with immunohistochemistry to detect tyrosine hydroxylase (TH)- or Fos-positive neurons, we used RNA sequencing to resolve the effects of IsH and 5-HT deficiency on the expression of serotonergic and catecholaminergic system genes in the VLM. 5-HT deficiency compromised autoresuscitation and survival. IsH significantly increased the number of identifiable TH-positive VLM neurons, an effect enhanced by 5-HT deficiency (P = 0.003). Contrary to our hypothesis, 5-HT-deficient pups had significantly more Fos-positive neurons following IsH (P = 0.008) and more activated TH-positive neurons following IsH or air (P = 0.04). In both groups the expression of the 5-HT transporter and TPH2 was increased following IsH. In 5-HT-deficient pups, the expression of the inhibitory 5-HT1A receptor was decreased following IsH, while the expression of DOPA decarboxylase was increased. These data show that the serotonergic and catecholaminergic systems in the VLM of the neonatal rat are dynamically upregulated by IsH, potentially adapting cardiorespiratory responses to severe hypoxia.

Discovery of a Novel Seminal Fluid Microbiome and Influence of Estrogen Receptor Alpha Genetic Status.

Bacteria harbored in the male reproductive system may influence reproductive function and health of the male and result in developmental origins of adult health and disease (DOHaD) effects in his offspring. Such effects could be due to the seminal fluid, which is slightly basic and enriched with carbohydrates; thereby, creating an ideal habitat for microbes or a potential seminal fluid microbiome (SFM). Using wild-type (WT) and estrogen receptor-alpha (ESR1) knockout (KO) male mice, we describe a unique SFM whose inhabitants differ from gut microbes. The bacterial composition of the SFM is influenced according to whether mice have functional Esr1 genes. Propionibacterium acnes, causative agent of chronic prostatitis possibly culminating in prostate cancer, is reduced in SFM of ESR1 KO compared to WT mice (P ≤ 0.0007). In certain genetic backgrounds, WT mice show a greater incidence of prostate cancer than ESR1 KO, which may be due to increased abundance of P. acnes. Additionally, select gut microbiome residents in ESR1 KO males, such as Lachnospiraceae and Christensenellaceae, might contribute to previously identified phenotypes, especially obesity, in these mutant mice. Understanding how genetics and environmental factors influence the SFM may provide the next frontier in male reproductive disorders and possibly paternal-based DOHaD diseases.

Analysis of genome sequences from plant pathogenic Rhodococcus reveals genetic novelties in virulence loci.

Members of Gram-positive Actinobacteria cause economically important diseases to plants. Within the Rhodococcus genus, some members can cause growth deformities and persist as pathogens on a wide range of host plants. The current model predicts that phytopathogenic isolates require a cluster of three loci present on a linear plasmid, with the fas operon central to virulence. The Fas proteins synthesize, modify, and activate a mixture of growth regulating cytokinins, which cause a hormonal imbalance in plants, resulting in abnormal growth. We sequenced and compared the genomes of 20 isolates of Rhodococcus to gain insights into the mechanisms and evolution of virulence in these bacteria. Horizontal gene transfer was identified as critical but limited in the scale of virulence evolution, as few loci are conserved and exclusive to phytopathogenic isolates. Although the fas operon is present in most phytopathogenic isolates, it is absent from phytopathogenic isolate A21d2. Instead, this isolate has a horizontally acquired gene chimera that encodes a novel fusion protein with isopentyltransferase and phosphoribohydrolase domains, predicted to be capable of catalyzing and activating cytokinins, respectively. Cytokinin profiling of the archetypal D188 isolate revealed only one activate cytokinin type that was specifically synthesized in a fas-dependent manner. These results suggest that only the isopentenyladenine cytokinin type is synthesized and necessary for Rhodococcus phytopathogenicity, which is not consistent with the extant model stating that a mixture of cytokinins is necessary for Rhodococcus to cause leafy gall symptoms. In all, data indicate that only four horizontally acquired functions are sufficient to confer the trait of phytopathogenicity to members of the genetically diverse clade of Rhodococcus.

Characterization of the rat developmental liver transcriptome.

Gene regulation and transcriptome studies have been enabled by the development of RNA-Seq applications for high-throughput sequencing platforms. Next generation sequencing is remarkably efficient and avoids many issues inherent in hybridization-based microarray methodologies including the exon-specific dependence of probe design. Biologically relevant transcripts including messenger and regulatory RNAs may now be quantified and annotated regardless of whether they have previously been observed. We used RNA-Seq to investigate global patterns of gene expression in early developing rat liver. Liver samples from timed-pregnant Lewis rats were collected at six fetal and neonatal stages [embryonic day (E)14, E16, E18, E20, postnatal day (P)1, P7], transcripts were sequenced using an Illumina HiSeq 2000, and data analysis was performed with the Tuxedo software suite. Genes and isoforms differing in abundance were queried for enrichment within functionally related gene groups using the Functional Annotation Tool of the DAVID Bioinformatics Database. While hematopoietic gene expression is initiated by E14, hepatocyte maturation is a gradual process involving clusters of genes responsible for response to nutrients and enzymes responsible for glycolysis and fatty acid catabolism. Following birth, a large cluster of differentially abundant genes was enriched for mitochondrial gene expression and cholesterol synthesis indicating that by 1 wk of age, the liver is engaged in lipid sensing and bile production. Clustering results for differentially abundant genes and isoforms were similar with the greatest difference for the E14/E16 comparison. Finally, a bioinformatic approach was used to annotate 1,307 novel liver transcripts assembled from sequences that aligned to intergenic regions of the rat genome.

Genome sequencing and comparative analysis of the carrot bacterial blight pathogen, Xanthomonas hortorum pv. carotae M081, for insights into pathogenicity and applications in molecular diagnostics.

Xanthomonas hortorum pv. carotae (Xhc) is an economically important pathogen of carrots. Its ability to epiphytically colonize foliar surfaces and infect seeds can result in bacterial blight of carrots when grown in warm and humid regions. We used high-throughput sequencing to determine the genome sequence of isolate M081 of Xhc. The short reads were de novo assembled and the resulting contigs were ordered using a syntenic reference genome sequence from X. campestris pv. campestris ATCC 33913. The improved, high-quality draft genome sequence of Xhc M081 is the first for its species. Despite its distance from other sequenced xanthomonads, Xhc M081 still shared a large inventory of orthologous genes, including many clusters of virulence genes common to other foliar pathogenic species of Xanthomonas. We also mined the genome sequence and identified at least 21 candidate type III effector genes. Two were members of the avrBs2 and xopQ families that demonstrably elicit effector-triggered immunity. We showed that expression in planta of these two type III effectors from Xhc M081 was sufficient to elicit resistance gene-mediated hypersensitive responses in heterologous plants, indicating a possibility for resistance gene-mediated control of Xhc. Finally, we identified regions unique to the Xhc M081 genome sequence, and demonstrated their potential in the design of molecular diagnostics for this pathogen.

Polyphyletic photosynthetic reaction centre genes in oligotrophic marine Gammaproteobacteria.

Ecological studies indicate that aerobic anoxygenic phototrophic bacteria (AAP) that use bacteriochlorophyll to support phototrophic electron transport are widely distributed in the oceans. All cultivated marine AAP are alpha-3 and alpha-4 Proteobacteria, but metagenomic evidence indicates that uncultured AAP Gammaproteobacteria are important members of ocean surface microbial communities. Here we report the description of obligately oligotrophic, marine Gammaproteobacteria that have genes for aerobic anoxygenic photosynthesis. Three strains belonging to the OM60 clade were isolated in autoclaved seawater media. Polymerase chain reaction assays for the pufM gene show that these strains contain photosynthetic reaction centre genes. DNA sequencing and phylogenetic analysis indicate that the pufM genes are polyphyletic, suggesting multiple instances of lateral gene transfer. Peptide sequences from six photosynthesis genes (pufL, pufM, pufC, pufB, pufA and puhA) were detected by proteomic analyses of strain HTCC2080 cells grown aerobically in seawater. They closely match predicted peptides from an environmental seawater bacterial artificial chromosome clone of gammaproteobacterial origin, thus identifying the OM60 clade as a significant source of gammaproteobacterial AAP genes in marine systems. The cell yield and rate of growth of HTCC2080 in autoclaved, aerobic seawater increased in the light. These findings identify the OM60 clade as a source of Gammaproteobacteria AAP genes in coastal oceans, and demonstrate that aerobic, anoxygenic photosynthetic metabolism can enhance the productivity of marine oligotrophic bacteria that also grow heterotrophically in darkness.

Toxicogenomic profiling of the hepatic tumor promoters indole-3-carbinol, 17beta-estradiol and beta-naphthoflavone in rainbow trout.

Indole-3-carbinol (I3C), from cruciferous vegetables, has been found to suppress or enhance tumors in several animal models. We previously reported that dietary I3C promotes hepatocarcinogenesis in rainbow trout (Oncorhynchus mykiss) at concentrations that differentially activated estrogen receptor (ER) or aryl hydrocarbon receptor (AhR)-mediated responses based on individual protein biomarkers. In this study, we evaluated the relative importance of these pathways as potential mechanisms for I3C on a global scale. Hepatic gene expression profiles were examined in trout after dietary exposure to 500 and 1500 ppm I3C and 3,3'-diindolylmethane (DIM), a major in vivo component of I3C, and were compared to the transcriptional signatures of two model hepatic tumor promoters: 17beta-estradiol (E2), an ER agonist, and beta-naphthoflavone, an AhR agonist. We demonstrate that I3C and DIM acted similar to E2 at the transcriptional level based on correlation analysis of expression profiles and clustering of gene responses. Of the genes regulated by E2 (fold change >or =2.0 or < or =0.50), most genes were regulated similarly by DIM (87-92%) and I3C (71%), suggesting a common mechanism of action. Of interest were upregulated genes associated with signaling pathways for cell growth and proliferation, vitellogenesis, and protein folding, stability, and transport. Other genes downregulated by E2, including those involved in acute-phase immune response, were also downregulated by DIM and I3C. Gene regulation was confirmed by qRT-PCR and Western blot. These data indicate I3C promotes hepatocarcinogenesis through estrogenic mechanisms in trout liver and suggest DIM may be an even more potent hepatic tumor promoter in this model.

Use of a rainbow trout oligonucleotide microarray to determine transcriptional patterns in aflatoxin B1-induced hepatocellular carcinoma compared to adjacent liver.

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide, and its occurrence is associated with a number of environmental factors including ingestion of the dietary contaminant aflatoxin B(1) (AFB(1)). Research over the last 40 years has revealed rainbow trout (Oncorhynchus mykiss) to be an excellent research model for study of AFB(1)-induced hepatocarcinogenesis; however, little is known about changes at the molecular level in trout tumors. We have developed a rainbow trout oligonucleotide array containing 1672 elements representing over 1400 genes of known or probable relevance to toxicology, comparative immunology, carcinogenesis, endocrinology, and stress physiology. In this study, we applied microarray technology to examine gene expression of AFB(1)-induced HCC in the rainbow trout tumor model. Carcinogenesis was initiated in trout embryos with 50 ppb AFB(1), and after 13 months control livers, tumors, and tumor-adjacent liver tissues were isolated from juvenile fish. Global gene expression was determined in histologically confirmed HCCs compared to noncancerous adjacent tissue and sham-initiated control liver. We observed distinct gene regulation patterns in HCC compared to noncancerous tissue including upregulation of genes important for cell cycle control, transcription, cytoskeletal formation, and the acute phase response and downregulation of genes involved in drug metabolism, lipid metabolism, and retinol metabolism. Interestingly, the expression profiles observed in trout HCC are similar to the transcriptional signatures found in human and rodent HCC, further supporting the validity of the model. Overall, these findings contribute to a better understanding of the mechanism of AFB(1)-induced hepatocarcinogenesis in trout and identify conserved genes important for carcinogenesis in species separated evolutionarily by more than 400 million years.

HIP14, a novel ankyrin domain-containing protein, links huntingtin to intracellular trafficking and endocytosis.

Huntington disease (HD) is caused by polyglutamine [poly(Q)] expansion in the protein huntingtin (htt). Although the exact mechanism of disease progression remains to be elucidated, altered interactions of mutant htt with its protein partners could contribute to the disease. Using the yeast two-hybrid system, we have isolated a novel htt interacting protein, HIP14. HIP14's interaction with htt is inversely correlated to the poly(Q) length in htt. mRNAs of 9 and 6 bp are transcribed from the HIP14 gene, with the 6 kb transcript being predominantly expressed in the brain. HIP14 protein is enriched in the brain, shows partial co-localization with htt in the striatum, and is found in medium spiny projection neurons, the subset of neurons affected in HD. HIP14 localizes to the Golgi, and to vesicles in the cytoplasm. The HIP14 protein has sequence similarity to Akr1p, a protein essential for endocytosis in Saccharomyces cerevisiae. Expression of human HIP14 results in rescue of the temperature-sensitive lethality in akr1 Delta yeast cells and, furthermore, restores their defect in endocytosis, demonstrating a role for HIP14 in intracellular trafficking. Our findings suggest that decreased interaction between htt and HIP14 could contribute to the neuronal dysfunction in HD by perturbing normal intracellular transport pathways in neurons.