Gene dispersion is the key determinant of the read count bias in differential expression analysis of RNA-seq data

Background: In differential expression analysis of RNA-sequencing (RNA-seq) read count data for two sample groups, it is known that highly expressed genes (or longer genes) are more likely to be differentially expressed which is called read count bias (or gene length bias). This bias had great effect on the downstream Gene Ontology over-representation analysis. However, such a bias has not been systematically analyzed for different replicate types of RNA-seq data. Results: We show that the dispersion coefficient of a gene in the negative binomial modeling of read counts is the critical determinant of the read count bias (and gene length bias) by mathematical inference and tests for a number of simulated and real RNA-seq datasets. We demonstrate that the read count bias is mostly confined to data with small gene dispersions (e.g., technical replicates and some of genetically identical replicates such as cell lines or inbred animals), and many biological replicate data from unrelated samples do not suffer from such a bias except for genes with some small counts. It is also shown that the sample-permuting GSEA method yields a considerable number of false positives caused by the read count bias, while the preranked method does not. Conclusion: We showed the small gene variance (similarly, dispersion) is the main cause of read count bias (and gene length bias) for the first time and analyzed the read count bias for different replicate types of RNA-seq data and its effect on gene-set enrichment analysis

A Dynamical Model of Oocyte Maturation Unveils Precisely Orchestrated Meiotic Decisions

Maturation of vertebrate oocytes into haploid gametes relies on two consecutive meioses without intervening DNA replication. The temporal sequence of cellular transitions driving eggs from G2 arrest to meiosis I (MI) and then to meiosis II (MII) is controlled by the interplay between cyclin-dependent and mitogen-activated protein kinases. In this paper, we propose a dynamical model of the molecular network that orchestrates maturation of Xenopus laevis oocytes. Our model reproduces the core features of maturation progression, including the characteristic non-monotonous time course of cyclin-Cdks, and unveils the network design principles underlying a precise sequence of meiotic decisions, as captured by bifurcation and sensitivity analyses. Firstly, a coherent and sharp meiotic resumption is triggered by the concerted action of positive feedback loops post-translationally activating cyclin-Cdks. Secondly, meiotic transition is driven by the dynamic antagonism between positive and negative feedback loops controlling cyclin turnover. Our findings reveal a highly modular network in which the coordination of distinct regulatory schemes ensures both reliable and flexible cell-cycle decisions

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation

Preterm low birthweight and the role of oral bacteria

Preterm and low birthweight (PTLBW) continues to be a major cause of mortality and morbidity across the world. In recent years, maternal periodontal disease has been implicated as a risk factor for adverse pregnancy outcomes. There is conflicting evidence to support such an outcome as illustrated by descriptive, case control and randomised controlled trials involving pregnant women from across the world, using different measurement tools to determine the level of periodontal disease. Whilst considering the literature, there is evidence for both arguments, based on the effect of periodontal inflammatory by products. Bacteria associated with periodontal disease are not dissimilar to those known to be associated with genito-urinary bacterial infections and adverse pregnancy outcomes. Several groups have demonstrated the apparent translocation of Fusobacterium nucleatum, Prevotella nigrescens, Prevotella intermedia, Porphyromonus gingivalis, Treponema denticola to the foetal placental unit whereby a maternal or foetal response has been detected resulting in premature birth or low birthweight. The normal process of parturition involves a cascade of events including a build-up of inflammatory mediators as linked to inflammation, whereby the maternal environment becomes hostile and threatens the well-being of the infant, and the foetus expelled. The question remains therefore, is there a greater risk of delivering a PTLBW infant when the mother has detectable periodontal disease, or is the release of inflammatory mediators and their translocation via the haematogenous route sufficient to induce a poor pregnancy outcome? The data investigated would suggest that there is a positive outcome when certain oral gram-negative bacteria create a cumulative effect sufficient to trigger early delivery, which represents the final straw to result in preterm or low birthweight delivery. There is equally sufficient epidemiological evidence that does not support this outcome, but it is agreed that maintaining oral health during pregnancy is beneficial to the mother and her infant

The Naturally Occurring YMDD Mutation among Patients Chronically Infected HBV and Untreated with Lamivudine: A Systematic Review and Meta-Analysis

Background: Several recent reports have demonstrated that tyrosine (Y)-methionine (M)-aspartic acid (D)-aspartic acid (D) (YMDD) motif mutations can naturally occur in chronic HBV patients without antiviral treatment such as lamivudine therapy. This paper aims to assess the overall spontaneous incidence and related risk factors of YMDD-motif mutations among lamivudine-naïve chronic HBV carriers, so as to provide some clue for clinical treatment of hepatitis B. Methodology/Principal Findings: Chinese and English literatures were searched for studies reporting natural YMDD mutations among untreated chronic HBV patients from 2001 to 2010. The incidence estimates were summarized and analyzed by meta-analyses. Forty-seven eligible articles from eight countries were selected in this review (13 in English and 34 in Chinese). The pooled incidence of YMDD-motif mutation among untreated chronic HBV patients from eight countries was 12.21 % (95 % CI: 9.69%–14.95%). China had an incidence of 13.38 % (95 % CI: 10.90%–16.07%) and seven other countries had an incidence of 9.90 % (95 % CI: 3.28%–19.55%), respectively. Lamivudine therapy would increase the risk of mutations 5.23 times higher than the untreated patients. A higher HBV DNA copy number was associated with increased incidence of natural YMDD mutation. No significant difference was found in YMDD mutation incidence between groups of different gender, age, HBeAg status, patients ’ ALT (alanine aminotransferase) level, and between the groups of HBV genotype B and C. Conclusions: The YMDD-motif mutations can occur spontaneously with a relatively high incidence in CHB patient

Pivotal Role of Inosine Triphosphate Pyrophosphatase in Maintaining Genome Stability and the Prevention of Apoptosis in Human Cells

Pure nucleotide precursor pools are a prerequisite for high-fidelity DNA replication and the suppression of mutagenesis and carcinogenesis. ITPases are nucleoside triphosphate pyrophosphatases that clean the precursor pools of the non-canonical triphosphates of inosine and xanthine. The precise role of the human ITPase, encoded by the ITPA gene, is not clearly defined. ITPA is clinically important because a widespread polymorphism, 94C>A, leads to null ITPase activity in erythrocytes and is associated with an adverse reaction to thiopurine drugs. We studied the cellular function of ITPA in HeLa cells using the purine analog 6-N hydroxylaminopurine (HAP), whose triphosphate is also a substrate for ITPA. In this study, we demonstrate that ITPA knockdown sensitizes HeLa cells to HAP-induced DNA breaks and apoptosis. The HAP-induced DNA damage and cytotoxicity observed in ITPA knockdown cells are rescued by an overexpression of the yeast ITPase encoded by the HAM1 gene. We further show that ITPA knockdown results in elevated mutagenesis in response to HAP treatment. Our studies reveal the significance of ITPA in preventing base analog-induced apoptosis, DNA damage and mutagenesis in human cells. This implies that individuals with defective ITPase are predisposed to genome damage by impurities in nucleotide pools, which is drastically augmented by therapy with purine analogs. They are also at an elevated risk for degenerative diseases and cancer

Mitogenomic phylogenetic analyses of the Delphinidae with an emphasis on the Globicephalinae

BACKGROUND: Previous DNA-based phylogenetic studies of the Delphinidae family suggest it has undergone rapid diversification, as characterised by unresolved and poorly supported taxonomic relationships (polytomies) for some of the species within this group. Using an increased amount of sequence data we test between alternative hypotheses of soft polytomies caused by rapid speciation, slow evolutionary rate and/or insufficient sequence data, and hard polytomies caused by simultaneous speciation within this family. Combining the mitogenome sequences of five new and 12 previously published species within the Delphinidae, we used Bayesian and maximum-likelihood methods to estimate the phylogeny from partitioned and unpartitioned mitogenome sequences. Further ad hoc tests were then conducted to estimate the support for alternative topologies. RESULTS: We found high support for all the relationships within our reconstructed phylogenies, and topologies were consistent between the Bayesian and maximum-likelihood trees inferred from partitioned and unpartitioned data. Resolved relationships included the placement of the killer whale (Orcinus orca) as sister taxon to the rest of the Globicephalinae subfamily, placement of the Risso's dolphin (Grampus griseus) within the Globicephalinae subfamily, removal of the white-beaked dolphin (Lagenorhynchus albirostris) from the Delphininae subfamily and the placement of the rough-toothed dolphin (Steno bredanensis) as sister taxon to the rest of the Delphininae subfamily rather than within the Globicephalinae subfamily. The additional testing of alternative topologies allowed us to reject all other putative relationships, with the exception that we were unable to reject the hypothesis that the relationship between L. albirostris and the Globicephalinae and Delphininae subfamilies was polytomic. CONCLUSION: Despite their rapid diversification, the increased sequence data yielded by mitogenomes enables the resolution of a strongly supported, bifurcating phylogeny, and a chronology of the divergences within the Delphinidae family. This highlights the benefits and potential application of large mitogenome datasets to resolve long-standing phylogenetic uncertainties

Genome-wide association analysis of insomnia complaints identifies risk genes and genetic overlap with psychiatric and metabolic traits.

To access publisher's full text version of this article, please click on the hyperlink in Additional Links field or click on the hyperlink at the top of the page marked FilesPersistent insomnia is among the most frequent complaints in general practice. To identify genetic factors for insomnia complaints, we performed a genome-wide association study (GWAS) and a genome-wide gene-based association study (GWGAS) in 113,006 individuals. We identify three loci and seven genes associated with insomnia complaints, with the associations for one locus and five genes supported by joint analysis with an independent sample (n = 7,565). Our top association (MEIS1, P < 5 × 10-8) has previously been implicated in restless legs syndrome (RLS). Additional analyses favor the hypothesis that MEIS1 exhibits pleiotropy for insomnia and RLS and show that the observed association with insomnia complaints cannot be explained only by the presence of an RLS subgroup within the cases. Sex-specific analyses suggest that there are different genetic architectures between the sexes in addition to shared genetic factors. We show substantial positive genetic correlation of insomnia complaints with internalizing personality traits and metabolic traits and negative correlation with subjective well-being and educational attainment. These findings provide new insight into the genetic architecture of insomnia.Netherlands Organization for Scientific Research NWO Brain & Cognition 433-09-228 European Research Council ERC-ADG-2014-671084 INSOMNIA Netherlands Scientific Organization (NWO) VU University (Amsterdam, the Netherlands) Dutch Brain Foundation Helmholtz Zentrum Munchen - German Federal Ministry of Education and Research state of Bavaria German Migraine & Headache Society (DMKG) Almirall AstraZeneca Berlin Chemie Boehringer Boots Health Care GlaxoSmithKline Janssen Cilag McNeil Pharma MSD Sharp Dohme Pfizer Institute of Epidemiology and Social Medicine at the University of Munster German Ministry of Education and Research (BMBF) German Restless Legs Patient Organisation (RLS Deutsche Restless Legs Vereinigung) Swiss RLS Patient Association (Schweizerische Restless Legs Selbsthilfegruppe

Epigenetic inactivation of TCF2 in ovarian cancer and various cancer cell lines

Transcription factor 2 gene (TCF2) encodes hepatocyte nuclear factor 1β (HNF1β), a transcription factor associated with development and metabolism. Mutation of TCF2 has been observed in renal cell cancer, and by screening aberrantly methylated genes, we have now identified TCF2 as a target for epigenetic inactivation in ovarian cancer. TCF2 was methylated in 53% of ovarian cancer cell lines and 26% of primary ovarian cancers, resulting in loss of the gene's expression. TCF2 expression was restored by treating cells with a methyltransferase inhibitor, 5-aza-2′deoxycitidine (5-aza-dC). In addition, chromatin immunoprecipitation showed deacetylation of histone H3 in methylated cells and, when combined with 5-aza-dC, the histone deacetylase inhibitor trichostatin A synergistically induced TCF2 expression. Epigenetic inactivation of TCF2 was also seen in colorectal, gastric and pancreatic cell lines, suggesting general involvement of epigenetic inactivation of TCF2 in tumorigenesis. Restoration of TCF2 expression induced expression of HNF4α, a transcriptional target of HNF1β, indicating that epigenetic silencing of TCF2 leads to alteration of the hepatocyte nuclear factor network in tumours. These results suggest that TCF2 is involved in the development of ovarian cancers and may represent a useful target for their detection and treatment

Developmental Patterns of Doublecortin Expression and White Matter Neuron Density in the Postnatal Primate Prefrontal Cortex and Schizophrenia

Postnatal neurogenesis occurs in the subventricular zone and dentate gyrus, and evidence suggests that new neurons may be present in additional regions of the mature primate brain, including the prefrontal cortex (PFC). Addition of new neurons to the PFC implies local generation of neurons or migration from areas such as the subventricular zone. We examined the putative contribution of new, migrating neurons to postnatal cortical development by determining the density of neurons in white matter subjacent to the cortex and measuring expression of doublecortin (DCX), a microtubule-associated protein involved in neuronal migration, in humans and rhesus macaques. We found a striking decline in DCX expression (human and macaque) and density of white matter neurons (humans) during infancy, consistent with the arrival of new neurons in the early postnatal cortex. Considering the expansion of the brain during this time, the decline in white matter neuron density does not necessarily indicate reduced total numbers of white matter neurons in early postnatal life. Furthermore, numerous cells in the white matter and deep grey matter were positive for the migration-associated glycoprotein polysialiated-neuronal cell adhesion molecule and GAD65/67, suggesting that immature migrating neurons in the adult may be GABAergic. We also examined DCX mRNA in the PFC of adult schizophrenia patients (n = 37) and matched controls (n = 37) and did not find any difference in DCX mRNA expression. However, we report a negative correlation between DCX mRNA expression and white matter neuron density in adult schizophrenia patients, in contrast to a positive correlation in human development where DCX mRNA and white matter neuron density are higher earlier in life. Accumulation of neurons in the white matter in schizophrenia would be congruent with a negative correlation between DCX mRNA and white matter neuron density and support the hypothesis of a migration deficit in schizophrenia