PANTHER version 7: improved phylogenetic trees, orthologs and collaboration with the Gene Ontology Consortium

Research into how carvacrol and metformin affect several human proteins in a hyperglycemic condition: A comparative study in silico and in vitro

Carvacrol (CV) is an organic compound found in the essential oils of many aromatic herbs. It is nearly unfeasible to analyze all the current human proteins for a query ligand using in vitro and in vivo methods. This study aimed to clarify whether CV possesses an anti-diabetic feature via Docking-based inverse docking and molecular dynamic (MD) simulation and in vitro characterization against a set of novel human protein targets. Herein, the best poses of CV docking simulations according to binding energy ranged from -7.9 to -3.5 (kcal/mol). After pathway analysis of the protein list through GeneMANIA and WebGestalt, eight interacting proteins (DPP4, FBP1, GCK, HSD11β1, INSR, PYGL, PPARA, and PPARG) with CV were determined, and these proteins exhibited stable structures during the MD process with CV. In vitro application, statistically significant results were achieved only in combined doses with CV or metformin. Considering all these findings, PPARG and INSR, among these target proteins of CV, are FDA-approved targets for treating diabetes. Therefore, CV may be on its way to becoming a promising therapeutic compound for treating Diabetes Mellitus (DM). Our outcomes expose formerly unexplored potential target human proteins, whose association with diabetic disorders might guide new potential treatments for DM.

The role of interleukin-10 receptor alpha (IL10Rα) in Mycobacterium avium subsp. paratuberculosis infection of a mammary epithelial cell line

BACKGROUND

Johne's disease is a chronic wasting disease caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease is highly contagious and MAP infection in dairy cattle can eventually lead to death. With no available treatment for Johne's disease, genetic selection and improvements in management practices could help reduce its prevalence. In a previous study, the gene coding interleukin-10 receptor subunit alpha (IL10Rα) was associated with Johne's disease in dairy cattle. Our objective was to determine how IL10Rα affects the pathogenesis of MAP by examining the effect of a live MAP challenge on a mammary epithelial cell line (MAC-T) that had IL10Rα knocked out using CRISPR/cas9. The wild type and the IL10Rα knockout MAC-T cell lines were exposed to live MAP bacteria for 72 h. Thereafter, mRNA was extracted from infected and uninfected cells. Differentially expressed genes were compared between the wild type and the IL10Rα knockout cell lines. Gene ontology was performed based on the differentially expressed genes to determine which biological pathways were involved.

RESULTS

Immune system processes pathways were targeted to determine the effect of IL10Rα on the response to MAP infection. There was a difference in immune response between the wild type and IL10Rα knockout MAC-T cell lines, and less difference in immune response between infected and not infected IL10Rα knockout MAC-T cells, indicating IL10Rα plays an important role in the progression of MAP infection. Additionally, these comparisons allowed us to identify other genes involved in inflammation-mediated chemokine and cytokine signalling, interleukin signalling and toll-like receptor pathways.

CONCLUSIONS

Identifying differentially expressed genes in wild type and ILR10α knockout MAC-T cells infected with live MAP bacteria provided further evidence that IL10Rα contributes to mounting an immune response to MAP infection and allowed us to identify additional potential candidate genes involved in this process. We found there was a complex immune response during MAP infection that is controlled by many genes.

An inducible genetic tool to track and manipulate specific microglial states reveals their plasticity and roles in remyelination

Recent single-cell RNA sequencing studies have revealed distinct microglial states in development and disease. These include proliferative-region-associated microglia (PAMs) in developing white matter and disease-associated microglia (DAMs) prevalent in various neurodegenerative conditions. PAMs and DAMs share a similar core gene signature. However, the extent of the dynamism and plasticity of these microglial states, as well as their functional significance, remains elusive, partly due to the lack of specific tools. Here, we generated an inducible Cre driver line, Clec7a-CreERT2, that targets PAMs and DAMs in the brain parenchyma. Utilizing this tool, we profiled labeled cells during development and in several disease models, uncovering convergence and context-dependent differences in PAM and DAM gene expression. Through long-term tracking, we demonstrated microglial state plasticity. Lastly, we specifically depleted DAMs in demyelination, revealing their roles in disease recovery. Together, we provide a versatile genetic tool to characterize microglial states in CNS development and disease.

2-Desaza-annomontine (C81) impedes angiogenesis through reduced VEGFR2 expression derived from inhibition of CDC2-like kinases

Angiogenesis is a crucial process in the progression of various pathologies, like solid tumors, wet age-related macular degeneration, and chronic inflammation. Current anti-angiogenic treatments still have major drawbacks like limited efficacy in diseases that also rely on inflammation. Therefore, new anti-angiogenic approaches are sorely needed, and simultaneous inhibition of angiogenesis and inflammation is desirable. Here, we show that 2-desaza-annomontine (C81), a derivative of the plant alkaloid annomontine previously shown to inhibit endothelial inflammation, impedes angiogenesis by inhibiting CDC2-like kinases (CLKs) and WNT/β-catenin signaling. C81 reduced choroidal neovascularization in a laser-induced murine in vivo model, inhibited sprouting from vascular endothelial growth factor A (VEGF-A)-activated murine aortic rings ex vivo, and reduced angiogenesis-related activities of endothelial cells in multiple functional assays. This was largely phenocopied by CLK inhibitors and knockdowns, but not by inhibitors of the other known targets of C81. Mechanistically, CLK inhibition reduced VEGF receptor 2 (VEGFR2) mRNA and protein expression as well as downstream signaling. This was partly caused by a reduction of WNT/β-catenin pathway activity, as activating the pathway induced, while β-catenin knockdown impeded VEGFR2 expression. Surprisingly, alternative splicing of VEGFR2 was not detected. In summary, C81 and other CLK inhibitors could be promising compounds in the treatment of diseases that depend on angiogenesis and inflammation due to their impairment of both processes.

Combination of computational new approach methodologies for enhancing evidence of biological pathway conservation across species

The ability to predict which chemicals are of concern for environmental safety is dependent, in part, on the ability to extrapolate chemical effects across many species. This work investigated the complementary use of two computational new approach methodologies to support cross-species predictions of chemical susceptibility: the US Environmental Protection Agency Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool and Unilever's recently developed Genes to Pathways - Species Conservation Analysis (G2P-SCAN) tool. These stand-alone tools rely on existing biological knowledge to help understand chemical susceptibility and biological pathway conservation across species. The utility and challenges of these combined computational approaches were demonstrated using case examples focused on chemical interactions with peroxisome proliferator activated receptor alpha (PPARα), estrogen receptor 1 (ESR1), and gamma-aminobutyric acid type A receptor subunit alpha (GABRA1). Overall, the biological pathway information enhanced the weight of evidence to support cross-species susceptibility predictions. Through comparisons of relevant molecular and functional data gleaned from adverse outcome pathways (AOPs) to mapped biological pathways, it was possible to gain a toxicological context for various chemical-protein interactions. The information gained through this computational approach could ultimately inform chemical safety assessments by enhancing cross-species predictions of chemical susceptibility. It could also help fulfill a core objective of the AOP framework by potentially expanding the biologically plausible taxonomic domain of applicability of relevant AOPs.

An inducible genetic tool for tracking and manipulating specific microglial states in development and disease

Recent single-cell RNA sequencing studies have revealed distinct microglial states in development and disease. These include proliferative region-associated microglia (PAM) in developing white matter and disease-associated microglia (DAM) prevalent in various neurodegenerative conditions. PAM and DAM share a similar core gene signature and other functional properties. However, the extent of the dynamism and plasticity of these microglial states, as well as their functional significance, remains elusive, partly due to the lack of specific tools. Here, we report the generation of an inducible Cre driver line, Clec7a-CreERT2, designed to target PAM and DAM in the brain parenchyma. Utilizing this tool, we profile labeled cells during development and in several disease models, uncovering convergence and context-dependent differences in PAM/DAM gene expression. Through long-term tracking, we demonstrate surprising levels of plasticity in these microglial states. Lastly, we specifically depleted DAM in cuprizone-induced demyelination, revealing their roles in disease progression and recovery.

Global Analysis of Post-Translational Side-Chain Arginylation Using Pan-Arginylation Antibodies

Arginylation is a post-translational modification mediated by the arginyltransferase 1 (ATE1), which transfers the amino acid arginine to a protein or peptide substrate from a tRNA molecule. Initially, arginylation was thought to occur only on N-terminally exposed acidic residues, and its function was thought to be limited to targeting proteins for degradation. However, more recent data have shown that ATE1 can arginylate side chains of internal acidic residues in a protein without necessarily affecting metabolic stability. This greatly expands the potential targets and functions of arginylation, but tools for studying this process have remained limited. Here, we report the first global screen specifically for side-chain arginylation. We generate and validate "pan-arginylation" antibodies, which are designed to detect side-chain arginylation in any amino acid sequence context. We use these antibodies for immunoaffinity enrichment of side-chain arginylated proteins from wildtype and Ate1 knockout cell lysates. In this way, we identify a limited set of proteins that likely undergo ATE1-dependent side-chain arginylation and that are enriched in specific cellular roles, including translation, splicing, and the cytoskeleton.

Novel genetic variant associated with globoid cell leukodystrophy in a family of mixed breed dogs

BACKGROUND

Globoid cell leukodystrophy (GCL) is a fatal autosomal recessive disease caused by variants in the galactosylceramidase (GALC) gene. Two dog breed-specific variants are reported.

OBJECTIVES

Characterize the putatively causative GALC variant for GCL in a family of dogs and determine population allele frequency.

ANIMALS

Four related mixed-breed puppies with signs of neurologic disease were evaluated. Subsequently, 33 related dogs were tested for genetic markers for parentage and the identified GALC variant. Additional GALC genotyping was performed on 278 banked samples from various breeds.

METHODS

The 4 affected puppies had neurological exams and necropsies. DNA was isolated from blood samples. Variants in GALC were identified via Sanger sequencing. Parentage testing was performed using short tandem repeat markers. Prevalence of the GALC variant of interest was investigated in other breeds.

RESULTS

GCL was confirmed histopathologically. A novel missense variant in GALC (NC_006590.4:g.58893972G>A) was homozygous in all affected animals (n = 4). A recessive mode of inheritance was confirmed by parentage testing as was variant linkage with the phenotype (LOD = 3.36). Among the related dogs (n = 33), 3 dogs were homozygous and 7 heterozygous. The variant allele was not detected in screening 278 dogs from 5 breeds. The novel variant is either unique to this family or has an extremely low allele frequency in the general population.

CONCLUSIONS AND CLINICAL IMPORTANCE

A novel GALC variant was identified that likely explains GCL in this cohort. The identification of multiple causal variants for GCL in dogs is consistent with findings in humans.

Ovarian transcriptional response to Wolbachia infection in D. melanogaster in the context of between-genotype variation in gene expression

Wolbachia is a maternally transmitted endosymbiotic bacteria that infects a wide variety of arthropod and nematode hosts. The effects of Wolbachia on host biology are far-reaching and include changes in host gene expression. However, previous work on the host transcriptional response has generally been investigated in the context of a single host genotype. Thus, the relative effect of Wolbachia infection versus vs. host genotype on gene expression is unknown. Here, we explicitly test the relative roles of Wolbachia infection and host genotype on host gene expression by comparing the ovarian transcriptomes of 4 strains of Drosophila melanogaster (D. melanogaster) infected and uninfected with Wolbachia. Our data suggest that infection explains a small amount of transcriptional variation, particularly in comparison to variation in gene expression among strains. However, infection specifically affects genes related to cell cycle, translation, and metabolism. We also find enrichment of cell division and recombination processes among genes with infection-associated differential expression. Broadly, the transcriptomic changes identified in this study provide novel understanding of the relative magnitude of the effect of Wolbachia infection on gene expression in the context of host genetic variation and also point to genes that are consistently differentially expressed in response to infection among multiple genotypes.

Genes-to-Pathways Species Conservation Analysis: Enabling the Exploration of Conservation of Biological Pathways and Processes Across Species

The last two decades have witnessed a strong momentum toward integration of cell-based and computational approaches in safety assessments. This is fueling a global regulatory paradigm shift toward reduction and replacement of the use of animals in toxicity tests while promoting the use of new approach methodologies. The understanding of conservation of molecular targets and pathways provides an opportunity to extrapolate effects across species and ultimately to determine the taxonomic applicability domain of assays and biological effects. Despite the wealth of genome-linked data available, there is a compelling need for improved accessibility, while ensuring that it reflects the underpinning biology. We present the novel pipeline Genes-to-Pathways Species Conservation Analysis (G2P-SCAN) to further support understanding on cross-species extrapolation of biological processes. This R package extracts, synthetizes, and structures the data available from different databases, that is, gene orthologs, protein families, entities, and reactions, linked to human genes and respective pathways across six relevant model species. The use of G2P-SCAN enables the overall analysis of orthology and functional families to substantiate the identification of conservation and susceptibility at the pathway level. In the present study we discuss five case studies, demonstrating the validity of the developed pipeline and its potential use as species extrapolation support. We foresee this pipeline will provide valuable biological insights and create space for the use of mechanistically based data to inform potential species susceptibility for research and safety decision purposes. Environ Toxicol Chem 2023;42:1152-1166. © 2023 UNILEVER GLOBAL IP LTD. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Genomic distribution and polymorphism of G-quadruplex motifs occupying ovine promoters and enhancers

G-quadruplexes are non-canonical DNA structures that are formed in regions with short runs of guanines. During the last decade they have gained considerable attention due to their involvement in basic cellular processes, linking them to several physiological processes and pathological conditions. Regulation of gene transcription is among the crucial roles that G-quadruplexes play in the cells. Several ways in which these structures affect transcription have been described, both negatively and positively. Recently, G-quadruplexes have been shown to be implicated in the three-dimensional rearrangement of the chromosomes that enables the interaction of enhancers and gene promoters during regulation of transcription. Sheep is a species for which almost no G-quadruplex-related studies have been conducted and thus research on this species is kept out from the progress that has been made in the G-quadruplex field. In this context, we investigated the DNA sequences with potential to form G-quadruplexes (G4-motifs) in the ovine enhancers and promoters. We describe the distribution of G4-motifs within the regulatory regions which is shown to be enriched in G4-motifs in a way similar to other mammals. Furthermore, our data suggest that G4-motifs promote promoter-enhancer interactions in sheep. The single nucleotide polymorphisms colocalizing with promoter- and enhancer-associated ovine G4-motifs constitute a considerable pool of polymorphism and given the crucial role of these specific G4-motifs on regulation of transcription, we suggest this polymorphism as an interesting target for ovine genetic studies.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Primate protein-ligand interfaces exhibit significant conservation and unveil human-specific evolutionary drivers

Despite the vast phenotypic differences observed across primates, their protein products are largely similar to each other at the sequence level. We hypothesized that, since proteins accomplish all their functions via interactions with other molecules, alterations in the sites that participate in these interactions may be of critical importance. To uncover the extent to which these sites evolve across primates, we built a structurally-derived dataset of ~4,200 one-to-one orthologous sequence groups across 18 primate species, consisting of ~68,000 ligand-binding sites that interact with DNA, RNA, small molecules, ions, or peptides. Using this dataset, we identify functionally important patterns of conservation and variation within the amino acid residues that facilitate protein-ligand interactions across the primate phylogeny. We uncover that interaction sites are significantly more conserved than other sites, and that sites binding DNA and RNA further exhibit the lowest levels of variation. We also show that the subset of ligand-binding sites that do vary are enriched in components of gene regulatory pathways and uncover several instances of human-specific ligand-binding site changes within transcription factors. Altogether, our results suggest that ligand-binding sites have experienced selective pressure in primates and propose that variation in these sites may have an outsized effect on phenotypic variation in primates through pleiotropic effects on gene regulation.

lncRNA‑miRNA‑mRNA network in female offspring born from obese dams

Maternal obesity is associated with disturbance of lipid metabolism and obesity in offspring; however, the pathogenesis is still unclear. The present study elucidated the role of potential lipid metabolism-associated long non-coding RNA (lncRNA) and identified the pathways involved in mice born to obese dams. In the present study, maternal obesity was induced by feeding a high-fat diet for 10 weeks in female C57/BL6 mice, whereas control mice were fed a standard diet. All female mice mated with healthy male mice and were allowed to deliver spontaneously. The results demonstrated that female offspring from obese dams presented a tendency to become overweight in the first 8 weeks after birth; however, maternal obesity did not significantly alter the body weight of male offspring. RNA-sequencing analysis was performed on female offspring liver at 3 weeks old. Significantly dysregulated lncRNAs and downstream targets in female offspring liver were identified using bioinformatics analysis. lncRNA, microRNA (miRNA or miR) and mRNA expression levels in liver and AML12 cells were assessed using reverse transcription-quantitative PCR. A total of 8 upregulated and 17 downregulated lncRNAs were demonstrated in offspring from obese dams and lncRNA Lockd was indicated to be a key dysregulated lncRNA. Competing endogenous RNA (ceRNA) models suggested that the lncRNA Lockd/miR-582-5p/Elovl5 pathway was key for lipid metabolism in the liver of offspring from obese dams. Finally, small interfering RNA and miRNA inhibitor transfection was used to evaluate the ceRNA models in AML12 cells. Taken together, the results of the present study indicated that lncRNA Lockd-miR-582-5p-Elovl5 network may be disrupted in lipid metabolism and lead to obesity in the offspring of obese dams. This research will provide new insights into the molecular mechanism of obesity and lipid metabolism disorder.

Genome-wide analysis of annexin gene family in Schrenkiella parvula and Eutrema salsugineum suggests their roles in salt stress response

Annexins (Anns) play an important role in plant development, growth and responses to various stresses. Although Ann genes have been characterized in some plants, their role in adaptation mechanisms and tolerance to environmental stresses have not been studied in extremophile plants. In this study, Ann genes in Schrenkiella parvula and Eutrema salsugineum were identified using a genome-wide method and phylogenetic relationships, subcellular distribution, gene structures, conserved residues and motifs and also promoter prediction have been studied through bioinformatics analysis. We identified ten and eight encoding putative Ann genes in S. parvula and E. salsugineum genome respectively, which were divided into six subfamilies according to phylogenetic relationships. By observing conservation in gene structures and protein motifs we found that the majority of Ann members in two extremophile plants are similar. Furthermore, promoter analysis revealed a greater number of GATA, Dof, bHLH and NAC transcription factor binding sites, as well as ABRE, ABRE3a, ABRE4, MYB and Myc cis-acting elements in compare to Arabidopsis thaliana. To gain additional insight into the putative roles of candidate Ann genes, the expression of SpAnn1, SpAnn2 and SpAnn6 in S. parvula was studied in response to salt stress, which indicated that their expression level in shoot increased. Similarly, salt stress induced expression of EsAnn1, 5 and 7, in roots and EsAnn1, 2 and 5 in leaves of E. salsugineum. Our comparative analysis implies that both halophytes have different regulatory mechanisms compared to A. thaliana and suggest SpAnn2 gene play important roles in mediating salt stress.

Cryptic lineages and standing genetic variation across independent cane toad introductions

Widespread introduced species can be leveraged to investigate the genetic, ecological and adaptive processes underlying rapid evolution and range expansion, particularly the contributions of genetic diversity to adaptation. Rhinella marina, the cane toad, has been a focus of invasion biology for decades in Australia. However, their introduction history in North America is less clear. Here, we investigated the roles of introduction history and genetic diversity in establishment success of cane toads across their introduced range. We used reduced representation sequencing (ddRAD) to obtain 34,000 SNPs from 247 toads in native (French Guiana, Guyana, Ecuador, Panama, Texas) and introduced (Bermuda, southern Florida, northern Florida, Hawai'i, Puerto Rico) populations. Unlike all other cane toad introductions, we found that Florida populations were more closely related to native Central American lineages (R. horribilis), than to native Southern American lineages (R. marina). Furthermore, we found high levels of diversity and population structure in the native range, corroborating suggestions that R. marina is a species complex. We also found that introduced populations exhibit only slightly lower genetic diversity than native populations. Together with demographic analyses, this indicates founding populations of toads in Florida were larger than previously reported. Lastly, within R. marina, only one of 245 putatively adaptive SNPs showed fixed differences between native and introduced ranges, suggesting that putative selection in these introduced populations is based upon existing genetic variation. Our findings highlight the importance of genetic sequencing in understanding biological introductions and hint at the role of standing genetic variation in range expansion.

Lactate-mediated neural plasticity genes emerged during the evolution of memory systems

The ability to record experiences and learning is present to different degrees in several species; however, the complexity and diversity of memory processes are cognitive function features that differentiate humans from other species. Lactate has recently been discovered to act as a signaling molecule for neuronal plasticity linked to long-term memory. Because lactate is not only an energy substrate for neurons but also a signaling molecule for plasticity (Magistretti and Allaman in Nat Rev Neurosci 19:235-249, 2018. https://doi.org/10.1038/nrn.2018.19 ), it is of particular interest to understand how and when memory-related genes and lactate-mediated neural plasticity (LMNP) genes emerged and evolved in humans. To understand the evolutionary origin and processes of memory and LMNP genes, we first collected information on genes related to memory and LMNP from the literature and then conducted a comparative analysis of these genes. We found that the memory and LMNP genes have different origins, suggesting that these genes may have become established gradually in evolutionarily and functional terms and not at the same time. We also found that memory and LMNP systems have a similar evolutionary history, having been formed with the gradual participation of newly emerging genes throughout their evolution. We propose that the function of LMNP as a signaling process may be evolutionarily associated with memory systems through an unidentified system that is linked by 13 common genes between memory and LMNP gene sets. This study provides evolutionary insight into the possible relationship between memory and the LMNP systems that deepens our understanding of the evolution of memory systems.

Differential expression of immune-regulatory proteins C5AR1, CLEC4A and NLRP3 on peripheral blood mononuclear cells in early-stage non-small cell lung cancer patients

Changes in gene expression profiling of peripheral blood mononuclear cells (PBMC) appear to represent the host's response to the cancer cells via paracrine signaling. We speculated that protein expression on circulating T-lymphocytes represent T-lymphocyte trafficking before infiltration into the tumor microenvironment. The possibility of using protein expression on circulating T-lymphocytes as a biomarker to discriminate early-stage non-small cell lung cancer (NSCLC) was explored. Four independent PBMC gene expression microarray datasets (GSE12771, GSE13255, GSE20189 and GSE3934) were analyzed. We selected C5AR1, CLEC4A and NLRP3 based on their significant protein expression in tumor-infiltrating lymphocytes, but not in normal lymphoid tissue. A validation study using automated flow cytometry was conducted in 141 study participants including 76 treatment-naive early-stage non-small cell lung cancer patients (NSCLC), 12 individuals with non-malignant pulmonary diseases, and 53 healthy individuals. Median ratios of C5AR1, CLEC4A and NLRP3 specific antibody staining to CD3 positive cells in early-stage NSCLC patients compared to healthy controls were 0.014 [0-0.37] vs. 0.01 [0-0.07, p = 0.13], 0.03 [0-0.87] vs. 0.02 [0-0.13, p = 0.10] and 0.19 [0-0.60] vs. 0.09 [0.02-0.31, p < 0.0001], respectively. Median fluorescence intensity (MFI) of CD3⁺C5AR1⁺, CD3⁺CLEC4A⁺ and CD3⁺NLRP3⁺ expression in early-stage NSCLC patients compared to healthy volunteers was 185 [64.2-4801] vs. 107.5 [27-229, p < 0.0001], 91.2 [42.4-2355] vs. 71.25 [46.2-103, p = 0.0005], and 1585 [478-5224] vs. 758.5 [318-1976, p < 0.0001], respectively. NLRP3:CD3 ratio, CD3⁺C5AR1⁺, CD3⁺CLEC4A⁺ and CD3⁺NLRP3⁺ MFI were significantly higher in early-stage NSCLC than healthy volunteers with an area under the ROC curve of 0.69-0.76. The CD3⁺NLRP3⁺ MFI provided the most distinguishable expression at 71.5% sensitivity and 70% specificity. Furthermore, CD3⁺NLRP3⁺ MFI potentially discriminated between early-stage NSCLC from malignant-mimic inflammation and infection pulmonary disease. Further validation in various pulmonary inflammatory disease might be warranted. Our proof-of-principle findings strengthen the hypothesis that malignancies generate distinctive protein expression fingerprints on circulating T-lymphocytes.

Investigating Grapevine Red Blotch Virus Infection in Vitis vinifera L. cv. Cabernet Sauvignon Grapes: A Multi-Omics Approach

Grapevine red blotch virus (GRBV) is a recently identified virus. Previous research indicates primarily a substantial impact on berry ripening in all varieties studied. The current study analyzed grapes’ primary and secondary metabolism across grapevine genotypes and seasons to reveal both conserved and variable impacts to GRBV infection. Vitis vinifera cv. Cabernet Sauvignon (CS) grapevines grafted on two different rootstocks (110R and 420A) were analyzed in 2016 and 2017. Metabolite profiling revealed a considerable impact on amino acid and malate acid levels, volatile aroma compounds derived from the lipoxygenase pathway, and anthocyanins synthesized in the phenylpropanoid pathway. Conserved transcriptional responses to GRBV showed induction of auxin-mediated pathways and photosynthesis with inhibition of transcription and translation processes mainly at harvest. There was an induction of plant-pathogen interactions at pre-veraison, for all genotypes and seasons, except for CS 110R in 2017. Lastly, differential co-expression analysis revealed a transcriptional shift from metabolic synthesis and energy metabolism to transcription and translation processes associated with a virus-induced gene silencing transcript. This plant-derived defense response transcript was only significantly upregulated at veraison for all genotypes and seasons, suggesting a phenological association with disease expression and plant immune responses.

Mining of transcriptome identifies CD109 and LRP12 as possible biomarkers and deregulation mechanism of T cell receptor pathway in Acute Myeloid Leukemia

Acute Myeloid Leukemia (AML) is a heterogeneous disease with highest mortality compared to other types of leukemia. There is a need to find the gene abnormalities and mechanisms behind them due to their heterogenic nature. The present study is aimed to understand genes, pathways and biomarker proteins influenced by transcriptomic deregulation due to AML. Differentially expressed gene (DEG), protein-protein interaction network, gene ontology, KEGG pathway, variant analysis and secretome analyses were performed using different AML RNAseq datasets. A total of 655 DEGs including 291 up-regulated and 364 down-regulated genes, which were satisfied with a fold change of 1.5 were identified. Top hub genes for AML were identified as TP53, PTPRC and AKT1. This integrative bioinformatics approach revealed the deregulation of T Cell Receptor (TCR) pathway and altered immune response related genes. The survival analysis revealed the associated deregulation of multiple TCR pathway related genes. Variant analysis identified the benign and likely benign nature of many important target genes and markers screened, which were found to have an important role in the progression of AML. DEGs and secretome analysis found out a set of seven molecules represents potential biomarkers for AML. In vitro analytical validation showed overexpression pattern of CD109 and LRP12 in AML cell line and HL-60 cells than the normal human bone marrow-derived stromal cell line HS-5. Here we report first time for CD109 and LRP12 as a possible biomarkers for the diagnostic significance. Amino acid substitutions detected by variant analysis and deregulation of immune checkpoint molecules revealed their role in reducing immune response and inability to fight cancer cells. In conclusion, this study highlights the possibility of new biomarkers for AML and the mechanism of decrease in immune response due to the downregulation of co-stimulatory immune molecules, which needs further clinical validation investigations.

•

Using RNA-seq data of AML patients, two biomarkers including CD109 and LRP12 for the diagnostic significance were identified based on DEGs, GO/KEGG, and PPI network analysis.

•

The transcriptome mining unmasked the complexity of gene alterations in AML by identifying immune response related genes deregulation and significance of TCR signalling.

•

Several genes were identified as AML hub genes by network analysis, variant analysis identified non-synonymous variants in co-stimulatory checkpoint targets and the co-inhibitory targets.

Regio-Specific N-Glycome and N-Glycoproteome Map of the Elderly Human Brain With and Without Alzheimer's Disease

The proteins in the cell membrane of the brain are modified by glycans in highly interactive regions. The glycans and glycoproteins are involved in cell-cell interactions that are of fundamental importance to the brain. In this study, the comprehensive N-glycome and N-glycoproteome of the brain were determined in 11 functional brain regions, some of them known to be affected with the progression of Alzheimer's disease. N-glycans throughout the regions were generally highly branched and highly sialofucosylated. Regional variations were also found with regard to the glycan types including high mannose and complex-type structures. Glycoproteomic analysis identified the proteins that differed in glycosylation in the various regions. To obtain the broader representation of glycan compositions, four subjects with two in their 70s and two in their 90s representing two Alzheimer's disease subjects, one hippocampal sclerosis subject, and one subject with no cognitive impairment were analyzed. The four subjects were all glycomically mapped across 11 brain regions. Marked differences in the glycomic and glycoproteomic profiles were observed between the samples.

Gene Expression and DNA Methylation in Human Papillomavirus Positive and Negative Head and Neck Squamous Cell Carcinomas

High-risk human papillomaviruses (HPV) are important agents, responsible for a large percentage of the 745,000 cases of head and neck squamous cell carcinomas (HNSCC), which were identified worldwide in 2020. In addition to being virally induced, tobacco and heavy alcohol consumption are believed to cause DNA damage contributing to the high number of HNSCC cases. Gene expression and DNA methylation differ between HNSCC based on HPV status. We used publicly available gene expression and DNA methylation profiles from the Cancer Genome Atlas and compared HPV positive and HPV negative HNSCC groups. We used differential gene expression analysis, differential methylation analysis, and a combination of these two analyses to identify the differences. Differential expression analysis identified 1854 differentially expressed genes, including PCNA, TNFRSF14, TRAF1, TRAF2, BCL2, and BIRC3. SYCP2 was identified as one of the top deregulated genes in the differential methylation analysis and in the combined differential expression and methylation analyses. Additionally, pathway and ontology analyses identified the extracellular matrix and receptor interaction pathway as the most altered between HPV negative and HPV positive HNSCC groups. Combining gene expression and DNA methylation can help in elucidating the genes involved in HPV positive HNSCC tumorigenesis, such as SYCP2 and TAF7L.

Quantitative proteomic analytic approaches to identify metabolic changes in the medial prefrontal cortex of rats exposed to space radiation

NASA’s planned mission to Mars will result in astronauts being exposed to ∼350 mSv/yr of Galactic Cosmic Radiation (GCR). A growing body of data from ground-based experiments indicates that exposure to space radiation doses (approximating those that astronauts will be exposed to on a mission to Mars) impairs a variety of cognitive processes, including cognitive flexibility tasks. Some studies report that 33% of individuals may experience severe cognitive impairment. Translating the results from ground-based rodent studies into tangible risk estimates for astronauts is an enormous challenge, but it would be germane for NASA to use the vast body of data from the rodent studies to start developing appropriate countermeasures, in the expectation that some level of space radiation (SR) -induced cognitive impairment could occur in astronauts. While some targeted studies have reported radiation-induced changes in the neurotransmission properties and/or increased neuroinflammation within space radiation exposed brains, there remains little information that can be used to start the development of a mechanism-based countermeasure strategy. In this study, we have employed a robust label-free mass spectrometry (MS) -based untargeted quantitative proteomic profiling approach to characterize the composition of the medial prefrontal cortex (mPFC) proteome in rats that have been exposed to 15 cGy of 600 MeV/n28Si ions. A variety of analytical techniques were used to mine the generated expression data, which in such studies is typically hampered by low and variable sample size. We have identified several pathways and proteins whose expression alters as a result of space radiation exposure, including decreased mitochondrial function, and a further subset of proteins differs in rats that have a high level of cognitive performance after SR exposure in comparison with those that have low performance levels. While this study has provided further insight into how SR impacts upon neurophysiology, and what adaptive responses can be invoked to prevent the emergence of SR-induced cognitive impairment, the main objective of this paper is to outline strategies that can be used by others to analyze sub-optimal data sets and to identify new information.

Genome-wide DNA Methylation Differences in Nonfunctioning Pituitary Adenomas With and Without Postsurgical Progression

CONTEXT

Tumor progression in surgically treated patients with nonfunctioning pituitary adenomas (NFPAs) is associated with excess mortality. Reliable biomarkers allowing early identification of tumor progression are missing.

OBJECTIVE

To explore DNA methylation patterns associated with tumor progression in NFPA patients.

METHODS

This case-controlled exploratory trial at a university hospital studied patients who underwent surgery for NFPA that had immunohistochemical characteristics of a gonadotropinoma. Cases included patients requiring reintervention due to tumor progression (reintervention group, n = 26) and controls who had a postoperative residual tumor without tumor progression for at least 5 years (radiologically stable group, n = 17). Genome-wide methylation data from each tumor sample were analyzed using the Infinium MethylationEPIC BeadChip platform.

RESULTS

The analysis showed that 605 CpG positions were significantly differently methylated (differently methylated positions, DMPs) between the patient groups (false discovery rate adjusted P value < 0.05, beta value > 0.2), mapping to 389 genes. The largest number of DMPs were detected in the genes NUP93 and LGALS1. The 3 hypomethylated DMPs and the 3 hypermethylated DMPs with the lowest P values were all significantly (P < 0.05) and individually associated with reintervention-free survival. One of the hypermethylated DMPs with the lowest P value was located in the gene GABRA1.

CONCLUSION

In this exploratory study, DNA methylation patterns in NFPA patients were associated with postoperative tumor progression requiring reintervention. The DMPs included genes that have been previously associated with tumor development. Our study is a step toward finding epigenetic signatures to predict tumor progression in patients with NFPA.

Exploring the expression patterns of palmitoylating and de-palmitoylating enzymes in the mouse brain using the curated RNA-seq database BrainPalmSeq

Protein S-palmitoylation is a reversible post-translational lipid modification that plays a critical role in neuronal development and plasticity, while dysregulated S-palmitoylation underlies a number of severe neurological disorders. Dynamic S-palmitoylation is regulated by a large family of ZDHHC palmitoylating enzymes, their accessory proteins, and a small number of known de-palmitoylating enzymes. Here, we curated and analyzed expression data for the proteins that regulate S-palmitoylation from publicly available RNAseq datasets, providing a comprehensive overview of their distribution in the mouse nervous system. We developed a web-tool that enables interactive visualization of the expression patterns for these proteins in the nervous system (http://brainpalmseq.med.ubc.ca/), and explored this resource to find region and cell-type specific expression patterns that give insight into the function of palmitoylating and de-palmitoylating enzymes in the brain and neurological disorders. We found coordinated expression of ZDHHC enzymes with their accessory proteins, de-palmitoylating enzymes and other brain-expressed genes that included an enrichment of S-palmitoylation substrates. Finally, we utilized ZDHHC expression patterns to predict and validate palmitoylating enzyme-substrate interactions.

Panzootic chytrid fungus exploits diverse amphibian host environments through plastic infection strategies

While some pathogens are limited to single species, others can colonize many hosts, likely contributing to the emergence of novel disease outbreaks. Despite this biodiversity threat, traits associated with host niche expansions are not well understood in multi-host pathogens. Here, we aimed to uncover functional machinery driving multi-host invasion by focusing on Batrachochytrium dendrobatidis (Bd), a pathogen that infects the skin of hundreds of amphibians worldwide. We performed a meta-analysis of Bd gene expression using data from published infection experiments and newly generated profiles. We analyzed Bd transcriptomic landscapes across the skin of 14 host species, reconstructed Bd isolates phylogenetic relationships, and inferred the origin and evolutionary history of differentially expressed genes under a phylogenetic framework comprising other 12 zoosporic fungi. Bd displayed plastic infection strategies when challenged by hosts with different disease susceptibility. Our analyses identified sets of differentially expressed genes under host environments with similar infection outcome. We stressed nutritional immunity and gene silencing as important processes required to overcome challenging skin environments in less susceptible hosts. Overall, Bd genes expressed during amphibian skin exploitation have arisen mainly via gene duplications with great family expansions, increasing the gene copy events previously described for this fungal species. Finally, we provide a comprehensive gene dataset that can be used to further examine eco-evolutionary hypotheses for this host-pathogen system. Our study supports the idea that host environments exert contrasting selective pressures, such that gene expression plasticity could be one of the evolutionary keys leading to the success of multi-host pathogens.

In vitro ve in silico analizi ile metforminin meme tümörü hücrelerinde protein profili üzerindeki etkinliği

Aim: This study aimed to uncover the varieties in protein profiles of Met in breast tumor (BT) cells by assessment of in vitro and in silico analysis. Materials and Methods: Here, the cells obtained from mastectomy patients were cultured, the effective Met-dose was determined as 25 mM through cell viability and BrdU tests. Protein identification in the breast tumor cells was implemented by employing LC-MS/MS technology. Results: The expression of SSR3, THAP3, FTH1, NEFM, ANP32A, ANP32B, KRT7 proteins was significantly decreased whereas the GARS protein increased in the 25 mM Met group compared to the Non-Met (0 mM) control group. In silico analysis, we analyzed the probable interactions of all these proteins with each other and other proteins, to evaluate the analysis of the larger protein network, and which metabolic pathway proteins are involved in. Conclusion: The stated proteomics analysis in our study proposes a better understanding of the prognosis of breast cancer and future studies to investigate the effect of metformin in this field on proteomic pathways in other sorts of cancer.

Dysregulated Gene Expression in Lymphoblasts from Parkinson’s Disease

Parkinson’s disease is the second largest neurodegenerative disease worldwide and is caused by a combination of genetics and environment. It is characterized by the death of neurons in the substantia nigra of the brain but is not solely a disease of the brain, as it affects multiple tissues and organs. Studying Parkinson’s disease in accessible tissues such as skin and blood has increased our understanding of the disease’s pathogenesis. Here, we used lymphoblast cell lines generated from Parkinson’s disease patient and healthy age- and sex-matched control groups and obtained their whole-cell transcriptomes and proteomes. Our analysis revealed, in both the transcriptomes and the proteomes of PD cells, a global downregulation of genes involved in protein synthesis, as well as the upregulation of immune processes and sphingolipid metabolism. In contrast, we discovered an uncoupling of mRNA and protein expression in processes associated with mitochondrial respiration in the form of a general downregulation in associated transcripts and an upregulation in proteins. Complex V was different to the other oxidative phosphorylation complexes in that the levels of its associated transcripts were also lower, but the levels of their encoded polypeptides were not elevated. This may suggest that further layers of regulation specific to Complex V are in play.

Multi-Omics Analysis of Glioblastoma and Glioblastoma Cell Line: Molecular Insights Into the Functional Role of GPR56 and TG2 in Mesenchymal Transition

G protein-coupled receptor 56 (GPR56/ADGRG1) is an adhesion GPCR with an essential role in brain development and cancer. Elevated expression of GPR56 was observed in the clinical specimens of Glioblastoma (GBM), a highly invasive primary brain tumor. However, we found the expression to be variable across the specimens, presumably due to the intratumor heterogeneity of GBM. Therefore, we re-examined GPR56 expression in public domain spatial gene expression data and single-cell expression data for GBM, which revealed that GPR56 expression was high in cellular tumors, infiltrating tumor cells, and proliferating cells, low in microvascular proliferation and peri-necrotic areas of the tumor, especially in hypoxic mesenchymal-like cells. To gain a better understanding of the consequences of GPR56 downregulation in tumor cells and other molecular changes associated with it, we generated a sh-RNA-mediated GPR56 knockdown in the GBM cell line U373 and performed transcriptomics, proteomics, and phospho-proteomics analysis. Our analysis revealed enrichment of gene signatures, pathways, and phosphorylation of proteins potentially associated with mesenchymal (MES) transition in the tumor and concurrent increase in cell invasion and migration behavior of the GPR56 knockdown GBM cells. Interestingly, our analysis also showed elevated expression of Transglutaminase 2 (TG2) - a known interactor of GPR56, in the knockdown cells. The inverse expression of GPR56 and TG2 was also observed in intratumoral, spatial gene expression data for GBM and in GBM cell lines cultured in vitro under hypoxic conditions. Integrating all these observations, we propose a putative functional link between the inverse expression of the two proteins, the hypoxic niche and the mesenchymal status in the tumor. Hypoxia-induced downregulation of GPR56 and activation of TG2 may result in a network of molecular events that contribute to the mesenchymal transition of GBM cells, and we propose a putative model to explain this functional and regulatory relationship of the two proteins.

Brucella activates the host RIDD pathway to subvert BLOS1-directed immune defense

The phagocytosis and destruction of pathogens in lysosomes constitute central elements of innate immune defense. Here, we show that Brucella, the causative agent of brucellosis, the most prevalent bacterial zoonosis globally, subverts this immune defense pathway by activating regulated IRE1α-dependent decay (RIDD) of Bloc1s1 mRNA encoding BLOS1, a protein that promotes endosome-lysosome fusion. RIDD-deficient cells and mice harboring a RIDD-incompetent variant of IRE1α were resistant to infection. Inactivation of the Bloc1s1 gene impaired the ability to assemble BLOC-1-related complex (BORC), resulting in differential recruitment of BORC-related lysosome trafficking components, perinuclear trafficking of Brucella-containing vacuoles (BCVs), and enhanced susceptibility to infection. The RIDD-resistant Bloc1s1 variant maintains the integrity of BORC and a higher-level association of BORC-related components that promote centrifugal lysosome trafficking, resulting in enhanced BCV peripheral trafficking and lysosomal destruction, and resistance to infection. These findings demonstrate that host RIDD activity on BLOS1 regulates Brucella intracellular parasitism by disrupting BORC-directed lysosomal trafficking. Notably, coronavirus murine hepatitis virus also subverted the RIDD-BLOS1 axis to promote intracellular replication. Our work establishes BLOS1 as a novel immune defense factor whose activity is hijacked by diverse pathogens.

Epicardium-derived cells organize through tight junctions to replenish cardiac muscle in salamanders

The contribution of the epicardium, the outermost layer of the heart, to cardiac regeneration has remained controversial due to a lack of suitable analytical tools. By combining genetic marker-independent lineage-tracing strategies with transcriptional profiling and loss-of-function methods, we report here that the epicardium of the highly regenerative salamander species Pleurodeles waltl has an intrinsic capacity to differentiate into cardiomyocytes. Following cryoinjury, CLDN6+ epicardium-derived cells appear at the lesion site, organize into honeycomb-like structures connected via focal tight junctions and undergo transcriptional reprogramming that results in concomitant differentiation into de novo cardiomyocytes. Ablation of CLDN6+ differentiation intermediates as well as disruption of their tight junctions impairs cardiac regeneration. Salamanders constitute the evolutionarily closest species to mammals with an extensive ability to regenerate heart muscle and our results highlight the epicardium and tight junctions as key targets in efforts to promote cardiac regeneration.

Genomic Resources for Salminus brasiliensis

The Neotropical region bears the most diverse freshwater fish fauna on the planet and is the stage for dramatic conservation struggles. Initiatives aiming for conservation of a single emblematic fish, a flagship species, to which different onlookers relate on a cultural/personal level, holds promise towards engagement and conservation actions benefiting whole biological communities and ecosystems. Here, we present the first comprehensive genomic resources for Salminus brasiliensis, a potential flagship Neotropical species. This fish faces pressing conservation issues, as well as taxonomic uncertainty, being a main species relevant to angling and commercial fisheries. We make available 178 million Illumina paired-end reads, 90 bases long, comprising 16 Gb (≈15X coverage) of filtered data, obtained from a primary genomic library of 500-bp fragments. We present the first de novo genomic assembly for S. brasiliensis, with ∼1 Gb (N₅₀ = 10,889), as well as the coding genome annotation of 12,962 putative genes from assembled genomic fragments over 10 kb, most of which could be identified from the Ostariophysi GenBank database. We also provide a genome-wide panel for more than 80,000 predicted microsatellite loci for low-cost, fast and abundant DNA marker development for this species. A total of 47, among 52 candidates, empirically assayed microsatellites were confirmed as polymorphic in this fish. All genomic data produced for S. brasiliensis is hereby made publicly accessible. With the disclosure of these results, we intend to foster general biology studies and to provide tools to be applied immediately in conservation and aquaculture in this candidate flagship Neotropical species.

SHOOT: phylogenetic gene search and ortholog inference

Determining the evolutionary relationships between genes is fundamental to comparative biological research. Here, we present SHOOT. SHOOT searches a user query sequence against a database of phylogenetic trees and returns a tree with the query sequence correctly placed within it. We show that SHOOT performs this analysis with comparable speed to a BLAST search. We demonstrate that SHOOT phylogenetic placements are as accurate as conventional tree inference, and it can identify orthologs with high accuracy. In summary, SHOOT is a fast and accurate tool for phylogenetic analyses of novel query sequences. It is available online at www.shoot.bio.

Autosomal Recessive Cerebellar Atrophy and Spastic Ataxia in Patients With Pathogenic Biallelic Variants in GEMIN5

The hereditary ataxias are a heterogenous group of disorders with an increasing number of causative genes being described. Due to the clinical and genetic heterogeneity seen in these conditions, the majority of such individuals endure a diagnostic odyssey or remain undiagnosed. Defining the molecular etiology can bring insights into the responsible molecular pathways and eventually the identification of therapeutic targets. Here, we describe the identification of biallelic variants in the GEMIN5 gene among seven unrelated families with nine affected individuals presenting with spastic ataxia and cerebellar atrophy. GEMIN5, an RNA-binding protein, has been shown to regulate transcription and translation machinery. GEMIN5 is a component of small nuclear ribonucleoprotein (snRNP) complexes and helps in the assembly of the spliceosome complexes. We found that biallelic GEMIN5 variants cause structural abnormalities in the encoded protein and reduce expression of snRNP complex proteins in patient cells compared with unaffected controls. Finally, knocking out endogenous Gemin5 in mice caused early embryonic lethality, suggesting that Gemin5 expression is crucial for normal development. Our work further expands on the phenotypic spectrum associated with GEMIN5-related disease and implicates the role of GEMIN5 among patients with spastic ataxia, cerebellar atrophy, and motor predominant developmental delay.

Tgfbr2 in Dental Pulp Cells Guides Neurite Outgrowth in Developing Teeth

Transforming growth factor β (TGFβ) plays an important role in tooth morphogenesis and mineralization. During postnatal development, the dental pulp (DP) mesenchyme secretes neurotrophic factors that guide trigeminal nerve fibers into and throughout the DP. This process is tightly linked with dentin formation and mineralization. Our laboratory established a mouse model in which Tgfbr2 was conditionally deleted in DP mesenchyme using an Osterix promoter-driven Cre recombinase (Tgfbr2 ^cko ). These mice survived postnatally with significant defects in bones and teeth, including reduced mineralization and short roots. Hematoxylin and eosin staining revealed reduced axon-like structures in the mutant mice. Reporter imaging demonstrated that Osterix-Cre activity within the tooth was active in the DP and derivatives, but not in neuronal afferents. Immunofluorescence staining for β3 tubulin (neuronal marker) was performed on serial cryosections from control and mutant molars on postnatal days 7 and 24 (P7, P24). Confocal imaging and pixel quantification demonstrated reduced innervation in Tgfbr2 ^cko first molars at both stages compared to controls, indicating that signals necessary to promote neurite outgrowth were disrupted by Tgfbr2 deletion. We performed mRNA-Sequence (RNA-Seq) and gene onotology analyses using RNA from the DP of P7 control and mutant mice to investigate the pathways involved in Tgfbr2-mediated tooth development. These analyses identified downregulation of several mineralization-related and neuronal genes in the Tgfbr2 ^cko DP compared to controls. Select gene expression patterns were confirmed by quantitative real-time PCR and immunofluorescence imaging. Lastly, trigeminal neurons were co-cultured atop Transwell filters overlying primary Tgfbr2 ^f/f DP cells. Tgfbr2 in the DP was deleted via Adenovirus-expressed Cre recombinase. Confocal imaging of axons through the filter pores showed increased axonal sprouting from neurons cultured with Tgfbr2-positive DP cells compared to neurons cultured alone. Axon sprouting was reduced when Tgfbr2 was knocked down in the DP cells. Immunofluorescence of dentin sialophosphoprotein in co-cultured DP cells confirmed reduced mineralization potential in cells with Tgfbr2 deletion. Both our proteomics and RNA-Seq analyses indicate that axonal guidance cues, particularly semaphorin signaling, were disrupted by Tgfbr2 deletion. Thus, Tgfbr2 in the DP mesenchyme appears to regulate differentiation and the cells' ability to guide neurite outgrowth during tooth mineralization and innervation.

Immune Checkpoint Blockade Augments Changes Within Oncolytic Virus-induced Cancer MHC-I Peptidome, Creating Novel Antitumor CD8 T Cell Reactivities

The combination cancer immunotherapies with oncolytic virus (OV) and immune checkpoint blockade (ICB) reinstate otherwise dysfunctional antitumor CD8 T cell responses. One major mechanism that aids such reinstatement of antitumor CD8 T cells involves the availability of new class I major histocompatibility complex (MHC-I)-bound tumor epitopes following therapeutic intervention. Thus, therapy-induced changes within the MHC-I peptidome hold the key to understanding the clinical implications for therapy-reinstated CD8 T cell responses. Here, using mass spectrometry-based immuno-affinity methods and tumor-bearing animals treated with OV and ICB (alone or in combination), we captured the therapy-induced alterations within the tumor MHC-I peptidome, which were then tested for their CD8 T cell response-stimulating activity. We found that the oncolytic reovirus monotherapy drives up- as well as downexpression of tumor MHC-I peptides in a cancer type and oncolysis susceptibility dependent manner. Interestingly, the combination of reovirus + ICB results in higher numbers of differentially expressed MHC-I-associated peptides (DEMHCPs) relative to either monotherapies. Most importantly, OV+ICB-driven DEMHCPs contain biologically active epitopes that stimulate interferon-gamma responses in cognate CD8 T cells, which may mediate clinically desired antitumor attack and cancer immunoediting. These findings highlight that the therapy-induced changes to the MHC-I peptidome contribute toward the reinstated antitumor CD8 T cell attack established following OV + ICB combination cancer immunotherapy.

Genomic Diversity Profiling and Breed-Specific Evolutionary Signatures of Selection in Arunachali Yak

Arunachali yak, the only registered yak breed of India, is crucial for the economic sustainability of pastoralist Monpa community. This study intended to determine the genomic diversity and to identify signatures of selection in the breed. Previously available double digest restriction-site associated DNA (ddRAD) sequencing data of Arunachali yak animals was processed and 99,919 SNPs were considered for further analysis. The genomic diversity profiled based on nucleotide diversity, π (π = 0.041 in 200 bp windows), effective population size, Ne (Ne = 83) and Runs of homozygosity (ROH) (predominance of shorter length ROHs) was found to be optimum. Subsequently, 207 regions were identified to be under selective sweeps through de-correlated composite of multiple signals (DCMS) statistic which combined three individual test statistics viz. π, Tajima’s D and |iHS| in non-overlapping 100 kb windows. Mapping of these regions revealed 611 protein-coding genes including KIT, KITLG, CDH12, FGG, FGA, FGB, PDGFRA, PEAR1, STXBP3, olfactory receptor genes (OR5K3, OR5H6 and OR1E1) and taste receptor genes (TAS2R1, TAS2R3 and TAS2R4). Functional annotation highlighted that biological processes like platelet aggregation and sensory perception were the most overrepresented and the associated regions could be considered as breed-specific signatures of selection in Arunachali yak. These findings point towards evolutionary role of natural selection in environmental adaptation of Arunachali yak population and provide useful insights for pursuing genome-wide association studies in future.

Genomic landscape, polymorphism and possible LINE-associated delivery of G-quadruplex motifs in the bovine genes

G-Quadruplex structures are non-B DNA structures that occur in regions carrying short runs of guanines. They are implicated in several biological processes including transcription, translation, replication and telomere maintenance as well as in several pathological conditions like cancer and thus they have gained the attention of the scientific community. The rise of the -omics era significantly affected the G-quadruplex research and the genome-wide characterization of G-Quadruplexes has been rendered a necessary first step towards applying genomics approaches for their study. While in human and several model organisms there is a considerable number of works studying genome-wide the DNA motifs with potential to form G-quadruplexes (G4-motifs), there is a total absence of any similar studies regarding livestock animals. The objectives of the present study were to provide a detailed characterization of the bovine genic G4-motifs' distribution and properties and to suggest a possible mechanism for the delivery of G4 motifs in the genes. Our data indicate that the distribution of G4-motifs within bovine genes and the annotation of said genes to Gene Ontology terms are similar to what is already shown for other organisms. By investigating their structural characteristics and polymorphism, it is obvious that the overall stability of the putative quadruplex structures is in line with the current notion in the G4 field. Similarly to human, the bovine G4-motifs are overrepresented in specific LINE repeat elements, the L1_BTs in the case of cattle. We highlight the potential role of these elements as vehicles for delivery of G4 motifs in the introns of the bovine genes. Lastly, it seems that a basis exists for connecting traits of agricultural importance to the genetic variation of G4 motifs, thus, the value of cattle as an interesting new model organism for G4-related genetic studies might be worth to be investigated.

Comparative proteomic analysis of nuclear and cytoplasmic compartments in human cardiac progenitor cells

Clinical trials evaluating cardiac progenitor cells (CPC) demonstrated feasibility and safety, but no clear functional benefits. Therefore a deeper understanding of CPC biology is warranted to inform strategies capable to enhance their therapeutic potential. Here we have defined, using a label-free proteomic approach, the differential cytoplasmic and nuclear compartments of human CPC (hCPC). Global analysis of cytoplasmic repertoire in hCPC suggested an important hypoxia response capacity and active collagen metabolism. In addition, comparative analysis of the nuclear protein compartment identified a significant regulation of a small number of proteins in hCPC versus human mesenchymal stem cells (hMSC). Two proteins significantly upregulated in the hCPC nuclear compartment, IL1A and IMP3, showed also a parallel increase in mRNA expression in hCPC versus hMSC, and were studied further. IL1A, subjected to an important post-transcriptional regulation, was demonstrated to act as a dual-function cytokine with a plausible role in apoptosis regulation. The knockdown of the mRNA binding protein (IMP3) did not negatively impact hCPC viability, but reduced their proliferation and migration capacity. Analysis of a panel of putative candidate genes identified HMGA2 and PTPRF as IMP3 targets in hCPC. Therefore, they are potentially involved in hCPC proliferation/migration regulation.

K2P18.1 translates T cell receptor signals into thymic regulatory T cell development

It remains largely unclear how thymocytes translate relative differences in T cell receptor (TCR) signal strength into distinct developmental programs that drive the cell fate decisions towards conventional (Tconv) or regulatory T cells (Treg). Following TCR activation, intracellular calcium (Ca²⁺) is the most important second messenger, for which the potassium channel K_2P18.1 is a relevant regulator. Here, we identify K_2P18.1 as a central translator of the TCR signal into the thymus-derived Treg (tTreg) selection process. TCR signal was coupled to NF-κB-mediated K_2P18.1 upregulation in tTreg progenitors. K_2P18.1 provided the driving force for sustained Ca²⁺ influx that facilitated NF-κB- and NFAT-dependent expression of FoxP3, the master transcription factor for Treg development and function. Loss of K_2P18.1 ion-current function induced a mild lymphoproliferative phenotype in mice, with reduced Treg numbers that led to aggravated experimental autoimmune encephalomyelitis, while a gain-of-function mutation in K_2P18.1 resulted in increased Treg numbers in mice. Our findings in human thymus, recent thymic emigrants and multiple sclerosis patients with a dominant-negative missense K_2P18.1 variant that is associated with poor clinical outcomes indicate that K_2P18.1 also plays a role in human Treg development. Pharmacological modulation of K_2P18.1 specifically modulated Treg numbers in vitro and in vivo. Finally, we identified nitroxoline as a K_2P18.1 activator that led to rapid and reversible Treg increase in patients with urinary tract infections. Conclusively, our findings reveal how K_2P18.1 translates TCR signals into thymic T cell fate decisions and Treg development, and provide a basis for the therapeutic utilization of Treg in several human disorders.

PANTHER: Making genome-scale phylogenetics accessible to all

Phylogenetics is a powerful tool for analyzing protein sequences, by inferring their evolutionary relationships to other proteins. However, phylogenetics analyses can be challenging: they are computationally expensive and must be performed carefully in order to avoid systematic errors and artifacts. Protein Analysis THrough Evolutionary Relationships (PANTHER; http://pantherdb.org) is a publicly available, user-focused knowledgebase that stores the results of an extensive phylogenetic reconstruction pipeline that includes computational and manual processes and quality control steps. First, fully reconciled phylogenetic trees (including ancestral protein sequences) are reconstructed for a set of "reference" protein sequences obtained from fully sequenced genomes of organisms across the tree of life. Second, the resulting phylogenetic trees are manually reviewed and annotated with function evolution events: inferred gains and losses of protein function along branches of the phylogenetic tree. Here, we describe in detail the current contents of PANTHER, how those contents are generated, and how they can be used in a variety of applications. The PANTHER knowledgebase can be downloaded or accessed via an extensive API. In addition, PANTHER provides software tools to facilitate the application of the knowledgebase to common protein sequence analysis tasks: exploring an annotated genome by gene function; performing "enrichment analysis" of lists of genes; annotating a single sequence or large batch of sequences by homology; and assessing the likelihood that a genetic variant at a particular site in a protein will have deleterious effects.

The novel anti-inflammatory activity of mcIRBP from Momordica charantia is associated with the improvement of diabetic nephropathy

Diabetic nephropathy is an inflammatory immune disorder accompanying diabetes.

The splicing factor DHX38/PRP16 is required for ovarian clear cell carcinoma tumorigenesis, as revealed by a CRISPR-Cas9 screen

Certain cancers, such as ovarian clear cell carcinoma (OCCC), display high levels of genetic variation between patients, making it difficult to develop effective therapies. In order to identify novel genes critical to OCCC growth, we carried out a comprehensive CRISPR‐Cas9 knockout screen against cell growth using an OCCC cell line and a normal ovarian surface epithelium cell line. We identified the gene encoding DHX38/PRP16, an ATP‐dependent RNA helicase involved in splicing, as critical for the growth and tumorigenesis of OCCC. DHX38/PRP16 knockdown in OCCC cells, but not normal cells, induces apoptosis and impairs OCCC tumorigenesis in a mouse model. Our results suggest that DHX38/PRP16 may play a role in OCCC tumorigenesis and could potentially be a promising therapeutic target.

Age-Related Changes of Gene Expression Profiles in Drosophila

An individual’s gene expression profile changes throughout their life. This change in gene expression is shaped by differences in physiological needs and functions between the younger and older organism. Despite intensive studies, the aging process is not fully understood, and several genes involved in this process may remain to be identified. Here we report a transcriptomic analysis of Drosophila melanogaster using microarrays. We compared the expression profiles of two-day-old female adult flies with those of 45-day-old flies. We identified 1184 genes with pronounced differences in expression level between young and old age groups. Most genes involved in muscle development/maintenance that display different levels of expression with age were downregulated in older flies. Many of these genes contributed to sarcomere formation and function. Several of these genes were functionally related to direct and indirect flight muscles; some of them were exclusively expressed in these muscles. Conversely, several genes involved in apoptosis processes were upregulated in aging flies. In addition, several genes involved in resistance to toxic chemicals were upregulated in aging flies, which is consistent with a global upregulation of the defense response system in aging flies. Finally, we randomly selected 12 genes among 232 genes with unknown function and generated transgenic flies expressing recombinant proteins fused with GFP protein to determine their subcellular expression. We also found that the knockdown of some of those 12 genes can affect the lifespan of flies.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

Arabidopsis bioinformatics: tools and strategies

The sequencing of the Arabidopsis thaliana genome 21 years ago ushered in the genomics era for plant research. Since then, an incredible variety of bioinformatic tools permit easy access to large repositories of genomic, transcriptomic, proteomic, epigenomic and other '-omic' data. In this review, we cover some more recent tools (and highlight the 'classics') for exploring such data in order to help formulate quality, testable hypotheses, often without having to generate new experimental data. We cover tools for examining gene expression and co-expression patterns, undertaking promoter analyses and gene set enrichment analyses, and exploring protein-protein and protein-DNA interactions. We will touch on tools that integrate different data sets at the end of the article.