Class-C SOX transcription factors control GnRH gene expression via the intronic transcriptional enhancer

SOX11 variants cause a neurodevelopmental disorder with infrequent ocular malformations and hypogonadotropic hypogonadism and with distinct DNA methylation profile

PURPOSE

This study aimed to undertake a multidisciplinary characterization of the phenotype associated with SOX11 variants.

METHODS

Individuals with protein altering variants in SOX11 were identified through exome and genome sequencing and international data sharing. Deep clinical phenotyping was undertaken by referring clinicians. Blood DNA methylation was assessed using Infinium MethylationEPIC array. The expression pattern of SOX11 in developing human brain was defined using RNAscope.

RESULTS

We reported 38 new patients with SOX11 variants. Idiopathic hypogonadotropic hypogonadism was confirmed as a feature of SOX11 syndrome. A distinctive pattern of blood DNA methylation was identified in SOX11 syndrome, separating SOX11 syndrome from other BAFopathies.

CONCLUSION

SOX11 syndrome is a distinct clinical entity with characteristic clinical features and episignature differentiating it from BAFopathies.

Pituitary Gonadotropin Gene Expression During Induced Onset of Postsmolt Maturation in Male Atlantic Salmon: In Vivo and Tissue Culture Studies

Precocious male maturation causes reduced welfare and increased production costs in Atlantic salmon (Salmo salar) aquaculture. The pituitary produces and releases follicle-stimulating hormone (Fsh), the gonadotropin triggering puberty in male salmonids. However, little is known about how Fsh production is regulated in Atlantic salmon. We examined, in vivo and ex vivo, transcriptional changes of gonadotropin-related genes accompanying the initial steps of testis maturation, in pituitaries of males exposed to photoperiod and temperature conditions promoting maturation (constant light and 16°C). Pituitary fshb, lhb and gnrhr2bba transcripts increased in vivo in maturing males (gonado-somatic index > 0.1%). RNA sequencing (RNAseq) analysis using pituitaries from genetically similar males carrying the same genetic predisposition to mature, but differing by responding or not responding to stimulatory environmental conditions, revealed 144 differentially expressed genes, ~2/3rds being up-regulated in responders, including fshb and other pituitary hormones, steroid-related and other puberty-associated transcripts. Functional enrichment analyses confirmed gene involvement in hormone/steroid production and gonad development. In ex vivo studies, whole pituitaries were exposed to a selection of hormones and growth factors. Gonadotropin-releasing hormone (Gnrh), 17β-estradiol (E₂) and 11-ketotestosterone (11-KT) up-regulated gnrhr2bba and lhb, while fshb was up-regulated by Gnrh but down-regulated by 11-KT in pituitaries from immature males. Also pituitaries from maturing males responded to Gnrh and sex steroids by increased gnrhr2bba and lhb transcript levels, but fshb expression remained unchanged. Growth factors (inhibin A, activin A and insulin-like growth factor 1) did not change gnrhr2bba, lhb or fshb transcript levels in pituitaries either from immature or maturing males. Additional pituitary ex vivo studies on candidates identified by RNAseq showed that these transcripts were preferentially regulated by Gnrh and sex steroids, but not by growth factors, and that Gnrh/sex steroids were less effective when incubating pituitaries from maturing males. Our results suggest that a yet to be characterized mechanism up-regulating fshb expression in the salmon pituitary is activated in response to stimulatory environmental conditions prior to morphological signs of testis maturation, and that the transcriptional program associated with this mechanism becomes unresponsive or less responsive to most stimulators ex vivo once males had entered pubertal developmental in vivo.

Clinical and molecular features of idiopathic hypogonadotropic hypogonadism in Taiwan: A single center experience

BACKGROUND

Idiopathic (isolated) hypogonadotropic hypogonadism (IHH) is a rare disease that can be classified as Kallmann syndrome (KS) or normosmic IHH (nIHH). This study investigated the phenotype and genotype of IHH in Taiwanese patients.

METHODS

Twenty-six unrelated IHH patients were included in this study and their clinical, hormonal, and radiological findings were analyzed retrospectively. Whole exome sequencing (WES) was performed to identify the etiology.

RESULTS

The 26 patients (M:F = 19:7) were divided into a KS group (n = 11) and a nIHH group (n = 15). The diagnosis was earlier in boys than in girls. Fifteen patients were found to have pathogenic/likely pathogenic (P/LP) variants of IHH-associated genes, and the mutation detection rate was 58%. CHD7, FGFR1, and ANOS1 were the most common genetic etiologies identified in this group. Two patients with nIHH were found to have de novo SOX11 mutations and Coffin-Siris syndrome features. After treatment, the height outcomes and secondary sexual characteristics were significantly improved. There were no obvious differences between the genetically resolved (GR), variants of uncertain significance (VUS) and genetically unresolved groups (GUR).

CONCLUSION

Whole exome sequencing is useful in patients with IHH, and we identified the SOX11 gene as a causal factor in this study. We described the clinical, hormonal, and molecular characteristics, and the treatment outcomes, of Taiwanese patients with IHH, which should aid therapeutic planning and further research.

Gonadotropin-releasing hormone neuron development in vertebrates

Gonadotropin-releasing hormone (GnRH) neurons are master regulators of the reproductive axis in vertebrates. During early mammalian embryogenesis, GnRH1 neurons emerge in the nasal/olfactory placode. These neurons undertake a long-distance migration, moving from the nose to the preoptic area and hypothalamus. While significant advances have been made in understanding the functional importance of the GnRH1 neurons in reproduction, where GnRH1 neurons come from and how are they specified during early development is still under debate. In addition to the GnRH1 gene, most vertebrate species including humans have one or two additional GnRH genes. Compared to the GnRH1 neurons, much less is known about the development and regulation of GnRH2 neuron and GnRH3 neurons. The objective of this article is to review what is currently known about GnRH neuron development. We will survey various cell autonomous and non-autonomous factors implicated in the regulation of GnRH neuron development. Finally, we will discuss emerging tools and new approaches to resolve open questions pertaining to GnRH neuron development.

Characterization of the human GnRH neuron developmental transcriptome using a GNRH1-TdTomato reporter line in human pluripotent stem cells

Gonadotropin-releasing hormone (GnRH) neurons provide a fundamental signal for the onset of puberty and subsequent reproductive functions by secretion of gonadotropin-releasing hormone. Their disrupted development or function leads to congenital hypogonadotropic hypogonadism (CHH). To model the development of human GnRH neurons, we generated a stable GNRH1-TdTomato reporter cell line in human pluripotent stem cells (hPSCs) using CRISPR-Cas9 genome editing. RNA-sequencing of the reporter clone, differentiated into GnRH neurons by dual SMAD inhibition and FGF8 treatment, revealed 6461 differentially expressed genes between progenitors and GnRH neurons. Expression of the transcription factor ISL1, one of the top 50 most upregulated genes in the TdTomato-expressing GnRH neurons, was confirmed in 10.5 gestational week-old human fetal GnRH neurons. Among the differentially expressed genes, we detected 15 genes that are implicated in CHH and several genes that are implicated in human puberty timing. Finally, FGF8 treatment in the neuronal progenitor pool led to upregulation of 37 genes expressed both in progenitors and in TdTomato-expressing GnRH neurons, which suggests upstream regulation of these genes by FGF8 signaling during GnRH neuron differentiation. These results illustrate how hPSC-derived human GnRH neuron transcriptomic analysis can be utilized to dissect signaling pathways and gene regulatory networks involved in human GnRH neuron development.This article has an associated First Person interview with the first author of the paper.

Nasal Placode Development, GnRH Neuronal Migration and Kallmann Syndrome

The development of Gonadotropin releasing hormone-1 (GnRH) neurons is important for a functional reproduction system in vertebrates. Disruption of GnRH results in hypogonadism and if accompanied by anosmia is termed Kallmann Syndrome (KS). From their origin in the nasal placode, GnRH neurons migrate along the olfactory-derived vomeronasal axons to the nasal forebrain junction and then turn caudally into the developing forebrain. Although research on the origin of GnRH neurons, their migration and genes associated with KS has identified multiple factors that influence development of this system, several aspects still remain unclear. This review discusses development of the olfactory system, factors that regulate GnRH neuron formation and development of the olfactory system, migration of the GnRH neurons from the nose into the brain, and mutations in humans with KS that result from disruption of normal GnRH/olfactory systems development.

Central Precocious Puberty Caused by Novel Mutations in the Promoter and 5'-UTR Region of the Imprinted MKRN3 Gene

Background: Central Precocious Puberty (CPP) is clinically defined by the development of secondary sexual characteristics before the age of 8 years in girls and 9 years in boys. To date, mutations in the coding region of KISS1, KISS1R, PROKR2, DLK1, and MKRN3 genes have been reported as causative for CPP. This study investigated the presence of causative mutations in both the promoter and the 5'-UTR regions of the MKRN3 gene. Methods: Sanger DNA sequencing was used for screening the proximal promoter and 5'-UTR region of the MKRN3 gene in a group of 73 index girls with CPP. Mutations identified were cloned in luciferase reporter gene vectors and transiently transfected in GN11 cells in order to check for changes in the activity of the MKRN3 promoter. GN11 cells were previously checked for Mkrn3 expression using lentivirus mediated knock-down. In silico analysis was implemented for the detection of changes in the mRNA secondary structure of the mutated MKRN3 5'-UTR. Results: Three novel heterozygous mutations (-166, -865, -886 nt upstream to the transcription start site) located in the proximal promoter region of the MKRN3 gene were identified in six non-related girls with CPP. Four of these girls shared the -865 mutation, one the -166, and another one the -886. A 5'-UTR (+13 nt downstream to the transcription start site) novel mutation was also identified in a girl with similar clinical phenotype. Gene reporter assay evaluated the identified promoter mutations and demonstrated a significant reduction of MKRN3 promoter activity in transfected GN11 cells. In silico analysis for the mutated 5'-UTR predicted a significant change of the mRNA secondary structure. The minimum free energy (MFE) of the mutated 5'-UTR was higher when compared to the corresponding wild-type indicating less stable RNA secondary structure. Conclusion: Our findings demonstrated novel genetic alterations in the promoter and 5'-UTR regulatory regions of the MKRN3 gene. These changes add to another region to check for the etiology of CPP.

Parallel adaptive evolution of geographically distant herring populations on both sides of the North Atlantic Ocean

Atlantic herring is an excellent species for studying the genetic basis of adaptation in geographically distant populations because of its characteristically large population sizes and low genetic drift. In this study we compared whole-genome resequencing data of Atlantic herring populations from both sides of the Atlantic Ocean. An important finding was the very low degree of genetic differentiation among geographically distant populations (fixation index = 0.026), suggesting lack of reproductive isolation across the ocean. This feature of the Atlantic herring facilitates the detection of genetic factors affecting adaptation because of the sharp contrast between loci showing genetic differentiation resulting from natural selection and the low background noise resulting from genetic drift. We show that genetic factors associated with timing of reproduction are shared between genetically distinct and geographically distant populations. The genes for thyroid-stimulating hormone receptor (TSHR), the SOX11 transcription factor (SOX11), calmodulin (CALM), and estrogen receptor 2 (ESR2A), all with a significant role in reproductive biology, were among the loci that showed the most consistent association with spawning time throughout the species range. In fact, the same two SNPs located at the 5' end of TSHR showed the most significant association with spawning time in both the east and west Atlantic. We also identified unexpected haplotype sharing between spring-spawning oceanic herring and autumn-spawning populations across the Atlantic Ocean and the Baltic Sea. The genomic regions showing this pattern are unlikely to control spawning time but may be involved in adaptation to ecological factor(s) shared among these populations.

The Arcuate Estrogen-Regulated Transcriptome: Estrogen Response Element-Dependent and -Independent Signaling of ERα in Female Mice

To influence energy homeostasis and reproduction, 17β-estradiol (E2) controls the arcuate nucleus (ARC) through multiple receptor-mediated mechanisms, but primarily via estrogen receptor (ER) α, which signals through both estrogen response element (ERE)-dependent and -independent mechanisms. To determine ERα-mediated, ERE-dependent, and ERE-independent E2 signaling in the ARC, we examined the differential regulation of the mouse arcuate transcriptome by E2 using three mice genotypes: (1) wild-type, (2) ERα knock-in/knockout (ERE-independent mechanisms), and (3) total ERα knockout (ERα-independent mechanisms). Females were ovariectomized and injected with oil or E2, and RNA sequencing on the ARC was used to identify E2-regulated genes in each genotype. Our results show that E2 regulates numerous genes involved in cell signaling, cytoskeleton structure, inflammation, neurotransmission, neuropeptide production, and transcription. Furthermore, ERE-independent signaling regulates ARC genes expressed in kisspeptin neurons and transcription factors that control the hypothalamic/pituitary/gonadal axis. Interestingly, a few genes involved in mitochondrial oxidative respiration were regulated by E2 through ERα-independent signaling. A comparison within oil- and E2-treated females across the three genotypes suggests that genes involved in cell growth and proliferation, extracellular matrices, neuropeptides, receptors, and transcription are differentially expressed across the genotypes. Comparing with previously published chromatin immunoprecipitation sequencing analysis, we found that ERE-independent regulation in the ARC is mainly mediated by tethering of ERα, which is consistent with previous findings. We conclude that the mouse arcuate estrogen-regulated transcriptome is regulated by multiple receptor-mediated mechanisms to modulate the central control of energy homeostasis and reproduction, including novel E2-responsive pathways.

Tumor Necrosis Factor-α Impairs Kisspeptin Signaling in Human Gonadotropin-Releasing Hormone Primary Neurons

CONTEXT

Inflammatory pathways may impair central regulatory networks involving gonadotropin-releasing hormone (GnRH) neuron activity. Studies in humans are limited by the lack of human GnRH neuron cell lines.

OBJECTIVE

To establish an in vitro model of human GnRH neurons and analyze the effects of proinflammatory cytokines.

DESIGN

The primary human fetal hypothalamic cells (hfHypo) were isolated from 12-week-old fetuses. Responsiveness to kisspeptin, the main GnRH neurons' physiological regulator, was evaluated for biological characterization. Tumor necrosis factor alpha (TNF-α) was used as a proinflammatory stimulus. Main Outcome Measures: Expression of specific GnRH neuron markers by quantitative reverse transcription-polymerase chain reaction, flow cytometry, and immunocytochemistry analyses; and GnRH-releasing ability and electrophysiological changes in response to kisspeptin.

RESULTS

The primary hfHypo showed a high percentage of GnRH-positive cells (80%), expressing a functional kisspeptin receptor (KISS1R) and able to release GnRH in response to kisspeptin. TNF-α exposure determined a specific inflammatory intracellular signaling and reduced GnRH secretion, KISS1R expression, and kisspeptin-induced depolarizing effect. Moreover, hfHypo possessed a primary cilium, whose assembly was inhibited by TNF-α.

CONCLUSION

The hfHypo cells represent a novel tool for investigations on human GnRH neuron biology. TNF-α directly affects GnRH neuron function by interfering with KISS1R expression and ciliogenesis, thereby impairing kisspeptin signaling.

Mice carrying a human GLUD2 gene recapitulate aspects of human transcriptome and metabolome development

Whereas all mammals have one glutamate dehydrogenase gene (GLUD1), humans and apes carry an additional gene (GLUD2), which encodes an enzyme with distinct biochemical properties. We inserted a bacterial artificial chromosome containing the human GLUD2 gene into mice and analyzed the resulting changes in the transcriptome and metabolome during postnatal brain development. Effects were most pronounced early postnatally, and predominantly genes involved in neuronal development were affected. Remarkably, the effects in the transgenic mice partially parallel the transcriptome and metabolome differences seen between humans and macaques analyzed. Notably, the introduction of GLUD2 did not affect glutamate levels in mice, consistent with observations in the primates. Instead, the metabolic effects of GLUD2 center on the tricarboxylic acid cycle, suggesting that GLUD2 affects carbon flux during early brain development, possibly supporting lipid biosynthesis.

The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.

Genomic structure and promoter functional analysis of GnRH3 gene in large yellow croaker (Larimichthys crocea)

Gonadotropin-releasing hormone III (GnRH3) is considered to be a key neurohormone in fish reproduction control. In the present study, the cDNA and genomic sequences of GnRH3 were cloned and characterized from large yellow croaker Larimichthys crocea. The cDNA encoded a protein of 99 amino acids with four functional motifs. The full-length genome sequence was composed of 3797 nucleotides, including four exons and three introns. Higher identities of amino acid sequences and conserved exon-intron organizations were found between LcGnRH3 and other GnRH3 genes. In addition, some special features of the sequences were detected in partial species. For example, two specific residues (V and A) were found in the family Sciaenidae, and the unique 75-72 bp type of the open reading frame 2 and 3 existed in the family Cyprinidae. Analysis of the 2576 bp promoter fragment of LcGnRH3 showed a number of transcription factor binding sites, such as AP1, CREB, GATA-1, HSF, FOXA2, and FOXL1. Promoter functional analysis using an EGFP reporter fusion in zebrafish larvae presented positive signals in the brain, including the olfactory region, the terminal nerve ganglion, the telencephalon, and the hypothalamus. The expression pattern was generally consistent with the endogenous GnRH3 GFP-expressing transgenic zebrafish lines, but the details were different. These results indicate that the structure and function of LcGnRH3 are generally similar to the other teleost GnRH3 genes, but there exist some distinctions among them.

Pleiotropic function of SRY-related HMG box transcription factor 4 in regulation of tumorigenesis

In addition to their critical roles in embryonic development, cell fate decision, and differentiation, members of Sox (Sry-related high-mobility group box) family of transcription factors including Sox4 have been implicated in various cancers. Multiple studies have revealed an increased expression along with specific oncogenic function of Sox4 in tumors, while others observed a reduced expression of Sox4 in different types of malignancies and suppression of tumor initiation or progression by this protein. More interestingly, the prognostic value of Sox4 is debated due to obvious differences between various reports as well as inconsistencies within specific studies. This review summarizes our current understanding of Sox4 expression pattern and its transcription-dependent, as well as transcription-independent, functions in tumor initiation or progression and its correlation with patient survival. We also discuss the existing discrepancies between different reports and their possible explanations.

Ontogenesis of gonadotropin-releasing hormone neurons: a model for hypothalamic neuroendocrine cell development

The vertebrate hypothalamo-pituitary-gonadal axis is the anatomical framework responsible for reproductive competence and species propagation. Essential to the coordinated actions of this three-tiered biological system is the fact that the regulatory inputs ultimately converge on the gonadotropin-releasing hormone (GnRH) neuronal system, which in rodents primarily resides in the preoptic/hypothalamic region. In this short review we will focus on: (1) the general embryonic temporal and spatial development of the rodent GnRH neuronal system, (2) the origin(s) of GnRH neurons, and (3) which transcription - and growth factors have been found to be critical for GnRH neuronal ontogenesis and cellular fate-specification. Moreover, we ask the question whether the molecular and cellular mechanisms involved in GnRH neuronal development may also play a role in the development of other hypophyseal secreting neuroendocrine cells in the hypothalamus.

The role of SRY-related HMG box transcription factor 4 (SOX4) in tumorigenesis and metastasis: friend or foe?

Development and progression of cancer are mediated by alterations in transcriptional networks, resulting in a disturbed balance between the activity of oncogenes and tumor suppressor genes. Transcription factors have the capacity to regulate global transcriptional profiles, and are consequently often found to be deregulated in their expression and function during tumorigenesis. Sex-determining region Y-related high-mobility-group box transcription factor 4 (SOX4) is a member of the group C subfamily of the SOX transcription factors and has a critical role during embryogenesis, where its expression is widespread and controls the development of numerous tissues. SOX4 expression is elevated in a wide variety of tumors, including leukemia, colorectal cancer, lung cancer and breast cancer, suggesting a fundamental role in the development of these malignancies. In many cancers, deregulated expression of this developmental factor has been correlated with increased cancer cell proliferation, cell survival, inhibition of apoptosis and tumor progression through the induction of an epithelial-to-mesenchymal transition and metastasis. However, in a limited subset of tumors, SOX4 has also been reported to act as a tumor suppressor. These opposing roles suggest that the outcome of SOX4 activation depends on the cellular context and the tumor origin. Indeed, SOX4 expression, transcriptional activity and target gene specificity can be controlled by signaling pathways, including the transforming growth factor-β and the WNT pathway, as well as at the post-translational level through regulation of protein stability and interaction with specific cofactors, such as TCF, syntenin-1 and p53. Here, we provide an overview of our current knowledge concerning the role of SOX4 in tumor development and progression.

SOX4 transcriptionally regulates multiple SEMA3/plexin family members and promotes tumor growth in pancreatic cancer

Semaphorin signaling through Plexin frequently participates in tumorigenesis and malignant progression in various types of cancer. In particular, the role of semaphorin signaling in pancreatic ductal adenocarcinoma (PDAC) remains unexplored, despite a high likelihood of metastasis and mortality. Unlike other epithelial malignancies that often express a small number of specific genes in the Semaphorin/Plexin family, five or more are often expressed in human PDAC. Such concomitant expression of these SEMA3/Plexin family members is not a result of gene amplification, but (at least partially) from increased gene transcription activated by SOX4 de novo expressed in PDAC. Via chromatin-immunoprecipitation, luciferase promoter activity assay and electrophoresis mobility shift assay, SOX4 is demonstrated to bind to the consensus site at the promoter of each SEMA3 and Plexin gene to enhance transcription activity. Conversely, RNAi-knockdown of SOX4 in PDAC cell lines results in decreased expression of SEMA3/Plexin family members and is associated with restricted tumor growth both in vitro and in SCID mice. We further demonstrate that SOX4 levels parallel with the summed expression of SEMA3/Plexin family members (P = 0.033, NPar Kruskal-Wallis one-way analysis), which also correlates with poor survival in human PDAC (P = 0.0409, Kaplan-Meier analysis). Intriguingly, miR-129-2 and miR-335, both of which target SOX4 for degradation, are co-repressed in human PDAC cases associated with up-regulated SOX4 in a statistically significant way. In conclusion, we disclose a miR-129-2(miR-335)/SOX4/Semaphorin-Plexin regulatory axis in the tumorigenesis of pancreatic cancer.