Next-generation sequencing identifies novel genes with rare variants in total anomalous pulmonary venous connection

Comprehensive analysis of miRNAs-lncRNAs-mRNAs modules and ceRNA network in acute liver failure: Hsa-miR3175 and C-reactive protein determination

Acute liver failure (ALF), also known as fulminant hepatitis, coagulation disorders, and worsening mental status. It has a poor prognosis and high mortality rate. Among these, the top 10 upregulated genes included GKA-DPA1, IGLL5, PLA2G7, CCL5, IGLJ, GUSBP11, RHOBT1, IGLL3P, CCL18, and ADRBK2. On the other hand, the top 10 downregulated genes were SLC6A1, PID1, AVPR1A, PP1R1A, ST3GAL6, TPST, ERO1LB, SLCO4C1, and KLF15. Furthermore, the DEGs were found to be enriched in processes related to LIAO metastasis and creighton endocrine therapy resistance. To explore the interactions among the DEGs, we constructed a PPI network. This network revealed 16 hub genes that play crucial roles in ALF pathogenesis. Within this network, hsa-mir-375 and hsa-mir-650 were identified as central nodes, indicating their potential importance in ALF. By identifying and analyzing the transcriptional-level ceRNA network, we have provided valuable insights into the etiology of ALF.

Genome-wide profiling of angiogenic cis-regulatory elements unravels cis-regulatory SNPs for vascular abnormality

Angiogenesis is extensively involved in embryonic development and requires complex regulation networks, whose defects can cause a variety of vascular abnormalities. Cis-regulatory elements control gene expression at all developmental stages, but they have not been studied or profiled in angiogenesis yet. In this study, we exploited public DNase-seq and RNA-seq datasets from a VEGFA-stimulated in vitro angiogenic model, and carried out an integrated analysis of the transcriptome and chromatin accessibility across the entire process. Totally, we generated a bank of 47,125 angiogenic cis-regulatory elements with promoter (marker by H3K4me3) and/or enhancer (marker by H3K27ac) activities. Motif enrichment analysis revealed that these angiogenic cis-regulatory elements interacted preferentially with ETS family TFs. With this tool, we performed an association study using our WES data of TAPVC and identified rs199530718 as a cis-regulatory SNP associated with disease risk. Altogether, this study generated a genome-wide bank of angiogenic cis-regulatory elements and illustrated its utility in identifying novel cis-regulatory SNPs for TAPVC, expanding new horizons of angiogenesis as well as vascular abnormality genetics.

m6A regulators are differently expressed and correlated with immune response of pancreatic adenocarcinoma

BACKGROUND

N6 methyladenosine (m6A) RNA methylation regulators play a key role in the occurrence and development of many tumors. However, the function of N6 methyladenosine (m6A) RNA methylation regulators in pancreatic adenocarcinoma (PAAD) has not been fully clarified.

METHODS

We used data set from GEPIA 2, UALCAN, TIMER, TISIDB, CBioPortal database to analyze the gene expression of 20 major m6A RNA methylation regulators.

RESULTS

Our study revealed that the irregularity of m6A regulators were associated with poor prognosis in PAAD. Meantime, 13 m6A regulators showed high expression in PAAD samples (ALKBH5, ELAVL1, FTO, HNRNPC, IGF2BP2, METTL14, METTL16 (METT10D), RBM15, VIRMA (KIAA1429), YTHDF1, YTHDF2, YTHDF3 and ZC3H13). In these regulators, we evaluated HNRNPC and IGF2BP2 were significantly correlated with worse outcomes and ALKBH5, IGF2BP2, METTL16 (METT10D), RBM15 were significantly correlated with PAAD in advanced stage. Moreover, we showed m6A regulators is correlated with Immuno-regulators' (Immunoinhibitors, Immunostimulators and MHC molecules) expression and levels of immune infiltration in PAAD. Bioinformatics further demonstrate m6A regulators were participated in revising in RNA processing.

CONCLUSIONS

Our study investigated that the m6A regulatory factors may serve as a biomarker and a potential target of immunotherapy for PAAD.

Identification of immune-related gene signature for predicting prognosis in uterine corpus endometrial carcinoma

The objective of this study is to develop a gene signature related to the immune system that can be used to create personalized immunotherapy for Uterine Corpus Endometrial Carcinoma (UCEC). To classify the UCEC samples into different immune clusters, we utilized consensus clustering analysis. Additionally, immune correlation algorithms were employed to investigate the tumor immune microenvironment (TIME) in diverse clusters. To explore the biological function, we conducted GSEA analysis. Next, we developed a Nomogram by integrating a prognostic model with clinical features. Finally, we performed experimental validation in vitro to verify our prognostic risk model. In our study, we classified UCEC patients into three clusters using consensus clustering. We hypothesized that cluster C1 represents the immune inflammation type, cluster C2 represents the immune rejection type, and cluster C3 represents the immune desert type. The hub genes identified in the training cohort were primarily enriched in the MAPK signaling pathway, as well as the PD-L1 expression and PD-1 checkpoint pathway in cancer, all of which are immune-related pathways. Cluster C1 may be a more suitable for immunotherapy. The prognostic risk model showed a strong predictive ability. Our constructed risk model demonstrated a high level of accuracy in predicting the prognosis of UCEC, while also effectively reflecting the state of TIME.

The 100 most cited papers on total anomalous pulmonary venous connection: a bibliometric analysis

BACKGROUND

The number of citations a paper receives reflects its impact on the scientific community. We aimed to identify and explore the characteristics of the most cited papers on total anomalous pulmonary venous connection (TAPVC).

METHODS

Web of Science Core Collection Expanded Science Citation Index (1900 to present) was searched and papers on TAPVC were reviewed. Articles were ranked by the number of citations and the 100 most cited papers were analyzed.

RESULTS

The 100 most cited papers were published between 1952 and 2018 with a mean number of citations of 52 (range 26 to 148). The 1990s was the most productive decade. All articles except one were written in English. The 100 most cited articles were published in 24 journals, led by Journal of Thoracic and Cardiovascular Surgery (21 articles), followed by Annals of Thoracic Surgery (20 articles), and Circulation (16 articles). The United States of America contributed most of the 100 most cited papers (60 articles). Hospital for Sick Children, Toronto led the list of citation classics with six papers. Christopher A. Caldarone, John W. Kirklin, and P. E. F. Daubeney were the most productive authors with 3 articles each. More than half of the papers were cohort studies (51 articles). Surgery, radiology and etiology were the main topics. Thirty-one articles were funded by public foundations, and none received support from commercial companies.

CONCLUSIONS

The bibliometric analysis gives a historical perspective on scientific progress in the field of TAPVC and lays the foundation for future research.

Genetic Screening of Targeted Region on the Chromosome 22q11.2 in Patients with Microtia and Congenital Heart Defect

Microtia is a congenital malformation characterized by a small, abnormally shaped auricle (pinna) ranging in severity. Congenital heart defect (CHD) is one of the comorbid anomalies with microtia. However, the genetic basis of the co-existence of microtia and CHD remains unclear. Copy number variations (CNVs) of 22q11.2 contribute significantly to microtia and CHD, respectively, thus suggesting a possible shared genetic cause embedded in this genomic region. In this study, 19 sporadic patients with microtia and CHD, as well as a nuclear family, were enrolled for genetic screening of single nucleotide variations (SNVs) and CNVs in 22q11.2 by target capture sequencing. We detected a total of 105 potential deleterious variations, which were enriched in ear- or heart-development-related genes, including TBX1 and DGCR8. The gene burden analysis also suggested that these genes carry more deleterious mutations in the patients, as well as several other genes associated with cardiac development, such as CLTCL1. Additionally, a microduplication harboring SUSD2 was validated in an independent cohort. This study provides new insights into the underlying mechanisms for the comorbidity of microtia and CHD focusing on chromosome 22q11.2, and suggests that a combination of genetic variations, including SNVs and CNVs, may play a crucial role instead of single gene mutation.

Integrated genomic analysis identifies novel low-frequency cis-regulatory variant rs2279658 associated with VSD risk in Chinese children

VSD combined with other cardiac or extracardiac malformations (defined as "complex VSD" by us) is one of the major causes of perinatal morbidity and mortality. Functional non-coding SNPs (cis-regulatory SNPs) have not been systematically studied in CHDs, including complex VSD. Here we report an exome-wide association analysis using WES data of 60 PA/VSD cases, 20 TOF cases and 100 controls in Chinese children. We identify 93 low-frequency non-coding SNPs associated with complex VSD risk. A functional genomics pipeline integrating ATAC-seq, ChIP-seq and promoter CHi-C recognizes the rs2279658 variant as a candidate cis-regulatory SNP. Specifically, rs2279658 resides in a cardiac-specific enhancer bound by FOXH1 and PITX2, and would abrogate binding of these two transcription factors to the identified enhancer during cardiac morphogenesis. COQ2 and FAM175A are predicted to be target genes for "rs2279658-FOXH1 or PITX2" pairs in the heart. These findings highlight the importance of cis-regulatory SNPs in the pathogenesis of complex VSD and broaden our understanding of this disease.

Identification of energy metabolism-related biomarkers for risk prediction of heart failure patients using random forest algorithm

Objective

Energy metabolism plays a crucial role in the improvement of heart dysfunction as well as the development of heart failure (HF). The current study is designed to identify energy metabolism-related diagnostic biomarkers for predicting the risk of HF due to myocardial infarction.

Methods

Transcriptome sequencing data of HF patients and non-heart failure (NF) people (GSE66360 and GSE59867) were obtained from gene expression omnibus (GEO) database. Energy metabolism-related differentially expressed genes (DEGs) were screened between HF and NF samples. The subtyping consistency analysis was performed to enable the samples to be grouped. The immune infiltration level among subtypes was assessed by single sample gene set enrichment analysis (ssGSEA). Random forest algorithm (RF) and support vector machine (SVM) were applied to identify diagnostic biomarkers, and the receiver operating characteristic curves (ROC) was plotted to validate the accuracy. Predictive nomogram was constructed and validated based on the result of the RF. Drug screening and gene-miRNA network were analyzed to predict the energy metabolism-related drugs and potential molecular mechanism.

Results

A total of 22 energy metabolism-related DEGs were identified between HF and NF patients. The clustering analysis showed that HF patients could be classified into two subtypes based on the energy metabolism-related genes, and functional analyses demonstrated that the identified DEGs among two clusters were mainly involved in immune response regulating signaling pathway and lipid and atherosclerosis. ssGSEA analysis revealed that there were significant differences in the infiltration levels of immune cells between two subtypes of HF patients. Random-forest and support vector machine algorithm eventually identified ten diagnostic markers (MEF2D, RXRA, PPARA, FOXO1, PPARD, PPP3CB, MAPK14, CREB1, MEF2A, PRMT1) for risk prediction of HF patients, and the proposed nomogram resulted in good predictive performance (GSE66360, AUC = 0.91; GSE59867, AUC = 0.84) and the clinical usefulness in HF patients. More importantly, 10 drugs and 15 miRNA were predicted as drug target and hub miRNA that associated with energy metabolism-related genes, providing further information on clinical HF treatment.

Conclusion

This study identified ten energy metabolism-related diagnostic markers using random forest algorithm, which may help optimize risk stratification and clinical treatment in HF patients.

HSPA5 Could Be a Prognostic Biomarker Correlated with Immune Infiltration in Breast Cancer

Background

Breast cancer (BC) is a frequent disease in females. The heat shock 70 kDa protein 5 (HSPA5) has recently been discovered to have an important function in tumor growth. However, the biological significance of HSPA5 in BC is unknown. Material and Method. Firstly, The Cancer Genome Atlas (TCGA) database was applied to analyze the expressions of HSPA5 in different cancer types, especially in BC. Then, the LinkedOmics database was used to screen genes coexpressed with HSPA5 in BC, presented by protein-protein interaction (PPI) and analyzed by functional enrichment analyses. Next, the Kaplan-Meier plotter was adopted to study the prognostic significance of HSPA5 and the relation between HSPA5 expression and different clinical factors in BC. Finally, the Tumor Immune Estimation Resource (TIMER) method was adopted to explore the relation between immune infiltration and HSPA5 in BC.

Result

HSPA5 was highly expressed in most cancers, including BC. Genes coexpressed with HSPA5 were mainly related to endoplasmic reticulum unfolded protein response, melanosome, thyroid hormone synthesis, N-glycan biosynthesis, and so on. In the survival analysis, high HSPA5 expression indicated a poor prognosis in BC, and the expression of HSPA5 in BC was elevated after the incidence of BC, changing with different clinical factors. In the immune infiltration, HSPA5 was positively correlated with most immune cells.

Conclusion

HSPA5 is an oncogene in BC progression, and it is connected with the prognosis and the immune infiltration in BC. Our findings suggest that HSPA5 could be an immunotherapy target and a prognostic biomarker in BC.

Peripheral Blood Transcripts Predict Preoperative Obstructive Total Anomalous Pulmonary Venous Connection

The lack of accessible noninvasive tools to examine the molecular alterations limits our understanding of the causes of total anomalous pulmonary venous connection (TAPVC), as well as the identification of effective operational strategies. Here, we consecutively enrolled peripheral leukocyte transcripts of 26 preoperative obstructive and 22 non-obstructive patients with TAPVC. Two-hundred and fifty six differentially expressed mRNA and 27 differentially expressed long noncoding RNA transcripts were dysregulated. The up-regulated mRNA was enriched in the hydrogen peroxide catabolic process, response to mechanical stimulus, neutrophil degranulation, hemostasis, response to bacterium, and the NABA CORE MATRISOME pathway, all of which are associated with the development of fibrosis. Furthermore, we constructed predictive models using multiple machine-learning algorithms and tested the performance in the validation set. The mRNA NR3C2 and lncRNA MEG3 were screened based on multiple iterations. The random forest prediction model can predict preoperative obstruction patients in the validation set with high accuracy (area under curve = 1; sensitivity = 1). These data highlight the potential of peripheral leukocyte transcripts to evaluate obstructive-related pathophysiological alterations, leading to precision healthcare solutions that could improve patient survival after surgery. It also provides a novel direction for the study of preoperative obstructive TAPVC.

Expression of ferroptosis-related gene correlates with immune microenvironment and predicts prognosis in gastric cancer

The study is to explore the role of ferroptosis-related genes (FRGs) in the occurrence and development of gastric cancer (GC), and to construct a new prognosis signature to predict the prognosis in GC. Clinical information and corresponding RNA data of GC patients were downloaded from TCGA and GEO databases. Consensus clustering was performed to identify new molecular subgroups. ESTIMATE, CIBERSORT, McpCounter and TIMER algorithm were used to analyze the infiltration of immune cells in two molecular subgroups. LASSO algorithm and multivariate Cox analysis were used to construct a prognostic risk signature. Functional analysis was conducted to elucidate the underlying mechanisms. Finally, the FRPGs were verified by Quantitative Real-Time PCR. We obtained 16 FRGs and divided GC patients into two subgroups by consistent clustering. Cluster C1 with a higher abundance of immune cell infiltration but lower probability in response to immunotherapy, it was reasonable to speculate that Cluster C1 was in accordance with the immune rejection type. Functional analysis showed that the biological process of DEGs in training cohort mainly included immune globulin, and human immune response mediated by circulating immune globulin. GSEA analysis showed that compared with Cluster C2, Cluster C1 showed lower expression in lipid metabolism. The nomogram combined with risk signature and clinical features can accurately predict the prognosis of GC patients. We identified two molecular subtypes, Clusters C1 and C2. In Cluster C1, patients with poor prognosis present with a hyperimmune status and low lipid metabolism, and we speculate that Cluster C1 was in accordance with the immune rejection type. The risk model based on FRPGs can accurately predict the prognosis of GC. These results indicated that ferroptosis is associated with TIME, and deserved considerable attention in determining immunotherapy treatment strategy for GC patients.

Identification of novel potential biomarkers in infantile hemangioma via weighted gene co-expression network analysis

Background

Infantile hemangioma (IH) is the most common benign tumor in children and is characterized by endothelial cells proliferation and angiogenesis. Some hub genes may play a critical role in angiogenesis. This study aimed to identify the hub genes and analyze their biological functions in IH.

Methods

Differentially expressed genes (DEGs) in hemangioma tissues, regardless of different stages, were identified by microarray analysis. The hub genes were selected through integrated weighted gene co-expression network analysis (WGCNA) and protein–protein interaction (PPI) network. Subsequently, detailed bioinformatics analysis of the hub genes was performed by gene set enrichment analysis (GSEA). Finally, quantitative real-time polymerase chain reaction (qRT-PCR) analysis was conducted to validate the hub genes expression in hemangioma-derived endothelial cells (HemECs) and human umbilical vein endothelial cells (HUVECs).

Results

In total, 1115 DEGs were identified between the hemangiomas and normal samples, including 754 upregulated genes and 361 downregulated genes. Two co-expression modules were identified by WGCNA and green module eigengenes were highly correlated with hemangioma (correlation coefficient = 0.87). Using module membership (MM) > 0.8 and gene significance (GS) > 0.8 as the cut-off criteria, 108 candidate genes were selected and put into the PPI network, and three most correlated genes (APLN, APLNR, TMEM132A) were identified as the hub genes. GSEA predicted that the hub genes would regulate endothelial cell proliferation and angiogenesis. The differential expression of these genes was validated by qRT-PCR.

Conclusions

This research suggested that the identified hub genes may be associated with the angiogenesis of IH. These genes may improve our understanding of the mechanism of IH and represent potential anti-angiogenesis therapeutic targets for IH.

Association of PLXND1 with a novel subtype of anomalous pulmonary venous return

Anomalous pulmonary venous return (APVR) is a potentially lethal congenital heart disease. Elucidating the genetic etiology is crucial for understanding its pathogenesis and improving clinical practice, while its genetic basis remains largely unknown due to complex genetic etiology. We thus performed whole-exome sequencing for 144 APVR patients and 1636 healthy controls and report a comprehensive atlas of APVR-related rare genetic variants. Novel singleton, loss-of-function and deleterious missense variants (DVars) were enriched in patients, particularly for genes highly-expressed in the developing human heart at the critical time point for pulmonary veins draining into the left atrium. Notably, PLXND1, encoding a receptor for semaphorins, represents a strong candidate gene of APVR (adjusted P = 1.1e-03, OR: 10.9-69.3), accounting for 4.17% of APVR. We further validated this finding in an independent cohort consisting of 82 case-control pairs. In these two cohorts, eight DVars were identified in different patients, which convergently disrupt the GTPase-activating protein-related domain of PLXND1. All variant carriers displayed strikingly similar clinical features, in that all anomalous drainage of pulmonary vein(s) occurred on the right side and incorrectly connected to the right atrium, may representing a novel subtype of APVR for molecular diagnosis. Studies in Plxnd1 knockout mice further revealed the effects of PLXND1 deficiency on severe heart and lung defects and cellular abnormalities related to APVR such as abnormal migration and vascular formation of vascular endothelial cells. These findings indicate the important role of PLXND1 in APVR pathogenesis, providing novel insights into the genetic etiology and molecular subtyping for APVR.

HMGA2 contributes to vascular development and sprouting angiogenesis by promoting IGFBP2 production

Angiogenesis is the process by which new blood vessels form from preexisting vessels and regulates the processes of embryonic development, wound healing and tumorigenesis. HMGA2 is involved in the occurrence of several cancers, but its biological role and the exact downstream genes involved in vascular development and sprouting angiogenesis remain largely unknown. Here, we first found that HMGA2 knockdown in zebrafish embryos resulted in defects of central artery formation. RNA sequencing revealed that IGFBP2 was significantly downregulated by interference with HMGA2, and IGFBP2 overexpression reversed the inhibition of brain vascular development caused by HMGA2 deficiency. In vitro, we further found that HMGA2 knockdown blocked the migration, tube formation and branching of HUVECs. Similarly, IGFBP2 protein overexpression attenuated the impairments induced by HMGA2 deficiency. Moreover, the promotion of angiogenesis by HMGA2 overexpression was verified in a Matrigel plug assay. We next found that HMGA2 bound directly to a region in the IGFBP2 promoter and positively regulated IGFBP2 expression. Interestingly, the mRNA expression levels of HMGA2 and IGFBP2 were increased significantly in the peripheral blood of hemangioma patients, indicating that overexpression of HMGA2 and IGFBP2 results in vessel formation, consistent with the results of the in vivo and in vitro experiments. In summary, our findings demonstrate that HMGA2 promotes central artery formation by modulating angiogenesis via IGFBP2 induction.

Identification of Hub mRNAs and lncRNAs in Atrial Fibrillation Using Weighted Co-expression Network Analysis With RNA-Seq Data

Atrial fibrillation (AF)/paroxysmal AF (PAF) is the main cause of cardiogenic embolism. In recent years, the progression from paroxysmal AF to persistent AF has attracted more and more attention. However, the molecular mechanism of the progression of AF is unclear. In this study, we performed RNA sequencing for normal samples, paroxysmal AF and persistent AF samples to identify differentially expressed gene (DEG) and explore the roles of these DEGs in AF. Totally, 272 differently expressed mRNAs (DEmRNAs) and 286 differentially expressed lncRNAs (DElncRNAs) were identified in paroxysmal AF compared to normal samples; 324 DEmRNAs and 258 DElncRNAs were found in persistent atrial fibrillation compared with normal samples; and 520 DEmRNAs and 414 DElncRNAs were identified in persistent AF compared to paroxysmal AF samples. Interestingly, among the DEGs, approximately 50% were coding genes and around 50% were non-coding RNAs, suggesting that lncRNAs may also have a crucial role in the progression of AF. Bioinformatics analysis demonstrated that these DEGs were significantly related to regulating multiple AF associated pathways, such as the regulation of vascular endothelial growth factor production and binding to the CXCR chemokine receptor. Furthermore, weighted gene co-expression network analysis (WGCNA) was conducted to identify key modules and hub RNAs and lncRNAs to determine their potential associations with AF. Five hub modules were identified in the progression of AF, including blue, brown, gray, turquoise and yellow modules. Interestingly, blue module and turquoise module were significantly negatively and positively correlated to the progression of AF respectively, indicating that they may have a more important role in the AF. Moreover, the hub protein-protein interaction (PPI) networks and lncRNA-mRNA regulatory network were constructed. Bioinformatics analysis on the hub PPI network in turquoise was involved in regulating immune response related signaling, such as leukocyte chemotaxis, macrophage activation, and positive regulation of α-β T cell activation. Our findings could clarify the underlying molecular changes associated fibrillation, and provide a useful resource for identifying AF marker.

The Functionalities and Clinical Significance of Tumor-Infiltrating Immune Cells in Esophageal Squamous Cell Carcinoma

Tumor-infiltrating immune cells have been implicated in the tumorigenesis and progression of esophageal squamous cell carcinoma (ESCC). However, the functionalities and clinical significance of immune cells remain largely unveiled. In this study, the gene expression data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were extracted. The relative infiltrating levels were estimated by single-sample gene set enrichment analysis. Some cytotoxic immune cells were attenuated, and resting cytotoxic immune cells were accumulated in ESCC. Remarkably, we also observed that infiltrating levels of macrophage M2 and resting natural killer (NK) cells were increased in nonresponders of CRT, and T cells that had anticancer activities such as activated memory CD4 and T helper 2 (Th2) cells were significantly reduced in ESCC tissues of the nonresponders. Moreover, the high infiltrations of the resting natural killer (NK) and dendritic cell (DC) were observed to result in a shorter overall survival in ESCC. Consistently, high expression of immune checkpoint genes, CTLA4 and HAVCR2, was associated with poor prognosis. Furthermore, STAT5B, a key transcription factor, as well as its target genes, involved in the regulation of T cells, was significantly downregulated in ESCC, especially subgroup I, indicating that downregulation of STAT5B might be associated with reduced T cell-mediated anticancer activity. In conclusion, the present study significantly improved our understanding of the regulatory roles of immune cells in ESCC.

Curcumin Improves Pulmonary Hypertension Rats by Regulating Mitochondrial Function

Objective

To investigate the role of curcumin in regulating pathogenesis of pulmonary arterial smooth muscle cells (PASMCs) derived from pulmonary arterial hypertension (PAH) model.

Methods

Male Sprague Dawley rats were injected with monocrotaline (MCT) to establish the PAH experimental model. The rats were divided into control group, MCT group, and curcumin group. At the end of the study, hemodynamic data were measured to determine pulmonary hypertension. Proliferation ability of PASMCs, a remodeling indicator of pulmonary artery and right ventricle, was detected. In addition, the morphology and function of mitochondria, antiglycolysis and antiproliferation pathways, and genes were also analyzed.

Results

Curcumin may function by reversing MCT-mediated pulmonary vascular remodeling in rats. Curcumin effectively improved pulmonary vascular remodeling, promoted PASMC apoptosis, and protected mitochondrial function. In addition, curcumin treatment suppressed the PI3K/AKT pathway in PASMCs and regulated the expression of antiproliferative genes.

Conclusion

Curcumin can improve energy metabolism and reverse the process of PAHS. However, there were side effects of curcumin in MCT-induced rats, suggesting that the dosage should be treated with caution and its toxicological mechanism should be further studied and evaluated.

Other Types of Glycosylation

O-Linked glycosylation such as O-fucose, O-glucose, and O-N-acetylglucosamine are considered to be unusual. As suggested by the high levels of evolutional conservation, these O-glycans are fundamentally important for life. In the last two decades, our understanding of the importance of these glycans has greatly advanced. In particular, identification of the glycosyltransferases responsible for the biosynthesis of these glycans has accelerated basic research on the functional significance and molecular mechanisms by which these O-glycans regulate protein functions as well as clinical research on human diseases due to changes in these types of O-glycosylation. Notably, Notch receptor signaling is modified with and regulated by these types of O-glycans. Here, we summarize the current view of the structures and the significance of these O-glycans mainly in the context of Notch signaling regulation and human diseases.

Upregulation of miR-29c-3p Hinders Melanoma Progression by Inhibiting CDCA4 Expression

OBJECTIVE

To investigate the expression and regulation mechanism of miR-29c-3p and cell division cycle associated 4 (CDCA4) in melanoma (MM). Data and Methods. Fifty-nine patients with MM admitted to our hospital were enrolled as the MM group. They were followed up for 3 years to analyze the prognostic factors; meanwhile, 51 healthy subjects were allocated into a normal group. MM cell lines (M21 and C8161) were transfected with miR-29c-3p-mimics, miR-29c-3p-inhibitor, miR-NC, si-CDCA4, and sh-CDCA4. The expression of miR-29c-3p, CDCA4, Bax, Caspase3, Bcl-2, N-cadherin, vimentin, and E-cadherin was quantified, and cell proliferation, migration, invasion, and apoptosis, as well as epithelial-mesenchymal transition (EMT), were determined.

RESULTS

Serum miR-29c-3p was lowly expressed and CDCA4 was highly expressed in the MM group. The area under the curve (AUC) of both for diagnosing MM was greater than 0.9. miR-29c-3p and CDCA4 were related to regional lymph node staging (N staging), distant metastasis (M staging), tumor diameter, and pathological differentiation. Low miR-29c-3p and high CDCA4 were associated with poor prognosis of MM. Overexpression of miR-29c-3p and suppression of CDCA4 hindered cell proliferation, migration, invasion, and expression of Bax, Caspase3, N-cadherin, and vimentin, but cell apoptosis and expression of Bcl-2 and E-cadherin were enhanced. Dual-luciferase reporter (DLR) assay confirmed the targeted relationship between miR-29c-3p and CDCA4. After miR-29c-3p-mimics+sh-CDCA4 was transfected into M21 and C8161 cells, the proliferation, invasion, and apoptosis were not different from those in the miR-NC group transfected with unrelated sequences.

CONCLUSION

Overexpression of miR-29c-3p suppresses CDCA4 expression and decreases proliferation, migration, invasion, apoptosis, and EMT of MM cells, thus hindering MM progression.

CircITGA7 Suppresses Gastric Cancer Progression Through miR-1471/MTDH Axis

In recent years, there have been reports about the involvement of circular RNAs (circRNAs) in the pathogenesis of gastric cancer (GC), but the molecular mechanism in cell proliferation, invasion, and migration is still unclear. Based on The Cancer Genome Atlas (TCGA) database, we analyzed differentially expressed circRNAs between GC and non-tumor tissues. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were used to clarify the functional role in GC. Here, we showed that circITGA7 was lowly expressed in GC tissues based on the TCGA database. In vitro, silencing the expression of circITGA7 increased cell proliferation and metastasis, whereas overexpression did the opposite. Mechanistically, miR-1471 has circITGA7 as a sponge, and miR-1471 has metadherin (MTDH) as a target gene. Consequently, functional analysis showed that the tumor suppressor effect of circITGA7 was the result of regulating the miR-1471/MTDH axis. Overall, the circITGA7/miR-1471/MTDH signaling pathway may play a crucial role in GC, providing a new potential mechanism involved in GC progression.

An in-depth analysis reveals two new genetic variants on 22q11.2 associated with vitiligo in the Chinese Han population

BACKGROUND

Vitiligo is a complex disease in which patchy depigmentation is the result of an autoimmune-induced loss of melanocytes in affected regions. On the basis of a genome-wide linkage analysis of vitiligo in the Chinese Han population, we previously showed significant evidence of a linkage between 22q12 and vitiligo. Our aim in the current study was to identify vitiligo susceptibility variants within an expanded region of the 22q12 locus.

METHODS AND RESULTS

An in-depth analysis of the expanded region of the 22q12 locus was performed by imputation using a large GWAS dataset consisting of 1117 cases and 1701 controls. Eight nominal SNPs were selected and genotyped in an independent cohort of Chinese Han individuals (2069 patients and 1370 control individuals) by using the Sequenom MassArray iPLEX1 system. The data were analyzed with PLINK 1.07 software. The C allele of rs730669 located in ZDHHC8/RTN4R showed a strong association with vitiligo (P = 3.25 × 10^-8, OR = 0.81). The C allele of rs4820338 located in VPREB1 and the A allele of rs2051582 (a SNP reported in our previous study) located in IL2RB showed a suggestive association with vitiligo (P = 1.04 × 10^-5, OR = 0.86; P = 1.78 × 10^-6, OR = 1.27). The three identified SNPs showed independent associations with vitiligo in a conditional logistic regression analysis (all P < 1.0 × 10^-5; all D' < 0.05 and r² < 1.0 × 10^-4).

CONCLUSIONS

The study reveals that two novel variants rs730669 (ZDHHC8/RTN4R) and rs4820338 (VPREB1) on 22q11.2 might confer susceptibility to vitiligo and affect disease subphenotypes. The presence of multiple independent variants emphasizes their important roles in the genetic pathogenesis of disease.

Integrative Analysis of Differently Expressed Genes Reveals a 17-Gene Prognosis Signature for Endometrial Carcinoma

Endometrial carcinoma (EC) is the fifth widely occurring malignant neoplasm among women all over the world. However, there is still lacking efficacy indicators for EC's prognosis. Here, we analyzed two databases including an RNA-sequencing-based TCGA dataset and a microarray-based GSE106191. After normalizing the raw data, we identified 114 common genes with upregulation and 308 common genes with downregulation in both the TCGA and GSE106191 databases. Bioinformatics analysis showed that the differently expressed genes in EC were related to the IL17 signaling pathway, PI3K-Akt signaling pathway, and cGMP-PKG signaling pathway. Furthermore, we performed the least absolute shrinkage and selection operator (LASSO) Cox regression analysis and generated a signature featuring 17 prognosis-related genes (MAL2, ANKRD22, METTL7B, IL32, ERFE, OAS1, TRPC1, SRPX, RAPGEF4, PSD3, SIMC1, TRPC6, WFS1, PGR, PAMR1, KCNK6, and FAM189A2) and found that it could predict OS in EC patients. The further analysis showed that OAS1, MAL2, ANKRD22, METTL7B, and IL32 were significantly upregulated in EC samples after comparison with normal samples. However, TRPC1, SRPX, RAPGEF4, PSD3, SIMC1, TRPC6, WFS1, PGR, PAMR1, KCNK6, and FAM189A2 were significantly downregulated in EC samples in comparison with normal samples. And correlation analysis showed that our results showed that the expressions of 17 prognosis-related hub genes were significantly correlated based on Pearson correlation. We here offer a newly genetic biomarker for the prediction of EC patients' prognosis.

Circular RNA circCCDC66 promotes glioma proliferation by acting as a ceRNA for miR-320a to regulate FOXM1 expression

BACKGROUND

In this study, we determine the potential roles and uncover the regulatory mechanisms of circCCDC66 in regulating cell growth and cell metastasis of glioma.

METHODS

qRT-PCR was used to detect the expressions of circCCDC66 in gliomas and tissues. The biological function of circCCDC66 in glioma cell lines was elucidated by functional experiments. Cell counting kit-8 and transwell were used to detect the effect of circCCDC66 on the proliferation, migration and invasion of glioma cells. Bioinformatics analysis was applied to reveal the targets of circCCDC66.

RESULTS

The results showed circCCDC66 was overexpressed in glioma and acted as an oncogene. CircCCDC66 knockdown suppressed the proliferation, migration, and invasion of glioma cells. We constructed a circCCDC66 regulating miRNA network and revealed miR-320a was a potential target of circCCDC66, which was down-regulated in high-grade gliomas compared to low-grade gliomas. Bioinformatics analysis showed circCCDC66-miR-320a/b axis was involved in regulating multiple cancer-related pathways. Furthermore, we identified FOXM1 as a key target of circCCDC66, which was involved in regulating DNA damage response pathways. In mechanism study, circCCDC66 could sponge miR-320a, thereby increasing the expression of FOXM1.

CONCLUSIONS

CircCCDC66 could facilitate glioma cells proliferation, invasion and migration by down-regulating miR-320a and up-regulating FOXM1.

Total Anomalous Pulmonary Venous Connection in Mother and Son with a Central 22q11.2 Microdeletion

In this clinical report, we describe a male infant and his mother, who had similar congenital heart defects. They were both diagnosed neonatally with total anomalous pulmonary venous connection (TAPVC) in combination with other heart defects. Neither of the two had any other organ malformations or dysmorphic facial features. SNP-array identified a central 22q11.2 microdeletion in the male infant and his mother as well as in the maternal grandmother and maternal aunt. The mother and the maternal aunt additionally harbored a 15q11.2 BP1-BP2 microdeletion. The maternal grandmother was unaffected by heart disease. However, heart computed tomography scan of the maternal aunt revealed a quadricuspid aortic valve. Additionally, the maternal grandmother and the maternal aunt both had significant learning disabilities. Rarely, TAPVC has been described in patients with the common 22q11.2 microdeletions. However, to the best of our knowledge, TAPVC has not previously been reported in patients with this small central 22q11.2 microdeletion. Haploinsufficiency of TBX1 was originally thought to be the main cause of the 22q11.2 microdeletion syndrome phenotype, but TBX1 is not included in the atypical central 22q11.2 microdeletion. Previous reports have suggested an association between TAPVC and the 15q11.2 BP1-BP2 microdeletion. Our report does not support this association as the maternal aunt, who harbors both microdeletions, is unaffected by TAPVC, and the male infant affected by TAPVC does not harbor the 15q11.2 BP1-BP2 microdeletion. Our findings support that genes located in the central 22q11.2 region are important for heart development and that haploinsufficiency of these genes plays a crucial role in the development of the rare heart defect TAPVC.

Long Non-coding RNA FIRRE Acts as a miR-520a-3p Sponge to Promote Gallbladder Cancer Progression via Mediating YOD1 Expression

OBJECTIVES

The role of lncRNAs in gallbladder cancer (GBC) remains poorly understood. In this study, we explored the function of functional intergenic repeating RNA element (FIRRE) in GBC.

MATERIALS AND METHODS

Whole transcriptome resequencing was performed in three pairs of GBC tissues and adjacent non-tumor tissues. lncRNA FIRRE expression was verified by real-time PCR. The function of FIRRE in GBC was evaluated by experiments in vitro and in vivo. The mechanism of FIRRE was investigated via fluorescent in situ hybridization, RNA pull-down, dual luciferase reporter assays, and RNA immunoprecipitation.

RESULTS

FIRRE level was dramatically increased in GBC tissues compared to that in the adjacent non-tumor tissues. High expression of FIRRE was closely related to clinical stage and poor prognosis in GBC patients. Moreover, FIRRE remarkably enhanced proliferation and migration, and inhibited apoptosis of GBC cells. Mechanistically, FIRRE modulated YOD1 expression by sponging miR-520a-3p, thus contributing to the development of GBC.

CONCLUSION

Our data revealed that FIRRE might act as a novel mediator in GBC progression by sponging miR-520a-3p and regulating YOD1. FIRRE might be regarded as a potential diagnostic marker or target for GBC treatment.

Pulmonary Vein Stenosis: Incremental Knowledge Gains to Improve Outcomes

Pulmonary vein stenosis remains a considerable clinical challenge, with high mortality still present in children with progressive disease. In this review, we discuss the clinical spectrum of pulmonary vein stenosis and what is known about the etiology and potential modifying and contributing factors in progressive pulmonary vein stenosis.

Upregulation Mitochondrial Carrier 1 (MTCH1) Is Associated with Cell Proliferation, Invasion, and Migration of Liver Hepatocellular Carcinoma

Among the primary causes of cancer-associated death in the world, liver hepatocellular carcinoma (LIHC) ranks the third. LIHC is defined as the sixth most frequently diagnosed carcinoma. The gene mitochondrial carrier 1 (MTCH1) is a protein-coding gene. Recent research suggests that MTCH1 may be associated with some diseases. Here, our study attempts to explore the role and implication of MTCH1 in LIHC. Kaplan Meier Plotter and GEPIA (Gene Expression Profiling Interactive Analysis) databases were employed to determine the expression of MTCH1 and its correlation with prognostic status in LIHC patients. For the first time, our results suggested that MTCH1 was aberrantly expressed in human pan-cancer and highly expressed in LIHC. Its high expression was closely associated with metastasis of tumor, stage of cancer, and poor survival of patients. Then, through enrichment analysis, MTCH1 was found to be closely related to RNA splicing in LIHC. Subsequently, we conducted a series of functional experiments. PCR data showed that LIHC cell lines and samples are highly expressed MTCH1. CCK-8 (Cell Counting Kit-8) assays and Transwell assays indicated that silencing MTCH1 certainly suppressed cell proliferation, migration, and invasion. These findings shed the clue that MTCH1 could be regarded as the potential prognosis biomarker of LIHC and a promising therapeutic target for LIHC.

Marker Genes Change of Synovial Fibroblasts in Rheumatoid Arthritis Patients

Background

Rheumatoid arthritis (RA) is a chronic condition that manifests as inflammation of synovial joints, leading to joint destruction and deformity.

Methods

We identified single-cell RNA-seq data of synovial fibroblasts from RA and osteoarthritis (OA) patients in GSE109449 dataset. RA- and OA-specific cellular subpopulations were identified, and enrichment analysis was performed. Further, key genes for RA and OA were obtained by combined analysis with differentially expressed genes (DEGs) between RA and OA in GSE56409 dataset. The diagnostic role of key genes for RA was predicted using receiver operating characteristic (ROC) curve. Finally, we identified differences in immune cell infiltration between RA and OA patients, and utilized flow cytometry, qRT-PCR, and Western blot were used to examine the immune cell and key genes in RA patients.

Results

The cluster 0 matched OA and cluster 3 matched RA and significantly enriched for neutrophil-mediated immunity and ECM receptor interaction, respectively. We identified 478 DEGs. In the top 20 degrees of connection in the PPI network, the key genes for RA were obtained by comparing with the gene markers of cluster 0 and cluster 3, respectively. ROC curve showed that CCL2 and MMP13 might be diagnostic markers for RA. We found aberrant levels of CD8+T, neutrophil, and B cells in RA fibroblasts, which were validated in clinical samples. Importantly, we also validated the differential expression of key genes between RA and OA.

Conclusion

High expression of CCL2 and MMP13 in RA may be a diagnostic and therapeutic target.

Comprehensive Analysis and Identification of Key Driver Genes for Distinguishing Between Esophageal Adenocarcinoma and Squamous Cell Carcinoma

Background: Esophageal cancer (EC) is one of the deadliest cancers in the world. However, the mechanism that drives the evolution of EC is still unclear. On this basis, we identified the key genes and molecular pathways that may be related to the progression of esophageal adenocarcinoma and squamous cell carcinoma to find potential markers or therapeutic targets. Methods: GSE26886 were obtained from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) among normal samples, EA, and squamous cell carcinoma were determined using R software. Then, potential functions of DEGs were determined using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The STRING software was used to identify the most important modules in the protein-protein interaction (PPI) network. The expression levels of hub genes were confirmed using UALCAN database. Kaplan-Meier plotters were used to confirm the correlation between hub genes and outcomes in EC. Results: In this study, we identified 1,098 genes induced in esophageal adenocarcinoma (EA) and esophageal squamous cell carcinoma (ESCC), and 669 genes were reduced in EA and ESCC, suggesting that these genes may play an important role in the occurrence and development of EC tumors. Bioinformatics analysis showed that these genes were involved in cell cycle regulation and p53 and phosphoinositide 3-kinase (PI3K)/Akt signaling pathway. In addition, we identified 147 induced genes and 130 reduced genes differentially expressed in EA and ESCC. The expression of ESCC in the EA group was different from that in the control group. By PPI network analysis, we identified 10 hub genes, including GNAQ, RGS5, MAPK1, ATP1B1, HADHA, HSDL2, SLC25A20, ACOX1, SCP2, and NLN. TCGA validation showed that these genes were present in the dysfunctional samples between EC and normal samples and between EA and ESCC. Kaplan-Meier analysis showed that MAPK1, ACOX1, SCP2, and NLN were associated with overall survival in patients with ESCC and EA. Conclusions: In this study, we identified a series of DEGs between EC and normal samples and between EA and ESCC samples. We also identified 10 key genes involved in the EC process. We believe that this study may provide a new biomarker for the prognosis of EA and ESCC.

Utility of a modified vascular corrosion casting technique in the diagnosis of fetal total anomalous pulmonary venous connection

Total anomalous pulmonary venous connection (TAPVC) is a rare congenital cardiac malformation, and prenatal detection of TAPVC malformation remains a challenging. TAPVC can be easily missed or misdiagnosed in prenatal examinations. This study was aimed to use the modified vascular corrosion casting technique to prepare fetal cardiovascular casts with TAPVC and investigate the utility of cardiovascular casting for the demonstration of fetal TAPVC. The retrospective study enrolled twenty fetuses (22 to 29 + 4 gestational weeks) with TAPVC diagnosed by prenatal echocardiography and casting technique from May 2015 to May 2020. Pre- and postnatal medical records, including results obtained by prenatal ultrasound, postpartum computed tomography angiography, as well as anatomic and cardiovascular casting findings were carefully reviewed and analyzed. In twenty cases, 80% (16/20) had intra- or extracardiac malformations. The TAPVC types were supracardiac (n = 8), cardiac (n = 6), infracardiac (n = 4), and mixed (n = 2). The diagnosis of 1 case each of supracardiac and cardiac TAPVC was modified to partial anomalous pulmonary venous connection; additionally, 4 malformations were missed and 2 were misdiagnosed, including an anomalous left brachiocephalic vein in supracardiac TAPVC, abnormal inflow of the hepatic vein and a double inferior vena cava in infracardiac TAPVC; and bilateral ductus arteriosus in infracardiac TAPVC; a tetralogy of Fallot in cardiac TAPVC that was corrected to right ventricular double outlet; and an absence of ductus arteriosus that was misdiagnosed as slim ductus arteriosus. Comparing with ultrasound, casting technique has its own superiority in exhibiting TAPVC abnormalities, especially in certain types such as course, origin and absence abnormalities of ductus. Postpartum cardiovascular casts can accurately depict the branch structure of the heart's larger vessels, and may be used as a clinical assessment and teaching method in complex cardiac malformations.

Construction of a Nine-MicroRNA-Based Signature to Predict the Overall Survival of Esophageal Cancer Patients

BACKGROUND

Esophageal cancer (EC) is a common malignant tumor. MicroRNAs (miRNAs) play a key role in the occurrence and metastasis and are closely related to the prognosis of EC. Therefore, it will provide a powerful tool to predict the overall survival (OS) of EC patients based on miRNAs expression in EC tissues and blood samples.

METHODS

Five independent databases, TCGA, GSE106817, GSE113486, GSE122497, and GSE112264, were used to construct nine-miRna signature and nomograms for prognosis. The bioinformatics analysis was used to predict the enrichment pathways of targets.

RESULTS

A total of 132 overexpressed miRNAs and 23 suppressed miRNAs showed significant differential expression in both EC serum and tissue samples compared with normal samples. We also showed that nine miRNAs were related to the prognosis of EC. Higher levels of miR-15a-5p, miR-92a-3p, miR-92a-1-5p, miR-590-5p, miR-324-5p, miR-25-3p, miR-181b-5p, miR-421, and miR-93-5p were correlated to the shorter survival time in patients with EC. In addition, we constructed a risk prediction model based on the levels of nine differentially expressed miRNAs (DEMs) and found that the OS time of EC patients with high-risk score was shorter than that of EC patients with low-risk score. Furthermore, our results showed that the risk prediction scores of EC samples were higher than those of normal samples. Finally, the area under the curve (AUC) was used to analyze the risk characteristics of EC and normal controls. By calculating the AUC and the calibration curve, the RNA signature showed a good performance. Bioinformatics analysis showed that nine DEMs were associated with several crucial signaling, including p53, FoxO, PI3K-Akt, HIF-1, and TORC1 signaling. Finally, 14 messenger RNAs (mRNAs) were identified as hub targets of nine miRNAs, including BTRC, SIAH1, RNF138, CDC27, NEDD4L, MKRN1, RLIM, FBXO11, RNF34, MYLIP, FBXW7, RNF4, UBE3C, and RNF111. TCGA dataset validation showed that these hub targets were significantly differently expressed in EC tissues compared with normal samples.

CONCLUSION

We have constructed maps and nomograms of nine-miRna risk signals associated with EC prognosis. Bioinformatics analysis revealed that nine DEMs were associated with several crucial signaling, including p53, FoxO, PI3K-Akt, HIF-1, and TORC1 signaling, in EC. We think that this study will provide clinicians with an effective decision-making tool.

Identification of Molecular Characteristics and New Prognostic Targets for Thymoma by Multiomics Analysis

BACKGROUND

Thymoma is a heterogeneous tumor originated from thymic epithelial cells. The molecular mechanism of thymoma remains unclear.

METHODS

The expression profile, methylation, and mutation data of thymoma were obtained from TCGA database. The coexpression network was constructed using the variance of gene expression through WGCNA. Enrichment analysis using clusterProfiler R package and overall survival (OS) analysis by Kaplan-Meier method were carried out for the intersection of differential expression genes (DEGs) screened by limma R package and important module genes. PPI network was constructed based on STRING database for genes with significant impact on survival. The impact of key genes on the prognosis of thymoma was evaluated by ROC curve and Cox regression model. Finally, the immune cell infiltration, methylation modification, and gene mutation were calculated.

RESULTS

We obtained eleven coexpression modules, and three of them were higher positively correlated with thymoma. DEGs in these three modules mainly involved in MAPK cascade and PPAR pathway. LIPE, MYH6, ACTG2, KLF4, SULT4A1, and TF were identified as key genes through the PPI network. AUC values of LIPE were the highest. Cox regression analysis showed that low expression of LIPE was a prognostic risk factor for thymoma. In addition, there was a high correlation between LIPE and T cells. Importantly, the expression of LIPE was modified by methylation. Among all the mutated genes, GTF2I had the highest mutation frequency.

CONCLUSION

These results suggested that the molecular mechanism of thymoma may be related to immune inflammation. LIPE may be the key genes affecting prognosis of thymoma. Our findings will help to elucidate the pathogenesis and therapeutic targets of thymoma.

The Interaction Analysis of SNP Variants and DNA Methylation Identifies Novel Methylated Pathogenesis Genes in Congenital Heart Diseases

Congenital heart defect (CHD) is a rare and complicated disease with a high mortality rate. Its etiology remains unclear and includes many aspects. DNA methylation has been indicated to be involved in heart development in the early stage of life, and aberrant methylation level was related to CHDs. This study provides the first evidence of the cross talk of SNP variants and DNA methylation in clarifying CHD underlying genomic cause. We gathered whole exome sequencing (WES) data for Group 1 consisting of patients with PA (n = 78), TOF (n = 20), TAPVC (n = 78), and PDA (n = 40), and 100 healthy children as control group. Rare non-synonymous mutations and novel genes were found and highlighted. Meanwhile, we carried out the second analysis of DNA methylation data from patients with PA (n = 3), TAPVC (n = 3), TOF (n = 3), and PDA (n = 2), and five healthy controls using 850 K array in Group 2. DNA methylation was linked to WES data, and we explored an obvious overlap of hyper/hypomethylated genes. Next, we identified some candidate genes by Fisher’s exact test and Burden analysis; then, those methylated genes were figured out by the criteria of the mutation located in the CpG islands of the genome, differential methylation sites (DMS), and DNA methylation quantitative trait loci (meQTLs) in the database, respectively. Also, the interaction of differentially methylated candidate genes with known CHD pathogenetic genes was depicted in a molecular network. Taken together, our findings show that nine novel genes (ANGPTL4, VEGFA, PAX3, MUC4, HLA-DRB1, TJP2, BCR, PKD1, and HK2) in methylation level are critical to CHD and reveal a new insight into the molecular pathogenesis of CHD.

Hippo Pathway Core Genes Based Prognostic Signature and Immune Infiltration Patterns in Lung Squamous Cell Carcinoma

BACKGROUND

We investigated the prognostic effects and their patterns of immune infiltration of hippo pathway core genes in lung squamous cell carcinoma, in order to find some clues for underlying mechanisms of LUSC tumorigenesis and help developing new therapeutic methods.

METHODS

The mutational data, transcriptome data and corresponding clinical medical information of LUSC patients were extracted from The Cancer Genome Atlas (TCGA) database. Differential expression genes (DEGs) and Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were explored. Survival analysis for the hippo core genes and the prognostic model were performed. Immune infiltration was estimated by CIBERSORT algorithm and some immune checkpoints-related genes were further investigated.

RESULTS

Overall, 551 LUSC samples were included in our study, consisting of 502 LUSC tumor samples and 49 adjacent normal samples, respectively. There were 1910 up-regulated DEGs and 2253 down-regulated DEGs were finally identified. The top five mutational hippo pathway core genes were LATS1 (4%), WWC1 (2%), TAOK1 (2%), TAOK3 (2%), and TAOK2 (2%), respectively. the mutation of LATS2 was highly associated with co-mutational NF2 (P <0.05) and TAOK1 (P <0.05). In survival analyses, we found only WWC1 (log-rank p = 0.046, HR = 1.32, 95% CI = 1-1.73) and LATS2 (log-rank p = 0.013, HR = 1.41, 95%CI = 1.08-1.86) had significant prognostic roles. After getting the three subgroups according to the subtyping results, we demonstrated that T cell gamma delta (p = 5.78e-6), B cell memory (p = 4.61e-4) and T cell CD4+ memory resting (p = 2.65e-5) had significant differences among the three groups. SIGLEC15 (P <0.01) and CD274 (P <0.05) also had statistical differences among the three subgroups.

CONCLUSIONS

Our study verified the prognostic roles of WWC1 and LATS2 in LUSC patients. Immune checkpoints-related genes SIGLEC15 and CD274 had statistical differences among the three subgroups, which may provide new perceptions on the molecular mechanisms in LUSC and maybe helpful for precisely selecting specific LUSC patients with potential immunotherapy benefits.

Identification of Molecular Regulatory Features and Markers for Acute Type A Aortic Dissection

BACKGROUND

Acute type A aortic dissection (ATAAD) is one of the most lethal cardiovascular diseases, and its molecular mechanism remains unclear.

METHODS

Differentially expressed genes (DEGs) between ATAAD and control were detected by limma R package in GSE52093, GSE153434, GSE98770, and GSE84827, respectively. The coexpression network of DEGs was identified by the WGCNA package. Enrichment analysis was performed for module genes that were positively correlated with ATAAD using clusterProfiler R package. In addition, differentially methylated markers between aortic dissection and control were identified by ChAMP package. After comparing with ATAAD-related genes, a protein-protein interaction (PPI) network was established based on the STRING database. The genes with the highest connectivity were identified as hub genes. Finally, differential immune cell infiltration between ATAAD and control was identified by ssGSEA.

RESULTS

From GSE52093 and GSE153434, 268 module genes were obtained with consistent direction of differential expression and high correlation with ATAAD. They were significantly enriched in T cell activation, HIF-1 signaling pathway, and cell cycle. In addition, 2060 differentially methylated markers were obtained from GSE84827. Among them, 77 methylation markers were ATAAD-related DEGs. Using the PPI network, we identified MYC, ITGA2, RND3, BCL2, and PHLPP2 as hub genes. Finally, we identified significantly differentially infiltrated immune cells in ATAAD.

CONCLUSION

The hub genes we identified may be regulated by methylation and participate in the development of ATAAD through immune inflammation and oxidative stress response. The findings may provide new insights into the molecular mechanisms and therapeutic targets for ATAAD.

Identification of Inflammatory Genes, Pathways, and Immune Cells in Necrotizing Enterocolitis of Preterm Infant by Bioinformatics Approaches

Background

Necrotizing enterocolitis (NEC) is one of the most serious gastrointestinal disease-causing high morbidity and mortality in premature infants. However, the underlying mechanism of the pathogenesis of NEC is still not fully understood.

Methods

RNA sequencing of intestinal specimens from 9 NEC and 5 controls was employed to quantify the gene expression levels. RNA sequencing was employed to quantify the gene expression levels. DESeq2 tool was used to identify the differentially expressed genes. The biological function, pathways, transcription factors, and immune cells dysregulated in NEC were characterized by gene set enrichment analysis.

Results

In the present study, we analyzed RNA sequencing data of NECs and controls and revealed that immune-related pathways were highly activated, while some cellular responses to external stimuli-related pathways were inactivated in NEC. Moreover, B cells, macrophages M1, and plasma cells were identified as the major cell types involved in NEC. Furthermore, we also found that inflammation-related transcription factor genes, such as STAT1, STAT2, and IRF2, were significantly activated in NEC, further suggesting that these TFs might play critical roles in NEC pathogenesis. In addition, NEC samples exhibited heterogeneity to some extent. Interestingly, two subgroups in the NEC samples were identified by hierarchical clustering analysis. Notably, B cells, T cells, Th1, and Tregs involved in adaptive immune were predicted to highly infiltrate into subgroup I, while subgroup II was significantly infiltrated by neutrophils. The heterogeneity of immune cells in NEC indicated that both innate and adaptive immunes might induce NEC-related inflammatory response.

Conclusions

In summary, we systematically analyzed inflammation-related genes, signaling pathways, and immune cells to characterize the NEC pathogenesis and samples, which greatly improved our understanding of the roles of inflammatory responses in NEC.

Gene Expression Signature of Traumatic Brain Injury

Background: Traumatic brain injury (TBI) is a brain function change caused by external forces, which is one of the main causes of death and disability worldwide. The aim of this study was to identify early diagnostic markers and potential therapeutic targets for TBI.

Methods: Differences between TBI and controls in GSE89866 and GSE104687 were analyzed. The two groups of differentially expressed genes (DEGs) were combined for coexpression analysis, and the modules of interest were performed using enrichment analysis. Hub genes were identified by calculating area under curve (AUC) values of module genes, PPI network analysis, and functional similarity. Finally, the difference in immune cell infiltration between TBI and control was calculated by ssGSEA.

Results: A total of 4,817 DEGs were identified in GSE89866 and 1,329 DEGs in GSE104687. They were clustered into nine modules. The genes of modules 1, 4, and 7 had the most crosstalk and were identified as important modules. Enrichment analysis revealed that they were mainly associated with neurodevelopment and immune inflammation. In the PPI network constructed by genes with top 50 AUC values in module genes, we identified the top 10 genes with the greatest connectivity. Among them, down-regulated RPL27, RPS4X, RPL23A, RPS15A, and RPL7A had similar functions and were identified as hub genes. In addition, DC and Tem were significantly up-regulated and down-regulated between TBI and control, respectively.

Conclusion: We found that hub genes may have a diagnostic role for TBI. Molecular dysregulation mechanisms of TBI are associated with neurological and immune inflammation. These results may provide new ideas for the diagnosis and treatment of TBI.

Comprehensive Analysis of the Prognostic Role and Mutational Characteristics of m6A-Related Genes in Lung Squamous Cell Carcinoma

BACKGROUND

There have been limited treatment therapies for lung squamous cell carcinoma (LUSC). M6A-related genes may be the next therapeutic targets for LUSC. In this study, we explored the prognostic role and mutational characteristics of m6A-related genes in LUSC.

METHODS

LUSC gene expression data, mutational data, and corresponding clinical information were extracted from The Cancer Genome Atlas database. Differentially expressed genes (DEGs) were identified, and the mutation characteristics of LUSC patients were explored. Then, m6A-related genes were extracted and the correlations among the genes were detected. Finally, the prognostic roles of the genes were investigated and the nomogram model was developed. Besides, the protein-protein interaction (PPI) network was used to explore the potential interactions among the genes.

RESULTS

In total, there are 551 LUSC samples enrolled in our study, containing 502 LUSC tumor samples and 49 adjacent normal LUSC samples, respectively. There were 2970 upregulated DEGs and 1806 downregulated DEGs were further explored. IGF2BP1 and RBM15 had significant co-occurrence frequency (p < 0.05). Besides, METTL14 and ZC3H13 or YTHDF3 also had significant co-occurrence frequency (p < 0.05). All the m6A-related genes represent the positive correlation. WTAP was identified as a prognostic gene in the TCGA database while YTHDC1 and YTHDF1 were identified as prognostic genes. In multivariate Cox analysis, YTHDF1, age, pN stage, pTNM stage, and smoking were all identified as significant prognostic factors for OS.

CONCLUSION

We investigated the expression patterns and mutational characteristics of LUSC patients and identified three potential independent prognostic m6A-related genes (WTAP, YTHDC1, and YTHDF1) for OS in LUSC patients.

Identification of Potential Therapeutic Targets and Molecular Regulatory Mechanisms for Osteoporosis by Bioinformatics Methods

BACKGROUND

Osteoporosis is characterized by low bone mass, deterioration of bone tissue structure, and susceptibility to fracture. New and more suitable therapeutic targets need to be discovered.

METHODS

We collected osteoporosis-related datasets (GSE56815, GSE99624, and GSE63446). The methylation markers were obtained by differential analysis. Degree, DMNC, MCC, and MNC plug-ins were used to screen the important methylation markers in PPI network, then enrichment analysis was performed. ROC curve was used to evaluate the diagnostic effect of osteoporosis. In addition, we evaluated the difference in immune cell infiltration between osteoporotic patients and control by ssGSEA. Finally, differential miRNAs in osteoporosis were used to predict the regulators of key methylation markers.

RESULTS

A total of 2351 differentially expressed genes and 5246 differentially methylated positions were obtained between osteoporotic patients and controls. We identified 19 methylation markers by PPI network. They were mainly involved in biological functions and signaling pathways such as apoptosis and immune inflammation. HIST1H3G, MAP3K5, NOP2, OXA1L, and ZFPM2 with higher AUC values were considered key methylation markers. There were significant differences in immune cell infiltration between osteoporotic patients and controls, especially dendritic cells and natural killer cells. The correlation between MAP3K5 and immune cells was high, and its differential expression was also validated by other two datasets. In addition, NOP2 was predicted to be regulated by differentially expressed hsa-miR-3130-5p.

CONCLUSION

Our efforts aim to provide new methylation markers as therapeutic targets for osteoporosis to better treat osteoporosis in the future.

m6A Regulators Is Differently Expressed and Correlated With Immune Response of Esophageal Cancer

N6 methyladenosine (m6A) RNA methylation regulators play an important role in the development of tumors. However, their function in esophageal cancer (EC) has not been fully elucidated. Here, we analyzed the gene expression data of 24 major m6A RNA methylation regulators from 775 patients with EC from TCGA dataset. The present study showed the aberrations of m6A regulators in genome were correlated to prognosis in human ECs. Meanwhile, 17 m6A regulators showed increased expression in EC samples, including YTHDC1, IGF2BP2, FTO, METTL14, YTHDF3, RBM15, WTAP, HNRNPA2B1, HNRNPC, ALKBH5, YTHDF2, METTL16, IGF2BP3, VIRMA, RBM15B, YTHDF1, KIAA1429, HAKAI, and ZC3H13. Among them, we found HNRNPC, YTHDC2, WTAP, VIRMA, IGF2BP3, and HNRNPA2B1 were significantly correlated to worse outcomes and advanced stage in EC. Furthermore, we showed levels of m6A regulators is correlated with the expression of Immuno-regulators (Immunoinhibitors, Immunostimulators, and MHC molecules) and immune infiltration levels in EC. Bioinformatics further confirm m6A regulators were involved in regulating RNA splicing, RNA stability, and cell proliferation. Our study showed m6A regulators are promising targets and biomarkers for cancer immunotherapy in EC.

Identification of Common Genes and Pathways in Eight Fibrosis Diseases

Acute and chronic inflammation often leads to fibrosis, which is also the common and final pathological outcome of chronic inflammatory diseases. To explore the common genes and pathogenic pathways among different fibrotic diseases, we collected all the reported genes of the eight fibrotic diseases: eye fibrosis, heart fibrosis, hepatic fibrosis, intestinal fibrosis, lung fibrosis, pancreas fibrosis, renal fibrosis, and skin fibrosis. We calculated the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment scores of all fibrotic disease genes. Each gene was encoded using KEGG and GO enrichment scores, which reflected how much a gene can affect this function. For each fibrotic disease, by comparing the KEGG and GO enrichment scores between reported disease genes and other genes using the Monte Carlo feature selection (MCFS) method, the key KEGG and GO features were identified. We compared the gene overlaps among eight fibrotic diseases and connective tissue growth factor (CTGF) was finally identified as the common key molecule. The key KEGG and GO features of the eight fibrotic diseases were all screened by MCFS method. Moreover, we interestingly found overlaps of pathways between renal fibrosis and skin fibrosis, such as GO:1901890-positive regulation of cell junction assembly, as well as common regulatory genes, such as CTGF, which is the key molecule regulating fibrogenesis. We hope to offer a new insight into the cellular and molecular mechanisms underlying fibrosis and therefore help leading to the development of new drugs, which specifically delay or even improve the symptoms of fibrosis.

Integrated WGCNA and PPI Network to Screen Hub Genes Signatures for Infantile Hemangioma

BACKGROUND

Infantile hemangioma (IH) is characterized by proliferation and regression.

METHODS

Based on the GSE127487 dataset, the differentially expressed genes (DEGs) between 6, 12, or 24 months and normal samples were screened, respectively. STEM software was used to screen the continued up-regulated or down-regulated in common genes. The modules were assessed by weighted gene co-expression network analysis (WGCNA). The enrichment analysis was performed to identified the biological function of important module genes. The area under curve (AUC) value and protein-protein interaction (PPI) network were used to identify hub genes. The differential expression of hub genes in IH and normal tissues was detected by qPCR.

RESULTS

There were 5,785, 4,712, and 2,149 DEGs between 6, 12, and 24 months and normal tissues. We found 1,218 DEGs were up-regulated or down-regulated expression simultaneously in common genes. They were identified as 10 co-expression modules. Module 3 and module 4 were positively or negatively correlated with the development of IH, respectively. These two module genes were significantly involved in immunity, cell cycle arrest and mTOR signaling pathway. The two module genes with AUC greater than 0.8 at different stages of IH were put into PPI network, and five genes with the highest degree were identified as hub genes. The differential expression of these genes was also verified by qRTPCR.

CONCLUSION

Five hub genes may distinguish for proliferative and regressive IH lesions. The WGCNA and PPI network analyses may help to clarify the molecular mechanism of IH at different stages.

Identification of driver genes with aberrantly alternative splicing events in pediatric patients with retinoblastoma

Retinoblastoma (RB) is one of the most common cancer in children. However, the specific mechanism about RB tumorigenesis has not been fully understood. In this study, to comprehensively characterize the splicing alterations in the tumorigenesis of RB, we analyzed the differential alternative splicing events in RB. Specifically, the isoforms of RB1 were downregulated in the RB samples, and a large proportion of differentially expressed genes had multiple differentially expressed transcripts (64%). We identified 1453 genes with differential alternative splicing, among which, SE accounted for the majority, followed by MXE, RI, A3SS, and A5SS. Furthermore, the biological function related to the normal function of eyes, and E2F family TFs were significantly enriched by the genes with differential alternative splicing. Among the genes associated with visual sense, ABCA4 was found to have two mutually exclusive exons, resulting in two isoforms with different functionalities. Notably, DAZAP1 was identified as one of the critical splicing factors in RB, which was potentially involved in E2F and RB pathways. Functionally, differential binding sites in DAZAP1 protein were significantly observed between RB and normal samples. Based on the comprehensive analysis of the differential alternative splicing events and splicing factors, we identified some driver genes with differential alternative splicing and critical splicing factors involved in RB, which would greatly improve our understanding of the alternative splicing process in the tumorigenesis of RB.

Recognition of Immune Microenvironment Landscape and Immune-Related Prognostic Genes in Breast Cancer

Background

Breast cancer (BC) is the most common malignant tumor in women. The immunophenotype of tumor microenvironment (TME) has shown great therapeutic potential in tumor.

Method

The transcriptome was obtained from TCGA and GEO data. Immune infiltration was analyzed by single-sample gene set enrichment (ssGSEA). The immune feature was constructed by Cox regression analysis. In addition, the coexpression of differential expression genes (DEGs) was identified. Through enrichment analysis, the function and pathway of module genes were identified. The somatic mutations related to immune characteristics were analyzed by Maftools. By using the consistency clustering algorithm, the molecular subtypes were constructed, and the overall survival time (OS) was predicted.

Results

Immune landscape can be divided into low immune infiltration and high immune infiltration. Cox regression analysis identified seven immune cells as protective factors of BC. In the coexpression modules for DEGs of high and low immune infiltration, module 1 was related to T cells and high immune infiltration. In particular, the area under the curve (AUC) value of hub gene SASH3 was the highest, and the correlation with T cells was stronger in the high immune infiltration. Enrichment analysis found that oxidative stress, T cell aggregation, and apoptosis were observed in high immune infiltration. In addition, TP53 was identified as the most important somatic gene mutation related to immune characteristics. Importantly, we also constructed seven immune cell-based breast cancer subtypes to predict OS.

Conclusion

We evaluated the immune landscape of BC and constructed the gene characteristics related to the immune landscape. The potential of T cells and SASH3 in immunotherapy of BC was revealed, which may guide the development of new clinical treatment strategies.

Next Generation Sequencing Identify Rare Copy Number Variants in Non-syndromic Patent Ductus Arteriosus

Patent ductus arteriosus (PDA) is a common congenital cardiovascular malformation with both inherited and acquired causes. Several genes have been reported to be related to PDA, but the molecular pathogenesis is still unclear. Here, we screened a population matched cohort of 39 patients with PDA and 100 healthy children using whole exome sequencing (WES). And identified 10 copy number variants (CNVs) and 20 candidate genes using Gene ontology (GO) functional enrichment analysis. In gene network analysis, we screened 7 pathogenic CNVs of 10 candidate genes (MAP3K1, MYC, VAV2, WDR5, RXRA, APLNR, TJP1, ERCC2, FOSB, CHRNA4). Further analysis of transcriptome array showed that 7 candidate genes (MAP3K1, MYC, VAV2, APLNR, TJP1, FOSB, CHRNA4) were indeed significantly expressed in human embryonic heart. Moreover, CHRNA4 was observed the most important genes. Our data provided rare CNVs as potential genetic cause of PDA in humans and also advance understanding of the genetic components of PDA.

Identification of Key Genes Involved in Acute Myocardial Infarction by Comparative Transcriptome Analysis

Background

Acute myocardial infarction (AMI) is regarded as an urgent clinical entity, and identification of differentially expressed genes, lncRNAs, and altered pathways shall provide new insight into the molecular mechanisms behind AMI.

Materials and Methods

Microarray data was collected to identify key genes and lncRNAs involved in AMI pathogenesis. The differential expression analysis and gene set enrichment analysis (GSEA) were employed to identify the upregulated and downregulated genes and pathways in AMI. The protein-protein interaction network and protein-RNA interaction analysis were utilized to reveal key long noncoding RNAs.

Results

In the present study, we utilized gene expression profiles of circulating endothelial cells (CEC) from 49 patients of AMI and 50 controls and identified a total of 552 differentially expressed genes (DEGs). Based on these DEGs, we also observed that inflammatory response-related genes and pathways were highly upregulated in AMI. Mapping the DEGs to the protein-protein interaction (PPI) network and identifying the subnetworks, we found that OMD and WDFY3 were the hub nodes of two subnetworks with the highest connectivity, which were found to be involved in circadian rhythm and organ- or tissue-specific immune response. Furthermore, 23 lncRNAs were differentially expressed between AMI and control groups. Specifically, we identified some functional lncRNAs, including XIST and its antisense RNA, TSIX, and three lncRNAs (LINC00528, LINC00936, and LINC01001), which were predicted to be interacting with TLR2 and participate in Toll-like receptor signaling pathway. In addition, we also employed the MMPC algorithm to identify six gene signatures for AMI diagnosis. Particularly, the multivariable SVM model based on the six genes has achieved a satisfying performance (AUC = 0.97).

Conclusion

In conclusion, we have identified key regulatory lncRNAs implicated in AMI, which not only deepens our understanding of the lncRNA-related molecular mechanism of AMI but also provides computationally predicted regulatory lncRNAs for AMI researchers.

Transcriptome Analysis Identifies Novel Prognostic Genes in Osteosarcoma

Osteosarcoma (OS), a malignant primary bone tumor often seen in young adults, is highly aggressive. The improvements in high-throughput technologies have accelerated the identification of various prognostic biomarkers for cancer survival prediction. However, only few studies focus on the prediction of prognosis in OS patients using gene expression data due to small sample size and the lack of public datasets. In the present study, the RNA-seq data of 82 OS samples, along with their clinical information, were collected from the TARGET database. To identify the prognostic genes for the OS survival prediction, we selected the top 50 genes of contribution as the initial candidate genes of the prognostic risk model, which were ranked by random forest model, and found that the prognostic model with five predictors including CD180, MYC, PROSER2, DNAI1, and FATE1 was the optimal multivariable Cox regression model. Moreover, based on a multivariable Cox regression model, we also developed a scoring method and stratified the OS patients into groups of different risks. The stratification for OS patients in the validation set further demonstrated that our model has a robust performance. In addition, we also investigated the biological function of differentially expressed genes between two risk groups and found that those genes were mainly involved with biological pathways and processes regarding immunity. In summary, the identification of novel prognostic biomarkers in OS would greatly assist the prediction of OS survival and development of molecularly targeted therapies, which in turn benefit patients' survival.

Network-Based Coexpression Analysis Identifies Functional and Prognostic Long Noncoding RNAs in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a primary liver cancer associated with a growing incidence and extremely high mortality. However, the pathogenic mechanism is still not fully understood. In the present study, we identified 1,631 upregulated and 1,515 downregulated genes and found that cell cycle and metabolism-related pathways or biological processes highly dysregulated in HCC. To assess the biological importance of these DEGs, we carried out weighted gene coexpression network analysis (WGCNA) to identify the functional modules potentially involved in HCC pathogenesis or progression. The five modules were detected with Dynamic Tree Cut algorithm, and GO enrichment analysis revealed that these modules exhibited different biological processes or signaling pathways, such as metabolism-related pathways, cell proliferation-related pathways, and molecules in tumor microenvironment. Moreover, we also observed two immune cells, namely, cytotoxic cells and macrophage enriched in modules grey and brown, respectively, while T helper cell-2 (Th2) was enriched in module turquoise. Among the WGCNA network, four hub long noncoding RNAs (lncRNAs) were identified to be associated with HCC prognostic outcomes, suggesting that coexpression network analysis could uncover lncRNAs with functional importance, which may be associated with prognostic outcomes of HCC patients. In summary, this study demonstrated that network-based analysis could identify some functional modules and some hub-lncRNAs, which may be critical for HCC pathogenesis or progression.

TEAD4 promotes tumor development in patients with lung adenocarcinoma via ERK signaling pathway

OBJECTIVE

Whether TEAD4 itself plays a vital role in the tumorigenesis and development of lung adenocarcinoma remains unclear. In our study, we aim to investigate the expression pattern and biological functions of TEAD4 and further investigate the potential mechanisms.

METHODS

Clinical tumor and paired normal samples were collected for preparing tissue microarray. Western blot and immunohistochemical (IHC) staining of TEAD4 expression in these tissues were conducted to explore the expression pattern. Moreover, A549 cell line was select for investigating the function of TEAD4 for lung adenocarcinoma in vitro and in vivo. RNA sequencing was finally performed to further detect the potential downstream genes.

RESULTS

The elevated TEAD4 expression level was observed in tumor tissues and the patients with higher TEAD4 expression tended to have worse overall survival. The knockdown of TEAD4 inhibits A549 cells proliferation ability and migration ability. A total of 431 differentially expressed genes (DEGs), including 239 down-regulated genes and 191 up-regulated genes, were finally identified and some of DEGs were validated. Moreover, knockdown of TEAD4 led to the down-regulation of pERK, which maybe the potential TEAD4-targeted signaling pathway to play the pro-tumorigenic function.

CONCLUSIONS

The expression level of TEAD4 is high in lung adenocarcinoma tumor tissues and positively associated with worse prognosis. Up-regulation of TEAD4 may lead to excessive transcription and phosphorylation of ERK proteins and therefore accelerates the process of tumor development. Our results demonstrate that overexpression of TEAD4 is a new mechanism of dysregulation of Hippo pathway.

Identification of novel candidate genes in heterotaxy syndrome patients with congenital heart diseases by whole exome sequencing

Heterotaxy syndrome (HS) involves dysfunction of multiple systems resulting from abnormal left-right (LR) body patterning. Most HS patients present with complex congenital heart diseases (CHD), the disability and mortality of HS patients are extremely high. HS has great heterogeneity in phenotypes and genotypes, which have rendered gene discovery challenging. The aim of this study was to identify novel genes that underlie pathogenesis of HS patients with CHD. Whole exome sequencing was performed in 25 unrelated HS cases and 100 healthy controls; 19 nonsynonymous variants in 6 novel candidate genes (FLNA, ITGA1, PCNT, KIF7, GLI1, KMT2D) were identified. The functions of candidate genes were further analyzed in zebrafish model by CRISPR/Cas9 technique. Genome-editing was successfully introduced into the gene loci of flna, kmt2d and kif7, but the phenotypes were heterogenous. Disruption of each gene disturbed normal cardiac looping while kif7 knockout had a more prominent effect on liver budding and pitx2 expression. Our results revealed three potential HS pathogenic genes with probably different molecular mechanisms.