A GWAS assessment of the contribution of genomic imprinting to the variation of body mass index in mice

Background

Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis.

Results

10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10⁻⁴ to 1.842 ×10⁻⁴ when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered.

Conclusions

When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).     

Genetics of complex disorders

The success stories of identifying genes in Mendelian disorders have stimulated research that aims at identifying the genetic determinants in complex disorders, in which both genetics, environment and chance affect the pathogenetic processes. This review summarizes the brief history and lessons learned from genetic analysis of complex disorders and outlines some landscapes ahead for medical research.     

Penalized Multimarker vs. Single-Marker Regression Methods for Genome-Wide Association Studies of Quantitative Traits

The data from genome-wide association studies (GWAS) in humans are still predominantly analyzed using single-marker association methods. As an alternative to single-marker analysis (SMA), all or subsets of markers can be tested simultaneously. This approach requires a form of penalized regression (PR) as the number of SNPs is much larger than the sample size. Here we review PR methods in the context of GWAS, extend them to perform penalty parameter and SNP selection by false discovery rate (FDR) control, and assess their performance in comparison with SMA. PR methods were compared with SMA, using realistically simulated GWAS data with a continuous phenotype and real data. Based on these comparisons our analytic FDR criterion may currently be the best approach to SNP selection using PR for GWAS. We found that PR with FDR control provides substantially more power than SMA with genome-wide type-I error control but somewhat less power than SMA with Benjamini–Hochberg FDR control (SMA-BH). PR with FDR-based penalty parameter selection controlled the FDR somewhat conservatively while SMA-BH may not achieve FDR control in all situations. Differences among PR methods seem quite small when the focus is on SNP selection with FDR control. Incorporating linkage disequilibrium into the penalization by adapting penalties developed for covariates measured on graphs can improve power but also generate more false positives or wider regions for follow-up. We recommend the elastic net with a mixing weight for the Lasso penalty near 0.5 as the best method.     

ezQTL: A Web Platform for Interactive Visualization and Colocalization of QTLs and GWAS Loci

Genome-wide association studies(GWAS) have identified thousands of genomic loci associated with complex diseases and traits, including cancer. The vast majority of common traitassociated variants identified via GWAS fall in non-coding regions of the genome, posing a challenge in elucidating the causal variants, genes, and mechanisms involved. Expression quantitative trait locus(e QTL) and other molecular QTL studies have been valuable resources in identifying candidate causal genes from GWAS loc...     

Pathway analysis of a genome-wide association study in schizophrenia

Objective

The aim of this study was to identify the candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms that contribute to schizophrenia susceptibility and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach.

Methods

We used schizophrenia GWAS data of the genotypes of 660,259 SNPs in 1378 controls and 1351 cases of European descent after quality control filtering. ICSNPathway (Identify candidate Causal SNPs and Pathways) analysis was applied to the schizophrenia GWAS dataset. The first stage involved the pre-selection of candidate SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate SNPs using improved-gene set enrichment analysis.

Results

ICSNPathway analysis identified fifteen candidate SNPs, ten candidate pathways, and nine hypothetical biological mechanisms. The most strongly associated potential pathways were as follows. First, rs1644731 and rs1644730 to RDH8 to estrogen biosynthetic process (p < 0.001, FDR < 0.001). The genes involved in this pathway are RDH8 and HSD3B1 (p < 0.05). All-trans-retinol dehydrogenase (RDH8) is a visual cycle enzyme that reduces all-trans-retinal to all-trans-retinol in the presence of NADPH. The chemical reactions and pathways involved result in the formation of estrogens, which are C18 steroid hormones that can stimulate the development of female sexual characteristics. Second, rs1146031 to ACVR1 to mesoderm formation and activin binding (p < 0.001, FDR = 0.032, 0.034). Two of 15 candidate genes are known genes associated with schizophrenia: KCNQ2 and APOL2. One of the 10 candidate pathways, estrogen biosynthetic process, is known to be associated with schizophrenia (p < 0.001, FDR < 0.001). However, 13 of candidate genes (RDH8, ACVR1, PSMD9, KCNAB1, SLC17A3, ARCN1, COG7, STAB2, LRPAP1, STAB1, CXCL16, COL4A4, EXOSC3) and 9 of candidate pathways were novel.

Conclusion

By applying ICSNPathway analysis to schizophrenia GWAS data, we identified candidate SNPs, genes like KCNQ2 and APOL2 and pathways involving the estrogen biosynthetic process may contribute to schizophrenia susceptibility. Further analyses are needed to validate the results of this analysis.     

中国精神分裂症的全基因组关联分析及其转化医学进展

我国在精神分裂症的遗传学和生命组学研究方面取得了很大进展,如在全基因组关联分析(genome-wide association study,GWAS)方面工作获得了一系列成果.随着我国对重大疾病转化医学的逐步关注和重视,利用在精神分裂症上已经获得的广泛和深入的研究结果,寻找精神分裂症各种临床应用的生物标记物研究,系统性地建立适合于类似精神分裂症这类复杂疾病的早期诊断、干预和预防的临床咨询和应用体系等将是该疾病转化医学方面可实施的方法和案例.精神分裂症的转化医学方面还涉及精神分裂症患者的个体化用药方案建立.药物疗效和药物不良反应的个体差异具有较复杂的环境和遗传背景,结合精神分裂症的遗传学病因和药物作用的遗传学差异,将有效发挥治疗药物的功效,并降低重大不良反应在敏感个体上的发生.对精神分裂症这类给国家和社会带来极其重大负担的重大疾病,积极推动我国在此类疾病上的基础研究成果转化和转化医学的实施具有重要的社会效应和积极的带动作用.     

Gene set analysis of genome-wide association studies: methodological issues and perspectives

Recent studies have demonstrated that gene set analysis, which tests disease association with genetic variants in a group of functionally related genes, is a promising approach for analyzing and interpreting genome-wide association studies (GWAS) data. These approaches aim to increase power by combining association signals from multiple genes in the same gene set. In addition, gene set analysis can also shed more light on the biological processes underlying complex diseases. However, current approaches for gene set analysis are still in an early stage of development in that analysis results are often prone to sources of bias, including gene set size and gene length, linkage disequilibrium patterns and the presence of overlapping genes. In this paper, we provide an in-depth review of the gene set analysis procedures, along with parameter choices and the particular methodology challenges at each stage. In addition to providing a survey of recently developed tools, we also classify the analysis methods into larger categories and discuss their strengths and limitations. In the last section, we outline several important areas for improving the analytical strategies in gene set analysis.     

GWAS在2型糖尿病相关基因研究中的应用

2型糖尿病(type 2 diabetes,T2D)是一种常见的复杂疾病,其发病受到遗传和环境因素的共同作用.全基因组关联研究(genome-wide association study,GWAS)是一种可在全基因组范围筛查疾病相关的序列变异的新型群体关联研究方法.近年来,采用GWAS以及在此基础上展开的meta分析,已分别在TCF7L2、HHEX-IDE、SLC30A8、CDKAL1、CDKN2A-CDKN2B、IGF2BP2、NOTCH2、CDC123-CAMK1D、ADAMTS9、THADA、TSPAN8-LGR5、JAZF1等12个基因区域鉴定出多个T2D相关的多态位点.已有的研究提示,上述多个基因可能在胰岛β细胞发育和功能维持方面扮演着重要角色.本文集中介绍了GWAS的原理及其在T2D研究中的优势;回顾了GWAS在T2D研究中的主要发现;并对运用GWAS在T2D研究中尚需解决的问题进行了总结和展望.     

A new method for post Genome-Wide Association Study (GWAS) analysis of colorectal cancer in Taiwan

Recently, single nucleotide polymorphisms (SNPs) located in specific loci or genes have been identified associated with susceptibility to colorectal cancer (CRC) in Genome-Wide Association Studies (GWAS). However, in different ethnicities and regions, the genetic variations and the environmental factors can widely vary. Therefore, here we propose a post-GWAS analysis method to investigate the CRC susceptibility SNPs in Taiwan by conducting a replication analysis and bioinformatics analysis. One hundred and forty-four significant SNPs from published GWAS results were collected by a literature survey, and two hundred and eighteen CRC samples and 385 normal samples were collected for post-GWAS analysis. Finally, twenty-six significant SNPs were identified and reported as associated with susceptibility to colorectal cancer, other cancers, obesity, and celiac disease in a previous GWAS study. Functional analysis results of 26 SNPs indicate that most biological processes identified are involved in regulating immune responses and apoptosis. In addition, an efficient prediction model was constructed by applying Jackknife feature selection and ANOVA testing. As compared to another risk prediction model of CRC for European Caucasians population, which performs 0.616 of AUC by using 54 SNPs, the proposed model shows good performance in predicting CRC risk within the Taiwanese population, i.e., 0.724 AUC by using 16 SNPs. We believe that the proposed risk prediction model is highly promising for predicting CRC risk within the Taiwanese population. In addition, the functional analysis results could be helpful to explore the potential associated regulatory mechanisms that may be involved in CRC development.     

A genome-wide association study of facial morphology identifies novel genetic loci in Han Chinese

The human face is a heritable surface with many complex sensory organs. In recent years, many genetic loci associated with facial features have been reported in different populations, yet there is a lack of studies on the Han Chinese population. Here, we report a genome-wide association study of 3 D normal human faces of 2,659 Han Chinese with autosegment phenotypes of facial morphology. We identify singlenucleotide polymorphisms(SNPs) encompassing four genomic regions showing significant associations with different facial regions, including SNPs in DENND1 B associated with the chin, SNPs among PISRT1 associated with eyes, SNPs between DCHS2 and SFRP2 associated with the nose, and SNPs in VPS13 B associated with the nose. We replicate 24 SNPs from previously reported genetic loci in different populations, whose candidate genes are DCHS2, SUPT3 H, HOXD1, SOX9, PAX3, and EDAR. These results provide a more comprehensive understanding of the genetic basis of variation in human facial morphology.     

Probability of detecting disease-associated single nucleotide polymorphisms in case-control genome-wide association studies

Some case-control genome-wide association studies (CCGWASs) select promising single nucleotide polymorphisms (SNPs) by ranking corresponding p-values, rather than by applying the same p-value threshold to each SNP. For such a study, we define the detection probability (DP) for a specific disease-associated SNP as the probability that the SNP will be "T-selected," namely have one of the top T largest chi-square values (or smallest p-values) for trend tests of association. The corresponding proportion positive (PP) is the fraction of selected SNPs that are true disease-associated SNPs. We study DP and PP analytically and via simulations, both for fixed and for random effects models of genetic risk, that allow for heterogeneity in genetic risk. DP increases with genetic effect size and case-control sample size and decreases with the number of nondisease-associated SNPs, mainly through the ratio of T to N, the total number of SNPs. We show that DP increases very slowly with T, and the increment in DP per unit increase in T declines rapidly with T. DP is also diminished if the number of true disease SNPs exceeds T. For a genetic odds ratio per minor disease allele of 1.2 or less, even a CCGWAS with 1000 cases and 1000 controls requires T to be impractically large to achieve an acceptable DP, leading to PP values so low as to make the study futile and misleading. We further calculate the sample size of the initial CCGWAS that is required to minimize the total cost of a research program that also includes follow-up studies to examine the T-selected SNPs. A large initial CCGWAS is desirable if genetic effects are small or if the cost of a follow-up study is large.     

The Statistical Power of Inclusive Composite Interval Mapping in Detecting Digenic Epistasis Showing Common F2 Segregation Ratios

Epistasis is a commonly observed genetic phenomenon and an important source of variation of complex traits,which could maintain additive variance and therefore assure the long-term genetic gain in breeding.Inclusive composite interval mapping（ICIM） is able to identify epistatic quantitative trait loci（QTLs） no matter whether the two interacting QTLs have any additive effects.In this article,we conducted a simulation study to evaluate detection power and false discovery rate（FDR） of ICIM epistatic mapping,by considering F₂ and doubled haploid（DH） populations,different F₂ segregation ratios and population sizes.Results indicated that estimations of QTL locations and effects were unbiased,and the detection power of epistatic mapping was largely affected by population size,heritability of epistasis,and the amount and distribution of genetic effects.When the same likelihood of odd（LOD） threshold was used,detection power of QTL was higher in F₂ population than power in DH population;meanwhile FDR in F₂ was also higher than that in DH.The increase of marker density from 10 cM to 5 cM led to similar detection power but higher FDR.In simulated populations,ICIM achieved better mapping results than multiple interval mapping（MIM） in estimation of QTL positions and effect.At the end,we gave epistatic mapping results of ICIM in one actual population in rice（Oryza sativa L.）.     

Genomic atlases of introgression and differentiation reveal breeding footprints in Chinese cultivated rice

The long history of cultivation and breeding has left a variety of footprints in the genomes of Asian cultivated rice (Oryza sativa L.). In this study, we focus on two types of genomic footprints, introgression and differentiation, in a population of more than 1200 Chinese rice accessions. We found that a Xian/indica and a temperate Geng/japonica accession respectively contained an average of 19.3-Mb and 6.8-Mb alien introgressed chromosomal segments, of which many contained functional sequence variants, quantitative trait loci, or genes controlling flowering, grain, and resistance traits. Notably, we found most introgressions, including the known heterotic loci Hd3a and TAC1, were distributed differentially between the female and male parents of three-line indica hybrid rice, indicating their potential contribution to heterosis. We also found many differentiated regions between subgroups within a subpopulation contained agronomically important loci, such as DTH7, Hd1 for heading date, and qCT7 for cold tolerance, providing new candidates for studying local adaptation or heterosis. Tracing these footprints allows us to better understand the genetic exchange or differentiation underlying agronomic traits in modern Chinese rice cultivars. These findings also provide potential targets for rice genetic research and breeding.     

Modelling the genetic architecture of flowering time control in barley through nested association mapping

Background

Barley, globally the fourth most important cereal, provides food and beverages for humans and feed for animal husbandry. Maximizing grain yield under varying climate conditions largely depends on the optimal timing of flowering. Therefore, regulation of flowering time is of extraordinary importance to meet future food and feed demands. We developed the first barley nested association mapping (NAM) population, HEB-25, by crossing 25 wild barleys with one elite barley cultivar, and used it to dissect the genetic architecture of flowering time.

Results

Upon cultivation of 1,420 lines in multi-field trials and applying a genome-wide association study, eight major quantitative trait loci (QTL) were identified as main determinants to control flowering time in barley. These QTL accounted for 64% of the cross-validated proportion of explained genotypic variance (p_G). The strongest single QTL effect corresponded to the known photoperiod response gene Ppd-H1. After sequencing the causative part of Ppd-H1, we differentiated twelve haplotypes in HEB-25, whereof the strongest exotic haplotype accelerated flowering time by 11 days compared to the elite barley haplotype. Applying a whole genome prediction model including main effects and epistatic interactions allowed predicting flowering time with an unmatched accuracy of 77% of cross-validated p_G.

Conclusions

The elaborated causal models represent a fundamental step to explain flowering time in barley. In addition, our study confirms that the exotic biodiversity present in HEB-25 is a valuable toolbox to dissect the genetic architecture of important agronomic traits and to replenish the elite barley breeding pool with favorable, trait-improving exotic alleles.

     

Analysis of secondary phenotypes in multigroup association studies

Although case-control association studies have been widely used, they are insufficient for many complex diseases, such as Alzheimer's disease and breast cancer, since these diseases may have multiple subtypes with distinct morphologies and clinical implications. Many multigroup studies, such as the Alzheimer's Disease Neuroimaging Initiative (ADNI), have been undertaken by recruiting subjects based on their multiclass primary disease status, while extensive secondary outcomes have been collected. The aim of this paper is to develop a general regression framework for the analysis of secondary phenotypes collected in multigroup association studies. Our regression framework is built on a conditional model for the secondary outcome given the multigroup status and covariates and its relationship with the population regression of interest of the secondary outcome given the covariates. Then, we develop generalized estimation equations to estimate the parameters of interest. We use both simulations and a large-scale imaging genetic data analysis from the ADNI to evaluate the effect of the multigroup sampling scheme on standard genome-wide association analyses based on linear regression methods, while comparing it with our statistical methods that appropriately adjust for the multigroup sampling scheme. Data used in preparation of this article were obtained from the ADNI database.     

A metagenomic approach to dissect the genetic composition of enterotypes in Han Chinese and two Muslim groups

Distinct enterotypes have been observed in the human gut but little is known about the genetic basis of the microbiome. Moreover, it is not clear how many genetic differences exist between enterotypes within or between populations. In this study, both the 16S rRNA gene and the metagenomes of the gut microbiota were sequenced from 48 Han Chinese, 48 Kazaks, and 96 Uyghurs, and taxonomies were assigned after de novo assembly. Single nucleotide polymorphisms were also identified by referring to data from the Human Microbiome Project. Systematic analysis of the gut communities in terms of their abundance and genetic composition was also performed, together with a genome-wide association study of the host genomes. The gut microbiota of 192 subjects was clearly classified into two enterotypes (Bacteroides and Prevotella). Interestingly, both enterotypes showed a clear genetic differentiation in terms of their functional catalogue of genes, especially for genes involved in amino acid and carbohydrate metabolism. In addition, several differentiated genera and genes were found among the three populations. Notably, one human variant (rs878394) was identified that showed significant association with the abundance of Prevotella, which is linked to LYPLAL1, a gene associated with body fat distribution, the waist-hip ratio and insulin sensitivity. Taken together, considerable differentiation was observed in gut microbes between enterotypes and among populations that was reflected in both the taxonomic composition and the genetic makeup of their functional genes, which could have been influenced by a variety of factors, such as diet and host genetic variation.     

Overexpression of human Argonaute2 inhibits cell and tumor growth

Background

Argonaute (Ago) proteins are essential for the biogenesis and function of ~ 20–30 nucleotide long RNAs such as microRNAs (miRNAs). Ago expression increases or decreases under various physiological conditions, although the functional consequences are unknown. In addition, while reduced global miRNA production was shown to enhance cellular transformation and tumorigenesis, how Ago proteins contribute to human diseases has not been reported.

Method

Ago2, an essential Ago isoform in mammals, was stably expressed in 293 T, the human embryonic kidney cell line, and H1299, the human lung adenocarcinoma cell line. miRNA and mRNA expression was investigated by quantitative PCR and microarray profiling. Cell proliferation and migration was examined by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay and scratch assay in the cell cultures, respectively. How Ago2 affected cell growth in vivo was determined by H1299 xenograft tumor growth in mice. Changes in Ago2 expression in human lung cancer samples were investigated by quantitative PCR and immunohistochemistry.

Results

Stable Ago2 overexpression elicited specific changes in miRNA and mRNA expression in both 293 T and H1299 cells. It also inhibited cell proliferation and migration in cell cultures as well as xenograft tumor growth in nude mice. Ago2 expression was lower in human lung adenocarcinomas than in the paired, non-cancerous tissues.

General significance

We concluded that changes in Ago2 expression might have significant physiological and pathological consequences in vivo.     

基因组时代林木抗病分子机理研究的新进展

林木的分子病理学研究长期以来落后于农业作物病理学。随着高通量测序技术的问世,林木的分子病理学研究迎来了一个崭新的时代。从2006年至今,杨树、云杉等重要森林树种的全基因组测序相继完成,这为全面解析林木的抗病过程提供了遗传背景。同时,转录组学和全基因组关联分析的应用使得人们能快速地积累大量的数据,从而为揭示林木和病原菌之间的分子互作机制奠定了基础。近两年来CRISPR/Cas9基因编辑等分子生物学技术创新不断。高效的分子生物学技术结合基因组学研究有利于林木育种的研究。以下阐述了林木对抗病原菌入侵的生理机制,综合论述了近十年来基因组学和转录组学研究在木本植物分子病理学方面所取得的成果,总结了分子生物学技术在林木抗病领域的研究成果,分析了存在的问题和未来发展的趋势,以期为林木抗病育种提供参考。