The adeno-associated virus 2 genome and Rep 68/78 proteins interact with cellular sites of DNA damage.

Nuclear DNA viruses simultaneously access cellular factors that aid their life cycle while evading inhibitory factors by localizing to distinct nuclear sites. Adeno-associated viruses (AAVs), which are Dependoviruses in the family Parvovirinae, are non-enveloped icosahedral viruses, which have been developed as recombinant AAV vectors to express transgenes. AAV2 expression and replication occur in nuclear viral replication centers (VRCs), which relies on cellular replication machinery as well as coinfection by helper viruses such as adenoviruses or herpesviruses, or exogenous DNA damage to host cells. AAV2 infection induces a complex cellular DNA damage response (DDR), in response to either viral DNA or viral proteins expressed in the host nucleus during infection, where VRCs co-localized with DDR proteins. We have previously developed a modified iteration of a viral chromosome conformation capture (V3C-seq) assay to show that the autonomous parvovirus minute virus of mice localizes to cellular sites of DNA damage to establish and amplify its replication. Similar V3C-seq assays to map AAV2 show that the AAV2 genome co-localized with cellular sites of DNA damage under both non-replicating and replicating conditions. The AAV2 non-structural protein Rep 68/78, also localized to cellular DDR sites during both non-replicating and replicating infections, and also when ectopically expressed. Ectopically expressed Rep could be efficiently re-localized to DDR sites induced by micro-irradiation. Recombinant AAV2 gene therapy vector genomes derived from AAV2 localized to sites of cellular DNA damage to a lesser degree, suggesting that the inverted terminal repeat origins of replication were insufficient for targeting.

NRIMD, a Web Server for Analyzing Protein Allosteric Interactions Based on Molecular Dynamics Simulation.

Long-range allosteric communication between distant sites and active sites in proteins is central to biological regulation but still poorly characterized, limiting the development of protein engineering and drug design. Addressing this gap, NRIMD is an open-access web server for analyzing long-range interactions in proteins from molecular dynamics (MD) simulations, such as the effect of mutations at distal sites or allosteric ligand binding at allosteric sites on the active center. Based on our recent works on neural relational inference using graph neural networks, this cloud-based web server accepts MD simulation data on any length of residues in the alpha-carbon skeleton format from mainstream MD software. The input trajectory data are validated at the frontend deployed on the cloud and then processed on the backend deployed on a high-performance computer system with a collection of complementary tools. The web server provides a one-stop-shop MD analysis platform to predict long-range interactions and their paths between distant sites and active sites. It provides a user-friendly interface for detailed analysis and visualization. To the best of our knowledge, NRIMD is the first-of-its-kind online service to provide comprehensive long-range interaction analysis on MD simulations, which significantly lowers the barrier of predictions on protein long-range interactions using deep learning. The NRIMD web server is publicly available at https://nrimd.luddy.indianapolis.iu.edu/.

The Allele Catalog Tool: a web-based interactive tool for allele discovery and analysis.

BACKGROUND: The advancement of sequencing technologies today has made a plethora of whole-genome re-sequenced (WGRS) data publicly available. However, research utilizing the WGRS data without further configuration is nearly impossible. To solve this problem, our research group has developed an interactive Allele Catalog Tool to enable researchers to explore the coding region allelic variation present in over 1,000 re-sequenced accessions each for soybean, Arabidopsis, and maize. RESULTS: The Allele Catalog Tool was designed originally with soybean genomic data and resources. The Allele Catalog datasets were generated using our variant calling pipeline (SnakyVC) and the Allele Catalog pipeline (AlleleCatalog). The variant calling pipeline is developed to parallelly process raw sequencing reads to generate the Variant Call Format (VCF) files, and the Allele Catalog pipeline takes VCF files to perform imputations, functional effect predictions, and assemble alleles for each gene to generate curated Allele Catalog datasets. Both pipelines were utilized to generate the data panels (VCF files and Allele Catalog files) in which the accessions of the WGRS datasets were collected from various sources, currently representing over 1,000 diverse accessions for soybean, Arabidopsis, and maize individually. The main features of the Allele Catalog Tool include data query, visualization of results, categorical filtering, and download functions. Queries are performed from user input, and results are a tabular format of summary results by categorical description and genotype results of the alleles for each gene. The categorical information is specific to each species; additionally, available detailed meta-information is provided in modal popups. The genotypic information contains the variant positions, reference or alternate genotypes, the functional effect classes, and the amino-acid changes of each accession. Besides that, the results can also be downloaded for other research purposes. CONCLUSIONS: The Allele Catalog Tool is a web-based tool that currently supports three species: soybean, Arabidopsis, and maize. The Soybean Allele Catalog Tool is hosted on the SoyKB website ( https://soykb.org/SoybeanAlleleCatalogTool/ ), while the Allele Catalog Tool for Arabidopsis and maize is hosted on the KBCommons website ( https://kbcommons.org/system/tools/AlleleCatalogTool/Zmays and https://kbcommons.org/system/tools/AlleleCatalogTool/Athaliana ). Researchers can use this tool to connect variant alleles of genes with meta-information of species.

The bioinformatics toolbox for circRNA discovery and analysis.

Circular RNAs (circRNAs) are a unique class of RNA molecule identified more than 40 years ago which are produced by a covalent linkage via back-splicing of linear RNA. Recent advances in sequencing technologies and bioinformatics tools have led directly to an ever-expanding field of types and biological functions of circRNAs. In parallel with technological developments, practical applications of circRNAs have arisen including their utilization as biomarkers of human disease. Currently, circRNA-associated bioinformatics tools can support projects including circRNA annotation, circRNA identification and network analysis of competing endogenous RNA (ceRNA). In this review, we collected about 100 circRNA-associated bioinformatics tools and summarized their current attributes and capabilities. We also performed network analysis and text mining on circRNA tool publications in order to reveal trends in their ongoing development.

scGNN 2.0: a graph neural network tool for imputation and clustering of single-cell RNA-Seq data.

MOTIVATION: Gene expression imputation has been an essential step of the single-cell RNA-Seq data analysis workflow. Among several deep-learning methods, the debut of scGNN gained substantial recognition in 2021 for its superior performance and the ability to produce a cell-cell graph. However, the implementation of scGNN was relatively time-consuming and its performance could still be optimized. RESULTS: The implementation of scGNN 2.0 is significantly faster than scGNN thanks to a simplified close-loop architecture. For all eight datasets, cell clustering performance was increased by 85.02% on average in terms of adjusted rand index, and the imputation Median L1 Error was reduced by 67.94% on average. With the built-in visualizations, users can quickly assess the imputation and cell clustering results, compare against benchmarks and interpret the cell-cell interaction. The expanded input and output formats also pave the way for custom workflows that integrate scGNN 2.0 with other scRNA-Seq toolkits on both Python and R platforms. AVAILABILITY AND IMPLEMENTATION: scGNN 2.0 is implemented in Python (as of version 3.8) with the source code available at https://github.com/OSU-BMBL/scGNN2.0. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Coxiella burnetii Virulent Phase I and Avirulent Phase II Variants Differentially Manipulate Autophagy Pathway in Neutrophils.

Coxiella burnetii is an obligate intracellular Gram-negative bacterium that causes Q fever in humans. The virulent C. burnetii Nine Mile phase I (NMI) strain causes disease in animal models, while the avirulent NM phase II (NMII) strain does not. In this study, we found that NMI infection induces severe splenomegaly and bacterial burden in the spleen in BALB/c mice, while NMII infection does not. A significantly higher number of CD11b+ Ly6G+ neutrophils accumulated in the liver, lung, and spleen of NMI-infected mice than in NMII-infected mice. Thus, neutrophil accumulation correlates with NMI and NMII infection-induced inflammatory responses. In vitro studies also demonstrated that although NMII exhibited a higher infection rate than NMI in mouse bone marrow neutrophils (BMNs), NMI-infected BMNs survived longer than NMII-infected BMNs. These results suggest that the differential interactions of NMI and NMII with neutrophils may be related to their ability to cause disease in animals. To understand the molecular mechanism underlying the differential interactions of NMI and NMII with neutrophils, global transcriptomic gene expressions were compared between NMI- and NMII-infected BMNs by RNA sequencing (RNA-seq) analysis. Interestingly, several genes involved in autophagy-related pathways, particularly membrane trafficking and lipid metabolism, are upregulated in NMII-infected BMNs but downregulated in NMI-infected BMNs. Immunofluorescence and immunoblot analyses indicate that compared to NMI-infected BMNs, vacuoles in NMII-infected-BMNs exhibit increased autophagic flux along with phosphatidylserine translocation in the cell membrane. Similar to neutrophils, NMII activated LC3-mediated autophagy in human macrophages. These findings suggest that the differential manipulation of autophagy of NMI and NMII may relate to their pathogenesis.

Astrocytic GABA transporter 1 deficit in novel SLC6A1 variants mediated epilepsy: Connected from protein destabilization to seizures in mice and humans.

OBJECTIVE: Mutations in γ-aminobutyric acid (GABA) transporter 1 (GAT-1)-encoding SLC6A1 have been associated with myoclonic atonic epilepsy and other phenotypes. We determined the patho-mechanisms of the mutant GAT-1, in order to identify treatment targets. METHODS: We conducted whole-exome sequencing of patients with myoclonic atonic epilepsy (MAE) and characterized the seizure phenotypes and EEG patterns. We studied the protein stability and structural changes with homology modeling and machine learning tools. We characterized the function and trafficking of the mutant GAT-1 with 3H radioactive GABA uptake assay and confocal microscopy. We utilized different models including a knockin mouse and human astrocytes derived from induced pluripotent stem cells (iPSCs). We focused on astrocytes because of their direct impact of astrocytic GAT-1 in seizures. RESULTS: We identified four novel SLC6A1 variants associated with MAE and 2 to 4 Hz spike-wave discharges as a common EEG feature. Machine learning tools predicted that the variant proteins are destabilized. The variant protein had reduced expression and reduced GABA uptake due to endoplasmic reticular retention. The consistent observation was made in cortical and thalamic astrocytes from variant-knockin mice and human iPSC-derived astrocytes. The Slc6a+/A288V mouse, representative of MAE, had increased 5-7 Hz spike-wave discharges and absence seizures. INTERPRETATION: SLC6A1 variants in various locations of the protein peptides can cause MAE with similar seizure phenotypes and EEG features. Reduced GABA uptake is due to decreased functional GAT-1, which, in thalamic astrocytes, could result in increased extracellular GABA accumulation and enhanced tonic inhibition, leading to seizures and abnormal EEGs.

The NS1 protein of the parvovirus MVM Aids in the localization of the viral genome to cellular sites of DNA damage.

The autonomous parvovirus Minute Virus of Mice (MVM) localizes to cellular DNA damage sites to establish and sustain viral replication centers, which can be visualized by focal deposition of the essential MVM non-structural phosphoprotein NS1. How such foci are established remains unknown. Here, we show that NS1 localized to cellular sites of DNA damage independently of its ability to covalently bind the 5' end of the viral genome, or its consensus DNA binding sequence. Many of these sites were identical to those occupied by virus during infection. However, localization of the MVM genome to DNA damage sites occurred only when wild-type NS1, but not its DNA-binding mutant was expressed. Additionally, wild-type NS1, but not its DNA binding mutant, could localize a heterologous DNA molecule containing the NS1 binding sequence to DNA damage sites. These findings suggest that NS1 may function as a bridging molecule, helping the MVM genome localize to cellular DNA damage sites to facilitate ongoing virus replication.

Bioinformatics for plant and agricultural discoveries in the age of multiomics: A review and case study of maize nodal root growth under water deficit.

Advances in next-generation sequencing and other high-throughput technologies have facilitated multiomics research, such as genomics, epigenomics, transcriptomics, proteomics, metabolomics, and phenomics. The resultant emerging multiomics data have brought new challenges as well as opportunities, as seen in the plant and agriculture science domains. We reviewed several bioinformatic and computational methods, models, and platforms, and we have highlighted some of our in-house developed efforts aimed at multiomics data analysis, integration, and management issues faced by the research community. A case study using multiomics datasets generated from our studies of maize nodal root growth under water deficit stress demonstrates the power of these datasets and some other publicly available tools. This analysis also sheds light on the landscape of such applied bioinformatic tools currently available for plant and crop science studies and introduces emerging trends and how they may affect the future.

A rectified factor network based biclustering method for detecting cancer-related coding genes and miRNAs, and their interactions.

Detecting cancer-related genes and their interactions is a crucial task in cancer research. For this purpose, we proposed an efficient method, to detect coding genes, microRNAs (miRNAs), and their interactions related to a particular cancer or a cancer subtype using their expression data from the same set of samples. Firstly, biclusters specific to a particular type of cancer are detected based on rectified factor networks and ranked according to their associations with general cancers. Secondly, coding genes and miRNAs in each bicluster are prioritized by considering their differential expression and differential correlation values, protein-protein interaction data, and potential cancer markers. Finally, a rank fusion process is used to obtain the final comprehensive rank by combining multiple ranking results. We applied our proposed method on breast cancer datasets. Results show that our method outperforms other methods in detecting breast cancer-related coding genes and miRNAs. Furthermore, our method is very efficient in computing time, which can handle tens of thousands genes/miRNAs and hundreds of patients in hours on a desktop. This work may aid researchers in studying the genetic architecture of complex diseases, and improving the accuracy of diagnosis.

G2S: a web-service for annotating genomic variants on 3D protein structures.

Motivation: Accurately mapping and annotating genomic locations on 3D protein structures is a key step in structure-based analysis of genomic variants detected by recent large-scale sequencing efforts. There are several mapping resources currently available, but none of them provides a web API (Application Programming Interface) that supports programmatic access. Results: We present G2S, a real-time web API that provides automated mapping of genomic variants on 3D protein structures. G2S can align genomic locations of variants, protein locations, or protein sequences to protein structures and retrieve the mapped residues from structures. G2S API uses REST-inspired design and it can be used by various clients such as web browsers, command terminals, programming languages and other bioinformatics tools for bringing 3D structures into genomic variant analysis. Availability and implementation: The webserver and source codes are freely available at https://g2s.genomenexus.org. Contact: g2s@genomenexus.org. Supplementary information: Supplementary data are available at Bioinformatics online.

PGen: large-scale genomic variations analysis workflow and browser in SoyKB.

BACKGROUND: With the advances in next-generation sequencing (NGS) technology and significant reductions in sequencing costs, it is now possible to sequence large collections of germplasm in crops for detecting genome-scale genetic variations and to apply the knowledge towards improvements in traits. To efficiently facilitate large-scale NGS resequencing data analysis of genomic variations, we have developed "PGen", an integrated and optimized workflow using the Extreme Science and Engineering Discovery Environment (XSEDE) high-performance computing (HPC) virtual system, iPlant cloud data storage resources and Pegasus workflow management system (Pegasus-WMS). The workflow allows users to identify single nucleotide polymorphisms (SNPs) and insertion-deletions (indels), perform SNP annotations and conduct copy number variation analyses on multiple resequencing datasets in a user-friendly and seamless way. RESULTS: We have developed both a Linux version in GitHub ( https://github.com/pegasus-isi/PGen-GenomicVariations-Workflow ) and a web-based implementation of the PGen workflow integrated within the Soybean Knowledge Base (SoyKB), ( http://soykb.org/Pegasus/index.php ). Using PGen, we identified 10,218,140 single-nucleotide polymorphisms (SNPs) and 1,398,982 indels from analysis of 106 soybean lines sequenced at 15X coverage. 297,245 non-synonymous SNPs and 3330 copy number variation (CNV) regions were identified from this analysis. SNPs identified using PGen from additional soybean resequencing projects adding to 500+ soybean germplasm lines in total have been integrated. These SNPs are being utilized for trait improvement using genotype to phenotype prediction approaches developed in-house. In order to browse and access NGS data easily, we have also developed an NGS resequencing data browser ( http://soykb.org/NGS_Resequence/NGS_index.php ) within SoyKB to provide easy access to SNP and downstream analysis results for soybean researchers. CONCLUSION: PGen workflow has been optimized for the most efficient analysis of soybean data using thorough testing and validation. This research serves as an example of best practices for development of genomics data analysis workflows by integrating remote HPC resources and efficient data management with ease of use for biological users. PGen workflow can also be easily customized for analysis of data in other species.

Evaluation of genetic variation among Brazilian soybean cultivars through genome resequencing.

BACKGROUND: Soybean [Glycine max (L.) Merrill] is one of the most important legumes cultivated worldwide, and Brazil is one of the main producers of this crop. Since the sequencing of its reference genome, interest in structural and allelic variations of cultivated and wild soybean germplasm has grown. To investigate the genetics of the Brazilian soybean germplasm, we selected soybean cultivars based on the year of commercialization, geographical region and maturity group and resequenced their genomes. RESULTS: We resequenced the genomes of 28 Brazilian soybean cultivars with an average genome coverage of 14.8X. A total of 5,835,185 single nucleotide polymorphisms (SNPs) and 1,329,844 InDels were identified across the 20 soybean chromosomes, with 541,762 SNPs, 98,922 InDels and 1,093 CNVs that were exclusive to the 28 Brazilian cultivars. In addition, 668 allelic variations of 327 genes were shared among all of the Brazilian cultivars, including genes related to DNA-dependent transcription-elongation, photosynthesis, ATP synthesis-coupled electron transport, cellular respiration, and precursors of metabolite generation and energy. A very homogeneous structure was also observed for the Brazilian soybean germplasm, and we observed 41 regions putatively influenced by positive selection. Finally, we detected 3,880 regions with copy-number variations (CNVs) that could help to explain the divergence among the accessions evaluated. CONCLUSIONS: The large number of allelic and structural variations identified in this study can be used in marker-assisted selection programs to detect unique SNPs for cultivar fingerprinting. The results presented here suggest that despite the diversification of modern Brazilian cultivars, the soybean germplasm remains very narrow because of the large number of genome regions that exhibit low diversity. These results emphasize the need to introduce new alleles to increase the genetic diversity of the Brazilian germplasm.

Exploring Human Diseases and Biological Mechanisms by Protein Structure Prediction and Modeling.

Protein structure prediction and modeling provide a tool for understanding protein functions by computationally constructing protein structures from amino acid sequences and analyzing them. With help from protein prediction tools and web servers, users can obtain the three-dimensional protein structure models and gain knowledge of functions from the proteins. In this chapter, we will provide several examples of such studies. As an example, structure modeling methods were used to investigate the relation between mutation-caused misfolding of protein and human diseases including epilepsy and leukemia. Protein structure prediction and modeling were also applied in nucleotide-gated channels and their interaction interfaces to investigate their roles in brain and heart cells. In molecular mechanism studies of plants, rice salinity tolerance mechanism was studied via structure modeling on crucial proteins identified by systems biology analysis; trait-associated protein-protein interactions were modeled, which sheds some light on the roles of mutations in soybean oil/protein content. In the age of precision medicine, we believe protein structure prediction and modeling will play more and more important roles in investigating biomedical mechanism of diseases and drug design.

A Bayesian model for detection of high-order interactions among genetic variants in genome-wide association studies.

BACKGROUND: A central question for disease studies and crop improvements is how genetics variants drive phenotypes. Genome Wide Association Study (GWAS) provides a powerful tool for characterizing the genotype-phenotype relationships in complex traits and diseases. Epistasis (gene-gene interaction), including high-order interaction among more than two genes, often plays important roles in complex traits and diseases, but current GWAS analysis usually just focuses on additive effects of single nucleotide polymorphisms (SNPs). The lack of effective computational modelling of high-order functional interactions often leads to significant under-utilization of GWAS data. RESULTS: We have developed a novel Bayesian computational method with a Markov Chain Monte Carlo (MCMC) search, and implemented the method as a Bayesian High-order Interaction Toolkit (BHIT) for detecting epistatic interactions among SNPs. BHIT first builds a Bayesian model on both continuous data and discrete data, which is capable of detecting high-order interactions in SNPs related to case--control or quantitative phenotypes. We also developed a pipeline that enables users to apply BHIT on different species in different use cases. CONCLUSIONS: Using both simulation data and soybean nutritional seed composition studies on oil content and protein content, BHIT effectively detected some high-order interactions associated with phenotypes, and it outperformed a number of other available tools. BHIT is freely available for academic users at http://digbio.missouri.edu/BHIT/.

Hydrogen Regulates Ulcerative Colitis by Affecting the Intestinal Redox Environment.

The redox balance in the intestine plays an important role in maintaining intestinal homeostasis, and it is closely related to the intestinal mucosal barrier, intestinal inflammation, and the gut microbiota. Current research on the treatment of ulcerative colitis has focused on immune disorders, excessive inflammation, and oxidative stress. However, an imbalance in intestinal redox reaction plays a particularly critical role. Hydrogen is produced by some anaerobic bacteria via hydrogenases in the intestine. Increasing evidence suggests that hydrogen, as an inert gas, is crucial for immunity, inflammation, and oxidative stress and plays a protective role in ulcerative colitis. Hydrogen maintains the redox state balance in the intestine in ulcerative colitis and reduces damage to intestinal epithelial cells by exerting its selective antioxidant ability. Hydrogen also regulates the intestinal flora, reduces the harmful effects of bacteria on the intestinal epithelial barrier, promotes the restoration of normal anaerobic bacteria in the intestines, and ultimately improves the integrity of the intestinal epithelial barrier. The present review focuses on the therapeutic mechanisms of hydrogen-targeting ulcerative colitis.

scBSP: A fast and accurate tool for identifying spatially variable genes from spatial transcriptomic data.

Spatially resolved transcriptomics have enabled the inference of gene expression patterns within two and three-dimensional space, while introducing computational challenges due to growing spatial resolutions and sparse expressions. Here, we introduce scBSP, an open-source, versatile, and user-friendly package designed for identifying spatially variable genes in large-scale spatial transcriptomics. scBSP implements sparse matrix operation to significantly increase the computational efficiency in both computational time and memory usage, processing the high-definition spatial transcriptomics data for 19,950 genes on 181,367 spots within 10 seconds. Applied to diverse sequencing data and simulations, scBSP efficiently identifies spatially variable genes, demonstrating fast computational speed and consistency across various sequencing techniques and spatial resolutions for both two and three-dimensional data with up to millions of cells. On a sample with hundreds of thousands of sports, scBSP identifies SVGs accurately in seconds to on a typical desktop computer.

Altered mucosal bacteria and metabolomics in patients with Peutz-Jeghers syndrome.

BACKGROUND: Peutz-Jeghers syndrome (PJS) is a rare genetic disorder characterized by the development of pigmented spots, gastrointestinal polyps and increased susceptibility to cancers. Currently, most studies have investigated intestinal microbiota through fecal microbiota, and there are few reports about mucosa-associated microbiota. It remains valuable to search for the key intestinal microbiota or abnormal metabolic pathways linked to PJS. AIM: This study aimed to assess the structure and composition of mucosa-associated microbiota in patients with PJS and to explore the potential influence of intestinal microbiota disorders and metabolite changes on PJS. METHODS: The bacterial composition was analyzed in 13 PJS patients and 12 controls using 16S rRNA gene sequencing (Illumina MiSeq) for bacteria. Differential analyses of the intestinal microbiota were performed from the phylum to species level. Liquid chromatography-tandem mass spectrometry (LC‒MS) was used to detect the differentially abundant metabolites of PJS patients and controls to identify different metabolites and metabolic biomarkers of small intestinal mucosa samples. RESULTS: High-throughput sequencing confirmed the special characteristics and biodiversity of the mucosa microflora in patients with PJS. They had lower bacterial biodiversity than controls. The abundance of intestinal mucosal microflora was significantly lower than that of fecal microflora. In addition, lipid metabolism, amino acid metabolism, carbohydrate metabolism, nucleotide metabolism and other pathways were significantly different from those of controls, which were associated with the development of the enteric nervous system, intestinal inflammation and development of tumors. CONCLUSION: This is the first report on the mucosa-associated microbiota and metabolite profile of subjects with PJS, which may be meaningful to provide a structural basis for further research on intestinal microecology in PJS.

Prediction of protein stability changes upon single-point variant using 3D structure profile.

Identifying protein thermodynamic stability changes upon single-point variants is crucial for studying mutation-induced alterations in protein biophysics, genomic variants, and mutation-related diseases. In the last decade, various computational methods have been developed to predict the effects of single-point variants, but the prediction accuracy is still far from satisfactory for practical applications. Herein, we review approaches and tools for predicting stability changes upon the single-point variant. Most of these methods require tertiary protein structure as input to achieve reliable predictions. However, the availability of protein structures limits the immediate application of these tools. To improve the performance of a computational prediction from a protein sequence without experimental structural information, we introduce a new computational framework: MU3DSP. This method assesses the effects of single-point variants on protein thermodynamic stability based on point mutated protein 3D structure profile. Given a protein sequence with a single variant as input, MU3DSP integrates both sequence-level features and averaged features of 3D structures obtained from sequence alignment to PDB to assess the change of thermodynamic stability induced by the substitution. MU3DSP outperforms existing methods on various benchmarks, making it a reliable tool to assess both somatic and germline substitution variants and assist in protein design. MU3DSP is available as an open-source tool at https://github.com/hurraygong/MU3DSP.

Dimension-agnostic and granularity-based spatially variable gene identification.

Identifying spatially variable genes (SVGs) is critical in linking molecular cell functions with tissue phenotypes. Spatially resolved transcriptomics captures cellular-level gene expression with corresponding spatial coordinates in two or three dimensions and can be used to infer SVGs effectively. However, current computational methods may not achieve reliable results and often cannot handle three-dimensional spatial transcriptomic data. Here we introduce BSP (big-small patch), a spatial granularity-guided and non-parametric model to identify SVGs from two or three-dimensional spatial transcriptomics data in a fast and robust manner. This new method has been extensively tested in simulations, demonstrating superior accuracy, robustness, and high efficiency. BSP is further validated by substantiated biological discoveries in cancer, neural science, rheumatoid arthritis, and kidney studies with various types of spatial transcriptomics technologies.