Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution.

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development.

Transient naive reprogramming corrects hiPS cells functionally and epigenetically.

Cells undergo a major epigenome reconfiguration when reprogrammed to human induced pluripotent stem cells (hiPS cells). However, the epigenomes of hiPS cells and human embryonic stem (hES) cells differ significantly, which affects hiPS cell function1-8. These differences include epigenetic memory and aberrations that emerge during reprogramming, for which the mechanisms remain unknown. Here we characterized the persistence and emergence of these epigenetic differences by performing genome-wide DNA methylation profiling throughout primed and naive reprogramming of human somatic cells to hiPS cells. We found that reprogramming-induced epigenetic aberrations emerge midway through primed reprogramming, whereas DNA demethylation begins early in naive reprogramming. Using this knowledge, we developed a transient-naive-treatment (TNT) reprogramming strategy that emulates the embryonic epigenetic reset. We show that the epigenetic memory in hiPS cells is concentrated in cell of origin-dependent repressive chromatin marked by H3K9me3, lamin-B1 and aberrant CpH methylation. TNT reprogramming reconfigures these domains to a hES cell-like state and does not disrupt genomic imprinting. Using an isogenic system, we demonstrate that TNT reprogramming can correct the transposable element overexpression and differential gene expression seen in conventional hiPS cells, and that TNT-reprogrammed hiPS and hES cells show similar differentiation efficiencies. Moreover, TNT reprogramming enhances the differentiation of hiPS cells derived from multiple cell types. Thus, TNT reprogramming corrects epigenetic memory and aberrations, producing hiPS cells that are molecularly and functionally more similar to hES cells than conventional hiPS cells. We foresee TNT reprogramming becoming a new standard for biomedical and therapeutic applications and providing a novel system for studying epigenetic memory.

Reprogramming roadmap reveals route to human induced trophoblast stem cells.

The reprogramming of human somatic cells to primed or naive induced pluripotent stem cells recapitulates the stages of early embryonic development1-6. The molecular mechanism that underpins these reprogramming processes remains largely unexplored, which impedes our understanding and limits rational improvements to reprogramming protocols. Here, to address these issues, we reconstruct molecular reprogramming trajectories of human dermal fibroblasts using single-cell transcriptomics. This revealed that reprogramming into primed and naive pluripotency follows diverging and distinct trajectories. Moreover, genome-wide analyses of accessible chromatin showed key changes in the regulatory elements of core pluripotency genes, and orchestrated global changes in chromatin accessibility over time. Integrated analysis of these datasets revealed a role for transcription factors associated with the trophectoderm lineage, and the existence of a subpopulation of cells that enter a trophectoderm-like state during reprogramming. Furthermore, this trophectoderm-like state could be captured, which enabled the derivation of induced trophoblast stem cells. Induced trophoblast stem cells are molecularly and functionally similar to trophoblast stem cells derived from human blastocysts or first-trimester placentas7. Our results provide a high-resolution roadmap for the transcription-factor-mediated reprogramming of human somatic cells, indicate a role for the trophectoderm-lineage-specific regulatory program during this process, and facilitate the direct reprogramming of somatic cells into induced trophoblast stem cells.

Transcriptome profiling reveals expression signatures of cranial neural crest cells arising from different axial levels.

BACKGROUND: Cranial neural crest cells (NCCs) are a unique embryonic cell type which give rise to a diverse array of derivatives extending from neurons and glia through to bone and cartilage. Depending on their point of origin along the antero-posterior axis cranial NCCs are rapidly sorted into distinct migratory streams that give rise to axial specific structures. These migratory streams mirror the underlying segmentation of the brain with NCCs exiting the diencephalon and midbrain following distinct paths compared to those exiting the hindbrain rhombomeres (r). The genetic landscape of cranial NCCs arising at different axial levels remains unknown. RESULTS: Here we have used RNA sequencing to uncover the transcriptional profiles of mouse cranial NCCs arising at different axial levels. Whole transcriptome analysis identified over 120 transcripts differentially expressed between NCCs arising anterior to r3 (referred to as r1-r2 migratory stream for simplicity) and the r4 migratory stream. Eight of the genes differentially expressed between these populations were validated by RT-PCR with 2 being further validated by in situ hybridisation. We also explored the expression of the Neuropilins (Nrp1 and Nrp2) and their co-receptors and show that the A-type Plexins are differentially expressed in different cranial NCC streams. CONCLUSIONS: Our analyses identify a large number of genes differentially regulated between cranial NCCs arising at different axial levels. This data provides a comprehensive description of the genetic landscape driving diversity of distinct cranial NCC streams and provides novel insight into the regulatory networks controlling the formation of specific skeletal elements and the mechanisms promoting migration along different paths.

Placental transcriptome co-expression analysis reveals conserved regulatory programs across gestation.

BACKGROUND: Mammalian development in utero is absolutely dependent on proper placental development, which is ultimately regulated by the placental genome. The regulation of the placental genome can be directly studied by exploring the underlying organisation of the placental transcriptome through a systematic analysis of gene-wise co-expression relationships. RESULTS: In this study, we performed a comprehensive analysis of human placental co-expression using RNA sequencing and intergrated multiple transcriptome datasets spanning human gestation. We identified modules of co-expressed genes that are preserved across human gestation, and also identifed modules conserved in the mouse indicating conserved molecular networks involved in placental development and gene expression patterns more specific to late gestation. Analysis of co-expressed gene flanking sequences indicated that conserved co-expression modules in the placenta are regulated by a core set of transcription factors, including ZNF423 and EBF1. Additionally, we identified a gene co-expression module enriched for genes implicated in the pregnancy pathology preeclampsia. By using an independnet transcriptome dataset, we show that these co-expressed genes are differentially expressed in preeclampsia. CONCLUSIONS: This study represents a comprehensive characterisation of placental co-expression and provides insight into potential transcriptional regulators that govern conserved molecular programs fundamental to placental development.

Recent progress towards understanding the role of DNA methylation in human placental development.

Epigenetic modifications, and particularly DNA methylation, have been studied in many tissues, both healthy and diseased, and across numerous developmental stages. The placenta is the only organ that has a transient life of 9 months and undergoes rapid growth and dynamic structural and functional changes across gestation. Additionally, the placenta is unique because although developing within the mother, its genome is identical to that of the foetus. Given these distinctive characteristics, it is not surprising that the epigenetic landscape affecting placental gene expression may be different to that in other healthy tissues. However, the role of epigenetic modifications, and particularly DNA methylation, in placental development remains largely unknown. Of particular interest is the fact that the placenta is the most hypomethylated human tissue and is characterized by the presence of large partially methylated domains (PMDs) containing silenced genes. Moreover, how and why the placenta is hypomethylated and what role DNA methylation plays in regulating placental gene expression across gestation are poorly understood. We review genome-wide DNA methylation studies in the human placenta and highlight that the different cell types that make up the placenta have very different DNA methylation profiles. Summarizing studies on DNA methylation in the placenta and its relationship with pregnancy complications are difficult due to the limited number of studies available for comparison. To understand the key steps in placental development and hence what may be perturbed in pregnancy complications requires large-scale genome-wide DNA methylation studies coupled with transcriptome analyses.

massiR: a method for predicting the sex of samples in gene expression microarray datasets.

UNLABELLED: High-throughput gene expression microarrays are currently the most efficient method for transcriptome-wide expression analyses. Consequently, gene expression data available through public repositories have largely been obtained from microarray experiments. However, the metadata associated with many publicly available expression microarray datasets often lacks sample sex information, therefore limiting the reuse of these data in new analyses or larger meta-analyses where the effect of sex is to be considered. Here, we present the massiR package, which provides a method for researchers to predict the sex of samples in microarray datasets. Using information from microarray probes representing Y chromosome genes, this package implements unsupervised clustering methods to classify samples into male and female groups, providing an efficient way to identify or confirm the sex of samples in mammalian microarray datasets. AVAILABILITY AND IMPLEMENTATION: massiR is implemented as a Bioconductor package in R. The package and the vignette can be downloaded at bioconductor.org and are provided under a GPL-2 license.

Integrative transcriptome meta-analysis reveals widespread sex-biased gene expression at the human fetal-maternal interface.

As males and females share highly similar genomes, the regulation of many sexually dimorphic traits is constrained to occur through sex-biased gene regulation. There is strong evidence that human males and females differ in terms of growth and development in utero and that these divergent growth strategies appear to place males at increased risk when in sub-optimal conditions. Since the placenta is the interface of maternal-fetal exchange throughout pregnancy, these developmental differences are most likely orchestrated by differential placental function. To date, progress in this field has been hampered by a lack of genome-wide information on sex differences in placental gene expression. Therefore, our motivation in this study was to characterize sex-biased gene expression in the human placenta. We obtained gene expression data for >300 non-pathological placenta samples from 11 microarray datasets and applied mapping-based array probe re-annotation and inverse-variance meta-analysis methods which showed that >140 genes (false discovery rate (FDR) <0.05) are differentially expressed between male and female placentae. A majority of these genes (>60%) are autosomal, many of which are involved in high-level regulatory processes such as gene transcription, cell growth and proliferation and hormonal function. Of particular interest, we detected higher female expression from all seven genes in the LHB-CGB cluster, which includes genes involved in placental development, the maintenance of pregnancy and maternal immune tolerance of the conceptus. These results demonstrate that sex-biased gene expression in the normal human placenta occurs across the genome and includes genes that are central to growth, development and the maintenance of pregnancy.

Approaches for the Analysis and Interpretation of Whole-Genome Bisulfite Sequencing Data.

DNA methylation is a covalent modification of DNA that plays important roles in processes such as the regulation of gene expression, transcription factor binding, and suppression of transposable elements. The use of whole-genome bisulfite sequencing (WGBS) enables the genome-wide identification and quantification of DNA methylation patterns at single-base resolution and is the gold standard for the analysis of DNA methylation. However, the computational analysis of WGBS data can be particularly challenging, as many computationally intensive steps are required. Here, we outline step-by-step an approach for the analysis and interpretation of WGBS data. First, sequencing reads must be trimmed, quality-checked, and aligned to the genome. Second, DNA methylation levels are estimated at each cytosine position using the aligned sequence reads of the bisulfite-treated DNA. Third, regions of differential cytosine methylation between samples can be identified. Finally, these data need to be visualized and interpreted in the context of the biological question at hand.

Future-proofing genomic data and consent management: a comprehensive review of technology innovations.

Genomic information is increasingly used to inform medical treatments and manage future disease risks. However, any personal and societal gains must be carefully balanced against the risk to individuals contributing their genomic data. Expanding our understanding of actionable genomic insights requires researchers to access large global datasets to capture the complexity of genomic contribution to diseases. Similarly, clinicians need efficient access to a patient's genome as well as population-representative historical records for evidence-based decisions. Both researchers and clinicians hence rely on participants to consent to the use of their genomic data, which in turn requires trust in the professional and ethical handling of this information. Here, we review existing and emerging solutions for secure and effective genomic information management, including storage, encryption, consent, and authorization that are needed to build participant trust. We discuss recent innovations in cloud computing, quantum-computing-proof encryption, and self-sovereign identity. These innovations can augment key developments from within the genomics community, notably GA4GH Passports and the Crypt4GH file container standard. We also explore how decentralized storage as well as the digital consenting process can offer culturally acceptable processes to encourage data contributions from ethnic minorities. We conclude that the individual and their right for self-determination needs to be put at the center of any genomics framework, because only on an individual level can the received benefits be accurately balanced against the risk of exposing private information.

Circulating epigenomic biomarkers correspond with kidney disease susceptibility in high-risk populations with type 2 diabetes mellitus.

AIMS: To investigate epigenomic indices of diabetic kidney disease (DKD) susceptibility among high-risk populations with type 2 diabetes mellitus. METHODS: KDIGO (Kidney Disease: Improving Global Outcomes) clinical guidelines were used to classify people living with or without DKD. Differential gene methylation of DKD was then assessed in a discovery Aboriginal Diabetes Study cohort (PROPHECY, 89 people) and an external independent study from Thailand (THEPTARIN, 128 people). Corresponding mRNA levels were also measured and linked to levels of albuminuria and eGFR. RESULTS: Increased DKD risk was associated with reduced methylation and elevated gene expression in the PROPHECY discovery cohort of Aboriginal Australians and these findings were externally validated in the THEPTARIN diabetes registry of Thai people living with type 2 diabetes mellitus. CONCLUSIONS: Novel epigenomic scores can improve diagnostic performance over clinical modelling using albuminuria and GFR alone and can distinguish DKD susceptibility.

Large-scale manipulation of promoter DNA methylation reveals context-specific transcriptional responses and stability.

BACKGROUND: Cytosine DNA methylation is widely described as a transcriptional repressive mark with the capacity to silence promoters. Epigenome engineering techniques enable direct testing of the effect of induced DNA methylation on endogenous promoters; however, the downstream effects have not yet been comprehensively assessed. RESULTS: Here, we simultaneously induce methylation at thousands of promoters in human cells using an engineered zinc finger-DNMT3A fusion protein, enabling us to test the effect of forced DNA methylation upon transcription, chromatin accessibility, histone modifications, and DNA methylation persistence after the removal of the fusion protein. We find that transcriptional responses to DNA methylation are highly context-specific, including lack of repression, as well as cases of increased gene expression, which appears to be driven by the eviction of methyl-sensitive transcriptional repressors. Furthermore, we find that some regulatory networks can override DNA methylation and that promoter methylation can cause alternative promoter usage. DNA methylation deposited at promoter and distal regulatory regions is rapidly erased after removal of the zinc finger-DNMT3A fusion protein, in a process combining passive and TET-mediated demethylation. Finally, we demonstrate that induced DNA methylation can exist simultaneously on promoter nucleosomes that possess the active histone modification H3K4me3, or DNA bound by the initiated form of RNA polymerase II. CONCLUSIONS: These findings have important implications for epigenome engineering and demonstrate that the response of promoters to DNA methylation is more complex than previously appreciated.

The emergence of the brain non-CpG methylation system in vertebrates.

Mammalian brains feature exceptionally high levels of non-CpG DNA methylation alongside the canonical form of CpG methylation. Non-CpG methylation plays a critical regulatory role in cognitive function, which is mediated by the binding of MeCP2, the transcriptional regulator that when mutated causes Rett syndrome. However, it is unclear whether the non-CpG neural methylation system is restricted to mammalian species with complex cognitive abilities or has deeper evolutionary origins. To test this, we investigated brain DNA methylation across 12 distantly related animal lineages, revealing that non-CpG methylation is restricted to vertebrates. We discovered that in vertebrates, non-CpG methylation is enriched within a highly conserved set of developmental genes transcriptionally repressed in adult brains, indicating that it demarcates a deeply conserved regulatory program. We also found that the writer of non-CpG methylation, DNMT3A, and the reader, MeCP2, originated at the onset of vertebrates as a result of the ancestral vertebrate whole-genome duplication. Together, we demonstrate how this novel layer of epigenetic information assembled at the root of vertebrates and gained new regulatory roles independent of the ancestral form of the canonical CpG methylation. This suggests that the emergence of non-CpG methylation may have fostered the evolution of sophisticated cognitive abilities found in the vertebrate lineage.

Transcriptional signature in microglia associated with Aß plaque phagocytosis.

The role of microglia cells in Alzheimer's disease (AD) is well recognized, however their molecular and functional diversity remain unclear. Here, we isolated amyloid plaque-containing (using labelling with methoxy-XO4, XO4+) and non-containing (XO4-) microglia from an AD mouse model. Transcriptomics analysis identified different transcriptional trajectories in ageing and AD mice. XO4+ microglial transcriptomes demonstrated dysregulated expression of genes associated with late onset AD. We further showed that the transcriptional program associated with XO4+ microglia from mice is present in a subset of human microglia isolated from brains of individuals with AD. XO4- microglia displayed transcriptional signatures associated with accelerated ageing and contained more intracellular post-synaptic material than XO4+ microglia, despite reduced active synaptosome phagocytosis. We identified HIF1α as potentially regulating synaptosome phagocytosis in vitro using primary human microglia, and BV2 mouse microglial cells. Together, these findings provide insight into molecular mechanisms underpinning the functional diversity of microglia in AD.

Convergent evolution of a vertebrate-like methylome in a marine sponge.

Vertebrates have highly methylated genomes at CpG positions, whereas invertebrates have sparsely methylated genomes. This increase in methylation content is considered a major regulatory innovation of vertebrate genomes. However, here we report that a sponge, proposed as the potential sister group to the rest of animals, has a highly methylated genome. Despite major differences in genome size and architecture, we find similarities between the independent acquisitions of the hypermethylated state. Both lineages show genome-wide CpG depletion, conserved strong transcription factor methyl-sensitivity and developmental methylation dynamics at 5-hydroxymethylcytosine enriched regions. Together, our findings trace back patterns associated with DNA methylation in vertebrates to the early steps of animal evolution. Thus, the sponge methylome challenges previous hypotheses concerning the uniqueness of vertebrate genome hypermethylation and its implications for regulatory complexity.

A single-cell atlas of entorhinal cortex from individuals with Alzheimer's disease reveals cell-type-specific gene expression regulation.

There is currently little information available about how individual cell types contribute to Alzheimer's disease. Here we applied single-nucleus RNA sequencing to entorhinal cortex samples from control and Alzheimer's disease brains (n = 6 per group), yielding a total of 13,214 high-quality nuclei. We detail cell-type-specific gene expression patterns, unveiling how transcriptional changes in specific cell subpopulations are associated with Alzheimer's disease. We report that the Alzheimer's disease risk gene APOE is specifically repressed in Alzheimer's disease oligodendrocyte progenitor cells and astrocyte subpopulations and upregulated in an Alzheimer's disease-specific microglial subopulation. Integrating transcription factor regulatory modules with Alzheimer's disease risk loci revealed drivers of cell-type-specific state transitions towards Alzheimer's disease. For example, transcription factor EB, a master regulator of lysosomal function, regulates multiple disease genes in a specific Alzheimer's disease astrocyte subpopulation. These results provide insights into the coordinated control of Alzheimer's disease risk genes and their cell-type-specific contribution to disease susceptibility. These results are available at http://adsn.ddnetbio.com.

Approaches for the Analysis and Interpretation of Whole Genome Bisulfite Sequencing Data.

DNA methylation is a covalent modification of DNA that plays important roles in processes such as the regulation of gene expression, transcription factor binding, and suppression of transposable elements. The use of whole genome bisulfite sequencing (WGBS) enables the genome-wide identification and quantification of DNA methylation patterns at single-base resolution and is the gold standard for analysis of DNA methylation. Computational analysis of WGBS data can be particularly challenging, as many computationally intensive steps are required. Here, we outline a step-by-step approach for the analysis and interpretation of WGBS data. First, sequencing reads must be trimmed, quality checked, and aligned to the genome. Second, DNA methylation levels are estimated at each cytosine position using the aligned sequence reads of the bisulfite treated DNA. Third, regions of differential cytosine methylation between samples can be identified. Finally, these data need to be visualized and interpreted in the context of the biological question at hand.

msgbsR: An R package for analysing methylation-sensitive restriction enzyme sequencing data.

Genotyping-by-sequencing (GBS) or restriction-site associated DNA marker sequencing (RAD-seq) is a practical and cost-effective method for analysing large genomes from high diversity species. This method of sequencing, coupled with methylation-sensitive enzymes (often referred to as methylation-sensitive restriction enzyme sequencing or MRE-seq), is an effective tool to study DNA methylation in parts of the genome that are inaccessible in other sequencing techniques or are not annotated in microarray technologies. Current software tools do not fulfil all methylation-sensitive restriction sequencing assays for determining differences in DNA methylation between samples. To fill this computational need, we present msgbsR, an R package that contains tools for the analysis of methylation-sensitive restriction enzyme sequencing experiments. msgbsR can be used to identify and quantify read counts at methylated sites directly from alignment files (BAM files) and enables verification of restriction enzyme cut sites with the correct recognition sequence of the individual enzyme. In addition, msgbsR assesses DNA methylation based on read coverage, similar to RNA sequencing experiments, rather than methylation proportion and is a useful tool in analysing differential methylation on large populations. The package is fully documented and available freely online as a Bioconductor package ( https://bioconductor.org/packages/release/bioc/html/msgbsR.html ).

Recurrent acquisition of cytosine methyltransferases into eukaryotic retrotransposons.

Transposable elements are in a constant arms race with the silencing mechanisms of their host genomes. One silencing mechanism commonly used by many eukaryotes is dependent on cytosine methylation, a covalent modification of DNA deposited by C5 cytosine methyltransferases (DNMTs). Here, we report how two distantly related eukaryotic lineages, dinoflagellates and charophytes, have independently incorporated DNMTs into the coding regions of distinct retrotransposon classes. Concomitantly, we show that dinoflagellates of the genus Symbiodinium have evolved cytosine methylation patterns unlike any other eukaryote, with most of the genome methylated at CG dinucleotides. Finally, we demonstrate the ability of retrotransposon DNMTs to methylate CGs de novo, suggesting that retrotransposons could self-methylate retrotranscribed DNA. Together, this is an example of how retrotransposons incorporate host-derived genes involved in DNA methylation. In some cases, this event could have implications for the composition and regulation of the host epigenomic environment.

Author Correction: Recurrent acquisition of cytosine methyltransferases into eukaryotic retrotransposons.

The original version of this Article contained an error in the spelling of the author Hongfei Li, which was incorrectly given as Fei Hong. This has now been corrected in both the PDF and HTML versions of the Article.