Differential expression analysis for sequence count data

*Motivation:* High-throughput nucleotide sequencing provides quantitative readouts in assays for RNA expression (RNA-Seq), protein-DNA binding (ChIP-Seq) or cell counting (barcode sequencing). Statistical inference of differential signal in such data requires estimation of their variability throughout the dynamic range. When the number of replicates is small, error modelling is needed to achieve statistical power.

*Results:* We propose an error model that uses the negative binomial distribution, with variance and mean linked by local regression, to model the null distribution of the count data. The method controls type-I error and provides good detection power. 

*Availability:* A free open-source R software package, _DESeq_, is available from the Bioconductor project and from "http://www-huber.embl.de/users/anders/DESeq":http://www-huber.embl.de/users/anders/DESeq

Expression of hereditary hemochromatosis C282Y HFE protein in HEK293 cells activates specific endoplasmic reticulum stress responses

<p>Abstract</p> <p>Background</p> <p>Hereditary Hemochromatosis (HH) is a genetic disease associated with iron overload, in which individuals homozygous for the mutant C282Y <it>HFE </it>associated allele are at risk for the development of a range of disorders particularly liver disease. Conformational diseases are a class of disorders associated with the expression of misfolded protein. HFE C282Y is a mutant protein that does not fold correctly and consequently is retained in the Endoplasmic Reticulum (ER). In this context, we sought to identify ER stress signals associated with mutant C282Y HFE protein expression, which may have a role in the molecular pathogenesis of HH.</p> <p>Results</p> <p>Vector constructs of Wild type HFE and Mutant C282Y HFE were made and transfected into HEK293 cell lines. We have shown that expression of C282Y HFE protein triggers both an unfolded protein response (UPR), as revealed by the increased GRP78, ATF6 and CHOP expression, and an ER overload response (EOR), as indicated by NF-κB activation. Furthermore, C282Y HFE protein induced apoptotic responses associated with activation of ER stress. Inhibition studies demonstrated that tauroursodeoxycholic acid, an endogenous bile acid, downregulates these events. Finally, we found that the co-existence of both C282Y HFE and Z alpha 1-antitrypsin protein (the protein associated with the liver disease of Z alpha 1-antitrypsin deficiency) expression on ER stress responses acted as potential disease modifiers with respect to each other.</p> <p>Conclusion</p> <p>Our novel observations suggest that both the ER overload response (EOR) and the unfolded protein response (UPR) are activated by mutant C282Y HFE protein.</p

Quantitative Fitness Analysis Shows That NMD Proteins and Many Other Protein Complexes Suppress or Enhance Distinct Telomere Cap Defects

To better understand telomere biology in budding yeast, we have performed systematic suppressor/enhancer analyses on yeast strains containing a point mutation in the essential telomere capping gene CDC13 (cdc13-1) or containing a null mutation in the DNA damage response and telomere capping gene YKU70 (yku70Δ). We performed Quantitative Fitness Analysis (QFA) on thousands of yeast strains containing mutations affecting telomere-capping proteins in combination with a library of systematic gene deletion mutations. To perform QFA, we typically inoculate 384 separate cultures onto solid agar plates and monitor growth of each culture by photography over time. The data are fitted to a logistic population growth model; and growth parameters, such as maximum growth rate and maximum doubling potential, are deduced. QFA reveals that as many as 5% of systematic gene deletions, affecting numerous functional classes, strongly interact with telomere capping defects. We show that, while Cdc13 and Yku70 perform complementary roles in telomere capping, their genetic interaction profiles differ significantly. At least 19 different classes of functionally or physically related proteins can be identified as interacting with cdc13-1, yku70Δ, or both. Each specific genetic interaction informs the roles of individual gene products in telomere biology. One striking example is with genes of the nonsense-mediated RNA decay (NMD) pathway which, when disabled, suppress the conditional cdc13-1 mutation but enhance the null yku70Δ mutation. We show that the suppressing/enhancing role of the NMD pathway at uncapped telomeres is mediated through the levels of Stn1, an essential telomere capping protein, which interacts with Cdc13 and recruitment of telomerase to telomeres. We show that increased Stn1 levels affect growth of cells with telomere capping defects due to cdc13-1 and yku70Δ. QFA is a sensitive, high-throughput method that will also be useful to understand other aspects of microbial cell biology

Dynamic changes in eIF4F-mRNA interactions revealed by global analyses of environmental stress responses

BACKGROUND: Translation factors eIF4E and eIF4G form eIF4F, which interacts with the messenger RNA (mRNA) 5' cap to promote ribosome recruitment and translation initiation. Variations in the association of eIF4F with individual mRNAs likely contribute to differences in translation initiation frequencies between mRNAs. As translation initiation is globally reprogrammed by environmental stresses, we were interested in determining whether eIF4F interactions with individual mRNAs are reprogrammed and how this may contribute to global environmental stress responses. RESULTS: Using a tagged-factor protein capture and RNA-sequencing (RNA-seq) approach, we have assessed how mRNA associations with eIF4E, eIF4G1 and eIF4G2 change globally in response to three defined stresses that each cause a rapid attenuation of protein synthesis: oxidative stress induced by hydrogen peroxide and nutrient stresses caused by amino acid or glucose withdrawal. We find that acute stress leads to dynamic and unexpected changes in eIF4F-mRNA interactions that are shared among each factor and across the stresses imposed. eIF4F-mRNA interactions stabilised by stress are predominantly associated with translational repression, while more actively initiating mRNAs become relatively depleted for eIF4F. Simultaneously, other mRNAs are insulated from these stress-induced changes in eIF4F association. CONCLUSION: Dynamic eIF4F-mRNA interaction changes are part of a coordinated early translational control response shared across environmental stresses. Our data are compatible with a model where multiple mRNA closed-loop complexes form with differing stability. Hence, unexpectedly, in the absence of other stabilising factors, rapid translation initiation on mRNAs correlates with less stable eIF4F interactions

Path and Ridge Regression Analysis of Seed Yield and Seed Yield Components of Russian Wildrye (Psathyrostachys juncea Nevski) under Field Conditions

The correlations among seed yield components, and their direct and indirect effects on the seed yield (Z) of Russina wildrye (Psathyrostachys juncea Nevski) were investigated. The seed yield components: fertile tillers m-2 (Y1), spikelets per fertile tillers (Y2), florets per spikelet- (Y3), seed numbers per spikelet (Y4) and seed weight (Y5) were counted and the Z were determined in field experiments from 2003 to 2006 via big sample size. Y1 was the most important seed yield component describing the Z and Y2 was the least. The total direct effects of the Y1, Y3 and Y5 to the Z were positive while Y4 and Y2 were weakly negative. The total effects (directs plus indirects) of the components were positively contributed to the Z by path analyses. The seed yield components Y1, Y2, Y4 and Y5 were significantly (P<0.001) correlated with the Z for 4 years totally, while in the individual years, Y2 were not significant correlated with Y3, Y4 and Y5 by Peason correlation analyses in the five components in the plant seed production. Therefore, selection for high seed yield through direct selection for large Y1, Y2 and Y3 would be effective for breeding programs in grasses. Furthermore, it is the most important that, via ridge regression, a steady algorithm model between Z and the five yield components was founded, which can be closely estimated the seed yield via the components

Darwinism, organizational evolution and survival: key challenges for future research

How do social organizations evolve? How do they adapt to environmental pressures? What resources and capabilities determine their survival within dynamic competition? Charles Darwin’s seminal work The Origin of Species (1859) has provided a significant impact on the development of the management and organization theory literatures on organizational evolution. This article introduces the JMG Special Issue focused on Darwinism, organizational evolution and survival. We discuss key themes in the organizational evolution research that have emerged in recent years. These include the increasing adoption of the co-evolutionary approach, with a particular focus on the definition of appropriate units of analysis, such as routines, and related challenges associated with exploring the relationship between co-evolution, re-use of knowledge, adaptation, and exaptation processes. We then introduce the three articles that we have finally accepted in this Special Issue after an extensive, multi-round, triple blind-review process. We briefly outline how each of these articles contributes to understanding among scholars, practitioners and policy makers of the continuous evolutionary processes within and among social organizations and systems

The Human Sweet Tooth

Humans love the taste of sugar and the word "sweet" is used to describe not only this basic taste quality but also something that is desirable or pleasurable, e.g., la dolce vita. Although sugar or sweetened foods are generally among the most preferred choices, not everyone likes sugar, especially at high concentrations. The focus of my group's research is to understand why some people have a sweet tooth and others do not. We have used genetic and molecular techniques in humans, rats, mice, cats and primates to understand the origins of sweet taste perception. Our studies demonstrate that there are two sweet receptor genes (TAS1R2 and TAS1R3), and alleles of one of the two genes predict the avidity with which some mammals drink sweet solutions. We also find a relationship between sweet and bitter perception. Children who are genetically more sensitive to bitter compounds report that very sweet solutions are more pleasant and they prefer sweet carbonated beverages more than milk, relative to less bitter-sensitive peers. Overall, people differ in their ability to perceive the basic tastes, and particular constellations of genes and experience may drive some people, but not others, toward a caries-inducing sweet diet. Future studies will be designed to understand how a genetic preference for sweet food and drink might contribute to the development of dental caries