Phenotypic correlations of carpal gland diverticular number with production traits and its genome-wide association analysis in multiple pig populations

Pig carpal glands play crucial roles in territorial recognition, reproductive behavior, and information exchange; however, their effects on production traits and underlying genetic mechanisms remain unclear. In this study, 1028 pigs from six populations were counted for the carpal gland diverticular numbers (CGDNs) on the left (CGDNL) and right (CGDNR) legs, and their carcass and meat quality traits were assessed. The CGDNs were significantly different among the populations, and Licha Black pigs had a lower CGDN than the Bama Xiang breed. It was also significantly different between sexes, with males having more diverticula than females (p ≤ 0.0391). Moreover, the number was asymmetric, with CGDNR being significantly higher than CGDNL. Notably, CGDNs was significantly correlated with each other in phenotype and genetics and with 24-h pH, 24-h meat color score, 24-h marbling score, fat content, moisture content, sodium salt content, and saturated fatty acid content in phenotype. Furthermore, genome-wide association analyses identified seven SNPs in association with CGDNs at a 5% genome-wide significance level, all of which were located in a 1.78-Mb (35.347-37.129 Mb) region on chromosome 1. CNC10010837 and CNC10010840 were the top SNPs: both had an additive effect of 0.789 ± 0.120 on CGDNR with p = 8.31E-10. These findings provide important insights into the functions and underlying genetic mechanisms of swine carpal glands.

Genome-wide association study identifies a novel candidate gene for egg production traits in chickens

Qingyuan partridge chicken is a renowned indigenous yellow broiler breed in China. Egg production traits are important economic traits for chickens. With the decreasing cost of whole genome resequencing, identifying candidate genes with more precision has become possible. In order to identify molecular markers and candidate genes associated with egg production traits, we conducted genome-wide association studies based on the resequencing data of 287 female Qingyuan partridge chickens. For each hen, age at first egg and egg laying rate were recorded and calculated, respectively. With a univariate linear mixed model, we detected one genome-wide significant single nucleotide polymorphism (SNP) and three chromosome-wide significant SNPs associated with egg laying rate. MTA2 is highly likely to be a functional gene for egg laying rate. Our study identifies MTA2 as the first time to be associated with egg laying rate. Findings in our study will advance our understanding of the genetic basis of egg production and have the potential to improve the efficiency of genomic selection in chickens.

Genetics in neovascular age-related macular degeneration susceptibility and treatment response to anti-VEGF intravitreal injection: A case series study

BACKGROUND

To identify genotypes associated with neovascular age-related macular degeneration (nAMD) and investigate the associations between genotype variations and anti-vascular endothelial growth factor (VEGF) treatment response.

METHODS

This observational, retrospective, case series study enrolled patients diagnosed with nAMD who received anti-VEGF treatment in National Taiwan University Hospital with at least one-year follow-up between 2012 and 2020. A genome-wide association study (GWAS) was conducted on enrolled patients and controls. Correlations between the genotypes identified from GWAS and the treatment response of functional/anatomical biomarkers, including visual acuity (VA), presence of intraretinal or subretinal fluid (SRF), serous or fibrovascular pigmented epithelium detachment (PED), and disruption of the ellipsoid zone (EZ), were analysed.

RESULTS

In total, 182 patients with nAMD and 1748 controls were enrolled. GWAS revealed 16 single nucleotide polymorphisms (SNPs) as risk loci for nAMD, including seven loci in CFH and ARMS2/HTRA1 and nine novel loci, including rs117517872 and rs79835234(COPB2-DT), rs7525578(RAP1A), rs2123738(LOC105376755), rs1374879(CNTN3), rs3812692(SAR1A), rs117501587(PRKCA), rs9965945(CNDP1), and rs189769231(MATK). Our study revealed rs800292(CFH), rs11200638(HTRA1), and rs2123738(LOC105376755) correlated with poor treatment response in VA (P = 0.005), SRF (P = 0.044), and fibrovascular PED (P = 0.007), respectively. Rs9965945(CNDP1) was correlated with poor response in disruption of EZ (P = 0.046) and serous PED (P = 0.049).

CONCLUSIONS

Among the 16 SNPs found in the GWAS, four loci-CFH, ARMS2/HTRA1, and two novel loci-were correlated with the susceptibility of nAMD and anatomical/functional responses after anti-VEGF treatment.

R2R3-MYB transcription factor CsMYB60 controls mature fruit skin color by regulating flavonoid accumulation in cucumber

Skin color is an important trait that determines the cosmetic appearance and quality of fruits. In cucumber, the skin color ranges from white to brown in mature fruits. However, the genetic basis for this important trait remains unclear. We conducted a genome-wide association study of natural cucumber populations, along with map-based cloning techniques, on an F2 population resulting from a cross between Pepino (with yellow-brown fruit skin) and Zaoer-N (with creamy fruit skin). We identified CsMYB60 as a candidate gene responsible for skin coloration in mature cucumber fruits. In cucumber accessions with white to pale yellow skin color, a premature stop mutation (C to T) was found in the second exon region of CsMYB60, whereas light yellow cucumber accessions exhibited splicing premature termination caused by an intronic mutator-like element insertion in CsMYB60. Transgenic CsMYB60c cucumber plants displayed a yellow-brown skin color by promoting accumulation of flavonoids, especially hyperoside, a yellow-colored flavonol. CsMYB60c encodes a nuclear protein that primarily acts as a transcriptional activator through its C-terminal activation motif. RNA sequencing and DNA affinity purification sequencing assays revealed that CsMYB60c promotes skin coloration by directly binding to the YYTACCTAMYT motif in the promoter regions of flavonoid biosynthetic genes, including CsF3'H, which encodes flavonoid 3'-hydroxylase. The findings of our study not only offer insight into the function of CsMYB60 as dominantly controlling fruit coloration, but also highlight that intronic DNA mutations can have a similar phenotypic impact as exonic mutations, which may be valuable in future cucumber breeding programs.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Genetically proxied intestinal microbiota and risk of erectile dysfunction

BACKGROUND

The interaction between intestinal microbiota and erectile dysfunction (ED) is less investigated. This study was performed to explore the association between intestinal microbiota and ED.

METHODS

In this two-sample Mendelian randomization (MR) study, genetic variants of gut microbiota were obtained from MiBioGen consortium containing 18,340 individuals. Six methods including inverse variance weighting (IVW), MR-Egger, weighted median, maximum likelihood, MR robust adjusted profile score, and MR pleiotropy residual sum and outlier were used to investigate the causal links between intestinal microbiota and ED. Furthermore, reverse MR analysis was performed to exclude the causal impact of ED on gut microbiota.

RESULTS

As revealed by the IVW estimator, the risks of ED were raised by genetically proxied Lachnospiraceae (OR: 1.27), Lachnospiraceae NC2004 group (OR: 1.17), Oscillibacter (OR: 1.20), Senegalimassilia (OR: 1.32) (All P < 0.05) and Tyzzerella-3 (OR: 1.14, P < 0.05). It was observed that Ruminococcaceae UCG013 exerted protective effect against ED (OR: 0.77, P < 0.05). These results were consistent with other estimators in sensitivity analyses. In reverse MR analyses, genetic liability to ED did not alter the abundances of Lachnospiraceae, Lachnospiraceae NC2004 group, Oscillibacter, Senegalimassilia, Tyzzerella-3, and Ruminococcaceae UCG013 (All P > 0.05). No heterogeneity and pleiotropy were detected by Cochran's Q-test, MR-Egger, and global test (All P > 0.05).

CONCLUSIONS

This study provided novel evidence that genetically proxied Lachnospiraceae, Lachnospiraceae NC2004 group, Oscillibacter, Senegalimassilia, Tyzzerella-3, and Ruminococcaceae UCG013 had potentially causal effects on ED. Further studies are needed to clarify the biological mechanisms linking intestinal microbiota to ED.

Flowering time regulator qFT13-3 involved in soybean adaptation to high latitudes

Soybean is a short-day plant that typically flowers earlier when exposed to short-day conditions. However, the identification of genes associated with earlier flowering time but without a yield penalty is rare. In this study, we conducted genome-wide association studies (GWAS) using two re-sequencing datasets that included 113 wild soybeans (G. soja) and 1192 cultivated soybeans (G. max), respectively, and simultaneously identified a candidate flowering gene, qFT13-3, which encodes a protein homologous to the pseudo-response regulator (PRR) transcription factor. We identified four major haplotypes of qFT13-3 in the natural population, with haplotype H4 (qFT13-3H4) being lost during domestication, while qFT13-3H1 underwent natural and artificial selection, increasing in proportion from 4.5% in G. soja to 43.8% in landrace and to 81.9% in improve cultivars. Notably, most cultivars harbouring qFT13-3H1 were located in high-latitude regions. Knockout of qFT13-3 accelerated flowering and maturity time under long-day conditions, indicating that qFT13-3 functions as a flowering inhibitor. Our results also showed that qFT13-3 directly downregulates the expression of GmELF3b-2 which is a component of the circadian clock evening complex. Field trials revealed that the qft13-3 mutants shorten the maturity period by 11 days without a concomitant penalty on yield. Collectively, qFT13-3 can be utilized for the breeding of high-yield cultivars with a short maturity time suitable for high latitudes.

Investigation of phenotypic, genetic and genomic background of Milk spectra in Sarda dairy sheep

Aim of this study was to analyse the genetic background of milk Fourier transform infrared (FTIR) spectra in dairy sheep. Individual milk FTIR spectra, with 1060 wavenumbers each, were available for 793 adult Sarda breed ewes genotyped at 45,813 SNP. The absorbance values of each wavenumber was analysed using a linear mixed model that included dim class, parity and lambing month as fixed effects and flock-test date and animal as random effects. The model was applied to estimate variance components and heritability and to perform a genome-wide association study for each wavenumber. Average h2 of wavenumbers absorbance was 0.13 ± 0.08, with the largest values observed in the regions associated with the characteristic bonds of carbonylic and methylenic groups of milk fat (h2 = 0.57 at 1724-1728 cm-1; and h2 = 0.34 at 2811-2834 cm-1, respectively). The absorbance values of wavenumbers were moderately correlated with the estimated heritabilities. After the Bonferroni correction, a total of nine markers were found to be significantly associated with 32 different wavenumbers. Of particular interest was the SNP s63269.1, mapped on chromosome 2, that was found to be associated with 27 wavenumbers. Genes previously found to be related to traits of interest (e.g. disease resistance, milk yield and quality, cheese firmness) are located close to the significant SNP. As expected, the heritability estimated for the absorbance of each wavenumbers seems to be associated with the related milk components.

A novel application of data-consistent inversion to overcome spurious inference in genome-wide association studies

The genome-wide association studies (GWAS) typically use linear or logistic regression models to identify associations between phenotypes (traits) and genotypes (genetic variants) of interest. However, the use of regression with the additive assumption has potential limitations. First, the normality assumption of residuals is the one that is rarely seen in practice, and deviation from normality increases the Type-I error rate. Second, building a model based on such an assumption ignores genetic structures, like, dominant, recessive, and protective-risk cases. Ignoring genetic variants may result in spurious conclusions about the associations between a variant and a trait. We propose an assumption-free model built upon data-consistent inversion (DCI), which is a recently developed measure-theoretic framework utilized for uncertainty quantification. This proposed DCI-derived model builds a nonparametric distribution on model inputs that propagates to the distribution of observed data without the required normality assumption of residuals in the regression model. This characteristic enables the proposed DCI-derived model to cover all genetic variants without emphasizing on additivity of the classic-GWAS model. Simulations and a replication GWAS with data from the COPDGene demonstrate the ability of this model to control the Type-I error rate at least as well as the classic-GWAS (additive linear model) approach while having similar or greater power to discover variants in different genetic modes of transmission.

Site-directed genotype screening for elimination of antinutritional saponins in quinoa seeds identifies TSARL1 as a master controller of saponin biosynthesis selectively in seeds

Climate change may result in a drier climate and increased salinization, threatening agricultural productivity worldwide. Quinoa (Chenopodium quinoa) produces highly nutritious seeds and tolerates abiotic stresses such as drought and high salinity, making it a promising future food source. However, the presence of antinutritional saponins in their seeds is an undesirable trait. We mapped genes controlling seed saponin content to a genomic region that includes TSARL1. We isolated desired genetic variation in this gene by producing a large mutant library of a commercial quinoa cultivar and screening the library for specific nucleotide substitutions using droplet digital PCR. We were able to rapidly isolate two independent tsarl1 mutants, which retained saponins in the leaves and roots for defence, but saponins were undetectable in the seed coat. We further could show that TSARL1 specifically controls seed saponin biosynthesis in the committed step after 2,3-oxidosqualene. Our work provides new important knowledge on the function of TSARL1 and represents a breakthrough for quinoa breeding.

Genome-wide association study reveals the genetic basis for petal-size formation in rapeseed (Brassica napus) and CRISPR-Cas9-mediated mutagenesis of BnFHY3 for petal-size reduction

Petals in rapeseed (Brassica napus) serve multiple functions, including protection of reproductive organs, nutrient acquisition, and attraction of pollinators. However, they also cluster densely at the top, forming a thick layer that absorbs and reflects a considerable amount of photosynthetically active radiation. Breeding genotypes with large, small, or even petal-less varieties, requires knowledge of primary genes for allelic selection and manipulation. However, our current understanding of petal-size regulation is limited, and the lack of markers and pre-breeding materials hinders targeted petal-size breeding. Here, we conducted a genome-wide association study on petal size using 295 diverse accessions. We identified 20 significant single nucleotide polymorphisms and 236 genes associated with petal-size variation. Through a cross-analysis of genomic and transcriptomic data, we focused on 14 specific genes, from which molecular markers for diverging petal-size features can be developed. Leveraging CRISPR-Cas9 technology, we successfully generated a quadruple mutant of Far-Red Elongated Hypocotyl 3 (q-bnfhy3), which exhibited smaller petals compared to the wild type. Our study provides insights into the genetic basis of petal-size regulation in rapeseed and offers abundant potential molecular markers for breeding. The q-bnfhy3 mutant unveiled a novel role of FHY3 orthologues in regulating petal size in addition to previously reported functions.

An Emerging Approach of Age-Related Hearing Loss Research: Application of Integrated Multi-Omics Analysis

As one of the most common otologic diseases in the elderly, age-related hearing loss (ARHL) usually characterized by hearing loss and cognitive disorders, which have a significant impact on the elderly's physical and mental health and quality of life. However, as a typical disease of aging, it is unclear why aging causes widespread hearing impairment in the elderly. As molecular biological experiments have been conducted for research recently, ARHL is gradually established at various levels with the application and development of integrated multi-omics analysis in the studies of ARHL. Here, the recent progress in the application of multi-omics analysis in the molecular mechanisms of ARHL development and therapeutic regimens, including the combined analysis of different omics, such as transcriptome, proteome, and metabolome, to screen for risk sites, risk genes, and differences in lipid metabolism, etc., is outlined and the integrated histological data further promote the profound understanding of the disease process as well as physiological mechanisms of ARHL. The advantages and disadvantages of multi-omics analysis in disease research are also discussed and the authors speculate on the future prospects and applications of this part-to-whole approach, which may provide more comprehensive guidance for ARHL and aging disease prevention and treatment.

An inverse causal relationship between serum 25-hydroxyvitamin D levels and pulmonary hypertension: A two-sample Mendelian randomization study

Observational studies have confirmed that 25-hydroxyvitamin D (25(OH)D) is associated with pulmonary hypertension (PH), but the causal association between each other is unclear. Therefore, Mendelian randomization (MR) method was performed to validate the causal association between PH and serum 25(OH)D levels. The summary data for 25(OH)D and PH were from the National Human Genome Research Institute-European Bioinformatics Institute. Catalog of human genome-wide association studies and FinnGen biobank consortium. MR analysis was utilized to explore the potential causal association between PH and 25(OH)D. To evaluate this association, inverse variance weighting was considered as the primary method. Cochran's Q test, MR-Egger intercept test, and "leave-one-out" sensitivity analyses were utilized to control the pleiotropy and heterogeneity in the study. Two-sample MR analysis revealed an inverse causal relationship between 25(OH)D and PH (odds ratio: 0.376, 95% confidence interval: 0.162-0.876, p = 2.334 × 10-2). There was no significant heterogeneity and pleiotropy. The present study confirmed the inverse causal relationship between 25(OH)D and PH. This pathway may provide another treatment pathway in PH. Further studies to elucidate this pathway is indicated.

Genome-wide meta-analysis and fine-mapping prioritize potential causal variants and genes related to leprosy

To date, genome-wide association studies (GWASs) have discovered 35 susceptible loci of leprosy; however, the cumulative effects of these loci can only partially explain the overall risk of leprosy, and the causal variants and genes within these loci remain unknown. Here, we conducted out new GWASs in two independent cohorts of 5007 cases and 4579 controls and then a meta-analysis in these newly generated and multiple previously published (2277 cases and 3159 controls) datasets were performed. Three novel and 15 previously reported risk loci were identified from these datasets, increasing the known leprosy risk loci of explained genetic heritability from 23.0 to 38.5%. A comprehensive fine-mapping analysis was conducted, and 19 causal variants and 14 causal genes were identified. Specifically, manual checking of epigenomic information from the Epimap database revealed that the causal variants were mainly located within the immune-relevant or immune-specific regulatory elements. Furthermore, by using gene-set, tissue, and cell-type enrichment analyses, we highlighted the key roles of immune-related tissues and cells and implicated the PD-1 signaling pathways in the pathogenetic mechanism of leprosy. Collectively, our study identified candidate causal variants and elucidated the potential regulatory and coding mechanisms for genes associated with leprosy.

Genome-Wide Association Study of Cardiovascular Resilience Identifies Protective Variation in the CETP Gene

Background The risk of atherosclerotic cardiovascular disease (ASCVD) increases sharply with age. Some older individuals, however, remain unaffected despite high predicted risk. These individuals may carry cardioprotective genetic variants that contribute to resilience. Our aim was to assess whether asymptomatic older individuals without prevalent ASCVD carry cardioprotective genetic variants that contribute to ASCVD resilience. Methods and Results We performed a genome-wide association study using a 10-year predicted ASCVD risk score as a quantitative trait, calculated only in asymptomatic older individuals aged ≥70 years without prevalent ASCVD. Our discovery genome-wide association study of N=12 031 ASCVD event-free individuals from the ASPREE (Aspirin in Reducing Events in the Elderly) trial identified 2 independent variants, rs9939224 (P<5×10-8) and rs56156922 (P<10-6), in the CETP (cholesteryl ester transfer protein) gene. The CETP gene is a regulator of plasma high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and lipoprotein(a) levels, and it is a therapeutic drug target. The associations were replicated in the UK Biobank (subpopulation of N=13 888 individuals aged ≥69 years without prevalent ASCVD). Carriers of the identified CETP variants (versus noncarriers) had higher plasma high-density lipoprotein cholesterol levels, lower plasma low-density lipoprotein cholesterol levels, and reduced risk of incident ASCVD events during follow-up. Expression quantitative trait loci analysis predicted the identified CETP variants reduce CETP gene expression across various tissues. Previously reported associations between genetic CETP inhibition and increased risk of age-related macular degeneration were not observed among the 3917 ASPREE trial participants with retinal imaging and genetic data available. Conclusions Common genetic variants in the CETP gene region are associated with cardiovascular resilience during aging. Registration URL: https://www.clinicaltrials.gov; Unique identifier: NCT01038583.

Genome-Wide Association Study of Pericardial Fat Area in 28 161 UK Biobank Participants

BACKGROUND Pericardial adipose tissue (PAT) is the visceral adipose tissue compartment surrounding the heart. Experimental and observational research has suggested that greater PAT deposition might mediate cardiovascular disease, independent of general or subcutaneous adiposity. We characterize the genetic architecture of adiposity-adjusted PAT and identify causal associations between PAT and adverse cardiac magnetic resonance imaging measures of cardiac structure and function in 28 161 UK Biobank participants. METHODS AND RESULTS The PAT phenotype was extracted from cardiac magnetic resonance images using an automated image analysis tool previously developed and validated in this cohort. A genome-wide association study was performed with PAT area set as the phenotype, adjusting for age, sex, and other measures of obesity. Functional mapping and Bayesian colocalization were used to understand the biologic role of identified variants. Mendelian randomization analysis was used to examine potential causal links between genetically determined PAT and cardiac magnetic resonance-derived measures of left ventricular structure and function. We discovered 12 genome-wide significant variants, with 2 independent sentinel variants (rs6428792, P=4.20×10-9 and rs11992444, P=1.30×10-12) at 2 distinct genomic loci, that were mapped to 3 potentially causal genes: T-box transcription factor 15 (TBX15), tryptophanyl tRNA synthetase 2, mitochondrial (WARS2) and early B-cell factor-2 (EBF2) through functional annotation. Bayesian colocalization additionally suggested a role of RP4-712E4.1. Genetically predicted differences in adiposity-adjusted PAT were causally associated with adverse left ventricular remodeling. CONCLUSIONS This study provides insights into the genetic architecture determining differential PAT deposition, identifies causal links with left structural and functional parameters, and provides novel data about the pathophysiological importance of adiposity distribution.

Single Nucleotide Polymorphism rs9277336 Controls the Nuclear Alpha Actinin 4-Human Leukocyte Antigen-DPA1 Axis and Pulmonary Endothelial Pathophenotypes in Pulmonary Arterial Hypertension

Background Pulmonary arterial hypertension (PAH) is a complex, fatal disease where disease severity has been associated with the single nucleotide polymorphism (SNP) rs2856830, located near the human leukocyte antigen DPA1 (HLA-DPA1) gene. We aimed to define the genetic architecture of functional variants associated with PAH disease severity by identifying allele-specific binding transcription factors and downstream targets that control endothelial pathophenotypes and PAH. Methods and Results Electrophoretic mobility shift assays of oligonucleotides containing SNP rs2856830 and 8 SNPs in linkage disequilibrium revealed functional SNPs via allele-imbalanced binding to human pulmonary arterial endothelial cell nuclear proteins. DNA pulldown proteomics identified SNP-binding proteins. SNP genotyping and clinical correlation analysis were performed in 84 patients with PAH at University of Pittsburgh Medical Center and in 679 patients with PAH in the All of Us database. SNP rs9277336 was identified as a functional SNP in linkage disequilibrium (r2>0.8) defined by rs2856830, and the minor allele was associated with decreased hospitalizations and improved cardiac output in patients with PAH, an index of disease severity. SNP pulldown proteomics showed allele-specific binding of nuclear ACTN4 (alpha actinin 4) protein to rs9277336 minor allele. Both ACTN4 and HLA-DPA1 were downregulated in pulmonary endothelium in human patients and rodent models of PAH. Via transcriptomic and phenotypic analyses, knockdown of HLA-DPA1 phenocopied knockdown of ACTN4, both similarly controlling cell structure pathways, immune pathways, and endothelial dysfunction. Conclusions We defined the pathogenic activity of functional SNP rs9277336, entailing the allele-specific binding of ACTN4 and controlling expression of the neighboring HLA-DPA1 gene. Through inflammatory or genetic means, downregulation of this ACTN4-HLA-DPA1 regulatory axis promotes endothelial pathophenotypes, providing a mechanistic explanation for the association between this SNP and PAH outcomes.

Genome‐wide association study of tuberculosis in the western Chinese Han and Tibetan population

Tuberculosis (TB) remains a serious global public health threat. Accumulated evidence has demonstrated that human susceptibility to TB has a strong genetic basis. And different susceptibility single nucleotide polymorphisms (SNP) have been reported in different studies. To gain greater insight into the host susceptibility to TB, we perform a two‐stage genome‐wide association study to identify the susceptible loci of TB. In the discovery stage, 3116 (1532 TB patients and 1584 healthy controls) and 439 (211 TB patients and 228 healthy controls) individuals were genome‐wide genotyped from a western Chinese Han and Tibetan population, respectively. Based on the additive genetic model, we discovered 14 and three independent loci that had potential associations with TB susceptibility in the Chinese Han and Tibetan populations, respectively (p < 1 × 10⁻⁵). Furthermore, we conducted an imputation‐based meta‐analysis on another two East Asia cohorts to replicate our findings. We identified one independent locus harbored by the human leukocyte antigen (HLA) class II genes that was genome‐wide significantly associated with TB (lead SNP rs111875628 with a p‐value of 2.20 × 10⁻⁹). Our findings suggest a novel mechanism of the interaction with the HLA class II genes and reinforce the importance of the HLA class II alleles in response to TB.

Genetic Variants Associated With Systolic Blood Pressure in Children and Adolescents

Background Genetics, along with lifestyle and behavioral characteristics, play an important role in hypertension in adults. Our aim was to identify genetic variants associated with blood pressure in childhood and adolescence. Methods and Results We conducted a candidate single-nucleotide polymorphism (SNP) analysis and genome-wide association study among 9778 participants aged <18 years in BioVU, the Vanderbilt University Medical Center biobank. The outcome was childhood blood pressure percentile from age 0 to 18 years. For the candidate SNP analysis, a total of 457 previously identified SNPs were examined. Linear regression was used to test the association between genetic variants and median systolic blood pressure (SBP) percentile. Adjusted models included median age, self-reported sex, race, the first 4 principal components of ancestry, and median body mass index Z score. Analyses were conducted in the overall cohort and stratified by age group. A polygenic risk score was calculated for each participant, and the association between polygenic risk score and median SBP percentile in childhood was examined using linear regression. In the overall candidate SNP analysis, 2 SNPs reached significance: rs1018148 (FBN1; P=1.0×10-4) and rs11105354 (ATP2B1; P=1.4×10-4). In the postpuberty age group, 1 SNP reached significance: rs1018148 (FBN1; P=2.2×10-5). In the genome-wide association study of all participants, no SNPs reached genome-wide significance. Higher polygenic risk score was associated with higher SBP percentile (β, 0.35 [95% CI, 0.10-0.60)], and there was a significant interaction with age (P for interaction<0.01). Conclusions These findings suggest that genetic variants play an important role in SBP in childhood and adolescence and provide evidence for age-specific genetic associations with SBP.

Circulating Soluble CD163, Associations With Cardiovascular Outcomes and Mortality, and Identification of Genetic Variants in Older Individuals: The Cardiovascular Health Study

Background Monocytes/macrophages participate in cardiovascular disease. CD163 (cluster of differentiation 163) is a monocyte/macrophage receptor, and the shed sCD163 (soluble CD163) reflects monocyte/macrophage activation. We examined the association of sCD163 with incident cardiovascular disease events and performed a genome-wide association study to identify sCD163-associated variants. Methods and Results We measured plasma sCD163 in 5214 adults (aged ≥65 years, 58.7% women, 16.2% Black) of the CHS (Cardiovascular Health Study). We used Cox regression models (associations of sCD163 with incident events and mortality); median follow-up was 26 years. Genome-wide association study analyses were stratified on race. Adjusted for age, sex, and race and ethnicity, sCD163 levels were associated with all-cause mortality (hazard ratio [HR], 1.08 [95% CI, 1.04-1.12] per SD increase), cardiovascular disease mortality (HR, 1.15 [95% CI, 1.09-1.21]), incident coronary heart disease (HR, 1.10 [95% CI, 1.04-1.16]), and incident heart failure (HR, 1.18 [95% CI, 1.12-1.25]). When further adjusted (eg, cardiovascular disease risk factors), only incident coronary heart disease lost significance. In European American individuals, genome-wide association studies identified 38 variants on chromosome 2 near MGAT5 (top result rs62165726, P=3.3×10^-18),19 variants near chromosome 17 gene ASGR1 (rs55714927, P=1.5×10^-14), and 18 variants near chromosome 11 gene ST3GAL4. These regions replicated in the European ancestry ADDITION-PRO cohort, a longitudinal cohort study nested in the Danish arm of the Anglo-Danish-Dutch study of Intensive Treatment Intensive Treatment In peOple with screeNdetcted Diabetes in Primary Care. In Black individuals, we identified 9 variants on chromosome 6 (rs3129781 P=7.1×10^-9) in the HLA region, and 3 variants (rs115391969 P=4.3×10^-8) near the chromosome 16 gene MYLK3_. Conclusions Monocyte function, as measured by sCD163, may be predictive of overall and cardiovascular-specific mortality and incident heart failure.

Genetic Architecture of Plasma Alpha-Aminoadipic Acid Reveals a Relationship With High-Density Lipoprotein Cholesterol

Background Elevated plasma levels of alpha-aminoadipic acid (2-AAA) have been associated with the development of type 2 diabetes and atherosclerosis. However, the nature of the association remains unknown. Methods and Results We identified genetic determinants of plasma 2-AAA through meta-analysis of genome-wide association study data in 5456 individuals of European, African, and Asian ancestry from the Framingham Heart Study, Diabetes Prevention Program, Jackson Heart Study, and Shanghai Women's and Men's Health Studies. No single nucleotide polymorphisms reached genome-wide significance across all samples. However, the top associations from the meta-analysis included single-nucleotide polymorphisms in the known 2-AAA pathway gene DHTKD1, and single-nucleotide polymorphisms in genes involved in mitochondrial respiration (NDUFS4) and macrophage function (MSR1). We used a Mendelian randomization instrumental variable approach to evaluate relationships between 2-AAA and cardiometabolic phenotypes in large disease genome-wide association studies. Mendelian randomization identified a suggestive inverse association between increased 2-AAA and lower high-density lipoprotein cholesterol (P=0.005). We further characterized the genetically predicted relationship through measurement of plasma 2-AAA and high-density lipoprotein cholesterol in 2 separate samples of individuals with and without cardiometabolic disease (N=98), and confirmed a significant negative correlation between 2-AAA and high-density lipoprotein (rs=-0.53, P<0.0001). Conclusions 2-AAA levels in plasma may be regulated, in part, by common variants in genes involved in mitochondrial and macrophage function. Elevated plasma 2-AAA associates with reduced levels of high-density lipoprotein cholesterol. Further mechanistic studies are required to probe this as a possible mechanism linking 2-AAA to future cardiometabolic risk.

Weak Association Between Genetic Markers of Hyperuricemia and Cardiorenal Outcomes: Insights From the STANISLAS Study Cohort With a 20-Year Follow-Up

Background

Hyperuricemia is associated with poor cardiovascular outcomes, although it is uncertain whether this relationship is causal in nature. This study aimed to: (1) assess the heritability of serum uric acid (SUA) levels, (2) conduct a genome‐wide association study on SUA levels, and (3) investigate the association between certain single‐nucleotide polymorphisms and target organ damage.

Methods and Results

The STANISLAS (Suivi Temporaire Annuel Non‐Invasif de la Santé des Lorrains Assurés Sociaux) study cohort is a single‐center longitudinal cohort recruited between 1993 and 1995 (visit 1), with a last visit (visit 4 [V4]) performed ≈20 years apart. Serum lipid profile, SUA, urinary albumin/creatinine ratio, estimated glomerular filtration rate, 24‐hour ambulatory blood pressure monitoring, transthoracic echocardiography, pulse wave velocity, and genotyping for each participant were assessed at V4. A total of 1573 participants were included at V4, among whom 1417 had available SUA data at visit 1. Genome‐wide association study results highlighted multiple single‐nucleotide polymorphisms on the SLC2A9 gene linked to SUA levels. Carriers of the most associated mutated SLC2A9 allele (rs16890979) had significantly lower SUA levels. Although SUA level at V4 was highly associated with diabetes, prediabetes, higher body mass index, CRP (C‐reactive protein) levels, estimated glomerular filtration rate variation (visit 1–V4), carotid intima‐media thickness, and pulse wave velocity, rs16890979 was only associated with higher carotid intima‐media thickness.

Conclusions

Our findings demonstrate that rs16890979, a genetic determinant of SUA levels located on the SLC2A9 gene, is associated with carotid intima‐media thickness despite significant associations between SUA levels and several clinical outcomes, thereby lending support to the hypothesis of a link between SUA and cardiovascular disease.

Genetic susceptibility to hepatocellular carcinoma in chromosome 22q13.31, findings of a genome-wide association study

BACKGROUND AND AIM

Chronic hepatitis C virus (HCV) infection, long-term alcohol use, cigarette smoking, and obesity are the major risk factors for hepatocellular carcinoma (HCC) in the United States, but the disease risk varies substantially among individuals with these factors, suggesting host susceptibility to and gene-environment interactions in HCC. To address genetic susceptibility to HCC, we conducted a genome-wide association study (GWAS).

METHODS

Two case-control studies on HCC were conducted in the United States. DNA samples were genotyped using the Illumian microarray chip with over 710 000 single nucleotide polymorphisms (SNPs). We compared these SNPs between 705 HCC cases and 1455 population controls for their associations with HCC and verified our findings in additional studies.

RESULTS

In this GWAS, we found that two SNPs were associated with HCC at P < 5E-8 and six SNPs at P < 5E-6 after adjusting for age, sex, and the top three principal components (PCs). Five of the SNPs in chromosome 22q13.31, three in PNPLA3 (rs2281135, rs2896019, and rs4823173) and two in SAMM50 (rs3761472, rs3827385), were replicated in a small US case-control study and a cohort study in Singapore. The associations remained significant after adjusting for body mass index and HCV infection. Meta-analysis of multiple datasets indicated that these SNPs were significantly associated with HCC.

CONCLUSIONS

SNPs in PNPLA3 and SAMM50 are known risk loci for nonalcoholic fatty liver disease (NAFLD) and are suspected to be associated with HCC. Our GWAS demonstrated the associations of these SNPs with HCC in a US population. Biological mechanisms underlying the relationship remain to be elucidated.

Genetic Profile of Endotoxemia Reveals an Association With Thromboembolism and Stroke

Background Translocation of lipopolysaccharide from gram-negative bacteria into the systemic circulation results in endotoxemia. In addition to acute infections, endotoxemia is detected in cardiometabolic disorders, such as cardiovascular diseases and obesity. Methods and Results We performed a genome-wide association study of serum lipopolysaccharide activity in 11 296 individuals from 6 different Finnish study cohorts. Endotoxemia was measured by limulus amebocyte lysate assay in the whole population and by 2 other techniques (Endolisa and high-performance liquid chromatography/tandem mass spectrometry) in subpopulations. The associations of the composed genetic risk score of endotoxemia and thrombosis-related clinical end points for 195 170 participants were analyzed in FinnGen. Lipopolysaccharide activity had a genome-wide significant association with 741 single-nucleotide polymorphisms in 5 independent loci, which were mainly located at genes affecting the contact activation of the coagulation cascade and lipoprotein metabolism and explained 1.5% to 9.2% of the variability in lipopolysaccharide activity levels. The closest genes included KNG1, KLKB1, F12, SLC34A1, YPEL4, CLP1, ZDHHC5, SERPING1, CBX5, and LIPC. The genetic risk score of endotoxemia was associated with deep vein thrombosis, pulmonary embolism, pulmonary heart disease, and venous thromboembolism. Conclusions The biological activity of lipopolysaccharide in the circulation (ie, endotoxemia) has a small but highly significant genetic component. Endotoxemia is associated with genetic variation in the contact activation pathway, vasoactivity, and lipoprotein metabolism, which play important roles in host defense, lipopolysaccharide neutralization, and thrombosis, and thereby thromboembolism and stroke.

Heterogeneity of a dwarf phenotype in Dutch traditional chicken breeds revealed by genomic analyses

The growth of animals is a complex trait, in chicken resulting in a diverse variety of forms, caused by a heterogeneous genetic basis. Bantam chicken, known as an exquisite form of dwarfism, has been used for crossbreeding to create corresponding dwarf counterparts for native fowls in the Dutch populations. Here, we demonstrate the heterogeneity of the bantam trait in Dutch chickens and reveal the underlying genetic causes, using whole-genome sequence data from matching pairs of bantam and normal-sized breeds. During the bantam-oriented crossbreeding, various bantam origins were used to introduce the bantam phenotype, and three major bantam sources were identified and clustered. The genome-wide association studies revealed multiple genetic variants and genes associated with bantam phenotype, including HMGA2 and PRDM16, genes involved in body growth and stature. The comparison of associated variants among studies illustrated differences related to divergent bantam origins, suggesting a clear heterogeneity among bantam breeds. We show that in neo-bantam breeds, the bantam-related regions underwent a strong haplotype introgression from the bantam source, outcompeting haplotypes from the normal-sized counterpart. The bantam heterogeneity is further confirmed by the presence of multiple haplotypes comprising associated alleles, which suggests the selection of the bantam phenotype is likely subject to a convergent direction across populations. Our study demonstrates that the diverse history of human-mediated crossbreeding has contributed to the complexity and heterogeneity of the bantam phenotype.

Steelhead (Oncorhynchus mykiss) lineages and sexes show variable patterns of association of adult migration timing and age-at-maturity traits with two genomic regions

As life history diversity plays a critical role in supporting the resilience of exploited populations, understanding the genetic basis of those life history variations is important for conservation management. However, effective application requires a robust understanding of the strength and universality of genetic associations. Here, we examine genetic variation of single nucleotide polymorphisms in genomic regions previously associated with migration phenology and age-at-maturity in steelhead (Oncorhynchus mykiss) from the Columbia River. We found chromosome 28 markers (GREB1L, ROCK1 genes) explained significant variance in migration timing in both coastal and inland steelhead. However, strength of association was much greater in coastal than inland steelhead (R 2 0.51 vs. 0.08), suggesting that genomic background and challenging inland migration pathways may act to moderate effects of this region. Further, we found that chromosome 25 candidate markers (SIX6 gene) were significantly associated with age and size at first return migration for inland steelhead, and this pattern was mediated by sex in a predictable pattern (males R 2 = 0.139-0.170; females R 2 = 0.096-0.111). While this encourages using these candidate regions in predicting life history characteristics, we suggest that stock-specific associations and haplotype frequencies will be useful in guiding implementation of genetic assays to inform management.

Natural variation of root lesion nematode antagonism in the biocontrol fungus Clonostachys rosea and identification of biocontrol factors through genome-wide association mapping

Biological control is a promising approach to reduce plant diseases caused by nematodes to ensure high productivity in agricultural production. Large-scale analyses of genetic variation in fungal species used for biocontrol can generate knowledge regarding interaction mechanisms that can improve efficacy of biocontrol applications. In this study, we performed a genome-wide association study (GWAS) for in vitro antagonism against the root lesion nematode Pratylenchus penetrans in 53 previously genome re-sequenced strains of the biocontrol fungus Clonostachys rosea. Nematode mortality in C. rosea potato dextrose broth (PDB) culture filtrates was highly variable and showed continuous variation (p < .001) between strains, indicating a polygenic inheritance. Twenty-one strains produced culture filtrates with higher (p ≤ .05) nematode mortality compared with the PDB control treatment, while ten strains lowered (p ≤ .05) the mortality. The difference in in vitro antagonism against P. penetrans correlated with antagonism against the soybean cyst nematode Heterodera glycines, indicating lack of host specificity in C. rosea. An empirical Bayesian multiple hypothesis testing approach identified 279 single nucleotide polymorphism markers significantly (local false sign rate < 10-10) associated with the trait. Genes present in the genomic regions associated with nematicidal activity included several membrane transporters, a chitinase and genes encoding proteins predicted to biosynthesize secondary metabolites. Gene deletion strains of the predicted nonribosomal peptide synthetase genes nps4 and nps5 were generated and showed increased (p ≤ .001) fungal growth and conidiation rates compared to the wild type. Deletion strains also exhibited reduced (p < .001) nematicidal activity and reduced (p ≤ .05) biocontrol efficacy against nematode root disease and against fusarium foot rot on wheat. In summary, we show that the GWAS approach can be used to identify biocontrol factors in C. rosea, specifically the putative nonribosomal peptide synthetases NPS4 and NPS5.

Developing Novel Therapeutic Strategies Against Ring Finger Protein 213 Vasculopathy: Is It an Achievable Goal?

See Article Choi et al.

Genetic Determinants of Circulating Glycine Levels and Risk of Coronary Artery Disease

Background Recent studies have revealed sexually dimorphic associations between the carbamoyl-phosphate synthase 1 locus, intermediates of the metabolic pathway leading from choline to urea, and risk of coronary artery disease ( CAD ) in women. Based on evidence from the literature, the atheroprotective association with carbamoyl-phosphate synthase 1 could be mediated by the strong genetic effect of this locus on increased circulating glycine levels. Methods and Results We sought to identify additional genetic determinants of circulating glycine levels by carrying out a meta-analysis of genome-wide association study data in up to 30 118 subjects of European ancestry. Mendelian randomization and other analytical approaches were used to determine whether glycine-associated variants were associated with CAD and traditional risk factors. Twelve loci were significantly associated with circulating glycine levels, 7 of which were not previously known to be involved in glycine metabolism ( ACADM , PHGDH , COX 18- ADAMTS 3, PSPH , TRIB 1, PTPRD , and ABO ). Glycine-raising alleles at several loci individually exhibited directionally consistent associations with decreased risk of CAD . However, these effects could not be attributed directly to glycine because of associations with other CAD -related traits. By comparison, genetic models that only included the 2 variants directly involved in glycine degradation and for which there were no other pleiotropic associations were not associated with risk of CAD or blood pressure, lipid levels, and obesity-related traits. Conclusions These results provide additional insight into the genetic architecture of glycine metabolism, but do not yield conclusive evidence for a causal relationship between circulating levels of this amino acid and risk of CAD in humans.

Variants at the APOE /C1/C2/C4 Locus Modulate Cholesterol Efflux Capacity Independently of High-Density Lipoprotein Cholesterol

Background Macrophage cholesterol efflux to high-density lipoproteins ( HDLs ) is the first step of reverse cholesterol transport. The cholesterol efflux capacity ( CEC ) of HDL particles is a protective risk factor for coronary artery disease independent of HDL cholesterol levels. Using a genome-wide association study approach, we aimed to identify pathways that regulate CEC in humans. Methods and Results We measured CEC in 5293 French Canadians. We tested the genetic association between 4 CEC measures and genotypes at >9 million common autosomal DNA sequence variants. These analyses yielded 10 genome-wide significant signals ( P<6.25×10-9) representing 7 loci. Five of these loci harbor genes with important roles in lipid biology ( CETP , LIPC , LPL , APOA 1/C3/A4/A5, and APOE /C1/C2/C4). Except for the APOE /C1/C2/C4 variant ( rs141622900, P nonadjusted=1.0×10-11; P adjusted=8.8×10-9), the association signals disappear when correcting for HDL cholesterol and triglyceride levels. The additional 2 significant signals were near the PPP 1 CB / PLB 1 and RBFOX 3/ ENPP 7 genes. In secondary analyses, we considered candidate functional variants for 58 genes implicated in HDL biology, as well as 239 variants associated with blood lipid levels and/or coronary artery disease risk by genome-wide association study . These analyses identified 27 significant CEC associations, implicating 5 additional loci ( GCKR , LIPG , PLTP , PPARA , and TRIB 1). Conclusions Our genome-wide association study identified common genetic variation at the APOE /C1/C2/C4 locus as a major determinant of CEC that acts largely independently of HDL cholesterol. We predict that HDL -based therapies aiming at increasing CEC will be modulated by changes in the expression of apolipoproteins in this gene cluster.

Exploring the genetics of lesion and nodal resistance in pea (Pisum sativum L.) to Sclerotinia sclerotiorum using genome-wide association studies and RNA-Seq

The disease white mold caused by the fungus Sclerotinia sclerotiorum is a significant threat to pea production, and improved resistance to this disease is needed. Nodal resistance in plants is a phenomenon where a fungal infection is prevented from passing through a node, and the infection is limited to an internode region. Nodal resistance has been observed in some pathosystems such as the pea (Pisum sativum L.)-S. sclerotiorum pathosystem. In addition to nodal resistance, different pea lines display different levels of stem lesion size restriction, referred to as lesion resistance. It is unclear whether the genetics of lesion resistance and nodal resistance are identical or different. This study applied genome-wide association studies (GWAS) and RNA-Seq to understand the genetic makeup of these two types of resistance. The time series RNA-Seq experiment consisted of two pea lines (the susceptible 'Lifter' and the partially resistant PI 240515), two treatments (mock inoculated samples and S. sclerotiorum-inoculated samples), and three time points (12, 24, and 48 hr post inoculation). Integrated results from GWAS and RNA-Seq analyses identified different redox-related transcripts for lesion and nodal resistances. A transcript encoding a glutathione S-transferase was the only shared resistance variant for both phenotypes. There were more leucine rich-repeat containing transcripts found for lesion resistance, while different candidate resistance transcripts such as a VQ motif-containing protein and a myo-inositol oxygenase were found for nodal resistance. This study demonstrated the robustness of combining GWAS and RNA-Seq for identifying white mold resistance in pea, and results suggest different genetics underlying lesion and nodal resistance.

Genome Wide Association Study Identifies the HMGCS2 Locus to be Associated With Chlorthalidone Induced Glucose Increase in Hypertensive Patients

Background

Thiazide and thiazide‐like diuretics are first‐line medications for treating uncomplicated hypertension. However, their use has been associated with adverse metabolic events, including hyperglycemia and incident diabetes mellitus, with incompletely understood mechanisms. Our goal was to identify genomic variants associated with thiazide‐like diuretic/chlorthalidone‐induced glucose change.

Methods and Results

Genome‐wide analysis of glucose change after treatment with chlorthalidone was performed by race among the white (n=175) and black (n=135) participants from the PEAR‐2 (Pharmacogenomic Evaluation of Antihypertensive Responses‐2). Single‐nucleotide polymorphisms with P<5×10⁻⁸ were further prioritized using in silico analysis based on their expression quantitative trait loci function. Among blacks, an intronic single‐nucleotide polymorphism (rs9943291) in the HMGCS2 was associated with increase in glucose levels following chlorthalidone treatment (ß=12.5; P=4.17×10⁻⁸). G‐allele carriers of HMGCS2 had higher glucose levels (glucose change=+16.29 mg/dL) post chlorthalidone treatment compared with noncarriers of G allele (glucose change=+2.80 mg/dL). This association was successfully replicated in an independent replication cohort of hydrochlorothiazide‐treated participants from the PEAR study (ß=5.54; P=0.023). A meta‐analysis of the 2 studies was performed by race in Meta‐Analysis Helper, where this single‐nucleotide polymorphism, rs9943291, was genome‐wide significant with a meta‐analysis P value of 3.71×10⁻⁸. HMGCS2, a part of the HMG‐CoA synthase family, is important for ketogenesis and cholesterol synthesis pathways that are essential in glucose homeostasis.

Conclusions

These results suggest that HMGCS2 is a promising candidate gene involved in chlorthalidone and Hydrochlorothiazide (HCTZ)‐induced glucose change. This may provide insights into the mechanisms involved in thiazide‐induced hyperglycemia that may ultimately facilitate personalized approaches to antihypertensive selection for hypertension treatment.

Clinical Trial Registration

URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00246519 and NCT01203852.

Causal Effect of Plasminogen Activator Inhibitor Type 1 on Coronary Heart Disease

Background

Plasminogen activator inhibitor type 1 (PAI‐1) plays an essential role in the fibrinolysis system and thrombosis. Population studies have reported that blood PAI‐1 levels are associated with increased risk of coronary heart disease (CHD). However, it is unclear whether the association reflects a causal influence of PAI‐1 on CHD risk.

Methods and Results

To evaluate the association between PAI‐1 and CHD, we applied a 3‐step strategy. First, we investigated the observational association between PAI‐1 and CHD incidence using a systematic review based on a literature search for PAI‐1 and CHD studies. Second, we explored the causal association between PAI‐1 and CHD using a Mendelian randomization approach using summary statistics from large genome‐wide association studies. Finally, we explored the causal effect of PAI‐1 on cardiovascular risk factors including metabolic and subclinical atherosclerosis measures. In the systematic meta‐analysis, the highest quantile of blood PAI‐1 level was associated with higher CHD risk comparing with the lowest quantile (odds ratio=2.17; 95% CI: 1.53, 3.07) in an age‐ and sex‐adjusted model. The effect size was reduced in studies using a multivariable‐adjusted model (odds ratio=1.46; 95% CI: 1.13, 1.88). The Mendelian randomization analyses suggested a causal effect of increased PAI‐1 level on CHD risk (odds ratio=1.22 per unit increase of log‐transformed PAI‐1; 95% CI: 1.01, 1.47). In addition, we also detected a causal effect of PAI‐1 on elevating blood glucose and high‐density lipoprotein cholesterol.

Conclusions

Our study indicates a causal effect of elevated PAI‐1 level on CHD risk, which may be mediated by glucose dysfunction.

Rare variants in optic disc area gene CARD10 enriched in primary open-angle glaucoma

BACKGROUND

Genome-wide association studies (GWAS) have identified association of common alleles with primary open-angle glaucoma (POAG) and its quantitative endophenotypes near numerous genes. This study aims to determine whether rare pathogenic variants in these disease-associated genes contribute to POAG.

METHODS

Participants fulfilled strict inclusion criteria of advanced POAG at a young age of diagnosis. Myocilin mutation carriers were excluded using direct sequencing. Whole exome sequencing was performed on 187 glaucoma cases and 103 local screened nonglaucoma controls then joint-called with exomes of 993 previously sequenced Australian controls. GWAS-associated genes were assessed for enrichment of rare predicted pathogenic variants in POAG. Significantly enriched genes were compared against Exome Aggregation Consortium (ExAC) public control.

RESULTS

Eighty-six GWAS disease or trait-associated glaucoma genes were captured and sequenced. CARD10 showed enrichment after Bonferroni correction for rare variants in glaucoma cases (OR = 13.2, P = 6.94 × 10^-5) with mutations identified in 4.28% of our POAG cohort compared to 0.27% in controls. CARD10 was significantly associated with optic disc parameters in previous GWAS. The whole GWAS gene set showed no enrichment in POAG overall (OR = 1.12, P = 0.51).

CONCLUSION

We report here an enrichment of rare predicted pathogenic coding variants within a GWAS-associated locus in POAG (CARD10). These findings indicate that both common and rare pathogenic coding variants in CARD10 may contribute to POAG pathogenesis.

Genome-wide time-to-event analysis on smoking progression stages in a family-based study

BACKGROUND

Various pivotal stages in smoking behavior can be identified, including initiation, conversion from experimenting to established use, development of tolerance, and cessation. Previous studies have shown high heritability for age of smoking initiation and cessation; however, time-to-event genome-wide association studies aiming to identify underpinning genes that accelerate or delay these transitions are missing to date.

METHODS

We investigated which single nucleotide polymorphisms (SNPs) across the whole genome contribute to the hazard ratio of transition between different stages of smoking behavior by performing time-to-event analyses within a large Finnish twin family cohort (N = 1962), and further conducted mediation analyses of plausible intermediate traits for significant SNPs.

RESULTS

Genome-wide significant signals were detected for three of the four transitions: (1) for smoking cessation on 10p14 (P = 4.47e-08 for rs72779075 flanked by RP11-575N15 and GATA3), (2) for tolerance on 11p13 (P = 1.29e-08 for rs11031684 in RP1-65P5.1), mediated by smoking quantity, and on 9q34.12 (P = 3.81e-08 for rs2304808 in FUBP3), independent of smoking quantity, and (3) for smoking initiation on 19q13.33 (P = 3.37e-08 for rs73050610 flanked by TRPM4 and SLC6A16) in analysis adjusted for first time sensations. Although our top SNPs did not replicate, another SNP in the TRPM4-SLC6A16 gene region showed statistically significant association after region-based multiple testing correction in an independent Australian twin family sample.

CONCLUSION

Our results suggest that the functional effect of the TRPM4-SLC6A16 gene region deserves further investigation, and that complex neurotransmitter networks including dopamine and glutamate may play a critical role in smoking initiation. Moreover, comparison of these results implies that genetic contributions to the complex smoking behavioral phenotypes vary among the transitions.

Identification of PTCSC3 as a Novel Locus for Large-Vessel Ischemic Stroke: A Genome-Wide Association Study

Background

Ischemic stroke is a major cause of death and disability in the world. A major ischemic stroke subtype, large‐vessel ischemic stroke (large artery atherosclerosis; LAA), has been shown to have some genetic components in individuals of European ancestry. However, it is not clear whether the genetic predisposition to LAA stroke varies among ethnicities. We sought to identify genetic factors that contribute to LAA stroke in 2 independent samples of Han Chinese individuals.

Methods and Results

Novel genetic variants that predispose individuals to LAA stroke were identified using a genome‐wide association study (GWAS) of 444 individuals with LAA stroke and 1727 controls in a Han Chinese population residing in Taiwan. The study was replicated in an independent Han Chinese population comprising an additional 319 cases and 1802 controls. We identified 5 single‐nucleotide polymorphisms, including rs2415317 (P=3.10×10⁻⁸), rs934075 (P=4.00×10⁻⁹), rs944289 (P=3.57×10⁻⁸), rs2787417 (P=1.76×10⁻⁸), and rs1952706 (P=2.92×10⁻⁸), at one novel locus on chromosome 14q13.3 within PTCSC3 (encoding papillary thyroid carcinoma susceptibility candidate 3) that were associated with LAA stroke at genome‐wide significance (P<5×10⁻⁸).

Conclusions

Our data provide strong support for future studies on the role of PTCSC3 in the pathogenesis of LAA stroke and the association between LAA stroke development and thyroid function. In addition, these findings provide insights into the genetic basis of LAA stroke and identify a novel pathway that might be applicable for future therapeutic intervention.

Genome-wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.

Genome-wide association study of pre-eclampsia detects novel maternal single nucleotide polymorphisms and copy-number variants in subsets of the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) study cohort

A genome-wide association study was undertaken to identify maternal single nucleotide polymorphisms (SNPs) and copy-number variants (CNVs) associated with pre-eclampsia. Case-control analysis was performed on 1070 Afro-Caribbean (n = 21 cases and 1049 controls) and 723 Hispanic (n = 62 cases and 661 controls) mothers and 1257 mothers of European ancestry (n = 50 cases and 1207 controls) from the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) study. European ancestry subjects were genotyped on Illumina Human610-Quad and Afro-Caribbean and Hispanic subjects were genotyped on Illumina Human1M-Duo BeadChip microarrays. Genome-wide SNP data were analyzed using PLINK. CNVs were called using three detection algorithms (GNOSIS, PennCNV, and QuantiSNP), merged using CNVision, and then screened using stringent criteria. SNP and CNV findings were compared to those of the Study of Pregnancy Hypertension in Iowa (SOPHIA), an independent pre-eclampsia case-control dataset of Caucasian mothers (n = 177 cases and 116 controls). A list of top SNPs were identified for each of the HAPO ethnic groups, but none reached Bonferroni-corrected significance. Novel candidate CNVs showing enrichment among pre-eclampsia cases were also identified in each of the three ethnic groups. Several variants were suggestively replicated in SOPHIA. The discovered SNPs and copy-number variable regions present interesting candidate genetic variants for pre-eclampsia that warrant further replication and investigation.