Measuring the exposome: a powerful basis for evaluating environmental exposures and cancer risk

Advances in laboratory sciences offer much in the challenge to unravel the complex etiology of cancer and to therefore provide an evidence-base for prevention. One area where improved measurements are particularly important to epidemiology is exposure assessment; this requirement has been highlighted through the concept of the exposome. In addition, the ability to observe genetic and epigenetic alterations in individuals exposed to putative risk factors also affords an opportunity to elucidate underlying mechanisms of carcinogenesis, which in turn may allow earlier detection and more refined molecular classification of disease. In this context the application of omics technologies to large population-based studies and their associated biobanks raise exciting new avenues of research. This review considers the areas of genomics, transcriptomics, epigenomics and metabolomics and the evidence to date that people exposed to well-defined factors (for example, tobacco, diet, occupational exposures, environmental pollutants) have specific omics profiles. Although in their early stages of development these approaches show promising evidence of distinct exposure-derived biological effects and indicate molecular pathways that may be particularly relevant to the carcinogenic process subsequent to environmental and lifestyle exposures. Such an interdisciplinary approach is vital if the full benefits of advances in laboratory sciences and investments in large-scale prospective cohort studies are to be realized in relation to cancer prevention.

Making sense of OMICS data in population-based environmental health studies

Although experience from the application of OMICS technologies in population-based environmental health studies is still relatively limited, the accumulated evidence shows that it can allow the identification of features (genes, proteins, and metabolites), or sets of such features, which are targeted by particular exposures or correlate with disease risk. Such features or profiles can therefore serve as biomarkers of exposure or disease risk. Blood-based OMIC profiles appear to reflect to some extent events occurring in target tissues and are associated with toxicity or disease and therefore have the potential to facilitate the elucidation of exposure-disease relationships. Further progress in this direction requires better understanding of the significance of exposure-induced network perturbations for disease initiation and progression and the development of a framework that combines agnostic searches with the utilization of prior knowledge, taking account of particular elements which characterize the structure and evolution of complex systems and brings in principles of systems biology.

Transcriptome analysis of the sulfate deficiency response in the marine microalga Emiliania huxleyi

The response to sulfate deficiency of plants and freshwater green algae has been extensively analysed by system biology approaches. By contrast, seawater sulfate concentration is high and very little is known about the sulfur metabolism of marine organisms. Here, we used a combination of metabolite analysis and transcriptomics to analyse the response of the marine microalga Emiliania huxleyi as it acclimated to sulfate limitation. Lowering sulfate availability in artificial seawater from 25 to 5 mM resulted in significant reduction in growth and intracellular concentrations of dimethylsulfoniopropionate and glutathione. Sulfate-limited E. huxleyi cells showed increased sulfate uptake but sulfate reduction to sulfite did not seem to be regulated. Sulfate limitation in E. huxleyi affected expression of 1718 genes. The vast majority of these genes were upregulated, including genes involved in carbohydrate and lipid metabolism, and genes involved in the general stress response. The acclimation response of E. huxleyi to sulfate deficiency shows several similarities to the well-described responses of Arabidopsis and Chlamydomonas, but also has many unique features. This dataset shows that even though E. huxleyi is adapted to constitutively high sulfate concentration, it retains the ability to re-program its gene expression in response to reduced sulfate availability.

Regulated expression of polysaccharide utilization and capsular biosynthesis loci in biofilm and planktonic Bacteroides thetaiotaomicron during growth in chemostats

Bacteroides thetaiotaomicron is a prominent member of the human distal gut microbiota that specializes in breaking down diet and host-derived polysaccharides. While polysaccharide utilization has been well studied in B. thetaiotaomicron, other aspects of its behavior are less well characterized, including the factors that allow it to maintain itself in the gut. Biofilm formation may be a mechanism for bacterial retention in the gut. Therefore, we used custom GeneChips to compare the transcriptomes of biofilm and planktonic B. thetaiotaomicron during growth in mono-colonized chemostats. We identified 1,154 genes with a fold-change greater than 2, with confidence greater than or equal to 95%. Among the prominent changes observed in biofilm populations were: (i) greater expression of genes in polysaccharide utilization loci that are involved in foraging of O-glycans normally found in the gut mucosa; and (ii) regulated expression of capsular polysaccharide biosynthesis loci. Hierarchical clustering of the data with different datasets, which were obtained during growth under a range of conditions in minimal media and in intestinal tracts of gnotobiotic mice, revealed that within this group of differentially expressed genes, biofilm communities were more similar to the in vivo samples than to planktonic cells and exhibited features of substrate limitation. The current study also validates the use of chemostats as an in vitro "gnotobiotic" model to study gene expression of attached populations of this bacterium. This is important to gut microbiota research, because bacterial attachment and the consequences of disruptions in attachment are difficult to study in vivo.

Reduction of photosynthetic sensitivity in response to abiotic stress in tomato is mediated by a new generation plant activator

BACKGROUND

Yield losses as a result of abiotic stress factors present a significant challenge for the future of global food production. While breeding technologies provide potential to combat negative stress-mediated outcomes over time, interventions which act to prime plant tolerance to stress, via the use of phytohormone-based elicitors for example, could act as a valuable tool for crop protection. However, the translation of fundamental biology into functioning solution is often constrained by knowledge-gaps.

RESULTS

Photosynthetic and transcriptomic responses were characterised in young tomato (Solanum lycopersicum L.) seedlings in response to pre-treatment with a new plant health activator technology, 'Alethea', followed by a subsequent 100 mM salinity stress. Alethea is a novel proprietary technology composed of three key constituent compounds; the hitherto unexplored compound potassium dihydrojasmonate, an analogue of jasmonic acid; sodium benzoate, a carboxylic acid precursor to salicylic acid, and the α-amino acid L-arginine. Salinity treatment led to a maximal 47% reduction in net photosynthetic rate 8 d following NaCl treatment, yet in Alethea pre-treated seedlings, sensitivity to salinity stress was markedly reduced during the experimental period. Microarray analysis of leaf transcriptional responses showed that while salinity stress and Alethea individually impacted on largely non-overlapping, distinct groups of genes, Alethea pre-treatment substantially modified the response to salinity. Alethea affected the expression of genes related to biotic stress, ethylene signalling, cell wall synthesis, redox signalling and photosynthetic processes. Since Alethea had clear effects on photosynthesis/chloroplastic function at the physiological and molecular levels, we also investigated the ability of Alethea to protect various crop species against methyl viologen, a potent generator of oxidative stress in chloroplasts. Alethea pre-treatment produced dramatic reductions in visible foliar necrosis caused by methyl viologen compared with non-primed controls.

CONCLUSIONS

'Alethea' technology mediates positive recovery of abiotic stress-induced photosynthetic and foliar loss of performance, which is accompanied by altered transcriptional responses to stress.

Specific mutations in H5N1 mainly impact the magnitude and velocity of the host response in mice

BACKGROUND

Influenza infection causes respiratory disease that can lead to death. The complex interplay between virus-encoded and host-specific pathogenicity regulators - and the relative contributions of each toward viral pathogenicity - is not well-understood.

RESULTS

By analyzing a collection of lung samples from mice infected by A/Vietnam/1203/2004 (H5N1; VN1203), we characterized a signature of transcripts and proteins associated with the kinetics of the host response. Using a new geometrical representation method and two criteria, we show that inoculation concentrations and four specific mutations in VN1203 mainly impact the magnitude and velocity of the host response kinetics, rather than specific sets of up- and down- regulated genes. We observed analogous kinetic effects using lung samples from mice infected with A/California/04/2009 (H1N1), and we show that these effects correlate with morbidity and viral titer.

CONCLUSIONS

We have demonstrated the importance of the kinetics of the host response to H5N1 pathogenesis and its relationship with clinical disease severity and virus replication. These kinetic properties imply that time-matched comparisons of 'omics profiles to viral infections give limited views to differentiate host-responses. Moreover, these results demonstrate that a fast activation of the host-response at the earliest time points post-infection is critical for protective mechanisms against fast replicating viruses.

Integrated OMICS guided engineering of biofuel butanol-tolerance in photosynthetic Synechocystis sp. PCC 6803

BACKGROUND

Photosynthetic cyanobacteria have been recently proposed as a 'microbial factory' to produce butanol due to their capability to utilize solar energy and CO2 as the sole energy and carbon sources, respectively. However, to improve the productivity, one key issue needed to be addressed is the low tolerance of the photosynthetic hosts to butanol.

RESULTS

In this study, we first applied a quantitative transcriptomics approach with a next-generation RNA sequencing technology to identify gene targets relevant to butanol tolerance in a model cyanobacterium Synechocystis sp. PCC 6803. The results showed that 278 genes were induced by the butanol exposure at all three sampling points through the growth time course. Genes encoding heat-shock proteins, oxidative stress related proteins, transporters and proteins involved in common stress responses, were induced by butanol exposure. We then applied GC-MS based metabolomics analysis to determine the metabolic changes associated with the butanol exposure. The results showed that 46 out of 73 chemically classified metabolites were differentially regulated by butanol treatment. Notably, 3-phosphoglycerate, glycine, serine and urea related to general stress responses were elevated in butanol-treated cells. To validate the potential targets, we constructed gene knockout mutants for three selected gene targets. The comparative phenotypic analysis confirmed that these genes were involved in the butanol tolerance.

CONCLUSION

The integrated OMICS analysis provided a comprehensive view of the complicated molecular mechanisms employed by Synechocystis sp. PCC 6803 against butanol stress, and allowed identification of a series of potential gene candidates for tolerance engineering in cyanobacterium Synechocystis sp. PCC 6803.

Bridging the clinical gaps: genetic, epigenetic and transcriptomic biomarkers for the early detection of lung cancer in the post-National Lung Screening Trial era

Lung cancer is the leading cause of cancer death worldwide in part due to our inability to identify which smokers are at highest risk and the lack of effective tools to detect the disease at its earliest and potentially curable stage. Recent results from the National Lung Screening Trial have shown that annual screening of high-risk smokers with low-dose helical computed tomography of the chest can reduce lung cancer mortality. However, molecular biomarkers are needed to identify which current and former smokers would benefit most from annual computed tomography scan screening in order to reduce the costs and morbidity associated with this procedure. Additionally, there is an urgent clinical need to develop biomarkers that can distinguish benign from malignant lesions found on computed tomography of the chest given its very high false positive rate. This review highlights recent genetic, transcriptomic and epigenomic biomarkers that are emerging as tools for the early detection of lung cancer both in the diagnostic and screening setting.

Bioinformatics and Gene Network Analyses of the Swine Mammary Gland Transcriptome during Late Gestation

We used the newly-developed Dynamic Impact Approach (DIA) and gene network analysis to study the sow mammary transcriptome at 80, 100, and 110 days of pregnancy. A swine oligoarray with 13,290 inserts was used for transcriptome profiling. An ANOVA with false discovery rate (FDR < 0.15) correction resulted in 1,409 genes with a significant time effect across time comparisons. The DIA uncovered that Fatty acid biosynthesis, Interleukin-4 receptor binding, Galactose metabolism, and mTOR signaling were among the most-impacted pathways. IL-4 receptor binding, ABC transporters, cytokine-cytokine receptor interaction, and Jak-STAT signaling were markedly activated at 110 days compared with 80 and 100 days. Epigenetic and transcription factor regulatory mechanisms appear important in coordinating the final stages of mammary development during pregnancy. Network analysis revealed a crucial role for TP53, ARNT2, E2F4, and PPARG. The bioinformatics analyses revealed a number of pathways and functions that perform an irreplaceable role during late gestation to farrowing.

The Saccharomyces cerevisiae transcriptome as a mirror of phytochemical variation in complex extracts of Equisetum arvense from America, China, Europe and India

BACKGROUND

Pattern-oriented chemical profiling is increasingly being used to characterize the phytochemical composition of herbal medicines for quality control purposes. Ideally, a fingerprint of the biological effects should complement the chemical fingerprint. For ethical and practical reasons it is not possible to test each herbal extract in laboratory animals or humans. What is needed is a test system consisting of an organism with relevant biology and complexity that can serve as a surrogate in vitro system. The purpose of this study was to test the hypothesis that the Saccharomyces cerevisiae transcriptome might be used as an indicator of phytochemical variation of closely-related yet distinctly different extracts prepared from a single species of a phytogeographically widely distributed medicinal plant. We combined phytochemical profiling using chromatographic methods (HPTLC, HPLC-PDA-MS/MS) and gene expression studies using Affymetrix Yeast 2.0 gene chip with principal component analysis and k-nearest neighbor clustering analysis to test this hypothesis using extracts prepared from the phytogeographically widely distributed medicinal plant Equisetum arvense as a test case.

RESULTS

We found that the Equisetum arvense extracts exhibited qualitative and quantitative differences in their phytochemical composition grouped along their phytogeographical origin. Exposure of yeast to the extracts led to changes in gene expression that reflected both the similarities and differences in the phytochemical composition of the extracts. The Equisetum arvense extracts elicited changes in the expression of genes involved in mRNA translation, drug transport, metabolism of energy reserves, phospholipid metabolism, and the cellular stress response.

CONCLUSIONS

Our data show that functional genomics in S. cerevisiae may be developed as a sensitive bioassay for the scientific investigation of the interplay between phytochemical composition and transcriptional effects of complex mixtures of chemical compounds. S. cerevisiae transcriptomics may also be developed for testing of mixtures of conventional drugs ("polypills") to discover novel antagonistic or synergistic effects of those drug combinations.

Gene expression profiling of blood in ruptured intracranial aneurysms: in search of biomarkers

The molecular mechanisms underlying the systemic response to subarachnoid hemorrhage (SAH) from ruptured intracranial aneurysms (RAs) are not fully understood. We investigated whether the analysis of gene expression in peripheral blood could provide clinically relevant information regarding the biologic consequences of SAH. Transcriptomics were performed using Illumina HumanHT-12v4 microarrays for 43 RA patients and 18 controls (C). Differentially expressed transcripts were analyzed for overrepresented functional groups and blood cell type-specific gene expression. The set of differentially expressed transcripts was validated using quantitative polymerase chain reaction in an independent group of subjects (15 RA patients and 14 C). There were 135 differentially expressed genes (false discovery rate 1%, absolute fold change 1.7): the abundant levels of 78 mRNAs increased and 57 mRNAs decreased. Among RA patients, transcripts specific to T lymphocyte subpopulations were downregulated, whereas those related to monocytes and neutrophils were upregulated. Expression profiles of a set of 16 genes and lymphocyte-to-monocyte-and-neutrophil gene expression ratios distinguished RA patients from C. These results indicate that SAH from RAs strongly influences the transcription profiles of blood cells. A specific pattern of these changes suggests suppression in lymphocyte response and enhancements in monocyte and neutrophil activities. This is probably related to the immunodepression observed in SAH.

Analysis of inflammatory and lipid metabolic networks across RAW264.7 and thioglycolate-elicited macrophages

Studies of macrophage biology have been significantly advanced by the availability of cell lines such as RAW264.7 cells. However, it is unclear how these cell lines differ from primary macrophages such as thioglycolate-elicited peritoneal macrophages (TGEMs). We used the inflammatory stimulus Kdo2-lipid A (KLA) to stimulate RAW264.7 and TGEM cells. Temporal changes of lipid and gene expression levels were concomitantly measured and a systems-level analysis was performed on the fold-change data. Here we present a comprehensive comparison between the two cell types. Upon KLA treatment, both RAW264.7 and TGEM cells show a strong inflammatory response. TGEM (primary) cells show a more rapid and intense inflammatory response relative to RAW264.7 cells. DNA levels (fold-change relative to control) are reduced in RAW264.7 cells, correlating with greater downregulation of cell cycle genes. The transcriptional response suggests that the cholesterol de novo synthesis increases considerably in RAW264.7 cells, but 25-hydroxycholesterol increases considerably in TGEM cells. Overall, while RAW264.7 cells behave similarly to TGEM cells in some ways and can be used as a good model for inflammation- and immune function-related kinetic studies, they behave differently than TGEM cells in other aspects of lipid metabolism and phenotypes used as models for various disorders such as atherosclerosis.

Role of sortase-dependent pili of Bifidobacterium bifidum PRL2010 in modulating bacterium-host interactions. Proc Natl Acad Sci U S A. 2013;110:11151-11156. [PMID: 23776216 DOI: 10.1073/pnas.1303897110]

Bifidobacteria represent one of the dominant groups of microorganisms colonizing the human infant intestine. Commensal bacteria that interact with a eukaryotic host are believed to express adhesive molecules on their cell surface that bind to specific host cell receptors or soluble macromolecules. Whole-genome transcription profiling of Bifidobacterium bifidum PRL2010, a strain isolated from infant stool, revealed a small number of commonly expressed extracellular proteins, among which were genes that specify sortase-dependent pili. Expression of the coding sequences of these B. bifidum PRL2010 appendages in nonpiliated Lactococcus lactis enhanced adherence to human enterocytes through extracellular matrix protein and bacterial aggregation. Furthermore, such piliated L. lactis cells evoked a higher TNF-α response during murine colonization compared with their nonpiliated parent, suggesting that bifidobacterial sortase-dependent pili not only contribute to adherence but also display immunomodulatory activity.

Biological aging alters circadian mechanisms in murine adipose tissue depots

Biological aging alters the metabolism and volume of adipose tissue depots. Recent evidence suggests that circadian mechanisms play a role in promoting adipogenesis, obesity, and lipodystrophy. The current study compared cohorts of younger (5-9 months) and older (24-28 months) C57BL/6 mice as a function of biological age and circadian time. Advanced age significantly reduced the weight of the brown, epididymal, inguinal, and retroperitoneal adipose depots but not total body weight. The older mice reduced their physical activity by >50% and delayed their activity initiation after light offset. The expressed transcriptome in brown and white adipose depots and liver of both cohorts displayed evidence of circadian rhythmicity; however, the oscillating mRNAs differed significantly between age groups and across tissues. The amplitude of Cry1, a component of the negative arm of the circadian apparatus, and downstream regulators such as Rev-erbα were elevated in the older relative to the younger cohorts as a function of circadian time. Overall, transcript levels differed significantly for 557 (inguinal adipose), 1,016 (liver), and 1,021 (brown adipose) expressed sequences between the cohorts as a function of age. These included transcripts encoding proteins within the canonical and non-canonical Wnt pathways. Since the Wnt pathway regulates adipose stem cell differentiation and shares a critical enzyme, glycogen synthase kinase 3β, with the circadian mechanism, the intersection between these two fundamental regulatory mechanisms merits further investigation with respect to biological aging of adipose tissues.

Expression quantitative trait loci analysis identifies associations between genotype and gene expression in human intestine

BACKGROUND & AIMS

Genome-wide association studies have greatly increased our understanding of intestinal disease. However, little is known about how genetic variations result in phenotypic changes. Some polymorphisms have been shown to modulate quantifiable phenotypic traits; these are called quantitative trait loci. Quantitative trait loci that affect levels of gene expression are called expression quantitative trait loci (eQTL), which can provide insight into the biological relevance of data from genome-wide association studies. We performed a comprehensive eQTL scan of intestinal tissue.

METHODS

Total RNA was extracted from ileal biopsy specimens and genomic DNA was obtained from whole-blood samples from the same cohort of individuals. Cis- and trans-eQTL analyses were performed using a custom software pipeline for samples from 173 subjects. The analyses determined the expression levels of 19,047 unique autosomal genes listed in the US National Center for Biotechnology Information database and more than 580,000 variants from the Single Nucleotide Polymorphism database.

RESULTS

The presence of more than 15,000 cis- and trans-eQTL was detected with statistical significance. eQTL associated with the same expression trait were in high linkage disequilibrium. Comparative analysis with previous eQTL studies showed that 30% to 40% of genes identified as eQTL in monocytes, liver tissue, lymphoblastoid cell lines, T cells, and fibroblasts are also eQTL in ileal tissue. Conversely, most of the significant eQTL have not been previously identified and could be tissue specific. These are involved in many cell functions, including division and antigen processing and presentation. Our analysis confirmed that previously published cis-eQTL are single nucleotide polymorphisms associated with inflammatory bowel disease: rs2298428/UBE2L3, rs1050152/SLC22A4, and SLC22A5. We identified many new associations between inflammatory bowel disease susceptibility loci and gene expression.

CONCLUSIONS

eQTL analysis of intestinal tissue supports findings that some eQTL remain stable across cell types, whereas others are specific to the sampled location. Our findings confirm and expand the number of known genotypes associated with expression and could help elucidate mechanisms of intestinal disease.

Characterizing vaccine responses using host genomic and transcriptomic analysis

Vaccines have had a profound influence on human health with no other health intervention rivaling their impact on the morbidity and mortality associated with infectious disease. However, the magnitude and persistence of vaccine immunity varies considerably between individuals, a phenomenon that is not well understood. Recent studies have used contemporary technologies to correlate variation in the genome and transcriptome to immunological measures of vaccine responsiveness. These approaches have provided fresh insight into the intrinsic factors determining the potency and duration of vaccine-induced immunity. The fundamental challenge will be to translate these findings into innovative and pragmatic strategies to develop new and more effective vaccines.

TIMPs of parasitic helminths - a large-scale analysis of high-throughput sequence datasets

BACKGROUND

Tissue inhibitors of metalloproteases (TIMPs) are a multifunctional family of proteins that orchestrate extracellular matrix turnover, tissue remodelling and other cellular processes. In parasitic helminths, such as hookworms, TIMPs have been proposed to play key roles in the host-parasite interplay, including invasion of and establishment in the vertebrate animal hosts. Currently, knowledge of helminth TIMPs is limited to a small number of studies on canine hookworms, whereas no information is available on the occurrence of TIMPs in other parasitic helminths causing neglected diseases.

METHODS

In the present study, we conducted a large-scale investigation of TIMP proteins of a range of neglected human parasites including the hookworm Necator americanus, the roundworm Ascaris suum, the liver flukes Clonorchis sinensis and Opisthorchis viverrini, as well as the schistosome blood flukes. This entailed mining available transcriptomic and/or genomic sequence datasets for the presence of homologues of known TIMPs, predicting secondary structures of defined protein sequences, systematic phylogenetic analyses and assessment of differential expression of genes encoding putative TIMPs in the developmental stages of A. suum, N. americanus and Schistosoma haematobium which infect the mammalian hosts.

RESULTS

A total of 15 protein sequences with high homology to known eukaryotic TIMPs were predicted from the complement of sequence data available for parasitic helminths and subjected to in-depth bioinformatic analyses.

CONCLUSIONS

Supported by the availability of gene manipulation technologies such as RNA interference and/or transgenesis, this work provides a basis for future functional explorations of helminth TIMPs and, in particular, of their role/s in fundamental biological pathways linked to long-term establishment in the vertebrate hosts, with a view towards the development of novel approaches for the control of neglected helminthiases.

Transcriptome profiling of mice testes following low dose irradiation

BACKGROUND

Radiotherapy is used routinely to treat testicular cancer. Testicular cells vary in radio-sensitivity and the aim of this study was to investigate cellular and molecular changes caused by low dose irradiation of mice testis and to identify transcripts from different cell types in the adult testis.

METHODS

Transcriptome profiling was performed on total RNA from testes sampled at various time points (n = 17) after 1 Gy of irradiation. Transcripts displaying large overall expression changes during the time series, but small expression changes between neighbouring time points were selected for further analysis. These transcripts were separated into clusters and their cellular origin was determined. Immunohistochemistry and in silico quantification was further used to study cellular changes post-irradiation (pi).

RESULTS

We identified a subset of transcripts (n = 988) where changes in expression pi can be explained by changes in cellularity. We separated the transcripts into five unique clusters that we associated with spermatogonia, spermatocytes, early spermatids, late spermatids and somatic cells, respectively. Transcripts in the somatic cell cluster showed large changes in expression pi, mainly caused by changes in cellularity. Further investigations revealed that the low dose irradiation seemed to cause Leydig cell hyperplasia, which contributed to the detected expression changes in the somatic cell cluster.

CONCLUSIONS

The five clusters represent gene expression in distinct cell types of the adult testis. We observed large expression changes in the somatic cell profile, which mainly could be attributed to changes in cellularity, but hyperplasia of Leydig cells may also play a role. We speculate that the possible hyperplasia may be caused by lower testosterone production and inadequate inhibin signalling due to missing germ cells.

Multitissue molecular, genomic, and developmental effects of the Deepwater Horizon oil spill on resident Gulf killifish (Fundulus grandis)

The Deepwater Horizon oil rig disaster resulted in crude oil contamination along the Gulf coast in sensitive estuaries. Toxicity from exposure to crude oil can affect populations of fish that live or breed in oiled habitats as seen following the Exxon Valdez oil spill. In an ongoing study of the effects of Deepwater Horizon crude oil on fish, Gulf killifish ( Fundulus grandis ) were collected from an oiled site (Grande Terre, LA) and two reference locations (coastal MS and AL) and monitored for measures of exposure to crude oil. Killifish collected from Grande Terre had divergent gene expression in the liver and gill tissue coincident with the arrival of contaminating oil and up-regulation of cytochrome P4501A (CYP1A) protein in gill, liver, intestine, and head kidney for over one year following peak landfall of oil (August 2011) compared to fish collected from reference sites. Furthermore, laboratory exposures of Gulf killifish embryos to field-collected sediments from Grande Terre and Barataria Bay, LA, also resulted in increased CYP1A and developmental abnormalities when exposed to sediments collected from oiled sites compared to exposure to sediments collected from a reference site. These data are predictive of population-level impacts in fish exposed to sediments from oiled locations along the Gulf of Mexico coast.

Blood meal-induced changes to antennal transcriptome profiles reveal shifts in odor sensitivities in Anopheles gambiae

Olfactory-driven behaviors are central to the lifecycle of the malaria vector mosquito Anopheles gambiae and are initiated by peripheral signaling in the antenna and other olfactory tissues. To continue gaining insight into the relationship between gene expression and olfaction, we have performed cohort comparisons of antennal transcript abundances at five time points after a blood meal, a key event in both reproduction and disease transmission cycles. We found that more than 5,000 transcripts displayed significant abundance differences, many of which were correlated by cluster analysis. Within the chemosensory gene families, we observed a general reduction in the level of chemosensory gene transcripts, although a subset of odorant receptors (AgOrs) was modestly enhanced in post-blood-fed samples. Integration of AgOr transcript abundance data with previously characterized AgOr excitatory odorant response profiles revealed potential changes in antennal odorant receptivity that coincided with the shift from host-seeking to oviposition behaviors in blood-fed female mosquitoes. Behavioral testing of ovipositing females to odorants highlighted by this synthetic analysis identified two unique, unitary oviposition cues for An. gambiae, 2-propylphenol and 4-methylcyclohexanol. We posit that modest, yet cumulative, alterations of AgOr transcript levels modulate peripheral odor coding resulting in biologically relevant behavioral effects. Moreover, these results demonstrate that highly quantitative, RNAseq transcript abundance data can be successfully integrated with functional data to generate testable hypotheses.

Molecular and cellular mechanisms of neutral lipid accumulation in diatom following nitrogen deprivation

BACKGROUND

Nitrogen limitation can induce neutral lipid accumulation in microalgae, as well as inhibiting their growth. Therefore, to obtain cultures with both high biomass and high lipid contents, and explore the lipid accumulation mechanisms, we implemented nitrogen deprivation in a model diatom Phaeodactylum tricornutum at late exponential phase.

RESULTS

Neutral lipid contents per cell subsequently increased 2.4-fold, both the number and total volume of oil bodies increased markedly, and cell density rose slightly. Transcriptional profile analyzed by RNA-Seq showed that expression levels of 1213 genes (including key carbon fixation, TCA cycle, glycerolipid metabolism and nitrogen assimilation genes) increased, with a false discovery rate cut-off of 0.001, under N deprivation. However, most light harvesting complex genes were down-regulated, extensive degradation of chloroplast membranes was observed under an electron microscope, and photosynthetic efficiency declined. Further identification of lipid classes showed that levels of MGDG and DGDG, the main lipid components of chloroplast membranes, dramatically decreased and triacylglycerol (TAG) levels significantly rose, indicating that intracellular membrane remodeling substantially contributed to the neutral lipid accumulation.

CONCLUSIONS

Our findings shed light on the molecular mechanisms of neutral lipid accumulation and the key genes involved in lipid metabolism in diatoms. They also provide indications of possible strategies for improving microalgal biodiesel production.

Sex-biased gene expression in the brown alga Fucus vesiculosus

BACKGROUND

The fucoid brown algae (Heterokontophyta, Phaeophyceae) are increasingly the focus of ecological genetics, biodiversity, biogeography and speciation research. The molecular genetics underlying mating system variation, where repeated dioecious - hermaphrodite switches during evolution are recognized, and the molecular evolution of sex-related genes are key questions currently hampered by a lack of genomic information. We therefore undertook a comparative analysis of male and female reproductive tissue transcriptomes against a vegetative background during natural reproductive cycles in Fucus vesiculosus.

RESULTS

Over 300 k reads were assembled and annotated against public protein databases including a brown alga. Compared with the vegetative tissue, photosynthetic and carbohydrate metabolism pathways were under-expressed, particularly in male tissue, while several pathways involved in genetic information processing and replication were over-expressed. Estimates of sex-biased gene (SBG) expression were higher for male (14% of annotated orthologues) than female tissue (9%) relative to the vegetative background. Mean expression levels and variance were also greater in male- than female-biased genes. Major female-biased genes were carbohydrate-modifying enzymes with likely roles in zygote cell wall biogenesis and/or modification. Male-biased genes reflected distinct sperm development and function, and orthologues for signal perception (a phototropin), transduction (several kinases), and putatively flagella-localized proteins (including candidate gamete-recognition proteins) were uniquely expressed in males. Overall, the results suggest constraint on female-biased genes (possible pleiotropy), and less constrained male-biased genes, mostly associated with sperm-specific functions.

CONCLUSIONS

Our results support the growing contention that males possess a large array of genes regulating male fitness, broadly supporting findings in evolutionarily distant heterogametic animal models. This work identifies an annotated set of F. vesiculosus gene products that potentially regulate sexual reproduction and may contribute to prezygotic isolation, one essential step towards developing tools for a functional understanding of species isolation and differentiation.

Beyond phylotyping: understanding the impact of gut microbiota on host biology

BACKGROUND

Microbial constituents of the gut microbiome interact with each other and the host to alter the luminal environment and impact development, motility, and homeostasis of the gut. Powerful methods are becoming available to investigate connections between the gut microbiome and human health. While high-throughput sequencing of 16S rRNA genes can be used to identify and enumerate microbes in the gut, advances in several techniques (e.g., metagenomics, metatranscriptomics, metabolomics, and metaproteomics) are providing a clearer view as to the specific activities of the microbiota in the context of functional host-microbe interactions. Testing emergent hypotheses regarding microbial effects on host biology, which arise from analyses of 'Big Data' generated from massive parallel high-throughput sequencing technology and spectroscopic techniques, to guide translational research is an important goal for the future. Insights regarding the fundamental operating principles of the gut microbiota should lay the foundation for rational manipulation of the microbiota to promote human health.

PURPOSE

In this review, we provide an overview of current research on the gut microbiome emphasizing current state-of-the-art technologies, approaches, and directions for improvement of our understanding of the impact of the gut microbiota with specific focus on gastrointestinal motility disorders.

Amniotic fluid: the use of high-dimensional biology to understand fetal well-being

Our aim was to review the use of high-dimensional biology techniques, specifically transcriptomics, proteomics, and metabolomics, in amniotic fluid to elucidate the mechanisms behind preterm birth or assessment of fetal development. We performed a comprehensive MEDLINE literature search on the use of transcriptomic, proteomic, and metabolomic technologies for amniotic fluid analysis. All abstracts were reviewed for pertinence to preterm birth or fetal maturation in human subjects. Nineteen articles qualified for inclusion. Most articles described the discovery of biomarker candidates, but few larger, multicenter replication or validation studies have been done. We conclude that the use of high-dimensional systems biology techniques to analyze amniotic fluid has significant potential to elucidate the mechanisms of preterm birth and fetal maturation. However, further multicenter collaborative efforts are needed to replicate and validate candidate biomarkers before they can become useful tools for clinical practice. Ideally, amniotic fluid biomarkers should be translated to a noninvasive test performed in maternal serum or urine.

Decreased hepatotoxic bile acid composition and altered synthesis in progressive human nonalcoholic fatty liver disease

Bile acids (BAs) have many physiological roles and exhibit both toxic and protective influences within the liver. Alterations in the BA profile may be the result of disease induced liver injury. Nonalcoholic fatty liver disease (NAFLD) is a prevalent form of chronic liver disease characterized by the pathophysiological progression from simple steatosis to nonalcoholic steatohepatitis (NASH). The hypothesis of this study is that the 'classical' (neutral) and 'alternative' (acidic) BA synthesis pathways are altered together with hepatic BA composition during progression of human NAFLD. This study employed the use of transcriptomic and metabolomic assays to study the hepatic toxicologic BA profile in progressive human NAFLD. Individual human liver samples diagnosed as normal, steatosis, and NASH were utilized in the assays. The transcriptomic analysis of 70 BA genes revealed an enrichment of downregulated BA metabolism and transcription factor/receptor genes in livers diagnosed as NASH. Increased mRNA expression of BAAT and CYP7B1 was observed in contrast to decreased CYP8B1 expression in NASH samples. The BA metabolomic profile of NASH livers exhibited an increase in taurine together with elevated levels of conjugated BA species, taurocholic acid (TCA) and taurodeoxycholic acid (TDCA). Conversely, cholic acid (CA) and glycodeoxycholic acid (GDCA) were decreased in NASH liver. These findings reveal a potential shift toward the alternative pathway of BA synthesis during NASH, mediated by increased mRNA and protein expression of CYP7B1. Overall, the transcriptomic changes of BA synthesis pathway enzymes together with altered hepatic BA composition signify an attempt by the liver to reduce hepatotoxicity during disease progression to NASH.

A combined blood based gene expression and plasma protein abundance signature for diagnosis of epithelial ovarian cancer--a study of the OVCAD consortium

BACKGROUND

The immune system is a key player in fighting cancer. Thus, we sought to identify a molecular 'immune response signature' indicating the presence of epithelial ovarian cancer (EOC) and to combine this with a serum protein biomarker panel to increase the specificity and sensitivity for earlier detection of EOC.

METHODS

Comparing the expression of 32,000 genes in a leukocytes fraction from 44 EOC patients and 19 controls, three uncorrelated shrunken centroid models were selected, comprised of 7, 14, and 6 genes. A second selection step using RT-qPCR data and significance analysis of microarrays yielded 13 genes (AP2A1, B4GALT1, C1orf63, CCR2, CFP, DIS3, NEAT1, NOXA1, OSM, PAPOLG, PRIC285, ZNF419, and BC037918) which were finally used in 343 samples (90 healthy, six cystadenoma, eight low malignant potential tumor, 19 FIGO I/II, and 220 FIGO III/IV EOC patients). Using new 65 controls and 224 EOC patients (thereof 14 FIGO I/II) the abundances of six plasma proteins (MIF, prolactin, CA125, leptin, osteopondin, and IGF2) was determined and used in combination with the expression values from the 13 genes for diagnosis of EOC.

RESULTS

Combined diagnostic models using either each five gene expression and plasma protein abundance values or 13 gene expression and six plasma protein abundance values can discriminate controls from patients with EOC with Receiver Operator Characteristics Area Under the Curve values of 0.998 and bootstrap .632+ validated classification errors of 3.1% and 2.8%, respectively. The sensitivities were 97.8% and 95.6%, respectively, at a set specificity of 99.6%.

CONCLUSIONS

The combination of gene expression and plasma protein based blood derived biomarkers in one diagnostic model increases the sensitivity and the specificity significantly. Such a diagnostic test may allow earlier diagnosis of epithelial ovarian cancer.

Catalyzing plant science research with RNA-seq

Next generation DNA sequencing technologies are driving increasingly rapid, affordable and high resolution analyses of plant transcriptomes through sequencing of their associated cDNA (complementary DNA) populations; an analytical platform commonly referred to as RNA-sequencing (RNA-seq). Since entering the arena of whole genome profiling technologies only a few years ago, RNA-seq has proven itself to be a powerful tool with a remarkably diverse range of applications, from detailed studies of biological processes at the cell type-specific level, to providing insights into fundamental questions in plant biology on an evolutionary time scale. Applications include generating genomic data for heretofore unsequenced species, thus expanding the boundaries of what had been considered "model organisms," elucidating structural and regulatory gene networks, revealing how plants respond to developmental cues and their environment, allowing a better understanding of the relationships between genes and their products, and uniting the "omics" fields of transcriptomics, proteomics, and metabolomics into a now common systems biology paradigm. We provide an overview of the breadth of such studies and summarize the range of RNA-seq protocols that have been developed to address questions spanning cell type-specific-based transcriptomics, transcript secondary structure and gene mapping.

Performance in omics analyses of blood samples in long-term storage: opportunities for the exploitation of existing biobanks in environmental health research

BACKGROUND

The suitability for omic analysis of biosamples collected in previous decades and currently stored in biobanks is unknown.

OBJECTIVES

We evaluated the influence of handling and storage conditions of blood-derived biosamples on transcriptomic, epigenomic (CpG methylation), plasma metabolomic [UPLC-ToFMS (ultra performance liquid chromatography-time-of-flight mass spectrometry)], and wide-target proteomic profiles.

METHODS

We collected fresh blood samples without RNA preservative in heparin, EDTA, or citrate and held them at room temperature for ≤ 24 hr before fractionating them into buffy coat, erythrocytes, and plasma and freezing the fractions at -80oC or in liquid nitrogen. We developed methodology for isolating RNA from the buffy coats and conducted omic analyses. Finally, we analyzed analogous samples from the EPIC-Italy and Northern Sweden Health and Disease Study biobanks.

RESULTS

Microarray-quality RNA could be isolated from buffy coats (including most biobank samples) that had been frozen within 8 hr of blood collection by thawing the samples in RNA preservative. Different anticoagulants influenced the metabolomic, proteomic, and to a lesser extent transcriptomic profiles. Transcriptomic profiles were most affected by the delay (as little as 2 hr) before blood fractionation, whereas storage temperature had minimal impact. Effects on metabolomic and proteomic profiles were noted in samples processed ≥ 8 hr after collection, but no effects were due to storage temperature. None of the variables examined significantly influenced the epigenomic profiles. No systematic influence of time-in-storage was observed in samples stored over a period of 13-17 years.

CONCLUSIONS

Most samples currently stored in biobanks are amenable to meaningful omics analysis, provided that they satisfy collection and storage criteria defined in this study.

Gene expression signatures of coronary heart disease

OBJECTIVE

To identify transcriptomic biomarkers of coronary heart disease (CHD) in 188 cases with CHD and 188 age- and sex-matched controls who were participants in the Framingham Heart Study.

APPROACH AND RESULTS

A total of 35 genes were differentially expressed in cases with CHD versus controls at false discovery rate<0.5, including GZMB, TMEM56, and GUK1. Cluster analysis revealed 3 gene clusters associated with CHD, 2 linked to increased erythrocyte production and a third to reduced natural killer and T cell activity in cases with CHD. Exon-level results corroborated and extended the gene-level results. Alternative splicing analysis suggested that GUK1 and 38 other genes were differentially spliced in cases with CHD versus controls. Gene Ontology analysis linked ubiquitination and T-cell-related pathways with CHD.

CONCLUSIONS

Two bioinformatically defined groups of genes show consistent associations with CHD. Our findings are consistent with the hypotheses that hematopoesis is upregulated in CHD, possibly reflecting a compensatory mechanism, and that innate immune activity is disrupted in CHD or altered by its treatment. Transcriptomic signatures may be useful in identifying pathways associated with CHD and point toward novel therapeutic targets for its treatment and prevention.

Blood transcriptomics: applications in toxicology

The number of new chemicals that are being synthesized each year has been steadily increasing. While chemicals are of immense benefit to mankind, many of them have a significant negative impact, primarily owing to their inherent chemistry and toxicity, on the environment as well as human health. In addition to chemical exposures, human exposures to numerous non-chemical toxic agents take place in the environment and workplace. Given that human exposure to toxic agents is often unavoidable and many of these agents are found to have detrimental human health effects, it is important to develop strategies to prevent the adverse health effects associated with toxic exposures. Early detection of adverse health effects as well as a clear understanding of the mechanisms, especially at the molecular level, underlying these effects are key elements in preventing the adverse health effects associated with human exposure to toxic agents. Recent developments in genomics, especially transcriptomics, have prompted investigations into this important area of toxicology. Previous studies conducted in our laboratory and elsewhere have demonstrated the potential application of blood gene expression profiling as a sensitive, mechanistically relevant and practical surrogate approach for the early detection of adverse health effects associated with exposure to toxic agents. The advantages of blood gene expression profiling as a surrogate approach to detect early target organ toxicity and the molecular mechanisms underlying the toxicity are illustrated and discussed using recent studies on hepatotoxicity and pulmonary toxicity. Furthermore, the important challenges this emerging field in toxicology faces are presented in this review article.

Current challenges and future potential of tomato breeding using omics approaches

As tomatoes are one of the most important vegetables in the world, improvements in the quality and yield of tomato are strongly required. For this purpose, omics approaches such as metabolomics and transcriptomics are used not only for basic research to understand relationships between important traits and metabolism but also for the development of next generation breeding strategies of tomato plants, because an increase in the knowledge improves the taste and quality, stress resistance and/or potentially health-beneficial metabolites and is connected to improvements in the biochemical composition of tomatoes. Such omics data can be applied to network analyses to potentially reveal unknown cellular regulatory networks in tomato plants. The high-quality tomato genome that was sequenced in 2012 will likely accelerate the application of omics strategies, including next generation sequencing for tomato breeding. In this review, we highlight the current studies of omics network analyses of tomatoes and other plant species, in particular, a gene coexpression network. Key applications of omics approaches are also presented as case examples to improve economically important traits for tomato breeding.

Next-generation sequencing: from understanding biology to personalized medicine

Within just a few years, the new methods for high-throughput next-generation sequencing have generated completely novel insights into the heritability and pathophysiology of human disease. In this review, we wish to highlight the benefits of the current state-of-the-art sequencing technologies for genetic and epigenetic research. We illustrate how these technologies help to constantly improve our understanding of genetic mechanisms in biological systems and summarize the progress made so far. This can be exemplified by the case of heritable heart muscle diseases, so-called cardiomyopathies. Here, next-generation sequencing is able to identify novel disease genes, and first clinical applications demonstrate the successful translation of this technology into personalized patient care.

Comparative transcriptomics reveals different strategies of Trichoderma mycoparasitism

BACKGROUND

Trichoderma is a genus of mycotrophic filamentous fungi (teleomorph Hypocrea) which possess a bright variety of biotrophic and saprotrophic lifestyles. The ability to parasitize and/or kill other fungi (mycoparasitism) is used in plant protection against soil-borne fungal diseases (biological control, or biocontrol). To investigate mechanisms of mycoparasitism, we compared the transcriptional responses of cosmopolitan opportunistic species and powerful biocontrol agents Trichoderma atroviride and T. virens with tropical ecologically restricted species T. reesei during confrontations with a plant pathogenic fungus Rhizoctonia solani.

RESULTS

The three Trichoderma spp. exhibited a strikingly different transcriptomic response already before physical contact with alien hyphae. T. atroviride expressed an array of genes involved in production of secondary metabolites, GH16 ß-glucanases, various proteases and small secreted cysteine rich proteins. T. virens, on the other hand, expressed mainly the genes for biosynthesis of gliotoxin, respective precursors and also glutathione, which is necessary for gliotoxin biosynthesis. In contrast, T. reesei increased the expression of genes encoding cellulases and hemicellulases, and of the genes involved in solute transport. The majority of differentially regulated genes were orthologues present in all three species or both in T. atroviride and T. virens, indicating that the regulation of expression of these genes is different in the three Trichoderma spp. The genes expressed in all three fungi exhibited a nonrandom genomic distribution, indicating a possibility for their regulation via chromatin modification.

CONCLUSION

This genome-wide expression study demonstrates that the initial Trichoderma mycotrophy has differentiated into several alternative ecological strategies ranging from parasitism to predation and saprotrophy. It provides first insights into the mechanisms of interactions between Trichoderma and other fungi that may be exploited for further development of biofungicides.

Biostatistical approaches for the reconstruction of gene co-expression networks based on transcriptomic data

Techniques in molecular biology have permitted the gathering of an extremely large amount of information relating organisms and their genes. The current challenge is assigning a putative function to thousands of genes that have been detected in different organisms. One of the most informative types of genomic data to achieve a better knowledge of protein function is gene expression data. Based on gene expression data and assuming that genes involved in the same function should have a similar or correlated expression pattern, a function can be attributed to those genes with unknown functions when they appear to be linked in a gene co-expression network (GCN). Several tools for the construction of GCNs have been proposed and applied to plant gene expression data. Here, we review recent methodologies used for plant gene expression data and compare the results, advantages and disadvantages in order to help researchers in their choice of a method for the construction of GCNs.

Depression pathogenesis and treatment: what can we learn from blood mRNA expression?

Alterations in several biological systems, including the neuroendocrine and immune systems, have been consistently demonstrated in patients with major depressive disorder. These alterations have been predominantly studied using easily accessible systems such as blood and saliva. In recent years there has been an increasing body of evidence supporting the use of peripheral blood gene expression to investigate the pathogenesis of depression, and to identify relevant biomarkers. In this paper we review the current literature on gene expression alterations in depression, focusing in particular on three important and interlinked biological domains: inflammation, glucocorticoid receptor functionality and neuroplasticity. We also briefly review the few existing transcriptomics studies. Our review summarizes data showing that patients with major depressive disorder exhibit an altered pattern of expression in several genes belonging to these three biological domains when compared with healthy controls. In particular, we show evidence for a pattern of 'state-related' gene expression changes that are normalized either by remission or by antidepressant treatment. Taken together, these findings highlight the use of peripheral blood gene expression as a clinically relevant biomarker approach.

Analysis of omics data with genome-scale models of metabolism

Over the past decade a massive amount of research has been dedicated to generating omics data to gain insight into a variety of biological phenomena, including cancer, obesity, biofuel production, and infection. Although most of these omics data are available publicly, there is a growing concern that much of these data sit in databases without being used or fully analyzed. Statistical inference methods have been widely applied to gain insight into which genes may influence the activities of others in a given omics data set, however, they do not provide information on the underlying mechanisms or whether the interactions are direct or distal. Biochemically, genetically, and genomically consistent knowledge bases are increasingly being used to extract deeper biological knowledge and understanding from these data sets than possible by inferential methods. This improvement is largely due to knowledge bases providing a validated biological context for interpreting the data.

A predictive signature gene set for discriminating active from latent tuberculosis in Warao Amerindian children

BACKGROUND

Tuberculosis (TB) continues to cause a high toll of disease and death among children worldwide. The diagnosis of childhood TB is challenged by the paucibacillary nature of the disease and the difficulties in obtaining specimens. Whereas scientific and clinical research efforts to develop novel diagnostic tools have focused on TB in adults, childhood TB has been relatively neglected. Blood transcriptional profiling has improved our understanding of disease pathogenesis of adult TB and may offer future leads for diagnosis and treatment. No studies applying gene expression profiling of children with TB have been published so far.

RESULTS

We identified a 116-gene signature set that showed an average prediction error of 11% for TB vs. latent TB infection (LTBI) and for TB vs. LTBI vs. healthy controls (HC) in our dataset. A minimal gene set of only 9 genes showed the same prediction error of 11% for TB vs. LTBI in our dataset. Furthermore, this minimal set showed a significant discriminatory value for TB vs. LTBI for all previously published adult studies using whole blood gene expression, with average prediction errors between 17% and 23%. In order to identify a robust representative gene set that would perform well in populations of different genetic backgrounds, we selected ten genes that were highly discriminative between TB, LTBI and HC in all literature datasets as well as in our dataset. Functional annotation of these genes highlights a possible role for genes involved in calcium signaling and calcium metabolism as biomarkers for active TB. These ten genes were validated by quantitative real-time polymerase chain reaction in an additional cohort of 54 Warao Amerindian children with LTBI, HC and non-TB pneumonia. Decision tree analysis indicated that five of the ten genes were sufficient to classify 78% of the TB cases correctly with no LTBI subjects wrongly classified as TB (100% specificity).

CONCLUSIONS

Our data justify the further exploration of our signature set as biomarkers for potential childhood TB diagnosis. We show that, as the identification of different biomarkers in ethnically distinct cohorts is apparent, it is important to cross-validate newly identified markers in all available cohorts.

Evolutionary conservation of the oocyte transcriptome among vertebrates and its implications for understanding human reproductive function

Cross-phylum and cross-species comparative transcriptomic analyses provide an evolutionary perspective on how specific tissues use genomic information. A significant mRNA subset present in the oocytes of most vertebrates is stabilized or stored for post-LH surge use. Since transcription is arrested in the oocyte before ovulation, this RNA is important for completing maturation and sustaining embryo development until zygotic genome activation. We compared the human oocyte transcriptome with an oocyte-enriched subset of mouse, bovine and frog (Xenopus laevis) genes in order to evaluate similarities between species. Graded temperature stringency hybridization on a multi-species oocyte cDNA array was used to measure the similarity of preferentially expressed sequences to the human oocyte library. Identity analysis of 679 human orthologs compared with each identified official gene symbol found in the subtractive (somatic-oocyte) libraries comprising our array revealed that bovine/human similarity was greater than mouse/human or frog/human similarity. However, based on protein sequence, mouse/human similarity was greater than bovine/human similarity. Among the genes over-expressed in oocytes relative to somatic tissue in Xenopus, Mus and Bos, a high level of conservation was found relative to humans, especially for genes involved in early embryonic development.

Systems biological approaches towards understanding cellulase production by Trichoderma reesei

Recent progress and improvement in "-omics" technologies has made it possible to study the physiology of organisms by integrated and genome-wide approaches. This bears the advantage that the global response, rather than isolated pathways and circuits within an organism, can be investigated ("systems biology"). The sequencing of the genome of Trichoderma reesei (teleomorph Hypocrea jecorina), a fungus that serves as a major producer of biomass-degrading enzymes for the use of renewable lignocellulosic material towards production of biofuels and biorefineries, has offered the possibility to study this organism and its enzyme production on a genome wide scale. In this review, I will highlight the use of genomics, transcriptomics, proteomics and metabolomics towards an improved and novel understanding of the biochemical processes that involve in the massive overproduction of secreted proteins.

De novo transcriptome characterization of Vitis vinifera cv. Corvina unveils varietal diversity

BACKGROUND

Plants such as grapevine (Vitis spp.) display significant inter-cultivar genetic and phenotypic variation. The genetic components underlying phenotypic diversity in grapevine must be understood in order to disentangle genetic and environmental factors.

RESULTS

We have shown that cDNA sequencing by RNA-seq is a robust approach for the characterization of varietal diversity between a local grapevine cultivar (Corvina) and the PN40024 reference genome. We detected 15,161 known genes including 9463 with novel splice isoforms, and identified 2321 potentially novel protein-coding genes in non-annotated or unassembled regions of the reference genome. We also discovered 180 apparent private genes in the Corvina genome which were missing from the reference genome.

CONCLUSIONS

The de novo assembly approach allowed a substantial amount of the Corvina transcriptome to be reconstructed, improving known gene annotations by robustly defining gene structures, annotating splice isoforms and detecting genes without annotations. The private genes we discovered are likely to be nonessential but could influence certain cultivar-specific characteristics. Therefore, the application of de novo transcriptome assembly should not be restricted to species lacking a reference genome because it can also improve existing reference genome annotations and identify novel, cultivar-specific genes.

Inter-Laboratory Variability in Array-Based RNA Quantification Methods

Ribonucleic acids (RNA) are hypothesized to have preceded their derivatives, deoxyribonucleic acids (DNA), as the molecular media of genetic information when life emerged on earth. Molecular biologists are accustomed to the dramatic effects a subtle variation in the ribose moiety composition between RNA and DNA can have on the stability of these molecules. While DNA is very stable after extraction from biological samples and subsequent treatment, RNA is notoriously labile. The short half-life property, inherent to RNA, benefits cells that do not need to express their entire repertoire of proteins. The cellular machinery turns off the production of a given protein by shutting down the transcription of its cognate coding gene and by either actively degrading the remaining mRNA or allowing it to decay on its own. The steady-state level of each mRNA in a given cell varies continuously and is specified by changing kinetics of synthesis and degradation. Because it is technically possible to simultaneously measure thousands of nucleic acid molecules, these quantities have been studied by the life sciences community to investigate a range of biological problems. Since the RNA abundance can change according to a wide range of perturbations, this makes it the molecule of choice for exploring biological systems; its instability, on the other hand, could be an underestimated source of technical variability. We found that a large fraction of the RNA abundance originally present in the biological system prior to extraction was masked by the RNA labeling and measurement procedure. The method used to extract RNA molecules from cells and to label them prior to hybridization operations on DNA arrays affects the original distribution of RNA. Only if RNA measurements are performed according to the same procedure can biological information be inferred from the assay read out.

Root-microbe systems: the effect and mode of interaction of Stress Protecting Agent (SPA) Stenotrophomonas rhizophila DSM14405(T.)

Stenotrophomonas rhizophila has great potential for applications in biotechnology and biological control due to its ability to both promote plant growth and protect roots against biotic and a-biotic stresses, yet little is known about the mode of interactions in the root-environment system. We studied mechanisms associated with osmotic stress using transcriptomic and microscopic approaches. In response to salt or root extracts, the transcriptome of S. rhizophila DSM14405(T) changed drastically. We found a notably similar response for several functional gene groups responsible for general stress protection, energy production, and cell motility. However, unique changes in the transcriptome were also observed: the negative regulation of flagella-coding genes together with the up-regulation of the genes responsible for biofilm formation and alginate biosynthesis were identified as a single mechanism of S. rhizophila DSM14405(T) against salt shock. However, production and excretion of glucosylglycerol (GG) were found as a remarkable mechanism for the stress protection of this Stenotrophomonas strain. For S. rhizophila treated with root exudates, the shift from the planktonic lifestyle to a sessile one was measured as expressed in the down-regulation of flagellar-driven motility. These findings fit well with the observed positive regulation of host colonization genes and microscopic images that show different colonization patterns of oilseed rape roots. Spermidine, described as a plant growth regulator, was also newly identified as a protector against stress. Overall, we identified mechanisms of Stenotrophomonas to protect roots against osmotic stress in the environment. In addition to both the changes in life style and energy metabolism, phytohormons, and osmoprotectants were also found to play a key role in stress protection.

OPTIMAS-DW: a comprehensive transcriptomics, metabolomics, ionomics, proteomics and phenomics data resource for maize

BACKGROUND

Maize is a major crop plant, grown for human and animal nutrition, as well as a renewable resource for bioenergy. When looking at the problems of limited fossil fuels, the growth of the world's population or the world's climate change, it is important to find ways to increase the yield and biomass of maize and to study how it reacts to specific abiotic and biotic stress situations. Within the OPTIMAS systems biology project maize plants were grown under a large set of controlled stress conditions, phenotypically characterised and plant material was harvested to analyse the effect of specific environmental conditions or developmental stages. Transcriptomic, metabolomic, ionomic and proteomic parameters were measured from the same plant material allowing the comparison of results across different omics domains. A data warehouse was developed to store experimental data as well as analysis results of the performed experiments.

DESCRIPTION

The OPTIMAS Data Warehouse (OPTIMAS-DW) is a comprehensive data collection for maize and integrates data from different data domains such as transcriptomics, metabolomics, ionomics, proteomics and phenomics. Within the OPTIMAS project, a 44K oligo chip was designed and annotated to describe the functions of the selected unigenes. Several treatment- and plant growth stage experiments were performed and measured data were filled into data templates and imported into the data warehouse by a Java based import tool. A web interface allows users to browse through all stored experiment data in OPTIMAS-DW including all data domains. Furthermore, the user can filter the data to extract information of particular interest. All data can be exported into different file formats for further data analysis and visualisation. The data analysis integrates data from different data domains and enables the user to find answers to different systems biology questions. Finally, maize specific pathway information is provided.

CONCLUSIONS

With OPTIMAS-DW a data warehouse for maize was established, which is able to handle different data domains, comprises several analysis results that will support researchers within their work and supports systems biological research in particular. The system is available at http://www.optimas-bioenergy.org/optimas_dw.

RNA-seq based identification and mutant validation of gene targets related to ethanol resistance in cyanobacterial Synechocystis sp. PCC 6803

BACKGROUND

Fermentation production of biofuel ethanol consumes agricultural crops, which will compete directly with the food supply. As an alternative, photosynthetic cyanobacteria have been proposed as microbial factories to produce ethanol directly from solar energy and CO2. However, the ethanol productivity from photoautotrophic cyanobacteria is still very low, mostly due to the low tolerance of cyanobacterial systems to ethanol stress.

RESULTS

To build a foundation necessary to engineer robust ethanol-producing cyanobacterial hosts, in this study we applied a quantitative transcriptomics approach with a next-generation sequencing technology, combined with quantitative reverse-transcript PCR (RT-PCR) analysis, to reveal the global metabolic responses to ethanol in model cyanobacterial Synechocystis sp. PCC 6803. The results showed that ethanol exposure induced genes involved in common stress responses, transporting and cell envelope modification. In addition, the cells can also utilize enhanced polyhydroxyalkanoates (PHA) accumulation and glyoxalase detoxication pathway as means against ethanol stress. The up-regulation of photosynthesis by ethanol was also further confirmed at transcriptional level. Finally, we used gene knockout strains to validate the potential target genes related to ethanol tolerance.

CONCLUSION

RNA-Seq based global transcriptomic analysis provided a comprehensive view of cellular response to ethanol exposure. The analysis provided a list of gene targets for engineering ethanol tolerance in cyanobacterium Synechocystis.

Regulators of gene expression as biomarkers for prostate cancer

Recent technological advancements in gene expression analysis have led to the discovery of a promising new group of prostate cancer (PCa) biomarkers that have the potential to influence diagnosis and the prediction of disease severity. The accumulation of deleterious changes in gene expression is a fundamental mechanism of prostate carcinogenesis. Aberrant gene expression can arise from changes in epigenetic regulation or mutation in the genome affecting either key regulatory elements or gene sequences themselves. At the epigenetic level, a myriad of abnormal histone modifications and changes in DNA methylation are found in PCa patients. In addition, many mutations in the genome have been associated with higher PCa risk. Finally, over- or underexpression of key genes involved in cell cycle regulation, apoptosis, cell adhesion and regulation of transcription has been observed. An interesting group of biomarkers are emerging from these studies which may prove more predictive than the standard prostate specific antigen (PSA) serum test. In this review, we discuss recent results in the field of gene expression analysis in PCa including the most promising biomarkers in the areas of epigenetics, genomics and the transcriptome, some of which are currently under investigation as clinical tests for early detection and better prognostic prediction of PCa.

Cellular Stress Following Water Deprivation in the Model Legume Lotus japonicus

Drought stress is one of the most important factors in the limitation of plant productivity worldwide. In order to cope with water deprivation, plants have adopted several strategies that produce major changes in gene expression. In this paper, the response to drought stress in the model legume Lotus japonicus was studied using a transcriptomic approach. Drought induced an extensive reprogramming of the transcriptome as related to various aspects of cellular metabolism, including genes involved in photosynthesis, amino acid metabolism and cell wall metabolism, among others. A particular focus was made on the genes involved in the cellular stress response. Key genes involved in the control of the cell cycle, antioxidant defense and stress signaling, were modulated as a consequence of water deprivation. Genes belonging to different families of transcription factors were also highly responsive to stress. Several of them were homologies to known stress-responsive genes from the model plant Arabidopsis thaliana, while some novel transcription factors were peculiar to the L. japonicus drought stress response.

Comparative Genomics and Systems Biology of Malaria Parasites Plasmodium

Malaria is a serious infectious disease that causes over one million deaths yearly. It is caused by a group of protozoan parasites in the genus Plasmodium. No effective vaccine is currently available and the elevated levels of resistance to drugs in use underscore the pressing need for novel antimalarial targets. In this review, we survey omics centered developments in Plasmodium biology, which have set the stage for a quantum leap in our understanding of the fundamental processes of the parasite life cycle and mechanisms of drug resistance and immune evasion.

Grape berry ripening delay induced by a pre-véraison NAA treatment is paralleled by a shift in the expression pattern of auxin- and ethylene-related genes

BACKGROUND

Auxins act as repressors of ripening inception in grape (véraison), while ethylene and abscisic acid (ABA) play a positive role as inducers of the syndrome. Despite the increasing amount of information made available on this topic, the complex network of interactions among these hormones remains elusive. In order to shed light on these aspects, a holistic approach was adopted to evaluate, at the transcriptomic level, the crosstalk between hormones in grape berries, whose ripening progression was delayed by applying naphtalenacetic acid (NAA) one week before véraison.

RESULTS

The NAA treatment caused significant changes in the transcription rate of about 1,500 genes, indicating that auxin delayed grape berry ripening also at the transcriptional level, along with the recovery of a steady state of its intracellular concentration. Hormone indices analysis carried out with the HORMONOMETER tool suggests that biologically active concentrations of auxins were achieved throughout a homeostatic recovery. This occurred within 7 days after the treatment, during which the physiological response was mainly unspecific and due to a likely pharmacological effect of NAA. This hypothesis is strongly supported by the up-regulation of genes involved in auxin conjugation (GH3-like) and action (IAA4- and IAA31-like). A strong antagonistic effect between auxin and ethylene was also observed, along with a substantial 'synergism' between auxins and ABA, although to a lesser extent.

CONCLUSIONS

This study suggests that, in presence of altered levels of auxins, the crosstalk between hormones involves diverse mechanisms, acting at both the hormone response and biosynthesis levels, creating a complex response network.

Transcriptome variation along bud development in grapevine (Vitis vinifera L.)

BACKGROUND

Vegetative buds provide plants in temperate environments the possibility for growth and reproduction when environmental conditions are favorable. In grapevine, crucial developmental events take place within buds during two growing seasons in consecutive years. The first season, the shoot apical meristem within the bud differentiates all the basic elements of the shoot including flowering transition in lateral primordia and development of inflorescence primordia. These events practically end with bud dormancy. The second season, buds resume shoot growth associated to flower formation and development. Gene expression has been previously monitored at specific stages of bud development but has never been followed along the two growing seasons.

RESULTS

Gene expression changes were analyzed along the bud annual cycle at eight different time points. Principal Components Analysis (PCA) revealed that the main factors explaining the global gene expression differences were the processes of bud dormancy and active growth as well as stress responses. Accordingly, non dormant buds showed an enrichment in functional categories typical of actively proliferating and growing cells together with the over abundance of transcripts belonging to stress response pathways. Differential expression analyses performed between consecutive time points indicated that major transcriptional changes were associated to para/endodormancy, endo/ecodormancy and ecodormancy/bud break transitions. Transcripts encoding key regulators of reproductive development were grouped in three major expression clusters corresponding to: (i) transcripts associated to flowering induction, (ii) transcripts associated to flower meristem specification and initiation and (iii) transcripts putatively involved in dormancy. Within this cluster, a MADS-box gene (VvFLC2) and other transcripts with similar expression patterns could participate in dormancy regulation.

CONCLUSIONS

This work provides a global view of major transcriptional changes taking place along bud development in grapevine, highlighting those molecular and biological functions involved in the main events of bud development. As reported in other woody species, the results suggest that genes regulating flowering could also be involved in dormancy regulatory pathways in grapevine.

Ecogenomics of microbial communities in bioremediation of chlorinated contaminated sites

Organohalide compounds such as chloroethenes, chloroethanes, and polychlorinated benzenes are among the most significant pollutants in the world. These compounds are often found in contamination plumes with other pollutants such as solvents, pesticides, and petroleum derivatives. Microbial bioremediation of contaminated sites, has become commonplace whereby key processes involved in bioremediation include anaerobic degradation and transformation of these organohalides by organohalide respiring bacteria and also via hydrolytic, oxygenic, and reductive mechanisms by aerobic bacteria. Microbial ecogenomics has enabled us to not only study the microbiology involved in these complex processes but also develop tools to better monitor and assess these sites during bioremediation. Microbial ecogenomics have capitalized on recent advances in high-throughput and -output genomics technologies in combination with microbial physiology studies to address these complex bioremediation problems at a system level. Advances in environmental metagenomics, transcriptomics, and proteomics have provided insights into key genes and their regulation in the environment. They have also given us clues into microbial community structures, dynamics, and functions at contaminated sites. These techniques have not only aided us in understanding the lifestyles of common organohalide respirers, for example Dehalococcoides, Dehalobacter, and Desulfitobacterium, but also provided insights into novel and yet uncultured microorganisms found in organohalide respiring consortia. In this paper, we look at how ecogenomic studies have aided us to understand the microbial structures and functions in response to environmental stimuli such as the presence of chlorinated pollutants.