Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD.

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

Unknown Metabolite Identification Using Machine Learning Collision Cross-Section Prediction and Tandem Mass Spectrometry.

Ion mobility (IM) spectrometry provides semiorthogonal data to mass spectrometry (MS), showing promise for identifying unknown metabolites in complex non-targeted metabolomics data sets. While current literature has showcased IM-MS for identifying unknowns under near ideal circumstances, less work has been conducted to evaluate the performance of this approach in metabolomics studies involving highly complex samples with difficult matrices. Here, we present a workflow incorporating de novo molecular formula annotation and MS/MS structure elucidation using SIRIUS 4 with experimental IM collision cross-section (CCS) measurements and machine learning CCS predictions to identify differential unknown metabolites in mutant strains of Caenorhabditis elegans. For many of those ion features, this workflow enabled the successful filtering of candidate structures generated by in silico MS/MS predictions, though in some cases, annotations were challenged by significant hurdles in instrumentation performance and data analysis. While for 37% of differential features we were able to successfully collect both MS/MS and CCS data, fewer than half of these features benefited from a reduction in the number of possible candidate structures using CCS filtering due to poor matching of the machine learning training sets, limited accuracy of experimental and predicted CCS values, and lack of candidate structures resulting from the MS/MS data. When using a CCS error cutoff of ±3%, on average, 28% of candidate structures could be successfully filtered. Herein, we identify and describe the bottlenecks and limitations associated with the identification of unknowns in non-targeted metabolomics using IM-MS to focus and provide insights into areas requiring further improvement.

An anchored experimental design and meta-analysis approach to address batch effects in large-scale metabolomics.

Untargeted metabolomics studies are unbiased but identifying the same feature across studies is complicated by environmental variation, batch effects, and instrument variability. Ideally, several studies that assay the same set of metabolic features would be used to select recurring features to pursue for identification. Here, we developed an anchored experimental design. This generalizable approach enabled us to integrate three genetic studies consisting of 14 test strains of Caenorhabditis elegans prior to the compound identification process. An anchor strain, PD1074, was included in every sample collection, resulting in a large set of biological replicates of a genetically identical strain that anchored each study. This enables us to estimate treatment effects within each batch and apply straightforward meta-analytic approaches to combine treatment effects across batches without the need for estimation of batch effects and complex normalization strategies. We collected 104 test samples for three genetic studies across six batches to produce five analytical datasets from two complementary technologies commonly used in untargeted metabolomics. Here, we use the model system C. elegans to demonstrate that an augmented design combined with experimental blocks and other metabolomic QC approaches can be used to anchor studies and enable comparisons of stable spectral features across time without the need for compound identification. This approach is generalizable to systems where the same genotype can be assayed in multiple environments and provides biologically relevant features for downstream compound identification efforts. All methods are included in the newest release of the publicly available SECIMTools based on the open-source Galaxy platform.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines.

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, and NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the nested association mapping pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long-term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

GAIT-GM integrative cross-omics analyses reveal cholinergic defects in a C. elegans model of Parkinson's disease.

Parkinson's disease (PD) is a disabling neurodegenerative disorder in which multiple cell types, including dopaminergic and cholinergic neurons, are affected. The mechanisms of neurodegeneration in PD are not fully understood, limiting the development of therapies directed at disease-relevant molecular targets. C. elegans is a genetically tractable model system that can be used to disentangle disease mechanisms in complex diseases such as PD. Such mechanisms can be studied combining high-throughput molecular profiling technologies such as transcriptomics and metabolomics. However, the integrative analysis of multi-omics data in order to unravel disease mechanisms is a challenging task without advanced bioinformatics training. Galaxy, a widely-used resource for enabling bioinformatics analysis by the broad scientific community, has poor representation of multi-omics integration pipelines. We present the integrative analysis of gene expression and metabolite levels of a C. elegans PD model using GAIT-GM, a new Galaxy tool for multi-omics data analysis. Using GAIT-GM, we discovered an association between branched-chain amino acid metabolism and cholinergic neurons in the C. elegans PD model. An independent follow-up experiment uncovered cholinergic neurodegeneration in the C. elegans model that is consistent with cholinergic cell loss observed in PD. GAIT-GM is an easy to use Galaxy-based tool for generating novel testable hypotheses of disease mechanisms involving gene-metabolite relationships.

Long-Term Metabolomics Reference Material.

The use of quality control samples in metabolomics ensures data quality, reproducibility, and comparability between studies, analytical platforms, and laboratories. Long-term, stable, and sustainable reference materials (RMs) are a critical component of the quality assurance/quality control (QA/QC) system; however, the limited selection of currently available matrix-matched RMs reduces their applicability for widespread use. To produce an RM in any context, for any matrix that is robust to changes over the course of time, we developed iterative batch averaging method (IBAT). To illustrate this method, we generated 11 independently grown Escherichia coli batches and made an RM over the course of 10 IBAT iterations. We measured the variance of these materials by nuclear magnetic resonance (NMR) and showed that IBAT produces a stable and sustainable RM over time. This E. coli RM was then used as a food source to produce a Caenorhabditis elegans RM for a metabolomics experiment. The metabolite extraction of this material, alongside 41 independently grown individual C. elegans samples of the same genotype, allowed us to estimate the proportion of sample variation in preanalytical steps. From the NMR data, we found that 40% of the metabolite variance is due to the metabolite extraction process and analysis and 60% is due to sample-to-sample variance. The availability of RMs in untargeted metabolomics is one of the predominant needs of the metabolomics community that reach beyond quality control practices. IBAT addresses this need by facilitating the production of biologically relevant RMs and increasing their widespread use.

Testcrosses are an efficient strategy for identifying cis-regulatory variation: Bayesian analysis of allele-specific expression (BayesASE).

Allelic imbalance (AI) occurs when alleles in a diploid individual are differentially expressed and indicates cis acting regulatory variation. What is the distribution of allelic effects in a natural population? Are all alleles the same? Are all alleles distinct? The approach described applies to any technology generating allele-specific sequence counts, for example for chromatin accessibility and can be applied generally including to comparisons between tissues or environments for the same genotype. Tests of allelic effect are generally performed by crossing individuals and comparing expression between alleles directly in the F1. However, a crossing scheme that compares alleles pairwise is a prohibitive cost for more than a handful of alleles as the number of crosses is at least (n2-n)/2 where n is the number of alleles. We show here that a testcross design followed by a hypothesis test of AI between testcrosses can be used to infer differences between nontester alleles, allowing n alleles to be compared with n crosses. Using a mouse data set where both testcrosses and direct comparisons have been performed, we show that the predicted differences between nontester alleles are validated at levels of over 90% when a parent-of-origin effect is present and of 60%-80% overall. Power considerations for a testcross, are similar to those in a reciprocal cross. In all applications, the testing for AI involves several complex bioinformatics steps. BayesASE is a complete bioinformatics pipeline that incorporates state-of-the-art error reduction techniques and a flexible Bayesian approach to estimating AI and formally comparing levels of AI between conditions. The modular structure of BayesASE has been packaged in Galaxy, made available in Nextflow and as a collection of scripts for the SLURM workload manager on github (https://github.com/McIntyre-Lab/BayesASE).

Risk Factors for Recurrent Staphylococcus aureus Bacteremia.

BACKGROUND: To understand the clinical, bacterial, and host characteristics associated with recurrent Staphylococcus aureus bacteremia (R-SAB), patients with R-SAB were compared to contemporaneous patients with a single episode of SAB (S-SAB). METHODS: All SAB isolates underwent spa genotyping. All isolates from R-SAB patients underwent pulsed-field gel electrophoresis (PFGE). PFGE-indistinguishable pairs from 40 patients underwent whole genome sequencing (WGS). Acute phase plasma from R-SAB and S-SAB patients was matched 1:1 for age, race, sex, and bacterial genotype, and underwent cytokine quantification using 25-analyte multiplex bead array. RESULTS: R-SAB occurred in 69 (9.1%) of the 756 study patients. Of the 69 patients, 30 experienced relapse (43.5%) and 39 reinfection (56.5%). Age, race, hemodialysis dependence, presence of foreign body, methicillin-resistant Staphyloccus aureus, and persistent bacteremia were individually associated with likelihood of recurrence. Multivariate risk modeling revealed that black hemodialysis patients were nearly 2 times more likely (odds ratio [OR] = 9.652 [95% confidence interval [CI], 5.402-17.418]) than white hemodialysis patients (OR = 4.53 [95% CI, 1.696-10.879]) to experience R-SAB. WGS confirmed PFGE interpretations in all cases. Median RANTES (regulated on activation, normal T cell expressed and secreted) levels in acute phase plasma from the initial episode of SAB were higher in R-SAB than in matched S-SAB controls (P = .0053, false discovery rate < 0.10). CONCLUSION: This study identified several risk factors for R-SAB. The largest risk for R-SAB is among black hemodialysis patients. Higher RANTES levels in R-SAB compared to matched controls warrants further study.

Uncovering hidden genetic variation in photosynthesis of field-grown maize under ozone pollution.

Ozone is the most damaging air pollutant to crops, currently reducing Midwest US maize production by up to 10%, yet there has been very little effort to adapt germplasm for ozone tolerance. Ozone enters plants through stomata, reacts to form reactive oxygen species in the apoplast and ultimately decreases photosynthetic C gain. In this study, 10 diverse inbred parents were crossed in a half-diallel design to create 45 F1 hybrids, which were tested for ozone response in the field using free air concentration enrichment (FACE). Ozone stress increased the heritability of photosynthetic traits and altered genetic correlations among traits. Hybrids from parents Hp301 and NC338 showed greater sensitivity to ozone stress, and disrupted relationships among photosynthetic traits. The physiological responses underlying sensitivity to ozone differed in hybrids from the two parents, suggesting multiple mechanisms of response to oxidative stress. FACE technology was essential to this evaluation because genetic variation in photosynthesis under elevated ozone was not predictable based on performance at ambient ozone. These findings suggest that selection under elevated ozone is needed to identify deleterious alleles in the world's largest commodity crop.

A Robust Methodology for Assessing Differential Homeolog Contributions to the Transcriptomes of Allopolyploids.

Polyploidy has played a pivotal and recurring role in angiosperm evolution. Allotetraploids arise from hybridization between species and possess duplicated gene copies (homeologs) that serve redundant roles immediately after polyploidization. Although polyploidization is a major contributor to plant evolution, it remains poorly understood. We describe an analytical approach for assessing homeolog-specific expression that begins with de novo assembly of parental transcriptomes and effectively (i) reduces redundancy in de novo assemblies, (ii) identifies putative orthologs, (iii) isolates common regions between orthologs, and (iv) assesses homeolog-specific expression using a robust Bayesian Poisson-Gamma model to account for sequence bias when mapping polyploid reads back to parental references. Using this novel methodology, we examine differential homeolog contributions to the transcriptome in the recently formed allopolyploids Tragopogon mirus and T. miscellus (Compositae). Notably, we assess a larger Tragopogon gene set than previous studies of this system. Using carefully identified orthologous regions and filtering biased orthologs, we find in both allopolyploids largely balanced expression with no strong parental bias. These new methods can be used to examine homeolog expression in any tetrapolyploid system without requiring a reference genome.

SECIMTools: a suite of metabolomics data analysis tools.

BACKGROUND: Metabolomics has the promise to transform the area of personalized medicine with the rapid development of high throughput technology for untargeted analysis of metabolites. Open access, easy to use, analytic tools that are broadly accessible to the biological community need to be developed. While technology used in metabolomics varies, most metabolomics studies have a set of features identified. Galaxy is an open access platform that enables scientists at all levels to interact with big data. Galaxy promotes reproducibility by saving histories and enabling the sharing workflows among scientists. RESULTS: SECIMTools (SouthEast Center for Integrated Metabolomics) is a set of Python applications that are available both as standalone tools and wrapped for use in Galaxy. The suite includes a comprehensive set of quality control metrics (retention time window evaluation and various peak evaluation tools), visualization techniques (hierarchical cluster heatmap, principal component analysis, modular modularity clustering), basic statistical analysis methods (partial least squares - discriminant analysis, analysis of variance, t-test, Kruskal-Wallis non-parametric test), advanced classification methods (random forest, support vector machines), and advanced variable selection tools (least absolute shrinkage and selection operator LASSO and Elastic Net). CONCLUSIONS: SECIMTools leverages the Galaxy platform and enables integrated workflows for metabolomics data analysis made from building blocks designed for easy use and interpretability. Standard data formats and a set of utilities allow arbitrary linkages between tools to encourage novel workflow designs. The Galaxy framework enables future data integration for metabolomics studies with other omics data.

Intergenerational environmental effects: functional signals in offspring transcriptomes and metabolomes after parental jasmonic acid treatment in apomictic dandelion.

Parental environments can influence offspring traits. However, the magnitude of the impact of parental environments on offspring molecular phenotypes is poorly understood. Here, we test the direct effects and intergenerational effects of jasmonic acid (JA) treatment, which is involved in herbivory-induced defense signaling, on transcriptomes and metabolomes in apomictic common dandelion (Taraxacum officinale). In a full factorial crossed design with parental and offspring JA and control treatments, we performed leaf RNA-seq gene expression analysis, LC-MS metabolomics and total phenolics assays in offspring plants. Expression analysis, leveraged by a de novo assembled transcriptome, revealed an induced response to JA exposure that is consistent with known JA effects. The intergenerational effect of treatment was considerable: 307 of 858 detected JA-responsive transcripts were affected by parental JA treatment. In terms of the numbers of metabolites affected, the magnitude of the chemical response to parental JA exposure was c. 10% of the direct JA treatment response. Transcriptome and metabolome analyses both identified the phosphatidylinositol signaling pathway as a target of intergenerational JA effects. Our results highlight that parental environments can have substantial effects in offspring generations. Transcriptome and metabolome assays provide a basis for zooming in on the potential mechanisms of inherited JA effects.

Reliable Detection of Herpes Simplex Virus Sequence Variation by High-Throughput Resequencing.

High-throughput sequencing (HTS) has resulted in data for a number of herpes simplex virus (HSV) laboratory strains and clinical isolates. The knowledge of these sequences has been critical for investigating viral pathogenicity. However, the assembly of complete herpesviral genomes, including HSV, is complicated due to the existence of large repeat regions and arrays of smaller reiterated sequences that are commonly found in these genomes. In addition, the inherent genetic variation in populations of isolates for viruses and other microorganisms presents an additional challenge to many existing HTS sequence assembly pipelines. Here, we evaluate two approaches for the identification of genetic variants in HSV1 strains using Illumina short read sequencing data. The first, a reference-based approach, identifies variants from reads aligned to a reference sequence and the second, a de novo assembly approach, identifies variants from reads aligned to de novo assembled consensus sequences. Of critical importance for both approaches is the reduction in the number of low complexity regions through the construction of a non-redundant reference genome. We compared variants identified in the two methods. Our results indicate that approximately 85% of variants are identified regardless of the approach. The reference-based approach to variant discovery captures an additional 15% representing variants divergent from the HSV1 reference possibly due to viral passage. Reference-based approaches are significantly less labor-intensive and identify variants across the genome where de novo assembly-based approaches are limited to regions where contigs have been successfully assembled. In addition, regions of poor quality assembly can lead to false variant identification in de novo consensus sequences. For viruses with a well-assembled reference genome, a reference-based approach is recommended.

High-Throughput Phenotyping of Maize Leaf Physiological and Biochemical Traits Using Hyperspectral Reflectance.

High-throughput, noninvasive field phenotyping has revealed genetic variation in crop morphological, developmental, and agronomic traits, but rapid measurements of the underlying physiological and biochemical traits are needed to fully understand genetic variation in plant-environment interactions. This study tested the application of leaf hyperspectral reflectance (λ = 500-2,400 nm) as a high-throughput phenotyping approach for rapid and accurate assessment of leaf photosynthetic and biochemical traits in maize (Zea mays). Leaf traits were measured with standard wet-laboratory and gas-exchange approaches alongside measurements of leaf reflectance. Partial least-squares regression was used to develop a measure of leaf chlorophyll content, nitrogen content, sucrose content, specific leaf area, maximum rate of phosphoenolpyruvate carboxylation, [CO2]-saturated rate of photosynthesis, and leaf oxygen radical absorbance capacity from leaf reflectance spectra. Partial least-squares regression models accurately predicted five out of seven traits and were more accurate than previously used simple spectral indices for leaf chlorophyll, nitrogen content, and specific leaf area. Correlations among leaf traits and statistical inferences about differences among genotypes and treatments were similar for measured and modeled data. The hyperspectral reflectance approach to phenotyping was dramatically faster than traditional measurements, enabling over 1,000 rows to be phenotyped during midday hours over just 2 to 4 d, and offers a nondestructive method to accurately assess physiological and biochemical trait responses to environmental stress.

Buffering of Genetic Regulatory Networks in Drosophila melanogaster.

Regulatory variation in gene expression can be described by cis- and trans-genetic components. Here we used RNA-seq data from a population panel of Drosophila melanogaster test crosses to compare allelic imbalance (AI) in female head tissue between mated and virgin flies, an environmental change known to affect transcription. Indeed, 3048 exons (1610 genes) are differentially expressed in this study. A Bayesian model for AI, with an intersection test, controls type I error. There are ∼200 genes with AI exclusively in mated or virgin flies, indicating an environmental component of expression regulation. On average 34% of genes within a cross and 54% of all genes show evidence for genetic regulation of transcription. Nearly all differentially regulated genes are affected in cis, with an average of 63% of expression variation explained by the cis-effects. Trans-effects explain 8% of the variance in AI on average and the interaction between cis and trans explains an average of 11% of the total variance in AI. In both environments cis- and trans-effects are compensatory in their overall effect, with a negative association between cis- and trans-effects in 85% of the exons examined. We hypothesize that the gene expression level perturbed by cis-regulatory mutations is compensated through trans-regulatory mechanisms, e.g., trans and cis by trans-factors buffering cis-mutations. In addition, when AI is detected in both environments, cis-mated, cis-virgin, and trans-mated-trans-virgin estimates are highly concordant with 99% of all exons positively correlated with a median correlation of 0.83 for cis and 0.95 for trans We conclude that the gene regulatory networks (GRNs) are robust and that trans-buffering explains robustness.

Presence of arachidonoyl-carnitine is associated with adverse cardiometabolic responses in hypertensive patients treated with atenolol.

INTRODUCTION: Atenolol, a commonly prescribed ß blocker for hypertension, is also associated with adverse cardiometabolic effects such as hyperglycemia and dyslipidemia. Knowledge of the mechanistic underpinnings of these adverse effects of atenolol is incomplete. OBJECTIVE: We sought to identify biomarkers associated with risk for these untoward effects of atenolol. We measured baseline blood serum levels of acylcarnitines (ACs) that are involved in a host of different metabolic pathways, to establish associations with adverse cardiometabolic responses after atenolol treatment. METHODS: Serum samples from Caucasian hypertensive patients (n = 224) who were treated with atenolol in the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) study were interrogated using a quantitative LC/MS assay for a large number of unique ACs in serum. For the 23 ACs that were detected in serum from ≥80 % of all patients, we conducted linear regression for changes in cardiometabolic factors with baseline AC levels, baseline cardiometabolic factors, age, sex, and BMI as covariates. For the 5 ACs that were detected in serum from 20 to 79 % of the patients, we similarly modeled changes in cardiometabolic factors, but with specifying the AC as present/absent in the regression. RESULTS: Among the 28 ACs, the presence (vs. absence) of arachidonoyl-carnitine (C20:4) was significantly associated with increased glucose (p = 0.0002), and was nominally associated with decreased plasma HDL-C (p = 0.017) and with less blood pressure (BP) lowering (p = 0.006 for systolic BP, p = 0.002 for diastolic BP), after adjustment. CONCLUSION: Serum level of C20:4 is a promising biomarker to predict adverse cardiometabolic responses including glucose and poor antihypertensive response to atenolol.

What the X has to do with it: differences in regulatory variability between the sexes in Drosophila simulans.

The mechanistic basis of regulatory variation and the prevailing evolutionary forces shaping that variation are known to differ between sexes and between chromosomes. Regulatory variation of gene expression can be due to functional changes within a gene itself (cis) or in other genes elsewhere in the genome (trans). The evolutionary properties of cis mutations are expected to differ from mutations affecting gene expression in trans. We analyze allele-specific expression across a set of X substitution lines in intact adult Drosophila simulans to evaluate whether regulatory variation differs for cis and trans, for males and females, and for X-linked and autosomal genes. Regulatory variation is common (56% of genes), and patterns of variation within D. simulans are consistent with previous observations in Drosophila that there is more cis than trans variation within species (47% vs. 25%, respectively). The relationship between sex-bias and sex-limited variation is remarkably consistent across sexes. However, there are differences between cis and trans effects: cis variants show evidence of purifying selection in the sex toward which expression is biased, while trans variants do not. For female-biased genes, the X is depleted for trans variation in a manner consistent with a female-dominated selection regime on the X. Surprisingly, there is no evidence for depletion of trans variation for male-biased genes on X. This is evidence for regulatory feminization of the X, trans-acting factors controlling male-biased genes are more likely to be found on the autosomes than those controlling female-biased genes.

Ago HITS-CLIP expands understanding of Kaposi's sarcoma-associated herpesvirus miRNA function in primary effusion lymphomas.

KSHV is the etiological agent of Kaposi's sarcoma (KS), primary effusion lymphoma (PEL), and a subset of multicentricCastleman's disease (MCD). The fact that KSHV-encoded miRNAs are readily detectable in all KSHV-associated tumors suggests a potential role in viral pathogenesis and tumorigenesis. MiRNA-mediated regulation of gene expression is a complex network with each miRNA having many potential targets, and to date only few KSHV miRNA targets have been experimentally determined. A detailed understanding of KSHV miRNA functions requires high-through putribonomics to globally analyze putative miRNA targets in a cell type-specific manner. We performed Ago HITS-CLIP to identify viral and cellular miRNAs and their cognate targets in two latently KSHV-infected PEL cell lines. Ago HITS-CLIP recovered 1170 and 950 cellular KSHV miRNA targets from BCBL-1 and BC-3, respectively. Importantly, enriched clusters contained KSHV miRNA seed matches in the 3'UTRs of numerous well characterized targets, among them THBS1, BACH1, and C/EBPß. KSHV miRNA targets were strongly enriched for genes involved in multiple pathways central for KSHV biology, such as apoptosis, cell cycle regulation, lymphocyte proliferation, and immune evasion, thus further supporting a role in KSHV pathogenesis and potentially tumorigenesis. A limited number of viral transcripts were also enriched by HITS-CLIP including vIL-6 expressed only in a subset of PEL cells during latency. Interestingly, Ago HITS-CLIP revealed extremely high levels of Ago-associated KSHV miRNAs especially in BC-3 cells where more than 70% of all miRNAs are of viral origin. This suggests that in addition to seed match-specific targeting of cellular genes, KSHV miRNAs may also function by hijacking RISCs, thereby contributing to a global de-repression of cellular gene expression due to the loss of regulation by human miRNAs. In summary, we provide an extensive list of cellular and viral miRNA targets representing an important resource to decipher KSHV miRNA function.

RNA-seq: technical variability and sampling.

BACKGROUND: RNA-seq is revolutionizing the way we study transcriptomes. mRNA can be surveyed without prior knowledge of gene transcripts. Alternative splicing of transcript isoforms and the identification of previously unknown exons are being reported. Initial reports of differences in exon usage, and splicing between samples as well as quantitative differences among samples are beginning to surface. Biological variation has been reported to be larger than technical variation. In addition, technical variation has been reported to be in line with expectations due to random sampling. However, strategies for dealing with technical variation will differ depending on the magnitude. The size of technical variance, and the role of sampling are examined in this manuscript. RESULTS: In this study three independent Solexa/Illumina experiments containing technical replicates are analyzed. When coverage is low, large disagreements between technical replicates are apparent. Exon detection between technical replicates is highly variable when the coverage is less than 5 reads per nucleotide and estimates of gene expression are more likely to disagree when coverage is low. Although large disagreements in the estimates of expression are observed at all levels of coverage. CONCLUSIONS: Technical variability is too high to ignore. Technical variability results in inconsistent detection of exons at low levels of coverage. Further, the estimate of the relative abundance of a transcript can substantially disagree, even when coverage levels are high. This may be due to the low sampling fraction and if so, it will persist as an issue needing to be addressed in experimental design even as the next wave of technology produces larger numbers of reads. We provide practical recommendations for dealing with the technical variability, without dramatic cost increases.

Nonflowering plants possess a unique folate-dependent phenylalanine hydroxylase that is localized in chloroplasts.

Tetrahydropterin-dependent aromatic amino acid hydroxylases (AAHs) are known from animals and microbes but not plants. A survey of genomes and ESTs revealed AAH-like sequences in gymnosperms, mosses, and algae. Analysis of full-length AAH cDNAs from Pinus taeda, Physcomitrella patens, and Chlamydomonas reinhardtii indicated that the encoded proteins form a distinct clade within the AAH family. These proteins were shown to have Phe hydroxylase activity by functional complementation of an Escherichia coli Tyr auxotroph and by enzyme assays. The P. taeda and P. patens AAHs were specific for Phe, required iron, showed Michaelian kinetics, and were active as monomers. Uniquely, they preferred 10-formyltetrahydrofolate to any physiological tetrahydropterin as cofactor and, consistent with preferring a folate cofactor, retained activity in complementation tests with tetrahydropterin-depleted E. coli host strains. Targeting assays in Arabidopsis thaliana mesophyll protoplasts using green fluorescent protein fusions, and import assays with purified Pisum sativum chloroplasts, indicated chloroplastic localization. Targeting assays further indicated that pterin-4a-carbinolamine dehydratase, which regenerates the AAH cofactor, is also chloroplastic. Ablating the single AAH gene in P. patens caused accumulation of Phe and caffeic acid esters. These data show that nonflowering plants have functional plastidial AAHs, establish an unprecedented electron donor role for a folate, and uncover a novel link between folate and aromatic metabolism.