DNA sequence-dependent compartmentalization and silencing of chromatin at the nuclear lamina.

A large fraction of the mammalian genome is organized into inactive chromosomal domains along the nuclear lamina. The mechanism by which these lamina associated domains (LADs) are established remains to be elucidated. Using genomic repositioning assays, we show that LADs, spanning the developmentally regulated IgH and Cyp3a loci contain discrete DNA regions that associate chromatin with the nuclear lamina and repress gene activity in fibroblasts. Lamina interaction is established during mitosis and likely involves the localized recruitment of Lamin B during late anaphase. Fine-scale mapping of LADs reveals numerous lamina-associating sequences (LASs), which are enriched for a GAGA motif. This repeated motif directs lamina association and is bound by the transcriptional repressor cKrox, in a complex with HDAC3 and Lap2ß. Knockdown of cKrox or HDAC3 results in dissociation of LASs/LADs from the nuclear lamina. These results reveal a mechanism that couples nuclear compartmentalization of chromatin domains with the control of gene activity.

Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes.

Methods to deconvolve single-cell RNA-sequencing (scRNA-seq) data are necessary for samples containing a mixture of genotypes, whether they are natural or experimentally combined. Multiplexing across donors is a popular experimental design that can avoid batch effects, reduce costs and improve doublet detection. By using variants detected in scRNA-seq reads, it is possible to assign cells to their donor of origin and identify cross-genotype doublets that may have highly similar transcriptional profiles, precluding detection by transcriptional profile. More subtle cross-genotype variant contamination can be used to estimate the amount of ambient RNA. Ambient RNA is caused by cell lysis before droplet partitioning and is an important confounder of scRNA-seq analysis. Here we develop souporcell, a method to cluster cells using the genetic variants detected within the scRNA-seq reads. We show that it achieves high accuracy on genotype clustering, doublet detection and ambient RNA estimation, as demonstrated across a range of challenging scenarios.

Cardelino: computational integration of somatic clonal substructure and single-cell transcriptomes.

Bulk and single-cell DNA sequencing has enabled reconstructing clonal substructures of somatic tissues from frequency and cooccurrence patterns of somatic variants. However, approaches to characterize phenotypic variations between clones are not established. Here we present cardelino (https://github.com/single-cell-genetics/cardelino), a computational method for inferring the clonal tree configuration and the clone of origin of individual cells assayed using single-cell RNA-seq (scRNA-seq). Cardelino flexibly integrates information from imperfect clonal trees inferred based on bulk exome-seq data, and sparse variant alleles expressed in scRNA-seq data. We apply cardelino to a published cancer dataset and to newly generated matched scRNA-seq and exome-seq data from 32 human dermal fibroblast lines, identifying hundreds of differentially expressed genes between cells from different somatic clones. These genes are frequently enriched for cell cycle and proliferation pathways, indicating a role for cell division genes in somatic evolution in healthy skin.

Common genetic variation drives molecular heterogeneity in human iPSCs.

Technology utilizing human induced pluripotent stem cells (iPS cells) has enormous potential to provide improved cellular models of human disease. However, variable genetic and phenotypic characterization of many existing iPS cell lines limits their potential use for research and therapy. Here we describe the systematic generation, genotyping and phenotyping of 711 iPS cell lines derived from 301 healthy individuals by the Human Induced Pluripotent Stem Cells Initiative. Our study outlines the major sources of genetic and phenotypic variation in iPS cells and establishes their suitability as models of complex human traits and cancer. Through genome-wide profiling we find that 5-46% of the variation in different iPS cell phenotypes, including differentiation capacity and cellular morphology, arises from differences between individuals. Additionally, we assess the phenotypic consequences of genomic copy-number alterations that are repeatedly observed in iPS cells. In addition, we present a comprehensive map of common regulatory variants affecting the transcriptome of human pluripotent cells.

Corrigendum: Common genetic variation drives molecular heterogeneity in human iPSCs.

This corrects the article DOI: 10.1038/nature22403.

Activin/nodal signaling and NANOG orchestrate human embryonic stem cell fate decisions by controlling the H3K4me3 chromatin mark.

Stem cells can self-renew and differentiate into multiple cell types. These characteristics are maintained by the combination of specific signaling pathways and transcription factors that cooperate to establish a unique epigenetic state. Despite the broad interest of these mechanisms, the precise molecular controls by which extracellular signals organize epigenetic marks to confer multipotency remain to be uncovered. Here, we use human embryonic stem cells (hESCs) to show that the Activin-SMAD2/3 signaling pathway cooperates with the core pluripotency factor NANOG to recruit the DPY30-COMPASS histone modifiers onto key developmental genes. Functional studies demonstrate the importance of these interactions for correct histone 3 Lys4 trimethylation and also self-renewal and differentiation. Finally, genetic studies in mice show that Dpy30 is also necessary to maintain pluripotency in the pregastrulation embryo, thereby confirming the existence of similar regulations in vivo during early embryonic development. Our results reveal the mechanisms by which extracellular factors coordinate chromatin status and cell fate decisions in hESCs.

Correction to: Niacin-mediated rejuvenation of macrophage/microglia enhances remyelination of the aging central nervous system.

The article Niacin­mediated rejuvenation of macrophage/microglia enhances remyelination of the aging central nervous system, written by Khalil S. Rawji, Adam M.H. Young, Tanay Ghosh, Nathan J. Michaels, Reza Mirzaei, Janson Kappen, Kathleen L. Kolehmainen, Nima Alaeiilkhchi, Brian Lozinski, Manoj K. Mishra, Annie Pu, Weiwen Tang, Salma Zein, Deepak K. Kaushik, Michael B. Keough, Jason R. Plemel, Fiona Calvert, Andrew J. Knights, Daniel J. Gaffney, Wolfram Tetzlaff, Robin J. M. Franklin and V. Wee Yong, was originally published electronically on the publisher's internet.

Niacin-mediated rejuvenation of macrophage/microglia enhances remyelination of the aging central nervous system.

Remyelination following CNS demyelination restores rapid signal propagation and protects axons; however, its efficiency declines with increasing age. Both intrinsic changes in the oligodendrocyte progenitor cell population and extrinsic factors in the lesion microenvironment of older subjects contribute to this decline. Microglia and monocyte-derived macrophages are critical for successful remyelination, releasing growth factors and clearing inhibitory myelin debris. Several studies have implicated delayed recruitment of macrophages/microglia into lesions as a key contributor to the decline in remyelination observed in older subjects. Here we show that the decreased expression of the scavenger receptor CD36 of aging mouse microglia and human microglia in culture underlies their reduced phagocytic activity. Overexpression of CD36 in cultured microglia rescues the deficit in phagocytosis of myelin debris. By screening for clinically approved agents that stimulate macrophages/microglia, we have found that niacin (vitamin B3) upregulates CD36 expression and enhances myelin phagocytosis by microglia in culture. This increase in myelin phagocytosis is mediated through the niacin receptor (hydroxycarboxylic acid receptor 2). Genetic fate mapping and multiphoton live imaging show that systemic treatment of 9-12-month-old demyelinated mice with therapeutically relevant doses of niacin promotes myelin debris clearance in lesions by both peripherally derived macrophages and microglia. This is accompanied by enhancement of oligodendrocyte progenitor cell numbers and by improved remyelination in the treated mice. Niacin represents a safe and translationally amenable regenerative therapy for chronic demyelinating diseases such as multiple sclerosis.

CRISPR/Cas9 editing in human pluripotent stem cell-cardiomyocytes highlights arrhythmias, hypocontractility, and energy depletion as potential therapeutic targets for hypertrophic cardiomyopathy.

Aims: Sarcomeric gene mutations frequently underlie hypertrophic cardiomyopathy (HCM), a prevalent and complex condition leading to left ventricle thickening and heart dysfunction. We evaluated isogenic genome-edited human pluripotent stem cell-cardiomyocytes (hPSC-CM) for their validity to model, and add clarity to, HCM. Methods and results: CRISPR/Cas9 editing produced 11 variants of the HCM-causing mutation c.C9123T-MYH7 [(p.R453C-ß-myosin heavy chain (MHC)] in 3 independent hPSC lines. Isogenic sets were differentiated to hPSC-CMs for high-throughput, non-subjective molecular and functional assessment using 12 approaches in 2D monolayers and/or 3D engineered heart tissues. Although immature, edited hPSC-CMs exhibited the main hallmarks of HCM (hypertrophy, multi-nucleation, hypertrophic marker expression, sarcomeric disarray). Functional evaluation supported the energy depletion model due to higher metabolic respiration activity, accompanied by abnormalities in calcium handling, arrhythmias, and contraction force. Partial phenotypic rescue was achieved with ranolazine but not omecamtiv mecarbil, while RNAseq highlighted potentially novel molecular targets. Conclusion: Our holistic and comprehensive approach showed that energy depletion affected core cardiomyocyte functionality. The engineered R453C-ßMHC-mutation triggered compensatory responses in hPSC-CMs, causing increased ATP production and αMHC to energy-efficient ßMHC switching. We showed that pharmacological rescue of arrhythmias was possible, while MHY7: MYH6 and mutant: wild-type MYH7 ratios may be diagnostic, and previously undescribed lncRNAs and gene modifiers are suggestive of new mechanisms.

DNase I sensitivity QTLs are a major determinant of human expression variation.

The mapping of expression quantitative trait loci (eQTLs) has emerged as an important tool for linking genetic variation to changes in gene regulation. However, it remains difficult to identify the causal variants underlying eQTLs, and little is known about the regulatory mechanisms by which they act. Here we show that genetic variants that modify chromatin accessibility and transcription factor binding are a major mechanism through which genetic variation leads to gene expression differences among humans. We used DNase I sequencing to measure chromatin accessibility in 70 Yoruba lymphoblastoid cell lines, for which genome-wide genotypes and estimates of gene expression levels are also available. We obtained a total of 2.7 billion uniquely mapped DNase I-sequencing (DNase-seq) reads, which allowed us to produce genome-wide maps of chromatin accessibility for each individual. We identified 8,902 locations at which the DNase-seq read depth correlated significantly with genotype at a nearby single nucleotide polymorphism or insertion/deletion (false discovery rate = 10%). We call such variants 'DNase I sensitivity quantitative trait loci' (dsQTLs). We found that dsQTLs are strongly enriched within inferred transcription factor binding sites and are frequently associated with allele-specific changes in transcription factor binding. A substantial fraction (16%) of dsQTLs are also associated with variation in the expression levels of nearby genes (that is, these loci are also classified as eQTLs). Conversely, we estimate that as many as 55% of eQTL single nucleotide polymorphisms are also dsQTLs. Our observations indicate that dsQTLs are highly abundant in the human genome and are likely to be important contributors to phenotypic variation.

Genetic background drives transcriptional variation in human induced pluripotent stem cells.

Human iPS cells have been generated using a diverse range of tissues from a variety of donors using different reprogramming vectors. However, these cell lines are heterogeneous, which presents a limitation for their use in disease modeling and personalized medicine. To explore the basis of this heterogeneity we generated 25 iPS cell lines under normalised conditions from the same set of somatic tissues across a number of donors. RNA-seq data sets from each cell line were compared to identify the majority contributors to transcriptional heterogeneity. We found that genetic differences between individual donors were the major cause of transcriptional variation between lines. In contrast, residual signatures from the somatic cell of origin, so called epigenetic memory, contributed relatively little to transcriptional variation. Thus, underlying genetic background variation is responsible for most heterogeneity between human iPS cell lines. We conclude that epigenetic effects in hIPSCs are minimal, and that hIPSCs are a stable, robust and powerful platform for large-scale studies of the function of genetic differences between individuals. Our data also suggest that future studies using hIPSCs as a model system should focus most effort on collection of large numbers of donors, rather than generating large numbers of lines from the same donor.

Early maturation and distinct tau pathology in induced pluripotent stem cell-derived neurons from patients with MAPT mutations.

Tauopathies, such as Alzheimer's disease, some cases of frontotemporal dementia, corticobasal degeneration and progressive supranuclear palsy, are characterized by aggregates of the microtubule-associated protein tau, which are linked to neuronal death and disease development and can be caused by mutations in the MAPT gene. Six tau isoforms are present in the adult human brain and they differ by the presence of 3(3R) or 4(4R) C-terminal repeats. Only the shortest 3R isoform is present in foetal brain. MAPT mutations found in human disease affect tau binding to microtubules or the 3R:4R isoform ratio by altering exon 10 splicing. We have differentiated neurons from induced pluripotent stem cells derived from fibroblasts of controls and patients with N279K and P301L MAPT mutations. Induced pluripotent stem cell-derived neurons recapitulate developmental tau expression, showing the adult brain tau isoforms after several months in culture. Both N279K and P301L neurons exhibit earlier electrophysiological maturation and altered mitochondrial transport compared to controls. Specifically, the N279K neurons show abnormally premature developmental 4R tau expression, including changes in the 3R:4R isoform ratio and AT100-hyperphosphorylated tau aggregates, while P301L neurons are characterized by contorted processes with varicosity-like structures, some containing both alpha-synuclein and 4R tau. The previously unreported faster maturation of MAPT mutant human neurons, the developmental expression of 4R tau and the morphological alterations may contribute to disease development.

Global properties and functional complexity of human gene regulatory variation.

Identification and functional interpretation of gene regulatory variants is a major focus of modern genomics. The application of genetic mapping to molecular and cellular traits has enabled the detection of regulatory variation on genome-wide scales and revealed an enormous diversity of regulatory architecture in humans and other species. In this review I summarise the insights gained and questions raised by a decade of genetic mapping of gene expression variation. I discuss recent extensions of this approach using alternative molecular phenotypes that have revealed some of the biological mechanisms that drive gene expression variation between individuals. Finally, I highlight outstanding problems and future directions for development.

Coloured sweat in two brothers: First report of familial chromhidrosis.

The uncommon diagnosis of chromhidrosis is most frequently made in young adults. This sweat gland disease, although benign, may impact significantly on the patient's quality of life. We describe the first report of familial chromhidrosis of pseudo-eccrine type (pseudochromhidrosis) occurring in two brothers aged 9 and 12 years. The classification and causality of chromhidrosis is described and approaches to assessment and management are outlined.

Controls of nucleosome positioning in the human genome.

Nucleosomes are important for gene regulation because their arrangement on the genome can control which proteins bind to DNA. Currently, few human nucleosomes are thought to be consistently positioned across cells; however, this has been difficult to assess due to the limited resolution of existing data. We performed paired-end sequencing of micrococcal nuclease-digested chromatin (MNase-seq) from seven lymphoblastoid cell lines and mapped over 3.6 billion MNase-seq fragments to the human genome to create the highest-resolution map of nucleosome occupancy to date in a human cell type. In contrast to previous results, we find that most nucleosomes have more consistent positioning than expected by chance and a substantial fraction (8.7%) of nucleosomes have moderate to strong positioning. In aggregate, nucleosome sequences have 10 bp periodic patterns in dinucleotide frequency and DNase I sensitivity; and, across cells, nucleosomes frequently have translational offsets that are multiples of 10 bp. We estimate that almost half of the genome contains regularly spaced arrays of nucleosomes, which are enriched in active chromatin domains. Single nucleotide polymorphisms that reduce DNase I sensitivity can disrupt the phasing of nucleosome arrays, which indicates that they often result from positioning against a barrier formed by other proteins. However, nucleosome arrays can also be created by DNA sequence alone. The most striking example is an array of over 400 nucleosomes on chromosome 12 that is created by tandem repetition of sequences with strong positioning properties. In summary, a large fraction of nucleosomes are consistently positioned--in some regions because they adopt favored sequence positions, and in other regions because they are forced into specific arrangements by chromatin remodeling or DNA binding proteins.

The contribution of RNA decay quantitative trait loci to inter-individual variation in steady-state gene expression levels.

Recent gene expression QTL (eQTL) mapping studies have provided considerable insight into the genetic basis for inter-individual regulatory variation. However, a limitation of all eQTL studies to date, which have used measurements of steady-state gene expression levels, is the inability to directly distinguish between variation in transcription and decay rates. To address this gap, we performed a genome-wide study of variation in gene-specific mRNA decay rates across individuals. Using a time-course study design, we estimated mRNA decay rates for over 16,000 genes in 70 Yoruban HapMap lymphoblastoid cell lines (LCLs), for which extensive genotyping data are available. Considering mRNA decay rates across genes, we found that: (i) as expected, highly expressed genes are generally associated with lower mRNA decay rates, (ii) genes with rapid mRNA decay rates are enriched with putative binding sites for miRNA and RNA binding proteins, and (iii) genes with similar functional roles tend to exhibit correlated rates of mRNA decay. Focusing on variation in mRNA decay across individuals, we estimate that steady-state expression levels are significantly correlated with variation in decay rates in 10% of genes. Somewhat counter-intuitively, for about half of these genes, higher expression is associated with faster decay rates, possibly due to a coupling of mRNA decay with transcriptional processes in genes involved in rapid cellular responses. Finally, we used these data to map genetic variation that is specifically associated with variation in mRNA decay rates across individuals. We found 195 such loci, which we named RNA decay quantitative trait loci ("rdQTLs"). All the observed rdQTLs are located near the regulated genes and therefore are assumed to act in cis. By analyzing our data within the context of known steady-state eQTLs, we estimate that a substantial fraction of eQTLs are associated with inter-individual variation in mRNA decay rates.

Randomized comparison of volar locking plates and intramedullary nails for unstable distal radius fractures.

PURPOSE: To compare minimally invasive intramedullary nails (IMN) and volar locking plates (VLP) for the treatment of unstable distal radius fractures by evaluating postoperative subjective, radiographic, and functional outcomes. The hypothesis was that IMN patients would have less pain and required less pain medication in the early postoperative period and returned to work earlier than VLP patients. METHODS: Sixty patients with closed, displaced, unstable, extra-articular, metaphyseal fractures of the distal radius were randomized to receive a VLP or an IMN for internal fixation. Functional outcomes (Quick Disabilities of the Arm, Shoulder, and Hand, Michigan Hand Questionnaire), radiographic measurements (ulnar variance, radial height, inclination, and volar tilt), and range of motion were assessed until final follow-up at 2 years after surgery. Narcotic pain medication use was documented for 5 weeks following surgery. RESULTS: There were 2 groups of 30 patients with IMN (mean age, 55 ± 14 y) or VLP (mean age, 55 ± 16 y) with similar demographics and comorbidities. Patients with IMN regained extension earlier but had similar range of motion to patients with VLP at final follow-up. There was similar improvement in Michigan Hand Questionnaire, Quick Disabilities of the Arm, Shoulder, and Hand, and strength between groups. Five weeks after surgery, fewer IMN patients required narcotic pain medication (13%) than VLP patients (33%). Radiographic outcomes were similar at final follow-up. There were 3 failures with IMN versus 1 failure with VLP. All 10 employed patients with IMN returned to previous work compared with 10 of the 12 employed patients with VLP. Time to return to work was similar for both groups. CONCLUSIONS: In a cohort of similar patients, IMN and VLP provided comparable improvement in functional and radiographic outcomes. Patients with IMN required less narcotic pain medication after surgery than VLP patients. TYPE OF STUDY/LEVEL OF EVIDENCE: Therapeutic II.

Accurate inference of transcription factor binding from DNA sequence and chromatin accessibility data.

Accurate functional annotation of regulatory elements is essential for understanding global gene regulation. Here, we report a genome-wide map of 827,000 transcription factor binding sites in human lymphoblastoid cell lines, which is comprised of sites corresponding to 239 position weight matrices of known transcription factor binding motifs, and 49 novel sequence motifs. To generate this map, we developed a probabilistic framework that integrates cell- or tissue-specific experimental data such as histone modifications and DNase I cleavage patterns with genomic information such as gene annotation and evolutionary conservation. Comparison to empirical ChIP-seq data suggests that our method is highly accurate yet has the advantage of targeting many factors in a single assay. We anticipate that this approach will be a valuable tool for genome-wide studies of gene regulation in a wide variety of cell types or tissues under diverse conditions.

Malignant change in seborrhoeic keratoses in a region with high solar ultraviolet levels.

Regions with high solar UV levels and high skin cancer rates may experience a greater incidence of malignancy in association with seborrhoeic keratoses (SebK) than in low UV regions. Previous reports have indicated that basal cell carcinoma is the most common neoplasm with reported rates of up to 4 per cent of excised SebK. The rates of such compound lesions occurring in our practice were reviewed, indicating that Bowen's disease was the most frequently observed neoplasm with a rate of 7 per cent. In total, 10 per cent of all excised specimens showed either frank malignancy or some degree of atypia.

The epidermal growth factor receptor in squamous cell carcinoma: An emerging drug target.

Squamous cell carcinoma (SCC) of the skin is an increasingly common diagnosis worldwide. Together with precursor lesions, SCC carry a significant morbidity, particularly in regions with high solar UV radiation levels. Advanced lesions are locally or sometimes widely metastatic and may be resistant to treatment. Drugs targeting the epidermal growth factor receptor (EGFR) are currently the only significant non-surgical option for advanced SCC beyond radiotherapy and conventional chemotherapy. The role of the EGFR in skin cancer is described and the outcomes of targeted anti-EGFR therapy published to date are summarised. The future of anti-EGFR targeted therapies in the treatment of skin cancer is discussed. Targeted molecular therapies are becoming increasingly widespread and an understanding of the evidence for their use as well as their side effect profile is important in order to offer patients informed and current advice.