Genome-wide SNP discovery and QTL mapping for fruit quality traits in inbred backcross lines (IBLs) of solanum pimpinellifolium using genotyping by sequencing

Deciphering salt stress responses in Solanum pimpinellifolium through high-throughput phenotyping

Soil salinity is a major environmental stressor affecting agricultural productivity worldwide. Understanding plant responses to salt stress is crucial for developing resilient crop varieties. Wild relatives of cultivated crops, such as wild tomato, Solanum pimpinellifolium, can serve as a useful resource to further expand the resilience potential of the cultivated germplasm, S. lycopersicum. In this study, we employed high-throughput phenotyping in the greenhouse and field conditions to explore salt stress responses of a S. pimpinellifolium diversity panel. Our study revealed extensive phenotypic variations in response to salt stress, with traits such as transpiration rate, shoot mass, and ion accumulation showing significant correlations with plant performance. We found that while transpiration was a key determinant of plant performance in the greenhouse, shoot mass strongly correlated with yield under field conditions. Conversely, ion accumulation was the least influential factor under greenhouse conditions. Through a Genome Wide Association Study, we identified candidate genes not previously associated with salt stress, highlighting the power of high-throughput phenotyping in uncovering novel aspects of plant stress responses. This study contributes to our understanding of salt stress tolerance in S. pimpinellifolium and lays the groundwork for further investigations into the genetic basis of these traits, ultimately informing breeding efforts for salinity tolerance in tomato and other crops.

Exploring the Global Trends of Bacillus, Trichoderma and Entomopathogenic Fungi for Pathogen and Pest Control in Chili Cultivation

Chili, renowned globally and deeply ingrained in various cultures. Regrettably, the onset of diseases instigated by pests and pathogens has inflicted substantial losses on chili crops, with some farms experiencing complete production decimation. Challenges confronting chili cultivation include threats from pathogenic microbes like Xanthomonas, Fusarium, Phytophthora, Verticillium, Rhizoctonia, Colletotrichium and Viruses, alongside pests such as whiteflies, mites, thrips, aphids, and fruit flies. While conventional farming practices often resort to chemical pesticides to combat these challenges, their utilization poses substantial risks to both human health and the environment. In response to this pressing issue, this review aims to evaluate the potential of microbe-based biological control as eco-friendly alternatives to chemical pesticides for chili cultivation. Biocontrol agents such as Bacillus spp., Trichoderma spp., and entomopathogenic fungi present safer and more environmentally sustainable alternatives to chemical pesticides. However, despite the recognised potential of biocontrol agents, research on their efficacy in controlling the array of pests and pathogens affecting chili farming remains limited. This review addresses this gap by evaluating the efficiency of biocontrol agents, drawing insights from existing studies conducted in other crop systems, regarding pest and pathogen management. Notably, an analysis of Scopus publications revealed fewer than 30 publications in 2023 focused on these three microbial agents. Intriguingly, India, as the world's largest chili producer, leads in the number of publications concerning Bacillus spp., Trichoderma spp., and entomopathogenic fungi in chili cultivation. Further research on microbial agents is imperative to mitigate infections and reduce reliance on chemical pesticides for sustainable chili production.

Genetic variation in NOTCH1 is associated with overweight and obesity in Brazilian elderly

Excessive weight (overweight and obesity) is a common disorder involving genetic and environmental factors, associated with cardiovascular diseases, type-2 diabetes, and others. NOTCH1 is critical for the maintenance of stem cells and adult tissues, being reported as a key player in metabolism and adipogenesis in animals. Thus, we test the hypothesis that NOTCH1 Single Nucleotide Polymorphisms (SNPs) are associated with excessive weight. Participants from the census-based cohort SABE (Saúde, Bem Estar e Envelhecimento-Health, Well-Being, and Aging), carried out in the city of São Paulo-Brazil, were stratified into cases and controls according to BMI. We filter the SNPs located at the start and end positions of NOTCH1 and 50 Kb on both sides. We selected SNPs with minor allelic frequency (MAF) greater than or equal to 0.01 and Hardy-Weinberg equilibrium (p > 0.05) and r2 ≥ 0.8. We performed an association study with genotypes and haplotypes, as well as in silico functional analysis of the identified SNPs. We observed an association of the SNP rs9411207 with the risk of excessive weight, under log-additive model, and the genotype distribution showed an increased frequency of homozygous TT (OR 1.50, CI 1.20-1.88; p = 0.0002). The haplotype GAT constructed from this and other SNPs in high Linkage Disequilibrium was more frequent in excessive-weight individuals (p = 0.003). In silico analyses suggested that these SNPs are likely to affect the transcription of NOTCH1 and other genes involved in adipogenesis and metabolism. This is the first study reporting association between NOTCH1 SNPs and the risk of excessive weight. Considering the possibility of NOTCH1 modulation, additional population studies are important to replicate these data and confirm the usefulness of risk genotypes for management strategies of excessive weight.

Diversification of FT-like genes in the PEBP family contributes to the variation of flowering traits in Sapindaceae species

Many species of Sapindaceae, such as lychee, longan, and rambutan, provide nutritious and delicious fruit. Understanding the molecular genetic mechanisms that underlie the regulation of flowering is essential for securing flower and fruit productivity. Most endogenous and exogenous flowering cues are integrated into the florigen encoded by FLOWERING LOCUS T. However, the regulatory mechanisms of flowering remain poorly understood in Sapindaceae. Here, we identified 60 phosphatidylethanolamine-binding protein-coding genes from six Sapindaceae plants. Gene duplication events led to the emergence of two or more paralogs of the FT gene that have evolved antagonistic functions in Sapindaceae. Among them, the FT1-like genes are functionally conserved and promote flowering, while the FT2-like genes likely serve as repressors that delay flowering. Importantly, we show here that the natural variation at nucleotide position - 1437 of the lychee FT1 promoter determined the binding affinity of the SVP protein (LcSVP9), which was a negative regulator of flowering, resulting in the differential expression of LcFT1, which in turn affected flowering time in lychee. This finding provides a potential molecular marker for breeding lychee. Taken together, our results reveal some crucial aspects of FT gene family genetics that underlie the regulation of flowering in Sapindaceae.

Deleterious Effects of Heat Stress on the Tomato, Its Innate Responses, and Potential Preventive Strategies in the Realm of Emerging Technologies

The tomato is a fruit vegetable rich in nutritional and medicinal value grown in greenhouses and fields worldwide. It is severely sensitive to heat stress, which frequently occurs with rising global warming. Predictions indicate a 0.2 °C increase in average surface temperatures per decade for the next three decades, which underlines the threat of austere heat stress in the future. Previous studies have reported that heat stress adversely affects tomato growth, limits nutrient availability, hammers photosynthesis, disrupts reproduction, denatures proteins, upsets signaling pathways, and damages cell membranes. The overproduction of reactive oxygen species in response to heat stress is toxic to tomato plants. The negative consequences of heat stress on the tomato have been the focus of much investigation, resulting in the emergence of several therapeutic interventions. However, a considerable distance remains to be covered to develop tomato varieties that are tolerant to current heat stress and durable in the perspective of increasing global warming. This current review provides a critical analysis of the heat stress consequences on the tomato in the context of global warming, its innate response to heat stress, and the elucidation of domains characterized by a scarcity of knowledge, along with potential avenues for enhancing sustainable tolerance against heat stress through the involvement of diverse advanced technologies. The particular mechanism underlying thermotolerance remains indeterminate and requires further elucidatory investigation. The precise roles and interplay of signaling pathways in response to heat stress remain unresolved. The etiology of tomato plants' physiological and molecular responses against heat stress remains unexplained. Utilizing modern functional genomics techniques, including transcriptomics, proteomics, and metabolomics, can assist in identifying potential candidate proteins, metabolites, genes, gene networks, and signaling pathways contributing to tomato stress tolerance. Improving tomato tolerance against heat stress urges a comprehensive and combined strategy including modern techniques, the latest apparatuses, speedy breeding, physiology, and molecular markers to regulate their physiological, molecular, and biochemical reactions.

Establishment of a Biorepository for Down Syndrome: Experience of the Inter-Institutional Multidisciplinary BioBank - BioBIM

BACKGROUND

Down syndrome, or Trisomy 21, is the leading genetic cause of cognitive disability in children and is associated with a high risk of several comorbidities, particularly congenital heart defects, early onset Alzheimer's disease, leukaemia, and autoimmune disorders.

OBJECTIVE

This study describes the design, methods, and operational procedures employed to establish a biobank dedicated to Down syndrome that can support research projects investigating the effects of various genetic and environmental factors on this complex disease.

METHODS

Blood was collected from all recruited subjects, processed, aliquoted and immediately frozen at -80 °C in the Interinstitutional Multidisciplinary BioBank (BioBIM) facilities. A small aliquot of the sample was used to perform blood tests for which analysis would not be feasible at a later date, such as blood cell counts. Each biological sample was coded, assigned a Standard PREanalytical Code, and registered in the oloBIOBANK software connected to a medical card containing all the donor's anamnestic data. All samples were stored under continuous real-time temperature recording using a freezer connected to a T-GUARD alarm system. In addition, a radiofrequency identification tracking system strictly monitored each cryopreservation operation performed throughout the sample lifecycle.

RESULTS

Biological samples were collected from 454 individuals with Down syndrome from 2007 to 2023. A total of 2233 biological samples were available for research purposes, including whole blood in different anticoagulants, serum, plasma, and frozen peripheral blood mononuclear cells. The quality of the nucleic acids obtained through specific standard operating procedures demonstrated that these samples were appropriate for clinical and basic research.

CONCLUSION

By establishing this biobank, we have gathered a significant number of biological samples and clinical data from individuals with Down syndrome, thereby fostering collaboration between different research groups in an open and transparent manner. Sharing expertise and resources among scientists will ultimately facilitate the transfer of knowledge to clinical practice, leading to the development of more effective therapeutic treatments to improve the outcomes and quality of life of patients with Down syndrome.

Whole genome discovery of regulatory genes responsible for the response of chicken to heat stress

Long noncoding RNAs (lncRNAs) are functional bridges connecting the genome with phenotypes by interacting with DNA, mRNA, and proteins. Using publically available acute heat stress (AHS)-related RNA-seq data, we discovered novel lncRNAs and tested their association with AHS along with ~ 8800 known lncRNAs and ~ 28,000 mRNA transcripts. Our pipeline discovered a total of 145 potentially novel-lncRNAs. One of them (Fishcomb_p-value = 0.06) along with another novel transcript (annotated as protein-coding; Fishcomb_p-value = 0.03) were identified as significantly associated with AHS. We found five known-lncRNAs and 134 mRNAs transcripts that were significantly associated with AHS. Four novel lncRNAs interact cis-regulated with 12 mRNA transcripts and are targeted by 11 miRNAs. Also six meta-lncRNAs associate with 134 meta-mRNAs through trans-acting co-expression, each targeted by 15 and 216 miRNAs, respectively. Three of the known-lncRNAs significantly co-expressed with almost 97 of the significant mRNAs (Pearson correlation p-value < 0.05). We report the mentioned three known-lncRNAs (ENSGALT00000099876, ENSGALT00000107573, and ENSGALT00000106323) as the most, significantly regulatory elements of AHS in chicken. It can be concluded that in order to alleviate the adverse effects of AHS on chicken, the manipulation of the three regulatory lncRNAs could lead to a more desirable result than the manipulation of the most significant mRNAs.

Pangenome Analysis of Helicobacter pylori Isolates from Selected Areas of Africa Indicated Diverse Antibiotic Resistance and Virulence Genes

The challenge facing Helicobacter pylori (H. pylori) infection management in some parts of Africa is the evolution of drug-resistant species, the lack of gold standard in diagnostic methods, and the ineffectiveness of current vaccines against the bacteria. It is being established that even though clinical consequences linked to the bacteria vary geographically, there is rather a generic approach to treatment. This situation has remained problematic in the successful fight against the bacteria in parts of Africa. As a result, this study compared the genomes of selected H. pylori isolates from selected areas of Africa and evaluated their virulence and antibiotic drug resistance, those that are highly pathogenic and are associated with specific clinical outcomes and those that are less virulent and rarely associated with clinical outcomes. 146 genomes of H. pylori isolated from selected locations of Africa were sampled, and bioinformatic tools such as Abricate, CARD RGI, MLST, Prokka, Roary, Phandango, Google Sheets, and iTOLS were used to compare the isolates and their antibiotic resistance or susceptibility. Over 20 k virulence and AMR genes were observed. About 95% of the isolates were genetically diverse, 90% of the isolates harbored shell genes, and 50% harbored cloud and core genes. Some isolates did not retain the cagA and vacA genes. Clarithromycin, metronidazole, amoxicillin, and tinidazole were resistant to most AMR genes (vacA, cagA, oip, and bab). Conclusion. This study found both virulence and AMR genes in all H. pylori strains in all the selected geographies around Africa with differing quantities. MLST, Pangenome, and ORF analyses showed disparities among the isolates. This in general could imply diversities in terms of genetics, evolution, and protein production. Therefore, generic administration of antibiotics such as clarithromycin, amoxicillin, and erythromycin as treatment methods in the African subregion could be contributing to the spread of the bacterium's antibiotic resistance.

Plant-derived exosomes: a green approach for cancer drug delivery

Plant-derived exosomes (PDEs) are natural extracellular vesicles (EVs). In the current decade, they have been highlighted for cancer therapeutic development. Cancer is a global health crisis and it requires an effective, affordable, and less side effect-based treatment. Emerging research based on PDEs suggests that they have immense potential to be considered as a therapeutic option. Research evidences indicate that PDEs' internal molecular cargos show impressive cancer prevention activity with less toxicity. PDEs-based drug delivery systems overcome several limitations of traditional drug delivery tools. Extraction of PDEs from plant sources employ diverse methodologies, encompassing ultracentrifugation, immunoaffinity, size-based isolation, and precipitation, each with distinct advantages and limitations. The core constituents of PDEs comprise of lipids, proteins, DNA, and RNA. Worldwide, a few clinical trials on plant-derived exosomes are underway, and regulatory affairs for their use as therapeutic agents are still not understood with clarity. This review aims to comprehensively analyze the current state of research on plant-derived exosomes as a promising avenue for drug delivery, highlighting anticancer activity, challenges, and future orientation in effective cancer therapeutic development.

Revitalizing agriculture: next-generation genotyping and -omics technologies enabling molecular prediction of resilient traits in the Solanaceae family

This review highlights -omics research in Solanaceae family, with a particular focus on resilient traits. Extensive research has enriched our understanding of Solanaceae genomics and genetics, with historical varietal development mainly focusing on disease resistance and cultivar improvement but shifting the emphasis towards unveiling resilience mechanisms in genebank-preserved germplasm is nowadays crucial. Collecting such information, might help researchers and breeders developing new experimental design, providing an overview of the state of the art of the most advanced approaches for the identification of the genetic elements laying behind resilience. Building this starting point, we aim at providing a useful tool for tackling the global agricultural resilience goals in these crops.

Assessment of Kaistella jeonii esterase conformational dynamics in response to poly(ethylene terephthalate) binding

The pervasive presence of plastic in the environment has reached a concerning scale, being identified in many ecosystems. Bioremediation is the cheapest and most eco-friendly alternative to remove this polymer from affected areas. Recent work described that a novel cold-active esterase enzyme extracted from the bacteria Kaistella jeonii could promiscuously degrade PET. Compared to the well-known PETase from Ideonella sakaiensis, this novel esterase presents a low sequence identity yet has a remarkably similar folding. However, enzymatic assays demonstrated a lower catalytic efficiency. In this work, we employed a strict computational approach to investigate the binding mechanism between the esterase and PET. Understanding the underlying mechanism of binding can shed light on the evolutive mechanism of how enzymes have been evolving to degrade these artificial molecules and help develop rational engineering approaches to improve PETase-like enzymes. Our results indicate that this esterase misses a disulfide bridge, keeping the catalytic residues closer and possibly influencing its catalytic efficiency. Moreover, we describe the structural response to the interaction between enzyme and PET, indicating local and global effects. Our results aid in deepening the knowledge behind the mechanism of biological catalysis of PET degradation and as a base for the engineering of novel PETases.

Genetic improvement in edible fish: status, constraints, and prospects on CRISPR-based genome engineering

Conventional selective breeding in aquaculture has been effective in genetically enhancing economic traits like growth and disease resistance. However, its advances are restricted by heritability, the extended period required to produce a strain with desirable traits, and the necessity to target multiple characteristics simultaneously in the breeding programs. Genome editing tools like zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) are promising for faster genetic improvement in fishes. CRISPR/Cas9 technology is the least expensive, most precise, and well compatible with multiplexing of all genome editing approaches, making it a productive and highly targeted approach for developing customized fish strains with specified characteristics. As a result, the use of CRISPR/Cas9 technology in aquaculture is rapidly growing, with the main traits researched being reproduction and development, growth, pigmentation, disease resistance, trans-GFP utilization, and omega-3 metabolism. However, technological obstacles, such as off-target effects, ancestral genome duplication, and mosaicism in founder population, need to be addressed to achieve sustainable fish production. Furthermore, present regulatory and risk assessment frameworks are inadequate to address the technical hurdles of CRISPR/Cas9, even though public and regulatory approval is critical to commercializing novel technology products. In this review, we examine the potential of CRISPR/Cas9 technology for the genetic improvement of edible fish, the technical, ethical, and socio-economic challenges to using it in fish species, and its future scope for sustainable fish production.

Different transcription of novel, functional long non-coding RNA genes by UV-B in green algae, Volvox carteri

Long non-coding RNAs (lncRNAs) are identified as important regulatory molecules related to diverse biological processes. In recent years, benefiting from the rapid development of high-throughput sequencing technology, RNA-seq, and analysis methods, more lncRNAs have been identified and discovered in various plant and algal species. However, so far, only limited studies related to algal lncRNAs are available. Volvox carteri f. nagariensis is the best multicellular model organism to study in developmental and evolutionary biology; therefore, studying and increasing information about this species is important. This study identified lncRNAs in the multicellular green algae Volvox carteri and 1457 lncRNAs were reported, using RNA-seq data and with the help of bioinformatics tools and software. This study investigated the effect of low-dose UV-B radiation on changes in the expression profile of lncRNAs in gonidial and somatic cells. The differential expression of lncRNAs was analyzed between the treatment (UV-B) and the control (WL) groups in gonidial and somatic cells. A total of 37 and 26 lncRNAs with significant differential expression in gonidial and somatic cells, respectively, were reported. Co-expression analysis between the lncRNAs and their neighbor protein-coding genes (in the interval of ± 10 Kb) was accomplished. In gonidial cells, 184 genes with a positive correlation and 13 genes with a negative correlation (greater than 0.95), and in somatic cells, 174 genes with a positive correlation, and 18 genes with a negative correlation were detected. Functional analysis of neighboring coding genes was also performed based on gene ontology. The results of the current work may help gain deeper insight into the regulation of gene expression in the studied model organism, Volvox carteri.

Reference genes for Eucalyptus spp. under Beauveria bassiana inoculation and subsequently infestation by the galling wasp Leptocybe invasa

Relative gene expression analysis through RT-qPCR is an important molecular technique that helps understanding different molecular mechanisms, such as the plant defense response to insect pests. However, the use of RT-qPCR for gene expression analysis can be affected by factors that directly affect the reliability of the results. Among these factors, the appropriate choice of reference genes is crucial and can strongly impact RT-qPCR relative gene expression analyses, highlighting the importance in correctly choosing the most suitable genes for the success of the analysis. Thus, this study aimed to select and validate reference genes for relative gene expression studies through RT-qPCR in hybrids of Eucalyptus tereticornis × Eucalyptus camaldulensis (drought tolerant and susceptible to Leptocybe invasa) under conditions of inoculation by the Beauveria bassiana fungus and subsequent infestation by L. invasa. The expression level and stability of eleven candidate genes were evaluated. Stability was analyzed using the RefFinder tool, which integrates the geNorm, NormFinder, BestKeeper, and Delta-Ct algorithms. The selected reference genes were validated through the expression analysis of the transcriptional factor EcDREB2 (dehydration-responsive element-binding protein 2). For all treatments evaluated, EcPTB, EcPP2A-1, and EcEUC12 were the best reference genes. The triplets EcPTB/EcEUC12/EcUBP6, EcPP2A-1/EcEUC12/EcPTB, EcIDH/EcSAND/Ecα-TUB, EcPP2A-1/Ecα-TUB/EcPTB, and EcPP2A-1/EcUPL7/EcSAND were the best reference genes for the control plants, mother plants, plants inoculated with B. bassiana, plants infested with L. invasa, and plants inoculated with B. bassiana and subsequently infested with L. invasa, respectively. The best determined reference genes were used to normalize the RT-qPCR expression data for each experimental condition evaluated. The results emphasize the importance of this type of study to ensure the reliability of relative gene expression analyses. Furthermore, the findings of this study can be used as a basis for future research, comprising gene expression analysis of different eucalyptus metabolic pathways.

Unveiling the Genetic Landscape of Feed Efficiency in Holstein Dairy Cows: Insights into Heritability, Genetic Markers, and Pathways via Meta-Analysis

Improving the feeding efficiency of dairy cows is a key component to improve the utilization of land resources and meet the demand for high-quality protein. Advances in genomic methods and omics techniques have made it possible to breed more efficient dairy cows through genomic selection. The aim of this review is to obtain a comprehensive understanding of the biological background of feed efficiency (FE) complex traits in purebred Holstein dairy cows including heritability estimate, and genetic markers, genes, and pathways participating in FE regulation mechanism. Through a literature search, we systematically reviewed the heritability estimation, molecular genetic markers, genes, biomarkers, and pathways of traits related to feeding efficiency in Holstein dairy cows. A meta-analysis based on a random-effects model was performed to combine reported heritability estimates of FE complex. The heritability of residual feed intake, dry matter intake, and energy balance was 0.20, 0.34, and 0.22, respectively, which proved that it was reasonable to include the related traits in the selection breeding program. For molecular genetic markers, a total of 13 single-nucleotide polymorphisms and copy number variance loci, associated genes, and functions were reported to be significant across populations. A total of 169 reported candidate genes were summarized on a large scale, using a higher threshold (adjusted P value < 0.05). Then, the subsequent pathway enrichment of these genes was performed. The important genes reported in the articles were included in a gene list and the gene list was enriched by gene ontology (GO):biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis. Three GO:BP terms and four KEGG terms were statistically significant, which mainly focused on adenosine triphosphate (ATP) synthesis, electron transport chain, and OXPHOS pathway. Among these pathways, involved genes such as ATP5MC2, NDUFA, COX7A2, UQCR, and MMP are particularly important as they were previously reported. Twenty-nine reported biological mechanisms along with involved genes were explained mainly by four biological pathways (insulin-like growth factor axis, lipid metabolism, oxidative phosphorylation pathways, tryptophan metabolism). The information from this study will be useful for future studies of genomic selection breeding and genetic structures influencing animal FE. A better understanding of the underlying biological mechanisms would be beneficial, particularly as it might address genetic antagonism.

Sodium butyrate administration improves intestinal development of suckling lambs

This study was conducted to investigate the effects of sodium butyrate (SB) supplementation on growth performance, intestinal barrier functions, and intestinal bacterial communities in sucking lambs. Forty lambs of 7 d old, with an average body weight (BW) of 4.46 ± 0.45 kg, were allocated into the control (CON) or SB group, with each group having five replicate pens (n = 5). Lambs were orally administered SB at 1.8 mL/kg BW in the SB group or the same volume of saline in the CON group. Treatments were administered from 7 to 35 d of age, when one lamb from each replicate was slaughtered to obtain intestinal tissues and contents. The results showed that supplementation with SB tended to increase the BW (P = 0.079) and the starter intake (P = 0.089) of lambs at 35 d of age. The average daily gain of lambs in the SB group was significantly greater than that in the CON group (P < 0.05). The villus height of jejunum in the SB group was markedly higher (P < 0.05) than that in the CON group. In ileum, lambs in the SB group had lower (P < 0.05) crypt depth and greater (P < 0.05) villus-to-crypt ratio than those in the CON group. Compared with the CON group, the mRNA and protein expressions of Claudin-1 and Occludin were increased (P < 0.05) in the SB group. Supplementation with SB decreased the relative abundances of pathogenic bacteria, including Clostridia_UCG-014 (P = 0.094) and Romboutsia (P < 0.05), which were negatively associated with the intestinal barrier function genes (P < 0.05). The relative abundance of Succiniclasticum (P < 0.05) was higher in the SB group, and it was positively correlated with the ratio of villi height to crypt depth in the jejunum (P < 0.05). Compared with the CON group, the function "Metabolism of Cofactors and Vitamins" was increased in the SB group lambs (P < 0.05). In conclusion, SB orally administration during suckling period could improve the small intestine development and growth performance of lambs by inhibiting the harmful bacteria (Clostridia_UCG-014, Romboutsia) colonization, and enhancing intestinal barrier functions.

Transcriptional profiling upon T cell stimulation reveals down-regulation of inflammatory pathways in T and B cells in SLE versus Sjögren's syndrome

Systemic lupus erythematosus (SLE) and primary Sjögren's syndrome (pSS) share clinical as well as pathogenic similarities. Although previous studies suggest various abnormalities in different immune cell compartments, dedicated cell-type specific transcriptomic signatures are often masked by patient heterogeneity. Here, we performed transcriptional profiling of isolated CD4, CD8, CD16 and CD19 lymphocytes from pSS and SLE patients upon T cell stimulation, in addition to a steady-state condition directly after blood drawing, in total comprising 581 sequencing samples. T cell stimulation, which induced a pronounced inflammatory response in all four cell types, gave rise to substantial re-modulation of lymphocyte subsets in the two autoimmune diseases compared to healthy controls, far exceeding the transcriptomic differences detected at steady-state. In particular, we detected cell-type and disease-specific down-regulation of a range of pro-inflammatory cytokine and chemokine pathways. Such differences between SLE and pSS patients are instrumental for selective immune targeting by future therapies.

PCMT1 knockdown attenuates malignant properties by globally regulating transcriptome profiles in triple-negative breast cancer cells

Background

As the most frequently diagnosed cancer in women, Breast cancer has high mortality and metastasis rate, especially triple-negative breast cancer (TNBC). As an oncogene, protein-L-isoaspartate (D-aspartate) O-methyltransferase (PCMT1) is a prognostic biomarker in breast cancer and is highly expressed, while its underlying functions remain unknown.

Methods

In this study, we silenced PCTM1 in TNBC MDA-MB-231 cells by short hairpin RNA (shPCMT1) to investigate its cellular functions using cell proliferation, apoptosis, migration, and invasion experiments. Following this, the transcriptome sequencing (RNA-seq) experiment was conducted to explore the molecular targets of PCMT1, including differentially expressed genes (DEGs) and regulated alternative splicing events (RASEs).

Results

The results showed that shPCMT1 inhibited the proliferation, migration, and invasion of MDA-MB-231 cells. We obtained 1,084 DEGs and 2,287 RASEs between shPCMT1 and negative control (NC) groups through RNA-seq. The DEGs were significantly enriched in immune or inflammation response and cell adhesion-associated pathways, pathways associated with PCMT1 cellular function in cell migration. The RASE genes were enriched in cell cycle-associated pathways and were associated with the altered cell proliferation rate. We finally validated the changed expression and splicing levels of DEGs and RASEs. We found that 34 RNA binding protein (RBP) genes were dysregulated by shPCMT1, including NQO1, S100A4, EEF1A2, and RBMS2. The dysregulated RBP genes could partially explain how PCMT1 regulates the global transcriptome profiles.

Conclusion

In conclusion, our study identified the molecular targets of PCMT1 in the TNBC cell line, expands our understanding of the regulatory mechanisms of PCMT1 in cancer progression, and provides novel insights into the progression of TNBC. The identified molecular targets are potential therapeutic targets for future TNBC treatment.

Genomics Reveals Complex Population History and Unexpected Diversity of Eurasian Otters (Lutra lutra) in Britain Relative to Genetic Methods

Conservation genetic analyses of many endangered species have been based on genotyping of microsatellite loci and sequencing of short fragments of mtDNA. The increase in power and resolution afforded by whole genome approaches may challenge conclusions made on limited numbers of loci and maternally inherited haploid markers. Here, we provide a matched comparison of whole genome sequencing versus microsatellite and control region (CR) genotyping for Eurasian otters (Lutra lutra). Previous work identified four genetically differentiated "stronghold" populations of otter in Britain, derived from regional populations that survived the population crash of the 1950s-1980s. Using whole genome resequencing data from 45 samples from across the British stronghold populations, we confirmed some aspects of population structure derived from previous marker-driven studies. Importantly, we showed that genomic signals of the population crash bottlenecks matched evidence from otter population surveys. Unexpectedly, two strongly divergent mitochondrial lineages were identified that were undetectable using CR fragments, and otters in the east of England were genetically distinct and surprisingly variable. We hypothesize that this previously unsuspected variability may derive from past releases of Eurasian otters from other, non-British source populations in England around the time of the population bottleneck. Our work highlights that even reasonably well-studied species may harbor genetic surprises, if studied using modern high-throughput sequencing methods.

Differential Gene Expression Responses to Salt and Drought Stress in Tall Fescue (Festuca arundinacea Schreb.)

Understanding gene expression kinetics and the underlying physiological mechanisms in stress combinations is a challenge for the purpose of stress resistance breeding. The novelty of this study is correlating the physiological mechanisms with the expression of key target genes in tall fescue under a combination of various salinity and osmotic stress treatments. Four drought- and salt-responsive genes belonging to different crucial pathways evaluated included one transcription factor FabZIP69, one for the cytosolic polyamine synthetase FaADC1, one for ABA signaling FaCYP707A1, and another one for the specific Na+/H+ plasma membrane antiporter FaSOS1 involve in osmotic homeostasis. FaSOS1, FaCYP707A1, and FabZIP69 were induced early at 6 h after NaCl treatment, while FaSOS1 and FaCYP707A1 were transcribed gradually after exposure to PEG. However, stress interactions showed a significantly increased expression in all genes. Expression of these genes was positively correlated to Pro, SSs, IL, DPPH, and antioxidant enzyme activity and negatively correlated with RWC, total Chl, and MSI. Chemical analyses showed that tall fescue plants exposed to the combination of stresses exhibited increased quantity of reactive oxygen species (H2O2), EL and DPPH, and higher levels of antioxidant enzyme activities (CAT, and SOD), Pro, and SSs content, compared with control seedlings. Under dual-stress conditions, the expression of FabZIP69 was effective in controlling the expression of FaSOS1 and FaADC1 genes differently.

Associating somatic mutation with clinical outcomes through kernel regression and optimal transport

Somatic mutations in cancer patients are inherently sparse and potentially high dimensional. Cancer patients may share the same set of deregulated biological processes perturbed by different sets of somatically mutated genes. Therefore, when assessing the associations between somatic mutations and clinical outcomes, gene-by-gene analysis is often under-powered because it does not capture the complex disease mechanisms shared across cancer patients. Rather than testing genes one by one, an intuitive approach is to aggregate somatic mutation data of multiple genes to assess their joint association with clinical outcomes. The challenge is how to aggregate such information. Building on the optimal transport method, we propose a principled approach to estimate the similarity of somatic mutation profiles of multiple genes between tumor samples, while accounting for gene-gene similarities defined by gene annotations or empirical mutational patterns. Using such similarities, we can assess the associations between somatic mutations and clinical outcomes by kernel regression. We have applied our method to analyze somatic mutation data of 17 cancer types and identified at least five cancer types, where somatic mutations are associated with overall survival, progression-free interval, or cytolytic activity.

Classification of Isatis indigotica Fortune and Isatis tinctoria Linnaeus via comparative analysis of chloroplast genomes

BACKGROUND

Isatis tinctoria Linnaeus and Isatis indigotica Fortune are very inconsistent in their morphological characteristics, but the Flora of China treats them as the same species. In this work, a new technology that differs from conventional barcodes is developed to prove that they are different species and to clarify their classification.

RESULTS AND METHODS

I. indigotica was indistinguishable from I. tinctoria when using ITS2. CPGAVAS2 was used to construct the chloroplast genomes. MAFFT and DnaSP were used to calculate nucleotide polymorphism, the chloroplast genomes of the two have high diversity in the rpl32 ~ trnL-UAG short region. When using this region as a mini barcode, it was found that there are obvious differences in the base numbers of I. tinctoria and different ploidy I. indigotica were found, but diploid and tetraploid I. indigotica had the same number of bases. Moreover, the reconstruction of the maximum likelihood (ML) tree, utilizing the mini-barcode, demonstrated that I. tinctoria and both diploid and tetraploid I. indigotica are located on distinct branches. The genome size of tetraploid I. indigotica was approximately 643.773 MB, the heterozygosity rate was approximately 0.98%, and the repeat sequence content was approximately 90.43%. This species has a highly heterozygous, extremely repetitive genome.

CONCLUSION

A new method was established to differentiate between I. indigotica and I. tinctoria. Furthermore, this approach provides a reference and basis for the directional breeding of Isatis.

Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes

Fruiting bodies (sporocarps, sporophores or basidiomata) of mushroom-forming fungi (Agaricomycetes) are among the most complex structures produced by fungi. Unlike vegetative hyphae, fruiting bodies grow determinately and follow a genetically encoded developmental program that orchestrates their growth, tissue differentiation and sexual sporulation. In spite of more than a century of research, our understanding of the molecular details of fruiting body morphogenesis is still limited and a general synthesis on the genetics of this complex process is lacking. In this paper, we aim at a comprehensive identification of conserved genes related to fruiting body morphogenesis and distil novel functional hypotheses for functionally poorly characterised ones. As a result of this analysis, we report 921 conserved developmentally expressed gene families, only a few dozens of which have previously been reported to be involved in fruiting body development. Based on literature data, conserved expression patterns and functional annotations, we provide hypotheses on the potential role of these gene families in fruiting body development, yielding the most complete description of molecular processes in fruiting body morphogenesis to date. We discuss genes related to the initiation of fruiting, differentiation, growth, cell surface and cell wall, defence, transcriptional regulation as well as signal transduction. Based on these data we derive a general model of fruiting body development, which includes an early, proliferative phase that is mostly concerned with laying out the mushroom body plan (via cell division and differentiation), and a second phase of growth via cell expansion as well as meiotic events and sporulation. Altogether, our discussions cover 1 480 genes of Coprinopsis cinerea, and their orthologs in Agaricus bisporus, Cyclocybe aegerita, Armillaria ostoyae, Auriculariopsis ampla, Laccaria bicolor, Lentinula edodes, Lentinus tigrinus, Mycena kentingensis, Phanerochaete chrysosporium, Pleurotus ostreatus, and Schizophyllum commune, providing functional hypotheses for ~10 % of genes in the genomes of these species. Although experimental evidence for the role of these genes will need to be established in the future, our data provide a roadmap for guiding functional analyses of fruiting related genes in the Agaricomycetes. We anticipate that the gene compendium presented here, combined with developments in functional genomics approaches will contribute to uncovering the genetic bases of one of the most spectacular multicellular developmental processes in fungi. Citation: Nagy LG, Vonk PJ, Künzler M, Földi C, Virágh M, Ohm RA, Hennicke F, Bálint B, Csernetics Á, Hegedüs B, Hou Z, Liu XB, Nan S, M. Pareek M, Sahu N, Szathmári B, Varga T, Wu W, Yang X, Merényi Z (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of Agaricomycetes. Studies in Mycology 104: 1-85. doi: 10.3114/sim.2022.104.01.

Identification of candidate lethal haplotypes and genomic association with post-natal mortality and reproductive traits in Nellore cattle

The wide use of genomic information has enabled the identification of lethal recessive alleles that are the major genetic causes of reduced conception rates, longer calving intervals, or lower survival for live-born animals. This study was carried out to screen the Nellore cattle genome for lethal recessive haplotypes based on deviation from the expected population homozygosity, and to test SNP markers surrounding the lethal haplotypes region for association with heifer rebreeding (HR), post-natal mortality (PNM) and stayability (STAY). This approach requires genotypes only from apparently normal individuals and not from affected embryos. A total of 62,022 animals were genotyped and imputed to a high-density panel (777,962 SNP markers). Expected numbers of homozygous individuals were calculated, and the probabilities of observing 0 homozygotes was obtained. Deregressed genomic breeding values [(G)EBVs] were used in a GWAS to identify candidate genes and biological mechanisms affecting HR, STAY and PNM. In the functional analyses, genes within 100 kb down and upstream of each significant SNP marker, were researched. Thirty haplotypes had high expected frequency, while no homozygotes were observed. Most of the alleles present in these haplotypes had a negative mean effect for PNM, HR and STAY. The GWAS revealed significant SNP markers involved in different physiological mechanisms, leading to harmful effect on the three traits. The functional analysis revealed 26 genes enriched for 19 GO terms. Most of the GO terms found for biological processes, molecular functions and pathways were related to tissue development and the immune system. More phenotypes underlying these putative regions in this population could be the subject of future investigation. Tests to find putative lethal haplotype carriers could help breeders to eliminate them from the population or manage matings in order to avoid homozygous.

The dissection of tomato flavor: biochemistry, genetics, and omics

Flavor and quality are the major drivers of fruit consumption in the US. However, the poor flavor of modern commercial tomato varieties is a major cause of consumer dissatisfaction. Studies in flavor research have informed the role of volatile organic compounds in improving overall liking and sweetness of tomatoes. These studies have utilized and applied the tools of molecular biology, genetics, biochemistry, omics, machine learning, and gene editing to elucidate the compounds and biochemical pathways essential for good tasting fruit. Here, we discuss the progress in identifying the biosynthetic pathways and chemical modifications of important tomato volatile compounds. We also summarize the advances in developing highly flavorful tomato varieties and future steps toward developing a "perfect tomato".

Mapping of quantitative trait loci for the nutritional value of fresh market tomato

The incidence of many diseases, such as cancer, cardiovascular diseases, and diabetes, is associated with malnutrition and an unbalanced daily diet. Vegetables are an important source of vitamins and essential compounds for human health. As a result, such metabolites have increasingly become the focus of breeding programs. Tomato is one of the most popular components of our daily diet. Therefore, the improvement of tomato's nutritional quality is an important goal. In the present study, we performed targeted metabolic profiling of an interspecific Solanum pimpinellifolium × S. lycopersicum inbred backcross line (IBL) population and identified quantitative trait loci responsible for the nutritional value of tomato. Transgressive segregation was apparent for many of the nutritional compounds such that some IBLs had extremely high levels of various amino acids and vitamins compared to their parents. A total of 117 QTLs for nutritional traits including 62 QTLs for amino acids, 18 QTLs for fatty acids, 12 QTLs for water-soluble vitamins, and 25 QTLs for fat-soluble vitamins were identified. Moreover, almost 24% of identified QTLs were confirmed in previous studies, and 40 possible gene candidates were found for 18 identified QTLs. These findings can help breeders to improve the nutritional value of tomato.

Investigating bone resorption in Atlantic herring fish intermuscular bones with solid-state NMR

Bones are connective tissues mainly made of collagen proteins with calcium phosphate deposits. They undergo constant remodeling, including destroying existing bones tissues (known as bone resorption) and rebuilding new ones. Bone remodeling has been well-described in mammals, but it is not the case in fish. Here, we focused on the mobile phase of the bone vascular system by carefully preserving moisture in adult Atlantic herring intermuscular bones. We detected pore water with high ionic strength and soluble degraded peptides whose ¹H-transverse relaxation times, T₂s, exceed 15 milliseconds. With favorable T₂s, we incorporated a solution state spinlock scheme into the INEPT techniques to unequivocally demonstrate collagen degradation. In addition, we detected a substantial amount of inorganic phosphate in solution with ³¹P-NMR in the considerable background of solid hydroxyapatite calcium phosphate by saturation recovery experiment. It is consistent with the idea that bone resorption degrades bone collagen and releases calcium ions and phosphate ions in the pore water with increased ionic strength. Our report is the first to probe the resorption process in the heterogenous bone microstructure with a rigorous characterization of ¹H and ¹³C relaxation behavior and direct assignments. In addition, we contribute to the fish bones literature by investigating fish bone remodeling using NMR for the first time.

Characterization and phylogenetic analysis of the chloroplast genome of galium spurium

Galium spurium is a farmland weed, with strong stress resistance. However, its chloroplast genome has never been reported. In this study, the complete sequence of the chloroplast genome of G. spurium was characterized, which is a circular molecule, 153,481 bp in length, and with a large single copy region of 84,334 bp, a small single copy region of 17,057 bp, and a pair of inverted repeat regions of 26,045 bp. The whole genome contained 127 genes, including 82 protein-coding genes, 37 transfer RNA genes, and eight ribosomal RNA genes. Phylogenetic analysis shows that it relates closely to G. aparine. This study provides a basis for the further phylogenic study of Galium.

QTL mapping and identification of candidate genes linked to red rot resistance in sugarcane

Sugarcane (Saccharum species hybrid) is one of the most important commercial crops cultivated worldwide for products like white sugar, bagasse, ethanol, etc. Red rot is a major sugarcane disease caused by a hemi-biotrophic fungus, Colletotrichum falcatum Went., which can potentially cause a reduction in yield up to 100%. Breeding for red rot-resistant sugarcane varieties has become cumbersome due to its complex genome and frequent generation of new pathotypes of red rot fungus. In the present study, a genetic linkage map was developed using a selfed population of a popular sugarcane variety CoS 96268. A QTL linked to red rot resistance (qREDROT) was identified, which explained 26% of the total phenotypic variation for the trait. A genotype-phenotype network analysis performed to account for epistatic interactions, identified the key markers involved in red rot resistance. The differential expression of the genes located in the genomic region between the two flanking markers of the qREDROT as well as in the vicinity of the markers identified through the genotype-phenotype network analysis in a set of contrasting genotypes for red rot infection further confirmed the mapping results. Further, the expression analysis revealed that the plant defense-related gene coding 26S protease regulatory subunit is strongly associated with the red rot resistance. The findings can help in the screening of disease resistant genotypes for developing red rot-resistant varieties of sugarcane.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03481-7.

Plasticity of maternal environment-dependent expression-QTLs of tomato seeds

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

Full-length transcriptome from different life stages of cobia (Rachycentron canadum, Rachycentridae)

Cobia (Rachycentron canadum, Rachycentridae) is one of the prospective species for mariculture. The transcriptome-based study on cobia was hampered by an inadequate reference genome and a lack of full-length cDNAs. We used a long-read based sequencing technology (PacBio Sequel II Iso-Seq3 SMRT) to obtain complete transcriptome sequences from larvae, juveniles, and various tissues of adult cobia, and a single SMRTcell generated 99 gigabytes of data and 51,205,946,694 bases. A total of 8609435, 7441673 and 9140164 subreads were generated from the larval, juvenile, and adult sample pools, with mean sub-read lengths of 2109.9, 1988.2 and 1996.2 bp, respectively. All samples were combined to increase transcript recovery and clustered into 35661 high-quality reads. This is the first report on a full-length transcriptome from R. canadum. Our results illustrate a significant increase in the identified amount of cobia LncRNAs and alternatively spliced transcripts, which will help improve genome annotation. Furthermore, this information will be beneficial for nutrigenomics and functional studies on cobia and other commercially important mariculture species.

Integrated single-molecule real-time sequencing and RNA sequencing reveal the molecular mechanisms of salt tolerance in a novel synthesized polyploid genetic bridge between maize and its wild relatives

BACKGROUND

Tripsacum dactyloides (2n = 4x = 72) and Zea perennis (2n = 4x = 40) are tertiary gene pools of Zea mays L. and exhibit many abiotic adaptations absent in modern maize, especially salt tolerance. A previously reported allopolyploid (hereafter referred to as MTP, 2n = 74) synthesized using Zea mays, Tripsacum dactyloides, and Zea perennis has even stronger salt tolerance than Z. perennis and T. dactyloides. This allopolyploid will be a powerful genetic bridge for the genetic improvement of maize. However, the molecular mechanisms underlying its salt tolerance, as well as the key genes involved in regulating its salt tolerance, remain unclear.

RESULTS

Single-molecule real-time sequencing and RNA sequencing were used to identify the genes involved in salt tolerance and reveal the underlying molecular mechanisms. Based on the SMRT-seq results, we obtained 227,375 reference unigenes with an average length of 2300 bp; most of the unigenes were annotated to Z. mays sequences (76.5%) in the NR database. Moreover, a total of 484 and 1053 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Functional enrichment analysis of DEGs revealed that multiple pathways responded to salt stress, including "Flavonoid biosynthesis," "Oxidoreductase activity," and "Plant hormone signal transduction" in the leaves and roots, and "Iron ion binding," "Acetyl-CoA carboxylase activity," and "Serine-type carboxypeptidase activity" in the roots. Transcription factors, such as those in the WRKY, B3-ARF, and bHLH families, and cytokinin negatively regulators negatively regulated the salt stress response. According to the results of the short time series-expression miner analysis, proteins involved in "Spliceosome" and "MAPK signal pathway" dynamically responded to salt stress as salinity changed. Protein-protein interaction analysis revealed that heat shock proteins play a role in the large interaction network regulating salt tolerance.

CONCLUSIONS

Our results reveal the molecular mechanism underlying the regulation of MTP in the response to salt stress and abundant salt-tolerance-related unigenes. These findings will aid the retrieval of lost alleles in modern maize and provide a new approach for using T. dactyloides and Z. perennis to improve maize.

Superoxide dismutase (SOD) family in durum wheat: promising candidates for improving crop resilience

The SOD family has been extensively analyzed at genome wide level in several crops. However, little is known about this family in durum wheat. In this study, a total of 14 TdSOD genes were identified in whole durum wheat genome including 8 TdCu-ZnSODs, 2 TdMnSODs, and 4 TdFeSODs. In silico analysis evinced that TdSOD family members displayed a closer evolutionary relationship, similar gene structure and protein features with their homologs from other plant species. Furthermore, the analysis of their promoter regions revealed the presence of a great number of cis-regulatory elements related to plant development, abiotic and biotic stresses, phytohormones, and several potential binding sites for transcription factors. Interestingly, 3D structure analysis revealed that TdCu-ZnSOD2A-2 and TdCu-ZnSOD2B-2, belonging to the Cu-Zn group, were modeled as copper chaperone for SOD like their homologs from rice and Arabidopsis. The expression profile of eight TdSOD candidate genes was investigated under salt, drought, cold, and ABA treatments. Notably, TdCu-ZnSOD2A-1, TdFeSOD4A-1, and TdFeSOD7A-1 were significantly up-regulated under all stress treatments. On the other hand, TdCu-ZnSOD7B and TdMnSOD2B were strongly expressed in roots and leaves under cold stress and TdCu-ZnSOD2B-2 was particularly up-regulated in leaves under ABA treatment. Ultimately, these findings provide valuable information for the identification of attractive candidate genes to improve wheat resilience.

Periodontitis: An Oral Disease with Severe Consequences

Periodontitis, being a multifactorial disorder is found to be the most common oral disease denoted by the inflammation of gingiva and resorption of tooth supporting alveolar bone. The disease being closely linked with fast life style and determined by unhygienic behavioural factors, the internal milieu of oral cavity and formation of plaque biofilm on the dental and gingival surfaces. Porphyromonas gingivalis, being the major keystone pathogen of the periodontal biofilm evokes host immune responses that causes damage of gingival tissues and resorption of bones. The biofilm associated microbial community progressively aggravates the condition resulting in chronic inflammation and finally tooth loss. The disease often maintains bidirectional relationship with different systemic, genetic, autoimmune, immunodeficiency diseases and even psychological disorders. The disease can be diagnosed and predicted by various genetic, radiographic and computer-aided design (CAD) & computer-aided engineering (CAE) and artificial neural network (ANN). The elucidation of genetic background explains the inheritance of the disease. The therapeutic approaches commonly followed include mechanical removal of dental plaque with the use of systemic antibiotics. Awareness generation amongst local people, adoption of good practice of timely tooth brushing preferably with fluoride paste or with nanoconjugate pastes will reduce the chance of periodontal plaque formation. Modern tissue engineering technology like 3D bioprinting of periodontal tissue may help in patient specific flawless regeneration of tooth structures and associated bones.

Whole-Genome Resequencing Deciphers New Insight Into Genetic Diversity and Signatures of Resistance in Cultivated Cotton Gossypium hirsutum

Cotton is an important crop that produces fiber and cottonseed oil for the textile and oil industry. However, cotton leaf curl virus disease (CLCuD) stress is limiting its yield in several Asian countries. In this study, we have sequenced Mac7 accession, a Gossypium hirsutum resistance source against several biotic stresses. By aligning with the Gossypium hirsutum (AD1) 'TM-1' genome, a total of 4.7 and 1.2 million SNPs and InDels were identified in the Mac7 genome. The gene ontology and metabolic pathway enrichment indicated SNPs and InDels role in nucleotide bindings, secondary metabolite synthesis, and plant-pathogen interaction pathways. Furthermore, the RNA-seq data in different tissues and qPCR expression profiling under CLCuD provided individual gene roles in resistant and susceptible accessions. Interestingly, the differential NLR genes demonstrated higher expression in resistant plants rather than in susceptible plants expression. The current resequencing results may provide primary data to identify DNA resistance markers which will be helpful in marker-assisted breeding for development of Mac7-derived resistance lines.

Strain-specific effect of Limosilactobacillus fermentum with distinct genetic lineages on loperamide-induced constipation in mice: attributing effects to certain genes

In 2013, Limosilactobacillus fermentum was regarded as a "generally recognized as safe" organism by the US Food and Drug Administration, and emerging evidence showed that it can exert beneficial health effects on humans. In this study, five L. fermentum strains from different phylogroups of a phylogenetic tree containing 224 L. fermentum strains were chosen, and their protective effects against loperamide-induced constipation in mice were studied. Animal experiments showed that L. fermentum YN54 significantly alleviated weight loss, increased fecal moisture, accelerated intestinal peristalsis, and increased the small intestinal transit rate in mice with constipation by regulating gastrointestinal peptides and increasing the amount of intestinal short-chain fatty acids. However, the other four L. fermentum strains (XJ61, CECT5716, WX115, and GD121) did not relieve constipation in mice treated with loperamide. A comparative genomic analysis of these strains was conducted and "L. fermentum YN54 only" genes were functionally annotated and validated with the other three L. fermentum strains (FJ12, GX51, and ZH1010) that had different functional genes. Finally, the genes involved in the synthesis of fatty acid hydrase, polysaccharides, and cell membranes were identified to be associated with the probiotic effect of L. fermentum on mice with constipation through preliminary experiments in this study.

Multiple instance neural networks based on sparse attention for cancer detection using T-cell receptor sequences

Early detection of cancers has been much explored due to its paramount importance in biomedical fields. Among different types of data used to answer this biological question, studies based on T cell receptors (TCRs) are under recent spotlight due to the growing appreciation of the roles of the host immunity system in tumor biology. However, the one-to-many correspondence between a patient and multiple TCR sequences hinders researchers from simply adopting classical statistical/machine learning methods. There were recent attempts to model this type of data in the context of multiple instance learning (MIL). Despite the novel application of MIL to cancer detection using TCR sequences and the demonstrated adequate performance in several tumor types, there is still room for improvement, especially for certain cancer types. Furthermore, explainable neural network models are not fully investigated for this application. In this article, we propose multiple instance neural networks based on sparse attention (MINN-SA) to enhance the performance in cancer detection and explainability. The sparse attention structure drops out uninformative instances in each bag, achieving both interpretability and better predictive performance in combination with the skip connection. Our experiments show that MINN-SA yields the highest area under the ROC curve scores on average measured across 10 different types of cancers, compared to existing MIL approaches. Moreover, we observe from the estimated attentions that MINN-SA can identify the TCRs that are specific for tumor antigens in the same T cell repertoire.

Towards a potential pan-cancer prognostic signature for gene expression based on probesets and ensemble machine learning

Cancer is one of the leading causes of death worldwide and can be caused by environmental aspects (for example, exposure to asbestos), by human behavior (such as smoking), or by genetic factors. To understand which genes might be involved in patients’ survival, researchers have invented prognostic genetic signatures: lists of genes that can be used in scientific analyses to predict if a patient will survive or not. In this study, we joined together five different prognostic signatures, each of them related to a specific cancer type, to generate a unique pan-cancer prognostic signature, that contains 207 unique probesets related to 187 unique gene symbols, with one particular probeset present in two cancer type-specific signatures (203072_at related to the MYO1E gene). We applied our proposed pan-cancer signature with the Random Forests machine learning method to 57 microarray gene expression datasets of 12 different cancer types, and analyzed the results. We also compared the performance of our pan-cancer signature with the performances of two alternative prognostic signatures, and with the performances of each cancer type-specific signature on their corresponding cancer type-specific datasets. Our results confirmed the effectiveness of our prognostic pan-cancer signature. Moreover, we performed a pathway enrichment analysis, which indicated an association between the signature genes and a protein-protein interaction analysis, that highlighted PIK3R2 and FN1 as key genes having a fundamental relevance in our signature, suggesting an important role in pan-cancer prognosis for both of them.

A comprehensive Bioconductor ecosystem for the design of CRISPR guide RNAs across nucleases and technologies

The success of CRISPR-mediated gene perturbation studies is highly dependent on the quality of gRNAs, and several tools have been developed to enable optimal gRNA design. However, these tools are not all adaptable to the latest CRISPR modalities or nucleases, nor do they offer comprehensive annotation methods for advanced CRISPR applications. Here, we present a new ecosystem of R packages, called crisprVerse, that enables efficient gRNA design and annotation for a multitude of CRISPR technologies. This includes CRISPR knockout (CRISPRko), CRISPR activation (CRISPRa), CRISPR interference (CRISPRi), CRISPR base editing (CRISPRbe) and CRISPR knockdown (CRISPRkd). The core package, crisprDesign, offers a user-friendly and unified interface to add off-target annotations, rich gene and SNP annotations, and on- and off-target activity scores. These functionalities are enabled for any RNA- or DNA-targeting nucleases, including Cas9, Cas12, and Cas13. The crisprVerse ecosystem is open-source and deployed through the Bioconductor project ( https://github.com/crisprVerse ).

Genetic Differentiation and Migration Fluxes of Viruses from Melon Crops and Crop Edge Weeds

Weeds surrounding crops may act as alternative hosts, playing important epidemiological roles as virus reservoirs and impacting virus evolution. We used high-throughput sequencing to identify viruses in Spanish melon crops and plants belonging to three pluriannual weed species, Ecballium elaterium, Malva sylvestris, and Solanum nigrum, sampled at the edges of the crops. Melon and E. elaterium, both belonging to the family Cucurbitaceae, shared three virus species, whereas there was no virus species overlap between melon and the other two weeds. The diversity of cucurbit aphid-borne yellows virus (CABYV) and tomato leaf curl New Delhi virus (ToLCNDV), both in melon and E. elaterium, was further studied by amplicon sequencing. Phylogenetic and population genetics analyses showed that the CABYV population was structured by the host, identifying three sites in the CABYV RNA-dependent RNA polymerase under positive selection, perhaps reflecting host adaptation. The ToLCNDV population was much less diverse than the CABYV one, likely as a consequence of the relatively recent introduction of ToLCNDV in Spain. In spite of its low diversity, we identified geographical but no host differentiation for ToLCNDV. Potential virus migration fluxes between E. elaterium and melon plants were also analyzed. For CABYV, no evidence of migration between the populations of the two hosts was found, whereas important fluxes were identified between geographically distant subpopulations for each host. For ToLCNDV, in contrast, evidence of migration from melon to E. elaterium was found, but not the other way around. IMPORTANCE It has been reported that about half of the emerging diseases affecting plants are caused by viruses. Alternative hosts often play critical roles in virus emergence as virus reservoirs, bridging host species that are otherwise unconnected and/or favoring virus diversification. In spite of this, the viromes of potential alternative hosts remain largely unexplored. In the case of crops, pluriannual weeds at the crop edges may play these roles. Here, we took advantage of the power of high-throughput sequencing to characterize the viromes of three weed species frequently found at the edges of melon crops. We identified three viruses shared by melon and the cucurbit weed, with two of them being epidemiologically relevant for melon crops. Further genetic analyses showed that these two viruses had contrasting patterns of diversification and migration, providing an interesting example on the role that weeds may play in the ecology and evolution of viruses affecting crops.

Host-pathogen coevolution drives innate immune response to Aphanomyces astaci infection in freshwater crayfish: transcriptomic evidence

BACKGROUND

For over a century, scientists have studied host-pathogen interactions between the crayfish plague disease agent Aphanomyces astaci and freshwater crayfish. It has been hypothesised that North American crayfish hosts are disease-resistant due to the long-lasting coevolution with the pathogen. Similarly, the increasing number of latent infections reported in the historically sensitive European crayfish hosts seems to indicate that similar coevolutionary processes are occurring between European crayfish and A. astaci. Our current understanding of these host-pathogen interactions is largely focused on the innate immunity processes in the crayfish haemolymph and cuticle, but the molecular basis of the observed disease-resistance and susceptibility remain unclear. To understand how coevolution is shaping the host's molecular response to the pathogen, susceptible native European noble crayfish and invasive disease-resistant marbled crayfish were challenged with two A. astaci strains of different origin: a haplogroup A strain (introduced to Europe at least 50 years ago, low virulence) and a haplogroup B strain (signal crayfish in lake Tahoe, USA, high virulence). Here, we compare the gene expression profiles of the hepatopancreas, an integrated organ of crayfish immunity and metabolism.

RESULTS

We characterised several novel innate immune-related gene groups in both crayfish species. Across all challenge groups, we detected 412 differentially expressed genes (DEGs) in the noble crayfish, and 257 DEGs in the marbled crayfish. In the noble crayfish, a clear immune response was detected to the haplogroup B strain, but not to the haplogroup A strain. In contrast, in the marbled crayfish we detected an immune response to the haplogroup A strain, but not to the haplogroup B strain.

CONCLUSIONS

We highlight the hepatopancreas as an important hub for the synthesis of immune molecules in the response to A. astaci. A clear distinction between the innate immune response in the marbled crayfish and the noble crayfish is the capability of the marbled crayfish to mobilise a higher variety of innate immune response effectors. With this study we outline that the type and strength of the host immune response to the pathogen is strongly influenced by the coevolutionary history of the crayfish with specific A. astaci strains.

BindVAE: Dirichlet variational autoencoders for de novo motif discovery from accessible chromatin

We present a novel unsupervised deep learning approach called BindVAE, based on Dirichlet variational autoencoders, for jointly decoding multiple TF binding signals from open chromatin regions. BindVAE can disentangle an input DNA sequence into distinct latent factors that encode cell-type specific in vivo binding signals for individual TFs, composite patterns for TFs involved in cooperative binding, and genomic context surrounding the binding sites. On the task of retrieving the motifs of expressed TFs in a given cell type, BindVAE is competitive with existing motif discovery approaches.

Metabolic alterations in alga Chlamydomonas reinhardtii exposed to nTiO2 materials

Nano-sized titanium dioxide (nTiO2) is one of the most commonly used materials, however the knowledge about the molecular basis for metabolic and physiological changes in phytoplankton is yet to be explored. In the present study we use a combination of targeted metabolomics, transcriptomics and physiological response studies to decipher the metabolic perturbation in green alga Chlamydomonas reinhardtii exposed for 72 h to increasing concentrations (2, 20, 100 and 200 mg L-1) of nTiO2 with primary sizes of 5, 15 and 20 nm. Results show that the exposure to all three nTiO2 materials induced perturbation of the metabolism of amino acids, nucleotides, fatty acids, tricarboxylic acids, antioxidants but not in the photosynthesis. The alterations of the most responsive metabolites were concentration and primary size-dependent despite the significant formation of micrometer-size aggregates and their sedimentation. The metabolic perturbations corroborate the observed physiological responses and transcriptomic results and confirmed the importance of oxidative stress as a major toxicity mechanism for nTiO2. Transcriptomics revealed also an important influence of nTiO2 treatments on the transport, adenosine triphosphate binding cassette transporters, and metal transporters, suggesting a perturbation in a global nutrition of the microalgal cell, which was most pronounced for exposure to 5 nm nTiO2. The present study provides for the first-time evidence for the main metabolic perturbations in green alga C. reinhardtii exposed to nTiO2 and helps to improve biological understanding of the molecular basis of these perturbations.

Interrelated Effects of Zinc Deficiency and the Microbiome on Group B Streptococcal Vaginal Colonization

Group B Streptococcus (GBS) in the vaginal tract is a risk factor for preterm birth and adverse pregnancy outcomes. GBS colonization is also transient in nature, which likely reflects the contributions of pathogen determinants, interactions with commensal flora, and host factors, making this environment particularly challenging to understand. Additionally, dietary zinc deficiency is a health concern on the global scale that is known to be associated with recurrent bacterial infection and increased rate of preterm birth or stillbirth. However, the impact of zinc deficiency on vaginal health has not yet been studied. Here we use a murine model to assess the role of dietary zinc on GBS burden and the impact of GBS colonization on the vaginal microbiome. We show that GBS vaginal colonization is increased in a zinc-deficient host and that the presence of GBS significantly alters the microbial community structure of the vagina. Using machine learning approaches, we show that vaginal community turnover during GBS colonization is driven by computationally predictable changes in key taxa, including several organisms not previously described in the context of the vaginal microbiota, such as Akkermansia muciniphila. We observed that A. muciniphila increases GBS vaginal persistence and, in a cohort of human vaginal microbiome samples collected throughout pregnancy, we observed an increased prevalence of codetection of GBS and A. muciniphila in patients who delivered preterm compared to those who delivered at full term. These findings reveal the importance and complexity of both host zinc availability and native microbiome to GBS vaginal persistence. IMPORTANCE The presence of group B Streptococcus (GBS) in the vaginal tract, perturbations in the vaginal microbiota, and dietary zinc deficiency are three factors that are independently known to be associated with increased risk of adverse pregnancy outcomes. Here, we developed an experimental mouse model to assess the impact of dietary zinc deficiency on GBS vaginal burden and persistence and to determine how changes in GBS colonization impact vaginal microbial structure. We have employed unique animal, in silica metabolic, and machine learning models, paired with analyses of human cohort data, to identify taxonomic biomarkers that contribute to host susceptibility to GBS vaginal persistence. Collectively, the data reported here identify that both dietary zinc deficiency and the presence of A. muciniphila could perpetuate an increased GBS burden and prolonged exposure in the vaginal tract, which potentiate the risk of invasive infection in utero and in the newborn.

The chloroplast genome of Farsetia hamiltonii Royle, phylogenetic analysis, and comparative study with other members of Clade C of Brassicaceae

BACKGROUND

Farsetia hamiltonii Royle is a medicinally important annual plant from the Cholistan desert that belongs to the tribe Anastaticeae and clade C of the Brassicaceae family. We provide the entire chloroplast sequence of F.hamiltonii, obtained using the Illumina HiSeq2500 and paired-end sequencing. We compared F. hamiltonii to nine other clade C species, including Farsetia occidentalis, Lobularia libyca, Notoceras bicorne, Parolinia ornata, Morettia canescens, Cochlearia borzaeana, Megacarpaea polyandra, Biscutella laevigata, and Iberis amara. We conducted phylogenetic research on the 22 Brassicaceae species, which included members from 17 tribes and six clades.

RESULTS

The chloroplast genome sequence of F.hamiltonii of 154,802 bp sizes with 36.30% GC content and have a typical structure comprised of a Large Single Copy (LSC) of 83,906 bp, a Small Single Copy (SSC) of 17,988 bp, and two copies of Inverted Repeats (IRs) of 26,454 bp. The genomes of F. hamiltonii and F. occidentalis show shared amino acid frequencies and codon use, RNA editing sites, simple sequence repeats, and oligonucleotide repeats. The maximum likelihood tree revealed Farsetia as a monophyletic genus, closely linked to Morettia, with a bootstrap score of 100. The rate of transversion substitutions (Tv) was higher than the rate of transition substitutions (Ts), resulting in Ts/Tv less than one in all comparisons with F. hamiltonii, indicating that the species are closely related. The rate of synonymous substitutions (Ks) was greater than non-synonymous substitutions (Ka) in all comparisons with F. hamiltonii, with a Ka/Ks ratio smaller than one, indicating that genes underwent purifying selection. Low nucleotide diversity values range from 0.00085 to 0.08516, and IR regions comprise comparable genes on junctions with minimal change, supporting the conserved status of the selected chloroplast genomes of the clade C of the Brassicaceae family. We identified ten polymorphic regions, including rps8-rpl14, rps15-ycf1, ndhG-ndhI, psbK-psbI, ccsA-ndhD, rpl36-rps8, petA-psbJ, ndhF-rpl32, psaJ-rpl3, and ycf1 that might be exploited to construct genuine and inexpensive to solve taxonomic discrepancy and understand phylogenetic relationship amongst Brassicaceae species.

CONCLUSION

The entire chloroplast sequencing of F. hamiltonii sheds light on the divergence of genic chloroplast sequences among members of the clade C. When other Farsetia species are sequenced in the future, the full F. hamiltonii chloroplast will be used as a source for comprehensive taxonomical investigations of the genus. The comparison of F. hamiltonii and other clade C species adds new information to the phylogenetic data and evolutionary processes of the clade. The results of this study will also provide further molecular uses of clade C chloroplasts for possible plant genetic modifications and will help recognise more Brassicaceae family species.

The IMD and Toll canonical immune pathways of Triatoma pallidipennis are preferentially activated by Gram-negative and Gram-positive bacteria, respectively, but cross-activation also occurs

BACKGROUND

Antimicrobial peptides (AMPs) participate in the humoral immune response of insects eliminating invasive microorganisms. The immune deficiency pathway (IMD) and Toll are the main pathways by which the synthesis of these molecules is regulated in response to Gram-negative (IMD pathway) or Gram-positive (Toll pathway) bacteria. Various pattern-recognition receptors (PRRs) participate in the recognition of microorganisms, such as pgrp-lc and toll, which trigger signaling cascades and activate NF-κB family transcription factors, such as relish, that translocate to the cell nucleus, mainly in the fat body, inducing AMP gene transcription.

METHODS

T. pallidipennis inhibited in Tppgrp-lc, Tptoll, and Tprelish were challenged with E. coli and M. luteus to analyze the expression of AMPs transcripts in the fat body and to execute survival assays.

RESULTS

In this work we investigated the participation of the pgrp-lc and toll receptor genes and the relish transcription factor (designated as Tppgrp-lc, Tptoll, and Tprelish), in the transcriptional regulation of defensin B, prolixicin, and lysozyme B in Triatoma pallidipennis, one of the main vectors of Chagas disease. AMP transcript abundance was higher in the fat body of blood-fed than non-fed bugs. Challenge with Escherichia coli or Micrococcus luteus induced differential increases in AMP transcripts. Additionally, silencing of Tppgrp-lc, Tptoll, and Tprelish resulted in reduced AMP transcription and survival of bugs after a bacterial challenge.

CONCLUSIONS

Our findings demonstrated that the IMD and Toll pathways in T. pallidipennis preferentially respond to Gram-negative and Gram-positive bacteria, respectively, by increasing the expression of AMP transcripts, but cross-induction also occurs.

Chromosome-level genome assembly of the diploid blueberry Vaccinium darrowii provides insights into its subtropical adaptation and cuticle synthesis

Vaccinium darrowii is a subtropical wild blueberry species that has been used to breed economically important southern highbush cultivars. The adaptive traits of V. darrowii to subtropical climates can provide valuable information for breeding blueberry and perhaps other plants, especially against the background of global warming. Here, we assembled the V. darrowii genome into 12 pseudochromosomes using Oxford Nanopore long reads complemented with Hi-C scaffolding technologies, and we predicted 41 815 genes using RNA-sequencing evidence. Syntenic analysis across three Vaccinium species revealed a highly conserved genome structure, with the highest collinearity between V. darrowii and Vaccinium corymbosum. This conserved genome structure may explain the high fertility observed during crossbreeding of V. darrowii with other blueberry cultivars. Analysis of gene expansion and tandem duplication indicated possible roles for defense- and flowering-associated genes in the adaptation of V. darrowii to the subtropics. Putative SOC1 genes in V. darrowii were identified based on phylogeny and expression analysis. Blueberries are covered in a thick cuticle layer and contain anthocyanins, which confer their powdery blue color. Using RNA sequencing, we delineated the cuticle biosynthesis pathways of Vaccinium species in V. darrowii. This result can serve as a reference for breeding berries whose colors are appealing to customers. The V. darrowii reference genome, together with the unique traits of this species, including its diploid genome, short vegetative phase, and high compatibility in hybridization with other blueberries, make V. darrowii a potential research model for blueberry species.

Aldehyde dehydrogenase superfamily in sorghum: genome-wide identification, evolution, and transcript profiling during development stages and stress conditions

BACKGROUND

Aldehyde dehydrogenases (ALDHs) are a family of NAD(P)⁺ dependent enzymes that detoxify aldehydes by promoting their oxidation to respective carboxylic acids. The role of ALDH enzymes in various plant species has been extensively studied, revealing their critical role in salinity, drought, heat, and heavy metal stress tolerance. Despite their physiological significance, ALDH genes in Sorghum bicolor have yet to be studied thoroughly.

RESULTS

In this study, a total of 19 ALDH genes have been identified that have been grouped into ten families based on the criteria of the ALDH gene nomenclature committee. Segmental duplication assisted more in the enhancement of SbALDH gene family members than tandem duplication. All the identified SbALDH members made a cluster with monocot rice and maize in the phylogenetic tree rather than dicot species, suggesting the pre-eudicot-monocot separation of the ALDH superfamily members. The gene structure and protein domain were found to be mostly conserved in separate phylogenetic classes, indicating that each family played an important role in evolution. Expression analysis revealed that several SbALDHs were expressed in various tissues, developmental stages, and in response to abiotic stresses, indicating that they can play roles in plant growth, development, or stress adaptation. Interestingly, the majority of the SbALDH genes were found to be highly responsive to drought stress, and the SbALDH18B1 transcript showed maximum enhancement in all the stress conditions. The presence of cis-acting elements (mainly ABRE and MBS) in the promoter region of these genes might have a significant role in drought tolerance.

CONCLUSIONS

Our findings add to the current understanding, evolutionary history, and contribution of SbALDHs in stress tolerance, and smooth the path of further functional validation of these genes.

Genome-wide association study reveals novel quantitative trait loci and candidate genes of lint percentage in upland cotton based on the CottonSNP80K array

Thirty-four SNPs corresponding with 22 QTLs for lint percentage, including 13 novel QTLs, was detected via GWAS. Two candidate genes underlying this trait were also identified. Cotton (Gossypium spp.) is an important natural textile fiber and oilseed crop cultivated worldwide. Lint percentage (LP, %) is one of the important yield components, and increasing LP is a core goal of cotton breeding improvement. However, the genetic and molecular mechanisms underlying LP in upland cotton remain unclear. Here, we performed a genome-wide association study (GWAS) for LP based on 254 upland cotton accessions in four environments as well as the best linear unbiased predictors using the high-density CottonSNP80K array. In total, 41,413 high-quality single-nucleotide polymorphisms (SNPs) were screened, and 34 SNPs within 22 quantitative trait loci (QTLs) were significantly associated with LP. In total, 175 candidate genes were identified from two major genomic loci (GR1 and GR2), and 50 hub genes were identified through GO enrichment and weighted gene co-expression network analysis. Two candidate genes (Gh_D01G0162 and Gh_D07G0463), which may participate in early fiber development to affect the number of fiber protrusions and LP, were also identified. Their genetic variation and expression were verified by linkage disequilibrium blocks, haplotypes, and quantitative real-time polymerase chain reaction, respectively. The weighted gene interaction network analysis showed that the expression of Gh_D07G0463 was significantly correlated with that of Gh_D01G0162. These identified SNPs, QTLs and candidate genes provide important insights into the genetic and molecular mechanisms underlying variations in LP and serve as a foundation for LP improvement via marker-assisted breeding.

Impact of 5HydroxyMethylCytosine (5hmC) on reverse/direct association of cell-cycle, apoptosis, and extracellular matrix pathways in gastrointestinal cancers

BACKGROUND

Aberrant levels of 5-hydroxymethylcytosine (5-hmC) can lead to cancer progression. Identification of 5-hmC-related biological pathways in cancer studies can produce better understanding of gastrointestinal (GI) cancers. We conducted a network-based analysis on 5-hmC levels extracted from circulating free DNAs (cfDNA) in GI cancers including colon, gastric, and pancreatic cancers, and from healthy donors. The co-5-hmC network was reconstructed using the weighted-gene co-expression network method. The cancer-related modules/subnetworks were detected. Preservation of three detected 5-hmC-related modules was assessed in an external dataset. The 5-hmC-related modules were functionally enriched, and biological pathways were identified. The relationship between modules was assessed using the Pearson correlation coefficient (p-value < 0.05). An elastic network classifier was used to assess the potential of the 5-hmC modules in distinguishing cancer patients from healthy individuals. To assess the efficiency of the model, the Area Under the Curve (AUC) was computed using five-fold cross-validation in an external dataset.

RESULTS

The main biological pathways were the cell cycle, apoptosis, and extracellular matrix (ECM) organization. Direct association between the cell cycle and apoptosis, inverse association between apoptosis and ECM organization, and inverse association between the cell cycle and ECM organization were detected for the 5-hmC modules in GI cancers. An AUC of 92% (0.73-1.00) was observed for the predictive model including 11 genes.

CONCLUSION

The intricate association between biological pathways of identified modules may reveal the hidden significance of 5-hmC in GI cancers. The identified predictive model and new biomarkers may be beneficial in cancer detection and precision medicine using liquid biopsy in the early stages.