The transcriptional and regulatory identity of erythropoietin producing cells.

Erythropoietin (Epo) is the master regulator of erythropoiesis and oxygen homeostasis. Despite its physiological importance, the molecular and genomic contexts of the cells responsible for renal Epo production remain unclear, limiting more-effective therapies for anemia. Here, we performed single-cell RNA and transposase-accessible chromatin (ATAC) sequencing of an Epo reporter mouse to molecularly identify Epo-producing cells under hypoxic conditions. Our data indicate that a distinct population of kidney stroma, which we term Norn cells, is the major source of endocrine Epo production in mice. We use these datasets to identify the markers, signaling pathways and transcriptional circuits characteristic of Norn cells. Using single-cell RNA sequencing and RNA in situ hybridization in human kidney tissues, we further provide evidence that this cell population is conserved in humans. These preliminary findings open new avenues to functionally dissect EPO gene regulation in health and disease and may serve as groundwork to improve erythropoiesis-stimulating therapies.

LGR5 expressing skin fibroblasts define a major cellular hub perturbed in scleroderma.

Systemic sclerosis (scleroderma, SSc) is an incurable autoimmune disease with high morbidity and mortality rates. Here, we conducted a population-scale single-cell genomic analysis of skin and blood samples of 56 healthy controls and 97 SSc patients at different stages of the disease. We found immune compartment dysfunction only in a specific subtype of diffuse SSc patients but global dysregulation of the stromal compartment, particularly in a previously undefined subset of LGR5+-scleroderma-associated fibroblasts (ScAFs). ScAFs are perturbed morphologically and molecularly in SSc patients. Single-cell multiome profiling of stromal cells revealed ScAF-specific markers, pathways, regulatory elements, and transcription factors underlining disease development. Systematic analysis of these molecular features with clinical metadata associates specific ScAF targets with disease pathogenesis and SSc clinical traits. Our high-resolution atlas of the sclerodermatous skin spectrum will enable a paradigm shift in the understanding of SSc disease and facilitate the development of biomarkers and therapeutic strategies.

Position effects at the FGF8 locus are associated with femoral hypoplasia.

Copy-number variations (CNVs) are a common cause of congenital limb malformations and are interpreted primarily on the basis of their effect on gene dosage. However, recent studies show that CNVs also influence the 3D genome chromatin organization. The functional interpretation of whether a phenotype is the result of gene dosage or a regulatory position effect remains challenging. Here, we report on two unrelated families with individuals affected by bilateral hypoplasia of the femoral bones, both harboring de novo duplications on chromosome 10q24.32. The ∼0.5 Mb duplications include FGF8, a key regulator of limb development and several limb enhancer elements. To functionally characterize these variants, we analyzed the local chromatin architecture in the affected individuals' cells and re-engineered the duplications in mice by using CRISPR-Cas9 genome editing. We found that the duplications were associated with ectopic chromatin contacts and increased FGF8 expression. Transgenic mice carrying the heterozygous tandem duplication including Fgf8 exhibited proximal shortening of the limbs, resembling the human phenotype. To evaluate whether the phenotype was a result of gene dosage, we generated another transgenic mice line, carrying the duplication on one allele and a concurrent Fgf8 deletion on the other allele, as a control. Surprisingly, the same malformations were observed. Capture Hi-C experiments revealed ectopic interaction with the duplicated region and Fgf8, indicating a position effect. In summary, we show that duplications at the FGF8 locus are associated with femoral hypoplasia and that the phenotype is most likely the result of position effects altering FGF8 expression rather than gene dosage effects.

Variants in KIAA0825 underlie autosomal recessive postaxial polydactyly.

Postaxial polydactyly (PAP) is a common limb malformation that often leads to cosmetic and functional complications. Molecular evaluation of polydactyly can serve as a tool to elucidate genetic and signaling pathways that regulate limb development, specifically, the anterior-posterior specification of the limb. To date, only five genes have been identified for nonsyndromic PAP: FAM92A, GLI1, GLI3, IQCE and ZNF141. In this study, two Pakistani multiplex consanguineous families with autosomal recessive nonsyndromic PAP were clinically and molecularly evaluated. From both pedigrees, a DNA sample from an affected member underwent exome sequencing. In each family, we identified a segregating frameshift (c.591dupA [p.(Q198Tfs*21)]) and nonsense variant (c.2173A > T [p.(K725*)]) in KIAA0825 (also known as C5orf36). Although KIAA0825 encodes a protein of unknown function, it has been demonstrated that its murine ortholog is expressed during limb development. Our data contribute to the establishment of a catalog of genes important in limb patterning, which can aid in diagnosis and obtaining a better understanding of the biology of polydactyly.

H2AFY promoter deletion causes PITX1 endoactivation and Liebenberg syndrome.

BACKGROUND: Structural variants (SVs) affecting non-coding cis-regulatory elements are a common cause of congenital limb malformation. Yet, the functional interpretation of these non-coding variants remains challenging. The human Liebenberg syndrome is characterised by a partial transformation of the arms into legs and has been shown to be caused by SVs at the PITX1 locus leading to its misregulation in the forelimb by its native enhancer element Pen. This study aims to elucidate the genetic cause of an unsolved family with a mild form of Liebenberg syndrome and investigate the role of promoters in long-range gene regulation. METHODS: Here, we identify SVs by whole genome sequencing (WGS) and use CRISPR-Cas9 genome editing in transgenic mice to assign pathogenicity to the SVs. RESULTS: In this study, we used WGS in a family with three mildly affected individuals with Liebenberg syndrome and identified the smallest deletion described so far including the first non-coding exon of H2AFY. To functionally characterise the variant, we re-engineered the 8.5 kb deletion using CRISPR-Cas9 technology in the mouse and showed that the promoter of the housekeeping gene H2afy insulates the Pen enhancer from Pitx1 in forelimbs; its loss leads to misexpression of Pitx1 by the pan-limb activity of the Pen enhancer causing Liebenberg syndrome. CONCLUSION: Our data indicate that housekeeping promoters may titrate promiscuous enhancer activity to ensure normal morphogenesis. The deletion of the H2AFY promoter as a cause of Liebenberg syndrome highlights this new mutational mechanism and its role in congenital disease.

Dynamic 3D chromatin architecture contributes to enhancer specificity and limb morphogenesis.

The regulatory specificity of enhancers and their interaction with gene promoters is thought to be controlled by their sequence and the binding of transcription factors. By studying Pitx1, a regulator of hindlimb development, we show that dynamic changes in chromatin conformation can restrict the activity of enhancers. Inconsistent with its hindlimb-restricted expression, Pitx1 is controlled by an enhancer (Pen) that shows activity in forelimbs and hindlimbs. By Capture Hi-C and three-dimensional modeling of the locus, we demonstrate that forelimbs and hindlimbs have fundamentally different chromatin configurations, whereby Pen and Pitx1 interact in hindlimbs and are physically separated in forelimbs. Structural variants can convert the inactive into the active conformation, thereby inducing Pitx1 misexpression in forelimbs, causing partial arm-to-leg transformation in mice and humans. Thus, tissue-specific three-dimensional chromatin conformation can contribute to enhancer activity and specificity in vivo and its disturbance can result in gene misexpression and disease.

Noncoding copy-number variations are associated with congenital limb malformation.

PurposeCopy-number variants (CNVs) are generally interpreted by linking the effects of gene dosage with phenotypes. The clinical interpretation of noncoding CNVs remains challenging. We investigated the percentage of disease-associated CNVs in patients with congenital limb malformations that affect noncoding cis-regulatory sequences versus genes sensitive to gene dosage effects.MethodsWe applied high-resolution copy-number analysis to 340 unrelated individuals with isolated limb malformation. To investigate novel candidate CNVs, we re-engineered human CNVs in mice using clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing.ResultsOf the individuals studied, 10% harbored CNVs segregating with the phenotype in the affected families. We identified 31 CNVs previously associated with congenital limb malformations and four novel candidate CNVs. Most of the disease-associated CNVs (57%) affected the noncoding cis-regulatory genome, while only 43% included a known disease gene and were likely to result from gene dosage effects. In transgenic mice harboring four novel candidate CNVs, we observed altered gene expression in all cases, indicating that the CNVs had a regulatory effect either by changing the enhancer dosage or altering the topological associating domain architecture of the genome.ConclusionOur findings suggest that CNVs affecting noncoding regulatory elements are a major cause of congenital limb malformations.

Mutations in MYO1H cause a recessive form of central hypoventilation with autonomic dysfunction.

BACKGROUND: Congenital central hypoventilation syndrome (CCHS) is a rare life-threatening disorder of respiratory and autonomic regulation. It is classically caused by dominant mutations in the transcription factor PHOX2B. The objective of the present study was to identify the molecular cause of a recessive form of central hypoventilation with autonomic dysfunction. METHODS: Here, we used homozygosity mapping and whole-genome sequencing in a consanguineous family with CCHS in combination with functional analyses in CRISPR/Cas9 engineered mice. RESULTS: We report on a consanguineous family with three affected children, all tested PHOX2B mutation negative, presenting with alveolar hypoventilation and symptoms of autonomic dysregulation. Whole-genome sequencing revealed a homozygous frameshift mutation in exon 25 of the MYO1H gene (c.2524_2524delA) segregating with the phenotype in the family. MYO1H encodes for the unconventional myosin IH, which is thought to function as a motor protein in intracellular transport and vesicle trafficking. We show that Myo1h is broadly expressed in the mouse lower medulla, including the CO2-sensitive Phox2b+ retrotrapezoid neurons. To test the pathogenicity of the variant, we engineered two Myo1h mutant mouse strains: the first strain (Myo1h*) resembling the human mutation and the second being a full knock-out (Myo1hFS ). Whole-body plethysmography studies in Myo1h* newborns with the re-engineered human mutation revealed hypoventilation and a blunted response to CO2, recapitulating the breathing phenotype observed in the kindred. CONCLUSIONS: Our results identify MYO1H as an important gene in CO2 sensitivity and respiratory control and as the cause of a rare recessive form of congenital central hypoventilation.

Deletions, Inversions, Duplications: Engineering of Structural Variants using CRISPR/Cas in Mice.

Structural variations (SVs) contribute to the variability of our genome and are often associated with disease. Their study in model systems was hampered until now by labor-intensive genetic targeting procedures and multiple mouse crossing steps. Here we present the use of CRISPR/Cas for the fast (10 weeks) and efficient generation of SVs in mice. We specifically produced deletions, inversions, and also duplications at six different genomic loci ranging from 1.1 kb to 1.6 Mb with efficiencies up to 42%. After PCR-based selection, clones were successfully used to create mice via aggregation. To test the practicability of the method, we reproduced a human 500 kb disease-associated deletion and were able to recapitulate the human phenotype in mice. Furthermore, we evaluated the regulatory potential of a large genomic interval by deleting a 1.5 Mb fragment. The method presented permits rapid in vivo modeling of genomic rearrangements.