Nearby Single Nucleotide Polymorphisms Research Articles

Structural variants contribute to phenotypic variation in maize.

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants, especially transposable elements are likely to affect phenotypic variation but we need better methods in maize for detecting polymorphic structural variants and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic structural variants and then genotyped a large maize diversity panel for these variants using short-read sequencing data. We characterized variation of SVs within the panel and identified SV polymorphisms that are associated with life history traits and genotype-by-environment interactions. While most of the SVs associated with traits contained TEs, only one of the SV's boundaries clearly matched TE breakpoints indicative of a TE insertion, whereas the other polymorphisms were likely caused by deletions. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that this method did not identify variants that would have been missed in a SNP association study.

A deep population reference panel of tandem repeat variation

Tandem repeats (TRs) represent one of the largest sources of genetic variation in humans and are implicated in a range of phenotypes. Here we present a deep characterization of TR variation based on high coverage whole genome sequencing from 3550 diverse individuals from the 1000 Genomes Project and H3Africa cohorts. We develop a method, EnsembleTR, to integrate genotypes from four separate methods resulting in high-quality genotypes at more than 1.7 million TR loci. Our catalog reveals novel sequence features influencing TR heterozygosity, identifies population-specific trinucleotide expansions, and finds hundreds of novel eQTL signals. Finally, we generate a phased haplotype panel which can be used to impute most TRs from nearby single nucleotide polymorphisms (SNPs) with high accuracy. Overall, the TR genotypes and reference haplotype panel generated here will serve as valuable resources for future genome-wide and population-wide studies of TRs and their role in human phenotypes.

Intergenic variants, rs1200114 and rs1200108 are genetically associated along with a decreased ATP1B1 expression in Fuchs Endothelial Corneal Dystrophy

Fuchs endothelial corneal dystrophy (FECD) is an age-related, bilateral corneal condition, characterized by apoptosis of the terminally differentiated endothelial cells. A genome-wide association study (GWAS) conducted in the European population in 2017, identified a new single nucleotide polymorphism (SNP), rs1200114 in the intergenic region between long intergenic non-protein coding RNA 970 (LINC00970) and ATPase Na+/K+ transporting subunit beta 1 (ATP1B1). The major focus of the current study is to understand the genetic association of this intergenic variant, rs1200114 with FECD in the Indian population. Sanger sequencing followed by statistical analysis indicated a significant difference in the allelic frequency between controls and cases (P = 0.01) with the minor allele ‘G’ of rs1200114 imparting a 1.64 fold increased risk for the disease. Luciferase reporter assay revealed no significant difference in the luciferase activity between allele ‘A’ and ‘G’ of rs1200114. However, quantitative RT-PCR assay revealed lower expression of ATP1B1 in FECD subjects compared with controls (P = 0.007). Therefore, to find whether another nearby SNP imparts regulatory effect, tag SNP association analysis was carried out; which revealed a significant association of another SNP, rs1200108, present in the intergenic region between LINC00970 and ATP1B1 with FECD (P = 0.009). The protective allele ‘A’ of rs1200108 displayed reduced reporter activity as opposed to the risk allele ‘G’ (P = 0.014). Furthermore, haplotype ‘A-A’ of rs1200108 - rs1200114 was present at a higher frequency in control subjects, suggesting it as a protective haplotype. Altogether, this study inferred the genetic association of rs1200114 and rs1200108 along with the decreased expression of ATP1B1 related to FECD pathogenesis in the Indian population.

Wavelet Screening: a novel approach to analyzing GWAS data

BackgroundTraditional methods for single-variant genome-wide association study (GWAS) incur a substantial multiple-testing burden because of the need to test for associations with a vast number of single-nucleotide polymorphisms (SNPs) simultaneously. Further, by ignoring more complex joint effects of nearby SNPs within a given region, these methods fail to consider the genomic context of an association with the outcome.ResultsTo address these shortcomings, we present a more powerful method for GWAS, coined ‘Wavelet Screening’ (WS), that greatly reduces the number of tests to be performed. This is achieved through the use of a sliding-window approach based on wavelets to sequentially screen the entire genome for associations. Wavelets are oscillatory functions that are useful for analyzing the local frequency and time behavior of signals. The signals can then be divided into different scale components and analyzed separately. In the current setting, we consider a sequence of SNPs as a genetic signal, and for each screened region, we transform the genetic signal into the wavelet space. The null and alternative hypotheses are modeled using the posterior distribution of the wavelet coefficients. WS is enhanced by using additional information from the regression coefficients and by taking advantage of the pyramidal structure of wavelets. When faced with more complex genetic signals than single-SNP associations, we show via simulations that WS provides a substantial gain in power compared to both the traditional GWAS modeling and another popular regional association test called SNP-set (Sequence) Kernel Association Test (SKAT). To demonstrate feasibility, we applied WS to a large Norwegian cohort (N=8006) with genotypes and information available on gestational duration.ConclusionsWS is a powerful and versatile approach to analyzing whole-genome data and lends itself easily to investigating various omics data types. Given its broader focus on the genomic context of an association, WS may provide additional insight into trait etiology by revealing genes and loci that might have been missed by previous efforts.

PIP-SNP: a pipeline for processing SNP data featured as linkage disequilibrium bin mapping, genotype imputing and marker synthesizing.

Genome-wide association study data analyses often face two significant challenges: (i) high dimensionality of single-nucleotide polymorphism (SNP) genotypes and (ii) imputation of missing values. SNPs are not independent due to physical linkage and natural selection. The correlation of nearby SNPs is known as linkage disequilibrium (LD), which can be used for LD conceptual SNP bin mapping, missing genotype inferencing and SNP dimension reduction. We used a stochastic process to describe the SNP signals and proposed two types of autocorrelations to measure nearby SNPs’ information redundancy. Based on the calculated autocorrelation coefficients, we constructed LD bins. We adopted a k-nearest neighbors algorithm (kNN) to impute the missing genotypes. We proposed several novel methods to find the optimal synthetic marker to represent the SNP bin. We also proposed methods to evaluate the information loss or information conservation between using the original genome-wide markers and using dimension-reduced synthetic markers. Our performance assessments on the real-life SNP data from a rice recombinant inbred line (RIL) population and a rice HapMap project show that the new methods produce satisfactory results. We implemented these functional modules in C/C++ and streamlined them into a web-based pipeline named PIP-SNP (https://bioinfo.noble.org/PIP_SNP/) for processing SNP data.

Genome-Wide Identification of Powdery Mildew Resistance in Common Bean (Phaseolus vulgaris L.).

Genome-wide association studies (GWAS) have been utilized to detect genetic variations related to several agronomic traits and disease resistance in common bean. However, its application in the powdery mildew (PM) disease to identify candidate genes and their location in the common bean genome has not been fully addressed. Single-nucleotide polymorphism (SNP) genotyping with a BeadChip containing 5398 SNPs was used to detect genetic variations related to PM disease resistance in a panel of 211 genotypes grown under two field conditions for two consecutive years. Significant SNPs identified on chromosomes Pv04 and Pv10 were repeatable, ensuring the phenotypic data’s reliability and the causal relationship. A cluster of resistance genes was revealed on the Pv04 of the common bean genome, coiled-coil-nucleotide-binding site–leucine-rich repeat (CC-NBS-LRR, CNL), and Toll/interleukin-1 receptor-nucleotide-binding site–leucine-rich repeat type (TIR-NBS-LRR, TNL)-like resistance genes were identified. Furthermore, two resistance genes, Phavu_010G1320001g and Phavu_010G136800g, were also identified on Pv10. Further sequence analysis showed that these genes were homologs to the disease-resistance protein (RLM1A-like) and the putative disease-resistance protein (At4g11170.1) in Arabidopsis. Significant SNPs related to two LRR receptor-like kinases (RLK) were only identified on Pv11 in 2018. Many genes encoding the auxin-responsive protein, TIFY10A protein, growth-regulating factor five-like, ubiquitin-like protein, and cell wall RBR3-like protein related to PM disease resistance were identified nearby significant SNPs. These results suggested that the resistance to PM pathogen involves a network of many genes constitutively co-expressed.

Kernel machine SNP set analysis finds the association of BUD13, ZPR1, and APOA5 variants with metabolic syndrome in Tehran Cardio-metabolic Genetics Study

Metabolic syndrome (MetS) is one of the most important risk factors for cardiovascular disease. The 11p23.3 chromosomal region plays a potential role in the pathogenesis of MetS. The present study aimed to assess the association between 18 single nucleotide polymorphisms (SNPs) located at the BUD13, ZPR1, and APOA5 genes with MetS in the Tehran Cardio-metabolic Genetics Study (TCGS). In 5421 MetS affected and non-affected participants, we analyzed the data using two models. The first model (MetS model) examined SNPs' association with MetS. The second model (HTg-MetS Model) examined the association of SNPs with MetS affection participants who had a high plasma triglyceride (TG). The four-gamete rules were used to make SNP sets from correlated nearby SNPs. The kernel machine regression models and single SNP regression evaluated the association between SNP sets and MetS. The kernel machine results showed two sets over three sets of correlated SNPs have a significant joint effect on both models (p < 0.0001). Also, single SNP regression results showed that the odds ratios (ORs) for both models are almost similar; however, the p-values had slightly higher significance levels in the HTg-MetS model. The strongest ORs in the HTg-MetS model belonged to the G allele in rs2266788 (MetS: OR = 1.3, p = 3.6 × 10–7; HTg-MetS: OR = 1.4, p = 2.3 × 10–11) and the T allele in rs651821 (MetS: OR = 1.3, p = 2.8 × 10–7; HTg-MetS: OR = 1.4, p = 3.6 × 10–11). In the present study, the kernel machine regression models could help assess the association between the BUD13, ZPR1, and APOA5 gene variants (11p23.3 region) with lipid-related traits in MetS and MetS affected with high TG.

A Genome-Wide Association Study Revealed Key SNPs/Genes Associated With Salinity Stress Tolerance In Upland Cotton.

Millions of hectares of land are too saline to produce economically valuable crop yields. Salt tolerance in cotton is an imperative approach for improvement in response to ever-increasing soil salinization. Little is known about the genetic basis of salt tolerance in cotton at the seedling stage. To address this issue, a genome-wide association study (GWAS) was conducted on a core collection of a genetically diverse population of upland cotton (Gossypium hirsutum L.) comprising of 419 accessions, representing various geographic origins, including China, USA, Pakistan, the former Soviet Union, Chad, Australia, Brazil, Mexico, Sudan, and Uganda. Phenotypic evaluation of 7 traits under control (0 mM) and treatment (150 mM) NaCl conditions depicted the presence of broad natural variation in the studied population. The association study was carried out with the efficient mixed-model association eXpedited software package. A total of 17,264 single-nucleotide polymorphisms (SNPs) associated with different salinity stress tolerance related traits were found. Twenty-three candidate SNPs related to salinity stress-related traits were selected. Final key SNPs were selected based on the r2 value with nearby SNPs in a linkage disequilibrium (LD) block. Twenty putative candidate genes surrounding SNPs, A10_95330133 and D10_61258588, associated with leaf relative water content, RWC_150, and leaf fresh weight, FW_150, were identified, respectively. We further validated the expression patterns of twelve candidate genes with qRT-PCR, which revealed different expression levels in salt-tolerant and salt-sensitive genotypes. The results of our GWAS provide useful knowledge about the genetic control of salt tolerance at the seedling stage, which could assist in elucidating the genetic and molecular mechanisms of salinity stress tolerance in cotton plants.

A novel association of rs13334070 in the RPGRIP1L gene with adiposity factors discovered by joint linkage and linkage disequilibrium analysis in Iranian pedigrees: Tehran Cardiometabolic Genetic Study (TCGS).

Understanding the genetic and metabolic bases of obesity is helpful in planning and developing health strategies. Therefore, the first family-based joint linkage and linkage disequilibrium study was conducted in Iranian pedigrees to assess the relationship between obesity and single-nucleotide polymorphisms (SNPs) located in the 16q12.2 region. In the present study, a total of 13,344 individuals were included, of whom 12,502 individuals were within 3,109 pedigrees and 842 were unrelated singletons. To investigate the relationship between obesity and genetic variants, a joint model of linkage and linkage disequilibrium was applied. Moreover, a sequence kernel association test (SKAT) was used to evaluate the association of the SNP set with body size and lipid profile measurements. The joint model showed that rs13334070, in the intron 4 of the RPGRIP1L gene, has a significant association with obesity. According to the 4-gamete rule, which is a procedure for constructing SNP sets by considering recombination occurrence between SNPs, this polymorphism has a high correlation with six nearby SNPs that make an SNP set. SKAT showed that this SNP set has a significant association with body size factors, but almost no association with most of the lipid profile measurements. In conclusion, from the result of this study, it might be reasonable to consider RPGRIP1L as an important gene whose variations could be associated with obesity risk factors.

Detection and analysis of CpG sites with multimodal DNA methylation level distributions and their relationships with SNPs

DNA methylation levels at cytosine-phosphate-guanine (CpG) sites with multimodal distributions among different samples have been reported recently. One possible explanation for such variability is that genetic variants might affect epigenetic variation. One obvious case is that mutations such as single-nucleotide polymorphisms (SNPs) interrupt CpG sites, resulting in different DNA methylation levels for different genotypes. However, the relationship between genetic variations and epigenetic differences has not been studied thoroughly, partially because of the lack of powerful and robust methods to survey genome-wide CpG sites with multimodal methylation level distributions (mmCpGs). In this article, we develop a Gaussian mixture-model clustering (GMMC)–based approach to systematically detect all mmCpGs across the genome based on the GAW20 data set. In total, 3785 and 3847 mmCpGs have been identified in pre- and posttreatment data sets, respectively. Result analysis shows that approximately 68 to 70% of mmCpGs detected from unrelated individuals either have direct overlaps with SNPs or have associations with nearby SNPs, suggesting a strong correlation between SNPs and mmCpGs. Comparison with an existing approach illustrates that our GMMC-based method is more consistent when the number of samples decreases. In conclusion, mmCpGs may reveal important connections between genetics and epigenetics and they should be carefully identified and evaluated.

Sodium valproate rescues expression of TRANK1 in iPSC-derived neural cells that carry a genetic variant associated with serious mental illness.

Biological characterization of genetic variants identified in genome-wide association studies (GWAS) remains a substantial challenge. Here we used human-induced pluripotent stem cells (iPSC) and their neural derivatives to characterize common variants on chromosome 3p22 that have been associated by GWAS with major mental illnesses. IPSC-derived neural progenitor cells carrying the risk allele of the single nucleotide polymorphism (SNP), rs9834970, displayed lower baseline TRANK1 expression that was rescued by chronic treatment with therapeutic dosages of valproic acid (VPA). VPA had the greatest effects on TRANK1 expression in iPSC, NPC, and astrocytes. Although rs9834970 has no known function, we demonstrated that a nearby SNP, rs906482, strongly affects binding by the transcription factor, CTCF, and that the high-affinity allele usually occurs on haplotypes carrying the rs9834970 risk allele. Decreased expression of TRANK1 perturbed expression of many genes involved in neural development and differentiation. These findings have important implications for the pathophysiology of major mental illnesses and the development of novel therapeutics.

Abstract 1221: Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic states in regulatory regions

Abstract Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic state in regulatory regions The vast majority of disease-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWASs) reside in noncoding genomic regions. Functional impact of noncoding casual variants on target gene expression often involves regulatory mechanisms, typically through altering local chromatin accessibility, histone modifications, transcription factor (TF) binding affinity and chromatin interactions. In chronic lymphocytic leukemia (CLL), we and others have identified 41 GWAS risk variants, whose mechanism of action in CLL etiology remains largely unknown. In this study, we aim to decipher the regulatory potential of these variants and to better understand the mechanisms contributing to dysregulation of genes implicated in CLL pathogenesis. We performed ChIP-seq for H3K4me3, H3K4me1, H3K27ac and H3K27me3 in CLL. In addition, we also analyzed over 250 public ChIP-seq, ATAC-seq, DNase-seq, FAIRE-seq, and chromatin interaction data in CLL, B cells and/or lymphoblastoid cell lines. Comparing the 41 index SNPs and the nearby SNPs they associated with in linkage disequilibrium (R2 at least 0.5 in the EUR ethnic group in 1000 Genomes Project) to the above epigenetic data sets revealed overlaps with regulatory regions in 38 of the risk loci, with 28 of the loci overlapping super-enhancers. We used two approaches to elucidate the mechanisms of potential SNP causality. First, we used GATK UnifiedGenotyper to infer genotype for the SNPs that overlapped open chromatin regions and ChIP-seq peaks, by which evidence of imbalance between reference and alternative allele was identified at heterozygous sites. Second, for SNPs located within open chromatin regions or ChIP-seq peaks, we used FIMO program to scan the 100bp sequences centered on the SNPs for significant matches (p=1e-05) to the position weight matrices in five TF binding motif databases; SNPs that had TF motifs altered by the alternative alleles were retained. Together, these two approaches revealed evidence for allele imbalance or TF binding motif alterations at 36 of the 38 CLL risk loci, suggesting the generality of epigenetic regulation in CLL pathogenesis. Further analysis of chromatin interaction, gene expression and expression quantitative trait loci data will likely link the casual SNPs to their target genes. Citation Format: Huihuang Yan, Shulan Tian, Geffen Kleinstern, Susan Slager. Chronic lymphocytic leukemia risk variants are associated with alteration of epigenetic states in regulatory regions [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1221.

Generation of two H1 hESC sublines carrying a heterozygous and homozygous knock-out of RB1

Retinoblastoma is a childhood cancer of the retina caused by biallelic inactivation of the tumor suppressor gene RB1. In heritable retinoblastoma, one allele is inherited in mutant form via one of the parental germ cells. To study molecular mechanisms in retinoblastoma, two sublines of H1 hESCs were generated, carrying a knock-out allele of RB1 in the heterozygous or homozygous state. Exon 3 of RB1 was targeted and modified by nucleotide deletions using the CRISPR/Cas9 nuclease system. Based on a nearby single nucleotide polymorphism, the modification could be assigned to one allele.

Empirical comparison of reduced representation bisulfite sequencing and Infinium BeadChip reproducibility and coverage of DNA methylation in humans

We empirically examined the strengths and weaknesses of two human genome-wide DNA methylation platforms: rapid multiplexed reduced representation bisulfite sequencing and Illumina’s Infinium BeadChip. Rapid multiplexed reduced representation bisulfite sequencing required less input DNA, offered more flexibility in coverage, and interrogated more CpG loci at a higher regional density. The Infinium covered slightly more protein coding, cancer-associated and mitochondrial-related genes, both platforms covered all known imprinting clusters, and rapid multiplexed reduced representation bisulfite sequencing covered more microRNA genes than the HumanMethylation450, but fewer than the MethylationEPIC. Rapid multiplexed reduced representation bisulfite sequencing did not always interrogate exactly the same CpG loci, but genomic tiling improved overlap between different libraries. Reproducibility of rapid multiplexed reduced representation bisulfite sequencing and concordance between the platforms increased with CpG density. Only rapid multiplexed reduced representation bisulfite sequencing could genotype samples and measure allele-specific methylation, and we confirmed that Infinium measurements are influenced by nearby single-nucleotide polymorphisms. The respective strengths and weaknesses of these two genome-wide DNA methylation platforms need to be considered when conducting human epigenetic studies.

Parkinson's disease-associated genetic variation is linked to quantitative expression of inflammatory genes.

Genome-wide association studies (GWAS) have linked dozens of single nucleotide polymorphisms (SNPs) with Parkinson’s disease (PD) risk. Ascertaining the functional and eventual causal mechanisms underlying these relationships has proven difficult. The majority of risk SNPs, and nearby SNPs in linkage disequilibrium (LD), are found in intergenic or intronic regions and confer risk through allele-dependent expression of multiple unknown target genes. Combining GWAS results with publicly available GTEx data, generated through eQTL (expression quantitative trait loci) identification studies, enables a direct association of SNPs to gene expression levels and aids in narrowing the large population of potential genetic targets for hypothesis-driven experimental cell biology. Separately, overlapping of SNPs with putative enhancer segmentations can strengthen target filtering. We report here the results of analyzing 7,607 PD risk SNPs along with an additional 23,759 high linkage disequilibrium-associated variants paired with eQTL gene expression. We found that enrichment analysis on the set of genes following target filtering pointed to a single large LD block at 6p21 that contained multiple HLA-MHC-II genes. These MHC-II genes remain associated with PD when the genes were filtered for correlation between GWAS significance and eQTL levels, strongly indicating a direct effect on PD etiology.

Selecting sequence variants to improve genomic predictions for dairy cattle

BackgroundMillions of genetic variants have been identified by population-scale sequencing projects, but subsets of these variants are needed for routine genomic predictions or genotyping arrays. Methods for selecting sequence variants were compared using simulated sequence genotypes and real July 2015 data from the 1000 Bull Genomes Project.MethodsCandidate sequence variants for 444 Holstein animals were combined with high-density (HD) imputed genotypes for 26,970 progeny-tested Holstein bulls. Test 1 included single nucleotide polymorphisms (SNPs) for 481,904 candidate sequence variants. Test 2 also included 249,966 insertions-deletions (InDels). After merging sequence variants with 312,614 HD SNPs and editing steps, Tests 1 and 2 included 762,588 and 1,003,453 variants, respectively. Imputation quality from findhap software was assessed with 404 of the sequenced animals in the reference population and 40 randomly chosen animals for validation. Their sequence genotypes were reduced to the subset of genotypes that were in common with HD genotypes and then imputed back to sequence. Predictions were tested for 33 traits using 2015 data of 3983 US validation bulls with daughters that were first phenotyped after August 2011.ResultsThe average percentage of correctly imputed variants across all chromosomes was 97.2 for Test 1 and 97.0 for Test 2. Total time required to prepare, edit, impute, and estimate the effects of sequence variants for 27,235 bulls was about 1 week using less than 33 threads. Many sequence variants had larger estimated effects than nearby HD SNPs, but prediction reliability improved only by 0.6 percentage points in Test 1 when sequence SNPs were added to HD SNPs and by 0.4 percentage points in Test 2 when sequence SNPs and InDels were included. However, selecting the 16,648 candidate SNPs with the largest estimated effects and adding them to the 60,671 SNPs used in routine evaluations improved reliabilities by 2.7 percentage points.ConclusionsReliabilities for genomic predictions improved when selected sequence variants were added; gains were similar for simulated and real data for the same population, and larger than previous gains obtained by adding HD SNPs. With many genotyped animals, many data sources, and millions of variants, computing strategies must efficiently balance costs of imputation, selection, and prediction to obtain subsets of markers that provide the highest accuracy.

Estrogen-dependent association of HDAC4 with fear in female mice and women with PTSD

Women are at increased risk of developing post-traumatic stress disorder (PTSD) following a traumatic event. Recent studies suggest that this may be mediated, in part, by circulating estrogen levels. This study evaluated the hypothesis that individual variation in response to estrogen levels contributes to fear regulation and PTSD risk in women. We evaluated DNA methylation from blood of female participants in the Grady Trauma Project and found that serum estradiol levels associates with DNA methylation across the genome. For genes expressed in blood, we examined the association between each CpG site and PTSD diagnosis using linear models that adjusted for cell proportions and age. After multiple test correction, PTSD associated with methylation of CpG sites in the HDAC4 gene, which encodes histone deacetylase 4, and is involved in long-term memory formation and behavior. DNA methylation of HDAC4 CpG sites were tagged by a nearby single-nucleotide polymorphism (rs7570903), which also associated with HDAC4 expression, fear-potentiated startle and resting-state functional connectivity of the amygdala in traumatized humans. Using auditory Pavlovian fear conditioning in a rodent model, we examined the regulation of Hdac4 in the amygdala of ovariectomized (OVX) female mice. Hdac4 messenger RNA levels were higher in the amygdala 2 h after tone-shock presentations, compared with OVX-homecage control females. In naturally cycling females, tone-shock presentations increased Hdac4 expression relative to homecage controls for metestrous (low estrogen) but not the proestrous (high estrogen) group. Together, these results support an estrogenic influence of HDAC4 regulation and expression that may contribute to PTSD in women.

Visual Assessment of Chest Computed Tomographic Images Is Independently Useful for Genetic Association Analysis in Studies of Chronic Obstructive Pulmonary Disease.

Automated analysis of computed tomographic (CT) lung images for epidemiologic and genetic association studies is increasingly common, but little is known about the utility of visual versus semiautomated emphysema and airway assessments for genetic association studies. Assess the relative utility of visual versus semiautomated emphysema and airway assessments for genetic association studies. A standardized inspection protocol was used to visually assess chest CT images for 1,540 non-Hispanic white subjects within the COPDGene Study for the presence and severity of radiographic features representing airway wall thickness and emphysema. A genome-wide association study (GWAS) was performed, and two sets of candidate single-nucleotide polymorphisms with a higher prior likelihood of association were specified a priori for separate analysis. For each visual CT examination feature, a corresponding semiautomated CT feature(s) was identified for comparison in the same subjects. GWAS for visual features of chest CT scans identified a genome-wide significant association with visual emphysema at the 15q25 locus (P = 6.3e-9). In the a priori-specified set of 19 previously identified GWAS loci, 7 and 8 loci were associated with airway measures or emphysema measures, respectively. In the a priori-specified candidate gene set, 13 of 196 candidate genes harbored a nearby single-nucleotide polymorphism significantly associated with an emphysema phenotype. Visual CT examination associations were robust to adjustment for semiautomated correlates in many cases. Standardized visual assessments of emphysema and airway disease are significantly associated with genetic loci previously associated with chronic obstructive pulmonary disease susceptibility or semiautomated CT examination phenotypes in GWAS. Visual CT measures of emphysema and airways disease offer independent information for genetic association studies in relation to standard semiautomated measures.

KCNK3 Variants Are Associated With Hyperaldosteronism and Hypertension.

Blood pressure (BP) is a complex trait that is the consequence of an interaction between genetic and environmental determinants. Previous studies have demonstrated increased BP in mice with global deletion of TASK-1 channels contemporaneous with diverse dysregulation of aldosterone production. In humans, genome-wide association studies in ≈100 000 individuals of European, East Asian, and South Asian ancestry identified a single nucleotide polymorphism (SNP) in KCNK3 (the gene encoding TASK-1) associated with mean arterial pressure. The current study was motivated by the hypotheses that (1) association of KCNK3 SNPs with BP and related traits extends to blacks and Hispanics, and (2) KCNK3 SNPs exhibit associations with plasma renin activity and aldosterone levels. We examined baseline BP measurements for 7840 participants from the Multi-Ethnic Study of Atherosclerosis (MESA), and aldosterone levels and plasma renin activity in a subset of 1653 MESA participants. We identified statistically significant association of the previously reported KCNK3 SNP (rs1275988) with mean arterial pressure in MESA blacks (P=0.024) and a nearby SNP (rs13394970) in MESA Hispanics (P=0.031). We discovered additional KCNK3 SNP associations with systolic BP, mean arterial pressure, and hypertension. We also identified statistically significant association of KCNK3 rs2586886 with plasma aldosterone level in MESA and demonstrated that global deletion of TASK-1 channels in mice produces a mild-hyperaldosteronism, not associated with a decrease in renin. Our results suggest that genetic variation in the KCNK3 gene may contribute to BP variation and less severe hypertensive disorders in which aldosterone may be one of several causative factors.

Blood and Intestine eQTLs from an Anti-TNF-Resistant Crohn's Disease Cohort Inform IBD Genetic Association Loci.

OBJECTIVES:Genome-wide association studies (GWAS) have identified loci reproducibly associated with inflammatory bowel disease (IBD) and other immune-mediated diseases; however, the molecular mechanisms underlying most of genetic susceptibility remain undefined. Expressional quantitative trait loci (eQTL) of disease-relevant tissue can be employed in order to elucidate the genes and pathways affected by disease-specific genetic variance.METHODS:In this study, we derived eQTLs for human whole blood and intestine tissues of anti-tumor necrosis factor-resistant Crohn's disease (CD) patients. We interpreted these eQTLs in the context of published IBD GWAS hits to inform on the disease process.RESULTS:At 10% false discovery rate, we discovered that 5,174 genes in blood and 2,063 genes in the intestine were controlled by a nearby single-nucleotide polymorphism (SNP) (i.e., cis-eQTL), among which 1,360 were shared between the two tissues. A large fraction of the identified eQTLs were supported by the regulomeDB database, showing that the eQTLs reside in regulatory elements (odds ratio; OR=3.44 and 3.24 for blood and intestine eQTLs, respectively) as opposed to protein-coding regions. Published IBD GWAS hits as a whole were enriched for blood and intestine eQTLs (OR=2.88 and 2.05; and P value=2.51E-9 and 0.013, respectively), thereby linking genetic susceptibility to control of gene expression in these tissues. Through a systematic search, we used eQTL data to inform 109 out of 372 IBD GWAS SNPs documented in National Human Genome Research Institute catalog, and we categorized the genes influenced by eQTLs according to their functions. Many of these genes have experimentally validated roles in specific cell types contributing to intestinal inflammation.CONCLUSIONS:The blood and intestine eQTLs described in this study represent a powerful tool to link GWAS loci to a regulatory function and thus elucidate the mechanisms underlying the genetic loci associated with IBD and related conditions. Overall, our eQTL discovery approach empirically identifies the disease-associated variants including their impact on the direction and extent of expression changes in the context of disease-relevant cellular pathways in order to infer the functional outcome of this aspect of genetic susceptibility.