Strongest Association Signal Research Articles

Enhancing portability of trans-ancestral polygenic risk scores through tissue-specific functional genomic data integration.

Portability of trans-ancestral polygenic risk scores is often confounded by differences in linkage disequilibrium and genetic architecture between ancestries. Recent literature has shown that prioritizing GWAS SNPs with functional genomic evidence over strong association signals can improve model portability. We leveraged three RegulomeDB-derived functional regulatory annotations-SURF, TURF, and TLand-to construct polygenic risk models across a set of quantitative and binary traits highlighting functional mutations tagged by trait-associated tissue annotations. Tissue-specific prioritization by TURF and TLand provide a significant improvement in model accuracy over standard polygenic risk score (PRS) models across all traits. We developed the Trans-ancestral Iterative Tissue Refinement (TITR) algorithm to construct PRS models that prioritize functional mutations across multiple trait-implicated tissues. TITR-constructed PRS models show increased predictive accuracy over single tissue prioritization. This indicates our TITR approach captures a more comprehensive view of regulatory systems across implicated tissues that contribute to variance in trait expression.

A shared genetic contribution to osteoarthritis and COVID-19 outcomes: a large-scale genome-wide cross-trait analysis.

Patients with osteoarthritis (OA) are exposed to an increased risk of adverse outcomes of COVID-19, and they tend to experience disruption in access to healthcare services and exercise facilities. However, a deep understanding of this comorbidity phenomenon and the underlying genetic architecture of the two diseases is still unclear. In this study, we aimed to untangle the relationship between OA and COVID-19 outcomes by conducting a large-scale genome-wide cross-trait analysis. Genetic correlation and causal relationships between OA and COVID-19 outcomes (critical COVID-19, COVID-19 hospitalization, and COVID-19 infection) were estimated by linkage disequilibrium score regression and Mendelian Randomization approaches. We further applied Multi-Trait Analysis of GWAS and colocalization analysis to identify putative functional genes associated with both OA and COVID-19 outcomes. Significant positive genetic correlations between OA susceptibility and both critical COVID-19 (rg=0.266, P=0.0097) and COVID-19 hospitalization (rg=0.361, P=0.0006) were detected. However, there was no evidence to support causal genetic relationships between OA and critical COVID-19 (OR=1.17[1.00-1.36], P=0.049) or OA and COVID-19 hospitalization OR=1.08[0.97-1.20], P=0.143). These results were robustly consistent after the removal of obesity-related single nucleotide polymorphisms (SNPs). Moreover, we identified a strong association signal located near the FYCO1 gene (lead SNPs: rs71325101 for critical COVID-19, Pmeta=1.02×10-34; rs13079478 for COVID-19 hospitalization, Pmeta=1.09×10-25). Our findings further confirmed the comorbidity of OA and COVID-19 severity, but indicate a non-causal impact of OA on COVID-19 outcomes. The study offers an instructive perspective that OA patients did not generate negative COVID-19 outcomes during the pandemic in a causal way. Further clinical guidance can be formulated to enhance the quality of self-management in vulnerable OA patients.

Analysis of genetically determined gene expression suggests role of inflammatory processes in exfoliation syndrome

BackgroundExfoliation syndrome (XFS) is an age-related systemic disorder characterized by excessive production and progressive accumulation of abnormal extracellular material, with pathognomonic ocular manifestations. It is the most common cause of secondary glaucoma, resulting in widespread global blindness. The largest global meta-analysis of XFS in 123,457 multi-ethnic individuals from 24 countries identified seven loci with the strongest association signal in chr15q22–25 region near LOXL1. Expression analysis have so far correlated coding and a few non-coding variants in the region with LOXL1 expression levels, but functional effects of these variants is unclear. We hypothesize that analysis of the contribution of the genetically determined component of gene expression to XFS risk can provide a powerful method to elucidate potential roles of additional genes and clarify biology that underlie XFS.ResultsTranscriptomic Wide Association Studies (TWAS) using PrediXcan models trained in 48 GTEx tissues leveraging on results from the multi-ethnic and European ancestry GWAS were performed. To eliminate the possibility of false-positive results due to Linkage Disequilibrium (LD) contamination, we i) performed PrediXcan analysis in reduced models removing variants in LD with LOXL1 missense variants associated with XFS, and variants in LOXL1 models in both multiethnic and European ancestry individuals, ii) conducted conditional analysis of the significant signals in European ancestry individuals, and iii) filtered signals based on correlated gene expression, LD and shared eQTLs, iv) conducted expression validation analysis in human iris tissues. We observed twenty-eight genes in chr15q22–25 region that showed statistically significant associations, which were whittled down to ten genes after statistical validations. In experimental analysis, mRNA transcript levels for ARID3B, CD276, LOXL1, NEO1, SCAMP2, and UBL7 were significantly decreased in iris tissues from XFS patients compared to control samples. TWAS genes for XFS were significantly enriched for genes associated with inflammatory conditions. We also observed a higher incidence of XFS comorbidity with inflammatory and connective tissue diseases.ConclusionOur results implicate a role for connective tissues and inflammation pathways in the etiology of XFS. Targeting the inflammatory pathway may be a potential therapeutic option to reduce progression in XFS.

PRDX6 AS1 gene polymorphisms and SLE susceptibility in Chinese populations.

BackgroundSystemic lupus erythematosus (SLE) is a complex, multisystem autoimmune disease that is characterized by the production of autoantibodies. Although accumulated evidence suggests that the dysregulation of long non-coding RNAs (lncRNAs) is involved in the pathogenesis of SLE, the genetic contributions of lncRNA coding genes to SLE susceptibility remain largely unknown. Here, we aimed to provide more evidence for the role of lncRNA coding genes to SLE susceptibility.MethodsThe genetic association analysis was first adopted from the previous genome-wide association studies (GWAS) and was then validated in an independent cohort. PRDX6-AS1 is located at chr1:173204199-173446294. It spans a region of approximately 240 kb, and 297 single nucleotide polymorphisms (SNPs) were covered by the previous GWAS. Differential expression at the mRNA level was analyzed based on the ArrayExpress Archive database.ResultsA total of 33 SNPs were associated with SLE susceptibility, with a P<1.68×10-4. The strongest association signal was detected at rs844649 (P=2.12×10-6), according to the previous GWAS. Combining the results from the GWAS Chinese cohort and our replication cohort, we pursued a meta-analysis approach and found a pronounced genetic association between PRDX6-AS1 rs844649 and SLE susceptibility (pmeta=1.24×10-13, OR 1.50, 95% CI: 1.34–1.67). The mRNA expression of PRDX6 was elevated in peripheral blood cells, peripheral blood mononuclear cells (PBMCs), and multiple cell subpopulations, such as B cells, CD4+ T cells, CD3+ cells, and monocytes in patients with SLE. The PRDX6 protein expression level was also increased in patients with SLE compared with healthy donors.ConclusionOur study provides new evidence that variants located in lncRNA coding genes are associated with SLE susceptibility.

Genome-wide analyses using multi-locus models revealed marker-trait associations for major agronomic traits in Sorghum bicolor.

Globally, sorghum is the fifth most important cereal crop, and it is a major crop in Ethiopia, where it has a high genetic diversity. The country’s sorghum gene pool contributes significantly to sorghum improvement worldwide. This study aimed to identify genomic regions and candidate genes associated with major agronomic traits in sorghum by using its genetic resources in Ethiopia for a genome-wide association study (GWAS). Phenotypic data of days to flowering (DTF), plant height (PH), panicle length (PALH), panicle width (PAWD), panicle weight (PAWT), and grain yield (GY) were collected from a GWAS panel comprising 324 sorghum accessions grown in three environments. SeqSNP, a targeted genotyping method, was used to genotype the panel using 5,000 gene-based single nucleotide polymorphism (SNP) markers. For marker-trait association (MTA) analyses, fixed and random model circulating probability unification (FarmCPU), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK) models were used. In all traits, high phenotypic variation was observed, with broad-sense heritability ranging from 0.32 (for GY) to 0.90 (for PALH). A population structure, principal component analysis, and kinship analysis revealed that the accessions could be divided into two groups. In total, 54 MTAs were identified, 11 of which were detected by both BLINK and farmCPU. MTAs identified for each trait ranged from five (PAWT and GY) to fourteen (PH) representing both novel and previously identified quantitative trait loci (QTLs). Three SNPs were associated with more than one trait, including a SNP within the Sobic.004G189200 gene that was associated with PH and PAWT. Major effect SNP loci, Sbi2393610 (PVE = 23.3%), Sbi10438246 (PVE = 35.2%), Sbi17789352 (PVE = 11.9%) and Sbi30169733 (PVE = 18.9%) on chromosomes 1, 3, 5 and 9 that showed strong association signals for PAWD, DTF, GY and PALH, respectively, were major findings of this study. The SNP markers and candidate genes identified in this study provide insights into the genetic control of grain yield and related agronomic traits, and once validated, the markers could be used in genomics-led breeding.

A natural uORF variant confers phosphorus acquisition diversity in soybean

Phosphorus (P) is an essential element for all organisms. Because P fertilizers are a non-renewable resource and high fixation in soils, sustainable agriculture requires researchers to improve crop P acquisition efficiency. Here, we report a strong association signal at a locus of CPU1 (component of phosphorus uptake 1), from a genome-wide association study of P acquisition efficiency in a soybean core collection grown in the field. A SEC12-like gene, GmPHF1, is identified as the causal gene for CPU1. GmPHF1 facilitates the ER (endoplasmic reticulum) exit of the phosphate transporter, GmPT4, to the plasma membrane of root epidermal cells. A common SNP in an upstream open reading frame (uORF) of GmPHF1, which alters the abundance of GmPHF1 in a tissue-specific manner, contributes to P acquisition diversity in soybean. A natural genetic variation conditions diversity in soybean P acquisition, which can be used to develop P-efficient soybean genotypes.

Genome-Wide Association Study Reveals That PvGUX1_1 Is Associated with Pod Stringlessness in Snap Bean (Phaseolus vulgaris L.).

Simple SummaryUsing 138 snap bean accessions as plant materials, we investigated their suture strings across two years. With the goal of identifying the gene(s) responsible for the formation of suture strings, we conducted a genome-wide association study. A strong association signal was found in a 266.19 kb region on Chr02. Within the region, 23 candidate genes were identified. Importantly, the sequence and gene expression of PvGUX1_1 differed significantly between sutured pods and non-sutured pods. In addition, PvGUX1_1 was also a domesticated locus that diverged from PvGUX1_2 during an early stage. The results obtained in this study can provide important information for the improvement of pod quality in snap beans.Suture strings are a particularly important pod trait that determine the quality and texture of snap beans (Phaseolus vulgaris L.). The St locus on chromosome 2 has been described as a major locus associated with suture strings. However, the gene and genetic basis underlying this locus remain unknown. Here, we investigated the suture strings of 138 snap bean accessions across two years. A total of 3.66 million single-nucleotide polymorphisms (SNPs) were obtained by deep resequencing. Based on these SNPs, we identified a strong association signal on Chr02 and a promising candidate gene, PvGUX1_1. Further analysis revealed that the 2 bp deletion in the exon of PvGUX1_1 was significantly associated with stringlessness. Comparative mapping indicated that PvGUX1_1 was a domesticated locus and diverged from PvGUX1_2 during an early stage. Our study provides important insights into the genetic mechanism of suture string formation and useful information for snap bean improvement.

GWAS in people of Middle Eastern descent reveals a locus protective of kidney function\u2014a cross-sectional study

BackgroundType 2 diabetes is one of the leading causes of chronic kidney failure, which increases globally and represents a significant threat to public health. People from the Middle East represent one of the largest immigrant groups in Europe today. Despite poor glucose regulation and high risk for early-onset insulin-deficient type 2 diabetes, they have better kidney function and lower rates of all-cause and cardiovascular-specific mortality compared with people of European ancestry. Here, we assessed the genetic basis of estimated glomerular filtration rate (eGFR) and other metabolic traits in people of Iraqi ancestry living in southern Sweden.MethodsGenome-wide association study (GWAS) analyses were performed in 1201 Iraqi-born residents of the city of Malmö for eGFR and ten other metabolic traits using linear mixed-models to account for family structure.ResultsThe strongest association signal was detected for eGFR in CST9 (rs13037490; P value = 2.4 × 10−13), a locus previously associated with cystatin C-based eGFR; importantly, the effect (major) allele here contrasts the effect (minor) allele in other populations, suggesting favorable selection at this locus. Additional novel genome-wide significant loci for eGFR (ERBB4), fasting glucose (CAMTA1, NDUFA10, TRIO, WWC1, TRAPPC9, SH3GL2, ABCC11), quantitative insulin-sensitivity check index (METTL16), and HbA1C (CAMTA1, ME1, PAK1, RORA) were identified.ConclusionsThe genetic effects discovered here may help explain why people from the Middle East have better kidney function than those of European descent. Genetic predisposition to preserved kidney function may also underlie the observed survival benefits in Middle Eastern immigrants with type 2 diabetes.

Genetic Analyses of Enamel Hypoplasia in Multiethnic Cohorts

Introduction: Enamel hypoplasia causes a reduction in the thickness of affected enamel and is one of the most common dental anomalies. This defect is caused by environmental and/or genetic factors that interfere with tooth formation, emphasizing the importance of investigating enamel hypoplasia on an epidemiological and genetic level. Methods: A genome-wide association of enamel hypoplasia was performed in multiple cohorts, overall comprising 7,159 individuals ranging in age from 7 to 82 years. Mixed models were used to test for genetic association while simultaneously accounting for relatedness and genetic population structure. Meta-analysis was then performed. More than 5 million single-nucleotide polymorphisms were tested in individual cohorts. Results: Analyses of the individual cohorts and meta-analysis identified association signals close to genome-wide significance (p < 5 × 10^–8), and many suggestive association signals (5 × 10^–8 < p < 5 × 10^–6) near genes with plausible roles in tooth/enamel development. Conclusion: The strongest association signal (p = 1.57 × 10^–9) was observed near BMP2K in one of the individual cohorts. Additional suggestive signals were observed near genes with plausible roles in tooth development in the meta-analysis, such as SLC4A4 which can influence enamel hypoplasia. Additional human genetic studies are needed to replicate these results and functional studies in model systems are needed to validate our findings.

Age-related differences of genetic susceptibility to patients with acute lymphoblastic leukemia.

Inherited predispositions to acute lymphoblastic leukemia have been well investigated in pediatric patients, but studies on adults, particularly Chinese patients, are limited. In this study, we conducted a genome-wide association study in 466 all-age Chinese patients with Acute lymphoblastic leukemia (ALL) and 1,466 non-ALL controls to estimate the impact of age on ALL susceptibility in the Chinese population. Among the 17 reported loci, 8 have been validated in pediatric and 1 in adult patients. The strongest association signal was identified at ARID5B locus and gradually decreased with age, while the signal at GATA3 exhibited the opposite trend and significantly impact on adult patients. With genome-wide approaches, germline variants at 2q14.3 rank as the top inherited predisposition to adult patients (e.g., rs73956024, P = 4.3 × 10-5) and separate the genetic risk of pediatric vs. adult patients (P = 3.6 × 10-6), whereas variants at 15q25.3 (e.g., rs11638062) have a similar impact on patients in different age groups (overall P = 2.9 × 10-7). Our analysis highlights the impact of age on genetic susceptibility to ALL in Chinese patients.

Association mapping of resistance to emerging stem rust pathogen races in spring wheat using genotyping-by-sequencing.

The identification and characterization of resistance genes should outpace the rapid emergence of new P. graminis f. sp. tritici races, such as TTRTF and TTKTT, to mitigate stem rust damage to wheat. The objective of the current study was to identify and characterize P. graminis f. sp. tritici race resistance association signals. A total of 250 North American spring wheat lines were evaluated at the seedling stage with a total of seven isolates including TKKTP, TKTTF, TKTTF, TRTTF, TTRTF, TTKSK, and TTKTT. The lines were genotyped by a GBS platform and 9,042 SNPs were used for identification of chromosome regions associated with resistance against the seven isolates. Strong association signals were detected on chromosomes 6BL (Sr11 gene region) and 4AL, likely Sr7a, for resistance against both TKKTP and TKTTF. Similarly, association signals were also detected on chromosomes 4AL (race TTRTF resistance) and 4BS (race TTKSK and TTKTT resistance). Association analysis based on mean phenotypic differences between closely related isolates identified QTL that were not elucidated by direct association mapping of the responses, individually. Overall, with the exception of race TRTTF, each race shared at least one association signal with another race. However, the number of race-specific association signals are larger than that of association signals common among races suggesting the need for identifying and characterizing QTL/genes for newly emerging stem rust pathogen races. There was also high concordance between PCA-based GWAS association signals and association signals from that of both single and multi-locus mixed models.

Genetic Analysis Reveals a Significant Contribution of CES1 to Prostate Cancer Progression in Taiwanese Men.

The genes that influence prostate cancer progression remain largely unknown. Since the carboxylesterase gene family plays a crucial role in xenobiotic metabolism and lipid/cholesterol homeostasis, we hypothesize that genetic variants in carboxylesterase genes may influence clinical outcomes for prostate cancer patients. A total of 478 (36 genotyped and 442 imputed) single nucleotide polymorphisms (SNPs) in five genes of the carboxylesterase family were assessed in terms of their associations with biochemical recurrence (BCR)-free survival in 643 Taiwanese patients with prostate cancer who underwent radical prostatectomy. The strongest association signal was shown in CES1 (P = 9.64 × 10−4 for genotyped SNP rs8192935 and P = 8.96 × 10−5 for imputed SNP rs8192950). After multiple test correction and adjustment for clinical covariates, CES1 rs8192935 (P = 9.67 × 10−4) and rs8192950 (P = 9.34 × 10−5) remained significant. These SNPs were correlated with CES1 expression levels, which in turn were associated with prostate cancer aggressiveness. Furthermore, our meta-analysis, including eight studies, indicated that a high CES1 expression predicted better outcomes among prostate cancer patients (hazard ratio 0.82, 95% confidence interval 0.70–0.97, P = 0.02). In conclusion, our findings suggest that CES1 rs8192935 and rs8192950 are associated with BCR and that CES1 plays a tumor suppressive role in prostate cancer.

Exome sequencing study revealed novel susceptibility loci in subarachnoid hemorrhage (SAH)

AimTo expand our current understanding of the genetic basis of subarachnoid hemorrhage (SAH), and reveal the susceptibility genes in SAH risk.MethodsWe conducted whole-exome sequencing (WES) in a cohort of 196 individuals, including 94 SAH patients and 94 controls, as well as 8 samples that belong to two pedigrees. Systematically examination for rare variations (through direct genotyping) and common variations (through genotyping and imputation) for SAHs were performed in this study.ResultsA total of 16,029 single-nucleotide polymorphisms (SNPs) and 108,999 short indels were detected in all samples, and among them, 30 SNPs distributed on 17 genes presented a strong association signal with SAH. Two novel pathogenic gene variants were identified as associated risk loci, including mutation in TPO and PALD1. The statistical analysis for rare, damaging variations in SAHs identified several susceptibility genes which were involved in degradation of the extracellular matrix and transcription factor signal pathways. And 25 putative pathogenic genes for SAH were also identified basic on functional interaction network analysis with the published SAH-associated genes. Additionally, pedigree analysis revealed autosomal dominant inheritance of pathogenic genes.ConclusionSystematical analysis revealed a key role for rare variations in SAH risk and discovered SNPs in new complex loci. Our study expanded the list of candidate genes associated with SAH risk, and will facilitate the investigation of disease-related mechanisms and potential clinical therapies.

QTL mapping and haplotype analysis revealed candidate genes for grain thickness in rice (Oryza sativa L.)

Grain size and shape are important for grain yield and quality in rice (Oryza sativa L.). Grain thickness (GT) is one of the constituent characteristics of grain shape and has a strong influence on grain weight in rice. In this study, a japonica cultivar (02428) with a thick grain phenotype was crossed with an indica cultivar (YZX) with a thin grain phenotype to construct a set of 192 recombinant inbred lines (RILs). Quantitative trait locus (QTL) analysis was performed with a high-density genetic map harbouring 2711 bin markers. We identified 6 QTLs for GT in three environments and two stable QTLs with high heritability residing on the known grain width (GW) genes GW5 and GL7/GW7, verifying the effectiveness and reliability of our study. We detected a novel QTL on chromosome 9 with strong association signals for both GW and GT. Using transcript profiles and quantitative RT-PCR, we identified a candidate gene, Os09g0535500, which encodes a zinc finger C3HC4-type domain, that showed the highest expression level and significant differential expression between parents in the development of young panicle and grain. Further haplotype-phenotype association analysis revealed that Os09g0535500 was strongly correlated with GW and GT. Thus, our results indicated that Os09g0535500 was the most promising candidate gene in WTG9 and is likely important for rice grain development, which lays a foundation for further functional validation and breeding utilization.

An integrative, genomic, transcriptomic and network-assisted study to identify genes associated with human cleft lip with or without cleft palate

BackgroundCleft lip with or without cleft palate (CL/P) is one of the most common congenital human birth defects. A combination of genetic and epidemiology studies has contributed to a better knowledge of CL/P-associated candidate genes and environmental risk factors. However, the etiology of CL/P remains not fully understood. In this study, to identify new CL/P-associated genes, we conducted an integrative analysis using our in-house network tools, dmGWAS [dense module search for Genome-Wide Association Studies (GWAS)] and EW_dmGWAS (Edge-Weighted dmGWAS), in a combination with GWAS data, the human protein-protein interaction (PPI) network, and differential gene expression profiles.ResultsA total of 87 genes were consistently detected in both European and Asian ancestries in dmGWAS. There were 31.0% (27/87) showed nominal significance with CL/P (gene-based p < 0.05), with three genes showing strong association signals, including KIAA1598, GPR183, and ZMYND11 (p < 1 × 10− 3). In EW_dmGWAS, we identified 253 and 245 module genes associated with CL/P for European ancestry and the Asian ancestry, respectively. Functional enrichment analysis demonstrated that these genes were involved in cell adhesion, protein localization to the plasma membrane, the regulation of the apoptotic signaling pathway, and other pathological conditions. A small proportion of genes (5.1% for European ancestry; 2.4% for Asian ancestry) had prior evidence in CL/P as annotated in CleftGeneDB database. Our analysis highlighted nine novel CL/P candidate genes (BRD1, CREBBP, CSK, DNM1L, LOR, PTPN18, SND1, TGS1, and VIM) and 17 previously reported genes in the top modules.ConclusionsThe genes identified through superimposing GWAS signals and differential gene expression profiles onto human PPI network, as well as their functional features, helped our understanding of the etiology of CL/P. Our multi-omics integrative analyses revealed nine novel candidate genes involved in CL/P.

Genome-Wide Association Study Identifies an Immune-Related Etiology for Severe Aplastic Anemia

Introduction. Acquired severe aplastic anemia (SAA) is a life-threatening disorder characterized by severe progressive pancytopenia and hypocellular bone marrow. The etiology of acquired SAA is not understood but believed to be related to abnormal immune responses to environmental exposures. We conducted a genome-wide association study (GWAS) to identify common germline variants associated with SAA. Methods. We identified 895 patients with SAA who underwent related or unrelated hematopoietic cell transplant (HCT) with clinical data and pre-HCT blood samples available in the Center for International Blood and Marrow Transplant Research (CIBMTR) database and biorepository. Pre-HCT DNA was extracted from blood of patients with SAA and genome-wide genotyping was conducted using the Illumina OmniExpress array. We excluded 93 inherited bone marrow failure cases. The SAA cases were grouped cases into discovery and validation sets based on time of batch sample receipt. Analyses were limited to patients of European ancestry based on principal component analyses to minimize the potential effect of population stratification. Controls were genomically matched and selected from previously scanned cancer-free subjects at the Cancer Genomics Research Laboratory, Division of Cancer Epidemiology and Genetics, NCI. The final analysis included 534 acquired SAA cases (359 in the discovery set and 175 in the validation set), and 2,455 controls (1,396 in the discovery set, and 1,059 in the validation set). Results. Patients with SAA in this study received HCT between 1989-2015 at a median age of 21 years, 56% were male, and the median time between SAA diagnosis and HCT was 11 months. Strong genome-wide association signals were identified across the human leukocyte antigen (HLA) genes encoded at the major histocompatibility complex (MHC) on chromosome 6p21 (Figure 1). The top SNP was located in the P4 binding pocket of the HLA class II gene HLA-DPB1(rs1042151A&gt;G, p.Met105Val, pooled-odds ratio [OR]=1.75, 95% confidence interval [CI]=1.50-2.03, p=1.94x10-13). The expression of HLA-DP in CD19+ cells from 175 healthy donors was significantly different by rs1042151 A&gt;G genotype (p=2.04x10-6) (Figure 2). A second SNP near HLA-B, rs28367832G&gt;A, also reached genome-wide significance (pooled-OR=1.49, 95% CI=1.22-1.78, p=7.27x10-9). Copy-number variant analysis and next generation sequencing also identified somatic, clonal copy-neutral loss-of-heterozygosity affecting class I HLA genes in 8.6% of the SAA cases and none of the controls. Conclusion. This SAA GWAS identified strong association signals between common germline genetic variants in HLA class I and II genes and SAA. The main SNP is associated with changes in HLA-DP expression suggesting a key role for this locus in SAA etiology. This study adds further evidence to the connection between SAA and immune dysregulation. Disclosures Cerhan: Janssen: Membership on an entity's Board of Directors or advisory committees; NanoString: Research Funding; Celgene: Research Funding. Lee:Incyte: Research Funding; Syndax: Research Funding; Amgen: Research Funding; Novartis: Research Funding; Takeda: Research Funding; Kadmon: Consultancy, Research Funding; Pfizer: Consultancy, Research Funding; AstraZeneca: Research Funding.

Single‐marker and haplotype‐based association analysis of anthracnose (Colletotrichum dematium)resistance in spinach (Spinacia oleracea)

AbstractAnthracnose (Colletotrichum dematium) is an important disease in spinach (Spinacia oleracea). Sources of resistance must be identified, and molecular tools must be developed to expedite cultivar development. In this study, a diverse collection of 276 spinach accessions was scored for anthracnose disease severity. We then evaluated marker identification approaches by testing how well haplotype‐based trait modelling compares to single markers in identifying strong association signals. Alleles in linkage disequilibrium were tagged in haplotype blocks, and anthracnose‐associated molecular markers were identified using single‐SNP (sSNP), pairwise haplotype (htP) and multi‐marker haplotype (htM) SNP tagging approaches. We identified 49 significantly associated markers distributed on several spinach chromosomes using all methods. The sSNP approach identified 13 markers, while htP identified 24 (~63% more) and htM 34 (~162% more). Of these markers, four were uniquely identified by the sSNP approach, nine by htP and nineteen by htM. The results indicate that resistance to anthracnose is polygenic and that haplotype‐based analysis may have more power than sSNP. Using a combination of these methods can improve the identification of molecular markers for spinach breeding.

Genome-wide association analysis of autism identified multiple loci that have been reported as strong signals for neuropsychiatric disorders.

Autism is a common neurodevelopmental disorder with a moderate to a high degree of heritability, but only a few common genetic variants that explain the heritability have been associated. We performed a genome-wide transmission disequilibrium test analysis of a newly genotyped autism case-parent triad samples (127 trios) in Han Chinese, identified top association signals at multiple single nucleotide polymorphisms (SNPs), including rs9839376 (OR = 2.59, P = 1.27 × 10-05 ) at KCNMB2, rs6044680 (OR = 0.319, P = 4.82 × 10-05 ) and rs7274133 (OR = 0.313, P = 3.22 × 10-05 ) at PCSK2, and rs310619 (OR = 2.40, P = 7.44 × 10-05 ) at EEF1A2. Furthermore, a genome-wide combined P-value of individual SNPs in two independent case-parent triad samples (total 402 triads, n = 1,206) identified SNPs at EGFLAM, ZDHHC2, AGBL1, and SNX29 as additional association signals for autism. While none of these signals achieved a genome-wide significance in the two samples of our study, they have been reported in a previous genome-wide association study of neuropsychiatric disorders, and the majority of these SNP have a significant cis-regulatory association with mRNA in human tissues (false discovery rate (FDR) < 0.05). Our study warrants further study or replication with additional sample for association with autism and other neuropsychiatric disorders. Autism Res 2020, 13: 382-396. © 2019 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: Autism is a common neurodevelopmental disorder, heritable, but only a few common genetic variants that explain the heritability have been associated. We conducted a genome-wide association study with two cohorts of autism case-parent triad samples in Han Chinese and identified multiple single nucleotide polymorphisms that were reported as strong association signals in a previous genome-wide association study of other neuropsychiatric disorders or related traits. Our study provides evidence for shared genetic variants among autism and other neuropsychiatric disorders.

Genomewide Association Study Confirming the Association of NAT2 with Susceptibility to Antituberculosis Drug-Induced Liver Injury in Thai Patients.

Antituberculosis drug-induced liver injury (ATDILI) is a common side effect leading to tuberculosis (TB) treatment disruption. The mechanism of the disease remains poorly understood. We conducted a genomewide association study (GWAS) to investigate all possible genetic factors of ATDILI in Thai patients. This study was carried out in Thai TB patients, including 79 ATDILI cases and 239 tolerant controls from our network hospitals in Thailand. Nearly 1 million single-nucleotide polymorphisms (SNPs) were genotyped across the whole genome using an Illumina OmniExpress Exome BeadChip array. In the discovery stage, we identified strong association signals on chromosome 8 originating from the N-acetyltransferase (NAT2) region. The A allele of rs1495741, the top SNP in the intergenic region of NAT2 and PSD3 (14 kb from NAT2), was significantly associated with ATDILI (recessive model, odds ratio of 6.01 [95% confidence interval, 3.42 to 10.57]; P = 6.86E-11). This particular SNP was reported as a tag SNP for NAT2 inferred phenotypes. The AA, AG, and GG genotypes represented NAT2 slow acetylators, intermediate acetylators, and fast acetylators, respectively. The tag SNP genotypes demonstrated a concordance rate of 94.98% with NAT2 acetylator phenotypes. This GWAS shows that NAT2 is the most important risk factor for ATDILI in the Thai population.

Serum Antibody Profiling Healthy Toddlers Reveals Low, Intermediate and High Response Cohorts

Abstract During their first two years of life, infant/toddlers produce their own antibodies in response to diverse environmental exposures, infections, and vaccinations. The specificity of these antibodies towards self-antigens, infections, and vaccine components and how they vary among individuals is poorly characterized. We developed a new antigen array comprising autoantigens, infectious agents, and vaccine antigens to assess the serum antibody specificities among a cohort of &amp;gt;150 health toddlers. An analysis of their responses reveals a stratification of healthy toddlers into 3 groups, low, intermediate, and high responders. Among those 16% that are high responders, their serum IgGs were able to bind a diverse array of self-antigens as well as infectious agents. Longitudinal follow-up suggests this pattern is relatively stable over time. Comparing clinical data reveals a significant correlation with the high IgG responder cohort and a family history of asthma and maternal gestational diabetes. Targeted DNA sequencing in the high responder group revealed a strong genetic association signal at the HLA locus, with polymorphisms at this locus associated with high ANA and IgG titers to many antigens. Findings from this study may provide insights into the natural history of human autoantibody formation.