Identifying reference genes with stable expression from high throughput sequence data

© The Author(s), 2012. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 3 (2012): 385, doi:10.3389/fmicb.2012.00385.Genes that are constitutively expressed across multiple environmental stimuli are crucial to quantifying differentially expressed genes, particularly when employing quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) assays. However, the identification of these potential reference genes in non-model organisms is challenging and is often guided by expression patterns in distantly related organisms. Here, transcriptome datasets from the diatom Thalassiosira pseudonana grown under replete, phosphorus-limited, iron-limited, and phosphorus and iron co-limited nutrient regimes were analyzed through literature-based searches for homologous reference genes, k-means clustering, and analysis of sequence counts (ASC) to identify putative reference genes. A total of 9759 genes were identified and screened for stable expression. Literature-based searches surveyed 18 generally accepted reference genes, revealing 101 homologs in T. pseudonana with variable expression and a wide range of mean tags per million. k-means analysis parsed the whole transcriptome into 15 clusters. The two most stable clusters contained 709 genes, but still had distinct patterns in expression. ASC analyses identified 179 genes that were stably expressed (posterior probability < 0.1 for 1.25 fold change). Genes known to have a stable expression pattern across the test treatments, like actin, were identified in this pool of 179 candidate genes. ASC can be employed on data without biological replicates and was more robust than the k-means approach in isolating genes with stable expression. The intersection of the genes identified through ASC with commonly used reference genes from the literature suggests that actin and ubiquitin ligase may be useful reference genes for T. pseudonana and potentially other diatoms. With the wealth of transcriptome sequence data becoming available, ASC can be easily applied to transcriptome datasets from other phytoplankton to identify reference genes.This research was funded by the National Science Foundation grant #OCE-0723667 (to Sonya T. Dyhrman, Mak A. Saito, Bethany D. Jenkins, and Tatiana A. Rynearson). Harriet Alexander is funded under a National Defense Science and Engineering Graduate (NDSEG) Fellowship

The transcriptome and proteome of the diatom Thalassiosira pseudonana reveal a diverse phosphorus stress response

© The Author(s), 2012. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in PLoS One 7 (2012): e33768, doi:10.1371/journal.pone.0033768.Phosphorus (P) is a critical driver of phytoplankton growth and ecosystem function in the ocean. Diatoms are an abundant class of marine phytoplankton that are responsible for significant amounts of primary production. With the control they exert on the oceanic carbon cycle, there have been a number of studies focused on how diatoms respond to limiting macro and micronutrients such as iron and nitrogen. However, diatom physiological responses to P deficiency are poorly understood. Here, we couple deep sequencing of transcript tags and quantitative proteomics to analyze the diatom Thalassiosira pseudonana grown under P-replete and P-deficient conditions. A total of 318 transcripts were differentially regulated with a false discovery rate of <0.05, and a total of 136 proteins were differentially abundant (p<0.05). Significant changes in the abundance of transcripts and proteins were observed and coordinated for multiple biochemical pathways, including glycolysis and translation. Patterns in transcript and protein abundance were also linked to physiological changes in cellular P distributions, and enzyme activities. These data demonstrate that diatom P deficiency results in changes in cellular P allocation through polyphosphate production, increased P transport, a switch to utilization of dissolved organic P through increased production of metalloenzymes, and a remodeling of the cell surface through production of sulfolipids. Together, these findings reveal that T. pseudonana has evolved a sophisticated response to P deficiency involving multiple biochemical strategies that are likely critical to its ability to respond to variations in environmental P availability.This research was supported by the National Science Foundation (NSF) Environmental Genomics and NSF Biological Oceanography Program through grant OCE-0723667 to Dr. Dyhrman, Dr. Jenkins, Dr. Saito, and Dr. Rynearson, the NSF Chemical Oceanography Program through grant OCE-0549794 to Dr. Dyhrman and OCE-0526800 to Dr. Jenkins, the G. B. Moore Foundation and OCE-0752291 to Dr. Saito, NSF-EPSCoR (NSF-0554548 & NSF-1004057) to the University of Rhode Island, the Center for Microbial Oceanography: Research and Education, and the Joint Genome Institute/DOE Community Sequencing Program (CSP795793) to Dr. Jenkins, Dr. Dyhrman, Dr. Rynearson and Dr. Saito

Empirical bayes analysis of sequencing-based transcriptional profiling without replicates

Background: Recent technological advancements have made high throughput sequencing an increasingly popular approach for transcriptome analysis. Advantages of sequencing-based transcriptional profiling over microarrays have been reported, including lower technical variability. However, advances in technology do not remove biological variation between replicates and this variation is often neglected in many analyses. Results: We propose an empirical Bayes method, titled Analysis of Sequence Counts (ASC), to detect differential expression based on sequencing technology. ASC borrows information across sequences to establish prior distribution of sample variation, so that biological variation can be accounted for even when replicates are not available. Compared to current approaches that simply tests for equality of proportions in two samples, ASC is less biased towards highly expressed sequences and can identify more genes with a greater log fold change at lower overall abundance. Conclusions: ASC unifies the biological and statistical significance of differential expression by estimating the posterior mean of log fold change and estimating false discovery rates based on the posterior mean. The implementation in R is available at http://www.stat.brown.edu/Zwu/research.aspx

A new geological map of the Lau Basin (southwestern Pacific Ocean) reveals crustal growth processes in arc-backarc systems

A 1:1,000,000-scale lithostratigraphic assemblage map of the Lau Basin (southwestern Pacific Ocean) has been created using remote predictive mapping (RPM) techniques developed by geological surveys on land. Formation-level geological units were identified in training sets at scales of 1:100,000–1:200,000 in different parts of the basin and then extrapolated to the areas where geological data are sparse. The final compilation is presented together with a quantitative analysis of assemblage-level crustal growth based on area-age relationships of the assigned units. The data sets used to develop mapping criteria and an internally consistent legend for the compilation included high-resolution ship-based multibeam, satellite- and ship-based gravity, magnetics, seafloor imaging, and sampling data. The correlation of units was informed by published geochronological information and kinematic models of basin opening. The map covers >1,000,000 km2 of the Lau-Tonga arc-backarc system, subdivided into nine assemblage types: forearc crust (9% by area), crust of the active volcanic arc (7%), backarc rifts and spreading centers (20%), transitional arc-backarc crust (13%), relict arc crust (38%), relict backarc crust (8%), and undivided arc-backarc assemblages (<5%), plus oceanic assemblages, intraplate volcanoes, and carbonate platforms. Major differences in the proportions of assemblage types compared to other intraoceanic subduction systems (e.g., Mariana backarc, North Fiji Basin) underscore the complex geological makeup of the Lau Basin. Backarc crust formed and is forming simultaneously at 12 different locations in the basin in response to widely distributed extension, and this is considered to be a dominant pattern of crustal accretion in large arc-backarc systems. Accelerated basin opening and a microplate breakout north of the Peggy Ridge has been accommodated by seven different spreading centers. The result is an intricate mosaic of small intact assemblages in the north of the basin, compared to fewer and larger assemblages in the south. Although the oldest rocks are Eocene (~40 m.y. old basement of the Lau and Tonga Ridges), half of the backarc crust in the map area formed within the last 3 m.y. and therefore represents some of the fastest growing crust on Earth, associated with prolific magmatic and hydro-thermal activity. These observations provide important clues to the geological evolution and makeup of ancient backarc basins and to processes of crustal growth that ultimately lead to the emergence of continents

Histone H3.3 beyond cancer: Germline mutations in Histone 3 Family 3A and 3B cause a previously unidentified neurodegenerative disorder in 46 patients

Although somatic mutations in Histone 3.3 (H3.3) are well-studied drivers of oncogenesis, the role of germline mutations remains unreported. We analyze 46 patients bearing de novo germline mutations in histone 3 family 3A (H3F3A) or H3F3B with progressive neurologic dysfunction and congenital anomalies without malignancies. Molecular modeling of all 37 variants demonstrated clear disruptions in interactions with DNA, other histones, and histone chaperone proteins. Patient histone posttranslational modifications (PTMs) analysis revealed notably aberrant local PTM patterns distinct from the somatic lysine mutations that cause global PTM dysregulation. RNA sequencing on patient cells demonstrated up-regulated gene expression related to mitosis and cell division, and cellular assays confirmed an increased proliferative capacity. A zebrafish model showed craniofacial anomalies and a defect in Foxd3-derived glia. These data suggest that the mechanism of germline mutations are distinct from cancer-associated somatic histone mutations but may converge on control of cell proliferation

Effects of eight neuropsychiatric copy number variants on human brain structure

Many copy number variants (CNVs) confer risk for the same range of neurodevelopmental symptoms and psychiatric conditions including autism and schizophrenia. Yet, to date neuroimaging studies have typically been carried out one mutation at a time, showing that CNVs have large effects on brain anatomy. Here, we aimed to characterize and quantify the distinct brain morphometry effects and latent dimensions across 8 neuropsychiatric CNVs. We analyzed T1-weighted MRI data from clinically and non-clinically ascertained CNV carriers (deletion/duplication) at the 1q21.1 (n = 39/28), 16p11.2 (n = 87/78), 22q11.2 (n = 75/30), and 15q11.2 (n = 72/76) loci as well as 1296 non-carriers (controls). Case-control contrasts of all examined genomic loci demonstrated effects on brain anatomy, with deletions and duplications showing mirror effects at the global and regional levels. Although CNVs mainly showed distinct brain patterns, principal component analysis (PCA) loaded subsets of CNVs on two latent brain dimensions, which explained 32 and 29% of the variance of the 8 Cohen’s d maps. The cingulate gyrus, insula, supplementary motor cortex, and cerebellum were identified by PCA and multi-view pattern learning as top regions contributing to latent dimension shared across subsets of CNVs. The large proportion of distinct CNV effects on brain morphology may explain the small neuroimaging effect sizes reported in polygenic psychiatric conditions. Nevertheless, latent gene brain morphology dimensions will help subgroup the rapidly expanding landscape of neuropsychiatric variants and dissect the heterogeneity of idiopathic conditions

Effect of angiotensin-converting enzyme inhibitor and angiotensin receptor blocker initiation on organ support-free days in patients hospitalized with COVID-19

IMPORTANCE Overactivation of the renin-angiotensin system (RAS) may contribute to poor clinical outcomes in patients with COVID-19. Objective To determine whether angiotensin-converting enzyme (ACE) inhibitor or angiotensin receptor blocker (ARB) initiation improves outcomes in patients hospitalized for COVID-19. DESIGN, SETTING, AND PARTICIPANTS In an ongoing, adaptive platform randomized clinical trial, 721 critically ill and 58 non–critically ill hospitalized adults were randomized to receive an RAS inhibitor or control between March 16, 2021, and February 25, 2022, at 69 sites in 7 countries (final follow-up on June 1, 2022). INTERVENTIONS Patients were randomized to receive open-label initiation of an ACE inhibitor (n = 257), ARB (n = 248), ARB in combination with DMX-200 (a chemokine receptor-2 inhibitor; n = 10), or no RAS inhibitor (control; n = 264) for up to 10 days. MAIN OUTCOMES AND MEASURES The primary outcome was organ support–free days, a composite of hospital survival and days alive without cardiovascular or respiratory organ support through 21 days. The primary analysis was a bayesian cumulative logistic model. Odds ratios (ORs) greater than 1 represent improved outcomes. RESULTS On February 25, 2022, enrollment was discontinued due to safety concerns. Among 679 critically ill patients with available primary outcome data, the median age was 56 years and 239 participants (35.2%) were women. Median (IQR) organ support–free days among critically ill patients was 10 (–1 to 16) in the ACE inhibitor group (n = 231), 8 (–1 to 17) in the ARB group (n = 217), and 12 (0 to 17) in the control group (n = 231) (median adjusted odds ratios of 0.77 [95% bayesian credible interval, 0.58-1.06] for improvement for ACE inhibitor and 0.76 [95% credible interval, 0.56-1.05] for ARB compared with control). The posterior probabilities that ACE inhibitors and ARBs worsened organ support–free days compared with control were 94.9% and 95.4%, respectively. Hospital survival occurred in 166 of 231 critically ill participants (71.9%) in the ACE inhibitor group, 152 of 217 (70.0%) in the ARB group, and 182 of 231 (78.8%) in the control group (posterior probabilities that ACE inhibitor and ARB worsened hospital survival compared with control were 95.3% and 98.1%, respectively). CONCLUSIONS AND RELEVANCE In this trial, among critically ill adults with COVID-19, initiation of an ACE inhibitor or ARB did not improve, and likely worsened, clinical outcomes. TRIAL REGISTRATION ClinicalTrials.gov Identifier: NCT0273570

HIV infection of dendritic cells subverts the IFN induction pathway via IRF-1 and inhibits type 1 IFN production

Many viruses have developed mechanisms to evade the IFN response. Here, HIV-1 was shown to induce a distinct subset of IFN-stimulated genes (ISGs) in monocyte-derived dendritic cells (DCs), without detectable type I or II IFN. These ISGs all contained an IFN regulatory factor 1 (IRF-1) binding site in their promoters, and their expression was shown to be driven by IRF-1, indicating this subset was induced directly by viral infection by IRF-1. IRF-1 and -7 protein expression was enriched in HIV p24 antigen-positive DCs. A HIV deletion mutant with the IRF-1 binding site deleted from the long terminal repeat showed reduced growth kinetics. Early and persistent induction of IRF-1 was coupled with sequential transient up-regulation of its 2 inhibitors, IRF-8, followed by IRF-2, suggesting a mechanism for IFN inhibition. HIV-1 mutants with Vpr deleted induced IFN, showing that Vpr is inhibitory. However, HIV IFN inhibition was mediated by failure of IRF-3 activation rather than by its degradation, as in T cells. In contrast, herpes simplex virus type 2 markedly induced IFN&beta; and a broader range of ISGs to higher levels, supporting the hypothesis that HIV-1 specifically manipulates the induction of IFN and ISGs to enhance its noncytopathic replication in DCs. <br /