An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases.

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Pharmacogenetic implementation for CYP2C19 and pharmacokinetics of voriconazole in children with malignancy or inborn errors of immunity.

BACKGROUND: Voriconazole pharmacokinetics (PK) are known to be affected by genetic polymorphisms of drug-metabolizing enzymes such as CYP2C19; however, such information is limited for the pediatric population. The primary aim of this study is to establish a voriconazole PK model incorporating CYP2C19 phenotypes in Japanese children with malignancy or inborn errors of immunity. METHODS: CYP2C19 genotypes were assessed by whole-genome genotyping and defined as follows: *17/*17: ultrarapid metabolizer (URM), *1/*17: rapid metabolizer (RM), *1/*1:normal metabolizer (NM), *1/*2, *1/*3, *2/*17:intermediate metabolizer (IM), and *2/*2, *2/*3, *3/*3: poor metabolizer (PM). Population PK analysis was performed. The voriconazole serum concentration profile was described by a two-compartment model with first-order absorption, mixed linear and nonlinear (Michaelis-Menten) elimination. RESULTS: Voriconazole concentration data were available from 60 patients with a median age of 5.3 years. The phenotypes predicted from CYP2C19 genotypes were RM in 1 (2 %), NM in 21 (35 %) patients, IM in 27 (45 %) patients, and PM in 11 (18 %) patients. Underlying diseases included 38 (63%) patients with hematological malignancy and 18 (30 %) patients with inborn errors of immunity. Among the CYP2C19 phenotypes, PM was predicted to show complete inhibition (the degree of Vmax inhibition [Vmax, inh] = 100 %; Vmax = 0). The estimated parameters of Vmax,inh were +0.8 higher in patients with gamma-glutamyl transpeptidase (γ-GTP) Grade 2 or higher and +2.7 higher when C-reactive protein (CRP) levels were 2.0 mg/dL or higher. CONCLUSION: CYP2C19 genetic polymorphisms, γ-GTP, and CRP affect Vmax,inh of voriconazole in children with malignancy or inborn errors of immunity.

A longitudinal population-based study identifies THBS2 as a susceptibility gene for intervertebral disc degeneration.

PURPOSE: Intervertebral disc degeneration (IDD) is a common degenerative disease associated with ageing. Additionally, IDD is recognized as one of the leading causes of low back pain and disability in the working-age population and is the first step in the process leading to degenerative spinal changes. However, the genetic factors and regulatory mechanisms of IDD remain unknown. Therefore, we selected eight single nucleotide polymorphisms of genes to reveal the progression of IDD in a 7-year longitudinal study of the general population in Japan. METHODS: IDD was evaluated in the Wakayama Spine Study (WSS), which is a population-based cohort study. Overall, 574 participants from the general population cohort who underwent whole spine magnetic resonance imaging and provided clinical information were included in this longitudinal survey. RESULTS: The progression of IDD was affected only by THBS2 at the lumbar region, T12-L1 (p = 0.0044) and L3-4 (p = 0.0045). The significant interaction between THBS2 and age with IDD negatively affected the thoracic spines and passively influenced both the thoracolumbar junction and thoracic spines. The higher progression per year of Pfirrmann's score was rapid in young people with age; however, this decelerated the IDD progression per year in different ages. CONCLUSION: Our longitudinal study found the genes associated with IDD progression and that genetic factors' impact on IDD differs depending on disc level and age.

RNF213 Variants, Vasospastic Angina, and Risk of Fatal Myocardial Infarction.

Importance: Vasospastic angina (VSA) is vasospasm of the coronary artery and is particularly prevalent in East Asian populations. However, the specific genetic architecture for VSA at genome-wide levels is not fully understood. Objective: To identify genetic factors associated with VSA. Design, Setting, and Participants: This was a case-control genome-wide association study of VSA. Data from Biobank Japan (BBJ; enrolled patients from 2002-2008 and 2013-2018) were used, and controls without coronary artery disease (CAD) were enrolled. Patients from the BBJ were genotyped using arrays or a set of arrays. Patients recruited between 2002 and 2005 were classified within the first dataset, and those recruited between 2006 and 2008 were classified within the second dataset. To replicate the genome-wide association study in the first and second datasets, VSA cases and control samples from the latest patients in the BBJ recruited between 2013 and 2018 were analyzed in a third dataset. Exposures: Single-nucleotide variants associated with VSA. Main Outcomes and Measures: Cases with VSA and controls without CAD. Results: A total of 5720 cases (mean [SD] age, 67 [10] years; 3672 male [64.2%]) and 153 864 controls (mean [SD] age, 62 [15] years; 77 362 male [50.3%]) in 3 datasets were included in this study. The variants at the RNF213 locus showed the strongest association with VSA across the 3 datasets (odds ratio [OR], 2.34; 95% CI, 1.99-2.74; P = 4.4 × 10-25). Additionally, rs112735431, an Asian-specific rare deleterious variant (p.Arg4810Lys) experimentally shown to be associated with reduced angiogenesis and a well-known causal risk for Moyamoya disease was the most promising candidate for a causal variant explaining the association. The effect size of rs112735431 on VSA was distinct from that of other CADs. Furthermore, homozygous carriers of rs112735431 showed an association with VSA characterized by a large effect estimate (OR, 18.34; 95% CI, 5.15-65.22; P = 7.0 × 10-6), deviating from the additive model (OR, 4.35; 95% CI, 1.18-16.05; P = .03). Stratified analyses revealed that rs112735431 exhibited a stronger association in males (χ21 = 7.24; P = .007) and a younger age group (OR, 3.06; 95% CI, 2.24-4.19), corresponding to the epidemiologic features of VSA. In the registry, carriers without CAD of the risk allele rs112735431 had a strikingly high mortality rate due to acute myocardial infarction during the follow-up period (hazard ratio, 2.71; 95% CI, 1.57-4.65; P = 3.3 × 10-4). As previously reported, a possible overlap between VSA and Moyamoya disease was not found. Conclusions and Relevance: Results of this study suggest that vascular cell dysfunction mediated by variants in the RNF213 locus may promote coronary vasospasm, and the presence of the risk allele could serve as a predictive factor for the prognosis.

Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

A genome-wide association study identifies a locus associated with knee extension strength in older Japanese individuals.

Sarcopenia is a common skeletal muscle disease in older people. Lower limb muscle strength is a good predictive value for sarcopenia; however, little is known about its genetic components. Here, we conducted a genome-wide association study (GWAS) for knee extension strength in a total of 3452 Japanese aged 60 years or older from two independent cohorts. We identified a significant locus, rs10749438 which is an intronic variant in TACC2 (transforming acidic coiled-coil-containing 2) (P = 4.2 × 10-8). TACC2, encoding a cytoskeleton-related protein, is highly expressed in skeletal muscle, and is reported as a target of myotonic dystrophy 1-associated splicing alterations. These suggest that changes in TACC2 expression are associated with variations in muscle strength in older people. The association was consistently observed in young and middle-aged subjects. Our findings would shed light on genetic components of lower limb muscle strength and indicate TACC2 as a potential therapeutic target for sarcopenia.

Anti-integrin αvß6 antibody in Takayasu arteritis patients with or without ulcerative colitis.

Background: It has been well documented that Takayasu arteritis (TAK) and ulcerative colitis (UC) coexist in the same patients. HLA-B*52 characterizes the co-occurrence, which is one of the common genetic features between these two diseases, indicating shared underlying pathologic mechanisms. Anti-integrin αvß6 antibody (Ab) is present in sera of UC patients in a highly specific manner. We investigated if there were any associations between anti-integrin αvß6 Ab and TAK, considering the risk HLA alleles. Methods: A total of 227 Japanese TAK patients were recruited in the current study and their serum samples were subjected to measurement of anti-integrin αvß6 Ab by ELISA. The clinical information, including the co-occurrence of UC, was collected. The HLA allele carrier status was determined by Luminex or genotype imputation. Results: The information about the presence of UC was available for 165 patients, among which eight (4.84%) patients had UC. Anti-integrin αvß6 antibody was identified in 7 out of 8 TAK subjects with UC (87.5%) while only 5 out of 157 (3.18%) TAK subjects without UC had the antibody (OR 121, p=7.46×10-8). A total of 99 out of 218 (45.4%) patients were HLA-B*52 carriers. There was no significant association between the presence of anti-integrin αvß6 Ab and HLA-B*52 carrier status in those without UC (OR 2.01, 95% CI 0.33-12.4, p = 0.189). Conclusions: The prevalence of anti-integrin αvß6 Ab was high in TAK patients with UC, but not in the absence of concomitant UC. The effect of HLA-B*52 on anti-integrin αvß6 Ab production would be minimal.

Annotation of nuclear lncRNAs based on chromatin interactions.

The human genome is pervasively transcribed and produces a wide variety of long non-coding RNAs (lncRNAs), constituting the majority of transcripts across human cell types. Some specific nuclear lncRNAs have been shown to be important regulatory components acting locally. As RNA-chromatin interaction and Hi-C chromatin conformation data showed that chromatin interactions of nuclear lncRNAs are determined by the local chromatin 3D conformation, we used Hi-C data to identify potential target genes of lncRNAs. RNA-protein interaction data suggested that nuclear lncRNAs act as scaffolds to recruit regulatory proteins to target promoters and enhancers. Nuclear lncRNAs may therefore play a role in directing regulatory factors to locations spatially close to the lncRNA gene. We provide the analysis results through an interactive visualization web portal at https://fantom.gsc.riken.jp/zenbu/reports/#F6_3D_lncRNA.

Fundamentals for predicting transcriptional regulations from DNA sequence patterns.

Cell-type-specific regulatory elements, cataloged through extensive experiments and bioinformatics in large-scale consortiums, have enabled enrichment analyses of genetic associations that primarily utilize positional information of the regulatory elements. These analyses have identified cell types and pathways genetically associated with human complex traits. However, our understanding of detailed allelic effects on these elements' activities and on-off states remains incomplete, hampering the interpretation of human genetic study results. This review introduces machine learning methods to learn sequence-dependent transcriptional regulation mechanisms from DNA sequences for predicting such allelic effects (not associations). We provide a concise history of machine-learning-based approaches, the requirements, and the key computational processes, focusing on primers in machine learning. Convolution and self-attention, pivotal in modern deep-learning models, are explained through geometrical interpretations using dot products. This facilitates understanding of the concept and why these have been used for machine learning for DNA sequences. These will inspire further research in this genetics and genomics field.

Relationship of Major Depressive Disorder and Schizophrenia Polygenic Risk Scores to Suicide: A Comparison Between European and Asian Ancestry Populations.

Psychiatric diagnosis rates in suicide decedents appear higher in European ancestry populations compared with East Asians. Shared genetic components exist between major depressive disorder (MDD)/schizophrenia (SCZ) and suicide, but no study has compared these shared polygenic architectures between Europeans and East Asians. We compared polygenic risk scores (PRSs) for MDD/SCZ determined from large data sets specific to each ancestry in European and East Asian suicide decedent samples. MDD/SCZ PRSs appeared more prominent in European suicides compared with Japanese suicides. A greater coexistence of psychiatric disorders in European suicide decedents than in East Asian suicide decedents may be partly explained by genetics. Our results are limited by the smaller sample size of our suicide decedents and sample size disparities between the European discovery data set and the East Asian data set for MDD/SCZ, resulting in less statistical power to detect robust difference between the two ancestries.

Psychiatric diagnosis in suicides appears more common in Europeans than in East Asians.This is the first comparison of suicide genome-wide association studies between Europeans and East Asians.Major depressive disorder/schizophrenia polygenic risk scores for suicide were more significant for Europeans than for East Asians.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing.

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Comparative evaluation of SNVs, indels, and structural variations detected with short- and long-read sequencing data.

Short- and long-read sequencing technologies are routinely used to detect DNA variants, including SNVs, indels, and structural variations (SVs). However, the differences in the quality and quantity of variants detected between short- and long-read data are not fully understood. In this study, we comprehensively evaluated the variant calling performance of short- and long-read-based SNV, indel, and SV detection algorithms (6 for SNVs, 12 for indels, and 13 for SVs) using a novel evaluation framework incorporating manual visual inspection. The results showed that indel-insertion calls greater than 10 bp were poorly detected by short-read-based detection algorithms compared to long-read-based algorithms; however, the recall and precision of SNV and indel-deletion detection were similar between short- and long-read data. The recall of SV detection with short-read-based algorithms was significantly lower in repetitive regions, especially for small- to intermediate-sized SVs, than that detected with long-read-based algorithms. In contrast, the recall and precision of SV detection in nonrepetitive regions were similar between short- and long-read data. These findings suggest the need for refined strategies, such as incorporating multiple variant detection algorithms, to generate a more complete set of variants using short-read data.

Disease-specific variant interpretation highlighted the genetic findings in 2325 Japanese patients with retinitis pigmentosa and allied diseases.

BACKGROUND: As gene-specific therapy for inherited retinal dystrophy (IRD) advances, unified variant interpretation across institutes is becoming increasingly important. This study aims to update the genetic findings of 86 retinitis pigmentosa (RP)-related genes in a large number of Japanese patients with RP by applying the standardised variant interpretation guidelines for Japanese patients with IRD (J-IRD-VI guidelines) built upon the American College of Medical Genetics and Genomics and the Association for Molecular Pathology rules, and assess the contribution of these genes in RP-allied diseases. METHODS: We assessed 2325 probands with RP (n=2155, including n=1204 sequenced previously with the same sequencing panel) and allied diseases (n=170, newly analysed), including Usher syndrome, Leber congenital amaurosis and cone-rod dystrophy (CRD). Target sequencing using a panel of 86 genes was performed. The variants were interpreted according to the J-IRD-VI guidelines. RESULTS: A total of 3564 variants were detected, of which 524 variants were interpreted as pathogenic or likely pathogenic. Among these 524 variants, 280 (53.4%) had been either undetected or interpreted as variants of unknown significance or benign variants in our earlier study of 1204 patients with RP. This led to a genetic diagnostic rate in 38.6% of patients with RP, with EYS accounting for 46.7% of the genetically solved patients, showing a 9% increase in diagnostic rate from our earlier study. The genetic diagnostic rate for patients with CRD was 28.2%, with RP-related genes significantly contributing over other allied diseases. CONCLUSION: A large-scale genetic analysis using the J-IRD-VI guidelines highlighted the population-specific genetic findings for Japanese patients with IRD; these findings serve as a foundation for the clinical application of gene-specific therapies.

Protein-altering variants at copy number-variable regions influence diverse human phenotypes.

Copy number variants (CNVs) are among the largest genetic variants, yet CNVs have not been effectively ascertained in most genetic association studies. Here we ascertained protein-altering CNVs from UK Biobank whole-exome sequencing data (n = 468,570) using haplotype-informed methods capable of detecting subexonic CNVs and variation within segmental duplications. Incorporating CNVs into analyses of rare variants predicted to cause gene loss of function (LOF) identified 100 associations of predicted LOF variants with 41 quantitative traits. A low-frequency partial deletion of RGL3 exon 6 conferred one of the strongest protective effects of gene LOF on hypertension risk (odds ratio = 0.86 (0.82-0.90)). Protein-coding variation in rapidly evolving gene families within segmental duplications-previously invisible to most analysis methods-generated some of the human genome's largest contributions to variation in type 2 diabetes risk, chronotype and blood cell traits. These results illustrate the potential for new genetic insights from genomic variation that has escaped large-scale analysis to date.

Fine-mapping genomic loci refines bipolar disorder risk genes.

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Attenuation of HOIL-1L ligase activity promotes systemic autoimmune disorders by augmenting linear ubiquitin signaling.

Linear ubiquitin chains, which are generated specifically by the linear ubiquitin assembly complex (LUBAC) ubiquitin ligase, play crucial roles in immune signaling, including NF-κB activation. LUBAC comprises catalytic large isoform of heme-oxidized iron regulatory protein 2 ubiquitin ligase 1 (HOIL-1L) interacting protein (HOIP), accessory HOIL-1L, and SHANK-associated RH domain-interacting protein (SHARPIN). Deletion of the ubiquitin ligase activity of HOIL-1L, an accessory ligase of LUBAC, augments LUBAC functions by enhancing LUBAC-mediated linear ubiquitination, which is catalyzed by HOIP. Here, we show that HOIL-1L ΔRING1 mice, which exhibit augmented LUBAC functions upon loss of the HOIL-1L ligase, developed systemic lupus erythematosus (SLE) and Sjögren's syndrome in a female-dominant fashion. Augmented LUBAC activity led to hyperactivation of both lymphoid and myeloid cells. In line with the findings in mice, we sought to identify missense single nucleotide polymorphisms/variations of the RBCK1/HOIL-1L gene in humans that attenuate HOIL-1L ligase activity. We found that the R464H variant, which is encoded by rs774507518 within the RBCK1/HOIL-1L gene, attenuated HOIL-1L ligase activity and augmented LUBAC-mediated immune signaling, including that mediated by Toll-like receptors. We also found that rs774507518 was enriched significantly in patients with SLE, strongly suggesting that RBCK1/HOIL-1L is an SLE susceptibility gene and that augmented linear ubiquitin signaling generated specifically by LUBAC underlies the pathogenesis of this prototype systemic autoimmune disease.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases.

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

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.

Adolescent idiopathic scoliosis (AIS) is a twisting deformity of the spine that occurs during periods of rapid growth in children worldwide. Children with severe cases of AIS require surgery to stop it from getting worse, presenting a significant financial burden to health systems and families. Although AIS is known to cluster in families, its genetic causes and its inheritance pattern have remained elusive. Additionally, AIS is known to be more prevalent in females, a bias that has not been explained. Advances in techniques to study the genetics underlying diseases have revealed that certain variations that increase the risk of AIS affect cartilage and connective tissue. In humans, one such variation is near a gene called Pax1, and it is female-specific. The extracellular matrix is a network of proteins and other molecules in the space between cells that help connect tissues together, and it is particularly important in cartilage and other connective tissues. One of the main components of the extracellular matrix is collagen. Yu, Kanshour, Ushiki et al. hypothesized that changes in the extracellular matrix could affect the cartilage and connective tissues of the spine, leading to AIS. To show this, the scientists screened over 100,000 individuals and found that AIS is associated with variants in two genes coding for extracellular matrix proteins. One of these variants was found in a gene called Col11a1, which codes for one of the proteins that makes up collagen. To understand the relationship between Pax1 and Col11a1, Yu, Kanshour, Ushiki et al. genetically modified mice so that they would lack the Pax1 gene. In these mice, the activation of Col11a1 was reduced in the mouse spine. They also found that the form of Col11a1 associated with AIS could not suppress the activation of a gene called Mmp3 in mouse cartilage cells as effectively as unmutated Col11a1. Going one step further, the researchers found that lowering the levels of an estrogen receptor altered the activation patterns of Pax1, Col11a1, and Mmp3 in mouse cartilage cells. These findings suggest a possible mechanism for AIS, particularly in females. The findings of Yu, Kanshour, Ushiki et al. highlight that cartilage cells in the spine are particularly relevant in AIS. The results also point to specific molecules within the extracellular matrix as important for maintaining proper alignment in the spine when children are growing rapidly. This information may guide future therapies aimed at maintaining healthy spinal cells in adolescent children, particularly girls.

Family History of Developmental Dysplasia of the Hip is a Risk Factor for the Progression of Hip Osteoarthritis.

BACKGROUND: Developmental dysplasia of the hip (DDH) is considered to have genetic predisposition and presents many intrafamilial occurrences. However, there is no report that evaluates the effect of DDH family history on the progression after the onset of hip osteoarthritis (OA). METHODS: Medical interviews about detailed clinical information including family history were conducted on 298 consecutive patients who had undergone surgery for OA due to DDH. Clinical or radiographic items that are associated with the severity of DDH (total hip arthroplasty [THA], involvement of bilateral DDH, onset age of hip pain, and three radiological indices of DDH: center-edge angle, sharp angle, and acetabular roof obliquity) were collected and evaluated in multivariate analyses for their associations with DDH family history in a qualitative or quantitative manner. Survival time analyses for THA as the endpoint was also performed to evaluate the effects of DDH family history on the progression of OA. RESULTS: The DDH family history showed significant associations with bilateral involvement of DDH (odds ratio = 2.09 [95% confidence interval {CI} 1.05 to 4.16]; P = .037), early onset of hip pain (P = .0065), and radiological severity of DDH (P = .016). The DDH family history showed a significant association with undergoing THA (odds ratio = 2.25 [95% CI 1.09 to 4.66]; P = .029), further supported by the Cox regression analyses (hazards ratio = 1.56 [95% CI 1.15 to 2.11]; P = .0044). CONCLUSION: A DDH family history is a risk factor for the progression of hip OA. Stronger genetic predisposition to DDH leads to faster onset and progression of hip OA.