Genome sequences of 10 new carnation mottle virus variants.

Here, we report the genome sequences of 10 Carnation mottle virus variants. Six variants originated from a single proprietary carnation cultivar, and four were derived from four different proprietary cultivars. All variants showed nucleotide differences, but the last four did not show any variation at the amino acid level.

Genome Sequence of a New Carnation Small Viroid-Like RNA, CarSV-1.

Here, we report the genome sequence of a new circular viroid-like RNA (CarSV-1) derived from Dianthus caryophyllus (carnation) leaves. The CarSV-1 genome has notable sequence similarity (62%) to the well-studied CarSV viroid-like RNA and comprises the complete hammerhead consensus sequences involved in self-cleavage. CarSV-1 co-occurs with carnation viruses, such as CarMV.

miRNA repertoire and host immune factor regulation upon avian coronavirus infection in eggs.

Avian infectious bronchitis virus (IBV) is a coronavirus with great economic impact on the poultry industry, causing an acute and highly contagious disease in chickens that primarily affects the respiratory and reproductive systems. The cellular regulation of IBV pathogenesis and the host immune responses involved remain to be fully elucidated. MicroRNAs (miRNAs) have emerged as a class of crucial regulators of numerous cellular processes, including responses to viral infections. Here, we employed a high-throughput sequencing approach to analyze the miRNA composition of the spleen and the lungs of chicken embryos upon IBV infection. Compared to healthy chicken embryos, 13 and six miRNAs were upregulated in the spleen and the lungs, respectively, all predicted to influence viral transcription, cytokine production, and lymphocyte functioning. Subsequent downregulation of NFATC3, NFAT5, SPPL3, and TGFB2 genes in particular was observed only in the spleen, demonstrating the biological functionality of the miRNAs in this lymphoid organ. This is the first study that describes the modulation of miRNAs and the related host immune factors by IBV in chicken embryos. Our data provide novel insight into complex virus-host interactions and specifically highlight components that could affect the host's immune response to IBV infection.

Identifying small RNAs derived from maternal- and somatic-type rRNAs in zebrafish development.

rRNAs are non-coding RNAs present in all prokaryotes and eukaryotes. In eukaryotes there are four rRNAs: 18S, 5.8S, 28S, originating from a common precursor (45S), and 5S. We have recently discovered the existence of two distinct developmental types of rRNA: a maternal-type, present in eggs and a somatic-type, expressed in adult tissues. Lately, next-generation sequencing has allowed the discovery of new small-RNAs deriving from longer non-coding RNAs, including small-RNAs from rRNAs (srRNAs). Here, we systemically investigated srRNAs of maternal- or somatic-type 18S, 5.8S, 28S, with small-RNAseq from many zebrafish developmental stages. We identified new srRNAs for each rRNA. For 5.8S, we found srRNA consisting of the 5' or 3' halves, with only the latter having different sequence for the maternal- and somatic-types. For 18S, we discovered 21 nt srRNA from the 5' end of the 18S rRNA with a striking resemblance to microRNAs; as it is likely processed from a stem-loop precursor and present in human and mouse Argonaute-complexed small-RNA. For 28S, an abundant 80 nt srRNA from the 3' end of the 28S rRNA was found. The expression levels during embryogenesis of these srRNA indicate they are not generated from rRNA degradation and might have a role in the zebrafish development.

Expression of distinct maternal and somatic 5.8S, 18S, and 28S rRNA types during zebrafish development.

There is mounting evidence that the ribosome is not a static translation machinery, but a cell-specific, adaptive system. Ribosomal variations have mostly been studied at the protein level, even though the essential transcriptional functions are primarily performed by rRNAs. At the RNA level, oocyte-specific 5S rRNAs are long known for Xenopus. Recently, we described for zebrafish a similar system in which the sole maternal-type 5S rRNA present in eggs is replaced completely during embryonic development by a somatic-type. Here, we report the discovery of an analogous system for the 45S rDNA elements: 5.8S, 18S, and 28S. The maternal-type 5.8S, 18S, and 28S rRNA sequences differ substantially from those of the somatic-type, plus the maternal-type rRNAs are also replaced by the somatic-type rRNAs during embryogenesis. We discuss the structural and functional implications of the observed sequence differences with respect to the translational functions of the 5.8S, 18S, and 28S rRNA elements. Finally, in silico evidence suggests that expansion segments (ES) in 18S rRNA, previously implicated in ribosome-mRNA interaction, may have a preference for interacting with specific mRNA genes. Taken together, our findings indicate that two distinct types of ribosomes exist in zebrafish during development, each likely conducting the translation machinery in a unique way.

Transcriptome dynamics in early zebrafish embryogenesis determined by high-resolution time course analysis of 180 successive, individual zebrafish embryos.

BACKGROUND: Recently, much progress has been made in the field of gene-expression in early embryogenesis. However, the dynamic behaviour of transcriptomes in individual embryos has hardly been studied yet and the time points at which pools of embryos are collected are usually still quite far apart. Here, we present a high-resolution gene-expression time series with 180 individual zebrafish embryos, obtained from nine different spawns, developmentally ordered and profiled from late blastula to mid-gastrula stage. On average one embryo per minute was analysed. The focus was on identification and description of the transcriptome dynamics of the expressed genes in this embryonic stage, rather than to biologically interpret profiles in cellular processes and pathways. RESULTS: In the late blastula to mid-gastrula stage, we found 6,734 genes being expressed with low variability and rather gradual changes. Ten types of dynamic behaviour were defined, such as genes with continuously increasing or decreasing expression, and all expressed genes were grouped into these types. Also, the exact expression starting and stopping points of several hundred genes during this developmental period could be pinpointed. Although the resolution of the experiment was so high, that we were able to clearly identify four known oscillating genes, no genes were observed with a peaking expression. Additionally, several genes showed expression at two or three distinct levels that strongly related to the spawn an embryo originated from. CONCLUSION: Our unique experimental set-up of whole-transcriptome analysis of 180 individual embryos, provided an unparalleled in-depth insight into the dynamics of early zebrafish embryogenesis. The existence of a tightly regulated embryonic transcriptome program, even between individuals from different spawns is shown. We have made the expression profile of all genes available for domain experts. The fact that we were able to separate the different spawns by their gene-expression variance over all expressed genes, underlines the importance of spawn specificity, as well as the unexpectedly tight gene-expression regulation in early zebrafish embryogenesis.

Genetic characterization of T-DNA insertions in the genome of the Arabidopsis thaliana sumo1/2 knock-down line.

Sumoylation is an essential post-translational modification in Arabidopsis thaliana, which entails the conjugation of the SUMO protein onto lysine residues in target proteins. In Arabidopsis, 2 closely related genes, SUMO1 and SUMO2, act redundantly and are in combination essential for plant development, i.e. the combined loss of SUMO1 and SUMO2 results in embryo-lethality. To circumvent this lethality, SUMO2 was previously knocked down in a sumo1 knockout background by expressing an artificial microRNA that targets SUMO2 (amiR-SUMO2). This sumo1/2KD line with low SUMO2 levels represents a valuable genetics tool to investigate SUMO function in planta. Here, we re-sequenced the whole-genome of this sumo1/2KD line and identified 2 amiR-SUMO2 insertions in this line, which were confirmed by PCR-genotyping. Identification of these 2 insertions enables genetics with this tool.

Linking maternal and somatic 5S rRNA types with different sequence-specific non-LTR retrotransposons.

5S rRNA is a ribosomal core component, transcribed from many gene copies organized in genomic repeats. Some eukaryotic species have two 5S rRNA types defined by their predominant expression in oogenesis or adult tissue. Our next-generation sequencing study on zebrafish egg, embryo, and adult tissue identified maternal-type 5S rRNA that is exclusively accumulated during oogenesis, replaced throughout the embryogenesis by a somatic-type, and thus virtually absent in adult somatic tissue. The maternal-type 5S rDNA contains several thousands of gene copies on chromosome 4 in tandem repeats with small intergenic regions, whereas the somatic-type is present in only 12 gene copies on chromosome 18 with large intergenic regions. The nine-nucleotide variation between the two 5S rRNA types likely affects TFIII binding and riboprotein L5 binding, probably leading to storage of maternal-type rRNA. Remarkably, these sequence differences are located exactly at the sequence-specific target site for genome integration by the 5S rRNA-specific Mutsu retrotransposon family. Thus, we could define maternal- and somatic-type MutsuDr subfamilies. Furthermore, we identified four additional maternal-type and two new somatic-type MutsuDr subfamilies, each with their own target sequence. This target-site specificity, frequently intact maternal-type retrotransposon elements, plus specific presence of Mutsu retrotransposon RNA and piRNA in egg and adult tissue, suggest an involvement of retrotransposons in achieving the differential copy number of the two types of 5S rDNA loci.

Transcriptome data on maternal RNA of 24 individual zebrafish eggs from five sibling mothers.

Maternal mRNA that is present in the mature oocyte plays an important role in the proper development of the early embryo. To elucidate the role of the maternal transcriptome we recently reported a microarray study on individual zebrafish eggs from five different clutches from sibling mothers and showed differences in maternal RNA abundance between and within clutches, "Mother-specific signature in the maternal transcriptome composition of mature, unfertilized Eggs" [1]. Here we provide in detail the applied preprocessing method as well as the R-code to identify expressed and non-expressed genes in the associated transcriptome dataset. Additionally, we provide a website that allows a researcher to search for the expression of their gene of interest in this experiment.

Confounding Factors in the Transcriptome Analysis of an In-Vivo Exposure Experiment.

CONFOUNDING FACTORS: In transcriptomics experimentation, confounding factors frequently exist alongside the intended experimental factors and can severely influence the outcome of a transcriptome analysis. Confounding factors are regularly discussed in methodological literature, but their actual, practical impact on the outcome and interpretation of transcriptomics experiments is, to our knowledge, not documented. For instance, in-vivo experimental factors; like Individual, Sample-Composition and Time-of-Day are potentially formidable confounding factors. To study these confounding factors, we designed an extensive in-vivo transcriptome experiment (n = 264) with UVR exposure of murine skin containing six consecutive samples from each individual mouse (n = 64). ANALYSIS APPROACH: Evaluation of the confounding factors: Sample-Composition, Time-of-Day, Handling-Stress, and Individual-Mouse resulted in the identification of many genes that were affected by them. These genes sometimes showed over 30-fold expression differences. The most prominent confounding factor was Sample-Composition caused by mouse-dependent skin composition differences, sampling variation and/or influx/efflux of mobile cells. Although we can only evaluate these effects for known cell type specifically expressed genes in our complex heterogeneous samples, it is clear that the observed variations also affect the cumulative expression levels of many other non-cell-type-specific genes. ANOVA: ANOVA analysis can only attempt to neutralize the effects of the well-defined confounding factors, such as Individual-Mouse, on the experimental factors UV-Dose and Recovery-Time. Also, by definition, ANOVA only yields reproducible gene-expression differences, but we found that these differences were very small compared to the fold changes induced by the confounding factors, questioning the biological relevance of these ANOVA-detected differences. Furthermore, it turned out that many of the differentially expressed genes found by ANOVA were also present in the gene clusters associated with the confounding factors. CONCLUSION: Hence our overall conclusion is that confounding factors have a major impact on the outcome of in-vivo transcriptomics experiments. Thus the set-up, analysis, and interpretation of such experiments should be approached with the utmost prudence.

Mother-Specific Signature in the Maternal Transcriptome Composition of Mature, Unfertilized Zebrafish Eggs.

Maternal mRNA present in mature oocytes plays an important role in the proper development of the early embryo. As the composition of the maternal transcriptome in general has been studied with pooled mature eggs, potential differences between individual eggs are unknown. Here we present a transcriptome study on individual zebrafish eggs from clutches of five mothers in which we focus on the differences in maternal mRNA abundance per gene between and within clutches. To minimize technical interference, we used mature, unfertilized eggs from siblings. About half of the number of analyzed genes was found to be expressed as maternal RNA. The expressed and non-expressed genes showed that maternal mRNA accumulation is a non-random process, as it is related to specific biological pathways and processes relevant in early embryogenesis. Moreover, it turned out that overall the composition of the maternal transcriptome is tightly regulated as about half of the expressed genes display a less than twofold expression range between the observed minimum and maximum expression values of a gene in the experiment. Even more, the maximum gene-expression difference within clutches is for 88% of the expressed genes lower than twofold. This means that expression differences observed in maternally expressed genes are primarily caused by differences between mothers, with only limited variability between eggs from the same mother. This was underlined by the fact that 99% of the expressed genes were found to be differentially expressed between any of the mothers in an ANOVA test. Furthermore, linking chromosome location, transcription factor binding sites, and miRNA target sites of the genes in clusters of distinct and unique mother-specific gene-expression, suggest biological relevance of the mother-specific signatures in the maternal transcriptome composition. Altogether, the maternal transcriptome composition of mature zebrafish oocytes seems to be tightly regulated with a distinct mother-specific signature.

WRKY40 and WRKY6 act downstream of the green leaf volatile E-2-hexenal in Arabidopsis.

Plants are known to be responsive to volatiles, but knowledge about the molecular players involved in transducing their perception remains scarce. We study the response of Arabidopsis thaliana to E-2-hexenal, one of the green leaf volatiles (GLV) that is produced upon wounding, herbivory or infection with pathogens. We have taken a transcriptomics approach to identify genes that are induced by E-2-hexenal, but not by defence hormones or other GLVs. Furthermore, by studying the promoters of early E-2-hexenal-induced genes we determined that the only statistically enriched cis-element was the W-box motif. Since members of the plant-specific family of WRKY transcription factors act in trans on this cis-element, we focused on WRKY6, 40 and 53 that were most strongly induced by E-2-hexenal. Root elongation of Arabidopsis seedlings of the wrky40 wrky6 double mutant was much less inhibited than in wt plants, similar to the E-2-hexenal-responsive mutant her1, which is perturbed in γ-amino butyric acid (GABA) metabolism. The induction of several of the E-2-hexenal-specific genes was much higher in the wrky40, wrky6 or wrky40 wrky6 mutants, including GAD4, a glutamate decarboxylase that catalyzes the formation of GABA from glutamate. In conclusion, WRKY6 and 40 seem to act as important players transducing E-2-hexenal perception.

Valuable lessons-learned in transcriptomics experimentation.

We have collected several valuable lessons that will help improve transcriptomics experimentation. These lessons relate to experiment design, execution, and analysis. The cautions, but also the pointers, may help biologists avoid common pitfalls in transcriptomics experimentation and achieve better results with their transcriptome studies.

Improving small RNA-seq by using a synthetic spike-in set for size-range quality control together with a set for data normalization.

There is an increasing interest in complementing RNA-seq experiments with small-RNA (sRNA) expression data to obtain a comprehensive view of a transcriptome. Currently, two main experimental challenges concerning sRNA-seq exist: how to check the size distribution of isolated sRNAs, given the sensitive size-selection steps in the protocol; and how to normalize data between samples, given the low complexity of sRNA types. We here present two separate sets of synthetic RNA spike-ins for monitoring size-selection and for performing data normalization in sRNA-seq. The size-range quality control (SRQC) spike-in set, consisting of 11 oligoribonucleotides (10-70 nucleotides), was tested by intentionally altering the size-selection protocol and verified via several comparative experiments. We demonstrate that the SRQC set is useful to reproducibly track down biases in the size-selection in sRNA-seq. The external reference for data-normalization (ERDN) spike-in set, consisting of 19 oligoribonucleotides, was developed for sample-to-sample normalization in differential-expression analysis of sRNA-seq data. Testing and applying the ERDN set showed that it can reproducibly detect differential expression over a dynamic range of 2(18). Hence, biological variation in sRNA composition and content between samples is preserved while technical variation is effectively minimized. Together, both spike-in sets can significantly improve the technical reproducibility of sRNA-seq.

A novel toxicogenomics-based approach to categorize (non-)genotoxic carcinogens.

Alternative methods to detect non-genotoxic carcinogens are urgently needed, as this class of carcinogens goes undetected in the current testing strategy for carcinogenicity under REACH. A complicating factor is that non-genotoxic carcinogens act through several distinctive modes of action, which makes prediction of their carcinogenic property difficult. We have recently demonstrated that gene expression profiling in primary mouse hepatocytes is a useful approach to categorize non-genotoxic carcinogens according to their modes of action. In the current study, we improved the methods used for analysis and added mouse embryonic stem cells as a second in vitro test system, because of their features complementary to hepatocytes. Our approach involved an unsupervised analysis based on the 30 most significantly up- and down-regulated genes per chemical. Mouse embryonic stem cells and primary mouse hepatocytes were exposed to a selected set of chemicals and subsequently subjected to gene expression profiling. We focused on non-genotoxic carcinogens, but also included genotoxic carcinogens and non-carcinogens to test the robustness of this approach. Application of the optimized comparison approach resulted in improved categorization of non-genotoxic carcinogens. Mouse embryonic stem cells were a useful addition, especially for genotoxic substances, but also for detection of non-genotoxic carcinogens that went undetected by primary hepatocytes. The approach presented here is an important step forward to categorize chemicals, especially those that are carcinogenic.

A day and night difference in the response of the hepatic transcriptome to cyclophosphamide treatment.

Application of omics-based technologies is a widely used approach in research aiming to improve testing strategies for human health risk assessment. In most of these studies, however, temporal variations in gene expression caused by the circadian clock are a commonly neglected pitfall. In the present study, we investigated the impact of the circadian clock on the response of the hepatic transcriptome after exposure of mice to the chemotherapeutic agent cyclophosphamide (CP). Analysis of the data without considering clock progression revealed common responses in terms of regulated pathways between light and dark phase exposure, including DNA damage, oxidative stress, and a general immune response. The overall response, however, was stronger in mice exposed during the day. Use of time-matched controls, thereby eliminating non-CP-responsive circadian clock-controlled genes, showed that this difference in response was actually even more pronounced: CP-related responses were only identified in mice exposed during the day. Only minor differences were found in acute toxicity pathways, namely lymphocyte counts and kidney weights, indicating that gene expression is subject to time of day effects. This study is the first to highlight the impact of the circadian clock on the identification of toxic responses by omics approaches.

Staphylococcus aureus ST398 gene expression profiling during ex vivo colonization of porcine nasal epithelium.

BACKGROUND: Staphylococcus aureus is a common human and animal opportunistic pathogen. In humans nasal carriage of S. aureus is a risk factor for various infections. Methicillin-resistant S. aureus ST398 is highly prevalent in pigs in Europe and North America. The mechanism of successful pig colonization by MRSA ST398 is poorly understood. Previously, we developed a nasal colonization model of porcine nasal mucosa explants to identify molecular traits involved in nasal MRSA colonization of pigs. RESULTS: We report the analysis of changes in the transcription of MRSA ST398 strain S0462 during colonization on the explant epithelium. Major regulated genes were encoding metabolic processes and regulation of these genes may represent metabolic adaptation to nasal mucosa explants. Colonization was not accompanied by significant changes in transcripts of the main virulence associated genes or known human colonization factors. Here, we documented regulation of two genes which have potential influence on S. aureus colonization; cysteine extracellular proteinase (scpA) and von Willebrand factor-binding protein (vWbp, encoded on SaPIbov5). Colonization with isogenic-deletion strains (Δvwbp and ΔscpA) did not alter the ex vivo nasal S. aureus colonization compared to wild type. CONCLUSIONS: Our results suggest that nasal colonization with MRSA ST398 is a complex event that is accompanied with changes in bacterial gene expression regulation and metabolic adaptation.

A range finding protocol to support design for transcriptomics experimentation: examples of in-vitro and in-vivo murine UV exposure.

In transcriptomics research, design for experimentation by carefully considering biological, technological, practical and statistical aspects is very important, because the experimental design space is essentially limitless. Usually, the ranges of variable biological parameters of the design space are based on common practices and in turn on phenotypic endpoints. However, specific sub-cellular processes might only be partially reflected by phenotypic endpoints or outside the associated parameter range. Here, we provide a generic protocol for range finding in design for transcriptomics experimentation based on small-scale gene-expression experiments to help in the search for the right location in the design space by analyzing the activity of already known genes of relevant molecular mechanisms. Two examples illustrate the applicability: in-vitro UV-C exposure of mouse embryonic fibroblasts and in-vivo UV-B exposure of mouse skin. Our pragmatic approach is based on: framing a specific biological question and associated gene-set, performing a wide-ranged experiment without replication, eliminating potentially non-relevant genes, and determining the experimental 'sweet spot' by gene-set enrichment plus dose-response correlation analysis. Examination of many cellular processes that are related to UV response, such as DNA repair and cell-cycle arrest, revealed that basically each cellular (sub-) process is active at its own specific spot(s) in the experimental design space. Hence, the use of range finding, based on an affordable protocol like this, enables researchers to conveniently identify the 'sweet spot' for their cellular process of interest in an experimental design space and might have far-reaching implications for experimental standardization.

Absence/presence calling in microarray-based CGH experiments with non-model organisms.

Structural variations in genomes are commonly studied by (micro)array-based comparative genomic hybridization. The data analysis methods to infer copy number variation in model organisms (human, mouse) are established. In principle, the procedures are based on signal ratios between test and reference samples and the order of the probe targets in the genome. These procedures are less applicable to experiments with non-model organisms, which frequently comprise non-sequenced genomes with an unknown order of probe targets. We therefore present an additional analysis approach, which does not depend on the structural information of a reference genome, and quantifies the presence or absence of a probe target in an unknown genome. The principle is that intensity values of target probes are compared with the intensities of negative-control probes and positive-control probes from a control hybridization, to determine if a probe target is absent or present. In a test, analyzing the genome content of a known bacterial strain: Staphylococcus aureus MRSA252, this approach proved to be successful, demonstrated by receiver operating characteristic area under the curve values larger than 0.9995. We show its usability in various applications, such as comparing genome content and validating next-generation sequencing reads from eukaryotic non-model organisms.

In vivo murine hepatic microRNA and mRNA expression signatures predicting the (non-)genotoxic carcinogenic potential of chemicals.

There is a high need to improve the assessment of, especially non-genotoxic, carcinogenic features of chemicals. We therefore explored a toxicogenomics-based approach using genome-wide microRNA and mRNA expression profiles upon short-term exposure in mice. For this, wild-type mice were exposed for seven days to three different classes of chemicals, i.e., four genotoxic carcinogens (GTXC), seven non-genotoxic carcinogens (NGTXC), and five toxic non-carcinogens. Hepatic expression patterns of mRNA and microRNA transcripts were determined after exposure and used to assess the discriminative power of the in vivo transcriptome for GTXC and NGTXC. A final classifier set, discriminative for GTXC and NGTXC, was generated from the transcriptomic data using a tiered approach. This appeared to be a valid approach, since the predictive power of the final classifier set in three different classifier algorithms was very high for the original training set of chemicals. Subsequent validation in an additional set of chemicals revealed that the predictive power for GTXC remained high, in contrast to NGTXC, which appeared to be more troublesome. Our study demonstrated that the in vivo microRNA-ome has less discriminative power to correctly identify (non-)genotoxic carcinogen classes. The results generally indicate that single mRNA transcripts do have the potential to be applied in risk assessment, but that additional (genomic) strategies are necessary to correctly predict the non-genotoxic carcinogenic potential of a chemical.