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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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).

The SUMO conjugation pathway in Populus: genomic analysis, tissue-specific and inducible SUMOylation and in vitro de-SUMOylation.

Covalent attachment of the small ubiquitin-like modifier (SUMO) to proteins in eukaryotic cells can regulate an assortment of cellular processes including transcription, and DNA-protein and protein-protein interactions. We identified gene models and found evidence for expression of genes involved in SUMOylation and SUMO deconjugation in Populus. We detected SUMOylated proteins in diverse organ and tissue types. SUMOylation was altered during responses to heat shock, desiccation, peroxide and irrigation of roots with high salt solution. SUMO deconjugation from substrates was sensitive to cysteine protease inhibitors. Product sizes and sensitivity to inhibitors are consistent with poly-SUMO chain formation as an intermediate step in SUMO redistribution to substrates in plant cells responding to treatments. The SUMOylation pathway is active in Populus and substrate conjugation to SUMO is a rapid response to multiple inducers.

Post-translational modification of an R2R3-MYB transcription factor by a MAP Kinase during xylem development.

Despite the pivotal role played by R2R3-MYB family members in the regulation of plant gene expression, little is known about post-translational regulation of these proteins. In animals, the MYB family member, c-MYB, is post-translationally modified by a mitogen-activated protein kinase (MAPK), p42(mapk). In order to test the hypothesis that R2R3-MYB proteins may be regulated by MAPK activity, interplay between a R2R3-MYB family member expressed in differentiating pine xylem (Pinus taeda MYB4, PtMYB4) and MAPK proteins expressed in the same tissue was examined. One of the MAPK proteins expressed in pine xylem, PtMAPK6, phosphorylated PtMYB4. Recombinant PtMAPK6 phosphorylated PtMYB4 on serine-236, located in the C-terminal activation domain of this transcription factor in a context that is found in other plant MYB proteins. Modification of the PtMAPK6 target serine in PtMYB4 did not appear to alter DNA binding in vitro but did alter the ability of PtMYB4 to promote transcriptional activation in yeast. PtMAPK6 activity was detected in developing xylem cells that had ceased cell division and formed secondary walls. Together, the data support a role for PtMAPK6 during early xylem development and suggest a function for this kinase in regulating gene expression through phosphorylation of PtMYB4.

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.

Salicylate and catechol levels are maintained in nahG transgenic poplar.

Metabolic profiling was used to investigate the molecular phenotypes of a transgenic Populus tremula x P. alba hybrid expressing the nahG transgene, a bacterial gene encoding salicylate hydroxylase that converts salicylic acid to catechol. Despite the efficacy of this transgenic approach to reduce salicylic acid levels in other model systems and thereby elucidate roles for salicylic acid in plant signaling, transgenic poplars had similar foliar levels of salicylic acid and catechol to that of non-transformed controls and exhibited no morphological phenotypes. To gain a deeper understanding of the basis for these observations, we analyzed metabolic profiles of leaves as influenced by transgene expression. Expression of nahG decreased quinic acid conjugates and increased catechol glucoside, while exerting little effect on levels of salicylic acid and catechol, the substrate and product, respectively, of the nahG enzyme. This suggests a biological role of elevated constitutive salicylic acid levels in Populus, in contrast to other plant systems in which nahG dramatically reduces salicylic acid levels.

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.