Homeodomain Proteins SIX3 and SIX6 Regulate Gonadotrope-specific Genes During Pituitary Development

LncRNA12097.1 contributes to endometrial cell growth by enhancing YES1 activating β-catenin via sponging miR-145-5p

Despite previous investigations elucidating the regulatory mechanisms of long non-coding RNAs (lncRNAs) in endometrial function and reproductive disorders, the precise pathways through which lncRNAs impact endometrial functions and fertility remain unclear. In this study, we performed an expression profile analysis of lncRNAs in the endometrial tissue of Hu sheep with different prolificacy, identifying 13,707 lncRNAs. We discovered a bidirectional lncRNA, designated lncRNA12097.1, exhibiting significant up-regulation exclusively in the endometrium of Hu sheep with high fecundity. Functional analyses revealed lncRNA12097.1 significantly enhanced proliferation and cell cycle progression in both endometrial epithelial cell (EEC) and stromal cells (ESC), while inhibiting apoptosis in these cell types. Mechanistically, we demonstrated a directly interaction between lncRNA12097.1 and miR-145-5p, with YES proto-oncogene 1 (YES1) being identified as a validated target of miR-145-5p. Interference with lncRNA12097.1 resulted in suppressed cell growth through down-regulation of YES1 expression, which could be rescued by miR-145-5p. Furthermore, lncRNA12097.1 functions as a competitive endogenous RNA (ceRNA) for miR-145-5p in ESCs, sequestering miR-145-5p and preventing its binding to the 3'UTR of YES1 mRNA. This interaction led to increased expression of YES1 and subsequent activation of downstream β-catenin signaling, thereby promoting ESC growth in Hu sheep. These findings provide novel molecular insights into the mechanisms underlying prolificacy in sheep.

Reduced nephron endowment in the common Six2-TGC tg mouse line is due to Six3 misexpression by aberrant enhancer-promoter interactions in the transgene

Lifelong kidney function relies on the complement of nephrons generated during mammalian development from a mesenchymal nephron progenitor cell (NPC) population. Low nephron endowment confers increased susceptibility to chronic kidney disease. We asked whether reduced nephron numbers in the popular Six2TGC transgenic mouse line 1 was due to disruption of a regulatory gene at the integration site or to ectopic expression of a gene(s) contained within the transgene. Targeted locus amplification identified integration of the Six2TGC transgene within an intron of Cntnap5a on chr1. We generated Hi-C datasets from NPCs isolated from the Six2TGC tg/+ mice, the Cited1 CreERT2/+ control mice, and the Six2TGC tg/+ ; Tsc1 +/Flox,2 mice that exhibited restored nephron number compared with Six2TGC tg/+ mice, and mapped the precise integration of Six2TGC and Cited1 CreERT2 transgenes to chr1 and chr14, respectively. No changes in topology, accessibility, or expression were observed within the 50-megabase region centered on Cntnap5a in Six2TGC tg/+ mice compared with control mice. By contrast, we identified an aberrant regulatory interaction between a Six2 distal enhancer and the Six3 promoter contained within the transgene. Increasing the Six2TGC tg to Six2 locus ratio or removing one Six2 allele in Six2TGC tg/+ mice, caused severe renal hypoplasia. Furthermore, CRISPR disruption of Six3 within the transgene ( Six2TGC Δ Six3CT ) restored nephron endowment to wildtype levels and abolished the stoichiometric effect. Data from genetic and biochemical studies together suggest that in Six2TGC, SIX3 interferes with SIX2 function in NPC renewal through its C-terminal domain.

Significance

Using high-resolution chromatin conformation and accessibility datasets we mapped the integration site of two popular transgenes used in studies of nephron progenitor cells and kidney development. Aberrant enhancer-promoter interactions drive ectopic expression of Six3 in the Six2TGC tg line which was correlated with disruption of nephrogenesis. Disruption of Six3 within the transgene restored nephron numbers to control levels; further genetic and biochemical studies suggest that Six3 interferes with Six2 -mediated regulation of NPC renewal.

Comprehensive analysis of differences in N6-methyladenosine RNA methylomes in the rat adenohypophysis after GnRH treatment

N6-methyladenosine is considered to be the most common and abundant internal chemical modification among the more than 150 identified chemical RNA modifications. It is involved in most biological processes and actively participates in the regulation of animal reproduction. However, the potential function of m⁶ A in the pituitaries of mammals is not yet clear. It is also unknown whether m⁶ A is involved in the secretion and regulation of FSH by GnRH, which in turn affects mammalian reproduction. In this study, rats were treated with gonadorelin to simulate physiological GnRH-mediated regulation of FSH synthesis and secretion, and m⁶ A-seq was used to analyze the differential m⁶ A modification of the rat pituitary after gonadorelin treatment. A whole-transcriptome map of m⁶ A in the rat pituitary gland before and after gonadorelin treatment was successfully created. A total of 6413 differential peaks were identified, of which 3764 m⁶ A peaks were upregulated and 2649 m⁶ A peaks were downregulated. Among the 709 differentially expressed genes, 250 genes were discovered with differential methylation modifications. Intriguingly, the altered m⁶ A peaks within mRNAs were enriched in steroid biosynthetic processes and responses to cAMP. The results of the study will lay a foundation for further exploration of the potential role of m⁶ A modification in the regulation of reproductive hormone secretion and provide a theoretical basis for the application of GnRH analogs in mammalian artificial reproduction.

Association of SIX1-SIX6 polymorphisms with peripapillary retinal nerve fibre layer thickness in children

PURPOSE

Association of SIX1-SIX6 variants with peripapillary retinal nerve fibre layer (p-RNFL) thickness had been reported in adults. This study aimed to investigate these associations in children, with further explorations by spatial, age and sex stratifications.

METHODS

2878 school children aged between 6 and 9 years were enrolled from the Hong Kong Children Eye Study. Three single-nucleotide polymorphisms (SNPs) at the SIX1-SIX6 locus were genotyped. The association of each SNP with p-RNFL thickness (including global and sectoral thickness) were evaluated using multiple linear regression.

RESULTS

SNPs rs33912345 (p=7.7×10^-4) and rs10483727 (p=0.0013) showed significant associations with temporal-inferior p-RNFL thickness. The C allele of rs33912345 was associated with a thinner temporal-inferior p-RNFL by an average of 2.44 µm, while rs10483727-T was associated with a thinner temporal-inferior p-RNFL by 2.32 µm. The association with temporal-inferior p-RNFL was the strongest in the 8-9 year-old group for rs33912345 (p=5.2×10^-4) and rs10483727 (p=3.3×10^-4). Both SNPs were significantly associated with temporal-inferior p-RNFL thickness in boys (p<0.0017), but not in girls (p>0.05). In contrast, rs12436579-C (β=1.66; p=0.0059), but not rs33912345-C (β=1.31; p=0.052) or rs10483727-T (β=1.19; p=0.078), was nominally associated with a thicker nasal-inferior p-RNFL.

CONCLUSIONS

Both rs33912345 and rs10483727 at SIX1-SIX6 were associated with p-RNFL thickness in children, especially at the temporal-inferior sector, with age-dependent and sex-specific effects. SNP rs12436579 was associated with nasal-inferior p-RNFL thickness. Our findings suggested a role of SIX1-SIX6 in RNFL variation during neural retina development in childhood.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

The SIX Family of Transcription Factors: Common Themes Integrating Developmental and Cancer Biology

The sine oculis (SIX) family of transcription factors are key regulators of developmental processes during embryogenesis. Members of this family control gene expression to promote self-renewal of progenitor cell populations and govern mechanisms of cell differentiation. When the function of SIX genes becomes disrupted, distinct congenital defects develops both in animal models and humans. In addition to the embryonic setting, members of the SIX family have been found to be critical regulators of tumorigenesis, promoting cell proliferation, epithelial-to-mesenchymal transition, and metastasis. Research in both the fields of developmental biology and cancer research have provided an extensive understanding of SIX family transcription factor functions. Here we review recent progress in elucidating the role of SIX family genes in congenital disease as well as in the promotion of cancer. Common themes arise when comparing SIX transcription factor function during embryonic and cancer development. We highlight the complementary nature of these two fields and how knowledge in one area can open new aspects of experimentation in the other.

Physiological Characterization and Transcriptomic Properties of GnRH Neurons Derived From Human Stem Cells

Gonadotropin-releasing hormone (GnRH) neurons in the hypothalamus play a key role in the regulation of reproductive function. In this study, we sought an efficient method for generating GnRH neurons from human embryonic and induced pluripotent stem cells (hESC and hiPSC, respectively). First, we found that exposure of primitive neuroepithelial cells, rather than neuroprogenitor cells, to fibroblast growth factor 8 (FGF8), was more effective in generating GnRH neurons. Second, addition of kisspeptin to FGF8 further increased the efficiency rates of GnRH neurogeneration. Third, we generated a fluorescent marker mCherry labeled human embryonic GnRH cell line (mCh-hESC) using a CRISPR-Cas9 targeting approach. Fourth, we examined physiological characteristics of GnRH (mCh-hESC) neurons: similar to GnRH neurons in vivo, they released the GnRH peptide in a pulsatile manner at ~60 min intervals; GnRH release increased in response to high potassium, kisspeptin, estradiol, and neurokinin B challenges; and injection of depolarizing current induced action potentials. Finally, we characterized developmental changes in transcriptomes of GnRH neurons using hESC, hiPSC, and mCh-hESC. The developmental pattern of transcriptomes was remarkably similar among the 3 cell lines. Collectively, human stem cell-derived GnRH neurons will be an important tool for establishing disease models to understand diseases, such as idiopathic hypothalamic hypogonadism, and testing contraceptive drugs.

Evidence That the Etiology of Congenital Hypopituitarism Has a Major Genetic Component but Is Infrequently Monogenic

Purpose

Congenital hypopituitarism usually occurs sporadically. In most patients, the etiology remains unknown.

Methods

We studied 13 children with sporadic congenital hypopituitarism. Children with non-endocrine, non-familial idiopathic short stature (NFSS) (n = 19) served as a control group. Exome sequencing was performed in probands and both unaffected parents. A burden testing approach was used to compare the number of candidate variants in the two groups.

Results

First, we assessed the frequency of rare, predicted-pathogenic variants in 42 genes previously reported to be associated with pituitary gland development. The average number of variants per individual was greater in probands with congenital hypopituitarism than those with NFSS (1.1 vs. 0.21, mean variants/proband, P = 0.03). The number of probands with at least 1 variant in a pituitary-associated gene was greater in congenital hypopituitarism than in NFSS (62% vs. 21%, P = 0.03). Second, we assessed the frequency of rare, predicted-pathogenic variants in the exome (to capture undiscovered causes) that were inherited in a fashion that could explain the sporadic occurrence of the proband's condition with a monogenic etiology (de novo mutation, autosomal recessive, or X-linked recessive) with complete penetrance. There were fewer monogenic candidates in the probands with congenital hypopituitarism than those with NFSS (1.3 vs. 2.5 candidate variants/proband, P = 0.024). We did not find any candidate variants (0 of 13 probands) in genes previously reported to explain the phenotype in congenital hypopituitarism, unlike NFSS (8 of 19 probands, P = 0.01).

Conclusion

Our findings provide evidence that the etiology of sporadic congenital hypopituitarism has a major genetic component but may be infrequently monogenic with full penetrance, suggesting a more complex etiology.

Androgen receptor positively regulates gonadotropin-releasing hormone receptor in pituitary gonadotropes

Within pituitary gonadotropes, the gonadotropin-releasing hormone receptor (GnRHR) receives hypothalamic input from GnRH neurons that is critical for reproduction. Previous studies have suggested that androgens may regulate GnRHR, although the mechanisms remain unknown. In this study, we demonstrated that androgens positively regulate Gnrhr mRNA in mice. We then investigated the effects of androgens and androgen receptor (AR) on Gnrhr promoter activity in immortalized mouse LβT2 cells, which represent mature gonadotropes. We found that AR positively regulates the Gnrhr proximal promoter, and that this effect requires a hormone response element (HRE) half site at -159/-153 relative to the transcription start site. We also identified nonconsensus, full-length HREs at -499/-484 and -159/-144, which are both positively regulated by androgens on a heterologous promoter. Furthermore, AR associates with the Gnrhr promoter in ChIP. Altogether, we report that GnRHR is positively regulated by androgens through recruitment of AR to the Gnrhr proximal promoter.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

SIX6 is a TAL1-regulated transcription factor in T-ALL and associated with inferior outcome

T-cell acute lymphoblastic leukemia (T-ALL) is a hematological malignancy driven by abnormal activity of transcription factors. Here we report an aberrant expression of the developmental transcription factor SIX6 in the TAL1-subtype of T-ALL. Our results demonstrate that the binding of TAL1 and GATA3 transcription factors into an upstream enhancer element directly regulates SIX6 expression. High expression of SIX6 was associated with inferior event-free survival within three independent patient cohorts. At a functional level, CRISPR-Cas9-mediated knockout of the SIX6 gene in TAL1 positive Jurkat cells induced changes in genes associated with the mTOR-, K-RAS-, and TNFα-related molecular signatures but did not impair cell proliferation or viability. There was also no acceleration of T-ALL development within a Myc driven zebrafish tumor model in vivo. Taken together, our results show that SIX6 belongs to the TAL1 regulatory gene network in T-ALL but is alone insufficient to influence the development or maintenance of T-ALL.

Copy Number Variants Contributing to Combined Pituitary Hormone Deficiency

Combined pituitary hormone deficiency represents a disorder with complex etiology. For many patients, causes of the disease remain unexplained, despite usage of advanced genetic testing. Although major and common transcription factors were identified two decades ago, we still struggle with identification of rare inborn factors contributing to pituitary function. In this report, we follow up genomic screening of CPHD patient cohort that were previously tested for changes in a coding sequences of genes with the use of the whole exome. We aimed to find contribution of rare copy number variations (CNVs). As a result, we identified genomic imbalances in 7 regions among 12 CPHD patients. Five out of seven regions showed copy gains whereas two presented losses of genomic fragment. Three regions with detected gains encompassed known CPHD genes namely LHX4, HESX1, and OTX2. Among new CPHD loci, the most interesting seem to be the region covering SIX3 gene, that is abundantly expressed in developing brain, and together with HESX1 contributes to pituitary organogenesis as it was evidenced before in functional studies. In conclusion, with the use of broadened genomic approach we identified copy number imbalances for 12 CPHD patients. Although further functional studies are required in order to estimate its true impact on expression pattern during pituitary organogenesis and CPHD etiology.

Augmentation of Nr4a3 and Suppression of Fshb Expression in the Pituitary Gland of Female Annexin A5 Null Mouse

GnRH enhances the expression of annexin A5 (ANXA5) in pituitary gonadotropes, and ANXA5 enhances gonadotropin secretion. However, the impact of ANXA5 regulation on the expression of pituitary hormone genes remains unclear. Here, using quantitative PCR, we demonstrated that ANXA5 deficiency in female mice reduced the expression of Fshb and Gh in their pituitary glands. Transcriptome analysis confirmed a specific increase in Nr4a3 mRNA expression in addition to lower levels of Fshb expression in ANXA5-deficient female pituitary glands. This gene was then found to be a GnRH-inducible immediate early gene, and its increased expression caused protein to accumulate in the nucleus after administration of a GnRH agonist in LβT2 cells, which are an in vitro pituitary gonadotrope model. The increase in ANXA5 protein levels in LβT2 cells clearly suppressed Nr4a3 expression. siRNA-mediated inhibition of Nr4a3 expression increased Fshb expression. The results revealed that GnRH stimulates Nr4a3 and Anxa5 sequentially. NR4A3 suppression of Fshb may be necessary for later massive secretion of FSH by GnRH in gonadotropes, and Nr4a3 would be negatively regulated by ANXA5 to increase FSH secretion.

Comprehensive Analysis of Differentially Expressed Profiles of mRNA, lncRNA, and circRNA in the Uterus of Seasonal Reproduction Sheep

Photoperiod is one of the important factors leading to seasonal reproduction of sheep. However, the molecular mechanisms underlying the photoperiod regulation of seasonal reproduction remain poorly understood. In this study, we compared the expression profiles of mRNAs, lncRNAs, and circRNAs in uterine tissues from Sunite sheep during three different photoperiods, namely, the short photoperiod (SP), short transfer to long photoperiod (SLP), and long photoperiod (LP). The results showed that 298, 403, and 378 differentially expressed (DE) mRNAs, 171, 491, and 499 DE lncRNAs, and 124, 270, and 400 DE circRNAs were identified between SP and LP, between SP and SLP, and between LP and SLP, respectively. Furthermore, functional enrichment analysis showed that the differentially expressed RNAs were mainly involved in the GnRH signaling pathway, thyroid hormone synthesis, and thyroid hormone signaling pathway. In addition, co-expression networks of lncRNA–mRNA were constructed based on the correlation analysis between the differentially expressed RNAs. Our study provides new insights into the expression changes of RNAs in different photoperiods, which might contribute to understanding the molecular mechanisms of seasonal reproduction in sheep.

Development of the Pituitary Gland

The development of the anterior pituitary gland occurs in distinct sequential developmental steps, leading to the formation of a complex organ containing five different cell types secreting six different hormones. During this process, the temporal and spatial expression of a cascade of signaling molecules and transcription factors plays a crucial role in organ commitment, cell proliferation, patterning, and terminal differentiation. The morphogenesis of the gland and the emergence of distinct cell types from a common primordium are governed by complex regulatory networks involving transcription factors and signaling molecules that may be either intrinsic to the developing pituitary or extrinsic, originating from the ventral diencephalon, the oral ectoderm, and the surrounding mesenchyme. Endocrine cells of the pituitary gland are organized into structural and functional networks that contribute to the coordinated response of endocrine cells to stimuli; these cellular networks are formed during embryonic development and are maintained or may be modified in adulthood, contributing to the plasticity of the gland. Abnormalities in any of the steps of pituitary development may lead to congenital hypopituitarism that includes a spectrum of disorders from isolated to combined hormone deficiencies including syndromic disorders such as septo-optic dysplasia. Over the past decade, the acceleration of next-generation sequencing has allowed for rapid analysis of the patient genome to identify novel mutations and novel candidate genes associated with hypothalmo-pituitary development. Subsequent functional analysis using patient fibroblast cells, and the generation of stem cells derived from patient cells, is fast replacing the need for animal models while providing a more physiologically relevant characterization of novel mutations. Furthermore, CRISPR-Cas9 as the method for gene editing is replacing previous laborious and time-consuming gene editing methods that were commonly used, thus yielding knockout cell lines in a fraction of the time. © 2020 American Physiological Society. Compr Physiol 10:389-413, 2020.

Functional Annotation of the Transcriptome of the Pig, Sus scrofa, Based Upon Network Analysis of an RNAseq Transcriptional Atlas

The domestic pig (Sus scrofa) is both an economically important livestock species and a model for biomedical research. Two highly contiguous pig reference genomes have recently been released. To support functional annotation of the pig genomes and comparative analysis with large human transcriptomic data sets, we aimed to create a pig gene expression atlas. To achieve this objective, we extended a previous approach developed for the chicken. We downloaded RNAseq data sets from public repositories, down-sampled to a common depth, and quantified expression against a reference transcriptome using the mRNA quantitation tool, Kallisto. We then used the network analysis tool Graphia to identify clusters of transcripts that were coexpressed across the merged data set. Consistent with the principle of guilt-by-association, we identified coexpression clusters that were highly tissue or cell-type restricted and contained transcription factors that have previously been implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes, such as the cell cycle and oxidative phosphorylation. The same approach was used to identify coexpression clusters within RNAseq data from multiple individual liver and brain samples, highlighting cell type, process, and region-specific gene expression. Evidence of conserved expression can add confidence to assignment of orthology between pig and human genes. Many transcripts currently identified as novel genes with ENSSSCG or LOC IDs were found to be coexpressed with annotated neighbouring transcripts in the same orientation, indicating they may be products of the same transcriptional unit. The meta-analytic approach to utilising public RNAseq data is extendable to include new data sets and new species and provides a framework to support the Functional Annotation of Animals Genomes (FAANG) initiative.

Dysregulation of Hypothalamic Gene Expression and the Oxytocinergic System by Soybean Oil Diets in Male Mice

Soybean oil consumption has increased greatly in the past half-century and is linked to obesity and diabetes. To test the hypothesis that soybean oil diet alters hypothalamic gene expression in conjunction with metabolic phenotype, we performed RNA sequencing analysis using male mice fed isocaloric, high-fat diets based on conventional soybean oil (high in linoleic acid, LA), a genetically modified, low-LA soybean oil (Plenish), and coconut oil (high in saturated fat, containing no LA). The 2 soybean oil diets had similar but nonidentical effects on the hypothalamic transcriptome, whereas the coconut oil diet had a negligible effect compared to a low-fat control diet. Dysregulated genes were associated with inflammation, neuroendocrine, neurochemical, and insulin signaling. Oxt was the only gene with metabolic, inflammation, and neurological relevance upregulated by both soybean oil diets compared to both control diets. Oxytocin immunoreactivity in the supraoptic and paraventricular nuclei of the hypothalamus was reduced, whereas plasma oxytocin and hypothalamic Oxt were increased. These central and peripheral effects of soybean oil diets were correlated with glucose intolerance but not body weight. Alterations in hypothalamic Oxt and plasma oxytocin were not observed in the coconut oil diet enriched in stigmasterol, a phytosterol found in soybean oil. We postulate that neither stigmasterol nor LA is responsible for effects of soybean oil diets on oxytocin and that Oxt messenger RNA levels could be associated with the diabetic state. Given the ubiquitous presence of soybean oil in the American diet, its observed effects on hypothalamic gene expression could have important public health ramifications.

Transcription Profiles of Age-at-Maturity-Associated Genes Suggest Cell Fate Commitment Regulation as a Key Factor in the Atlantic Salmon Maturation Process

Despite recent taxonomic diversification in studies linking genotype with phenotype, follow-up studies aimed at understanding the molecular processes of such genotype-phenotype associations remain rare. The age at which an individual reaches sexual maturity is an important fitness trait in many wild species. However, the molecular mechanisms regulating maturation timing processes remain obscure. A recent genome-wide association study in Atlantic salmon (Salmo salar) identified large-effect age-at-maturity-associated chromosomal regions including genes vgll3, akap11 and six6, which have roles in adipogenesis, spermatogenesis and the hypothalamic-pituitary-gonadal (HPG) axis, respectively. Here, we determine expression patterns of these genes during salmon development and their potential molecular partners and pathways. Using Nanostring transcription profiling technology, we show development- and tissue-specific mRNA expression patterns for vgll3, akap11 and six6. Correlated expression levels of vgll3 and akap11, which have adjacent chromosomal location, suggests they may have shared regulation. Further, vgll3 correlating with arhgap6 and yap1, and akap11 with lats1 and yap1 suggests that Vgll3 and Akap11 take part in actin cytoskeleton regulation. Tissue-specific expression results indicate that vgll3 and akap11 paralogs have sex-dependent expression patterns in gonads. Moreover, six6 correlating with slc38a6 and rtn1, and Hippo signaling genes suggests that Six6 could have a broader role in the HPG neuroendrocrine and cell fate commitment regulation, respectively. We conclude that Vgll3, Akap11 and Six6 may influence Atlantic salmon maturation timing via affecting adipogenesis and gametogenesis by regulating cell fate commitment and the HPG axis. These results may help to unravel general molecular mechanisms behind maturation.

Cell Type- and Sex-Dependent Transcriptome Profiles of Rat Anterior Pituitary Cells

Understanding the physiology and pathology of an organ composed of a variety of cell populations depends critically on genome-wide information on each cell type. Here, we report single-cell transcriptome profiling of over 6,800 freshly dispersed anterior pituitary cells from postpubertal male and female rats. Six pituitary-specific cell types were identified based on known marker genes and characterized: folliculostellate cells and hormone-producing corticotrophs, gonadotrophs, thyrotrophs, somatotrophs, and lactotrophs. Also identified were endothelial and blood cells from the pituitary capillary network. The expression of numerous developmental and neuroendocrine marker genes in both folliculostellate and hormone-producing cells supports that they have a common origin. For several genes, the validity of transcriptome analysis was confirmed by qRT-PCR and single cell immunocytochemistry. Folliculostellate cells exhibit impressive transcriptome diversity, indicating their major roles in production of endogenous ligands and detoxification enzymes, and organization of extracellular matrix. Transcriptome profiles of hormone-producing cells also indicate contributions toward those functions, while also clearly demonstrating their endocrine function. This survey highlights many novel genetic markers contributing to pituitary cell type identity, sexual dimorphism, and function, and points to relationships between hormone-producing and folliculostellate cells.

Deletion of the Homeodomain Protein Six6 From GnRH Neurons Decreases GnRH Gene Expression, Resulting in Infertility

Hypothalamic GnRH (luteinizing hormone-releasing hormone) neurons are crucial for the hypothalamic-pituitary-gonadal (HPG) axis, which regulates mammalian fertility. Insufficient GnRH disrupts the HPG axis and is often associated with the genetic condition idiopathic hypogonadotropic hypogonadism (IHH). The homeodomain protein sine oculis-related homeobox 6 (Six6) is required for the development of GnRH neurons. Although it is known that Six6 is specifically expressed within a more mature GnRH neuronal cell line and that overexpression of Six6 induces GnRH transcription in these cells, the direct role of Six6 within the GnRH neuron in vivo is unknown. Here we find that global Six6 knockout (KO) embryos show apoptosis of GnRH neurons beginning at embryonic day 14.5 with 90% loss of GnRH neurons by postnatal day 1. We sought to determine whether the hypogonadism and infertility reported in the Six6KO mice are generated via actions within the GnRH neuron in vivo by creating a Six6-flox mouse and crossing it with the LHRHcre mouse. Loss of Six6 specifically within the GnRH neuron abolished GnRH expression in ∼0% of GnRH neurons. We further demonstrated that deletion of Six6 only within the GnRH neuron leads to infertility, hypogonadism, hypogonadotropism, and delayed puberty. We conclude that Six6 plays distinct roles in maintaining fertility in the GnRH neuron vs in the migratory environment of the GnRH neuron by maintaining expression of GnRH and survival of GnRH neurons, respectively. These results increase knowledge of the role of Six6 in the brain and may offer insight into the mechanism of IHH.

Identification of a pituitary ERα-activated enhancer triggering the expression of Nr5a1, the earliest gonadotrope lineage-specific transcription factor

Background

Gonadotrope lineage differentiation is a stepwise process taking place during pituitary development. The early step of gonadotrope lineage specification is characterized by the expression of the Nr5a1 transcription factor, a crucial factor for gonadotrope cell fate determination. Abnormalities affecting Nr5a1 expression lead to hypogonadotropic hypogonadism and infertility. Although significant knowledge has been gained on the signaling and transcriptional events controlling gonadotrope differentiation, epigenetic mechanisms regulating Nr5a1 expression during early gonadotrope lineage specification are still poorly understood.

Results

Using ATAC chromatin accessibility analyses on three cell lines recapitulating gradual stages of gonadotrope differentiation and in vivo on developing pituitaries, we demonstrate that a yet undescribed enhancer is transiently recruited during gonadotrope specification. Using CRISPR/Cas9, we show that this enhancer is mandatory for the emergence of Nr5a1 during gonadotrope specification. Furthermore, we identify a highly conserved estrogen-binding element and demonstrate that the enhancer activation is dependent upon estrogen acting through ERα. Lastly, we provide evidence that binding of ERα is crucial for chromatin remodeling of Nr5a1 enhancer and promoter, leading to RNA polymerase recruitment and transcription.

Conclusion

This study identifies the earliest regulatory sequence involved in gonadotrope lineage specification and highlights the key epigenetic role played by ERα in this differentiation process.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0291-8) contains supplementary material, which is available to authorized users.

Molecular regulation of follicle-stimulating hormone synthesis, secretion and action

Follicle-stimulating hormone (FSH) plays fundamental roles in male and female fertility. FSH is a heterodimeric glycoprotein expressed by gonadotrophs in the anterior pituitary. The hormone-specific FSHβ-subunit is non-covalently associated with the common α-subunit that is also present in the luteinizing hormone (LH), another gonadotrophic hormone secreted by gonadotrophs and thyroid-stimulating hormone (TSH) secreted by thyrotrophs. Several decades of research led to the purification, structural characterization and physiological regulation of FSH in a variety of species including humans. With the advent of molecular tools, availability of immortalized gonadotroph cell lines and genetically modified mouse models, our knowledge on molecular mechanisms of FSH regulation has tremendously expanded. Several key players that regulate FSH synthesis, sorting, secretion and action in gonads and extragonadal tissues have been identified in a physiological setting. Novel post-transcriptional and post-translational regulatory mechanisms have also been identified that provide additional layers of regulation mediating FSH homeostasis. Recombinant human FSH analogs hold promise for a variety of clinical applications, whereas blocking antibodies against FSH may prove efficacious for preventing age-dependent bone loss and adiposity. It is anticipated that several exciting new discoveries uncovering all aspects of FSH biology will soon be forthcoming.

WDR11-mediated Hedgehog signalling defects underlie a new ciliopathy related to Kallmann syndrome

WDR11 has been implicated in congenital hypogonadotropic hypogonadism (CHH) and Kallmann syndrome (KS), human developmental genetic disorders defined by delayed puberty and infertility. However, WDR11's role in development is poorly understood. Here, we report that WDR11 modulates the Hedgehog (Hh) signalling pathway and is essential for ciliogenesis. Disruption of WDR11 expression in mouse and zebrafish results in phenotypic characteristics associated with defective Hh signalling, accompanied by dysgenesis of ciliated tissues. Wdr11-null mice also exhibit early-onset obesity. We find that WDR11 shuttles from the cilium to the nucleus in response to Hh signalling. WDR11 regulates the proteolytic processing of GLI3 and cooperates with the transcription factor EMX1 in the induction of downstream Hh pathway gene expression and gonadotrophin-releasing hormone production. The CHH/KS-associated human mutations result in loss of function of WDR11. Treatment with the Hh agonist purmorphamine partially rescues the WDR11 haploinsufficiency phenotypes. Our study reveals a novel class of ciliopathy caused by WDR11 mutations and suggests that CHH/KS may be a part of the human ciliopathy spectrum.

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.

Chromatin status and transcription factor binding to gonadotropin promoters in gonadotrope cell lines

BACKGROUND

Proper expression of key reproductive hormones from gonadotrope cells of the pituitary is required for pubertal onset and reproduction. To further our understanding of the molecular events taking place during embryonic development, leading to expression of the glycoproteins luteinizing hormone (LH) and follicle-stimulating hormone (FSH), we characterized chromatin structure changes, imparted mainly by histone modifications, in model gonadotrope cell lines.

METHODS

We evaluated chromatin status and gene expression profiles by chromatin immunoprecipitation assays, DNase sensitivity assay, and RNA sequencing in three developmentally staged gonadotrope cell lines, αT1-1 (progenitor, expressing Cga), αT3-1 (immature, expressing Cga and Gnrhr), and LβT2 (mature, expressing Cga, Gnrhr, Lhb, and Fshb), to assess changes in chromatin status and transcription factor access of gonadotrope-specific genes.

RESULTS

We found the common mRNA α-subunit of LH and FSH, called Cga, to have an open chromatin conformation in all three cell lines. In contrast, chromatin status of Gnrhr is open only in αT3-1 and LβT2 cells. Lhb begins to open in LβT2 cells and was further opened by activin treatment. Histone H3 modifications associated with active chromatin were high on Gnrhr in αT3-1 and LβT2, and Lhb in LβT2 cells, while H3 modifications associated with repressed chromatin were low on Gnrhr, Lhb, and Fshb in LβT2 cells. Finally, chromatin status correlates with the progressive access of LHX3 to Cga and Gnrhr, followed by PITX1 binding to the Lhb promoter.

CONCLUSION

Our data show the gonadotrope-specific genes Cga, Gnrhr, Lhb, and Fshb are not only controlled by developmental transcription factors, but also by epigenetic mechanisms that include the modulation of chromatin structure, and histone modifications.

An ancient selective sweep linked to reproductive life history evolution in sockeye salmon

Study of parallel (or convergent) phenotypic evolution can provide important insights into processes driving sympatric, ecologically-mediated divergence and speciation, as ecotype pairs may provide a biological replicate of the underlying signals and mechanisms. Here, we provide evidence for a selective sweep creating an island of divergence associated with reproductive behavior in sockeye salmon (Oncorhynchus nerka), identifying a series of linked single nucleotide polymorphisms across a ~22,733 basepair region spanning the leucine-rich repeat-containing protein 9 gene exhibiting signatures of divergent selection associated with stream- and shore-spawning in both anadromous and resident forms across their pan-Pacific distribution. This divergence likely occurred ~3.8 Mya (95% HPD = 2.1–6.03 Mya), after sockeye separated from pink (O. gorbuscha) and chum (O. keta) salmon, but prior to the Pleistocene glaciations. Our results suggest recurrent evolution of reproductive ecotypes across the native range of O. nerka is at least partially associated with divergent selection of pre-existing genetic variation within or linked to this region. As sockeye salmon are unique among Pacific salmonids in their flexibility to spawn in lake-shore benthic environments, this region provides great promise for continued investigation of the genomic basis of O. nerka life history evolution, and, more broadly, for increasing our understanding of the heritable basis of adaptation of complex traits in novel environments.

Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post- versus prepubertal Brahman heifers

Puberty onset is a developmental process influenced by genetic determinants, environment, and nutrition. Mutations and regulatory gene networks constitute the molecular basis for the genetic determinants of puberty onset. The emerging knowledge of these genetic determinants presents opportunities for innovation in the breeding of early pubertal cattle. This paper presents new data on hypothalamic gene expression related to puberty in (Brahman) in age- and weight-matched heifers. Six postpubertal heifers were compared with 6 prepubertal heifers using whole-genome RNA sequencing methodology for quantification of global gene expression in the hypothalamus. Five transcription factors (TF) with potential regulatory roles in the hypothalamus were identified in this experiment: , , , , and . These TF genes were significantly differentially expressed in the hypothalamus of postpubertal versus prepubertal heifers and were also identified as significant according to the applied regulatory impact factor metric ( < 0.05). Two of these 5 TF, and , were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The gene belongs to the family of homologues of Drosophila sine oculis () genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as and are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in were associated with puberty in humans. Mutations in these TF, together with other genetic determinants previously discovered, could be used in genomic selection to predict the genetic merit of cattle (i.e., the likelihood of the offspring presenting earlier than average puberty for Brahman). Knowledge of key mutations involved in genetic traits is an advantage for genomic prediction because it can increase its accuracy.

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.

Epigenetic regulation of alternative promoters and enhancers in progenitor, immature, and mature gonadotrope cell lines

Gonadotrope cell identity genes emerge in a stepwise process during mouse pituitary development. Cga, encoding for the α-subunit of TSH, LH, and FSH, is initially detected at E11.5 followed by Gnrhr and steroidogenic factor Sf1 at E13.5, specifying cells engaged in a gonadotrope cell fate. Lhb and Fshb appear at E16.5 and 17.5, respectively, typifying differentiated gonadotrope cells. Using the αT1-1, αT3-1 and LβT2 cell lines recapitulating these stages of gonadotrope differentiation, DNA methylation at Gnrhr and Sf1 was investigated. Regulatory regions were found hypermethylated in progenitor αT1-1 cells and hypomethylated in differentiated LβT2 cells. Abundance of RNA polymerase II together with active histone modifications including H3K4me1, H3K4me3, and H3K27ac were strictly correlated with DNA hypomethylation. Analyses of epigenomic modifications and chromatin accessibility were further extended to Isl1, Lhx3, Gata2, and Pitx2, highlighting alternative usages of specific regulatory gene domains in progenitor αT1-1, immature αT3-1, and mature LβT2 gonadotrope cells.

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.

Polymorphisms and genes associated with puberty in heifers

Puberty onset is a multifactorial process influenced by genetic determinants and environmental conditions, especially nutritional status. Genes, genetic variations, and regulatory networks compose the molecular basis of achieving puberty. In this article, we reviewed the discovery of multiple polymorphisms and genes associated with heifer puberty phenotypes and discuss the opportunities to use this evolving knowledge of genetic determinants for breeding early pubertal Bos indicus-influenced cattle. The discovery of polymorphisms and genes was mainly achieved through candidate gene studies, quantitative trait loci analyses, genome-wide association studies, and recently, global gene expression studies (transcriptome). These studies are recapitulated and summarized in the current review.

Fshb-iCre mice are efficient and specific Cre deleters for the gonadotrope lineage

Genetic analysis of development and function of the gonadotrope cell lineage within mouse anterior pituitary has been greatly facilitated by at least three currently available Cre strains in which Cre was either knocked into the Gnrhr locus or expressed as a transgene from Cga and Lhb promoters. However, in each case there are some limitations including CRE expression in thyrotropes within pituitary or ectopic expression outside of pituitary, for example in some populations of neurons or gonads. Hence, these Cre strains often pose problems with regard to undesirable deletion of alleles in non-gonadotrope cells, fertility and germline transmission of mutant alleles. Here, we describe generation and characterization of a new Fshb-iCre deleter strain using 4.7 kb of ovine Fshb promoter regulatory sequences driving iCre expression exclusively in the gonadotrope lineage within anterior pituitary. Fshb-iCre mice develop normally, display no ectopic CRE expression in gonads and are fertile. When crossed onto a loxP recombination-mediated red to green color switch reporter mouse genetic background, in vivo CRE recombinase activity is detectable in gonadotropes at more than 95% efficiency and the GFP-tagged gonadotropes readily purified by fluorescence activated cell sorting. We demonstrate the applicability of this Fshb-iCre deleter strain in a mouse model in which Dicer is efficiently and selectively deleted in gonadotropes. We further show that loss of DICER-dependent miRNAs in gonadotropes leads to profound suppression of gonadotropins resulting in male and female infertility. Thus, Fshb-iCre mice serve as a new genetic tool to efficiently manipulate gonadotrope-specific gene expression in vivo.

Molecular insights into the aetiology of female reproductive ageing

As age at pubertal onset declines and age at first pregnancy increases, the mechanisms that regulate female reproductive lifespan become increasingly relevant to population health. The timing of menarche and menopause can have profound effects not only on fertility but also on the risk of diseases such as type 2 diabetes mellitus, cardiovascular disease and breast cancer. Genetic studies have identified dozens of highly penetrant rare mutations associated with reproductive disorders, and also ∼175 common genetic variants associated with the timing of puberty or menopause. These findings, alongside other functional studies, have highlighted a diverse range of mechanisms involved in reproductive ageing, implicating core biological processes such as cell cycle regulation and energy homeostasis. The aim of this article is to review the contribution of such genetic findings to our understanding of the molecular regulation of reproductive timing, as well as the biological basis of the epidemiological links between reproductive ageing and disease risk.