Identification of EGFLAM, SPATC1L and RNASE13 as novel susceptibility loci for aortic aneurysm in Japanese individuals by exome-wide association studies

EGFLAM exhibits oncogenic activity and shows promise as a prognostic biomarker and therapeutic target in glioblastoma

Glioblastoma (GBM) remains the most lethal primary brain tumor, characterized by dismal survival rates. Novel molecular targets are urgently required to enhance therapeutic outcomes. A combination of bioinformatics analysis and experimental validation was employed to investigate the role of EGFLAM in GBM. The Chinese Glioma Genome Atlas provided a platform for gene expression profiling, while siRNA-mediated knockdown and overexpression assays in GBM cell lines, alongside in vivo tumorigenesis models, facilitated functional validation. EGFLAM was found to be significantly overexpressed in GBM tissues, correlating with adverse prognostic factors and higher tumor grades, particularly in patients over the age of 41. Functional assays indicated that EGFLAM is vital for maintaining GBM cell proliferation, viability, and invasiveness. Knockdown of EGFLAM expression led to a marked decrease in tumorigenic capabilities. Proteomic interactions involving EGFLAM, such as with NUP205, were implicated in cell cycle regulation, providing insight into its oncogenic mechanism. In vivo studies further demonstrated that silencing EGFLAM expression could inhibit tumor growth, underscoring its therapeutic potential. The study identifies EGFLAM as a pivotal oncogenic factor in GBM, serving as both a prognostic biomarker and a viable therapeutic target. These findings lay the groundwork for future research into EGFLAM-targeted therapies, aiming to improve clinical outcomes for GBM patients.

Role of Pericytes in Cardiomyopathy-Associated Myocardial Infarction Revealed by Multiple Single-Cell Sequencing Analysis

Acute myocardial infarction (AMI) is one of the leading causes of cardiovascular death worldwide. AMI with cardiomyopathy is accompanied by a poor long-term prognosis. However, limited studies have focused on the mechanism of cardiomyopathy associated with AMI. Pericytes are important to the microvascular function in the heart, yet little attention has been paid to their function in myocardial infarction until now. In this study, we integrated single-cell data from individuals with cardiomyopathy and myocardial infarction (MI) GWAS data to reveal the potential function of pericytes in cardiomyopathy-associated MI. We found that pericytes were concentrated in the left atrium and left ventricle tissues. DLC1/GUCY1A2/EGFLAM were the top three uniquely expressed genes in pericytes (p < 0.05). The marker genes of pericytes were enriched in renin secretion, vascular smooth muscle contraction, gap junction, purine metabolism, and diabetic cardiomyopathy pathways (p < 0.05). Among these pathways, the renin secretion and purine metabolism pathways were also found in the process of MI. In cardiomyopathy patients, the biosynthesis of collagen, modulating enzymes, and collagen formation were uniquely negatively regulated in pericytes compared to other cell types (p < 0.05). COL4A2/COL4A1/SMAD3 were the hub genes in pericyte function involved in cardiomyopathy and AMI. In conclusion, this study provides new evidence about the importance of pericytes in the pathogenesis of cardiomyopathy-associated MI. DLC1/GUCY1A2/EGFLAM were highly expressed in pericytes. The hub genes COL4A2/COL4A1/SMAD3 may be potential research targets for cardiomyopathy-associated MI.

Genome-Wide Association Study of Alopecia Areata in Taiwan: The Conflict Between Individuals and Hair Follicles

Purpose

Alopecia areata (AA) is one of the most prevalent autoimmune diseases affecting humans. Given that hair follicles are immune-privileged, autoimmunity can result in disfiguring hair loss. However, the genetic basis for AA in the Taiwanese population remains unknown.

Materials and Methods

A genome-wide association study was conducted using a cohort of 408 AA cases and 8167 controls. To link variants to gene relationships, we used 882 SNPs (P<1E-05) within 74 genes that were associated with AA group to build the biological networks by IPA software. HLA diplotypes and haplotypes were analyzed using Attribute Bagging (HIBAG)-R package and chi-square analysis.

Results

Seven single nucleotide polymorphisms (SNPs) including LINC02006 (rs531166736, rs187306735), APC (rs112800832_C_CAT), SRP19 (rs139948960, rs144784670), EGFLAM (rs16903975) and LDLRAD3 (rs79874564) were closely associated with the AA phenotype (P<5E-08). Examination of biological networks revealed that these genomic areas are associated with antigen presentation signaling, B cell and T cell development, Th1 and Th2 activation pathways, Notch signaling, crosstalk signaling between dendritic cells and natural killer cells, and phagosome maturation. Based on human leukocyte antigen (HLA) genotype analysis, four HLA genotypes (HLA-B*15:01-*40:01, HLA-DQA1*01:02-*03:03, HLA-DQA1*01:02, and HLA-DQB1*02:01) were found to be associated with AA (adjusted p-value<0.05). HLA-DQA1*01:02 is the most significantly related gene in the Taiwanese population (adjusted p-value = 2.09E-05).

Conclusion

This study successfully identified susceptibility loci associated with AA in the Taiwanese population. These findings not only shed light on the origins of AA within the Taiwanese context but also contribute to a comprehensive understanding of the genetic factors influencing AA susceptibility.

Whole Exome Sequencing of a Patient with a Milder Phenotype of Xeroderma Pigmentosum Group C

A 17-year-old female Korean patient (XP115KO) was previously diagnosed with Xeroderma pigmentosum group C (XPC) by Direct Sanger sequencing, which revealed a homozygous nonsense mutation in the XPC gene (rs121965088: c.1735C > T, p.Arg579Ter). While rs121965088 is associated with a poor prognosis, our patient presented with a milder phenotype. Hence, we conducted whole-exome sequencing in the patient and her family members to detect coexisting mutations that may have resulted in a milder phenotype of rs121965088 through genetic interaction. Materials and Methods: the whole-exome sequencing analysis of samples obtained from the patient and her family members (father, mother, and brother) was performed. To identify the underlying genetic cause of XPC, the extracted DNA was analyzed using Agilent's SureSelect XT Human All Exon v5. The functional effects of the resultant variants were predicted using the SNPinfo web server, and structural changes in the XPC protein using the 3D protein modeling program SWISS-MODEL. Results: Eight biallelic variants, homozygous in the patient and heterozygous in her parents, were detected. Four were found in the XPC gene: one nonsense variant (rs121965088: c.1735C > T, p.Arg579Ter) and three silent variants (rs2227998: c.2061G > A, p. Arg687Arg; rs2279017: c.2251-6A > C, intron; rs2607775: c.-27G > C, 5'UTR). The remaining four variants were found in non-XP genes, including one frameshift variant [rs72452004 of olfactory receptor family 2 subfamily T member 35 (OR2T35)], three missense variants [rs202089462 of ALF transcription elongation factor 3 (AFF3), rs138027161 of TCR gamma alternate reading frame protein (TARP), and rs3750575 of annexin A7 (ANXA7)]. Conclusions: potential candidates for genetic interactions with rs121965088 were found. The rs2279017 and rs2607775 of XPC involved mutations in the intron region, which affected RNA splicing and protein translation. The genetic variants of AFF3, TARP, and ANXA7 are all frameshift or missense mutations, inevitably disturbing the translation and function of the resultant proteins. Further research on their functions in DNA repair pathways may reveal undiscovered cellular relationships within xeroderma pigmentosum.

Genome-wide association study reveals new QTL and functional candidate genes for the number of ribs and carcass length in pigs

The number of ribs (NR) and carcass length (CL) are important economic traits in pig breeding programs. Pigs with a higher NR and longer CL produce greater pork yields. In the present study, Suhuai pigs with NR and CL phenotypes were genotyped using the Neogen^® GGP Porcine 80 K SNP array to identify the QTL affecting NR and CL and dissect the candidate genes for the two traits. The SNP-chip data was imputed to the whole-genome sequence (iWGS) to increase the probability of identifying causal variants. Through genome-wide association studies (GWAS) based on both chip and iWGS data, significant SNPs were detected on Sus scrofa chromosome (SSC) 1, SSC4 and SSC7 for NR and on SSC5, SSC16 and SSC17 for CL. Moreover, two SNPs (H3GA0022644 and WU_10.2_7_103460706) on SSC7 detected in chip-based GWAS were significantly associated with both NR and CL. Through Bayes fine mapping, one reported QTL for NR on SSC7 and two reported QTL for CL on SSC17 were verified, and two new QTL (SSC1: 14.05-15.84 Mb and SSC4: 64.83-66.59 Mb) affecting NR and two new QTL (SSC5: 58.31-59.84 Mb and SSC16: 22.98-23.43 Mb) affecting CL were detected. According to the biological functions of genes, MTHFD1L on SSC1 and SULF1 on SSC4 are novel functional candidate genes for NR, and EMP1 on SSC5 and EGFLAM on SSC16 are novel functional candidate genes for CL. Overall, our findings provide a basis for identifying new causal genes and mutations affecting NR and CL.

Marginal Zinc Deficiency in Mice Increased the Number of Abnormal Sperm and Altered the Expression Level of Spermatogenesis-Related Genes

Marginal zinc deficiency is more common than severe zinc deficiency, and the effect of marginal zinc deficiency on male reproduction is unknown. This study investigated the effect of marginal zinc deficiency on spermatogenesis and its mechanism. Male ICR mice were fed normal zinc (30 mg/kg) and marginal zinc deficiency (10 mg/kg) diets for 35 days. Zinc-dependent proteins and enzymes were measured as biomarkers of zinc levels in the body. Metallothionein and Cu-Zn SOD levels in the control group were higher than those in the marginal zinc deficiency group. Hematoxylin and eosin staining showed that the marginal zinc deficiency diet caused histopathological changes in the testis and destruction of the sperm head under scanning electron microscopy. Sperm parameters and sex hormone levels were also affected by marginal zinc deficiency. In addition, marginal zinc deficiency led to alter expression level of several important spermatogenesis-related genes in the epididymis and testes. These results indicate that although zinc intake in marginal zinc deficiency is close to the recommended reference value, low zinc intake interferes with the expression of genes related to spermatogenesis and may lead to sperm abnormalities in mice.

Differential Expression of Super-Enhancer-Associated Long Non-coding RNAs in Uterine Leiomyomas

Super-enhancer-associated long non-coding RNAs (SE-lncRNAs) are a specific set of lncRNAs transcribed from super-enhancer (SE) genomic regions. Recent studies have revealed that SE-lncRNAs play essential roles in tumorigenesis through the regulation of oncogenes. The objective of this study was to elucidate the expression profile of SE-lncRNAs with concurrent assessment of associated mRNAs in leiomyomas and paired myometrium. Arraystar SE-lncRNAs arrays were used to systematically profile the differentially expressed SE-lncRNAs along with the corresponding SE-regulated protein coding genes in eight leiomyomas and paired myometrium. The analysis indicated 7680 SE-lncRNAs were expressed, of which 721 SE-lncRNAs were overexpressed, while 247 SE-lncRNAs were underexpressed by 1.5-fold or greater in leiomyoma. Thirteen novel SE-lncRNAs and their corresponding protein coding genes were selected, and their expression was confirmed in eighty-one paired leiomyoma tissues by quantitative real-time PCR. The thirteen pairs of SE-lncRNAs and their corresponding protein coding genes included RP11-353N14.2/CBX4, SOCS2-AS1/SOCS2, RP1-170O19.14/HOXA11, CASC15/PRL, EGFLAM-AS1/EGFLAM, RP11-225H22/NEURL1, RP5-1086K13.1/CD58, AC092839.3/SPTBN1, RP11-69I8.3/CTGF, TM4SF1-AS1/TM4SF1, RP11-373D23/FOSL2, RP11-399K21.11/COMTD1, and CTB-113P19.1/SPARC. Among these SE-lncRNAs, the expression of SOCS2-AS1/SOCS2, RP11-353N14.2/CBX4, RP1-170O19.14/HOXA11, and RP11-225H22/NEURL1 was significantly higher in African Americans as compared with Caucasians. The expression of RP11-353N14.2/CBX4, SOCS2-AS1/SOCS2, CASC15/PRL, and CTB-113P19.1/SPARC was significantly higher in tumors with MED12-mutation-positive as compared with MED12-mutation-negative tumors. Collectively, our results indicate that the differential expression of SE in leiomyomas is another mechanism contributing to dysregulation of protein coding genes in leiomyomas and that race and MED12 mutation can influence the expression of a select group of SE.

Genetic Variants Associated With Resilience in Human and Animal Studies

Resilience is broadly defined as the ability to maintain or regain functioning in the face of adversity and is influenced by both environmental and genetic factors. The identification of specific genetic factors and their biological pathways underpinning resilient functioning can help in the identification of common key factors, but heterogeneities in the operationalisation of resilience have hampered advances. We conducted a systematic review of genetic variants associated with resilience to enable the identification of general resilience mechanisms. We adopted broad inclusion criteria for the definition of resilience to capture both human and animal model studies, which use a wide range of resilience definitions and measure very different outcomes. Analyzing 158 studies, we found 71 candidate genes associated with resilience. OPRM1 (Opioid receptor mu 1), NPY (neuropeptide Y), CACNA1C (calcium voltage-gated channel subunit alpha1 C), DCC (deleted in colorectal carcinoma), and FKBP5 (FKBP prolyl isomerase 5) had both animal and human variants associated with resilience, supporting the idea of shared biological pathways. Further, for OPRM1, OXTR (oxytocin receptor), CRHR1 (corticotropin-releasing hormone receptor 1), COMT (catechol-O-methyltransferase), BDNF (brain-derived neurotrophic factor), APOE (apolipoprotein E), and SLC6A4 (solute carrier family 6 member 4), the same allele was associated with resilience across divergent resilience definitions, which suggests these genes may therefore provide a starting point for further research examining commonality in resilience pathways.

Integrative analysis of transcriptome-wide association study and mRNA expression profile identified candidate genes and pathways associated with aortic aneurysm and dissection

BACKGROUND

Aortic aneurysm and dissection (AAD) are a set of life-threatening diseases. This study aimed to investigate the genetic mechanisms of AAD by integrating transcriptome-wide association study (TWAS) and mRNA expression profile.

METHODS

The genome-wide association study (GWAS) summary data of AAD was obtained from the UK Biobank, which contains 452,264 White British individuals, including 1470 AAD patients. The TWAS analysis was performed by integrating expression quantitative trait loci (eQTL) data of aorta and the GWAS dataset of AAD using the FUSION software. The TWAS significant genes and differentially expressed genes (DEGs) identified by mRNA expression profile of aortic dissection were integrated to find common genes and biological process. For TWAS significant genes, protein-protein interaction (PPI) network analysis was further conducted based on STRING database.

RESULTS

TWAS identified 423 genes with P < 0.05. After comparing the results of TWAS and mRNA expression profile, 11 overlapping genes (PDE8B, IKBKE, HMGA1, PKM, CHST1, DUS3L, S100A16, PTGS1, RAB38, PDLIM5, NOL6) and 15 common gene ontology (GO) terms (including extracellular matrix organization, external encapsulating structure organization, cell-substrate adhesion, actin filament-based process, focal adhesion, protein kinase activity) were identified. 9 hub genes of the TWAS results were identified via PPI network analysis, including RPS9, RPS18, RSRC1, DNAJC3, HBS1L, PRKCA, NCAM1, ITGB3, FTSJ3.

CONCLUSION

Multiple candidate genes and biological processes associated with AAD were identified by the present integrative study of TWAS and mRNA expression profile. Further studies are needed to elucidate the genetic mechanisms of AAD.

Integrated bioinformatics analysis reveals novel key biomarkers and potential candidate small molecule drugs in gestational diabetes mellitus

Gestational diabetes mellitus (GDM) is the metabolic disorder that appears during pregnancy. The current investigation aimed to identify central differentially expressed genes (DEGs) in GDM. The transcription profiling by array data (E-MTAB-6418) was obtained from the ArrayExpress database. The DEGs between GDM samples and non-GDM samples were analyzed. Functional enrichment analysis were performed using ToppGene. Then we constructed the protein–protein interaction (PPI) network of DEGs by the Search Tool for the Retrieval of Interacting Genes database (STRING) and module analysis was performed. Subsequently, we constructed the miRNA–hub gene network and TF–hub gene regulatory network. The validation of hub genes was performed through receiver operating characteristic curve (ROC). Finally, the candidate small molecules as potential drugs to treat GDM were predicted by using molecular docking. Through transcription profiling by array data, a total of 869 DEGs were detected including 439 up-regulated and 430 down-regulated genes. Functional enrichment analysis showed these DEGs were mainly enriched in reproduction, cell adhesion, cell surface interactions at the vascular wall and extracellular matrix organization. Ten genes, HSP90AA1, EGFR, RPS13, RBX1, PAK1, FYN, ABL1, SMAD3, STAT3 and PRKCA were associated with GDM, according to ROC analysis. Finally, the most significant small molecules were predicted based on molecular docking. This investigation identified hub genes, signal pathways and therapeutic agents, which might help us, enhance our understanding of the mechanisms of GDM and find some novel therapeutic agents for GDM.

EGFLAM correlates with cell proliferation, migration, invasion and poor prognosis in glioblastoma

EGFLAM as a novel gene biomarker has been reported in some cancers but not glioblastoma (GBM) yet. To clarify the functional role of EGFLAM in GBM, we performed this study. Firstly, based on TCGA and Oncomine database, EGFLAM expression and clinical significance in GBM patients was analyzed. Furthermore, the biological effect of EGFLAM in GBM cells was determined by qRT-PCR, CCK-8 assay, colony formation assay, wound healing assay, transwell assays and western blot analysis. The databases analysis showed that EGFLAM expression was at higher levels in GBM patients with poor prognosis. The results indicated that EGFLAM silence inhibited the proliferation, migration and invasion of U87 cells, which was regulated through repression of PI3K/AKT pathway. Accordingly, the data from our work shed some light on EGFLAM might be a prognostic biomarker and therapeutic target for GBM.

Identification of nine genes as novel susceptibility loci for early-onset ischemic stroke, intracerebral hemorrhage, or subarachnoid hemorrhage

Given that substantial genetic components have been shown in ischemic stroke, intracerebral hemorrhage (ICH), and subarachnoid hemorrhage (SAH), heritability may be higher in early-onset than late-onset individuals with these conditions. Although genome-wide association studies (GWASs) have identified various genes and loci significantly associated with ischemic stroke, ICH, or intracranial aneurysm mainly in European ancestry populations, genetic variants that contribute to susceptibility to these disorders remain to be identified definitively. We performed exome-wide association studies (EWASs) to identify genetic variants that confer susceptibility to ischemic stroke, ICH, or SAH in early-onset subjects with these conditions. A total of 6,649 individuals aged ≤65 years were examined. For the EWAS of ischemic or hemorrhagic stroke, 6,224 individuals (450 subjects with ischemic stroke, 5,774 controls) or 6,179 individuals (261 subjects with ICH, 176 subjects with SAH, 5,742 controls), respectively, were examined. EWASs were performed with the use of Illumina Human Exome-12 v1.2 DNA Analysis BeadChip or Infinium Exome-24 v1.0 BeadChip. To compensate for multiple comparisons of allele frequencies with ischemic stroke, ICH, or SAH, we applied a false discovery rate (FDR) of <0.05 for statistical significance of association. The association of allele frequencies of 31,245 single nucleotide polymorphisms (SNPs) that passed quality control to ischemic stroke was examined with Fisher's exact test, and 31 SNPs were significantly (FDR <0.05) associated with ischemic stroke. The association of allele frequencies of 31,253 or 30,970 SNPs to ICH or SAH, respectively, was examined with Fisher's exact test, and six or two SNPs were significantly associated with ICH or SAH, respectively. Multivariable logistic regression analysis with adjustment for age, sex, and the prevalence of hypertension and diabetes mellitus revealed that 12 SNPs were significantly [P<0.0004 (0.05/124)] related to ischemic stroke. Similar analysis with adjustment for age, sex, and the prevalence of hypertension revealed that six or two SNPs were significantly [P<0.0016 (0.05/32)] related to ICH or SAH, respectively. After examination of linkage disequilibrium of identified SNPs and results of previous GWASs, we identified HHIPL2, CTNNA3, LOC643770, UTP20, and TRIB3 as susceptibility loci for ischemic stroke, DNTTIP2 and FAM205A as susceptibility loci for ICH, and FAM160A1 and OR52E4 as such loci for SAH. Therefore, to the best of our knowledge, we have newly identified nine genes that confer susceptibility to early-onset ischemic stroke, ICH, or SAH. Determination of genotypes for the SNPs in these genes may prove informative for assessment of the genetic risk for ischemic stroke, ICH, or SAH in Japanese.