Functional Heatmap: an automated and interactive pattern recognition tool to integrate time with multi-omics assays

The Role of Alternative Splicing Factors hnRNP G and Fox-2 in the Progression and Prognosis of Esophageal Cancer

Aim. Alternative splicing (AS) has been widely demonstrated in the occurrence and progression of many cancers. Nevertheless, the involvement of cancer-associated splicing factors in the development of esophageal carcinoma (ESCA) remains to be explored. Method. RNA-Seq data and the corresponding clinical information of the ESCA cohort were downloaded from The Cancer Genome Atlas database. Bioinformatics methods were used to further analyzed the differently expressed AS (DEAS) events and their splicing network. Kaplan–Meier, Cox regression, and unsupervised cluster analyses were used to assess the association between AS events and clinical characteristics of ESCA patients. The splicing factors screened out were verified in vitro at the cellular level. Results. A total of 50,342 AS events were identified, of which 3,988 were DEAS events and 46 of these were associated with overall survival (OS) of ESCA patients, with a 5-year OS rate of 0.941. By constructing a network of AS events with survival-related splicing factors, the AS factors related to prognosis can be further identified. In vitro experiments and database analysis confirmed that the high expression of hnRNP G in ESCA is related to the high invasion ability of ESCA cells and the poor prognosis of ESCA patients. In contrast, the low expression of fox-2 in esophageal cancer is related to a better prognosis. Conclusion. ESCA-associated AS factors hnRNP G and Fox-2 are of great value in deciphering the underlying mechanisms of AS in ESCA and providing clues for therapeutic goals for further validation.

Genome-wide characterization of C2H2 zinc-finger gene family provides insight into the mechanisms and evolution of the dehydration-rehydration responses in Physcomitrium and Arabidopsis

Dehydration tolerance is a vital factor for land plant evolution and world agricultural production. Numerous studies enlightened that the plant-specific C2H2-type zinc-finger proteins (C2H2-ZFPs) as master regulators played pivotal roles in the abiotic stress responses of plants. However, a comprehensive understanding of the evolution of C2H2-ZFPs in terrestrial plants and its regulatory mechanism in dehydration and rehydration response remains a mystery. In this study, the genome-wide identification of C2H2-ZFP genes revealed 549 homologs in the representatives of terrestrial plant lineages from liverwort to angiosperms. Based on the characteristics of the conserved C2H2-ZF domains, four major C2H2-ZF types (M-, Z-, Q-, and D-type) were identified in the C2H2-ZFPs, with the dominants of M-type in all selected species and followed by Z-type in non-seed plants and Q-type in seed plants, respectively. Phylogenetic analyses of the identified C2H2-ZFPs supported four major groups in the land plant representatives, among which the members from the desiccation-tolerant Physcomitrium patens and the dehydration-sensitive Arabidopsis thaliana displayed different topological relationships in the phylogenies reconstructed for a single species. C2H2-ZFPs clustered in the same subclades shared similar features in their conserved domains and gene structures. Approximately, 81% of the C2H2-ZFP promoters of all 549 identified C2H2-ZFPs harbored the conserved ABA-responsive elements (ABREs) and/or dehydration-responsive elements (DREs). Comparative transcriptomic analyses showed that 50 PpZFPs and 56 AtZFPs significantly changed their transcripts abundance. Interestingly, most of the dehydration- and rehydration-responsive PpZPFs and AtZFPs had been predicted to contain the ABRE and DRE elements in their promoter regions and with over half of which phylogenetically belonging to group III. The differences in the expression patterns of C2H2-ZFPs in responses to dehydration and rehydration between P. patens and A. thaliana reflected their different strategies to adapt to dehydration. The identified candidate PpZFPs were specifically induced by moderate dehydration and reached the peak transcript abundance in severe dehydration. Our study lays the foundations for further functional investigation of C2H2-ZFPs in dehydration responses from an evolutionary perspective in land plants. The findings will provide us with genetic resources and potential targets for drought tolerance breeding in crops and beyond.

Identification of unique transcriptomic signatures and key genes through RNA sequencing and integrated WGCNA and PPI network analysis in HIV infected lung cancer

With the widespread use of highly active antiretroviral therapy (HARRT), the survival time of AIDS patients has been greatly extended. However, the incidence of lung cancer in HIV-infected patients is increasing and has become a major problem threatening the survival of AIDS patients. The aim of this study is to use Weighted Gene Co-expression Network Analysis (WGCNA) and differential gene analysis to find possible key genes involved in HIV-infected lung cancer. In this study, using lung tissue samples from five pairs of HIV-infected lung cancer patients, second-generation sequencing was performed and transcriptomic data were obtained. A total of 132 HIV-infected lung cancer-related genes were screened out by WGCNA and differential gene expression analysis methods. Based on gene annotation analysis, these genes were mainly enriched in mitosis-related functions and pathways. In addition, in protein-protein interaction (PPI) analysis, a total of 39 hub genes were identified. Among them, five genes (ASPM, CDCA8, CENPF, CEP55, and PLK1) were present in both three hub gene lists (intersection gene, DEGs, and WCGNA module) suggesting that these five genes may become key genes involved in HIV-infected lung cancer.

MOVIS: A multi-omics software solution for multi-modal time-series clustering, embedding, and visualizing tasks

Thanks to recent advances in sequencing and computational technologies, many researchers with biological and/or medical backgrounds are now producing multiple data sets with an embedded temporal dimension. Multi-modalities enable researchers to explore and investigate different biological and physico-chemical processes with various technologies. Motivated to explore multi-omics data and time-series multi-omics specifically, the exploration process has been hindered by the separation introduced by each omics-type. To effectively explore such temporal data sets, discover anomalies, find patterns, and better understand their intricacies, expertise in computer science and bioinformatics is required. Here we present MOVIS, a modular time-series multi-omics exploration tool with a user-friendly web interface that facilitates the data exploration of such data. It brings into equal participation each time-series omic-type for analysis and visualization. As of the time of writing, two time-series multi-omics data sets have been integrated and successfully reproduced. The resulting visualizations are task-specific, reproducible, and publication-ready. MOVIS is built on open-source software and is easily extendable to accommodate different analytical tasks. An online version of MOVIS is available under https://movis.mathematik.uni-marburg.de/ and on Docker Hub (https://hub.docker.com/r/aanzel/movis).

Graphene-triphenyl phosphate (TPP) co-exposure in the marine environment: Interference with metabolism and immune regulation in mussel Mytilus galloprovincialis

Both immune regulation and endocrine systems are great challenges to marine organisms, and effective protocols for determining these adverse outcome pathways are limited, especially in vivo. The increasing usage of graphene nanomaterials can lead to the frequent exposure to marine organisms. Triphenyl phosphate (TPP), an organophosphate flame retardant, is frequently detected in natural environments. In this study, the combined toxic effects of co-exposure to graphene and TPP was investigated in Mytilus galloprovincialis using computational toxicology and multi-omics technology. Noticeably, graphene could disturb the membrane stability and increase the tissue accumulation of TPP. The adsorption behavior of TPP on graphene could inhibit the surface activity of graphene. In the digestive gland, transcriptomics analysis revealed the down-regulated genes in graphene + TPP treatment, including glyceraldehyde-3-phosphate dehydrogenase (GAPDH), sorbitol dehydrogenase (SORD), glutathione s-transferase mu 3 (GSTM3) and 4-aminobutyrate aminotransferase (ABAT), were mainly associated with oxidative stress and energy metabolism. Moreover, metabolic responses indicated that graphene + TPP could cause disturbances in energy metabolism and osmotic regulation marked by differentially altered ATP, glucose and taurine in mussels. These data underline the need for further knowledge on the potential interactions of nanomaterials with existing contaminants in marine organisms.

Sensopeptidomic Kinetic Approach Combined with Decision Trees and Random Forests to Study the Bitterness during Enzymatic Hydrolysis Kinetics of Micellar Caseins

Protein hydrolysates are, in general, mixtures of amino acids and small peptides able to supply the body with the constituent elements of proteins in a directly assimilable form. They are therefore characterised as products with high nutritional value. However, hydrolysed proteins display an unpleasant bitter taste and possible off-flavours which limit the field of their nutrition applications. The successful identification and characterisation of bitter protein hydrolysates and, more precisely, the peptides responsible for this unpleasant taste are essential for nutritional research. Due to the large number of peptides generated during hydrolysis, there is an urgent need to develop methods in order to rapidly characterise the bitterness of protein hydrolysates. In this article, two enzymatic hydrolysis kinetics of micellar milk caseins were performed for 9 h. For both kinetics, the optimal time to obtain a hydrolysate with appreciable organoleptic qualities is 5 h. Then, the influence of the presence or absence of peptides and their intensity over time compared to the different sensory characteristics of hydrolysates was studied using heat maps, random forests and regression trees. A total of 22 peptides formed during the enzymatic proteolysis of micellar caseins and influencing the bitterness the most were identified. These methods represent simple and efficient tools to identify the peptides susceptibly responsible for bitterness intensity and predict the main sensory feature of micellar casein enzymatic hydrolysates.

Time-dependent changes in placental mRNA expression after delivery due to delayed specimen collection

PROBLEM

Current scientific guidelines recommend collecting placental specimens within two hours of delivery for gene expression analysis. However, collecting samples in a narrow time window is a challenge in the dynamic and unpredictable clinical setting, so delays in placental specimen collection are possible. The purpose of our analysis was to investigate temporal changes in placental gene expression by longitudinally sampling placentas over a 24 h period.

METHOD OF STUDY

Eight placentas from individuals with uncomplicated, term pregnancies delivered by scheduled cesarean section were collected and sampled following the placental delivery and again at 1, 2, 4, 6, and 24 h post-delivery. At each time point, biopsies of chorionic villous tissue were taken from 3 cotyledons to account for intra-placental heterogeneity. The 3 biopsies from each time point were pooled prior to RNA extraction. Expression of 382 mRNA transcripts was quantified using the NanoString nCounter System. Fold change values were calculated for each time point relative to delivery, and a fold change threshold of 1.25 was used to determine a meaningful change from delivery.

RESULTS

Based on a fold change threshold of 1.25, 84.3% of transcripts were stable for at least 1 h, 80.2% were stable for at least two hours, and 20.6% of transcripts were stable through the collection at 24 h.

CONCLUSION

Our results suggest that for some mRNA transcripts, expression changes as time to sample collection increases. We have developed a Web application to allow investigators to explore transcripts relevant to their research interests and to set appropriate thresholds to aid in determining whether placentas with delayed sample collection can be included in analyses (https://placentaexpression.foundationsofhealth.org/).

The nasal viromes of cattle on arrival at western Canadian feedlots and their relationship to development of bovine respiratory disease

Bovine respiratory disease (BRD) has a complex pathogenesis and aetiology, being the costliest disease affecting the cattle industry in North America. In this study, we applied Nanopore-based viral metagenomic sequencing to explore the nasal virome of cattle upon arrival at feedlot and related the findings to the development of BRD. Deep nasal swabs (DNS) from 310 cattle for which BRD outcomes were known (155 cattle developed BRD within 40 days and 155 remained healthy) were included. The most prevalent virus in on-arrival samples was bovine coronavirus (BCV) (45.2%, 140/310), followed by bovine rhinitis virus B (BRBV) (21.9%, 68/310), enterovirus E (EVE) (19.6%, 60/310), bovine parainfluenza virus 3 (BPIV3) (10.3%, 32/310), ungulate tetraparvovirus 1 (UTPV1) (9.7%, 30/310) and influenza D virus (7.1%, 22/310). No relationship was found between BRD development and the number of viruses detected, the presence of any specific individual virus or combination of viruses. Bovine kobuvirus (BKV) was detected in 2.6% of animals (8/310), being the first report of this virus in Canada. Results of this study demonstrate the diversity of viruses in bovine DNS collected upon arrival at feedlot and highlights the need for further research into prediction of BRD development in the context of mixed infections.

Circulating MicroRNAs from Serum Exosomes May Serve as a Putative Biomarker in the Diagnosis and Treatment of Patients with Focal Cortical Dysplasia

Focal cortical dysplasia (FCD) is a congenital malformation of cortical development where the cortical neurons located in the brain area fail to migrate in the proper formation. Epilepsy, particularly medically refractory epilepsy, is the most common clinical presentation for all types of FCD. This study aimed to explore the expression change of circulating miRNAs in patients with FCD from serum exosomes. A total of nine patients with FCD and four healthy volunteers were enrolled in this study. The serum exosomes were isolated from the peripheral blood of the subjects. Transmission electron microscopy (TEM) was used to identify the exosomes. Both exosomal markers and neuronal markers were detected by Western blotting analysis to prove that we could obtain central nervous system-derived exosomes from the circulation. The expression profiles of circulating exosomal miRNAs were assessed using next-generation sequencing analysis (NGS). We obtained a total of 107 miRNAs with dominant fold change (>2-fold) from both the annotated 5p-arm and 3p-arm of 2780 mature miRNAs. Based on the integrated platform of HMDD v3.2, miRway DB and DIANA-miRPath v3.0 online tools, and confirmed by MiRBase analysis, four potentially predicted miRNAs from serum exosomes in patients with FCD were identified, including miR194-2-5p, miR15a-5p, miR-132-3p, and miR-145-5p. All four miRNAs presented upregulated expression in patients with FCD compared with controls. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and pathway category of four target miRNAs, we found eight possible signaling pathways that may be related to FCD. Among them, we suggest that the mTOR signaling pathway, PI3K-Akt signaling pathway, p53 signaling pathway, and cell cycle regulation and TGF-beta signaling pathway are high-risk pathways that play a crucial role in the pathogenesis of FCD and refractory epilepsy. Our results suggest that the circulating miRNAs from exosomes may provide a potential biomarker for diagnostic, prognostic, and therapeutic adjuncts in patients with FCD and refractory epilepsy.

Eugenol-loaded chitosan emulsion holds the texture of chilled hairtail (Trichiurus lepturus) better: mechanism exploration by proteomic analysis

Eugenol-loaded chitosan emulsion (ELCE) has been proved to have an excellent antibacterial property.

The pulmonary virome, bacteriological and histopathological findings in bovine respiratory disease from western Canada

The aetiology and pathogenesis of bovine respiratory disease (BRD) are complex and involve the interplay of infectious agents, management and environmental factors. Previous studies of BRD focused on ante‐mortem samples from the upper respiratory tract and identified several unconventional viruses. The lung, however, is the primary location where significant BRD lesions are usually found and is a common post‐mortem diagnostic specimen. In this study, results of high‐throughput virome sequencing, bacterial culture, targeted real‐time PCR and histological examination of 130 bovine pneumonic lungs from western Canadian cattle were combined to explore associations of microorganisms with different types of pneumonia. Fibrinous bronchopneumonia (FBP) was the predominant type of pneumonia (46.2%, 60/130) and was associated with the detection of Mannheimia haemolytica. Detection of Histophilus somni and Pasteurella multocida was associated with suppurative bronchopneumonia (SBP) and concurrent bronchopneumonia and bronchointerstitial pneumonia (BP&BIP), respectively. Sixteen viruses were identified, of which bovine parvovirus 2 (BPV2) was the most prevalent (11.5%, 15/130) followed by ungulate tetraparvovirus 1 (UTPV1, 8.5%, 11/130) and bovine respiratory syncytial virus (BRSV, 8.5%, 11/130). None of these viruses, however, were significantly associated with a particular type of pneumonia. Unconventional viruses such as influenza D virus (IDV) and bovine rhinitis B virus (BRBV) were detected, although sparsely, consistent with our previous findings in upper respiratory tract samples. Taken together, our results show that while virus detection in post‐mortem lung samples is of relatively little diagnostic value, the strong associations of H. somni and M. haemolytica with SBP and FBP, respectively, indicate that histopathology can be useful in differentiating bacterial aetiologies.

A-Lister: a tool for analysis of differentially expressed omics entities across multiple pairwise comparisons

BACKGROUND

Researchers commonly analyze lists of differentially expressed entities (DEEs), such as differentially expressed genes (DEGs), differentially expressed proteins (DEPs), and differentially methylated positions/regions (DMPs/DMRs), across multiple pairwise comparisons. Large biological studies can involve multiple conditions, tissues, and timepoints that result in dozens of pairwise comparisons. Manually filtering and comparing lists of DEEs across multiple pairwise comparisons, typically done by writing custom code, is a cumbersome task that can be streamlined and standardized.

RESULTS

A-Lister is a lightweight command line and graphical user interface tool written in Python. It can be executed in a differential expression mode or generic name list mode. In differential expression mode, A-Lister accepts as input delimited text files that are output by differential expression tools such as DESeq2, edgeR, Cuffdiff, and limma. To allow for the most flexibility in input ID types, to avoid database installation requirements, and to allow for secure offline use, A-Lister does not validate or impose restrictions on entity ID names. Users can specify thresholds to filter the input file(s) by column(s) such as p-value, q-value, and fold change. Additionally, users can filter the pairwise comparisons within the input files by fold change direction (sign). Queries composed of intersection, fuzzy intersection, difference, and union set operations can also be performed on any number of pairwise comparisons. Thus, the user can filter and compare any number of pairwise comparisons within a single A-Lister differential expression command. In generic name list mode, A-Lister accepts delimited text files containing lists of names as input. Queries composed of intersection, fuzzy intersection, difference, and union set operations can then be performed across these lists of names.

CONCLUSIONS

A-Lister is a flexible tool that enables the user to rapidly narrow down large lists of DEEs to a small number of most significant entities. These entities can then be further analyzed using visualization, pathway analysis, and other bioinformatics tools.