Transcriptome profiling of kisspeptin neurons from the mouse arcuate nucleus reveals new mechanisms in estrogenic control of fertility.

Kisspeptin neurons in the mediobasal hypothalamus (MBH) are critical targets of ovarian estrogen feedback regulating mammalian fertility. To reveal molecular mechanisms underlying this signaling, we thoroughly characterized the estrogen-regulated transcriptome of kisspeptin cells from ovariectomized transgenic mice substituted with 17ß-estradiol or vehicle. MBH kisspeptin neurons were harvested using laser-capture microdissection, pooled, and subjected to RNA sequencing. Estrogen treatment significantly (p.adj. < 0.05) up-regulated 1,190 and down-regulated 1,139 transcripts, including transcription factors, neuropeptides, ribosomal and mitochondrial proteins, ion channels, transporters, receptors, and regulatory RNAs. Reduced expression of the excitatory serotonin receptor-4 transcript (Htr4) diminished kisspeptin neuron responsiveness to serotonergic stimulation. Many estrogen-regulated transcripts have been implicated in puberty/fertility disorders. Patients (n = 337) with congenital hypogonadotropic hypogonadism (CHH) showed enrichment of rare variants in putative CHH-candidate genes (e.g., LRP1B, CACNA1G, FNDC3A). Comprehensive characterization of the estrogen-dependent kisspeptin neuron transcriptome sheds light on the molecular mechanisms of ovary-brain communication and informs genetic research on human fertility disorders.

CoverageMaster: comprehensive CNV detection and visualization from NGS short reads for genetic medicine applications.

CoverageMaster (CoM) is a copy number variation (CNV) calling algorithm based on depth-of-coverage maps designed to detect CNVs of any size in exome [whole exome sequencing (WES)] and genome [whole genome sequencing (WGS)] data. The core of the algorithm is the compression of sequencing coverage data in a multiscale Wavelet space and the analysis through an iterative Hidden Markov Model. CoM processes WES and WGS data at nucleotide scale resolution and accurately detects and visualizes full size range CNVs, including single or partial exon deletions and duplications. The results obtained with this approach support the possibility for coverage-based CNV callers to replace probe-based methods such as array comparative genomic hybridization and multiplex ligation-dependent probe amplification in the near future.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.