Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome.

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

FOXI3 pathogenic variants cause one form of craniofacial microsomia.

Craniofacial microsomia (CFM; also known as Goldenhar syndrome), is a craniofacial developmental disorder of variable expressivity and severity with a recognizable set of abnormalities. These birth defects are associated with structures derived from the first and second pharyngeal arches, can occur unilaterally and include ear dysplasia, microtia, preauricular tags and pits, facial asymmetry and other malformations. The inheritance pattern is controversial, and the molecular etiology of this syndrome is largely unknown. A total of 670 patients belonging to unrelated pedigrees with European and Chinese ancestry with CFM, are investigated. We identify 18 likely pathogenic variants in 21 probands (3.1%) in FOXI3. Biochemical experiments on transcriptional activity and subcellular localization of the likely pathogenic FOXI3 variants, and knock-in mouse studies strongly support the involvement of FOXI3 in CFM. Our findings indicate autosomal dominant inheritance with reduced penetrance, and/or autosomal recessive inheritance. The phenotypic expression of the FOXI3 variants is variable. The penetrance of the likely pathogenic variants in the seemingly dominant form is reduced, since a considerable number of such variants in affected individuals were inherited from non-affected parents. Here we provide suggestive evidence that common variation in the FOXI3 allele in trans with the pathogenic variant could modify the phenotypic severity and accounts for the incomplete penetrance.

Functional characteristics of a broad spectrum of TBX6 variants in Mayer-Rokitansky-Küster-Hauser syndrome.

PURPOSE: Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) is characterized by congenital absence of the uterus, cervix, and the upper part of the vagina in females. Whole-gene deletion and loss-of-function variants in TBX6 have been identified in association with MRKHS. We aimed to expand the spectrum of TBX6 variants in MRKHS and explore the biological effect of the variant alleles. METHODS: Rare variants in TBX6 were called from a combined multiethnic cohort of 622 probands with MRKHS who underwent exome sequencing or genome sequencing. Multiple in vitro functional experiments were performed, including messenger RNA analysis, western blotting, transcriptional activity assay, and immunofluorescence staining. RESULTS: We identified 16 rare variants in TBX6 from the combined cohort, including 1 protein-truncating variant reported in our previous study and 15 variants with unknown effects. By comparing the prevalence of TBX6 variants in the Chinese MRKHS cohort vs 1038 female controls, we observed a significant mutational burden of TBX6 in affected individuals (P = .0004, odds ratio = 5.25), suggesting a causal role of TBX6 variants in MRKHS. Of the 15 variants with uncertain effects, 7 were shown to induce a loss-of-function effect through various mechanisms. The c.423G>A (p.Leu141=) and c.839+5G>A variants impaired the normal splicing of TBX6 messenger RNA, c.422T>C (p.Leu141Pro) and c.745G>A (p.Val249Met) led to decreased protein expression, c.10C>T (p.Pro4Ser) and c.400G>A (p.Glu134Lys) resulted in perturbed transcriptional activity, and c.356G>A (p.Arg119His) caused protein mislocalization. We observed incomplete penetrance and variable expressivity in families carrying deleterious variants, which indicates a more complex genetic mechanism than classical Mendelian inheritance. CONCLUSION: Our study expands the mutational spectrum of TBX6 in MRKHS and delineates the molecular pathogenesis of TBX6 variants, supporting the association between deleterious variants in TBX6 and MRKHS.

Dominant monoallelic variant in the PAK2 gene causes Knobloch syndrome type 2.

Knobloch syndrome is an autosomal recessive phenotype mainly characterized by retinal detachment and encephalocele caused by biallelic pathogenic variants in the COL18A1 gene. However, there are patients clinically diagnosed as Knobloch syndrome with unknown molecular etiology not linked to COL18A1. We studied an historical pedigree (published in 1998) designated as KNO2 (Knobloch type 2 syndrome with intellectual disability, autistic behavior, retinal degeneration, encephalocele). Whole exome sequencing of the two affected siblings and the normal parents resulted in the identification of a PAK2 non-synonymous substitution p.(Glu435Lys) as a causative variant. The variant was monoallelic and apparently de novo in both siblings indicating a likely germ-line mosaicism in one of the parents; the mosaicism, however, could not be observed after deep sequencing of blood parental DNA. PAK2 encodes a member of a small group of serine/threonine kinases; these P21-activating kinases (PAKs) are essential in signal transduction and cellular regulation (cytoskeletal dynamics, cell motility, death and survival signaling and cell cycle progression). Structural analysis of the PAK2 p.(Glu435Lys) variant that is located in the kinase domain of the protein predicts a possible compromise in the kinase activity. Functional analysis of the p.(Glu435Lys) PAK2 variant in transfected HEK293T cells results in a partial loss of the kinase activity. PAK2 has been previously suggested as an autism-related gene. Our results show that PAK2-induced phenotypic spectrum is broad and not fully understood. We conclude that the KNO2 syndrome in the studied family is dominant and caused by a deleterious variant in the PAK2 gene.

Perturbations of genes essential for Müllerian duct and Wölffian duct development in Mayer-Rokitansky-Küster-Hauser syndrome.

Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) is associated with congenital absence of the uterus, cervix, and the upper part of the vagina; it is a sex-limited trait. Disrupted development of the Müllerian ducts (MD)/Wölffian ducts (WD) through multifactorial mechanisms has been proposed to underlie MRKHS. In this study, exome sequencing (ES) was performed on a Chinese discovery cohort (442 affected subjects and 941 female control subjects) and a replication MRKHS cohort (150 affected subjects of mixed ethnicity from North America, South America, and Europe). Phenotypic follow-up of the female reproductive system was performed on an additional cohort of PAX8-associated congenital hypothyroidism (CH) (n = 5, Chinese). By analyzing 19 candidate genes essential for MD/WD development, we identified 12 likely gene-disrupting (LGD) variants in 7 genes: PAX8 (n = 4), BMP4 (n = 2), BMP7 (n = 2), TBX6 (n = 1), HOXA10 (n = 1), EMX2 (n = 1), and WNT9B (n = 1), while LGD variants in these genes were not detected in control samples (p = 1.27E-06). Interestingly, a sex-limited penetrance with paternal inheritance was observed in multiple families. One additional PAX8 LGD variant from the replication cohort and two missense variants from both cohorts were revealed to cause loss-of-function of the protein. From the PAX8-associated CH cohort, we identified one individual presenting a syndromic condition characterized by CH and MRKHS (CH-MRKHS). Our study demonstrates the comprehensive utilization of knowledge from developmental biology toward elucidating genetic perturbations, i.e., rare pathogenic alleles involving the same loci, contributing to human birth defects.

Single cell transcriptome in aneuploidies reveals mechanisms of gene dosage imbalance.

Aneuploidy is a major source of gene dosage imbalance due to copy number alterations (CNA), and viable human trisomies are model disorders of altered gene expression. We study gene and allele-specific expression (ASE) of 9668 single-cell fibroblasts from trisomy 21 (T21) discordant twins and from mosaic T21, T18, T13 and T8. We examine 928 single cells with deep scRNAseq. Expected and observed overexpression of trisomic genes in trisomic vs. diploid bulk RNAseq is not detectable in trisomic vs. diploid single cells. Instead, for trisomic genes with low-to-average expression, their altered gene dosage is mainly due to the higher fraction of trisomic cells simultaneously expressing these genes, in agreement with a stochastic 2-state burst-like model of transcription. These results, confirmed in a further analysis of 8740 single fibroblasts with shallow scRNAseq, suggest that the specific transcriptional profile of each gene contributes to the phenotypic variability of trisomies. We propose an improved model to understand the effects of CNA and, generally, of gene regulation on gene dosage imbalance.

Multi-omic measurements of heterogeneity in HeLa cells across laboratories.

Reproducibility in research can be compromised by both biological and technical variation, but most of the focus is on removing the latter. Here we investigate the effects of biological variation in HeLa cell lines using a systems-wide approach. We determine the degree of molecular and phenotypic variability across 14 stock HeLa samples from 13 international laboratories. We cultured cells in uniform conditions and profiled genome-wide copy numbers, mRNAs, proteins and protein turnover rates in each cell line. We discovered substantial heterogeneity between HeLa variants, especially between lines of the CCL2 and Kyoto varieties, and observed progressive divergence within a specific cell line over 50 successive passages. Genomic variability has a complex, nonlinear effect on transcriptome, proteome and protein turnover profiles, and proteotype patterns explain the varying phenotypic response of different cell lines to Salmonella infection. These findings have implications for the interpretation and reproducibility of research results obtained from human cultured cells.

Extensive cellular heterogeneity of X inactivation revealed by single-cell allele-specific expression in human fibroblasts.

X-chromosome inactivation (XCI) provides a dosage compensation mechanism where, in each female cell, one of the two X chromosomes is randomly silenced. However, some genes on the inactive X chromosome and outside the pseudoautosomal regions escape from XCI and are expressed from both alleles (escapees). We investigated XCI at single-cell resolution combining deep single-cell RNA sequencing with whole-genome sequencing to examine allelic-specific expression in 935 primary fibroblast and 48 lymphoblastoid single cells from five female individuals. In this framework we integrated an original method to identify and exclude doublets of cells. In fibroblast cells, we have identified 55 genes as escapees including five undescribed escapee genes. Moreover, we observed that all genes exhibit a variable propensity to escape XCI in each cell and cell type and that each cell displays a distinct expression profile of the escapee genes. A metric, the Inactivation Score-defined as the mean of the allelic expression profiles of the escapees per cell-enables us to discover a heterogeneous and continuous degree of cellular XCI with extremes represented by "inactive" cells, i.e., cells exclusively expressing the escaping genes from the active X chromosome and "escaping" cells expressing the escapees from both alleles. We found that this effect is associated with cell-cycle phases and, independently, with the XIST expression level, which is higher in the quiescent phase (G0). Single-cell allele-specific expression is a powerful tool to identify novel escapees in different tissues and provide evidence of an unexpected cellular heterogeneity of XCI.

Deciphering Cell Lineage Specification during Male Sex Determination with Single-Cell RNA Sequencing.

The gonad is a unique biological system for studying cell-fate decisions. However, major questions remain regarding the identity of somatic progenitor cells and the transcriptional events driving cell differentiation. Using time-series single-cell RNA sequencing on XY mouse gonads during sex determination, we identified a single population of somatic progenitor cells prior to sex determination. A subset of these progenitors differentiates into Sertoli cells, a process characterized by a highly dynamic genetic program consisting of sequential waves of gene expression. Another subset of multipotent cells maintains their progenitor state but undergoes significant transcriptional changes restricting their competence toward a steroidogenic fate required for the differentiation of fetal Leydig cells. Our findings confirm the presence of a unique multipotent progenitor population in the gonadal primordium that gives rise to both supporting and interstitial lineages. These also provide the most granular analysis of the transcriptional events occurring during testicular cell-fate commitment.

Systematic proteome and proteostasis profiling in human Trisomy 21 fibroblast cells.

Down syndrome (DS) is mostly caused by a trisomy of the entire Chromosome 21 (Trisomy 21, T21). Here, we use SWATH mass spectrometry to quantify protein abundance and protein turnover in fibroblasts from a monozygotic twin pair discordant for T21, and to profile protein expression in 11 unrelated DS individuals and matched controls. The integration of the steady-state and turnover proteomic data indicates that protein-specific degradation of members of stoichiometric complexes is a major determinant of T21 gene dosage outcome, both within and between individuals. This effect is not apparent from genomic and transcriptomic data. The data also reveal that T21 results in extensive proteome remodeling, affecting proteins encoded by all chromosomes. Finally, we find broad, organelle-specific post-transcriptional effects such as significant downregulation of the mitochondrial proteome contributing to T21 hallmarks. Overall, we provide a valuable proteomic resource to understand the origin of DS phenotypic manifestations.

The genomic landscape of human cellular circadian variation points to a novel role for the signalosome.

The importance of natural gene expression variation for human behavior is undisputed, but its impact on circadian physiology remains mostly unexplored. Using umbilical cord fibroblasts, we have determined by genome-wide association how common genetic variation impacts upon cellular circadian function. Gene set enrichment points to differences in protein catabolism as one major source of clock variation in humans. The two most significant alleles regulated expression of COPS7B, a subunit of the COP9 signalosome. We further show that the signalosome complex is imported into the nucleus in timed fashion to stabilize the essential circadian protein BMAL1, a novel mechanism to oppose its proteasome-mediated degradation. Thus, circadian clock properties depend in part upon a genetically-encoded competition between stabilizing and destabilizing forces, and genetic alterations in these mechanisms provide one explanation for human chronotype.

SPINK2 deficiency causes infertility by inducing sperm defects in heterozygotes and azoospermia in homozygotes.

Azoospermia, characterized by the absence of spermatozoa in the ejaculate, is a common cause of male infertility with a poorly characterized etiology. Exome sequencing analysis of two azoospermic brothers allowed the identification of a homozygous splice mutation in SPINK2, encoding a serine protease inhibitor believed to target acrosin, the main sperm acrosomal protease. In accord with these findings, we observed that homozygous Spink2 KO male mice had azoospermia. Moreover, despite normal fertility, heterozygous male mice had a high rate of morphologically abnormal spermatozoa and a reduced sperm motility. Further analysis demonstrated that in the absence of Spink2, protease-induced stress initiates Golgi fragmentation and prevents acrosome biogenesis leading to spermatid differentiation arrest. We also observed a deleterious effect of acrosin overexpression in HEK cells, effect that was alleviated by SPINK2 coexpression confirming its role as acrosin inhibitor. These results demonstrate that SPINK2 is necessary to neutralize proteases during their cellular transit toward the acrosome and that its deficiency induces a pathological continuum ranging from oligoasthenoteratozoospermia in heterozygotes to azoospermia in homozygotes.

Detection of Imprinted Genes by Single-Cell Allele-Specific Gene Expression.

Genomic imprinting results in parental-specific gene expression. Imprinted genes are involved in the etiology of rare syndromes and have been associated with common diseases such as diabetes and cancer. Standard RNA bulk cell sequencing applied to whole-tissue samples has been used to detect imprinted genes in human and mouse models. However, lowly expressed genes cannot be detected by using RNA bulk approaches. Here, we report an original and robust method that combines single-cell RNA-seq and whole-genome sequencing into an optimized statistical framework to analyze genomic imprinting in specific cell types and in different individuals. Using samples from the probands of 2 family trios and 3 unrelated individuals, 1,084 individual primary fibroblasts were RNA sequenced and more than 700,000 informative heterozygous single-nucleotide variations (SNVs) were genotyped. The allele-specific coverage per gene of each SNV in each single cell was used to fit a beta-binomial distribution to model the likelihood of a gene being expressed from one and the same allele. Genes presenting a significant aggregate allelic ratio (between 0.9 and 1) were retained to identify of the allelic parent of origin. Our approach allowed us to validate the imprinting status of all of the known imprinted genes expressed in fibroblasts and the discovery of nine putative imprinted genes, thereby demonstrating the advantages of single-cell over bulk RNA-seq to identify imprinted genes. The proposed single-cell methodology is a powerful tool for establishing a cell type-specific map of genomic imprinting.

A Case of Wiedemann-Steiner Syndrome Associated with a 46,XY Disorder of Sexual Development and Gonadal Dysgenesis.

We report the case of a female patient suffering from a 46,XY disorder of sexual development (DSD) with complete gonadal dysgenesis and Wiedemann-Steiner Syndrome (WDSTS). The coexistence of these 2 conditions has not yet been reported. Using whole exome sequencing and comparative genome hybridization array, we identified a de novo MLL/KMT2A gene nonsense mutation which explains the WDSTS phenotype. In addition, we discovered novel genetic variants, which could explain the testicular dysgenesis observed in the patient, a maternally inherited 167-kb duplication of DAAM2 and MOCS1 genes and a de novo LRRC33/NRROS gene mutation. These genes, some of which are expressed during mouse gonadal development, could be considered as potentially new candidate genes for DSD.

DNA-Methylation Patterns in Trisomy 21 Using Cells from Monozygotic Twins.

DNA methylation is essential in mammalian development. We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in Down syndrome (DS). In order to determine the methylation differences in T21 without interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. Reduced Representation Bisulfite Sequencing (RRBS) revealed 35 differentially methylated promoter regions (DMRs) (Absolute methylation differences = 25%, FDR < 0.001) in MZ twins discordant for T21 that have also been observed in comparison between unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis (FDR = 1.60 e -03) and include genes of the HOXB and HOXD clusters. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) observed in the iPSC of the T21 versus normal twin. Altogether, the results of this study highlight the epigenetic effects of the extra chromosome 21 in T21 on loci outside of this chromosome that are relevant to DS associated phenotypes.

HSA21 Single-Minded 2 (Sim2) Binding Sites Co-Localize with Super-Enhancers and Pioneer Transcription Factors in Pluripotent Mouse ES Cells.

The HSA21 encoded Single-minded 2 (SIM2) transcription factor has key neurological functions and is a good candidate to be involved in the cognitive impairment of Down syndrome. We aimed to explore the functional capacity of SIM2 by mapping its DNA binding sites in mouse embryonic stem cells. ChIP-sequencing revealed 1229 high-confidence SIM2-binding sites. Analysis of the SIM2 target genes confirmed the importance of SIM2 in developmental and neuronal processes and indicated that SIM2 may be a master transcription regulator. Indeed, SIM2 DNA binding sites share sequence specificity and overlapping domains of occupancy with master transcription factors such as SOX2, OCT4 (Pou5f1), NANOG or KLF4. The association between SIM2 and these pioneer factors is supported by co-immunoprecipitation of SIM2 with SOX2, OCT4, NANOG or KLF4. Furthermore, the binding of SIM2 marks a particular sub-category of enhancers known as super-enhancers. These regions are characterized by typical DNA modifications and Mediator co-occupancy (MED1 and MED12). Altogether, we provide evidence that SIM2 binds a specific set of enhancer elements thus explaining how SIM2 can regulate its gene network in neuronal features.

Opposite phenotypes of muscle strength and locomotor function in mouse models of partial trisomy and monosomy 21 for the proximal Hspa13-App region.

The trisomy of human chromosome 21 (Hsa21), which causes Down syndrome (DS), is the most common viable human aneuploidy. In contrast to trisomy, the complete monosomy (M21) of Hsa21 is lethal, and only partial monosomy or mosaic monosomy of Hsa21 is seen. Both conditions lead to variable physiological abnormalities with constant intellectual disability, locomotor deficits, and altered muscle tone. To search for dosage-sensitive genes involved in DS and M21 phenotypes, we created two new mouse models: the Ts3Yah carrying a tandem duplication and the Ms3Yah carrying a deletion of the Hspa13-App interval syntenic with 21q11.2-q21.3. Here we report that the trisomy and the monosomy of this region alter locomotion, muscle strength, mass, and energetic balance. The expression profiling of skeletal muscles revealed global changes in the regulation of genes implicated in energetic metabolism, mitochondrial activity, and biogenesis. These genes are downregulated in Ts3Yah mice and upregulated in Ms3Yah mice. The shift in skeletal muscle metabolism correlates with a change in mitochondrial proliferation without an alteration in the respiratory function. However, the reactive oxygen species (ROS) production from mitochondrial complex I decreased in Ms3Yah mice, while the membrane permeability of Ts3Yah mitochondria slightly increased. Thus, we demonstrated how the Hspa13-App interval controls metabolic and mitochondrial phenotypes in muscles certainly as a consequence of change in dose of Gabpa, Nrip1, and Atp5j. Our results indicate that the copy number variation in the Hspa13-App region has a peripheral impact on locomotor activity by altering muscle function.

Research resource: the dynamic transcriptional profile of sertoli cells during the progression of spermatogenesis.

Sertoli cells (SCs), the only somatic cells within seminiferous tubules, associate intimately with developing germ cells. They not only provide physical and nutritional support but also secrete factors essential to the complex developmental processes of germ cell proliferation and differentiation. The SC transcriptome must therefore adapt rapidly during the different stages of spermatogenesis. We report comprehensive genome-wide expression profiles of pure populations of SCs isolated at 5 distinct stages of the first wave of mouse spermatogenesis, using RNA sequencing technology. We were able to reconstruct about 13 901 high-confidence, nonredundant coding and noncoding transcripts, characterized by complex alternative splicing patterns with more than 45% comprising novel isoforms of known genes. Interestingly, roughly one-fifth (2939) of these genes exhibited a dynamic expression profile reflecting the evolving role of SCs during the progression of spermatogenesis, with stage-specific expression of genes involved in biological processes such as cell cycle regulation, metabolism and energy production, retinoic acid synthesis, and blood-testis barrier biogenesis. Finally, regulatory network analysis identified the transcription factors endothelial PAS domain-containing protein 1 (EPAS1/Hif2α), aryl hydrocarbon receptor nuclear translocator (ARNT/Hif1ß), and signal transducer and activator of transcription 1 (STAT1) as potential master regulators driving the SC transcriptional program. Our results highlight the plastic transcriptional landscape of SCs during the progression of spermatogenesis and provide valuable resources to better understand SC function and spermatogenesis and its related disorders, such as male infertility.