Lavinia Gordon

SWAN: Subset-quantile Within Array Normalization for Illumina Infinium HumanMethylation450 BeadChips

DNA methylation is the most widely studied epigenetic mark and is known to be essential to normal development and frequently disrupted in disease. The Illumina HumanMethylation450 BeadChip assays the methylation status of CpGs at 485,577 sites across the genome. Here we present Subset-quantile Within Array Normalization (SWAN), a new method that substantially improves the results from this platform by reducing technical variation within and between arrays. SWAN is available in the minfi Bioconductor package.

Neonatal DNA methylation profile in human twins is specified by a complex interplay between intrauterine environmental and genetic factors, subject to tissue-specific influence

Comparison between groups of monozygotic (MZ) and dizygotic (DZ) twins enables an estimation of the relative contribution of genetic and shared and nonshared environmental factors to phenotypic variability. Using DNA methylation profiling of ∼20,000 CpG sites as a phenotype, we have examined discordance levels in three neonatal tissues from 22 MZ and 12 DZ twin pairs. MZ twins exhibit a wide range of within-pair differences at birth, but show discordance levels generally lower than DZ pairs. Within-pair methylation discordance was lowest in CpG islands in all twins and increased as a function of distance from islands. Variance component decomposition analysis of DNA methylation in MZ and DZ pairs revealed a low mean heritability across all tissues, although a wide range of heritabilities was detected for specific genomic CpG sites. The largest component of variation was attributed to the combined effects of nonshared intrauterine environment and stochastic factors. Regression analysis of methylation on birth weight revealed a general association between methylation of genes involved in metabolism and biosynthesis, providing further support for epigenetic change in the previously described link between low birth weight and increasing risk for cardiovascular, metabolic, and other complex diseases. Finally, comparison of our data with that of several older twins revealed little evidence for genome-wide epigenetic drift with increasing age. This is the first study to analyze DNA methylation on a genome scale in twins at birth, further highlighting the importance of the intrauterine environment on shaping the neonatal epigenome.

Longitudinal, genome-scale analysis of DNA methylation in twins from birth to 18 months of age reveals rapid epigenetic change in early life and pair-specific effects of discordance

BackgroundThe extent to which development- and age-associated epigenetic changes are influenced by genetic, environmental and stochastic factors remains to be discovered. Twins provide an ideal model with which to investigate these influences but previous cross-sectional twin studies provide contradictory evidence of within-pair epigenetic drift over time. Longitudinal twin studies can potentially address this discrepancy.ResultsIn a pilot, genome-scale study of DNA from buccal epithelium, a relatively homogeneous tissue, we show that one-third of the CpGs assayed show dynamic methylation between birth and 18 months. Although all classes of annotated genomic regions assessed show an increase in DNA methylation over time, probes located in intragenic regions, enhancers and low-density CpG promoters are significantly over-represented, while CpG islands and high-CpG density promoters are depleted among the most dynamic probes. Comparison of co-twins demonstrated that within-pair drift in DNA methylation in our cohort is specific to a subset of pairs, who show more differences at 18 months. The rest of the pairs show either minimal change in methylation discordance, or more similar, converging methylation profiles at 18 months. As with age-associated regions, sites that change in their level of within-pair discordance between birth and 18 months are enriched in genes involved in development, but the average magnitude of change is smaller than for longitudinal change.ConclusionsOur findings suggest that DNA methylation in buccal epithelium is influenced by non-shared stochastic and environmental factors that could reflect a degree of epigenetic plasticity within an otherwise constrained developmental program.

Evidence for widespread changes in promoter methylation profile in human placenta in response to increasing gestational age and environmental/stochastic factors

BackgroundThe human placenta facilitates the exchange of nutrients, gas and waste between the fetal and maternal circulations. It also protects the fetus from the maternal immune response. Due to its role at the feto-maternal interface, the placenta is subject to many environmental exposures that can potentially alter its epigenetic profile. Previous studies have reported gene expression differences in placenta over gestation, as well as inter-individual variation in expression of some genes. However, the factors contributing to this variation in gene expression remain poorly understood.ResultsIn this study, we performed a genome-wide DNA methylation analysis of gene promoters in placenta tissue from three pregnancy trimesters. We identified large-scale differences in DNA methylation levels between first, second and third trimesters, with an overall progressive increase in average methylation from first to third trimester. The most differentially methylated genes included many immune regulators, reflecting the change in placental immuno-modulation as pregnancy progresses. We also detected increased inter-individual variation in the third trimester relative to first and second, supporting an accumulation of environmentally induced (or stochastic) changes in DNA methylation pattern. These highly variable genes were enriched for those involved in amino acid and other metabolic pathways, potentially reflecting the adaptation of the human placenta to different environments.ConclusionsThe identification of cellular pathways subject to drift in response to environmental influences provide a basis for future studies examining the role of specific environmental factors on DNA methylation pattern and placenta-associated adverse pregnancy outcomes.

Copy Number Variation in Patients with Disorders of Sex Development Due to 46,XY Gonadal Dysgenesis

Disorders of sex development (DSD), ranging in severity from mild genital abnormalities to complete sex reversal, represent a major concern for patients and their families. DSD are often due to disruption of the genetic programs that regulate gonad development. Although some genes have been identified in these developmental pathways, the causative mutations have not been identified in more than 50% 46,XY DSD cases. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to analyse copy number variation in 23 individuals with unexplained 46,XY DSD due to gonadal dysgenesis (GD). Here we describe three discrete changes in copy number that are the likely cause of the GD. Firstly, we identified a large duplication on the X chromosome that included DAX1 (NR0B1). Secondly, we identified a rearrangement that appears to affect a novel gonad-specific regulatory region in a known testis gene, SOX9. Surprisingly this patient lacked any signs of campomelic dysplasia, suggesting that the deletion affected expression of SOX9 only in the gonad. Functional analysis of potential SRY binding sites within this deleted region identified five putative enhancers, suggesting that sequences additional to the known SRY-binding TES enhancer influence human testis-specific SOX9 expression. Thirdly, we identified a small deletion immediately downstream of GATA4, supporting a role for GATA4 in gonad development in humans. These CNV analyses give new insights into the pathways involved in human gonad development and dysfunction, and suggest that rearrangements of non-coding sequences disturbing gene regulation may account for significant proportion of DSD cases.

Evidence for age-related and individual-specific changes in DNA methylation profile of mononuclear cells during early immune development in humans

Environment induced epigenetic effects on gene expression in early life are likely to play important roles in mediating the risk of several immune-related diseases. In order to investigate this fully, it is essential to first document temporal changes in epigenetic profile in disease-free individuals as a prelude to defining environmentally mediated changes. Mononuclear cells (MC) were collected longitudinally from a small number of females at birth, 1 year, 2.5 years and 5 years of age and examined for changes in genome-scale DNA methylation profiles using the Illumina Infinium HumanMethylation27 BeadChip array platform. MC from two males were included for comparative purposes. Flow cytometry was used to define MC cell populations in each sample in order to exclude this as the major driver of epigenetic change. The data underwent quality control and normalization within the R programming environment. Unsupervised hierarchical clustering of samples clearly delineated neonatal MC from all other ages. A further clear distinction was observed between 1 year and 5 year samples, with 2.5 year samples showing a mixed distribution between the 1 and 5 year groups. Gene ontology of probes significantly variable over the neonatal period revealed methylation changes in genes associated with cell surface receptor and signal transduction events. In the postnatal period, methylation changes were mostly associated with the development of effector immune responses and homeostasis. Unlike all other chromosomes tested, a predominantly genetic effect was identified as controlling maintenance of X-chromosome methylation profile in females, largely refractory to change over time. This data suggests that the primary driver of neonatal epigenome is determined in utero, whilst postnatally, multiple genetic and environmental factors are implicated in the development of MC epigenetic profile, particularly between the ages of 1–5 years, when the highest level of inter individual variation is apparent. This supports a model for differential sensitivity of specific individuals to disruption in the developing epigenome during the first years of life. Further studies are now needed to examine evolving epigenetic variations in specific cell populations in relation to environmental exposures, immune phenotype and subsequent disease susceptibility.

Disruption of the histone acetyltransferase MYST4 leads to a Noonan syndrome–like phenotype and hyperactivated MAPK signaling in humans and mice

Epigenetic regulation of gene expression, through covalent modification of histones, is a key process controlling growth and development. Accordingly, the transcription factors regulating these processes are important targets of genetic diseases. However, surprisingly little is known about the relationship between aberrant epigenetic states, the cellular process affected, and their phenotypic consequences. By chromosomal breakpoint mapping in a patient with a Noonan syndrome-like phenotype that encompassed short stature, blepharoptosis, and attention deficit hyperactivity disorder, we identified haploinsufficiency of the histone acetyltransferase gene MYST histone acetyltransferase (monocytic leukemia) 4 (MYST4), as the underlying cause of the phenotype. Using acetylation, whole genome expression, and ChIP studies in cells from the patient, cell lines in which MYST4 expression was knocked down using siRNA, and the Myst4 querkopf mouse, we found that H3 acetylation is important for neural, craniofacial, and skeletal morphogenesis, mainly through its ability to specifically regulating the MAPK signaling pathway. This finding further elucidates the complex role of histone modifications in mammalian development and adds what we believe to be a new mechanism to the pathogenic phenotypes resulting from misregulation of the RAS signaling pathway.

Changes in the Chondrocyte and Extracellular Matrix Proteome during Post-natal Mouse Cartilage Development

Skeletal growth by endochondral ossification involves tightly coordinated chondrocyte differentiation that creates reserve, proliferating, prehypertrophic, and hypertrophic cartilage zones in the growth plate. Many human skeletal disorders result from mutations in cartilage extracellular matrix (ECM) components that compromise both ECM architecture and chondrocyte function. Understanding normal cartilage development, composition, and structure is therefore vital to unravel these disease mechanisms. To study this intricate process in vivo by proteomics, we analyzed mouse femoral head cartilage at developmental stages enriched in either immature chondrocytes or maturing/hypertrophic chondrocytes (post-natal days 3 and 21, respectively). Using LTQ-Orbitrap tandem mass spectrometry, we identified 703 cartilage proteins. Differentially abundant proteins (q < 0.01) included prototypic markers for both early and late chondrocyte differentiation (epiphycan and collagen X, respectively) and novel ECM and cell adhesion proteins with no previously described roles in cartilage development (tenascin X, vitrin, Urb, emilin-1, and the sushi repeat-containing proteins SRPX and SRPX2). Meta-analysis of cartilage development in vivo and an in vitro chondrocyte culture model (Wilson, R., Diseberg, A. F., Gordon, L., Zivkovic, S., Tatarczuch, L., Mackie, E. J., Gorman, J. J., and Bateman, J. F. (2010) Comprehensive profiling of cartilage extracellular matrix formation and maturation using sequential extraction and label-free quantitative proteomics. Mol. Cell. Proteomics 9, 1296–1313) identified components involved in both systems, such as Urb, and components with specific roles in vivo, including vitrin and CILP-2 (cartilage intermediate layer protein-2). Immunolocalization of Urb, vitrin, and CILP-2 indicated specific roles at different maturation stages. In addition to ECM-related changes, we provide the first biochemical evidence of changing endoplasmic reticulum function during cartilage development. Although the multifunctional chaperone BiP was not differentially expressed, enzymes and chaperones required specifically for collagen biosynthesis, such as the prolyl 3-hydroxylase 1, cartilage-associated protein, and peptidyl prolyl cis-trans isomerase B complex, were down-regulated during maturation. Conversely, the lumenal proteins calumenin, reticulocalbin-1, and reticulocalbin-2 were significantly increased, signifying a shift toward calcium binding functions. This first proteomic analysis of cartilage development in vivo reveals the breadth of protein expression changes during chondrocyte maturation and ECM remodeling in the mouse femoral head.

Expression discordance of monozygotic twins at birth: Effect of intrauterine environment and a possible mechanism for fetal programming

Within-pair comparison of monozygotic (MZ) twins provides an ideal model for studying factors that regulate epigenetic profile, by controlling for genetic variation. Previous reports have demonstrated epigenetic variability within MZ pairs, but the contribution of early life exposures to this variation remains unclear. As epigenetic marks govern gene expression, we have used gene expression discordance as a proxy measure of epigenetic discordance in MZ twins at birth in two cell types. We found strong evidence of expression discordance at birth in both cell types and some evidence for higher discordance in twin pairs with separate placentas. Genes previously defined as being involved in response to the external environment showed the most variable expression within pairs, independent of cell type, supporting the idea that even slight differences in intrauterine environment can influence expression profile. Focusing on birthweight, previously identified as a predisposing factor for cardiovascular, metabolic and other complex diseases, and using a statistical model that estimated association based on within-pair variation of expression and birthweight, we found some association between birthweight and expression of genes involved in metabolism and cardiovascular function. This study is the first to examine expression discordance in newborn twins. It provides evidence of a link between birthweight and activity of specific cellular pathways and, as evidence points to gene expression profiles being maintained through cell division by epigenetic factors, provides a plausible biological mechanism for the previously described link between low birthweight and increased risk of later complex disease.

Comprehensive Profiling of Cartilage Extracellular Matrix Formation and Maturation Using Sequential Extraction and Label-free Quantitative Proteomics

Articular cartilage is indispensable for joint function but has limited capacity for self-repair. Engineering of neocartilage in vitro is therefore a major target for autologous cartilage repair in arthritis. Previous analysis of neocartilage has targeted cellular organization and specific molecular components. However, the complexity of extracellular matrix (ECM) development in neocartilage has not been investigated by proteomics. To redress this, we developed a mouse neocartilage culture system that produces a cartilaginous ECM. Differential analysis of the tissue proteome of 3-week neocartilage and 3-day postnatal mouse cartilage using solubility-based protein fractionation targeted components involved in neocartilage development, including ECM maturation. Initially, SDS-PAGE analysis of sequential extracts revealed the transition in protein solubility from a high proportion of readily soluble (NaCl-extracted) proteins in juvenile cartilage to a high proportion of poorly soluble (guanidine hydrochloride-extracted) proteins in neocartilage. Label-free quantitative mass spectrometry (LTQ-Orbitrap) and statistical analysis were then used to filter three significant protein groups: proteins enriched according to extraction condition, proteins differentially abundant between juvenile cartilage and neocartilage, and proteins with differential solubility properties between the two tissue types. Classification of proteins differentially abundant between NaCl and guanidine hydrochloride extracts (n = 403) using bioinformatics revealed effective partitioning of readily soluble components from subunits of larger protein complexes. Proteins significantly enriched in neocartilage (n = 78) included proteins previously not reported or with unknown function in cartilage (integrin-binding protein DEL1; coiled-coil domain-containing protein 80; emilin-1 and pigment epithelium derived factor). Proteins with differential extractability between juvenile cartilage and neocartilage included ECM components (nidogen-2, perlecan, collagen VI, matrilin-3, tenascin and thrombospondin-1), and the relationship between protein extractability and ECM ultrastructural organization was supported by electron microscopy. Additionally, one guanidine extract-specific neocartilage protein, protease nexin-1, was confirmed by immunohistochemistry as a novel component of developing articular cartilage in vivo. The extraction profile and matrix-associated immunostaining implicates protease nexin-1 in cartilage development in vitro and in vivo.