Human adult spermatogenesis balances spermatogonial stem cell (SSC) self-renewal and differentiation, alongside complex germ cell-niche interactions, to ensure long-term fertility and faithful genome propagation. Here, we performed single-cell RNA sequencing of ~6500 testicular cells from young adults. We found five niche/somatic cell types (Leydig, myoid, Sertoli, endothelial, macrophage), and observed germline-niche interactions and key human-mouse differences. Spermatogenesis, including meiosis, was reconstructed computationally, revealing sequential coding, non-coding, and repeat-element transcriptional signatures. Interestingly, we identified five discrete transcriptional/developmental spermatogonial states, including a novel early SSC state, termed State 0. Epigenetic features and nascent transcription analyses suggested developmental plasticity within spermatogonial States. To understand the origin of State 0, we profiled testicular cells from infants, and identified distinct similarities between adult State 0 and infant SSCs. Overall, our datasets describe key transcriptional and epigenetic signatures of the normal adult human testis, and provide new insights into germ cell developmental transitions and plasticity.

SummaryHuman adult spermatogonial stem cells (hSSCs) must balance self-renewal and differentiation. To understand how this is achieved, we profiled DNA methylation and open chromatin (ATAC-seq) in SSEA4+ hSSCs, analyzed bulk and single-cell RNA transcriptomes (RNA-seq) in SSEA4+ hSSCs and differentiating c-KIT+ spermatogonia, and performed validation studies via immunofluorescence. First, DNA hypomethylation at embryonic developmental genes supports their epigenetic “poising” in hSSCs for future/embryonic expression, while core pluripotency genes (OCT4 and NANOG) were transcriptionally and epigenetically repressed. Interestingly, open chromatin in hSSCs was strikingly enriched in binding sites for pioneer factors (NFYA/B, DMRT1, and hormone receptors). Remarkably, single-cell RNA-seq clustering analysis identified four cellular/developmental states during hSSC differentiation, involving major transitions in cell-cycle and transcriptional regulators, splicing and signaling factors, and glucose/mitochondria regulators. Overall, our results outline the dynamic chromatin/transcription landscape operating in hSSCs and identify crucial molecular pathways that accompany the transition from quiescence to proliferation and differentiation.

Bestrophin-1 is preferentially expressed at the basolateral membrane of the retinal pigmented epithelium (RPE) of the retina. Mutations in the BEST1 gene cause the retinal dystrophies vitelliform macular dystrophy, autosomal-dominant vitreochoroidopathy, and autosomal-recessive bestrophinopathy. Here, we describe four missense mutations in bestrophin-1, three that we believe are previously unreported, in patients diagnosed with autosomal-dominant and -recessive forms of retinitis pigmentosa (RP). The physiological function of bestrophin-1 remains poorly understood although its heterologous expression induces a Cl--specific current. We tested the effect of RP-causing variants on Cl- channel activity and cellular localization of bestrophin-1. Two (p.L140V and p.I205T) produced significantly decreased chloride-selective whole-cell currents in comparison to those of wild-type protein. In a model system of a polarized epithelium, two of three mutations (p.L140V and p.D228N) caused mislocalization of bestrophin-1 from the basolateral membrane to the cytoplasm. Mutations in bestrophin-1 are increasingly recognized as an important cause of inherited retinal dystrophy.

The RAS proto-oncogene Harvey rat sarcoma viral oncogene homolog (HRAS) encodes a small GTPase that transduces signals from cell surface receptors to intracellular effectors to control cellular behavior. Although somatic HRAS mutations have been described in many cancers, germline mutations cause Costello syndrome (CS), a congenital disorder associated with predisposition to malignancy. Based on the epidemiology of CS and the occurrence of HRAS mutations in spermatocytic seminoma, we proposed that activating HRAS mutations become enriched in sperm through a process akin to tumorigenesis, termed selfish spermatogonial selection. To test this hypothesis, we quantified the levels, in blood and sperm samples, of HRAS mutations at the p.G12 codon and compared the results to changes at the p.A11 codon, at which activating mutations do not occur. The data strongly support the role of selection in determining HRAS mutation levels in sperm, and hence the occurrence of CS, but we also found differences from the mutation pattern in tumorigenesis. First, the relative prevalence of mutations in sperm correlates weakly with their in vitro activating properties and occurrence in cancers. Second, specific tandem base substitutions (predominantly GC>TT/AA) occur in sperm but not in cancers; genomewide analysis showed that this same mutation is also overrepresented in constitutional pathogenic and polymorphic variants, suggesting a heightened vulnerability to these mutations in the germline. We developed a statistical model to show how both intrinsic mutation rate and selfish selection contribute to the mutational burden borne by the paternal germline.paternal age effect | male mutation bias | RASopathy | testis

Congenital hereditary endothelial dystrophy 1 (CHED1) and posterior polymorphous corneal dystrophy 1 (PPCD1) are autosomal-dominant corneal endothelial dystrophies that have been genetically mapped to overlapping loci on the short arm of chromosome 20. We combined genetic and genomic approaches to identify the cause of disease in extensive pedigrees comprising over 100 affected individuals. After exclusion of pathogenic coding, splice-site, and copy-number variations, a parallel approach using targeted and whole-genome sequencing facilitated the identification of pathogenic variants in a conserved region of the OVOL2 proximal promoter sequence in the index families (c.-339_361dup for CHED1 and c.-370T>C for PPCD1). Direct sequencing of the OVOL2 promoter in other unrelated affected individuals identified two additional mutations within the conserved proximal promoter sequence (c.-274T>G and c.-307T>C). OVOL2 encodes ovo-like zinc finger 2, a C2H2 zinc-finger transcription factor that regulates mesenchymal-to-epithelial transition and acts as a direct transcriptional repressor of the established PPCD-associated gene ZEB1. Interestingly, we did not detect OVOL2 expression in the normal corneal endothelium. Our in vitro data demonstrate that all four mutated OVOL2 promoters exhibited more transcriptional activity than the corresponding wild-type promoter, and we postulate that the mutations identified create cryptic cis-acting regulatory sequence binding sites that drive aberrant OVOL2 expression during endothelial cell development. Our data establish CHED1 and PPCD1 as allelic conditions and show that CHED1 represents the extreme of what can be considered a disease spectrum. They also implicate transcriptional dysregulation of OVOL2 as a common cause of dominantly inherited corneal endothelial dystrophies.

PurposeTo demonstrate the utility of an amplification-free long-read sequencing method to characterize the Fuchs endothelial corneal dystrophy (FECD)-associated intronic TCF4 triplet repeat (CTG18.1).MethodsWe applied an amplification-free method, utilizing the CRISPR/Cas9 system, in combination with PacBio single-molecule real-time (SMRT) long-read sequencing, to study CTG18.1. FECD patient samples displaying a diverse range of CTG18.1 allele lengths and zygosity status (n = 11) were analyzed. A robust data analysis pipeline was developed to effectively filter, align, and interrogate CTG18.1-specific reads. All results were compared with conventional polymerase chain reaction (PCR)-based fragment analysis.ResultsCRISPR-guided SMRT sequencing of CTG18.1 provided accurate genotyping information for all samples and phasing was possible for 18/22 alleles sequenced. Repeat length instability was observed for all expanded (≥50 repeats) phased CTG18.1 alleles analyzed. Furthermore, higher levels of repeat instability were associated with increased CTG18.1 allele length (mode length ≥91 repeats) indicating that expanded alleles behave dynamically.ConclusionCRISPR-guided SMRT sequencing of CTG18.1 has revealed novel insights into CTG18.1 length instability. Furthermore, this study provides a framework to improve the molecular diagnostic accuracy for CTG18.1-mediated FECD, which we anticipate will become increasingly important as gene-directed therapies are developed for this common age-related and sight threatening disease.

De novo point mutations arise predominantly in the male germline and increase in frequency with age, but it has not previously been possible to locate specific, identifiable mutations directly within the seminiferous tubules of human testes. Using microdissection of tubules exhibiting altered expression of the spermatogonial markers MAGEA4, FGFR3, and phospho-AKT, whole genome amplification, and DNA sequencing, we establish an in situ strategy for discovery and analysis of pathogenic de novo mutations. In 14 testes from men aged 39-90 y, we identified 11 distinct gain-of-function mutations in five genes (fibroblast growth factor receptors FGFR2 and FGFR3, tyrosine phosphatase PTPN11, and RAS oncogene homologs HRAS and KRAS) from 16 of 22 tubules analyzed; all mutations have known associations with severe diseases, ranging from congenital or perinatal lethal disorders to somatically acquired cancers. These results support proposed selfish selection of spermatogonial mutations affecting growth factor receptor-RAS signaling, highlight its prevalence in older men, and enable direct visualization of the microscopic anatomy of elongated mutant clones.iscerning the source of spontaneous germline mutations is fundamental to understanding the causes of many diseases, including monogenic developmental disorders (1) and complex conditions such as autism (2, 3) and schizophrenia (4). Recent whole genome sequencing studies of parent-child trios show that most mutations (such as nucleotide substitutions) originate from the paternal germline and increase in frequency with the father's age (5, 6), an issue of particular significance given the demographic shift to delayed reproduction in many populations (7). The deduction that the testes of older men harbor a greater burden of mutations, compared with younger men, is consistent with indirect measures of genetic decline, ranging from high indices of arrested germ cell divisions to complete involution of the seminiferous tubules (7-9). Surprisingly, however, it has not previously been possible to trace specific mutations back to their origins within individual germ cells (spermatogonia) of human testes.One mechanism proposed to contribute to the age-related increase in male mutations is selfish spermatogonial selection, a process equivalent to neoplasia but occurring in the unique context of the germ cell (10). In this process, specific point mutations that confer gain-of-function to components of the growth factor receptor-RAS signaling pathway occur rarely in spermatogonial stem cells of the adult testis but show a steep increase in prevalence with age, attributed to clonal expansion of mutant spermatogonia over time (11)(12)(13)(14)(15)(16). Fertilization of the egg by a mutant sperm leads to serious congenital disorders in the next generation, characterized by multiple malformations and, in some cases, a predisposition to malignancy. These disorders include Apert, Crouzon, and (16, 26). Based on the unexpectedly high birth prevalence of several of the associated congenital d...

The dominant congenital disorders Apert syndrome, achondroplasia and multiple endocrine neoplasia–caused by specific missense mutations in the FGFR2, FGFR3 and RET proteins respectively–represent classical examples of paternal age-effect mutation, a class that arises at particularly high frequencies in the sperm of older men. Previous analyses of DNA from randomly selected cadaveric testes showed that the levels of the corresponding FGFR2, FGFR3 and RET mutations exhibit very uneven spatial distributions, with localised hotspots surrounded by large mutation-negative areas. These studies imply that normal testes are mosaic for clusters of mutant cells: these clusters are predicted to have altered growth and signalling properties leading to their clonal expansion (selfish spermatogonial selection), but DNA extraction eliminates the possibility to study such processes at a tissue level. Using a panel of antibodies optimised for the detection of spermatocytic seminoma, a rare tumour of spermatogonial origin, we demonstrate that putative clonal events are frequent within normal testes of elderly men (mean age: 73.3 yrs) and can be classed into two broad categories. We found numerous small (less than 200 cells) cellular aggregations with distinct immunohistochemical characteristics, localised to a portion of the seminiferous tubule, which are of uncertain significance. However more infrequently we identified additional regions where entire seminiferous tubules had a circumferentially altered immunohistochemical appearance that extended through multiple serial sections that were physically contiguous (up to 1 mm in length), and exhibited enhanced staining for antibodies both to FGFR3 and a marker of downstream signal activation, pAKT. These findings support the concept that populations of spermatogonia in individual seminiferous tubules in the testes of older men are clonal mosaics with regard to their signalling properties and activation, thus fulfilling one of the specific predictions of selfish spermatogonial selection.