Genome-Wide Mendelian Randomization Study Reveals Druggable Genes for Cerebral Small Vessel Disease.

BACKGROUND: Cerebral small vessel disease (CSVD) is a group of neurological disorders that affect the small blood vessels within the brain, for which no effective treatments are currently available. We conducted a Mendelian randomization (MR) study to identify candidate therapeutic genes for CSVD. METHODS: We retrieved genome-wide association study data from 6 recently conducted, extensive investigations focusing on CSVD magnetic resonance imaging markers and performed a 2-sample MR analysis to assess the potential causal effects of gene expression and protein level within druggable genes on CSVD in blood and brain tissues. Colocalization analyses and repeat studies were undertaken to verify the relationship. Additionally, mediation analysis was conducted to explore the potential mechanisms involving druggable genes and known risk factors for CSVD. Finally, phenome-wide MR analyses were applied to evaluate the potential adverse effects related to the identified druggable genes for CSVD treatment. RESULTS: Overall, 5 druggable genes consistently showed associations with CSVD in MR analyses across both the discovery and validation cohorts. Notably, the ALDH2 and KLHL24 genes were identified as associated with CSVD in both blood and brain tissues, whereas the genes ADRB1, BTN3A2, and EFEMP1 were exclusively detected in brain tissue. Moreover, mediation analysis elucidated the proportion of the total effects mediated by CSVD risk factors through candidate druggable genes, which ranged from 5.5% to 18.5%, and offered potential explanations for the observed results. A comprehensive phenome-wide MR analysis further emphasized both the therapeutic benefits and potential side effects of targeting these candidate druggable genes. CONCLUSIONS: This study provides genetic evidence supporting the potential therapeutic benefits of targeting druggable genes for treating CSVD, which will be useful for prioritizing CSVD drug development.

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.

Correlation between Circle of Willis configuration and intracranial arterial dolichoectasia, and genetic contributions.

OBJECTIVES: Intracranial arterial dolichoectasia (IADE) is characterized by the dilation, elongation, and tortuosity of intracranial arteries. We aimed to investigate the association between variations of the Circle of Willis (COW) and IADE in the general population, as well as estimate the genetic correlation between COW variations and IADE. METHODS: A total of 981 individuals from a population-based cohort were included. Brain magnetic resonance angiography was performed to assess COW variants and measure the diameters of intracranial arteries. IADE was defined as a total intracranial volume-adjusted diameter ≥ 2 standard deviations. Logistic regression models were used to analyze the association between COW variations and IADE. The heritability and genetic correlation were estimated using genome-wide complex trait analysis (GCTA) based on single nucleotide polymorphism (SNP) array data. RESULTS: The prevalence of IADE was 6.2 %. Hypoplastic/absent A1 segments were associated with an increase in contralateral ICA diameter (ß ± SE, 0.279 ± 0.049; p = 0.001) and a decrease in ipsilateral ICA diameter (ß ± SE, -0.300 ± 0.050; p = 0.001). Fetal-type posterior cerebral artery (FTP) was associated with a larger ICA diameter (ß ± SE, 0.326 ± 0.048; p = 0.001) and a smaller BA diameter (ß ± SE, -0.662 ± 0.043; p = 0.001). FTP revealed a positive genetic correlation with ICA dilation (rG = 0.259 ± 0.175; p = 0.0009) and a negative genetic correlation with BA dilation (rG = -0.192 ± 0.153, p = 0.015). CONCLUSIONS: There was an association between COW variations and larger intracranial arterial diameters in the general population. Genetic factors may play a role in the development of intracranial arterial dilation and the formation of COW variants.

Rare RNF213 variants is related to early-onset intracranial atherosclerosis: A Chinese community-based study.

BACKGROUND: The relationship between rare variants in Ring finger protein 213 (RNF213) and intracranial atherosclerosis (ICAS) remained unelucidated. Using whole-exome sequencing (WES) and high-resolution magnetic resonance imaging (HR-MRI), this study aimed at investigating the association between rare RNF213 variants and ICAS within a Chinese community-dwelling population. METHODS: The present study included 821 participants from Shunyi cohort. Genetic data of rare RNF213 variants were acquired by WES and were categorized by functional domains. Intracranial and extracranial atherosclerosis were assessed by brain HR-MRI and carotid ultrasound, respectively. Logistic regression and generalized linear regression were applied to evaluate the effects of rare RNF213 variants on atherosclerosis. Stratification by age were conducted with 50 years old set as the cutoff value. RESULTS: Ninety-five participants were identified as carriers of rare RNF213 variants. Carotid plaques were observed in 367 (44.7 %) participants, while ICAS was identified in 306 (37.3 %). Rare variants of RNF213 was not associated with ECAS. Employing HR-MRI, both the presence of rare variants (ß = 0.150, P = 0.025) and numerical count of variants (ß = 0.182, P = 0.003) were significantly correlated with ICAS within the group of age ≤50 years. Both variant existence (ß = 0.154, P = 0.014) and variant count (ß = 0.188, P = 0.003) were significantly associated with plaques in middle cerebral arteries within younger subgroup, rather than basilar arteries. Furthermore, a significant association was observed between variants that located outside the N-arm domain and ICAS in the younger subgroup (OR = 2.522, P = 0.030). Statistical results remained robust after adjusted for age, gender, and cardiovascular risk factors. CONCLUSIONS: Rare variants of RNF213 is associated with age-related ICAS in general Chinese population, highlighting the potential role of RNF213 as a genetic contributor to early-onset ICAS.

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.

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.

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.

The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease.

With the recent advance in genome-wide association studies (GWAS), disease-associated single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) have been extensively reported. Accordingly, the issue of incorrect identification of recombination events that can induce the distortion of multi-allelic or hemizygous variants has received more attention. However, the potential distorted calculation bias or significance of a detected association in a GWAS due to the coexistence of CNVs and SNPs in the same genomic region may remain under-recognized. Here we performed the association study within a congenital scoliosis (CS) cohort whose genetic etiology was recently elucidated as a compound inheritance model, including mostly one rare variant deletion CNV null allele and one common variant non-coding hypomorphic haplotype of the TBX6 gene. We demonstrated that the existence of a deletion in TBX6 led to an overestimation of the contribution of the SNPs on the hypomorphic allele. Furthermore, we generalized a model to explain the calculation bias, or distorted significance calculation for an association study, that can be 'induced' by CNVs at a locus. Meanwhile, overlapping between the disease-associated SNPs from published GWAS and common CNVs (overlap 10%) and pathogenic/likely pathogenic CNVs (overlap 99.69%) was significantly higher than the random distribution (p < 1 × 10-6 and p = 0.034, respectively), indicating that such co-existence of CNV and SNV alleles might generally influence data interpretation and potential outcomes of a GWAS. We also verified and assessed the influence of colocalizing CNVs to the detection sensitivity of disease-associated SNP variant alleles in another adolescent idiopathic scoliosis (AIS) genome-wide association study. We proposed that detecting co-existent CNVs when evaluating the association signals between SNPs and disease traits could improve genetic model analyses and better integrate GWAS with robust Mendelian principles.

Fat Mass and Obesity-Associated (FTO) Gene Polymorphisms Are Associated with Risk of Intervertebral Disc Degeneration in Chinese Han Population: A Case Control Study.

BACKGROUND The present study aimed to evaluate whether the fat mass and obesity-associated (FTO) gene polymorphisms are associated with risk of intervertebral disc degeneration (IDD) in a largest Chinese Han population. MATERIAL AND METHODS There were 502 IDD patients and 497 healthy controls enrolled in this study. Nineteen single nucleotide polymorphisms (SNPs) in the FTO gene were tested using the Sequenom MassARRAY platform. The Hardy-Weinberg equilibrium test, followed by allelic, genotypic, haplotypic association, and SNP interaction analyses were used for SNP evaluation. The Genotype-Tissue Expression (GTEx) database was used to evaluate expression quantitative trait loci (eQTL) value of polymorphism. Spearman rank correlation and logistic regression analyses were used for assessing the internal relation between genotypic changes and the risk of IDD. RESULTS Seventeen SNPs survived the Hardy-Weinberg equilibrium test. Allelic analysis showed that allele T of SNP rs1121980 was a risk allele. Haplotypic and SNP interaction analyses suggested that 2 haplotypes and 5 SNP combinations were associated with the predisposition of IDD respectively. GTEx database revealed that the SNP rs1121980 might interfere with the expression of the FTO gene in the muscle-skeletal system. Through clinical statistics analysis, the different genotypes of rs1121980 can present different disease severity of IDD. CONCLUSIONS Our study suggests that rs1121980 can become a biomarker for the screening and prognosis of IDD. The 2 haplotype blocks and 5 SNP-SNP combinations that we discovered might be indicative of the onset of IDD. Therefore, our study might serve as evidence for future IDD molecular diagnosis.

RNA editome in rhesus macaque shaped by purifying selection.

Understanding of the RNA editing process has been broadened considerably by the next generation sequencing technology; however, several issues regarding this regulatory step remain unresolved--the strategies to accurately delineate the editome, the mechanism by which its profile is maintained, and its evolutionary and functional relevance. Here we report an accurate and quantitative profile of the RNA editome for rhesus macaque, a close relative of human. By combining genome and transcriptome sequencing of multiple tissues from the same animal, we identified 31,250 editing sites, of which 99.8% are A-to-G transitions. We verified 96.6% of editing sites in coding regions and 97.5% of randomly selected sites in non-coding regions, as well as the corresponding levels of editing by multiple independent means, demonstrating the feasibility of our experimental paradigm. Several lines of evidence supported the notion that the adenosine deamination is associated with the macaque editome--A-to-G editing sites were flanked by sequences with the attributes of ADAR substrates, and both the sequence context and the expression profile of ADARs are relevant factors in determining the quantitative variance of RNA editing across different sites and tissue types. In support of the functional relevance of some of these editing sites, substitution valley of decreased divergence was detected around the editing site, suggesting the evolutionary constraint in maintaining some of these editing substrates with their double-stranded structure. These findings thus complement the "continuous probing" model that postulates tinkering-based origination of a small proportion of functional editing sites. In conclusion, the macaque editome reported here highlights RNA editing as a widespread functional regulation in primate evolution, and provides an informative framework for further understanding RNA editing in human.

Selectively Constrained RNA Editing Regulation Crosstalks with piRNA Biogenesis in Primates.

Although millions of RNA editing events have been reported to modify hereditary information across the primate transcriptome, evidence for their functional significance remains largely elusive, particularly for the vast majority of editing sites in noncoding regions. Here, we report a new mechanism for the functionality of RNA editing-a crosstalk with PIWI-interacting RNA (piRNA) biogenesis. Exploiting rhesus macaque as an emerging model organism closely related to human, in combination with extensive genome and transcriptome sequencing in seven tissues of the same animal, we deciphered accurate RNA editome across both long transcripts and the piRNA species. Superimposing and comparing these two distinct RNA editome profiles revealed 4,170 editing-bearing piRNA variants, or epiRNAs, that primarily derived from edited long transcripts. These epiRNAs represent distinct entities that evidence an intersection between RNA editing regulations and piRNA biogenesis. Population genetics analyses in a macaque population of 31 independent animals further demonstrated that the epiRNA-associated RNA editing is maintained by purifying selection, lending support to the functional significance of this crosstalk in rhesus macaque. Correspondingly, these findings are consistent in human, supporting the conservation of this mechanism during the primate evolution. Overall, our study reports the earliest lines of evidence for a crosstalk between selectively constrained RNA editing regulation and piRNA biogenesis, and further illustrates that such an interaction may contribute substantially to the diversification of the piRNA repertoire in primates.

Evolutionary interrogation of human biology in well-annotated genomic framework of rhesus macaque.

With genome sequence and composition highly analogous to human, rhesus macaque represents a unique reference for evolutionary studies of human biology. Here, we developed a comprehensive genomic framework of rhesus macaque, the RhesusBase2, for evolutionary interrogation of human genes and the associated regulations. A total of 1,667 next-generation sequencing (NGS) data sets were processed, integrated, and evaluated, generating 51.2 million new functional annotation records. With extensive NGS annotations, RhesusBase2 refined the fine-scale structures in 30% of the macaque Ensembl transcripts, reporting an accurate, up-to-date set of macaque gene models. On the basis of these annotations and accurate macaque gene models, we further developed an NGS-oriented Molecular Evolution Gateway to access and visualize macaque annotations in reference to human orthologous genes and associated regulations (www.rhesusbase.org/molEvo). We highlighted the application of this well-annotated genomic framework in generating hypothetical link of human-biased regulations to human-specific traits, by using mechanistic characterization of the DIEXF gene as an example that provides novel clues to the understanding of digestive system reduction in human evolution. On a global scale, we also identified a catalog of 9,295 human-biased regulatory events, which may represent novel elements that have a substantial impact on shaping human transcriptome and possibly underpin recent human phenotypic evolution. Taken together, we provide an NGS data-driven, information-rich framework that will broadly benefit genomics research in general and serves as an important resource for in-depth evolutionary studies of human biology.

RhesusBase: a knowledgebase for the monkey research community.

Although the rhesus macaque is a unique model for the translational study of human diseases, currently its use in biomedical research is still in its infant stage due to error-prone gene structures and limited annotations. Here, we present RhesusBase for the monkey research community (http://www.rhesusbase.org). We performed strand-specific RNA-Seq studies in 10 macaque tissues and generated 1.2 billion 90-bp paired-end reads, covering >97.4% of the putative exon in macaque transcripts annotated by Ensembl. We found that at least 28.7% of the macaque transcripts were previously mis-annotated, mainly due to incorrect exon-intron boundaries, incomplete untranslated regions (UTRs) and missed exons. Compared with the previous gene models, the revised transcripts show clearer sequence motifs near splicing junctions and the end of UTRs, as well as cleaner patterns of exon-intron distribution for expression tags and cross-species conservation scores. Strikingly, 1292 exon-intron boundary revisions between coding exons corrected the previously mis-annotated open reading frames. The revised gene models were experimentally verified in randomly selected cases. We further integrated functional genomics annotations from >60 categories of public and in-house resources and developed an online accessible database. User-friendly interfaces were developed to update, retrieve, visualize and download the RhesusBase meta-data, providing a 'one-stop' resource for the monkey research community.

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.

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.

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.

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.

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.

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.