ggmsa: a visual exploration tool for multiple sequence alignment and associated data.

The identification of the conserved and variable regions in the multiple sequence alignment (MSA) is critical to accelerating the process of understanding the function of genes. MSA visualizations allow us to transform sequence features into understandable visual representations. As the sequence-structure-function relationship gains increasing attention in molecular biology studies, the simple display of nucleotide or protein sequence alignment is not satisfied. A more scalable visualization is required to broaden the scope of sequence investigation. Here we present ggmsa, an R package for mining comprehensive sequence features and integrating the associated data of MSA by a variety of display methods. To uncover sequence conservation patterns, variations and recombination at the site level, sequence bundles, sequence logos, stacked sequence alignment and comparative plots are implemented. ggmsa supports integrating the correlation of MSA sequences and their phenotypes, as well as other traits such as ancestral sequences, molecular structures, molecular functions and expression levels. We also design a new visualization method for genome alignments in multiple alignment format to explore the pattern of within and between species variation. Combining these visual representations with prime knowledge, ggmsa assists researchers in discovering MSA and making decisions. The ggmsa package is open-source software released under the Artistic-2.0 license, and it is freely available on Bioconductor (https://bioconductor.org/packages/ggmsa) and Github (https://github.com/YuLab-SMU/ggmsa).

Chemotherapy plus therapeutic plasmapheresis with 4% human albumin solution in multiple myeloma patients with acute kidney injury: a prospective, open-label, proof-of-concept study.

As no unified treatment protocol or evidence yet exists for plasmapheresis without plasma, this study explored the outcomes of using 4% human albumin (ALB) solution as a replacement solution in patients undergoing plasma exchange for multiple myeloma (MM) patients with acute kidney injury (AKI). This study was prospectively registered (ChiCTR2000030640 and NCT05251896). Bortezomib-based chemotherapy plus therapeutic plasmapheresis (TPP) with 4% human ALB solution was assessed for three years in patients with MM aged >18 years, with AKI according to the Kidney Disease Improving Global Outcomes criteria, and without previous renal impairment from other causes. The primary endpoints were changes in renal function over 18 weeks and survival outcomes at 36 months. The secondary endpoints were the incidence of adverse reactions and symptom improvement. Among the 119 patients included in the analysis, 108 experienced renal reactions. The M protein (absolute changes: median -12.12%, interquartile ranges (IQRs) -18.62 to -5.626) and creatine (median -46.91 µmol/L, IQR -64.70 to -29.12) levels decreased, whereas the estimated glomerular filtration rate (eGFR) increased (median 20.66 mL/(min·1.73 m2), IQR 16.03-25.29). Regarding patient survival, 68.1% and 35.3% of patients survived for >12 and >36 months, respectively. The three symptoms with the greatest relief were urine foam, poor appetite, and blurred vision. All 11 patients (7.6%) who experienced mild adverse reactions achieved remission. In conclusion, in MM patients with AKI, plasma-free plasmapheresis with 4% human ALB solution and bortezomib-based chemotherapy effectively alleviated light chain damage to kidney function while improving patient quality of life.

ggtreeExtra: Compact Visualization of Richly Annotated Phylogenetic Data.

We present the ggtreeExtra package for visualizing heterogeneous data with a phylogenetic tree in a circular or rectangular layout (https://www.bioconductor.org/packages/ggtreeExtra). The package supports more data types and visualization methods than other tools. It supports using the grammar of graphics syntax to present data on a tree with richly annotated layers and allows evolutionary statistics inferred by commonly used software to be integrated and visualized with external data. GgtreeExtra is a universal tool for tree data visualization. It extends the applications of the phylogenetic tree in different disciplines by making more domain-specific data to be available to visualize and interpret in the evolutionary context.

Using clusterProfiler to characterize multiomics data.

With the advent of multiomics, software capable of multidimensional enrichment analysis has become increasingly crucial for uncovering gene set variations in biological processes and disease pathways. This is essential for elucidating disease mechanisms and identifying potential therapeutic targets. clusterProfiler stands out for its comprehensive utilization of databases and advanced visualization features. Importantly, clusterProfiler supports various biological knowledge, including Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, through performing over-representation and gene set enrichment analyses. A key feature is that clusterProfiler allows users to choose from various graphical outputs to visualize results, enhancing interpretability. This protocol describes innovative ways in which clusterProfiler has been used for integrating metabolomics and metagenomics analyses, identifying and characterizing transcription factors under stress conditions, and annotating cells in single-cell studies. In all cases, the computational steps can be completed within ~2 min. clusterProfiler is released through the Bioconductor project and can be accessed via https://bioconductor.org/packages/clusterProfiler/ .

MicrobiotaProcess: A comprehensive R package for deep mining microbiome.

The data output from microbiome research is growing at an accelerating rate, yet mining the data quickly and efficiently remains difficult. There is still a lack of an effective data structure to represent and manage data, as well as flexible and composable analysis methods. In response to these two issues, we designed and developed the MicrobiotaProcess package. It provides a comprehensive data structure, MPSE, to better integrate the primary and intermediate data, which improves the integration and exploration of the downstream data. Around this data structure, the downstream analysis tasks are decomposed and a set of functions are designed under a tidy framework. These functions independently perform simple tasks and can be combined to perform complex tasks. This gives users the ability to explore data, conduct personalized analyses, and develop analysis workflows. Moreover, MicrobiotaProcess can interoperate with other packages in the R community, which further expands its analytical capabilities. This article demonstrates the MicrobiotaProcess for analyzing microbiome data as well as other ecological data through several examples. It connects upstream data, provides flexible downstream analysis components, and provides visualization methods to assist in presenting and interpreting results.

Characterization of toxin-antitoxin systems from public sequencing data: A case study in Pseudomonas aeruginosa.

The toxin-antitoxin (TA) system is a widely distributed group of genetic modules that play important roles in the life of prokaryotes, with mobile genetic elements (MGEs) contributing to the dissemination of antibiotic resistance gene (ARG). The diversity and richness of TA systems in Pseudomonas aeruginosa, as one of the bacterial species with ARGs, have not yet been completely demonstrated. In this study, we explored the TA systems from the public genomic sequencing data and genome sequences. A small scale of genomic sequencing data in 281 isolates was selected from the NCBI SRA database, reassembling the genomes of these isolates led to the findings of abundant TA homologs. Furthermore, remapping these identified TA modules on 5,437 genome/draft genomes uncovers a great diversity of TA modules in P. aeruginosa. Moreover, manual inspection revealed several TA systems that were not yet reported in P. aeruginosa including the hok-sok, cptA-cptB, cbeA-cbtA, tomB-hha, and ryeA-sdsR. Additional annotation revealed that a large number of MGEs were closely distributed with TA. Also, 16% of ARGs are located relatively close to TA. Our work confirmed a wealth of TA genes in the unexplored P. aeruginosa pan-genomes, expanded the knowledge on P. aeruginosa, and provided methodological tips on large-scale data mining for future studies. The co-occurrence of MGE, ARG, and TA may indicate a potential interaction in their dissemination.

Exploring Epigenomic Datasets by ChIPseeker.

In many aspects of life, epigenetics, or the altering of phenotype without changes in sequences, play an essential role in biological function. A vast number of epigenomic datasets are emerging as a result of the advent of next-generation sequencing. Annotation, comparison, visualization, and interpretation of epigenomic datasets remain key aspects of computational biology. ChIPseeker is a Bioconductor package for performing these analyses among variable epigenomic datasets. The fundamental functions of ChIPseeker, including data preparation, annotation, comparison, and visualization, are explained in this article. ChIPseeker is a freely available open-source package that may be found at https://www.bioconductor.org/packages/ChIPseeker. © 2022 Wiley Periodicals LLC. Basic Protocol 1: ChIPseeker and epigenomic dataset preparation Basic Protocol 2: Annotation of epigenomic datasets Basic Protocol 3: Comparison of epigenomic datasets Basic Protocol 4: Visualization of annotated results Basic Protocol 5: Functional analysis of epigenomic datasets Basic Protocol 6: Genome-wide and locus-specific distribution of epigenomic datasets Basic Protocol 7: Heatmaps and metaplots of epigenomic datasets.

Stemness Analysis Uncovers That The Peroxisome Proliferator-Activated Receptor Signaling Pathway Can Mediate Fatty Acid Homeostasis In Sorafenib-Resistant Hepatocellular Carcinoma Cells.

Hepatocellular carcinoma (HCC) stem cells are regarded as an important part of individualized HCC treatment and sorafenib resistance. However, there is lacking systematic assessment of stem-like indices and associations with a response of sorafenib in HCC. Our study thus aimed to evaluate the status of tumor dedifferentiation for HCC and further identify the regulatory mechanisms under the condition of resistance to sorafenib. Datasets of HCC, including messenger RNAs (mRNAs) expression, somatic mutation, and clinical information were collected. The mRNA expression-based stemness index (mRNAsi), which can represent degrees of dedifferentiation of HCC samples, was calculated to predict drug response of sorafenib therapy and prognosis. Next, unsupervised cluster analysis was conducted to distinguish mRNAsi-based subgroups, and gene/geneset functional enrichment analysis was employed to identify key sorafenib resistance-related pathways. In addition, we analyzed and confirmed the regulation of key genes discovered in this study by combining other omics data. Finally, Luciferase reporter assays were performed to validate their regulation. Our study demonstrated that the stemness index obtained from transcriptomic is a promising biomarker to predict the response of sorafenib therapy and the prognosis in HCC. We revealed the peroxisome proliferator-activated receptor signaling pathway (the PPAR signaling pathway), related to fatty acid biosynthesis, that was a potential sorafenib resistance pathway that had not been reported before. By analyzing the core regulatory genes of the PPAR signaling pathway, we identified four candidate target genes, retinoid X receptor beta (RXRB), nuclear receptor subfamily 1 group H member 3 (NR1H3), cytochrome P450 family 8 subfamily B member 1 (CYP8B1) and stearoyl-CoA desaturase (SCD), as a signature to distinguish the response of sorafenib. We proposed and validated that the RXRB and NR1H3 could directly regulate NR1H3 and SCD, respectively. Our results suggest that the combined use of SCD inhibitors and sorafenib may be a promising therapeutic approach.

An enhancer variant at 16q22.1 predisposes to hepatocellular carcinoma via regulating PRMT7 expression.

Most cancer causal variants are found in gene regulatory elements, e.g., enhancers. However, enhancer variants predisposing to hepatocellular carcinoma (HCC) remain unreported. Here we conduct a genome-wide survey of HCC-susceptible enhancer variants through a three-stage association study in 11,958 individuals and identify rs73613962 (T > G) within the intronic region of PRMT7 at 16q22.1 as a susceptibility locus of HCC (OR = 1.41, P = 6.02 × 10-10). An enhancer dual-luciferase assay indicates that the rs73613962-harboring region has allele-specific enhancer activity. CRISPR-Cas9/dCas9 experiments further support the enhancer activity of this region to regulate PRMT7 expression. Mechanistically, transcription factor HNF4A binds to this enhancer region, with preference to the risk allele G, to promote PRMT7 expression. PRMT7 upregulation contributes to in vitro, in vivo, and clinical HCC-associated phenotypes, possibly by affecting the p53 signaling pathway. This concept of HCC pathogenesis may open a promising window for HCC prevention/treatment.

nCov2019: an R package for studying the COVID-19 coronavirus pandemic.

BACKGROUND: The global spreading of the COVID-19 coronavirus is still a serious public health challenge. Although there are a large number of public resources that provide statistics data, tools for retrospective historical data and convenient visualization are still valuable. To provide convenient access to data and visualization on the pandemic we developed an R package, nCov2019 (https://github.com/YuLab-SMU/nCov2019). METHODS: We collect stable and reliable data of COVID-19 cases from multiple authoritative and up-to-date sources, and aggregate the most recent and historical data for each country or even province. Medical progress information, including global vaccine development and therapeutics candidates, were also collected and can be directly accessed in our package. The nCov2019 package provides an R language interfaces and designed functions for data operation and presentation, a set of interfaces to fetch data subset intuitively, visualization methods, and a dashboard with no extra coding requirement for data exploration and interactive analysis. RESULTS: As of January 14, 2021, the global health crisis is still serious. The number of confirmed cases worldwide has reached 91,268,983. Following the USA, India has reached 10 million confirmed cases. Multiple peaks are observed in many countries. Under the efforts of researchers, 51 vaccines and 54 drugs are under development and 14 of these vaccines are already in the pre-clinical phase. DISCUSSION: The nCov2019 package provides detailed statistics data, visualization functions and the Shiny web application, which allows researchers to keep abreast of the latest epidemic spread overview.

clusterProfiler 4.0: A universal enrichment tool for interpreting omics data.

Functional enrichment analysis is pivotal for interpreting high-throughput omics data in life science. It is crucial for this type of tool to use the latest annotation databases for as many organisms as possible. To meet these requirements, we present here an updated version of our popular Bioconductor package, clusterProfiler 4.0. This package has been enhanced considerably compared with its original version published 9 years ago. The new version provides a universal interface for functional enrichment analysis in thousands of organisms based on internally supported ontologies and pathways as well as annotation data provided by users or derived from online databases. It also extends the dplyr and ggplot2 packages to offer tidy interfaces for data operation and visualization. Other new features include gene set enrichment analysis and comparison of enrichment results from multiple gene lists. We anticipate that clusterProfiler 4.0 will be applied to a wide range of scenarios across diverse organisms.

Ethylene trimerization with a half-sandwich titanium complex bearing a pendant thienyl group.

A half-sandwich titanium complex with a pendant thienyl group, activated by methylalumoxane (MAO), can trimerize ethylene to 1-hexene with considerable activity and high selectivity; the coodination of the sulfur atom to the titanium center is proposed to be responsible for the selective trimerization of ethylene.

A micro-nano porous oxide hybrid for efficient oxygen reduction in reduced-temperature solid oxide fuel cells.

Tremendous efforts to develop high-efficiency reduced-temperature (≤ 600°C) solid oxide fuel cells are motivated by their potentials for reduced materials cost, less engineering challenge, and better performance durability. A key obstacle to such fuel cells arises from sluggish oxygen reduction reaction kinetics on the cathodes. Here we reported that an oxide hybrid, featuring a nanoporous Sm(0.5)Sr(0.5)CoO(3-δ) (SSC) catalyst coating bonded onto the internal surface of a high-porosity La(0.9)Sr(0.1)Ga(0.8)Mg(0.2)O(3-δ) (LSGM) backbone, exhibited superior catalytic activity for oxygen reduction reactions and thereby yielded low interfacial resistances in air, e.g., 0.021 Ω cm(2) at 650°C and 0.043 Ω cm(2) at 600°C. We further demonstrated that such a micro-nano porous hybrid, adopted as the cathode in a thin LSGM electrolyte fuel cell, produced impressive power densities of 2.02 W cm(-2) at 650°C and 1.46 W cm(-2) at 600°C when operated on humidified hydrogen fuel and air oxidant.