Ryan E. Mills

Mapping copy number variation by population-scale genome sequencing

Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.

An integrated map of structural variation in 2,504 human genomes

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

Natural Mutagenesis of Human Genomes by Endogenous Retrotransposons

Two abundant classes of mobile elements, namely Alu and L1 elements, continue to generate new retrotransposon insertions in human genomes. Estimates suggest that these elements have generated millions of new germline insertions in individual human genomes worldwide. Unfortunately, current technologies are not capable of detecting most of these young insertions, and the true extent of germline mutagenesis by endogenous human retrotransposons has been difficult to examine. Here, we describe technologies for detecting these young retrotransposon insertions and demonstrate that such insertions indeed are abundant in human populations. We also found that new somatic L1 insertions occur at high frequencies in human lung cancer genomes. Genome-wide analysis suggests that altered DNA methylation may be responsible for the high levels of L1 mobilization observed in these tumors. Our data indicate that transposon-mediated mutagenesis is extensive in human genomes and is likely to have a major impact on human biology and diseases.

Comprehensive assessment of array-based platforms and calling algorithms for detection of copy number variants.

We have systematically compared copy number variant (CNV) detection on eleven microarrays to evaluate data quality and CNV calling, reproducibility, concordance across array platforms and laboratory sites, breakpoint accuracy and analysis tool variability. Different analytic tools applied to the same raw data typically yield CNV calls with <50% concordance. Moreover, reproducibility in replicate experiments is <70% for most platforms. Nevertheless, these findings should not preclude detection of large CNVs for clinical diagnostic purposes because large CNVs with poor reproducibility are found primarily in complex genomic regions and would typically be removed by standard clinical data curation. The striking differences between CNV calls from different platforms and analytic tools highlight the importance of careful assessment of experimental design in discovery and association studies and of strict data curation and filtering in diagnostics. The CNV resource presented here allows independent data evaluation and provides a means to benchmark new algorithms.

Discovery of common Asian copy number variants using integrated high-resolution array CGH and massively parallel DNA sequencing

Copy number variants (CNVs) account for the majority of human genomic diversity in terms of base coverage. Here, we have developed and applied a new method to combine high-resolution array comparative genomic hybridization (CGH) data with whole-genome DNA sequencing data to obtain a comprehensive catalog of common CNVs in Asian individuals. The genomes of 30 individuals from three Asian populations (Korean, Chinese and Japanese) were interrogated with an ultra-high-resolution array CGH platform containing 24 million probes. Whole-genome sequencing data from a reference genome (NA10851, with 28.3× coverage) and two Asian genomes (AK1, with 27.8× coverage and AK2, with 32.0× coverage) were used to transform the relative copy number information obtained from array CGH experiments into absolute copy number values. We discovered 5,177 CNVs, of which 3,547 were putative Asian-specific CNVs. These common CNVs in Asian populations will be a useful resource for subsequent genetic studies in these populations, and the new method of calling absolute CNVs will be essential for applying CNV data to personalized medicine.

Natural genetic variation caused by small insertions and deletions in the human genome

Human genetic variation is expected to play a central role in personalized medicine. Yet only a fraction of the natural genetic variation that is harbored by humans has been discovered to date. Here we report almost 2 million small insertions and deletions (INDELs) that range from 1 bp to 10,000 bp in length in the genomes of 79 diverse humans. These variants include 819,363 small INDELs that map to human genes. Small INDELs frequently were found in the coding exons of these genes, and several lines of evidence indicate that such variation is a major determinant of human biological diversity. Microarray-based genotyping experiments revealed several interesting observations regarding the population genetics of small INDEL variation. For example, we found that many of our INDELs had high levels of linkage disequilibrium (LD) with both HapMap SNPs and with high-scoring SNPs from genome-wide association studies. Overall, our study indicates that small INDEL variation is likely to be a key factor underlying inherited traits and diseases in humans.

Recently Mobilized Transposons in the Human and Chimpanzee Genomes

Transposable genetic elements are abundant in the genomes of most organisms, including humans. These endogenous mutagens can alter genes, promote genomic rearrangements, and may help to drive the speciation of organisms. In this study, we identified almost 11,000 transposon copies that are differentially present in the human and chimpanzee genomes. Most of these transposon copies were mobilized after the existence of a common ancestor of humans and chimpanzees, approximately 6 million years ago. Alu, L1, and SVA insertions accounted for >95% of the insertions in both species. Our data indicate that humans have supported higher levels of transposition than have chimpanzees during the past several million years and have amplified different transposon subfamilies. In both species, approximately 34% of the insertions were located within known genes. These insertions represent a form of species-specific genetic variation that may have contributed to the differential evolution of humans and chimpanzees. In addition to providing an initial overview of recently mobilized elements, our collections will be useful for assessing the impact of these insertions on their hosts and for studying the transposition mechanisms of these elements.

Extensive Genetic Diversity and Substructuring Among Zebrafish Strains Revealed through Copy Number Variant Analysis

Copy number variants (CNVs) represent a substantial source of genomic variation in vertebrates and have been associated with numerous human diseases. Despite this, the extent of CNVs in the zebrafish, an important model for human disease, remains unknown. Using 80 zebrafish genomes, representing three commonly used laboratory strains and one native population, we constructed a genome-wide, high-resolution CNV map for the zebrafish comprising 6,080 CNV elements and encompassing 14.6% of the zebrafish reference genome. This amount of copy number variation is four times that previously observed in other vertebrates, including humans. Moreover, 69% of the CNV elements exhibited strain specificity, with the highest number observed for Tubingen. This variation likely arose, in part, from Tubingens large founding size and composite population origin. Additional population genetic studies also provided important insight into the origins and substructure of these commonly used laboratory strains. This extensive variation among and within zebrafish strains may have functional effects that impact phenotype and, if not properly addressed, such extensive levels of germ-line variation and population substructure in this commonly used model organism can potentially confound studies intended for translation to human diseases.

Intersection of diverse neuronal genomes and neuropsychiatric disease: The Brain Somatic Mosaicism Network

Single-cell diversity in the brain The cells that make up an organism may all start from one genome, but somatic mutations mean that somewhere along the line of development, an organisms individual cellular genomes diverge. McConnell et al. review the implications and causes of single-cell genomic diversity for brain function. Somatic mutations caused by mobile genetic elements or errors in DNA repair may underlie certain neuropsychiatric disorders. Science, this issue p. eaal1641 BACKGROUND Elucidating the genetic architecture of neuropsychiatric disorders remains a major scientific and medical challenge. Emerging genomic technologies now permit the analysis of somatic mosaicism in human tissues. The measured frequencies of single-nucleotide variants (SNVs), small insertion/deletion (indel) mutations, structural variants [including copy number variants (CNVs), inversions, translocations, and whole-chromosome gains or losses], and mobile genetic element insertions (MEIs) indicate that each neuron may harbor hundreds of somatic mutations. Given the long life span of neurons and their central role in neural circuits and behavior, somatic mosaicism represents a potential mechanism that may contribute to neuronal diversity and the etiology of numerous neuropsychiatric disorders. ADVANCES Somatic mutations that confer cellular proliferative or cellular survival phenotypes have been identified in patients with cortical malformations. These data have led to the hypothesis that somatic mutations may also confer phenotypes to subsets of neurons, which could increase the risk of developing certain neuropsychiatric disorders. Genomic technologies, including advances in long-read, next-generation DNA sequencing technologies, single-cell genomics, and cutting-edge bioinformatics, can now make it possible to determine the types and frequencies of somatic mutations within the human brain. However, a comprehensive understanding of the contribution of somatic mosaicism to neurotypical brain development and neuropsychiatric disease requires a coordinated, multi-institutional effort. The National Institute of Mental Health (NIMH) has formed a network of 18 investigative teams representing 15 institutions called the Brain Somatic Mosaicism Network (BSMN). Each research team will use an array of genomic technologies to exploit well-curated human tissue repositories in an effort to define the frequency and pattern of somatic mutations in neurotypical individuals and in schizophrenia, autism spectrum disorder, bipolar disorder, Tourette syndrome, and epilepsy patient populations. Collectively, these efforts are estimated to generate a community resource of more than 10,000 DNA-sequencing data sets and will enable a cross-platform integrated analysis with other NIMH initiatives, such as the PsychENCODE project and the CommonMind Consortium. OUTLOOK A fundamental open question in neurodevelopmental genetics is whether and how somatic mosaicism may contribute to neuronal diversity within the neurotypical spectrum and in diseased brains. Healthy individuals may harbor known pathogenic somatic mutations at subclinical frequencies, and the local composition of neural cell types may be altered by mutations conferring prosurvival phenotypes in subsets of neurons. By extension, the neurotypical architecture of somatic mutations may confer circuit-level differences that would not be present if every neuron had an identical genome. Given the apparent abundance of somatic mutations within neurons, an in-depth understanding of how different types of somatic mosaicism affect neural function could yield mechanistic insight into the etiology of neurodevelopmental and neuropsychiatric disorders. The BSMN will examine large collections of postmortem brain tissue from neurotypical individuals and patients with neuropsychiatric disorders. By sequencing brain DNA and single neuronal genomes directly, rather than genomic DNA derived from peripheral blood or other somatic tissues, the BSMN will test the hypothesis that brain somatic variants contribute to neuropsychiatric disease. Notably, it is also possible that some inherited germline variants confer susceptibility to disease, which is later exacerbated by somatic mutations. Confirming such a scenario could increase our understanding of the genetic risk architecture of neuropsychiatric disease and may, in part, explain discordant neuropsychiatric phenotypes between identical twins. Results from these studies may lead to the discovery of biomarkers and genetic targets to improve the treatment of neuropsychiatric disease and may offer hope for improving the lives of patients and their families. Collectively, somatic SNVs, indels, structural variants (e.g., CNVs), and MEIs (e.g., L1 retrotransposition events) shape the genomic landscape of individual neurons. The Brain Somatic Mosaicism Network aims to systematically generate pioneering data on the types and frequencies of brain somatic mutations in both neurotypical individuals and those with neuropsychiatric disease. The resulting data will be shared as a large community resource. Neuropsychiatric disorders have a complex genetic architecture. Human genetic population-based studies have identified numerous heritable sequence and structural genomic variants associated with susceptibility to neuropsychiatric disease. However, these germline variants do not fully account for disease risk. During brain development, progenitor cells undergo billions of cell divisions to generate the ~80 billion neurons in the brain. The failure to accurately repair DNA damage arising during replication, transcription, and cellular metabolism amid this dramatic cellular expansion can lead to somatic mutations. Somatic mutations that alter subsets of neuronal transcriptomes and proteomes can, in turn, affect cell proliferation and survival and lead to neurodevelopmental disorders. The long life span of individual neurons and the direct relationship between neural circuits and behavior suggest that somatic mutations in small populations of neurons can significantly affect individual neurodevelopment. The Brain Somatic Mosaicism Network has been founded to study somatic mosaicism both in neurotypical human brains and in the context of complex neuropsychiatric disorders.

The genomic landscape of polymorphic human nuclear mitochondrial insertions

The transfer of mitochondrial genetic material into the nuclear genomes of eukaryotes is a well-established phenomenon that has been previously limited to the study of static reference genomes. The recent advancement of high throughput sequencing has enabled an expanded exploration into the diversity of polymorphic nuclear mitochondrial insertions (NumtS) within human populations. We have developed an approach to discover and genotype novel Numt insertions using whole genome, paired-end sequencing data. We have applied this method to a thousand individuals in 20 populations from the 1000 Genomes Project and other datasets and identified 141 new sites of Numt insertions, extending our current knowledge of existing NumtS by almost 20%. We find that recent Numt insertions are derived from throughout the mitochondrial genome, including the D-loop, and have integration biases that differ in some respects from previous studies on older, fixed NumtS in the reference genome. We determined the complete inserted sequence for a subset of these events and have identified a number of nearly full-length mitochondrial genome insertions into nuclear chromosomes. We further define their age and origin of insertion and present an analysis of their potential impact to ongoing studies of mitochondrial heteroplasmy and disease.