Search

Abstract The intricate nature of microbiota sequencing data—high dimensionality and sparsity—presents a challenge in identifying informative and reproducible microbial features for both research and clinical applications. Addressing this, we introduce PreLect, an innovative feature selection framework that harnesses microbes’ prevalence to facilitate consistent selection in sparse microbiota data. Upon rigorous benchmarking against established feature selection methodologies across 42 microbiome datasets, PreLect demonstrated superior classification capabilities compared to statistical methods and outperformed machine learning-based methods by selecting features with greater prevalence and abundance. A significant strength of PreLect lies in its ability to reliably identify reproducible microbial features across varied cohorts. Applied to colorectal cancer, PreLect identifies key microbes and highlights crucial pathways, such as lipopolysaccharide and glycerophospholipid biosynthesis, in cancer progression. This case study exemplifies PreLect’s utility in discerning clinically relevant microbial signatures. In summary, PreLect’s accuracy and robustness make it a significant advancement in the analysis of complex microbiota data.

Abstract Background The nuclear factor kappa B (NFκB) regulatory pathways downstream of tumor necrosis factor (TNF) play a critical role in carcinogenesis. However, the widespread influence of NFκB in cells can result in off-target effects, making it a challenging therapeutic target. Ensemble learning is a machine learning technique where multiple models are combined to improve the performance and robustness of the prediction. Accordingly, an ensemble learning model could uncover more precise targets within the NFκB/TNF signaling pathway for cancer therapy. Methods In this study, we trained an ensemble learning model on the transcriptome profiles from 16 cancer types in the TCGA database to identify a robust set of genes that are consistently associated with the NFκB/TNF pathway in cancer. Our model uses cancer patients as features to predict the genes involved in the NFκB/TNF signaling pathway and can be adapted to predict the genes for different cancer types by switching the cancer type of patients. We also performed functional analysis, survival analysis, and a case study of triple-negative breast cancer to demonstrate our model's potential in translational cancer medicine. Results Our model accurately identified genes regulated by NFκB in response to TNF in cancer patients. The downstream analysis showed that the identified genes are typically involved in the canonical NFκB-regulated pathways, particularly in adaptive immunity, anti-apoptosis, and cellular response to cytokine stimuli. These genes were found to have oncogenic properties and detrimental effects on patient survival. Our model also could distinguish patients with a specific cancer subtype, triple-negative breast cancer (TNBC), which is known to be influenced by NFκB-regulated pathways downstream of TNF. Furthermore, a functional module known as mononuclear cell differentiation was identified that accurately predicts TNBC patients and poor short-term survival in non-TNBC patients, providing a potential avenue for developing precision medicine for cancer subtypes. Conclusions In conclusion, our approach enables the discovery of genes in NFκB-regulated pathways in response to TNF and their relevance to carcinogenesis. We successfully categorized these genes into functional groups, providing valuable insights for discovering more precise and targeted cancer therapeutics.

Abstract The 15q13.3 microdeletion has pleiotropic effects ranging from apparently healthy to severely affected individuals. The underlying basis of the variable phenotype remains elusive. We analyzed gene expression using blood from three individuals with 15q13.3 microdeletion and brain cortex tissue from ten mice Df[h15q13]/+. We assessed differentially expressed genes (DEGs), protein–protein interaction (PPI) functional modules, and gene expression in brain developmental stages. The deleted genes’ haploinsufficiency was not transcriptionally compensated, suggesting a dosage effect may contribute to the pathomechanism. DEGs shared between tested individuals and a corresponding mouse model show a significant overlap including genes involved in monogenic neurodevelopmental disorders. Yet, network-wide dysregulatory effects suggest the phenotype is not caused by a single critical gene. A significant proportion of blood DEGs, silenced in adult brain, have maximum expression during the prenatal brain development. Based on DEGs and their PPI partners we identified altered functional modules related to developmental processes, including nervous system development. We show that the 15q13.3 microdeletion has a ubiquitous impact on the transcriptome pattern, especially dysregulation of genes involved in brain development. The high phenotypic variability seen in 15q13.3 microdeletion could stem from an increased vulnerability during brain development, instead of a specific pathomechanism.

Abstract Background Recently, non-coding RNAs are of growing interest, and more scientists attach importance to research on their functions. Long non-coding RNAs (lncRNAs) are defined as non-protein coding transcripts longer than 200 nucleotides. We already knew that lncRNAs are related to cancers and will be dysregulated in them. But most of their functions are still left to further study. A mechanism of RNA regulation, known as competing endogenous RNAs (ceRNAs), has been proposed to explain the complex relationships among mRNAs and lncRNAs by competing for binding with shared microRNAs (miRNAs). Methods We proposed an analysis framework to construct the association networks among lncRNA, mRNA, and miRNAs based on their expression patterns and decipher their network modules. Results We collected a large-scale gene expression dataset of 1,061 samples from breast invasive carcinoma (BRCA) patients, each consisted of the expression profiles of 4,359 lncRNAs, 16,517 mRNAs, and 534 miRNAs, and applied the proposed analysis approach to interrogate them. We have uncovered the underlying ceRNA modules and the key modulatory lncRNAs for different subtypes of breast cancer. Conclusions We proposed a modulatory analysis to infer the ceRNA effects among mRNAs and lncRNAs and performed functional analysis to reveal the plausible mechanisms of lncRNA modulation in the four breast cancer subtypes. Our results might provide new directions for breast cancer therapeutics and the proposed method could be readily applied to other diseases.

IntroductionThe microbiota-gut-brain axis is implicated in Alzheimer’s disease. Gut microbiota alterations in mild cognitive impairment (MCI) are inconsistent and remain to be understood. This study aims to investigate the gut microbial composition associated with MCI, cognitive functions, and structural brain differences.MethodsA nested case-control study was conducted in a community-based prospective cohort where detailed cognitive functions and structural brain images were collected. Thirty-one individuals with MCI were matched to sixty-five cognitively normal controls by age strata, gender, and urban/rural area. Fecal samples were examined using 16S ribosomal RNA (rRNA) V3–V4 sequencing. Compositional differences between the two groups were identified and correlated with the cognitive functions and volumes/thickness of brain structures.ResultsThere was no significant difference in alpha and beta diversity between MCIs and cognitively normal older adults. However, the abundance of the genus Ruminococcus, Butyricimonas, and Oxalobacter decreased in MCI patients, while an increased abundance of nine other genera, such as Flavonifractor, were found in MCIs. Altered genera discriminated MCI patients well from controls (AUC = 84.0%) and were associated with attention and executive function.ConclusionThis study provides insights into the role of gut microbiota in the neurodegenerative process.

Abstract Background RNA-seq emerges as a valuable method for clinical genetics. The transcriptome is “dynamic” and tissue-specific, but typically the probed tissues to analyze (TA) are different from the tissue of interest (TI) based on pathophysiology. Results We developed Phenotype-Tissue Expression and Exploration (PTEE), a tool to facilitate the decision about the most suitable TA for RNA-seq. We integrated phenotype-annotated genes, used 54 tissues from GTEx to perform correlation analyses and identify expressed genes and transcripts between TAs and TIs. We identified skeletal muscle as the most appropriate TA to inquire for cardiac arrhythmia genes and skin as a good proxy to study neurodevelopmental disorders. We also explored RNA-seq limitations and show that on-off switching of gene expression during ontogenesis or circadian rhythm can cause blind spots for RNA-seq-based analyses. Conclusions PTEE aids the identification of tissues suitable for RNA-seq for a given pathology to increase the success rate of diagnosis and gene discovery. PTEE is freely available at https://bioinf.eva.mpg.de/PTEE/

MicroRNAs (miRNAs) are approximately 20–22 nucleotides in length, which are well known to participate in the post-transcriptional modification. The mature miRNAs were observed to be varied on 5′ or 3′ that raise another term—the isoforms of mature miRNAs (isomiRs), which have been proven not the artifacts and discussed widely recently. In our research, we focused on studying the 5′ isomiRs in lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA). We characterized 75 isomiRs significantly associated with better prognosis and 43 isomiRs with poor prognosis. The 75 protective isomiRs can successfully distinguish tumors from normal samples and are expressed differently between patients of early and late stages. We also found that most of the protective isomiRs tend to be with downstream shift and upregulated compared with those with upstream shift, implying that a possible selection occurs during cancer development. Among these protective isomiRs, we observed a highly positive and significant correlation, as well as in harmful isomiRs, suggesting cooperation within the group. However, between protective and harmful, there is no such a concordance but conversely more negative correlation, suggesting the possible antagonistic effect between protective and harmful isomiRs. We also identified that two isomiRs miR-181a-3p|-3 and miR-181a-3p|2, respectively, belong to the harmful and protective groups, suggesting a bidirectional regulation of their originated archetype—miR-181a-3p. Additionally, we found that the protective isomiRs of miR-21-5p, which is an oncomiR, may be evolved as the tumor suppressors through producing isomiRs to hinder metastasis. In summary, these results displayed the characteristics of the protective isomiRs and their potential for developing the treatment of lung cancer.

Cancer development and progression result from somatic evolution by an accumulation of genomic alterations. The effects of those alterations on the fitness of somatic cells lead to evolutionary adaptations such as increased cell proliferation, angiogenesis, and altered anticancer drug responses. However, there are few general mathematical models to quantitatively examine how perturbations of a single gene shape subsequent evolution of the cancer genome. In this study, we proposed the gene gravity model to study the evolution of cancer genomes by incorporating the genome-wide transcription and somatic mutation profiles of ~3,000 tumors across 9 cancer types from The Cancer Genome Atlas into a broad gene network. We found that somatic mutations of a cancer driver gene may drive cancer genome evolution by inducing mutations in other genes. This functional consequence is often generated by the combined effect of genetic and epigenetic (e.g., chromatin regulation) alterations. By quantifying cancer genome evolution using the gene gravity model, we identified six putative cancer genes (AHNAK, COL11A1, DDX3X, FAT4, STAG2, and SYNE1). The tumor genomes harboring the nonsynonymous somatic mutations in these genes had a higher mutation density at the genome level compared to the wild-type groups. Furthermore, we provided statistical evidence that hypermutation of cancer driver genes on inactive X chromosomes is a general feature in female cancer genomes. In summary, this study sheds light on the functional consequences and evolutionary characteristics of somatic mutations during tumorigenesis by propelling adaptive cancer genome evolution, which would provide new perspectives for cancer research and therapeutics.

MicroRNAs are small non-coding RNAs that can regulate expressions of their target genes at the post-transcriptional level. In this study, we propose a tri-component strategy that combines the conservation of microRNAs, homology of mRNA coding regions, and conserved microRNA binding sites in the 3' untranslated regions to discover conserved microRNA-mRNA interactions. To validate the performance of our conservation strategy, we collected the experimentally validated microRNA-mRNA interactions from three databases as the golden standard. We found that the proposed strategy can improve the performance of existing target prediction algorithms by approximately 2-4 fold. In addition, we demonstrated that the proposed strategy could efficiently retain highly confident interactions from the intersection results of the existing algorithms and filter out the possible false positive predictions in the union one. Furthermore, this strategy can facilitate our ability to trace the homologues in different species that are targeted by the same miRNA family because it combines these three features to identify the conserved miRNA-mRNA interactions during evolution. Through an extensive application of the proposed conservation strategy to a study of the miR-1/206 regulatory network, we demonstrate that the target mRNA recruiting process could be associated with expansion of miRNA family during its evolution. We also uncovered the functional evolution of the miR-1/206 regulatory network. In this network, the early targeted genes tend to participate in more general and development-related functions. In summary, the conservation strategy is capable of helping to highlight the highly confident miRNA-mRNA interactions and can be further applied to reveal the evolutionary features of miRNA regulatory network and functions.

Alternative polyadenylation (APA) could result in mRNA isoforms with variable lengths of 3' UTRs. Gain of microRNA target sites in the 3' UTR of a long mRNA isoform may cause different regulation from the corresponding short isoform. It has been known that cancer cells globally exhibit a lower ratio of long and short isoforms (LSR); that is, they tend to express larger amounts of short isoforms. The objective of this study is to illustrate the relationship between microRNA differential regulation and LSR. We retrieved public APA annotations and isoform expression profiles of breast cancer and normal cells from a high-throughput sequencing method study specific for the mRNA 3' end. Combining microRNA expression profiles, we performed statistical analysis to reveal and estimate microRNA regulation on APA patterns in a global scale. First, we found that the amount of microRNA target sites in the alternative UTR (aUTR), the region only present in long isoforms, could affect the LSR of the target genes. Second, we observed that the genes whose aUTRs were targeted by up-regulated microRNAs in cancer cells had an overall lower LSR. Furthermore, the target sites of up-regulated microRNAs tended to appear in aUTRs. Finally, we demonstrated that the amount of target sites for up-regulated microRNAs in aUTRs correlated with the LSR change between cancer and normal cells. The results indicate that up-regulation of microRNAs might cause lower LSRs of target genes in cancer cells through degradation of their long isoforms. Our findings provide evidence of how microRNAs might play a crucial role in APA pattern shifts from normal to cancerous or proliferative states.

Many biological networks are signed molecular networks which consist of positive and negative links. To reveal the distinct features between links with different signs, we proposed signed link-clustering coefficients that assess the similarity of inter-action profiles between linked molecules. We found that positive links tended to cluster together, while negative links usually behaved like bridges between positive clusters. Positive links with higher adhesiveness tended to share protein domains, be associated with protein-protein interactions and make intra-connections within protein complexes. Negative links that were more bridge-like tended to make interconnections between protein complexes. Utilizing the proposed measures to group positive links, we observed hierarchical modules that could be well characterized by functional annotations or known protein complexes. Our results imply that the proposed sign-specific measures can help reveal the network structural characteristics and the embedded biological contexts of signed links, as well as the functional organization of signed molecular networks.

Ultraviolet (UV) irradiation has emerged as a viable alternative for water/wastewater disinfection. Laboratory dose-response data from collimated-beam tests are commonly used as a basis for determining the necessary delivered UV dose for full-scale UV systems as measured by UV intensity and exposure time. While researchers often think that germicidal dose can be estimated confidently in a collimated beam system, the reported dose-response relationships vary considerably. Numerous factors may affect the test results to some extent. They include apparatus setup, column dimensions, UV lamp type and output, intensity measurement, shutter type and operation, petri dish specifications, sample volume and depth of the liquid, mixing condition, laboratory settings, microbial organism preparation and testing, and water quality. The methodology used to calculate the UV dose for collimated beam tests is also a critical factor. To ensure reproducibility or to have a meaningful comparison of results from different collimated beam tests, a standardized protocol for collimated beam testing and its dose calculation is necessary. The standardized collimated beam test can then serve as a tool to 'calibrate' the dose calculation models for field-scale UV reactors and, consequently, to meaningfully compare the dose values derived from these models. A comprehensive standardized collimated beam testing protocol is warranted to advance applications of UV disinfection for water/ wastewater. However, development of such a protocol is a significant research effort that will require input from researchers of various technical and disciplinary backgrounds. In addition, some technical issues still need to be clarified by further research. This paper presents important issues associated with collimated beam testing in the context of developing a standard protocol as well as future research needs. The main objective is to serve as a starting point for development of a standardized collimated beam testing protocol that is universally acceptable by researchers, industries, and regulatory agencies. CE Database subject headings: Disinfection; Ultraviolet radiation; Tests.

Abstract Background Orbital infiltration in patients with multiple myeloma is a rare condition, with less than 50 cases reported in the medical literature. Most patients undergo conservative treatment because multiple myeloma is a disseminated systemic disease. Case presentation A 43-year-old male subject with multiple myeloma and long-term survival presented with orbital involvement. The subject lacked the typical features and poor prognostic factors associated with multiple myeloma, such as renal failure, hypercalcemia, and paraprotein in the serum and urine. The orbital computed tomographic scan revealed the tumor encasing the optic nerve, but without prominent bony destruction. Therefore, a frontal skull craniotomy with an epidural entrance to the orbital space was performed, to completely extirpate the orbital mass. The surgical procedure was followed by moderate-dose radiation therapy. After 32 months of follow-up care, the subject is doing well with excellent local control. Conclusion Although the effectiveness and applicability of this approach remains to be determined, this case report demonstrates that accurate and early detection combined with local surgical treatment and appropriate radio/chemotherapy, can be applied to effectively extend an orbital multiple myeloma patient's life.