The Bgee suite: integrated curated expression atlas and comparative transcriptomics in animals

Effect of deletion of the protein kinase PRKD1 on development of the mouse embryonic heart

Congenital heart disease (CHD) is the most common congenital anomaly, with an overall incidence of approximately 1% in the United Kingdom. Exome sequencing in large CHD cohorts has been performed to provide insights into the genetic aetiology of CHD. This includes a study of 1891 probands by our group in collaboration with others, which identified three novel genes-CDK13, PRKD1, and CHD4, in patients with syndromic CHD. PRKD1 encodes a serine/threonine protein kinase, which is important in a variety of fundamental cellular functions. Individuals with a heterozygous mutation in PRKD1 may have facial dysmorphism, ectodermal dysplasia and may have CHDs such as pulmonary stenosis, atrioventricular septal defects, coarctation of the aorta and bicuspid aortic valve. To obtain a greater appreciation for the role that this essential protein kinase plays in cardiogenesis and CHD, we have analysed a Prkd1 transgenic mouse model (Prkd1em1) carrying deletion of exon 2, causing loss of function. High-resolution episcopic microscopy affords detailed morphological 3D analysis of the developing heart and provides evidence for an essential role of Prkd1 in both normal cardiac development and CHD. We show that homozygous deletion of Prkd1 is associated with complex forms of CHD such as atrioventricular septal defects, and bicuspid aortic and pulmonary valves, and is lethal. Even in heterozygotes, cardiac differences occur. However, given that 97% of Prkd1 heterozygous mice display normal heart development, it is likely that one normal allele is sufficient, with the defects seen most likely to represent sporadic events. Moreover, mRNA and protein expression levels were investigated by RT-qPCR and western immunoblotting, respectively. A significant reduction in Prkd1 mRNA levels was seen in homozygotes, but not heterozygotes, compared to WT littermates. While a trend towards lower PRKD1 protein expression was seen in the heterozygotes, the difference was only significant in the homozygotes. There was no compensation by the related Prkd2 and Prkd3 at transcript level, as evidenced by RT-qPCR. Overall, we demonstrate a vital role of Prkd1 in heart development and the aetiology of CHD.

Procyanidin C1 inhibits tumor growth and metastasis in colon cancer via modulating miR-501-3p/HIGD1A axis

INTRODUCTION

Although colon (COAD) and rectal adenocarcinoma (READ) combined to refer to colorectal cancer (CRC), substantial clinical evidence urged that CRC should be treated as two different cancers due to compared with READ, COAD showed higher morbidity and worse 5-year survival.

OBJECTIVES

This study has tried to screen for the crucial gene that caused the worse prognosis and investigate its mechanism for mediating tumor growth and metastases in COAD. Meanwhile, the potential anti-COAD compound implicated in this mechanism was identified and testified from 1,855 food-borne chemical kits. This study aims to bring a new perspective to the development of new anti-COAD drugs and personalized medicine for patients with COAD.

METHODS AND RESULTS

The survival-related hub genes in COAD and READ were screened out from The Cancer Genome Atlas (TCGA) database and the results showed that HIGD1A, lower expressed in COAD than in READ, was associated with poor prognosis in COAD patients, but not in READ. Over-expressed HIGD1A suppressed CRC cell proliferation, invasion, and migration in vitro and in vivo. Meanwhile, the different expressed microRNA profiles between COAD and READ showed that miR-501-3p was highly expressed in COAD and inhibited HIGD1A expression by targeting 3'UTR of HIGD1A. MiR-501-3p mimics promoted cell proliferation and metastasis in CRC cells. In addition, Procyanidin C1 (PCC1), a kind of natural polyphenol has been verified as a potential miR-501-3p inhibitor. In vitro and in vivo, PCC1 promoted HIGD1A expression by suppressing miR-501-3p and resulted in inhibited tumor growth and metastasis.

CONCLUSION

The present study verified that miR-501-3p/HIGD1A axis mediated tumor growth and metastasis in COAD. PCC1, a flavonoid that riched in food exerts anti-COAD effects by inhibiting miR-501-3p and results in the latter losing the ability to suppress HIGD1A expression. Subsequently, unfettered HIGD1A inhibited tumor growth and metastasis in COAD.

Decoding bee cleptoparasitism through comparative transcriptomics of Coelioxoides waltheriae and its host Tetrapedia diversipes

Cleptoparasitism, also known as brood parasitism, is a widespread strategy among bee species in which the parasite lays eggs into the nests of the host species. Even though this behavior has significant ecological implications for the dynamics of several species, little is known about the molecular pathways associated with cleptoparasitism. To shed some light on this issue, we used gene expression data to perform a comparative analysis between two solitary neotropical bees: Coelioxoides waltheriae, an obligate parasite, and their specific host Tetrapedia diversipes. We found that ortholog genes involved in signal transduction, sensory perception, learning, and memory formation were differentially expressed between the cleptoparasite and the host. We hypothesize that these genes and their associated molecular pathways are engaged in cleptoparasitism-related processes and, hence, are appealing subjects for further investigation into functional and evolutionary aspects of cleptoparasitism in bees.

CXCR3 participates in asymmetric division of mouse oocytes by modulating actin dynamics

Extensive research has been conducted on the role of CXCR3 in immune responses and inflammation. However, the role of CXCR3 in the reproductive system, particularly in oocyte development, remains unknown. In this study, we present findings on the involvement of CXCR3 in the meiotic division process of mouse oocytes. We found CXCR3 was expressed consistently throughout the entire maturation process of mouse oocyte. Inhibition of CXCR3 impaired the asymmetric division of oocyte, while the injection of Cxcr3 mRNA was capable of restoring these defects. Further study showed that inhibition of CXCR3 perturbed spindle migration by affecting LIMK/cofilin pathway-mediated actin remodeling. Knockout of CXCR3 led to an upregulation of actin-binding protein and an increased ATP level in GV-stage oocytes, while maintaining normal actin dynamics during the process of meiosis. Additionally, we noticed the expression level of DYNLT1 is markedly elevated in CXCR3-null oocytes. DYNLT1 bound with the Arp2/3 complex, and knockdown of DYNLT1 in CXCR3-null oocytes impaired the organization of cytoplasmic actin, suggesting the regulatory role of DYNLT1 in actin organization, and the compensatory expression of DYNLT1 may contribute to maintain normal actin dynamics in CXCR3-knockout oocytes. In summary, our findings provide insights into the intricate network of actin dynamics associated with CXCR3 during oocyte meiosis.

Molecular mechanisms reconstruction from single-cell multi-omics data with HuMMuS

MOTIVATION

The molecular identity of a cell results from a complex interplay between heterogeneous molecular layers. Recent advances in single-cell sequencing technologies have opened the possibility to measure such molecular layers of regulation.

RESULTS

Here, we present HuMMuS, a new method for inferring regulatory mechanisms from single-cell multi-omics data. Differently from the state-of-the-art, HuMMuS captures cooperation between biological macromolecules and can easily include additional layers of molecular regulation. We benchmarked HuMMuS with respect to the state-of-the-art on both paired and unpaired multi-omics datasets. Our results proved the improvements provided by HuMMuS in terms of transcription factor (TF) targets, TF binding motifs and regulatory regions prediction. Finally, once applied to snmC-seq, scATAC-seq and scRNA-seq data from mouse brain cortex, HuMMuS enabled to accurately cluster scRNA profiles and to identify potential driver TFs.

AVAILABILITY AND IMPLEMENTATION

HuMMuS is available at https://github.com/cantinilab/HuMMuS.

An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Chapare virus infection and current perspectives on dentistry

OBJECTIVES

This narrative review addresses relevant points about Chapare virus (CHAV) entry in oral cells, CHAV transmission, and preventive strategies in dental clinical settings. It is critical in dentistry due to the frequent presence of gingival hemorrhage occurred in CHAV-infected patients.

MATERIALS AND METHODS

Studies related to CHAV were searched in MEDLINE/PubMed, Scopus, EMBASE, and Web-of-Science databases without language restriction or year of publication.

RESULTS

Recently, the PAHO/WHO and CDC indicate a presence of human-to-human transmission of CHAV associated with direct contact with saliva, blood, or urine, and also through droplets or aerosols created in healthcare procedures. CHAV was detected in human oropharyngeal saliva and gingival bleeding was confirmed in all cases of CHAV hemorrhagic fever, including evidence of nosocomial CHAV transmission in healthcare workers. We revisited the human transferrin receptor 1 (TfR1) expression in oral, nasal, and salivary glands tissues, as well as, we firstly identified the critical residues in the pre-glycoprotein (GP) complex of CHAV that interacts with human TfR1 using cutting-edge in silico bioinformatics platforms associated with molecular dynamic analysis.

CONCLUSIONS

In this multidisciplinary view, we also point out critical elements to provide perspectives on the preventive strategies for dentists and frontline healthcare workers against CHAV, and in the implementation of salivary diagnostic platforms for virus detection, which can be critical to an urgent plan to prevent human-to-human transmission based on current evidence.

CLINICAL RELEVANCE

The preventive strategies in dental clinical settings are pivotal due to the aerosol-generating procedures in dentistry with infected patients or suspected cases of CHAV infection.

Secreted novel AID/APOBEC-like deaminase 1 (SNAD1) - a new important player in fish immunology

The AID/APOBECs are a group of zinc-dependent cytidine deaminases that catalyse the deamination of bases in nucleic acids, resulting in a cytidine to uridine transition. Secreted novel AID/APOBEC-like deaminases (SNADs), characterized by the presence of a signal peptide are unique among all of intracellular classical AID/APOBECs, which are the central part of antibody diversity and antiviral defense. To date, there is no available knowledge on SNADs including protein characterization, biochemical characteristics and catalytic activity. We used various in silico approaches to define the phylogeny of SNADs, their common structural features, and their potential structural variations in fish species. Our analysis provides strong evidence of the universal presence of SNAD1 proteins/transcripts in fish, in which expression commences after hatching and is highest in anatomical organs linked to the immune system. Moreover, we searched published fish data and identified previously, "uncharacterized proteins" and transcripts as SNAD1 sequences. Our review into immunological research suggests SNAD1 role in immune response to infection or immunization, and interactions with the intestinal microbiota. We also noted SNAD1 association with temperature acclimation, environmental pollution and sex-based expression differences, with females showing higher level. To validate in silico predictions we performed expression studies of several SNAD1 gene variants in carp, which revealed distinct patterns of responses under different conditions. Dual sensitivity to environmental and pathogenic stress highlights its importance in the fish and potentially enhancing thermotolerance and immune defense. Revealing the biological roles of SNADs represents an exciting new area of research related to the role of DNA and/or RNA editing in fish biology.

In vivo CRISPR screening directly targeting testicular cells

CRISPR-Cas9 short guide RNA (sgRNA) library screening is a powerful approach to understand the molecular mechanisms of biological phenomena. However, its in vivo application is currently limited. Here, we developed our previously established in vitro revival screening method into an in vivo one to identify factors involved in spermatogenesis integrity by utilizing sperm capacitation as an indicator. By introducing an sgRNA library into testicular cells, we successfully pinpointed the retinal degeneration 3 (Rd3) gene as a significant factor in spermatogenesis. Single-cell RNA sequencing (scRNA-seq) analysis highlighted the high expression of Rd3 in round spermatids, and proteomics analysis indicated that Rd3 interacts with mitochondria. To search for cell-type-specific signaling pathways based on scRNA-seq and proteomics analyses, we developed a computational tool, Hub-Explorer. Through this, we discovered that Rd3 modulates oxidative stress by regulating mitochondrial distribution upon ciliogenesis induction. Collectively, our screening system provides a valuable in vivo approach to decipher molecular mechanisms in biological processes.

Where Are the Formerly Y-linked Genes in the Ryukyu Spiny Rat that has Lost its Y Chromosome?

It has been predicted that the highly degenerate mammalian Y chromosome will be lost eventually. Indeed, Y was lost in the Ryukyu spiny rat Tokudaia osimensis, but the fate of the formerly Y-linked genes is not completely known. We looked for all 12 ancestrally Y-linked genes in a draft T. osimensis genome sequence. Zfy1, Zfy2, Kdm5d, Eif2s3y, Usp9y, Uty, and Ddx3y are putatively functional and are now located on the X chromosome, whereas Rbmy, Uba1y, Ssty1, Ssty2, and Sry are missing or pseudogenized. Tissue expressions of the mouse orthologs of the retained genes are significantly broader/higher than those of the lost genes, suggesting that the destinies of the formerly Y-linked genes are related to their original expressions. Interestingly, patterns of gene retention/loss are significantly more similar than by chance across four rodent lineages where Y has been independently lost, indicating a level of certainty in the fate of Y-linked genes even when the chromosome is gone.

IGF2BPs as novel m6A readers: Diverse roles in regulating cancer cell biological functions, hypoxia adaptation, metabolism, and immunosuppressive tumor microenvironment

m6A methylation is the most frequent modification of mRNA in eukaryotes and plays a crucial role in cancer progression by regulating biological functions. Insulin-like growth factor 2 mRNA-binding proteins (IGF2BP) are newly identified m6A 'readers'. They belong to a family of RNA-binding proteins, which bind to the m6A sites on different RNA sequences and stabilize them to promote cancer progression. In this review, we summarize the mechanisms by which different upstream factors regulate IGF2BP in cancer. The current literature analyzed here reveals that the IGF2BP family proteins promote cancer cell proliferation, survival, and chemoresistance, inhibit apoptosis, and are also associated with cancer glycolysis, angiogenesis, and the immune response in the tumor microenvironment. Therefore, with the discovery of their role as 'readers' of m6A and the characteristic re-expression of IGF2BPs in cancers, it is important to elucidate their mechanism of action in the immunosuppressive tumor microenvironment. We also describe in detail the regulatory and interaction network of the IGF2BP family in downstream target RNAs and discuss their potential clinical applications as diagnostic and prognostic markers, as well as recent advances in IGF2BP biology and associated therapeutic value.

Parathyroid hormone-PTH1R signaling in cardiovascular disease and homeostasis

Primary hyperparathyroidism (pHPT) afflicts our aging population with an incidence approaching 50 per 100 000 patient-years at a female:male ratio of ~3:1. Decisions surrounding surgical management are currently driven by age, hypercalcemia severity, presence of osteoporosis, renal insufficiency, or hypercalciuria with or without nephrolithiasis. Cardiovascular (CV) disease (CVD) is not systematically considered. This is notable since the parathyroid hormone (PTH) 1 receptor (PTH1R) is biologically active in the vasculature, and adjusted CV mortality risk is increased almost threefold in individuals with pHPT who do not meet contemporary recommendations for surgical cure. We provide an overview of epidemiology, pharmacology, and physiology that highlights the need to: (i) identify biomarkers that establish a healthy 'set point' for CV PTH1R signaling tone; (ii) better understand the pharmacokinetic-pharmacodynamic (PK-PD) relationships of PTH1R ligands in CV homeostasis; and (iii) incorporate CVD risk assessment into the management of hyperparathyroidism.

DrosOMA: the Drosophila Orthologous Matrix browser

Background

Comparative genomic analyses to delineate gene evolutionary histories inform the understanding of organismal biology by characterising gene and gene family origins, trajectories, and dynamics, as well as enabling the tracing of speciation, duplication, and loss events, and facilitating the transfer of gene functional information across species. Genomic data are available for an increasing number of species from the genus Drosophila, however, a dedicated resource exploiting these data to provide the research community with browsable results from genus-wide orthology delineation has been lacking.

Methods

Using the OMA Orthologous Matrix orthology inference approach and browser deployment framework, we catalogued orthologues across a selected set of Drosophila species with high-quality annotated genomes. We developed and deployed a dedicated instance of the OMA browser to facilitate intuitive exploration, visualisation, and downloading of the genus-wide orthology delineation results.

Results

DrosOMA - the Drosophila Orthologous Matrix browser, accessible from https://drosoma.dcsr.unil.ch/ - presents the results of orthology delineation for 36 drosophilids from across the genus and four outgroup dipterans. It enables querying and browsing of the orthology data through a feature-rich web interface, with gene-view, orthologous group-view, and genome-view pages, including comprehensive gene name and identifier cross-references together with available functional annotations and protein domain architectures, as well as tools to visualise local and global synteny conservation.

Conclusions

The DrosOMA browser demonstrates the deployability of the OMA browser framework for building user-friendly orthology databases with dense sampling of a selected taxonomic group. It provides the Drosophila research community with a tailored resource of browsable results from genus-wide orthology delineation.

The Monarch Initiative in 2024: an analytic platform integrating phenotypes, genes and diseases across species

Bridging the gap between genetic variations, environmental determinants, and phenotypic outcomes is critical for supporting clinical diagnosis and understanding mechanisms of diseases. It requires integrating open data at a global scale. The Monarch Initiative advances these goals by developing open ontologies, semantic data models, and knowledge graphs for translational research. The Monarch App is an integrated platform combining data about genes, phenotypes, and diseases across species. Monarch's APIs enable access to carefully curated datasets and advanced analysis tools that support the understanding and diagnosis of disease for diverse applications such as variant prioritization, deep phenotyping, and patient profile-matching. We have migrated our system into a scalable, cloud-based infrastructure; simplified Monarch's data ingestion and knowledge graph integration systems; enhanced data mapping and integration standards; and developed a new user interface with novel search and graph navigation features. Furthermore, we advanced Monarch's analytic tools by developing a customized plugin for OpenAI's ChatGPT to increase the reliability of its responses about phenotypic data, allowing us to interrogate the knowledge in the Monarch graph using state-of-the-art Large Language Models. The resources of the Monarch Initiative can be found at monarchinitiative.org and its corresponding code repository at github.com/monarch-initiative/monarch-app.

The SIB Swiss Institute of Bioinformatics Semantic Web of data

The SIB Swiss Institute of Bioinformatics (https://www.sib.swiss/) is a federation of bioinformatics research and service groups. The international life science community in academia and industry has been accessing the freely available databases provided by SIB since its inception in 1998. In this paper we present the 11 databases which currently offer semantically enriched data in accordance with the FAIR principles (Findable, Accessible, Interoperable, Reusable), as well as the Swiss Personalized Health Network initiative (SPHN) which also employs this enrichment. The semantic enrichment facilitates the manipulation of large data sets from public databases and private data sets. Examples are provided to illustrate that the data from the SIB databases can not only be queried using precise criteria individually, but also across multiple databases, including a variety of non-SIB databases. Data manipulation, be it exploration, extraction, annotation, combination, and publication, is possible using the SPARQL query language. Providing documentation, tutorials and sample queries makes it easier to navigate this web of semantic data. Through this paper, the reader will discover how the existing SIB knowledge graphs can be leveraged to tackle the complex biological or clinical questions that are being addressed today.

Enhancing immune responses of ESC-based TAA cancer vaccines with a novel OMV delivery system

Embryonic stem cell (ESC)-derived epitopes can act as therapeutic tumor vaccines against different types of tumors Jin (Adv Healthc Mater 2023). However, these epitopes have poor immunogenicity and stimulate insufficient CD8+ T cell responses, which motivated us to develop a new method to deliver and enhance their effectiveness. Bacterial outer membrane vesicles (OMVs) can serve as immunoadjuvants and act as a delivery vector for tumor antigens. In the current study, we engineered a new OMV platform for the co-delivery of ESC-derived tumor antigens and immune checkpoint inhibitors (PD-L1 antibody). An engineered Staphylococcal Protein A (SpA) was created to non-specifically bind to anti-PD-L1 antibody. SpyCatcher (SpC) and SpA were fused into the cell outer membrane protein OmpA to capture SpyTag-attached peptides and PD-L1 antibody, respectively. The modified OMV was able to efficiently conjugate with ESC-derived TAAs and PD-L1 antibody (SpC-OMVs + SpT-peptides + anti-PD-L1), increasing the residence time of TAAs in the body. The results showed that the combination therapy of ESC-based TAAs and PD-L1 antibody delivered by OMV had significant inhibitory effects in mouse tumor model. Specifically, it was effective in reducing tumor growth by enhancing IFN-γ-CD8+ T cell responses and increasing the number of CD8+ memory cells and antigen-specific T cells. Overall, the new OMV delivery system is a versatile platform that can enhance the immune responses of ESC-based TAA cancer vaccines.

Collision-Induced Unfolding, Tandem MS, Bottom-up Proteomics, and Interactomics for Identification of Protein Complexes in Native Surface Mass Spectrometry

Endogenously occurring salts and nonvolatile matrix components in untreated biological surfaces can suppress protein ionization and promote adduct formation, challenging protein identification. Characterization of labile proteins within biological specimens is particularly demanding because additional purification or sample treatment steps can be time-intensive and can disrupt noncovalent interactions. It is demonstrated that the combined use of collision-induced unfolding, tandem mass spectrometry, and bottom-up proteomics improves protein characterization in native surface mass spectrometry (NSMS). This multiprong analysis is achieved by acquiring NSMS, MS/MS, ion mobility (IM), and bottom-up proteomics data from a single surface extracted sample. The validity of this multiprong approach was confirmed by the successful characterization of nine surface-deposited proteins, with molecular weights ranging from 8 to 147 kDa, in two separate mixtures. Bottom-up proteomics provided a list of proteins to match against observed proteins in NSMS and their detected subunits in tandem MS. The method was applied to characterize endogenous proteins from untreated chicken liver samples. The subcapsular liver sampling for NSMS analysis allowed for the detection of endogenous proteins with molecular weights of up to ∼220 kDa. Moreover, using IM-MS, collision cross sections and collision-induced unfolding pathways of enzymatic proteins and protein complexes of up to 145 kDa were obtained.

AutophagyNet: high-resolution data source for the analysis of autophagy and its regulation

Macroautophagy/autophagy is a highly-conserved catabolic procss eliminating dysfunctional cellular components and invading pathogens. Autophagy malfunction contributes to disorders such as cancer, neurodegenerative and inflammatory diseases. Understanding autophagy regulation in health and disease has been the focus of the last decades. We previously provided an integrated database for autophagy research, the Autophagy Regulatory Network (ARN). For the last eight years, this resource has been used by thousands of users. Here, we present a new and upgraded resource, AutophagyNet. It builds on the previous database but contains major improvements to address user feedback and novel needs due to the advancement in omics data availability. AutophagyNet contains updated interaction curation and integration of over 280,000 experimentally verified interactions between core autophagy proteins and their protein, transcriptional and post-transcriptional regulators as well as their potential upstream pathway connections. AutophagyNet provides annotations for each core protein about their role: 1) in different types of autophagy (mitophagy, xenophagy, etc.); 2) in distinct stages of autophagy (initiation, expansion, termination, etc.); 3) with subcellular and tissue-specific localization. These annotations can be used to filter the dataset, providing customizable download options tailored to the user's needs. The resource is available in various file formats (e.g. CSV, BioPAX and PSI-MI), and data can be analyzed and visualized directly in Cytoscape. The multi-layered regulation of autophagy can be analyzed by combining AutophagyNet with tissue- or cell type-specific (multi-)omics datasets (e.g. transcriptomic or proteomic data). The resource is publicly accessible at http://autophagynet.org.Abbreviations: ARN: Autophagy Regulatory Network; ATG: autophagy related; BCR: B cell receptor pathway; BECN1: beclin 1; GABARAP: GABA type A receptor-associated protein; IIP: innate immune pathway; LIR: LC3-interacting region; lncRNA: long non-coding RNA; MAP1LC3B: microtubule associated protein 1 light chain 3 beta; miRNA: microRNA; NHR: nuclear hormone receptor; PTM: post-translational modification; RTK: receptor tyrosine kinase; TCR: T cell receptor; TLR: toll like receptor.

Wnt16 Promotes Vascular Smooth Muscle Contractile Phenotype and Function via Taz (Wwtr1) Activation in Male LDLR-/- Mice

Wnt16 is expressed in bone and arteries, and maintains bone mass in mice and humans, but its role in cardiovascular physiology is unknown. We show that Wnt16 protein accumulates in murine and human vascular smooth muscle (VSM). WNT16 genotypes that convey risk for bone frailty also convey risk for cardiovascular events in the Dallas Heart Study. Murine Wnt16 deficiency, which causes postnatal bone loss, also reduced systolic blood pressure. Electron microscopy demonstrated abnormal VSM mitochondrial morphology in Wnt16-null mice, with reductions in mitochondrial respiration. Following angiotensin-II (AngII) infusion, thoracic ascending aorta (TAA) dilatation was greater in Wnt16-/- vs Wnt16+/+ mice (LDLR-/- background). Acta2 (vascular smooth muscle alpha actin) deficiency has been shown to impair contractile phenotype and worsen TAA aneurysm with concomitant reductions in blood pressure. Wnt16 deficiency reduced expression of Acta2, SM22 (transgelin), and other contractile genes, and reduced VSM contraction induced by TGFβ. Acta2 and SM22 proteins were reduced in Wnt16-/- VSM as was Ankrd1, a prototypic contractile target of Yap1 and Taz activation via TEA domain (TEAD)-directed transcription. Wnt16-/- VSM exhibited reduced nuclear Taz and Yap1 protein accumulation. SiRNA targeting Wnt16 or Taz, but not Yap1, phenocopied Wnt16 deficiency, and Taz siRNA inhibited contractile gene upregulation by Wnt16. Wnt16 incubation stimulated mitochondrial respiration and contraction (reversed by verteporfin, a Yap/Taz inhibitor). SiRNA targeting Taz inhibitors Ccm2 and Lats1/2 mimicked Wnt16 treatment. Wnt16 stimulated Taz binding to Acta2 chromatin and H3K4me3 methylation. TEAD cognates in the Acta2 promoter conveyed transcriptional responses to Wnt16 and Taz. Wnt16 regulates cardiovascular physiology and VSM contractile phenotype, mediated via Taz signaling.

Introgression of pigs in Taihu Lake region possibly contributed to the improvement of fertility in Danish Large White pigs

BACKGROUND

Eurasian pigs have undergone lineage admixture throughout history. It has been confirmed that the genes of indigenous pig breeds in China have been introduced into Western commercial pigs, providing genetic materials for breeding Western pigs. Pigs in Taihu Lake region (TL), such as the Meishan pig and Erhualian pig, serve as typical representatives of indigenous pig breeds in China due to their high reproductive performances. These pigs have also been imported into European countries in 1970 and 1980 s. They have played a positive role in improving the reproductive performances in European commercial pigs such as French Large White pigs (FLW). However, it is currently unclear if the lineage of TL pigs have been introgressed into the Danish Large White pigs (DLW), which are also known for their high reproductive performances in European pigs. To systematically identify genomic regions in which TL pigs have introgressed into DLW pigs and their physiological functions, we collected the re-sequencing data from 304 Eurasian pigs, to identify shared haplotypes between DLW and TL pigs.

RESULTS

The findings revealed the presence of introgressed genomic regions from TL pigs in the genome of DLW pigs indeed. The genes annotated within these regions were found to be mainly enriched in neurodevelopmental pathways. Furthermore, we found that the 115 kb region located in SSC16 exhibited highly shared haplotypes between TL and DLW pigs. The major haplotype of TL pigs in this region could significantly improve reproductive performances in various pig populations. Around this genomic region, NDUFS4 gene was highly expressed and showed differential expression in multiple reproductive tissues between extremely high and low farrowing Erhualian pigs. This suggested that NDUFS4 gene could be an important candidate causal gene responsible for affecting the reproductive performances of DLW pigs.

CONCLUSIONS

Our study has furthered our knowledge of the pattern of introgression from TL into DLW pigs and the potential effects on the fertility of DLW pigs.

CoSIA: an R Bioconductor package for CrOss Species Investigation and Analysis

SUMMARY

High-throughput sequencing technologies have enabled cross-species comparative transcriptomic studies; however, there are numerous challenges for these studies due to biological and technical factors. We developed CoSIA (Cross-Species Investigation and Analysis), a Bioconductor R package and Shiny app that provides an alternative framework for cross-species transcriptomic comparison of non-diseased wild-type RNA sequencing gene expression data from Bgee across tissues and species (human, mouse, rat, zebrafish, fly, and nematode) through visualization of variability, diversity, and specificity metrics.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lasseignelab/CoSIA.

Genetic variation in solute carrier family 5 member 2 mimicking sodium-glucose co-transporter 2-inhibition and risk of cardiovascular disease and all-cause mortality: reduced risk not explained by lower plasma glucose

AIMS

Treatment with sodium-glucose co-transporter 2 (SGLT2)-inhibitors reduces the risk of cardiovascular disease and mortality, but the mechanism is unclear. We hypothesized that a functional genetic variant in solute carrier family 5 member 2 (SLC5A2), known to be associated with familial renal glucosuria, would mimic pharmacological SGLT2-inhibition, and thus provide an opportunity to examine potential mediators of the effects on lower risk of cardiovascular disease and mortality.

METHODS AND RESULTS

We examined 112 712 individuals from the Copenhagen City Heart Study and Copenhagen General Population Study (CCHS + CGPS), 488 687 from the UK Biobank, and 342 499 from FinnGen, genotyped for SLC5A2 rs61742739, c.1961A > G; p.(Asn654Ser). The 2.0% heterozygotes and 0.01% homozygotes were pooled as carriers and compared with the 98% non-carriers. First, we examined the risk of cardiovascular disease and mortality; second, whether carrying the variant was associated with potential mediators of the effect; and third, whether identified potential mediators could explain the observed reduced risk of cardiovascular disease and mortality. In the CCHS + CGPS, carriers vs. non-carries had a 31% lower risk of heart failure, 21% lower risk of myocardial infarction, 16% lower risk of ischaemic heart disease, and 22% lower risk of all-cause mortality. Corresponding values in meta-analyses of the three studies combined were lower risk by 10%, 6%, 6%, and 10%, respectively. The SLC5A2 rs61742739 variant was not associated with a risk of ischaemic stroke or cardiovascular mortality. Of the lower risks observed in CCHS + CGPS, lower plasma glucose mediated 2.0%(P = 0.004) on heart failure, 3.1%(P = 0.09) on myocardial infarction, 4.1%(P = 0.02) on ischaemic heart disease, and 6.0%(P = 0.39) on all-cause mortality; corresponding values in the UK Biobank were 2.9%(P = 0.70), 1.5%(P = 0.77), 4.1%(P = 0.23), and 3.1%(P = 0.21), respectively.

CONCLUSION

A functional genetic variant in SLC5A2, mimicking SGLT2-inhibition, was associated with a lower risk of heart failure, myocardial infarction, ischaemic heart disease, and all-cause mortality. These effects were at most minimally mediated through lower plasma glucose.

Identification of Neurotransmission and Synaptic Biological Processes Disrupted in Autism Spectrum Disorder Using Interaction Networks and Community Detection Analysis

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder characterized by communication deficits and repetitive behavioral patterns. Hundreds of candidate genes have been implicated in ASD, including neurotransmission and synaptic (NS) genes; however, the genetic architecture of this disease is far from clear. In this study, we seek to clarify the biological processes affected by NS gene variants identified in individuals with ASD and the global networks that link those processes together. For a curated list of 1216 NS candidate genes, identified in multiple databases and the literature, we searched for ultra-rare (UR) loss-of-function (LoF) variants in the whole-exome sequencing dataset from the Autism Sequencing Consortium (N = 3938 cases). Filtering for population frequency was carried out using gnomAD (N = 60,146 controls). NS genes with UR LoF variants were used to construct a network of protein-protein interactions, and the network's biological communities were identified by applying the Leiden algorithm. We further explored the expression enrichment of network genes in specific brain regions. We identified 356 variants in 208 genes, with a preponderance of UR LoF variants in the PDE11A and SYTL3 genes. Expression enrichment analysis highlighted several subcortical structures, particularly the basal ganglia. The interaction network defined seven network communities, clustering synaptic and neurotransmitter pathways with several ubiquitous processes that occur in multiple organs and systems. This approach also uncovered biological pathways that are not usually associated with ASD, such as brain cytochromes P450 and brain mitochondrial metabolism. Overall, the community analysis suggests that ASD involves the disruption of synaptic and neurotransmitter pathways but also ubiquitous, but less frequently implicated, biological processes.

Characterization of an Osmr Conditional Knockout Mouse Model

Oncostatin M (OSM) is a member of the interleukin-6 (IL-6) family of cytokines and has been found to have distinct anti-inflammatory and pro-inflammatory properties in various cellular and disease contexts. OSM signals through two receptor complexes, one of which includes OSMRβ. To investigate OSM-OSMRβ signaling in adult hematopoiesis, we utilized the readily available conditional Osmrfl/fl mouse model B6;129-Osmrtm1.1Nat/J, which is poorly characterized in the literature. This model contains loxP sites flanking exon 2 of the Osmr gene. We crossed Osmrfl/fl mice to interferon-inducible Mx1-Cre, which is robustly induced in adult hematopoietic cells. We observed complete recombination of the Osmrfl allele and loss of exon 2 in hematopoietic (bone marrow) as well as non-hematopoietic (liver, lung, kidney) tissues. Using a TaqMan assay with probes downstream of exon 2, Osmr transcript was lower in the kidney but equivalent in bone marrow, lung, and liver from Osmrfl/fl Mx1-Cre versus Mx1-Cre control mice, suggesting that transcript is being produced despite loss of this exon. Western blots show that liver cells from Osmrfl/fl Mx1-Cre mice had complete loss of OSMR protein, while bone marrow, kidney, and lung cells had reduced OSMR protein at varying levels. RNA-seq analysis of a subpopulation of bone marrow cells (hematopoietic stem cells) finds that some OSM-stimulated genes, but not all, are suppressed in Osmrfl/fl Mx1-Cre cells. Together, our data suggest that the B6;129-Osmrtm1.1Nat/J model should be utilized with caution as loss of Osmr exon 2 has variable and tissue-dependent impact on mRNA and protein expression.

Genomic Variation, Population History, and Long-Term Genetic Adaptation to High Altitudes in Tibetan Partridge (Perdix hodgsoniae)

Species residing across elevational gradients display adaptations in response to environmental changes such as oxygen availability, ultraviolet radiation, and temperature. Here, we study genomic variation, gene expression, and long-term adaptation in Tibetan Partridge (Perdix hodgsoniae) populations residing across the elevational gradient of the Tibetan Plateau. We generated a high-quality draft genome and used it to carry out downstream population genomic and transcriptomic analysis. The P. hodgsoniae populations residing across various elevations were genetically distinct, and their phylogenetic clustering was consistent with their geographic distribution. We identified possible evidence of gene flow between populations residing in <3,000 and >4,200 m elevation that is consistent with known habitat expansion of high-altitude populations of P. hodgsoniae to a lower elevation. We identified a 60 kb haplotype encompassing the Estrogen Receptor 1 (ESR1) gene, showing strong genetic divergence between populations of P. hodgsoniae. We identified six single nucleotide polymorphisms within the ESR1 gene fixed for derived alleles in high-altitude populations that are strongly conserved across vertebrates. We also compared blood transcriptome profiles and identified differentially expressed genes (such as GAPDH, LDHA, and ALDOC) that correlated with differences in altitude among populations of P. hodgsoniae. These candidate genes from population genomics and transcriptomics analysis were enriched for neutrophil degranulation and glycolysis pathways, which are known to respond to hypoxia and hence may contribute to long-term adaptation to high altitudes in P. hodgsoniae. Our results highlight Tibetan Partridges as a useful model to study molecular mechanisms underlying long-term adaptation to high altitudes.

Evolution and divergence of teleost adrenergic receptors: why sometimes 'the drugs don't work' in fish

Adrenaline and noradrenaline, released as hormones and/or neurotransmitters, exert diverse physiological functions in vertebrates, and teleost fishes are widely used as model organisms to study adrenergic regulation; however, such investigations often rely on receptor subtype-specific pharmacological agents (agonists and antagonists; see Glossary) developed and validated in mammals. Meanwhile, evolutionary (phylogenetic and comparative genomic) studies have begun to unravel the diversification of adrenergic receptors (ARs) and reveal that whole-genome duplications and pseudogenization events in fishes results in notable distinctions from mammals in their genomic repertoire of ARs, while lineage-specific gene losses within teleosts have generated significant interspecific variability. In this Review, we visit the evolutionary history of ARs (including α1-, α2- and β-ARs) to highlight the prominent interspecific differences in teleosts, as well as between teleosts and other vertebrates. We also show that structural modelling of teleost ARs predicts differences in ligand binding affinity compared with mammalian orthologs. To emphasize the difficulty of studying the roles of different AR subtypes in fish, we collate examples from the literature of fish ARs behaving atypically compared with standard mammalian pharmacology. Thereafter, we focus on specific case studies of the liver, heart and red blood cells, where our understanding of AR expression has benefited from combining pharmacological approaches with molecular genetics. Finally, we briefly discuss the ongoing advances in 'omics' technologies that, alongside classical pharmacology, will provide abundant opportunities to further explore adrenergic signalling in teleosts.

Oral Infection, Oral Pathology and Salivary Diagnostics of Mpox Disease: Relevance in Dentistry and OMICs Perspectives

In this narrative review, we aim to point out the close relationship between mpox virus (MPXV) infection and the role of saliva as a diagnostic tool for mpox, considering the current molecular approach and in the perspective of OMICs application. The MPXV uses the host cell's rough endoplasmic reticulum, ribosomes, and cytoplasmic proteins to replicate its genome and synthesize virions for cellular exit. The presence of oral mucosa lesions associated with mpox infection is one of the first signs of infection; however, current diagnostic tools find it difficult to detect the virus before the rashes begin. MPXV transmission occurs through direct contact with an infected lesion and infected body fluids, including saliva, presenting a potential use of this fluid for diagnostic purposes. Currently available diagnostic tests for MPXV detection are performed either by real-time quantitative PCR (RT-qPCR) or ELISA, which presents several limitations since they are invasive tests. Despite current clinical trials with restricted sample size, MPXV DNA was detected in saliva with a sensitivity of 85%-100%. In this context, the application of transcriptomics, metabolomics, lipidomics, or proteomics analyses coupled with saliva can identify novel disease biomarkers. Thus, it is important to note that the identification and quantification of salivary DNA, RNA, lipid, protein, and metabolite can provide novel non-invasive biomarkers through the use of OMICs platforms aiding in the early detection and diagnosis of MPXV infection. Untargeted mass spectrometry (MS)-based proteomics reveals that some proteins also expressed in saliva were detected with greater expression differences in blood plasma when comparing mpox patients and healthy subjects, suggesting a promising alternative to be applied in screening or diagnostic platforms for mpox salivary diagnostics coupled to OMICs.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Molecular evolution of male reproduction across species with highly divergent sperm morphology in diverse murine rodents

Sperm competition can drive rapid evolution of male reproductive traits, but it remains unclear how variation in sperm competition intensity shapes phenotypic and molecular diversity across clades. Old World mice and rats (subfamily Murinae) comprise a rapid radiation and exhibit incredible diversity in sperm morphology and production. We combined phenotype and sequence data to model the evolution of reproductive traits and genes across 78 murine species. We identified several shifts towards smaller relative testes mass, a trait reflective of reduced sperm competition. Several sperm traits were associated with relative testes mass, suggesting that mating system evolution likely selects for convergent traits related to sperm competitive ability. Molecular evolutionary rates of spermatogenesis proteins also correlated with relative testes mass, but in an unexpected direction. We predicted that sperm competition would result in rapid divergence among species with large relative testes mass, but instead found that many spermatogenesis genes evolve more rapidly in species with smaller relative testes mass due to relaxed purifying selection. While some reproductive genes evolved under positive selection, relaxed selection played a greater role underlying rapid evolution in small testes species. Our work demonstrates that sexual selection can impose strong purifying selection shaping the evolution of male reproduction.

Analysis and Prediction of Pathogen Nucleic Acid Specificity for Toll-like Receptors in Vertebrates

Identification of key sequence, expression and function related features of nucleic acid-sensing host proteins is of fundamental importance to understand the dynamics of pathogen-specific host responses. To meet this objective, we considered toll-like receptors (TLRs), a representative class of membrane-bound sensor proteins, from 17 vertebrate species covering mammals, birds, reptiles, amphibians, and fishes in this comparative study. We identified the molecular signatures of host TLRs that are responsible for sensing pathogen nucleic acids or other pathogen-associated molecular patterns (PAMPs), and potentially play important roles in host defence mechanism. Interestingly, our findings reveal that such host-specific features are directly related to the strand (single or double) specificity of nucleic acid from pathogens. However, during host-pathogen interactions, such features were unable to explain the pathogenic PAMP (i.e., DNA, RNA or other) selectivity, suggesting a more complex mechanism. Using these features, we developed a number of machine learning models, of which Random Forest achieved a high performance (94.57% accuracy) to predict strand specificity of TLRs from protein-derived features. We applied the trained model to propose strand specificity of some previously uncharacterized distinct fish-specific novel TLRs (TLR18, TLR23, TLR24, TLR25, TLR27).

CoSIA: an R Bioconductor package for CrOss Species Investigation and Analysis

High throughput sequencing technologies have enabled cross-species comparative transcriptomic studies; however, there are numerous challenges for these studies due to biological and technical factors. We developed CoSIA (Cross-Species Investigation and Analysis), an Bioconductor R package and Shiny app that provides an alternative framework for cross-species transcriptomic comparison of non-diseased wild-type RNA sequencing gene expression data from Bgee across tissues and species (human, mouse, rat, zebrafish, fly, and nematode) through visualization of variability, diversity, and specificity metrics.

One-carbon metabolism is required for epigenetic stability in the mouse placenta

One-carbon metabolism, including the folate cycle, has a crucial role in fetal development though its molecular function is complex and unclear. The hypomorphic Mtrr ^gt allele is known to disrupt one-carbon metabolism, and thus methyl group availability, leading to several developmental phenotypes (e.g., neural tube closure defects, fetal growth anomalies). Remarkably, previous studies showed that some of the phenotypes were transgenerationally inherited. Here, we explored the genome-wide epigenetic impact of one-carbon metabolism in placentas associated with fetal growth phenotypes and determined whether specific DNA methylation changes were inherited. Firstly, methylome analysis of Mtrr ^gt/gt homozygous placentas revealed genome-wide epigenetic instability. Several differentially methylated regions (DMRs) were identified including at the Cxcl1 gene promoter and at the En2 gene locus, which may have phenotypic implications. Importantly, we discovered hypomethylation and ectopic expression of a subset of ERV elements throughout the genome of Mtrr ^gt/gt placentas with broad implications for genomic stability. Next, we determined that known spermatozoan DMRs in Mtrr ^gt/gt males were reprogrammed in the placenta with little evidence of direct or transgenerational germline DMR inheritance. However, some spermatozoan DMRs were associated with placental gene misexpression despite normalisation of DNA methylation, suggesting the inheritance of an alternative epigenetic mechanism. Integration of published wildtype histone ChIP-seq datasets with Mtrr ^gt/gt spermatozoan methylome and placental transcriptome datasets point towards H3K4me3 deposition at key loci. These data suggest that histone modifications might play a role in epigenetic inheritance in this context. Overall, this study sheds light on the mechanistic complexities of one-carbon metabolism in development and epigenetic inheritance.

Integrating Multi-Omics Analysis for Enhanced Diagnosis and Treatment of Glioblastoma: A Comprehensive Data-Driven Approach

The most aggressive primary malignant brain tumor in adults is glioblastoma (GBM), which has poor overall survival (OS). There is a high relapse rate among patients with GBM despite maximally safe surgery, radiation therapy, temozolomide (TMZ), and aggressive treatment. Hence, there is an urgent and unmet clinical need for new approaches to managing GBM. The current study identified modules (MYC, EGFR, PIK3CA, SUZ12, and SPRK2) involved in GBM disease through the NeDRex plugin. Furthermore, hub genes were identified in a comprehensive interaction network containing 7560 proteins related to GBM disease and 3860 proteins associated with signaling pathways involved in GBM. By integrating the results of the analyses mentioned above and again performing centrality analysis, eleven key genes involved in GBM disease were identified. ProteomicsDB and Gliovis databases were used for determining the gene expression in normal and tumor brain tissue. The NetworkAnalyst and the mGWAS-Explorer tools identified miRNAs, SNPs, and metabolites associated with these 11 genes. Moreover, a literature review of recent studies revealed other lists of metabolites related to GBM disease. The enrichment analysis of identified genes, miRNAs, and metabolites associated with GBM disease was performed using ExpressAnalyst, miEAA, and MetaboAnalyst tools. Further investigation of metabolite roles in GBM was performed using pathway, joint pathway, and network analyses. The results of this study allowed us to identify 11 genes (UBC, HDAC1, CTNNB1, TRIM28, CSNK2A1, RBBP4, TP53, APP, DAB1, PINK1, and RELN), five miRNAs (hsa-mir-221-3p, hsa-mir-30a-5p, hsa-mir-15a-5p, hsa-mir-130a-3p, and hsa-let-7b-5p), six metabolites (HDL, N6-acetyl-L-lysine, cholesterol, formate, N, N-dimethylglycine/xylose, and X2. piperidinone) and 15 distinct signaling pathways that play an indispensable role in GBM disease development. The identified top genes, miRNAs, and metabolite signatures can be targeted to establish early diagnostic methods and plan personalized GBM treatment strategies.

Dicodon-based measures for modeling gene expression

MOTIVATION

Codon usage preference patterns have been associated with modulation of translation efficiency, protein folding, and mRNA decay. However, new studies support that codon pair usage has also a remarkable effect at the gene expression level. Here, we expand the concept of CAI to answer if codon pair usage patterns can be understood in terms of codon usage bias, or if they offer new information regarding coding translation efficiency.

RESULTS

Through the implementation of a weighting strategy to consider the dicodon contributions, we observe that the dicodon-based measure has greater correlations with gene expression level than CAI. Interestingly, we have noted that dicodons associated with a low value of adaptiveness are related to dicodons which mediate strong translational inhibition in yeast. We have also noticed that some codon-pairs have a smaller dicodon contribution than estimated by the product of the respective codon contributions.

AVAILABILITY AND IMPLEMENTATION

Scripts, implemented in Python, are freely available for download at https://zenodo.org/record/7738276#.ZBIDBtLMIdU.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Biomedical discovery through the integrative biomedical knowledge hub (iBKH)

The abundance of biomedical knowledge gained from biological experiments and clinical practices is an invaluable resource for biomedicine. The emerging biomedical knowledge graphs (BKGs) provide an efficient and effective way to manage the abundant knowledge in biomedical and life science. In this study, we created a comprehensive BKG called the integrative Biomedical Knowledge Hub (iBKH) by harmonizing and integrating information from diverse biomedical resources. To make iBKH easily accessible for biomedical research, we developed a web-based, user-friendly graphical portal that allows fast and interactive knowledge retrieval. Additionally, we also implemented an efficient and scalable graph learning pipeline for discovering novel biomedical knowledge in iBKH. As a proof of concept, we performed our iBKH-based method for computational in-silico drug repurposing for Alzheimer's disease. The iBKH is publicly available.

Nuclear SUN1 stabilizes endothelial cell junctions via microtubules to regulate blood vessel formation

Endothelial cells line all blood vessels, where they coordinate blood vessel formation and the blood-tissue barrier via regulation of cell-cell junctions. The nucleus also regulates endothelial cell behaviors, but it is unclear how the nucleus contributes to endothelial cell activities at the cell periphery. Here, we show that the nuclear-localized linker of the nucleoskeleton and cytoskeleton (LINC) complex protein SUN1 regulates vascular sprouting and endothelial cell-cell junction morphology and function. Loss of murine endothelial Sun1 impaired blood vessel formation and destabilized junctions, angiogenic sprouts formed but retracted in SUN1-depleted sprouts, and zebrafish vessels lacking Sun1b had aberrant junctions and defective cell-cell connections. At the cellular level, SUN1 stabilized endothelial cell-cell junctions, promoted junction function, and regulated contractility. Mechanistically, SUN1 depletion altered cell behaviors via the cytoskeleton without changing transcriptional profiles. Reduced peripheral microtubule density, fewer junction contacts, and increased catastrophes accompanied SUN1 loss, and microtubule depolymerization phenocopied effects on junctions. Depletion of GEF-H1, a microtubule-regulated Rho activator, or the LINC complex protein nesprin-1 rescued defective junctions of SUN1-depleted endothelial cells. Thus, endothelial SUN1 regulates peripheral cell-cell junctions from the nucleus via LINC complex-based microtubule interactions that affect peripheral microtubule dynamics and Rho-regulated contractility, and this long-range regulation is important for proper blood vessel sprouting and junction integrity.

Natural selection and convergent evolution of the HOX gene family in Carnivora

HOX genes play a central role in the development and regulation of limb patterns. For mammals in the order Carnivora, limbs have evolved in different forms, and there are interesting cases of phenotypic convergence, such as the pseudothumb of the giant and red pandas, and the flippers or specialized limbs of the pinnipeds and sea otter. However, the molecular bases of limb development remain largely unclear. Here, we studied the molecular evolution of the HOX9 ~ 13 genes of 14 representative species in Carnivora and explored the molecular evolution of other HOX genes. We found that only one limb development gene, HOXC10, underwent convergent evolution between giant and red pandas and was thus an important candidate gene related to the development of pseudothumbs. No signals of amino acid convergence and natural selection were found in HOX9 ~ 13 genes between pinnipeds and sea otter, but there was evidence of positive selection and rapid evolution in four pinniped species. Overall, few HOX genes evolve via natural selection or convergent evolution, and these could be important candidate genes for further functional validation. Our findings provide insights into potential molecular mechanisms of the development of specialized pseudothumbs and flippers (or specialized limbs).

Gene expression databases for physiologically based pharmacokinetic modeling of humans and animal species

In drug research, developing a sound understanding of the key mechanistic drivers of pharmacokinetics (PK) for new molecular entities is essential for human PK and dose predictions. Here, characterizing the absorption, distribution, metabolism, and excretion (ADME) processes is crucial for a mechanistic understanding of the drug-target and drug-body interactions. Sufficient knowledge on ADME processes enables reliable interspecies and human PK estimations beyond allometric scaling. The physiologically based PK (PBPK) modeling framework allows the explicit consideration of organ-specific ADME processes. The sum of all passive and active ADME processes results in the observed plasma PK. Gene expression information can be used as surrogate for protein abundance and activity within PBPK models. The absolute and relative expression of ADME genes can differ between species and strains. This is affecting both, the PK and pharmacodynamics and is therefore posing a challenge for the extrapolation from preclinical findings to humans. We developed an automated workflow that generates whole-body gene expression databases for humans and other species relevant in drug development, animal health, nutritional sciences, and toxicology. Solely, bulk RNA-seq data curated and provided by the Swiss Institute of Bioinformatics from healthy, normal, and untreated primary tissue samples were considered as an unbiased reference of normal gene expression. The databases are interoperable with the Open Systems Pharmacology Suite (PK-Sim and MoBi) and enable seamless access to a central source of curated cross-species gene expression data. This will increase data transparency, increase reliability and reproducibility of PBPK model simulations, and accelerate mechanistic PBPK model development in the future.

Dietary Sugar Shifts Mitochondrial Metabolism and Small RNA Biogenesis in Sperm

Aims: Increasing concentrations of dietary sugar results in a linear accumulation of triglycerides in male Drosophila, while inducing a U-shaped obesity response in their offspring. Here, using a combination of proteomics and small RNA (sRNA) sequencing, we aimed at understanding the molecular underpinning in sperm for such plasticity. Results: Proteomic analysis of seminal vesicles revealed that increasing concentrations of dietary sugar resulted in a bell-shaped induction of proteins involved in metabolic/redox regulation. Using stains and in vivo redox reporter flies, this pattern could be explained by changes in sperm production of reactive oxygen species (ROS), more exactly mitochondria-derived H₂O₂. By quenching ROS with the antioxidant N-acetyl cysteine and performing sRNA-seq on sperm, we found that sperm miRNA is increased in response to ROS. Moreover, we found sperm mitosRNA to be increased in high-sugar diet conditions (independent of ROS). Reanalyzing our previously published data revealed a similar global upregulation of human sperm mitosRNA in response to a high-sugar diet, suggesting evolutionary conserved mechanisms. Innovation: This work highlights a fast response to dietary sugar in mitochondria-produced H₂O₂ in Drosophila sperm and identifies redox-sensitive miRNA downstream of this event. Conclusions: Our data support a model where changes in the sperm mitochondria in response to dietary sugar are the primary event, and changes in redox homoeostasis are secondary to mitochondrial ROS production. These data provide multiple candidates for paternal intergenerational metabolic responses as well as potential biomarkers for human male fertility.

The scalable precision medicine open knowledge engine (SPOKE): a massive knowledge graph of biomedical information

MOTIVATION

Knowledge graphs (KGs) are being adopted in industry, commerce and academia. Biomedical KG presents a challenge due to the complexity, size and heterogeneity of the underlying information.

RESULTS

In this work, we present the Scalable Precision Medicine Open Knowledge Engine (SPOKE), a biomedical KG connecting millions of concepts via semantically meaningful relationships. SPOKE contains 27 million nodes of 21 different types and 53 million edges of 55 types downloaded from 41 databases. The graph is built on the framework of 11 ontologies that maintain its structure, enable mappings and facilitate navigation. SPOKE is built weekly by python scripts which download each resource, check for integrity and completeness, and then create a 'parent table' of nodes and edges. Graph queries are translated by a REST API and users can submit searches directly via an API or a graphical user interface. Conclusions/Significance: SPOKE enables the integration of seemingly disparate information to support precision medicine efforts.

AVAILABILITY AND IMPLEMENTATION

The SPOKE neighborhood explorer is available at https://spoke.rbvi.ucsf.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Building a knowledge graph to enable precision medicine

Developing personalized diagnostic strategies and targeted treatments requires a deep understanding of disease biology and the ability to dissect the relationship between molecular and genetic factors and their phenotypic consequences. However, such knowledge is fragmented across publications, non-standardized repositories, and evolving ontologies describing various scales of biological organization between genotypes and clinical phenotypes. Here, we present PrimeKG, a multimodal knowledge graph for precision medicine analyses. PrimeKG integrates 20 high-quality resources to describe 17,080 diseases with 4,050,249 relationships representing ten major biological scales, including disease-associated protein perturbations, biological processes and pathways, anatomical and phenotypic scales, and the entire range of approved drugs with their therapeutic action, considerably expanding previous efforts in disease-rooted knowledge graphs. PrimeKG contains an abundance of 'indications', 'contradictions', and 'off-label use' drug-disease edges that lack in other knowledge graphs and can support AI analyses of how drugs affect disease-associated networks. We supplement PrimeKG's graph structure with language descriptions of clinical guidelines to enable multimodal analyses and provide instructions for continual updates of PrimeKG as new data become available.

Proteome-Wide Detection and Annotation of Receptor Tyrosine Kinases (RTKs): RTK-PRED and the TyReK Database

Receptor tyrosine kinases (RTKs) form a highly important group of protein receptors of the eukaryotic cell membrane. They control many vital cellular functions and are involved in the regulation of complex signaling networks. Mutations in RTKs have been associated with different types of cancers and other diseases. Although they are very important for proper cell function, they have been experimentally studied in a limited range of eukaryotic species. Currently, there is no available database for RTKs providing information about their function, expression, and interactions. Therefore, the identification of RTKs in multiple organisms, the documentation of their characteristics, and the collection of related information would be very useful. In this paper, we present a novel RTK detection pipeline (RTK-PRED) and the Receptor Tyrosine Kinases Database (TyReK-DB). RTK-PRED combines profile HMMs with transmembrane topology prediction to identify and classify potential RTKs. Proteins of all eukaryotic reference proteomes of the UniProt database were used as input in RTK-PRED leading to a filtered dataset of 20,478 RTKs. Based on the information collected for these RTKs from multiple databases, the relational TyReK database was created.

A Computational Inter-Species Study on Safrole Phase I Metabolism-Dependent Bioactivation: A Mechanistic Insight into the Study of Possible Differences among Species

Safrole, a 162.2 Da natural compound belonging to the alkenylbenzenes class, is classified as a possible carcinogen to humans by IARC (group IIB) and has proven to be genotoxic and carcinogenic to rodents. Despite its use as a food or feed additive, it is forbidden in many countries due to its documented toxicity; yet, it is still broadly present within food and feed and is particularly abundant in spices, herbs and essential oils. Specifically, safrole may exert its toxicity upon bioactivation to its proximate carcinogen 1'-hydroxy-safrole via specific members of the cytochrome P450 protein family with a certain inter/intra-species variability. To investigate this variability, an in-silico workflow based on molecular modelling, docking and molecular dynamics has been successfully applied. This work highlighted the mechanistic basis underpinning differences among humans, cats, chickens, goats, sheep, dogs, mice, pigs, rats and rabbits. The chosen metric to estimate the likeliness of formation of 1'-hydroxy-safrole by the species-specific cytochrome P450 under investigation allowed for the provision of a knowledge-based ground to rationally design and prioritise further experiments and deepen the current understanding of alkenylbenzenes bioactivation and CYPs mechanics. Both are crucial for a more informed framework of analysis for safrole toxicity.

High glucose causes developmental abnormalities in neuroepithelial cysts with actin and HK1 distribution changes

It has been reported that the offspring of diabetic pregnant women have an increased risk for neural tube defects. Previous studies in animal models suggested that high glucose induces cell apoptosis and epigenetic changes in the developing neural tube. However, effects on other cellular aspects such as the cell shape changes were not fully investigated. Actin dynamics plays essential roles in cell shape change. Disruption on actin dynamics is known to cause neural tube defects. In the present study, we used a 3D neuroepithelial cyst model and a rosette model, both cultured from human embryonic stem cells, to study the cellular effects caused by high glucose. By using these models, we observed couple of new changes besides increased apoptosis. First, we observed that high glucose disturbed the distribution of pH3 positive cells in the neuroepithelial cysts. Secondly, we found that high glucose exposure caused a relatively smaller actin inner boundary enclosed area, which was unlikely due to osmolarity changes. We further investigated key glucose metabolic enzymes in our models and the results showed that the distribution of hexokinase1 (HK1) was affected by high glucose. We observed that hexokinase1 has an apical-basal polarized distribution and is highest next to actin at the boundaries. hexokinase1 was more diffused and distributed less polarized under high glucose condition. Together, our observations broadened the cellular effects that may be caused by high glucose in the developing neural tube, especially in the secondary neurulation process.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Ten simple rules for using public biological data for your research

With an increasing amount of biological data available publicly, there is a need for a guide on how to successfully download and use this data. The 10 simple rules for using public biological data are: (1) use public data purposefully in your research; (2) evaluate data for your use case; (3) check data reuse requirements and embargoes; (4) be aware of ethics for data reuse; (5) plan for data storage and compute requirements; (6) know what you are downloading; (7) download programmatically and verify integrity; (8) properly cite data; (9) make reprocessed data and models Findable, Accessible, Interoperable, and Reusable (FAIR) and share; and (10) make pipelines and code FAIR and share. These rules are intended as a guide for researchers wanting to make use of available data and to increase data reuse and reproducibility.