RNA-Seq of amniotic fluid cell-free RNA: a discovery phase study of the pathophysiology of congenital cytomegalovirus infection.

BACKGROUND: Congenital cytomegalovirus infection is the most common perinatal infection and a significant cause of sensorineural hearing loss, cerebral palsy, and neurodevelopmental disability. There is a paucity of human gene expression studies examining the pathophysiology of cytomegalovirus infection. OBJECTIVE: This study aimed to perform a whole transcriptomic assessment of amniotic fluid from pregnancies with live fetuses to identify differentially expressed genes and enriched Gene Ontology categories associated with congenital cytomegalovirus infection. STUDY DESIGN: Amniotic fluid supernatant was prospectively collected from pregnant women undergoing amniocentesis for suspected congenital cytomegalovirus infection because of first-trimester maternal primary infection or ultrasound features suggestive of fetal infection. Women who had received therapy to prevent fetal infection were excluded. Congenital cytomegalovirus infection was diagnosed via viral polymerase chain reaction of amniotic fluid; cytomegalovirus-infected fetuses were paired with noninfected controls, matched for gestational age and fetal sex. Paired-end RNA sequencing was performed on amniotic fluid cell-free RNA with the Novaseq 6000 at a depth of 30 million reads per sample. Following quality control and filtering, reads were mapped to the human genome and counts summarized across genes. Differentially expressed genes were identified using 2 approaches: voomWithQualityWeights in conjunction with limma and RUVSeq with edgeR. Genes with a false discovery rate <0.05 were considered statistically significant. Differential exon use was analyzed using DEXSeq. Functional analysis was performed using gene set enrichment analysis and Ingenuity Pathway Analysis. Manual curation of differentially regulated genes was also performed. RESULTS: Amniotic fluid samples were collected from 50 women; 16 (32%) had congenital cytomegalovirus infection confirmed by polymerase chain reaction. After excluding 3 samples without matched controls, 13 cytomegalovirus-infected samples collected at 18 to 23 weeks and 13 cytomegalovirus-negative gestation-matched controls were submitted for RNA sequencing and analysis (N=26). Ten of the 13 pregnancies with cytomegalovirus-infected fetuses had amniocentesis because of serologic evidence of maternal primary infection with normal fetal ultrasound, and 3 had amniocentesis because of ultrasound abnormality suggestive of cytomegalovirus infection. Four cytomegalovirus-infected pregnancies ended in termination (n=3) or fetal death (n=1), and 9 resulted in live births. Pregnancy outcomes were available for 11 of the 13 cytomegalovirus-negative controls; all resulted in live births of clinically-well infants. Differential gene expression analysis revealed 309 up-regulated and 32 down-regulated genes in the cytomegalovirus-infected group compared with the cytomegalovirus-negative group. Gene set enrichment analysis showed significant enrichment of multiple Gene Ontology categories involving the innate immune response to viral infection and interferon signaling. Of the 32 significantly down-regulated genes, 8 were known to be involved in neurodevelopment and preferentially expressed by the brain. Six specific cellular restriction factors involved in host defense to cytomegalovirus infection were up-regulated in the cytomegalovirus-infected group. Ingenuity Pathway Analysis predicted the activation of pathways involved in progressive neurologic disease and inflammatory neurologic disease. CONCLUSION: In this next-generation sequencing study, we revealed new insights into the pathophysiology of congenital cytomegalovirus infection. These data on the up-regulation of the intraamniotic innate immune response to cytomegalovirus infection and the dysregulation of neurodevelopmental genes may inform future approaches to developing prognostic markers and assessing fetal responses to in utero therapy.

Patient-iPSC-Derived Kidney Organoids Show Functional Validation of a Ciliopathic Renal Phenotype and Reveal Underlying Pathogenetic Mechanisms.

Despite the increasing diagnostic rate of genomic sequencing, the genetic basis of more than 50% of heritable kidney disease remains unresolved. Kidney organoids differentiated from induced pluripotent stem cells (iPSCs) of individuals affected by inherited renal disease represent a potential, but unvalidated, platform for the functional validation of novel gene variants and investigation of underlying pathogenetic mechanisms. In this study, trio whole-exome sequencing of a prospectively identified nephronophthisis (NPHP) proband and her parents identified compound-heterozygous variants in IFT140, a gene previously associated with NPHP-related ciliopathies. IFT140 plays a key role in retrograde intraflagellar transport, but the precise downstream cellular mechanisms responsible for disease presentation remain unknown. A one-step reprogramming and gene-editing protocol was used to derive both uncorrected proband iPSCs and isogenic gene-corrected iPSCs, which were differentiated to kidney organoids. Proband organoid tubules demonstrated shortened, club-shaped primary cilia, whereas gene correction rescued this phenotype. Differential expression analysis of epithelial cells isolated from organoids suggested downregulation of genes associated with apicobasal polarity, cell-cell junctions, and dynein motor assembly in proband epithelial cells. Matrigel cyst cultures confirmed a polarization defect in proband versus gene-corrected renal epithelium. As such, this study represents a "proof of concept" for using proband-derived iPSCs to model renal disease and illustrates dysfunctional cellular pathways beyond the primary cilium in the setting of IFT140 mutations, which are established for other NPHP genotypes.

Direct reprogramming to human nephron progenitor-like cells using inducible piggyBac transposon expression of SNAI2-EYA1-SIX1.

All nephrons in the mammalian kidney arise from a transient nephron progenitor population that is lost close to the time of birth. The generation of new nephron progenitors and their maintenance in culture are central to the success of kidney regenerative strategies. Using a lentiviral screening approach, we previously generated a human induced nephron progenitor-like state in vitro using a pool of six transcription factors. Here, we sought to develop a more efficient approach for direct reprogramming of human cells that could be applied in vivo. PiggyBac transposons are a non-viral integrating gene delivery system that is suitable for in vivo use and allows for simultaneous delivery of multiple genes. Using an inducible piggyBac transposon system, we optimized a protocol for the direct reprogramming of HK2 cells to induced nephron progenitor-like cells with expression of only 3 transcription factors (SNAI2, EYA1, and SIX1). Culture in conditions supportive of the nephron progenitor state further increased the expression of nephron progenitor genes. The refined protocol was then applied to primary human renal epithelial cells, which integrated into developing nephron structures in vitro and in vivo. Such inducible reprogramming to nephron progenitor-like cells could facilitate direct cellular reprogramming for kidney regeneration.

missMethyl: an R package for analyzing data from Illumina's HumanMethylation450 platform.

UNLABELLED: DNA methylation is one of the most commonly studied epigenetic modifications due to its role in both disease and development. The Illumina HumanMethylation450 BeadChip is a cost-effective way to profile >450 000 CpGs across the human genome, making it a popular platform for profiling DNA methylation. Here we introduce missMethyl, an R package with a suite of tools for performing normalization, removal of unwanted variation in differential methylation analysis, differential variability testing and gene set analysis for the 450K array. AVAILABILITY AND IMPLEMENTATION: missMethyl is an R package available from the Bioconductor project at www.bioconductor.org. CONTACT: alicia.oshlack@mcri.edu.au SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Removing unwanted variation in a differential methylation analysis of Illumina HumanMethylation450 array data.

Due to their relatively low-cost per sample and broad, gene-centric coverage of CpGs across the human genome, Illumina's 450k arrays are widely used in large scale differential methylation studies. However, by their very nature, large studies are particularly susceptible to the effects of unwanted variation. The effects of unwanted variation have been extensively documented in gene expression array studies and numerous methods have been developed to mitigate these effects. However, there has been much less research focused on the appropriate methodology to use for accounting for unwanted variation in methylation array studies. Here we present a novel 2-stage approach using RUV-inverse in a differential methylation analysis of 450k data and show that it outperforms existing methods.

The polycomb repressive complex 2 governs life and death of peripheral T cells.

Differentiation of naïve CD4(+) T cells into effector (Th1, Th2, and Th17) and induced regulatory (iTreg) T cells requires lineage-specifying transcription factors and epigenetic modifications that allow appropriate repression or activation of gene transcription. The epigenetic silencing of cytokine genes is associated with the repressive H3K27 trimethylation mark, mediated by the Ezh2 or Ezh1 methyltransferase components of the polycomb repressive complex 2 (PRC2). Here we show that silencing of the Ifng, Gata3, and Il10 loci in naïve CD4(+) T cells is dependent on Ezh2. Naïve CD4(+) T cells lacking Ezh2 were epigenetically primed for overproduction of IFN-γ in Th2 and iTreg and IL-10 in Th2 cells. In addition, deficiency of Ezh2 accelerated effector Th cell death via death receptor-mediated extrinsic and intrinsic apoptotic pathways, confirmed in vivo for Ezh2-null IFN-γ-producing CD4(+) and CD8(+) T cells responding to Listeria monocytogenes infection. These findings demonstrate the key role of PRC2/Ezh2 in differentiation and survival of peripheral T cells and reveal potential immunotherapeutic targets.

Genome-wide DNA methylation analysis identifies hypomethylated genes regulated by FOXP3 in human regulatory T cells.

Regulatory T cells (Treg) prevent the emergence of autoimmune disease. Prototypic natural Treg (nTreg) can be reliably identified by demethylation at the Forkhead-box P3 (FOXP3) locus. To explore the methylation landscape of nTreg, we analyzed genome-wide methylation in human naive nTreg (rTreg) and conventional naive CD4(+) T cells (Naive). We detected 2315 differentially methylated cytosine-guanosine dinucleotides (CpGs) between these 2 cell types, many of which clustered into 127 regions of differential methylation (RDMs). Activation changed the methylation status of 466 CpGs and 18 RDMs in Naive but did not alter DNA methylation in rTreg. Gene-set testing of the 127 RDMs showed that promoter methylation and gene expression were reciprocally related. RDMs were enriched for putative FOXP3-binding motifs. Moreover, CpGs within known FOXP3-binding regions in the genome were hypomethylated. In support of the view that methylation limits access of FOXP3 to its DNA targets, we showed that increased expression of the immune suppressive receptor T-cell immunoglobulin and immunoreceptor tyrosine-based inhibitory motif domain (TIGIT), which delineated Treg from activated effector T cells, was associated with hypomethylation and FOXP3 binding at the TIGIT locus. Differential methylation analysis provides insight into previously undefined human Treg signature genes and their mode of regulation.

Stability of gene expression and epigenetic profiles highlights the utility of patient-derived paediatric acute lymphoblastic leukaemia xenografts for investigating molecular mechanisms of drug resistance.

BACKGROUND: Patient-derived tumour xenografts are an attractive model for preclinical testing of anti-cancer drugs. Insights into tumour biology and biomarkers predictive of responses to chemotherapeutic drugs can also be gained from investigating xenograft models. As a first step towards examining the equivalence of epigenetic profiles between xenografts and primary tumours in paediatric leukaemia, we performed genome-scale DNA methylation and gene expression profiling on a panel of 10 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) tumours that were stratified by prednisolone response. RESULTS: We found high correlations in DNA methylation and gene expression profiles between matching primary and xenograft tumour samples with Pearson's correlation coefficients ranging between 0.85 and 0.98. In order to demonstrate the potential utility of epigenetic analyses in BCP-ALL xenografts, we identified DNA methylation biomarkers that correlated with prednisolone responsiveness of the original tumour samples. Differential methylation of CAPS2, ARHGAP21, ARX and HOXB6 were confirmed by locus specific analysis. We identified 20 genes showing an inverse relationship between DNA methylation and gene expression in association with prednisolone response. Pathway analysis of these genes implicated apoptosis, cell signalling and cell structure networks in prednisolone responsiveness. CONCLUSIONS: The findings of this study confirm the stability of epigenetic and gene expression profiles of paediatric BCP-ALL propagated in mouse xenograft models. Further, our preliminary investigation of prednisolone sensitivity highlights the utility of mouse xenograft models for preclinical development of novel drug regimens with parallel investigation of underlying gene expression and epigenetic responses associated with novel drug responses.

Benchmarking single-cell hashtag oligo demultiplexing methods.

Sample multiplexing is often used to reduce cost and limit batch effects in single-cell RNA sequencing (scRNA-seq) experiments. A commonly used multiplexing technique involves tagging cells prior to pooling with a hashtag oligo (HTO) that can be sequenced along with the cells' RNA to determine their sample of origin. Several tools have been developed to demultiplex HTO sequencing data and assign cells to samples. In this study, we critically assess the performance of seven HTO demultiplexing tools: hashedDrops, HTODemux, GMM-Demux, demuxmix, deMULTIplex, BFF (bimodal flexible fitting) and HashSolo. The comparison uses data sets where each sample has also been demultiplexed using genetic variants from the RNA, enabling comparison of HTO demultiplexing techniques against complementary data from the genetic 'ground truth'. We find that all methods perform similarly where HTO labelling is of high quality, but methods that assume a bimodal count distribution perform poorly on lower quality data. We also suggest heuristic approaches for assessing the quality of HTO counts in an scRNA-seq experiment.

Conservation of the ST6Gal I gene and its expression in the mammary gland.

Milk sialoglycoconjugates can protect the gastrointestinal tract of the suckling neonate by competitively binding to invading pathogens and promoting growth of beneficial flora, and their potential role in postnatal brain development is of particular interest in human infant nutrition. Although the concentration and the distribution of sialoglycoconjugates have been extensively studied in the milk of various species, the investigation of sialyltransferase gene expression in the mammary gland, in the context of lactation, has been limited. The sialyltransferase enzyme ST6Gal I transfers sialic acid from CMP-sialic acid to type 2 (Galß1,4GlcNAc) free disaccharides or the termini of N- or O-linked oligosaccharides using an α2,6-linkage. Expression of the ST6Gal I gene is primarily regulated at the level of transcription through the use of several cell and development-specific promoters, producing transcripts with divergent 5' untranslated regions (UTR). In the mouse mammary gland, the novel 5'UTR exon (L) appears to be associated with a drastic increase in ST6Gal I gene expression during lactation. We find that rats also possess an exon (L), suggesting conservation of this regulatory mechanism in rodents. In contrast, an exon (L)-containing transcript was not detected in the lactating bovine or human mammary gland. We also observed a trend of increasing ST6Gal I gene expression in the bovine mammary gland, culminating in involution. This is in contrast to species such as mice where the greatest change in ST6Gal I gene expression occurs between pregnancy and lactation, suggesting different roles in rodents vs. other mammals for α2,6-sialylated oligosaccharides present in milk.

Gene set enrichment analysis for genome-wide DNA methylation data.

DNA methylation is one of the most commonly studied epigenetic marks, due to its role in disease and development. Illumina methylation arrays have been extensively used to measure methylation across the human genome. Methylation array analysis has primarily focused on preprocessing, normalization, and identification of differentially methylated CpGs and regions. GOmeth and GOregion are new methods for performing unbiased gene set testing following differential methylation analysis. Benchmarking analyses demonstrate GOmeth outperforms other approaches, and GOregion is the first method for gene set testing of differentially methylated regions. Both methods are publicly available in the missMethyl Bioconductor R package.

DNA Methylation Profiles of Purified Cell Types in Bronchoalveolar Lavage: Applications for Mixed Cell Paediatric Pulmonary Studies.

In epigenome-wide association studies analysing DNA methylation from samples containing multiple cell types, it is essential to adjust the analysis for cell type composition. One well established strategy for achieving this is reference-based cell type deconvolution, which relies on knowledge of the DNA methylation profiles of purified constituent cell types. These are then used to estimate the cell type proportions of each sample, which can then be incorporated to adjust the association analysis. Bronchoalveolar lavage is commonly used to sample the lung in clinical practice and contains a mixture of different cell types that can vary in proportion across samples, affecting the overall methylation profile. A current barrier to the use of bronchoalveolar lavage in DNA methylation-based research is the lack of reference DNA methylation profiles for each of the constituent cell types, thus making reference-based cell composition estimation difficult. Herein, we use bronchoalveolar lavage samples collected from children with cystic fibrosis to define DNA methylation profiles for the four most common and clinically relevant cell types: alveolar macrophages, granulocytes, lymphocytes and alveolar epithelial cells. We then demonstrate the use of these methylation profiles in conjunction with an established reference-based methylation deconvolution method to estimate the cell type composition of two different tissue types; a publicly available dataset derived from artificial blood-based cell mixtures and further bronchoalveolar lavage samples. The reference DNA methylation profiles developed in this work can be used for future reference-based cell type composition estimation of bronchoalveolar lavage. This will facilitate the use of this tissue in studies examining the role of DNA methylation in lung health and disease.

DNA methylation biomarkers of future health outcomes in children.

Biomarkers which predict future health outcomes are key to the goals of precision health. Such biomarkers do not have to be involved in the causal pathway of a disease, and their performance is best assessed using statistical tests of clinical performance and evaluation of net health impact. DNA methylation is the most commonly studied epigenetic process and represents a potential biomarker of future health outcomes. We review 25 studies in non-oncological paediatric conditions where DNA methylation biomarkers of future health outcomes are assessed. Whilst a number of positive findings have been described, the body of evidence is severely limited by issues with outcome measures, tissue-specific samples, accounting for sample cell type heterogeneity, lack of appropriate statistical testing, small effect sizes, limited validation, and no assessment of net health impact. Future studies should concentrate on careful study design to overcome these issues, and integration of DNA methylation data with other 'omic', clinical, and environmental data to generate the most clinically useful biomarkers of paediatric disease.

Knockdown of a disintegrin A metalloprotease 12 (ADAM12) during adipogenesis reduces cell numbers, delays differentiation, and increases lipid accumulation in 3T3-L1 cells.

Mouse models have shown that a disintegrin A metalloprotease 12 (ADAM12) is implicated during adipogenesis; the molecular pathways are not well understood. Stealth RNA interference was used to knock down ADAM12 in 3T3-L1 cells. Using gene profiling and metabolic enzymatic markers, we have identified signaling pathways ADAM12 impacts upon during proliferation, differentiation, and maturation of adipocytes. ADAM12 reduced cell numbers in proliferating preadipocytes, delayed differentiation of preadipocytes to adipocytes, and increased lipid accumulation in mature adipocytes. The pathway most affected by ADAM12 knockdown was regulation of insulin-like growth factor (IGF) activity by insulin-like growth factor binding proteins (IGFBPs); ADAM12 is known to cleave IGFBP3 and IGFBP5. The IGF/mTOR signaling pathway was down-regulated, supporting a role for ADAM12 in the IGFBP/IGF/mTOR-growth pathway. PPARγ signaling was also down-regulated by ADAM12 knockdown. Gene ontology (GO) analysis revealed that the extracellular matrix was the cellular compartment most impacted. Filtering for matrisome genes, connective tissue growth factor ( Ctgf) was up-regulated. CTGF and IGBP3 can interact with PPARγ to hinder its regulation. Increased expression of these molecules could have influenced PPARγ signaling reducing differentiation and an imbalance of lipids. We believe ADAM12 regulates cell proliferation of preadipocytes through IGFBP/IGF/mTOR signaling and delays differentiation through altered PPAR signaling to cause an imbalance of lipids within mature adipocytes.

Gestational diabetes and maternal obesity are associated with epigenome-wide methylation changes in children.

Offspring of women with gestational diabetes mellitus (GDM) are at increased risk of developing metabolic disease, potentially mediated by epigenetic mechanisms. We recruited 608 GDM and 626 control offspring from the Danish National Birth Cohort, aged between 9 and 16 years. DNA methylation profiles were measured in peripheral blood of 93 GDM offspring and 95 controls using the Illumina HumanMethylation450 BeadChip. Pyrosequencing was performed for validation/replication of putative GDM-associated, differentially methylated CpGs in additional 905 offspring (462 GDM, 444 control offspring). We identified 76 differentially methylated CpGs in GDM offspring compared with controls in the discovery cohort (FDR, P < 0.05). Adjusting for offspring BMI did not affect the association between methylation levels and GDM status for any of the 76 CpGs. Most of these epigenetic changes were due to confounding by maternal prepregnancy BMI; however, 13 methylation changes were independently associated with maternal GDM. Three prepregnancy BMI-associated CpGs (cg00992687 and cg09452568 of ESM1 and cg14328641 of MS4A3) were validated in the replication cohort, while cg09109411 (PDE6A) was found to be associated with GDM status. The identified methylation changes may reflect developmental programming of organ disease mechanisms and/or may serve as disease biomarkers.

Epigenome-wide analysis in newborn blood spots from monozygotic twins discordant for cerebral palsy reveals consistent regional differences in DNA methylation.

Background: Cerebral palsy (CP) is a clinical description for a group of motor disorders that are heterogeneous with respect to causes, symptoms and severity. A diagnosis of CP cannot usually be made at birth and in some cases may be delayed until 2-3 years of age. This limits opportunities for early intervention that could otherwise improve long-term outcomes. CP has been recorded in monozygotic twins discordant for the disorder, indicating a potential role of non-genetic factors such as intrauterine infection, hypoxia-ischaemia, haemorrhage and thrombosis. The aim of this exploratory study was to utilise the discordant monozygotic twin model to understand and measure epigenetic changes associated with the development of CP. Methods: We performed a genome-wide analysis of DNA methylation using the Illumina Infinium Human Methylation 450 BeadChip array with DNA from newborn blood spots of 15 monozygotic twin pairs who later became discordant for CP. Quality control and data preprocessing were undertaken using the minfi R package. Differential methylation analysis was performed using the remove unwanted variation (RUVm) method, taking twin pairing into account in order to identify CP-specific differentially methylated probes (DMPs), and bumphunter was performed to identify differentially methylated regions (DMRs). Results: We identified 33 top-ranked DMPs based on a nominal p value cut-off of p < 1 × 10-4 and two DMRs (p < 1 × 10-3) associated with CP. The top-ranked probes related to 25 genes including HNRNPL, RASSF5, CD3D and KALRN involved in immune signalling pathways, in addition to TBC1D24, FBXO9 and VIPR2 previously linked to epileptic encephalopathy. Gene ontology and pathway analysis of top-ranked DMP-associated genes revealed enrichment of inflammatory signalling pathways, regulation of cytokine secretion and regulation of leukocyte-mediated immunity. We also identified two top-ranked DMRs including one on chromosome 6 within the promoter region of LTA gene encoding tumour necrosis factor-beta (TNF-ß), an important regulator of inflammation and brain development. The second was within the transcription start site of the LIME1 gene, which plays a key role in inflammatory pathways such as MAPK signalling. CP-specific differential DNA methylation within one of our two top DMRs was validated using an independent platform, MassArray EpiTyper. Conclusions: Ours is the first epigenome-wide association study of CP in disease-discordant monozygotic twin pairs and suggests a potential role for immune dysfunction in this condition.

Cord Blood CD8+ T Cells Have a Natural Propensity to Express IL-4 in a Fatty Acid Metabolism and Caspase Activation-Dependent Manner.

How T cells differentiate in the neonate may critically determine the ability of the infant to cope with infections, respond to vaccines and avert allergies. Previously, we found that naïve cord blood CD4+ T cells differentiated toward an IL-4-expressing phenotype when activated in the presence of TGF-ß and monocyte-derived inflammatory cytokines, the latter are more highly secreted by infants who developed food allergy. Here, we show that in the absence of IL-2 or IL-12, naïve cord blood CD8+ T cells have a natural propensity to differentiate into IL-4-producing non-classic TC2 cells when they are activated alone, or in the presence of TGF-ß and/or inflammatory cytokines. Mechanistically, non-classic TC2 development is associated with decreased expression of IL-2 receptor alpha (CD25) and glycolysis, and increased fatty acid metabolism and caspase-dependent cell death. Consequently, the short chain fatty acid, sodium propionate (NaPo), enhanced IL-4 expression, but exogenous IL-2 or pan-caspase inhibition prevented IL-4 expression. In children with endoscopically and histologically confirmed non-inflammatory bowel disease and non-infectious pediatric idiopathic colitis, the presence of TGF-ß, NaPo, and IL-1ß or TNF-α promoted TC2 differentiation in vitro. In vivo, colonic mucosa of children with colitis had significantly increased expression of IL-4 in CD8+ T cells compared with controls. In addition, activated caspase-3 and IL-4 were co-expressed in CD8+ T cells in the colonic mucosa of children with colitis. Thus, in the context of colonic inflammation and limited IL-2 signaling, CD8+ T cells differentiate into non-classic TC2 that may contribute to the pathology of inflammatory/allergic diseases in children.

Epigenetic interplay at the ß-globin locus.

During development, the α- and ß-globin genes exhibit a highly conserved pattern of expression, giving rise to several developmental stage-specific hemoglobin variants. Networks of regulatory proteins interact with epigenetic complexes to regulate DNA accessibility and histone modifications, thereby determining appropriate patterns of globin gene expression. In this review, we focus on recent advances in the understanding of the molecular mechanisms that underpin globin gene expression, focusing on multi-subunit regulatory complexes that bind to specific regions of DNA to orchestrate globin gene transcription throughout development.

A cross-package Bioconductor workflow for analysing methylation array data.

Methylation in the human genome is known to be associated with development and disease. The Illumina Infinium methylation arrays are by far the most common way to interrogate methylation across the human genome. This paper provides a Bioconductor workflow using multiple packages for the analysis of methylation array data. Specifically, we demonstrate the steps involved in a typical differential methylation analysis pipeline including: quality control, filtering, normalization, data exploration and statistical testing for probe-wise differential methylation. We further outline other analyses such as differential methylation of regions, differential variability analysis, estimating cell type composition and gene ontology testing. Finally, we provide some examples of how to visualise methylation array data.

Differentiation of human embryonic stem cells to HOXA+ hemogenic vasculature that resembles the aorta-gonad-mesonephros.

The ability to generate hematopoietic stem cells from human pluripotent cells would enable many biomedical applications. We find that hematopoietic CD34+ cells in spin embryoid bodies derived from human embryonic stem cells (hESCs) lack HOXA expression compared with repopulation-competent human cord blood CD34+ cells, indicating incorrect mesoderm patterning. Using reporter hESC lines to track the endothelial (SOX17) to hematopoietic (RUNX1C) transition that occurs in development, we show that simultaneous modulation of WNT and ACTIVIN signaling yields CD34+ hematopoietic cells with HOXA expression that more closely resembles that of cord blood. The cultures generate a network of aorta-like SOX17+ vessels from which RUNX1C+ blood cells emerge, similar to hematopoiesis in the aorta-gonad-mesonephros (AGM). Nascent CD34+ hematopoietic cells and corresponding cells sorted from human AGM show similar expression of cell surface receptors, signaling molecules and transcription factors. Our findings provide an approach to mimic in vitro a key early stage in human hematopoiesis for the generation of AGM-derived hematopoietic lineages from hESCs.