Exome sequencing and the genetic basis of complex traits. Nat Genet. 2012;44:623-630. [PMID: 22641211 DOI: 10.1038/ng.2303]

Rare Functional Variant in TM2D3 is Associated with Late-Onset Alzheimer's Disease

We performed an exome-wide association analysis in 1393 late-onset Alzheimer's disease (LOAD) cases and 8141 controls from the CHARGE consortium. We found that a rare variant (P155L) in TM2D3 was enriched in Icelanders (~0.5% versus <0.05% in other European populations). In 433 LOAD cases and 3903 controls from the Icelandic AGES sub-study, P155L was associated with increased risk and earlier onset of LOAD [odds ratio (95% CI) = 7.5 (3.5-15.9), p = 6.6x10-9]. Mutation in the Drosophila TM2D3 homolog, almondex, causes a phenotype similar to loss of Notch/Presenilin signaling. Human TM2D3 is capable of rescuing these phenotypes, but this activity is abolished by P155L, establishing it as a functionally damaging allele. Our results establish a rare TM2D3 variant in association with LOAD susceptibility, and together with prior work suggests possible links to the β-amyloid cascade.

Investigation of rare and low-frequency variants using high-throughput sequencing with pooled DNA samples

High-throughput sequencing using pooled DNA samples can facilitate genome-wide studies on rare and low-frequency variants in a large population. Some major questions concerning the pooling sequencing strategy are whether rare and low-frequency variants can be detected reliably, and whether estimated minor allele frequencies (MAFs) can represent the actual values obtained from individually genotyped samples. In this study, we evaluated MAF estimates using three variant detection tools with two sets of pooled whole exome sequencing (WES) and one set of pooled whole genome sequencing (WGS) data. Both GATK and Freebayes displayed high sensitivity, specificity and accuracy when detecting rare or low-frequency variants. For the WGS study, 56% of the low-frequency variants in Illumina array have identical MAFs and 26% have one allele difference between sequencing and individual genotyping data. The MAF estimates from WGS correlated well (r = 0.94) with those from Illumina arrays. The MAFs from the pooled WES data also showed high concordance (r = 0.88) with those from the individual genotyping data. In conclusion, the MAFs estimated from pooled DNA sequencing data reflect the MAFs in individually genotyped samples well. The pooling strategy can thus be a rapid and cost-effective approach for the initial screening in large-scale association studies.

Distinct genetic architectures for syndromic and nonsyndromic congenital heart defects identified by exome sequencing

Congenital heart defects (CHDs) have a neonatal incidence of 0.8-1% (refs. 1,2). Despite abundant examples of monogenic CHD in humans and mice, CHD has a low absolute sibling recurrence risk (∼2.7%), suggesting a considerable role for de novo mutations (DNMs) and/or incomplete penetrance. De novo protein-truncating variants (PTVs) have been shown to be enriched among the 10% of 'syndromic' patients with extra-cardiac manifestations. We exome sequenced 1,891 probands, including both syndromic CHD (S-CHD, n = 610) and nonsyndromic CHD (NS-CHD, n = 1,281). In S-CHD, we confirmed a significant enrichment of de novo PTVs but not inherited PTVs in known CHD-associated genes, consistent with recent findings. Conversely, in NS-CHD we observed significant enrichment of PTVs inherited from unaffected parents in CHD-associated genes. We identified three genome-wide significant S-CHD disorders caused by DNMs in CHD4, CDK13 and PRKD1. Our study finds evidence for distinct genetic architectures underlying the low sibling recurrence risk in S-CHD and NS-CHD.

Inferring Crohn's disease association from exome sequences by integrating biological knowledge

Background

Exome sequencing has been emerged as a primary method to identify detailed sequence variants associated with complex diseases including Crohn’s disease in the protein-coding regions of human genome. However, constructing an interpretable model for exome sequencing data is challenging because of the huge diversity of genomic variation. In addition, it has been known that utilizing biologically relevant information in a rigorous manner is essential for effectively extracting disease-associated information.

Results

In this paper, we incorporate three different types of biological knowledge such as predicted pathogenicity, disease gene annotation, and functional interaction network of human genes, and integrate them with exome sequence data in non-negative matrix tri-factorization framework. Based on the proposed method, we successfully identified Crohn’s disease patients from exome sequencing data and achieved the area under the receiver operating characteristics curve (AUC) of 0.816, while other clustering methods not using biological information achieved the AUC of 0.786. Moreover, the disease association score derived from our method showed higher correlation with Crohn’s disease genes than other unrelated genes.

Conclusions

As a consequence, by integrating biological information across multiple levels such as variant, gene, and systems, our method could be useful for identifying disease susceptibility and its associated genes from exome sequencing data.

Using proteomic and genomic methods to understand JDM

Juvenile dermatomyositis is a complex illness characterized by vascular/perivascular inflammation, primarily in the skin and muscles. In this review, we discuss how proteomic and genomic technologies have expanded our understanding of the immune pathogenesis of this disease. We will also discuss further directions that the field may take to use existing and developing technologies to further our understanding of this often-perplexing disease.

Variation in PTCHD2, CRISP3, NAP1L4, FSCB, and AP3B2 associated with spherical equivalent

PURPOSE

Ocular refraction is measured in spherical equivalent as the power of the external lens required to focus images on the retina. Myopia (nearsightedness) and hyperopia (farsightedness) are the most common refractive errors, and the leading causes of visual impairment and blindness in the world. The goal of this study is to identify rare and low-frequency variants that influence spherical equivalent.

METHODS

We conducted variant-level and gene-level quantitative trait association analyses for mean spherical equivalent, using data from 1,560 individuals in the Beaver Dam Eye Study. Genotyping was conducted using the Illumina exome array. We analyzed 34,976 single nucleotide variants and 11,571 autosomal genes across the genome, using single-variant tests as well as gene-based tests.

RESULTS

Spherical equivalent was significantly associated with five genes in gene-based analysis: PTCHD2 at 1p36.22 (p = 3.6 × 10(-7)), CRISP3 at 6p12.3 (p = 4.3 × 10(-6)), NAP1L4 at 11p15.5 (p = 3.6 × 10(-6)), FSCB at 14q21.2 (p = 1.5 × 10(-7)), and AP3B2 at 15q25.2 (p = 1.6 × 10(-7)). The variant-based tests identified evidence suggestive of association with two novel variants in linkage disequilibrium (pairwise r(2) = 0.80) in the TCTE1 gene region at 6p21.1 (rs2297336, minor allele frequency (MAF) = 14.1%, β = -0.62 p = 3.7 × 10(-6); rs324146, MAF = 16.9%, β = -0.55, p = 1.4 × 10(-5)). In addition to these novel findings, we successfully replicated a previously reported association with rs634990 near GJD2 at 15q14 (MAF = 47%, β = -0.29, p=1.8 × 10(-3)). We also found evidence of association with spherical equivalent on 2q37.1 in PRSS56 at rs1550094 (MAF = 31%, β = -0.33, p = 1.7 × 10(-3)), a region previously associated with myopia.

CONCLUSIONS

We identified several novel candidate genes that may play a role in the control of spherical equivalent. However, further studies are needed to replicate these findings. In addition, our results contribute to the increasing evidence that variation in the GJD2 and PRSS56 genes influence the development of refractive errors. Identifying that variation in these genes is associated with spherical equivalent may provide further insight into the etiology of myopia and consequent vision loss.

Identification of genetic variants associated with susceptibility to West Nile virus neuroinvasive disease

West Nile virus (WNV) infection results in a diverse spectrum of outcomes, and host genetics are likely to influence susceptibility to neuroinvasive disease (West Nile neuroinvasive disease (WNND)). We performed whole-exome sequencing of 44 individuals with WNND and identified alleles associated with severe disease by variant filtration in cases, kernel association testing in cases and controls and single-nucleotide polymorphism (SNP) imputation into a larger cohort of WNND cases and seropositive controls followed by genome-wide association analysis. Variant filtration prioritized genes based on the enrichment of otherwise rare variants, but did not unambiguously implicate variants shared by a majority of cases. Kernel association demonstrated enrichment for risk and protective alleles in the human leukocyte antigen (HLA)-A and HLA-DQB1 loci that have well understood roles in antiviral immunity. Two loci, HERC5 and an intergenic region between CD83 and JARID2, were implicated by multiple imputed SNPs and exceeded genome-wide significance in a discovery cohort (n=862). SNPs at two additional loci, TFCP2L1 and CACNA1H, achieved genome-wide significance after association testing of directly genotyped and imputed SNPs in a discovery cohort (n=862) and a separate replication cohort (n=1387). The context of these loci suggests that immunoregulatory, ion channel and endothelial barrier functions may be important elements of the host response to WNV.

A Review of Whole-Exome Sequencing Efforts Toward Hereditary Breast Cancer Susceptibility Gene Discovery

Inherited genetic risk factors contribute toward breast cancer (BC) onset. BC risk variants can be divided into three categories of penetrance (high, moderate, and low) that reflect the probability of developing the disease. Traditional BC susceptibility gene discovery approaches that searched for high- and moderate-risk variants in familial BC cases have had limited success; to date, these risk variants explain only ∼30% of familial BC cases. Next-generation sequencing technologies can be used to search for novel high and moderate BC risk variants, and this manuscript reviews 12 familial BC whole-exome sequencing efforts. Study design, filtering strategies, and segregation and validation analyses are discussed. Overall, only a modest number of novel BC risk genes were identified, and 90% and 97% of the exome-sequenced families and cases, respectively, had no BC risk variants reported. It is important to learn from these studies and consider alternate strategies in order to make further advances. The discovery of new BC susceptibility genes is critical for improved risk assessment and to provide insight toward disease mechanisms for the development of more effective therapies.

Pharmacogenetics of Drug-Induced QT Interval Prolongation: An Update

A prolonged QT interval is an important risk factor for ventricular arrhythmias and sudden cardiac death. QT prolongation can be caused by drugs. There are multiple risk factors for drug-induced QT prolongation, including genetic variation. QT prolongation is one of the most common reasons for withdrawal of drugs from the market, despite the fact that these drugs may be beneficial for certain patients and not harmful in every patient. Identifying genetic variants associated with drug-induced QT prolongation might add to tailored pharmacotherapy and prevent beneficial drugs from being withdrawn unnecessarily. In this review, our objective was to provide an overview of the genetic background of drug-induced QT prolongation, distinguishing pharmacokinetic and pharmacodynamic pathways. Pharmacokinetic-mediated genetic susceptibility is mainly characterized by variation in genes encoding drug-metabolizing cytochrome P450 enzymes or drug transporters. For instance, the P-glycoprotein drug transporter plays a role in the pharmacokinetic susceptibility of drug-induced QT prolongation. The pharmacodynamic component of genetic susceptibility is mainly characterized by genes known to be associated with QT interval duration in the general population and genes in which the causal mutations of congenital long QT syndromes are located. Ethnicity influences susceptibility to drug-induced QT interval prolongation, with Caucasians being more sensitive than other ethnicities. Research on the association between pharmacogenetic interactions and clinical endpoints such as sudden cardiac death is still limited. Future studies in this area could enable us to determine the risk of arrhythmias more adequately in clinical practice.

Rare variants in known and novel candidate genes predisposing to statin-associated myopathy

AIM

Genetic variants affecting statin uptake, metabolism or predisposing to muscular diseases may confer susceptibility to statin-induced myopathy. Besides the SLCO1B1 rs4149056 genotype, common genetic variants do not seem to determine statin-associated myopathy. Here we aimed to address the potential role of rare variants.

METHODS

We performed whole exome sequencing in 88 individuals suffering from statin-associated myopathy and assessed the burden of rare variants using candidate-gene and exome-wide association analysis.

RESULTS

In the novel candidate gene CLCN1, we identified a heterozygote truncating mutation p.R894* in four patients. In addition, we detected predictably pathogenic case-specific variants in MYOT, CYP3A5, SH3TC2, FBXO32 and RBM20.

CONCLUSION

These findings support the role of rare variants and nominate loci for follow-up studies.

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.

Discovery of rare variants for complex phenotypes

With the rise of sequencing technologies, it is now feasible to assess the role rare variants play in the genetic contribution to complex trait variation. While some of the earlier targeted sequencing studies successfully identified rare variants of large effect, unbiased gene discovery using exome sequencing has experienced limited success for complex traits. Nevertheless, rare variant association studies have demonstrated that rare variants do contribute to phenotypic variability, but sample sizes will likely have to be even larger than those of common variant association studies to be powered for the detection of genes and loci. Large-scale sequencing efforts of tens of thousands of individuals, such as the UK10K Project and aggregation efforts such as the Exome Aggregation Consortium, have made great strides in advancing our knowledge of the landscape of rare variation, but there remain many considerations when studying rare variation in the context of complex traits. We discuss these considerations in this review, presenting a broad range of topics at a high level as an introduction to rare variant analysis in complex traits including the issues of power, study design, sample ascertainment, de novo variation, and statistical testing approaches. Ultimately, as sequencing costs continue to decline, larger sequencing studies will yield clearer insights into the biological consequence of rare mutations and may reveal which genes play a role in the etiology of complex traits.

Global Carrier Rates of Rare Inherited Disorders Using Population Exome Sequences

Exome sequencing has revealed the causative mutations behind numerous rare, inherited disorders, but it is challenging to find reliable epidemiological values for rare disorders. Here, I provide a genetic epidemiology method to identify the causative mutations behind rare, inherited disorders using two population exome sequences (1000 Genomes and NHLBI). I created global maps of carrier rate distribution for 18 recessive disorders in 16 diverse ethnic populations. Out of a total of 161 mutations associated with 18 recessive disorders, I detected 24 mutations in either or both exome studies. The genetic mapping revealed strong international spatial heterogeneities in the carrier patterns of the inherited disorders. I next validated this methodology by statistically evaluating the carrier rate of one well-understood disorder, sickle cell anemia (SCA). The population exome-based epidemiology of SCA [African (allele frequency (AF) = 0.0454, N = 2447), Asian (AF = 0, N = 286), European (AF = 0.000214, N = 4677), and Hispanic (AF = 0.0111, N = 362)] was not significantly different from that obtained from a clinical prevalence survey. A pair-wise proportion test revealed no significant differences between the two exome projects in terms of AF (46/48 cases; P > 0.05). I conclude that population exome-based carrier rates can form the foundation for a prospectively maintained database of use to clinical geneticists. Similar modeling methods can be applied to many inherited disorders.

Emerging Genomic Tools for Legume Breeding: Current Status and Future Prospects

Legumes play a vital role in ensuring global nutritional food security and improving soil quality through nitrogen fixation. Accelerated higher genetic gains is required to meet the demand of ever increasing global population. In recent years, speedy developments have been witnessed in legume genomics due to advancements in next-generation sequencing (NGS) and high-throughput genotyping technologies. Reference genome sequences for many legume crops have been reported in the last 5 years. The availability of the draft genome sequences and re-sequencing of elite genotypes for several important legume crops have made it possible to identify structural variations at large scale. Availability of large-scale genomic resources and low-cost and high-throughput genotyping technologies are enhancing the efficiency and resolution of genetic mapping and marker-trait association studies. Most importantly, deployment of molecular breeding approaches has resulted in development of improved lines in some legume crops such as chickpea and groundnut. In order to support genomics-driven crop improvement at a fast pace, the deployment of breeder-friendly genomics and decision support tools seems appear to be critical in breeding programs in developing countries. This review provides an overview of emerging genomics and informatics tools/approaches that will be the key driving force for accelerating genomics-assisted breeding and ultimately ensuring nutritional and food security in developing countries.

Amplicon-based semiconductor sequencing of human exomes: performance evaluation and optimization strategies

The Ion Proton platform allows to perform whole exome sequencing (WES) at low cost, providing rapid turnaround time and great flexibility. Products for WES on Ion Proton system include the AmpliSeq Exome kit and the recently introduced HiQ sequencing chemistry. Here, we used gold standard variants from GIAB consortium to assess the performances in variants identification, characterize the erroneous calls and develop a filtering strategy to reduce false positives. The AmpliSeq Exome kit captures a large fraction of bases (>94 %) in human CDS, ClinVar genes and ACMG genes, but with 2,041 (7 %), 449 (13 %) and 11 (19 %) genes not fully represented, respectively. Overall, 515 protein coding genes contain hard-to-sequence regions, including 90 genes from ClinVar. Performance in variants detection was maximum at mean coverage >120×, while at 90× and 70× we measured a loss of variants of 3.2 and 4.5 %, respectively. WES using HiQ chemistry showed ~71/97.5 % sensitivity, ~37/2 % FDR and ~0.66/0.98 F1 score for indels and SNPs, respectively. The proposed low, medium or high-stringency filters reduced the amount of false positives by 10.2, 21.2 and 40.4 % for indels and 21.2, 41.9 and 68.2 % for SNP, respectively. Amplicon-based WES on Ion Proton platform using HiQ chemistry emerged as a competitive approach, with improved accuracy in variants identification. False-positive variants remain an issue for the Ion Torrent technology, but our filtering strategy can be applied to reduce erroneous variants.

Mutations and altered expression of SERPINF1 in patients with familial otosclerosis

Otosclerosis is a relatively common heterogenous condition, characterized by abnormal bone remodelling in the otic capsule leading to fixation of the stapedial footplate and an associated conductive hearing loss. Although familial linkage and candidate gene association studies have been performed in recent years, little progress has been made in identifying disease-causing genes. Here, we used whole-exome sequencing in four families exhibiting dominantly inherited otosclerosis to identify 23 candidate variants (reduced to 9 after segregation analysis) for further investigation in a secondary cohort of 84 familial cases. Multiple mutations were found in the SERPINF1 (Serpin Peptidase Inhibitor, Clade F) gene which encodes PEDF (pigment epithelium-derived factor), a potent inhibitor of angiogenesis and known regulator of bone density. Six rare heterozygous SERPINF1 variants were found in seven patients in our familial otosclerosis cohort; three are missense mutations predicted to be deleterious to protein function. The other three variants are all located in the 5′-untranslated region (UTR) of an alternative spliced transcript SERPINF1-012. RNA-seq analysis demonstrated that this is the major SERPINF1 transcript in human stapes bone. Analysis of stapes from two patients with the 5′-UTR mutations showed that they had reduced expression of SERPINF1-012. All three 5′-UTR mutations are predicted to occur within transcription factor binding sites and reporter gene assays confirmed that they affect gene expression levels. Furthermore, RT-qPCR analysis of stapes bone cDNA showed that SERPINF1-012 expression is reduced in otosclerosis patients with and without SERPINF1 mutations, suggesting that it may be a common pathogenic pathway in the disease.

Society for Pediatric Research 2015 Young Investigator Award: genetics of human hematopoiesis-what patients can teach us about blood cell production

Blood cell production or hematopoiesis is one of the most well-understood paradigms of cell differentiation in the body. The majority of work on hematopoiesis comes from studies that have primarily been conducted in mice, zebrafish, or other valuable model systems. However, it is clear that such model organisms may not consistently and faithfully mimic what is observed in humans with blood disorders. Moreover, there is significant divergence between species that is increasingly being appreciated at the genomic level. As a result, there is an opportunity to use observations in humans to provide a refined view of hematopoiesis. Here, we discuss vignettes from our work that illustrate how insight from human genetics can improve our understanding of blood cell production and identify promising therapeutic approaches for blood disorders.

Genetics of Bipolar Disorder: Recent Update and Future Directions

Although genetic studies of Bipolar Disorder have been pursued for decades, it has only been in the last several years that clearly replicated findings have emerged. These findings, typically of modest effects, point to a polygenic genetic architecture consisting of multiple common and rare susceptibility variants. While larger genome-wide association studies are ongoing, the advent of whole exome and genome sequencing should lead to the identification of rare, and potentially more penetrant, variants. Progress along both fronts will provide novel insights into the biology of Bipolar Disorder and help usher in a new era of personalized medicine and improved treatments.

Next Generation Sequencing in Alzheimer's Disease

For the first time in the history of human genetics research, it is now both technically feasible and economically affordable to screen individual genomes for novel disease-causing mutations at base-pair resolution using "next-generation sequencing" (NGS). One popular aim in many of today's NGS studies is genome resequencing (in part or whole) to identify DNA variants potentially accounting for the "missing heritability" problem observed in many genetically complex traits. Thus far, only relatively few projects have applied these powerful new technologies to search for novel Alzheimer's disease (AD) related sequence variants. In this review, I summarize the findings from the first NGS-based resequencing studies in AD and discuss their potential implications and limitations. Notable recent discoveries using NGS include the identification of rare susceptibility modifying alleles in APP, TREM2, and PLD3. Several other large-scale NGS projects are currently underway so that additional discoveries can be expected over the coming years.

Meta-Analysis of Rare Variant Association Tests in Multiethnic Populations

Several methods have been proposed to increase power in rare variant association testing by aggregating information from individual rare variants (MAF < 0.005). However, how to best combine rare variants across multiple ethnicities and the relative performance of designs using different ethnic sampling fractions remains unknown. In this study, we compare the performance of several statistical approaches for assessing rare variant associations across multiple ethnicities. We also explore how different ethnic sampling fractions perform, including single-ethnicity studies and studies that sample up to four ethnicities. We conducted simulations based on targeted sequencing data from 4,611 women in four ethnicities (African, European, Japanese American, and Latina). As with single-ethnicity studies, burden tests had greater power when all causal rare variants were deleterious, and variance component-based tests had greater power when some causal rare variants were deleterious and some were protective. Multiethnic studies had greater power than single-ethnicity studies at many loci, with inclusion of African Americans providing the largest impact. On average, studies including African Americans had as much as 20% greater power than equivalently sized studies without African Americans. This suggests that association studies between rare variants and complex disease should consider including subjects from multiple ethnicities, with preference given to genetically diverse groups.

Histone Modifier Genes Alter Conotruncal Heart Phenotypes in 22q11.2 Deletion Syndrome

We performed whole exome sequence (WES) to identify genetic modifiers on 184 individuals with 22q11.2 deletion syndrome (22q11DS), of whom 89 case subjects had severe congenital heart disease (CHD) and 95 control subjects had normal hearts. Three genes including JMJD1C (jumonji domain containing 1C), RREB1 (Ras responsive element binding protein 1), and SEC24C (SEC24 family member C) had rare (MAF < 0.001) predicted deleterious single-nucleotide variations (rdSNVs) in seven case subjects and no control subjects (p = 0.005; Fisher exact and permutation tests). Because JMJD1C and RREB1 are involved in chromatin modification, we investigated other histone modification genes. Eighteen case subjects (20%) had rdSNVs in four genes (JMJD1C, RREB1, MINA, KDM7A) all involved in demethylation of histones (H3K9, H3K27). Overall, rdSNVs were enriched in histone modifier genes that activate transcription (Fisher exact p = 0.0004, permutations, p = 0.0003, OR = 5.16); however, rdSNVs in control subjects were not enriched. This implicates histone modification genes as influencing risk for CHD in presence of the deletion.

Exome sequencing reveals a high genetic heterogeneity on familial Hirschsprung disease

Hirschsprung disease (HSCR; OMIM 142623) is a developmental disorder characterized by aganglionosis along variable lengths of the distal gastrointestinal tract, which results in intestinal obstruction. Interactions among known HSCR genes and/or unknown disease susceptibility loci lead to variable severity of phenotype. Neither linkage nor genome-wide association studies have efficiently contributed to completely dissect the genetic pathways underlying this complex genetic disorder. We have performed whole exome sequencing of 16 HSCR patients from 8 unrelated families with SOLID platform. Variants shared by affected relatives were validated by Sanger sequencing. We searched for genes recurrently mutated across families. Only variations in the FAT3 gene were significantly enriched in five families. Within-family analysis identified compound heterozygotes for AHNAK and several genes (N = 23) with heterozygous variants that co-segregated with the phenotype. Network and pathway analyses facilitated the discovery of polygenic inheritance involving FAT3, HSCR known genes and their gene partners. Altogether, our approach has facilitated the detection of more than one damaging variant in biologically plausible genes that could jointly contribute to the phenotype. Our data may contribute to the understanding of the complex interactions that occur during enteric nervous system development and the etiopathology of familial HSCR.

Prioritization Of Nonsynonymous Single Nucleotide Variants For Exome Sequencing Studies Via Integrative Learning On Multiple Genomic Data

The rapid advancement of next generation sequencing technology has greatly accelerated the progress for understanding human inherited diseases via such innovations as exome sequencing. Nevertheless, the identification of causative variants from sequencing data remains a great challenge. Traditional statistical genetics approaches such as linkage analysis and association studies have limited power in analyzing exome sequencing data, while relying on simply filtration strategies and predicted functional implications of mutations to pinpoint pathogenic variants are prone to produce false positives. To overcome these limitations, we herein propose a supervised learning approach, termed snvForest, to prioritize candidate nonsynonymous single nucleotide variants for a specific type of disease by integrating 11 functional scores at the variant level and 8 association scores at the gene level. We conduct a series of large-scale in silico validation experiments, demonstrating the effectiveness of snvForest across 2,511 diseases of different inheritance styles and the superiority of our approach over two state-of-the-art methods. We further apply snvForest to three real exome sequencing data sets of epileptic encephalophathies and intellectual disability to show the ability of our approach to identify causative de novo mutations for these complex diseases. The online service and standalone software of snvForest are found at http://bioinfo.au.tsinghua.edu.cn/jianglab/snvforest.

The dichotomy between disease phenotype databases and the implications for understanding complex diseases involving the major histocompatibility complex

Many genes related to innate and adaptive immunity reside within the major histocompatibility complex (MHC) and have been associated with a multitude of complex, immune-related disorders. Despite years of genetic study, this region has seen few causative determinants discovered for immune-mediated diseases. Reported associations have been curated in various databases including the Genetic Association Database, NCBI database of clinically relevant variants (ClinVar) and the Human Gene Mutation Database and together capture genetic associations and annotated pathogenic loci within the MHC and across the genome for a variety of complex, immune-mediated diseases. A review of these three distinct databases reveals disparate annotations between associated genes and pathogenic loci, alluding to the polygenic, multifactorial nature of immune-mediated diseases and the pleiotropic character of genes within the MHC. The technical limitations and inherent biases imposed by current approaches and technologies in studying the MHC create a strong case for the need to perform targeted deep sequencing of the MHC and other immunologically relevant loci in order to fully elucidate and study the causative elements of complex immune-mediated diseases.

Regulatory domain or CpG site variation in SLC12A5, encoding the chloride transporter KCC2, in human autism and schizophrenia

Many encoded gene products responsible for neurodevelopmental disorders (NDs) like autism spectrum disorders (ASD), schizophrenia (SCZ), intellectual disability (ID), and idiopathic generalized epilepsy (IGE) converge on networks controlling synaptic function. An increase in KCC2 (SLC12A5) Cl⁻ transporter activity drives the developmental GABA excitatory-inhibitory sequence, but the role of KCC2 in human NDs is essentially unknown. Here, we report two rare, non-synonymous (NS), functionally-impairing variants in the KCC2 C-terminal regulatory domain (CTRD) in human ASD (R952H and R1049C) and SCZ (R952H) previously linked with IGE and familial febrile seizures, and another novel NS KCC2 variant in ASD (R1048W) with highly-predicted pathogenicity. Exome data from 2517 simplex families in the ASD Simon Simplex Collection (SSC) revealed significantly more KCC2 CTRD variants in ASD cases than controls, and interestingly, these were more often synonymous and predicted to disrupt or introduce a CpG site. Furthermore, full gene analysis showed ASD cases are more likely to contain rare KCC2 variants affecting CpG sites than controls. These data suggest genetically-encoded dysregulation of KCC2-dependent GABA signaling may contribute to multiple human NDs.

Connectivity Homology Enables Inter-Species Network Models of Synthetic Lethality

Synthetic lethality is a genetic interaction wherein two otherwise nonessential genes cause cellular inviability when knocked out simultaneously. Drugs can mimic genetic knock-out effects; therefore, our understanding of promiscuous drugs, polypharmacology-related adverse drug reactions, and multi-drug therapies, especially cancer combination therapy, may be informed by a deeper understanding of synthetic lethality. However, the colossal experimental burden in humans necessitates in silico methods to guide the identification of synthetic lethal pairs. Here, we present SINaTRA (Species-INdependent TRAnslation), a network-based methodology that discovers genome-wide synthetic lethality in translation between species. SINaTRA uses connectivity homology, defined as biological connectivity patterns that persist across species, to identify synthetic lethal pairs. Importantly, our approach does not rely on genetic homology or structural and functional similarity, and it significantly outperforms models utilizing these data. We validate SINaTRA by predicting synthetic lethality in S. pombe using S. cerevisiae data, then identify over one million putative human synthetic lethal pairs to guide experimental approaches. We highlight the translational applications of our algorithm for drug discovery by identifying clusters of genes significantly enriched for single- and multi-drug cancer therapies.

Synthetic lethality is a genetic interaction that has promising implications for informing novel cancer therapies. Over 200 million pairwise tests would be required to identify all pairwise synthetic lethal interactions in humans–currently, an impossibly large experimental burden. To simplify the process, we have developed a method to predict human synthetic lethal pairs in translation from a well-studied species to one in which synthetic lethality is understudied using both species’ protein-protein interaction networks. Here, we explore the model’s success in translation from S. cerevisiae to S. pombe. We then predict human synthetic lethality and suggest novel areas of inquiry for cancer therapies.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Likelihood-based complex trait association testing for arbitrary depth sequencing data

In next generation sequencing (NGS)-based genetic studies, researchers typically perform genotype calling first and then apply standard genotype-based methods for association testing. However, such a two-step approach ignores genotype calling uncertainty in the association testing step and may incur power loss and/or inflated type-I error. In the recent literature, a few robust and efficient likelihood based methods including both likelihood ratio test (LRT) and score test have been proposed to carry out association testing without intermediate genotype calling. These methods take genotype calling uncertainty into account by directly incorporating genotype likelihood function (GLF) of NGS data into association analysis. However, existing LRT methods are computationally demanding or do not allow covariate adjustment; while existing score tests are not applicable to markers with low minor allele frequency (MAF). We provide an LRT allowing flexible covariate adjustment, develop a statistically more powerful score test and propose a combination strategy (UNC combo) to leverage the advantages of both tests. We have carried out extensive simulations to evaluate the performance of our proposed LRT and score test. Simulations and real data analysis demonstrate the advantages of our proposed combination strategy: it offers a satisfactory trade-off in terms of computational efficiency, applicability (accommodating both common variants and variants with low MAF) and statistical power, particularly for the analysis of quantitative trait where the power gain can be up to ∼60% when the causal variant is of low frequency (MAF < 0.01).

AVAILABILITY AND IMPLEMENTATION

UNC combo and the associated R files, including documentation, examples, are available at http://www.unc.edu/∼yunmli/UNCcombo/

CONTACT

yunli@med.unc.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

An efficient resampling method for calibrating single and gene-based rare variant association analysis in case-control studies

For aggregation tests of genes or regions, the set of included variants often have small total minor allele counts (MACs), and this is particularly true when the most deleterious sets of variants are considered. When MAC is low, commonly used asymptotic tests are not well calibrated for binary phenotypes and can have conservative or anti-conservative results and potential power loss. Empirical p-values obtained via resampling methods are computationally costly for highly significant p-values and the results can be conservative due to the discrete nature of resampling tests. Based on the observation that only the individuals containing minor alleles contribute to the score statistics, we develop an efficient resampling method for single and multiple variant score-based tests that can adjust for covariates. Our method can improve computational efficiency >1000-fold over conventional resampling for low MAC variant sets. We ameliorate the conservativeness of results through the use of mid-p-values. Using the estimated minimum achievable p-value for each test, we calibrate QQ plots and provide an effective number of tests. In analysis of a case-control study with deep exome sequence, we demonstrate that our methods are both well calibrated and also reduce computation time significantly compared with resampling methods.

Phenotypic extremes in rare variant study designs

Currently, next-generation sequencing studies aim to identify rare and low-frequency variation that may contribute to disease. For a given effect size, as the allele frequency decreases, the power to detect genes or variants of interest also decreases. Although many methods have been proposed for the analysis of such data, study design and analytic issues still persist in data interpretation. In this study we present sequencing data for ABCA1 that has known rare variants associated with high-density lipoprotein cholesterol (HDL-C). We contrast empirical findings from two study designs: a phenotypic extreme sample and a population-based random sample. We found differing strengths of association with HDL-C across the two study designs (P=0.0006 with n=701 phenotypic extremes vs P=0.03 with n=1600 randomly sampled individuals). To explore this apparent difference in evidence for association, we performed a simulation study focused on the impact of phenotypic selection on power. We demonstrate that the power gain for an extreme phenotypic selection study design is much greater in rare variant studies than for studies of common variants. Our study confirms that studying phenotypic extremes is critical in rare variant studies because it boosts power in two ways: the typical increases from extreme sampling and increasing the proportion of relevant functional variants ascertained and thereby tested for association. Furthermore, we show that when combining statistical evidence through meta-analysis from an extreme-selected sample and a second separate population-based random sample, power is lower when a traditional sample size weighting is used compared with weighting by the noncentrality parameter.

Challenges in exome analysis by LifeScope and its alternative computational pipelines

Background

Every next generation sequencing (NGS) platform relies on proprietary and open source computational tools to analyze sequencing data. NGS tools for Illumina platforms are well documented which is not the case with AB SOLiD systems. We applied several computational and variant calling pipelines to analyse targeted exome sequencing data obtained using AB SOLiD 5500 system. Our investigated tools comprised proprietary LifeScope’s pipeline in combination with open source color-space competent mapping programs and a variant caller. We present instrumental details of the pipelines that were used and quantitative comparative analysis of variant lists generated by LifeScope’s pipeline versus open source tools.

Results

Sufficient coverage of targeted regions was achieved by all investigated pipelines. High variability was observed in identities of variants across the mapping programs. We observed less than 50 % concordance of variant lists produced by approaches based on different mapping algorithms. We summarized different approaches with regards to coverage (DP) and quality (QUAL) properties of the variants provided by GATK and found that LifeScope’s computational pipeline is superior. Fusion of information on mapping profiles (pileup) at genomic positions of variants in several different alignments proved to be a useful strategy to assess questionable singleton variants.

Conclusions

We quantitatively supported a conclusion that Lifescope’s pipeline is superior for processing sequencing data obtained by AB SOLiD 5500 system. Nevertheless the use of alternative pipelines is encouraged because aggregation of information from other mapping and variant calling approaches helps to resolve questionable calls and increases the confidence of the call. It was noted that a coverage threshold for variant to be considered for further analysis has to be chosen in data-driven way to prevent a loss of important information.

Common and Rare Variant Association Study for Plasma Lipids and Coronary Artery Disease

Blood lipid levels are highly heritable and modifiable risk factors for coronary artery disease (CAD), and are the leading cause of death worldwide. These facts have motivated human genetic association studies that have the substantial potential to define the risk factors that are causal and to identify pathways and therapeutic targets for lipids and CAD.The success of the HapMap project that provided an extensive catalog of human genetic variations and the development of microarray based genotyping chips (typically containing variations with allele frequencies ＞ 5%) facilitated common variant association study (CVAS; formerly termed genome-wide association study, GWAS) identifying disease-associated variants in a genome-wide manner. To date, 157 loci associated with blood lipids and 46 loci with CAD have been successfully identified, accounting for approximately 12%-14% of heritability for lipids and 10% of heritability for CAD. However, there is yet a major challenge termed "missing heritability problem," namely the observation that loci detected by CVAS explain only a small fraction of the inferred genetic variations. To explain such missing portions, focuses in genetic association studies have shifted from common to rare variants. However, it is challenging to apply rare variant association study (RVAS) in an unbiased manner because such variants typically lack the sufficient number to be identified statistically.In this review, we provide a current understanding of the genetic architecture mostly derived from CVAS, and several updates on the progress and limitations of RVAS for lipids and CAD.

Aetiological coding sequence variants in non-syndromic premature ovarian failure: From genetic linkage analysis to next generation sequencing

Premature ovarian failure (POF) is a frequent pathology affecting 1-1.5% of women under 40 years old. Despite advances in diagnosing and treating human infertility, POF is still classified as being idiopathic in 50-80% of cases, strongly suggesting a genetic origin for the disease. Different types of autosomal and X-linked genetic anomalies can originate the phenotype in syndromic and non-syndromic POF cases. Particular interest has been focused on research into non-syndromic POF causative coding variants during the past two decades. This has been based on the assumption that amino acid substitutions might modify the intrinsic physicochemical properties of functional proteins, thereby inducing pathological phenotypes. In this case, a restricted number of mutations might originate the disease. However, like other complex pathologies, POF might result from synergistic/compensatory effects caused by several low-to-mildly drastic mutations which have frequently been classified as non-functional SNPs. Indeed, reproductive phenotypes can be considered as quantitative traits resulting from the subtle interaction of many genes. Although numerous sequencing projects have involved candidate genes, only a few coding mutations explaining a low percentage of cases have been described. Such apparent failure to identify aetiological coding sequence variations might have been due to the inherent molecular complexity of mammalian reproduction and to the difficulty of simultaneously analysing large genomic regions by Sanger sequencing. The purpose of this review is to present the molecular and cellular effects caused by non-synonymous mutations which have been formally associated, by functional tests, with the aetiology of hypergonadotropic non-syndromic POF. Considerations have also been included regarding the polygenic nature of reproduction and POF, as well as future approaches for identifying novel aetiological genes based on next generation sequencing (NGS).

Performance comparison of four commercial human whole-exome capture platforms

Whole exome sequencing (WXS) is widely used to identify causative genetic mutations of diseases. However, not only have several commercial human exome capture platforms been developed, but substantial updates have been released in the past few years. We report a performance comparison for the latest release of four commercial platforms, Roche/NimbleGen’s SeqCap EZ Human Exome Library v3.0, Illumina’s Nextera Rapid Capture Exome (v1.2), Agilent’s SureSelect XT Human All Exon v5 and Agilent’s SureSelect QXT, using the same DNA samples. Agilent XT showed the highest target enrichment efficiency and the best SNV and short indel detection sensitivity in coding regions with the least amount of sequencing. Agilent QXT had slightly inferior target enrichment than Agilent XT. Illumina, with additional sequencing, detected SNVs and short indels at the same quality as Agilent XT, and showed the best performance in coverage of medically interesting mutations. NimbleGen detected more SNVs and indels in untranslated regions than the others. We also found that the platforms, which enzymatically fragment the genomic DNA (gDNA), detected more homozygous SNVs than those using sonicated gDNA. We believe that our analysis will help investigators when selecting a suitable exome capture platform for their particular research.

Linkage Analysis of Urine Arsenic Species Patterns in the Strong Heart Family Study

Arsenic toxicokinetics are important for disease risks in exposed populations, but genetic determinants are not fully understood. We examined urine arsenic species patterns measured by HPLC-ICPMS among 2189 Strong Heart Study participants 18 years of age and older with data on ~400 genome-wide microsatellite markers spaced ~10 cM and arsenic speciation (683 participants from Arizona, 684 from Oklahoma, and 822 from North and South Dakota). We logit-transformed % arsenic species (% inorganic arsenic, %MMA, and %DMA) and also conducted principal component analyses of the logit % arsenic species. We used inverse-normalized residuals from multivariable-adjusted polygenic heritability analysis for multipoint variance components linkage analysis. We also examined the contribution of polymorphisms in the arsenic metabolism gene AS3MT via conditional linkage analysis. We localized a quantitative trait locus (QTL) on chromosome 10 (LOD 4.12 for %MMA, 4.65 for %DMA, and 4.84 for the first principal component of logit % arsenic species). This peak was partially but not fully explained by measured AS3MT variants. We also localized a QTL for the second principal component of logit % arsenic species on chromosome 5 (LOD 4.21) that was not evident from considering % arsenic species individually. Some other loci were suggestive or significant for 1 geographical area but not overall across all areas, indicating possible locus heterogeneity. This genome-wide linkage scan suggests genetic determinants of arsenic toxicokinetics to be identified by future fine-mapping, and illustrates the utility of principal component analysis as a novel approach that considers % arsenic species jointly.

Genome-wide Polygenic Burden of Rare Deleterious Variants in Sudden Unexpected Death in Epilepsy

Sudden unexpected death in epilepsy (SUDEP) represents the most severe degree of the spectrum of epilepsy severity and is the commonest cause of epilepsy-related premature mortality. The precise pathophysiology and the genetic architecture of SUDEP remain elusive. Aiming to elucidate the genetic basis of SUDEP, we analysed rare, protein-changing variants from whole-exome sequences of 18 people who died of SUDEP, 87 living people with epilepsy and 1479 non-epilepsy disease controls. Association analysis revealed a significantly increased genome-wide polygenic burden per individual in the SUDEP cohort when compared to epilepsy (P = 5.7 × 10(- 3)) and non-epilepsy disease controls (P = 1.2 × 10(- 3)). The polygenic burden was driven both by the number of variants per individual, and over-representation of variants likely to be deleterious in the SUDEP cohort. As determined by this study, more than a thousand genes contribute to the observed polygenic burden within the framework of this study. Subsequent gene-based association analysis revealed five possible candidate genes significantly associated with SUDEP or epilepsy, but no one single gene emerges as common to the SUDEP cases. Our findings provide further evidence for a genetic susceptibility to SUDEP, and suggest an extensive polygenic contribution to SUDEP causation. Thus, an overall increased burden of deleterious variants in a highly polygenic background might be important in rendering a given individual more susceptible to SUDEP. Our findings suggest that exome sequencing in people with epilepsy might eventually contribute to generating SUDEP risk estimates, promoting stratified medicine in epilepsy, with the eventual aim of reducing an individual patient's risk of SUDEP.

The Use of Non-Variant Sites to Improve the Clinical Assessment of Whole-Genome Sequence Data

Genetic testing, which is now a routine part of clinical practice and disease management protocols, is often based on the assessment of small panels of variants or genes. On the other hand, continuous improvements in the speed and per-base costs of sequencing have now made whole exome sequencing (WES) and whole genome sequencing (WGS) viable strategies for targeted or complete genetic analysis, respectively. Standard WGS/WES data analytical workflows generally rely on calling of sequence variants respect to the reference genome sequence. However, the reference genome sequence contains a large number of sites represented by rare alleles, by known pathogenic alleles and by alleles strongly associated to disease by GWAS. It’s thus critical, for clinical applications of WGS and WES, to interpret whether non-variant sites are homozygous for the reference allele or if the corresponding genotype cannot be reliably called. Here we show that an alternative analytical approach based on the analysis of both variant and non-variant sites from WGS data allows to genotype more than 92% of sites corresponding to known SNPs compared to 6% genotyped by standard variant analysis. These include homozygous reference sites of clinical interest, thus leading to a broad and comprehensive characterization of variation necessary to an accurate evaluation of disease risk. Altogether, our findings indicate that characterization of both variant and non-variant clinically informative sites in the genome is necessary to allow an accurate clinical assessment of a personal genome. Finally, we propose a highly efficient extended VCF (eVCF) file format which allows to store genotype calls for sites of clinical interest while remaining compatible with current variant interpretation software.

Investigation of the role of TCF4 rare sequence variants in schizophrenia

Transcription factor 4 (TCF4) is one of the most robust of all reported schizophrenia risk loci and is supported by several genetic and functional lines of evidence. While numerous studies have implicated common genetic variation at TCF4 in schizophrenia risk, the role of rare, small-sized variants at this locus-such as single nucleotide variants and short indels which are below the resolution of chip-based arrays requires further exploration. The aim of the present study was to investigate the association between rare TCF4 sequence variants and schizophrenia. Exon-targeted resequencing was performed in 190 German schizophrenia patients. Six rare variants at the coding exons and flanking sequences of the TCF4 gene were identified, including two missense variants and one splice site variant. These six variants were then pooled with nine additional rare variants identified in 379 European participants of the 1000 Genomes Project, and all 15 variants were genotyped in an independent German sample (n = 1,808 patients; n = 2,261 controls). These data were then analyzed using six statistical methods developed for the association analysis of rare variants. No significant association (P < 0.05) was found. However, the results from our association and power analyses suggest that further research into the possible involvement of rare TCF4 sequence variants in schizophrenia risk is warranted by the assessment of larger cohorts with higher statistical power to identify rare variant associations.

Exome Analyses of Long QT Syndrome Reveal Candidate Pathogenic Mutations in Calmodulin-Interacting Genes

Long QT syndrome (LQTS) is an arrhythmogenic disorder that can lead to sudden death. To date, mutations in 15 LQTS-susceptibility genes have been implicated. However, the genetic cause for approximately 20% of LQTS patients remains elusive. Here, we performed whole-exome sequencing analyses on 59 LQTS and 61 unaffected individuals in 35 families and 138 unrelated LQTS cases, after genetic screening of known LQTS genes. Our systematic analysis of familial cases and subsequent verification by Sanger sequencing identified 92 candidate mutations in 88 genes for 23 of the 35 families (65.7%): these included eleven de novo, five recessive (two homozygous and three compound heterozygous) and seventy-three dominant mutations. Although no novel commonly mutated gene was identified other than known LQTS genes, protein-protein interaction (PPI) network analyses revealed ten new pathogenic candidates that directly or indirectly interact with proteins encoded by known LQTS genes. Furthermore, candidate gene based association studies using an independent set of 138 unrelated LQTS cases and 587 controls identified an additional novel candidate. Together, mutations in these new candidates and known genes explained 37.1% of the LQTS families (13 in 35). Moreover, half of the newly identified candidates directly interact with calmodulin (5 in 11; comparison with all genes; p=0.042). Subsequent variant analysis in the independent set of 138 cases identified 16 variants in the 11 genes, of which 14 were in calmodulin-interacting genes (87.5%). These results suggest an important role of calmodulin and its interacting proteins in the pathogenesis of LQTS.

The impact of host genetic variation on infection with HIV-1

The outcome after infection with the human immunodeficiency virus type 1 (HIV-1) is a complex phenotype determined by interactions among the pathogen, the human host and the surrounding environment. An impact of host genetic variation on HIV-1 susceptibility was identified early in the pandemic, with a major role attributed to the genes encoding class I human leukocyte antigens (HLA) and the chemokine receptor CCR5. Studies using genome-wide data sets have underscored the strength of these associations relative to variants located throughout the rest of the genome. However, the extent to which additional polymorphisms influence HIV-1 disease progression, and how much of the variability in outcome can be attributed to host genetics, remain largely unclear. Here we discuss findings concerning the functional impact of associated variants, outline methods for quantifying the host genetic component and examine how available genome-wide data sets may be leveraged to discover gene variants that affect the outcome of HIV-1 infection.

Clinical massively parallel sequencing for the diagnosis of myopathies

Massively parallel sequencing, otherwise known as high-throughput or next-generation sequencing, is rapidly gaining wide use in clinical practice due to possibility of simultaneous exploration of multiple genomic regions. More than 300 genes have been implicated in neuromuscular disorders, meaning that many genes need to be considered in a differential diagnosis for a patient affected with myopathy. By providing sequencing information for numerous genes at the same time, massively parallel sequencing greatly accelerates the diagnostic processes of myopathies compared to the classical "gene-after-gene" approach by Sanger sequencing. In this review, we describe multiple advantages of this powerful sequencing method for applications in myopathy diagnosis. We also outline recent studies that used this approach to discover new myopathy-causing genes and to diagnose cohorts of patients with muscular disorders. Finally, we highlight the key aspects and limitations of massively parallel sequencing that a neurologist considering this test needs to know in order to interpret the results of the test and to deal with other issues concerning the test.

Next generation sequencing in women affected by nonsyndromic premature ovarian failure displays new potential causative genes and mutations

OBJECTIVE

To identify new molecular actors involved in nonsyndromic premature ovarian failure (POF) etiology.

DESIGN

This is a retrospective case-control cohort study.

SETTING

University research group and IVF medical center.

PATIENT(S)

Twelve women affected by nonsyndromic POF. The control group included 176 women whose menopause had occurred after age 50 and had no antecedents regarding gynecological disease. A further 345 women from the same ethnic origin (general population group) were also recruited to assess allele frequency for potentially deleterious sequence variants.

INTERVENTION(S)

Next generation sequencing (NGS), Sanger sequencing, and bioinformatics analysis.

MAIN OUTCOME MEASURE(S)

The complete coding regions of 70 candidate genes were massively sequenced, via NGS, in POF patients. Bioinformatics and genetics were used to confirm NGS results and to identify potential sequence variants related to the disease pathogenesis.

RESULT(S)

We have identified mutations in two novel genes, ADAMTS19 and BMPR2, that are potentially related to POF origin. LHCGR mutations, which might have contributed to the phenotype, were also detected.

CONCLUSION(S)

We thus recommend NGS as a powerful tool for identifying new molecular actors in POF and for future diagnostic/prognostic purposes.

Pleiotropy analysis of quantitative traits at gene level by multivariate functional linear models

In genetics, pleiotropy describes the genetic effect of a single gene on multiple phenotypic traits. A common approach is to analyze the phenotypic traits separately using univariate analyses and combine the test results through multiple comparisons. This approach may lead to low power. Multivariate functional linear models are developed to connect genetic variant data to multiple quantitative traits adjusting for covariates for a unified analysis. Three types of approximate F-distribution tests based on Pillai-Bartlett trace, Hotelling-Lawley trace, and Wilks's Lambda are introduced to test for association between multiple quantitative traits and multiple genetic variants in one genetic region. The approximate F-distribution tests provide much more significant results than those of F-tests of univariate analysis and optimal sequence kernel association test (SKAT-O). Extensive simulations were performed to evaluate the false positive rates and power performance of the proposed models and tests. We show that the approximate F-distribution tests control the type I error rates very well. Overall, simultaneous analysis of multiple traits can increase power performance compared to an individual test of each trait. The proposed methods were applied to analyze (1) four lipid traits in eight European cohorts, and (2) three biochemical traits in the Trinity Students Study. The approximate F-distribution tests provide much more significant results than those of F-tests of univariate analysis and SKAT-O for the three biochemical traits. The approximate F-distribution tests of the proposed functional linear models are more sensitive than those of the traditional multivariate linear models that in turn are more sensitive than SKAT-O in the univariate case. The analysis of the four lipid traits and the three biochemical traits detects more association than SKAT-O in the univariate case.

The power of gene-based rare variant methods to detect disease-associated variation and test hypotheses about complex disease

Genome and exome sequencing in large cohorts enables characterization of the role of rare variation in complex diseases. Success in this endeavor, however, requires investigators to test a diverse array of genetic hypotheses which differ in the number, frequency and effect sizes of underlying causal variants. In this study, we evaluated the power of gene-based association methods to interrogate such hypotheses, and examined the implications for study design. We developed a flexible simulation approach, using 1000 Genomes data, to (a) generate sequence variation at human genes in up to 10K case-control samples, and (b) quantify the statistical power of a panel of widely used gene-based association tests under a variety of allelic architectures, locus effect sizes, and significance thresholds. For loci explaining ~1% of phenotypic variance underlying a common dichotomous trait, we find that all methods have low absolute power to achieve exome-wide significance (~5-20% power at α=2.5×10^-6) in 3K individuals; even in 10K samples, power is modest (~60%). The combined application of multiple methods increases sensitivity, but does so at the expense of a higher false positive rate. MiST, SKAT-O, and KBAC have the highest individual mean power across simulated datasets, but we observe wide architecture-dependent variability in the individual loci detected by each test, suggesting that inferences about disease architecture from analysis of sequencing studies can differ depending on which methods are used. Our results imply that tens of thousands of individuals, extensive functional annotation, or highly targeted hypothesis testing will be required to confidently detect or exclude rare variant signals at complex disease loci.

Re-sequencing technologies allow for a more complete interrogation of the role of human variation in complex disease. The inadequate power of single variant methods to assess the role of less common variation has led to the development of numerous statistical methods for testing aggregate groups of variants for association with disease. Such endeavors pose substantial analytical challenges, however, due to the diverse array of genetic hypotheses that need to be considered. In this work, we systematically quantify and compare the performance of a panel of commonly used gene-based association methods under a range of allelic architectures, significance thresholds, locus effect sizes, sample sizes, and filters for neutral variation. We find that MiST, SKAT-O, and KBAC have the highest mean power across simulated datasets. Across all methods, however, the power to detect even loci of relatively large effect is very low at exome-wide significance thresholds for sample sizes comparable with those of ongoing sequencing studies; as such, the absence of signal in studies of a few thousand individuals does not exclude a role for rare variation in complex traits. Finally, we directly compare the results reported by different gene-based methods in order to identify their comparative advantages and disadvantages under distinct locus architectures. Our findings have implications for meaningful interpretation of both positive and negative findings in ongoing and future sequencing studies.

Rare variants in MYD88, IRAK4 and IKBKG and susceptibility to invasive pneumococcal disease: a population-based case-control study

Although rare variants within the Toll-like receptor signalling pathway genes have been found to underlie human primary immunodeficiencies associated with selective predisposition to invasive pneumococcal disease (IPD), the contribution of variants in these genes to IPD susceptibility at the population level remains unknown. Complete re-sequencing of IRAK4, MYD88 and IKBKG genes was undertaken in 164 IPD cases from the UK and 164 geographically-matched population-based controls. 233 single-nucleotide variants (SNVs) were identified, of which ten were in coding regions. Four rare coding variants were predicted to be deleterious, two variants in MYD88 and two in IRAK4. The predicted deleterious variants in MYD88 were observed as two heterozygote cases but not seen in controls. Frequencies of predicted deleterious IRAK4 SNVs were the same in cases and controls. Our findings suggest that rare, functional variants in MYD88, IRAK4 or IKBKG do not significantly contribute to IPD susceptibility in adults at the population level.

Whole-exome SNP array identifies 15 new susceptibility loci for psoriasis

Genome-wide association studies (GWASs) have reproducibly associated ∼40 susceptibility loci with psoriasis. However, the missing heritability is evident and the contributions of coding variants have not yet been systematically evaluated. Here, we present a large-scale whole-exome array analysis for psoriasis consisting of 42,760 individuals. We discover 16 SNPs within 15 new genes/loci associated with psoriasis, including C1orf141, ZNF683, TMC6, AIM2, IL1RL1, CASR, SON, ZFYVE16, MTHFR, CCDC129, ZNF143, AP5B1, SYNE2, IFNGR2 and 3q26.2-q27 (P<5.00 × 10(-08)). In addition, we also replicate four known susceptibility loci TNIP1, NFKBIA, IL12B and LCE3D-LCE3E. These susceptibility variants identified in the current study collectively account for 1.9% of the psoriasis heritability. The variant within AIM2 is predicted to impact protein structure. Our findings increase the number of genetic risk factors for psoriasis and highlight new and plausible biological pathways in psoriasis.