DNA replication stress restricts ribosomal DNA copy number

Ribosomal RNAs (rRNAs) in budding yeast are encoded by ~100–200 repeats of a 9.1kb sequence arranged in tandem on chromosome XII, the ribosomal DNA (rDNA) locus. Copy number of rDNA repeat units in eukaryotic cells is maintained far in excess of the requirement for ribosome biogenesis. Despite the importance of the repeats for both ribosomal and non-ribosomal functions, it is currently not known how “normal” copy number is determined or maintained. To identify essential genes involved in the maintenance of rDNA copy number, we developed a droplet digital PCR based assay to measure rDNA copy number in yeast and used it to screen a yeast conditional temperature-sensitive mutant collection of essential genes. Our screen revealed that low rDNA copy number is associated with compromised DNA replication. Further, subculturing yeast under two separate conditions of DNA replication stress selected for a contraction of the rDNA array independent of the replication fork blocking protein, Fob1. Interestingly, cells with a contracted array grew better than their counterparts with normal copy number under conditions of DNA replication stress. Our data indicate that DNA replication stresses select for a smaller rDNA array. We speculate that this liberates scarce replication factors for use by the rest of the genome, which in turn helps cells complete DNA replication and continue to propagate. Interestingly, tumors from mini chromosome maintenance 2 (MCM2)-deficient mice also show a loss of rDNA repeats. Our data suggest that a reduction in rDNA copy number may indicate a history of DNA replication stress, and that rDNA array size could serve as a diagnostic marker for replication stress. Taken together, these data begin to suggest the selective pressures that combine to yield a “normal” rDNA copy number

What is the correct cost functional for variational data assimilation?

Variational approaches to data assimilation, and weakly constrained four dimensional variation (WC-4DVar) in particular, are important in the geosciences but also in other communities (often under different names). The cost functions and the resulting optimal trajectories may have a probabilistic interpretation, for instance by linking data assimilation with maximum aposteriori (MAP) estimation. This is possible in particular if the unknown trajectory is modelled as the solution of a stochastic differential equation (SDE), as is increasingly the case in weather forecasting and climate modelling. In this situation, the MAP estimator (or “most probable path” of the SDE) is obtained by minimising the Onsager–Machlup functional. Although this fact is well known, there seems to be some confusion in the literature, with the energy (or “least squares”) functional sometimes been claimed to yield the most probable path. The first aim of this paper is to address this confusion and show that the energy functional does not, in general, provide the most probable path. The second aim is to discuss the implications in practice. Although the mentioned results pertain to stochastic models in continuous time, they do have consequences in practice where SDE’s are approximated by discrete time schemes. It turns out that using an approximation to the SDE and calculating its most probable path does not necessarily yield a good approximation to the most probable path of the SDE proper. This suggest that even in discrete time, a version of the Onsager–Machlup functional should be used, rather than the energy functional, at least if the solution is to be interpreted as a MAP estimator

Non-contrast cardiac computed tomography can accurately detect chronic myocardial infarction: Validation study

BackgroundThis study evaluates whether non-contrast cardiac computed tomography (CCT) can detect chronic myocardial infarction (MI) in patients with irreversible perfusion defects on nuclear myocardial perfusion imaging (MPI).MethodsOne hundred twenty-two symptomatic patients with irreversible perfusion defect (N = 62) or normal MPI (N = 60) underwent coronary artery calcium (CAC) scanning. MI on these non-contrast CCTs was visually detected based on the hypo-attenuation areas (dark) in the myocardium and corresponding Hounsfield units (HU) were measured.ResultsNon-contrast CCT accurately detected MI in 57 patients with irreversible perfusion defect on MPI, yielding a sensitivity of 92%, specificity of 72%, negative predictive value (NPV) of 90%, and a positive predictive value (PPV) of 77%. On a per myocardial region analysis, non-contrast CT showed a sensitivity of 70%, specificity of 85%, NPV of 91%, and a PPV of 57%. The ROC curve showed that the optimal cutoff value of LV myocardium HU to predict MI on non-contrast CCT was 21.7 with a sensitivity of 97.4% and specificity of 99.7%.ConclusionNon-contrast CCT has an excellent agreement with MPI in detecting chronic MI. This study highlights a novel clinical utility of non-contrast CCT in addition to assessment of overall burden of atherosclerosis measured by CAC

Cell-specific deletion of C1qa identifies microglia as the dominant source of C1q in mouse brain

BACKGROUND: The complement cascade not only provides protection from infection but can also mediate destructive inflammation. Complement is also involved in elimination of neuronal synapses which is essential for proper development, but can be detrimental during aging and disease. C1q, required for several of these complement-mediated activities, is present in the neuropil, microglia, and a subset of interneurons in the brain. METHODS: To identify the source(s) of C1q in the brain, the C1qa gene was selectively inactivated in the microglia or Thy-1(+) neurons in both wild type mice and a mouse model of Alzheimer’s disease (AD), and C1q synthesis assessed by immunohistochemistry, QPCR, and western blot analysis. RESULTS: While C1q expression in the brain was unaffected after inactivation of C1qa in Thy-1(+) neurons, the brains of C1qa (FL/FL) :Cx3cr1 (CreERT2) mice in which C1qa was ablated in microglia were devoid of C1q with the exception of limited C1q in subsets of interneurons. Surprisingly, this loss of C1q occurred even in the absence of tamoxifen by 1 month of age, demonstrating that Cre activity is tamoxifen-independent in microglia in Cx3cr1 (CreERT2/WganJ) mice. C1q expression in C1qa (FL/FL) : Cx3cr1 (CreERT2/WganJ) mice continued to decline and remained almost completely absent through aging and in AD model mice. No difference in C1q was detected in the liver or kidney from C1qa (FL/FL) : Cx3cr1 (CreERT2/WganJ) mice relative to controls, and C1qa (FL/FL) : Cx3cr1 (CreERT2/WganJ) mice had minimal, if any, reduction in plasma C1q. CONCLUSIONS: Thus, microglia, but not neurons or peripheral sources, are the dominant source of C1q in the brain. While demonstrating that the Cx3cr1 (CreERT2/WganJ) deleter cannot be used for adult-induced deletion of genes in microglia, the model described here enables further investigation of physiological roles of C1q in the brain and identification of therapeutic targets for the selective control of complement-mediated activities contributing to neurodegenerative disorders. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12974-017-0814-9) contains supplementary material, which is available to authorized users

Association study of functional genetic variants of innate immunity related genes in celiac disease

BACKGROUND: Recent evidence suggest that the innate immune system is implicated in the early events of celiac disease (CD) pathogenesis. In this work for the first time we have assessed the relevance of different proinflammatory mediators typically related to innate immunity in CD predisposition. METHODS: We performed a familial study in which 105 celiac families characterized by the presence of an affected child with CD were genotyped for functional polymorphisms located at regulatory regions of IL-1α, IL-1β, IL-1RN, IL-18, RANTES and MCP-1 genes. Familial data was analysed with a transmission disequilibrium test (TDT) that revealed no statistically significant differences in the transmission pattern of the different genetic markers considered. RESULTS: The TDT analysis for IL-1α, IL-1β, IL-1RN, IL-18, and MCP-1 genes genetic variants did not reveal biased transmission to the affected offspring. Only a borderline association of RANTES promoter genetic variants with CD predisposition was observed. CONCLUSION: Our results suggest that the analysed polymorphisms of IL-1α, IL-1β, IL-1RN, IL-18, RANTES and MCP-1 genes do not seem to play a major role in CD genetic predisposition in our population

Ribosomal DNA Deletions Modulate Genome-Wide Gene Expression: “rDNA–Sensitive” Genes and Natural Variation

The ribosomal rDNA gene array is an epigenetically-regulated repeated gene locus. While rDNA copy number varies widely between and within species, the functional consequences of subtle copy number polymorphisms have been largely unknown. Deletions in the Drosophila Y-linked rDNA modifies heterochromatin-induced position effect variegation (PEV), but it has been unknown if the euchromatic component of the genome is affected by rDNA copy number. Polymorphisms of naturally occurring Y chromosomes affect both euchromatin and heterochromatin, although the elements responsible for these effects are unknown. Here we show that copy number of the Y-linked rDNA array is a source of genome-wide variation in gene expression. Induced deletions in the rDNA affect the expression of hundreds to thousands of euchromatic genes throughout the genome of males and females. Although the affected genes are not physically clustered, we observed functional enrichments for genes whose protein products are located in the mitochondria and are involved in electron transport. The affected genes significantly overlap with genes affected by natural polymorphisms on Y chromosomes, suggesting that polymorphic rDNA copy number is an important determinant of gene expression diversity in natural populations. Altogether, our results indicate that subtle changes to rDNA copy number between individuals may contribute to biologically relevant phenotypic variation

Quantitative utilization of prior biological knowledge in the Bayesian network modeling of gene expression data

<p>Abstract</p> <p>Background</p> <p>Bayesian Network (BN) is a powerful approach to reconstructing genetic regulatory networks from gene expression data. However, expression data by itself suffers from high noise and lack of power. Incorporating prior biological knowledge can improve the performance. As each type of prior knowledge on its own may be incomplete or limited by quality issues, integrating multiple sources of prior knowledge to utilize their consensus is desirable.</p> <p>Results</p> <p>We introduce a new method to incorporate the quantitative information from multiple sources of prior knowledge. It first uses the Naïve Bayesian classifier to assess the likelihood of functional linkage between gene pairs based on prior knowledge. In this study we included cocitation in PubMed and schematic similarity in Gene Ontology annotation. A candidate network edge reservoir is then created in which the copy number of each edge is proportional to the estimated likelihood of linkage between the two corresponding genes. In network simulation the Markov Chain Monte Carlo sampling algorithm is adopted, and samples from this reservoir at each iteration to generate new candidate networks. We evaluated the new algorithm using both simulated and real gene expression data including that from a yeast cell cycle and a mouse pancreas development/growth study. Incorporating prior knowledge led to a ~2 fold increase in the number of known transcription regulations recovered, without significant change in false positive rate. In contrast, without the prior knowledge BN modeling is not always better than a random selection, demonstrating the necessity in network modeling to supplement the gene expression data with additional information.</p> <p>Conclusion</p> <p>our new development provides a statistical means to utilize the quantitative information in prior biological knowledge in the BN modeling of gene expression data, which significantly improves the performance.</p

Large Scale Gene Expression Profiles of Regenerating Inner Ear Sensory Epithelia

Loss of inner ear sensory hair cells (HC) is a leading cause of human hearing loss and balance disorders. Unlike mammals, many lower vertebrates can regenerate these cells. We used cross-species microarrays to examine this process in the avian inner ear. Specifically, changes in expression of over 1700 transcription factor (TF) genes were investigated in hair cells of auditory and vestibular organs following treatment with two different damaging agents and regeneration in vitro. Multiple components of seven distinct known signaling pathways were clearly identifiable: TGFβ, PAX, NOTCH, WNT, NFKappaB, INSULIN/IGF1 and AP1. Numerous components of apoptotic and cell cycle control pathways were differentially expressed, including p27KIP and TFs that regulate its expression. A comparison of expression trends across tissues and treatments revealed identical patterns of expression that occurred at identical times during regenerative proliferation. Network analysis of the patterns of gene expression in this large dataset also revealed the additional presence of many components (and possible network interactions) of estrogen receptor signaling, circadian rhythm genes and parts of the polycomb complex (among others). Equal numbers of differentially expressed genes were identified that have not yet been placed into any known pathway. Specific time points and tissues also exhibited interesting differences: For example, 45 zinc finger genes were specifically up-regulated at later stages of cochlear regeneration. These results are the first of their kind and should provide the starting point for more detailed investigations of the role of these many pathways in HC recovery, and for a description of their possible interactions

Physical mapping and BAC-end sequence analysis provide initial insights into the flax (Linum usitatissimum L.) genome

<p>Abstract</p> <p>Background</p> <p>Flax (<it>Linum usitatissimum </it>L.) is an important source of oil rich in omega-3 fatty acids, which have proven health benefits and utility as an industrial raw material. Flax seeds also contain lignans which are associated with reducing the risk of certain types of cancer. Its bast fibres have broad industrial applications. However, genomic tools needed for molecular breeding were non existent. Hence a project, Total Utilization Flax GENomics (TUFGEN) was initiated. We report here the first genome-wide physical map of flax and the generation and analysis of BAC-end sequences (BES) from 43,776 clones, providing initial insights into the genome.</p> <p>Results</p> <p>The physical map consists of 416 contigs spanning ~368 Mb, assembled from 32,025 fingerprints, representing roughly 54.5% to 99.4% of the estimated haploid genome (370-675 Mb). The N50 size of the contigs was estimated to be ~1,494 kb. The longest contig was ~5,562 kb comprising 437 clones. There were 96 contigs containing more than 100 clones. Approximately 54.6 Mb representing 8-14.8% of the genome was obtained from 80,337 BES. Annotation revealed that a large part of the genome consists of ribosomal DNA (~13.8%), followed by known transposable elements at 6.1%. Furthermore, ~7.4% of sequence was identified to harbour novel repeat elements. Homology searches against flax-ESTs and NCBI-ESTs suggested that ~5.6% of the transcriptome is unique to flax. A total of 4064 putative genomic SSRs were identified and are being developed as novel markers for their use in molecular breeding.</p> <p>Conclusion</p> <p>The first genome-wide physical map of flax constructed with BAC clones provides a framework for accessing target loci with economic importance for marker development and positional cloning. Analysis of the BES has provided insights into the uniqueness of the flax genome. Compared to other plant genomes, the proportion of rDNA was found to be very high whereas the proportion of known transposable elements was low. The SSRs identified from BES will be valuable in saturating existing linkage maps and for anchoring physical and genetic maps. The physical map and paired-end reads from BAC clones will also serve as scaffolds to build and validate the whole genome shotgun assembly.</p