¡Descarga Family Whole Genome Sequencing: Precise Inheritance Analysis & Rare Variant Identification y más Guías, Proyectos, Investigaciones en PDF de Biología solo en Docsity!
Analysis of Genetic Inheritance in a Family Quartet by Whole
Genome Sequencing
Jared C. Roach 1,, Gustavo Glusman 1,, Arian F.A. Smit 1,, Chad D. Huff 1,2,, Robert Hubley^1 , Paul T. Shannon^1 , Lee Rowen^1 , Krishna P. Pant^3 , Nathan Goodman^1 , Michael Bamshad^4 , Jay Shendure^5 , Radoje Drmanac^3 , Lynn B. Jorde^2 , Leroy Hood 1,+, and David J. Galas 1,+ (^1) Institute for Systems Biology, Seattle, WA 98103 (^2) Department of Human Genetics, Eccles Institute of Human Genetics, University of Utah, Salt Lake City, Utah 84112, USA (^3) Complete Genomics, Inc. (CGI), Mountain View, CA 94043 (^4) Department of Pediatrics, University of Washington, Seattle WA 98195 (^5) Department of Genome Sciences, University of Washington, Seattle WA 98195
Abstract
We analyzed the whole genome sequences of a family of four, consisting of two siblings and their parents. Family-based sequencing allowed us to delineate recombination sites precisely, identify 70% of the sequencing errors, and identify very rare SNVs. We also directly estimated a human intergeneration mutation rate of ∼1.1×10-8^ per position per haploid genome. Both offspring in this family have two recessive disorders--Miller syndrome, for which the gene was concurrently identified, and primary ciliary dyskinesia, for which causative genes have been previously identified. Family-based genome analysis enabled us to narrow the candidate genes for both of these Mendelian disorders to only four. Our results demonstrate the unique value of complete genome sequencing in families.
Keywords whole genome sequencing; rare genetic disease; inheritance analysis; recessive models; de novo mutations; recombination hotspot; crossover; haploidentity; haploidentical block; inheritance state; inheritance vector; HMM; haplotype; Miller syndrome; POADS; DHODH; DNAH5; KIAA0556; CES
Whole-genome sequences from four members of a family represent a qualitatively different type of genetic data than whole-genome sequences from individual or sets of unrelated genomes. They enable inheritance analyses that detect errors and permit the identification of precise locations of recombination events. This leads in turn to near-complete knowledge of inheritance states by precisely determining the parental chromosomal origins of sequence blocks in offspring. Confident predictions of inheritance states and haplotypes power analyses including identification of genomic features with non-classical inheritance patterns such as hemizygous deletions or copy number variants (CNVs). Identification of inheritance patterns in the pedigree permits the detection of ∼70% of sequencing errors and sharply
- (^) Corresponding authors: dgalas@systemsbiology.org; lhood@systemsbiology.org. *These authors contributed equally.
NIH Public Access
Author Manuscript
Science. Author manuscript; available in PMC 2011 February 10.
Published in final edited form as:
Science. 2010 April 30; 328(5978): 636–639. doi:10.1126/science.1186802.
NIH-PA Author Manuscript
NIH-PA Author Manuscript
NIH-PA Author Manuscript
reduces the search space for disease-causing variants. These analyses would be far less powerful in studies that had fewer markers (e.g., standard genotype or exome datasets) or that had sequence from fewer family members.
DNA from each family member was extracted from peripheral blood cells and sequenced by Complete Genomics Inc. with a nanoarray-based short-read sequencing-by-ligation technology (1) including an adaptation of the pairwise end-sequencing strategy (2). Reads were mapped to the NCBI reference genome (coverage averaged 40×; Fig. S1, Tables S1 & S2). Polymorphic markers employed for this analysis were SNVs with at least two variants among the four genotypes of the family, averaging 802 bp between markers. We observed 4,471,510 positions at which at least one family member had an allele that varied from the reference genome. This corresponds to a Watterson's theta ( θW ) of 9.5×10-4^ per site for the two parents and the reference sequence (3), given the fraction of the genome successfully genotyped in each parent (Fig. S1). This is a close match to the estimate of θW = 9.3×10- we obtained by combining two previously published European genomes and the reference sequence (7). Of the 4.5 million variant positions, 3,665,772 were variable within the family
- the rest were homozygous and identical in all four members. Comparisons to known SNVs show that 323,255 of these 3.7 million SNVs are novel.
For each meiosis in a pedigree, each base position in a resulting gamete will have inherited one of two parental alleles. The number of inheritance patterns of the segregation of alleles in gametes is therefore 2 n , where n is the number of meioses in a pedigree. In a nuclear family of four, the Mendelian inheritance patterns can be grouped into four inheritance states for each variant position, with children receiving: 1) the same allele from both the mother and the father (identical), 2) the same allele from the mother but opposites from the father (haploidentical maternal), 3) the same allele from the father, but opposites from the mother (haploidentical paternal), and 4) opposites from both parents (nonidentical) (Fig. S2). Adjacent variant base pairs in alignments of the family genomes have the same inheritance state unless a recombination has occurred between these bases in one of the meioses. This delineates inheritance blocks.
Many algorithms can identify the boundaries of blocks, and theory-driven implementations are in wide use (4-6). For our complete genome sequence data, we developed an algorithm to identify all states, including non-Mendelian states. One non-Mendelian state will occur in regions where highly similar sequences are inadvertently compressed computationally (e.g., during sequence assembly of CNVs). In such a “compression block,” many positions will appear to be heterozygous in all individuals, regardless of the inheritance patterns of the positions contributing to the compression. Other non-Mendelian patterns are seen in regions prone to errors in sequence calling or assembly, or that have inherited hemizygous deletions. For both of these patterns, many positions will be observed as Mendelian inheritance errors (MIEs). Our algorithm identified six states: one for each of the four Mendelian inheritance states, one for a compression state, and one for an MIE-prone state (7). We identified 1.5% of the genome in this pedigree as 409 compression blocks and 1.7% as 126 error-prone blocks. Since these blocks are a source of false positives for recombination predictions, SNVs, and disease candidate alleles, it is important to identify them (Fig. 1). The power to precisely determine inheritance state boundaries is striking in families of at least four, and would be reduced had we sequenced fewer individuals (Fig. 2). Meiotic gene conversions could in principle be recognized in the same way as inheritance blocks: they would be indistinguishable from a short region flanked by meiotic recombinations in the same meiosis. We found that the great majority of candidate gene conversions regions were caused by reads mismapped to repetitive DNA, like CNVs or satellites, and did not conclusively identify gene conversion regions.
NIH-PA Author Manuscript
NIH-PA Author Manuscript
NIH-PA Author Manuscript
Our estimate is within one standard deviation (SD) of an earlier estimate of 1.7×10-8^ (SD: 9×10-9) based on 20 disease-causing loci (20). The rate we report is for autosomes, and should be lower than that of the Y chromosome, as in the male germline more cell divisions occur per generation. Though our rate differs approximately as expected from the recently reported estimate of 3.0×10-8^ (95% CI: 8.9×10-9^ – 7.0×10-8) for the Y chromosome, the error rates make this difference not significant (21).
Genomic inheritance analysis facilitates the identification of alleles that cause genetic disorders. Because genome sequences from a family of four provide near-exact determination of inheritance state boundaries, the number of false-positive disease gene candidates is greatly reduced compared to analyses lacking the context of a pedigree or complete genome sequence (Tables S3 & S4, Fig. 3). Two disorders in this family, Miller syndrome and primary ciliary dyskinesia, which affect both offspring but neither parent, provided an opportunity to test this application. A parsimonious explanation is that each phenotype arises from defects in a single gene, or a site regulating a single gene. The inheritance mode is undetermined, but a recessive mode is more consistent with observed data. We therefore examined each candidate variant by testing each of three inheritance modes: dominant, simple recessive, or compound heterozygote (a subcategory of recessive).
The two recessive modes require that both offspring have identical, dysfunctional variants for which the parents are heterozygous, and which may come either from the same position (simple recessive), or occur at distinct positions within the same gene (compound heterozygote). Genes consistent with these two recessive modes must lie in “identical” inheritance blocks, since both offspring are affected, limiting the search space to the 22% of the genome in these blocks. Since the phenotypes are rare, they are likely to be encoded by rare variants, which further limits the possibilities. Only two missense SNVs in the CES gene matched the simple recessive mode (Table S4), while three genes fit the compound heterozygote mode: DHODH , DNAH5 , and KIAA0556 (Fig. 1). A small number of possibly detrimental variants outside exons also matched the simple recessive mode: two in highly conserved regions, one in an intronic sequence near a splice site, five in non-protein coding transcripts, and one in a UTR. Concurrent with this study, the core exomes of the two affected offspring were sequenced along with those of two unrelated individuals with Miller syndrome (22). Compared to that study of only affected individuals, our analysis of just two affected offspring and their unaffected parents reduced the number of gene candidates in the core exome from nine to four; had we not sequenced the parents, we would have had 34 rather than four candidates (Fig. 3 & Table S5). The exome study supported DHODH as the primary gene for Miller syndrome. DNAH5 had been previously identified as a cause of primary ciliary dyskinesia, and so was likely the cause in these offspring as well (23).
Family genome analysis can clearly be effective for finding candidate genes encoding Mendelian traits because sequence accuracy is enhanced. In addition, delineation of recombination sites identifies inherited chromosome segments precisely, and reduces the chromosomal search space for candidate genes (in this case to 22% of the genome). The ability to identify large effects of very rare alleles in small pedigrees can complement the power of genome-wide association studies in identifying weak effects of common alleles in large populations. An unknown fraction of important phenotypes in humans are encoded by non-exonic variants identified only by whole-genome sequencing. When the cost of recruiting additional families is expensive relative to sequencing costs, sequencing genomes of families will be an economical strategy for the identification of many disease-causing genes. Constraining searches to very rare variants can provide considerable power, as recently demonstrated for Freeman-Sheldon syndrome and congenital chloride diarrhea (24,25). D e novo mutations can be assayed, either as we have reported here, or through family sequencing of more than two generations. As our knowledge of gene function
NIH-PA Author Manuscript
NIH-PA Author Manuscript
NIH-PA Author Manuscript
increases, we will be able to use the power of family genome analysis rapidly to identify disease-gene candidates. These data, along with relevant environmental and medical information, will characterize the integrated medical records of the future.
Supplementary Material
Refer to Web version on PubMed Central for supplementary material.
Acknowledgments
This study was supported by the University of Luxembourg Institute for Systems Biology Program and by these NIH grants: Center for Systems Biology GM076547 (L.H. and L.R.), RO1GM081083 (A.F.S. and G.G.), R01HL094976 and R21HG004749 (J.S.), RC2HG005608 (M.B. and J.S.), and R01HD048895 (M.B.). H. Tabor assisted with ethical review. J. Xing performed the PCA analysis. H. Mefford performed CNV analysis. A. Bigham and K. Buckingham evaluated candidate genes in unrelated individuals. D. Ballinger, A. Sparks, A. Halpern, and G. Nilsen assisted with sequencing and analysis. R. Bressler, S. Dee, and D. Mauldin assisted with bioinformatics. S. Ng and R. Qiu performed the capture array. S. Bloom obtained the resequencing data on the Illumina Genome Analyzer. M. Janer and S. Li performed Sequenom analysis. D. Cox commented on an early version of the manuscript. R. Durbin and D. Altshuler granted permission for our use of 1000genomes SNP data. The dbGAP accessions can be found at www.ncbi.nlm.nih.gov/projects/gap/cgibin/study.cgi?study_id=phs000244.v1.p (accession phs000244.v1.p1).
References and Notes
- Drmanac R, et al. Science November 5;2009 :1181498.
- Roach JC, Boysen C, Wang K, Hood L. Genomics Mar 20;1995 26:345. [PubMed: 7601461]
- Watterson GA. Theor Popul Biol Apr;1975 7:256. [PubMed: 1145509]
- Donnelly KP. Theor Popul Biol Feb;1983 23:34. [PubMed: 6857549]
- Kruglyak L, Daly MJ, Reeve-Daly MP, Lander ES. Am J Hum Genet Jun;1996 58:1347. [PubMed: 8651312]
- Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Nat Genet Jan;2002 30:97. [PubMed: 11731797]
- Supplemental Online Material.
- Petkov PM, Broman KW, Szatkiewicz JP, Paigen K. Trends Genet Nov;2007 23:539. [PubMed: 17964681]
- Chimpanzee Sequencing and Analysis Consortium. Nature Sep 1;2005 437:69. [PubMed: 16136131]
- Nachman MW, Crowell SL. Genetics Sep;2000 156:297. [PubMed: 10978293]
- Haile-Selassie Y. Nature Jul 12;2001 412:178. [PubMed: 11449272]
- Haile-Selassie Y, Asfaw B, White TD. Am J Phys Anthropol Jan;2004 123:1. [PubMed: 14669231]
- Haile-Selassie Y, Suwa G, White TD. Science Mar 5;2004 303:1503. [PubMed: 15001775]
- Deino AL, Tauxe L, Monaghan M, Hill A. J Hum Evol Jan-Feb;2002 42:117. [PubMed: 11795971]
- Brunet M, et al. Nature Jul 11;2002 418:145. [PubMed: 12110880]
- Chen FC, Li WH. Am J Hum Genet Feb;2001 68:444. [PubMed: 11170892]
- Burgess R, Yang Z. Mol Biol Evol Sep;2008 25:1979. [PubMed: 18603620]
- Takahata N. Jpn J Genet Dec;1993 68:539. [PubMed: 8031575]
- Wall JD. Genetics Jan;2003 163:395. [PubMed: 12586724]
- Kondrashov AS. Hum Mutat Jan;2003 21:12. [PubMed: 12497628]
- Xue Y, et al. Curr Biol Sep 15;2009 19:1453. [PubMed: 19716302]
- Ng SB, et al. Nat Genet Jan;2010 42:30. [PubMed: 19915526]
- Olbrich H, et al. Nat Genet Feb;2002 30:143. [PubMed: 11788826]
- Ng SB, et al. Nature Sep 10;2009 461:272. [PubMed: 19684571]
NIH-PA Author Manuscript
NIH-PA Author Manuscript
NIH-PA Author Manuscript
Figure 1. The landscape of recombination. Each chromosome in this schematic karyotype is used to represent information abstracted from the four corresponding chromosomes of the two children in the pedigree. It is vertically split to indicate the inheritance state from the father (left half) and mother (right half) as shown in the key. The three compound heterozygous ( DHODH, DNAH5, KIAA0556 ) and one recessive ( CES1 ) candidate gene, depicted by red bands, lie in “identical” blocks. Inset: Scatterplot of HapMap recombination rates (in centimorgans per megabase, cM/Mb) within the predicted crossover regions. The maximum value of cM/Mb found in each window is shown in red. The left hand histogram shows the size distribution of recombination windows (log 10 value: -0.58 ± 0.92). The upper graph shows the cM/Mb distribution for the observed maximal values (red), for similarly sized windows shifted by 6 kb (orange), and for similarly sized windows randomly chosen from the entire genome (blue). Note that a shift of 6 kb from the observed locations eliminates the correlation with hotspots. Of 155 recombination windows, 92 contained a HapMap site with
10 cM/Mb. Only five randomly picked windows are expected to contain such high recombination rates.
NIH-PA Author Manuscript
NIH-PA Author Manuscript
NIH-PA Author Manuscript
Figure 2. Power of four. Inheritance states illustrated for a single chromosome in six scenarios representing restrictions of the dataset to the exome (for two siblings only or for the full family) or to subsets of the family (parents and one child, two siblings, siblings and one parent), compared to inheritance state consistency analysis (ISCA) with full data from all four family members. The most supported state for each bin is shown as a color; the height of each histogram bar is proportional to the number of informative markers supporting that state. The father has two regions of homozygosity (thin red lines, bottom panel) on the short arm of the chromosome, where it is not possible to distinguish the haploidentical maternal from identical states (Fig. S2A, panel b). These regions are undetected when the mother's genotypes are missing, because all markers positions in the region are uninformative (second to bottom panel). A pedigree of two parents and one child has only one inheritance state, and so provides no information on recombination. Red, identical; blue, nonidentical; green, haploidentical maternal; yellow, haploidentical paternal. Chromosome structure is annotated as in Fig. 1.
NIH-PA Author Manuscript
NIH-PA Author Manuscript
NIH-PA Author Manuscript
