Progress towards an OECD reporting framework for transcriptomics and metabolomics in regulatory toxicology.

Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.

Molecular Assessment of Kidney Allografts: Are We Closer to a Daily Routine?

Kidney allograft pathology assessment has been traditionally based on clinical and histological criteria. Despite improvements in Banff histological classification, the diagnostics in particular cases is problematic reflecting a complex pathogenesis of graft injuries. With the advent of molecular techniques, polymerase-chain reaction, oligo- and microarray technologies allowed to study molecular phenotypes of graft injuries, especially acute and chronic rejections. Moreover, development of the molecular microscope diagnostic system (MMDx) to assess kidney graft biopsies, represents the first clinical application of a microarray-based method in transplantation. Whether MMDx may replace conventional pathology is the subject of ongoing research, however this platform is particularly useful in complex histological findings and may help clinicians to guide the therapy.

Concentration-dependent transcriptome of zebrafish embryo for environmental chemical assessment.

Mechanistic information is essential to screen and predict the adverse effects of a large number of chemicals during early-life exposure. Concentration-dependent omics can capture the extent of perturbations of biological pathways or processes and provide information on the mechanism of toxicity. However, the application of concentration-dependent transcriptome to assess the developmental toxicity of environmental chemicals is still limited. Here, twelve chemicals representing five different modes of action (MOAs) were tested by the concentration-dependent reduced zebrafish transcriptome approach (CRZT) in combination with a phenotype-based high content screen (PHCS). The responsiveness, sensitivity and mechanistic differentiation of CRZT were validated in comparison with PHCS. First, PHCS identified 10 chemicals with obvious embryotoxicity (LD50 range: 2.11-70.68 µM), while the potencies of the biological pathways perturbed by 12 chemicals (PODpath20 range: 0.002-2.1 µM) were demonstrated by CRZT. Second, although the potency of the transcriptome perturbations was positively correlated with lethality (LD50) (R2 = 0.64, P-value < 0.05) for most tested chemicals, BbF was non-embryotoxic but was the most potent on the perturbance of biological pathways. Finally, the profiles of the perturbed biological processes and the transcriptome potency (PODpath20) captured by CRZT could effectively classify most chemicals corresponding to their known MOAs. In summary, CRZT could significantly improve testing the developmental toxicity of environmental chemicals.

Age and poverty status alter the coding and noncoding transcriptome.

Emerging evidence indicates that noncoding RNAs play regulatory roles in aging and disease. The functional roles of long noncoding RNAs (lncRNAs) in physiology and disease are not completely understood. Little is known about lncRNAs in the context of human aging and socio-environmental conditions. Microarray profiling of lncRNAs and mRNAs from peripheral blood mononuclear cells from young and old white (n=16) and African American (AA) males (n=16) living above or below poverty from the Healthy Aging in Neighborhoods of Diversity across the Life Span study revealed changes in both lncRNAs and mRNAs with age and poverty status in white males, but not in AA males. We validated lncRNA changes in an expanded cohort (n=40); CTD-3247F14.2, GAS5, H19, TERC and MEG3 changed significantly with age, whereas AK022914,GAS5, KB-1047C11.2, MEG3 and XLOC_003262 changed with poverty. Mitochondrial function and response to DNA damage and stress were pathways enriched in younger individuals. Response to stress, viral infection, and immune signals were pathways enriched in individuals living above poverty. These data show that both human age and a marker of social adversity influence lncRNA expression, which may provide insight about molecular pathways underlying aging and social factors that affect disparities in aging and disease.

A transcriptome-based assessment of the astrocytic dystrophin-associated complex in the developing human brain.

Astrocytes play a critical role in regulating the interface between the cerebral vasculature and the central nervous system. Contributing to this is the astrocytic endfoot domain, a specialized structure that ensheathes the entirety of the vasculature and mediates signaling between endothelial cells, pericytes, and neurons. The astrocytic endfoot has been implicated as a critical element of the glymphatic pathway, and changes in protein expression profiles in this cellular domain are linked to Alzheimer's disease pathology. Despite this, basic physiological properties of this structure remain poorly understood including the developmental timing of its formation, and the protein components that localize there to mediate its functions. Here we use human transcriptome data from male and female subjects across several developmental stages and brain regions to characterize the gene expression profile of the dystrophin-associated complex (DAC), a known structural component of the astrocytic endfoot that supports perivascular localization of the astroglial water channel aquaporin-4. Transcriptomic profiling is also used to define genes exhibiting parallel expression profiles to DAC elements, generating a pool of candidate genes that encode gene products that may contribute to the physiological function of the perivascular astrocytic endfoot domain. We found that several genes encoding transporter proteins are transcriptionally associated with DAC genes.

Assessment of estradiol-induced gene regulation and proliferation in an immortalized mouse immature Sertoli cell line.

AIMS: The number of Sertoli cells during proliferative phase determines the fate of the germ cells in male reproductive system. A well-characterized cell line may help in better understanding of Sertoli cell biology. Hence, the present study assessed estradiol signaling in a mouse immature Sertoli cell line (MSC-1) as an alternative model in place of primary culture of Sertoli cells. MAIN METHODS: In this study, we used MSC-1 cell line, derived from 10-day old mice. The cell cycle parameters were assessed, and the expression and regulation of Sertoli cell-specific secretory genes (ABP; androgen-binding protein) and tight junction genes (claudin-5, occludin, and vimentin) in response to estradiol was studied. KEY FINDINGS: The results obtained suggested the presence of both estrogen receptors (ERα and ERß) in MSC-1 cells. In vitro scratch assay and cell-cycle analysis suggested the proliferative effects of estradiol in both time- and dose-dependent manner. The gene expression profiles of ABP, claudin-5, and occludin showed biphasic regulation at low and high doses of estradiol. Analysis of signaling pathways suggested the activation of extracellular signal-regulated kinase (ERK) pathway with significantly increased pERK/ERK ratio (p<0.05). The results also suggested down regulation in the expression of mir-17 family members (mir-17, mir-20b, and mir-106a) (p<0.05). SIGNIFICANCE: Considering the limited number of Sertoli cell lines and long-term survival inability of primary culture of Sertoli cells, MSC-1 cells could be a potential cell line for understanding the mechanisms of various cellular events in Sertoli cells.

Transcriptional assessment by microarray analysis and large-scale meta-analysis of the metabolic capacity of cardiac and skeletal muscle tissues to cope with reduced nutrient availability in Gilthead Sea Bream (Sparus aurata L.).

The effects of nutrient availability on the transcriptome of cardiac and skeletal muscle tissues were assessed in juvenile gilthead sea bream fed with a standard diet at two feeding levels: (1) full ration size and (2) 70 % satiation followed by a finishing phase at the maintenance ration. Microarray analysis evidenced a characteristic transcriptomic profile for each muscle tissue following changes in oxidative capacity (heart > red skeletal muscle > white skeletal muscle). The transcriptome of heart and secondly that of red skeletal muscle were highly responsive to nutritional changes, whereas that of glycolytic white skeletal muscle showed less ability to respond. The highly expressed and nutritionally regulated genes of heart were mainly related to signal transduction and transcriptional regulation. In contrast, those of white muscle were enriched in gene ontology (GO) terms related to proteolysis and protein ubiquitination. Microarray meta-analysis using the bioinformatic tool Fish and Chips ( http://fishandchips.genouest.org/index.php ) showed the close association of a representative cluster of white skeletal muscle with some of cardiac and red skeletal muscle, and many GO terms related to mitochondrial function appeared to be common links between them. A second round of cluster comparisons revealed that mitochondria-related GOs also linked differentially expressed genes of heart with those of liver from cortisol-treated gilthead sea bream. These results show that mitochondria are among the first responders to environmental and nutritional stress stimuli in gilthead sea bream, and functional phenotyping of this cellular organelle is highly promising to obtain reliable markers of growth performance and well-being in this fish species.

Assessment of myometrial transcriptome changes associated with spontaneous human labour by high-throughput RNA-seq.

The transition of the human uterus from a quiescent to a contractile state takes place over a number of weeks. On such biological time scales, cellular phenotype is modified by changes in the transcriptome, which in turn is under the control of the underlying endocrine, paracrine, and biophysical processes resulting from the ongoing pregnancy. In this study, we characterize the transition of the human myometrial transcriptome at term from not in labour (NIL) to in labour (LAB) using high throughput RNA sequencing (RNA-seq). RNA was isolated from the myometrium of uterine biopsies from patients at term who were not in labour (n = 5) and at term in spontaneous labour (n = 5) without augmentation. A total of 143.6 million separate reads were sequenced, achieving, on average, ∼13 times coverage of the expressed human transcriptome per sample. Principal component analysis indicated that the NIL and LAB transcriptomes could be distinguished as two distinct clusters. A comparison of the NIL and LAB groups, using three different statistical approaches (baySeq, edgeR, and DESeq), demonstrated an overlap of 764 differentially expressed genes. A comparison with currently available microarray data revealed only a partial overlap in differentially expressed genes. We conclude that the described RNA-seq data sets represent the first fully annotated catalogue of expressed mRNAs in human myometrium. When considered together, the full expression repertoire and the differentially expressed gene sets should provide an excellent resource for formulating new hypotheses of physiological function, as well as the discovery of novel therapeutic targets.

Salivaomics.

BACKGROUND: The ability to monitor health and wellness, as well as detect oral and systemic illnesses early through noninvasive means, are highly desirable goals in health care promotion and delivery. Saliva is an emerging medium to be explored for health and disease surveillance, as well as for personalized medicine. A major mandate is to demonstrate clinicians' ability to use saliva to detect and monitor systemic diseases. METHODS: To realize the translational and clinical vision of salivary diagnostics, two prerequisites are essential. The first is the need to develop and optimize diagnostic tools tailored to saliva. The second is the need to substantiate the scientific underpinnings of salivary biomarkers reflecting systemic diseases. RESULTS: The author describes five diagnostic alphabets (proteome, transcriptome, microRNA, metabolome and microbiome) and point-of-care technology platforms that are in place to advance the translational and clinical path. For mechanistic studies (that is, basic science studies), animal models are in place to elucidate the scientific mechanisms of systemic diseases reflected in saliva. CONCLUSIONS: Significant advancements have been made in the development of salivary diagnostic tools. The translation of the scientific mechanisms of systemic diseases reflected in saliva is in progress. CLINICAL IMPLICATIONS: On the scientific credentialing of salivary biomarkers for the detection of systemic diseases, salivary diagnostics will have an effect on access to care, health disparities and global health. Dentistry can advance into the realm of primary health care with integration of chairside screening for medical conditions.

An automated microfluidic device for assessment of mammalian cell genetic stability.

Single-cell transcriptome contains reliable gene regulatory relationships because gene-gene interactions only happen within a mammalian cell. While the study of gene-gene interactions enables us to understand the molecular mechanism of cellular events and evaluate molecular characteristics of a mammalian cell population, its complexity requires an analysis of a large number of single-cells at various stages. However, many existing microfluidic platforms cannot process single-cells effectively for routine molecular analysis. To address these challenges, we develop an integrated system with individual controller for effective single-cell transcriptome analysis. In this paper, we report an integrated microfluidic approach to rapidly measure gene expression in individual cells for genetic stability assessment of a cell population. Inside this integrated microfluidic device, the cells are individually manipulated and isolated in an array using micro sieve structures, then transferred into different nanoliter reaction chambers for parallel processing of single-cell transcriptome analysis. This device enables us to manipulate individual single-cells into nanoliter reactor with high recovery rate. We have performed gene expression analysis for a large number of HeLa cells and 293T cells expanded from a single-cell. Our data shows that even the house-keeping genes are expressed at heterogeneous levels within a clone of cells. The heterogeneity of actin expression reflects the genetic stability, and the expression distribution is different between cancer cells (HeLa) and immortalized 293T cells. The result demonstrates that this platform has the potential for assessment of genetic stability in cancer diagnosis.

IQSeq: integrated isoform quantification analysis based on next-generation sequencing.

With the recent advances in high-throughput RNA sequencing (RNA-Seq), biologists are able to measure transcription with unprecedented precision. One problem that can now be tackled is that of isoform quantification: here one tries to reconstruct the abundances of isoforms of a gene. We have developed a statistical solution for this problem, based on analyzing a set of RNA-Seq reads, and a practical implementation, available from archive.gersteinlab.org/proj/rnaseq/IQSeq, in a tool we call IQSeq (Isoform Quantification in next-generation Sequencing). Here, we present theoretical results which IQSeq is based on, and then use both simulated and real datasets to illustrate various applications of the tool. In order to measure the accuracy of an isoform-quantification result, one would try to estimate the average variance of the estimated isoform abundances for each gene (based on resampling the RNA-seq reads), and IQSeq has a particularly fast algorithm (based on the Fisher Information Matrix) for calculating this, achieving a speedup of ~ 500 times compared to brute-force resampling. IQSeq also calculates an information theoretic measure of overall transcriptome complexity to describe isoform abundance for a whole experiment. IQSeq has many features that are particularly useful in RNA-Seq experimental design, allowing one to optimally model the integration of different sequencing technologies in a cost-effective way. In particular, the IQSeq formalism integrates the analysis of different sample (i.e. read) sets generated from different technologies within the same statistical framework. It also supports a generalized statistical partial-sample-generation function to model the sequencing process. This allows one to have a modular, "plugin-able" read-generation function to support the particularities of the many evolving sequencing technologies.

Molecular profiling of marine fauna: integration of omics with environmental assessment of the world's oceans.

Many species that contribute to the commercial and ecological richness of our marine ecosystems are harbingers of environmental change. The ability of organisms to rapidly detect and respond to changes in the surrounding environment represents the foundation for application of molecular profiling technologies towards marine sentinel species in an attempt to identify signature profiles that may reside within the transcriptome, proteome, or metabolome and that are indicative of a particular environmental exposure event. The current review highlights recent examples of the biological information obtained for marine sentinel teleosts, mammals, and invertebrates. While in its infancy, such basal information can provide a systems biology framework in the detection and evaluation of environmental chemical contaminant effects on marine fauna. Repeated evaluation across different seasons and local marine environs will lead to discrimination between signature profiles representing normal variation within the complex milieu of environmental factors that trigger biological response in a given sentinel species and permit a greater understanding of normal versus anthropogenic-associated modulation of biological pathways, which prove detrimental to marine fauna. It is anticipated that incorporation of contaminant-specific molecular signatures into current risk assessment paradigms will lead to enhanced wildlife management strategies that minimize the impacts of our industrialized society on marine ecosystems.

Relating CNVs to transcriptome data at fine resolution: assessment of the effect of variant size, type, and overlap with functional regions.

Copy-number variants (CNVs) form an abundant class of genetic variation with a presumed widespread impact on individual traits. While recent advances, such as the population-scale sequencing of human genomes, facilitated the fine-scale mapping of CNVs, the phenotypic impact of most of these CNVs remains unclear. By relating copy-number genotypes to transcriptome sequencing data, we have evaluated the impact of CNVs, mapped at fine scale, on gene expression. Based on data from 129 individuals with ancestry from two populations, we identified CNVs associated with the expression of 110 genes, with 13% of the associations involving complex, multiallelic CNVs. Categorization of CNVs according to variant type, size, and gene overlap enabled us to examine the impact of different CNV classes on expression variation. While many small (<4 kb) CNVs were associated with expression variation, overall we observed an enrichment of large duplications and deletions, including large intergenic CNVs, relative to the entire set of expression-associated CNVs. Furthermore, the copy number of genes intersecting with CNVs typically correlated positively with the genes' expression, and also was more strongly correlated with expression than nearby single nucleotide polymorphisms, suggesting a frequent causal role of CNVs in expression quantitative trait loci (eQTLs). We also elucidated unexpected cases of negative correlations between copy number and expression by assessing the CNVs' effects on the structure and regulation of genes. Finally, we examined dosage compensation of transcript levels. Our results suggest that association studies can gain in resolution and power by including fine-scale CNV information, such as those obtained from population-scale sequencing.

Combining resources to obtain a comprehensive survey of the bovine embryo transcriptome through deep sequencing and microarrays.

While most assisted reproductive technologies (ART) are considered routine for the reproduction of species of economical importance, such as the bovine, the impact of these manipulations on the developing embryo remains largely unknown. In an effort to obtain a comprehensive survey of the bovine embryo transcriptome and how it is modified by ART, resources were combined to design an embryo-specific microarray. Close to one million high-quality reads were produced from subtracted bovine embryo libraries using Roche 454 Titanium deep sequencing technology, which enabled the creation of an augmented bovine genome catalog. This catalog was enriched with bovine embryo transcripts, and included newly discovered indel type and 3'UTR variants. Using this augmented bovine genome catalog, the EmbryoGENE Bovine Microarray was designed and is composed of a total of 42,242 probes, including 21,139 known reference genes; 9,322 probes for novel transcribed regions (NTRs); 3,677 alternatively spliced exons; 3,353 3'-tiling probes; and 3,723 controls. A suite of bioinformatics tools was also developed to facilitate microrarray data analysis and database creation; it includes a quality control module, a Laboratory Information Management System (LIMS) and microarray analysis software. Results obtained during this study have already led to the identification of differentially expressed blastocyst targets, NTRs, splice variants of the indel type, and 3'UTR variants. We were able to confirm microarray results by real-time PCR, indicating that the EmbryoGENE bovine microarray has the power to detect physiologically relevant changes in gene expression.