Detection of genes regulated by Lmx1b during limb dorsalization

Versican controlled by Lmx1b regulates hyaluronate density and hydration for semicircular canal morphogenesis

During inner ear semicircular canal morphogenesis in zebrafish, patterned canal-genesis zones express genes for extracellular matrix component synthesis. These include hyaluronan and the hyaluronan-binding chondroitin sulfate proteoglycan Versican, which are abundant in the matrices of many developing organs. Charged hyaluronate polymers play a key role in canal morphogenesis through osmotic swelling. However, the developmental factor(s) that pattern the synthesis of the matrix components and regulation of hyaluronate density and swelling are unknown. Here, we identify the transcription factor Lmx1b as a positive transcriptional regulator of hyaluronan, Versican, and chondroitin synthesis genes crucial for canal morphogenesis. We show that Versican regulates hyaluronan density through its protein core, whereas the charged chondroitin side chains contribute to the hydration of hyaluronate-containing extracellular matrices. Versican-tuned properties of hyaluronate matrices may be a broadly used mechanism in morphogenesis with important implications for understanding diseases in which these matrices are impaired, and for hydrogel engineering for tissue regeneration.

Versican controlled by Lmx1b regulates hyaluronate density and hydration for semicircular canal morphogenesis

During inner ear semicircular canal morphogenesis in zebrafish, patterned canal-genesis zones express genes for extracellular matrix component synthesis. These include hyaluronan and the hyaluronan-binding chondroitin sulfate proteoglycan Versican, which are abundant in the matrices of many developing organs. Charged hyaluronate polymers play a key role in canal morphogenesis through osmotic swelling. However, the developmental factor(s) that control the synthesis of the matrix components and regulation of hyaluronate density and swelling are unknown. Here, we identify the transcription factor, Lmx1b, as a positive transcriptional regulator of hyaluronan, Versican, and chondroitin synthesis genes crucial for canal morphogenesis. We show that Versican regulates hyaluronan density through its protein core, whereas the charged chondroitin side chains contribute to the osmotic swelling of hyaluronate. Versican-tuned properties of hyaluronate matrices may be a broadly used mechanism in morphogenesis with important implications for understanding diseases where these matrices are impaired, and for hydrogel engineering for tissue regeneration.

Sox, Fox, and Lmx1b binding sites differentially regulate a Gdf5-Associated regulatory region during elbow development

Introduction: The articulating ends of limb bones have precise morphology and asymmetry that ensures proper joint function. Growth differentiation factor 5 (Gdf5) is a secreted morphogen involved in cartilage and bone development that contributes to the architecture of developing joints. Dysregulation of Gdf5 results in joint dysmorphogenesis often leading to progressive joint degeneration or osteoarthritis (OA). The transcription factors and cis-regulatory modules (CRMs) that regulate Gdf5 expression are not well characterized. We previously identified a Gdf5-associated regulatory region (GARR) that contains predicted binding sites for Lmx1b, Osr2, Fox, and the Sox transcription factors. These transcription factors are recognized factors involved in joint morphogenesis and skeletal development. Methods: We used in situ hybridization to Gdf5, Col2A1, and the transcription factors of interest in developing chicken limbs to determine potential overlap in expression. We further analyzed scRNA-seq data derived from limbs and knees in published mouse and chicken datasets, identifying cells with coexpression of Gdf5 and the transcription factors of interest. We also performed site-directed mutatgenesis of the predicted transcription factor binding sites in a GARR-reporter construct and determined any change in activity using targeted regional electroporation (TREP) in micromass and embryonic chicken wing bioassays. Results: Gdf5 expression overlapped the expression of these transcription factors during joint development both by in situ hybridization (ISH) and scRNA-seq analyses. Within the GARR CRM, mutation of two binding sites common to Fox and Sox transcripstion factors reduced enhancer activity to background levels in micromass cultures and in ovo embryonic chicken wing bioassays, whereas mutation of two Sox-only binding sites caused a significant increase in activity. These results indicate that the Fox/Sox binding sites are required for activity, while the Sox-only sites are involved in repression of activity. Mutation of Lmx1b binding sites in GARR caused an overall reduction in enhancer activity in vitro and a dorsal reduction in ovo. Despite a recognized role for Osr2 in joint development, disruption of the predicted Osr2 site did not alter GARR activity. Conclusion: Taken together, our data indicates that GARR integrates positive, repressive, and asymmetrical inputs to fine-tune the expression of Gdf5 during elbow joint development.

A Genome-Wide Association Study Meta-Analysis of Alpha Angle Suggests Cam-Type Morphology May Be a Specific Feature of Hip Osteoarthritis in Older Adults

OBJECTIVE

To examine the genetic architecture of cam morphology using alpha angle (AA) as a proxy measure and conduct an AA genome-wide association study (GWAS) followed by Mendelian randomization (MR) to evaluate its causal relationship with hip osteoarthritis (OA).

METHODS

Observational analyses examined associations between AA measurements derived from hip dual x-ray absorptiometry (DXA) scans from the UK Biobank study and radiographic hip OA outcomes and subsequent total hip replacement. Following these analyses, an AA GWAS meta-analysis was performed (N = 44,214) using AA measurements previously derived in the Rotterdam Study. Linkage disequilibrium score regression assessed the genetic correlation between AA and hip OA. Genetic associations considered significant (P < 5 × 10-8 ) were used as AA genetic instrument for 2-sample MR analysis.

RESULTS

DXA-derived AA showed expected associations between AA and radiographic hip OA (adjusted odds ratio [OR] 1.63 [95% confidence interval (95% CI) 1.58, 1.67]) and between AA and total hip replacement (adjusted hazard ratio 1.45 [95% CI 1.33, 1.59]) in the UK Biobank study cohort. The heritability of AA was 10%, and AA had a moderate genetic correlation with hip OA (rg  = 0.26 [95% CI 0.10, 0.43]). Eight independent genetic signals were associated with AA. Two-sample MR provided weak evidence of causal effects of AA on hip OA risk (inverse variance weighted OR 1.84 [95% CI 1.14, 2.96], P = 0.01). In contrast, genetic predisposition for hip OA had stronger evidence of a causal effect on increased AA (inverse variance weighted β = 0.09 [95% CI 0.04, 0.13], P = 4.58 × 10-5 ).

CONCLUSION

Expected observational associations between AA and related clinical outcomes provided face validity for the DXA-derived AA measurements. Evidence of bidirectional associations between AA and hip OA, particularly for risk of hip OA on AA, suggests that hip shape modeling secondary to a genetic predisposition to hip OA contributes to the well-established relationship between hip OA and cam morphology in older adults.

Identification of limb-specific Lmx1b auto-regulatory modules with Nail-patella syndrome pathogenicity

LMX1B haploinsufficiency causes Nail-patella syndrome (NPS; MIM 161200), characterized by nail dysplasia, absent/hypoplastic patellae, chronic kidney disease, and glaucoma. Accordingly in mice, Lmx1b has been shown to play crucial roles in the development of the limb, kidney and eye. Although one functional allele of Lmx1b appears adequate for development, Lmx1b null mice display ventral-ventral distal limbs with abnormal kidney, eye and cerebellar development, more disruptive, but fully concordant with NPS. In Lmx1b functional knockouts (KOs), Lmx1b transcription in the limb is decreased nearly 6-fold, indicating autoregulation. Herein, we report on two conserved Lmx1b-associated cis-regulatory modules (LARM1 and LARM2) that are bound by Lmx1b, amplify Lmx1b expression with unique spatial modularity in the limb, and are necessary for Lmx1b-mediated limb dorsalization. These enhancers, being conserved across vertebrates (including coelacanth, but not other fish species), and required for normal locomotion, provide a unique opportunity to study the role of dorsalization in the fin to limb transition. We also report on two NPS patient families with normal LMX1B coding sequence, but with loss-of-function variations in the LARM1/2 region, stressing the role of regulatory modules in disease pathogenesis.

Nail-patella syndrome (NPS) is characterized by nail dysplasia, absent/hypoplastic patellae, chronic kidney disease, and glaucoma and can be caused by haploinsufficiency of LMX1B; however, not all patients harbor pathogenic LMX1B mutations. Here the authors show that loss-of-function variations in upstream enhancer sequences are responsible for a limb specific form of human NPS.

Hereditary Hearing Impairment with Cutaneous Abnormalities

Syndromic hereditary hearing impairment (HHI) is a clinically and etiologically diverse condition that has a profound influence on affected individuals and their families. As cutaneous findings are more apparent than hearing-related symptoms to clinicians and, more importantly, to caregivers of affected infants and young individuals, establishing a correlation map of skin manifestations and their underlying genetic causes is key to early identification and diagnosis of syndromic HHI. In this article, we performed a comprehensive PubMed database search on syndromic HHI with cutaneous abnormalities, and reviewed a total of 260 relevant publications. Our in-depth analyses revealed that the cutaneous manifestations associated with HHI could be classified into three categories: pigment, hyperkeratosis/nail, and connective tissue disorders, with each category involving distinct molecular pathogenesis mechanisms. This outline could help clinicians and researchers build a clear atlas regarding the phenotypic features and pathogenetic mechanisms of syndromic HHI with cutaneous abnormalities, and facilitate clinical and molecular diagnoses of these conditions.

LMX1B mRNA expression and its gene body CpG methylation are valuable prognostic biomarkers for laryngeal squamous cell carcinoma

This study aimed to explore the prognostic value of LMX1B mRNA expression and the methylation of its CpG sites in patients with laryngeal squamous cell carcinoma (LSCC). An in-silicon analysis was performed using data from the cancer genome atlas (TCGA)-Head and Neck Squamous Carcinoma (HNSC). After screening, 112 LSCC and 10 adjacent normal tissues were identified as eligible samples for analysis. Results showed that LMX1B expression was significantly upregulated in the cancer tissues (p < 0.01) and was an independent prognostic indicator in terms of OS (HR: 1.233, 95%CI: 1.082-1.405, p = 0.002) and RFS (HR: 1.200, 95%CI: 1.002-1.438, p = 0.048). By examining the methylation profile of 55 CpG sites in LMX1B locus, we found that the promoter methylation status was irrelevant to LMX1B expression. In comparison, LMX1B expression was generally positively correlated with gene body methylation. Among the gene body CpG sites, cg13600622 methylation showed a better predictive value than LMX1B expression in terms of OS (HR: 12.363, 95%CI: 1.076-142.033, p = 0.043), while cg14204784 methylation was a better marker of shorter RFS (HR: 12.363, 95%CI: 1.076-142.033, p = 0.043). Among the known downstream genes of LMX1B, only NR4A2 expression showed a moderately negative correlation (Pearson's r = -0.54) with it in LSCC tissues. However, this correlation was inconsistent with previous publications those reported a positive correlation between them. Based on these findings, we infer that upregulated LMX1B mRNA expression had an independent prognostic value in LSCC patients. Increased gene body methylation might be an important mechanism of its upregulation. Among the gene body CpG sites, cg13600622 and cg14204784 methylation level might be better prognostic markers than LMX1B mRNA expression in terms of OS and RFS respectively.

Common and distinct transcriptional signatures of mammalian embryonic lethality

The Deciphering the Mechanisms of Developmental Disorders programme has analysed the morphological and molecular phenotypes of embryonic and perinatal lethal mouse mutant lines in order to investigate the causes of embryonic lethality. Here we show that individual whole-embryo RNA-seq of 73 mouse mutant lines (>1000 transcriptomes) identifies transcriptional events underlying embryonic lethality and associates previously uncharacterised genes with specific pathways and tissues. For example, our data suggest that Hmgxb3 is involved in DNA-damage repair and cell-cycle regulation. Further, we separate embryonic delay signatures from mutant line-specific transcriptional changes by developing a baseline mRNA expression catalogue of wild-type mice during early embryogenesis (4-36 somites). Analysis of transcription outside coding sequence identifies deregulation of repetitive elements in Morc2a mutants and a gene involved in gene-specific splicing. Collectively, this work provides a large scale resource to further our understanding of early embryonic developmental disorders.

Pigeon foot feathering reveals conserved limb identity networks

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

LHX2 Mediates the FGF-to-SHH Regulatory Loop during Limb Development

During limb development, fibroblast growth factors (Fgfs) govern proximal⁻distal outgrowth and patterning. FGFs also synchronize developmental patterning between the proximal⁻distal and anterior⁻posterior axes by maintaining Sonic hedgehog (Shh) expression in cells of the zone of polarizing activity (ZPA) in the distal posterior mesoderm. Shh, in turn, maintains Fgfs in the apical ectodermal ridge (AER) that caps the distal tip of the limb bud. Crosstalk between Fgf and Shh signaling is critical for patterned limb development, but the mechanisms underlying this feedback loop are not well-characterized. Implantation of Fgf beads in the proximal posterior limb bud can maintain SHH expression in the former ZPA domain (evident 3 h after application), while prolonged exposure (24 h) can induce SHH outside of this domain. Although temporally and spatially disparate, comparative analysis of transcriptome data from these different populations accentuated genes involved in SHH regulation. Comparative analysis identified 25 candidates common to both treatments, with eight linked to SHH expression or function. Furthermore, we demonstrated that LHX2, a LIM Homeodomain transcription factor, is an intermediate in the FGF-mediated regulation of SHH. Our data suggest that LHX2 acts as a competency factor maintaining distal posterior SHH expression subjacent to the AER.

Genome-wide association study in Japanese females identifies fifteen novel skin-related trait associations

Skin trait variation impacts quality-of-life, especially for females from the viewpoint of beauty. To investigate genetic variation related to these traits, we conducted a GWAS of various skin phenotypes in 11,311 Japanese women and identified associations for age-spots, freckles, double eyelids, straight/curly hair, eyebrow thickness, hairiness, and sweating. In silico annotation with RoadMap Epigenomics epigenetic state maps and colocalization analysis of GWAS and GTEx Project eQTL signals provided information about tissue specificity, candidate causal variants, and functional target genes. Novel signals for skin-spot traits neighboured AKAP1/MSI2 (rs17833789; P = 2.2 × 10-9), BNC2 (rs10810635; P = 2.1 × 10-22), HSPA12A (rs12259842; P = 7.1 × 10-11), PPARGC1B (rs251468; P = 1.3 × 10-21), and RAB11FIP2 (rs10444039; P = 5.6 × 10-21). HSPA12A SNPs were the only protein-coding gene eQTLs identified across skin-spot loci. Double edged eyelid analysis identified that a signal around EMX2 (rs12570134; P = 8.2 × 10-15) was also associated with expression of EMX2 and the antisense-RNA gene EMX2OS in brain putamen basal ganglia tissue. A known hair morphology signal in EDAR was associated with both eyebrow thickness (rs3827760; P = 1.7 × 10-9) and straight/curly hair (rs260643; P = 1.6 × 10-103). Excessive hairiness signals' top SNPs were also eQTLs for TBX15 (rs984225; P = 1.6 × 10-8), BCL2 (rs7226979; P = 7.3 × 10-11), and GCC2 and LIMS1 (rs6542772; P = 2.2 × 10-9). For excessive sweating, top variants in two signals in chr2:28.82-29.05 Mb (rs56089836; P = 1.7 × 10-11) were eQTLs for either PPP1CB or PLB1, while a top chr16:48.26-48.45 Mb locus SNP was a known ABCC11 missense variant (rs6500380; P = 6.8 × 10-10). In total, we identified twelve loci containing sixteen association signals, of which fifteen were novel. These findings will help dermatologic researchers better understand the genetic underpinnings of skin-related phenotypic variation in human populations.

Steroid-resistant nephrotic syndrome as the initial presentation of nail-patella syndrome: a case of a de novo LMX1B mutation

BACKGROUND

Nail-patella syndrome (NPS) is an autosomal dominant disorder caused by mutations in the LMX1B gene and is characterized by nail dysplasia, skeletal abnormalities, and nephropathy. We herein report a case of steroid-resistant nephrotic syndrome (SRNS) prior to overt orthopedic symptoms in a patient with NPS.

CASE PRESENTATION

A 24-year-old woman presented to our hospital with knee pain. She had poorly developed nails, hypoplastic patellas, dislocation of the elbows, and iliac horns in the pelvis. At the age of 7, she developed nephrotic syndrome and was diagnosed with primary focal segmental glomerulosclerosis by renal biopsy. She received long-term corticosteroid therapy with no obvious response. Her clinical course and orthopedic manifestations indicated NPS, and a genetic analysis showed a de novo mutation in the LMX1B gene (c.819 + 1G > A). Nephropathy in this case was considered to be associated with NPS. Therefore, we discontinued corticosteroids without the exacerbation of nephrotic syndrome.

CONCLUSIONS

Patients with NPS may develop nephrotic syndrome prior to overt orthopedic symptoms and only show non-specific findings in renal biopsy at an early stage of NPS nephropathy. Hereditary nephrotic syndrome, often presenting as childhood-onset SRNS, may also be difficult to diagnose in patients with the following conditions: renal symptoms prior to overt extrarenal symptoms, de novo mutations, and non-specific findings in renal biopsy. Therefore, in the management of SRNS in children, we need to reconsider the possibility of hereditary diseases such as NPS even without a family history.

Lmx1b-targeted cis-regulatory modules involved in limb dorsalization

Lmx1b is a homeodomain transcription factor responsible for limb dorsalization. Despite striking double-ventral (loss-of-function) and double-dorsal (gain-of-function) limb phenotypes, no direct gene targets in the limb have been confirmed. To determine direct targets, we performed a chromatin immunoprecipitation against Lmx1b at E12.5 followed by next generation sequencing (ChIP-seq). Nearly 84% (n=617) of the Lmx1b-bound genomic intervals (LBIs) identified overlap with chromatin regulatory marks indicative of potential cis-regulatory modules (PCRMs). In addition, 73 LBIs mapped to known CRMs active during limb development. We compared Lmx1b-bound PCRMs to genes differentially expressed by Lmx1b and found 292 PCRMs within 1 Mb of 254 Lmx1b-regulated genes. Gene ontologic analysis suggests that Lmx1b targets extracellular matrix production, bone/joint formation, axonal guidance, vascular development, cell proliferation and cell movement. We validated the functional activity of a PCRM associated with joint-related Gdf5 that provides a mechanism for Lmx1b-mediated joint modification and a PCRM associated with Lmx1b that suggests a role in autoregulation. This is the first report to describe genome-wide Lmx1b binding during limb development, directly linking Lmx1b to targets that accomplish limb dorsalization.

Genetics of syndromic and non-syndromic hereditary nail disorders

The nail is a unique epithelial skin appendage made up of a fully keratinized nail plate. The nail can be affected in several systemic illnesses, dermatological diseases, and inherited nail disorders. Nail dystrophies can present as isolated disorders or as a part of syndromes. Substantial progress has been achieved in the management and diagnosis of nail diseases; however, not much is known about the underlying molecular controls of nail growth. The homeostasis and development of the nail appendage depend on the intricate interactions between the epidermis and underlying mesenchyme, and comprise different signaling pathways such as the WNT signaling pathway. Digit-tip regeneration in mice and humans has been a known fact for the past six decades; however, only recently the underlying biological mechanisms by which the nail organ achieves digit regeneration have been elucidated. Moreover, significant progress has been made in identifying nail stem cells and localizing stem cell niches in the nail unit. More fascinating, however, is the role they play in orchestrating the processes that lead to the regeneration of the digit. Further elucidating the role of nail stem cells and the signaling pathways driving epithelial-mesenchymal interactions in the nail unit might contribute to the development of novel therapeutic tools for amputees.

Muscle-type Identity of Proprioceptors Specified by Spatially Restricted Signals from Limb Mesenchyme

The selectivity with which proprioceptive sensory neurons innervate their central and peripheral targets implies that they exhibit distinctions in muscle-type identity. The molecular correlates of proprioceptor identity and its origins remain largely unknown, however. In screens to define muscle-type proprioceptor character, we find all-or-none differences in gene expression for proprioceptors that control antagonistic muscles at a single hindlimb joint. Analysis of three of these genes, cadherin13 (cdh13), semaphorin5a (sema5a), and cartilage-acidic protein-1 (crtac1), reveals expression in proprioceptor subsets that supply muscle groups located at restricted dorsoventral and proximodistal domains of the limb. Genetically altering the dorsoventral character of the limb mesenchyme elicits a change in the profile of proprioceptor cdh13, sema5a, and crtac1 expression. These findings indicate that proprioceptors acquire aspects of their muscle-type identity in response to mesenchymal signals expressed in restricted proximodistal and dorsoventral domains of the developing limb.

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis.

The signaling pathways of LMX1B and its role in glomerulosclerosis

LMX1B, a developmental LIM-homeodomain transcription factor, is widely expressed in vertebrate embryos, and it takes part in the development of diverse structures such as limbs, kidneys, eyes, brains, etc. LMX1B contributes to transcriptional regulation of glomerular basement membrane (GBM) collagen expression by podocytes. The normal function of podocytes and the normal morphology of GBM are very important to maintain the healthy renal filtration barrier. Recent discoveries find that the LMX1B gene is pivotal in glomus development and it is implicated in the dysfunction of the podocytes. Here, we review the signal transduction pathways of LMX1B and its role in the pathogenesis of glomerulosclerosis.