Opportunities for basic, clinical, and bioethics research at the intersection of machine learning and genomics.

The data-intensive fields of genomics and machine learning (ML) are in an early stage of convergence. Genomics researchers increasingly seek to harness the power of ML methods to extract knowledge from their data; conversely, ML scientists recognize that genomics offers a wealth of large, complex, and well-annotated datasets that can be used as a substrate for developing biologically relevant algorithms and applications. The National Human Genome Research Institute (NHGRI) inquired with researchers working in these two fields to identify common challenges and receive recommendations to better support genomic research efforts using ML approaches. Those included increasing the amount and variety of training datasets by integrating genomic with multiomics, context-specific (e.g., by cell type), and social determinants of health datasets; reducing the inherent biases of training datasets; prioritizing transparency and interpretability of ML methods; and developing privacy-preserving technologies for research participants' data.

Conventional Co-Housing Modulates Murine Gut Microbiota and Hematopoietic Gene Expression.

Specific-pathogen-free (SPF) mice have improved hematopoietic characteristics relative to germ-free mice, however, it is not clear whether improvements in hematopoietic traits will continue when the level of microorganism exposure is further increased. We co-housed SPF C57BL/6 mice in a conventional facility (CVT) and found a significant increase in gut microbiota diversity along with increased levels of myeloid cells and T cells, especially effector memory T cells. Through single cell RNA sequencing of sorted KL (c-Kit+Lin-) cells, we imputed a decline in long-term hematopoietic stem cells and an increase in granulocyte-monocyte progenitors in CVT mice with up-regulation of genes associated with cell survival. Bone marrow transplantation through competitive repopulation revealed a significant increase in KSL (c-Kit+Sca-1+Lin-) cell reconstitution in recipients of CVT donor cells which occurred when donors were co-housed for both one and twelve months. However, there was minimal to no gain in mature blood cell engraftment in recipients of CVT donor cells relative to those receiving SPF donor cells. We conclude that co-housing SPF mice with mice born in a conventional facility increased gut microbiota diversity, augmented myeloid cell production and T cell activation, stimulated KSL cell reconstitution, and altered hematopoietic gene expression.

TREML4 Promotes Inflammatory Programs in Human and Murine Macrophages and Alters Atherosclerosis Lesion Composition in the Apolipoprotein E Deficient Mouse.

The Triggering Receptor Expressed on Myeloid cells-like 4 (TREML4) is a member of the TREM receptor family, known modulators of inflammatory responses. We have previously found that TREML4 expression positively correlates with human coronary arterial calcification (CAC). However, the role of TREML4 in the pathogenesis of cardiovascular disease remains incompletely defined. Since macrophages play a key role in inflammatory conditions, we investigated if activated macrophages selectively expressed TREML4 and found that carriage of either one of the eQTL SNP's previously associated with increased TREML4 expression conferred higher expression in human inflammatory macrophages (M1) compared to alternatively activated macrophages (M2). Furthermore, we found that TREML4 expression in human M1 dysregulated several inflammatory pathways related to leukocyte activation, apoptosis and extracellular matrix degradation. Similarly, murine M1 expressed substantial levels of Treml4, as did oxLDL treated macrophages. Transcriptome analysis confirmed that murine Treml4 controls the expression of genes related to inflammation and lipid regulation pathways, suggesting a possible role in atherosclerosis. Analysis of Apoe-/-/Treml4-/- mice showed reduced plaque burden and lesion complexity as indicated by decreased stage scores, macrophage content and collagen deposition. Finally, transcriptome analysis of oxLDL-loaded murine macrophages showed that Treml4 represses a specific set of genes related to carbohydrate, ion and amino acid membrane transport. Metabolomic analysis confirmed that Treml4 deficiency may promote a beneficial relationship between iron homeostasis and glucose metabolism. Together, our results suggest that Treml4 plays a role in the development of cardiovascular disease, as indicated by Treml4-dependent dysregulation of macrophage inflammatory pathways, macrophage metabolism and promotion of vulnerability features in advanced lesions.

Attenuation of immune-mediated bone marrow damage in conventionally housed mice.

In humans, bone marrow (BM) failure syndromes, both constitutional and acquired, predispose to myeloid malignancies. We have modeled acquired immune aplastic anemia, the paradigmatic disease of these syndromes, in the mouse by infusing lymph node cells from specific pathogen-free (SPF) CD45.1 congenic C57BL/6 (B6) donors into hybrid CByB6F1 recipients housed either in conventional (CVB) or SPF facilities. The severity of BM damage was reduced in CVB recipients; they also had reduced levels of CD44+ CD62L- effector memory T cells, reduced numbers of donor-type CD44+ T cells, and reduced expansion of donor-type CD8 T cells carrying T-cell receptor ß-variable regions 07, 11, and 17. Analyses of fecal samples through 16S ribosomal RNA amplicon sequencing revealed greater gut microbial alpha diversity in CVB mice relative to that of SPF mice. Thus, the presence of a broader spectrum of gut microorganisms in CVB-housed CByB6F1 could have primed recipient animal's immune system leading to suppression of allogeneic donor T-cell activation and expansion and attenuation of host BM destruction. These results suggest the potential benefit of diverse gut microbiota in patients receiving BM transplants.

Laboratory mice born to wild mice have natural microbiota and model human immune responses.

Laboratory mouse studies are paramount for understanding basic biological phenomena but also have limitations. These include conflicting results caused by divergent microbiota and limited translational research value. To address both shortcomings, we transferred C57BL/6 embryos into wild mice, creating "wildlings." These mice have a natural microbiota and pathogens at all body sites and the tractable genetics of C57BL/6 mice. The bacterial microbiome, mycobiome, and virome of wildlings affect the immune landscape of multiple organs. Their gut microbiota outcompete laboratory microbiota and demonstrate resilience to environmental challenges. Wildlings, but not conventional laboratory mice, phenocopied human immune responses in two preclinical studies. A combined natural microbiota- and pathogen-based model may enhance the reproducibility of biomedical studies and increase the bench-to-bedside safety and success of immunological studies.

NAA10 polyadenylation signal variants cause syndromic microphthalmia.

BACKGROUND: A single variant in NAA10 (c.471+2T>A), the gene encoding N-acetyltransferase 10, has been associated with Lenz microphthalmia syndrome. In this study, we aimed to identify causative variants in families with syndromic X-linked microphthalmia. METHODS: Three families, including 15 affected individuals with syndromic X-linked microphthalmia, underwent analyses including linkage analysis, exome sequencing and targeted gene sequencing. The consequences of two identified variants in NAA10 were evaluated using quantitative PCR and RNAseq. RESULTS: Genetic linkage analysis in family 1 supported a candidate region on Xq27-q28, which included NAA10. Exome sequencing identified a hemizygous NAA10 polyadenylation signal (PAS) variant, chrX:153,195,397T>C, c.*43A>G, which segregated with the disease. Targeted sequencing of affected males from families 2 and 3 identified distinct NAA10 PAS variants, chrX:g.153,195,401T>C, c.*39A>G and chrX:g.153,195,400T>C, c.*40A>G. All three variants were absent from gnomAD. Quantitative PCR and RNAseq showed reduced NAA10 mRNA levels and abnormal 3' UTRs in affected individuals. Targeted sequencing of NAA10 in 376 additional affected individuals failed to identify variants in the PAS. CONCLUSION: These data show that PAS variants are the most common variant type in NAA10-associated syndromic microphthalmia, suggesting reduced RNA is the molecular mechanism by which these alterations cause microphthalmia/anophthalmia. We reviewed recognised variants in PAS associated with Mendelian disorders and identified only 23 others, indicating that NAA10 harbours more than 10% of all known PAS variants. We hypothesise that PAS in other genes harbour unrecognised pathogenic variants associated with Mendelian disorders. The systematic interrogation of PAS could improve genetic testing yields.

Commensal-specific T cell plasticity promotes rapid tissue adaptation to injury.

Barrier tissues are primary targets of environmental stressors and are home to the largest number of antigen-experienced lymphocytes in the body, including commensal-specific T cells. We found that skin-resident commensal-specific T cells harbor a paradoxical program characterized by a type 17 program associated with a poised type 2 state. Thus, in the context of injury and exposure to inflammatory mediators such as interleukin-18, these cells rapidly release type 2 cytokines, thereby acquiring contextual functions. Such acquisition of a type 2 effector program promotes tissue repair. Aberrant type 2 responses can also be unleashed in the context of local defects in immunoregulation. Thus, commensal-specific T cells co-opt tissue residency and cell-intrinsic flexibility as a means to promote both local immunity and tissue adaptation to injury.

Gut microbiome-mediated bile acid metabolism regulates liver cancer via NKT cells.

Primary liver tumors and liver metastasis currently represent the leading cause of cancer-related death. Commensal bacteria are important regulators of antitumor immunity, and although the liver is exposed to gut bacteria, their role in antitumor surveillance of liver tumors is poorly understood. We found that altering commensal gut bacteria in mice induced a liver-selective antitumor effect, with an increase of hepatic CXCR6+ natural killer T (NKT) cells and heightened interferon-γ production upon antigen stimulation. In vivo functional studies showed that NKT cells mediated liver-selective tumor inhibition. NKT cell accumulation was regulated by CXCL16 expression of liver sinusoidal endothelial cells, which was controlled by gut microbiome-mediated primary-to-secondary bile acid conversion. Our study suggests a link between gut bacteria-controlled bile acid metabolism and liver antitumor immunosurveillance.

Non-classical Immunity Controls Microbiota Impact on Skin Immunity and Tissue Repair.

Mammalian barrier surfaces are constitutively colonized by numerous microorganisms. We explored how the microbiota was sensed by the immune system and the defining properties of such responses. Here, we show that a skin commensal can induce T cell responses in a manner that is restricted to non-classical MHC class I molecules. These responses are uncoupled from inflammation and highly distinct from pathogen-induced cells. Commensal-specific T cells express a defined gene signature that is characterized by expression of effector genes together with immunoregulatory and tissue-repair signatures. As such, non-classical MHCI-restricted commensal-specific immune responses not only promoted protection to pathogens, but also accelerated skin wound closure. Thus, the microbiota can induce a highly physiological and pleiotropic form of adaptive immunity that couples antimicrobial function with tissue repair. Our work also reveals that non-classical MHC class I molecules, an evolutionarily ancient arm of the immune system, can promote homeostatic immunity to the microbiota.

Genotype-driven identification of a molecular network predictive of advanced coronary calcium in ClinSeq® and Framingham Heart Study cohorts.

BACKGROUND: One goal of personalized medicine is leveraging the emerging tools of data science to guide medical decision-making. Achieving this using disparate data sources is most daunting for polygenic traits. To this end, we employed random forests (RFs) and neural networks (NNs) for predictive modeling of coronary artery calcium (CAC), which is an intermediate endo-phenotype of coronary artery disease (CAD). METHODS: Model inputs were derived from advanced cases in the ClinSeq®; discovery cohort (n=16) and the FHS replication cohort (n=36) from 89 th -99 th CAC score percentile range, and age-matched controls (ClinSeq®; n=16, FHS n=36) with no detectable CAC (all subjects were Caucasian males). These inputs included clinical variables and genotypes of 56 single nucleotide polymorphisms (SNPs) ranked highest in terms of their nominal correlation with the advanced CAC state in the discovery cohort. Predictive performance was assessed by computing the areas under receiver operating characteristic curves (ROC-AUC). RESULTS: RF models trained and tested with clinical variables generated ROC-AUC values of 0.69 and 0.61 in the discovery and replication cohorts, respectively. In contrast, in both cohorts, the set of SNPs derived from the discovery cohort were highly predictive (ROC-AUC ≥0.85) with no significant change in predictive performance upon integration of clinical and genotype variables. Using the 21 SNPs that produced optimal predictive performance in both cohorts, we developed NN models trained with ClinSeq®; data and tested with FHS data and obtained high predictive accuracy (ROC-AUC=0.80-0.85) with several topologies. Several CAD and "vascular aging" related biological processes were enriched in the network of genes constructed from the predictive SNPs. CONCLUSIONS: We identified a molecular network predictive of advanced coronary calcium using genotype data from ClinSeq®; and FHS cohorts. Our results illustrate that machine learning tools, which utilize complex interactions between disease predictors intrinsic to the pathogenesis of polygenic disorders, hold promise for deriving predictive disease models and networks.

Pendrin localizes to the adrenal medulla and modulates catecholamine release.

Pendrin (Slc26a4) is a Cl(-)/HCO3 (-) exchanger expressed in renal intercalated cells and mediates renal Cl(-) absorption. With pendrin gene ablation, blood pressure and vascular volume fall, which increases plasma renin concentration. However, serum aldosterone does not significantly increase in pendrin-null mice, suggesting that pendrin regulates adrenal zona glomerulosa aldosterone production. Therefore, we examined pendrin expression in the adrenal gland using PCR, immunoblots, and immunohistochemistry. Pendrin protein was detected in adrenal lysates from wild-type but not pendrin-null mice. However, immunohistochemistry and qPCR of microdissected adrenal zones showed that pendrin was expressed in the adrenal medulla, rather than in cortex. Within the adrenal medulla, pendrin localizes to both epinephrine- and norepinephrine-producing chromaffin cells. Therefore, we examined plasma catecholamine concentration and blood pressure in wild-type and pendrin-null mice under basal conditions and then after 5 and 20 min of immobilization stress. Under basal conditions, blood pressure was lower in the mutant than in the wild-type mice, although epinephrine and norepinephrine concentrations were similar. Catecholamine concentration and blood pressure increased markedly in both groups with stress. With 20 min of immobilization stress, epinephrine and norepinephrine concentrations increased more in pendrin-null than in wild-type mice, although stress produced a similar increase in blood pressure in both groups. We conclude that pendrin is expressed in the adrenal medulla, where it blunts stress-induced catecholamine release.

Integrative DNA, RNA, and protein evidence connects TREML4 to coronary artery calcification.

Coronary artery calcification (CAC) is a heritable and definitive morphologic marker of atherosclerosis that strongly predicts risk for future cardiovascular events. To search for genes involved in CAC, we used an integrative transcriptomic, genomic, and protein expression strategy by using next-generation DNA sequencing in the discovery phase with follow-up studies using traditional molecular biology and histopathology techniques. RNA sequencing of peripheral blood from a discovery set of CAC cases and controls was used to identify dysregulated genes, which were validated by ClinSeq and Framingham Heart Study data. Only a single gene, TREML4, was upregulated in CAC cases in both studies. Further examination showed that rs2803496 was a TREML4 cis-eQTL and that the minor allele at this locus conferred up to a 6.5-fold increased relative risk of CAC. We characterized human TREML4 and demonstrated by immunohistochemical techniques that it is localized in macrophages surrounding the necrotic core of coronary plaques complicated by calcification (but not in arteries with less advanced disease). Finally, we determined by von Kossa staining that TREML4 colocalizes with areas of microcalcification within coronary plaques. Overall, we present integrative RNA, DNA, and protein evidence implicating TREML4 in coronary artery calcification. Our findings connect multimodal genomics data with a commonly used clinical marker of cardiovascular disease.

Identification of candidate genes involved in coronary artery calcification by transcriptome sequencing of cell lines.

BACKGROUND: Massively-parallel cDNA sequencing (RNA-Seq) is a new technique that holds great promise for cardiovascular genomics. Here, we used RNA-Seq to study the transcriptomes of matched coronary artery disease cases and controls in the ClinSeq® study, using cell lines as tissue surrogates. RESULTS: Lymphoblastoid cell lines (LCLs) from 16 cases and controls representing phenotypic extremes for coronary calcification were cultured and analyzed using RNA-Seq. All cell lines were then independently re-cultured and along with another set of 16 independent cases and controls, were profiled with Affymetrix microarrays to perform a technical validation of the RNA-Seq results. Statistically significant changes (p < 0.05) were detected in 186 transcripts, many of which are expressed at extremely low levels (5-10 copies/cell), which we confirmed through a separate spike-in control RNA-Seq experiment. Next, by fitting a linear model to exon-level RNA-Seq read counts, we detected signals of alternative splicing in 18 transcripts. Finally, we used the RNA-Seq data to identify differential expression (p < 0.0001) in eight previously unannotated regions that may represent novel transcripts. Overall, differentially expressed genes showed strong enrichment (p = 0.0002) for prior association with cardiovascular disease. At the network level, we found evidence for perturbation in pathways involving both cardiovascular system development and function as well as lipid metabolism. CONCLUSIONS: We present a pilot study for transcriptome involvement in coronary artery calcification and demonstrate how RNA-Seq analyses using LCLs as a tissue surrogate may yield fruitful results in a clinical sequencing project. In addition to canonical gene expression, we present candidate variants from alternative splicing and novel transcript detection, which have been unexplored in the context of this disease.

Massively parallel sequencing, aCGH, and RNA-Seq technologies provide a comprehensive molecular diagnosis of Fanconi anemia.

Current methods for detecting mutations in Fanconi anemia (FA)-suspected patients are inefficient and often miss mutations. We have applied recent advances in DNA sequencing and genomic capture to the diagnosis of FA. Specifically, we used custom molecular inversion probes or TruSeq-enrichment oligos to capture and sequence FA and related genes, including introns, from 27 samples from the International Fanconi Anemia Registry at The Rockefeller University. DNA sequencing was complemented with custom array comparative genomic hybridization (aCGH) and RNA sequencing (RNA-seq) analysis. aCGH identified deletions/duplications in 4 different FA genes. RNA-seq analysis revealed lack of allele specific expression associated with a deletion and splicing defects caused by missense, synonymous, and deep-in-intron variants. The combination of TruSeq-targeted capture, aCGH, and RNA-seq enabled us to identify the complementation group and biallelic germline mutations in all 27 families: FANCA (7), FANCB (3), FANCC (3), FANCD1 (1), FANCD2 (3), FANCF (2), FANCG (2), FANCI (1), FANCJ (2), and FANCL (3). FANCC mutations are often the cause of FA in patients of Ashkenazi Jewish (AJ) ancestry, and we identified 2 novel FANCC mutations in 2 patients of AJ ancestry. We describe here a strategy for efficient molecular diagnosis of FA.

Gene expression signatures of coronary heart disease.

OBJECTIVE: To identify transcriptomic biomarkers of coronary heart disease (CHD) in 188 cases with CHD and 188 age- and sex-matched controls who were participants in the Framingham Heart Study. APPROACH AND RESULTS: A total of 35 genes were differentially expressed in cases with CHD versus controls at false discovery rate<0.5, including GZMB, TMEM56, and GUK1. Cluster analysis revealed 3 gene clusters associated with CHD, 2 linked to increased erythrocyte production and a third to reduced natural killer and T cell activity in cases with CHD. Exon-level results corroborated and extended the gene-level results. Alternative splicing analysis suggested that GUK1 and 38 other genes were differentially spliced in cases with CHD versus controls. Gene Ontology analysis linked ubiquitination and T-cell-related pathways with CHD. CONCLUSIONS: Two bioinformatically defined groups of genes show consistent associations with CHD. Our findings are consistent with the hypotheses that hematopoesis is upregulated in CHD, possibly reflecting a compensatory mechanism, and that innate immune activity is disrupted in CHD or altered by its treatment. Transcriptomic signatures may be useful in identifying pathways associated with CHD and point toward novel therapeutic targets for its treatment and prevention.

Synthetic spike-in standards improve run-specific systematic error analysis for DNA and RNA sequencing.

While the importance of random sequencing errors decreases at higher DNA or RNA sequencing depths, systematic sequencing errors (SSEs) dominate at high sequencing depths and can be difficult to distinguish from biological variants. These SSEs can cause base quality scores to underestimate the probability of error at certain genomic positions, resulting in false positive variant calls, particularly in mixtures such as samples with RNA editing, tumors, circulating tumor cells, bacteria, mitochondrial heteroplasmy, or pooled DNA. Most algorithms proposed for correction of SSEs require a data set used to calculate association of SSEs with various features in the reads and sequence context. This data set is typically either from a part of the data set being "recalibrated" (Genome Analysis ToolKit, or GATK) or from a separate data set with special characteristics (SysCall). Here, we combine the advantages of these approaches by adding synthetic RNA spike-in standards to human RNA, and use GATK to recalibrate base quality scores with reads mapped to the spike-in standards. Compared to conventional GATK recalibration that uses reads mapped to the genome, spike-ins improve the accuracy of Illumina base quality scores by a mean of 5 Phred-scaled quality score units, and by as much as 13 units at CpG sites. In addition, since the spike-in data used for recalibration are independent of the genome being sequenced, our method allows run-specific recalibration even for the many species without a comprehensive and accurate SNP database. We also use GATK with the spike-in standards to demonstrate that the Illumina RNA sequencing runs overestimate quality scores for AC, CC, GC, GG, and TC dinucleotides, while SOLiD has less dinucleotide SSEs but more SSEs for certain cycles. We conclude that using these DNA and RNA spike-in standards with GATK improves base quality score recalibration.

Computational methods for the analysis of primate mobile elements.

Transposable elements (TE), defined as discrete pieces of DNA that can move from one site to another site in genomes, represent significant components of eukaryotic genomes, including primates. Comparative genome-wide analyses have revealed the considerable structural and functional impact of TE families on primate genomes. Insights into these questions have come in part from the development of computational methods that allow detailed and reliable identification, annotation, and evolutionary analyses of the many TE families that populate primate genomes. Here, we present an overview of these computational methods and describe efficient data mining strategies for providing a comprehensive picture of TE biology in newly available genome sequences.

Internal priming: an opportunistic pathway for L1 and Alu retrotransposition in hominins.

Retrotransposons, specifically Alu and L1 elements, have been especially successful in their expansion throughout primate genomes. While most of these elements integrate through an endonuclease-mediated process termed target primed reverse transcription, a minority integrate using alternative methods. Here we present evidence for one such mechanism, which we term internal priming and demonstrate that loci integrating through this mechanism are qualitatively different from "classical" insertions. Previous examples of this mechanism are limited to cell culture assays, which show that reverse transcription can initiate upstream of the 3' poly-A tail during retrotransposon integration. To detect whether this mechanism occurs in vivo as well as in cell culture, we have analyzed the human genome for internal priming events using recently integrated L1 and Alu elements. Our examination of the human genome resulted in the recovery of twenty events involving internal priming insertions, which are structurally distinct from both classical TPRT-mediated insertions and non-classical insertions. We suggest two possible mechanisms by which these internal priming loci are created and provide evidence supporting a role in staggered DNA double-strand break repair. Also, we demonstrate that the internal priming process is associated with inter-chromosomal duplications and the insertion of filler DNA.

Repair-mediated duplication by capture of proximal chromosomal DNA has shaped vertebrate genome evolution.

DNA double-strand breaks (DSBs) are a common form of cellular damage that can lead to cell death if not repaired promptly. Experimental systems have shown that DSB repair in eukaryotic cells is often imperfect and may result in the insertion of extra chromosomal DNA or the duplication of existing DNA at the breakpoint. These events are thought to be a source of genomic instability and human diseases, but it is unclear whether they have contributed significantly to genome evolution. Here we developed an innovative computational pipeline that takes advantage of the repetitive structure of genomes to detect repair-mediated duplication events (RDs) that occurred in the germline and created insertions of at least 50 bp of genomic DNA. Using this pipeline we identified over 1,000 probable RDs in the human genome. Of these, 824 were intra-chromosomal, closely linked duplications of up to 619 bp bearing the hallmarks of the synthesis-dependent strand-annealing repair pathway. This mechanism has duplicated hundreds of sequences predicted to be functional in the human genome, including exons, UTRs, intron splice sites and transcription factor binding sites. Dating of the duplication events using comparative genomics and experimental validation revealed that the mechanism has operated continuously but with decreasing intensity throughout primate evolution. The mechanism has produced species-specific duplications in all primate species surveyed and is contributing to genomic variation among humans. Finally, we show that RDs have also occurred, albeit at a lower frequency, in non-primate mammals and other vertebrates, indicating that this mechanism has been an important force shaping vertebrate genome evolution.

Mobile elements create structural variation: analysis of a complete human genome.

Structural variants (SVs) are common in the human genome. Because approximately half of the human genome consists of repetitive, transposable DNA sequences, it is plausible that these elements play an important role in generating SVs in humans. Sequencing of the diploid genome of one individual human (HuRef) affords us the opportunity to assess, for the first time, the impact of mobile elements on SVs in an individual in a thorough and unbiased fashion. In this study, we systematically evaluated more than 8000 SVs to identify mobile element-associated SVs as small as 100 bp and specific to the HuRef genome. Combining computational and experimental analyses, we identified and validated 706 mobile element insertion events (including Alu, L1, SVA elements, and nonclassical insertions), which added more than 305 kb of new DNA sequence to the HuRef genome compared with the Human Genome Project (HGP) reference sequence (hg18). We also identified 140 mobile element-associated deletions, which removed approximately 126 kb of sequence from the HuRef genome. Overall, approximately 10% of the HuRef-specific indels larger than 100 bp are caused by mobile element-associated events. More than one-third of the insertion/deletion events occurred in genic regions, and new Alu insertions occurred in exons of three human genes. Based on the number of insertions and the estimated time to the most recent common ancestor of HuRef and the HGP reference genome, we estimated the Alu, L1, and SVA retrotransposition rates to be one in 21 births, 212 births, and 916 births, respectively. This study presents the first comprehensive analysis of mobile element-related structural variants in the complete DNA sequence of an individual and demonstrates that mobile elements play an important role in generating inter-individual structural variation.