A new approach for rare variation collapsing on functional protein domains implicates specific genic regions in ALS

The Current Potential Pathogenesis of Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease mainly characterized by the accumulation of ubiquitinated proteins in the affected motor neurons. At present, the accurate pathogenesis of ALS remains unclear and there are still no effective treatment measures for ALS. The potential pathogenesis of ALS mainly includes the misfolding of some pathogenic proteins, the genetic variation, mitochondrial dysfunction, autophagy disorders, neuroinflammation, the misregulation of RNA, the altered axonal transport, and gut microbial dysbiosis. Exploring the pathogenesis of ALS is a critical step in searching for the effective therapeutic approaches. The current studies suggested that the genetic variation, gut microbial dysbiosis, the activation of glial cells, and the transportation disorder of extracellular vesicles may play some important roles in the pathogenesis of ALS. This review conducts a systematic review of these current potential promising topics closely related to the pathogenesis of ALS; it aims to provide some new evidences and clues for searching the novel treatment measures of ALS.

DYNATE: Localizing rare-variant association regions via multiple testing embedded in an aggregation tree

Rare-variants (RVs) genetic association studies enable researchers to uncover the variation in phenotypic traits left unexplained by common variation. Traditional single-variant analysis lacks power; thus, researchers have developed various methods to aggregate the effects of RVs across genomic regions to study their collective impact. Some existing methods utilize a static delineation of genomic regions, often resulting in suboptimal effect aggregation, as neutral subregions within the test region will result in an attenuation of signal. Other methods use varying windows to search for signals but often result in long regions containing many neutral RVs. To pinpoint short genomic regions enriched for disease-associated RVs, we developed a novel method, DYNamic Aggregation TEsting (DYNATE). DYNATE dynamically and hierarchically aggregates smaller genomic regions into larger ones and performs multiple testing for disease associations with a controlled weighted false discovery rate. DYNATE's main advantage lies in its strong ability to identify short genomic regions highly enriched for disease-associated RVs. Extensive numerical simulations demonstrate the superior performance of DYNATE under various scenarios compared with existing methods. We applied DYNATE to an amyotrophic lateral sclerosis study and identified a new gene, EPG5, harboring possibly pathogenic mutations.

Microtubule association of TRIM3 revealed by differential extraction proteomics

The microtubule network is formed from polymerised tubulin subunits and associating proteins, which govern microtubule dynamics and a diverse array of functions. To identify novel microtubule-binding proteins, we have developed an unbiased biochemical assay, which relies on the selective extraction of cytosolic proteins from U2OS cells, while leaving behind the microtubule network. Candidate proteins are linked to microtubules by their sensitivities to the depolymerising drug nocodazole or the microtubule-stabilising drug taxol, which is quantitated by mass spectrometry. Our approach is benchmarked by co-segregation of tubulin and previously established microtubule-binding proteins. We then identify several novel candidate microtubule-binding proteins, from which we have selected the ubiquitin E3 ligase tripartite motif-containing protein 3 (TRIM3) for further characterisation. We map TRIM3 microtubule binding to its C-terminal NHL-repeat region. We show that TRIM3 is required for the accumulation of acetylated tubulin, following treatment with taxol. Furthermore, loss of TRIM3 partially recapitulates the reduction in nocodazole-resistant microtubules characteristic of α-tubulin acetyltransferase 1 (ATAT1) depletion. These results can be explained by a decrease in ATAT1 following depletion of TRIM3 that is independent of transcription.

Rare variant analyses validate known ALS genes in a multi-ethnic population and identifies ANTXR2 as a candidate in PLS

Background

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease affecting over 30,000 people in the United States. It is characterized by the progressive decline of the nervous system that leads to the weakening of muscles which impacts physical function. Approximately, 15% of individuals diagnosed with ALS have a known genetic variant that contributes to their disease. As therapies that slow or prevent symptoms, such as antisense oligonucleotides, continue to develop, it is important to discover novel genes that could be targets for treatment. Additionally, as cohorts continue to grow, performing analyses in ALS subtypes, such as primary lateral sclerosis (PLS), becomes possible due to an increase in power. These analyses could highlight novel pathways in disease manifestation.

Methods

Building on our previous discoveries using rare variant association analyses, we conducted rare variant burden testing on a substantially larger cohort of 6,970 ALS patients from a large multi-ethnic cohort as well as 166 PLS patients, and 22,524 controls. We used intolerant domain percentiles based on sub-region Residual Variation Intolerance Score (subRVIS) that have been described previously in conjunction with gene based collapsing approaches to conduct burden testing to identify genes that associate with ALS and PLS.

Results

A gene based collapsing model showed significant associations with SOD1, TARDBP, and TBK1 (OR=19.18, p = 3.67 × 10-39; OR=4.73, p = 2 × 10-10; OR=2.3, p = 7.49 × 10-9, respectively). These genes have been previously associated with ALS. Additionally, a significant novel control enriched gene, ALKBH3 (p = 4.88 × 10-7), was protective for ALS in this model. An intolerant domain based collapsing model showed a significant improvement in identifying regions in TARDBP that associated with ALS (OR=10.08, p = 3.62 × 10-16). Our PLS protein truncating variant collapsing analysis demonstrated significant case enrichment in ANTXR2 (p=8.38 × 10-6).

Conclusions

In a large multi-ethnic cohort of 6,970 ALS patients, rare variant burden testing validated known ALS genes and identified a novel potentially protective gene, ALKBH3. A first-ever analysis in 166 patients with PLS found a candidate association with loss-of-function mutations in ANTXR2.

Rare variant analyses validate known ALS genes in a multi-ethnic population and identifies ANTXR2 as a candidate in PLS

Background

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease affecting over 30,000 people in the United States. It is characterized by the progressive decline of the nervous system that leads to the weakening of muscles which impacts physical function. Approximately, 15% of individuals diagnosed with ALS have a known genetic variant that contributes to their disease. As therapies that slow or prevent symptoms, such as antisense oligonucleotides, continue to develop, it is important to discover novel genes that could be targets for treatment. Additionally, as cohorts continue to grow, performing analyses in ALS subtypes, such as primary lateral sclerosis (PLS), becomes possible due to an increase in power. These analyses could highlight novel pathways in disease manifestation.

Methods

Building on our previous discoveries using rare variant association analyses, we conducted rare variant burden testing on a substantially larger cohort of 6,970 ALS patients from a large multi-ethnic cohort as well as 166 PLS patients, and 22,524 controls. We used intolerant domain percentiles based on sub-region Residual Variation Intolerance Score (subRVIS) that have been described previously in conjunction with gene based collapsing approaches to conduct burden testing to identify genes that associate with ALS and PLS.

Results

A gene based collapsing model showed significant associations with SOD1, TARDBP, and TBK1 (OR=19.18, p = 3.67 × 10-39; OR=4.73, p = 2 × 10-10; OR=2.3, p = 7.49 × 10-9, respectively). These genes have been previously associated with ALS. Additionally, a significant novel control enriched gene, ALKBH3 (p = 4.88 × 10-7), was protective for ALS in this model. An intolerant domain based collapsing model showed a significant improvement in identifying regions in TARDBP that associated with ALS (OR=10.08, p = 3.62 × 10-16). Our PLS protein truncating variant collapsing analysis demonstrated significant case enrichment in ANTXR2 (p=8.38 × 10-6).

Conclusions

In a large multi-ethnic cohort of 6,970 ALS patients, rare variant burden testing validated known ALS genes and identified a novel potentially protective gene, ALKBH3. A first-ever analysis in 166 patients with PLS found a candidate association with loss-of-function mutations in ANTXR2.

Amyotrophic Lateral Sclerosis Risk Genes and Suppressor

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease that leads to death by progressive paralysis and respiratory failure within 2-4 years of onset. About 90-95% of ALS cases are sporadic (sALS), and 5-10% are inherited through family (fALS). Though the mechanisms of the disease are still poorly understood, so far, approximately 40 genes have been reported as ALS causative genes. The mutations in some crucial genes, like SOD1, C9ORF72, FUS, and TDP-43, are majorly associated with ALS, resulting in ROS-associated oxidative stress, excitotoxicity, protein aggregation, altered RNA processing, axonal and vesicular trafficking dysregulation, and mitochondrial dysfunction. Recent studies show that dysfunctional cellular pathways get restored as a result of the repair of a single pathway in ALS. In this review article, our aim is to identify putative targets for therapeutic development and the importance of a single suppressor to reduce multiple symptoms by focusing on important mutations and the phenotypic suppressors of dysfunctional cellular pathways in crucial genes as reported by other studies.

Whole exome sequencing in Alopecia Areata identifies rare variants in KRT82

Alopecia areata is a complex genetic disease that results in hair loss due to the autoimmune-mediated attack of the hair follicle. We previously defined a role for both rare and common variants in our earlier GWAS and linkage studies. Here, we identify rare variants contributing to Alopecia Areata using a whole exome sequencing and gene-level burden analyses approach on 849 Alopecia Areata patients compared to 15,640 controls. KRT82 is identified as an Alopecia Areata risk gene with rare damaging variants in 51 heterozygous Alopecia Areata individuals (6.01%), achieving genome-wide significance (p = 2.18E−07). KRT82 encodes a hair-specific type II keratin that is exclusively expressed in the hair shaft cuticle during anagen phase, and its expression is decreased in Alopecia Areata patient skin and hair follicles. Finally, we find that cases with an identified damaging KRT82 variant and reduced KRT82 expression have elevated perifollicular CD8 infiltrates. In this work, we utilize whole exome sequencing to successfully identify a significant Alopecia Areata disease-relevant gene, KRT82, and reveal a proposed mechanism for rare variant predisposition leading to disrupted hair shaft integrity.

Common variants have been discovered to be associated with Alopecia Areata; however, rare variants have been less well studied. Here, the authors use whole-exome sequencing to identify associated rare variants in the hair keratin gene KRT82. Further, they find that individuals with Alopecia Areata have reduced expression of KRT82 in the skin and hair follicle.

Constitutively active SARM1 variants that induce neuropathy are enriched in ALS patients

Background

In response to injury, neurons activate a program of organized axon self-destruction initiated by the NAD⁺ hydrolase, SARM1. In healthy neurons SARM1 is autoinhibited, but single amino acid changes can abolish autoinhibition leading to constitutively active SARM1 enzymes that promote degeneration when expressed in cultured neurons.

Methods

To investigate whether naturally occurring human variants might disrupt SARM1 autoinhibition and potentially contribute to risk for neurodegenerative disease, we assayed the enzymatic activity of all 42 rare SARM1 alleles identified among 8507 amyotrophic lateral sclerosis (ALS) patients and 9671 controls. We then intrathecally injected mice with virus expressing SARM1 constructs to test the capacity of an ALS-associated constitutively active SARM1 variant to promote neurodegeneration in vivo.

Results

Twelve out of 42 SARM1 missense variants or small in-frame deletions assayed exhibit constitutive NADase activity, including more than half of those that are unique to the ALS patients or that occur in multiple patients. There is a > 5-fold enrichment of constitutively active variants among patients compared to controls. Expression of constitutively active ALS-associated SARM1 alleles in cultured dorsal root ganglion (DRG) neurons is pro-degenerative and cytotoxic. Intrathecal injection of an AAV expressing the common SARM1 reference allele is innocuous to mice, but a construct harboring SARM1^V184G, the constitutively active variant found most frequently among the ALS patients, causes axon loss, motor dysfunction, and sustained neuroinflammation.

Conclusions

These results implicate rare hypermorphic SARM1 alleles as candidate genetic risk factors for ALS and other neurodegenerative conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13024-021-00511-x.

Motoneuron Diseases

Motoneuron diseases (MNDs) represent a heterogeneous group of progressive paralytic disorders, mainly characterized by the loss of upper (corticospinal) motoneurons, lower (spinal) motoneurons or, often both. MNDs can occur from birth to adulthood and have a highly variable clinical presentation, even within gene-positive forms, suggesting the existence of environmental and genetic modifiers. A combination of cell autonomous and non-cell autonomous mechanisms contributes to motoneuron degeneration in MNDs, suggesting multifactorial pathogenic processes.

Protein network analysis to prioritize key genes in amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis (ALS) is a fatal disease, progressive nature characterizes by loss of both upper and lower motor neuron functions. One of the major challenge is to understand the mechanism of ALS multifactorial nature. We aimed to explore some key genes related to ALS through bioinformatics methods for its therapeutic intervention. Here, we applied a systems biology approach involving experimentally validated 148 ALS-associated proteins and construct ALS protein-protein interaction network (ALS-PPIN). The network was further statistically analysed and identified bottleneck-hubs. The network is also subjected to identify modules which could have similar functions. The interaction between the modules and bottleneck-hubs provides the functional regulatory role of the ALS mechanism. The ALS-PPIN demonstrated a hierarchical scale-free nature. We identified 17 bottleneck-hubs, in which CDC5L, SNW1, TP53, SOD1, and VCP were the high degree nodes (hubs) in ALS-PPIN. CDC5L was found to control highly cluster modules and play a vital role in the stability of the overall network followed by SNW1, TP53, SOD1, and VCP. HSPA5 and HSPA8 acting as a common connector for CDC5L and TP53 bottleneck-hubs. The functional and disease association analysis showed ALS has a strong correlation with mRNA processing, protein deubiquitination, and neoplasms, nervous system, immune system disease classes. In the future, biochemical investigation of the observed bottleneck-hubs and their interacting partners could provide a further understanding of their role in the pathophysiology of ALS.

•

Amyotrophic Lateral Sclerosis protein-protein interaction network (ALS-PPIN) followed a hierarchical scale-free nature.

•

We identified 17 bottleneck-hubs in the ALS-PPIN.

•

Among bottleneck-hubs we found CDC5L, SNW1, TP53, SOD1, and VCP were the high degree nodes (hubs) in the ALS-PPIN.

•

CDC5L is the effective communicator with all five modules in the ALS-PPIN and followed by SNW1 and TP53.

•

Modules are highly associated with various disease classes like neoplasms, nervous systems and others.

End-Truncated LAMB1 Causes a Hippocampal Memory Defect and a Leukoencephalopathy

OBJECTIVE

The majority of patients with a familial cerebral small vessel disease (CSVD) referred for molecular screening do not show pathogenic variants in known genes. In this study, we aimed to identify novel CSVD causal genes.

METHODS

We performed a gene-based collapsing test of rare protein-truncating variants identified in exome data of 258 unrelated CSVD patients of an ethnically matched control cohort and of 2 publicly available large-scale databases, gnomAD and TOPMed. Western blotting was used to investigate the functional consequences of variants. Clinical and magnetic resonance imaging features of mutated patients were characterized.

RESULTS

We showed that LAMB1 truncating variants escaping nonsense-mediated messenger RNA decay are strongly overrepresented in CSVD patients, reaching genome-wide significance (p < 5 × 10^-8 ). Using 2 antibodies recognizing the N- and C-terminal parts of LAMB1, we showed that truncated forms of LAMB1 are expressed in the endogenous fibroblasts of patients and trapped in the cytosol. These variants are associated with a novel phenotype characterized by the association of a hippocampal type episodic memory defect and a diffuse vascular leukoencephalopathy.

INTERPRETATION

These findings are important for diagnosis and clinical care, to avoid unnecessary and sometimes invasive investigations, and also from a mechanistic point of view to understand the role of extracellular matrix proteins in neuronal homeostasis. ANN NEUROL 2021;90:962-975.

Sub-genic intolerance, ClinVar, and the epilepsies: A whole-exome sequencing study of 29,165 individuals

Both mild and severe epilepsies are influenced by variants in the same genes, yet an explanation for the resulting phenotypic variation is unknown. As part of the ongoing Epi25 Collaboration, we performed a whole-exome sequencing analysis of 13,487 epilepsy-affected individuals and 15,678 control individuals. While prior Epi25 studies focused on gene-based collapsing analyses, we asked how the pattern of variation within genes differs by epilepsy type. Specifically, we compared the genetic architectures of severe developmental and epileptic encephalopathies (DEEs) and two generally less severe epilepsies, genetic generalized epilepsy and non-acquired focal epilepsy (NAFE). Our gene-based rare variant collapsing analysis used geographic ancestry-based clustering that included broader ancestries than previously possible and revealed novel associations. Using the missense intolerance ratio (MTR), we found that variants in DEE-affected individuals are in significantly more intolerant genic sub-regions than those in NAFE-affected individuals. Only previously reported pathogenic variants absent in available genomic datasets showed a significant burden in epilepsy-affected individuals compared with control individuals, and the ultra-rare pathogenic variants associated with DEE were located in more intolerant genic sub-regions than variants associated with non-DEE epilepsies. MTR filtering improved the yield of ultra-rare pathogenic variants in affected individuals compared with control individuals. Finally, analysis of variants in genes without a disease association revealed a significant burden of loss-of-function variants in the genes most intolerant to such variation, indicating additional epilepsy-risk genes yet to be discovered. Taken together, our study suggests that genic and sub-genic intolerance are critical characteristics for interpreting the effects of variation in genes that influence epilepsy.

ATAV: a comprehensive platform for population-scale genomic analyses

Background

A common approach for sequencing studies is to do joint-calling and store variants of all samples in a single file. If new samples are continually added or controls are re-used for several studies, the cost and time required to perform joint-calling for each analysis can become prohibitive.

Results

We present ATAV, an analysis platform for large-scale whole-exome and whole-genome sequencing projects. ATAV stores variant and per site coverage data for all samples in a centralized database, which is efficiently queried by ATAV to support diagnostic analyses for trios and singletons, as well as rare-variant collapsing analyses for finding disease associations in complex diseases. Runtime logs ensure full reproducibility and the modularized ATAV framework makes it extensible to continuous development. Besides helping with the identification of disease-causing variants for a range of diseases, ATAV has also enabled the discovery of disease-genes by rare-variant collapsing on datasets containing more than 20,000 samples. Analyses to date have been performed on data of more than 110,000 individuals demonstrating the scalability of the framework.

To allow users to easily access variant-level data directly from the database, we provide a web-based interface, the ATAV data browser (http://atavdb.org/). Through this browser, summary-level data for more than 40,000 samples can be queried by the general public representing a mix of cases and controls of diverse ancestries. Users have access to phenotype categories of variant carriers, as well as predicted ancestry, gender, and quality metrics. In contrast to many other platforms, the data browser is able to show data of newly-added samples in real-time and therefore evolves rapidly as more and more samples are sequenced.

Conclusions

Through ATAV, users have public access to one of the largest variant databases for patients sequenced at a tertiary care center and can look up any genes or variants of interest. Additionally, since the entire code is freely available on GitHub, ATAV can easily be deployed by other groups that wish to build their own platform, database, and user interface.

Serine biosynthesis defect due to haploinsufficiency of PHGDH causes retinal disease

Macular telangiectasia type 2 (MacTel) is a progressive, late-onset retinal degenerative disease linked to decreased serum levels of serine that elevate circulating levels of a toxic ceramide species, deoxysphingolipids (deoxySLs); however, causal genetic variants that reduce serine levels in patients have not been identified. Here we identify rare, functional variants in the gene encoding the rate-limiting serine biosynthetic enzyme, phosphoglycerate dehydrogenase (PHGDH), as the single locus accounting for a significant fraction of MacTel. Under a dominant collapsing analysis model of a genome-wide enrichment analysis of rare variants predicted to impact protein function in 793 MacTel cases and 17,610 matched controls, the PHGDH gene achieves genome-wide significance (P = 1.2 × 10^-13) with variants explaining ~3.2% of affected individuals. We further show that the resulting functional defects in PHGDH cause decreased serine biosynthesis and accumulation of deoxySLs in retinal pigmented epithelial cells. PHGDH is a significant locus for MacTel that explains the typical disease phenotype and suggests a number of potential treatment options.

Unique roles of rare variants in the genetics of complex diseases in humans

Genome-wide association studies have identified >10,000 genetic variants associated with various phenotypes and diseases. Although the majority are common variants, rare variants with >0.1% of minor allele frequency have been investigated by imputation and using disease-specific custom SNP arrays. Rare variants sequencing analysis mainly revealed have played unique roles in the genetics of complex diseases in humans due to their distinctive features, in contrast to common variants. Unique roles are hypothesis-free evidence for gene causality, a precise target of functional analysis for understanding disease mechanisms, a new favorable target for drug development, and a genetic marker with high disease risk for personalized medicine. As whole-genome sequencing continues to identify more rare variants, the roles associated with rare variants will also increase. However, a better estimation of the functional impact of rare variants across whole genome is needed to enhance their contribution to improvements in human health.

Genetics of Amyotrophic Lateral Sclerosis

Purpose of Review

Amyotrophic lateral sclerosis and frontotemporal dementia (ALS-FTD) spectrum disorder is a rare fatal disease with strong genetic influences. The implementation of short-read sequencing methodologies in increasingly large patient cohorts has rapidly expanded our knowledge of the complex genetic architecture of the disease. We aim to convey the broad history of ALS gene discovery as context for a focused review of 11 ALS gene associations reported over the last 5 years. We also summarize the current level of genetic evidence for all previously reported genes.

Recent Findings

The history of ALS gene discovery has occurred in at least four identifiable phases, each powered by different technologies and scale of investigation. The most recent epoch, benefitting from population-scale genome data, large international consortia, and low-cost sequencing, has yielded 11 new gene associations. We summarize the current level of genetic evidence supporting these ALS genes, highlighting any genotype-phenotype or genotype-pathology correlations, and discussing preliminary understanding of molecular pathogenesis. This era has also raised uncertainty around prior ALS-associated genes and clarified the role of others.

Summary

Our understanding of the genetic underpinning of ALS has expanded rapidly over the last 25 years and has led directly to the clinical application of molecularly driven therapies. Ongoing sequencing efforts in ALS will identify new causative and risk factor genes while clarifying the status of genes reported in prior eras of research.

Causal Genetic Variants in Stillbirth

BACKGROUND

In the majority of cases, the cause of stillbirth remains unknown despite detailed clinical and laboratory evaluation. Approximately 10 to 20% of stillbirths are attributed to chromosomal abnormalities. However, the causal nature of single-nucleotide variants and small insertions and deletions in exomes has been understudied.

METHODS

We generated exome sequencing data for 246 stillborn cases and followed established guidelines to identify causal variants in disease-associated genes. These genes included those that have been associated with stillbirth and strong candidate genes. We also evaluated the contribution of 18,653 genes in case-control analyses stratified according to the degree of depletion of functional variation (described here as "intolerance" to variation).

RESULTS

We identified molecular diagnoses in 15 of 246 cases of stillbirth (6.1%) involving seven genes that have been implicated in stillbirth and six disease genes that are good candidates for phenotypic expansion. Among the cases we evaluated, we also found an enrichment of loss-of-function variants in genes that are intolerant to such variation in the human population (odds ratio, 2.15; 95% confidence interval [CI], 1.46 to 3.06). Loss-of-function variants in intolerant genes were concentrated in genes that have not been associated with human disease (odds ratio, 2.22; 95% CI, 1.41 to 3.34), findings that differ from those in two postnatal clinical populations that were also evaluated in this study.

CONCLUSIONS

Our findings establish the diagnostic utility of clinical exome sequencing to evaluate the role of small genomic changes in stillbirth. The strength of the novel risk signal (as generated through the stratified analysis) was similar to that in known disease genes, which indicates that the genetic cause of stillbirth remains largely unknown. (Funded by the Institute for Genomic Medicine.).

Exome sequencing in genetic disease: recent advances and considerations

Over the past decade, exome sequencing (ES) has allowed significant advancements to the field of disease research. By targeting the protein-coding regions of the genome, ES combines the depth of knowledge on protein-altering variants with high-throughput data generation and ease of analysis. New discoveries continue to be made using ES, and medical science has benefitted both theoretically and clinically from its continued use. In this review, we describe recent advances and successes of ES in disease research. Through selected examples of recent publications, we explore how ES continues to be a valuable tool to find variants that might explain disease etiology or provide insight into the biology underlying the disease. We then discuss shortcomings of ES in terms of variant discoveries made by other sequencing technologies that would be missed because of the scope and techniques of ES. We conclude with a brief outlook on the future of ES, suggesting that although newer and more thorough sequencing methods will soon supplant ES, its results will continue to be useful for disease research.

Rare-variant collapsing analyses for complex traits: guidelines and applications

The first phase of genome-wide association studies (GWAS) assessed the role of common variation in human disease. Advances optimizing and economizing high-throughput sequencing have enabled a second phase of association studies that assess the contribution of rare variation to complex disease in all protein-coding genes. Unlike the early microarray-based studies, sequencing-based studies catalogue the full range of genetic variation, including the evolutionarily youngest forms. Although the experience with common variants helped establish relevant standards for genome-wide studies, the analysis of rare variation introduces several challenges that require novel analysis approaches.