N. Alice Yamada

Mapping and sequencing of structural variation from eight human genomes

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale—particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation—a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

The Fine-Scale and Complex Architecture of Human Copy-Number Variation

Despite considerable excitement over the potential functional significance of copy-number variants (CNVs), we still lack knowledge of the fine-scale architecture of the large majority of CNV regions in the human genome. In this study, we used a high-resolution array-based comparative genomic hybridization (aCGH) platform that targeted known CNV regions of the human genome at approximately 1 kb resolution to interrogate the genomic DNAs of 30 individuals from four HapMap populations. Our results revealed that 1020 of 1153 CNV loci (88%) were actually smaller in size than what is recorded in the Database of Genomic Variants based on previously published studies. A reduction in size of more than 50% was observed for 876 CNV regions (76%). We conclude that the total genomic content of currently known common human CNVs is likely smaller than previously thought. In addition, approximately 8% of the CNV regions observed in multiple individuals exhibited genomic architectural complexity in the form of smaller CNVs within larger ones and CNVs with interindividual variation in breakpoints. Future association studies that aim to capture the potential influences of CNVs on disease phenotypes will need to consider how to best ascertain this previously uncharacterized complexity.

Characterization of missing human genome sequences and copy-number polymorphic insertions

The extent of human genomic structural variation suggests that there must be portions of the genome yet to be discovered, annotated and characterized at the sequence level. We present a resource and analysis of 2,363 new insertion sequences corresponding to 720 genomic loci. We found that a substantial fraction of these sequences are either missing, fragmented or misassigned when compared to recent de novo sequence assemblies from short-read next-generation sequence data. We determined that 18–37% of these new insertions are copy-number polymorphic, including loci that show extensive population stratification among Europeans, Asians and Africans. Complete sequencing of 156 of these insertions identified new exons and conserved noncoding sequences not yet represented in the reference genome. We developed a method to accurately genotype these new insertions by mapping next-generation sequencing datasets to the breakpoint, thereby providing a means to characterize copy-number status for regions previously inaccessible to single-nucleotide polymorphism microarrays.

A Prader–Willi locus lncRNA cloud modulates diurnal genes and energy expenditure

Prader–Willi syndrome (PWS), a genetic disorder of obesity, intellectual disability and sleep abnormalities, is caused by loss of non-coding RNAs on paternal chromosome 15q11-q13. The imprinted minimal PWS locus encompasses a long non-coding RNA (lncRNA) transcript processed into multiple SNORD116 small nucleolar RNAs and the spliced exons of the host gene, 116HG. However, both the molecular function and the disease relevance of the spliced lncRNA 116HG are unknown. Here, we show that 116HG forms a subnuclear RNA cloud that co-purifies with the transcriptional activator RBBP5 and active metabolic genes, remains tethered to the site of its transcription and increases in size in post-natal neurons and during sleep. Snord116del mice lacking 116HG exhibited increased energy expenditure corresponding to the dysregulation of diurnally expressed Mtor and circadian genes Clock, Cry1 and Per2. These combined genomic and metabolic analyses demonstrate that 116HG regulates the diurnal energy expenditure of the brain. These novel molecular insights into the energy imbalance in PWS should lead to improved therapies and understanding of lncRNA roles in complex neurodevelopmental and metabolic disorders.

Chromosome Conformation Paints Reveal The Role Of Lamina Association In Genome Organization And Regulation

Non-random, dynamic three-dimensional organization of the nucleus is important for regulation of gene expression. Numerous studies using chromosome conformation capture strategies have uncovered ensemble organizational principles of individual chromosomes, including organization into active (A) and inactive (B) compartments. In addition, large inactive regions of the genome appear to be associated with the nuclear lamina, the so-called Lamina Associated Domains (LADs). However, the interrelationship between overall chromosome conformation and association of domains with the nuclear lamina remains unclear. In particular, the 3D organization of LADs within the context of the entire chromosome has not been investigated. In this study, we describe “chromosome conformation paints” to determine the relationship in situ between LAD and non-LAD regions of the genome in single cells. We find that LADs organize into constrained and compact regions at the nuclear lamina, and these findings are supported by an integrated analysis of both DamID and Hi-C data. Using a refined algorithm to identify active (A) and inactive (B) compartments from Hi-C data, we demonstrate that the LADs correspond to the B compartment. We demonstrate that in situ single cell chromosome organization is strikingly predicted by integrating both Hi-C and DamID data into a chromosome conformation model. In addition, using the chromosome conformation paints, we demonstrate that LAD (and B-compartment) organization is dependent upon both chromatin state and Lamin A/C. Finally, we demonstrate that small regions within LADs escape the repressive regime at the peripheral zone to interact with the A-compartment and are enriched for both transcription start sites (TSSs) and active enhancers.