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.

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.

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.

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.

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.

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 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.

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.

Hydrogen peroxide primes heart regeneration with a derepression mechanism.

While the adult human heart has very limited regenerative potential, the adult zebrafish heart can fully regenerate after 20% ventricular resection. Although previous reports suggest that developmental signaling pathways such as FGF and PDGF are reused in adult heart regeneration, the underlying intracellular mechanisms remain largely unknown. Here we show that H2O2 acts as a novel epicardial and myocardial signal to prime the heart for regeneration in adult zebrafish. Live imaging of intact hearts revealed highly localized H2O2 (~30 µM) production in the epicardium and adjacent compact myocardium at the resection site. Decreasing H2O2 formation with the Duox inhibitors diphenyleneiodonium (DPI) or apocynin, or scavenging H2O2 by catalase overexpression markedly impaired cardiac regeneration while exogenous H2O2 rescued the inhibitory effects of DPI on cardiac regeneration, indicating that H2O2 is an essential and sufficient signal in this process. Mechanistically, elevated H2O2 destabilized the redox-sensitive phosphatase Dusp6 and hence increased the phosphorylation of Erk1/2. The Dusp6 inhibitor BCI achieved similar pro-regenerative effects while transgenic overexpression of dusp6 impaired cardiac regeneration. H2O2 plays a dual role in recruiting immune cells and promoting heart regeneration through two relatively independent pathways. We conclude that H2O2 potentially generated from Duox/Nox2 promotes heart regeneration in zebrafish by unleashing MAP kinase signaling through a derepression mechanism involving Dusp6.