Low-frequency Nonsynonymous Variants Research Articles

Characterizing rare and low-frequency height-associated variants in the Japanese population

Human height is a representative phenotype to elucidate genetic architecture. However, the majority of large studies have been performed in European population. To investigate the rare and low-frequency variants associated with height, we construct a reference panel (N = 3,541) for genotype imputation by integrating the whole-genome sequence data from 1,037 Japanese with that of the 1000 Genomes Project, and perform a genome-wide association study in 191,787 Japanese. We report 573 height-associated variants, including 22 rare and 42 low-frequency variants. These 64 variants explain 1.7% of the phenotypic variance. Furthermore, a gene-based analysis identifies two genes with multiple height-increasing rare and low-frequency nonsynonymous variants (SLC27A3 and CYP26B1; PSKAT-O < 2.5 × 10−6). Our analysis shows a general tendency of the effect sizes of rare variants towards increasing height, which is contrary to findings among Europeans, suggesting that height-associated rare variants are under different selection pressure in Japanese and European populations.

Exome Chip Analysis Identifies Low-Frequency and Rare Variants in MRPL38 for White Matter Hyperintensities on Brain Magnetic Resonance Imaging.

Background and Purpose- White matter hyperintensities (WMH) on brain magnetic resonance imaging are typical signs of cerebral small vessel disease and may indicate various preclinical, age-related neurological disorders, such as stroke. Though WMH are highly heritable, known common variants explain a small proportion of the WMH variance. The contribution of low-frequency/rare coding variants to WMH burden has not been explored. Methods- In the discovery sample we recruited 20 719 stroke/dementia-free adults from 13 population-based cohort studies within the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, among which 17 790 were of European ancestry and 2929 of African ancestry. We genotyped these participants at ≈250 000 mostly exonic variants with Illumina HumanExome BeadChip arrays. We performed ethnicity-specific linear regression on rank-normalized WMH in each study separately, which were then combined in meta-analyses to test for association with single variants and genes aggregating the effects of putatively functional low-frequency/rare variants. We then sought replication of the top findings in 1192 adults (European ancestry) with whole exome/genome sequencing data from 2 independent studies. Results- At 17q25, we confirmed the association of multiple common variants in TRIM65, FBF1, and ACOX1 ( P<6×10-7). We also identified a novel association with 2 low-frequency nonsynonymous variants in MRPL38 (lead, rs34136221; PEA=4.5×10-8) partially independent of known common signal ( PEA(conditional)=1.4×10-3). We further identified a locus at 2q33 containing common variants in NBEAL1, CARF, and WDR12 (lead, rs2351524; Pall=1.9×10-10). Although our novel findings were not replicated because of limited power and possible differences in study design, meta-analysis of the discovery and replication samples yielded stronger association for the 2 low-frequency MRPL38 variants ( Prs34136221=2.8×10-8). Conclusions- Both common and low-frequency/rare functional variants influence WMH. Larger replication and experimental follow-up are essential to confirm our findings and uncover the biological causal mechanisms of age-related WMH.

Targeted Sequencing Reveals Low-Frequency Variants in EPHA Genes as Markers of Paclitaxel-Induced Peripheral Neuropathy.

Purpose: Neuropathy is the dose-limiting toxicity of paclitaxel and a major cause for decreased quality of life. Genetic factors have been shown to contribute to paclitaxel neuropathy susceptibility; however, the major causes for interindividual differences remain unexplained. In this study, we identified genetic markers associated with paclitaxel-induced neuropathy through massive sequencing of candidate genes.Experimental Design: We sequenced the coding region of 4 EPHA genes, 5 genes involved in paclitaxel pharmacokinetics, and 30 Charcot-Marie-Tooth genes, in 228 cancer patients with no/low neuropathy or high-grade neuropathy during paclitaxel treatment. An independent validation series included 202 paclitaxel-treated patients. Variation-/gene-based analyses were used to compare variant frequencies among neuropathy groups, and Cox regression models were used to analyze neuropathy along treatment.Results: Gene-based analysis identified EPHA6 as the gene most significantly associated with paclitaxel-induced neuropathy. Low-frequency nonsynonymous variants in EPHA6 were present exclusively in patients with high neuropathy, and all affected the ligand-binding domain of the protein. Accumulated dose analysis in the discovery series showed a significantly higher neuropathy risk for EPHA5/6/8 low-frequency nonsynonymous variant carriers [HR, 14.60; 95% confidence interval (CI), 2.33-91.62; P = 0.0042], and an independent cohort confirmed an increased neuropathy risk (HR, 2.07; 95% CI, 1.14-3.77; P = 0.017). Combining the series gave an estimated 2.5-fold higher risk of neuropathy (95% CI, 1.46-4.31; P = 9.1 × 10-4).Conclusions: This first study sequencing EPHA genes revealed that low-frequency variants in EPHA6, EPHA5, and EPHA8 contribute to the susceptibility to paclitaxel-induced neuropathy. Furthermore, EPHA's neuronal injury repair function suggests that these genes might constitute important neuropathy markers for many neurotoxic drugs. Clin Cancer Res; 23(5); 1227-35. ©2016 AACR.

Strong genetic structure revealed by multilocus patterns of variation in Giardia duodenalis isolates of patients from Galicia (NW-Iberian Peninsula)

We report a survey of genetic variation at three coding loci in Giardia duodenalis of assemblages A and B obtained from stool samples of patients from Santiago de Compostela (Galicia, NW-Iberian Peninsula). The mean pooled synonymous diversity for assemblage A was nearly five times lower than for assemblage B (0.77%±0.30% and 4.14%±1.65%, respectively). Synonymous variation in both assemblages was in mutation-drift equilibrium and an excess of low-frequency nonsynonymous variants suggested the action of purifying selection at the three loci. Differences between isolates contributed to 40% and 60% of total genetic variance in assemblages A and B, respectively, which revealed a significant genetic structure. These results, together with the lack of evidence for recombination, support that (i) Giardia assemblages A and B are in demographic equilibrium and behave as two genetically isolated populations, (ii) infections are initiated by a reduced number of individuals, which may be genetically diverse and even belong to different assemblages, and (iii) parasites reproduce clonally within the host. However, the observation of invariant loci in some isolates means that mechanisms for the homogenization of the genetic content of the two diploid nuclei in each individual must exist.

Low‐frequency nonsynonymous variants in FKBPL and ARPC1B genes are associated with breast cancer risk in Chinese women

Genome-wide association studies have reported more than 100 independent common loci associated with breast cancer risk. The contribution of low-frequency or rare variants to breast cancer susceptibility has not been well explored. Thus, we applied exome chip to genotype >200 000 low-frequency and rare variants in 1064 breast cancer cases and 1125 cancer-free controls and subsequently validated promising associations in another 1040 breast cancer cases and 1240 controls. We identified two low-frequency nonsynonymous variants at FKBPL (rs200847762, OR = 0.34, 95% CI = 0.20-0.57, P = 4.31 × 10-5 ) and ARPC1B (rs1045012, OR = 0.56, 95% CI = 0.43-0.74, P = 4.30 × 10-5 ) associated with breast cancer risk. In stratification analyses, we found that the protective effect of rs200847762 was stronger in ER-positive breast cancer (OR = 0.18, 95% CI = 0.06-0.42) than that in ER-negative one (OR = 0.59, 95% CI = 0.31-1.05). Our findings indicate that low-frequency variants may also contribute to breast cancer susceptibility and genetic variants in 6p21.33 and 7q22.1 are important in breast carcinogenesis. © 2016 Wiley Periodicals, Inc.

Europeans have a higher proportion of high‑frequency deleterious variants than Africans.

Recent studies have shown that a high proportion of rare variants in European and African populations are deleterious in nature. However, the deleterious fraction of high-frequency variants is unclear. Using more than 6500 exomes we show a much higher fraction (11 %) of relatively high-frequency nonsynonymous (amino acid altering) variants (DAF 0.1–10 %) in European Americans (EA) compared to those from African Americans (AA). In contrast, this difference was only marginal (<2 %) for low-frequency nonsynonymous variants (DAF <0.1 %). Our results also revealed that the proportion of high-frequency deleterious nonsynonymous variants in EA was much higher (24 %) than that in AA and this difference was very small (4 %) for the low-frequency deleterious amino acid altering variants. We also show that EA have significantly more number of high-frequency deleterious nonsynonymous variants per genome than AA. The high proportion of high-frequency deleterious variants in EA could be the result of the well-known bottleneck experienced by European populations in which harmful mutations may have drifted to high frequencies. The estimated ages of deleterious variants support this prediction. Our results suggest that high-frequency variants could be relatively more likely to be associated with diseases in Europeans than in Africans and hence emphasize the need for population-specific strategies in genomic medicine studies.

Abstract 946: Exome genotyping array identifies rare and low-frequency variants that may be associated with ovarian cancer risk

Abstract Background: Common genetic variants identified by genome-wide association studies and follow-up genotyping initiatives only explain a small proportion of the heritability of epithelial ovarian cancer (EOC). We hypothesized that some of the missing heritability may be explained by rare (minor allele frequency (MAF) &amp;lt;0.5%) and low-frequency (MAF 0.5-5%) protein-coding variants that are not well captured by genotyping arrays that focus on common variants. Methods: To examine the association between uncommon variants and EOC susceptibility, we genotyped 7,060 EOC cases and 6,712 controls from the Ovarian Cancer Association Consortium using the Affymetrix Axiom Exome Genotyping Array. This platform contains 300,785 putative functional coding variants selected from exome sequences of 12,000+ individuals representing multiple ethnicities and complex traits. After performing quality control (QC), we used unconditional logistic regression treating the number of minor alleles carried as an ordinal variable (log-additive model) to evaluate SNP-EOC risk associations. Genotype cluster plots were visually inspected for top-ranking variants to determine assay success, and allele frequencies of samples from the 1000 Genomes Project and other cohorts were evaluated for QC purposes. Results: A total of 13,318 samples (97%) passed QC, and 12,779 eligible samples (6,293 invasive serous and endometrioid cases and 6,486 controls) were included in the analyses. 280,737 (93%) of the 300,785 attempted markers were successfully genotyped and had call rates &amp;gt;90%. After removing 31,436 monomorphic markers, the 249,301 remaining markers were included in association analyses. We confirmed variants at previously reported EOC susceptibility loci (9p22, 3q25, 17cen-q21.3, and 8q24), with p-values between 10-5 and 10-11. Several correlated low-frequency non-synonymous variants (MAF=5%) residing in a novel EOC susceptibility gene at 3q25.31 were risk-associated. The most significant variant at 3q25.31 had an OR=1.36 and p = 5.58x10-8. Also noteworthy is a rare variant at 15q15.1-q21.1 (MAF=0.5%) that introduces a stop codon in a gene belonging to the tyrosine kinase family (OR=0.39, CI=0.26-0.58, p = 4.52x10-6). We also identified a more common synonymous coding variant in a phosphatase at 6q21.3 that was associated with increased risk (OR=1.22, CI=1.13-1.33, p = 1.84x10-6, MAF=9.6%). Population stratification and gene-based analyses are underway to further examine the impact of the newly-identified variants on EOC risk. A meta-analysis to include data from an additional 2104 EOC cases and 2516 controls genotyped with Illumina's HumanExome Beadchip rare variant array is also underway. Functional analysis of the most promising variants will follow. Conclusions: Preliminary data from this large-scale study reveals novel rare and low-frequency variants that may influence susceptibility to epithelial ovarian cancer. Citation Format: Jennifer Permuth-Wey, Y. Ann Chen, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Simon Gayther, Ellen L. Goode, Edwin Iversen, Alvaro N.A. Monteiro, Leigh Pearce, Paul D.P. Pharoah, Catherine M. Phelan, Ailith Pirie, Susan Ramus, Mary Ann Rossing, Joellen M. Schildkraut, Thomas A. Sellers, on behalf of the Ovarian Cancer Association Consortium. Exome genotyping array identifies rare and low-frequency variants that may be associated with ovarian cancer risk. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 946. doi:10.1158/1538-7445.AM2014-946

SP0009 The Expanding Spectrum of Autoinflammation

The term autoinflammation was coined by McDermott and colleagues in 1999 to describe inflammatory conditions with apparently unprovoked inflammation without high titer autoantibodies or antigen-specific T lymphocytes. Unlike autoimmune diseases,...

Dense Genotyping of Candidate Gene Loci Identifies Variants Associated With High-Density Lipoprotein Cholesterol

Plasma levels of high-density lipoprotein cholesterol (HDL-C) are known to be heritable, but only a fraction of the heritability is explained. We used a high-density genotyping array containing single-nucleotide polymorphisms (SNPs) from HDL-C candidate genes selected on known biology of HDL-C metabolism, mouse genetic studies, and human genetic association studies. SNP selection was based on tagging SNPs and included low-frequency nonsynonymous SNPs. Association analysis in a cohort containing extremes of HDL-C (case-control, n=1733) provided a discovery phase, with replication in 3 additional populations for a total meta-analysis in 7857 individuals. We replicated the majority of loci identified through genome-wide association studies and present on the array (including ABCA1, APOA1/C3/A4/A5, APOB, APOE/C1/C2, CETP, CTCF-PRMT8, FADS1/2/3, GALNT2, LCAT, LILRA3, LIPC, LIPG, LPL, LRP4, SCARB1, TRIB1, ZNF664) and provide evidence that suggests an association in several previously unreported candidate gene loci (including ABCG1, GPR109A/B/81, NFKB1, PON1/2/3/4). There was evidence for multiple, independent association signals in 5 loci, including association with low-frequency nonsynonymous variants. Genetic loci associated with HDL-C are likely to harbor multiple, independent causative variants, frequently with opposite effects on the HDL-C phenotype. Cohorts comprising subjects at the extremes of the HDL-C distribution may be efficiently used in a case-control discovery of quantitative traits.