Genotype imputation

LinkImputeR: user-guided genotype calling and imputation for non-model organisms

Background

Genomic studies such as genome-wide association and genomic selection require genome-wide genotype data. All existing technologies used to create these data result in missing genotypes, which are often then inferred using genotype imputation software. However, existing imputation methods most often make use only of genotypes that are successfully inferred after having passed a certain read depth threshold. Because of this, any read information for genotypes that did not pass the threshold, and were thus set to missing, is ignored. Most genomic studies also choose read depth thresholds and quality filters without investigating their effects on the size and quality of the resulting genotype data. Moreover, almost all genotype imputation methods require ordered markers and are therefore of limited utility in non-model organisms.

Results

Here we introduce LinkImputeR, a software program that exploits the read count information that is normally ignored, and makes use of all available DNA sequence information for the purposes of genotype calling and imputation. It is specifically designed for non-model organisms since it requires neither ordered markers nor a reference panel of genotypes. Using next-generation DNA sequence (NGS) data from apple, cannabis and grape, we quantify the effect of varying read count and missingness thresholds on the quantity and quality of genotypes generated from LinkImputeR. We demonstrate that LinkImputeR can increase the number of genotype calls by more than an order of magnitude, can improve genotyping accuracy by several percent and can thus improve the power of downstream analyses. Moreover, we show that the effects of quality and read depth filters can differ substantially between data sets and should therefore be investigated on a per-study basis.

Conclusions

By exploiting DNA sequence data that is normally ignored during genotype calling and imputation, LinkImputeR can significantly improve both the quantity and quality of genotype data generated from NGS technologies. It enables the user to quickly and easily examine the effects of varying thresholds and filters on the number and quality of the resulting genotype calls. In this manner, users can decide on thresholds that are most suitable for their purposes. We show that LinkImputeR can significantly augment the value and utility of NGS data sets, especially in non-model organisms with poor genomic resources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3873-5) contains supplementary material, which is available to authorized users.

Influence of common and rare genetic variation on warfarin dose among African-Americans and European-Americans using the exome array

AIM

We conducted a genome-wide association study using the Illumina Exome Array to identify coding SNPs that may explain additional warfarin dose variability.

PATIENTS & METHODS

Analysis was performed after adjustment for clinical variables and genetic factors known to influence warfarin dose among 1680 warfarin users (838 European-Americans and 842 African-Americans). Replication was performed in an independent sample.

RESULTS

We confirmed the influence of known genetic variants on warfarin dose variability. Our study is the first to show the association between rs12772169 and warfarin dose in African-Americans. In addition, genes COX15 and FGF5 showed significant association in European-Americans.

CONCLUSION

We identified some novel genes/SNPs that underpin warfarin dose response. Further replication is needed to confirm our findings.

SEQSpark: A Complete Analysis Tool for Large-Scale Rare Variant Association Studies Using Whole-Genome and Exome Sequence Data

Massively parallel sequencing technologies provide great opportunities for discovering rare susceptibility variants involved in complex disease etiology via large-scale imputation and exome and whole-genome sequence-based association studies. Due to modest effect sizes, large sample sizes of tens to hundreds of thousands of individuals are required for adequately powered studies. Current analytical tools are obsolete when it comes to handling these large datasets. To facilitate the analysis of large-scale sequence-based studies, we developed SEQSpark which implements parallel processing based on Spark to increase the speed and efficiency of performing data quality control, annotation, and association analysis. To demonstrate the versatility and speed of SEQSpark, we analyzed whole-genome sequence data from the UK10K, testing for associations with waist-to-hip ratios. The analysis, which was completed in 1.5 hr, included loading data, annotation, principal component analysis, and single variant and rare variant aggregate association analysis of >9 million variants. For rare variant aggregate analysis, an exome-wide significant association (p < 2.5 × 10^-6) was observed with CCDC62 (SKAT-O [p = 6.89 × 10^-7], combined multivariate collapsing [p = 1.48 × 10^-6], and burden of rare variants [p = 1.48 × 10^-6]). SEQSpark was also used to analyze 50,000 simulated exomes and it required 1.75 hr for the analysis of a quantitative trait using several rare variant aggregate association methods. Additionally, the performance of SEQSpark was compared to Variant Association Tools and PLINK/SEQ. SEQSpark was always faster and in some situations computation was reduced to a hundredth of the time. SEQSpark will empower large sequence-based epidemiological studies to quickly elucidate genetic variation involved in the etiology of complex traits.

Taller height as a risk factor for venous thromboembolism: a Mendelian randomization meta-analysis

Essentials Observational data suggest taller people have a higher risk of venous thromboembolism (VTE). We used Mendelian randomization techniques to further explore this association in three studies. Risk of VTE increased by 30-40% for each 10 cm increment in height. Height was more strongly associated with deep vein thrombosis than with pulmonary embolism.

SUMMARY

Background Taller height is associated with a greater risk of venous thromboembolism (VTE). Objectives To use instrumental variable (IV) techniques (Mendelian randomization) to further explore this relationship. Methods Participants of European ancestry were included from two cohort studies (Atherosclerosis Risk in Communities [ARIC] study and Cardiovascular Health Study [CHS]) and one case-control study (Mayo Clinic VTE Study [Mayo]). We created two weighted genetic risk scores (GRSs) for height; the full GRS included 668 single-nucleotide polymorphisms (SNPs) from a previously published meta-analysis, and the restricted GRS included a subset of 362 SNPs not associated with weight independently of height. Standard logistic regression and IV models were used to estimate odds ratios (ORs) for VTE per 10-cm increment in height. ORs were pooled across the three studies by the use of inverse variance-weighted random effects meta-analysis. Results Among 9143 ARIC and 3180 CHS participants free of VTE at baseline, there were 367 and 109 incident VTE events. There were 1143 VTE cases and 1292 controls included from Mayo. The pooled ORs from non-IV models and models using the full and restricted GRSs as IVs were 1.27 (95% confidence interval [CI] 1.11-1.46), 1.34 (95% CI 1.04-1.73) and 1.45 (95% CI 1.04-2.01) per 10-cm greater height, respectively. Conclusions Taller height is associated with an increased risk of VTE in adults of European ancestry. Possible explanations for this association, including taller people having a greater venous surface area, a higher number of venous valves, or greater hydrostatic pressure, need to be explored further.

Polymorphisms in Renal Ammonia Metabolism Genes Correlate With 24-Hour Urine pH

Introduction

Urine pH is critical for net acid and solute excretion, but the genetic factors that contribute to its regulation are incompletely understood.

Methods

We tested the association of single nucleotide polymorphisms (SNPs) from 16 genes related to ammonia (NH₃) metabolism (15 biological candidates selected a priori, 1 selected from a previous genome-wide association study analysis) to that of 24-hour urine pH in 2493 individuals of European descent across 2 different cohorts using linear regression, adjusting for age, sex, and body mass index.

Results

Of 2871 total SNPs in these genes, 13 SNPs in ATP6V0A4 (a4 subunit of hydrogen− adenosine triphosphatase), SLC9A3 (sodium/hydrogen exchanger, isoform 3), and RHCG (Rhesus C glycoprotein), and 12 SNPs from insulin-like growth factor binding protein 7 (IGFBP7) had a meta-analysis P value <0.01 in the joint analysis plus a consistent direction of effect and at a least suggestive association (P < 0.1) in both cohorts. The maximal effect size (in pH units) for each additional minor allele of the identified SNPs was −0.13 for IGFBP7, −0.08 for ATP6V0A4, 0.06 for RHCG, and −0.06 for SLC9A3; SNP rs34447434 in IGFBP7 had the lowest meta-analysis P value (P = 7.1 × 10⁻⁸). After adjusting for net alkali absorption, urine pH remained suggestively associated with multiple SNPs in IGFBP, 1 SNP in ATP6V0A4, and a new SNP in GLS (phosphate-dependent glutaminase).

Discussion

Overall, these findings suggest that variants in common genes involved in ammonia metabolism may substantively contribute to basal urine pH regulation. These variations might influence the likelihood of developing disease conditions associated with altered urine pH, such as uric acid or calcium phosphate kidney stones.

A Bayesian mathematical model of motor and cognitive outcomes in Parkinson's disease

BACKGROUND

There are few established predictors of the clinical course of PD. Prognostic markers would be useful for clinical care and research.

OBJECTIVE

To identify predictors of long-term motor and cognitive outcomes and rate of progression in PD.

METHODS

Newly diagnosed PD participants were followed for 7 years in a prospective study, conducted at 55 centers in the United States and Canada. Analyses were conducted in 244 participants with complete demographic, clinical, genetic, and dopamine transporter imaging data. Machine learning dynamic Bayesian graphical models were used to identify and simulate predictors and outcomes. The outcomes rate of cognition changes are assessed by the Montreal Cognitive Assessment scores, and rate of motor changes are assessed by UPDRS part-III.

RESULTS

The most robust and consistent longitudinal predictors of cognitive function included older age, baseline Unified Parkinson's Disease Rating Scale (UPDRS) parts I and II, Schwab and England activities of daily living scale, striatal dopamine transporter binding, and SNP rs11724635 in the gene BST1. The most consistent predictor of UPDRS part III was baseline level of activities of daily living (part II). Key findings were replicated using long-term data from an independent cohort study.

CONCLUSIONS

Baseline function near the time of Parkinson's disease diagnosis, as measured by activities of daily living, is a consistent predictor of long-term motor and cognitive outcomes. Additional predictors identified may further characterize the expected course of Parkinson's disease and suggest mechanisms underlying disease progression. The prognostic model developed in this study can be used to simulate the effects of the prognostic variables on motor and cognitive outcomes, and can be replicated and refined with data from independent longitudinal studies.

Inbred Strain Variant Database (ISVdb): A Repository for Probabilistically Informed Sequence Differences Among the Collaborative Cross Strains and Their Founders

The Collaborative Cross (CC) is a panel of recently established multiparental recombinant inbred mouse strains. For the CC, as for any multiparental population (MPP), effective experimental design and analysis benefit from detailed knowledge of the genetic differences between strains. Such differences can be directly determined by sequencing, but until now whole-genome sequencing was not publicly available for individual CC strains. An alternative and complementary approach is to infer genetic differences by combining two pieces of information: probabilistic estimates of the CC haplotype mosaic from a custom genotyping array, and probabilistic variant calls from sequencing of the CC founders. The computation for this inference, especially when performed genome-wide, can be intricate and time-consuming, requiring the researcher to generate nontrivial and potentially error-prone scripts. To provide standardized, easy-to-access CC sequence information, we have developed the Inbred Strain Variant Database (ISVdb). The ISVdb provides, for all the exonic variants from the Sanger Institute mouse sequencing dataset, direct sequence information for CC founders and, critically, the imputed sequence information for CC strains. Notably, the ISVdb also: (1) provides predicted variant consequence metadata; (2) allows rapid simulation of F1 populations; and (3) preserves imputation uncertainty, which will allow imputed data to be refined in the future as additional sequencing and genotyping data are collected. The ISVdb information is housed in an SQL database and is easily accessible through a custom online interface (http://isvdb.unc.edu), reducing the analytic burden on any researcher using the CC.

Identification of the functional variant driving ORMDL3 and GSDMB expression in human chromosome 17q12-21 in primary biliary cholangitis

Numerous genome-wide association studies (GWAS) have been performed to identify susceptibility genes to various human complex diseases. However, in many cases, neither a functional variant nor a disease susceptibility gene have been clarified. Here, we show an efficient approach for identification of a functional variant in a primary biliary cholangitis (PBC)-susceptible region, chromosome 17q12-21 (ORMDL3-GSDMB-ZPBP2-IKZF3). High-density association mapping was carried out based on SNP imputation analysis by using the whole-genome sequence data from a reference panel of 1,070 Japanese individuals (1KJPN), together with genotype data from our previous GWAS (PBC patients: n = 1,389; healthy controls: n = 1,508). Among 23 single nucleotide polymorphisms (SNPs) with P < 1.0 × 10^-8, rs12946510 was identified as the functional variant that influences gene expression via alteration of Forkhead box protein O1 (FOXO1) binding affinity in vitro. Moreover, expression-quantitative trait locus (e-QTL) analyses showed that the PBC susceptibility allele of rs12946510 was significantly associated with lower endogenous expression of ORMDL3 and GSDMB in whole blood and spleen. This study not only identified the functional variant in chr.17q12-21 and its molecular mechanism through which it conferred susceptibility to PBC, but it also illustrated an efficient systematic approach for post-GWAS analysis that is applicable to other complex diseases.

Associations between RNA splicing regulatory variants of stemness-related genes and racial disparities in susceptibility to prostate cancer

Evidence suggests that cells with a stemness phenotype play a pivotal role in oncogenesis, and prostate cells exhibiting this phenotype have been identified. We used two genome-wide association study (GWAS) datasets of African descendants, from the Multiethnic/Minority Cohort Study of Diet and Cancer (MEC) and the Ghana Prostate Study, and two GWAS datasets of non-Hispanic whites, from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial and the Breast and Prostate Cancer Cohort Consortium (BPC3), to analyze the associations between genetic variants of stemness-related genes and racial disparities in susceptibility to prostate cancer. We evaluated associations of single-nucleotide polymorphisms (SNPs) in 25 stemness-related genes with prostate cancer risk in 1,609 cases and 2,550 controls of non-Hispanic whites (4,934 SNPs) and 1,144 cases and 1,116 controls of African descendants (5,448 SNPs) with correction by false discovery rate ≤0.2. We identified 32 SNPs in five genes (TP63, ALDH1A1, WNT1, MET and EGFR) that were significantly associated with prostate cancer risk, of which six SNPs in three genes (TP63, ALDH1A1 and WNT1) and eight EGFR SNPs showed heterogeneity in susceptibility between these two racial groups. In addition, 13 SNPs in MET and one in ALDH1A1 were found only in African descendants. The in silico bioinformatics analyses revealed that EGFR rs2072454 and SNPs in linkage with the identified SNPs in MET and ALDH1A1 (r²  > 0.6) were predicted to regulate RNA splicing. These variants may serve as novel biomarkers for racial disparities in prostate cancer risk.

Meta-analysis identifies five novel loci associated with endometriosis highlighting key genes involved in hormone metabolism

Endometriosis is a heritable hormone-dependent gynecological disorder, associated with severe pelvic pain and reduced fertility; however, its molecular mechanisms remain largely unknown. Here we perform a meta-analysis of 11 genome-wide association case-control data sets, totalling 17,045 endometriosis cases and 191,596 controls. In addition to replicating previously reported loci, we identify five novel loci significantly associated with endometriosis risk (P<5 × 10⁻⁸), implicating genes involved in sex steroid hormone pathways (FN1, CCDC170, ESR1, SYNE1 and FSHB). Conditional analysis identified five secondary association signals, including two at the ESR1 locus, resulting in 19 independent single nucleotide polymorphisms (SNPs) robustly associated with endometriosis, which together explain up to 5.19% of variance in endometriosis. These results highlight novel variants in or near specific genes with important roles in sex steroid hormone signalling and function, and offer unique opportunities for more targeted functional research efforts.

Endometriosis is a major cause of infertility. Molecular mechanisms underlying the disease involve genetic and environmental risk factors. In a meta-analysis of eleven GWA studies, Sapkota and colleagues identify five novel risk loci, implicating steroid sex hormone pathways in the pathogenesis.

Linkage disequilibrium matches forensic genetic records to disjoint genomic marker sets

Combining genotypes across datasets is central in facilitating advances in genetics. Data aggregation efforts often face the challenge of record matching-the identification of dataset entries that represent the same individual. We show that records can be matched across genotype datasets that have no shared markers based on linkage disequilibrium between loci appearing in different datasets. Using two datasets for the same 872 people-one with 642,563 genome-wide SNPs and the other with 13 short tandem repeats (STRs) used in forensic applications-we find that 90-98% of forensic STR records can be connected to corresponding SNP records and vice versa. Accuracy increases to 99-100% when ∼30 STRs are used. Our method expands the potential of data aggregation, but it also suggests privacy risks intrinsic in maintenance of databases containing even small numbers of markers-including databases of forensic significance.

Serum magnesium and the risk of prediabetes: a population-based cohort study

AIMS/HYPOTHESIS

Previous studies have found an association between serum magnesium and incident diabetes; however, this association may be due to reverse causation, whereby diabetes may induce urinary magnesium loss. In contrast, in prediabetes (defined as impaired fasting glucose), serum glucose levels are below the threshold for urinary magnesium wasting and, hence, unlikely to influence serum magnesium levels. Thus, to study the directionality of the association between serum magnesium levels and diabetes, we investigated its association with prediabetes. We also investigated whether magnesium-regulating genes influence diabetes risk through serum magnesium levels. Additionally, we quantified the effect of insulin resistance in the association between serum magnesium levels and diabetes risk.

METHODS

Within the population-based Rotterdam Study, we used Cox models, adjusted for age, sex, lifestyle factors, comorbidities, kidney function, serum levels of electrolytes and diuretic use, to study the association between serum magnesium and prediabetes/diabetes. In addition, we performed two mediation analyses: (1) to study if common genetic variation in eight magnesium-regulating genes influence diabetes risk through serum magnesium levels; and (2) to quantify the proportion of the effect of serum magnesium levels on diabetes that is mediated through insulin resistance (quantified by HOMA-IR).

RESULTS

A total of 8555 participants (mean age, 64.7 years; median follow-up, 5.7 years) with normal glucose levels (mean ± SD: 5.46 ± 0.58 mmol/l) at baseline were included. A 0.1 mmol/l decrease in serum magnesium level was associated with an increase in diabetes risk (HR 1.18 [95% CI 1.04, 1.33]), confirming findings from previous studies. Of interest, a similar association was found between serum magnesium levels and prediabetes risk (HR 1.12 [95% CI 1.01, 1.25]). Genetic variation in CLDN19, CNNM2, FXYD2, SLC41A2, and TRPM6 significantly influenced diabetes risk (p < 0.05), and for CNNM2, FXYD2, SLC41A2 and TRPM6 this risk was completely mediated by serum magnesium levels. We found that 29.1% of the effect of serum magnesium levels on diabetes was mediated through insulin resistance, whereas for prediabetes 13.4% was mediated through insulin resistance.

CONCLUSIONS/INTERPRETATION

Low serum magnesium levels are associated with an increased risk of prediabetes and this increased risk is similar to that of diabetes. Furthermore, common variants in magnesium-regulating genes modify diabetes risk through serum magnesium levels. Both findings support a potential causal role of magnesium in the development of diabetes, where the hypothesised pathway is partly mediated through insulin resistance.

Leveraging cell type specific regulatory regions to detect SNPs associated with tissue factor pathway inhibitor plasma levels

Tissue factor pathway inhibitor (TFPI) regulates the formation of intravascular blood clots, which manifest clinically as ischemic heart disease, ischemic stroke, and venous thromboembolism (VTE). TFPI plasma levels are heritable, but the genetics underlying TFPI plasma level variability are poorly understood. Herein we report the first genome-wide association scan (GWAS) of TFPI plasma levels, conducted in 251 individuals from five extended French-Canadian Families ascertained on VTE. To improve discovery, we also applied a hypothesis-driven (HD) GWAS approach that prioritized single nucleotide polymorphisms (SNPs) in (1) hemostasis pathway genes, and (2) vascular endothelial cell (EC) regulatory regions, which are among the highest expressers of TFPI. Our GWAS identified 131 SNPs with suggestive evidence of association (P-value < 5 × 10^-8 ), but no SNPs reached the genome-wide threshold for statistical significance. Hemostasis pathway genes were not enriched for TFPI plasma level associated SNPs (global hypothesis test P-value = 0.147), but EC regulatory regions contained more TFPI plasma level associated SNPs than expected by chance (global hypothesis test P-value = 0.046). We therefore stratified our genome-wide SNPs, prioritizing those in EC regulatory regions via stratified false discovery rate (sFDR) control, and reranked the SNPs by q-value. The minimum q-value was 0.27, and the top-ranked SNPs did not show association evidence in the MARTHA replication sample of 1,033 unrelated VTE cases. Although this study did not result in new loci for TFPI, our work lays out a strategy to utilize epigenomic data in prioritization schemes for future GWAS studies.

Improved imputation accuracy of rare and low-frequency variants using population-specific high-coverage WGS-based imputation reference panel

Genetic imputation is a cost-efficient way to improve the power and resolution of genome-wide association (GWA) studies. Current publicly accessible imputation reference panels accurately predict genotypes for common variants with minor allele frequency (MAF)≥5% and low-frequency variants (0.5≤MAF<5%) across diverse populations, but the imputation of rare variation (MAF<0.5%) is still rather limited. In the current study, we evaluate imputation accuracy achieved with reference panels from diverse populations with a population-specific high-coverage (30 ×) whole-genome sequencing (WGS) based reference panel, comprising of 2244 Estonian individuals (0.25% of adult Estonians). Although the Estonian-specific panel contains fewer haplotypes and variants, the imputation confidence and accuracy of imputed low-frequency and rare variants was significantly higher. The results indicate the utility of population-specific reference panels for human genetic studies.

Comprehensive analysis of the association between UBAC2 polymorphisms and Behçet's disease in a Japanese population

Behçet’s disease (BD) is reportedly associated with polymorphisms of the ubiquitin-associated domain containing 2 (UBAC2) gene in Turkish, Italian, and Chinese populations. Here we investigated whether UBAC2 polymorphisms were associated with BD in a Japanese population. Using data from 611 Japanese BD patients and 737 Japanese controls who participated in our previous genome-wide association study, we analyzed the 58 genotyped single-nucleotide polymorphisms (SNPs) in the region 100 kb upstream and downstream of UBAC2. We also performed imputation analysis in the region, with 562 imputed SNPs included in the statistical analyses. Association testing revealed that the T allele of rs9517723 in the lncRNA LOC107984558 was significantly associated with ocular and central nervous system (CNS) lesions and showed the strongest association under the recessive model (TT vs. CT+CC: ocular lesion, Pc = 0.0099, OR = 1.56; CNS lesion, Pc = 0.0052, OR = 3.42). Expression analysis revealed that rs9517723 TT homozygotes showed significantly increased UBAC2 expression (P < 0.05). Our findings suggest that enhanced UBAC2 expression associated with the homozygous risk allele (TT) of rs9517723 could induce overactivation of ubiquitination-related pathway, resulting in the development of ocular and CNS lesions in BD.

Genome-wide analysis identifies an african-specific variant in SEMA4D associated with body mass index

OBJECTIVE

The prevalence of obesity varies between ethnic groups. No genome-wide association study (GWAS) for body mass index (BMI) has been conducted in continental Africans.

METHODS

We performed a GWAS for BMI in 1,570 West Africans (WA). Replication was conducted in independent samples of WA (n = 1,411) and African Americans (AA) (n = 9,020).

RESULTS

We identified a novel genome-wide significant African-specific locus for BMI (SEMA4D, rs80068415; minor allele frequency = 0.008, P = 2.10 × 10-8 ). This finding was replicated in independent samples of WA (P = 0.013) and AA (P = 0.017). Individuals with obesity had higher serum SEMA4D levels compared to those without obesity (P < 0.0001), and elevated levels of serum SEMA4D were associated with increased obesity risk (OR = 4.2, P < 1 × 10-4 ). The prevalence of obesity was higher in individuals with the CT versus TT genotypes (55.6% vs. 22.9%).

CONCLUSIONS

A novel variant in SEMA4D was significantly associated with BMI. Carriers of the C allele were 4.6 BMI units heavier than carriers of the T allele (P = 0.0007). This variant is monomorphic in Europeans and Asians, highlighting the importance of studying diverse populations. While there is evidence for the involvement of SEMA4D in inflammatory processes, this study is the first to implicate SEMA4D in obesity pathophysiology.

Influence of genetic variants associated with body mass index on eating behavior in childhood

OBJECTIVE

Childhood eating behaviors are associated with body mass index (BMI). Recent genome-wide association studies have identified many single-nucleotide polymorphisms (SNPs) associated with adult and childhood BMI. This study hypothesized that these SNPs also influence eating behavior.

METHODS

In a population-based prospective cohort study among 3,031 children (mean age [standard deviation]: 4.0 [0.1] years), two weighted genetic risk scores, based on 15 childhood and 97 adult BMI SNPs, and ten individual appetite- and/or satiety-related SNPs were tested for association with food fussiness, food responsiveness, enjoyment of food, satiety responsiveness, and slowness in eating.

RESULTS

The 15 SNP-based childhood BMI genetic risk score was not associated with the eating behavior subscales. The 97 SNP-based adult BMI genetic risk score was nominally associated with satiety responsiveness (β: -0.007 standard deviation, 95% confidence interval [CI] -0.013, 0.000). Of the 10 individual SNPs, rs11030104 in BDNF and rs10733682 in LMX1B were nominally associated with satiety responsiveness (β: -0.057 standard deviation, 95% CI -0.112, -0.002).

CONCLUSIONS

These findings do not strongly support the hypothesis that BMI-associated SNPs also influence eating behavior at this age. A potential role for BMI SNPs in satiety responsiveness during childhood was observed; however, no associations with the other eating behavior subscales were found.

Risk prediction of pulmonary tuberculosis using genetic and conventional risk factors in adult Korean population

A complex interplay among host, pathogen, and environmental factors is believed to contribute to the risk of developing pulmonary tuberculosis (PTB). The lack of replication of published genome-wide association study (GWAS) findings limits the clinical utility of reported single nucleotide polymorphisms (SNPs). We conducted a GWAS using 467 PTB cases and 1,313 healthy controls obtained from two community-based cohorts in Korea. We evaluated the performance of PTB risk models based on different combinations of genetic and nongenetic factors and validated the results in an independent Korean population comprised of 179 PTB cases and 500 healthy controls. We demonstrated the polygenic nature of PTB and nongenetic factors such as age, sex, and body mass index (BMI) were strongly associated with PTB risk. None of the SNPs achieved genome-wide significance; instead, we were able to replicate the associations between PTB and ten SNPs near or in the genes, CDCA7, GBE1, GADL1, SPATA16, C6orf118, KIAA1432, DMRT2, CTR9, CCDC67, and CDH13, which may play roles in the immune and inflammatory pathways. Among the replicated SNPs, an intergenic SNP, rs9365798, located downstream of the C6orf118 gene showed the most significant association under the dominant model (OR = 1.59, 95% CI 1.32–1.92, P = 2.1×10⁻⁶). The performance of a risk model combining the effects of ten replicated SNPs and six nongenetic factors (i.e., age, sex, BMI, cigarette smoking, systolic blood pressure, and hemoglobin) were validated in the replication set (AUC = 0.80, 95% CI 0.76–0.84). The strategy of combining genetic and nongenetic risk factors ultimately resulted in better risk prediction for PTB in the adult Korean population.

Scanning indels in the 5q22.1 region and identification of the TMEM232 susceptibility gene that is associated with atopic dermatitis in the Chinese Han population

BACKGROUND

Atopic dermatitis (AD) is a chronic inflammatory skin disease. The 5q22.1 region was found to have an association with AD in our previous genome-wide association study (GWAS).

OBJECTIVE

To identify the AD susceptibility gene in 5q22.1 and observe its expression in AD tissues.

METHODS

Suggestive indels from the GWAS data were genotyped in 3013 AD patients and 5075 controls from the Chinese Han population with the SequenomMassArray system. Association, Bayesian and bioinformatics analyses were used to identify possible causal indels and genes in the 5q22.1 region. Immunohistochemistry (IHC) was performed to observe protein expression in the tissues. PLINK 1.07 software was used for all statistical analyses.

RESULTS

The genotyping and association analysis showed that six deletions and four SNPs were associated with AD (P<0.005). The rs11357450 (P_combined=7.79E-04, OR=1.39, logBayes Factor=1.29) deletion located in TMEM232 was identified to be the strongest variant. Analysis of the genetic model revealed that the dominant model best described rs11357450 (P=1.96E-03, OR=1.22; 95% CI=1.07-1.37). IHC showed that the expression of TMEM232 decreased gradually from the granular layer to the basal layer in AD, but in normal tissues, this trend was reversed. Additionally, positive cytoplasm staining was found in lymphocytes around the blood vessels in AD.

CONCLUSIONS

The study indicates that TMEM232 in the 5q22.1 region is the causal gene for AD in the Chinese Han population.

Gene-based genome-wide association study identified 19p13.3 for lean body mass

Lean body mass (LBM) is a complex trait for human health. To identify genomic loci underlying LBM, we performed a gene-based genome-wide association study of lean mass index (LMI) in 1000 unrelated Caucasian subjects, and replicated in 2283 unrelated Caucasians subjects. Gene-based association analyses highlighted the significant associations of three genes UQCR, TCF3 and MBD3 in one single locus 19p13.3 (discovery p = 6.10 × 10-5, 1.65 × 10-4 and 1.10 × 10-4; replication p = 2.21 × 10-3, 1.84 × 10-3 and 6.95 × 10-3; combined p = 2.26 × 10-6, 4.86 × 10-6 and 1.15 × 10-5, respectively). These results, together with the known functional relevance of the three genes to LMI, suggested that the 19p13.3 region containing UQCR, TCF3 and MBD3 genes was a novel locus underlying lean mass variation.

Joint genome-wide prediction in several populations accounting for randomness of genotypes: A hierarchical Bayes approach. I: Multivariate Gaussian priors for marker effects and derivation of the joint probability mass function of genotypes

It is important to consider heterogeneity of marker effects and allelic frequencies in across population genome-wide prediction studies. Moreover, all regression models used in genome-wide prediction overlook randomness of genotypes. In this study, a family of hierarchical Bayesian models to perform across population genome-wide prediction modeling genotypes as random variables and allowing population-specific effects for each marker was developed. Models shared a common structure and differed in the priors used and the assumption about residual variances (homogeneous or heterogeneous). Randomness of genotypes was accounted for by deriving the joint probability mass function of marker genotypes conditional on allelic frequencies and pedigree information. As a consequence, these models incorporated kinship and genotypic information that not only permitted to account for heterogeneity of allelic frequencies, but also to include individuals with missing genotypes at some or all loci without the need for previous imputation. This was possible because the non-observed fraction of the design matrix was treated as an unknown model parameter. For each model, a simpler version ignoring population structure, but still accounting for randomness of genotypes was proposed. Implementation of these models and computation of some criteria for model comparison were illustrated using two simulated datasets. Theoretical and computational issues along with possible applications, extensions and refinements were discussed. Some features of the models developed in this study make them promising for genome-wide prediction, the use of information contained in the probability distribution of genotypes is perhaps the most appealing. Further studies to assess the performance of the models proposed here and also to compare them with conventional models used in genome-wide prediction are needed.

No prognostic value added by vitamin D pathway SNPs to current prognostic system for melanoma survival

The prognostic improvement attributed to genetic markers over current prognostic system has not been well studied for melanoma. The goal of this study is to evaluate the added prognostic value of Vitamin D Pathway (VitD) SNPs to currently known clinical and demographic factors such as age, sex, Breslow thickness, mitosis and ulceration (CDF). We utilized two large independent well-characterized melanoma studies: the Genes, Environment, and Melanoma (GEM) and MD Anderson studies, and performed variable selection of VitD pathway SNPs and CDF using Random Survival Forest (RSF) method in addition to Cox proportional hazards models. The Harrell's C-index was used to compare the performance of model predictability. The population-based GEM study enrolled 3,578 incident cases of cutaneous melanoma (CM), and the hospital-based MD Anderson study consisted of 1,804 CM patients. Including both VitD SNPs and CDF yielded C-index of 0.85, which provided slight but not significant improvement by CDF alone (C-index = 0.83) in the GEM study. Similar results were observed in the independent MD Anderson study (C-index = 0.84 and 0.83, respectively). The Cox model identified no significant associations after adjusting for multiplicity. Our results do not support clinically significant prognostic improvements attributable to VitD pathway SNPs over current prognostic system for melanoma survival.

Association of genetic variation in the tachykinin receptor 3 locus with hot flashes and night sweats in the Women's Health Initiative Study

OBJECTIVE

Vasomotor symptoms (VMS, ie, hot flashes or night sweats) are reported by many, but not all, women. The extent to which VMS are genetically determined is unknown. We evaluated the relationship of genetic variation and VMS.

METHODS

In this observational study, we accessed data from three genome-wide association studies (GWAS) (SNP Health Association Resource cohort [SHARe], WHI Memory Study cohort [WHIMS+], and Genome-Wide Association Studies of Treatment Response in Randomized Clinical Trials [GARNET] studies, total n = 17,695) of European American, African American, and Hispanic American postmenopausal women aged 50 to 79 years at baseline in the Women's Health Initiative Study. We examined genetic variation in relation to VMS (yes/no) in each study and using trans-ethnic inverse variance fixed-effects meta-analysis. A total of 11,078,977 single-nucleotide polymorphisms (SNPs) met the quality criteria.

RESULTS

After adjustment for covariates and population structure, three SNPs (on chromosomes 3 and 11) were associated with VMS at the genome-wide threshold of 5 × 10 in the African American SHARe GWAS, but were not associated in the other cohorts. In the meta-analysis, 14 SNPs, all located on chromosome 4 in the tachykinin receptor 3 (TACR3) locus, however, had P < 5 × 10. These SNPs' effect sizes were similar across studies/participants' ancestry (odds ratio ∼1.5).

CONCLUSIONS

Genetic variation in TACR3 may contribute to the risk of VMS. To our knowledge, this is the first GWAS to examine SNPs associated with VMS. These results support the biological hypothesis of a role for TACR3 in VMS, which was previously hypothesized from animal and human studies. Further study of these variants may lead to new insights into the biological pathways involved in VMS, which are poorly understood.

Scavenger Receptor A1 Prevents Metastasis of Non-Small Cell Lung Cancer via Suppression of Macrophage Serum Amyloid A1

Mechanisms of cross-talk between tumor cells and tumor-associated macrophages (TAM), which drive metastasis, are not fully understood. Scavenger receptor A1 (SR-A1) expressed primarily in macrophages has been associated with lung tumorigenesis. In this study, we used population genetics, transcriptomics, and functional analyses to uncover how SR-A1 is involved in lung cancer and its prognosis. SR-A1 genetic variants were investigated for possible association with survival of advanced stage NSCLC patients in the Harvard Lung Cancer Study cohort. Two SNPs (rs17484273, rs1484751) in SR-A1 were associated significantly with poor overall survival in this cohort. Data from The Cancer Genome Atlas showed considerable downregulation of SR-A1 in lung tumor tissues. The association of SR-A1 with prognosis was validated in animal models in the context of lung cancer metastasis. Macrophages derived from mice genetically deficient for SR-A1 exhibited accelerated metastasis in a model of lung cancer. On the other hand, tumor cell seeding, migration, and invasion, as well as macrophage accumulation in lung cancer tissue, were enhanced in SR-A1-deficient mice. SR-A1 deletion upregulated serum amyloid A1 (SAA1) in macrophages via MAPK/IκB/NFκB signaling. SAA1 promoted tumor cell invasion and macrophage migration in vitro and in vivo, but these effects were blocked by administration of an anti-SAA1 antibody. Overall, our findings show how SR-A1 suppresses lung cancer metastasis by downregulating SAA1 production in TAMs. Cancer Res; 77(7); 1586-98. ©2017 AACR.

Quantifying Interdependent Risks in Genomic Privacy

The rapid progress in human-genome sequencing is leading to a high availability of genomic data. These data is notoriously very sensitive and stable in time, and highly correlated among relatives. In this article, we study the implications of these familial correlations on kin genomic privacy. We formalize the problem and detail efficient reconstruction attacks based on graphical models and belief propagation. With our approach, an attacker can infer the genomes of the relatives of an individual whose genome or phenotype are observed by notably relying on Mendel’s Laws, statistical relationships between the genomic variants, and between the genome and the phenotype. We evaluate the effect of these dependencies on privacy with respect to the amount of observed variants and the relatives sharing them. We also study how the algorithmic performance evolves when we take these various relationships into account. Furthermore, to quantify the level of genomic privacy as a result of the proposed inference attack, we discuss possible definitions of genomic privacy metrics, and compare their values and evolution. Genomic data reveals Mendelian disorders and the likelihood of developing severe diseases, such as Alzheimer’s. We also introduce the quantification of health privacy , specifically, the measure of how well the predisposition to a disease is concealed from an attacker. We evaluate our approach on actual genomic data from a pedigree and show the threat extent by combining data gathered from a genome-sharing website as well as an online social network.

The Rare-Variant Generalized Disequilibrium Test for Association Analysis of Nuclear and Extended Pedigrees with Application to Alzheimer Disease WGS Data

Whole-genome and exome sequence data can be cost-effectively generated for the detection of rare-variant (RV) associations in families. Causal variants that aggregate in families usually have larger effect sizes than those found in sporadic cases, so family-based designs can be a more powerful approach than population-based designs. Moreover, some family-based designs are robust to confounding due to population admixture or substructure. We developed a RV extension of the generalized disequilibrium test (GDT) to analyze sequence data obtained from nuclear and extended families. The GDT utilizes genotype differences of all discordant relative pairs to assess associations within a family, and the RV extension combines the single-variant GDT statistic over a genomic region of interest. The RV-GDT has increased power by efficiently incorporating information beyond first-degree relatives and allows for the inclusion of covariates. Using simulated genetic data, we demonstrated that the RV-GDT method has well-controlled type I error rates, even when applied to admixed populations and populations with substructure. It is more powerful than existing family-based RV association methods, particularly for the analysis of extended pedigrees and pedigrees with missing data. We analyzed whole-genome sequence data from families affected by Alzheimer disease to illustrate the application of the RV-GDT. Given the capability of the RV-GDT to adequately control for population admixture or substructure and analyze pedigrees with missing genotype data and its superior power over other family-based methods, it is an effective tool for elucidating the involvement of RVs in the etiology of complex traits.

Securing the use of existing sample collections for future human genetic research

Hundreds of thousands of individuals have been genotyped in the past decades using genotyping arrays, representing both a valuable data resource for future biomedical research and a substantial investment in human genetic research. However, novel chip designs and their altered sets of single-nucleotide polymorphisms (SNPs) pose the question of how well established data resources, such as large samples of healthy controls genotyped on legacy arrays, can be combined with newer samples genotyped on those novel arrays using genotype imputation. We exemplarily investigated this question based on genotype data of 30 European and 30 African unrelated samples from the 1000 Genomes project and on markers present on two legacy SNP arrays, namely Affymetrix's Human SNP 6.0 and Illumina's 550k array, and three newer arrays, namely two Axiom arrays from Affymetrix and an OmniExpress array from Illumina. We cross-compared the imputation accuracy as well as efficacy and assessed genotype concordance among these arrays. Although the accuracy of genotype prediction was uniformly high across all arrays, the imputation efficacy, that is, the proportion of successfully imputed markers, differed considerably between array combinations in both sample sets, with legacy arrays showing a trend towards lower efficacy values compared with newer arrays when serving as imputation basis. We conclude that, given the substantial losses of markers covered by the legacy arrays, the re-genotyping of existing samples sets, in particular those of healthy population controls, would be a worthwhile endeavor to secure their continued use in the future.

A comparison of different algorithms for phasing haplotypes using Holstein cattle genotypes and pedigree data

Phasing genotypes to haplotypes is becoming increasingly important due to its applications in the study of diseases, population and evolutionary genetics, imputation, and so on. Several studies have focused on the development of computational methods that infer haplotype phase from population genotype data. The aim of this study was to compare phasing algorithms implemented in Beagle, Findhap, FImpute, Impute2, and ShapeIt2 software using 50k and 777k (HD) genotyping data. Six scenarios were considered: no-parents, sire-progeny pairs, sire-dam-progeny trios, each with and without pedigree information in Holstein cattle. Algorithms were compared with respect to their phasing accuracy and computational efficiency. In the studied population, Beagle and FImpute were more accurate than other phasing algorithms. Across scenarios, phasing accuracies for Beagle and FImpute were 99.49-99.90% and 99.44-99.99% for 50k, respectively, and 99.90-99.99% and 99.87-99.99% for HD, respectively. Generally, FImpute resulted in higher accuracy when genotypic information of at least one parent was available. In the absence of parental genotypes and pedigree information, Beagle and Impute2 (with double the default number of states) were slightly more accurate than FImpute. Findhap gave high phasing accuracy when parents' genotypes and pedigree information were available. In terms of computing time, Findhap was the fastest algorithm followed by FImpute. FImpute was 30 to 131, 87 to 786, and 353 to 1,400 times faster across scenarios than Beagle, ShapeIt2, and Impute2, respectively. In summary, FImpute and Beagle were the most accurate phasing algorithms. Moreover, the low computational requirement of FImpute makes it an attractive algorithm for phasing genotypes of large livestock populations.

Genetic Basis of Chronotype in Humans: Insights From Three Landmark GWAS

Study Objectives

Chronotype, or diurnal preference, refers to behavioral manifestations of the endogenous circadian system that governs preferred timing of sleep and wake. As variations in circadian timing and system perturbations are linked to disease development, the fundamental biology of chronotype has received attention for its role in the regulation and dysregulation of sleep and related illnesses. Family studies indicate that chronotype is a heritable trait, thus directing attention toward its genetic basis. Although discoveries from molecular studies of candidate genes have shed light onto its genetic architecture, the contribution of genetic variation to chronotype has remained unclear with few related variants identified. In the advent of large-scale genome-wide association studies (GWAS), scientists now have the ability to discover novel common genetic variants associated with complex phenotypes. Three recent large-scale GWASs of chronotype were conducted on subjects of European ancestry from the 23andMe cohort and the UK Biobank. This review discusses the findings of these landmark GWASs in the context of prior research.

Methods

We systematically reviewed and compared methodological and analytical approaches and results across the three GWASs of chronotype.

Results

A good deal of consistency was observed across studies with 9 genes identified in 2 of the 3 GWASs. Several genes previously unknown to influence chronotype were identified.

Conclusions

GWAS is an important tool in identifying common variants associated with the complex chronotype phenotype, the findings of which can supplement and guide molecular science. Future directions in model systems and discovery of rare variants are discussed.

Comparison of HapMap and 1000 Genomes Reference Panels in a Large-Scale Genome-Wide Association Study

An increasing number of genome-wide association (GWA) studies are now using the higher resolution 1000 Genomes Project reference panel (1000G) for imputation, with the expectation that 1000G imputation will lead to the discovery of additional associated loci when compared to HapMap imputation. In order to assess the improvement of 1000G over HapMap imputation in identifying associated loci, we compared the results of GWA studies of circulating fibrinogen based on the two reference panels. Using both HapMap and 1000G imputation we performed a meta-analysis of 22 studies comprising the same 91,953 individuals. We identified six additional signals using 1000G imputation, while 29 loci were associated using both HapMap and 1000G imputation. One locus identified using HapMap imputation was not significant using 1000G imputation. The genome-wide significance threshold of 5×10-8 is based on the number of independent statistical tests using HapMap imputation, and 1000G imputation may lead to further independent tests that should be corrected for. When using a stricter Bonferroni correction for the 1000G GWA study (P-value < 2.5×10-8), the number of loci significant only using HapMap imputation increased to 4 while the number of loci significant only using 1000G decreased to 5. In conclusion, 1000G imputation enabled the identification of 20% more loci than HapMap imputation, although the advantage of 1000G imputation became less clear when a stricter Bonferroni correction was used. More generally, our results provide insights that are applicable to the implementation of other dense reference panels that are under development.

Genome-wide association test of multiple continuous traits using imputed SNPs

More and more large cohort studies have conducted or are conducting genome-wide association studies (GWAS) to reveal the genetic components of many complex human diseases. These large cohort studies often collected a broad array of correlated phenotypes that reflect common physiological processes. By jointly analyzing these correlated traits, we can gain more power by aggregating multiple weak effects and shed light on the mechanisms underlying complex human diseases. The majority of existing multi-trait association test methods are based on jointly modeling the multivariate traits conditional on the genotype as covariate, and can readily accommodate the imputed SNPs by using their imputed dosage as a covariate. An alternative class of multi-trait association tests is based on the inverted regression, which models the distribution of genotypes conditional on the covariate and multivariate traits, and has been shown to have competitive performance. To our knowledge, all existing inverted regression approaches have implicitly used the "best-guess" genotypes, which is not efficient and known to lead to dramatic power loss, and there have not been any proposed methods of incorporating imputation uncertainty into inverted regressions. In this work, we propose a general and efficient framework that can account for the imputation uncertainty to further improve the association test power of inverted regression models for imputed SNPs. We demonstrate through extensive numerical studies that the proposed method has competitive performance. We further illustrate its usefulness by application to association test of diabetes-related glycemic traits in the Atherosclerosis Risk in Communities (ARIC) Study.

The Analysis of Ethnic Mixtures

Population of ethnic mixtures can be useful in genetic studies. Admixture mapping, or mapping by admixture linkage disequilibrium (MALD), is specially developed for admixed populations and can supplement traditional genome-wide association analyses in the search for genetic variants underlying complex traits. Admixture mapping tests the association between a trait and locus-specific ancestries. The locus-specific ancestries are in linkage disequilibrium (LD), which is generated by an admixture process between genetically distinct ancestral populations. Because of the highly correlated-locus specific ancestries, admixture mapping performs many fewer independent tests across the genome than current genome-wide association analysis. Therefore, admixture mapping can be more powerful because it reduces the penalty due to multiple tests. In this chapter, we introduce the theory behind admixture mapping and explain how to conduct the analysis in practice.

Using Genetic Variation to Explore the Causal Effect of Maternal Pregnancy Adiposity on Future Offspring Adiposity: A Mendelian Randomisation Study

BACKGROUND

It has been suggested that greater maternal adiposity during pregnancy affects lifelong risk of offspring fatness via intrauterine mechanisms. Our aim was to use Mendelian randomisation (MR) to investigate the causal effect of intrauterine exposure to greater maternal body mass index (BMI) on offspring BMI and fat mass from childhood to early adulthood.

METHODS AND FINDINGS

We used maternal genetic variants as instrumental variables (IVs) to test the causal effect of maternal BMI in pregnancy on offspring fatness (BMI and dual-energy X-ray absorptiometry [DXA] determined fat mass index [FMI]) in a MR approach. This was investigated, with repeat measurements, from ages 7 to 18 in the Avon Longitudinal Study of Parents and Children (ALSPAC; n = 2,521 to 3,720 for different ages). We then sought to replicate findings with results for BMI at age 6 in Generation R (n = 2,337 for replication sample; n = 6,057 for total pooled sample). In confounder-adjusted multivariable regression in ALSPAC, a 1 standard deviation (SD, equivalent of 3.7 kg/m2) increase in maternal BMI was associated with a 0.25 SD (95% CI 0.21-0.29) increase in offspring BMI at age 7, with similar results at later ages and when FMI was used as the outcome. A weighted genetic risk score was generated from 32 genetic variants robustly associated with BMI (minimum F-statistic = 45 in ALSPAC). The MR results using this genetic risk score as an IV in ALSPAC were close to the null at all ages (e.g., 0.04 SD (95% CI -0.21-0.30) at age 7 and 0.03 SD (95% CI -0.26-0.32) at age 18 per SD increase in maternal BMI), which was similar when a 97 variant generic risk score was used in ALSPAC. When findings from age 7 in ALSPAC were meta-analysed with those from age 6 in Generation R, the pooled confounder-adjusted multivariable regression association was 0.22 SD (95% CI 0.19-0.25) per SD increase in maternal BMI and the pooled MR effect (pooling the 97 variant score results from ALSPAC with the 32 variant score results from Generation R) was 0.05 SD (95%CI -0.11-0.21) per SD increase in maternal BMI (p-value for difference between the two results = 0.05). A number of sensitivity analyses exploring violation of the MR results supported our main findings. However, power was limited for some of the sensitivity tests and further studies with relevant data on maternal, offspring, and paternal genotype are required to obtain more precise (and unbiased) causal estimates.

CONCLUSIONS

Our findings provide little evidence to support a strong causal intrauterine effect of incrementally greater maternal BMI resulting in greater offspring adiposity.

Integration of Population-Level Genotype Data with Functional Annotation Reveals Over-Representation of Long Noncoding RNAs at Ovarian Cancer Susceptibility Loci

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple loci associated with epithelial ovarian cancer (EOC) susceptibility, but further progress requires integration of epidemiology and biology to illuminate true risk loci below genome-wide significance levels (P < 5 × 10-8). Most risk SNPs lie within non-protein-encoding regions, and we hypothesize that long noncoding RNA (lncRNA) genes are enriched at EOC risk regions and represent biologically relevant functional targets.

METHODS

Using imputed GWAS data from about 18,000 invasive EOC cases and 34,000 controls of European ancestry, the GENCODE (v19) lncRNA database was used to annotate SNPs from 13,442 lncRNAs for permutation-based enrichment analysis. Tumor expression quantitative trait locus (eQTL) analysis was performed for sub-genome-wide regions (1 × 10-5 > P > 5 × 10-8) overlapping lncRNAs.

RESULTS

Of 5,294 EOC-associated SNPs (P < 1.0 × 10-5), 1,464 (28%) mapped within 53 unique lncRNAs and an additional 3,484 (66%) SNPs were correlated (r2 > 0.2) with SNPs within 115 lncRNAs. EOC-associated SNPs comprised 130 independent regions, of which 72 (55%) overlapped with lncRNAs, representing a significant enrichment (P = 5.0 × 10-4) that was more pronounced among a subset of 5,401 lncRNAs with active epigenetic regulation in normal ovarian tissue. EOC-associated lncRNAs and their putative promoters and transcription factors were enriched for biologically relevant pathways and eQTL analysis identified five novel putative risk regions with allele-specific effects on lncRNA gene expression.

CONCLUSIONS

lncRNAs are significantly enriched at EOC risk regions, suggesting a mechanistic role for lncRNAs in driving predisposition to EOC.

IMPACT

lncRNAs represent key candidates for integrative epidemiologic and functional studies. Further research on their biologic role in ovarian cancer is indicated. Cancer Epidemiol Biomarkers Prev; 26(1); 116-25. ©2016 AACR.

Age at menarche and age at natural menopause in East Asian women: a genome-wide association study

Age at menarche (AM) and age at natural menopause (ANM) are complex traits with a high heritability. Abnormal timing of menarche or menopause is associated with a reduced span of fertility and risk for several age-related diseases including breast, endometrial and ovarian cancer, cardiovascular disease, and osteoporosis. To identify novel genetic loci for AM or ANM in East Asian women and to replicate previously identified loci primarily in women of European ancestry by genome-wide association studies (GWASs), we conducted a two-stage GWAS. Stage I aimed to discover promising novel AM and ANM loci using GWAS data of 8073 women from Shanghai, China. The Stage II replication study used the data from another Chinese GWAS (n = 1230 for AM and n = 1458 for ANM), a Korean GWAS (n = 4215 for AM and n = 1739 for ANM), and de novo genotyping of 2877 additional Chinese women. Previous GWAS-identified loci for AM and ANM were also evaluated. We identified two suggestive menarcheal age loci tagged by rs79195475 at 10q21.3 (beta = -0.118 years, P = 3.4 × 10^-6) and rs1023935 at 4p15.1 (beta = -0.145 years, P = 4.9 × 10^-6) and one menopausal age locus tagged by rs3818134 at 22q12.2 (beta = -0.276 years, P = 8.8 × 10^-6). These suggestive loci warrant a further validation in independent populations. Although limited by low statistical power, we replicated 19 of the 98 menarche loci and 5 of the 20 menopause loci previously identified in women of European ancestry in East Asian women, suggesting a shared genetic architecture for these two traits across populations.

The Future is The Past: Methylation QTLs in Schizophrenia

Genome-wide association studies (GWAS) have remarkably advanced insight into the genetic basis of schizophrenia (SCZ). Still, most of the functional variance in disease risk remains unexplained. Hence, there is a growing need to map genetic variability-to-genes-to-functions for understanding the pathophysiology of SCZ and the development of better treatments. Genetic variation can regulate various cellular functions including DNA methylation, an epigenetic mark with important roles in transcription and the mediation of environmental influences. Methylation quantitative trait loci (meQTLs) are derived by mapping levels of DNA methylation in genetically different, genotyped individuals and define loci at which DNA methylation is influenced by genetic variation. Recent evidence points to an abundance of meQTLs in brain tissues whose functional contributions to development and mental diseases are still poorly understood. Interestingly, fetal meQTLs reside in regulatory domains affecting methylome reconfiguration during early brain development and are enriched in loci identified by GWAS for SCZ. Moreover, fetal meQTLs are preserved in the adult brain and could trace early epigenomic deregulation during vulnerable periods. Overall, these findings highlight the role of fetal meQTLs in the genetic risk for and in the possible neurodevelopmental origin of SCZ.

Radiogenomics: A systems biology approach to understanding genetic risk factors for radiotherapy toxicity?

Adverse reactions in normal tissue after radiotherapy (RT) limit the dose that can be given to tumour cells. Since 80% of individual variation in clinical response is estimated to be caused by patient-related factors, identifying these factors might allow prediction of patients with increased risk of developing severe reactions. While inactivation of cell renewal is considered a major cause of toxicity in early-reacting normal tissues, complex interactions involving multiple cell types, cytokines, and hypoxia seem important for late reactions. Here, we review 'omics' approaches such as screening of genetic polymorphisms or gene expression analysis, and assess the potential of epigenetic factors, posttranslational modification, signal transduction, and metabolism. Furthermore, functional assays have suggested possible associations with clinical risk of adverse reaction. Pathway analysis incorporating different 'omics' approaches may be more efficient in identifying critical pathways than pathway analysis based on single 'omics' data sets. Integrating these pathways with functional assays may be powerful in identifying multiple subgroups of RT patients characterised by different mechanisms. Thus 'omics' and functional approaches may synergise if they are integrated into radiogenomics 'systems biology' to facilitate the goal of individualised radiotherapy.

Associations of genetic variants for adult lipid levels with lipid levels in children. The Generation R Study

Lipid concentrations are heritable traits. Recently, the number of known genetic loci associated with lipid levels in adults increased from 95 to 157. The effects of these 157 loci have not been tested in children. Considering that lipid levels track from childhood to adulthood, we studied to determine whether these variants already affected lipid concentrations in a large group of 2,645 children with a median age of 6.0 years (95% range 5.7-7.3 years) from the population-based Generation R Study. Twenty-eight SNPs associated with TGs, 39 SNPs associated with total cholesterol (TC), 28 SNPs associated with LDL cholesterol (LDL-C), and 56 SNPs associated with HDL cholesterol (HDL-C) were analyzed individually and combined into genetic risk scores (GRSs). All risk scores were associated with their specific outcomes. The differences in mean absolute lipid and lipoprotein values between the 10% of children with the highest lipid or lipoprotein GRS versus the 10% with the lowest score were 0.28, 0.25, 0.32, and 0.30 mmol/l for TGs, TC, LDL-C, and HDL-C, respectively. In conclusion, we show for the first time that GRSs based on 157 SNPs associated with adult lipid concentrations are associated with lipid levels in children. The genetic background of these phenotypes at least partly overlaps between children and adults.

Phenotypically Enriched Genotypic Imputation in Genetic Association Tests

BACKGROUND

In longitudinal epidemiological studies there may be individuals with rich phenotype data who die or are lost to follow-up before providing DNA for genetic studies. Often, the genotypic and phenotypic data of the relatives are available. Two strategies for analyzing the incomplete data are to exclude ungenotyped subjects from analysis (the complete-case method, CC) and to include phenotyped but ungenotyped individuals in analysis by using relatives' genotypes for genotype imputation (GI). In both strategies, the information in the phenotypic data was not used to handle the missing-genotype problem.

METHODS

We propose a phenotypically enriched genotypic imputation (PEGI) method that uses the EM (expectation-maximization)-based maximum likelihood method to incorporate observed phenotypes into genotype imputation.

RESULTS

Our simulations with genotypes missing completely at random show that, for a single-nucleotide polymorphism (SNP) with moderate to strong effect on a phenotype, PEGI improves power more than GI without excess type I errors. Using the Framingham Heart Study data set, we compare the ability of the PEGI, GI, and CC to detect the associations between 5 SNPs and age at natural menopause.

CONCLUSION

The PEGI method may improve power to detect an association over both CC and GI under many circumstances.

A genome-wide association study in multiple system atrophy

OBJECTIVE

To identify genetic variants that play a role in the pathogenesis of multiple system atrophy (MSA), we undertook a genome-wide association study (GWAS).

METHODS

We performed a GWAS with >5 million genotyped and imputed single nucleotide polymorphisms (SNPs) in 918 patients with MSA of European ancestry and 3,864 controls. MSA cases were collected from North American and European centers, one third of which were neuropathologically confirmed.

RESULTS

We found no significant loci after stringent multiple testing correction. A number of regions emerged as potentially interesting for follow-up at p < 1 × 10-6, including SNPs in the genes FBXO47, ELOVL7, EDN1, and MAPT. Contrary to previous reports, we found no association of the genes SNCA and COQ2 with MSA.

CONCLUSIONS

We present a GWAS in MSA. We have identified several potentially interesting gene loci, including the MAPT locus, whose significance will have to be evaluated in a larger sample set. Common genetic variation in SNCA and COQ2 does not seem to be associated with MSA. In the future, additional samples of well-characterized patients with MSA will need to be collected to perform a larger MSA GWAS, but this initial study forms the basis for these next steps.

Comparing performance of modern genotype imputation methods in different ethnicities

A variety of modern software packages are available for genotype imputation relying on advanced concepts such as pre-phasing of the target dataset or utilization of admixed reference panels. In this study, we performed a comprehensive evaluation of the accuracy of modern imputation methods on the basis of the publicly available POPRES samples. Good quality genotypes were masked and re-imputed by different imputation frameworks: namely MaCH, IMPUTE2, MaCH-Minimac, SHAPEIT-IMPUTE2 and MaCH-Admix. Results were compared to evaluate the relative merit of pre-phasing and the usage of admixed references. We showed that the pre-phasing framework SHAPEIT-IMPUTE2 can overestimate the certainty of genotype distributions resulting in the lowest percentage of correctly imputed genotypes in our case. MaCH-Minimac performed better than SHAPEIT-IMPUTE2. Pre-phasing always reduced imputation accuracy. IMPUTE2 and MaCH-Admix, both relying on admixed-reference panels, showed comparable results. MaCH showed superior results if well-matched references were available (Nei’s G_ST ≤ 0.010). For small to medium datasets, frameworks using genetically closest reference panel are recommended if the genetic distance between target and reference data set is small. Our results are valid for small to medium data sets. As shown on a larger data set of population based German samples, the disadvantage of pre-phasing decreases for larger sample sizes.

An Adaptive Fisher's Combination Method for Joint Analysis of Multiple Phenotypes in Association Studies

Currently, the analyses of most genome-wide association studies (GWAS) have been performed on a single phenotype. There is increasing evidence showing that pleiotropy is a widespread phenomenon in complex diseases. Therefore, using only one single phenotype may lose statistical power to identify the underlying genetic mechanism. There is an increasing need to develop and apply powerful statistical tests to detect association between multiple phenotypes and a genetic variant. In this paper, we develop an Adaptive Fisher’s Combination (AFC) method for joint analysis of multiple phenotypes in association studies. The AFC method combines p-values obtained in standard univariate GWAS by using the optimal number of p-values which is determined by the data. We perform extensive simulations to evaluate the performance of the AFC method and compare the power of our method with the powers of TATES, Tippett’s method, Fisher’s combination test, MANOVA, MultiPhen, and SUMSCORE. Our simulation studies show that the proposed method has correct type I error rates and is either the most powerful test or comparable with the most powerful test. Finally, we illustrate our proposed methodology by analyzing whole-genome genotyping data from a lung function study.

Genetic variants in RBFOX3 are associated with sleep latency

Time to fall asleep (sleep latency) is a major determinant of sleep quality. Chronic, long sleep latency is a major characteristic of sleep-onset insomnia and/or delayed sleep phase syndrome. In this study we aimed to discover common polymorphisms that contribute to the genetics of sleep latency. We performed a meta-analysis of genome-wide association studies (GWAS) including 2 572 737 single nucleotide polymorphisms (SNPs) established in seven European cohorts including 4242 individuals. We found a cluster of three highly correlated variants (rs9900428, rs9907432 and rs7211029) in the RNA-binding protein fox-1 homolog 3 gene (RBFOX3) associated with sleep latency (P-values=5.77 × 10(-08), 6.59 × 10(-)(08) and 9.17 × 10(-)(08)). These SNPs were replicated in up to 12 independent populations including 30 377 individuals (P-values=1.5 × 10(-)(02), 7.0 × 10(-)(03) and 2.5 × 10(-)(03); combined meta-analysis P-values=5.5 × 10(-07), 5.4 × 10(-07) and 1.0 × 10(-07)). A functional prediction of RBFOX3 based on co-expression with other genes shows that this gene is predominantly expressed in brain (P-value=1.4 × 10(-316)) and the central nervous system (P-value=7.5 × 10(-)(321)). The predicted function of RBFOX3 based on co-expression analysis with other genes shows that this gene is significantly involved in the release cycle of neurotransmitters including gamma-aminobutyric acid and various monoamines (P-values<2.9 × 10(-11)) that are crucial in triggering the onset of sleep. To conclude, in this first large-scale GWAS of sleep latency we report a novel association of variants in RBFOX3 gene. Further, a functional prediction of RBFOX3 supports the involvement of RBFOX3 with sleep latency.

Single nucleotide polymorphisms in an intergenic chromosome 2q region associated with tissue factor pathway inhibitor plasma levels and venous thromboembolism

Essentials Tissue factor pathway inhibitor (TFPI) regulates the blood coagulation cascade. We replicated previously reported linkage of TFPI plasma levels to the chromosome 2q region. The putative causal locus, rs62187992, was associated with TFPI plasma levels and thrombosis. rs62187992 was marginally associated with TFPI expression in human aortic endothelial cells. Click to hear Ann Gil's presentation on new insights into thrombin activatable fibrinolysis inhibitor SUMMARY: Background Tissue factor pathway inhibitor (TFPI) regulates fibrin clot formation, and low TFPI plasma levels increase the risk of arterial thromboembolism and venous thromboembolism (VTE). TFPI plasma levels are also heritable, and a previous linkage scan implicated the chromosome 2q region, but no specific genes. Objectives To replicate the finding of the linkage region in an independent sample, and to identify the causal locus. Methods We first performed a linkage analysis of microsatellite markers and TFPI plasma levels in 251 individuals from the F5L Family Study, and replicated the finding of the linkage peak on chromosome 2q (LOD = 3.06). We next defined a follow-up region that included 112 603 single nucleotide polymorphisms (SNPs) under the linkage peak, and meta-analyzed associations between these SNPs and TFPI plasma levels across the F5L Family Study and the Marseille Thrombosis Association (MARTHA) Study, a study of 1033 unrelated VTE patients. SNPs with false discovery rate q-values of < 0.10 were tested for association with TFPI plasma levels in 892 patients with coronary artery disease in the AtheroGene Study. Results and Conclusions One SNP, rs62187992, was associated with TFPI plasma levels in all three samples (β = + 0.14 and P = 4.23 × 10-6 combined; β = + 0.16 and P = 0.02 in the F5L Family Study; β = + 0.13 and P = 6.3 × 10-4 in the MARTHA Study; β = + 0.17 and P = 0.03 in the AtheroGene Study), and contributed to the linkage peak in the F5L Family Study. rs62187992 was also associated with clinical VTE (odds ratio 0.90, P = 0.03) in the INVENT Consortium of > 7000 cases and their controls, and was marginally associated with TFPI expression (β = + 0.19, P = 0.08) in human aortic endothelial cells, a primary site of TFPI synthesis. The biological mechanisms underlying these associations remain to be elucidated.

Generalized reduced rank latent factor regression for high dimensional tensor fields, and neuroimaging-genetic applications

We propose a generalized reduced rank latent factor regression model (GRRLF) for the analysis of tensor field responses and high dimensional covariates. The model is motivated by the need from imaging-genetic studies to identify genetic variants that are associated with brain imaging phenotypes, often in the form of high dimensional tensor fields. GRRLF identifies from the structure in the data the effective dimensionality of the data, and then jointly performs dimension reduction of the covariates, dynamic identification of latent factors, and nonparametric estimation of both covariate and latent response fields. After accounting for the latent and covariate effects, GRLLF performs a nonparametric test on the remaining factor of interest. GRRLF provides a better factorization of the signals compared with common solutions, and is less susceptible to overfitting because it exploits the effective dimensionality. The generality and the flexibility of GRRLF also allow various statistical models to be handled in a unified framework and solutions can be efficiently computed. Within the field of neuroimaging, it improves the sensitivity for weak signals and is a promising alternative to existing approaches. The operation of the framework is demonstrated with both synthetic datasets and a real-world neuroimaging example in which the effects of a set of genes on the structure of the brain at the voxel level were measured, and the results compared favorably with those from existing approaches.

Next-generation genotype imputation service and methods

Genotype imputation is a key component of genetic association studies, where it increases power, facilitates meta-analysis, and aids interpretation of signals. Genotype imputation is computationally demanding and, with current tools, typically requires access to a high-performance computing cluster and to a reference panel of sequenced genomes. Here we describe improvements to imputation machinery that reduce computational requirements by more than an order of magnitude with no loss of accuracy in comparison to standard imputation tools. We also describe a new web-based service for imputation that facilitates access to new reference panels and greatly improves user experience and productivity.

Cardiometabolic risk loci share downstream cis- and trans-gene regulation across tissues and diseases

Genome-wide association studies (GWAS) have identified hundreds of cardiometabolic disease (CMD) risk loci. However, they contribute little to genetic variance, and most downstream gene-regulatory mechanisms are unknown. We genotyped and RNA-sequenced vascular and metabolic tissues from 600 coronary artery disease patients in the Stockholm-Tartu Atherosclerosis Reverse Networks Engineering Task study (STARNET). Gene expression traits associated with CMD risk single-nucleotide polymorphism (SNPs) identified by GWAS were more extensively found in STARNET than in tissue- and disease-unspecific gene-tissue expression studies, indicating sharing of downstream cis-/trans-gene regulation across tissues and CMDs. In contrast, the regulatory effects of other GWAS risk SNPs were tissue-specific; abdominal fat emerged as an important gene-regulatory site for blood lipids, such as for the low-density lipoprotein cholesterol and coronary artery disease risk gene PCSK9 STARNET provides insights into gene-regulatory mechanisms for CMD risk loci, facilitating their translation into opportunities for diagnosis, therapy, and prevention.

Associations of genetic risk scores based on adult adiposity pathways with childhood growth and adiposity measures

BACKGROUND

Results from genome-wide association studies (GWAS) identified many loci and biological pathways that influence adult body mass index (BMI). We aimed to identify if biological pathways related to adult BMI also affect infant growth and childhood adiposity measures.

METHODS

We used data from a population-based prospective cohort study among 3,975 children with a mean age of 6 years. Genetic risk scores were constructed based on the 97 SNPs associated with adult BMI previously identified with GWAS and on 28 BMI related biological pathways based on subsets of these 97 SNPs. Outcomes were infant peak weight velocity, BMI at adiposity peak and age at adiposity peak, and childhood BMI, total fat mass percentage, android/gynoid fat ratio, and preperitoneal fat area. Analyses were performed using linear regression models.

RESULTS

A higher overall adult BMI risk score was associated with infant BMI at adiposity peak and childhood BMI, total fat mass, android/gynoid fat ratio, and preperitoneal fat area (all p-values < 0.05). Analyses focused on specific biological pathways showed that the membrane proteins genetic risk score was associated with infant peak weight velocity, and the genetic risk scores related to neuronal developmental processes, hypothalamic processes, cyclicAMP, WNT-signaling, membrane proteins, monogenic obesity and/or energy homeostasis, glucose homeostasis, cell cycle, and muscle biology pathways were associated with childhood adiposity measures (all p-values <0.05). None of the pathways were associated with childhood preperitoneal fat area.

CONCLUSIONS

A genetic risk score based on 97 SNPs related to adult BMI was associated with peak weight velocity during infancy and general and abdominal fat measurements at the age of 6 years. Risk scores based on genetic variants linked to specific biological pathways, including central nervous system and hypothalamic processes, influence body fat development from early life onwards.