Understanding the consequences of regulatory variation in the human genome remains a major challenge, with important implications for understanding gene regulation and interpreting the many disease-risk variants that fall outside of protein-coding regions. Here, we provide a direct window into the regulatory consequences of genetic variation by sequencing RNA from 922 genotyped individuals. We present a comprehensive description of the distribution of regulatory variation--by the specific expression phenotypes altered, the properties of affected genes, and the genomic characteristics of regulatory variants. We detect variants influencing expression of over ten thousand genes, and through the enhanced resolution offered by RNA-sequencing, for the first time we identify thousands of variants associated with specific phenotypes including splicing and allelic expression. Evaluating the effects of both long-range intra-chromosomal and trans (cross-chromosomal) regulation, we observe modularity in the regulatory network, with three-dimensional chromosomal configuration playing a particular role in regulatory modules within each chromosome. We also observe a significant depletion of regulatory variants affecting central and critical genes, along with a trend of reduced effect sizes as variant frequency increases, providing evidence that purifying selection and buffering have limited the deleterious impact of regulatory variation on the cell. Further, generalizing beyond observed variants, we have analyzed the genomic properties of variants associated with expression and splicing and developed a Bayesian model to predict regulatory consequences of genetic variants, applicable to the interpretation of individual genomes and disease studies. Together, these results represent a critical step toward characterizing the complete landscape of human regulatory variation.

Upon vessel injury, platelets become activated and rapidly reorganize their actin cytoskeleton to adhere to the site of endothelial damage, triggering the formation of a fibrin-rich plug to prevent further blood loss. Inactivation of Pdlim7 provides the new perspective that regulation of actin cytoskeletal changes in platelets is dependent on the encoded PDZ-LIM protein. Loss-of-function of Pdlim7 triggers hypercoagulopathy and causes significant perinatal lethality in mice. Our in vivo and in vitro studies reveal that Pdlim7 is dynamically distributed along actin fibers, and lack of Pdlim7 leads to a marked inability to rearrange the actin cytoskeleton. Specifically, the absence of Pdlim7 prevents platelets from bundling actin fibers into a concentric ring that defines the round spread shape of activated platelets. Similarly, in mouse embryonic fibroblasts, loss of Pdlim7 abolishes the formation of stress fibers needed to adopt the typical elongated fibroblast shape. In addition to revealing a fundamental cell biological role in actin cytoskeletal organization, we also demonstrate a function of Pdlim7 in regulating the cycling between the GTP/GDP-bound states of Arf6. The small GTPase Arf6 is an essential factor required for actin dynamics, cytoskeletal rearrangements, and platelet activation. Consistent with our findings of significantly elevated initial F-actin ratios and subsequent morphological aberrations, loss of Pdlim7 causes a shift in balance towards an increased Arf6-GTP level in resting platelets. These findings identify a new Pdlim7-Arf6 axis controlling actin dynamics and implicate Pdlim7 as a primary endogenous regulator of platelet-dependent hemostasis.

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

As a consequence of the accumulation of insertion events over evolutionary time, mobile elements now comprise nearly half of the human genome. The Alu, L1, and SVA mobile element families are still duplicating, generating variation between individual genomes. Mobile element insertions (MEI) have been identified as causes for genetic diseases, including hemophilia, neurofibromatosis, and various cancers. Here we present a comprehensive map of 7,380 MEI polymorphisms from the 1000 Genomes Project whole-genome sequencing data of 185 samples in three major populations detected with two detection methods. This catalog enables us to systematically study mutation rates, population segregation, genomic distribution, and functional properties of MEI polymorphisms and to compare MEI to SNP variation from the same individuals. Population allele frequencies of MEI and SNPs are described, broadly, by the same neutral ancestral processes despite vastly different mutation mechanisms and rates, except in coding regions where MEI are virtually absent, presumably due to strong negative selection. A direct comparison of MEI and SNP diversity levels suggests a differential mobile element insertion rate among populations.

Genome amplification and cellular senescence are commonly associated with pathological processes. While physiological roles for polyploidization and senescence have been described in mouse development, controversy exists over their significance in humans. Here, we describe tetraploidization and senescence as phenomena of normal human placenta development. During pregnancy, placental extravillous trophoblasts (EVTs) invade the pregnant endometrium, termed decidua, to establish an adapted microenvironment required for the developing embryo. This process is critically dependent on continuous cell proliferation and differentiation, which is thought to follow the classical model of cell cycle arrest prior to terminal differentiation. Strikingly, flow cytometry and DNAseq revealed that EVT formation is accompanied with a genome-wide polyploidization, independent of mitotic cycles. DNA replication in these cells was analysed by a fluorescent cell-cycle indicator reporter system, cell cycle marker expression and EdU incorporation. Upon invasion into the decidua, EVTs widely lose their replicative potential and enter a senescent state characterized by high senescence-associated (SA) β-galactosidase activity, induction of a SA secretory phenotype as well as typical metabolic alterations. Furthermore, we show that the shift from endocycle-dependent genome amplification to growth arrest is disturbed in androgenic complete hydatidiform moles (CHM), a hyperplastic pregnancy disorder associated with increased risk of developing choriocarinoma. Senescence is decreased in CHM-EVTs, accompanied by exacerbated endoreduplication and hyperploidy. We propose induction of cellular senescence as a ploidy-limiting mechanism during normal human placentation and unravel a link between excessive polyploidization and reduced senescence in CHM.

DNA sequence information underpins genetic research, enabling discoveries of important biological or medical benefit. Sequencing projects have traditionally used long (400-800 base pair) reads, but the existence of reference sequences for the human and many other genomes makes it possible to develop new, fast approaches to re-sequencing, whereby shorter reads are compared to a reference to identify intraspecies genetic variation. Here we report an approach that generates several billion bases of accurate nucleotide sequence per experiment at low cost. Single molecules of DNA are attached to a flat surface, amplified in situ and used as templates for synthetic sequencing with fluorescent reversible terminator deoxyribonucleotides. Images of the surface are analysed to generate high-quality sequence. We demonstrate application of this approach to human genome sequencing on flow-sorted X chromosomes and then scale the approach to determine the genome sequence of a male Yoruba from Ibadan, Nigeria. We build an accurate consensus sequence from >30x average depth of paired 35-base reads. We characterize four million single-nucleotide polymorphisms and four hundred thousand structural variants, many of which were previously unknown. Our approach is effective for accurate, rapid and economical whole-genome re-sequencing and many other biomedical applications.

Although DNA methylation is a key regulator of gene expression, the comprehensive methylation landscape of metastatic cancer has never been defined. Through whole-genome bisulfite sequencing paired with deep whole-genome and transcriptome sequencing of 100 castration-resistant prostate metastases, we discovered alterations affecting driver genes that were detectable only with integrated whole-genome approaches. Notably, we observed that 22% of tumors exhibited a novel epigenomic subtype associated with hypermethylation and somatic mutations in TET2, DNMT3B, IDH1 and BRAF. We also identified intergenic regions where methylation is associated with RNA expression of the oncogenic driver genes AR, MYC and ERG. Finally, we showed that differential methylation during progression preferentially occurs at somatic mutational hotspots and putative regulatory regions. This study is a large integrated study of whole-genome, whole-methylome and whole-transcriptome sequencing in metastatic cancer that provides a comprehensive overview of the important regulatory role of methylation in metastatic castration-resistant prostate cancer.

The 15q13.3 microdeletion is associated with several neuropsychiatric disorders, including autism and schizophrenia. Previous association and functional studies have investigated the potential role of several genes within the deletion in neuronal dysfunction, but the molecular effects of the deletion as a whole remain largely unknown.

The early environment, including maternal characteristics, provides many cues to young organisms that shape their long-term physical and mental health. Identifying the earliest molecular events that precede observable developmental outcomes could help identify children in need of support prior to the onset of physical and mental health difficulties. In this study, we examined whether mothers' attachment insecurity, maltreatment history, and depressive symptoms were associated with alterations in DNA methylation patterns in their infants, and whether these correlates in the infant epigenome were associated with socioemotional and behavioral functioning in toddlerhood. We recruited 156 women oversampled for histories of depression, who completed psychiatric interviews and depression screening during pregnancy, then provided follow-up behavioral data on their children at 18 months. Buccal cell DNA was obtained from 32 of their infants for a large-scale analysis of methylation patterns across 5 × 106 individual CpG dinucleotides, using clustering-based significance criteria to control for multiple comparisons. We found that tens of thousands of individual infant CpGs were alternatively methylated in association with maternal attachment insecurity, maltreatment in childhood, and antenatal and postpartum depressive symptoms, including genes implicated in developmental patterning, cell-cell communication, hormonal regulation, immune function/inflammatory response, and neurotransmission. Density of DNA methylation at selected genes from the result set was also significantly associated with toddler socioemotional and behavioral problems. This is the first report to identify novel regions of the human infant genome at which DNA methylation patterns are associated longitudinally both with maternal characteristics and with offspring socioemotional and behavioral problems in toddlerhood.

Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells.

Structural variation in the complement 4 gene (C4) confers genetic risk for schizophrenia. The variation includes numbers of the increased C4A copy number, which predicts increased C4A mRNA expression. C4-anaphylatoxin (C4-ana) is a C4 protein fragment released upon C4 protein activation that has the potential to change the blood-brain barrier (BBB). We hypothesized that elevated plasma levels of C4-ana occur in individuals with schizophrenia (iSCZ). Blood was collected from 15 iSCZ with illness duration < 5 years and from 14 healthy controls (HC). Plasma C4-ana was measured by radioimmunoassay. Other complement activation products C3-ana, C5-ana, and terminal complement complex (TCC) were also measured. Digital-droplet PCR was used to determine C4 gene structural variation state. Recombinant C4-ana was added to primary brain endothelial cells (BEC) and permeability was measured in vitro. C4-ana concentration was elevated in plasma from iSCZ compared to HC (mean = 654 ± 16 ng/mL, 557 ± 94 respectively, p = 0.01). The patients also carried more copies of the C4AL gene and demonstrated a positive correlation between plasma C4-ana concentrations and C4A gene copy number. Furthermore, C4-ana increased the permeability of a monolayer of BEC in vitro. Our findings are consistent with a specific role for C4A protein in schizophrenia and raise the possibility that its activation product, C4-ana, increases BBB permeability. Exploratory analyses suggest the novel hypothesis that the relationship between C4-ana levels and C4A gene copy number could also be altered in iSCZ, suggesting an interaction with unknown genetic and/or environmental risk factors.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

DNA capture technologies combined with high-throughput sequencing now enable cost-effective, deep-coverage, targeted sequencing of complete exomes. This is well suited for SNP discovery and genotyping. However there has been little attention devoted to Copy Number Variation (CNV) detection from exome capture datasets despite the potentially high impact of CNVs in exonic regions on protein function.

Few studies have been conducted to understand post-zygotic accumulation of mutations in cells of the healthy human body. We reprogrammed 32 skin fibroblast cells from families of donors into human induced pluripotent stem cell (hiPSC) lines. The clonal nature of hiPSC lines allows a high-resolution analysis of the genomes of the founder fibroblast cells without being confounded by the artifacts of single-cell whole-genome amplification. We estimate that on average a fibroblast cell in children has 1035 mostly benign mosaic SNVs. On average, 235 SNVs could be directly confirmed in the original fibroblast population by ultradeep sequencing, down to an allele frequency (AF) of 0.1%. More sensitive droplet digital PCR experiments confirmed more SNVs as mosaic with AF as low as 0.01%, suggesting that 1035 mosaic SNVs per fibroblast cell is the true average. Similar analyses in adults revealed no significant increase in the number of SNVs per cell, suggesting that a major fraction of mosaic SNVs in fibroblasts arises during development. Mosaic SNVs were distributed uniformly across the genome and were enriched in a mutational signature previously observed in cancers and in de novo variants and which, we hypothesize, is a hallmark of normal cell proliferation. Finally, AF distribution of mosaic SNVs had distinct narrow peaks, which could be a characteristic of clonal cell selection, clonal expansion, or both. These findings reveal a large degree of somatic mosaicism in healthy human tissues, link de novo and cancer mutations to somatic mosaicism, and couple somatic mosaicism with cell proliferation.

Early life adversity and insecure attachment style are known risk factors for perinatal depression. The biological pathways linking these experiences, however, have not yet been elucidated. We hypothesized that overlap in patterns of DNA methylation in association with each of these phenomena could identify genes and pathways of importance. Specifically, we wished to distinguish between allostatic-load and role-transition hypotheses of perinatal depression. We conducted a large-scale analysis of methylation patterns across 5 × 106 individual CG dinucleotides in 54 women participating in a longitudinal prospective study of perinatal depression, using clustering-based criteria for significance to control for multiple comparisons. We identified 1580 regions in which methylation density was associated with childhood adversity, 3 in which methylation density was associated with insecure attachment style, and 6 in which methylation density was associated with perinatal depression. Shorter telomeres were observed in association with childhood trauma but not with perinatal depression or attachment insecurity. A detailed analysis of methylation density in the oxytocin receptor gene revealed similar patterns of DNA methylation in association with perinatal depression and with insecure attachment style, while childhood trauma was associated with a distinct methylation pattern in this gene. Clinically, attachment style was strongly associated with depression only in pregnancy and the early postpartum, whereas the association of childhood adversity with depression was time-invariant. We concluded that the broad DNA methylation signature and reduced telomere length associated with childhood adversity could indicate increased allostatic load across multiple body systems, whereas perinatal depression and attachment insecurity may be narrower phenotypes with more limited DNA methylation signatures outside the CNS, and no apparent association with telomere length or, by extension, allostatic load. In contrast, the finding of matching DNA methylation patterns within the oxytocin receptor gene for perinatal depression and attachment insecurity is consistent with the theory that the perinatal period is a time of activation of existing attachment schemas for the purpose of structuring the mother-child relationship, and that such activation may occur in part through specific patterns of methylation of the oxytocin receptor gene.

HepG2 is one of the most widely used human cancer cell lines in biomedical research and one of the main cell lines of ENCODE. Although the functional genomic and epigenomic characteristics of HepG2 are extensively studied, its genome sequence has never been comprehensively analyzed and higher order genomic structural features are largely unknown. The high degree of aneuploidy in HepG2 renders traditional genome variant analysis methods challenging and partially ineffective. Correct and complete interpretation of the extensive functional genomics data from HepG2 requires an understanding of the cell line's genome sequence and genome structure. Using a variety of sequencing and analysis methods, we identified a wide spectrum of genome characteristics in HepG2: copy numbers of chromosomal segments at high resolution, SNVs and Indels (corrected for aneuploidy), regions with loss of heterozygosity, phased haplotypes extending to entire chromosome arms, retrotransposon insertions and structural variants (SVs) including complex and somatic genomic rearrangements. A large number of SVs were phased, sequence assembled and experimentally validated. We re-analyzed published HepG2 datasets for allele-specific expression and DNA methylation and assembled an allele-specific CRISPR/Cas9 targeting map. We demonstrate how deeper insights into genomic regulatory complexity are gained by adopting a genome-integrated framework.

Somatic mosaicism, manifesting as single nucleotide variants (SNVs), mobile element insertions, and structural changes in the DNA, is a common phenomenon in human brain cells, with potential functional consequences. Using a clonal approach, we previously detected 200-400 mosaic SNVs per cell in three human fetal brains (15-21 wk postconception). However, structural variation in the human fetal brain has not yet been investigated. Here, we discover and validate four mosaic structural variants (SVs) in the same brains and resolve their precise breakpoints. The SVs were of kilobase scale and complex, consisting of deletion(s) and rearranged genomic fragments, which sometimes originated from different chromosomes. Sequences at the breakpoints of these rearrangements had microhomologies, suggesting their origin from replication errors. One SV was found in two clones, and we timed its origin to ∼14 wk postconception. No large scale mosaic copy number variants (CNVs) were detectable in normal fetal human brains, suggesting that previously reported megabase-scale CNVs in neurons arise at later stages of development. By reanalysis of public single nuclei data from adult brain neurons, we detected an extrachromosomal circular DNA event. Our study reveals the existence of mosaic SVs in the developing human brain, likely arising from cell proliferation during mid-neurogenesis. Although relatively rare compared to SNVs and present in ∼10% of neurons, SVs in developing human brain affect a comparable number of bases in the genome (∼6200 vs. ∼4000 bp), implying that they may have similar functional consequences.

The human pangenome, a new reference sequence, addresses many limitations of the current GRCh38 reference. The first release is based on 94 high-quality haploid assemblies from individuals with diverse backgrounds. We employed a k-mer indexing strategy for comparative analysis across multiple assemblies, including the pangenome reference, GRCh38, and CHM13, a telomere-to-telomere reference assembly. Our k-mer indexing approach enabled us to identify a valuable collection of universally conserved sequences across all assemblies, referred to as "pan-conserved segment tags" (PSTs). By examining intervals between these segments, we discerned highly conserved genomic segments and those with structurally related polymorphisms. We found 60,764 polymorphic intervals with unique geo-ethnic features in the pangenome reference. In this study, we utilized ultra-conserved sequences (PSTs) to forge a link between human pangenome assemblies and reference genomes. This methodology enables the examination of any sequence of interest within the pangenome, using the reference genome as a comparative framework.

Accurate and efficient genome-wide detection of copy number variants (CNVs) is essential for understanding human genomic variation, genome-wide CNV association type studies, cytogenetics research and diagnostics, and independent validation of CNVs identified from sequencing based technologies. Numerous, array-based platforms for CNV detection exist utilizing array Comparative Genome Hybridization (aCGH), Single Nucleotide Polymorphism (SNP) genotyping or both. We have quantitatively assessed the abilities of twelve leading genome-wide CNV detection platforms to accurately detect Gold Standard sets of CNVs in the genome of HapMap CEU sample NA12878, and found significant differences in performance. The technologies analyzed were the NimbleGen 4.2 M, 2.1 M and 3×720 K Whole Genome and CNV focused arrays, the Agilent 1×1 M CGH and High Resolution and 2×400 K CNV and SNP+CGH arrays, the Illumina Human Omni1Quad array and the Affymetrix SNP 6.0 array. The Gold Standards used were a 1000 Genomes Project sequencing-based set of 3997 validated CNVs and an ultra high-resolution aCGH-based set of 756 validated CNVs. We found that sensitivity, total number, size range and breakpoint resolution of CNV calls were highest for CNV focused arrays. Our results are important for cost effective CNV detection and validation for both basic and clinical applications.

Gene transcription mediates many vital aspects of mammalian embryonic development. A comprehensive characterization and analysis of the dynamics of gene transcription in the embryo is therefore likely to provide significant insights into the basic mechanisms of this process. We used microarrays to map transcription in the mouse embryo in the important period from embryonic day 8 (e8.0) to postnatal day 1 (p1) during which the bulk of the differentiation and development of organ systems takes place. Analysis of these expression profiles revealed distinct patterns of gene expression which correlate with the differentiation of organs including the nervous system, liver, skin, lungs, and digestive system, among others. Statistical analysis of the data based on Gene Ontology (GO) group annotation showed that specific temporal sequence patterns in gene class utilization across development are very similar to patterns seen during the embryonic development of Drosophila, suggesting conservation of the temporal progression of these processes across 550 million years of evolution. The temporal profiles of gene expression and activation of processes revealed here provide intriguing insights into the mechanisms of mammalian development, embryogenesis, and organogenesis, as well as into the evolution of developmental processes.