Protocol for high-throughput DNA methylation profiling in rat tissues using automated reduced representation bisulfite sequencing.

Although reduced representation bisulfite sequencing (RRBS) measures DNA methylation (DNAme) across CpG-rich genomic regions with high sensitivity, the assay can be time-consuming and prone to batch effects. Here, we present a high-throughput, automated RRBS protocol starting with DNA extraction from frozen rat tissues. We describe steps for RRBS library preparation, library quality control, and sequencing. We also detail an optimized pipeline for sequencing data processing. This protocol has been applied successfully to DNAme profiling across multiple rat tissues. For complete details on the use and execution of this protocol, please refer to Nair et al.1.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

A methylation clock model of mild SARS-CoV-2 infection provides insight into immune dysregulation.

DNA methylation comprises a cumulative record of lifetime exposures superimposed on genetically determined markers. Little is known about methylation dynamics in humans following an acute perturbation, such as infection. We characterized the temporal trajectory of blood epigenetic remodeling in 133 participants in a prospective study of young adults before, during, and after asymptomatic and mildly symptomatic SARS-CoV-2 infection. The differential methylation caused by asymptomatic or mildly symptomatic infections was indistinguishable. While differential gene expression largely returned to baseline levels after the virus became undetectable, some differentially methylated sites persisted for months of follow-up, with a pattern resembling autoimmune or inflammatory disease. We leveraged these responses to construct methylation-based machine learning models that distinguished samples from pre-, during-, and postinfection time periods, and quantitatively predicted the time since infection. The clinical trajectory in the young adults and in a diverse cohort with more severe outcomes was predicted by the similarity of methylation before or early after SARS-CoV-2 infection to the model-defined postinfection state. Unlike the phenomenon of trained immunity, the postacute SARS-CoV-2 epigenetic landscape we identify is antiprotective.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training.

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

SARS-CoV-2 Outbreak Dynamics in an Isolated US Military Recruit Training Center With Rigorous Prevention Measures.

BACKGROUND: Marine recruits training at Parris Island experienced an unexpectedly high rate of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, despite preventive measures including a supervised, 2-week, pre-entry quarantine. We characterize SARS-CoV-2 transmission in this cohort. METHODS: Between May and November 2020, we monitored 2,469 unvaccinated, mostly male, Marine recruits prospectively during basic training. If participants tested negative for SARS-CoV-2 by quantitative polymerase chain reaction (qPCR) at the end of quarantine, they were transferred to the training site in segregated companies and underwent biweekly testing for 6 weeks. We assessed the effects of coronavirus disease 2019 (COVID-19) prevention measures on other respiratory infections with passive surveillance data, performed phylogenetic analysis, and modeled transmission dynamics and testing regimens. RESULTS: Preventive measures were associated with drastically lower rates of other respiratory illnesses. However, among the trainees, 1,107 (44.8%) tested SARS-CoV-2-positive, with either mild or no symptoms. Phylogenetic analysis of viral genomes from 580 participants revealed that all cases but one were linked to five independent introductions, each characterized by accumulation of mutations across and within companies, and similar viral isolates in individuals from the same company. Variation in company transmission rates (mean reproduction number R 0 ; 5.5 [95% confidence interval [CI], 5.0, 6.1]) could be accounted for by multiple initial cases within a company and superspreader events. Simulations indicate that frequent rapid-report testing with case isolation may minimize outbreaks. CONCLUSIONS: Transmission of wild-type SARS-CoV-2 among Marine recruits was approximately twice that seen in the community. Insights from SARS-CoV-2 outbreak dynamics and mutations spread in a remote, congregate setting may inform effective mitigation strategies.

Multi-omics profiling of single nuclei from frozen archived postmortem human pituitary tissue.

Concomitant profiling of transcriptome and chromatin accessibility in isolated nuclei can reveal gene regulatory control mechanisms in health and disease. We report a single nucleus multi-omics analysis protocol optimized for frozen archived postmortem human pituitaries that is also effective for frozen ovine and murine pituitaries and human skeletal muscle biopsies. Its main advantages are that (1) it is not limited to fresh tissue, (2) it avoids tissue dissociation-induced transcriptional changes, and (3) it includes a novel, automated quality control pipeline. For complete details on the use and execution of this protocol, please refer to Ruf-Zamojski et al. (2021) and Zhang et al. (2022).

Optimization of the Omni-ATAC protocol to chromatin accessibility profiling in snap-frozen rat adipose and muscle tissues.

ATAC-seq is a fast and sensitive method for the epigenomic profiling of open chromatin and for mapping of transcription factor binding sites [1]. Despite the development of the Omni-ATAC protocol for the profiling of chromatin accessibility in frozen tissues [2], studies in adipose tissue have been restricted due to technical challenges including the high lipid content of adipocytes and reproducibility issues between replicates. Here, we provide a modified Omni-ATAC protocol that achieves high data reproducibility in various tissue types from rat, including adipose and muscle tissues [3].•This protocol describes a methodology that enables chromatin accessibility profiling from snap-frozen rat adipose and muscle tissues.•The technique comprises an optimized bead-based tissue homogenization process that substitutes to Dounce homogenization, reduces variability in the experimental procedure, and is adaptable to various tissue types.•In comparison with the Omni-ATAC protocol, the method described here results in improved ATAC-seq data quality that complies with ENCODE quality standards.

Single nucleus transcriptome and chromatin accessibility of postmortem human pituitaries reveal diverse stem cell regulatory mechanisms.

Despite their importance in tissue homeostasis and renewal, human pituitary stem cells (PSCs) are incompletely characterized. We describe a human single nucleus RNA-seq and ATAC-seq resource from pediatric, adult, and aged postmortem pituitaries (snpituitaryatlas.princeton.edu) and characterize cell-type-specific gene expression and chromatin accessibility programs for all major pituitary cell lineages. We identify uncommitted PSCs, committing progenitor cells, and sex differences. Pseudotime trajectory analysis indicates that early-life PSCs are distinct from the other age groups. Linear modeling of same-cell multiome data identifies regulatory domain accessibility sites and transcription factors that are significantly associated with gene expression in PSCs compared with other cell types and within PSCs. We identify distinct deterministic mechanisms that contribute to heterogeneous marker expression within PSCs. These findings characterize human stem cell lineages and reveal diverse mechanisms regulating key PSC genes and cell type identity.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle.

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Single nucleus multi-omics regulatory landscape of the murine pituitary.

To provide a multi-omics resource and investigate transcriptional regulatory mechanisms, we profile the transcriptome, chromatin accessibility, and methylation status of over 70,000 single nuclei (sn) from adult mouse pituitaries. Paired snRNAseq and snATACseq datasets from individual animals highlight a continuum between developmental epigenetically-encoded cell types and transcriptionally-determined transient cell states. Co-accessibility analysis-based identification of a putative Fshb cis-regulatory domain that overlaps the fertility-linked rs11031006 human polymorphism, followed by experimental validation illustrate the use of this resource for hypothesis generation. We also identify transcriptional and chromatin accessibility programs distinguishing each major cell type. Regulons, which are co-regulated gene sets sharing binding sites for a common transcription factor driver, recapitulate cell type clustering. We identify both cell type-specific and sex-specific regulons that are highly correlated with promoter accessibility, but not with methylation state, supporting the centrality of chromatin accessibility in shaping cell-defining transcriptional programs. The sn multi-omics atlas is accessible at snpituitaryatlas.princeton.edu.

Sedentary and Trained Older Men Have Distinct Circulating Exosomal microRNA Profiles at Baseline and in Response to Acute Exercise.

Exercise has multi-systemic benefits and attenuates the physiological impairments associated with aging. Emerging evidence suggests that circulating exosomes mediate some of the beneficial effects of exercise via the transfer of microRNAs between tissues. However, the impact of regular exercise and acute exercise on circulating exosomal microRNAs (exomiRs) in older populations remains unknown. In the present study, we analyzed circulating exomiR expression in endurance-trained elderly men (n = 5) and age-matched sedentary males (n = 5) at baseline (Pre), immediately after a forty minute bout of aerobic exercise on a cycle ergometer (Post), and three hours after this acute exercise (3hPost). Following the isolation and enrichment of exosomes from plasma, exosome-enriched preparations were characterized and exomiR levels were determined by sequencing. The effect of regular exercise on circulating exomiRs was assessed by comparing the baseline expression levels in the trained and sedentary groups. The effect of acute exercise was determined by comparing baseline and post-training expression levels in each group. Regular exercise resulted in significantly increased baseline expression of three exomiRs (miR-486-5p, miR-215-5p, miR-941) and decreased expression of one exomiR (miR-151b). Acute exercise altered circulating exomiR expression in both groups. However, exomiRs regulated by acute exercise in the trained group (7 miRNAs at Post and 8 at 3hPost) were distinct from those in the sedentary group (9 at Post and 4 at 3hPost). Pathway analysis prediction and reported target validation experiments revealed that the majority of exercise-regulated exomiRs are targeting genes that are related to IGF-1 signaling, a pathway involved in exercise-induced muscle and cardiac hypertrophy. The immediately post-acute exercise exomiR signature in the trained group correlates with activation of IGF-1 signaling, whereas in the sedentary group it is associated with inhibition of IGF-1 signaling. While further validation is needed, including measurements of IGF-1/IGF-1 signaling in blood or skeletal muscle, our results suggest that training status may counteract age-related anabolic resistance by modulating circulating exomiR profiles both at baseline and in response to acute exercise.

Cytogenetic, Genomic, and Functional Characterization of Pituitary Gonadotrope Cell Lines.

LßT2 and αT3-1 are important, widely studied cell line models for the pituitary gonadotropes that were generated by targeted tumorigenesis in transgenic mice. LßT2 cells are more mature gonadotrope precursors than αT3-1 cells. Microsatellite authentication patterns, chromosomal characteristics, and their intercellular variation have not been reported. We performed microsatellite and cytogenetic analysis of both cell types at early passage numbers. Short tandem repeat (STR) profiling was consistent with a mixed C57BL/6J × BALB/cJ genetic background, with distinct patterns for each cell type. Spectral karyotyping in αT3-1 cells revealed cell-to-cell variation in chromosome composition and pseudodiploidy. In LßT2 cells, chromosome counting and karyotyping demonstrated pseudotriploidy and high chromosomal variation among cells. Chromosome copy number variation was confirmed by single-cell DNA sequencing. Chromosomal compositions were consistent with a male sex for αT3-1 and a female sex for LßT2 cells. Among LßT2 stocks used in multiple laboratories, we detected two genetically similar but distinguishable lines via STR authentication, LßT2a and LßT2b. The two lines differed in morphological appearance, with LßT2a having significantly smaller cell and nucleus areas. Analysis of immediate early gene and gonadotropin subunit gene expression revealed variations in basal expression and responses to continuous and pulsatile GnRH stimulation. LßT2a showed higher basal levels of Egr1, Fos, and Lhb but lower Fos induction. Fshb induction reached significance only in LßT2b cells. Our study highlights the heterogeneity in gonadotrope cell line genomes and provides reference STR authentication patterns that can be monitored to improve experimental reproducibility and facilitate comparisons of results within and across laboratories.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Regulatory Architecture of the LßT2 Gonadotrope Cell Underlying the Response to Gonadotropin-Releasing Hormone.

The LßT2 mouse pituitary cell line has many characteristics of a mature gonadotrope and is a widely used model system for studying the developmental processes and the response to gonadotropin-releasing hormone (GnRH). The global epigenetic landscape, which contributes to cell-specific gene regulatory mechanisms, and the single-cell transcriptome response variation of LßT2 cells have not been previously investigated. Here, we integrate the transcriptome and genome-wide chromatin accessibility state of LßT2 cells during GnRH stimulation. In addition, we examine cell-to-cell variability in the transcriptional response to GnRH using Gel bead-in-Emulsion Drop-seq technology. Analysis of a bulk RNA-seq data set obtained 45 min after exposure to either GnRH or vehicle identified 112 transcripts that were regulated >4-fold by GnRH (FDR < 0.05). The top regulated transcripts constitute, as determined by Bayesian massive public data integration analysis, a human pituitary-relevant coordinated gene program. Chromatin accessibility [assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq)] data sets generated from GnRH-treated LßT2 cells identified more than 58,000 open chromatin regions, some containing notches consistent with bound transcription factor footprints. The study of the most prominent open regions showed that 75% were in transcriptionally active promoters or introns, supporting their involvement in active transcription. Lhb, Cga, and Egr1 showed significantly open chromatin over their promoters. While Fshb was closed over its promoter, several discrete significantly open regions were found at -40 to -90 kb, which may represent novel upstream enhancers. Chromatin accessibility determined by ATAC-seq was associated with high levels of gene expression determined by RNA-seq. We obtained high-quality single-cell Gel bead-in-Emulsion Drop-seq transcriptome data, with an average of >4,000 expressed genes/cell, from 1,992 vehicle- and 1,889 GnRH-treated cells. While the individual cell expression patterns showed high cell-to-cell variation, representing both biological and measurement variation, the average expression patterns correlated well with bulk RNA-seq data. Computational assignment of each cell to its precise cell cycle phase showed that the response to GnRH was unaffected by cell cycle. To our knowledge, this study represents the first genome-wide epigenetic and single-cell transcriptomic characterization of this important gonadotrope model. The data have been deposited publicly and should provide a resource for hypothesis generation and further study.

Pandemic H1N1 influenza A viruses suppress immunogenic RIPK3-driven dendritic cell death.

The risk of emerging pandemic influenza A viruses (IAVs) that approach the devastating 1918 strain motivates finding strain-specific host-pathogen mechanisms. During infection, dendritic cells (DC) mature into antigen-presenting cells that activate T cells, linking innate to adaptive immunity. DC infection with seasonal IAVs, but not with the 1918 and 2009 pandemic strains, induces global RNA degradation. Here, we show that DC infection with seasonal IAV causes immunogenic RIPK3-mediated cell death. Pandemic IAV suppresses this immunogenic DC cell death. Only DC infected with seasonal IAV, but not with pandemic IAV, enhance maturation of uninfected DC and T cell proliferation. In vivo, circulating T cell levels are reduced after pandemic, but not seasonal, IAV infection. Using recombinant viruses, we identify the HA genomic segment as the mediator of cell death inhibition. These results show how pandemic influenza viruses subvert the immune response.

Modeling and high-throughput experimental data uncover the mechanisms underlying Fshb gene sensitivity to gonadotropin-releasing hormone pulse frequency.

Neuroendocrine control of reproduction by brain-secreted pulses of gonadotropin-releasing hormone (GnRH) represents a longstanding puzzle about extracellular signal decoding mechanisms. GnRH regulates the pituitary gonadotropin's follicle-stimulating hormone (FSH) and luteinizing hormone (LH), both of which are heterodimers specified by unique ß subunits (FSHß/LHß). Contrary to Lhb, Fshb gene induction has a preference for low-frequency GnRH pulses. To clarify the underlying regulatory mechanisms, we developed three biologically anchored mathematical models: 1) parallel activation of Fshb inhibitory factors (e.g. inhibin α and VGF nerve growth factor-inducible), 2) activation of a signaling component with a refractory period (e.g. G protein), and 3) inactivation of a factor needed for Fshb induction (e.g. growth differentiation factor 9). Simulations with all three models recapitulated the Fshb expression levels obtained in pituitary gonadotrope cells perifused with varying GnRH pulse frequencies. Notably, simulations altering average concentration, pulse duration, and pulse frequency revealed that the apparent frequency-dependent pattern of Fshb expression in model 1 actually resulted from variations in average GnRH concentration. In contrast, models 2 and 3 showed "true" pulse frequency sensing. To resolve which components of this GnRH signal induce Fshb, we developed a high-throughput parallel experimental system. We analyzed over 4,000 samples in experiments with varying near-physiological GnRH concentrations and pulse patterns. Whereas Egr1 and Fos genes responded only to variations in average GnRH concentration, Fshb levels were sensitive to both average concentration and true pulse frequency. These results provide a foundation for understanding the role of multiple regulatory factors in modulating Fshb gene activity.

Characterization of Gonadotrope Secretoproteome Identifies Neurosecretory Protein VGF-derived Peptide Suppression of Follicle-stimulating Hormone Gene Expression.

Reproductive function is controlled by the pulsatile release of hypothalamic gonadotropin-releasing hormone (GnRH), which regulates the expression of the gonadotropins luteinizing hormone and FSH in pituitary gonadotropes. Paradoxically, Fshb gene expression is maximally induced at lower frequency GnRH pulses, which provide a very low average concentration of GnRH stimulation. We studied the role of secreted factors in modulating gonadotropin gene expression. Inhibition of secretion specifically disrupted gonadotropin subunit gene regulation but left early gene induction intact. We characterized the gonadotrope secretoproteome and global mRNA expression at baseline and after Gαs knockdown, which has been found to increase Fshb gene expression (1). We identified 1077 secreted proteins or peptides, 19 of which showed mRNA regulation by GnRH or/and Gαs knockdown. Among several novel secreted factors implicated in Fshb gene regulation, we focused on the neurosecretory protein VGF. Vgf mRNA, whose gene has been implicated in fertility (2), exhibited high induction by GnRH and depended on Gαs In contrast with Fshb induction, Vgf induction occurred preferentially at high GnRH pulse frequency. We hypothesized that a VGF-derived peptide might regulate Fshb gene induction. siRNA knockdown or extracellular immunoneutralization of VGF augmented Fshb mRNA induction by GnRH. GnRH stimulated the secretion of the VGF-derived peptide NERP1. NERP1 caused a concentration-dependent decrease in Fshb gene induction. These findings implicate a VGF-derived peptide in selective regulation of the Fshb gene. Our results support the concept that signaling specificity from the cell membrane GnRH receptor to the nuclear Fshb gene involves integration of intracellular signaling and exosignaling regulatory motifs.

Low-variance RNAs identify Parkinson's disease molecular signature in blood.

The diagnosis of Parkinson's disease (PD) is usually not established until advanced neurodegeneration leads to clinically detectable symptoms. Previous blood PD transcriptome studies show low concordance, possibly resulting from the use of microarray technology, which has high measurement variation. The Leucine-rich repeat kinase 2 (LRRK2) G2019S mutation predisposes to PD. Using preclinical and clinical studies, we sought to develop a novel statistically motivated transcriptomic-based approach to identify a molecular signature in the blood of Ashkenazi Jewish PD patients, including LRRK2 mutation carriers. Using a digital gene expression platform to quantify 175 messenger RNA (mRNA) markers with low coefficients of variation (CV), we first compared whole-blood transcript levels in mouse models (1) overexpressing wild-type (WT) LRRK2, (2) overexpressing G2019S LRRK2, (3) lacking LRRK2 (knockout), and (4) and in WT controls. We then studied an Ashkenazi Jewish cohort of 34 symptomatic PD patients (both WT LRRK2 and G2019S LRRK2) and 32 asymptomatic controls. The expression profiles distinguished the four mouse groups with different genetic background. In patients, we detected significant differences in blood transcript levels both between individuals differing in LRRK2 genotype and between PD patients and controls. Discriminatory PD markers included genes associated with innate and adaptive immunity and inflammatory disease. Notably, gene expression patterns in levodopa-treated PD patients were significantly closer to those of healthy controls in a dose-dependent manner. We identify whole-blood mRNA signatures correlating with LRRK2 genotype and with PD disease state. This approach may provide insight into pathogenesis and a route to early disease detection.

High incidence of LRAT promoter hypermethylation in colorectal cancer correlates with tumor stage.

Lecithin:retinol acyltransferase (LRAT) is a major enzyme involved in vitamin A/retinol metabolism, which regulates various physiological processes like cell proliferation and differentiation. LRAT expression is reduced in numerous cancers, yet the underlying mechanisms have remained undefined. We hypothesized that methylation silencing may contribute to decreased LRAT gene expression in colorectal cancer (CRC). LRAT hypermethylation status was analyzed in five CRC cell lines, 167 colorectal tumors, and 69 adjacent normal colonic mucosae, using a quantitative bisulfite/PCR/LDR/Universal Array assay. LRAT transcription levels were determined by real-time RT-PCR in a subset of tumors and matched normal tissues and in CRC cell lines that were treated with a demethylating agent, 5-aza-2'-deoxycytidine. The incidence of LRAT hypermethylation was significantly higher in colorectal tumors than in adjacent normal mucosae (p = 0.0025). Aberrant methylation occurred in 51 % of microsatellite-stable CRCs, in 84 % of microsatellite-unstable CRCs, and in 12 out of 13 colonic polyps. The number of hypermethylated LRAT events was inversely correlated with CRC stage (p < 0.0001). Importantly, LRAT hypermethylation was associated with decreased mRNA level in CRC clinical specimens, and demethylation treatment resulted in LRAT transcriptional reactivation. Our data support the idea that LRAT promoter hypermethylation associates with LRAT gene expression in CRC. The higher frequency of LRAT hypermethylation in colonic polyps and early-stage CRCs indicates that it may occur early in malignant progression.

Growth differentiation factor 9 (GDF9) forms an incoherent feed-forward loop modulating follicle-stimulating hormone ß-subunit (FSHß) gene expression.

Gonadotropin-releasing hormone (GnRH) is secreted in brief pulses from the hypothalamus and regulates follicle-stimulating hormone ß-subunit (FSHß) gene expression in pituitary gonadotropes in a frequency-sensitive manner. The mechanisms underlying its preferential and paradoxical induction of FSHß by low frequency GnRH pulses are incompletely understood. Here, we identify growth differentiation factor 9 (GDF9) as a GnRH-suppressed autocrine inducer of FSHß gene expression. GDF9 gene transcription and expression were preferentially decreased by high frequency GnRH pulses. GnRH regulation of GDF9 was concentration-dependent and involved ERK and PKA. GDF9 knockdown or immunoneutralization reduced FSHß mRNA expression. Conversely, exogenous GDF9 induced FSHß expression in immortalized gonadotropes and in mouse primary pituitary cells. GDF9 exposure increased FSH secretion in rat primary pituitary cells. GDF9 induced Smad2/3 phosphorylation, which was impeded by ALK5 knockdown and by activin receptor-like kinase (ALK) receptor inhibitor SB-505124, which also suppressed FSHß expression. Smad2/3 knockdown indicated that FSHß induction by GDF9 involved Smad2 and Smad3. FSHß mRNA induction by GDF9 and GnRH was synergistic. We hypothesized that GDF9 contributes to a regulatory loop that tunes the GnRH frequency-response characteristics of the FSHß gene. To test this, we determined the effects of GDF9 knockdown on FSHß induction at different GnRH pulse frequencies using a parallel perifusion system. Reduction of GDF9 shifted the characteristic pattern of GnRH pulse frequency sensitivity. These results identify GDF9 as contributing to an incoherent feed-forward loop, comprising both intracellular and secreted components, that regulates FSHß expression in response to activation of cell surface GnRH receptors.