Programming and inheritance of parental DNA methylomes in mammals.

The reprogramming of parental methylomes is essential for embryonic development. In mammals, paternal 5-methylcytosines (5mCs) have been proposed to be actively converted to oxidized bases. These paternal oxidized bases and maternal 5mCs are believed to be passively diluted by cell divisions. By generating single-base resolution, allele-specific DNA methylomes from mouse gametes, early embryos, and primordial germ cell (PGC), as well as single-base-resolution maps of oxidized cytosine bases for early embryos, we report the existence of 5hmC and 5fC in both maternal and paternal genomes and find that 5mC or its oxidized derivatives, at the majority of demethylated CpGs, are converted to unmodified cytosines independent of passive dilution from gametes to four-cell embryos. Therefore, we conclude that paternal methylome and at least a significant proportion of maternal methylome go through active demethylation during embryonic development. Additionally, all the known imprinting control regions (ICRs) were classified into germ-line or somatic ICRs.

Blastocyst-like structures generated from human pluripotent stem cells.

Limited access to embryos has hampered the study of human embryogenesis and disorders that occur during early pregnancy. Human pluripotent stem cells provide an alternative means to study human development in a dish1-7. Recent advances in partial embryo models derived from human pluripotent stem cells have enabled human development to be examined at early post-implantation stages8-14. However, models of the pre-implantation human blastocyst are lacking. Starting from naive human pluripotent stem cells, here we developed an effective three-dimensional culture strategy with successive lineage differentiation and self-organization to generate blastocyst-like structures in vitro. These structures-which we term 'human blastoids'-resemble human blastocysts in terms of their morphology, size, cell number, and composition and allocation of different cell lineages. Single-cell RNA-sequencing analyses also reveal the transcriptomic similarity of blastoids to blastocysts. Human blastoids are amenable to embryonic and extra-embryonic stem cell derivation and can further develop into peri-implantation embryo-like structures in vitro. Using chemical perturbations, we show that specific isozymes of protein kinase C have a critical function in the formation of the blastoid cavity. Human blastoids provide a readily accessible, scalable, versatile and perturbable alternative to blastocysts for studying early human development, understanding early pregnancy loss and gaining insights into early developmental defects.

Multiplexed Engineering and Analysis of Combinatorial Enhancer Activity in Single Cells.

The study of enhancers has been hampered by the scarcity of methods to systematically quantify their endogenous activity. We develop Mosaic-seq to systematically perturb enhancers and measure their endogenous activities at single-cell resolution. Mosaic-seq uses a CRISPR barcoding system to jointly measure a cell's transcriptome and its sgRNA modulators, thus quantifying the effects of dCas9-KRAB-mediated enhancer repression in single cells. Applying Mosaic-seq to 71 constituent enhancers from 15 super-enhancers, our analysis of 51,448 sgRNA-induced transcriptomes finds that only a small number of constituents are major effectors of target gene expression. Binding of p300 and RNAPII are key features of these constituents. We determine two key parameters of enhancer activity in single cells: their penetrance in a population and their contribution to expression in these cells. Through combinatorial interrogation, we find that simultaneous repression of multiple weak constituents can alter super-enhancer activity in a manner greatly exceeding repression of individual constituents.

FBA: feature barcoding analysis for single cell RNA-Seq.

MOTIVATION: Single cell RNA-Seq (scRNA-Seq) has broadened our understanding of cellular heterogeneity and provided valuable insights into cellular functions. Recent experimental strategies extend scRNA-Seq readouts to include additional features, including cell surface proteins and genomic perturbations. These 'feature barcoding' strategies rely on converting molecular and cellular features to unique sequence barcodes, which are then detected with the transcriptome. RESULTS: Here, we introduce FBA, a flexible and streamlined package to perform quality control, quantification, demultiplexing, multiplet detection, clustering and visualization of feature barcoding assays. AVAILABILITYAND IMPLEMENTATION: FBA is available on PyPi at https://pypi.org/project/fba and on GitHub at https://github.com/jlduan/fba. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Global chromatin landscapes identify candidate noncoding modifiers of cardiac rhythm.

Comprehensive cis-regulatory landscapes are essential for accurate enhancer prediction and disease variant mapping. Although cis-regulatory element (CRE) resources exist for most tissues and organs, many rare - yet functionally important - cell types remain overlooked. Despite representing only a small fraction of the heart's cellular biomass, the cardiac conduction system (CCS) unfailingly coordinates every life-sustaining heartbeat. To globally profile the mouse CCS cis-regulatory landscape, we genetically tagged CCS component-specific nuclei for comprehensive assay for transposase-accessible chromatin-sequencing (ATAC-Seq) analysis. Thus, we established a global CCS-enriched CRE database, referred to as CCS-ATAC, as a key resource for studying CCS-wide and component-specific regulatory functions. Using transcription factor (TF) motifs to construct CCS component-specific gene regulatory networks (GRNs), we identified and independently confirmed several specific TF sub-networks. Highlighting the functional importance of CCS-ATAC, we also validated numerous CCS-enriched enhancer elements and suggested gene targets based on CCS single-cell RNA-Seq data. Furthermore, we leveraged CCS-ATAC to improve annotation of existing human variants related to cardiac rhythm and nominated a potential enhancer-target pair that was dysregulated by a specific SNP. Collectively, our results established a CCS-regulatory compendium, identified novel CCS enhancer elements, and illuminated potential functional associations between human genomic variants and CCS component-specific CREs.

Large-scale production of human blastoids amenable to modeling blastocyst development and maternal-fetal cross talk.

Recent advances in human blastoids have opened new avenues for modeling early human development and implantation. One limitation of our first protocol for human blastoid generation was relatively low efficiency. We now report an optimized protocol for the efficient generation of large quantities of high-fidelity human blastoids from naive pluripotent stem cells. This enabled proteomics analysis that identified phosphosite-specific signatures potentially involved in the derivation and/or maintenance of the signaling states in human blastoids. Additionally, we uncovered endometrial stromal effects in promoting trophoblast cell survival, proliferation, and syncytialization during co-culture with blastoids and blastocysts. Side-by-side single-cell RNA sequencing revealed similarities and differences in transcriptome profiles between pre-implantation blastoids and blastocysts, as well as post-implantation cultures, and uncovered a population resembling early migratory trophoblasts during co-culture with endometrial stromal cells. Our optimized protocol will facilitate broader use of human blastoids as an accessible, perturbable, scalable, and tractable model for human blastocysts.

Optimizing de novo common wheat transcriptome assembly using short-read RNA-Seq data.

BACKGROUND: Rapid advances in next-generation sequencing methods have provided new opportunities for transcriptome sequencing (RNA-Seq). The unprecedented sequencing depth provided by RNA-Seq makes it a powerful and cost-efficient method for transcriptome study, and it has been widely used in model organisms and non-model organisms to identify and quantify RNA. For non-model organisms lacking well-defined genomes, de novo assembly is typically required for downstream RNA-Seq analyses, including SNP discovery and identification of genes differentially expressed by phenotypes. Although RNA-Seq has been successfully used to sequence many non-model organisms, the results of de novo assembly from short reads can still be improved by using recent bioinformatic developments. RESULTS: In this study, we used 212.6 million pair-end reads, which accounted for 16.2 Gb, to assemble the hexaploid wheat transcriptome. Two state-of-the-art assemblers, Trinity and Trans-ABySS, which use the single and multiple k-mer methods, respectively, were used, and the whole de novo assembly process was divided into the following four steps: pre-assembly, merging different samples, removal of redundancy and scaffolding. We documented every detail of these steps and how these steps influenced assembly performance to gain insight into transcriptome assembly from short reads. After optimization, the assembled transcripts were comparable to Sanger-derived ESTs in terms of both continuity and accuracy. We also provided considerable new wheat transcript data to the community. CONCLUSIONS: It is feasible to assemble the hexaploid wheat transcriptome from short reads. Special attention should be paid to dealing with multiple samples to balance the spectrum of expression levels and redundancy. To obtain an accurate overview of RNA profiling, removal of redundancy may be crucial in de novo assembly.

De novo assembly and characterisation of the transcriptome during seed development, and generation of genic-SSR markers in peanut (Arachis hypogaea L.).

BACKGROUND: The peanut (Arachis hypogaea L.) is an important oilseed crop in tropical and subtropical regions of the world. However, little about the molecular biology of the peanut is currently known. Recently, next-generation sequencing technology, termed RNA-seq, has provided a powerful approach for analysing the transcriptome, and for shedding light on the molecular biology of peanut. RESULTS: In this study, we employed RNA-seq to analyse the transcriptomes of the immature seeds of three different peanut varieties with different oil contents. A total of 26.1-27.2 million paired-end reads with lengths of 100 bp were generated from the three varieties and 59,077 unigenes were assembled with N50 of 823 bp. Based on sequence similarity search with known proteins, a total of 40,100 genes were identified. Among these unigenes, only 8,252 unigenes were annotated with 42 gene ontology (GO) functional categories. And 18,028 unigenes mapped to 125 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, 3,919 microsatellite markers were developed in the unigene library, and 160 PCR primers of SSR loci were used for validation of the amplification and the polymorphism. CONCLUSION: We completed a successful global analysis of the peanut transcriptome using RNA-seq, a large number of unigenes were assembled, and almost four thousand SSR primers were developed. These data will facilitate gene discovery and functional genomic studies of the peanut plant. In addition, this study provides insight into the complex transcriptome of the peanut and established a biotechnological platform for future research.

New cis-regulatory elements in the Rht-D1b locus region of wheat.

Fifteen gene-containing BACs with accumulated length of 1.82-Mb from the Rht-D1b locus region were sequenced and compared in detail with the orthologous regions of rice, sorghum, and maize. Our results show that Rht-D1b represents a conserved genomic region as implied by high gene sequence identity, good maintenance of gene colinearity, and the presence of multiple conserved noncoding sequences (CNSs) that are shared by other grass species. Eight cis-regulatory elements in these CNSs around grass DELLA genes were detected.

A tandem segmental duplication (TSD) in green revolution gene Rht-D1b region underlies plant height variation.

• Rht-D1c (Rht10) carried by Chinese wheat (Triticum aestivum) line Aibian 1 is an allele at the Rht-D1 locus. Among the Rht-1 alleles, little is known about Rht-D1c although it determines an extreme dwarf phenotype in wheat. • Here, we cloned and functionally characterized Rht-D1c using a combination of Southern blotting, target region sequencing, gene expression analysis and transgenic experiments. • We found that the Rht-D1c allele was generated through a tandem segmental duplication (TSD) of a > 1 Mb region, resulting in two copies of the Rht-D1b. Two copies of Rht-D1b in the TSD were three-fold more effective in reducing plant height than a single copy, and transformation with a segment containing the tandemly duplicated copy of Rht-D1b resulted in the same level of reduction of plant height as the original copy in Aibian 1. • Our results suggest that changes in gene copy number are one of the important sources of genetic diversity and some of these changes could be directly associated with important traits in crops.

Author Correction: Using Gjd3-CreEGFP mice to examine atrioventricular node morphology and composition.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Global Analysis of Enhancer Targets Reveals Convergent Enhancer-Driven Regulatory Modules.

Single-cell screens enable high-throughput functional assessment of enhancers in their endogenous genomic context. However, the design of current studies limits their application to identifying the primary gene targets of enhancers. Here, we improve the experimental and computational parameters of single-cell enhancer screens to identify the secondary gene targets of enhancers. Our analysis of >500 putative enhancers in K562 cells reveals an interwoven enhancer-driven gene regulatory network. We find that enhancers from distinct genomic loci converge to modulate the expression of common sub-modules, including the α- and ß-globin loci, by directly regulating transcription factors. Our analysis suggests that several genetic variants associated with myeloid blood cell traits alter the activity of a distal enhancer of MYB (∼140 kb away), with downstream consequences on hemoglobin genes expression and cell state. These data have implications for the understanding of enhancer-associated traits and emphasize the flexibility of controlling transcriptional systems by modifying enhancer activity.

Rational Reprogramming of Cellular States by Combinatorial Perturbation.

Ectopic expression of transcription factors (TFs) can reprogram cell state. However, because of the large combinatorial space of possible TF cocktails, it remains difficult to identify TFs that reprogram specific cell types. Here, we develop Reprogram-Seq to experimentally screen thousands of TF cocktails for reprogramming performance. Reprogram-Seq leverages organ-specific cell-atlas data with single-cell perturbation and computational analysis to predict, evaluate, and optimize TF combinations that reprogram a cell type of interest. Focusing on the cardiac system, we perform Reprogram-Seq on MEFs using an undirected library of 48 cardiac factors and, separately, a directed library of 10 epicardial-related TFs. We identify a combination of three TFs, which efficiently reprogram MEFs to epicardial-like cells that are transcriptionally, molecularly, morphologically, and functionally similar to primary epicardial cells. Reprogram-Seq holds promise to accelerate the generation of specific cell types for regenerative medicine.

Using Gjd3-CreEGFP mice to examine atrioventricular node morphology and composition.

The atrioventricular node (AVN) coordinates the timing of atrial and ventricular contraction to optimize cardiac performance. To study this critical function using mouse genetics, however, new reagents are needed that allow AVN-specific manipulation. Here we describe a novel Gjd3-CreEGFP mouse line that successfully recombines floxed alleles within the AVN beginning at E12.5. These mice have been engineered to express CreEGFP under the control of endogenous Gjd3 regulatory elements without perturbing native protein expression. Detailed histological analysis of Gjd3-CreEGFP mice reveals specific labeling of AVN cardiomyocytes and a subset of cardiac endothelial cells. Importantly, we show that Gjd3-CreEGFP mice have preserved cardiac mechanical and electrical function. In one application of our newly described mouse line, we provide a three-dimensional (3D) view of the AVN using tissue clearing combined with confocal microscopy. With this 3D model as a reference, we identify specific AVN sub-structures based on marker staining characteristics. In addition, we use our Gjd3-CreEGFP mice to guide microdissection of the AVN and construction of a single-cell atlas. Thus, our results establish a new transgenic tool for AVN-specific recombination, provide an updated model of AVN morphology, and describe a roadmap for exploring AVN cellular heterogeneity.

DNA methylation reprogramming of functional elements during mammalian embryonic development.

DNA methylation plays important roles during development. However, the DNA methylation reprogramming of functional elements has not been fully investigated during mammalian embryonic development. Herein, using our modified MethylC-Seq library generation method and published post-bisulphite adapter-tagging (PBAT) method, we generated genome-wide DNA methylomes of human gametes and early embryos at single-base resolution and compared them with mouse methylomes. We showed that the dynamics of DNA methylation in functional elements are conserved between humans and mice during early embryogenesis, except for satellite repeats. We further found that oocyte-specific hypomethylated promoters usually exhibit low CpG densities. Genes with oocyte-specific hypomethylated promoters generally show oocyte-specific hypomethylated genic and intergenic regions, and these hypomethylated regions contribute to the hypomethylation pattern of mammalian oocytes. Furthermore, hypomethylated genic regions with low CG densities correlate with gene silencing in oocytes, whereas hypomethylated genic regions with high CG densities correspond to high gene expression. We further show that methylation reprogramming of enhancers during early embryogenesis is highly associated with the development of almost all human organs. Our data support the hypothesis that DNA methylation plays important roles during mammalian development.

A TRIM insertion in the promoter of Ms2 causes male sterility in wheat.

The male-sterile ms2 mutant has been known for 40 years and has become extremely important in the commercial production of wheat. However, the gene responsible for this phenotype has remained unknown. Here we report the map-based cloning of the Ms2 gene. The Ms2 locus is remarkable in several ways that have implications in basic biology. Beyond having no functional annotation, barely detectable transcription in fertile wild-type wheat plants, and accumulated destructive mutations in Ms2 orthologs, the Ms2 allele in the ms2 mutant has acquired a terminal-repeat retrotransposon in miniature (TRIM) element in its promoter. This TRIM element is responsible for the anther-specific Ms2 activation that confers male sterility. The identification of Ms2 not only unravels the genetic basis of a historically important breeding trait, but also shows an example of how a TRIM element insertion near a gene can contribute to genetic novelty and phenotypic plasticity.

Genome wide abnormal DNA methylome of human blastocyst in assisted reproductive technology.

Proper reprogramming of parental DNA methylomes is essential for mammalian embryonic development. However, it is unknown whether abnormal methylome reprogramming occurs and is associated with the failure of embryonic development. Here we analyzed the DNA methylomes of 57 blastocysts and 29 trophectoderm samples with different morphological grades during assisted reproductive technology (ART) practices. Our data reveal that the global methylation levels of high-quality blastocysts are similar (0.30 ± 0.02, mean ± SD), while the methylation levels of low-quality blastocysts are divergent and away from those of high-quality blastocysts. The proportion of blastocysts with a methylation level falling within the range of 0.30 ± 0.02 in different grades correlates with the live birth rate for that grade. Moreover, abnormal methylated regions are associated with the failure of embryonic development. Furthermore, we can use the methylation data of cells biopsied from trophectoderm to predict the blastocyst methylation level as well as to detect the aneuploidy of the blastocysts. Our data indicate that global abnormal methylome reprogramming often occurs in human embryos, and suggest that DNA methylome is a potential biomarker in blastocyst selection in ART.

Conversion of Human Gastric Epithelial Cells to Multipotent Endodermal Progenitors using Defined Small Molecules.

Endodermal stem/progenitor cells have diverse potential applications in research and regenerative medicine, so a readily available source could have widespread uses. Here we describe derivation of human induced endodermal progenitor cells (hiEndoPCs) from gastrointestinal epithelial cells using a cocktail of defined small molecules along with support from tissue-specific mesenchymal feeders. The hiEndoPCs show clonal expansion in culture and give rise to hepatocytes, pancreatic endocrine cells, and intestinal epithelial cells when treated with defined soluble molecules directing differentiation. The hiEndoPC-derived hepatocytes are able to rescue liver failure in Fah-/-Rag2-/- mice after transplantation, and, unlike hESCs, transplanted hiEndoPCs do not give rise to teratomas. Since human gastric epithelial cells are readily available from donors of many ages, this conversion strategy can generate clonally expandable cell populations with a variety of potential applications, including personalized drug screening and therapeutic strategies for liver failure and diabetes.