Implementation and Feasibility of Clinical Genome Sequencing Embedded Into the Outpatient Nephrology Care for Patients With Proteinuric Kidney Disease.

Introduction: The diagnosis and management of proteinuric kidney diseases such as focal segmental glomerulosclerosis (FSGS) are challenging. Genetics holds the promise to improve clinical decision making for these diseases; however, it is often performed too late to enable timely clinical action and it is not implemented within routine outpatient nephrology visits. Methods: We sought to test the implementation and feasibility of clinical rapid genome sequencing (GS) in guiding decision making in patients with proteinuric kidney disease in real-time and embedded in the outpatient nephrology setting. Results: We enrolled 10 children or young adults with biopsy-proven FSGS (9 cases) or minimal change disease (1 case). The mean age at enrollment was 16.2 years (range 2-30). The workflow did not require referral to external genetics clinics but was conducted entirely during the nephrology standard-of-care appointments. The total turn-around-time from enrollment to return-of-results and clinical decision averaged 21.8 days (12.4 for GS), which is well within a time frame that allows clinically relevant treatment decisions. A monogenic or APOL1-related form of kidney disease was diagnosed in 5 of 10 patients. The genetic findings resulted in a rectified diagnosis in 6 patients. Both positive and negative GS findings determined a change in pharmacological treatment. In 3 patients, the results were instrumental for transplant evaluation, donor selection, and the immunosuppressive treatment. All patients and families received genetic counseling. Conclusion: Clinical GS is feasible and can be implemented in real-time in the outpatient care to help guiding clinical management. Additional studies are needed to confirm the cost-effectiveness and broader utility of clinical GS across the phenotypic and demographic spectrum of kidney diseases.

Implementation of Rapid Genome Sequencing for Critically Ill Infants With Complex Congenital Heart Disease.

BACKGROUND: Rapid genome sequencing (rGS) has been shown to improve care of critically ill infants. Congenital heart disease (CHD) is a leading cause of infant mortality and is often caused by genetic disorders, yet the utility of rGS has not been prospectively studied in this population. METHODS: We conducted a prospective evaluation of rGS to improve the care of infants with complex CHD in our cardiac neonatal intensive care unit. RESULTS: In a cohort of 48 infants with complex CHD, rGS diagnosed 14 genetic disorders in 13 (27%) individuals and led to changes in clinical management in 8 (62%) cases with diagnostic results. These included 2 cases in whom genetic diagnoses helped avert intensive, futile interventions before cardiac neonatal intensive care unit discharge, and 3 cases in whom eye disease was diagnosed and treated in early childhood. CONCLUSIONS: Our study provides the first prospective evaluation of rGS for infants with complex CHD to our knowledge. We found that rGS diagnosed genetic disorders in 27% of cases and led to changes in management in 62% of cases with diagnostic results. Our model of care depended on coordination between neonatologists, cardiologists, surgeons, geneticists, and genetic counselors. These findings highlight the important role of rGS in CHD and demonstrate the need for expanded study of how to implement this resource to a broader population of infants with CHD.

Identification of copy number variants with genome sequencing: Clinical experiences from the NYCKidSeq program.

Copy number variations (CNVs) play a significant role in human disease. While chromosomal microarray has traditionally been the first-tier test for CNV detection, use of genome sequencing (GS) is increasing. We report the frequency of CNVs detected with GS in a diverse pediatric cohort from the NYCKidSeq program and highlight specific examples of its clinical impact. A total of 1052 children (0-21 years) with neurodevelopmental, cardiac, and/or immunodeficiency phenotypes received GS. Phenotype-driven analysis was used, resulting in 183 (17.4%) participants with a diagnostic result. CNVs accounted for 20.2% of participants with a diagnostic result (37/183) and ranged from 0.5 kb to 16 Mb. Of participants with a diagnostic result (n = 183) and phenotypes in more than one category, 5/17 (29.4%) were solved by a CNV finding, suggesting a high prevalence of diagnostic CNVs in participants with complex phenotypes. Thirteen participants with a diagnostic CNV (35.1%) had previously uninformative genetic testing, of which nine included a chromosomal microarray. This study demonstrates the benefits of GS for reliable detection of CNVs in a pediatric cohort with variable phenotypes.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

PURPOSE: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. METHODS: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. RESULTS: A total of 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses (P < .001). Yield was greater for GS vs TGPs in Hispanic/Latino(a) (17.2% vs 9.5%, P < .001) and White/European American (19.8% vs 7.9%, P < .001) but not in Black/African American (11.5% vs 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. CONCLUSION: GS may yield up to twice as many diagnoses in pediatric patients compared with TGP testing but not yet across all population groups.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

Purpose: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. Methods: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. Results: 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses ( P < .001). Yield was greater for GS vs . TGPs in Hispanic/Latino(a) (17.2% vs . 9.5%, P < .001) and White/European American (19.8% vs . 7.9%, P < .001), but not in Black/African American (11.5% vs . 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs . White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. Conclusion: GS may yield up to twice as many diagnoses in pediatric patients compared to TGP testing, but not yet across all population groups.

Embryonic lethal genetic variants and chromosomally normal pregnancy loss.

OBJECTIVE: To examine whether rare damaging genetic variants are associated with chromosomally normal pregnancy loss and estimate the magnitude of the association. DESIGN: Case-control. SETTING: Cases were derived from a consecutive series of karyotyped losses at one New Jersey hospital. Controls were derived from the National Database for Autism Research. PATIENT(S): Cases comprised 19 chromosomally normal loss conceptus-parent trios. Controls comprised 547 unaffected siblings of autism case-parent trios. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): The rate of damaging variants in the exome (loss of function and missense-damaging) and the proportions of probands with at least one such variant among cases vs. controls. RESULTS: The proportions of probands with at least one rare damaging variant were 36.8% among cases and 22.9% among controls (odds ratio, 2.0; 99% confidence interval, 0.5-7.3). No case had a variant in a known fetal anomaly gene. The proportion with variants in possibly embryonic lethal genes increased in case probands (odds ratio, 14.5; 99% confidence interval, 1.5-89.7); variants occurred in BAZ1A, FBN2, and TIMP2. CONCLUSION(S): Rare genetic variants in the conceptus may be a cause of chromosomally normal pregnancy loss. A larger sample is needed to estimate the magnitude of the association with precision and identify relevant biologic pathways.

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism.

BACKGROUND: Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737. RESULTS: We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10-3), and combined dataset (p = 1.1 × 10-4). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10-35, loss-of-function p = 2.26 × 10-13) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10-6, OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia. CONCLUSIONS: In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

The NYCKidSeq project: study protocol for a randomized controlled trial incorporating genomics into the clinical care of diverse New York City children.

BACKGROUND: Increasingly, genomics is informing clinical practice, but challenges remain for medical professionals lacking genetics expertise, and in access to and clinical utility of genomic testing for minority and underrepresented populations. The latter is a particularly pernicious problem due to the historical lack of inclusion of racially and ethnically diverse populations in genomic research and genomic medicine. A further challenge is the rapidly changing landscape of genetic tests and considerations of cost, interpretation, and diagnostic yield for emerging modalities like whole-genome sequencing. METHODS: The NYCKidSeq project is a randomized controlled trial recruiting 1130 children and young adults predominantly from Harlem and the Bronx with suspected genetic disorders in three disease categories: neurologic, cardiovascular, and immunologic. Two clinical genetic tests will be performed for each participant, either proband, duo, or trio whole-genome sequencing (depending on sample availability) and proband targeted gene panels. Clinical utility, cost, and diagnostic yield of both testing modalities will be assessed. This study will evaluate the use of a novel, digital platform (GUÍA) to digitize the return of genomic results experience and improve participant understanding for English- and Spanish-speaking families. Surveys will collect data at three study visits: baseline (0 months), result disclosure visit (ROR1, + 3 months), and follow-up visit (ROR2, + 9 months). Outcomes will assess parental understanding of and attitudes toward receiving genomic results for their child and behavioral, psychological, and social impact of results. We will also conduct a pilot study to assess a digital tool called GenomeDiver designed to enhance communication between clinicians and genetic testing labs. We will evaluate GenomeDiver's ability to increase the diagnostic yield compared to standard practices, improve clinician's ability to perform targeted reverse phenotyping, and increase the efficiency of genetic testing lab personnel. DISCUSSION: The NYCKidSeq project will contribute to the innovations and best practices in communicating genomic test results to diverse populations. This work will inform strategies for implementing genomic medicine in health systems serving diverse populations using methods that are clinically useful, technologically savvy, culturally sensitive, and ethically sound. TRIAL REGISTRATION: ClinicalTrials.gov NCT03738098 . Registered on November 13, 2018 Trial Sponsor: Icahn School of Medicine at Mount Sinai Contact Name: Eimear Kenny, PhD (Principal Investigator) Address: Icahn School of Medicine at Mount Sinai, One Gustave L. Levy Pl., Box 1003, New York, NY 10029 Email: eimear.kenny@mssm.edu.

Detection of Copy Number Variants by Short Multiply Aggregated Sequence Homologies.

Chromosomal microarray testing is indicated for patients with diagnoses including unexplained developmental delay or intellectual disability, autism spectrum disorders, and multiple congenital anomalies. The short multiply aggregated sequence homologies (SMASH) genomic assay is a novel next-generation sequencing technology that performs copy number analysis at resolution similar to high-coverage whole genome sequencing but requires far less capacity. We benchmarked the performance of SMASH on a panel of genomic DNAs containing known copy number variants (CNVs). SMASH was able to detect pathogenic copy number variants of ≥10 kb in 77 of 77 samples. No pathogenic events were seen in 32 of 32 controls, indicating 100% sensitivity and specificity for detecting pathogenic CNVs >10 kb. Repeatability (interassay precision) and reproducibility (intra-assay precision) were assessed with 13 samples and showed perfect concordance. We also established that SMASH had a limit of detection of 20% for detection of large mosaic CNVs. Finally, we analyzed seven blinded specimens by SMASH analysis and successfully identified all pathogenic events. These results establish the efficacy of the SMASH genomic assay as a clinical test for the detection of pathogenic copy number variants at a resolution comparable to chromosomal microarray analysis.

Variant Classification Concordance using the ACMG-AMP Variant Interpretation Guidelines across Nine Genomic Implementation Research Studies.

Harmonization of variant pathogenicity classification across laboratories is important for advancing clinical genomics. The two CLIA-accredited Electronic Medical Record and Genomics Network sequencing centers and the six CLIA-accredited laboratories and one research laboratory performing genome or exome sequencing in the Clinical Sequencing Evidence-Generating Research Consortium collaborated to explore current sources of discordance in classification. Eight laboratories each submitted 20 classified variants in the ACMG secondary finding v.2.0 genes. After removing duplicates, each of the 158 variants was annotated and independently classified by two additional laboratories using the ACMG-AMP guidelines. Overall concordance across three laboratories was assessed and discordant variants were reviewed via teleconference and email. The submitted variant set included 28 P/LP variants, 96 VUS, and 34 LB/B variants, mostly in cancer (40%) and cardiac (27%) risk genes. Eighty-six (54%) variants reached complete five-category (i.e., P, LP, VUS, LB, B) concordance, and 17 (11%) had a discordance that could affect clinical recommendations (P/LP versus VUS/LB/B). 21% and 63% of variants submitted as P and LP, respectively, were discordant with VUS. Of the 54 originally discordant variants that underwent further review, 32 reached agreement, for a post-review concordance rate of 84% (118/140 variants). This project provides an updated estimate of variant concordance, identifies considerations for LP classified variants, and highlights ongoing sources of discordance. Continued and increased sharing of variant classifications and evidence across laboratories, and the ongoing work of ClinGen to provide general as well as gene- and disease-specific guidance, will lead to continued increases in concordance.

Distinct roles of BRCA2 in replication fork protection in response to hydroxyurea and DNA interstrand cross-links.

DNA interstrand cross-links (ICLs) are a form of DNA damage that requires the interplay of a number of repair proteins including those of the Fanconi anemia (FA) and the homologous recombination (HR) pathways. Pathogenic variants in the essential gene BRCA2/FANCD1, when monoallelic, predispose to breast and ovarian cancer, and when biallelic, result in a severe subtype of Fanconi anemia. BRCA2 function in the FA pathway is attributed to its role as a mediator of the RAD51 recombinase in HR repair of programmed DNA double-strand breaks (DSB). BRCA2 and RAD51 functions are also required to protect stalled replication forks from nucleolytic degradation during response to hydroxyurea (HU). While RAD51 has been shown to be necessary in the early steps of ICL repair to prevent aberrant nuclease resection, the role of BRCA2 in this process has not been described. Here, based on the analysis of BRCA2 DNA-binding domain (DBD) mutants (c.8488-1G>A and c.8524C>T) discovered in FA patients presenting with atypical FA-like phenotypes, we establish that BRCA2 is necessary for the protection of DNA at ICLs. Cells carrying BRCA2 DBD mutations are sensitive to ICL-inducing agents but resistant to HU treatment consistent with relatively high HR repair in these cells. BRCA2 function at an ICL protects against DNA2-WRN nuclease-helicase complex and not the MRE11 nuclease that is implicated in the resection of HU-induced stalled replication forks. Our results also indicate that unlike the processing at HU-induced stalled forks, the function of the SNF2 translocases (SMARCAL1, ZRANB3, or HLTF), implicated in fork reversal, are not an integral component of the ICL repair, pointing to a different mechanism of fork protection at different DNA lesions.

Prenatal Diagnosis by Whole Exome Sequencing in Fetuses with Ultrasound Abnormalities.

Whole-exome sequencing (WES) has been used as a standard of care for postnatal diagnosis in the clinical setting in the past few years for children and adults with undiagnosed disease. Many rare disorders have been diagnosed through WES, which is less expensive than the traditional serial genetic testing where patients had previously spent years on an uninformative diagnostic odyssey. Seeking a diagnosis often entails enduring time consuming, and sometimes invasive procedures which may be associated with medical risks that are stressful for families and impose a heavy burden on the health-care system. However, the use of WES is considered impractical in the prenatal and neonatal testing period because of the technical and computational challenges of performing genomic sequencing from small amounts of genetic material, and the need for faster turnaround time (TAT) than the current 6-8 weeks TAT provided by most clinical labs offering postnatal testing. With the rapidly evolving methods of sequence analysis, there are clinical challenges such as the constantly increasing number of genes being identified which are not yet fully phenotypically characterized, especially when ascertained prenatally or neonatally before all the clinical features may be evident. Despite these challenges, there are many clinical benefits to acquiring genomic information in the prenatal and neonatal period. These include superior prognostic information which allows for prenatal planning of mode of delivery and hospital for delivery and optimized neonatal management. We have developed a clinical WES assay using small amounts of DNA with a TAT of 10 days for use in the prenatal or neonatal setting. This test is used to detect small nucleotide variants and indels in fetuses and neonates with structural abnormalities.

Analytical Validation of Clinical Whole-Genome and Transcriptome Sequencing of Patient-Derived Tumors for Reporting Targetable Variants in Cancer.

We developed and validated a clinical whole-genome and transcriptome sequencing (WGTS) assay that provides a comprehensive genomic profile of a patient's tumor. The ability to fully capture the mappable genome with sufficient sequencing coverage to precisely call DNA somatic single nucleotide variants, insertions/deletions, copy number variants, structural variants, and RNA gene fusions was analyzed. New York State's Department of Health next-generation DNA sequencing guidelines were expanded for establishing performance validation applicable to whole-genome and transcriptome sequencing. Whole-genome sequencing laboratory protocols were validated for the Illumina HiSeq X Ten platform and RNA sequencing for Illumina HiSeq2500 platform for fresh or frozen and formalin-fixed, paraffin-embedded tumor samples. Various bioinformatics tools were also tested, and CIs for sensitivity and specificity thresholds in calling clinically significant somatic aberrations were determined. The validation was performed on a set of 125 tumor normal pairs. RNA sequencing was performed to call fusions and to confirm the DNA variants or exonic alterations. Here, we present our results and WGTS standards for variant allele frequency, reproducibility, analytical sensitivity, and present limit of detection analysis for single nucleotide variant calling, copy number identification, and structural variants. We show that The New York Genome Center WGTS clinical assay can provide a comprehensive patient variant discovery approach suitable for directed oncologic therapeutic applications.

CDG: An Online Server for Detecting Biologically Closest Disease-Causing Genes and its Application to Primary Immunodeficiency.

High-throughput genomic technologies yield about 20,000 variants in the protein-coding exome of each individual. A commonly used approach to select candidate disease-causing variants is to test whether the associated gene has been previously reported to be disease-causing. In the absence of known disease-causing genes, it can be challenging to associate candidate genes with specific genetic diseases. To facilitate the discovery of novel gene-disease associations, we determined the putative biologically closest known genes and their associated diseases for 13,005 human genes not currently reported to be disease-associated. We used these data to construct the closest disease-causing genes (CDG) server, which can be used to infer the closest genes with an associated disease for a user-defined list of genes or diseases. We demonstrate the utility of the CDG server in five immunodeficiency patient exomes across different diseases and modes of inheritance, where CDG dramatically reduced the number of candidate genes to be evaluated. This resource will be a considerable asset for ascertaining the potential relevance of genetic variants found in patient exomes to specific diseases of interest. The CDG database and online server are freely available to non-commercial users at: http://lab.rockefeller.edu/casanova/CDG.

Clinical whole exome sequencing from dried blood spot identifies novel genetic defect underlying asparagine synthetase deficiency.

We add two novel variants to the existing mutation spectrum of ASNS gene. Loss of ASNS function should be suspected in newborns presenting with congenital microcephaly, intellectual disability, progressive cerebral atrophy, and intractable seizures. Acquisition and sequencing of stored newborn blood spot can be a valuable option when no biological samples are available from a deceased child.

Human Adaptive Immunity Rescues an Inborn Error of Innate Immunity.

The molecular basis of the incomplete penetrance of monogenic disorders is unclear. We describe here eight related individuals with autosomal recessive TIRAP deficiency. Life-threatening staphylococcal disease occurred during childhood in the proband, but not in the other seven homozygotes. Responses to all Toll-like receptor 1/2 (TLR1/2), TLR2/6, and TLR4 agonists were impaired in the fibroblasts and leukocytes of all TIRAP-deficient individuals. However, the whole-blood response to the TLR2/6 agonist staphylococcal lipoteichoic acid (LTA) was abolished only in the index case individual, the only family member lacking LTA-specific antibodies (Abs). This defective response was reversed in the patient, but not in interleukin-1 receptor-associated kinase 4 (IRAK-4)-deficient individuals, by anti-LTA monoclonal antibody (mAb). Anti-LTA mAb also rescued the macrophage response in mice lacking TIRAP, but not TLR2 or MyD88. Thus, acquired anti-LTA Abs rescue TLR2-dependent immunity to staphylococcal LTA in individuals with inherited TIRAP deficiency, accounting for incomplete penetrance. Combined TIRAP and anti-LTA Ab deficiencies underlie staphylococcal disease in this patient.

Differential burden of rare protein truncating variants in Alzheimer's disease patients compared to centenarians.

We compared coding region variants of 53 cognitively healthy centenarians and 45 patients with Alzheimer's disease (AD), all of Ashkenazi Jewish (AJ) ancestry. Despite the small sample size, the known AD risk variant APOE4 reached genome-wide significance, indicating the advantage of utilizing 'super-controls'. We restricted our subsequent analysis to rare variants observed at most once in the 1000 Genomes database and having a minor allele frequency below 2% in our AJ sample. We compared the burden of predicted protein altering variants between cases and controls as normalized by the level of rare synonymous variants. We observed an increased burden among AD subjects for predicted loss-of-function (LoFs) variants defined as stop-gain, frame shift, initiation codon (INIT) and splice site mutations (n = 930, OR = 1.3, P = 1.5×E-5). There was no enrichment across all rare protein altering variants defined as missense plus LoFs, in frame indels and stop-loss variants (n = 13 014, OR = 0.97, P = 0.47). Among LoFs, the strongest burden was observed for INIT (OR = 2.16, P = 0.0097) and premature stop variants predicted to cause non-sense-mediated decay in the majority of transcripts (NMD) (OR = 1.98, P = 0.02). Notably, this increased burden of NMD, INIT and splice variants was more pronounced in a set of 1397 innate immune genes (OR = 4.55, P = 0.0043). Further comparison to additional exomes indicates that the difference in LoF burden originated both from the AD and centenarian sample. In summary, we observed an overall increased burden of rare LoFs in AD subjects as compared to centenarians, and this enrichment is more pronounced for innate immune genes.

Whole-exome sequencing to analyze population structure, parental inbreeding, and familial linkage.

Principal component analysis (PCA), homozygosity rate estimations, and linkage studies in humans are classically conducted through genome-wide single-nucleotide variant arrays (GWSA). We compared whole-exome sequencing (WES) and GWSA for this purpose. We analyzed 110 subjects originating from different regions of the world, including North Africa and the Middle East, which are poorly covered by public databases and have high consanguinity rates. We tested and applied a number of quality control (QC) filters. Compared with GWSA, we found that WES provided an accurate prediction of population substructure using variants with a minor allele frequency > 2% (correlation = 0.89 with the PCA coordinates obtained by GWSA). WES also yielded highly reliable estimates of homozygosity rates using runs of homozygosity with a 1,000-kb window (correlation = 0.94 with the estimates provided by GWSA). Finally, homozygosity mapping analyses in 15 families including a single offspring with high homozygosity rates showed that WES provided 51% less genome-wide linkage information than GWSA overall but 97% more information for the coding regions. At the genome-wide scale, 76.3% of linked regions were found by both GWSA and WES, 17.7% were found by GWSA only, and 6.0% were found by WES only. For coding regions, the corresponding percentages were 83.5%, 7.4%, and 9.1%, respectively. With appropriate QC filters, WES can be used for PCA and adjustment for population substructure, estimating homozygosity rates in individuals, and powerful linkage analyses, particularly in coding regions.