The genetic risk factors for cerebral venous thrombosis: a case-control study in a Chinese national comprehensive hospital.

BACKGROUND: About 13-25% of cerebral venous thrombosis (CVT) cases lack clear etiology, which may be associated with underlying genetic factors. This study aims to investigate genetic factors in CVT patients using whole exome sequencing (WES). METHODS: Thirty-eight CVT patients hospitalized underwent WES. 977 subjects with WES data from a community cohort study --the Shunyi cohort were as the control group. Using bioinformatics analysis, differential genes with rare damaging variants between two groups were filtered (P < 0.05). KEGG enrichment analysis was performed on the screened genes to identify pathways associated with CVT. RESULTS: Through analysis of medical history, routine tests, and imaging examinations, the etiology of 38 patients: 8 cases of antiphospholipid syndrome, 6 cases with hematologic diseases, 3 cases of protein C deficiency, and 2 cases of protein S deficiency. Five cases occurred during pregnancy or puerperium, and 3 cases had a history of oral contraceptive use, and so on. The etiology was unknown in 12 cases (31.6%), and the etiology of 4 patients were further clarified through WES: F9 c.838 + 1_838 + 16del, Hemizygote: F9 EX1-EX7 Dup; CBS c.430G > A, CBS c.949 A > G; F2 c.1787G > A; SERPINC1 c.409-11G > T. Comparing the WES data of two groups, a total of 179 different genes with rare damaging variants were screened (P < 0.05), with 5 genes of interest (JAK2, C3, PROC, PROZ, SERPIND1). Enrichment analysis of the 179 different genes revealed the complement and coagulation pathway and the mitogen activated protein kinases (MAPK) pathway were associated with CVT. CONCLUSION: For CVT patients with unknown etiology, WES could help identify the cause of CVT early, which is of great significance for treatment decisions and prognosis. In addition to the complement and coagulation pathway, MAPK pathway is associated with CVT, potentially related to platelet regulation and inflammatory response.

Transcriptomic and network analysis identifies shared and unique pathways and immune changes across fibrotic interstitial lung diseases.

BACKGROUND: Interstitial lung disease (ILD) encompasses a diverse group of disorders characterized by chronic inflammation and fibrosis of the pulmonary interstitium. Three ILDs, namely idiopathic pulmonary fibrosis (IPF), fibrotic hypersensitivity pneumonitis (fHP), and connective tissue disease-associated ILD (CTD-ILD), exhibit similar progressive fibrosis phenotypes, yet possess distinct etiologies, encouraging us to explore their different underlying mechanisms. METHODS: Transcriptome data of fibrotic lung tissues from patients with IPF, fHP, and CTD-ILD were subjected to functional annotation, network, and pathway analyses. Additionally, we employed the xCell deconvolution algorithm to predict immune cell infiltration in patients with fibrotic ILDs and healthy controls. RESULTS: We identified a shared progressive fibrosis-related module in these diseases which was related to extracellular matrix (ECM) degradation and production and potentially regulated by the p53 family transcription factors. In IPF, neuron-related processes emerged as a critical specific mechanism in functional enrichment. In fHP, we observed that B cell signaling and immunoglobulin A (IgA) production may act as predominant processes, which was further verified by B cell infiltration and the central role of CD19 gene. In CTD-ILD, active chemokine processes were enriched, and active dendritic cells (aDCs) were predicted to infiltrate the lung tissues. CONCLUSIONS: This study revealed shared and specific molecular and cellular pathways among IPF, fHP, and CTD-ILD, providing a basis for understanding their pathogenesis and identifying potential therapeutic targets.

RNF213 in moyamoya disease: Genotype-phenotype association and the underlying mechanism.

ABSTRACT: Moyamoya disease (MMD) is a cerebrovascular disorder characterized by a steno-occlusive internal carotid artery and compensatory vascular network formation. Although the precise pathogenic mechanism remains elusive, genetic association studies have identified RNF213 as the principal susceptibility gene for MMD, with the single nucleotide polymorphism p.R4810K recognized as the founder variant predominantly in the Asian populations. Distinct genotype-phenotype correlations are observable in RNF213-related MMD. The clinical manifestations linked to p.R4810K bear commonalities within Asian cohort, including familial predisposition, earlier age of onset, ischemic episodes, and involvement of the posterior cerebral artery (PCA). However, despite these shared phenotypic characteristics, there is significant heterogeneity in RNF213-related MMD presentations. This diversity manifests as variations across ethnic groups, inconsistent clinical symptoms and prognosis, and occurrence of other vasculopathies involving RNF213. This heterogeneity, in conjunction with the observed low disease penetrance of RNF213 mutations, suggests that the presence of these mutations may not be sufficient to cause MMD, underscoring the potential influence of other genetic or environmental factors. Although the current research might not have fully identified these additional contributors, experimental evidence points toward the involvement of RNF213 in angiogenesis, lipid metabolism, and the immune response. Future research is required to unveil the molecular mechanisms and identify the factors that synergize with RNF213 in the pathogenesis of MMD.

A large pedigree study confirmed the CGG repeat expansion of RILPL1 Is associated with oculopharyngodistal myopathy.

BACKGROUND: Oculopharyngodistal myopathy (OPDM) is an autosomal dominant adult-onset degenerative muscle disorder characterized by ptosis, ophthalmoplegia and weakness of the facial, pharyngeal and limb muscles. Trinucleotide repeat expansions in non-coding regions of LRP12, G1PC1, NOTCH2NLC and RILPL1 were reported to be the etiologies for OPDM. RESULTS: In this study, we performed long-read whole-genome sequencing in a large five-generation family of 156 individuals, including 21 patients diagnosed with typical OPDM. We identified CGG repeat expansions in 5'UTR of RILPL1 gene in all patients we tested while no CGG expansion in unaffected family members. Repeat-primed PCR and fluorescence amplicon length analysis PCR were further confirmed the segregation of CGG expansions in other family members and 1000 normal Chinese controls. Methylation analysis indicated that methylation levels of the RILPL1 gene were unaltered in OPDM patients, which was consistent with previous studies. Our findings provide evidence that RILPL1 is associated OPDM in this large pedigree. CONCLUSIONS: Our results identified RILPL1 is the associated the disease in this large pedigree.

Association Study Identified HLA-DQA1 as a Novel Genetic Risk of Systemic Lupus Erythematosus-Associated Pulmonary Arterial Hypertension.

OBJECTIVE: Pulmonary arterial hypertension (PAH) is a severe complication of systemic lupus erythematosus (SLE). However, the genetic signatures of SLE-associated PAH have not been well studied. We aimed to identify genetic variants implicated in SLE-associated PAH susceptibility within the major histocompatibility complex (MHC) region and assess the contribution to clinical outcomes. METHODS: A total of 172 patients with SLE-associated PAH confirmed by right heart catheterization, 1,303 patients with SLE without PAH, and 9,906 healthy controls were included. Deep sequencing of the MHC region was performed to identify alleles, single-nucleotide polymorphisms, and amino acids. We compared patients with SLE-associated PAH with patients with SLE without PAH and healthy controls. Clinical association study was conducted to explore the contribution to phenotypes. RESULTS: A total of 19,881 genetic variants were identified within the MHC region. HLA-DQA1*03:02 was identified as a novel genetic variant associated with SLE-associated PAH in the discovery cohort (P = 5.68 × 10-12 ) and authenticated in an independent replication cohort (P = 1.30 × 10-9 ). The strongest associated amino acid position was mapped to HLA-DQα1 in the region affecting MHC/peptide-CD4+ T cell receptor affinity and antigen binding. Clinical association study demonstrated that patients with SLE-associated PAH with HLA-DQA1*03:02 had significantly lower rates of target role achievement (P = 0.005) and survival (P = 0.04). CONCLUSION: This study, based on the largest cohort of SLE-associated PAH, is the first to investigate how MHC region genetic variants contribute to SLE-associated PAH susceptibility. HLA-DQA1*03:02 is a novel genetic risk factor and a prognostic factor in SLE-associated PAH. Patients with SLE with this allele require regular monitoring and careful follow-up for early diagnosis and interventions for potential PAH.

A-to-I RNA Editing in Klebsiella pneumoniae Regulates Quorum Sensing and Affects Cell Growth and Virulence.

Millions of adenosine (A) to inosine (I) RNA editing events are reported and well-studied in eukaryotes; however, many features and functions remain unclear in prokaryotes. By combining PacBio Sequel, Illumina whole-genome sequencing, and RNA Sequencing data of two Klebsiella pneumoniae strains with different virulence, a total of 13 RNA editing events are identified. The RNA editing event of badR is focused, which shows a significant difference in editing levels in the two K. pneumoniae strains and is predicted to be a transcription factor. A hard-coded Cys is mutated on DNA to simulate the effect of complete editing of badR. Transcriptome analysis identifies the cellular quorum sensing (QS) pathway as the most dramatic change, demonstrating the dynamic regulation of RNA editing on badR related to coordinated collective behavior. Indeed, a significant difference in autoinducer 2 activity and cell growth is detected when the cells reach the stationary phase. Additionally, the mutant strain shows significantly lower virulence than the WT strain in the Galleria mellonella infection model. Furthermore, RNA editing regulation of badR is highly conserved across K. pneumoniae strains. Overall, this work provides new insights into posttranscriptional regulation in bacteria.

Integrated bioinformatics analysis identifies shared immune changes between ischemic stroke and COVID 19.

Although COVID-19 is primarily a respiratory disease, its neurological complications, such as ischemic stroke (IS), have aroused growing concerns and reports. However, the molecular mechanisms that underlie IS and COVID-19 are not well understood. Therefore, we implemented transcriptomic analysis from eight GEO datasets consist of 1191 samples to detect common pathways and molecular biomarkers in IS and COVID-19 that help understand the linkage between them. Differentially expressed genes (DEGs) were detected for IS and COVID-19 separately for finding shared mechanisms and we found that immune-related pathways were outlined with statistical significance. JAK2, which was identified as a hub gene, was supposed to be a potential therapeutic gene targets during the immunological process of COVID-19 and IS. Besides, we found a decrease in the proportion of CD8+ T and T helper 2 cells in the peripheral circulation of both COVID and IS patients, and NCR3 expression was significantly correlated with this change. In conclusion, we demonstrated that transcriptomic analyses reported in this study could make a deeper understanding of the common mechanism and might be promising for effective therapeutic for IS and COVID-19.

Transcriptomic analysis identifies shared biological foundations between ischemic stroke and Alzheimer's disease.

Aim: Alzheimer's disease (AD) and ischemic stroke (IS), two major neurological diseases, are suggested to be associated in clinical and pathophysiological levels. Previous studies have provided some insights into the possible genetic mechanisms behind the correlation between AD and IS, but this issue is still not clear. We implemented transcriptomic analysis to detect common hub genes and pathways to help promote the understanding of this issue. Materials and methods: Four gene expression profiling datasets (GSE16561, GSE58294, GSE63060, and GSE63061) of peripheral whole blood, which contain 108 IS samples, 284 AD samples, and 285 matched controls, were employed to detect differentially expressed genes (DEGs) for AD and IS, which were further analyzed for shared biological pathways, candidate drugs, and transcription factors. Protein-protein interaction (PPI) network and drug-target interaction analysis were applied to identify hub genes and drug targets, respectively. Result verification was done with other independent datasets (GSE37587, GSE46480, and GSE140829). The difference in proportions of various immune cells in the peripheral blood of AD and IS patients were evaluated using CIBERSORT. Results: We identified 74 DEGs and 18 biological processes with statistical significance shared by AD and IS, 9 of which were immune-related pathways. Five hub genes scored high in the topological analysis of the PPI network, and we also found eight drug target genes and candidate drugs which were associated with AD and IS. As for immunological changes, an increase in the proportion of M0 macrophages was found in the peripheral circulation of both AD and IS patients, and SOD1 expression was significantly correlated with this change. Conclusion: Collectively, the common DEGs and shared pathways found in this study suggest a potential shared etiology between AD and IS, behind which immune system, particularly the M0 macrophage elevation, might have important roles. While, the shared hub genes, potential therapeutic gene targets and drugs reported in this study provide promising treatment strategies for AD and IS.

Genotype-Phenotype Correlation Analysis and Identification of a Novel SRD5A2 Mutation in Four Unrelated Chinese Patients with 5α-Reductase Deficiency.

Objective: The 5α-reductase type 2 deficiency is mainly caused by mutations in the SRD5A2 gene. Our study aims to investigate the SRD5A2 gene mutations and their corresponding manifestations. Methods: Four unrelated Chinese patients with 46, XY ambiguous genitalia were studied. Molecular genetic alterations and clinical presentations were analyzed. Results: Five variants in the SRD5A2 gene were identified, all highly conserved in vertebrate orthologs. The p.P251A was a novel variant, predicted to "Affect protein function" and to be "probably damaging". Combining patients' gene mutations with their external genitalia and male sexual characteristics, we found that three variants, p.Q6X, p.N193S, and p.H90Y, were associated with severe undervirilization of external genitalia, and the other two, p.G203S and p.P251A, probably retained part of the enzyme activity. Conclusion: Mutation analysis of SRD5A2 gene is crucial for differential diagnosis in patients with 5α-reductase type 2 deficiency. Patients' variable manifestations depend on the mutation type and residual enzyme activity. The novel variant p.P251A enlarges the spectrum of SRD5A2 mutations.

Differential Regulation of the Immune System in Peripheral Blood Following Ischemic Stroke.

Method: 108 IS samples and 47 matched controls were obtained from the GEO database. Immune-related genes (IRGs) and their associated drugs were collected from the ImmPort and PharmGBK databases, respectively. Random forest (RF) regression and least absolute shrinkage and selection operator (LASSO) logistic regression were applied to identify immune-related genetic biomarkers (IRGBs) of IS, and accuracy was verified using neural network models. Finally, proportion changes of various immune cells in peripheral blood of IS patients were evaluated using CIBERSORT and xCell and correlation analyses were performed between IRGBs and differentially distributed immune cells. Results: A total of 537 genes were differentially expressed between IS and control samples. Four immune-related differential expressed genes identified by regression analysis presented strong predictive power (AUC = 0.909) which we suggeseted them as immune-related genetic biomarkers (IRGBs). We also demonstrated six immune-related genes targeted by known drugs. In addition, post-IS immune system presented an increase in the proportion of innate immune cells and a decrease in adaptive immune cells in the peripheral circulation, and IRGBs showing significance were associated with this process. Conclusion: The study identified CARD11, ICAM2, VIM, and CD19 as immune-related genetic biomarkers of IS. Six immune-related DEGs targeted by known drugs were found and provide new candidate drug targets for modulating the post-IS immune system. The innate immune cells and adaptive immune cells are diversified in the post-IS immune system, and IRGBs might play important role during this process.

Investigating the Genetic Characteristics of Hippocampal Volume and Plasma ß-Amyloid in a Chinese Community-Dwelling Population.

BACKGROUND AND OBJECTIVES: The genetic characteristics and correlations of hippocampal volume and plasm Aß, probable endophenotypes for dementia, remain to be explored in Chinese community cohort. Using whole-exome sequencing and SNP-array genotyping, we sought to identify rare and common variants and genes influencing these two endophenotypes, and calculate their heritability and genetic correlation. METHODS: Association analyses with both whole-exome sequencing and SNP-array genotyping data were performed for hippocampal volumes and plasm Aß with mixed-effect linear regression model adjusted for sex, age, and total intracranial volume or APOE ε4 while considering familial relatedness. We also performed gene-level analysis for common and gene-burden analysis for rare variants. Heritability and genetic correlation were further examined. RESULTS: Totally 1,261 participants from a Chinese community cohort were included and we identified one gene, PTPRT, for hippocampal volume, with the top significant SNPs by whole genome-wide association study. rs6030076 (P=5.48×10-8, ß=-0.092, SE=0.017) from whole-exome sequencing and rs6030088 (P=8.24×10-9, ß=-105.22 SE=18.09) from SNP-array data, both located in this gene. Gene-burden analysis based on rare mutations detected 6 genes to be significantly associated with Aß. The SNP-based heritability was 0.43±0.13 for hippocampal volume and 0.2-0.3 for plasma Aß. The SNP-based genetic correlation between hippocampal volume and plasma Aß were negative. DISCUSSION: In this study, we identified several SNPs and one gene, PTPRT, which were not reported in previous GWASs, associated with hippocampal volume. Besides, the heritability and the genetic correlation gave an overview of hippocampal volume and plasma Aß. Our findings provide insights into the mechanisms behind the individual variances in these endophenotypes.

Exome sequencing reveals genetic architecture in patients with isolated or syndromic short stature.

Short stature is among the most common endocrinological disease phenotypes of childhood and may occur as an isolated finding or in conjunction with other clinical manifestations. Although the diagnostic utility of clinical genetic testing in short stature has been implicated, the genetic architecture and the utility of genomic studies such as exome sequencing (ES) in a sizable cohort of patients with short stature have not been investigated systematically. In this study, we recruited 561 individuals with short stature from two centers in China during a 4-year period. We performed ES for all patients and available parents. All patients were retrospectively divided into two groups: an isolated short stature group (group I, n = 257) and an apparently syndromic short stature group (group II, n = 304). Causal variants were identified in 135 of 561 (24.1%) patients. In group I, 29 of 257 (11.3%) of the patients were solved by variants in 24 genes. In group II, 106 of 304 (34.9%) patients were solved by variants in 57 genes. Genes involved in fundamental cellular process played an important role in the genetic architecture of syndromic short stature. Distinct genetic architectures and pathophysiological processes underlie isolated and syndromic short stature.

Phenotypic and genetic spectrum of isolated macrodactyly: somatic mosaicism of PIK3CA and AKT1 oncogenic variants.

BACKGROUND: Isolated macrodactyly is a severe congenital hand anomaly with functional and physiological impact. Known causative genes include PIK3CA, AKT1 and PTEN. The aim of this study is to gain insights into the genetics basis of isolated macrodactyly. RESULTS: We enrolled 24 patients with isolated macrodactyly. Four of them were diagnosed with Proteus syndrome based on skin presentations characteristic to this disease. Targeted next-generation sequencing was performed using patients' blood and affected tissues. Overall, 20 patients carry mosaic PIK3CA pathogenic variants, i.e. p.His1047Arg (N = 7), p.Glu542Lys (N = 6), p.Glu545Lys (N = 2), p.His1047Leu (N = 2), p.Glu453Lys (N = 1), p.Gln546Lys (N = 1) and p.His1047Tyr (N = 1). Four patients who met the diagnostic criteria of Proteus syndrome carry mosaic AKT1 p.Glu17Lys variant. Variant allele frequencies of these mosaic variants obtained through next-generation sequencing range from 10 to 33%. In genotype-phenotype correlation analysis of patients with PIK3CA variant, we found that patients with the macrodactyly of the lower limbs tend to carry PIK3CA variants located in the helical domain (P = 0.005). CONCLUSIONS: Mosaic PIK3CA and AKT1 variants can be found in all of our samples with isolated macrodactyly. Insights into phenotypic and genetic spectrum of isolated macrodactyly may be helpful in perusing a more precise and effective management of isolated macrodactyly.

Mutational landscape and genetic signatures of cell-free DNA in tumour-induced osteomalacia.

Tumour-induced osteomalacia (TIO) is a very rare paraneoplastic syndrome with bone pain, fractures and muscle weakness, which is mostly caused by phosphaturic mesenchymal tumours (PMTs). Cell-free DNA (cfDNA) has been regarded as a non-invasive liquid biopsy for many malignant tumours. However, it has not been studied in benign tumours, which prompted us to adopt the targeted next-generation sequencing approach to compare cfDNAs of 4 TIO patients, four patients with bone metastasis (BM) and 10 healthy controls. The mutational landscapes of cfDNA in TIO and BM groups were similar in the spectrum of allele frequencies and mutation types. Markedly, deleterious missense mutations in FGFR1 and loss-of-function mutations in MED12 were found in 3/4 TIO patients but none of BM patients. The gene ontology analysis strongly supported that these mutated genes found in TIOs would play a potential role in PMTs' process. The genetic signatures and corresponding change in expression of FGFR1 and FGF23 were further validated in PMT tissues from a test cohort of another three TIO patients. In summary, we reported the first study of the mutational landscape and genetic signatures of cfDNA in TIO/PMTs.

TBX6 missense variants expand the mutational spectrum in a non-Mendelian inheritance disease.

Congenital scoliosis (CS) is a birth defect with variable clinical and anatomical manifestations due to spinal malformation. The genetic etiology underlying about 10% of CS cases in the Chinese population is compound inheritance by which the gene dosage is reduced below that of haploinsufficiency. In this genetic model, the trait manifests as a result of the combined effect of a rare variant and common pathogenic variant allele at a locus. From exome sequencing (ES) data of 523 patients in Asia and two patients in Texas, we identified six TBX6 gene-disruptive variants from 11 unrelated CS patients via ES and in vitro functional testing. The in trans mild hypomorphic allele was identified in 10 of the 11 subjects; as anticipated these 10 shared a similar spinal deformity of hemivertebrae. The remaining case has a homozygous variant in TBX6 (c.418C>T) and presents a more severe spinal deformity phenotype. We found decreased transcriptional activity and abnormal cellular localization as the molecular mechanisms for TBX6 missense loss-of-function alleles. Expanding the mutational spectrum of TBX6 pathogenic alleles enabled an increased molecular diagnostic detection rate, provided further evidence for the gene dosage-dependent genetic model underlying CS, and refined clinical classification.

Exome sequencing reveals a novel variant in NFX1 causing intracranial aneurysm in a Chinese family.

BACKGROUND: Genetic risk factors play an important role in the pathogenesis of familial intracranial aneurysms (FIAs); however, the molecular mechanisms remain largely unknown. OBJECTIVE: To investigate potential FIA-causing genetic variants by rare variant interrogation and a family-based genomics approach in a large family with an extensive multigenerational pedigree with FIAs. METHOD: Exome sequencing (ES) was performed in a dominant likely family with intracranial aneurysms (IAs). Variants were analyzed by an in-house developed pipeline and prioritized using various filtering strategies, including population frequency, variant type, and predicted variant pathogenicity. Sanger sequencing was also performed to evaluate the segregation of the variants with the phenotype. RESULTS: Based on the ES data obtained from five individuals from a family with 7/21 living members affected with IAs, a total of 14 variants were prioritized as candidate variants. Familial segregation analysis revealed that NFX1 c.2519T>C (p.Leu840Pro) segregated in accordance with Mendelian expectations with the phenotype within the family-that is, present in all IA-affected cases and absent from all unaffected members of the second generation. This missense variant is absent from public databases (1000genome, ExAC, gnomAD, ESP5400), and has damaging predictions by bioinformatics tools (Gerp ++ score = 5.88, CADD score = 16.43, MutationTaster score = 1, LRT score = 0). In addition, 840Leu in NFX1 is robustly conserved in mammals and maps in a region before the RING-type zinc finger domain. CONCLUSION: NFX1 c.2519T>C (p.Leu840Pro) may contribute to the pathogenetics of a subset of FIAs.

Genetic and molecular mechanism for distinct clinical phenotypes conveyed by allelic truncating mutations implicated in FBN1.

BACKGROUND: The molecular and genetic mechanisms by which different single nucleotide variant alleles in specific genes, or at the same genetic locus, cause distinct disease phenotypes often remain unclear. Allelic truncating mutations of FBN1 could cause either classical Marfan syndrome (MFS) or a more complicated phenotype associated with Marfanoid-progeroid-lipodystrophy syndrome (MPLS). METHODS: We investigated a small cohort, encompassing two classical MFS and one MPLS subjects from China, whose clinical presentation included scoliosis potentially requiring surgical intervention. Targeted next generation sequencing was performed on all the participants. We analyzed the molecular diagnosis, clinical features, and the potential molecular mechanism involved in the MPLS subject in our cohort. RESULTS: We report a novel de novo FBN1 mutation for the first Chinese subject with MPLS, a more complicated fibrillinopathy, and two subjects with more classical MFS. We further predict that the MPLS truncating mutation, and others previously reported, is prone to escape the nonsense-mediated decay (NMD), while MFS mutations are predicted to be subjected to NMD. Also, the MPLS mutation occurs within the glucogenic hormone asprosin domain of FBN1. In vitro experiments showed that the single MPLS mutation p.Glu2759Cysfs*9 appears to perturb proper FBN1 protein aggregation as compared with the classical MFS mutation p.Tyr2596Thrfs*86. Both mutations appear to upregulate SMAD2 phosphorylation in vitro. CONCLUSION: We provide direct evidence that a dominant-negative interaction of FBN1 potentially explains the complex MPLS phenotypes through genetic and functional analysis. Our study expands the mutation spectrum of FBN1 and highlights the potential molecular mechanism for MPLS.

Identification of novel FBN1 variations implicated in congenital scoliosis.

Congenital scoliosis (CS) is a form of scoliosis caused by congenital vertebral malformations. Genetic predisposition has been demonstrated in CS. We previously reported that TBX6 loss-of-function causes CS in a compound heterozygous model; however, this model can explain only 10% of CS. Many monogenic and polygenic CS genes remain to be elucidated. In this study, we analyzed exome sequencing (ES) data of 615 Chinese CS from the Deciphering Disorders Involving Scoliosis and COmorbidities (DISCO) project. Cosegregation studies for 103 familial CS identified a novel heterozygous nonsense variant, c.2649G>A (p.Trp883Ter) in FBN1. The association between FBN1 and CS was then analyzed by extracting FBN1 variants from ES data of 574 sporadic CS and 828 controls; 30 novel variants were identified and prioritized for further analyses. A mutational burden test showed that the deleterious FBN1 variants were significantly enriched in CS subjects (OR = 3.9, P = 0.03 by Fisher's exact test). One missense variant, c.2613A>C (p.Leu871Phe) was recurrent in two unrelated CS subjects, and in vitro functional experiments for the variant suggest that FBN1 may contribute to CS by upregulating the transforming growth factor beta (TGF-ß) signaling. Our study expanded the phenotypic spectrum of FBN1, and provided nove insights into the genetic etiology of CS.

Evolutionarily significant A-to-I RNA editing events originated through G-to-A mutations in primates.

BACKGROUND: Recent studies have revealed thousands of A-to-I RNA editing events in primates, but the origination and general functions of these events are not well addressed. RESULTS: Here, we perform a comparative editome study in human and rhesus macaque and uncover a substantial proportion of macaque A-to-I editing sites that are genomically polymorphic in some animals or encoded as non-editable nucleotides in human. The occurrence of these recent gain and loss of RNA editing through DNA point mutation is significantly more prevalent than that expected for the nearby regions. Ancestral state analyses further demonstrate that an increase in recent gain of editing events contribute to the over-representation, with G-to-A mutation site as a favorable location for the origination of robust A-to-I editing events. Population genetics analyses of the focal editing sites further reveal that a portion of these young editing events are evolutionarily significant, indicating general functional relevance for at least a fraction of these sites. CONCLUSIONS: Overall, we report a list of A-to-I editing events that recently originated through G-to-A mutations in primates, representing a valuable resource to investigate the features and evolutionary significance of A-to-I editing events at the population and species levels. The unique subset of primate editome also illuminates the general functions of RNA editing by connecting it to particular gene regulatory processes, based on the characterized outcome of a gene regulatory level in different individuals or primate species with or without these editing events.

Whole-Genome Methylation Analysis of Phenotype Discordant Monozygotic Twins Reveals Novel Epigenetic Perturbation Contributing to the Pathogenesis of Adolescent Idiopathic Scoliosis.

Background: Adolescent idiopathic scoliosis (AIS) is a complex disease affecting a large number of teenagers, especially in female. This study reveals novel epigenetic perturbation to the pathogenesis of AIS. Methods: A female monozygotic (MZ) twin pair discordant for AIS were examined for whole-exome sequencing and epigenome difference. Sets of differentially methylated regions (DMRs) were validated using MethylTarget™ method in 20 AIS female patients and 20 healthy female controls. Results: Few exome difference but several potential DMRs were found between the MZ twins. We identified 313 hypermethylated DMRs and 397 hypomethylated DMRs, respectively. Most of them were enriched in the MAPK and PI3K-Akt signaling pathway, which may contribute to the discordance of AIS. Several DMRs related to scoliosis genes were tested, and the NDN: TSS-DMR (chr15:23932133-23932304, hg19) was confirmed in additional samples. The methylation level of this DMR was significantly higher in the AIS group than in the control group (p = 0.04). Conclusions: We described the epigenome difference in an AIS female discordant MZ twin pair using Whole Genome Bisulfite Sequencing (WGBS). The NDN: TSS-DMR had higher methylation level in female AIS, which can help elucidate the potential etiology of AIS.