Similar Transcriptomes Research Articles

Homozygous iMycCα transgenic mice as a model of plasma B-cell lymphomas

The 3' regulatory region (3'RR) located downstream from the Cα gene is the conductor of transcription, accessibility, and remodeling of the IgH locus at mature B-cell stages. Convincing demonstrations of the essential contributions of the 3'RR in B-cell lymphomagenesis have been provided by mouse models which bring the oncogene c-Myc under the 3'RR transcriptional control. In this study, we developed a mouse model of CD138+ plasma B-cell lymphomas. If the KI of c-myc directly into Cα just 5' to the 3'RR in iMycCα mice produced B-cell lymphomas with low kinetics, we enforced c-myc production in iMycCα mice by the generation of homozygous c-myc transgenic mice. Our results show that homozygous iMycCα mice lead to a mouse model of plasma CD138+ B-cell lymphomas with interesting and wide transcriptomic similarities to human multiple myeloma and appropriated emergence kinetics that can be used to test new experimental therapeutic approaches.

Medullary Thyroid Carcinoma Mutational Spectrum Update and Signaling-Type Inference by Transcriptional Profiles: Literature Meta-Analysis and Study of Tumor Samples

Simple SummaryMedullary thyroid carcinoma (MTC) is a rare but clinically relevant tumor based on its aggressiveness and the limited therapeutic opportunities currently available for advanced cases. A better understanding of the mechanisms of MTC development is crucial to identify more effective means of intervention and therapies. Several studies have shown that RET and RAS genes play a central role in MTC. However, little is known about the signaling processes operating downstream of these genes. Here, we report mutation and gene expression profiles in proprietary sporadic MTCs, including both primary and metastatic tumors. We show that tumors derived from the same patient display similar expression profiles and that the latter can be used to obtain information about specific downstream signaling, identifying distinct molecular subtypes. Furthermore, by reviewing the relevant literature, we highlight that, along with RET and RAS, other less frequent genes are emerging as possible new players in MTC.Medullary thyroid carcinoma (MTC) is a rare but aggressive tumor. Although RET and RAS genes are recognized drivers in MTC, associated downstream signaling pathways are largely unknown. In this study, we report 17 sporadic MTCs, collected at our institution, comprising patient-matched primary and lymph node metastatic tumors investigated for mutational and transcriptional profiles. As we identified two uncommon RET deletions (D898_E901del and E632_L633del), we also performed a literature review and meta-analysis to assess the occurrence of unconventional alterations in MTC, focusing on next-generation sequencing studies. We found that new gene alterations are emerging, along with the known RET/RAS drivers, involving not only RET by multiple concurrent mutations or deletions but also other previously underestimated cancer-related genes, especially in sporadic MTCs. In our MTC gene profiles, we found transcriptome similarity between patient-matched tissues and expression of immune genes only by a few samples. Furthermore, we defined a gene signature able to stratify samples into two distinct signaling types, termed MEN2B-like and MEN2A-like. We provide an updated overview of the MTC mutational spectrum and describe how transcriptional profiles can be used to define distinct MTC signaling subtypes that appear to be shared by various gene drivers, including the unconventional ones.

Elevated Kir2.1/nuclear N2ICD defines a highly malignant subtype of non-WNT/SHH medulloblastomas

Medulloblastoma (MB) is one of the most common childhood malignant brain tumors (WHO grade IV), traditionally divided into WNT, SHH, Group 3, and Group 4 subgroups based on the transcription profiles, somatic DNA alterations, and clinical outcomes. Unlike WNT and SHH subgroup MBs, Group 3 and Group 4 MBs have similar transcriptomes and lack clearly specific drivers and targeted therapeutic options. The recently revised WHO Classification of CNS Tumors has assigned Group 3 and 4 to a provisional non-WNT/SHH entity. In the present study, we demonstrate that Kir2.1, an inwardly-rectifying potassium channel, is highly expressed in non-WNT/SHH MBs, which promotes tumor cell invasion and metastasis by recruiting Adam10 to enhance S2 cleavage of Notch2 thereby activating the Notch2 signaling pathway. Disruption of the Notch2 pathway markedly inhibited the growth and metastasis of Kir2.1-overexpressing MB cell-derived xenograft tumors in mice. Moreover, Kir2.1high/nuclear N2ICDhigh MBs are associated with the significantly shorter lifespan of the patients. Thus, Kir2.1high/nuclear N2ICDhigh can be used as a biomarker to define a novel subtype of non-WNT/SHH MBs. Our findings are important for the modification of treatment regimens and the development of novel-targeted therapies for non-WNT/SHH MBs.

Uncovering the molecular identity of cardiosphere-derived cells (CDCs) by single-cell RNA sequencing

Cardiosphere-derived cells (CDCs) generated from human cardiac biopsies have been shown to have disease-modifying bioactivity in clinical trials. Paradoxically, CDCs’ cellular origin in the heart remains elusive. We studied the molecular identity of CDCs using single-cell RNA sequencing (sc-RNAseq) in comparison to cardiac non-myocyte and non-hematopoietic cells (cardiac fibroblasts/CFs, smooth muscle cells/SMCs and endothelial cells/ECs). We identified CDCs as a distinct and mitochondria-rich cell type that shared biological similarities with non-myocyte cells but not with cardiac progenitor cells derived from human-induced pluripotent stem cells. CXCL6 emerged as a new specific marker for CDCs. By analysis of sc-RNAseq data from human right atrial biopsies in comparison with CDCs we uncovered transcriptomic similarities between CDCs and CFs. By direct comparison of infant and adult CDC sc-RNAseq data, infant CDCs revealed GO-terms associated with cardiac development. To analyze the beneficial effects of CDCs (pro-angiogenic, anti-fibrotic, anti-apoptotic), we performed functional in vitro assays with CDC-derived extracellular vesicles (EVs). CDC EVs augmented in vitro angiogenesis and did not stimulate scarring. They also reduced the expression of pro-apoptotic Bax in NRCMs. In conclusion, CDCs were disclosed as mitochondria-rich cells with unique properties but also with similarities to right atrial CFs. CDCs displayed highly proliferative, secretory and immunomodulatory properties, characteristics that can also be found in activated or inflammatory cell types. By special culture conditions, CDCs earn some bioactivities, including angiogenic potential, which might modify disease in certain disorders.

Basal-Like Cell-Conditioned Medium Exerts Anti-Fibrotic Effects In Vitro and In Vivo.

In idiopathic pulmonary fibrosis (IPF), basal-like cells are atypically present in the alveolar region, where they may affect adjacent stromal cells by paracrine mechanisms. We here aimed to confirm the presence of basal-like cells in peripheral IPF lung tissue in vivo, to culture and characterize the cells in vitro, and to investigate their paracrine effects on IPF fibroblasts in vitro and in bleomycin-injured rats in vivo. Basal-like cells are mainly localized in areas of pathological bronchiolization or honeycomb cysts in peripheral IPF lung tissue. Single-cell RNA sequencing (scRNA-seq) demonstrated an overall homogeneity, the expression of the basal cell markers cytokeratin KRT5 and KRT17, and close transcriptomic similarities to basal cells in the majority of cells cultured in vitro. Basal-like cells secreted significant levels of prostaglandin E2 (PGE2), and their conditioned medium (CM) inhibited alpha-smooth muscle actin (α-SMA) and collagen 1A1 (Col1A1) and upregulated matrix metalloproteinase-1 (MMP-1) and hepatocyte growth factor (HGF) by IPF fibroblasts in vitro. The instillation of CM in bleomycin-injured rat lungs resulted in reduced collagen content, improved lung architecture, and reduced α-SMA-positive cells. Our data suggested that basal-like cells may limit aberrant fibroblast activation and differentiation in IPF through paracrine mechanisms.

The Application of Consensus Weighted Gene Co-expression Network Analysis to Comparative Transcriptome Meta-Datasets of Multiple Sclerosis in Gray and White Matter.

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system characterized by demyelination, which leads to the formation of white matter lesions (WMLs) and gray matter lesions (GMLs). Recently, a large amount of transcriptomics or proteomics research works explored MS, but few studies focused on the differences and similarities between GMLs and WMLs in transcriptomics. Furthermore, there are astonishing pathological differences between WMLs and GMLs, for example, there are differences in the type and abundance of infiltrating immune cells between WMLs and GMLs. Here, we used consensus weighted gene co-expression network analysis (WGCNA), single-sample gene set enrichment analysis (ssGSEA), and machine learning methods to identify the transcriptomic differences and similarities of the MS between GMLs and WMLs, and to find the co-expression modules with significant differences or similarities between them. Through weighted co-expression network analysis and ssGSEA analysis, CD56 bright natural killer cell was identified as the key immune infiltration factor in MS, whether in GM or WM. We also found that the co-expression networks between the two groups are quite similar (density = 0.79), and 28 differentially expressed genes (DEGs) are distributed in the midnightblue module, which is most related to CD56 bright natural killer cell in GM. Simultaneously, we also found that there are huge disparities between the modules, such as divergences between darkred module and lightyellow module, and these divergences may be relevant to the functions of the genes in the modules.

Single-Cell Transcriptome Analysis Reveals Embryonic Endothelial Heterogeneity at Spatiotemporal Level and Multifunctions of MicroRNA-126 in Mice.

Endothelial cells (ECs) play a critical role in angiogenesis and vascular remodeling. The heterogeneity of ECs has been reported at adult stages, yet it has not been fully investigated. This study aims to assess the transcriptional heterogeneity of developmental ECs at spatiotemporal level and to reveal the changes of embryonic ECs clustering when endothelium-enriched microRNA-126 (miR-126) was specifically knocked out. C57BL/6J mice embryos at day 14.5 were harvested and digested, followed by fluorescence-activated cell sorting to enrich ECs. Then, single-cell RNA sequencing was applied to enriched embryonic ECs. Tie2 (Tek receptor tyrosine kinase)-cre-mediated ECs-specific miR-126 knockout mice were constructed, and ECs from Tie2-cre-mediated ECs-specific miR-126 knockout embryos were subjected to single-cell RNA sequencing. Embryonic ECs were clustered into 11 groups corresponding to anatomic characteristics. The vascular bed (arteries, capillaries, veins, lymphatics) exhibited transcriptomic similarity across the developmental stage. Embryonic ECs had higher proliferative potential than adult ECs. Integrating analysis showed that 3 ECs populations (hepatic, mesenchymal transition, and pulmonary ECs) were apparently disorganized after miR-126 being knocked out. Gene ontology analysis revealed that disrupted ECs were mainly related to hypoxia, glycometabolism, and vascular calcification. Additionally, in vivo experiment showed that Tie2-cre-mediated ECs-specific miR-126 knockout mice exhibited excessive intussusceptive angiogenesis; reductive glucose and pyruvate tolerance; and excessive accumulation of calcium. Agonist miR-126-3p agomir significantly rescued the phenotype of glucose metabolic dysfunction in Tie2-cre-mediated ECs-specific miR-126 knockout mice. The heterogeneity of ECs is established as early as the embryonic stage. The deficiency of miR-126 disrupts the differentiation and diversification of embryonic ECs, suggesting that miR-126 plays an essential role in the maintenance of ECs heterogeneity.

Human Trophectoderm Spheroid Derived from Human Embryonic Stem Cells.

The use of human embryos for studying the early implantation processes and trophoblast is restricted by ethical concerns. The development of models mimicking the peri-implantation embryos is critical for understanding the physiology of human embryos and many pathophysiological disorders including recurrent implantation failure and miscarriage. Three-dimensional (3D) models of trophoblastic spheroids have been successfully derived from human embryonic stem cells (hESC). Simultaneous activation of the BMP pathway and blockage of the Activin/Nodal pathway favor the differentiation of hESC into trophoblast. Here we describe a 3D trophectoderm differentiation protocol with the use of BAP (BMP4, A83-01, and PD173074) to generate hESC-derived trophectoderm spheroids (BAP-EB). BAP-EB is highly reproducible and exhibits morphological and transcriptomic similarities to human early blastocysts.

Transcriptome analysis reveals the contribution of oligodendrocyte and radial glia-derived cues for maintenance of microglia identity.

Microglia are increasingly being recognized as druggable targets in neurodegenerative disorders, and good in vitro models are crucial to address cell biological questions. Major challenges are to recapitulate the complex microglial morphology and their in vivo transcriptome. We have therefore exposed primary microglia from adult rhesus macaques to a variety of different culture conditions including exposure to soluble factors as M-CSF, IL-34, and TGF-β as well as serum replacement approaches, and compared their morphologies and transcriptomes to those of mature, homeostatic in vivo microglia. This enabled us to develop a new, partially serum-free, monoculture protocol, that yields high numbers of ramified cells. We also demonstrate that exposure of adult microglia to M-CSF or IL-34 induces similar transcriptomes, and that exposure to TGF-β has much less pronounced effects than it does on rodent microglia. However, regardless of culture conditions, the transcriptomes of in vitro and in vivo microglia remained substantially different. Analysis of differentially expressed genes inspired us to perform 3D-spherical coculture experiments of microglia with oligodendrocytes and radial glia. In such spheres, microglia signature genes were strongly induced, even in the absence of neurons and astrocytes. These data reveal a novel role for oligodendrocyte and radial glia-derived cues in the maintenance of microglial identity, providing new anchor points to study microglia in health and disease.

Shared transcriptomic signature between Huntington's disease and early-onset Alzheimer's disease.

Neurodegenerative diseases share common pathological features, including progressive loss of neurons and protein misfolding, which ultimately culminates in dementia. Interestingly, dementia is not restricted to elderly people. In fact, many diseases have been identified as causes of early-onset dementia (< 65 years of age) such as Huntington's disease (HD) and early-onset Alzheimer's disease (EOAD). Importantly, disease-specific genetic mutations have already been identified for EOAD and HD. Thus, one could suggest that the molecular link between EOAD and HD may arise from alterations at the transcriptomic level, which is yet to be determined. In keeping with this, we hypothesized that EOAD and HD share a transcriptomic signature. In this work, we aimed at identifying transcriptome similarities between both diseases. We collected data sample of postmortem cerebral cortex from 6 AD (GSE33000, GSE118553, GSE122063, GSE39420, GSE44770, GSE44771) and 1 HD (GSE33000) microarray studies in the Gene Expression Omnibus, a public transcriptomic data repository. Of note, only subjects with age at death under 65 were selected (EOAD: n = 65, controls: n = 266; HD: n = 158, controls: n = 158). Differential expression and functional enrichment analyses were performed for each group versus their respective age-matched controls. We identified 1,260 differentially expressed genes (476 up- and 784 down-regulated genes) and 675 enriched gene ontology (GO) terms between EOAD and HD. The GO terms were grouped in 11 clusters of biological processes (Table 1). Cluster observation suggests a link in the apoptotic signaling, immune response and synaptic transmission processes in both diseases. Our results demonstrate a transcriptomic signature shared by EOAD and HD. Unveiling the similarities between EOAD and HD at the transcriptomic level could advance our knowledge about the pathogenesis of these diseases and may help developing therapeutic strategies targeting early-onset neurodegenerative diseases.

Hippocampal transcriptomic profiling of Alzheimer’s disease patients and mouse models: Differences and similarities

Rodent models of Alzheimer's disease (AD) have provide invaluable insights on AD pathophysiology in the past decades. Therefore, the evaluation of similarities between models and human pathology can improve AD research. Here, we compared hippocampal transcriptomic profiles of human late onset-AD (LOAD) and early-onset AD (EOAD) individuals with three mouse models (hAβ-KI, APP/PS1 and 5xFAD) to evaluate transcriptomic similarities between these animal models and human pathology. While APP/PS1 and 5xFAD are well established EOAD models, hAβ-KI was recently introduced as a potential LOAD model. We hypothesized that APP/PS1 and 5xFAD will have more similiatires to EOAD human AD and hAβ-KI with human LOAD. Five publicly available human AD/cognitively unimpaired transcriptomic profiles were collected from GEO (https://www.ncbi.nlm.nih.gov/geo/), merged and submitted to differential expression analysis (DEA) and functional enrichment analysis (FEA) using R. APP/PS1 mouse model RNAseq data from two datasets (GSE149661 and GSE145907) were also collected. Finally, RNAseq data of 5xFAD and hAβ-KI mouse models were obtained from AMP-AD Knowledge Portal (https://www.synapse.org/) using synapser and synapserutils packages. All animal models were submitted to DEA and FEA. DEA identified 1164, 3261 and 1782 differentially expressed genes (DEGs) for the comparisons between hAβ-KI, 5xFAD and APP/PS1 mutants versus wild type controls, respectively. The hAβ-KI mice model exhibited more DEGs in common with LOAD than with EOAD patients. Interestingly, APP/PS1 and 5xFAD mice also presented more DEGs in common with LOAD than with EOAD patients. However, 5xFAD showed significantly higher DEGs intersection with EOAD than with hAβ-KI and APP/PS1. FEA of Gene Ontology biological processes revealed 92% overlap of terms enriched in hAβ-KI with LOAD, but only 32% with EOAD. Unexpectedly, both 5xFAD and APP/PS1 presented more than 75% of their enriched terms in common with LOAD patients compared to ∼45% in EOAD. New animal models are being introduced constantly to better simulate AD pathology. Hence, the evaluation of their overlap with human pathology is paramount. We found that the new hAβ-KI shows a good fit towards LOAD, as was its original purpose. APP/PS1 and 5xFAD on the other hand, were less optimal than expected to model EOAD.

Single Cell RNA Sequencing Driven Characterization of Rare B/Myeloid and T/Myeloid Mixed Phenotype Acute Leukemia

Introduction: Pediatric mixed phenotype acute leukemia (MPAL), a rare subgroup of leukemia, contains features of both myeloid and lymphoid lineage blasts, which makes the disease more difficult to diagnose/treat. More information is needed to understand the origins of the major pediatric MPAL subtypes, B/Myeloid (B-MPAL) and T/Myeloid (T-MPAL), and how they relate to other leukemias. Single-cell RNA sequencing (scRNA-seq) analysis of bone marrow (BM) can provide in-depth information about the leukemia microenvironment and reveal differences/similarities between MPAL subtypes and other types of leukemia that could be exploited to develop novel diagnostics/therapies.Methods: We analyzed ~16,000 cells from five pediatric MPAL BM samples to generate a transcriptomic landscape of B-MPAL and T-MPAL blasts and associated microenvironment cells. Samples collected at the time of diagnosis (Dx) were used to generate scRNA-Seq data using a droplet-based barcoding technique (Panigrahy et al. JCI 2019, Tellechea et al. JID, 2020). After data normalization, cell clusters were identified using principal component analysis (PCA) and Uniform Manifold Approximation and Projection (UMAP) approach (Becht et al. Nat. Biotech 2018). Meta-analysis was performed using single cell samples from ongoing studies in the Bhasin lab (Bhasin, et al. Blood 2020 (ASH), Thomas et al. Blood 2020 (ASH)) and publicly available single cell data from GEO biorepository. Unsupervised analysis using UMAP and PCA was performed to determine the overall relationship among B-MPAL, T-MPAL and other leukemias (acute myeloid leukemia (AML), B-cell acute lymphoblastic leukemia (B-ALL), T-cell ALL (T-ALL)). Supervised differentially expressed gene (DEG) analysis was performed to identify B- and T- MPAL blast cell signatures (P value < 0.001 and log2 FC > 0.5). Transcriptomic profiles in MPAL samples and normal BM stem and immune cells were compared using data from the Human Cell Atlas Data Portal (humancellatlas.org). Gene set enrichment analysis (GSEA) was performed, and significantly enriched pathways were compared in MPAL subtypes (P value < 0.001).Results: PCA analysis showed transcriptome similarity between B-MPAL and both B-ALL and AML, while T-MPAL transcriptome correlated with T-ALL and AML (Fig. 1A). B- and T-MPAL subtype blasts clustered separately from each other in UMAP analyses, with T-MPAL blasts clustering with T-ALL blasts, and B-MPAL somewhat overlapping with B-ALL blasts. Subtype DEG analysis of leukemia blasts and healthy BM revealed distinct significantly upregulated gene signatures in B-MPAL (YBX3, SOCS2, BCL11A, and HIST1H1C) and T-MPAL (ITM2A, HPGD, PDLIM1, and TRDC) blasts (Fig. 1B). Pathway analysis showed upregulated gene activity related to TNFA signaling via NFKB, B-cell survival, and the AP1, FRA, and NGF transcription factors in B-MPAL blasts. In contrast, IL-17 and IL-12, T-cell apoptosis, and Stathmin pathways were upregulated in T-MPAL blasts (Fig. 1C). T-MPAL T-cells also expressed higher levels of T-cell exhaustion markers compared to T-cells in B-MPAL samples and healthy bone marrow.After filtering out genes that are significantly expressed in immune cells, we identified genes that are differentially expressed at diagnosis in MPAL blasts from patients that relapsed after treatment (Dx1) versus patients that achieved remission (Dx2). These genes are potential prognostic markers for B-MPAL and T-MPAL relapse/remission. These include MDM2 and NEIL1 from Dx1 and FOSL2 and CDKN1A in Dx2 B-MPAL blasts. In T-MPAL, expression of HES4 and SPINK2 is associated with Dx1 blasts and GNAQ and ITGA4 with Dx2 blasts. Pathway enrichment analysis on B-MPAL blasts revealed upregulation of interferon gamma and PD-1 signaling in Dx1 samples and increased HSP27 and Cell Cycle pathways in the Dx2 subset. T-MPAL Dx1 associated pathways included prostaglandin synthesis and IL-17, while cell-cell junction and extracellular matrix interactions were increased in T-MPAL Dx2 samples (Fig. 1D).Conclusion: Single-cell profiling was used to characterize the molecular landscapes of MPAL blasts and the bone marrow microenvironment and identified gene signatures and pathways that are specifically enriched in B- and T-MPAL subtypes. [Display omitted] DisclosuresDeRyckere: Meryx: Other: Equity ownership. Graham: Meryx: Membership on an entity's Board of Directors or advisory committees, Other: Equity ownership.

Direct reprogramming of human Sertoli cells into male germline stem cells with the self-renewal and differentiation potentials via overexpressing DAZL/DAZ2/BOULE genes.

SummaryWe propose a new concept that human somatic cells can be converted to become male germline stem cells by the defined factors. Here, we demonstrated that the overexpression of DAZL, DAZ2, and BOULE could directly reprogram human Sertoli cells into cells with the characteristics of human spermatogonial stem cells (SSCs), as shown by their similar transcriptomes and proteomics with human SSCs. Significantly, human SSCs derived from human Sertoli cells colonized and proliferated in vivo, and they could differentiate into spermatocytes and haploid spermatids in vitro. Human Sertoli cell-derived SSCs excluded Y chromosome microdeletions and assumed normal chromosomes. Collectively, human somatic cells could be converted directly to human SSCs with the self-renewal and differentiation potentials and high safety. This study is of unusual significance, because it provides an effective approach for reprogramming human somatic cells into male germ cells and offers invaluable male gametes for treating male infertility.

Single-cell and single-nucleus RNA-seq uncovers shared and distinct axes of variation in dorsal LGN neurons in mice, non-human primates, and humans.

Abundant evidence supports the presence of at least three distinct types of thalamocortical (TC) neurons in the primate dorsal lateral geniculate nucleus (dLGN) of the thalamus, the brain region that conveys visual information from the retina to the primary visual cortex (V1). Different types of TC neurons in mice, humans, and macaques have distinct morphologies, distinct connectivity patterns, and convey different aspects of visual information to the cortex. To investigate the molecular underpinnings of these cell types, and how these relate to differences in dLGN between human, macaque, and mice, we profiled gene expression in single nuclei and cells using RNA-sequencing. These efforts identified four distinct types of TC neurons in the primate dLGN: magnocellular (M) neurons, parvocellular (P) neurons, and two types of koniocellular (K) neurons. Despite extensively documented morphological and physiological differences between M and P neurons, we identified few genes with significant differential expression between transcriptomic cell types corresponding to these two neuronal populations. Likewise, the dominant feature of TC neurons of the adult mouse dLGN is high transcriptomic similarity, with an axis of heterogeneity that aligns with core vs. shell portions of mouse dLGN. Together, these data show that transcriptomic differences between principal cell types in the mature mammalian dLGN are subtle relative to the observed differences in morphology and cortical projection targets. Finally, alignment of transcriptome profiles across species highlights expanded diversity of GABAergic neurons in primate versus mouse dLGN and homologous types of TC neurons in primates that are distinct from TC neurons in mouse.

Multicellular Human Cardiac Organoids Transcriptomically Model Distinct Tissue-Level Features of Adult Myocardium.

Human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have been widely used for disease modeling and drug cardiotoxicity screening. To this end, we recently developed human cardiac organoids (hCOs) for modeling human myocardium. Here, we perform a transcriptomic analysis of various in vitro hiPSC-CM platforms (2D iPSC-CM, 3D iPSC-CM and hCOs) to deduce the strengths and limitations of these in vitro models. We further compared iPSC-CM models to human myocardium samples. Our data show that the 3D in vitro environment of 3D hiPSC-CMs and hCOs stimulates the expression of genes associated with tissue formation. The hCOs demonstrated diverse physiologically relevant cellular functions compared to the hiPSC-CM only models. Including other cardiac cell types within hCOs led to more transcriptomic similarities to adult myocardium. hCOs lack matured cardiomyocytes and immune cells, which limits a complete replication of human adult myocardium. In conclusion, 3D hCOs are transcriptomically similar to myocardium, and future developments of engineered 3D cardiac models would benefit from diversifying cell populations, especially immune cells.

A targeted transcriptomics approach for the determination of mixture effects of pesticides

While real-life exposure occurs to complex chemical mixtures, toxicological risk assessment mostly focuses on individual compounds. There is an increasing demand for in vitro tools and strategies for mixture toxicity analysis. Based on a previously established set of hepatotoxicity marker genes, we analyzed mixture effects of non-cytotoxic concentrations of different pesticides in exposure-relevant binary mixtures in human HepaRG hepatocarcinoma cells using targeted transcriptomics. An approach for mixture analysis at the level of a complex endpoint such as a transcript pattern is presented, including mixture design based on relative transcriptomic potencies and similarities. From a mechanistic point of view, goal of the study was to evaluate combinations of chemicals with varying degrees of similarity in order to determine whether differences in mechanisms of action lead to different mixtures effects. Using a model deviation ratio-based approach for assessing mixture effects, it was revealed that most data points are consistent with the assumption of dose addition. A tendency for synergistic effects was only observed at high concentrations of some combinations of the test compounds azoxystrobin, cyproconazole, difenoconazole, propiconazole and thiacloprid, which may not be representative of human real-life exposure. In summary, the findings of our study suggest that, for the pesticide mixtures investigated, risk assessment based on the general assumption of dose addition can be considered sufficiently protective for consumers. The way of data analysis presented in this paper can pave the way for a more comprehensive use of multi-gene expression data in experimental studies related to mixture toxicity.

Comparative Transcriptomics Reveals Distinct Patterns of Gene Expression Conservation through Vertebrate Embryogenesis.

Despite life’s diversity, studies of variation often remind us of our shared evolutionary past. Abundant genome sequencing and analyses of gene regulatory networks illustrate that genes and entire pathways are conserved, reused, and elaborated in the evolution of diversity. Predating these discoveries, 19th-century embryologists observed that though morphology at birth varies tremendously, certain stages of vertebrate embryogenesis appear remarkably similar across vertebrates. In the mid to late 20th century, anatomical variability of early and late-stage embryos and conservation of mid-stages embryos (the “phylotypic” stage) was named the hourglass model of diversification. This model has found mixed support in recent analyses comparing gene expression across species possibly owing to differences in species, embryonic stages, and gene sets compared. We compare 186 microarray and RNA-seq data sets covering embryogenesis in six vertebrate species. We use an unbiased clustering approach to group stages of embryogenesis by transcriptomic similarity and ask whether gene expression similarity of clustered embryonic stages deviates from a null expectation. We characterize expression conservation patterns of each gene at each evolutionary node after correcting for phylogenetic nonindependence. We find significant enrichment of genes exhibiting early conservation, hourglass, late conservation patterns in both microarray and RNA-seq data sets. Enrichment of genes showing patterned conservation through embryogenesis indicates diversification of embryogenesis may be temporally constrained. However, the circumstances under which each pattern emerges remain unknown and require both broad evolutionary sampling and systematic examination of embryogenesis across species.

Time-course transcriptome landscape of achene development in lettuce

Lettuce (Lactuca sativa L.), which belongs to the large Asteraceae (Compositae) family, breeds by sexual reproduction and produces seeds. Actually, lettuce seeds are achenes, which are defined as fruits. However, few studies have described the morphological characteristics of the lettuce achenes, and genes essential for achene development are largely unknown in lettuce. To investigate the gene activity during achene development and determine the possible mechanisms that influence achene development in lettuce, we performed a time-course transcriptome analysis of lettuce achenes. A total of 27,390 expressed genes were detected at the five achene development stages. We investigated the gene expression patterns during achene development and identified the enriched biological processes at the corresponding stages. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses revealed a variety of transcriptomic similarities and differences at different achene development stages. Further, transcription factors and phytohormones were found to play important roles during achene development. Finally, we proposed a working model to illustrate the gene expression modules and possible molecular mechanisms underlying achene development. Our time-course transcriptome data also provide a foundation for future functional studies to reveal the genetic control of achene development in lettuce.

Highly conserved sperm function-related transcripts across three species: human, rat and mouse

Assessing male reproductive toxicity of environmental and therapeutic agents relies on the histopathology of the testis and epididymis in a pre-clinical setting. Animal histopathology poorly correlates with human sperm parameters, and none of these current methods are strong indicators of sperm health or reproductive potential. Therefore, there is an urgent need to identify a translatable, non-invasive and reliable approach to monitor environmental and therapeutic agents’ effects on male reproductive health. mRNA sequences were analyzed in mouse, rat and human sperm samples to identify sperm transcriptomic similarities across species that could be used as biomarkers to predict male reproductive toxicity in animal models. Semen specimens were collected from men aged 18 to 55 years with proven fertility. Rat and mouse semen specimens were collected via needle punctures of the cauda epididymides. Sperm RNAs were extracted using an optimized sperm RNA isolation protocol and subjected to polyA-purified mRNA-sequencing. Bioinformatics analyses, including differential abundance and gene set enrichment analysis, were used to investigate the biological and molecular functions of all shared and differentially abundant transcripts across species. Transcriptome profiling identified 6,684 similarly expressed transcripts within the three species of which 1,579 transcripts were found to be involved in spermatogenic functions. Our findings have shown that sperm transcriptome is highly species dependent, however, there are some key similarities among transcripts that are required for fertility. Based on these similarities, sperm mRNA biomarker may be developed to monitor male reproductive toxicity where rodent models would make suitable laboratory substitutes for human.

Missing Nurse Bees-Early Transcriptomic Switch From Nurse Bee to Forager Induced by Sublethal Imidacloprid.

The environmental residue/sublethal doses of neonicotinoid insecticides are believed to generate a negative impact on pollinators, including honey bees. Here we report our recent investigation on how imidacloprid, one of the major neonicotinoids, affects worker bees by profiling the transcriptomes of various ages of bees exposed to different doses of imidacloprid during the larval stage. The results show that imidacloprid treatments during the larval stage severely altered the gene expression profiles and may induce precocious foraging. Differential expression of foraging regulators was found in 14-day-old treated adults. A high transcriptome similarity between larvae-treated 14-day-old adults and 20-day-old controls was also observed, and the similarity was positively correlated with the dose of imidacloprid. One parts per billion (ppb) of imidacloprid was sufficient to generate a long-term impact on the bee’s gene expression as severe as with 50 ppb imidacloprid. The disappearance of nurse bees may be driven not only by the hive member constitution but also by the neonicotinoid-induced precocious foraging behavior.