An engineered multicellular stem cell niche for the 3D derivation of human myogenic progenitors from iPSCs.

Fate decisions in the embryo are controlled by a plethora of microenvironmental interactions in a three-dimensional niche. To investigate whether aspects of this microenvironmental complexity can be engineered to direct myogenic human-induced pluripotent stem cell (hiPSC) differentiation, we here screened murine cell types present in the developmental or adult stem cell niche in heterotypic suspension embryoids. We identified embryonic endothelial cells and fibroblasts as highly permissive for myogenic specification of hiPSCs. After two weeks of sequential Wnt and FGF pathway induction, these three-component embryoids are enriched in Pax7-positive embryonic-like myogenic progenitors that can be isolated by flow cytometry. Myogenic differentiation of hiPSCs in heterotypic embryoids relies on a specialized structural microenvironment and depends on MAPK, PI3K/AKT, and Notch signaling. After transplantation in a mouse model of Duchenne muscular dystrophy, embryonic-like myogenic progenitors repopulate the stem cell niche, reactivate after repeated injury, and, compared to adult human myoblasts, display enhanced fusion and lead to increased muscle function. Altogether, we provide a two-week protocol for efficient and scalable suspension-based 3D derivation of Pax7-positive myogenic progenitors from hiPSCs.

Statistical analysis of high-dimensional biomedical data: a gentle introduction to analytical goals, common approaches and challenges.

BACKGROUND: In high-dimensional data (HDD) settings, the number of variables associated with each observation is very large. Prominent examples of HDD in biomedical research include omics data with a large number of variables such as many measurements across the genome, proteome, or metabolome, as well as electronic health records data that have large numbers of variables recorded for each patient. The statistical analysis of such data requires knowledge and experience, sometimes of complex methods adapted to the respective research questions. METHODS: Advances in statistical methodology and machine learning methods offer new opportunities for innovative analyses of HDD, but at the same time require a deeper understanding of some fundamental statistical concepts. Topic group TG9 "High-dimensional data" of the STRATOS (STRengthening Analytical Thinking for Observational Studies) initiative provides guidance for the analysis of observational studies, addressing particular statistical challenges and opportunities for the analysis of studies involving HDD. In this overview, we discuss key aspects of HDD analysis to provide a gentle introduction for non-statisticians and for classically trained statisticians with little experience specific to HDD. RESULTS: The paper is organized with respect to subtopics that are most relevant for the analysis of HDD, in particular initial data analysis, exploratory data analysis, multiple testing, and prediction. For each subtopic, main analytical goals in HDD settings are outlined. For each of these goals, basic explanations for some commonly used analysis methods are provided. Situations are identified where traditional statistical methods cannot, or should not, be used in the HDD setting, or where adequate analytic tools are still lacking. Many key references are provided. CONCLUSIONS: This review aims to provide a solid statistical foundation for researchers, including statisticians and non-statisticians, who are new to research with HDD or simply want to better evaluate and understand the results of HDD analyses.

Association of amino acid metabolites with osteoporosis, a metabolomic approach: Bushehr elderly health program.

INTRODUCTION AND OBJECTIVES: Amino acids are the most frequently reported metabolites associated with low bone mineral density (BMD) in metabolomics studies. We aimed to evaluate the association between amino acid metabolic profile and bone indices in the elderly population. METHODS: 400 individuals were randomly selected from 2384 elderly men and women over 60 years participating in the second stage of the Bushehr elderly health (BEH) program, a population-based prospective cohort study that is being conducted in Bushehr, a southern province of Iran. Frozen plasma samples were used to measure 29 amino acid and derivatives metabolites using the UPLC-MS/MS-based targeted metabolomics platform. We conducted Elastic net regression analysis to detect the metabolites associated with BMD of different sites and lumbar spine trabecular bone score, and also to examine the ability of the measured metabolites to differentiate osteoporosis. RESULTS: We adjusted the analysis for possible confounders (age, BMI, diabetes, smoking, physical activity, vitamin D level, and sex). Valine, leucine, isoleucine, and alanine in women and tryptophan in men were the most important amino acids inversely associated with osteoporosis (OR range from 0.77 to 0.89). Sarcosine, followed by tyrosine, asparagine, alpha aminobutyric acid, and ADMA in women and glutamine in men and when both women and men were considered together were the most discriminating amino acids detected in individuals with osteoporosis (OR range from 1.15 to 1.31). CONCLUSION: We found several amino acid metabolites associated with possible bone status in elderly individuals. Further studies are required to evaluate the utility of these metabolites as clinical biomarkers for osteoporosis prediction and their effect on bone health as dietary supplements.

Diurnal regulation of RNA polymerase III transcription is under the control of both the feeding-fasting response and the circadian clock.

RNA polymerase III (Pol III) synthesizes short noncoding RNAs, many of which are essential for translation. Accordingly, Pol III activity is tightly regulated with cell growth and proliferation by factors such as MYC, RB1, TRP53, and MAF1. MAF1 is a repressor of Pol III transcription whose activity is controlled by phosphorylation; in particular, it is inactivated through phosphorylation by the TORC1 kinase complex, a sensor of nutrient availability. Pol III regulation is thus sensitive to environmental cues, yet a diurnal profile of Pol III transcription activity is so far lacking. Here, we first use gene expression arrays to measure mRNA accumulation during the diurnal cycle in the livers of (1) wild-type mice, (2) arrhythmic Arntl knockout mice, (3) mice fed at regular intervals during both night and day, and (4) mice lacking the Maf1 gene, and so provide a comprehensive view of the changes in cyclic mRNA accumulation occurring in these different systems. We then show that Pol III occupancy of its target genes rises before the onset of the night, stays high during the night, when mice normally ingest food and when translation is known to be increased, and decreases in daytime. Whereas higher Pol III occupancy during the night reflects a MAF1-dependent response to feeding, the rise of Pol III occupancy before the onset of the night reflects a circadian clock-dependent response. Thus, Pol III transcription during the diurnal cycle is regulated both in response to nutrients and by the circadian clock, which allows anticipatory Pol III transcription.

Apolipoprotein E4 Expression Causes Gain of Toxic Function in Isogenic Human Induced Pluripotent Stem Cell-Derived Endothelial Cells.

OBJECTIVE: The ApoE (apolipoprotein) allele epsilon 4 is a major genetic risk factor for Alzheimer disease, cardiovascular disorders, and stroke, indicating that it significantly impacts cerebral and vascular systems. However, very little is known about how APOE genotype affects brain endothelial cells, which form a network of tight junctions to regulate communication between the brain and circulating blood factors. Approach and Results: Here, we present a novel model of endothelial dysfunction using isogenic human induced pluripotent stem cell-derived cells harboring different alleles of the APOE gene, specifically ApoE 3/3, 3/4, and 4/4. We show for the first time that ApoE4 expression by endothelial cells is sufficient to cause a toxic gain of cellular dysfunction. Using RNAseq, we found significant effects of ApoE4 on signaling pathways involved in blood coagulation and barrier function. These changes were associated with altered cell function, including increased binding of platelets to ECs with the 3/4 or 4/4 genotype. ApoE4-positive cells exhibited a proinflammatory state and prothrombotic state, evidenced by higher secretion of Aß (amyloid-ß) 40 and 42, increased release of cytokines, and overexpression of the platelet-binding protein VWF (vonWillebrand factor). Immunohistochemistry of human brain Alzheimer disease brains also showed increased VWF expression with the ApoE4/4 genotype. Finally, pharmacological inhibition of inflammation in ECs by celastrol rescued overexpression of VWF in cells expressing ApoE4. CONCLUSIONS: These cells provide novel insight into ApoE4-mediated endothelial dysfunction and provide a new platform to test potential therapies for vascular disorders.

Domains of genome-wide gene expression dysregulation in Down's syndrome.

Trisomy 21 is the most frequent genetic cause of cognitive impairment. To assess the perturbations of gene expression in trisomy 21, and to eliminate the noise of genomic variability, we studied the transcriptome of fetal fibroblasts from a pair of monozygotic twins discordant for trisomy 21. Here we show that the differential expression between the twins is organized in domains along all chromosomes that are either upregulated or downregulated. These gene expression dysregulation domains (GEDDs) can be defined by the expression level of their gene content, and are well conserved in induced pluripotent stem cells derived from the twins' fibroblasts. Comparison of the transcriptome of the Ts65Dn mouse model of Down's syndrome and normal littermate mouse fibroblasts also showed GEDDs along the mouse chromosomes that were syntenic in human. The GEDDs correlate with the lamina-associated (LADs) and replication domains of mammalian cells. The overall position of LADs was not altered in trisomic cells; however, the H3K4me3 profile of the trisomic fibroblasts was modified and accurately followed the GEDD pattern. These results indicate that the nuclear compartments of trisomic cells undergo modifications of the chromatin environment influencing the overall transcriptome, and that GEDDs may therefore contribute to some trisomy 21 phenotypes.

Proteomes of Paired Human Cerebrospinal Fluid and Plasma: Relation to Blood-Brain Barrier Permeability in Older Adults.

The systems-level relationship between the proteomes of cerebrospinal fluid (CSF) and plasma has not been comprehensively described so far. Recently developed shotgun proteomic workflows allow for deeper characterization of the proteomes from body fluids in much larger sample size. We deployed state-of-the-art mass spectrometry-based proteomics in paired CSF and plasma samples volunteered by 120 elders with and without cognitive impairment to comprehensively characterize and examine compartmental proteome differences and relationships between both body fluids. We further assessed the influence of blood-brain barrier (BBB) integrity and tested the hypothesis that BBB breakdown can be identified from CSF and plasma proteome alterations in nondemented elders. We quantified 790 proteins in CSF and 422 proteins in plasma, and 255 of the proteins were identified in both compartments. Pearson's statistics determined 28 proteins with associated levels between CSF and plasma. BBB integrity as defined with the CSF/serum albumin index influenced 76 CSF/plasma protein ratios. In least absolute shrinkage and selection operator models, CSF and plasma proteins improved identification of BBB impairment. In conclusion, we provide here a first comprehensive draft map of interacting human CSF and plasma proteomes, in view of their complex and dynamic compositions, and influence of the BBB.

Reciprocal Effects on Neurocognitive and Metabolic Phenotypes in Mouse Models of 16p11.2 Deletion and Duplication Syndromes.

The 16p11.2 600 kb BP4-BP5 deletion and duplication syndromes have been associated with developmental delay; autism spectrum disorders; and reciprocal effects on the body mass index, head circumference and brain volumes. Here, we explored these relationships using novel engineered mouse models carrying a deletion (Del/+) or a duplication (Dup/+) of the Sult1a1-Spn region homologous to the human 16p11.2 BP4-BP5 locus. On a C57BL/6N inbred genetic background, Del/+ mice exhibited reduced weight and impaired adipogenesis, hyperactivity, repetitive behaviors, and recognition memory deficits. In contrast, Dup/+ mice showed largely opposite phenotypes. On a F1 C57BL/6N × C3B hybrid genetic background, we also observed alterations in social interaction in the Del/+ and the Dup/+ animals, with other robust phenotypes affecting recognition memory and weight. To explore the dosage effect of the 16p11.2 genes on metabolism, Del/+ and Dup/+ models were challenged with high fat and high sugar diet, which revealed opposite energy imbalance. Transcriptomic analysis revealed that the majority of the genes located in the Sult1a1-Spn region were sensitive to dosage with a major effect on several pathways associated with neurocognitive and metabolic phenotypes. Whereas the behavioral consequence of the 16p11 region genetic dosage was similar in mice and humans with activity and memory alterations, the metabolic defects were opposite: adult Del/+ mice are lean in comparison to the human obese phenotype and the Dup/+ mice are overweight in comparison to the human underweight phenotype. Together, these data indicate that the dosage imbalance at the 16p11.2 locus perturbs the expression of modifiers outside the CNV that can modulate the penetrance, expressivity and direction of effects in both humans and mice.

A Potential Contributory Role for Ciliary Dysfunction in the 16p11.2 600 kb BP4-BP5 Pathology.

The 16p11.2 600 kb copy-number variants (CNVs) are associated with mirror phenotypes on BMI, head circumference, and brain volume and represent frequent genetic lesions in autism spectrum disorders (ASDs) and schizophrenia. Here we interrogated the transcriptome of individuals carrying reciprocal 16p11.2 CNVs. Transcript perturbations correlated with clinical endophenotypes and were enriched for genes associated with ASDs, abnormalities of head size, and ciliopathies. Ciliary gene expression was also perturbed in orthologous mouse models, raising the possibility that ciliary dysfunction contributes to 16p11.2 pathologies. In support of this hypothesis, we found structural ciliary defects in the CA1 hippocampal region of 16p11.2 duplication mice. Moreover, by using an established zebrafish model, we show genetic interaction between KCTD13, a key driver of the mirrored neuroanatomical phenotypes of the 16p11.2 CNV, and ciliopathy-associated genes. Overexpression of BBS7 rescues head size and neuroanatomical defects of kctd13 morphants, whereas suppression or overexpression of CEP290 rescues phenotypes induced by KCTD13 under- or overexpression, respectively. Our data suggest that dysregulation of ciliopathy genes contributes to the clinical phenotypes of these CNVs.

Blood-brain barrier breakdown, neuroinflammation, and cognitive decline in older adults.

INTRODUCTION: Blood-brain barrier (BBB) breakdown is observed in older versus younger adults and in late-onset Alzheimer's disease versus age-matched controls, but its causes and consequences in aging are unclear. We tested the hypothesis that BBB breakdown is associated with cognitive decline and inflammation in nondemented elders. METHODS: Cerebrospinal fluid and serum inflammatory markers were measured using sandwich immunoassays in 120 subjects. Least Absolute Shrinkage and Selection Operator-logistic regression selected cerebrospinal fluid and serum signatures that best classified BBB impairment defined by the cerebrospinal fluid albumin index ≥9. Linear regression examined changes in Clinical Dementia Rating sum of boxes as a function of BBB integrity at baseline. RESULTS: Mean age was 70 years, mean Mini­Mental State Examination was 27, and BBB impairment was recorded in 13.5%. BBB breakdown was associated with cognitive decline (P = .015). Cerebrospinal fluid intercellular adhesion molecule-1, vascular endothelial growth factor, interleukin-8, serum amyloid A, macrophage derived chemokine, and gender generated an area under the curve of 0.95 for BBB impairment, and serum IL-16, VEGF-D, IL-15, and other variables generated an AUC of 0.92 for BBB impairment. DISCUSSION: BBB breakdown is associated with more rapid cognitive decline. Inflammatory mechanisms, including cell adhesion, neutrophil migration, lipid metabolism, and angiogenesis may be implicated.

PPARγ Controls Ectopic Adipogenesis and Cross-Talks with Myogenesis During Skeletal Muscle Regeneration.

Skeletal muscle is a regenerative tissue which can repair damaged myofibers through the activation of tissue-resident muscle stem cells (MuSCs). Many muscle diseases with impaired regeneration cause excessive adipose tissue accumulation in muscle, alter the myogenic fate of MuSCs, and deregulate the cross-talk between MuSCs and fibro/adipogenic progenitors (FAPs), a bi-potent cell population which supports myogenesis and controls intra-muscular fibrosis and adipocyte formation. In order to better characterize the interaction between adipogenesis and myogenesis, we studied muscle regeneration and MuSC function in whole body Pparg null mice generated by epiblast-specific Cre/lox deletion (PpargΔ/Δ). We demonstrate that deletion of PPARγ completely abolishes ectopic muscle adipogenesis during regeneration and impairs MuSC expansion and myogenesis after injury. Ex vivo assays revealed that perturbed myogenesis in PpargΔ/Δ mice does not primarily result from intrinsic defects of MuSCs or from perturbed myogenic support from FAPs. The immune transition from a pro- to anti-inflammatory MuSC niche during regeneration is perturbed in PpargΔ/Δ mice and suggests that PPARγ signaling in macrophages can interact with ectopic adipogenesis and influence muscle regeneration. Altogether, our study demonstrates that a PPARγ-dependent adipogenic response regulates muscle fat infiltration during regeneration and that PPARγ is required for MuSC function and efficient muscle repair.

ICP-MS/MS-Based Ionomics: A Validated Methodology to Investigate the Biological Variability of the Human Ionome.

We here describe the development, validation and application of a quantitative methodology for the simultaneous determination of 29 elements in human serum using state-of-the-art inductively coupled plasma triple quadrupole mass spectrometry (ICP-MS/MS). This new methodology offers high-throughput elemental profiling using simple dilution of minimal quantity of serum samples. We report the outcomes of the validation procedure including limits of detection/quantification, linearity of calibration curves, precision, recovery and measurement uncertainty. ICP-MS/MS-based ionomics was used to analyze human serum of 120 older adults. Following a metabolomic data mining approach, the generated ionome profiles were subjected to principal component analysis revealing gender and age-specific differences. The ionome of female individuals was marked by higher levels of calcium, phosphorus, copper and copper to zinc ratio, while iron concentration was lower with respect to male subjects. Age was associated with lower concentrations of zinc. These findings were complemented with additional readouts to interpret micronutrient status including ceruloplasmin, ferritin and inorganic phosphate. Our data supports a gender-specific compartmentalization of the ionome that may reflect different bone remodelling in female individuals. Our ICP-MS/MS methodology enriches the panel of validated "Omics" approaches to study molecular relationships between the exposome and the ionome in relation with nutrition and health.

Novel H3K4me3 marks are enriched at human- and chimpanzee-specific cytogenetic structures.

Human and chimpanzee genomes are 98.8% identical within comparable sequences. However, they differ structurally in nine pericentric inversions, one fusion that originated human chromosome 2, and content and localization of heterochromatin and lineage-specific segmental duplications. The possible functional consequences of these cytogenetic and structural differences are not fully understood and their possible involvement in speciation remains unclear. We show that subtelomeric regions--regions that have a species-specific organization, are more divergent in sequence, and are enriched in genes and recombination hotspots--are significantly enriched for species-specific histone modifications that decorate transcription start sites in different tissues in both human and chimpanzee. The human lineage-specific chromosome 2 fusion point and ancestral centromere locus as well as chromosome 1 and 18 pericentric inversion breakpoints showed enrichment of human-specific H3K4me3 peaks in the prefrontal cortex. Our results reveal an association between plastic regions and potential novel regulatory elements.

Markers of neuroinflammation associated with Alzheimer's disease pathology in older adults.

BACKGROUND: In vitro and animal studies have linked neuroinflammation to Alzheimer's disease (AD) pathology. Studies on markers of inflammation in subjects with mild cognitive impairment or AD dementia provided inconsistent results. We hypothesized that distinct blood and cerebrospinal fluid (CSF) inflammatory markers are associated with biomarkers of amyloid and tau pathology in older adults without cognitive impairment or with beginning cognitive decline. OBJECTIVE: To identify blood-based and CSF neuroinflammation marker signatures associated with AD pathology (i.e. an AD CSF biomarker profile) and to investigate associations of inflammation markers with CSF biomarkers of amyloid, tau pathology, and neuronal injury. DESIGN/METHODS: Cross-sectional analysis was performed on data from 120 older community-dwelling adults with normal cognition (n=48) or with cognitive impairment (n=72). CSF Aß1-42, tau and p-tau181, and a panel of 37 neuroinflammatory markers in both CSF and serum were quantified. Least absolute shrinkage and selection operator (LASSO) regression was applied to determine a reference model that best predicts an AD CSF biomarker profile defined a priori as p-tau181/Aß1-42 ratio >0.0779. It was then compared to a second model that included the inflammatory markers from either serum or CSF. In addition, the correlations between inflammatory markers and CSF Aß1-42, tau and p-tau181 levels were assessed. RESULTS: Forty-two subjects met criteria for having an AD CSF biomarker profile. The best predictive models included 8 serum or 3 CSF neuroinflammatory markers related to cytokine mediated inflammation, vascular injury, and angiogenesis. Both models improved the accuracy to predict an AD biomarker profile when compared to the reference model. In analyses separately performed in the subgroup of participants with cognitive impairment, adding the serum or the CSF neuroinflammation markers also improved the accuracy of the diagnosis of AD pathology. None of the inflammatory markers correlated with the CSF Aß1-42 levels. Six CSF markers (IL-15, MCP-1, VEGFR-1, sICAM1, sVCAM-1, and VEGF-D) correlated with the CSF tau and p-tau181 levels, and these associations remained significant after controlling for age, sex, cognitive impairment, and APOEε4 status. CONCLUSIONS: The identified serum and CSF neuroinflammation biomarker signatures improve the accuracy of classification for AD pathology in older adults. Our results suggest that inflammation, vascular injury, and angiogenesis as reflected by CSF markers are closely related to cerebral tau pathology.

Genome-wide analysis of SREBP1 activity around the clock reveals its combined dependency on nutrient and circadian signals.

In mammals, the circadian clock allows them to anticipate and adapt physiology around the 24 hours. Conversely, metabolism and food consumption regulate the internal clock, pointing the existence of an intricate relationship between nutrient state and circadian homeostasis that is far from being understood. The Sterol Regulatory Element Binding Protein 1 (SREBP1) is a key regulator of lipid homeostasis. Hepatic SREBP1 function is influenced by the nutrient-response cycle, but also by the circadian machinery. To systematically understand how the interplay of circadian clock and nutrient-driven rhythm regulates SREBP1 activity, we evaluated the genome-wide binding of SREBP1 to its targets throughout the day in C57BL/6 mice. The recruitment of SREBP1 to the DNA showed a highly circadian behaviour, with a maximum during the fed status. However, the temporal expression of SREBP1 targets was not always synchronized with its binding pattern. In particular, different expression phases were observed for SREBP1 target genes depending on their function, suggesting the involvement of other transcription factors in their regulation. Binding sites for Hepatocyte Nuclear Factor 4 (HNF4) were specifically enriched in the close proximity of SREBP1 peaks of genes, whose expression was shifted by about 8 hours with respect to SREBP1 binding. Thus, the cross-talk between hepatic HNF4 and SREBP1 may underlie the expression timing of this subgroup of SREBP1 targets. Interestingly, the proper temporal expression profile of these genes was dramatically changed in Bmal1-/- mice upon time-restricted feeding, for which a rhythmic, but slightly delayed, binding of SREBP1 was maintained. Collectively, our results show that besides the nutrient-driven regulation of SREBP1 nuclear translocation, a second layer of modulation of SREBP1 transcriptional activity, strongly dependent from the circadian clock, exists. This system allows us to fine tune the expression timing of SREBP1 target genes, thus helping to temporally separate the different physiological processes in which these genes are involved.

The complex SNP and CNV genetic architecture of the increased risk of congenital heart defects in Down syndrome.

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21-specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. RESULTS: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays. AVAILABILITY: The R package abs filter for library clonality simulations and detection of amplification-biased sites is available from http://updepla1srv1.epfl.ch/waszaks/absfilter

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.