Genome-wide Significance Level Research Articles

QTL mapping in outbred tetraploid (and diploid) diallel populations.

Over the last decade, multiparental populations have become a mainstay of genetics research in diploid species. Our goal was to extend this paradigm to autotetraploids by developing software for quantitative trait locus (QTL) mapping in connected F1 populations derived from a set of shared parents. For QTL discovery, phenotypes are regressed on the dosage of parental haplotypes to estimate additive effects. Statistical properties of the model were explored by simulating half-diallel diploid and tetraploid populations with different population sizes and numbers of parents. Across scenarios, the number of progeny per parental haplotype (pph) largely determined the statistical power for QTL detection and accuracy of the estimated haplotype effects. Multiallelic QTL with heritability 0.2 were detected with 90% probability at 25 pph and genome-wide significance level 0.05, and the additive haplotype effects were estimated with over 90% accuracy. Following QTL discovery, the software enables a comparison of models with multiple QTL and nonadditive effects. To illustrate, we analyzed potato tuber shape in a half-diallel population with three tetraploid parents. A well-known QTL on chromosome 10 was detected, for which the inclusion of digenic dominance lowered the Deviance Information Criterion (DIC) by 17 points compared to the additive model. The final model also contained a minor QTL on chromosome 1, but higher-order dominance and epistatic effects were excluded based on the DIC. In terms of practical impacts, the software is already being used to select offspring based on the effect and dosage of particular haplotypes in breeding programs.

Association of Variants Near the Bradykinin Receptor B2 Gene With Angioedema in Patients Taking ACE Inhibitors

BackgroundAngioedema is a rare but potentially life-threatening adverse reaction associated with angiotensin-converting enzyme (ACE) inhibitors. Identification of potential genetic factors related to this adverse event may help identify at-risk patients. ObjectivesThe aim of this study was to identify genetic factors associated with ACE inhibitor–associated angioedema. MethodsA genomewide association study involving patients of European descent, all taking ACE inhibitors, was conducted in a discovery cohort (Copenhagen Hospital Biobank), and associations were confirmed in a replication cohort (Swedegene). Cases were defined as subjects with angioedema events and filled prescriptions for ACE inhibitors ≤180 days before the events. Control subjects were defined as those with continuous treatment with ACE inhibitors without any history of angioedema. Odds ratios (ORs) and 95% confidence intervals (CIs) were computed for angioedema risk using logistic mixed model regression analysis. Summary statistics from the discovery and replication cohorts were analyzed using a fixed-effects meta-analysis model. ResultsThe discovery cohort consisted of 462 cases and 53,391 ACE inhibitor–treated control subjects. The replication cohort consisted of 142 cases and 1,345 ACE inhibitor–treated control subjects. In the discovery cohort, 1 locus, residing at chromosome 14q32.2, was identified that associated with angioedema at the genomewide significance level of P <5 × 10–8. The lead variant at this locus, rs34485356, is an intergenic variant located 60 kb upstream of BDKRB2 (OR: 1.62; 95% CI: 1.38 to 1.90; P = 4.3 × 10−9). This variant was validated in our replication cohort with a similar direction and effect size (OR: 1.60; 95% CI: 1.13 to 2.25; P = 7.2 × 10−3). We found that carriers of the risk allele had significantly lower systolic (−0.46 mm Hg per T allele; 95% CI: −0.83 to −0.10; P = 0.013) and diastolic (−0.26 mm Hg per T allele; 95% CI: −0.46 to −0.05; P = 0.013) blood pressure. ConclusionsIn this genomewide association study involving individuals treated with ACE inhibitors, we found that common variants located in close proximity to the bradykinin receptor B2 gene were associated with increased risk for ACE inhibitor–related angioedema.

Assessing the role of cortisol in cancer: a wide-ranged Mendelian randomisation study.

Cortisol's immunosuppressive, obesogenic, and hyperglycaemic effects suggest that it may play a role in cancer development. However, whether cortisol increases cancer risk is not known. We investigated the potential causal association between plasma cortisol and risk of overall and common site-specific cancers using Mendelian randomisation. Three genetic variants associated with morning plasma cortisol levels at the genome-wide significance level (P < 5 × 10-8) in the Cortisol Network consortium were used as genetic instruments. Summary-level genome-wide association study data for the cancer outcomes were obtained from large-scale cancer consortia, the UK Biobank, and the FinnGen consortium. Two-sample Mendelian randomisation analyses were performed using the fixed-effects inverse-variance weighted method. Estimates across data sources were combined using meta-analysis. A standard deviation increase in genetically predicted plasma cortisol was associated with increased risk of endometrial cancer (odds ratio 1.50, 95% confidence interval 1.13-1.99; P = 0.005). There was no significant association between genetically predicted plasma cortisol and risk of other common site-specific cancers, including breast, ovarian, prostate, colorectal, lung, or malignant skin cancer, or overall cancer. These results indicate that elevated plasma cortisol levels may increase the risk of endometrial cancer but not other cancers. The mechanism by which this occurs remains to be investigated.

A Genome-Wide Association Study of Age-Related Hearing Impairment in Middle- and Old-Aged Chinese Twins.

Background Age-related hearing impairment (ARHI) is considered an unpreventable disorder. We aimed to detect specific genetic variants that are potentially related to ARHI via genome-wide association study (GWAS). Methods A sample of 131 dizygotic twins was genotyped for single-nucleotide polymorphism- (SNP-) based GWAS. Gene-based test was performed using VEGAS2. Pathway enrichment analysis was conducted by PASCAL. Results The twins are with a median age of 49 years, of which 128 were females and 134 were males. rs6633657 was the only SNP that reached the genome-wide significance level for better ear hearing level (BEHL) at 2.0 kHz (P = 1.19 × 10−8). Totally, 9, 10, 42, 7, 17, and 5 SNPs were suggestive evidence level for (P < 1 × 10−5) BEHLs at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz and pure tone average (PTA), respectively. Several promising genetic regions in chromosomes (near the C20orf196, AQPEP, UBQLN3, OR51B5, OR51I2, OR52D1, GLTP, GIT2, and PARK2) nominally associated with ARHI were identified. Gene-based analysis revealed 165, 173, 77, 178, 170, and 145 genes nominally associated with BEHLs at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz and PTA, respectively (P < 0.05). For BEHLs at 0.5, 1.0, and 2.0 kHz, the main enriched pathways were phosphatidylinositol signaling system, regulation of ornithine decarboxylase, eukaryotic translation initiation factor (EIF) pathway, amine compound solute carrier (SLC) transporters, synthesis of phosphoinositides (PIPS) at the plasma membrane, and phosphatidylinositols (PI) metabolism. Conclusions The genetic variations reported herein are significantly involved in functional genes and regulatory domains that mediate ARHI pathogenesis. These findings provide clues for the further unraveling of the molecular physiology of hearing functions and identifying novel diagnostic biomarkers and therapeutic targets of ARHI.

Genome-wide association study in the Taiwan Biobank identifies four novel genes for human height: NABP2, RASA2, RNF41 and SLC39A5.

Numerous genome-wide association studies (GWASs) have been conducted for the identification of genetic variants involved with human height. The vast majority of these studies, however, have been conducted in populations of European ancestry. Here, we report the first GWAS of adult height in the Taiwan Biobank using a discovery sample of 14 571 individuals and an independent replication sample of 20 506 individuals. From our analysis, we generalize to the Taiwanese population genome-wide significant associations with height and 18 previously identified genes in European and non-Taiwanese East Asian populations. We also identify and replicate, at the genome-wide significance level, associated variants for height in four novel genes at two loci that have not previously been reported: RASA2 on chromosome 3 and NABP2, RNF41 and SLC39A5 at 12q13.3 on chromosome 12. RASA2 and RNF41 are strong candidates for having a role in height with copy number and loss of function variants in RASA2 previously found to be associated with short stature disorders, and decreased expression of the RNF41 gene resulting in insulin resistance in skeletal muscle. The results from our analysis of the Taiwan Biobank underscore the potential for the identification of novel genetic discoveries in underrepresented worldwide populations, even for traits, such as height, that have been extensively investigated in large-scale studies of European ancestry populations.

A further survey of the quantitative trait loci affecting swine body size and carcass traits in five related pig populations.

Breeding for good meat quality performance while maintaining large body size and desirable carcass traits has been the major challenge for modern swine selective breeding. To address this goal, in the present work we studied five related populations produced by two commercial breeds (Berkshire and Duroc) and two Chinese breeds (Licha black pig and Lulai black pig). A single-trait GWAS performed on 20 body size and carcass traits using a self-developed China Chip-1 porcine SNP50K BeadChip identified 11 genome-wide significant QTL on nine chromosomes and 22 suggestive QTL on 15 chromosomes. For the 11 genome-wide significant QTL, eight were detected in at least two populations, and the rest were population-specific and only mapped in Shanxia black pig. Most of the genome-wide significant QTL were pleiotropic; for example, the QTL around 75.65Mb on SSC4 was associated with four traits at genome-wide significance level. After screening the genes within 50kb of the top SNP for each genome-wide significant QTL, NR6A1 and VRTN were chosen as candidate genes for vertebrae number; PLAG1 and BMP2 were identified as candidate genes for body size; and MC4R was the strong candidate gene for body weight. The four genes have been reported as candidates for thoracic vertebrae number, lumbar vertebrae number, carcass length and body weight respectively in previous studies. The effects of VRTN on thoracic vertebrae number, carcass length and body length have been verified in Shanxia black pig. Therefore, the VRTN genotype could be used in gene-assisted selection, and this could accelerate genetic improvement of body size and carcass traits in Shanxia black pig.

Clinical and Germline Risk Factors for the Occurrence of Multiple Treatment Toxicities during Childhood Acute Lymphoblastic Leukemia (ALL) Therapy

Treatment-related toxicities can be severe, life-threatening events that impact on delivery of ALL therapy and impair quality of life. The occurrence of multiple (≥ 2) toxicities in any one individual is of particular clinical concern, due to cumulative short and long-term impacts on health and functioning. Few studies have examined the co-occurrence of and risk factors for multiple non-infectious, non hematological toxicities in children with ALL.We evaluated clinical factors and germline single nucleotide polymorphisms (SNPs), to determine the incidence of multiple toxicities and associated risk factors; and to provide insights into their pathogenesis and prevention. We hypothesized that susceptibility factors associated with a two-fold increased risk of multiple toxicity could be identified, and that these may serve as predictive biomarkers.A systematic retrospective review of 1251 Australian children (1-18 years) treated on consecutive ALL protocols, between 1998-2013, was performed. Selected toxicities were venous thromboembolism (VTE), insulin-requirement, bone toxicity (osteonecrosis and fractures); neurotoxicity (central and peripheral). Toxicities chosen were those with a reported incidence >5%, robustness of retrospective data capture and clinical significance. Multiple toxicities were defined as ≥2 selected toxicities occurring in the one patient during first-line ALL therapy. Multivariate logistic regression, Cox regression and Kaplan-Meier survival analyses were performed. Two-tailed P values <0.05 were considered significant. A genome-wide association study (GWAS) was conducted on 702 Caucasian individuals within the cohort, which included 64 cases (≥2 toxicities) and 638 controls (≤ 1 toxicity). Samples were genotyped on the Illumina Oncoarray platform and imputed to the 1000G reference set. Genome-wide significance levels for the GWAS were P < 5E-08.There were 1033 children treated on BFM-based protocols and 218 children on COG-based protocols (Table 1, “Baseline characteristics”). The median age of the cohort was 59 months (range 9-218 months) and median follow-up time was 78 months (range 3-186 months).Five-year overall survival (OS), event-free survival (EFS) and leukemia-free survival (LFS) were 92% ±0.8%, 83.8% ±1.1% and 85.6% ±1.1% respectively. Multiple (≥ 2) toxicities occurred in 133/1240 evaluable children (10.73%). There was no difference in survival for children who experienced multiple toxicities (n=133) compared to the rest of the cohort (n=1107) (P =0.138 for OS, P= 0.595 for EFS, P = 0.276 for LFS).Independent clinical risk factors in multivariate analysis were age ≥ 10 years, female sex, and high-risk ALL. Additional laboratory test risk factors in induction/consolidation were an elevated gamma glutamyl-transferase (GGT) > 5 x upper limit normal (ULN) and serum albumin <20g/L (Table 2). The analysis was adjusted for treatment platform (BFM vs COG), gender, age and immunophenotype (T vs B).An intronic SNP rs78342591 in the MUC16 locus reached genome-wide significance for multiple toxicity (P =2.42E-08, odds ratio 5.89, 95% confidence interval (CI) 3.23-10.74). This was an imputed SNP.Clinical and germline factors could identify children who may benefit from targeted supportive care and surveillance for acute toxicity following ALL therapy. Future research directions include replication, direct sequencing and subsequent functional validation, focusing on MUC16 expression and modulation of drug toxicity. Translation of this research may involve analysis of both risk and protective mechanisms associated with development of multiple toxicities in childhood ALL. [Display omitted] DisclosuresNo relevant conflicts of interest to declare.

Identification of a Locus Near ULK1 Associated With Progression-Free Survival in Ovarian Cancer.

Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10-8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10-8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. This finding provides insight into genetic markers associated with EOC outcome and potential treatment options.See related commentary by Peres and Monteiro, p. 1604.

Identifying novel genetic variants for brain amyloid deposition: a genome-wide association study in the Korean population

BackgroundGenome-wide association studies (GWAS) have identified a number of genetic variants for Alzheimer’s disease (AD). However, most GWAS were conducted in individuals of European ancestry, and non-European populations are still underrepresented in genetic discovery efforts. Here, we performed GWAS to identify single nucleotide polymorphisms (SNPs) associated with amyloid β (Aβ) positivity using a large sample of Korean population.MethodsOne thousand four hundred seventy-four participants of Korean ancestry were recruited from multicenters in South Korea. Discovery dataset consisted of 1190 participants (383 with cognitively unimpaired [CU], 330 with amnestic mild cognitive impairment [aMCI], and 477 with AD dementia [ADD]) and replication dataset consisted of 284 participants (46 with CU, 167 with aMCI, and 71 with ADD). GWAS was conducted to identify SNPs associated with Aβ positivity (measured by amyloid positron emission tomography). Aβ prediction models were developed using the identified SNPs. Furthermore, bioinformatics analysis was conducted for the identified SNPs.ResultsIn addition to APOE, we identified nine SNPs on chromosome 7, which were associated with a decreased risk of Aβ positivity at a genome-wide suggestive level. Of these nine SNPs, four novel SNPs (rs73375428, rs2903923, rs3828947, and rs11983537) were associated with a decreased risk of Aβ positivity (p < 0.05) in the replication dataset. In a meta-analysis, two SNPs (rs7337542 and rs2903923) reached a genome-wide significant level (p < 5.0 × 10−8). Prediction performance for Aβ positivity increased when rs73375428 were incorporated (area under curve = 0.75; 95% CI = 0.74–0.76) in addition to clinical factors and APOE genotype. Cis-eQTL analysis demonstrated that the rs73375428 was associated with decreased expression levels of FGL2 in the brain.ConclusionThe novel genetic variants associated with FGL2 decreased risk of Aβ positivity in the Korean population. This finding may provide a candidate therapeutic target for AD, highlighting the importance of genetic studies in diverse populations.

Genetic risk of clozapine-induced leukopenia and neutropenia: a genome-wide association study

BackgroundClozapine is considered to be the most effective antipsychotic medication for schizophrenia. However, it is associated with several adverse effects such as leukopenia, and the underlying mechanism has not yet been fully elucidated. The authors performed a genome-wide association study (GWAS) in a Chinese population to identify genetic markers for clozapine-induced leukopenia (CIL) and clozapine-induced neutropenia (CIN).MethodsA total of 1879 patients (225 CIL cases, including 43 CIN cases, and 1,654 controls) of Chinese descent were included. Data from common and rare single nucleotide polymorphisms (SNPs) were tested for association. The authors also performed a trans-ancestry meta-analysis with GWAS results of European individuals from the Clozapine-Induced Agranulocytosis Consortium (CIAC).ResultsThe authors identified several novel loci reaching the threshold of genome-wide significance level (P < 5 × 10−8). Three novel loci were associated with CIL while six were associated with CIN, and two T cell related genes (TRAC and TRAT1) were implicated. The authors also observed that one locus with evidence close to genome-wide significance (P = 5.08 × 10−8) was near the HLA-B gene in the major histocompatibility complex region in the trans-ancestry meta-analysis.ConclusionsThe associations provide novel and valuable understanding of the genetic and immune causes of CIL and CIN, which is useful for improving clinical management of clozapine related treatment for schizophrenia. Causal variants and related underlying molecular mechanisms need to be understood in future developments.

Genetically predicted circulating vitamin C in relation to cardiovascular disease.

We conducted a two-sample Mendelian randomization (MR) study to assess the associations of genetically predicted circulating vitamin C levels with cardiovascular diseases (CVDs). Ten lead single-nucleotide polymorphisms associated with plasma vitamin C levels at the genome-wide significance level were used as instrumental variables. Summary-level data for 15 CVDs were obtained from corresponding genetic consortia, the UK Biobank study, and the FinnGen consortium. The inverse-variance-weighted method was the primary analysis method, supplemented by the weighted median and MR-Egger methods. Estimates for each CVD from different sources were combined. Genetically predicted vitamin C levels were not associated with any CVD after accounting for multiple testing. However, there were suggestive associations of higher genetically predicted vitamin C levels (per 1 standard deviation increase) with lower risk of cardioembolic stroke [odds ratio, 0.79; 95% confidence interval (CI), 0.64, 0.99; P = 0.038] and higher risk of atrial fibrillation (odds ratio, 1.09; 95% CI, 1.00, 1.18; P = 0.049) in the inverse-variance-weighted method and with lower risk of peripheral artery disease (odds ratio, 0.76, 95% CI, 0.62, 0.93; P = 0.009) in the weighted median method. We found limited evidence with MR techniques for an overall protective role of vitamin C in the primary prevention of CVD. The associations of vitamin C levels with cardioembolic stroke, atrial fibrillation, and peripheral artery disease need further study.

Joint Genome-Wide Association Analyses Identified 49 Novel Loci For Age at Natural Menopause.

Age at natural menopause (ANM) is an important index for women's health. Either early or late ANM is associated with a series of adverse outcomes later in life. Despite being an inheritable trait, its genetic determinant has not yet been fully understood. Aiming to better characterize the genetic architecture of ANM, we conducted genome-wide association study (GWAS) meta-analyses in European-specific as well as trans-ancestry samples by using GWAS summary statistics from the following 3 large studies: the Reproductive Genetics Consortium (ReproGen; N = 69 626), the UK Biobank cohort (UKBB; N = 111 593) and the BioBank Japan Project (BBJ; N = 43 861), followed by a series of bioinformatical assessments and functional annotations. By integrating the summary statistics from the 3 GWAS of up to 225 200 participants, this largest meta-analysis identified 49 novel loci and 3 secondary signals that were associated with ANM at the genome-wide significance level (P < 5 × 10-8). No population specificity or heterogeneity was observed at most of the associated loci. Functional annotations prioritized 90 candidate genes at the newly identified loci. Among the 26 traits that were genetically correlated with ANM, hormone replacement therapy (HRT) exerted a causal relationship, implying a causal pattern by which HRT was determined by ANM. Our findings improved our understanding of the etiology of female menopause, as well as shed light on potential new therapies for abnormal menopause.

Are there monogenic hereditary forms of bladder cancer or only genetic susceptibilities?

Bladder cancer (BC) is the most common cancer involving the urinary system and the ninth most common cancer worldwide. Tobacco smoking is the most important environmental risk factor of BC. Several single nucleotide polymorphisms have been validated by genome-wide association studies as genetic risk factors for BC. However, the identification of DNA mismatch-repair genes, including MSH2 in Lynch syndrome and MUTYH in MUTYH-associated polyposis, raises the possibility of monogenic hereditary forms of BC. Moreover, other genetic mutations may play a key role in familial and hereditary transmissions of BC. Therefore, the aim of this review is to focus on the major hereditary syndromes involved in the development of BC and to report BC genetic susceptibilities established with genome-wide significance level.

Causal effect of renal function on venous thromboembolism: a two-sample Mendelian randomization investigation

Whether renal function is causally associated with venous thromboembolism (VTE) is not yet fully elucidated. We conducted a two-sample Mendelian randomization (MR) study to determine the causal effect of renal function, measured as estimated glomerular filtration rate (eGFR), on VTE. Single-nucleotide polymorphisms associated with eGFR were selected as instrumental variables at the genome-wide significance level (p < 5 × 10−8) from a meta-analysis of 122 genome-wide association studies including up to 1,046,070 individuals. Summary-level data for VTE were obtained from the FinnGen consortium (6913 VTE cases and 169,986 non-cases) and UK Biobank study (4620 VTE cases and 356,574 non-cases). MR estimates were calculated using the random-effects inverse-variance weighted method and combined using fixed-effects meta-analysis. Genetically predicted decreased eGFR was significantly associated with an increased risk of VTE in both FinnGen and UK Biobank. For one-unit decrease in log-transformed eGFR, the odds ratios of VTE were 2.93 (95% confidence interval (CI) 1.25, 6.84) and 4.46 (95% CI 1.59, 12.5) when using data from FinnGen and UK Biobank, respectively. The combined odds ratio was 3.47 (95% CI 1.80, 6.68). Results were consistent in all sensitivity analyses and no horizontal pleiotropy was detected. This MR-study supported a casual role of impaired renal function in VTE.

Genome wide association study of response to interval and continuous exercise training: the Predict-HIIT study

BackgroundLow cardiorespiratory fitness (V̇O2peak) is highly associated with chronic disease and mortality from all causes. Whilst exercise training is recommended in health guidelines to improve V̇O2peak, there is considerable inter-individual variability in the V̇O2peak response to the same dose of exercise. Understanding how genetic factors contribute to V̇O2peak training response may improve personalisation of exercise programs. The aim of this study was to identify genetic variants that are associated with the magnitude of V̇O2peak response following exercise training.MethodsParticipant change in objectively measured V̇O2peak from 18 different interventions was obtained from a multi-centre study (Predict-HIIT). A genome-wide association study was completed (n = 507), and a polygenic predictor score (PPS) was developed using alleles from single nucleotide polymorphisms (SNPs) significantly associated (P < 1 × 10–5) with the magnitude of V̇O2peak response. Findings were tested in an independent validation study (n = 39) and compared to previous research.ResultsNo variants at the genome-wide significance level were found after adjusting for key covariates (baseline V̇O2peak, individual study, principal components which were significantly associated with the trait). A Quantile–Quantile plot indicates there was minor inflation in the study. Twelve novel loci showed a trend of association with V̇O2peak response that reached suggestive significance (P < 1 × 10–5). The strongest association was found near the membrane associated guanylate kinase, WW and PDZ domain containing 2 (MAGI2) gene (rs6959961, P = 2.61 × 10–7). A PPS created from the 12 lead SNPs was unable to predict V̇O2peak response in a tenfold cross validation, or in an independent (n = 39) validation study (P > 0.1). Significant correlations were found for beta coefficients of variants in the Predict-HIIT (P < 1 × 10–4) and the validation study (P < × 10–6), indicating that general effects of the loci exist, and that with a higher statistical power, more significant genetic associations may become apparent.ConclusionsOngoing research and validation of current and previous findings is needed to determine if genetics does play a large role in V̇O2peak response variance, and whether genomic predictors for V̇O2peak response trainability can inform evidence-based clinical practice.Trial registration Australian New Zealand Clinical Trials Registry (ANZCTR), Trial Id: ACTRN12618000501246, Date Registered: 06/04/2018, http://www.anzctr.org.au/Trial/Registration/TrialReview.aspx?id=374601&isReview=true.

Genome-wide association study reveals novel loci associated with feeding behavior in Pekin ducks

BackgroundFeeding behavior traits are an essential part of livestock production. However, the genetic base of feeding behavior traits remains unclear in Pekin ducks. This study aimed to determine novel loci related to feeding behavior in Pekin ducks.ResultsIn this study, the feeding information of 540 Pekin ducks was recorded, and individual genotype was evaluated using genotyping-by-sequencing methods. Genome-wide association analysis (GWAS) was conducted for feeding behavior traits. Overall, thirty significant (P-value < 4.74E-06) SNPs for feeding behavior traits were discovered, and four of them reached the genome-wide significance level (P-value < 2.37E-07). One genome-wide significance locus associated with daily meal times was located in a 122.25 Mb region on chromosome 2, which was within the intron of gene ubiquitin-conjugating enzyme E2 E2 (UBE2E2), and could explain 2.64% of the phenotypic variation. This locus was also significantly associated with meal feed intake, and explained 2.72% of this phenotypic variation.ConclusionsThis study is the first GWAS for feeding behavior traits in ducks. Our results provide a list of candidate genes associated with feeding behavior, and also help to better understand the genetic mechanisms of feeding behavior patterns in ducks.

DNA methylome analysis identifies BMI-related epigenetic changes associated with non-small cell lung cancer susceptibility.

BackgroundBody mass index (BMI) has been reported to be inversely associated with incident risk of non‐small cell lung cancer (NSCLC). However, the underlying mechanism is still unclear. This study aimed to investigate the role of DNA methylation in the relationship between BMI and NSCLC.MethodsWe carried out a genome‐wide DNA methylation study of BMI in peripheral blood among 2266 Chinese participants by using Illumina Methylation arrays. For the BMI‐related DNA methylation changes, their associations with NSCLC risk were further analyzed and their mediation effects on BMI‐NSCLC association were also evaluated.ResultsThe methylation levels of four CpGs (cg12593793, cg17061862, cg11024682, and cg06500161, annotated to LMNA, ZNF143, SREBF1, and ABCG1, respectively) were found to be significantly associated with BMI. Methylation levels of cg12593793, cg11024682, and cg06500161 were observed to be inversely associated with NSCLC risk [OR (95%CI) =0.22 (0.16, 0.31), 0.39 (0.30, 0.50), and 0.66 (0.53, 0.82), respectively]. Additionally, cg11024682 in SREBF1 and cg06500161 in ABCG1 mediated 45.3% and 19.5% of the association between BMI and decreased NSCLC risk, respectively.ConclusionsIn this study, we identified four DNA methylation sites associated with BMI in the Chinese populations at the genome‐wide significant level. We also found that the BMI‐related methylations of SREBF1 and ABCG1 could mediate about a quintile‐to‐half of the effect of BMI on reduced NSCLC risk, which adds a potential mechanism underlying this association.

Homocysteine, B vitamins, and cardiovascular disease: a Mendelian randomization study

BackgroundWhether a modestly elevated homocysteine level is causally associated with an increased risk of cardiovascular disease remains unestablished. We conducted a Mendelian randomization study to assess the associations of circulating total homocysteine (tHcy) and B vitamin levels with cardiovascular diseases in the general population.MethodsIndependent single nucleotide polymorphisms associated with tHcy (n = 14), folate (n = 2), vitamin B6 (n = 1), and vitamin B12 (n = 14) at the genome-wide significance level were selected as instrumental variables. Summary-level data for 12 cardiovascular endpoints were obtained from genetic consortia, the UK Biobank study, and the FinnGen consortium.ResultsHigher genetically predicted circulating tHcy levels were associated with an increased risk of stroke. For each one standard deviation (SD) increase in genetically predicted tHcy levels, the odds ratio (OR) was 1.11 (95% confidence interval (CI), 1.03, 1.21; p = 0.008) for any stroke, 1.26 (95% CI, 1.05, 1.51; p = 0.013) for subarachnoid hemorrhage, and 1.11 (95% CI, 1.03, 1.21; p = 0.011) for ischemic stroke. Higher genetically predicted folate levels were associated with decreased risk of coronary artery disease (ORSD, 0.88; 95% CI, 0.78, 1.00, p = 0.049) and any stroke (ORSD, 0.86; 95% CI, 0.76, 0.97, p = 0.012). Genetically predicted increased vitamin B6 levels were associated with a reduced risk of ischemic stroke (ORSD, 0.88; 95% CI, 0.81, 0.97, p = 0.009). None of these associations persisted after multiple testing correction. There was no association between genetically predicted vitamin B12 and cardiovascular disease.ConclusionsThis study reveals suggestive evidence that B vitamin therapy and lowering of tHcy may reduce the risk of stroke, particularly subarachnoid hemorrhage and ischemic stroke.

Functional Genomic Analyses of the 21q22.3 Locus Identifying Functional Variants and Candidate Gene YBEY for Breast Cancer Risk.

Simple SummaryPrevious research has revealed a genetic predisposition to breast carcinogenesis. Thus, identifying causal genetic variants and their associated gene networks may improve breast cancer diagnostics and risk assessment. Our study investigated YBEY, an uncharacterized gene in humans, and its functions in breast cancer risk and progression. We identified two genetic variants associated with YBEY expression that may have causal functions in breast cancer risk. We performed in vitro functional assays using MCF-7, T47D, and MDA-MB-231 breast cancer cell lines and showed that knockdown of YBEY expression significantly inhibited proliferation, colony formation, and invasion/migration. We utilized RNA sequencing to identify gene networks associated with YBEY knockdown including inflammation and metabolic pathways. Further, we used data available in The Cancer Genome Atlas to explore trends in YBEY expression patterns in normal and tumor tissues. Our study provides a role for YBEY in breast carcinogenesis, and further studies investigating its mechanistic functions are warranted.We previously identified a locus at 21q22.3, tagged by the single nucleotide polymorphism (SNP) rs35418111, being associated with breast cancer risk at a genome-wide significance level; however, the underlying causal functional variants and gene(s) responsible for this association are unknown. We performed functional genomic analyses to identify potential functional variants and target genes that may mediate this association. Functional annotation for SNPs in high linkage disequilibrium (LD, r2 > 0.8) with rs35418111 in Asians showed evidence of promoter and/or enhancer activities, including rs35418111, rs2078203, rs8134832, rs57385578, and rs8126917. These five variants were assessed for interactions with nuclear proteins by electrophoretic mobility shift assays. Our results showed that the risk alleles for rs2078203 and rs35418111 altered DNA-protein interaction patterns. Cis-expression quantitative loci (cis-eQTL) analysis, using data from the Genotype-Tissue Expression database (GTEx) European-ancestry female normal breast tissue, indicated that the risk allele of rs35418111 was associated with a decreased expression of the YBEY gene, a relatively uncharacterized endoribonuclease in humans. We investigated the biological effects of YBEY on breast cancer cell lines by transient knock-down of YBEY expression in MCF-7, T47D, and MDA-MB-231 cell lines. Knockdown of YBEY mRNA in breast cancer cell lines consistently decreased cell proliferation, colony formation, and migration/invasion, regardless of estrogen receptor status. We performed RNA sequencing in MDA-MB-231 cells transfected with siRNA targeting YBEY and subsequent gene set enrichment analysis to identify gene networks associated with YBEY knockdown. These data indicated YBEY was involved in networks associated with inflammation and metabolism. Finally, we showed trends in YBEY expression patterns in breast tissues from The Cancer Genome Atlas (TCGA); early-stage breast cancers had elevated YBEY expression compared with normal tissue, but significantly decreased expression in late-stage disease. Our study provides evidence of a significant role for the human YBEY gene in breast cancer pathogenesis and the association between the rs35418111/21q22.3 locus and breast cancer risk, which may be mediated through functional SNPs, rs35418111 and rs2078203, that regulate expression of YBEY.

Genome-Wide Association Study Demonstrates the Role Played by the CD226 Gene in Rasa Aragonesa Sheep Reproductive Seasonality.

Simple SummaryTo elucidate the genetic basis of reproductive seasonality in Rasa Aragonesa sheep breed, we performed a genome-wide association study (GWAS) in order to detect single nucleotide polymorphisms (SNPs) or regions associated with traits related to ovarian function and behavioural signs of estrous. The GWAS included 205 ewes with genotypes for 583882 SNPs. Only one SNP overcame the genome-wide significance level. Nine potential SNPs overcame the chromosome-wise significance level (FDR 10%). Gene annotation demonstrated that CD226 molecule (CD226) and neuropeptide Y (NPY) genes that could be involved in reproductive seasonality were close to the significant SNPs. To validate the results, we sequenced the entire coding region of the NPY gene and four exons of the CD226 gene to search for polymorphisms that could be involved in the phenotypes studied. Two synonymous and two nonsynonymous SNPs in the NPY and CD226 genes, respectively, were genotyped in the whole population. We demonstrated that the AA genotype of the SNP rs404360094 located in exon 3 of the CD226 gene was associated with higher and lower total days of anoestrus and oestrous cycling months, respectively. Therefore, this SNP could be utilized as a genetic marker for assisted selection marker to reduce seasonality.A genome-wide association study (GWAS) was used to identify genomic regions influencing seasonality reproduction traits in Rasa Aragonesa sheep. Three traits associated with either ovarian function based on blood progesterone levels (total days of anoestrus and progesterone cycling months) or behavioral signs of oestrous (oestrous cycling months) were studied. The GWAS included 205 ewes genotyped using the 50k and 680k Illumina Ovine Beadchips. Only one SNP associated with the progesterone cycling months overcame the genome-wide significance level (rs404991855). Nine SNPs exhibited significant associations at the chromosome level, being the SNPs rs404991855 and rs418191944, that are located in the CD226 molecule (CD226) gene, associated with the three traits. This gene is related to reproductive diseases. Two other SNPs were located close to the neuropeptide Y (NPY) gene, which is involved in circadian rhythms. To validate the GWAS, partial characterization of both genes by Sanger sequencing, and genotyping of two synonymous and two nonsynonymous SNPs in the NPY and CD226 genes, respectively, were performed. SNP association analysis showed that only SNP rs404360094 in the exon 3 of the CD226 gene, which produces an amino acid substitution from asparagine (uncharged polar) to aspartic acid (acidic), was associated with the three seasonality traits. Our results suggest that the CD226 gene may be involved in the reproductive seasonality in Rasa Aragonesa.