Compound heterozygotes occur when different variants at the same locus on both maternal and paternal chromosomes produce a recessive trait. Here we present the tool VarCount for the quantification of variants at the individual level. We used VarCount to characterize compound heterozygous coding variants in patients with epileptic encephalopathy and in the 1000 Genomes Project participants. The Epi4k data contains variants identified by whole exome sequencing in patients with either Lennox-Gastaut Syndrome (LGS) or infantile spasms (IS), as well as their parents. We queried the Epi4k dataset (264 trios) and the phased 1000 Genomes Project data (2504 participants) for recessive variants. To assess enrichment, transcript counts were compared between the Epi4k and 1000 Genomes Project participants using minor allele frequency (MAF) cutoffs of 0.5 and 1.0%, and including all ancestries or only probands of European ancestry. In the Epi4k participants, we found enrichment for rare, compound heterozygous variants in six genes, including three involved in neuronal growth and development - PRTG (p = 0.00086, 1% MAF, combined ancestries), TNC (p = 0.022, 1% MAF, combined ancestries) and MACF1 (p = 0.0245, 0.5% MAF, EU ancestry). Due to the total number of transcripts considered in these analyses, the enrichment detected was not significant after correction for multiple testing and higher powered or prospective studies are necessary to validate the candidacy of these genes. However, PRTG, TNC and MACF1 are potential novel recessive epilepsy genes and our results highlight that compound heterozygous variants should be considered in sporadic epilepsy.

Noonan syndrome and related disorders are a group of clinically and genetically heterogeneous conditions caused by mutations in genes of the RAS/MAPK pathway. Noonan syndrome causes multiple congenital anomalies, which are frequently accompanied by hypertrophic cardiomyopathy (HCM). We report here a Tunisian patient with a severe phenotype of Noonan syndrome including neonatal HCM, facial dysmorphism, severe failure to thrive, cutaneous abnormalities, pectus excavatum and severe stunted growth, who died in her eighth month of life. Using whole exome sequencing, we identified a de novo mutation in exon 7 of the RAF1 gene: c.776C > A (p.Ser259Tyr). This mutation affects a highly conserved serine residue, a main mediator of Raf-1 inhibition via phosphorylation. To our knowledge the c.776C > A mutation has been previously reported in only one case with prenatally diagnosed Noonan syndrome. Our study further supports the striking correlation of RAF1 mutations with HCM and highlights the clinical severity of Noonan syndrome associated with a RAF1 p.Ser259Tyr mutation.

Chronic periodontitis (CP) is a common inflammatory illness affecting a large proportion of humans. Genetic factors are thought to play important roles in its onset and development. A functional polymorphism (rs1800587) in the promoter of the interleukin-1α gene (-889 C/T) has been found to confer risk of CP primarily in Europeans, but the association between this variant and CP in the Chinese population remains less conclusive. In the current study, we aimed to investigate the association between rs1800587 and CP in case-control samples of Han Chinese origin. A total of 1,777 study subjects, including 884 CP patients and 893 healthy controls, were collected. Genotyping of rs1800587 was performed using the SNAPSHOT method, and statistical analyses were conducted to evaluate the association between rs1800587 and CP. In our sample, rs1800587 was significantly associated with the onset of CP (additive model, T-allele vs. C-allele, p = 0.00359, odds ratio = 1.446, 95% confidence intervals (CIs) = 1.127-1.855; dominant model, (TT + TC) vs. CC, p = 0.00250, odds ratio = 1.502, 95% CIs = 1.152-1.957; overdominant model, TC vs. (TT + CC), p = 0.00264, odds ratio = 1.508, 95% CIs = 1.152-1.976). The T-allele and [TC] genotypes of rs1800587 were significantly overrepresented in CP patients compared with controls. Our data suggest that rs1800587 of IL-1α is significantly associated with the risk of CP in Han Chinese subjects, further confirming its possible involvement in the disease.

Primary aldosteronism is the most common form of secondary hypertension, and aldosteronoma makes up a significant proportion of primary aldosteronism cases. Aldosteronoma is also called aldosterone-producing adenoma (APA). Although there have been many studies about APA, the pathogenesis of this disease is not yet fully understood. In this study, we aimed to find out the difference of gene expression patterns between APA and nonfunctional adrenocortical adenoma (NFAA) using a weighted gene coexpression network (WGCNA) and differentially expressed gene (DEG) analysis; only the genes that meet the corresponding standards of both methods were defined as real hub genes and then used for further analysis. Twenty-nine real hub genes were found out, most of which were enriched in the phospholipid metabolic process. WISP2, S100A10, SSTR5-AS1, SLC29A1, APOC1, and SLITRK4 are six real hub genes with the same gene expression pattern between the combined and validation datasets, three of which indirectly or directly participate in lipid metabolism including WISP2, S100A10, and APOC1. According to the gene expression pattern of DEGs, we speculated five candidate drugs with potential therapeutic value for APA, one of which is cycloheximide, an inhibitor for phospholipid biosynthesis. All the evidence suggests that phospholipid metabolism may be an important pathophysiological mechanism for APA. Our study provides a new perspective regarding the pathophysiological mechanism of APA and offers some small molecules that may possibly be effective drugs against APA.

Pancreatic adenocarcinoma (PAAD) has been a huge challenge to public health due to its increasing incidence, frequent early metastasis, and poor outcome. The molecular basis of tumorigenesis and metastasis in PAAD is largely unclear. Here, we identified a novel tumor-suppressor long noncoding RNA (lncRNA) MBNL1-AS1, in PAAD and revealed its downstream mechanism. Quantitative real-time PCR (qRT-PCR) data showed that MBNL1-AS1 expression was significantly downregulated in PAAD tissues and cells, which was closely associated with metastasis and poor prognosis. Cell counting kit-8 (CCK-8) assay, transwell assay, and western blot verified that overexpression of MBNL1-AS1 suppressed cell proliferation, migration, and epithelial mesenchymal transformation (EMT) behavior in PAAD cells. By using a dual luciferase reporter gene system, we confirmed that miR-301b-3p was a direct target of MBNL1-AS1. Further mechanismic study revealed that upregulation of miR-301b-3p abolished the inhibitory effect of MBNL1-AS1 overexpression on cell proliferation, tumorigenesis, migration and EMT. Our results demonstrate that MBNL1-AS1 plays a tumor-suppressive role in PAAD mainly by downregulating miR-301b-3p, providing a novel therapeutic target for PAAD.

Because of the ability to metabolize a large number of electron acceptors such as nitrate, nitrite, fumarate, and metal oxides, Shewanella species have attracted much attention in recent years. Generally, the use of these electron acceptors is mainly achieved through electron transfer proteins and their interactions which will dynamically change across different environmental conditions in cells. Therefore, functional analysis of condition-specific molecular networks can reveal biological information on electron transfer processes. By integrating expression data and molecular networks, we constructed condition-specific molecular networks for Shewanella piezotolerans WP3. We then identified condition-specific key genes and studied their potential functions with an emphasis on their roles in electron transfer processes. Functional module analysis showed that different flagellar assembly modules appeared under these conditions and suggested that flagellar proteins are important for these conditions. We also identified the electron transfer modules underlying these various environmental conditions. The present results could help with screening electron transfer genes and understanding electron transfer processes under various environmental conditions for the Shewanella species.

This study aimed at investigating the effects of wrist-ankle acupuncture-moxibustion and hot compression with Chinese herbal medicine on pain symptoms, endothelin-1 (ET-1), and calcitonin gene-related peptide (CGRP) concentrations of patients with cervical radiculopathy (CR).

Colon cancer is the most common malignant tumor of the gastrointestinal tract, and approximately 80%-90% of colon cancers are colon adenocarcinomas (COADs). This study aimed to screen key microRNAs (miRNAs) associated with COAD. Differentially expressed (DE) miRNAs were screened between COAD and adjacent cancer samples based on the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas obtained from datasets. The miRNAs of interest were validated using quantitative real-time polymerase chain reaction. Moreover, the effects of hsa-miR-135b-5p on the biological behavior of COAD cells were observed. To obtain the target genes of hsa-miR-135b-5p, transcriptome sequencing of the SW480 cells was performed, followed by protein-protein interaction (PPI) network and hsa-miR-135b-5p-target gene regulatory network construction and prognostic analysis. Downregulation of hsa-miR-135b-5p significantly inhibited SW480 cell proliferation, migration, and invasion and significantly facilitated apoptosis (P < 0.05). A total of 3384 DEmRNAs were screened, and enrichment analysis showed that the upregulated mRNAs were enriched in 25 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 326 Gene Ontology Biological Processes (GO-BPs) while the downregulated mRNAs were enriched in 20 KEGG pathways and 276 GO-BPs. A PPI network was then constructed, and H2BC14, H2BC3, and H4C11 had a higher degree. In addition, a total of 352 hsa-miR-135b-5p-gene regulatory relationships were identified. Prognostic analysis showed that FOXN2, NSA2, MYCBP, DIRAS2, DESI1, and RAB33B had prognostic significance (P < 0.05). In addition, the validation analysis results showed that FOXN2, NSA2, and DESI1 were significantly expressed between the miR-135b-5p-inhibitor and negative control groups (P < 0.05). Therefore, downregulation of hsa-miR-135b-5p inhibits cell proliferation, migration, and invasion in COAD, and carcinogenesis may function by targeting FOXN2, NSA2, MYCBP, DIRAS2, DESI1, and RAB33B.

Assays of transposase accessible chromatin sequencing (ATAC-seq) is an efficient assay to investigate chromatin accessibility, which depends on the activity of a robust Tn5 transposase to fragment the genome while cutting in the sequencing adapters.

Sepsis is a severe disease with high mortality, and liver injury is an independent risk factor for sepsis morbidity and mortality. We analyzed co-differentially expressed genes (co-DEGs) to explore potential biomarkers and therapeutic targets for sepsis-related liver injury. Three gene expression datasets (GSE60088, GSE23767, and GSE71530) were downloaded from the Gene Expression Omnibus (GEO). DEGs were screened between sepsis and control samples using GEO2R. The association of these DEGs with infection and liver disease was analyzed by using the CTD database. GO functional analysis, KEGG pathway enrichment analysis, and protein-protein interaction (PPI) network analysis were performed to elucidate the potential molecular mechanism of DEGs. DEGs of different tissues in GSE60088 were analyzed again to obtain specific markers of septic liver injury. Mouse model of sepsis was also established by cecal ligation and puncture (CLP), and the expression of specific markers in liver, lung, and kidney tissues was analyzed using Western blot. Here, we identified 21 DEGs in three datasets with 8 hub genes, all of which showed higher inference scores in liver diseases than bacterial infections. Among them, only TNFRSF1A had a liver-specific differential expression. TNFRSF1A was also confirmed to be specifically reduced in septic liver tissues in mice. Therefore, TNFRSF1A may serve as a potential biomarker for septic liver injury.

Methylmalonic acidemia (MMA) is an autosomal recessive metabolic disorder mainly caused by mutations in the methylmalonyl coenzyme A mutase (MCM) gene (MMUT) and leads to the reduced activity of MCM. In this study, a 3-year-old girl was diagnosed with carnitine deficiency secondary to methylmalonic acidemia by tandem mass spectrometry (MS/MS) and gas chromatography/mass spectrometry (GS/MS). Whole-exome sequencing (WES) was performed on the patient and identified two compound heterozygous mutations in MMUT: c.554C>T (p. S185F) and c.729-730insTT (p. D244Lfs ∗ 39). Bioinformatics analysis predicted that the rare missense mutation of c.554C>T would be damaging. Moreover, this rare mutation resulted in the reduced levels of MMUT mRNA and MMUT protein. Collectively, our findings provide a greater understanding of the effects of MMUT variants and will facilitate the diagnosis and treatment of patients with MMA.

Glioblastoma (GBM) is a highly prevalent brain tumor characterized by high rates of morbidity, recurrence, and mortality. While temozolomide (TMZ) is commonly used as a first-line treatment for this cancer, the emergence of TMZ resistance limits its utility. The long noncoding RNA HOXA-AS2 reportedly drives GBM progression, but whether it can influence therapeutic resistance to TMZ has yet to be established.

Despite much attention given to the history of goat evolution in Kenya, information on the origin, demographic history, dispersal route, and genetic diversity of Galla goats remains unclear. Here, we examined the genetic background, diversity, demographic history, and population genetic variation of Galla goats using mtDNA D-loop and HSP70 single-nucleotide polymorphism markers. The results revealed 90 segregating sites and 68 haplotypes in a 600-bp mtDNA D-loop sequence. The overall mean mitochondrial haplotype diversity was 0.993. The haplotype diversities ranged between 0.8939 ± 0.0777 and 1.0000 ± 0.0221 in all populations supporting high genetic diversity. Mitochondrial phylogenetic analysis revealed three Galla goat haplogroups (A, G, and D), supporting multiple maternal ancestries, of which haplogroup A was the most predominant. Analysis of molecular variance (AMOVA) showed considerable variation within populations at 94.39%, evidence of high genetic diversity. Bimodal mismatch distribution patterns were observed while most populations recorded negative results for Tajima and Fu's Fs neutrality tests supporting population expansion. Genetic variation among populations was also confirmed using HSP70 gene fragment sequences, where six polymorphic sites which defined 21 haplotypes were discovered. Analysis of molecular variance revealed a significant FST index value of 0.134 and a high FIS index value of 0.746, an indication of inbreeding. This information will pave the way for conservation strategies and informed breeding to improve Galla or other goat breeds for climate-smart agriculture.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a worldwide pandemic, activates signaling cascades and leads to innate immune responses and secretion of multiple chemokines and cytokines. Long noncoding RNAs (lncRNAs) have a crucial role in inflammatory pathways. Through our search on the PubMed database, we discovered that existing research has primarily focused on examining the regulatory impacts of five lncRNAs in the context of viral infections. However, their role in regulating other conditions, including SARS-CoV-2, has not been explored. Therefore, this study aimed to investigate the expression pattern of lncRNAs in the peripheral blood mononuclear cells (PBMC) and their potential roles in SARS-CoV-2 infection. Potentially significant competing endogenous RNA (ceRNA) networks of these five lncRNAs were found using online in-silico techniques.

DNA methylation is an epigenetic marker that has been shown to vary significantly across different tissues. Taking advantage of the methylation differences between placenta-derived cell-free DNA and maternal blood, several groups employed different approaches for the discovery of fetal-specific biomarkers. The aim of this study was to analyse whole-genome fetal and maternal methylomes in order to identify and confirm the presence of differentially methylated regions (DMRs). We have initially utilized methylated DNA immunoprecipitation (MeDIP) and next-generation sequencing (NGS) to identify genome-wide DMRs between chorionic villus sampling (CVS) and female non-pregnant plasma (PL) and peripheral blood (WBF) samples. Next, using specific criteria, 331 fetal-specific DMRs were selected and confirmed in eight CVS, eight WBF and eight PL samples by combining MeDIP and in-solution targeted enrichment followed by NGS. Results showed higher enrichment in CVS samples as compared to both WBF and PL samples, confirming the distinct methylation levels between fetal and maternal DNA for the selected DMRs. We have successfully implemented a novel approach for the discovery and confirmation of a significant number of fetal-specific DMRs by combining for the first time MeDIP and in-solution targeted enrichment followed by NGS. The implementation of this double-enrichment approach is highly efficient and enables the detailed analysis of multiple DMRs by targeted NGS. Also, this is, to our knowledge, the first reported application of MeDIP on plasma samples, which leverages the implementation of our enrichment methodology in the detection of fetal abnormalities in maternal plasma.

Postmenopausal women experiences osteoporotic structural damage and bone fragility resulting from reduced bone formation and increased bone resorption. Osteoporosis frequently affects the vertebral column and causes compression fractures. This study aims to characterize roles of miRNAs in osteoporosis and subsequent incidence risk of vertebral fractures for postmenopausal women. Methods. Differentially expressed miRNAs between osteoporotic patients with vertebral fractures and osteoporotic patients without fracture were identified. This retrospective study included 78 osteoporotic patients with vertebral fractures and 82 osteoporotic patients without vertebral fractures. The plasma levels of bone metabolic markers, 25-hydroxyvitamin D (25-(OH)VitD), propeptide of type I procollagen (PINP), and β-Carboxyl terminal peptide (β-CTx), were detected using the patented electro-chemiluminescence (ECLIA) method. The expression levels of miR-491-5p and miR-485-3p were determined by qRT-PCR. Pearson correlation analysis was carried out to assess the relationship between miR-491-5p, miR-485-3p, and bone metabolic markers. Receiver operating characteristic (ROC) curves and the area under the ROC curve (AUC) were used to evaluate the performance of miR-491-5p and miR-485-3p in diagnosing the occurrence of vertebral fractures in osteoporotic patients.Results: The plasma levels of PINP and β-CTx were elevated but the plasma level of 25-(OH)VitD was declined in osteoporotic patients with vertebral fractures when comparable to those without (< 0.05). The plasma expression levels of miR-491-5p and miR-485-3p were declined osteoporotic patients with vertebral fractures when comparable to those without (< 0.001). Pearson correlation analysis revealed that the relative expression level of miR-491-5p was negatively correlated with the level of 25-(OH)VitD (r = -0.518, < 0.001) but positively correlated with the levels of PINP (r = 0.547, < 0.001) and β-CTx (r = 0.380, < 0.001). We also observed a negative correlation between the relative expression level of miR-485-3p and 25-(OH)VitD (r = -0.388, < 0.001), a positive correlation between miR-485-3p and PINP (r = 0.422,< 0.001). ROC curves for prediction of vertebral fracture following osteoporosis in postmenopausal women by miR-491-5p expression yielded 0.866 AUC and by miR-485-3p expression produced 0.848 AUC.

Psoriasis is an immune and inflammation-related skin disease. Triptolide with immunosuppressive and anti-inflammatory properties has been utilized for psoriasis treatment. However, the potential immunological mechanisms of triptolide have not been fully elucidated.

Osteopontin (OPN) is a phosphorylated glycoprotein. There is increasing evidence that the OPN gene played a major role in the progression of solid organ tumors. However, few studies have clarified how OPN regulated the functional role of human esophageal squamous cell carcinoma (ESCC). This study was designed to investigate the effect of OPN in esophageal squamous cell carcinoma.

Among primary liver carcinoma cases, the proportion of liver hepatocellular carcinoma (LIHC) cases is 75%-85%. Current treatments for LIHC include chemotherapy, surgical excision, and liver transplantation, which are effective for early LIHC treatment. Nevertheless, the early symptoms of liver carcinoma are atypical, so a large proportion of LIHC patients are diagnosed at an advanced stage. Histocompatibility minor 13 (HM13), located in the endoplasmic reticulum, is responsible for catalysing the hydrolysis of some signal peptides after cleavage from the precursor protein. Here, we studied the role of HM13 in LIHC development through bioinformatics analysis. Database analysis showed that HM13 was of great significance for LIHC tumorigenesis. Compared to normal liver tissues, HM13 expression was increased to a greater extent in LIHC tissues. After analysis of Kaplan‒Meier plotter and Gene Expression Profiling Interactive Analysis (GEPIA) datasets, we discovered that highly expressed HM13 exhibited an association with shorter overall survival (OS), disease-free survival (DFS), and disease-specific survival (DSS). We conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to analyse HM13-related genes, and the data indicated that these genes obviously participated in rRNA processing, ribosome biogenesis, spliceosome, Huntington's disease, and ATP-dependent helicase activity. The Cell Counting Kit-8 (CCK-8) assay and Transwell assay showed that reducing HM13 expression hindered LIHC cell proliferation, migration, and invasion. In conclusion, these findings indicate that HM13 is a biomarker and is related to the poor prognosis of LIHC. Our results are conducive to discovering new targets for LIHC treatment.

Object. β-Elemene is an emerging antitumor Chinese medicine, but the exact mechanism of action of β-elemene in colorectal cancer (CRC) remains unclear. This study aimed to explore the mechanism of the lncRNA-miRNA-mRNA network in the process of β-elemene inhibiting CRC. Methods. RNA sequencing was performed on CRC cells from the control group (untreated) and the case group (β-elemene-treated). According to the sequencing data, we screened the differentially expressed (DE) lncRNAs, miRNAs, and mRNAs and then analyzed them by functional enrichment analyses. Through the lncRNA-miRNA-mRNA network, the key miRNAs and mRNAs involved in the process of β-elemene inhibiting CRC were further identified. Results. Totally, 607 upregulated and 599 downregulated DElncRNAs, 12 downregulated and 24 upregulated DEmiRNAs, and 3153 downregulated and 3248 upregulated DEmRNAs were identified. Through the lncRNA-miRNA-mRNA network, 3 miRNAs (miR-7109-3p, miR-4506, and miR-3182), 7 prognostic mRNAs (ALPG, DTX1, HOXD13, RIMS3, SLC16A8, SYT1, and TNNT1), and 2 key mRNAs (RIMS3 and SLC16A8) were determined to participate in the inhibitory mechanism of β-elemene in CRC. Conclusion. This study revealed for the first time that the lncRNA-miRNA-mRNA network is involved in the regulation of β-elemene in CRC, and these identified miRNAs and mRNAs could be new clinical prognostic biomarkers and therapeutic targets for CRC patients.