Using GOstats to test gene lists for GO term association

CHD8 haploinsufficiency links autism to transient alterations in excitatory and inhibitory trajectories

Mutations in the chromodomain helicase DNA-binding 8 (CHD8) gene are a frequent cause of autism spectrum disorder (ASD). While its phenotypic spectrum often encompasses macrocephaly, implicating cortical abnormalities, how CHD8 haploinsufficiency affects neurodevelopmental is unclear. Here, employing human cerebral organoids, we find that CHD8 haploinsufficiency disrupted neurodevelopmental trajectories with an accelerated and delayed generation of, respectively, inhibitory and excitatory neurons that yields, at days 60 and 120, symmetrically opposite expansions in their proportions. This imbalance is consistent with an enlargement of cerebral organoids as an in vitro correlate of patients' macrocephaly. Through an isogenic design of patient-specific mutations and mosaic organoids, we define genotype-phenotype relationships and uncover their cell-autonomous nature. Our results define cell-type-specific CHD8-dependent molecular defects related to an abnormal program of proliferation and alternative splicing. By identifying cell-type-specific effects of CHD8 mutations, our study uncovers reproducible developmental alterations that may be employed for neurodevelopmental disease modeling.

Feature selection and classification over the network with missing node observations

Jointly analyzing transcriptomic data and the existing biological networks can yield more robust and informative feature selection results, as well as better understanding of the biological mechanisms. Selecting and classifying node features over genome-scale networks has become increasingly important in genomic biology and genomic medicine. Existing methods have some critical drawbacks. The first is they do not allow flexible modeling of different subtypes of selected nodes. The second is they ignore nodes with missing values, very likely to increase bias in estimation. To address these limitations, we propose a general modeling framework for Bayesian node classification (BNC) with missing values. A new prior model is developed for the class indicators incorporating the network structure. For posterior computation, we resort to the Swendsen-Wang algorithm for efficiently updating class indicators. BNC can naturally handle missing values in the Bayesian modeling framework, which improves the node classification accuracy and reduces the bias in estimating gene effects. We demonstrate the advantages of our methods via extensive simulation studies and the analysis of the cutaneous melanoma dataset from The Cancer Genome Atlas.

Altered intercellular communication and extracellular matrix signaling as a potential disease mechanism in human hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is considered a primary disorder of the sarcomere resulting in unexplained left ventricular hypertrophy but the paradoxical association of nonmyocyte phenotypes such as fibrosis, mitral valve anomalies and microvascular occlusion is unexplained. To understand the interplay between cardiomyocyte and nonmyocyte cell types in human HCM, single nuclei RNA-sequencing was performed on myectomy specimens from HCM patients with left ventricular outflow tract obstruction and control samples from donor hearts free of cardiovascular disease. Clustering analysis based on gene expression patterns identified a total of 34 distinct cell populations, which were classified into 10 different cell types based on marker gene expression. Differential gene expression analysis comparing HCM to Normal datasets revealed differences in sarcomere and extracellular matrix gene expression. Analysis of expressed ligand-receptor pairs across multiple cell types indicated profound alteration in HCM intercellular communication, particularly between cardiomyocytes and fibroblasts, fibroblasts and lymphocytes and involving integrin β1 and its multiple extracellular matrix (ECM) cognate ligands. These findings provide a paradigm for how sarcomere dysfunction is associated with reduced cardiomyocyte secretion of ECM ligands, altered fibroblast ligand-receptor interactions with other cell types and increased fibroblast to lymphocyte signaling, which can further alter the ECM composition and promote nonmyocyte phenotypes.

Development of a Model System to Study Expression Profile of RAC2 Gene in Breast Cancer MDA-MB-231 Cell Line

The RAC2 gene encoding GTPases involve cellular signaling of actin polymerization, cell migration, and formation of the phagocytic NADPH oxidase complex. Oncogenic mutations in the RAC2 gene have been identified in various cancers, and extensive research is in progress to delineate its signaling pathways and identify potential therapeutic targets in breast cancers. This paper explored developing a bioinformatics model system to understand the RAC2 gene expression pattern concerning estrogenic receptor status in breast cancers. We have used the MDA-MB-231 breast cancer cell line to identify RAC2 gene expression. To simplify the development of model system with one dataset, we retrieved the microarray dataset GSE27515 from the Gene Expression Omnibus (GEO) for the differential gene expression analysis. Then, network analysis, pathway enrichment analysis, volcano plot, ORA, and the up/downregulated genes were used to highlight genes involved in signaling network pathways. We observed that the RAC2 gene is upregulated in the GSM679722, GSM676923, and GSM679724 downregulated in the samples GSM676925, GSM676926, and GSM676927 from the GEO dataset. Our observation found that the RAC2 gene is upregulated in the estrogen receptor (ER) negative breast cancers and downregulated in ER-positive breast cancer, involving pathways such as focal adhesion, MAPK signaling, axon guidance, and VEGF signaling pathway.

Identification and characterisation of sPEPs in Cryptococcus neoformans

Short open reading frame (sORF)-encoded peptides (sPEPs) have been found across a wide range of genomic locations in a variety of species. To date, their identification, validation, and characterisation in the human fungal pathogen Cryptococcus neoformans has been limited due to a lack of standardised protocols. We have developed an enrichment process that enables sPEP detection within a protein sample from this polysaccharide-encapsulated yeast, and implemented proteogenomics to provide insights into the validity of predicted and hypothetical sORFs annotated in the C. neoformans genome. Novel sORFs were discovered within the 5' and 3' UTRs of known transcripts as well as in "non-coding" RNAs. One novel candidate, dubbed NPB1, that resided in an RNA annotated as "non-coding", was chosen for characterisation. Through the creation of both specific point mutations and a full deletion allele, the function of the new sPEP, Npb1, was shown to resemble that of the bacterial trans-translation protein SmpB.

miR-378 affects metabolic disturbances in the mdx model of Duchenne muscular dystrophy

Although Duchenne muscular dystrophy (DMD) primarily affects muscle tissues, the alterations to systemic metabolism manifested in DMD patients contribute to the severe phenotype of this fatal disorder. We propose that microRNA-378a (miR-378) alters carbohydrate and lipid metabolism in dystrophic mdx mice. In our study, we utilized double knockout animals which lacked both dystrophin and miR-378 (mdx/miR-378^-/-). RNA sequencing of the liver identified 561 and 194 differentially expressed genes that distinguished mdx versus wild-type (WT) and mdx/miR-378^-/- versus mdx counterparts, respectively. Bioinformatics analysis predicted, among others, carbohydrate metabolism disorder in dystrophic mice, as functionally proven by impaired glucose tolerance and insulin sensitivity. The lack of miR-378 in mdx animals mitigated those effects with a faster glucose clearance in a glucose tolerance test (GTT) and normalization of liver glycogen levels. The absence of miR-378 also restored the expression of genes regulating lipid homeostasis, such as Acly, Fasn, Gpam, Pnpla3, and Scd1. In conclusion, we report for the first time that miR-378 loss results in increased systemic metabolism of mdx mice. Together with our previous finding, demonstrating alleviation of the muscle-related symptoms of DMD, we propose that the inhibition of miR-378 may represent a new strategy to attenuate the multifaceted symptoms of DMD.

Plasma metabolomic profile changes in females with phenylketonuria following a camp intervention

BACKGROUND

There remains a limited understanding of the metabolic perturbations, beyond phenylalanine (Phe) metabolism, that contribute to phenotypic variability in phenylketonuria (PKU).

OBJECTIVES

This study aimed to characterize changes in the PKU plasma metabolome following a 5-d metabolic camp intervention and to compare PKU profiles with those of matched healthy controls.

METHODS

In 28 females (aged 12-57 y), fasting plasma samples were collected on the first (day 1) and final (day 5) days of camp to measure metabolic control and to complete untargeted metabolomic profiling. Three-day dietary records were collected to assess changes in dietary adherence and composition. Univariate (Wilcoxon signed-rank and Mann-Whitney U test) and multivariate (random forest, hierarchical clustering) analyses were performed to identify clinical and metabolic features that were associated with the intervention and disease state.

RESULTS

Relative to healthy controls, Phe catabolites, ketones, and carnitine- and glycine-conjugated fatty acids were elevated in females with PKU at baseline, whereas fatty acylcholine metabolites were substantially lower. After the camp intervention, plasma Phe concentrations decreased [median change: -173 µmol/L (IQR: -325, -28 µmol/L)] and 70% of PKU participants demonstrated improved dietary adherence by decreasing Phe intake and/or increasing medical food consumption. This was accompanied by a shift in abundance for 223 metabolites (q < 0.05). Compounds associated with the metabolism of Phe, fatty acids, and choline contributed most to profile differences between camp days 1 and 5.

CONCLUSIONS

In females with PKU, untargeted metabolomics identified prominent perturbations in amino acid and lipid metabolites associated with bioenergetic impairment and oxidative stress. Choline-conjugated lipids could have fundamental roles in these pathways and they have not been previously evaluated in PKU. A short-term camp intervention was effective for improving or fully normalizing the abundance of the identified discriminatory metabolites.

Row and Column Structure-Based Biclustering for Gene Expression Data

Due to the development of high-throughput technologies for gene analysis, the biclustering method has attracted much attention. However, existing methods have problems with high time and space complexity. This paper proposes a biclustering method, called Row and Column Structure-based Biclustering (RCSBC), with low time and space complexity to find checkerboard patterns within microarray data. First, the paper describes the structure of bicluster by using the structure of rows and columns. Second, the paper chooses the representative rows and columns with two algorithms. Finally, the gene expression data are biclustered on the space spanned by representative rows and columns. To the best of our knowledge, this paper is the first to exploit the relationship between the row/column structure of a gene expression matrix and the structure of biclusters. Both the synthetic datasets and the real-life gene expression datasets are used to validate the effectiveness of our method. It can be seen from the experiment results that the RCSBC outperforms the state-of-the-art algorithms both on clustering accuracy and time/space complexity. This study offers new insights into biclustering the large-scale gene expression data without loading the whole data into memory.

Network analysis of TCGA and GTEx gene expression datasets for identification of trait-associated biomarkers in human cancer

Advances in high-throughput sequencing technologies now yield unprecedented volumes of OMICs data with opportunities to conduct systematic data analyses and derive novel biological insights. Here, we provide protocols to perform differential-expressed gene analysis of TCGA and GTEx RNA-Seq data from human cancers, complete integrative GO and network analyses with focus on clinical and survival data, and identify differential correlation of trait-associated biomarkers.

For complete details on the use and execution of this protocol, please refer to Chen and MacDonald (2021).

•

Protocols for the identification of trait-associated molecular correlates in cancer

•

Differentially-expressed gene (DEG) analysis of TCGA and GTEx transcriptomic data

•

Protocols for integrative network analysis of RNA-seq, clinical, and survival data

•

Differential correlation of trait-associated biomarkers for hypothesis testing

Combined use of Oxford Nanopore and Illumina sequencing yields insights into soybean structural variation biology

BACKGROUND

Structural variants (SVs), including deletions, insertions, duplications, and inversions, are relatively long genomic variations implicated in a diverse range of processes from human disease to ecology and evolution. Given their complex signatures, tendency to occur in repeated regions, and large size, discovering SVs based on short reads is challenging compared to single-nucleotide variants. The increasing availability of long-read technologies has greatly facilitated SV discovery; however, these technologies remain too costly to apply routinely to population-level studies. Here, we combined short-read and long-read sequencing technologies to provide a comprehensive population-scale assessment of structural variation in a panel of Canadian soybean cultivars.

RESULTS

We used Oxford Nanopore long-read sequencing data (~12× mean coverage) for 17 samples to both benchmark SV calls made from Illumina short-read data and predict SVs that were subsequently genotyped in a population of 102 samples using Illumina data. Benchmarking results show that variants discovered using Oxford Nanopore can be accurately genotyped from the Illumina data. We first use the genotyped deletions and insertions for population genetics analyses and show that results are comparable to those based on single-nucleotide variants. We observe that the population frequency and distribution within the genome of deletions and insertions are constrained by the location of genes. Gene Ontology and PFAM domain enrichment analyses also confirm previous reports that genes harboring high-frequency deletions and insertions are enriched for functions in defense response. Finally, we discover polymorphic transposable elements from the deletions and insertions and report evidence of the recent activity of a Stowaway MITE.

CONCLUSIONS

We show that structural variants discovered using Oxford Nanopore data can be genotyped with high accuracy from Illumina data. Our results demonstrate that long-read and short-read sequencing technologies can be efficiently combined to enhance SV analysis in large populations, providing a reusable framework for their study in a wider range of samples and non-model species.

Physiological and transcriptome analysis reveals the differences in nitrate content between lamina and midrib of flue-cured tobacco

Nitrate is an important precursor of tobacco-specific nitrosamines (TSNAs) and a remarkable difference in nitrate accumulation between lamina and midrib of flue-cured tobacco has long been observed. However, the physiological and molecular mechanisms underpinning this difference remain poorly understood. In this study, physiological and genetic factors impacting nitrate accumulation were identified in pot experiments using flue-cured tobacco K326 with contrasting nitrate content between lamina and midrib. The results showed that three times higher NO3-N content was observed in midrib than that in the lamina, along with lower pigment, NH4-N content, nitrate reductase activity (NRA), sucrose synthetase activity (SSA), and glutamine synthetase activity (GSA) in midrib. Transcriptome analysis revealed that expression of genes involved in porphyrin and chlorophyll metabolism, carotenoid biosynthesis, photosynthesis-antenna proteins, photosynthesis, carbon fixation in photosynthetic organisms, starch and sucrose metabolism, nitrogen metabolism, and biosynthesis of amino acids were significantly lower in midrib than in lamina. qRT-PCR results showed that the expression level of nitrate transporter genes LOC107782967, LOC107806749, LOC107775674, LOC107829632, LOC107799198, LOC107768465 decreased by 2.74, 1.81, 49.5, 3.5, 2.64 and 2.96-folds while LOC107789301 increased by 8.23-folds in midrib but not in lamina. Reduced chlorophyll content might result in low carbohydrate formation which is the source of energy and carbon skeleton supply, then the low capacity of nitrogen reduction, assimilation and transportation, and the poor ability of nitrate reallocation but the high capacity of accumulation might lead to nitrate accumulation in midrib. The results laid the foundation for reducing nitrate content and TSNA formation in tobacco midribs and their products.

The Nephrotoxin Puromycin Aminonucleoside Induces Injury in Kidney Organoids Differentiated from Induced Pluripotent Stem Cells

Kidney diseases, including acute kidney injury (AKI) and chronic kidney disease (CKD), which can progress to end stage renal disease (ESRD), are a worldwide health burden. Organ transplantation or kidney dialysis are the only effective available therapeutic tools. Therefore, in vitro models of kidney diseases and the development of prospective therapeutic options are urgently needed. Within the kidney, the glomeruli are involved in blood filtration and waste excretion and are easily affected by changing cellular conditions. Puromycin aminonucleoside (PAN) is a nephrotoxin, which can be employed to induce acute glomerular damage and to model glomerular disease. For this reason, we generated kidney organoids from three iPSC lines and treated these with PAN in order to induce kidney injury. Morphological observations revealed the disruption of glomerular and tubular structures within the kidney organoids upon PAN treatment, which were confirmed by transcriptome analyses. Subsequent analyses revealed an upregulation of immune response as well as inflammatory and cell-death-related processes. We conclude that the treatment of iPSC-derived kidney organoids with PAN induces kidney injury mediated by an intertwined network of inflammation, cytoskeletal re-arrangement, DNA damage, apoptosis and cell death. Furthermore, urine-stem-cell-derived kidney organoids can be used to model kidney-associated diseases and drug discovery.

Changes in Skeletal Muscle PAK1 Levels Regulate Tissue Crosstalk to Impact Whole Body Glucose Homeostasis

Skeletal muscle accounts for ~80% of insulin-stimulated glucose uptake. The Group I p21-activated kinase 1 (PAK1) is required for the non-canonical insulin-stimulated GLUT4 vesicle translocation in skeletal muscle cells. We found that the abundances of PAK1 protein and its downstream effector in muscle, ARPC1B, are significantly reduced in the skeletal muscle of humans with type 2 diabetes, compared to the non-diabetic controls, making skeletal muscle PAK1 a candidate regulator of glucose homeostasis. Although whole-body PAK1 knockout mice exhibit glucose intolerance and are insulin resistant, the contribution of skeletal muscle PAK1 in particular was unknown. As such, we developed inducible skeletal muscle-specific PAK1 knockout (skmPAK1-iKO) and overexpression (skmPAK1-iOE) mouse models to evaluate the role of PAK1 in skeletal muscle insulin sensitivity and glucose homeostasis. Using intraperitoneal glucose tolerance and insulin tolerance testing, we found that skeletal muscle PAK1 is required for maintaining whole body glucose homeostasis. Moreover, PAK1 enrichment in GLUT4-myc-L6 myoblasts preserves normal insulin-stimulated GLUT4 translocation under insulin resistance conditions. Unexpectedly, skmPAK1-iKO also showed aberrant plasma insulin levels following a glucose challenge. By applying conditioned media from PAK1-enriched myotubes or myoblasts to β-cells in culture, we established that a muscle-derived circulating factor(s) could enhance β-cell function. Taken together, these data suggest that PAK1 levels in the skeletal muscle can regulate not only skeletal muscle insulin sensitivity, but can also engage in tissue crosstalk with pancreatic β-cells, unveiling a new molecular mechanism by which PAK1 regulates whole-body glucose homeostasis.

Genomic signatures and insights into host niche adaptation of the entomopathogenic fungus Metarhizium humberi

The genus Metarhizium is composed of species used in biological control programs of agricultural pests worldwide. This genus includes common fungal pathogen of many insects and mites and endophytes that can increase plant growth. Metarhizium humberi was recently described as a new species. This species is highly virulent against some insect pests and promotes growth in sugarcane, strawberry, and soybean crops. In this study, we sequenced the genome of M. humberi, isolate ESALQ1638, and performed a functional analysis to determine its genomic signatures and highlight the genes and biological processes associated with its lifestyle. The genome annotation predicted 10633 genes in M. humberi, of which 92.0% are assigned putative functions, and ∼17% of the genome was annotated as repetitive sequences. We found that 18.5% of the M. humberi genome is similar to experimentally validated proteins associated with pathogen-host interaction. Compared to the genomes of eight Metarhizium species, the M. humberi ESALQ1638 genome revealed some unique traits that stood out, e.g., more genes functionally annotated as polyketide synthases (PKSs), overrepresended GO-terms associated to transport of ions, organic and amino acid, a higher percentage of repetitive elements, and higher levels of RIP-induced point mutations. The M. humberi genome will serve as a resource for promoting studies on genome structure and evolution that can contribute to research on biological control and plant biostimulation. Thus, the genomic data supported the broad host range of this species within the generalist PARB clade and suggested that M. humberi ESALQ1638 might be particularly good at producing secondary metabolites and might be more efficient in transporting amino acids and organic compounds.

Single-Cell Ribonucleic Acid Sequencing Clarifies Cold Tolerance Mechanisms in the Pacific White Shrimp (Litopenaeus Vannamei)

To characterize the cold tolerance mechanism of the Pacific white shrimp (Litopenaeus vannamei), we performed single-cell RNA sequencing (scRNA-seq) of ∼5185 hepatopancreas cells from cold-tolerant (Lv-T) and common (Lv-C) L. vannamei at preferred and low temperatures (28°C and 10°C, respectively). The cells fell into 10 clusters and 4 cell types: embryonic, resorptive, blister-like, and fibrillar. We identified differentially expressed genes between Lv-T and Lv-C, which were mainly associated with the terms “immune system,” “cytoskeleton,” “antioxidant system,” “digestive enzyme,” and “detoxification,” as well as the pathways “metabolic pathways of oxidative phosphorylation,” “metabolism of xenobiotics by cytochrome P450,” “chemical carcinogenesis,” “drug metabolism-cytochrome P450,” and “fatty acid metabolism.” Reconstruction of fibrillar cell trajectories showed that, under low temperature stress, hepatopancreas cells had two distinct fates, cell fate 1 and cell fate 2. Cell fate 1 was mainly involved in signal transduction and sensory organ development. Cell fate 2 was mainly involved in metabolic processes. This study preliminarily clarifies the molecular mechanisms underlying cold tolerance in L. vannamei, which will be useful for the breeding of shrimp with greater cold tolerance.

Zika virus disrupts gene expression in human myoblasts and myotubes: Relationship with susceptibility to infection

The tropism of Zika virus (ZIKV) has been described in the nervous system, blood, placenta, thymus, and skeletal muscle. We investigated the mechanisms of skeletal muscle susceptibility to ZIKV using an in vitro model of human skeletal muscle myogenesis, in which myoblasts differentiate into myotubes. Myoblasts were permissive to ZIKV infection, generating productive viral particles, while myotubes controlled ZIKV replication. To investigate the underlying mechanisms, we used gene expression profiling. First, we assessed gene changes in myotubes compared with myoblasts in the model without infection. As expected, we observed an increase in genes and pathways related to the contractile muscle system in the myotubes, a reduction in processes linked to proliferation, migration and cytokine production, among others, confirming the myogenic capacity of our system in vitro. A comparison between non-infected and infected myoblasts revealed more than 500 differentially expressed genes (DEGs). In contrast, infected myotubes showed almost 2,000 DEGs, among which we detected genes and pathways highly or exclusively expressed in myotubes, including those related to antiviral and innate immune responses. Such gene modulation could explain our findings showing that ZIKV also invades myotubes but does not replicate in these differentiated cells. In conclusion, we showed that ZIKV largely (but differentially) disrupts gene expression in human myoblasts and myotubes. Identifying genes involved in myotube resistance can shed light on potential antiviral mechanisms against ZIKV infection.

Spatiotemporal analysis identifies ABF2 and ABF3 as key hubs of endodermal response to nitrate

Nitrate is a nutrient and a potent signal that impacts global gene expression in plants. However, the regulatory factors controlling temporal and cell type-specific nitrate responses remain largely unknown. We assayed nitrate-responsive transcriptome changes in five major root cell types of the Arabidopsis thaliana root as a function of time. We found that gene-expression response to nitrate is dynamic and highly localized and predicted cell type-specific transcription factor (TF)-target interactions. Among cell types, the endodermis stands out as having the largest and most connected nitrate-regulatory gene network. ABF2 and ABF3 are major hubs for transcriptional responses in the endodermis cell layer. We experimentally validated TF-target interactions for ABF2 and ABF3 by chromatin immunoprecipitation followed by sequencing and a cell-based system to detect TF regulation genome-wide. Validated targets of ABF2 and ABF3 account for more than 50% of the nitrate-responsive transcriptome in the endodermis. Moreover, ABF2 and ABF3 are involved in nitrate-induced lateral root growth. Our approach offers an unprecedented spatiotemporal resolution of the root response to nitrate and identifies important components of cell-specific gene regulatory networks.

A novel gene functional similarity calculation model by utilizing the specificity of terms and relationships in gene ontology

Background

Recently, with the foundation and development of gene ontology (GO) resources, numerous works have been proposed to compute functional similarity of genes and achieved series of successes in some research fields. Focusing on the calculation of the information content (IC) of terms is the main idea of these methods, which is essential for measuring functional similarity of genes. However, most approaches have some deficiencies, especially when measuring the IC of both GO terms and their corresponding annotated term sets. To this end, measuring functional similarity of genes accurately is still challenging.

Results

In this article, we proposed a novel gene functional similarity calculation method, which especially encapsulates the specificity of terms and edges (STE). The proposed method mainly contains three steps. Firstly, a novel computing model is put forward to compute the IC of terms. This model has the ability to exploit the specific structural information of GO terms. Secondly, the IC of term sets are computed by capturing the genetic structure between the terms contained in the set. Lastly, we measure the gene functional similarity according to the IC overlap ratio of the corresponding annotated genes sets. The proposed method accurately measures the IC of not only GO terms but also the annotated term sets by leveraging the specificity of edges in the GO graph.

Conclusions

We conduct experiments on gene functional classification in biological pathways, gene expression datasets, and protein-protein interaction datasets. Extensive experimental results show the better performances of our proposed STE against several baseline methods.

KRAB zinc finger protein ZNF676 controls the transcriptional influence of LTR12-related endogenous retrovirus sequences

BACKGROUND

Transposable element-embedded regulatory sequences (TEeRS) and their KRAB-containing zinc finger protein (KZFP) controllers are increasingly recognized as modulators of gene expression. We aim to characterize the contribution of this system to gene regulation in early human development and germ cells.

RESULTS

Here, after studying genes driven by the long terminal repeat (LTR) of endogenous retroviruses, we identify the ape-restricted ZNF676 as the sequence-specific repressor of a subset of contemporary LTR12 integrants responsible for a large fraction of transpochimeric gene transcripts (TcGTs) generated during human early embryogenesis. We go on to reveal that the binding of this KZFP correlates with the epigenetic marking of these TEeRS in the germline, and is crucial to the control of genes involved in ciliogenesis/flagellogenesis, a biological process that dates back to the last common ancestor of eukaryotes.

CONCLUSION

These results illustrate how KZFPs and their TE targets contribute to the evolutionary turnover of transcription networks and participate in the transgenerational inheritance of epigenetic traits.

Genome-wide perturbations of Alu expression and Alu-associated post-transcriptional regulations distinguish oligodendroglioma from other gliomas

Alu is a primate-specific repeat element in the human genome and has been increasingly appreciated as a regulatory element in many biological processes. But the appreciation of Alu has been limited in tumorigenesis, especially for brain tumor. To investigate the relevance of Alu to the gliomagenesis, we studied Alu element-associated post-transcriptional processes and the RNA expression of the element by performing RNA-seq for a total of 41 pairs of neurotypical and diverse glioma brain tissues. We find that A-to-I editing and circular RNA levels, as well as Alu RNA expression, are decreased overall in gliomas, compared to normal tissue. Interestingly, grade 2 oligodendrogliomas are least affected in A-to-I editing and circular RNA levels among gliomas, whereas they have a higher proportion of down-regulated Alu subfamilies, compared to the other gliomas. These findings collectively imply a unique pattern of Alu-associated transcriptomes in grade 2 oligodendroglioma, providing an insight to gliomagenesis from the perspective of an evolutionary genetic element.

Common and Distinctive Intercellular Communication Patterns in Human Obstructive and Nonobstructive Hypertrophic Cardiomyopathy

Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder characterized by unexplained left ventricular hypertrophy with or without left ventricular outflow tract (LVOT) obstruction. Single-nuclei RNA-sequencing (snRNA-seq) of both obstructive and nonobstructive HCM patient samples has revealed alterations in communication between various cell types, but no direct and integrated comparison between the two HCM phenotypes has been reported. We performed a bioinformatic analysis of HCM snRNA-seq datasets from obstructive and nonobstructive patient samples to identify differentially expressed genes and distinctive patterns of intercellular communication. Differential gene expression analysis revealed 37 differentially expressed genes, predominantly in cardiomyocytes but also in other cell types, relevant to aging, muscle contraction, cell motility, and the extracellular matrix. Intercellular communication was generally reduced in HCM, affecting the extracellular matrix, growth factor binding, integrin binding, PDGF binding, and SMAD binding, but with increases in adenylate cyclase binding, calcium channel inhibitor activity, and serine-threonine kinase activity in nonobstructive HCM. Increases in neuron to leukocyte and dendritic cell communication, in fibroblast to leukocyte and dendritic cell communication, and in endothelial cell communication to other cell types, largely through changes in the expression of integrin-β1 and its cognate ligands, were also noted. These findings indicate both common and distinct physiological mechanisms affecting the pathogenesis of obstructive and nonobstructive HCM and provide opportunities for the personalized management of different HCM phenotypes.

Quantitative Proteome and Transcriptome Dynamics Analysis Reveals Iron Deficiency Response Networks and Signature in Neuronal Cells

Iron and oxygen deficiencies are common features in pathophysiological conditions, such as ischemia, neurological diseases, and cancer. Cellular adaptive responses to such deficiencies include repression of mitochondrial respiration, promotion of angiogenesis, and cell cycle control. We applied a systematic proteomics analysis to determine the global proteomic changes caused by acute hypoxia and chronic and acute iron deficiency (ID) in hippocampal neuronal cells. Our analysis identified over 8600 proteins, revealing similar and differential effects of each treatment on activation and inhibition of pathways regulating neuronal development. In addition, comparative analysis of ID-induced proteomics changes in cultured cells and transcriptomic changes in the rat hippocampus identified common altered pathways, indicating specific neuronal effects. Transcription factor enrichment and correlation analysis identified key transcription factors that were activated in both cultured cells and tissue by iron deficiency, including those implicated in iron regulation, such as HIF1, NFY, and NRF1. We further identified MEF2 as a novel transcription factor whose activity was induced by ID in both HT22 proteome and rat hippocampal transcriptome, thus linking iron deficiency to MEF2-dependent cellular signaling pathways in neuronal development. Taken together, our study results identified diverse signaling networks that were differentially regulated by hypoxia and ID in neuronal cells.

Quantitative metaproteomics and activity-based protein profiling of patient fecal microbiome identifies host and microbial serine-type endopeptidase activity associated with ulcerative colitis

The gut microbiota plays an important yet incompletely understood role in the induction and propagation of ulcerative colitis (UC). Organism-level efforts to identify UC-associated microbes have revealed the importance of community structure, but less is known about the molecular effectors of disease. We performed 16S rRNA gene sequencing in parallel with label-free data-dependent LC-MS/MS proteomics to characterize the stool microbiomes of healthy (n = 8) and UC (n = 10) patients. Comparisons of taxonomic composition between techniques revealed major differences in community structure partially attributable to the additional detection of host, fungal, viral, and food peptides by metaproteomics. Differential expression analysis of metaproteomic data identified 176 significantly enriched protein groups between healthy and UC patients. Gene ontology analysis revealed several enriched functions with serine-type endopeptidase activity overrepresented in UC patients. Using a biotinylated fluorophosphonate probe and streptavidin-based enrichment, we show that serine endopeptidases are active in patient fecal samples and that additional putative serine hydrolases are detectable by this approach compared with unenriched profiling. Finally, as metaproteomic databases expand, they are expected to asymptotically approach completeness. Using ComPIL and de novo peptide sequencing, we estimate the size of the probable peptide space unidentified (“dark peptidome”) by our large database approach to establish a rough benchmark for database sufficiency. Despite high variability inherent in patient samples, our analysis yielded a catalog of differentially enriched proteins between healthy and UC fecal proteomes. This catalog provides a clinically relevant jumping-off point for further molecular-level studies aimed at identifying the microbial underpinnings of UC.

•

Identified 176 significantly altered protein groups between healthy and UC patients.

•

Serine-type endopeptidase activity is overrepresented in UC patients.

•

Fluorophosphonate ABPP shows that endopeptidases are active in fecal samples.

•

ABPP enrichment helps identify additional putative serine hydrolases in samples.

•

De novo sequencing used to estimate number of MS2 spectra unidentified by ComPIL.

Genomic loci mispositioning in Tmem120a knockout mice yields latent lipodystrophy

Little is known about how the observed fat-specific pattern of 3D-spatial genome organisation is established. Here we report that adipocyte-specific knockout of the gene encoding nuclear envelope transmembrane protein Tmem120a disrupts fat genome organisation, thus causing a lipodystrophy syndrome. Tmem120a deficiency broadly suppresses lipid metabolism pathway gene expression and induces myogenic gene expression by repositioning genes, enhancers and miRNA-encoding loci between the nuclear periphery and interior. Tmem120a-/- mice, particularly females, exhibit a lipodystrophy syndrome similar to human familial partial lipodystrophy FPLD2, with profound insulin resistance and metabolic defects that manifest upon exposure to an obesogenic diet. Interestingly, similar genome organisation defects occurred in cells from FPLD2 patients that harbour nuclear envelope protein encoding LMNA mutations. Our data indicate TMEM120A genome organisation functions affect many adipose functions and its loss may yield adiposity spectrum disorders, including a miRNA-based mechanism that could explain muscle hypertrophy in human lipodystrophy.

The amniotic fluid proteome changes with gestational age in normal pregnancy: a cross-sectional study

The cell-free transcriptome in amniotic fluid (AF) has been shown to be informative of physiologic and pathologic processes in pregnancy; however, the change in AF proteome with gestational age has mostly been studied by targeted approaches. The objective of this study was to describe the gestational age-dependent changes in the AF proteome during normal pregnancy by using an omics platform. The abundance of 1310 proteins was measured on a high-throughput aptamer-based proteomics platform in AF samples collected from women during midtrimester (16-24 weeks of gestation, n = 15) and at term without labor (37-42 weeks of gestation, n = 13). Only pregnancies without obstetrical complications were included in the study. Almost 25% (320) of AF proteins significantly changed in abundance between the midtrimester and term gestation. Of these, 154 (48.1%) proteins increased, and 166 (51.9%) decreased in abundance at term compared to midtrimester. Tissue-specific signatures of the trachea, salivary glands, brain regions, and immune system were increased while those of the gestational tissues (uterus, placenta, and ovary), cardiac myocytes, and fetal liver were decreased at term compared to midtrimester. The changes in AF protein abundance were correlated with those previously reported in the cell-free AF transcriptome. Intersecting gestational age-modulated AF proteins and their corresponding mRNAs previously reported in the maternal blood identified neutrophil-related protein/mRNA pairs that were modulated in the same direction. The first study to utilize an aptamer-based assay to profile the AF proteome modulation with gestational age, it reveals that almost one-quarter of the proteins are modulated as gestation advances, which is more than twice the fraction of altered plasma proteins (~ 10%). The results reported herein have implications for future studies focused on discovering biomarkers to predict, monitor, and diagnose obstetrical diseases.

Endothelial and systemic upregulation of miR-34a-5p fine-tunes senescence in progeria

Endothelial defects significantly contribute to cardiovascular pathology in the premature aging disease Hutchinson-Gilford progeria syndrome (HGPS). Using an endothelium-specific progeria mouse model, we identify a novel, endothelium-specific microRNA (miR) signature linked to the p53-senescence pathway and a senescence-associated secretory phenotype (SASP). Progerin-expressing endothelial cells exert profound cell-non-autonomous effects initiating senescence in non-endothelial cell populations and causing immune cell infiltrates around blood vessels. Comparative miR expression analyses revealed unique upregulation of senescence-associated miR34a-5p in endothelial cells with strong accumulation at atheroprone aortic arch regions but also, in whole cardiac- and lung tissues as well as in the circulation of progeria mice. Mechanistically, miR34a-5p knockdown reduced not only p53 levels but also late-stage senescence regulator p16 with no effect on p21 levels, while p53 knockdown reduced miR34a-5p and partially rescued p21-mediated cell cycle inhibition with a moderate effect on SASP. These data demonstrate that miR34a-5p reinforces two separate senescence regulating branches in progerin-expressing endothelial cells, the p53- and p16-associated pathways, which synergistically maintain a senescence phenotype that contributes to cardiovascular pathology. Thus, the key function of circulatory miR34a-5p in endothelial dysfunction-linked cardiovascular pathology offers novel routes for diagnosis, prognosis and treatment for cardiovascular aging in HGPS and potentially geriatric patients.

Exercise-Induced Extracellular Vesicles Delay the Progression of Prostate Cancer

Increasing evidence suggests that regular physical exercise not only reduces the risk of cancer but also improves functional capacity, treatment efficacy and disease outcome in cancer patients. At least partially, these effects are mediated by the secretome of the tissues responding to exercise. The secreted molecules can be released in a carrier-free form or enclosed into extracellular vesicles (EVs). Several recent studies have shown that EVs are actively released into circulation during physical exercise. Here, we for the first time investigated the effects of exercise-induced EVs on the progression of cancer in an F344 rat model of metastatic prostate cancer. Although we did not observe a consistent increase in the circulating EV levels, RNA sequencing analysis demonstrated substantial changes in the RNA content of EVs collected before and immediately after forced wheel running exercise as well as differences between EVs from runners at resting state and sedentary rats. The major RNA biotype in EVs was mRNA, followed by miRNA and rRNA. Molecular functions of differentially expressed RNAs reflected various physiological processes including protein folding, metabolism and regulation of immune responses triggered by the exercise in the parental cells. Intravenous administration of exercise-induced EVs into F344 rats with orthotopically injected syngeneic prostate cancer cells PLS10, demonstrated reduction of the primary tumor volume by 35% and possibly—attenuation of lung metastases. Hence, our data provide the first evidence that exercise-induced EVs may modulate tumor physiology and delay the progression of cancer.

Grainyhead 1 acts as a drug-inducible conserved transcriptional regulator linked to insulin signaling and lifespan

Aging is impacted by interventions across species, often converging on metabolic pathways. Transcription factors regulate longevity yet approaches for their pharmacological modulation to exert geroprotection remain sparse. We show that increased expression of the transcription factor Grainyhead 1 (GRH-1) promotes lifespan and pathogen resistance in Caenorhabditis elegans. A compound screen identifies FDA-approved drugs able to activate human GRHL1 and promote nematodal GRH-1-dependent longevity. GRHL1 activity is regulated by post-translational lysine methylation and the phosphoinositide (PI) 3-kinase C2A. Consistently, nematodal longevity following impairment of the PI 3-kinase or insulin/IGF-1 receptor requires grh-1. In BXD mice, Grhl1 expression is positively correlated with lifespan and insulin sensitivity. In humans, GRHL1 expression positively correlates with insulin receptor signaling and also with lifespan. Fasting blood glucose levels, including in individuals with type 2 diabetes, are negatively correlated with GRHL1 expression. Thereby, GRH-1/GRHL1 is identified as a pharmacologically malleable transcription factor impacting insulin signaling and lifespan.

Life- and healthspan of organisms can be modulated by dietary, genetic, or pharmacological interventions, which often affect metabolic pathways. Here the authors report that Grainyhead 1 is an evolutionarily conserved, drug-inducible transcription factor that promotes longevity in C. elegans, and thus a potential target for the development of geroprotective drugs.

Cadmium interference with iron sensing reveals transcriptional programs sensitive and insensitive to reactive oxygen species

Iron (Fe) is an essential micronutrient whose uptake is tightly regulated to prevent either deficiency or toxicity. Cadmium (Cd) is a non-essential element that induces both Fe deficiency and toxicity; however, the mechanisms behind these Fe/Cd-induced responses are still elusive. Here we explored Cd- and Fe-associated responses in wild-type Arabidopsis and in a mutant that overaccumulates Fe (opt3-2). Gene expression profiling revealed a large overlap between transcripts induced by Fe deficiency and Cd exposure. Interestingly, the use of opt3-2 allowed us to identify additional gene clusters originally induced by Cd in the wild type but repressed in the opt3-2 background. Based on the high levels of H2O2 found in opt3-2, we propose a model where reactive oxygen species prevent the induction of genes that are induced in the wild type by either Fe deficiency or Cd. Interestingly, a defined cluster of Fe-responsive genes was found to be insensitive to this negative feedback, suggesting that their induction by Cd is more likely to be the result of an impaired Fe sensing. Overall, our data suggest that Fe deficiency responses are governed by multiple inputs and that a hierarchical regulation of Fe homeostasis prevents the induction of specific networks when Fe and H2O2 levels are elevated.

Nearby transposable elements impact plant stress gene regulatory networks: a meta-analysis in A. thaliana and S. lycopersicum

BACKGROUND

Transposable elements (TE) make up a large portion of many plant genomes and are playing innovative roles in genome evolution. Several TEs can contribute to gene regulation by influencing expression of nearby genes as stress-responsive regulatory motifs. To delineate TE-mediated plant stress regulatory networks, we took a 2-step computational approach consisting of identifying TEs in the proximity of stress-responsive genes, followed by searching for cis-regulatory motifs in these TE sequences and linking them to known regulatory factors. Through a systematic meta-analysis of RNA-seq expression profiles and genome annotations, we investigated the relation between the presence of TE superfamilies upstream, downstream or within introns of nearby genes and the differential expression of these genes in various stress conditions in the TE-poor Arabidopsis thaliana and the TE-rich Solanum lycopersicum.

RESULTS

We found that stress conditions frequently expressed genes having members of various TE superfamilies in their genomic proximity, such as SINE upon proteotoxic stress and Copia and Gypsy upon heat stress in A. thaliana, and EPRV and hAT upon infection, and Harbinger, LINE and Retrotransposon upon light stress in S. lycopersicum. These stress-specific gene-proximal TEs were mostly located within introns and more detected near upregulated than downregulated genes. Similar stress conditions were often related to the same TE superfamily. Additionally, we detected both novel and known motifs in the sequences of those TEs pointing to regulatory cooption of these TEs upon stress. Next, we constructed the regulatory network of TFs that act through binding these TEs to their target genes upon stress and discovered TE-mediated regulons targeted by TFs such as BRB/BPC, HD, HSF, GATA, NAC, DREB/CBF and MYB factors in Arabidopsis and AP2/ERF/B3, NAC, NF-Y, MYB, CXC and HD factors in tomato.

CONCLUSIONS

Overall, we map TE-mediated plant stress regulatory networks using numerous stress expression profile studies for two contrasting plant species to study the regulatory role TEs play in the response to stress. As TE-mediated gene regulation allows plants to adapt more rapidly to new environmental conditions, this study contributes to the future development of climate-resilient plants.

OUP accepted manuscript

Motivation

There is a plethora of measures to evaluate functional similarity (FS) of genes based on their co-expression, protein–protein interactions and sequence similarity. These measures are typically derived from hand-engineered and application-specific metrics to quantify the degree of shared information between two genes using their Gene Ontology (GO) annotations.

Results

We introduce deepSimDEF, a deep learning method to automatically learn FS estimation of gene pairs given a set of genes and their GO annotations. deepSimDEF’s key novelty is its ability to learn low-dimensional embedding vector representations of GO terms and gene products and then calculate FS using these learned vectors. We show that deepSimDEF can predict the FS of new genes using their annotations: it outperformed all other FS measures by >5–10% on yeast and human reference datasets on protein–protein interactions, gene co-expression and sequence homology tasks. Thus, deepSimDEF offers a powerful and adaptable deep neural architecture that can benefit a wide range of problems in genomics and proteomics, and its architecture is flexible enough to support its extension to any organism.

Availability and implementation

Source code and data are available at https://github.com/ahmadpgh/deepSimDEF

Supplementary information

Supplementary data are available at Bioinformatics online.

Loss of Planar Cell Polarity Effector Fuzzy Causes Renal Hypoplasia by Disrupting Several Signaling Pathways

In vertebrates, the planar cell polarity (PCP) pathway regulates tissue morphogenesis during organogenesis, including the kidney. Mutations in human PCP effector proteins have been associated with severe syndromic ciliopathies. Importantly, renal hypoplasia has been reported in some patients. However, the developmental disturbance that causes renal hypoplasia is unknown. Here, we describe the early onset of profound renal hypoplasia in mice homozygous for null mutation of the PCP effector gene, Fuzzy. We found that this phenotype is caused by defective branching morphogenesis of the ureteric bud (UB) in the absence of defects in nephron progenitor specification or in early steps of nephrogenesis. By using various experimental approaches, we show that the loss of Fuzzy affects multiple signaling pathways. Specifically, we found mild involvement of GDNF/c-Ret pathway that drives UB branching. We noted the deficient expression of molecules belonging to the Bmp, Fgf and Shh pathways. Analysis of the primary cilia in the UB structures revealed a significant decrease in ciliary length. We conclude that renal hypoplasia in the mouse Fuzzy mutants is caused by defective UB branching associated with dysregulation of ciliary and non-ciliary signaling pathways. Our work suggests a PCP effector-dependent pathogenetic mechanism that contributes to renal hypoplasia in mice and humans.

Honey Bee Larval and Adult Microbiome Life Stages Are Effectively Decoupled with Vertical Transmission Overcoming Early Life Perturbations

Microbiomes provide a range of benefits to their hosts which can lead to the coevolution of a joint ecological niche. However, holometabolous insects, some of the most successful organisms on Earth, occupy different niches throughout development, with larvae and adults being physiologically and morphologically highly distinct. Furthermore, transition between the stages usually involves the loss of the gut microbiome since the gut is remodeled during pupation. Most eusocial organisms appear to have evolved a workaround to this problem by sharing their communal microbiome across generations. However, whether this vertical microbiome transmission can overcome perturbations of the larval microbiome remains untested. Honey bees have a relatively simple, conserved, coevolved adult microbiome which is socially transmitted and affects many aspects of their biology. In contrast, larval microbiomes are more variable, with less clear roles. Here, we manipulated the gut microbiome of in vitro-reared larvae, and after pupation of the larvae, we inoculated the emerged bees with adult microbiome to test whether adult and larval microbiome stages may be coupled (e.g., through immune priming). Larval treatments differed in bacterial composition and abundance, depending on diet, which also drove larval gene expression. Nonetheless, adults converged on the typical core taxa and showed limited gene expression variation. This work demonstrates that honey bee adult and larval stages are effectively microbiologically decoupled, and the core adult microbiome is remarkably stable to early developmental perturbations. Combined with the transmission of the microbiome in early adulthood, this allows the formation of long-term host-microbiome associations. IMPORTANCE This work investigated host-microbiome interactions during a crucial developmental stage-the transition from larvae to adults, which is a challenge to both, the insect host and its microbiome. Using the honey bee as a tractable model system, we showed that microbiome transfer after emergence overrides any variation in the larvae, indicating that larval and adult microbiome stages are effectively decoupled. Together with the reliable vertical transfer in the eusocial system, this decoupling ensures that the adults are colonized with a consistent and derived microbiome after eclosion. Taken all together, our data provide additional support that the evolution of sociality, at least in the honey bee system tested here, is linked with host-microbiome relationships.

A Sensitization-Free Dimethyl Fumarate Prodrug, Isosorbide Di-(Methyl Fumarate), Provides a Topical Treatment Candidate for Psoriasis

Dimethyl fumarate (DMF) is an effective oral treatment for psoriasis administered in Europe for nearly 60 years. However, its potential has been limited by contact dermatitis that prohibits topical application. This paper characterizes a DMF derivative, isosorbide DMF (IDMF), which was designed to have antipsoriatic effects without skin-sensitizing properties. We show that IDMF exhibits neither genotoxicity nor radiation sensitivity in skin fibroblasts and is nonirritating and nonsensitizing in animal models (rat, rabbit, guinea pig). Microarray analysis of cytokine-stimulated keratinocytes showed that IDMF represses the expression of genes specifically upregulated in psoriatic skin lesions but not those of other skin diseases. IDMF also downregulated genes induced by IL-17A and TNF in keratinocytes as well as predicted targets of NF-κB and the antidifferentiation noncoding RNA (i.e., ANCR). IDMF further stimulated the transcription of oxidative stress response genes (NQO1, GPX2, GSR) with stronger NRF2/ARE activation compared to DMF. Finally, IDMF reduced erythema and scaling while repressing the expression of immune response genes in psoriasiform lesions elicited by topical application of imiquimod in mice. These data show that IDMF exhibits antipsoriatic activity that is similar or improved compared with that exhibited by DMF, without the harsh skin-sensitizing effects that have prevented topical delivery of the parent molecule.

Exploiting a Y chromosome-linked Cas9 for sex selection and gene drive

CRISPR-based genetic engineering tools aimed to bias sex ratios, or drive effector genes into animal populations, often integrate the transgenes into autosomal chromosomes. However, in species with heterogametic sex chromsomes (e.g. XY, ZW), sex linkage of endonucleases could be beneficial to drive the expression in a sex-specific manner to produce genetic sexing systems, sex ratio distorters, or even sex-specific gene drives, for example. To explore this possibility, here we develop a transgenic line of Drosophila melanogaster expressing Cas9 from the Y chromosome. We functionally characterize the utility of this strain for both sex selection and gene drive finding it to be quite effective. To explore its utility for population control, we built mathematical models illustrating its dynamics as compared to other state-of-the-art systems designed for both population modification and suppression. Taken together, our results contribute to the development of current CRISPR genetic control tools and demonstrate the utility of using sex-linked Cas9 strains for genetic control of animals.

The neutrophil protein CD177 is a novel PDPN receptor that regulates human cancer-associated fibroblast physiology

The cancer-associated fibroblast (CAF) marker podoplanin (PDPN) is generally correlated with poor clinical outcomes in cancer patients and thus represents a promising therapeutic target. Despite its biomedical relevance, basic aspects of PDPN biology such as its cellular functions and cell surface ligands remain poorly uncharacterized, thus challenging drug development. Here, we utilize a high throughput platform to elucidate the PDPN cell surface interactome, and uncover the neutrophil protein CD177 as a new binding partner. Quantitative proteomics analysis of the CAF phosphoproteome reveals a role for PDPN in cell signaling, growth and actomyosin contractility, among other processes. Moreover, cellular assays demonstrate that CD177 is a functional antagonist, recapitulating the phenotype observed in PDPN-deficient CAFs. In sum, starting from the unbiased elucidation of the PDPN co-receptome, our work provides insights into PDPN functions and reveals the PDPN/CD177 axis as a possible modulator of fibroblast physiology in the tumor microenvironment.

YAP/TAZ and ATF4 drive resistance to Sorafenib in hepatocellular carcinoma by preventing ferroptosis

Understanding the mechanisms underlying evasive resistance in cancer is an unmet medical need to improve the efficacy of current therapies. In this study, a combination of shRNA-mediated synthetic lethality screening and transcriptomic analysis revealed the transcription factors YAP/TAZ as key drivers of Sorafenib resistance in hepatocellular carcinoma (HCC) by repressing Sorafenib-induced ferroptosis. Mechanistically, in a TEAD-dependent manner, YAP/TAZ induce the expression of SLC7A11, a key transporter maintaining intracellular glutathione homeostasis, thus enabling HCC cells to overcome Sorafenib-induced ferroptosis. At the same time, YAP/TAZ sustain the protein stability, nuclear localization, and transcriptional activity of ATF4 which in turn cooperates to induce SLC7A11 expression. Our study uncovers a critical role of YAP/TAZ in the repression of ferroptosis and thus in the establishment of Sorafenib resistance in HCC, highlighting YAP/TAZ-based rewiring strategies as potential approaches to overcome HCC therapy resistance.

Error-prone, stress-induced 3' flap-based Okazaki fragment maturation supports cell survival

How cells with DNA replication defects acquire mutations that allow them to escape apoptosis under environmental stress is a long-standing question. Here, we report that an error-prone Okazaki fragment maturation (OFM) pathway is activated at restrictive temperatures in rad27Δ yeast cells. Restrictive temperature stress activated Dun1, facilitating transformation of unprocessed 5′ flaps into 3′ flaps, which were removed by 3′ nucleases, including DNA polymerase δ (Polδ). However, at certain regions, 3′ flaps formed secondary structures that facilitated 3′ end extension rather than degradation, producing alternative duplications with short spacer sequences, such as pol3 internal tandem duplications. Consequently, little 5′ flap was formed, suppressing rad27Δ-induced lethality at restrictive temperatures. We define a stress-induced, error-prone OFM pathway that generates mutations that counteract replication defects and drive cellular evolution and survival.

Tumor evolution selectively inactivates the core microRNA machinery for immune evasion

Cancer cells acquire genetic heterogeneity to escape from immune surveillance during tumor evolution, but a systematic approach to distinguish driver from passenger mutations is lacking. Here we investigate the impact of different immune pressure on tumor clonal dynamics and immune evasion mechanism, by combining massive parallel sequencing of immune edited tumors and CRISPR library screens in syngeneic mouse tumor model and co-culture system. We find that the core microRNA (miRNA) biogenesis and targeting machinery maintains the sensitivity of cancer cells to PD-1-independent T cell-mediated cytotoxicity. Genetic inactivation of the machinery or re-introduction of ANKRD52 frequent patient mutations dampens the JAK-STAT-interferon-γ signaling and antigen presentation in cancer cells, largely by abolishing miR-155-targeted silencing of suppressor of cytokine signaling 1 (SOCS1). Expression of each miRNA machinery component strongly correlates with intratumoral T cell infiltration in nearly all human cancer types. Our data indicate that the evolutionarily conserved miRNA pathway can be exploited by cancer cells to escape from T cell-mediated elimination and immunotherapy.

MicroRNA-100-5p and microRNA-298-5p released from apoptotic cortical neurons are endogenous Toll-like receptor 7/8 ligands that contribute to neurodegeneration

Background

MicroRNA (miRNA) expression in the brain is altered in neurodegenerative diseases. Recent studies demonstrated that selected miRNAs conventionally regulating gene expression at the post-transcriptional level can act extracellularly as signaling molecules. The identity of miRNA species serving as membrane receptor ligands involved in neuronal apoptosis in the central nervous system (CNS), as well as the miRNAs’ sequence and structure required for this mode of action remained largely unresolved.

Methods

Using a microarray-based screening approach we analyzed apoptotic cortical neurons of C56BL/6 mice and their supernatant with respect to alterations in miRNA expression/presence. HEK-Blue Toll-like receptor (TLR) 7/8 reporter cells, primary microglia and macrophages derived from human and mouse were employed to test the potential of the identified miRNAs released from apoptotic neurons to serve as signaling molecules for the RNA-sensing receptors. Biophysical and bioinformatical approaches, as well as immunoassays and sequential microscopy were used to analyze the interaction between candidate miRNA and TLR. Immunocytochemical and -histochemical analyses of murine CNS cultures and adult mice intrathecally injected with miRNAs, respectively, were performed to evaluate the impact of miRNA-induced TLR activation on neuronal survival and microglial activation.

Results

We identified a specific pattern of miRNAs released from apoptotic cortical neurons that activate TLR7 and/or TLR8, depending on sequence and species. Exposure of microglia and macrophages to certain miRNA classes released from apoptotic neurons resulted in the sequence-specific production of distinct cytokines/chemokines and increased phagocytic activity. Out of those miRNAs miR-100-5p and miR-298-5p, which have consistently been linked to neurodegenerative diseases, entered microglia, located to their endosomes, and directly bound to human TLR8. The miRNA-TLR interaction required novel sequence features, but no specific structure formation of mature miRNA. As a consequence of miR-100-5p- and miR-298-5p-induced TLR activation, cortical neurons underwent cell-autonomous apoptosis. Presence of miR-100-5p and miR-298-5p in cerebrospinal fluid led to neurodegeneration and microglial accumulation in the murine cerebral cortex through TLR7 signaling.

Conclusion

Our data demonstrate that specific miRNAs are released from apoptotic cortical neurons, serve as endogenous TLR7/8 ligands, and thereby trigger further neuronal apoptosis in the CNS. Our findings underline the recently discovered role of miRNAs as extracellular signaling molecules, particularly in the context of neurodegeneration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13024-021-00498-5.

A Novel Calibration Step in Gene Co-Expression Network Construction

High-throughput technologies such as DNA microarrays and RNA-sequencing are used to measure the expression levels of large numbers of genes simultaneously. To support the extraction of biological knowledge, individual gene expression levels are transformed to Gene Co-expression Networks (GCNs). In a GCN, nodes correspond to genes, and the weight of the connection between two nodes is a measure of similarity in the expression behavior of the two genes. In general, GCN construction and analysis includes three steps; 1) calculating a similarity value for each pair of genes 2) using these similarity values to construct a fully connected weighted network 3) finding clusters of genes in the network, commonly called modules. The specific implementation of these three steps can significantly impact the final output and the downstream biological analysis. GCN construction is a well-studied topic. Existing algorithms rely on relatively simple statistical and mathematical tools to implement these steps. Currently, software package WGCNA appears to be the most widely accepted standard. We hypothesize that the raw features provided by sequencing data can be leveraged to extract modules of higher quality. A novel preprocessing step of the gene expression data set is introduced that in effect calibrates the expression levels of individual genes, before computing pairwise similarities. Further, the similarity is computed as an inner-product of positive vectors. In experiments, this provides a significant improvement over WGCNA, as measured by aggregate p-values of the gene ontology term enrichment of the computed modules.

Bioinformatic Analysis and Integration of Transcriptome and Proteome Results Identify Key Coding and Noncoding Genes Predicting Malignancy in Intraductal Papillary Mucinous Neoplasms of the Pancreas

BACKGROUND

Intraductal papillary mucinous neoplasms (IPMNs) are precursor lesions of pancreatic ductal adenocarcinoma (PDAC). IPMNs are generally associated with high risk of developing malignancy and therefore need to be diagnosed and assessed accurately, once detected. Existing diagnostic methods are inadequate, and identification of efficient biomarker capable of detecting high-risk IPMNs is necessitated. Moreover, the mechanism of development of malignancy in IPMNs is also elusive.

METHODS

Gene expression meta-analysis conducted using 12 low-risk IPMN and 23 high-risk IPMN tissue samples. We have also listed all the altered miRNAs and long noncoding RNAs (lncRNAs), identified their target genes, and performed pathway analysis. We further enlisted cyst fluid proteins detected to be altered in high-risk or malignant IPMNs and compared them with fraction of differentially expressed genes secreted into cyst fluid.

RESULTS

Our meta-analysis identified 270 upregulated and 161 downregulated genes characteristically altered in high-risk IPMNs. We further identified 61 miRNAs and 14 lncRNAs and their target genes and key pathways contributing towards understanding of the gene regulation during the progression of the disease. Most importantly, we have detected 12 genes altered significantly both in cystic lesions and cyst fluid.

CONCLUSION

Our study reports, for the first time, a meta-analysis identifying key changes in gene expression between low-risk and high-risk IPMNs and also explains the regulatory aspect through construction of a miRNA-lncRNA-mRNA interaction network. The 12-gene-signature could function as potential biomarker in cyst fluid for detection of IPMN with a high risk of developing malignancy.

Gene Expression Signature Associated with Clinical Outcome in ALK-Positive Anaplastic Large Cell Lymphoma

Simple Summary

Anaplastic large cell lymphomas associated with ALK translocation have a good outcome after CHOP treatment; however, the 2-year relapse rate remains at 30%. Microarray gene-expression profiling, high throughput RT-qPCR, and RNA sequencing of 48 ALK-positive anaplastic large cell lymphoma (ALK⁺ ALCL) samples obtained at diagnosis enable the identification of genes associated with clinical outcome. More particularly, our molecular signatures indicate that the FN1 gene, a matrix key regulator, might also be involved in the prognosis and the therapeutic response in anaplastic lymphomas.

Abstract

Anaplastic large cell lymphomas associated with ALK translocation have a good outcome after CHOP treatment; however, the 2-year relapse rate remains at 30%. Microarray gene-expression profiling of 48 samples obtained at diagnosis was used to identify 47 genes that were differentially expressed between patients with early relapse/progression and no relapse. In the relapsing group, the most significant overrepresented genes were related to the regulation of the immune response and T-cell activation while those in the non-relapsing group were involved in the extracellular matrix. Fluidigm technology gave concordant results for 29 genes, of which FN1, FAM179A, and SLC40A1 had the strongest predictive power after logistic regression and two classification algorithms. In parallel with 39 samples, we used a Kallisto/Sleuth pipeline to analyze RNA sequencing data and identified 20 genes common to the 28 genes validated by Fluidigm technology—notably, the FAM179A and FN1 genes. Interestingly, FN1 also belongs to the gene signature predicting longer survival in diffuse large B-cell lymphomas treated with CHOP. Thus, our molecular signatures indicate that the FN1 gene, a matrix key regulator, might also be involved in the prognosis and the therapeutic response in anaplastic lymphomas.

The transcription factor PagLBD3 contributes to the regulation of secondary growth in Populus

LATERAL ORGAN BOUNDARIES DOMAIN (LBD) genes encode plant-specific transcription factors that participate in regulating various developmental processes. In this study, we genetically characterized PagLBD3 encoding an important regulator of secondary growth in poplar (Populus alba × Populus glandulosa). Overexpression of PagLBD3 increased stem secondary growth in Populus with a significantly higher rate of cambial cell differentiation into phloem, while dominant repression of PagLBD3 significantly decreased the rate of cambial cell differentiation into phloem. Furthermore, we identified 1756 PagLBD3 genome-wide putative direct target genes (DTGs) through RNA sequencing (RNA-seq)-coupled DNA affinity purification followed by sequencing (DAP-seq) assays. Gene Ontology analysis revealed that genes regulated by PagLBD3 were enriched in biological pathways regulating meristem development, xylem development, and auxin transport. Several central regulator genes for vascular development, including PHLOEM INTERCALATED WITH XYLEM (PXY), WUSCHEL RELATED HOMEOBOX4 (WOX4), Secondary Wall-Associated NAC Domain 1s (SND1-B2), and Vascular-Related NAC-Domain 6s (VND6-B1), were identified as PagLBD3 DTGs. Together, our results indicate that PagLBD3 and its DTGs form a complex transcriptional network to modulate cambium activity and phloem/xylem differentiation.

Transcriptome changes in maternal peripheral blood during term parturition mimic perturbations preceding spontaneous Preterm birth†

The complex physiologic process of parturition includes the onset of labor, which requires the orchestrated stimulation of a common pathway involving uterine contractility, cervical ripening, and chorioamniotic membrane activation. However, the labor-specific processes taking place in these tissues have limited use as predictive biomarkers unless they can be probed in non-invasive samples, such as the peripheral blood. Herein, we utilized a transcriptomic dataset to assess labor-specific changes in the peripheral blood of women who delivered at term. We identified a set of genes that were differentially expressed with labor and enriched for immunological processes, and these gene expression changes were strongly correlated with results from prior studies, providing in silico validation of our findings. We then identified significant correlations between labor-specific transcriptomic changes in the maternal circulation and those reported in the chorioamniotic membranes, myometrium, and cervix of women at term, demonstrating that tissue-specific labor signatures are partly mirrored in the peripheral blood. Last, we demonstrated a significant overlap between the peripheral blood transcriptomic changes in term parturition and those observed in asymptomatic women prior to the diagnosis of preterm prelabor rupture of membranes who delivered preterm. Collectively, we provide evidence that the normal process of labor at term is characterized by a unique immunological expression signature, which may serve as a useful tool for assessing labor status and potentially identifying women at risk for preterm birth.

Genome and transcriptome analysis of the beet armyworm Spodoptera exigua reveals targets for pest control

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies.

Human Induced Pluripotent Stem Cell-Derived Mesenchymal Stem Cells Acquire Rejuvenation and Reduced Heterogeneity

Despite the uniform selection criteria for the isolation of human mesenchymal stem cells (MSCs), considerable heterogeneity exists which reflects the distinct tissue origins and differences between individuals with respect to their genetic background and age. This heterogeneity is manifested by the variabilities seen in the transcriptomes, proteomes, secretomes, and epigenomes of tissue-specific MSCs. Here, we review literature on different aspects of MSC heterogeneity including the role of epigenetics and the impact of MSC heterogeneity on therapies. We then combine this with a meta-analysis of transcriptome data from distinct MSC subpopulations derived from bone marrow, adipose tissue, cruciate, tonsil, kidney, umbilical cord, fetus, and induced pluripotent stem cells derived MSCs (iMSCs). Beyond that, we investigate transcriptome differences between tissue-specific MSCs and pluripotent stem cells. Our meta-analysis of numerous MSC-related data sets revealed markers and associated biological processes characterizing the heterogeneity and the common features of MSCs from various tissues. We found that this heterogeneity is mainly related to the origin of the MSCs and infer that microenvironment and epigenetics are key drivers. The epigenomes of MSCs alter with age and this has a profound impact on their differentiation capabilities. Epigenetic modifications of MSCs are propagated during cell divisions and manifest in differentiated cells, thus contributing to diseased or healthy phenotypes of the respective tissue. An approach used to reduce heterogeneity caused by age- and tissue-related epigenetic and microenvironmental patterns is the iMSC concept: iMSCs are MSCs generated from induced pluripotent stem cells (iPSCs). During iMSC generation epigenetic and chromatin remodeling result in a gene expression pattern associated with rejuvenation thus allowing to overcome age-related shortcomings (e.g., limited differentiation and proliferation capacity). The importance of the iMSC concept is underlined by multiple clinical trials. In conclusion, we propose the use of rejuvenated iMSCs to bypass tissue- and age-related heterogeneity which are associated with native MSCs.

An integrated network analysis approach to identify potential key genes, transcription factors, and microRNAs regulating human hematopoietic stem cell aging

Hematopoietic stem cells (HSCs) undergo functional deterioration with increasing age that causes loss of their self-renewal and regenerative potential. Despite various efforts, significant success in identifying molecular regulators of HSC aging has not been achieved, one prime reason being the non-availability of appropriate human HSC samples. To demonstrate the scope of integrating and re-analyzing the HSC transcriptomics data available, we used existing tools and databases to structure a sequential data analysis pipeline to predict potential candidate genes, transcription factors, and microRNAs simultaneously. This sequential approach comprises (i) collecting matched young and aged mice HSC sample datasets, (ii) identifying differentially expressed genes, (iii) identifying human homologs of differentially expressed genes, (iv) inferring gene co-expression network modules, and (v) inferring the microRNA-transcription factor-gene regulatory network. Systems-level analyses of HSC interaction networks provided various insights based on which several candidates were predicted. For example, 16 HSC aging-related candidate genes were predicted (e.g., CD38, BRCA1, AGTR1, GSTM1, etc.) from GCN analysis. Following this, the shortest path distance-based analyses of the regulatory network predicted several novel candidate miRNAs and TFs. Among these, miR-124-3p was a common regulator in candidate gene modules, while TFs MYC and SP1 were identified to regulate various candidate genes. Based on the regulatory interactions among candidate genes, TFs, and miRNAs, a potential regulation model of biological processes in each of the candidate modules was predicted, which provided systems-level insights into the molecular complexity of each module to regulate HSC aging.

Substances for regenerative wound healing during antler renewal stimulated scar-less restoration of rat cutaneous wounds

Scarification is the outcome of cutaneous wound healing under normal conditions. Although considerable effort has been expended in this field, scar-less healing has not been achieved satisfactorily. The lack of a good model of scar-free healing has contributed to this undesirable situation. However, the annual regeneration of deer antlers, which starts from regenerative wound healing over the top of the pedicles (permanent bony protuberances), may provide such a model. Therefore, in this study, we investigated the process of pedicle wound healing at the organ, tissue, cell, and molecular levels. Our results convincingly demonstrate that wounds over the pedicle preceded a regenerative healing process including regeneration of skin appendages, such as hair follicles. Compared to the scar healing in rats, regenerative healing of the pedicle wound exhibited a weaker inflammatory response, lack of myofibroblast induction, and higher ratios of Col III/Col I, TGF-β3/TGF-β1, and MMP/TIMP. Importantly, our periosteal transplantation experiments in vivo revealed that this regenerative healing process was achieved through induction of antler stem cells (ASCs). Further study showed that this effect of ASCs on regenerative healing was not species-specific but more generic and could be applied to other mammalian species, as injection of ASCs stimulated regenerative healing of full-thickness excisional cutaneous wounds in rats. Overall, our findings show that ASCs may have therapeutic potential in enhancing the quality of wound healing and preventing scar formation in clinical settings.

Deciphering two rounds of cell lineage segregations during bovine preimplantation development

Blastocyst formation gives rise to the inner cell mass (ICM) and trophectoderm (TE) and is followed by the differentiation of the epiblast (Epi) and primitive endoderm (PrE) within the ICM. Although these two-round cell lineage differentiations underpin proper embryogenesis in every mammal, their spatiotemporal dynamics are quite diverse among species. Here, molecular details of the blastocyst stage in cattle were dissected using an optimized in vitro culture method. Blastocyst embryos were placed on agarose gel filled with nutrient-rich media to expose embryos to both gaseous and liquid phases. Embryos derived from this "on-gel" culture were transferred to surrogate mothers on day (D) 10 after fertilization and successfully implanted. Immunofluorescent studies using on-gel-cultured embryos revealed that the proportion of TE cells expressing the pluripotent ICM marker, OCT4, which was beyond 80% on D8, was rapidly reduced after D9 and reached 0% on D9.5. This first lineage segregation process was temporally parallel with the second one, identified by the spatial separation of Epi cells expressing SOX2 and PrE cells expressing SOX17. RNA-seq comparison of TE cells from D8 in vitro fertilized embryos and D14 in vivo embryos revealed that besides drastic reduction of pluripotency-related genes, TE cells highly expressed Wnt, FGF, and VEGF signaling pathways-related genes to facilitate the functional maturation required for feto-maternal interaction. Quantitative PCR analysis of TE cells derived from on-gel culture further confirmed time-dependent increments in the expression of key TE markers. Altogether, the present study provides platforms to understand species-specific strategies for mammalian preimplantation development.