bluuebunny/crossref_metadata_2025 huggingface dataset:
Allometric scaling and allocation patterns: Implications for predicting productivity across plant communities
Authors: Sileshi Gudeta W. , Nath Arun Jyoti , Kuyah Shem
Published: January 2023

As the application of allometry continues to expand, the variability in the allometry exponent has generated a great deal of debate in forest ecology. Some studies have reported counterintuitive values of the exponent, but the sources of such values have remained both unexplored and unexplained. Therefore, the objectives of our analyses were to: (1) uncover the global patterns of allometric variation in stem height with stem diameter, crown radius with stem diameter or stem height, crown depth with stem diameter, crown volume with stem diameter, crown depth with crown diameter, aboveground biomass with stem diameter or height, and belowground biomass with aboveground biomass; (2) assess variations in allometry parameters with taxonomic levels, climate zones, biomes and historical disturbance regimes; and (3) identify the sources of counterintuitive values of the allometry exponents. Here, we provide novel insights into the tight allometric co-variations between stem and crown dimensions and tree biomass. We also show a striking similarity in scaling across climate zones, biomes and disturbance regimes consistent with the allometry constraint hypothesis. We show that the central tendency of the exponent is toward 2/3 for the scaling of stem height with diameter, crown dimensions with stem diameter and height, 5/2–8/3 for the scaling of aboveground biomass with stem diameter, and 1 for the scaling of belowground biomass with aboveground biomass. This is indicative of an integrated growth regulation acting in tandem on growth in stem diameter, height, crown dimensions and biomass allocation. We also demonstrate that counterintuitive values of the exponent arise as artifacts of small sample sizes (N < 60), measurement errors, sampling biases and inappropriate regression techniques. We strongly recommend the use of larger sample sizes (N > 60) and representative samples of the target population when testing hypothesis about allometric variation. We also caution against conflation of statistical artifacts with violations of theoretical predictions.
Association of methylenetetrahydrofolate reductase (MTHFR) rs1801133 (677C>T) gene polymorphism with ischemic stroke risk in different populations: An updated meta-analysis
Authors: Zhao Lili , Li Tao , Dang Meijuan , Li Ye , Fan Hong , Hao Qian , Song Dingli , Lu Jialiang , Lu Ziwei , Jian Yating , Wang Heying , Wang Xiaoya , Wu Yulun , Zhang Guilian
Published: January 2023

Background: Recently, increasing evidence has implicated methylenetetrahydrofolate reductase (MTHFR) gene mutation as a risk factor for ischemic stroke (IS) in the general population. However, studies have been inconclusive and lack evidence on specific populations. We aim to determine whether the rs1801133 (NC_000001.11 (MTHFR):g. 677C>T (p.Ala222Val) variant, we termed as MTHFR rs1801133 (677 C>T), is linked to an increased risk of IS in different age groups and ancestry groups.

Methods: The literature relevant to our study was found by searching the PubMed, Cochrane Library, Web of Science, EMBASE, and CNKI databases. A random effect model analysis was used to calculate the pooled odds ratio (OR) and 95% confidence interval (CI) to evaluate any possible association. We conducted a subgroup analysis based on the age and ancestry groups of the included populations.

Results: As of March 2022, 1,925 citations had been identified in electronic databases, of which 96 studies involving 34,814 subjects met our eligibility criteria. A strong link was found between IS and the MTHFR gene rs1801133 (677C>T) polymorphism in all genetic models [dominant genetic model (OR = 1.47; 95%CI = 1.33–1.61; p < 0.001), recessive genetic model (OR = 1.52; 95%CI = 1.36–1.71; p < 0.001), heterozygous model (OR = 1.36; 95%CI = 1.24–1.48; p < 0.001), homozygous model (OR = 1.82; 95%CI = 1.58–2.11; p < 0.001), and T allelic genetic model (OR = 1.37; 95%CI = 1.27–1.48; p < 0.001)]. Further subgroup analyses indicated that the MTHFR rs1801133 (677C>T) variant may increase the risk of IS in Asian, Hispanic, or Latin population, middle-aged, and elderly populations (p < 0.001).

Conclusion: Our results implied that mutation of the T allele of MTHFR rs1801133 (677C>T) could be a risk factor for IS. A significant association was found among Asian, Hispanic, or Latin population, middle-aged, and elderly people.
Growth performance, survivability and profitability of improved smallholder chicken genetics in Nigeria: A COVID-19 intervention study
Authors: Bamidele Oladeji , Akinsola Oludayo Micheal , Yakubu Abdulmojeed , Hassan Waheed Akinola , Ogundu Uduak Emmanuel , Amole Tunde Adegoke
Published: January 2023

The impact of COVID-19 pandemic on smallholder farming households (SFH) includes increased poverty, and loss of livelihoods. Provision of livestock to SFH is a helpful intervention to mitigate this impact. This study provided a total of 150 smallholder poultry farmers, randomly selected from three states (Kebbi, Nasarawa, and Imo) in Nigeria, with ten 5-week-old chickens (mixed sexes) each, of either FUNAAB Alpha or Noiler chicken genetics. The improved, dual-purpose chickens were evaluated for growth performance (GP), survivability and profitability. The birds were managed under semi-scavenging production system. Body weight, mortality, and cost of production (COP) were recorded every 4 weeks until 21 weeks of age. Profitability was a function of the COP, and the selling price for live-birds (cocks). Body weight of Noiler (1,927 g) birds was not significantly (p > .05) higher than FUNAAB Alpha (1,792 g) at 21 weeks. Agroecology and genetics had significant (p < .05) effects on GP and survivability. Survivability of FUNAAB Alpha was higher (p < .05) than Noiler, with Nasarawa (81%–96%), having the highest (p < .0001) survival rate compared to Imo (62%–81%), and Kebbi (58%–75%). At 21 weeks, the number of cocks and hens differed significantly (p < .05) within the states (Imo: 2.4 ± .2 and 5.4 ± .3; Kebbi: 2.6 ± .2 and 5.5 ± .3; and Nasarawa: 2.9 ± .2 and 5.8 ± .3). Nasarawa (NGN 7,808; USD 19) ranked best for profitability, followed by Kebbi (NGN 6,545; USD 16) and Imo (NGN 5,875; USD 14). Overall, this study demonstrates that provision of improved chickens to vulnerable SFH in Nigeria holds great potential for economic growth, and resilience during emergencies, such as the COVID-19 pandemic.
Experimental verification and comprehensive analysis of m7G methylation regulators in the subcluster classification of ischemic stroke
Authors: Tian Yunze , Yu Beibei , Lv Boqiang , Zhang Yongfeng , Fu Longhui , Yang Shijie , Li Jianzhong , Gong Shouping
Published: January 2023

Background: Ischemic stroke (IS) is a fatal cerebrovascular disease involving several pathological mechanisms. Modification of 7-methylguanosine (m7G) has multiple regulatory functions. However, the expression pattern and mechanism of m7G in IS remain unknown. Herein, we aimed to explore the effect of m7G modification on IS.

Methods: We screened significantly different m7G-regulated genes in Gene Expression Omnibus datasets, GSE58294 and GSE22255. The random forest (RF) algorithm was selected to identify key m7G-regulated genes that were subsequently validated using the middle cerebral artery occlusion (MCAO) model and quantitative polymerase chain reaction (qPCR). A risk model was subsequently generated using key m7G-regulated genes. Then, “ConsensusClusterPlus” package was used to distinguish different m7G clusters of patients with IS. Simultaneously, between two m7G clusters, differentially expressed genes (DEGs) and immune infiltration differences were also explored. Finally, we investigated functional enrichment and the mRNA–miRNA–transcription factor network of DEGs.

Results: RF and qPCR confirmed that EIF3D, CYFIP2, NCBP2, DCPS, and NUDT1 were key m7G-related genes in IS that could accurately predict clinical risk (area under the curve = 0.967). NCBP2 was the most significantly associated gene with immune infiltration. Based on the expression profiles of these key m7G-related genes, the IS group could be divided into two clusters. According to the single-sample gene set enrichment analysis algorithm, four types of immune cells (immature dendritic cells, macrophages, natural killer T cells, and TH1 cells) were significantly different in the two m7G clusters. The functional enrichment of 282 DEGs between the two clusters was mainly concentrated in the “regulation of apoptotic signaling pathway,” “cellular response to DNA damage stimulus,” “adaptive immune system,” and “pyroptosis.” The miR-214–LTF–FOXJ1 axis may be a key regulatory pathway for IS.

Conclusion: Our findings suggest that EIF3D, CYFIP2, NCBP2, DCPS, and NUDT1 may serve as potential diagnostic biomarkers for IS and that the m7G clusters developed by these genes provide more evidence for the regulation of m7G in IS.
Prognostic implication and immunotherapy response prediction of a ubiquitination-related gene signature in breast cancer
Authors: Guo Yangyang , Chen Qiaoqiao , Zhang Yingjue , Cheng Xu , Cen Kenan , Dai Ying , Mai Yifeng , Hong Kai
Published: January 2023

Breast cancer (BC) is one of the most common tumor types and has poor outcomes. In this study, a ubiquitination-related prognostic signature was constructed, and its association with immunotherapy response in BC was explored. A list of ubiquitination-related genes was obtained from the molecular signatures database, and a ubiquitination-related gene signature was obtained by least absolute shrinkage and selection operator Cox regression. The genes, TCN1, DIRAS3, and IZUMO4, had significant influence on BC outcomes. Patients were categorized into two clusters—a high-risk group with poor survival and a low-risk group with greater chances of controlling BC progression. Univariate and multivariate Cox regression analyses revealed that the risk signature was an independent prognostic factor for BC. Gene set enrichment analysis suggested that the high-risk group was enriched in cell cycle and DNA replication pathways. The risk score was positively linked to the tumor microenvironment and negatively correlated with the immunotherapy response. The IC50 values for rapamycin were higher in the low-risk group, whereas those for axitinib, AZD6244, erlotinib, GDC0941, GSK650394, GSK269962A, lapatinib, and PD0325901 were higher in the high-risk group. Therefore, the ubiquitination-related signature is considered a promising tool for predicting a BC patient’s immunotherapy response.
Prognostic prediction of clear cell renal cell carcinoma based on lipid metabolism-related lncRNA risk coefficient model
Authors: Qu GenYi , Liu Lu , Yi Lai , Tang Cheng , Yang Guang , Chen Dan , Xu Yong
Published: January 2023

Objective: In order to predict the prognosis in patients with clear cell renal cell carcinoma (ccRCC) so as to understand cancer lipid metabolism and sensitivity to immune-targeting drugs, model algorithms were used to establish a risk coefficient model of long non-coding RNAs (lncRNAs) associated with lipid metabolism.

Methods: The transcriptome data were retrieved from TCGA, and lncRNAs associated with lipid metabolism were obtained through Pearson correlation and differential expression analyses. Differentially expressed lipid metabolism-related lncRNAs and lipid metabolism-related lncRNA pairs were obtained using the R language software. The minimum absolute shrinkage method and the selector operation regression method were used to construct the model and draw the receiver operator characteristic curve. High-risk patients were differentiated from low-risk patients through the cut-off value, and the correlation analyses of the high-risk subgroup and low-risk subgroup were performed.

Results: This research discovered that 25 pairs of lncRNAs were associated with the lipid metabolism of ccRCC, and 12 of these pairs were utilized to build the model. In combination with clinical data, the areas under the 1-, 3- and 5-year survival curves of ccRCC patients were 0.809, 0.764 and 0.792, separately. The cut-off value was used to perform subgroup analysis. The results showed that high-risk patients had poor prognosis. The results of Cox multivariate regressive analyses revealed that age and risk score were independent prediction factors of ccRCC prognosis. In addition, immune cell infiltration, the levels of gene expression at immune checkpoints, and high-risk patients more susceptible to sunitinib-targeted treatment were assessed by the risk model.

Conclusion: Our team identified new prognostic markers of ccRCC and established risk models that could assess the prognosis of ccRCC patients and help determine which type of patients were more susceptible to sunitinib. These discoveries are vital for the optimization of risk stratification and personalized management.
Association of environmental exposure to chromium with differential DNA methylation: An epigenome-wide study
Authors: Zhao Meiduo , Wu Jingtao , Xu Jing , Li Ang , Mei Yayuan , Ge Xiaoyu , Yin Guohuan , Liu Xiaolin , Wei Lanping , Xu Qun
Published: January 2023

Introduction: Previous studies have reported that chromium (Cr)-induced epigenetic alterations and DNA methylation play a vital role in the pathogenesis of diseases induced by chromium exposure. Epigenomic analyses have been limited and mainly focused on occupational chromium exposure; their findings are not generalizable to populations with environmental Cr exposure.

Methods: We identified the differential methylation of genes and regions to elucidate the mechanisms of toxicity related to environmental chromium exposure. DNA methylation was measured in blood samples collected from individuals in Cr-contaminated (n = 10) and unexposed areas (n = 10) by using the Illumina Infinium HumanMethylation850K array. To evaluate the relationship between chromium levels in urine and CpG methylation at 850 thousand sites, we investigated differentially methylated positions (DMPs) and differentially methylated regions (DMRs) by using linear models and DMRcate method, respectively. The model was adjusted for biologically relevant variables and estimated cell-type compositions.

Results: At the epigenome-wide level, we identified five CpGs [cg20690919 (pFDR =0.006), cg00704664 (pFDR =0.024), cg10809143 (pFDR =0.043), cg27057652 (pFDR =0.047), cg05390480 (pFDR =0.024)] and one DMR (chr17: 19,648,718-19,648,972), annotated to ALDH3A1 genes (p < 0.05) as being significantly associated with log2 transformed urinary chromium levels.

Discussion: Environmental chromium exposure is associated with DNA methylation, and the significant DMPs and DMR being annotated to cause DNA damage and genomic instability were found in this work. Research involving larger samples is required to further explore the epigenetic effect of environmental chromium exposure on health outcomes through DNA methylation.
Integrative approach for validation of six important fish species inhabiting River Poonch of north-west Himalayan region (India)
Authors: Awas Mohd , Ahmed Imtiaz , Ahmad Syed Mudasir , Al-Anazi Khalid Mashay , Farah Mohammad Abul , Bhat Basharat
Published: January 2023

Traditionally, species of fish are identified based on morphological characteristics. Although these taxonomic descriptions are essential, there are cases where the morphological characters distinguishing these species show marginal differences. For instance, in the Poonch River in the Himalayas, there are 21 species, out of which some are morphologically similar, and the taxonomic distinction between these species is unclear. Therefore, in this study, we used sequences from two mitochondrial genes, Cytochrome b (Cyt b) and a larger ribosomal subunit (16S rRNA), as well as the morphological analysis to address any taxonomic ambiguities among the six fish species. Maximum Likelihood results revealed that all the species were clustered according to their families and genera. The phenotypic analysis also supported this statement, as all the species of different genera likeSchizothorax,Tor,Garra,Traqilabeo, andGlyptothoraxare grouped in their particular cluster, it shows that species of a separate class share a mutual morphological characteristic. While genetic analyses of these species suggest nucleotide diversity (p) and haplotype diversity, with Hd values as 0.644 forCyt band 0.899 for 16S rRNA, confirming the rich genetic diversity in the river. Overall, we recommend that the integrative approach in delimiting the fish species is more effective than the individual one and can be used to rapidly diagnose a species and understand the evolutionary relationship between the species.
Comparative transcriptomic analysis of long noncoding RNAs in Leishmania-infected human macrophages
Authors: Fernandes Juliane C. R. , Gonçalves André N. A. , Floeter-Winter Lucile M. , Nakaya Helder I. , Muxel Sandra M.
Published: January 2023

It is well established that infection with Leishmania alters the host cell’s transcriptome. Since mammalian cells have multiple mechanisms to control gene expression, different molecules, such as noncoding RNAs, can be involved in this process. MicroRNAs have been extensively studied upon Leishmania infection, but whether long noncoding RNAs (lncRNAs) are also altered in macrophages is still unexplored. We performed RNA-seq from THP-1-derived macrophages infected with Leishmania amazonensis (La), L. braziliensis (Lb), and L. infantum (Li), investigating a previously unappreciated fraction of macrophage transcriptome. We found that more than 24% of the total annotated transcripts and 30% of differentially expressed (DE) RNAs in Leishmania-infected macrophage correspond to lncRNAs. LncRNAs and protein coding RNAs with altered expression are similar among macrophages infected with the Leishmania species. Still, some species-specific alterations could occur due to distinct pathophysiology in which Li infection led to a more significant number of exclusively DE RNAs. The most represented classes among DE lncRNAs were intergenic and antisense lncRNAs. We also found enrichment for immune response-related pathways in the DE protein coding RNAs, as well as putative targets of the lncRNAs. We performed a coexpression analysis to explore potential cis regulation of coding and antisense noncoding transcripts. We identified that antisense lncRNAs are similarly regulated as its neighbor protein coding genes, such as the BAALC/BAALC-AS1, BAALC/BAALC-AS2, HIF1A/HIF1A-AS1, HIF1A/HIF1A-AS3 and IRF1/IRF1-AS1 pairs, which can occur as a species-specific modulation. These findings are a novelty in the field because, to date, no study has focused on analyzing lncRNAs in Leishmania-infected macrophage. Our results suggest that lncRNAs may account for a novel mechanism by which Leishmania can control macrophage function. Further research must validate putative lncRNA targets and provide additional prospects in lncRNA function during Leishmania infection.
Forensic efficiencies of individual identification, kinship testing and ancestral inference in three Yunnan groups based on a self-developed multiple DIP panel
Authors: Chen Man , Lan Qiong , Nie Shengjie , Hu Liping , Fang Yating , Cui Wei , Bai Xiaole , Liu Liu , Zhu Bofeng
Published: January 2023

Deletion/insertion polymorphism (DIP), as a short insertion/deletion sequence polymorphic genetic marker, has attracted the attention of forensic genetic scientist due to its lack of stutter, short amplicon and abundant ancestral information. In this study, based on a self-developed 43 autosomal deletion/insertion polymorphism (A-DIP) loci panel which could meet the forensic application purposes of individual identification, kinship testing and ancestral inference to some extent, we evaluated the forensic efficiencies of the above three forensic objectives in Chinese Yi, Hani and Miao groups of Yunnan province. The cumulative match probability (CPM) and combined probability of exclusion (CPE) of these three groups were 1.11433E-18, 8.24299E-19, 4.21721E-18; 0.999610217, 0.999629285 and 0.999582084, respectively. Average 96.65% full sibling pairs could be identified from unrelated individual pairs (as likelihood ratios > 1) using this DIP panel, whereas the average false positive rate was 3.69% in three target Yunnan groups. With the biogeographical ancestor prediction models constructed by extreme gradient boosting (XGBoost) and support vector machine (SVM) algorithms, 0.8239 (95% CI 0.7984, 0.8474) of the unrelated individuals could be correctly divided according to the continental origins based on the 43 A-DIPs which were large frequency distribution differentiations among different continental populations. The present results of principal component analysis (PCA), multidimensional scaling (MDS), neighbor joining (NJ) and maximum likelihood (ML) phylogenetic trees and STRUCTURE analyses indicated that these three Yunnan groups had relatively close genetic distances with East Asian populations.
Differential DNA methylation associated with multiple sclerosis and disease modifying treatments in an underrepresented minority population
Authors: Bingen Jeremy M. , Clark Lindsay V. , Band Mark R. , Munzir Ilyas , Carrithers Michael D.
Published: January 2023

Black and Hispanic American patients frequently develop earlier onset of multiple sclerosis (MS) and a more severe disease course that can be resistant to disease modifying treatments. The objectives were to identify differential methylation of genomic DNA (gDNA) associated with disease susceptibility and treatment responses in a cohort of MS patients from underrepresented minority populations. Patients with MS and controls with non-inflammatory neurologic conditions were consented and enrolled under an IRB-approved protocol. Approximately 64% of donors identified as Black or African American and 30% as White, Hispanic-Latino. Infinium MethylationEPIC bead arrays were utilized to measure epigenome-wide gDNA methylation of whole blood. Data were analyzed in the presence and absence of adjustments for unknown covariates in the dataset, some of which corresponded to disease modifying treatments. Global patterns of differential methylation associated with MS were strongest for those probes that showed relative demethylation of loci with lower M values. Pathway analysis revealed unexpected associations with shigellosis and amoebiasis. Enrichment analysis revealed an over-representation of probes in enhancer regions and an under-representation in promoters. In the presence of adjustments for covariates that included disease modifying treatments, analysis revealed 10 differentially methylated regions (DMR’s) with an FDR <1E-77. Five of these genes (ARID5B, BAZ2B, RABGAP1, SFRP2, WBP1L) are associated with cancer risk and cellular differentiation and have not been previously identified in MS studies. Hierarchical cluster and multi-dimensional scaling analysis of differential DNA methylation at 147 loci within those DMR’s was sufficient to differentiate MS donors from controls. In the absence of corrections for disease modifying treatments, differential methylation in patients treated with dimethyl fumarate was associated with immune regulatory pathways that regulate cytokine and chemokine signaling, axon guidance, and adherens junctions. These results demonstrate possible associations of gastrointestinal pathogens and regulation of cellular differentiation with MS susceptibility in our patient cohort. This work further suggests that analyses can be performed in the presence and absence of corrections for immune therapies. Because of their high representation in our patient cohort, these results may be of specific relevance in the regulation of disease susceptibility and treatment responses in Black and Hispanic Americans.
Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata)
Authors: Hagadorn Mallory A. , Hunter Frances K. , DeLory Tim , Johnson Makenna M. , Pitts-Singer Theresa L. , Kapheim Karen M.
Published: January 2023

Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.
Construction of HBV gene-related prognostic and diagnostic models for hepatocellular carcinoma
Authors: Ma Keqiang , Wu Hongsheng , Ji Lei
Published: January 2023

Background: Hepatocellular carcinoma (HCC) is a main cause of malignancy-related death all over the world with a poor prognosis. The current research is focused on developing novel prognostic and diagnostic models of Hepatocellular carcinoma from the perspective of hepatitis B virus (HBV)-related genes, and predicting its prognostic characteristics and potential reliable biomarkers for Hepatocellular carcinoma diagnosis.

Methods: As per the information related to Hepatocellular carcinoma expression profile and the clinical data in multiple public databases, we utilized limma for assessing the differentially expressed genes (DEGs) in HBV vs non- hepatitis B virus groups, and the gene set was enriched, analyzed and annotated by WebGestaltR package. Then, STRING was employed to investigate the protein interactions. A risk model for evaluating Hepatocellular carcinoma prognosis was built with Lasso Cox regression analysis. The effect patients receiving immunotherapy was predicted using Tumor Immune Dysfunction and Exclusion (TIDE). Additionally, pRRophetic was used to investigate the drug sensitivity. Lastly, the Support Vector Machine (SVM) approach was utilized for building the diagnostic model.

Results: The Hepatocellular Carcinoma Molecular Atlas 18 (HCCDB18) data set was utilized for the identification of 1344 HBV-related differentially expressed genes, mainly associated with cell division activities. Five functional modules were established and then we built a prognostic model in accordance with the protein-protein interaction (PPI) network. Five HBV-related genes affecting prognosis were identified for constructing a prognostic model. Then, the samples were assigned into RS-high and -low groups as per their relevant prognostic risk score (RS). High-risk group showed worse prognosis, higher mutation rate of TP53, lower sensitivity to immunotherapy but higher response to chemotherapeutic drugs than low-risk group. Finally, the hepatitis B virus diagnostic model of Hepatocellular carcinoma was established.

Conclusion: In conclusion, the prognostic and diagnostic models of hepatitis B virus gene-related Hepatocellular carcinoma were constructed. ABCB6, IPO7, TIMM9, FZD7, and ACAT1, the five HBV-related genes that affect the prognosis, can work as reliable biomarkers for the diagnosis of Hepatocellular carcinoma, giving a new insight for improving the prognosis, diagnosis, and treatment outcomes of HBV-type Hepatocellular carcinoma.
Whole-genome characterization of large-cell lung carcinoma: A comparative analysis based on the histological classification
Authors: Wu Xiaowei , Yin Jin , Deng Yu , Zu Yukun
Published: January 2023

Background: According to the 2015 World Health Organization classification, large cell neuroendocrine carcinoma (LCNEC) was isolated from Large-cell lung cancer (LCLC) tumors, which constitutes 2%–3% of non-small cell lung cancer (NSCLC). However, LCLC tumors are still fairly vaguely defined at the molecular level compared to other subgroups.

Materials and Methods: In this study, whole-genome sequencing (WGS) was performed on 23 LCLC and 15 LCNEC tumor specimens. Meanwhile, data from the TCGA (586 LUADs and 511 LUSCs) and U Cologne (120 SCLCs) were analyzed and compared.

Results: The most common driver mutations were found in TP53 (13/23, 57%), FAM135B (8/23, 35%) and FAT3 (7/23, 30%) in LCLC, while their counterparts in LCNEC were TP53 (13/15, 87%), LRP1B (6/15, 40%) and FAT1 (6/15, 40%). Notably, FAM135B mutations only occurred in LCLC (P = 0.013). Cosmic signature analysis revealed widespread defective DNA mismatch repair and tobacco-induced mutations in both LCLC and LCNEC. Additionally, LCNEC had a higher incidence of chromosomal copy number variations (CNVs) and structural variations (SVs) compared with LCLC, although the differences were not statistically significant. Particularly, chromothripsis SVs was significantly associated with CNVs. Furthermore, mutational landscape of different subtypes indicated differences between subtypes, and there seems to be more commonalty between our cohort and SCLC than with other subtypes. SMARCA4 mutations may be specific driver gene alteration in our cohort.

Conclusion: Our results support that LCLC and LCNEC tumors follow distinct tumorigenic pathways. To our knowledge, this is the first genome-wide profiling comparison of LCLC and LCNEC.
Biological characteristics of γδT cells and application in tumor immunotherapy
Authors: Zhu Renhong , Yan Qian , Wang Yashu , Wang Keqiang
Published: January 2023

Human γδT cells are a special immune cell type which exist in small quantities in the body, do not require processing and presentation for antigen recognition, and have non-major histocompatibility complex (MHC)-restricted immune response. They play an important role in the body’s anti-tumor, anti-infection, immune regulation, immune surveillance and maintenance of immune tolerance. This article reviews the generation and development of human γδT cells, genetic characteristics, classification, recognition and role of antigens, and research progress in tumor immunotherapy.
CRISPR/Cas9-editing of KISS1 to generate pigs with hypogonadotropic hypogonadism as a castration free trait
Authors: Flórez Julio M. , Martins Kyra , Solin Staci , Bostrom Jonathan R. , Rodríguez-Villamil Paula , Ongaratto Felipe , Larson Sabreena A. , Ganbaatar Uyanga , Coutts Alexander W. , Kern Doug , Murphy Thomas W. , Kim Eui-Soo , Carlson Daniel F. , Huisman Abe , Sonstegard Tad S. , Lents Clay A.
Published: January 2023

Introduction: Most male pigs are surgically castrated to avoid puberty-derived boar taint and aggressiveness. However, this surgical intervention represents a welfare concern in swine production. Disrupting porcine KISS1 is hypothesized to delay or abolish puberty by inducing variable hypogonadotropism and thus preventing the need for castration.

Methods: To test this hypothesis, we generated the first KISS1-edited large animal using CRISPR/Cas9-ribonucleoproteins and single-stranded donor oligonucleotides. The targeted region preceded the sequence encoding a conserved core motif of kisspeptin. Genome editors were intracytoplasmically injected into 684 swine zygotes and transferred to 19 hormonally synchronized surrogate sows. In nine litters, 49 American Yorkshire and 20 Duroc liveborn piglets were naturally farrowed.

Results: Thirty-five of these pigs bore KISS1-disruptive alleles ranging in frequency from 5% to 97% and did not phenotypically differ from their wild-type counterparts. In contrast, four KISS1-edited pigs (two boars and two gilts) with disruptive allele frequencies of 96% and 100% demonstrated full hypogonadotropism, infantile reproductive tracts, and failed to reach sexual maturity. Change in body weight during development was unaffected by editing KISS1. Founder pigs partially carrying KISS1-disruptive alleles were bred resulting in a total of 53 KISS1+/+, 60 KISS1+/−, and 34 KISS1−/− F1 liveborn piglets, confirming germline transmission.

Discussion: Results demonstrate that a high proportion of KISS1 alleles in pigs must be disrupted before variation in gonadotropin secretion is observed, suggesting that even a small amount of kisspeptin ligand is sufficient to confer proper sexual development and puberty in pigs. Follow-on studies will evaluate fertility restoration in KISS1 KO breeding stock to fully realize the potential of KISS1 gene edits to eliminate the need for surgical castration.
Genome-wide identification, evolution and transcriptome analysis of GRAS gene family in Chinese chestnut (Castanea mollissima)
Authors: Yu Liyang , Hui Cai , Huang Ruimin , Wang Dongsheng , Fei Cao , Guo Chunlei , Zhang Jingzheng
Published: January 2023

GRAS transcription factors play an important role in regulating various biological processes in plant growth and development. However, their characterization and potential function are still vague in Chinese chestnut (Castanea mollissima), an important nut with rich nutrition and high economic value. In this study, 48 CmGRAS genes were identified in Chinese chestnut genome and phylogenetic analysis divided CmGRAS genes into nine subfamilies, and each of them has distinct conserved structure domain and features. Genomic organization revealed that CmGRAS tend to have a representative GRAS domain and fewer introns. Tandem duplication had the greatest contribution for the CmGRAS expansion based on the comparative genome analysis, and CmGRAS genes experienced strong purifying selection pressure based on the Ka/Ks. Gene expression analysis revealed some CmGRAS members with potential functions in bud development and ovule fertility. CmGRAS genes with more homologous relationships with reference species had more cis-acting elements and higher expression levels. Notably, the lack of DELLA domain in members of the DELLA subfamily may cause de functionalization, and the differences between the three-dimensional structures of them were exhibited. This comprehensive study provides theoretical and practical basis for future research on the evolution and function of GRAS gene family.
A novel cuproptosis-related lncRNA signature predicts prognosis and therapeutic response in bladder cancer
Authors: Cai Jinming , Xie Haoran , Yan Yilin , Huang Zhengnan , Tang Pengfei , Cao Xiangqian , Wang Zeyi , Yang Chenkai , Wen Jiling , Tan Mingyue , Zhang Fang , Shen Bing
Published: January 2023

Bladder cancer (BC) ranks the tenth in the incidence of global tumor epidemiology. LncRNAs and cuproptosis were discovered to regulate the cell death. Herein, we downloaded transcriptome profiling, mutational data, and clinical data on patients from The Cancer Genome Atlas (TCGA). High- and low-risk BC patients were categorized. Three CRLs (AL590428.1, AL138756.1 and GUSBP11) were taken into prognostic signature through least absolute shrinkage and selection operator (LASSO) Cox regression. Worse OS and PFS were shown in high-risk group (p < 0.05). ROC, independent prognostic analyses, nomogram and C-index were predicted via CRLs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated IncRNAs play a biological role in BC progression. Immune-related functions showed the high-risk group received more benefit from immunotherapy and had stronger immune responses, and the overall survival was better (p < 0.05). Finally, a more effective outcome (p < 0.05) was found from clinical immunotherapy via the TIDE algorithm and many potential anti-tumor drugs were identified. In our study, the cuproptosis-related signature provided a novel tool to predict the prognosis in BC patients accurately and provided a novel strategy for clinical immunotherapy and clinical applications.
Recent advances in epigenetic anticancer therapeutics and future perspectives
Authors: Ren Liwen , Yang Yihui , Li Wan , Yang Hong , Zhang Yizhi , Ge Binbin , Zhang Sen , Du Guanhua , Wang Jinhua
Published: January 2023

Tumor development is frequently accompanied by abnormal expression of multiple genomic genes, which can be broadly viewed as decreased expression of tumor suppressor genes and upregulated expression of oncogenes. In this process, epigenetic regulation plays an essential role in the regulation of gene expression without alteration of DNA or RNA sequence, including DNA methylation, RNA methylation, histone modifications and non-coding RNAs. Therefore, drugs developed for the above epigenetic modulation have entered clinical use or preclinical and clinical research stages, contributing to the development of antitumor drugs greatly. Despite the efficacy of epigenetic drugs in hematologic caners, their therapeutic effects in solid tumors have been less favorable. A growing body of research suggests that epigenetic drugs can be applied in combination with other therapies to increase efficacy and overcome tumor resistance. In this review, the progress of epigenetics in tumor progression and oncology drug development is systematically summarized, as well as its synergy with other oncology therapies. The future directions of epigenetic drug development are described in detail.
Construction of an original anoikis-related prognostic model closely related to immune infiltration in gastric cancer
Authors: Zhao Zhihong , Li Cun , Peng Ye , Liu Rui , Li Qian
Published: January 2023

Background: Anoikis is considered as a particular type of programmed cell death, the weakness or resistance of which contributes greatly to the development and progression of most malignant solid tumors. However, the latent impact of anoikis-related genes (ARGs) on gastric cancer (GC) is still ambiguous. Based on these, this study established an anoikis-related prognostic model of GC to identify the prognosis of patients and provide more effective treatment in clinical practice.

Methods: First, we extracted four public datasets containing the gene expression and clinicopathological information of GC, which were worked as the training and validating sets, separately. Then, an anoikis-related survival-predicted model of GC was developed via Lasso and COX regression analyses and verified by using the Kaplan-Meier (KM) curve and receiver operating characteristic (ROC) curve analyses. Next, we assigned GC patients to two groups characterized by the risk score calculated and analyzed somatic mutation, functional pathways, and immune infiltration between the different two groups. Finally, a unique nomogram was offered to clinicians to forecast the personal survival probability of GC patients.

Results: Based on seven anoikis-related markers screened and identified, a carcinogenic model of risk score was produced. Patients placed in the high-score group suffered significantly worse overall survival (OS) in four cohorts. Additionally, the model revealed a high sensitivity and specificity to prognosticate the prognoses of GC patients [area under the ROC curve (AUC) at 5-year = 0.713; GSE84437, AUC at 5-year = 0.639; GSE15459, AUC at 5-year = 0.672; GSE62254, AUC at 5-year = 0.616]. Apart from the excellent predictive performance, the model was also identified as an independent prediction factor from other clinicopathological characteristics. Combining anoikis-related prognostic model with GC clinical features, we built a more comprehensive nomogram to foresee the likelihood of survival of GC patients in a given year, showing a well-accurate prediction performance.

Conclusion: In summary, this study created a new anoikis-related signature for GC, which has potentially provided new critical insights into survival prediction and individualized therapy development.
A machine learning-based approach to ERα bioactivity and drug ADMET prediction
Authors: An Tianbo , Chen Yueren , Chen Yefeng , Ma Leyu , Wang Jingrui , Zhao Jian
Published: January 2023

By predicting ERα bioactivity and mining the potential relationship between Absorption, Distribution, Metabolism, Excretion, Toxicity (ADMET) attributes in drug research and development, the development efficiency of specific drugs for breast cancer will be effectively improved and the misjudgment rate of R&D personnel will be reduced. The quantitative prediction model of ERα bioactivity and classification prediction model of Absorption, Distribution, Metabolism, Excretion, Toxicity properties were constructed. The prediction results of ERα bioactivity were compared by XGBoot, Light GBM, Random Forest and MLP neural network. Two models with high prediction accuracy were selected and fused to obtain ERα bioactivity prediction model from Mean absolute error (MAE), mean squared error (MSE) and R2. The data were further subjected to model-based feature selection and FDR/FPR-based feature selection, respectively, and the results were placed in a voting machine to obtain Absorption, Distribution, Metabolism, Excretion, Toxicity classification prediction model. In this study, 430 molecular descriptors were removed, and finally 20 molecular descriptors with the most significant effect on biological activity obtained by the dual feature screening combined optimization method were used to establish a compound molecular descriptor prediction model for ERα biological activity, and further classification and prediction of the Absorption, Distribution, Metabolism, Excretion, Toxicity properties of the drugs were made. Eighty variables were selected by the model ExtraTreesClassifier Classifie, and 40 variables were selected by the model GradientBoostingClassifier to complete the model-based feature selection. At the same time, the feature selection method based on FDR/FPR is also selected, and the three classification models obtained by the two methods are placed into the voting machine to obtain the final model. The experimental results showed that the model‘s evaluation indexes and roc diagram were excellent and could accurately predict ERα bioactivity and Absorption, Distribution, Metabolism, Excretion, Toxicity properties. The model constructed in this study has high accuracy, fast convergence and robustness, has a very high accuracy for Absorption, Distribution, Metabolism, Excretion, Toxicity and ERα classification prediction, has bright prospects in the biopharmaceutical field, and is an important method for energy conservation and yield increase in the future.
Differential methylation of microRNA encoding genes may contribute to high myopia
Authors: Swierkowska Joanna , Vishweswaraiah Sangeetha , Mrugacz Malgorzata , Radhakrishna Uppala , Gajecka Marzena
Published: January 2023

Introduction: High myopia (HM), an eye disorder with a refractive error ≤−6.0 diopters, has multifactorial etiology with environmental and genetic factors involved. Recent studies confirm the impact of alterations in DNA methylation and microRNAs (miRNAs) on myopia. Here, we studied the combined aspects evaluating to the role of methylation of miRNA encoding genes in HM.

Materials and Methods: From the genome-wide DNA methylation data of 18 Polish children with HM and 18 matched controls, we retrieved differentially methylated CG dinucleotides localized in miRNA encoding genes. Putative target genes of the highest-ranked miRNAs were obtained from the miRDB and included in overrepresentation analyses in the ConsensusPathDB. Expression of target genes was assessed using the RNA sequencing data of retinal ARPE-19 cell line.

Results: We identified differential methylation of CG dinucleotides in promoter regions of MIR3621, MIR34C, MIR423 (increased methylation level), and MIR1178, MIRLET7A2, MIR885, MIR548I3, MIR6854, MIR675, MIRLET7C, MIR99A (decreased methylation level) genes. Several targets of these miRNAs, e.g. GNAS, TRAM1, CTNNB1, EIF4B, TENM3 and RUNX were previously associated with myopia/HM/refractive error in Europeans in genome-wide association studies. Overrepresentation analyses of miRNAs’ targets revealed enrichment in pathways/processes related to eye structure/function, such as axon guidance, transcription, focal adhesion, and signaling pathways of TGF-β, insulin, MAPK and EGF-EGFR.

Conclusion: Differential methylation of indicated miRNA encoding genes might influence their expression and contribute to HM pathogenesis via disrupted regulation of transcription of miRNAs’ target genes. Methylation of genes encoding miRNAs may be a new direction in research on both the mechanisms determining HM and non-invasive indicators in diagnostics.
A microarray data analysis investigating the pathogenesis and potential biomarkers of autophagy and ferroptosis in intervertebral disc degeneration
Authors: Kuang Wenhao , Jiang Cong , Yu Cheng , Hu Jinwei , Duan Yang , Chen Zhong
Published: January 2023

Background: Intervertebral disc degeneration (IDD) entails complex pathological changes and causes lower back pain (LBP). However, there is still a lack of understanding of the mechanisms involved in IDD, particularly regarding the roles of autophagy and ferroptosis. The current study used microarray data to investigate the pathogenesis of IDD and potential biomarkers related to autophagy and ferroptosis in IDD.

Methods: Differentially expressed genes (DEGs) were identified by analyzing the mRNA and miRNA expression profiles of IDD patients from the Gene Expression Omnibus (GEO). The protein-protein interaction network, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and gene set enrichment analysis (GSEA) were utilized. The Human Autophagy Database (HADb) and Ferroptosis Database were used in conjunction with hub genes to identify autophagy- and ferroptosis-related genes. The Transcription Factor -hub gene-miRNA network was constructed. Lastly, the expression of DEGs in normal and degenerated nucleus pulposus cells (NPCs) was investigated via the quantitative reverse transcription polymerase chain reaction (qRT-PCR).

Results: A total of 362 DEGs associated with IDD were identified. GO and KEGG analyses indicated that oxidative stress, extracellular matrix, PI3K-AKT signaling pathway, and ferroptosis were key factors in IDD occurrence. GSEA indicated that IDD was associated with changes in autophagy, iron ion homeostasis, extracellular matrix, and oxidative stress. Eighty-nine hub genes were obtained, including five that were autophagy-related and three that were ferroptosis-related. Of these, TP53 and SESN2 were the intersections of autophagy- and ferroptosis-related genes. In qRT-PCR analysis, CANX, SLC38A1, and TP53 were downregulated in degenerative NPCs, whereas GNAI3, SESN2, and VAMP3 were upregulated.

Conclusion: The current study revealed aspects of autophagy- and ferroptosis-related genes involved in IDD pathogenesis, warranting further investigation.