Low-frequency Coding Variants Research Articles

Genome scan reveals several loci associated with torus palatinus.

Torus Palatinus (TP) is a common trait with an unclear aetiology. Although prior studies suggest a hereditary component, the genetic factors that influence TP risk remain unknown. The purpose of this study is to identify genetic variants associated with TP. We assessed the TP status of 829 individuals from various ancestral backgrounds using 3D palate scans. We then carried out a genome-wide association study (GWAS) to identify common variants associated with TP. We also performed gene-based tests across the exome to investigate the role of low-frequency coding variants. Our GWAS did not identify any genome-wide significant signals but identified suggestive associations including hits on chromosomes 2, 5 and 17 with p-values less than 5 × 10-6. Candidate genes at these suggestive loci have been implicated in normal-range craniofacial features, syndromes with facial and oral malformations, and bone density. We did not find evidence that low-frequency coding variants influence TP risk. In addition, we failed to replicate associations identified in prior genetic studies of TP. These findings suggest that multiple genes likely influence the development of TP. Independent replication will be required to confirm our suggestive associations.

Autoimmune hypothyroidism GWAS reveals independent autoimmune and thyroid-specific contributions and an inverse relation with cancer risk.

The high prevalence of autoimmune hypothyroidism (AIHT) - more than 5% in human populations - provides a unique opportunity to unlock the most complete picture to date of genetic loci that underlie systemic and organ-specific autoimmunity. Using a meta-analysis of 81,718 AIHT cases in FinnGen and the UK Biobank, we dissect associations along axes of thyroid dysfunction and autoimmunity. This largest-to-date scan of hypothyroidism identifies 418 independent associations (p < 5x10 - 8 ), more than half of which have not previously been documented in thyroid disease. In 48 of these, a protein-coding variant is the lead SNP or is highly correlated (r 2 > 0.95) with the lead SNP at the locus, including low-frequency coding variants at LAG3 , ZAP70 , TG, TNFSF11, IRF3, S1PR4, HABP2, ZNF429 as well as established variants at ADCY7, IFIH1 and TYK2 . The variants at LAG3 (P67T), ZAP70 (T155M), and TG (Q655X) are highly enriched in Finland and functional experiments in T-cells demonstrate that the ZAP70 :T155M allele reduces T-cell activation. By employing a large-scale scan of non-thyroid autoimmunity and a published meta-analysis of TSH levels, we use a Bayesian classifier to dissect the associated loci into distinct groupings and from this estimate, a significant proportion are involved in systemic (i.e., general to multiple autoimmune conditions) autoimmunity (34%) and another subset in thyroid-specific dysfunction (17%). By comparing these association results further to other common disease endpoints, we identify a noteworthy overlap with skin cancer, with 10% of AIHT loci showing a consistent but opposite pattern of association where alleles that increase the risk of hypothyroidism have protective effects for skin cancer. The association results, including genes encoding checkpoint inhibitors and other genes affecting protein levels of PD1, bolster the causal role of natural variation in autoimmunity influencing cancer outcomes.

Low-frequency Coding Variants Associated With Body Mass Index Affect the Success of Bariatric Surgery.

A recent study identified 14 low-frequency coding variants associated with body mass index (BMI) in 718 734 individuals predominantly of European ancestry. We investigated the association of 2 genetic scores (GS) with i) the risk of severe/morbid obesity, ii) BMI variation before weight-loss intervention, iii) BMI change in response to an 18-month lifestyle/behavioral intervention program, and iv) BMI change up to 24 months after bariatric surgery. The 14 low-frequency coding variants were genotyped or sequenced in 342 French adults with severe/morbid obesity and 574 French adult controls from the general population. We built risk and protective GS based on 6 BMI-increasing and 5 BMI-decreasing low-frequency coding variants that were polymorphic in our study. While the risk GS was not associated with severe/morbid obesity status, BMI-decreasing low-frequency coding variants were significantly less frequent in patients with severe/morbid obesity than in French adults from the general population. Neither the risk nor the protective GS was associated with BMI before intervention in patients with severe/morbid obesity, nor did they affect BMI change in response to a lifestyle/behavioral modification program. The protective GS was associated with a greater BMI decrease following bariatric surgery. The risk and protective GS were associated with a higher and lower risk of BMI regain after bariatric surgery. Our data indicate that in populations of European descent, low-frequency coding variants associated with BMI in the general population also affect the outcomes of bariatric surgery in patients with severe/morbid obesity.

Gene expression atlas of energy balance brain regions.

Energy balance is controlled by interconnected brain regions in the hypothalamus, brainstem, cortex, and limbic system. Gene expression signatures of these regions can help elucidate the pathophysiology underlying obesity. RNA sequencing was conducted on P56 C57BL/6NTac male mice and E14.5 C57BL/6NTac embryo punch biopsies in 16 obesity-relevant brain regions. The expression of 190 known obesity-associated genes (monogenic, rare, and low-frequency coding variants; GWAS; syndromic) was analyzed in each anatomical region. Genes associated with these genetic categories of obesity had localized expression patterns across brain regions. Known monogenic obesity causal genes were highly enriched in the arcuate nucleus of the hypothalamus and developing hypothalamus. The obesity-associated genes clustered into distinct “modules” of similar expression profile, and these were distinct from expression modules formed by similar analysis with genes known to be associated with other disease phenotypes (type 1 and type 2 diabetes, autism, breast cancer) in the same energy balance–relevant brain regions.

778 EXOME SEQUENCING IN 30,000 CASES DEFINES NOVEL RISK FACTORS FOR CROHN'S DISEASE

Background Genome-wide association studies (GWASs) have identified >200 genomic regions associated with IBD. These signals are mainly driven by common variants in noncoding regions, making it a challenge to extract biological insight. Targeted sequencing of genes in these regions has successfully identified putative causal variants, often rarer and independent of the initial GWAS hit. These coding variants have led to more direct functional experiments demonstrating causal mechanisms for at least ten genes. Methods To further rare variant-based discovery, we pursued large-scale exome sequencing of more than 30,000 Crohn's disease (CD) cases and population controls from more than 20 centers in the International IBD Genetics Consortium. Sequencing was performed at the Broad Institute using both Nextera (11,125 CD cases) and TWIST (7,442 CD cases) exome captures - and at the Sanger Institute (whole genomes 6,404 CD cases and Agilent exome 3,848 CD cases). Further replication was performed with a sample of 4,359 cases from Kiel sequenced at Regeneron. In each technical experiment, a comparable or larger number of sequenced controls were available for genetic association analysis. Results We conducted meta-analysis of two exome captures at the Broad Institute. To expand on GWAS, we focused on rare and low-frequency coding variants between 0.01% and 10% and estimated roughly 85% of all protein coding variants in this frequency range are reliably analyzed in both exome captures. 119 variants (figure below) were identified with p Conclusion These findings provide novel launch points for mechanistic studies that will further our understanding of disease pathogenesis. For example, SDF2L1 (R161H) flags a gene not previously implicated in IBD but which is reported to regulate feeding-induced ER stress, with a series of CRISPR screens identifying it as an essential regulator of ER homeostasis. Download : Download high-res image (141KB) Download : Download full-size image Effect size versus frequency for significant low frequency coding variants (credit: Hailiang Huang)

Multi-trait association studies discover pleiotropic loci between Alzheimer\u2019s disease and cardiometabolic traits

BackgroundIdentification of genetic risk factors that are shared between Alzheimer’s disease (AD) and other traits, i.e., pleiotropy, can help improve our understanding of the etiology of AD and potentially detect new therapeutic targets. Previous epidemiological correlations observed between cardiometabolic traits and AD led us to assess the pleiotropy between these traits.MethodsWe performed a set of bivariate genome-wide association studies coupled with colocalization analysis to identify loci that are shared between AD and eleven cardiometabolic traits. For each of these loci, we performed colocalization with Genotype-Tissue Expression (GTEx) project expression quantitative trait loci (eQTL) to identify candidate causal genes.ResultsWe identified three previously unreported pleiotropic trait associations at known AD loci as well as four novel pleiotropic loci. One associated locus was tagged by a low-frequency coding variant in the gene DOCK4 and is potentially implicated in its alternative splicing. Colocalization with GTEx eQTL data identified additional candidate genes for the loci we detected, including ACE, the target of the hypertensive drug class of ACE inhibitors. We found that the allele associated with decreased ACE expression in brain tissue was also associated with increased risk of AD, providing human genetic evidence of a potential increase in AD risk from use of an established anti-hypertensive therapeutic.ConclusionOur results support a complex genetic relationship between AD and these cardiometabolic traits, and the candidate causal genes identified suggest that blood pressure and immune response play a role in the pleiotropy between these traits.

Impact of low-frequency coding variants on human facial shape

The contribution of low-frequency variants to the genetic architecture of normal-range facial traits is unknown. We studied the influence of low-frequency coding variants (MAF < 1%) in 8091 genes on multi-dimensional facial shape phenotypes in a European cohort of 2329 healthy individuals. Using three-dimensional images, we partitioned the full face into 31 hierarchically arranged segments to model facial morphology at multiple levels, and generated multi-dimensional phenotypes representing the shape variation within each segment. We used MultiSKAT, a multivariate kernel regression approach to scan the exome for face-associated low-frequency variants in a gene-based manner. After accounting for multiple tests, seven genes (AR, CARS2, FTSJ1, HFE, LTB4R, TELO2, NECTIN1) were significantly associated with shape variation of the cheek, chin, nose and mouth areas. These genes displayed a wide range of phenotypic effects, with some impacting the full face and others affecting localized regions. The missense variant rs142863092 in NECTIN1 had a significant effect on chin morphology and was predicted bioinformatically to have a deleterious effect on protein function. Notably, NECTIN1 is an established craniofacial gene that underlies a human syndrome that includes a mandibular phenotype. We further showed that nectin1a mutations can affect zebrafish craniofacial development, with the size and shape of the mandibular cartilage altered in mutant animals. Findings from this study expanded our understanding of the genetic basis of normal-range facial shape by highlighting the role of low-frequency coding variants in several novel genes.

Proprotein convertase 7 (PCSK7) reduces apoA-V levels.

The locus of the human proprotein convertase subtilisin-kexin type-7 (PC7) gene (PCSK7) is on chromosome 11q23.3 close to the gene cluster APOA5/APOA4/APOC3/APOA1, a region implicated in the regulation of lipoprotein metabolism. A GWAS reported the association of PCSK7 SNPs with plasma triglyceride(TG), and exome sequencing of African Americans revealed the association of a low-frequency coding variant of PC7 (R504H; SNP rs142953140) with a ~30% TG reduction. Another PCSK7 SNP rs508487 is in linkage disequilibrium with a promoter variant of the liver-derived apolipoprotein A-V (apoA-V), an indirect activator of the lipoprotein lipase (LpL), and is associated with elevated TG levels. We thus hypothesized that PC7 regulates the levels/activity of apoA-V. Studies in the human hepatic cell line HuH7 revealed that wild-type (WT) PC7 and its endoplasmic reticulum (ER)-retained forms bind to and enhance the degradation of human apoA-V in acidic lysosomes in a nonenzymatic fashion. PC7-induced degradation of apoA-V is inhibited by bafilomycin A1 and the alkalinizing agents: chloroquine and NH4 Cl. Thus, the PC7-induced apoA-V degradation implicates an ER-lysosomal communication inhibited by bafilomycin A1. Invitro, the natural R504H mutant enhances PC7 Ser505 phosphorylation at the structurally exposed Ser-X-Glu507 motif recognized by the secretory kinase Fam20C. Co-expression of the phosphomimetic PC7-S505E with apoA-V resulted in lower degradation compared to WT, suggesting that Ser505 phosphorylation of PC7 lowers TG levels via reduced apoA-V degradation. In agreement, in Pcsk7-/- mice fed high-fat diet, plasma apoA-V levels and adipocyte LpL activity are increased, providing an invivo mechanistic link for a role of liver PC7 in enhanced TG storage in adipocytes.

Sequence variants with large effects on cardiac electrophysiology and disease

Features of the QRS complex of the electrocardiogram, reflecting ventricular depolarisation, associate with various physiologic functions and several pathologic conditions. We test 32.5 million variants for association with ten measures of the QRS complex in 12 leads, using 405,732 electrocardiograms from 81,192 Icelanders. We identify 190 associations at 130 loci, the majority of which have not been reported before, including associations with 21 rare or low-frequency coding variants. Assessment of genes expressed in the heart yields an additional 13 rare QRS coding variants at 12 loci. We find 51 unreported associations between the QRS variants and echocardiographic traits and cardiovascular diseases, including atrial fibrillation, complete AV block, heart failure and supraventricular tachycardia. We demonstrate the advantage of in-depth analysis of the QRS complex in conjunction with other cardiovascular phenotypes to enhance our understanding of the genetic basis of myocardial mass, cardiac conduction and disease.

Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis Risk

(Cell 175, 1679–1687.e1–e7; November 29, 2018) We the editors of Cell have been informed by the Lead Contact and corresponding author of this paper, Dr. Chris Cotsapas, that further analyses following the publication of the paper have indicated that two variants (rs61999302 encoding PRKRA.D33G and rs62176112 encoding PRKRA.P11L) reported in this paper to be associated with MS risk were in fact spurious associations due to the presence of a dispersed duplication event, in which this region of chromosome 2 is duplicated in the MHC region of chromosome 6. The potential for this dispersed duplication to lead to spurious associations has been previously documented (Wellcome Trust Case Control Consortium et al., 2010Craddock N. Hurles M.E. Cardin N. Pearson R.D. Plagnol V. Robson S. Vukcevic D. Barnes C. Conrad D.F. Giannoulatou E. et al.Wellcome Trust Case Control ConsortiumGenome-wide association study of CNVs in 16,000 cases of eight common diseases and 3,000 shared controls.Nature. 2010; 464: 713-720Google Scholar). Upon learning of the potential problems associated with signals in this region, the authors used the data reported in this paper to impute classical HLA alleles and amino acid variants using SNP2HLA (Jia et al., 2013Jia X. Han B. Onengut-Gumuscu S. Chen W.M. Concannon P.J. Rich S.S. Raychaudhuri S. de Bakker P.I.W. Imputing amino acid polymorphisms in human leukocyte antigens.PLoS ONE. 2013; 8: e64683Google Scholar) and found that the two reported variants in PRKRA showed extensive linkage disequilibrium with genotyped and imputed HLA variants, including known MS risk alleles. Conditioning on any of these known variants was sufficient to explain the reported association, leading to the conclusion that these two signals are driven by the duplication and are not true associations of the reported variants with MS risk. With this note, we and the authors would like to alert readers to this issue. The remaining findings in the paper are unchanged, as are the main conclusions of the paper. The authors would also like to thank Dr. Matthew Hurles (Wellcome Sanger Institute) for bringing this issue to their attention. Furthermore, in the printed version of the experimental model and subject details section of the STAR methods, Köln is listed as one of the patient recruitment centers, but Düsseldorf should have been listed instead. One of the authors, Dr. Clemens Warnke, although currently based in Köln, in fact contributed patient samples collected in Düsseldorf to the study. This has been corrected in the article online, and the authors apologize for any confusion these issues may have caused. Low-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis RiskMitrovič et al.CellOctober 18, 2018In BriefIn a large multi-cohort study, unexplained heritability for multiple sclerosis is detected in low-frequency coding variants that are missed by GWAS analyses, further underscoring the role of immune genes in MS pathology. Full-Text PDF Open AccessLow-Frequency and Rare-Coding Variation Contributes to Multiple Sclerosis RiskMitrovič et al.CellJanuary 23, 2020In Brief(Cell 175, 1679–1687.e1–e7; November 29, 2018) Full-Text PDF Open Access

Exploring rare and low-frequency variants in the Saguenay\u2013Lac-Saint-Jean population identified genes associated with asthma and allergy traits

The Saguenay–Lac-Saint-Jean (SLSJ) region is located in northeastern Quebec and is known for its unique demographic history and founder effect. As founder populations are enriched with population-specific variants, we characterized the variants distribution in SLSJ and compared it with four European populations (Finnish, Sweden, United Kingdom and France), of which the Finnish population is another founder population. Targeted sequencing of the coding and non-coding immune regulatory regions of the SLSJ asthma familial cohort and the four European populations were performed. Rare and low-frequency coding and non-coding regulatory variants identified in the SLSJ population were then investigated for variant- and gene-level associations with asthma and allergy-related traits (eosinophil percentage, immunoglobulin (Ig) E levels and lung function). Our data showed that (1) rare or deleterious variants were not enriched in the two founder populations as compared with the three non-founder European populations; (2) a larger proportion of founder population-specific variants occurred with higher frequencies; and (3) low-frequency variants appeared to be more deleterious. Furthermore, a rare variant, rs1386931, located in the 3ʹ-UTR of CXCR6 and intron of FYCO1 was found to be associated with eosinophil percentage. Gene-based analyses identified NRP2, MRPL44 and SERPINE2 to be associated with various asthma and allergy-related traits. Our study demonstrated the usefulness of using a founder population to identify new genes associated with asthma and allergy-related traits; thus better understand the genes and pathways implicated in pathophysiology.

Identification of a missense variant in the WFS1 gene that causes a mild form of Wolfram syndrome and is associated with risk for type 2 diabetes in Ashkenazi Jewish individuals.

Wolfram syndrome is a rare, autosomal recessive syndrome characterised by juvenile-onset diabetes and optic atrophy and is caused by bi-allelic mutations in the WFS1 gene. In a recent sequencing study, an individual with juvenile-onset diabetes was observed to be homozygous for a rare missense variant (c.1672C>T, p.R558C) in the WFS1 gene. The aim of this study was to perform the genetic characterisation of this variant and to determine whether it is causal for young-onset diabetes and Wolfram syndrome. We analysed the allele frequency of the missense variant in multiple variant databases. We genotyped the variant in 475 individuals with type 1 diabetes and 2237 control individuals of Ashkenazi Jewish ancestry and analysed the phenotypes of homozygotes. We also investigated the association of this variant with risk for type 2 diabetes using genotype and sequence data for type 2 diabetes cases and controls. The missense variant demonstrated an allele frequency of 1.4% in individuals of Ashkenazi Jewish ancestry, 60-fold higher than in other populations. Genotyping of this variant in 475 individuals diagnosed with type 1 diabetes identified eight homozygotes compared with none in 2237 control individuals (genotype relative risk 135.3, p = 3.4 × 10-15). The age at diagnosis of diabetes for these eight individuals (17.8 ± 8.3years) was several times greater than for typical Wolfram syndrome (5 ± 4years). Further, optic atrophy was observed in only one of the eight individuals, while another individual had the Wolfram syndrome-relevant phenotype of neurogenic bladder. Analysis of sequence and genotype data in two case-control cohorts of Ashkenazi ancestry demonstrated that this variant is also associated with an increased risk of type 2 diabetes in heterozygotes (OR 1.81, p = 0.004). We have identified a low-frequency coding variant in the WFS1 gene that is enriched in Ashkenazi Jewish individuals and causes a mild form of Wolfram syndrome characterised by young-onset diabetes and reduced penetrance for optic atrophy. This variant should be considered for genetic testing in individuals of Ashkenazi ancestry diagnosed with young-onset non-autoimmune diabetes and should be included in Ashkenazi carrier screening panels.

Coding variants in RPL3L and MYZAP increase risk of atrial fibrillation

Most sequence variants identified hitherto in genome-wide association studies (GWAS) of atrial fibrillation are common, non-coding variants associated with risk through unknown mechanisms. We performed a meta-analysis of GWAS of atrial fibrillation among 29,502 cases and 767,760 controls from Iceland and the UK Biobank with follow-up in samples from Norway and the US, focusing on low-frequency coding and splice variants aiming to identify causal genes. We observe associations with one missense (OR = 1.20) and one splice-donor variant (OR = 1.50) in RPL3L, the first ribosomal gene implicated in atrial fibrillation to our knowledge. Analysis of 167 RNA samples from the right atrium reveals that the splice-donor variant in RPL3L results in exon skipping. We also observe an association with a missense variant in MYZAP (OR = 1.38), encoding a component of the intercalated discs of cardiomyocytes. Both discoveries emphasize the close relationship between the mechanical and electrical function of the heart.

Recurrent, low-frequency coding variants contributing to colorectal cancer in the Swedish population.

Genome-wide association studies (GWAS) have identified dozens of common genetic variants associated with risk of colorectal cancer (CRC). However, the majority of CRC heritability remains unclear. In order to discover low-frequency, high-risk CRC susceptibility variants in Swedish population, we genotyped 1 515 CRC patients enriched for familial cases, and 12 108 controls. Case/control association analysis suggested eight novel variants associated with CRC risk (OR 2.0–17.6, p-value < 2.0E-07), comprised of seven coding variants in genes RAB11FIP5, POTEA, COL27A1, MUC5B, PSMA8, MYH7B, and PABPC1L as well as one variant downstream of NEU1 gene. We also confirmed 27 out of 30 risk variants previously reported from GWAS in CRC with a mixed European population background. This study identified rare, coding sequence variants associated with CRC risk through analysis in a relatively homogeneous population. The segregation data suggest a complex mode of inheritance in seemingly dominant pedigrees.

P164 USING EXOME SEQUENCING TO EXPAND THE GENETIC ARCHITECTURE OF INFLAMMATORY BOWEL DISEASE

IBD studies over the past decade have confirmed association to 250 gene loci. A handful of associations have led to specific validated functional variants highlighting intracellular response to microbes and regulation of adaptive immunity in IBD pathogenesis. For the vast majority, however, the specific implicated gene and causal functional variants remain unknown. In partnership with IBD researchers around the world, we launched an exome sequencing initiative to define the full allelic spectrum of protein-altering variation in genes associated to CD and/or UC, assess their role in clinical course and response to therapy, and to determine whether truncating variants confer risk or protection in each IBD gene in order to highlight opportune therapeutic targets. We have already completed exomes of 13,000 IBD cases providing a high-resolution view of coding variation at each GWAS hit and demonstrated a convincing excess of rare exome signal collectively. The cases are drawn from individual substudies that focus on isolated populations (Ashkenazi, Finnish, French-Canadian), admixed populations and clinical extremes, each providing unique opportunities for discovery. Early findings include novel protective truncating variants, the complete allelic series (including unique founder population alleles), non-additive inheritance models at known genes such as NOD2, overlooked low-frequency coding variants that explain GWAS hits (PRDM1, NRIP1). The integration of low frequency and rare functional variants with GWAS is moving us much closer towards completion of the genetic architecture of IBD. Early analyses of the data have revealed novel non-coding associations in known IBD genes as well as coding variants that implicate genes at loci in which the ‘causative’ gene was unknown. Furthermore, initial results have suggested variation in enrichment of different processes across different ethnicities.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1–5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

The GBA variant E326K is associated with Parkinson's disease and explains a genome-wide association signal.

Coding variants in the GBA gene have been identified as the numerically most important genetic risk factors for Parkinson's disease (PD). In addition, genome-wide association studies (GWAS) have identified associations with PD in the SYT11-GBA region on chromosome 1q22, but the relationship to GBA coding variants have remained unclear. The aim of this study was to sequence the complete GBA gene in a clinical cohort and to investigate whether coding variants within the GBA gene may be driving reported association signals. We analyzed high-throughput sequencing data of all coding exons of GBA in 366 patients with PD. The identified low-frequency coding variants were genotyped in three Scandinavian case-controls series (786 patients and 713 controls). Previously reported risk variants from two independent association signals within the SYT11-GBA locus on chromosome 1 were also genotyped in the same samples. We performed association analyses and evaluated linkage disequilibrium (LD) between the variants. We identified six rare mutations (1.6%) and two low-frequency coding variants in GBA. E326K (rs2230288) was significantly more frequent in PD patients compared to controls (OR 1.65, p=0.03). There was no clear association of T369M (rs75548401) with disease (OR 1.43, p=0.24). Genotyping the two GWAS hits rs35749011 and rs114138760 in the same sample set, we replicated the association between rs35749011 and disease status (OR 1.67, p=0.03), while rs114138760 was found to have similar allele frequencies in patients and controls. Analyses revealed that E326K and rs35749011 are in very high LD (r2 0.95). Our results confirm that the GBA variant E326K is a susceptibility allele for PD. The results suggest that E326K may fully account for the primary association signal observed at chromosome 1q22 in previous GWAS of PD.

Whole-genome sequencing study of serum peptide levels: the Atherosclerosis Risk in Communities study.

Oligopeptides are important markers of protein metabolism, as they are cleaved from larger polypeptides and proteins. Genetic association studies may help elucidate their origin and function. In 1,552 European Americans and 1,872 African Americans of the Atherosclerosis Risk in Communities study, we performed whole-genome and whole-exome sequencing and measured serum levels of 25 peptides. Common variants (minor allele frequency > 5%) were analysed individually. We grouped low-frequency variants (minor allele frequency ≤ 5%) by a genome-wide sliding window using region-based aggregate tests. Furthermore, low-frequency regulatory variants were grouped by gene, as were functional coding variants. All analyses were performed separately in each ancestry group and then meta-analysed. We identified 22 common variant associations with peptide levels (P-value < 4.2 × 10-10), including 16 novel gene-peptide pairs. Notably, variants in kinin-kallikrein genes KNG1, F12, KLKB1, and ACE were associated with several different peptides. Variants in KLKB1 and ACE were associated with a fragment of complement component 3f. Both common variants and low-frequency coding variants in CPN1 were associated with a fibrinogen cleavage peptide. Four sliding windows were significantly associated with peptide levels (P-value < 4.2 × 10-10). Our results highlight the importance of the kinin-kallikrein system in the regulation of serum peptide levels, strengthen the evidence for a broad link between the kinin-kallikrein and complement systems, and suggest a role of CPN1 in the conversion of fibrinogen to fibrin.

Association studies of low-frequency coding variants in nonsyndromic cleft lip with or without cleft palate.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a group of common human birth defects with complex etiology. Although genome-wide association studies have successfully identified a number of risk loci, these loci only account for about 20% of the heritability of orofacial clefts. The "missing" heritability may be found in rare variants, copy number variants, or interactions. In this study, we investigated the role of low-frequency variants genotyped in 1995 cases and 1626 controls on the Illumina HumanCore + Exome chip. We performed two statistical tests, Sequence Kernel Association Test (SKAT) and Combined Multivariate and Collapsing (CMC) method using two minor allele frequency cutoffs (1% and 5%). We found that a burden of low-frequency coding variants in N4BP2, CDSN, PRTG, and AHRR were associated with increased risk of NSCL/P. Low-frequency variants in other genes were associated with decreased risk of NSCL/P. These results demonstrate that low-frequency variants contribute to the genetic etiology of NSCL/P.

A low-frequency variant in SMAD7 modulates TGF-β signaling and confers risk for colorectal cancer in Chinese population.

The TGF-β pathway plays an essential role in regulating cell proliferation and differentiation. GWASs and candidate approaches have identified a battery of genetic variants in the TGF-β pathway contributing to colorectal cancer (CRC). However, most of the significant variants are common variants and their functions remain ambiguous. To identify causal variants with low-frequency in the TGF-β pathway contributing to CRC susceptibility in Chinese population, we performed targeted sequencing of 12 key genes in TGF-β signaling in CRC patients followed by a two-stage case-control study with a total of 5109 cases and 5169 controls. Bioinformatic annotations and biochemical experiments were applied to reveal the potential functions of significant variants. Seven low-frequency genetic variants were captured through targeted sequencing. The two stage association studies showed that missense variant rs3764482 (c. 83C>T; p. S28F) in the gene SMAD7 was consistently and significantly associated with CRC risk. Compared with the wild type, the ORs for variant allele were 1.37 (95%CI: 1.10-1.70, P = 0.005), 1.55 (95%CI: 1.30-1.86, P = 1.15 × 106 ), and 1.48 (1.29-1.70, P = 2.44 × 10;8 ) in stage 1, stage 2, and the combined analyses, respectively. Functional annotations revealed that the minor allele T of rs3764482 was more effective than the major allele C in blocking the TGF-β signaling and inhibiting the phosphorylation of receptor-regulated SMADs (R-SMADs). In conclusion, low-frequency coding variant rs3764482 in SMAD7 is associated with CRC risk in Chinese population. The rs3764482 variant may block the TGF-β signaling via impeding the activation of downstream genes, leading to cancer cell proliferation, thus contributing to CRC pathogenesis.