Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A genome-wide association study identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we use mouse enhancer assays, three mouse enhancer knockouts, and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterize a sequence, PEC7, which overlaps the AIS-associated variant, and find it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, or deletion of both sequences lead to a kinky tail phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicates several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

TAD evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find 14% of all human TAD boundaries to be shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons compared to species-specific boundaries. CRISPR-Cas9 knockouts of an ultraconserved boundary in a mouse model lead to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in the upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations and showcases the functional importance of TAD evolution.

Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A GWAS identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we used mouse enhancer assays, three mouse enhancer knockouts and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterized a sequence, PEC7, that overlaps the AIS-associated variant, and found it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, and deletion of both sequences led to a kinky phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicated several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

TAD Evolutionary and functional characterization reveals diversity in mammalian TAD boundary properties and function

Topological associating domains (TADs) are self-interacting genomic units crucial for shaping gene regulation patterns. Despite their importance, the extent of their evolutionary conservation and its functional implications remain largely unknown. In this study, we generate Hi-C and ChIP-seq data and compare TAD organization across four primate and four rodent species, and characterize the genetic and epigenetic properties of TAD boundaries in correspondence to their evolutionary conservation. We find that only 14% of all human TAD boundaries are shared among all eight species (ultraconserved), while 15% are human-specific. Ultraconserved TAD boundaries have stronger insulation strength, CTCF binding, and enrichment of older retrotransposons, compared to species-specific boundaries. CRISPR-Cas9 knockouts of two ultraconserved boundaries in mouse models leads to tissue-specific gene expression changes and morphological phenotypes. Deletion of a human-specific boundary near the autism-related AUTS2 gene results in upregulation of this gene in neurons. Overall, our study provides pertinent TAD boundary evolutionary conservation annotations, and showcase the functional importance of TAD evolution.

A role of the frontotemporal lobar degeneration risk factor TMEM106B in myelination

TMEM106B encodes a lysosomal membrane protein and was initially identified as a risk factor for frontotemporal lobar degeneration. Recently, a dominant D252N mutation in TMEM106B was shown to cause hypomyelinating leukodystrophy. However, how TMEM106B regulates myelination is still unclear. Here we show that TMEM106B is expressed and localized to the lysosome compartment in oligodendrocytes. TMEM106B deficiency in mice results in myelination defects with a significant reduction of protein levels of proteolipid protein (PLP) and myelin oligodendrocyte glycoprotein (MOG), the membrane proteins found in the myelin sheath. The levels of many lysosome proteins are significantly decreased in the TMEM106B-deficient Oli-neu oligodendroglial precursor cell line. TMEM106B physically interacts with the lysosomal protease cathepsin D and is required to maintain proper cathepsin D levels in oligodendrocytes. Furthermore, we found that TMEM106B deficiency results in lysosome clustering in the perinuclear region and a decrease in lysosome exocytosis and cell surface PLP levels. Moreover, we found that the D252N mutation abolished lysosome enlargement and lysosome acidification induced by wild-type TMEM106B overexpression. Instead, it stimulates lysosome clustering near the nucleus as seen in TMEM106B-deficient cells. Our results support that TMEM106B regulates myelination through modulation of lysosome function in oligodendrocytes.

Loss of TMEM106B and PGRN leads to severe lysosomal abnormalities and neurodegeneration in mice

Haploinsufficiency of progranulin (PGRN) is a leading cause of frontotemporal lobar degeneration (FTLD). Loss of PGRN leads to lysosome dysfunction during aging. TMEM106B, a gene encoding a lysosomal membrane protein, is the main risk factor for FTLD with PGRN haploinsufficiency. But how TMEM106B affects FTLD disease progression remains to be determined. Here, we report that TMEM106B deficiency in mice leads to accumulation of lysosome vacuoles at the distal end of the axon initial segment in motor neurons and the development of FTLD‐related pathology during aging. Ablation of both PGRN and TMEM106B in mice results in severe neuronal loss and glial activation in the spinal cord, retina, and brain. Enlarged lysosomes are frequently found in both microglia and astrocytes. Loss of both PGRN and TMEM106B results in an increased accumulation of lysosomal vacuoles in the axon initial segment of motor neurons and enhances the manifestation of FTLD phenotypes with a much earlier onset. These results provide novel insights into the role of TMEM106B in the lysosome, in brain aging, and in FTLD pathogenesis.