Search

BackgroundEsophageal cancer (EC) is a prevalent malignancy characterized by a low 5-year survival rate, primarily attributed to delayed diagnosis and limited therapeutic options. Currently, early detection of EC heavily relies on endoscopy and pathological examination, which pose challenges due to their invasiveness and high costs, leading to low patient compliance. The detection of DNA methylation offers a non-endoscopic, cost-effective, and secure approach that holds promising prospects for early EC detection.MethodsTo identify improved methylation markers for early EC detection, we conducted a comprehensive review of relevant literature, summarized the performance of DNA methylation markers based on different input samples and analytical methods in EC early detection and screening.FindingsThis review reveals that blood cell free DNA methylation-based method is an effective non-invasive method for early detection of EC, although there is still a need to improve its sensitivity and specificity. Another highly sensitive and specific non-endoscopic approach for early detection of EC is the esophageal exfoliated cells based-DNA methylation analysis. However, while there are substantial studies in esophageal adenocarcinoma, further more validation is required in esophageal squamous cell carcinoma.ConclusionIn conclusion, DNA methylation detection holds significant potential as an early detection and screening technology for EC.

Background: Spontaneous abortion is the most common complication of early pregnancy. In this study, we aim to investigate the clinical application value of genetic diagnosis using single nucleotide polymorphism (SNP) microarray analysis on the products of conception and to characterize the types of genetic abnormalities and their prevalence in pregnancy loss in Northwest China.Methods: Over 48 months, we selected 652 products of conception, which included chorionic villi, fetal tissues, germ cell samples, amniotic fluid samples, cord blood samples, and a cardiac blood sample. We analyzed the distribution of chromosomal abnormalities leading to fetal arrest or abortion using SNP array. The patients were then categorized divided into groups based on maternal age, gestational age, number of miscarriages, and maternal ethnic background. The incidences of various chromosomal abnormalities in each group were compared.Results: Of the 652 cases, 314 (48.16%) exhibited chromosomal abnormalities. These included 286 cases with numerical chromosomal abnormalities, 24 cases with copy number variation, and four cases with loss of heterozygosity. Among them, there were 203 trisomy cases, 55 monosomy cases, and 28 polyploidy cases. In the subgroup analysis, significant differences were found in the frequency of numerical chromosomal abnormalities and copy number variation between the advanced and younger maternal age group as well as between the early and late abortion groups. Furthermore, we identified significant differences in the frequency of numerical chromosomal abnormalities between the first spontaneous abortion and recurrent miscarriage groups. However, there were no significant differences in the frequency of numerical chromosomal abnormalities between the Han and Uighur groups.Conclusion: Our research highlights chromosomal abnormalities as the primary cause of spontaneous abortion, with a higher incidence in early pregnancy and among women of advanced age. The use of SNP array analysis emerges as an effective and reliable technique for chromosome analysis in aborted fetuses. This method offers a comprehensive and dependable genetic investigation into the etiology of miscarriage, establishing itself as a valuable routine selection for genetic analysis in cases of natural abortions.

Background:PRDM12 is a newly discovered gene responsible for congenital insensitivity to pain (CIP). Its clinical manifestations are various and not widely known.Methods: The clinical data of two infants diagnosed with CIP associated with PRDM12 mutation were collected. A literature review was performed, and the clinical characteristics of 20 cases diagnosed with a mutation of PRDM12 were summarized and analyzed.Results: Two patients had pain insensitivity, tongue and lip defects, and corneal ulcers. The genomic analysis results showed that variants of PRDM12 were detected in the two families. The case 1 patient carried heterozygous variations of c.682+1G > A and c.502C > T (p.R168C), which were inherited from her father and mother, respectively. We enrolled 22 patients diagnosed with CIP through a literature review together with our cases. There were 16 male (72.7%) and 6 female (27.3%) patients. The age of onset ranged from 6 months to 57 years. The prevalence of clinic manifestation was 14 cases with insensitivity to pain (63.6%), 19 cases with self-mutilation behaviors (86.4%), 11 cases with tongue and lip defects (50%), 5 cases with mid-facial lesions (22.7%), 6 cases with distal phalanx injury (27.3%), 11 cases of recurrent infection (50%), 3 cases (13.6%) with anhidrosis, and 5 cases (22.7%) with global developmental delay. The prevalence of ocular symptoms was 11 cases (50%) with reduced tear secretion, 6 cases (27.3%) with decreased corneal sensitivity, 7 cases (31.8%) with disappeared corneal reflexes, 5.5 cases (25%, 0.5 indicated a single eye) with corneal opacity, 5 cases (22.7%) with corneal ulceration, and 1 case (4.5%) with a corneal scar.Conclusion: The syndrome caused by PRDM12 mutation is a clinically distinct and diagnosable disease that requires joint multidisciplinary management to control the development of the disease and minimize the occurrence of complications.

Premature ovarian insufficiency (POI) is characterized by early loss of ovarian function before the age of 40 years. It is confirmed to have a strong and indispensable genetic component. Caseinolytic mitochondrial matrix peptidase proteolytic subunit (CLPP) is a key inducer of mitochondrial protein quality control for the clearance of misfolded or damaged proteins, which is necessary to maintain mitochondrial function. Previous findings have shown that the variation in CLPP is closely related to the occurrence of POI, which is consistent with our findings. This study identified a novel CLPP missense variant (c.628G > A) in a woman with POI who presented with secondary amenorrhea, ovarian dysfunction, and primary infertility. The variant was located in exon 5 and resulted in a change from alanine to threonine (p.Ala210Thr). Importantly, Clpp was mainly localized in the cytoplasm of mouse ovarian granulosa cells and oocytes, and was relatively highly expressed in granulosa cells. Moreover, the overexpression of c.628G > A variant in human ovarian granulosa cells decreased the proliferative capacity. Functional experiments revealed that the inhibition of CLPP decreased the content and activity of oxidative respiratory chain complex IV by affecting the degradation of aggregated or misfolded COX5A, leading to the accumulation of reactive oxygen species and reduction of mitochondrial membrane potential, ultimately activating the intrinsic apoptotic pathways. The present study demonstrated that CLPP affected the apoptosis of granulosa cells, which might be one of the mechanisms by which CLPP aberrations led to the development of POI.

Genome-wide association study (GWAS) have identified over 1,000 loci associated with blood pressure. However, these loci only explain 6% of heritability. Transcriptome-wide association studies (TWAS) combine GWAS summary data with expression quantitative trait loci (eQTL) to provide a better approach to finding genes associated with complex traits. GWAS summary data (N = 450,584) for essential hypertension originating from European samples were subjected to Post-GWAS analysis using FUMA software and then combined with eQTL data from Genotype-Tissues Expression Project (GTEx) v8 for TWAS analysis using UTMOST, FUSION software, and then validated the results with SMR. FUMA identified 346 significant genes associated with hypertension, FUSION identified 461, and UTMOST cross-tissue analysis identified 34, of which 5 were common. SMR validation identified 3 key genes: ENPEP, USP38, and KCNK3. In previous GWAS studies on blood pressure regulation, the association of ENPEP and KCNK3 with hypertension has been established, and the association between USP38 and blood pressure regulation still needs further validation.

Introduction: Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer with low 5-year survival rate. Cellular senescence, characterized by permanent and irreversible cell proliferation arrest, plays an important role in tumorigenesis and development. This study aims to develop a cellular senescence-based stratified model, and a multivariable-based nomogram for guiding clinical therapy for HCC.Materials and methods: The mRNAs expression data of HCC patients and cellular senescence-related genes were obtained from TCGA and CellAge database, respectively. Through multiple analysis, a four cellular senescence-related genes-based prognostic stratified model was constructed and its predictive performance was validated through various methods. Then, a nomogram based on the model was constructed and HCC patients stratified by the model were analyzed for tumor mutation burden, tumor microenvironment, immune infiltration, drug sensitivity and immune checkpoint. Functional enrichment analysis was performed to explore potential biological pathways. Finally, we verified this model by siRNA transfection, scratch assay and Transwell Assay.Results: We established an cellular senescence-related genes-based stratified model, and a multivariable-based nomogram, which could accurately predict the prognosis of HCC patients in the ICGC database. The low and high risk score HCC patients stratified by the model showed different tumor mutation burden, tumor microenvironment, immune infiltration, drug sensitivity and immune checkpoint expressions. Functional enrichment analysis suggested several biological pathways related to the process and prognosis of HCC. Scratch assay and transwell assay indicated the promotion effects of the four cellular senescence-related genes (EZH2, G6PD, CBX8, and NDRG1) on the migraiton and invasion of HCC.Conclusion: We established a cellular senescence-based stratified model, and a multivariable-based nomogram, which could predict the survival of HCC patients and guide clinical treatment.

It is well known that histone modifications play an important part in various chromatin-dependent processes such as DNA replication, repair, and transcription. Using computational models to predict gene expression based on histone modifications has been intensively studied. However, the accuracy of the proposed models still has room for improvement, especially in cross-cell lines gene expression prediction. In the work, we proposed a new model TransferChrome to predict gene expression from histone modifications based on deep learning. The model uses a densely connected convolutional network to capture the features of histone modifications data and uses self-attention layers to aggregate global features of the data. For cross-cell lines gene expression prediction, TransferChrome adopts transfer learning to improve prediction accuracy. We trained and tested our model on 56 different cell lines from the REMC database. The experimental results show that our model achieved an average Area Under the Curve (AUC) score of 84.79%. Compared to three state-of-the-art models, TransferChrome improves the prediction performance on most cell lines. The experiments of cross-cell lines gene expression prediction show that TransferChrome performs best and is an efficient model for predicting cross-cell lines gene expression.

Background: Numerous studies have revealed that the long non-coding RNA LINC00662 is irregularly expressed in various cancers, as well as is correlated with cancer development and progression. Nevertheless, the clinical value of LINC00662 remains controversial. Hence, we explored the correlation of LINC00662 with cancer prognosis through meta-analysis and bioinformatics analysis.Methods: From the beginning through 12 March 2022, we searched for correlational studies on Web of Science, Embase, PubMed and The Cochrane Library. We used pooled hazard ratios (HRs) and odds ratios (ORs) with 95% confidence intervals (CIs) to determine the significance of studies on survival outcomes and clinicopathological aspects in human cancers. Additionally, the Gene Expression Profiling Interactive Analysis (GEPIA) database was employed to confirm our findings.Results: Our meta-analysis of 14 studies comprising a total of 960 cancer patients revealed that LINC00662 overexpression was correlated with poor overall survival (HR = 1.91, 95% CI 1.49–2.45, p < 0.001) in cancer patients and relapse-free survival (HR = 2.12, 95% CI 1.19–3.76, p = 0.010) in hepatocellular carcinoma patients. The correlation between LINC00662 and OS was further supported by the results of subgroup analyses according to cancer type, follow-up time, HR availability, and NOS score. In addition, LINC00662 overexpression predicted advanced tumor stage (OR = 4.23, 95% CI 2.50–7.17, p < 0.001), larger tumor size (OR = 1.49, 95% CI 1.11–1.99, p = 0.008), earlier lymph node metastasis (OR = 2.40, 95% CI 1.25–4.59, p = 0.008), and earlier distant metastasis (OR = 4.78, 95% CI 2.57–8.88, p < 0.001). However, there were no statistically significant differences in age (OR = 1.16, 95% CI 0.90–1.51, p = 0.246), gender (OR = 1.10, 95% CI 0.79–1.53, p = 0.578), or differentiation grade (OR = 1.53, 95% CI 0.71–3.33, p = 0.280).Conclusion: LINC00662 expression upregulation is associated with poor prognosis and advanced clinicopathological features in patients with multiple tumors. LINC00662 may serve as a biomarker for the diagnosis and treatment of patients with tumors.

High-fat diet could lead to a series of metabolic diseases, including obesity, and its mechanism is not clear. In this study, the rabbit individuals were fed with high-fat diet, the liver tissues were collected, high-throughput sequencing technology was used to reveal the expression of lncRNA and miRNA difference, and the molecular regulation mechanism of lncRNA-miRNA. A total of 24,615 DE lncRNAs and 52 DE miRNAs were identified, including 15 novel discovered DE miRNAs (5 upregulated and 10 downregulated). Furthermore, five miRNAs and three mRNAs were verified by qRT-PCR, and the results showed that the expression of the DE miRNAs and DE lncRNAs in the two groups was consistent with our sequencing results. GO and KEGG analyzed 7,57,139 target genes respectively, enriching the pathways related to lipid metabolism, including mucin O-glycan biosynthesis pathway, insulin resistance and glucagon signaling pathway. Moreover, 65 targeting relationships were obtained. Among them, LOC103348122/miR-450a-5p, LOC103350359/miR-450a-3p and LOC103350429/miR-148a-5p were proposed the first time. Significantly, LOC103348122/miR-450a-5p and LOC103350429/miR-148a-5p were related to lipid metabolism in the liver. This study is of great significance to the CeRNA regulatory network related to lipid metabolism in the liver of rabbits, and provides a basis for understanding hepatic steatosis in rabbits.

Nosema ceranae is a widespread fungal parasite for honey bees, causing bee nosemosis. Based on deep sequencing and bioinformatics, identification of circular RNAs (circRNAs) in Apis cerana workers’ midguts and circRNA-regulated immune response of host to N. ceranae invasion were conducted in this current work, followed by molecular verification of back-splicing sites and expression trends of circRNAs. Here, 10185 and 7405 circRNAs were identified in the midguts of workers at 7 days (AcT1) and 10 days (AcT2) post inoculation days post-inoculation with N. ceranae. PCR amplification result verified the back-splicing sites within three specific circRNAs (novel_circ_005123, novel_circ_007177, and novel_circ_015140) expressed in N. ceranae-inoculated midgut. In combination with transcriptome data from corresponding un-inoculated midguts (AcCK1 and AcCK2), 2266 circRNAs were found to be shared by the aforementioned four groups, whereas the numbers of specific ones were 2618, 1917, 5691, and 3723 respectively. Further, 83 52) differentially expressed circRNAs (DEcircRNAs) were identified in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group. Source genes of DEcircRNAs in workers’ midgut at seven dpi were involved in two cellular immune-related pathways such as endocytosis and ubiquitin mediated proteolysis. Additionally, competing endogenous RNA (ceRNA) network analysis showed that 23 13) DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group could target 18 14) miRNAs and further link to 1111 (1093) mRNAs. These target mRNAs were annotated to six cellular immunity pathways including endocytosis, lysosome, phagosome, ubiquitin mediated proteolysis, metabolism of xenobiotics by cytochrome P450, and insect hormone biosynthesis. Moreover, 284 164) internal ribosome entry site and 54 26) ORFs were identified from DEcircRNAs in AcCK1 vs. AcT1 (AcCK2 vs. AcT2) comparison group; additionally, ORFs in DEcircRNAs in midgut at seven dpi with N. ceranae were associated with several cellular immune pathways including endocytosis and ubiquitin-mediated proteolysis. Ultimately, RT-qPCR results showed that the expression trends of six DEcircRNAs were consistent with those in transcriptome data. These results demonstrated that N. ceranae altered the expression pattern of circRNAs in A. c. cerana workers’ midguts, and DEcircRNAs were likely to regulate host cellular and humoral immune response to microsporidian infection. Our findings lay a foundation for clarifying the mechanism underlying host immune response to N. ceranae infection and provide a new insight into interaction between Asian honey bee and microsporidian.

Background: As a new style of cell death, necroptosis plays a crucial role in tumor immune microenvironment. LncRNAs have been identified to act as competitive RNAs to influence genes involved in necroptosis. Therefore, we aim to create a signature based on necroptosis-related lncRNAs to predict the prognosis and immune landscape of lung adenocarcinoma (LUAD) patients in this study.Methods: TCGA database was used to acquire RNA sequencing (RNA-Seq) data and clinical information for 59 lung normal samples and 535 lung adenocarcinoma samples. The Pearson correlation analysis, univariate cox regression analysis and least absolute shrinkage and selection operator (LASSO) cox regression were performed to construct the prognostic NRlncRNAs signature. Then we used Kaplan-Meier (K-M) analysis, time-dependent ROC curves, univariate and multivariate cox regression analysis, and nomogram to validate this signature. In addition, GO, KEGG, and GSVA were analyzed to investigate the potential molecular mechanism. Moreover, we analyzed the relationship between our identified signature and immune microenvironment, TMB, and some clinical characteristics. Finally, we detected the expression of the six necroptosis-related lncRNAs in cells and tissues.Results: We constructed a NRlncRNAs signature consisting of six lncRNAs (FRMD6-AS1, LINC01480, FAM83A-AS1, FRMD6-AS1, MED4-AS1, and LINC01415) in LUAD. LUAD patients with high risk scores had lower chance of survival with an AUC of 0.739, 0.709, and 0.733 for 1-year, 3-year, and 5-year respectively. The results based on GO, KEGG, and GSVA enrichment analysis demonstrated that NRlncRNAs signature-related genes were mainly correlated with immune pathways, metabolic-and cell growth-related pathways, cell cycle, and apoptosis. Moreover, the risk score was correlated with the immune status of LUAD patients. Patients with higher risk scores had lower ESTIMATE scores and higher TIDE scores. The risk score was positively correlated with TMB. LINC01415, FRMD6-AS1 and FAM83A-AS1 were significantly overexpressed in lung adenocarcinoma, while the expression levels of MED4-AS1 and LINC01480 were lower in lung adenocarcinoma.Conclusion: Overall, an innovative prognostic signature based on NRlncRNAs was developed for LUAD through comprehensive bioinformatics analysis, which can act as a predictor of immunotherapy and may provide guidance for clinicians.

Background: Previous studies revealed that the gene signatures are associated with the modulation and pathogenesis of pulmonary arterial hypertension (PAH). However, identifying critical transcriptional signatures in the blood of PAH patients remains lacking.Methods: The differentially expressed transcriptional signatures in the blood of PAH patients were identified by a meta-analysis from four microarray datasets. Then we investigated the enrichment of gene ontology and KEGG pathways and identified top hub genes. Besides, we investigated the correlation of crucial hub genes with immune infiltrations, hallmark gene sets, and blood vessel remodeling genes. Furthermore, we investigated the diagnostic efficacy of essential hub genes and their expression validation in an independent cohort of PAH, and we validate the expression level of hub genes in monocrotaline (MCT) induced PAH rats’ model. Finally, we have identified the FDA-approved drugs that target the hub genes and their molecular docking.Results: We found 1,216 differentially expressed genes (DEGs), including 521 up-regulated and 695 down-regulated genes, in the blood of the PAH patients. The up-regulated DEGs are significantly associated with the enrichment of KEGG pathways mainly involved with immune regulation, cellular signaling, and metabolisms. We identified 13 master transcriptional regulators targeting the dysregulated genes in PAH. The STRING-based investigation identified the function of hub genes associated with multiple immune-related pathways in PAH. The expression levels of RPS27A, MAPK1, STAT1, RPS6, FBL, RPS3, RPS2, and GART are positively correlated with ssGSEA scores of various immune cells as positively correlated with the hallmark of oxidative stress. Besides, we found that these hub genes also regulate the vascular remodeling in PAH. Furthermore, the expression levels of identified hub genes showed good diagnostic efficacy in the blood of PAH, and we validated most of the hub genes are consistently dysregulated in an independent PAH cohort. Validation of hub genes expression level in the monocrotaline (MCT)-induced lung tissue of rats with PAH revealed that 5 screened hub genes (MAPK1, STAT1, TLR4, TLR2, GART) are significantly highly expressed in PAH rats, and 4 screened hub genes (RPS6, FBL, RPS3, and RPS2) are substantially lowly expressed in rats with PAH. Finally, we analyzed the interaction of hub proteins and FDA-approved drugs and revealed their molecular docking, and the results showed that MAPK1, TLR4, and GART interact with various drugs with appropriate binding affinity.Conclusion: The identified blood-derived key transcriptional signatures significantly correlate with immune infiltrations, hypoxia, glycolysis, and blood vessel remodeling genes. These findings may provide new insight into the diagnosis and treatment of PAH patients.

Background: Glioblastoma multiforme (GBM) is the most aggressive primary nervous system brain tumor. There is still a lack of effective methods to control its progression and recurrence in clinical treatment. It is clinically found that Xiaoliu Decoction (XLD) has the effect of treating brain tumors and preventing tumor recurrence. However, its mechanism is still unclear.Methods: Search the Traditional Chinese Medicine System Pharmacology Database (TCSMP) for efficient substances for the treatment of XLD in the treatment of GBM, and target the targeted genes of the effective ingredients to construct a network. At the same time, download GBM-related gene expression data from the TCGA and GTEX databases, screen differential expression bases, and establish a drug target disease network. Through bioinformatics analysis, the target genes and shared genes of the selected Chinese medicines are analyzed. Finally, molecular docking was performed to further clarify the possibility of XLD in multiple GBMs.Results: We screened 894 differentially expressed genes in GBM, 230 XLD active ingredients and 169 predicted targets of its active compounds, of which 19 target genes are related to the differential expression of GBM. Bioinformatics analysis shows that these targets are closely related to cell proliferation, cell cycle regulation, and DNA synthesis. Finally, through molecular docking, it was further confirmed that Tanshinone IIA, the active ingredient of XLD, was tightly bound to key proteins.Conclusion: To sum up, the results of this study suggest that the mechanism of XLD in the treatment of GBM involves multiple targets and signal pathways related to tumorigenesis and development. This study not only provides a new theoretical basis for the treatment of glioblastoma multiforme with traditional Chinese medicine, but also provides a new idea for the research and development of targeted drugs for the treatment of glioblastoma multiforme.

Background: Primary central nervous system lymphoma (PCNSL) is a rare B-cell lymphoma of central nervous system, which is often found in immunocompromised patients. The common clinical treatment of PCNSL is methotrexate (MTX) and whole brain radiation therapy. With the development of tumour immunology research, the tumour microenvironment of PCNSL is characterised by abnormal expression of different immune signature molecules and patients with PCNSL may benefit from tumour immunotherapy.Methods: In our research, RNA-seq data from 82 PCNSL patients were collated by mining the microarray data from the GEO database. All samples were classified into three types related to tumour immune response by the Cibersort algorithm and consistent clustering. Differential analysis of genes was used to uncover 2 sets of differential genes associated with tumour immunity. The ICI scores of each sample were obtained by PCA algorithm, and the relationship between ICI scores and immune checkpoint expression, immunotherapy and drug sensitivity was investigated. Genes associated with ICI scores and their functional characteristics were investigated by WGCNA analysis and PPI analysis, based on the ICI scores of each sample.Results: The tumour microenvironment in PCNSL has a greater relationship with the tumour immune response. ICI scores obtained from 375 differential genes were associated with multiple immune responses in PCNSL. PCNSL patients with higher ICI scores had a better tumour microenvironment and were sensitive to immunotherapy and some small molecule drug. This study also identified 64 genes associated with ICI scores, which may serve as important therapeutic and prognostic targets for PCNSL.Conclusion: The presence of multiple immunosuppressive responses in the tumour microenvironment of PCNSL which suggested that improving the immune function of PCNSL patients through immunotherapy and targeted therapies can be an effective treatment for PCNSL. And the ICI score and associated genes may also provide a better predictor of the clinical use of immunotherapy.

Aldehydes are primary volatile organic compounds (VOCs) in local Chinese chicken meat and contribute green grass, fatty, citrus, and bitter almond aromas to chicken meat. To understand the genetic basis of these aldehyde VOC aromas, we used approximately 500 Chinese Jingxing Yellow (JXY) chickens to conduct genome-wide association studies (GWAS) on the flavor traits with the data of single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs). In total, 501 association variants (253 SNPs and 248 INDELs) were found to be suggestively (SNPs: p-value < 2.77e-06 and INDELs: p-value < 3.78e-05) associated with total aldehydes (the sum of nine aldehydes), hexanal, heptanal, benzaldehyde, (E,E)-2,4-nonadienal, octanal, (E)-2-decenal, nonanal, decanal, and octadecanal. Of them, six SNPs and 23 INDELs reached a genome-wide significance level (SNPs: p-value < 1.38e-07 and INDELs: p-value < 1.89e-06). Potential candidate aldehyde genes were functionally annotated for lipid metabolism, especially fatty acid-related pathways and phospholipid-related gene ontology (GO) terms. Moreover, the GWAS analysis of total aldehydes, hexanal, and nonanal generated the most significant signals, and phenotypic content differed between different genotypes at candidate gene-related loci. For total aldehydes and hexanal traits, candidate genes were annotated based on the significant and suggestive variants on chromosomes 3 and 8 with highly polymorphic linkage blocks. The following candidate genes were also identified: GALM, MAP4K3, GPCPD1, RPS6KA2, CRLS1, ASAP1, TRMT6, SDC1, PUM2, ALDH9A1, MGST3, GMEB1, MECR, LDLRAP1, GPAM and ACSL5. We also found that polyunsaturated fatty acids (PUFAs) (C18:2n6c linoleic acid and C18:3n3 linolenic acid) were significantly correlated with total aldehydes and hexanal contents. PUFAs are important aldehyde precursors, and consistently, our results suggested that candidate genes involved in fatty acid pathways and phospholipid GO terms were identified in association loci. This work provides an understanding of the genetic basis of aldehyde formation, which is a key flavor-forming compound.

To protect the germplasm resources of Schizothorax biddulphi, we developed and used 20 pairs of polymorphic microsatellite primers to analyze the genetic diversity and structure of populations. A total of 126 samples were collected from the Qarqan River (CEC), Kizil River (KZL), and Aksu River (AKS) in Xinjiang, China. The results showed that 380 alleles were detected in 20 pairs of primers and the average number of alleles was 19.0. The effective allele numbers and Nei’s gene diversity ranged from 1.1499 to 1.1630 and 0.0962 to 0.1136, respectively. The Shannon index range suggested low levels of genetic diversity in all populations. The genetic distance between the CEC and AKS populations was the largest, and the genetic similarity was the smallest. There was a significant genetic differentiation between CEC and the other two populations. The UPGMA clustering tree was constructed based on population genetic distance, and the clustering tree constructed by individuals showed that the AKS population and KZL population were clustered together, and the CEC population was clustered separately. Also, the group structure analysis also got the same result. It can be seen that although the three populations of S. biddulphi do not have high genetic diversity, the differentiation between the populations was high and the gene flow was limited, especially the differentiation between the CEC population and the other two populations. This study not only provided genetic markers for the research of S. biddulphi but the results of this study also suggested the need for enhanced management of S. biddulphi populations.

Liver cancer is the sixth most frequently diagnosed primary malignancy and ranks as the third leading cause of cancer-related death worldwide in 2020. ER stress also plays a vital role in the pathogenesis of malignancies. In the current study, we aimed to construct an endoplasmic reticulum stress-related genes (ERGs) signature to predict the overall survival (OS) of patients with HCC. Differentially expressed ERGs (DE-ERGs) were analyzed using The Cancer Genome Atlas (TCGA-LIHC cohort) and International Cancer Genome Consortium (ICGC-LIRI-JP cohort) databases. The prognostic gene signature was identified by the univariate Cox regression and Least Absolute Shrinkage and Selection Operator (LASSO)-penalized Cox proportional hazards regression analysis. The predictive ability of the model was evaluated by utilizing Kaplan–Meier curves and time-dependent receiver operating characteristic (ROC) curves. Gene set variant analysis (GSVA) was performed to explore the underlying biological processes and signaling pathways. CIBERPORT and single-sample Gene Set Enrichment Analysis (ssGSEA) were implemented to estimate the immune status between the different risk groups. A total of 113 DE-ERGs were identified between 50 normal samples and 365 HCC samples in the TCGA-LIHC cohort, and 48 DE-ERGs were associated with OS through the univariate Cox regression. A six DE-ERGs (PPARGC1A, SQSTM1, SGK1, PON1, CDK1, and G6PD) signature was constructed and classified patients into high-risk and low-risk groups. The risk score was an independent prognostic indicator for OS (HR > 1, p < 0.001). The function enrichment analysis indicated that cell cycle, RNA degradation, protein localization, and cell division were the main biological processes. The high-risk group had higher immune cell infiltration levels than those of the low-risk group. We predicted the response to targeted therapy in high- and low-risk patients with HCC and found that the high-risk patients were more sensitive to pazopanib. At last, we verified the expression of the six gene patterns in HCC tissues by qRT-PCR and immunohistochemistry. This signature may be a potential tool to provide a choice for prognosis prediction and personal management of patients with HCC.

Background: Previous studies have verified that Inscuteable Spindle Orientation Adaptor Protein (INSC) can regulate cell proliferation and differentiation in the developing nervous system. It also plays an important role in spindle orientation during mitosis and asymmetric division of fibroblasts and participates in the process of stratification of the squamous epithelium. The role and potential mechanism of INSC in the development of colonic adenocarcinoma (COAD) have not been fully understood. This study aimed at exploring the prognostic value of INSC in COAD and the correlation of its expression with immune infiltration.Methods: The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx) project, Gene Expression Profiling Interactive Analysis (GEPIA), and Gene Expression Omnibus (GEO) database were used to analyze the expression of INSC in COAD. The INSC protein expression level was analyzed by immunohistochemistry staining and the Human Protein Atlas (HPA) database. The diagnostic and prognostic values of INSC in COAD patients were analyzed using receiver operating characteristic (ROC) and Kaplan–Meier (KM) survival curves. In order to understand whether INSC is an independent prognostic factor, we used univariable and multivariate Cox analyses to analyze INSC expression and several clinical characteristics with survival. We use STRING analysis to find INSC-related proteins and related biological events analyzed by Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. At last, GEPIA and the Tumor Immune Estimation Resource (TIMER) were employed to explore the relationship between INSC and immune infiltrates and its marker gene set.Results: INSC was lower expressed in COAD tissues than in normal colon tissues, which was correlated with tumor stage. Patients with lower expression of INSC had shorter overall survival (OS). Moreover, univariable Cox analysis demonstrated that high expression of INSC was an independent prognostic factor for COAD. ROC analysis showed INSC was an accurate marker for identifying tumors from normal colon tissue, and the AUC of the curve was 0.923. Significant GO term analysis by GSEA showed that genes correlated with INSC were found to be enriched in several immune-related pathways. Specifically, INSC expression showed significant negative correlations with infiltration levels of B cells, CD4+ T cells, macrophages, DCs, and their marker sets in COAD.Conclusion: INSC was provided with prognostic value in COAD and related to immune invasion.

Aims: Current idiopathic pulmonary fibrosis (IPF) therapies usually show a poor outcome or treatment efficacy. The search for new risk factors has significant implications in preventing, delaying, and treating IPF. The association between obesity and the risk of IPF is not clear. This study aimed to investigate the role of different obesity types in IPF risk, which provides the possibility of weight loss as a new approach for IPF prevention.Methods: We conducted a two-sample Mendelian randomization (MR) analysis to assess the causal effect of obesity on IPF risk. We collected summary data of genetically determined obesity-related traits, including body mass index (BMI), waist circumference (WC), and waist-to-hip ratio (WHR) from large-scale consortia (the sample size ranging from 232,101 to 681,275), and genetic association with IPF from one of the largest meta-analyses including 2,668 cases. A total of 35–469 single nucleotide polymorphisms were selected as instrumental variables for obesity-related traits. We further performed multivariable MR to estimate the independent effect of BMI and WC on the risk of IPF.Results: Increased BMI and WC were associated with higher risk of IPF [odds ratio (OR) = 1.51, 95% confidence interval (CI) (1.22–1.87), p = 1.27 × 10–4, and OR = 1.71, 95% CI (1.08–2.72), p = 2.33 × 10–2, respectively]. Similar results for the BMI and WC were obtained in the replicated analysis. Subsequently, only the result for BMI survived following the multiple testing correction and showed good consistency with the weighted median estimator. Sensitivity analyses indicated that there was no heterogeneity or horizontal pleiotropy for MR estimations. Further multivariable MR suggested that the BMI showed the same direction and similar magnitude with that in the univariable MR analysis. There was little evidence to support the causal role of WHR on the risk of IPF in this study.Conclusion: Genetically determined BMI demonstrates a causal risk for IPF, which offers a novel insight into probing potential mechanisms. Meanwhile, these results also suggest that weight loss may be beneficial to IPF prevention.

Background: Some studies have shown that the base excision repair (BER) pathway has an effect on HIV-1 replication. APEX1 and XRCC1 as key BER genes may affect DNA repair capacity. However, the roles of single nucleotide polymorphisms (SNPs) in APEX1 and XRCC1 and their impact on HIV-1 infection and AIDS progression remain unclear.Methods: A custom-designed 48-Plex SNPscan Kit was used for detection of single nucleotide polymorphisms. 601 HIV-1-infected men who have sex with men (MSM) and 624 age-matched healthy individuals were recruited in northern China. Four SNPs (rs1130409, rs1760944, rs2307486 and rs3136817) in APEX1 gene and three SNPs (rs1001581, rs25487 and rs25489) in XRCC1 gene were genotyped. The generalized multifactor dimension reduction (GMDR) method was used to identify the SNP-SNP interactions.Results: In this study, rs1130409 G allele, rs1001581 C allele and rs25487 C allele were associated with a higher risk of HIV-1 infection susceptibility (p = 0.020, p = 0.007 and p = 0.032, respectively). The frequencies of APEX1 haplotype TT and XRCC1 haplotype CT showed significant differences between cases and controls (p = 0.0372 and p = 0.0189, respectively). Interestingly, stratified analysis showed that the frequency of rs1001581 C allele was significantly higher in AIDS patients with the CD4+ T-lymphocyte count 200 cells/μl (p = 0.022). Moreover, significant gene-gene interactions among rs1130409, rs1001581 and rs25487 were identified by GMDR (p = 0.0107). Specially, individuals with five to six risk alleles have a higher susceptibility to HIV-1 infection than those with zero to two risk alleles (p < 0.001).Conclusion:APEX1 and XRCC1 gene polymorphisms were associated with the susceptibility to HIV-1 infection and AIDS progression in MSM populations in northern China.

Salt-alkali tolerance is one of the important breeding traits of Portunus trituberculatus. Identification of molecular markers linked to salt-alkali tolerance is prerequisite to develop such molecular marker-assisted breeding. In this study, Bulked Segregant Analysis (BSA) was used to screen molecular markers associated with salt-alkali tolerance trait in P. trituberculatus. Two DNA mixing pools with significant difference in salt-alkali tolerance were prepared and 94.83G of high-quality sequencing data was obtained. 855 SNPs and 1051 Indels were firstly selected as candidate markers by BSA analysis, out of which, 20 markers were further selected via △index value (close to 0 or 1) and eight of those were successfully verified. In addition, based on the located information of the markers in genome, eight candidate genes related to salt-alkali tolerance were anchored including ubiquitin-conjugating enzyme, aspartate–tRNA ligase, vesicle-trafficking protein, and so on. qPCR results showed that the expression patterns of all these genes changed significantly after salt-alkali stress, suggesting that they play certain roles in salt-alkali adaptation. Our results will provide applicable markers for molecular marker-assisted breeding and help to clarify the mechanisms of salt-alkali adaptation of P. trituberculatus.

Background: Mainstream application of cancer immunotherapy is hampered by the low response rate of most cancer patients. A novel immunotherapeutic target or a biomarker predicting response to immunotherapy needs to be developed. Guanylate-binding protein 1 (GBP1) is an interferon (IFN)-inducible guanosine triphosphatases (GTPases) involving inflammation and infection. However, the immunological effects of GBP1 in pan-cancer patients are still obscure.Methods: Using large-scale public data, we delineated the landscape of GBP1 across 33 cancer types. The correlation between GBP1 expression or mutation and immune cell infiltration was estimated by ESTIMATE, TIMER, xCell, and quanTIseq algorithms. GBP1-related genes and proteins were subjected to function enrichment analysis. Clustering analysis explored the relationship between GBP1 expression and anti-tumor immune phenotypes. We assessed the patient’s response to immunotherapy using the tumor immune dysfunction and exclusion (TIDE) score and immunophenoscore (IPS). Furthermore, we validated the predictive power of GBP1 expression in four independent immunotherapy cohorts.Results: GBP1 was differentially expressed in tumors and normal tissues in multiple cancer types. Distinct correlations existed between GBP1 expression and prognosis in cancer patients. GBP1 expression and mutation were positively associated with immune cell infiltration. Function enrichment analysis showed that GBP1-related genes were enriched in immune-related pathways. Positive correlations were also observed between GBP1 expression and the expression of immune checkpoints, as well as tumor mutation burden (TMB). Pan-cancer patients with higher GBP1 expression were more inclined to display “hot” anti-tumor immune phenotypes and had lower TIDE scores and higher immunophenoscore, suggesting that these patients had better responses to immunotherapy. Patients with higher GBP1 expression exhibited improved overall survival and clinical benefits in immunotherapy cohorts, including the Gide et al. cohort [area under the curve (AUC): 0.813], the IMvigor210 cohort (AUC: 0.607), the Lauss et al. cohort (AUC: 0.740), and the Kim et al. cohort (AUC: 0.793).Conclusion: This study provides comprehensive insights into the role of GBP1 in a pan-cancer manner. We identify GBP1 expression as a predictive biomarker for immunotherapy, potentially enabling more precise and personalized immunotherapeutic strategies in the future.

Coronary heart disease (CHD) is a global health concern with high morbidity and mortality rates. This study aimed to identify the possible long non-coding RNA (lncRNA) biomarkers of CHD. The lncRNA- and mRNA-related data of patients with CHD were downloaded from the Gene Expression Omnibus database (GSE113079). The limma package was used to identify differentially expressed lncRNAs and mRNAs (DElncRNAs and DEmRNAs, respectively). Then, miRcode, TargetScan, miRDB, and miRTarBase databases were used to form the competing endogenous RNA (ceRNA) network. Furthermore, SPSS Modeler 18.0 was used to construct a logistic stepwise regression prediction model for CHD diagnosis based on DElncRNAs. Of the microarray data, 70% was used as a training set and 30% as a test set. Moreover, a validation cohort including 30 patients with CHD and 30 healthy controls was used to verify the hub lncRNA expression through real-time reverse transcription-quantitative PCR (RT-qPCR). A total of 185 DElncRNAs (114 upregulated and 71 downregulated) and 382 DEmRNAs (162 upregulated and 220 downregulated) between CHD and healthy controls were identified from the microarray data. Furthermore, through bioinformatics prediction, a 38 lncRNA-21miRNA-40 mRNA ceRNA network was constructed. Next, by constructing a logistic stepwise regression prediction model for 38 DElncRNAs, we screened two hub lncRNAs AC010082.1 and AC011443.1 (p < 0.05). The sensitivity, specificity, and area under the curve were 98.41%, 100%, and 0.995, respectively, for the training set and 93.33%, 91.67%, and 0.983, respectively, for the test set. We further verified the significant upregulation of AC010082.1 (p < 0.01) and AC011443.1 (p < 0.05) in patients with CHD using RT-qPCR in the validation cohort. Our results suggest that lncRNA AC010082.1 and AC011443.1 are potential biomarkers of CHD. Their pathological mechanism in CHD requires further validation.

Hearing loss (HL) is one of the most common disabilities in the world. In industrialized countries, HL occurs in 1–2/1,000 newborns, and approximately 60% of HL is caused by genetic factors. Next generation sequencing (NGS) has been widely used to identify many candidate genes and variants in patients with HL, but the data are scattered in multitudinous studies. It is a challenge for scientists, clinicians, and biologists to easily obtain and analyze HL genes and variant data from these studies. Thus, we developed a one-stop database of HL-related genes and variants, Gene4HL (http://www.genemed.tech/gene4hl/), making it easy to catalog, search, browse and analyze the genetic data. Gene4HL integrates the detailed genetic and clinical data of 326 HL-related genes from 1,608 published studies, along with 62 popular genetic data sources to provide comprehensive knowledge of candidate genes and variants associated with HL. Additionally, Gene4HL supports the users to analyze their own genetic engineering network data, performs comprehensive annotation, and prioritizes candidate genes and variations using custom parameters. Thus, Gene4HL can help users explain the function of HL genes and the clinical significance of variants by correlating the genotypes and phenotypes in humans.

Background: Acute Myeloid Leukemia (AML) is a complex and heterogeneous hematologic malignancy. However, the function of prognosis-related signature genes in AML remains unclear.Methods: In the current study, transcriptome sequencing was performed on 15 clinical samples, differentially expressed RNAs were identified using R software. The potential interactions network was constructed by using the common genes between target genes of differentially expressed miRNAs with transcriptome sequencing results. Functional and pathway enrichment analysis was performed to identify candidate gene-mediated aberrant signaling pathways. Hub genes were identified by the cytohubba plugin in Cytoscape software, which then expanded the potential interactions regulatory module for hub genes. TCGA-LAML clinical data were used for the prognostic analysis of the hub genes in the regulatory network, and GVSA analysis was used to identify the immune signature of prognosis-related hub genes. qRT-PCR was used to verify the expression of hub genes in independent clinical samples.Results: We obtained 1,610 differentially expressed lncRNAs, 233 differentially expressed miRNAs, and 2,217 differentially expressed mRNAs from transcriptome sequencing. The potential interactions network is constructed by 12 lncRNAs, 25 miRNAs, and 692 mRNAs. Subsequently, a sub-network including 15 miRNAs as well as 12 lncRNAs was created based on the expanded regulatory modules of 25 key genes. The prognostic analysis results show that CCL5 and lncRNA UCA1 was a significant impact on the prognosis of AML. Besides, we found three potential interactions networks such as lncRNA UCA1/hsa-miR-16-5p/COL4A5, lncRNA UCA1/hsa-miR-16-5p/SPARC, and lncRNA SNORA27/hsa-miR-17-5p/CCL5 may play an important role in AML. Furthermore, the evaluation of the immune infiltration shows that CCL5 is positively correlated with various immune signatures, and lncRNA UCA1 is negatively correlated with the immune signatures. Finally, the result of qRT-PCR showed that CCL5 is down-regulated and lncRNA UCA1 is up-regulated in AML samples separately.Conclusions: In conclusion, we propose that CCL5 and lncRNA UCA1 could be recognized biomarkers for predicting survival prognosis based on constructing competing endogenous RNAs in AML, which will provide us novel insight into developing novel prognostic, diagnostic, and therapeutic for AML.

Background: The liver is the central metabolic organ of animals. In chicken, knowledge on the relationship between gene expression in the liver and fat deposition during development is still limited. A time-course transcriptomic study from the embryonic (day 12) to the egg-producing period (day 180 after hatch) was performed to profile slow-growing meat type chicken liver gene expression and to investigate its correlation with abdominal fat deposition.Results: The transcriptome profiles showed a separation of the different developmental stages. In total, 13,096 genes were ubiquitously expressed at all the tested developmental stages. The analysis of differentially expressed genes between adjacent developmental stages showed that biosynthesis of unsaturated fatty acids pathway was enriched from day 21 to day 140 after hatch. The correlation between liver gene expression and the trait abdominal fat weight (AFW) was analyzed by weighted gene co-expression network analysis. The genes MFGE8, HHLA1, CKAP2, and ACSBG2 were identified as hub genes in AFW positively correlated modules, which suggested important roles of these genes in the lipid metabolism in chicken liver.Conclusion: Our results provided a resource of developmental transcriptome profiles in chicken liver and suggested that the gene ACSBG2 among other detected genes can be used as a candidate gene for selecting low AFW chickens.

BackgroundPrevious studies have highlighted the important role of vitamin D and calcium pathway genes in immune modulation, cell differentiation and proliferation, and inflammation regulation, all closely implicated in the pathogenesis of non-alcoholic fatty liver disease (NAFLD).ObjectiveThis study aims to investigate whether 11 candidate single nucleotide polymorphisms (SNPs) in vitamin D and calcium pathway genes (CYP2R1, CYP24A1, and CYP27B1) are associated with the risk of NAFLD.MethodsIn this case-control study, a total of 3,023 subjects were enrolled, including 1,114 NAFLD cases and 1,909 controls. Eleven genetic variants in CYP2R1, CYP24A1, and CYP27B1 genes were genotyped. Logistic regression analysis was used to assess the effects of these variants on NAFLD risk. The functional annotations of positive SNPs were further evaluated by bioinformatics analysis.ResultsAfter adjusting for age, gender, and metabolic measures, we identified that CYP24A1 rs2296241 variant genotypes (recessive model: OR, 1.316; 95% CI, 1.048–1.653; p = 0.018), rs2248359 variant genotypes (recessive model: OR, 1.315; 95% CI, 1.033–1.674; p = 0.026), and CYP27B1 rs4646536 variant genotypes (additive model: OR, 1.147; 95% CI, 1.005–1.310; p = 0.042) were associated with an elevated risk of NAFLD. In combined effects analysis, we found that NAFLD risk significantly increased among patients carrying more rs2296241-A, rs2248359-T, and rs4646536-T alleles (ptrend = 0.049). Multivariate stepwise analysis indicated that age, visceral obesity, ALT, γ-GT, hypertriglyceridemia, hypertension, low HDL-C, hyperglycemia, and unfavorable alleles were independent predictors of NAFLD (all p < 0.05). The area under the receiver operating characteristic curve was 0.789 for all the above factors.ConclusionThe polymorphisms of vitamin family genes CYP24A1 (rs2296241, CYP24A1, and rs2248359) and CYP27B1 (rs4646536) were associated with NAFLD risk in Chinese Han population, which might provide new insight into NAFLD pathogenesis and tools for screening high-risk population.

The immune response within the tumor microenvironment plays a key role in tumorigenesis and determines the clinical outcomes of head and neck squamous cell carcinoma (HNSCC). However, to date, very limited robust and reliable immunological biomarkers have been developed that are capable of estimating prognosis in HNSCC patients. In this study, we aimed to identify the effects of novel immune-related gene signatures (IRGs) that can predict HNSCC prognosis. Based on gene expression profiles and clinical data of HNSCC patient cohorts from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, a total of 439 highly variable expressed immune-related genes (including 239 upregulated and 200 downregulated genes) were identified by using differential gene expression analysis. Pathway enrichment analysis indicated that these immune-related differentially expressed genes were enriched in inflammatory functions. After process screening in the training TCGA cohort, six immune-related genes (PLAU, STC2, TNFRSF4, PDGFA, DKK1, and CHGB) were significantly associated with overall survival (OS) based on the LASSO Cox regression model. Integrating these genes with clinicopathological features, a multivariable model was built and suggested better performance in determining patients’ OS in the testing cohort, and the independent validation cohort. In conclusion, a well-established model encompassing both immune-related gene signatures and clinicopathological factors would serve as a promising tool for the prognostic prediction of HNSCC.

Circular RNA (circRNA) is an important factor for regulating the progression of many cardiovascular diseases, including acute myocardial infarction (AMI). However, the role of circ_0124644 in AMI progression remains unclear. Hypoxia was used to induce cardiomyocytes injury. The expression of circ_0124644, microRNA (miR)-590-3p, and SRY-box transcription factor 4 (SOX4) mRNA was measured by qRT-PCR. Cell counting kit 8 (CCK8) assay and flow cytometry were utilized to detect cell viability, cell cycle progression, and apoptosis. The protein levels of apoptosis markers and SOX4 were determined by western blot (WB) analysis, and the levels of oxidative stress markers were assessed using commercial Assay Kits. Dual-luciferase reporter assay, RIP assay, and RNA pull-down assay were employed to confirm the interaction between miR-590-3p and circ_0124644 or SOX4. Circ_0124644 was upregulated in AMI patients and hypoxia-induced cardiomyocytes. Hypoxia could inhibit cardiomyocytes viability, cell cycle process, and promote apoptosis and oxidative stress, while silencing circ_0124644 could alleviate hypoxia-induced cardiomyocytes injury. In terms of mechanism, circ_0124644 could target miR-590-3p. MiR-590-3p overexpression could relieve hypoxia-induced cardiomyocytes injury. Also, the suppressive effect of circ_0124644 knockdown on hypoxia-induced cardiomyocytes injury could be reversed by miR-590-3p inhibitor. Moreover, SOX4 was found to be a target of miR-590-3p, and its overexpression also could reverse the regulation of miR-590-3p on hypoxia-induced cardiomyocytes injury. Circ_0124644 silencing could alleviate hypoxia-induced cardiomyocytes injury by regulating the miR-590-3p/SOX4 axis, suggesting that it might be a target for alleviating AMI.

Genotypic data provide deep insights into the population history and medical genetics. The local ancestry inference (LAI) (also termed local ancestry deconvolution) method uses the hidden Markov model (HMM) to solve the mathematical problem of ancestry reconstruction based on genomic data. HMM is combined with other statistical models and machine learning techniques for particular genetic tasks in a series of computer tools. In this article, we surveyed the mathematical structure, application characteristics, historical development, and benchmark analysis of the LAI method in detail, which will help researchers better understand and further develop LAI methods. Firstly, we extensively explore the mathematical structure of each model and its characteristic applications. Next, we use bibliometrics to show detailed model application fields and list articles to elaborate on the historical development. LAI publications had experienced a peak period during 2006–2016 and had kept on moving in the following years. The efficiency, accuracy, and stability of the existing models were evaluated by the benchmark. We find that phased data had higher accuracy in comparison with unphased data. We summarize these models with their distinct advantages and disadvantages. The Loter model uses dynamic programming to obtain a globally optimal solution with its parameter-free advantage. Aligned bases can be used directly in the Seqmix model if the genotype is hard to call. This research may help model developers to realize current challenges, develop more advanced models, and enable scholars to select appropriate models according to given populations and datasets.

YWHAG, which encodes an adapter protein 14-3-3γ, is highly expressed in the brain and regulates a diverse range of cell signaling pathways. Previously, eight YWHAG mutations have been identified in patients with epileptic encephalopathy (EE). In this study, using trios-based whole exome sequencing, we identified two novel YWHAG mutations in two unrelated families with childhood myoclonic epilepsy and/or febrile seizures (FS). The identified mutations included a heterozygous truncating mutation (c.124C>T/p.Arg42Ter) and a de novo missense mutation (c.373A>G/p.Lys125Glu). The two probands experienced daily myoclonic seizures that were recorded with ictal generalized polyspike-slow waves, but became seizure-free with simple valproate treatment. The other affected individuals presented FS. The truncating mutation was identified in the family with six individuals of mild phenotype, suggesting that YWHAG mutations of haploinsufficiency are relatively less pathogenic. Analysis on all missense mutations showed that nine mutations were located within 14-3-3γ binding groove and another mutation was located at residues critical for dimerization, indicating a molecular sub-regional effect. Mutation Arg132Cys, which was identified recurrently in five patients with EE, would have the strongest influence on binding affinity. 14-3-3γ dimers supports target proteins activity. Thus, a heterozygous missense mutation would lead to majority dimers being mutants; whereas a heterozygous truncating mutation would lead to only decreasing the number of wild-type dimer, being one of the explanations for phenotypical variation. This study suggests that YWHAG is potentially a candidate pathogenic gene of childhood myoclonic epilepsy and FS. The spectrum of epilepsy caused by YWHAG mutations potentially range from mild myoclonic epilepsy and FS to severe EE.

BackgroundCarcass traits are crucial characteristics of broilers. However, the underlying genetic mechanisms are not well understood. In the current study, significant loci and major-effect candidate genes affecting nine carcass traits related to meat production were analyzed in 873 purebred broilers using an imputation-based genome-wide association study.ResultsThe heritability estimates of nine carcass traits, including carcass weight, thigh muscle weight, and thigh muscle percentage, were moderate to high and ranged from 0.21 to 0.39. Twelve genome-wide significant SNPs and 118 suggestively significant SNPs of 546,656 autosomal variants were associated with carcass traits. All SNPs for six weight traits (body weight at 42 days of age, carcass weight, eviscerated weight, whole thigh weight, thigh weight, and thigh muscle weight) were clustered around the 24.08 Kb region (GGA24: 5.73–5.75 Mb) and contained only one candidate gene (DRD2). The most significant SNP, rs15226023, accounted for 4.85–7.71% of the estimated genetic variance of the six weight traits. The remaining SNPs for carcass composition traits (whole thigh percentage and thigh percentage) were clustered around the 42.52 Kb region (GGA3: 53.03–53.08 Mb) and contained only one candidate gene (ADGRG6). The most significant SNP in this region, rs13571431, accounted for 11.89–13.56% of the estimated genetic variance of two carcass composition traits. Some degree of genetic differentiation in ADGRG6 between large and small breeds was observed.ConclusionWe identified one 24.08 Kb region for weight traits and one 42.52 Kb region for thigh-related carcass traits. DRD2 was the major-effect candidate gene for weight traits, and ADGRG6 was the major-effect candidate gene for carcass composition traits. Our results supply essential information for causative mutation identification of carcass traits in broilers.

For lung adenocarcinoma (LUAD), patients of different stages have strong heterogeneity, and their overall prognosis varies greatly. Thus, exploration of novel biomarkers to better clarify the characteristics of LUAD is urgent. Multi-omics information of LUAD patients were collected form TCGA. Three independent LUAD cohorts were obtained from gene expression omnibus (GEO). A multi-omics correlation analysis of METTL5 was performed in TCGA dataset. To build a METTL5-associated prognostic score (MAPS). Spathial and random forest methods were first applied for feature selection. Then, LASSO was implemented to develop the model in TCGA cohort. The prognostic value of MAPS was validated in three independent GEO datasets. Finally, functional annotation was conducted using gene set enrichment analysis (GSEA) and the abundances of infiltrated immune cells were estimated by ImmuCellAI algorithm. A total of 901 LUAD patients were included. The expression of METTL5 in LUAD was significantly higher than that in normal lung tissue. And high expression of METTL5 indicated poor prognosis in all different stages (P < 0.001, HR = 1.81). Five genes (RAC1, C11of24, METTL5, RCCD1, and SLC7A5) were used to construct MAPS and MAPS was significantly correlated with poor prognosis (P < 0.001, HR = 2.15). Furthermore, multivariate Cox regression analysis suggested MAPS as an independent prognostic factor. Functional enrichment revealed significant association between MAPS and several immune components and pathways. This study provides insights into the potential significance of METTL5 in LUAD and MAPS can serve as a promising biomarker for LUAD.

Ketosis is a common metabolic disease in dairy cows during early lactation. However, information about the metabolomic and proteomic profiles associated with the incidence and progression of ketosis is still limited. In this study, an integrated metabolomics and proteomics approach was performed on blood serum sampled from cows diagnosed with clinical ketosis (case, ≥ 2.60 mmol/L plasma β-hydroxybutyrate; BHBA) and healthy controls (control, < 1.0 mmol/L BHBA). Samples were taken 2 weeks before parturition and 2 weeks after parturition from 19 animals (nine cases, 10 controls). All serum samples (n = 38) were subjected to Liquid Chromatography-Mass Spectrometry (LC-MS) based metabolomic analysis, and 20 samples underwent Data-Independent Acquisition (DIA) LC-MS based proteomic analysis. A total of 97 metabolites and 540 proteins were successfully identified, and multivariate analysis revealed significant differences in both metabolomic and proteomic profiles between cases and controls. We investigated clinical ketosis-associated metabolomic and proteomic changes using statistical analyses. Correlation analysis of statistically significant metabolites and proteins showed 78 strong correlations (correlation coefficient, R ≥ 0.7) between 38 metabolites and 25 proteins, which were then mapped to pathways using IMPaLA. Results showed that ketosis altered a wide range of metabolic pathways, such as metabolism, metabolism of proteins, gene expression and post-translational protein modification, vitamin metabolism, signaling, and disease related pathways. Findings presented here are relevant for identifying molecular targets for ketosis and biomarkers for ketosis detection during the transition period.

IntroductionIdiopathic focal epilepsy (IFE) is a group of self-limited epilepsies. The etiology for the majority of the patients with IFE remains elusive. We thus screened disease-causing variants in the patients with IFE.MethodsWhole-exome sequencing was performed in a cohort of 323 patients with IFE. Protein modeling was performed to predict the effects of missense variants. The genotype–phenotype correlation of the newly defined causative gene was analyzed.ResultsFour novel heterozygous variants in PGM3, including two de novo variants, were identified in four unrelated individuals with IFE. The variants included one truncating variant (c.1432C > T/p.Q478X) and three missense variants (c.478C > T/p.P160S, c.1239C > G/p.N413K, and c.1659T > A/p.N553K), which had no allele frequency in the gnomAD database. The missense variants were predicted to be damaging and affect hydrogen bonds with surrounding amino acids. Mutations Q478X, P160S, and N413K were associated with benign childhood epilepsy with centrotemporal electroencephalograph (EEG) spikes. P160S and N413K were located in the inner side of the enzyme active center. Mutation N553K was associated with benign occipital epilepsy with incomplete penetrance, located in the C-terminal of Domain 4. Further analysis demonstrated that previously reported biallelic PGM3 mutations were associated with severe immunodeficiency and/or congenital disorder of glycosylation, commonly accompanied by neurodevelopmental abnormalities, while monoallelic mutations were associated with milder symptoms like IFE.ConclusionThe genetic and molecular evidence from the present study implies that the PGM3 variants identified in IFE patients lead to defects of the PGM3 gene, suggesting that the PGM3 gene is potentially associated with epilepsy. The genotype–phenotype relationship of PGM3 mutations suggested a quantitative correlation between genetic impairment and phenotypic severity, which helps explain the mild symptoms and incomplete penetrance in individuals with IFE.

Mung bean (Vigna radiata L.) is an important but understudied food legume in Asia and now worldwide. Genetic studies may help to accelerate the exploitation of new genes for breeding in this crop. Here, we used a recombination inbred line population to construct an SNP genetic linkage map by genome sequencing technology. We obtained 21,508 high-quality SNP markers integrated into 1,946 bin markers that were mapped onto 11 linkage groups (LGs) with 99–258 bin markers per LG. The total genetic length of the map was 1060.2 cM (38.76–168.03 cM per LG), with an average distance between markers of 0.54 cM. However, there were 18 gaps >5 cM, distribution on LG1, 3, 5, 7, and 9. Gene mapping for lobed and indented leaflets was conducted using the map. A major quantitative trait locus (QTL) associated with indented leaflets was detected on chromosome 10, with phenotypic variation explained (PVE) values of 39.7% and 45.4% under two different environments. Several QTLs for lobed leaflets were detected and most of them were tightly linked together on Chromosome 3. However, only one major QTL, which explained the largest phenotypic variation (27.7–69.5%), was stably detected under two different environments using both R and Q methods. In the two main stable QTLs regions on chromosomes 3 and 10, candidate genes for regulating the molecular mechanism of different leaflet shapes were detected by functional annotation. The overlap of major QTLs under different environments indicated that the present map would be good enough for precisely mapping genes, and both the QTL analysis and gene prediction were useful for investigating the mechanism of leaf development in mung bean or legumes.

Fat traits are important in the chicken industry where there is a desire for high intramuscular fat (IMF) and low abdominal fat. However, there is limited knowledge on the relationship between the dynamic status of gene expression and the body fat deposition in chicken. Transcriptome data were obtained from breast muscle and abdominal fat of female chickens from nine developmental stages (from embryonic day 12 to hatched day 180). In total, 8,545 genes in breast muscle and 6,824 genes in abdominal fat were identified as developmentally dynamic genes. Weighted correlation network analysis was used to identify gene modules and the hub genes. Twenty-one hub genes were identified, e.g., ENSGALG00000041996, which represents a candidate for high IMF, and CREB3L1, which relates to low abdominal fat weight. The transcript factor L3MBTL1 and the transcript factor cofactors TNIP1, HAT1, and BEND6 related to both high breast muscle IMF and low abdominal fat weight. Our results provide a resource of developmental transcriptome profiles in chicken breast muscle and abdominal fat. The candidate genes can be used in the selection for increased IMF content and/or a decrease in abdominal fat weight which would contribute to the improvement of these traits.

Bladder cancer is the most common malignant tumor of the urinary system, and it has high incidence, high degree of malignancy, and easy recurrence after surgery. The etiology and pathogenesis of bladder cancer are not fully understood, but more and more studies have shown that its development may be regulated by some core molecules. To identify key molecules in bladder cancer, we constructed a three-layer network by merging lncRNA-miRNA regulatory network, miRNA-mRNA regulatory network, and lncRNA-mRNA coexpression network, and further analyzed the topology attributes of the network including the degree, betweenness centrality and closeness centrality of nodes. We found that miRNA-93 and miRNA-195 are controllers for a three-layer network and regulators of numerous target genes associated with bladder cancer. Functional enrichment analysis of their target mRNAs revealed that miRNA-93 and miRNA-195 may be closely related to bladder cancer by disturbing the homeostasis of the cell cycle or HTLV-I infection. In addition, since E2F1 and E2F2 are enriched in various KEGG signaling pathways, we conclude that they are important target genes of miRNA-93, and participate in the apoptotic process by forming a complex with a certain protein or transcription factor activity, sequence-specific DNA binding in bladder cancer. Similarly, AKT3 is an important target gene of miRNA-195, its expression is associated with PI3K-Akt-mTOR signaling pathway and AMPK-mTOR signaling pathway. Therefore, we speculate that AKT3 may participate in proliferation and apoptosis of bladder cancer cells through these pathways, and ultimately affect the biological behavior of tumor cells. Furthermore, through survival analysis, we found that miRNA-195 and miRNA-93 are associated with poor prognosis of bladder cancer. And the Kaplan-Meier curve showed that 24 mRNAs and nine lncRNAs are closely related to overall survival of bladder cancer.

Autosomal recessive non-syndromic hearing loss (ARNSHL) is a highly heterogeneous disease involving more than 70 pathogenic genes. However, most ARNSHL families have small-sized pedigrees with limited genetic information, rendering challenges for the molecular diagnosis of these patients. Therefore, we attempted to establish a strategy for identifying deleterious variants associated with ARNSHL by applying proband whole-exome sequencing (proband-WES). Aside from desiring to improve molecular diagnostic rates, we also aimed to search for novel deafness genes shared by patients with similar phenotype, making up for the deficiency of small ARNSHL families. In this study, 48.5% (16/33) families were detected the pathogenic variants in eight known deafness genes, including 10 novel variants identified in TMPRSS3 (MIM 605551), MYO15A (MIM 602666), TMC1 (MIM 606706), ADGRV1 (MIM 602851), and PTPRQ (MIM 603317). Apart from six novel variants with a truncating effect (nonsense, deletion, insertion, and splice-site), four novel missense variants were not found in 200 unrelated control population by using Sanger sequencing. It is important to note that none of novel genes were shared across different pedigrees, indicating that a larger sample size might be needed. Proband-WES is a cost-effective and precise way of identifying causative variants in nuclear families with ARNSHL. This economical strategy may be appropriated as a clinical application to provide molecular diagnostics, genetic counseling, and individualized health maintenance measures for patients with ARNSHL at hearing clinics.

The Xichuan black-bone chicken, which is a rare local chicken species in China, is an important genetic resource of black-bone chickens. Tyrosine can affect melanin production, but the molecular mechanism underlying tyrosine-induced melanin deposition in Xichuan black-bone chickens is poorly understood. Here, the blackness degree and melanin content of the breast muscle of Xichuan black-bone chickens fed a basic diet with five levels of added tyrosine (i.e., 0.2, 0.4, 0.6, 0.8, and 1.0%; these groups were denoted test groups I-V, respectively) were assessed, and the results showed that 0.8% tyrosine was the optimal level of added tyrosine. Moreover, the effects of tyrosine supplementation on the proliferation and tyrosinase content of melanocytes in Xichuan black-bone chickens were evaluated. The results revealed a dose-dependent relationship between tyrosine supplementation and melanocyte proliferation. In addition, 417 differentially expressed genes (DEGs), including 160 upregulated genes and 257 downregulated genes, were identified in a comparative analysis of the transcriptome profiles constructed using the pooled total RNA from breast muscle tissues of the control group and test group IV, respectively (fold change ≥2.0, P < 0.05). These DEGs were mainly involved in melanogenesis, the calcium signaling pathway, the Wnt signaling pathway, the mTOR signaling pathway, and vascular smooth muscle contraction. The pathway analysis of the DEGs identified some key genes associated with pigmentation, such as DCT and EDNRB2. In summary, the melanin content of breast muscle could be markedly enhanced by adding an appropriate amount of tyrosine to the diet of Xichuan black-bone chickens, and the EDNRB2-mediated molecular regulatory network could play a key role in the biological process of tyrosine-induced melanin deposition. These results have deepened the understanding of the molecular regulatory mechanism of melanin deposition in black-bone chickens and provide a basis for the regulation of nutrition and genetic breeding associated with melanin deposition in Xichuan black-bone chickens.

Salmonella enteritidis (SE) is a foodborne pathogen that negatively affects both animal and human health. Controlling poultry SE infection will have great practical significance for human public health, as poultry are considered to be important sources and carriers of the disease. In this study, the splenic transcriptomes of challenged-susceptible (S), challenged-resistant (R) and non-challenged (C) chicks (3-days old, specific-pathogen-free White Leghorn) were characterized in order to identify the immune-related gene markers and pathways. A total of 934 significant differentially expressed genes (DEGs) were identified in comparisons among the C, R and S birds. First reported here, the DEGs involved in the Forkhead box O (FoxO) signaling pathway, especially FoxO3, were identified as potential markers for host resistance to SE infection. The challenged-susceptible birds exhibited strong activation of the FoxO signaling pathway, which may be a major defect causing immune cell apoptosis as part of SE-induced pathology; these S birds also showed weak activation of mitogen-activated protein kinase (MAPK)-related genes, contrasting with strong splenic activation in the R birds. Interestingly, suppression of several pathways in the immune response against Salmonella, including cytokine-cytokine receptor interaction and Jak-STAT, was only found in S birds and there was evidence of cross-talk among these pathways, perhaps contributing to susceptibility to Salmonella infection. These findings will help facilitate understanding resistance and susceptibility to SE infection in the earliest phases of the host immune response through Salmonella-induced pathways, provide new approaches to develop strategies for SE prevention and treatment, and may enhance innate resistance by genetic selection in animals.

Introduction: LMNA splicing mutations occur in 9.1% of cases with cardiac involvement cases, but the phenotype and severity of disease they cause have not yet been systematically studied. The aim of this study was to understand the clinical and pathogenic characteristics of the LMNA splice-site mutation phenotype in patients with LMNA -related dilated cardiomyopathy (DCM) and sudden cardiac death (SCD). Methods and Results: First, we reported a novel family with LMNA -related DCM and SCD, and the clinical characteristics of all current patients with LMNA splicing mutations were further summarized through the ClinVar database. Seventeen families with a total of 134 individuals, containing a total of 15 LMNA splicing mutation sites, were enrolled. A total of 42 subjects (31.3%) had SCD. Compared without with the non-DCM group ( n = 56), the patients within the DCM group ( n = 78) presented a lower incidence of atrioventricular block (AVB) ( p = 0.015) and a higher incidence rates of non-sustained ventricular tachycardia ( p = 0.004),) and implantable cardioverter defibrillator (ICD) implantation ( p = 0.005). Kaplan‒Meier survival analysis showed that the patients with pacemaker (PM) implantation had a significantly reduced the occurrence of SCD compared to patientswith those without PM implantation (log-rank p < 0.001), while there was no significant difference in ICD implantation between the two groups (log-rank p = 0.73). Second, we identified the family that we reported with a mutation in an LMNA c.513+1 G>A mutation in the reported family, and pathogenic prediction analysis showed that the mutation site was extremely harmful. Next, we conducted gene expression levels and cardiac pathological biopsy studies on the proband of this family. We found that the expression of normal LMNA mRNA from the proband was significantly downregulated in peripheral blood mononuclear cells than incompared with healthy individuals. Finally, we comprehensively summarized the pathological characteristics of LMNA -related DCM, including hypertrophy, atrophy, fibrosis, white blood cell infiltration, intercalated disc remodeling, and downregulation of desmin and connexin 43 (Cx43) expression. Discussion: Above all, Cardiaccardiac involvement in patients with LMNA splice-site mutation presented with a high rate of SCD. Implanting a pacemaker significantly reduced the SCD rate in non-DCM patients with AVB. The pathogenic characterization was not only haveinvolved suppressed the expression of the healthy LMNA allele, but was also associated with abnormal expression and distribution of desmin and Cx43., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Ling, Hou, Jia, Lan, Wu, Wu, Jie, Liu, Huang, Wan, Li, Guo and Liang.)