Genome-wide Linkage Disequilibrium Research Articles

Shared genomic architectures of COVID-19 and antisocial behavior

Little is known about the genetics of norm violation and aggression in relation to coronavirus disease 2019 (COVID-19). To investigate this, we used summary statistics from genome-wide association studies and linkage disequilibrium score regression to calculate a matrix of genetic correlations (rgs) for antisocial behavior (ASB), COVID-19, and various health and behavioral traits. After false-discovery rate correction, ASB was genetically correlated with COVID-19 (rg = 0.51; P = 1.54E-02) and 19 other traits. ASB and COVID-19 were both positively genetically correlated with having a noisy workplace, doing heavy manual labor, chronic obstructive pulmonary disease, and genitourinary diseases. ASB and COVID-19 were both inversely genetically correlated with average income, education years, healthspan, verbal reasoning, lifespan, cheese intake, and being breastfed as a baby. But keep in mind that rgs are not necessarily causal. And, if causal, their prevailing directions of effect (which causes which) are indiscernible from rgs alone. Moreover, the SNP-heritability (h_g^2) estimates for two measures of COVID-19 were very small, restricting the overlap of genetic variance in absolute terms between ASB and COVID-19. Nonetheless, our findings suggest that those with antisocial tendencies possibly have a higher risk of exposure to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) than those without antisocial tendencies. This may have been especially true early in the pandemic before vaccines against SARS-CoV-2 were available and before the emergence of the highly transmissible Omicron variant.

Genotyping-by-Sequencing Based Molecular Genetic Diversity of Pakistani Bread Wheat (Triticum aestivum L.) Accessions.

Spring wheat (Triticum aestivum L.) is one of the most imperative staple food crops, with an annual production of 765 million tons globally to feed ∼40% world population. Genetic diversity in available germplasm is crucial for sustainable wheat improvement to ensure global food security. A diversity panel of 184 Pakistani wheat accessions was genotyped using 123,596 high-quality single nucleotide polymorphism (SNP) markers generated by genotyping-by-sequencing with 42% of the SNPs mapped on B, 36% on A, and 22% on D sub-genomes of wheat. Chromosome 2B contains the most SNPs (9,126), whereas 4D has the least (2,660) markers. The mean polymorphic information content, genetic diversity, and major allele frequency of the population were 0.157, 0.1844, and 0.87, respectively. Analysis of molecular variance revealed a higher genetic diversity (80%) within the sub-population than among the sub-populations (20%). The genome-wide linkage disequilibrium was 0.34 Mbp for the whole wheat genome. Among the three subgenomes, A has the highest LD decay value (0.29 Mbp), followed by B (0.2 Mbp) and D (0.07 Mbp) genomes, respectively. The results of population structure, principal coordinate analysis, phylogenetic tree, and kinship analysis also divided the whole population into three clusters comprising 31, 33, and 120 accessions in group 1, group 2, and group 3, respectively. All groups were dominated by the local wheat accessions. Estimation of genetic diversity will be a baseline for the selection of breeding parents for mutations and the genome-wide association and marker-assisted selection studies.

Genome-wide association studies in non-anxiety individuals identified novel risk loci for depression.

Depression is a debilitating mental disorder that often coexists with anxiety. The genetic mechanisms of depression and anxiety have considerable overlap, and studying depression in non-anxiety samples could help to discover novel gene. We assess the genetic variation of depression in non-anxiety samples, using genome-wide association studies (GWAS) and linkage disequilibrium score regression (LDSC). The GWAS of depression score and self-reported depression were conducted using the UK Biobank samples, comprising 99,178 non-anxiety participants with anxiety score <5 and 86,503 non-anxiety participants without self-reported anxiety, respectively. Replication analysis was then performed using two large-scale GWAS summary data of depression from Psychiatric Genomics Consortium (PGC). LDSC was finally used to evaluate genetic correlations with 855 health-related traits based on the primary GWAS. Two genome-wide significant loci for non-anxiety depression were identified: rs139702470 (p=1.54×10-8, OR=0.29) locate in PIEZO2, and rs6046722 (p=2.52×10-8, OR=1.09) locate in CFAP61. These associated genes were replicated in two GWAS of depression from PGC, such as rs1040582 (preplication GWAS1=0.02, preplication GWAS2=2.71×10-3) in CFAP61, and rs11661122 (preplication GWAS1=8.16×10-3, preplication GWAS2=8.08×10-3) in PIEZO2. LDSC identified 19 traits genetically associated with non-anxiety depression (p<0.001), such as marital separation/divorce (rg=0.45, SE=0.15). Our findings provide novel clues for understanding of the complex genetic architecture of depression.

Genetic diversity and population structure of rice (Oryza sativa L.) landraces from Kerala, India analyzed through genotyping-by-sequencing.

Researchers stand at the vanguard of advancement and application of next-generation sequencing technology for developing dominant strategies for the sustainable management of genetically diverse crops. We attempt to fill the existing research lacuna in the molecular characterization of potent rice landraces in Kerala. Genotyping-by-sequencing (GBS) was performed on 96 Kerala rice accessions to identify single-nucleotide polymorphisms (SNPs), to examine the genetic diversity, population structure, and to delineate linkage disequilibrium (LD) pattern. GBS identified 5856 high-quality SNPs. The structure analysis indicated three subpopulations with the highest probability for population clustering with significant genetic differentiation, confirmed by principal component analysis. The genome-wide LD decay distance was 772kb, at which the r2 dropped to half its maximum value. The analysis of genetic properties of the identified SNP panel with an average polymorphism information content (PIC) value of 0.22 and a minor allele frequency (MAF) > 0.1 unveiled their efficacy in genome-wide association studies (GWAS). High FST (0.266) and low Nm (0.692) portray a strong genetic differentiation among the rice landraces, complementing the genetic structuring observed in the studied population. Slow LD decay in the rice landraces reflects their self-pollinating behavior and the indirect selection of desired traits by domestication. Moreover, the high LD entails only a minimum number of SNP markers for detecting marker-trait association. The diverse germplasm utilized in this study can be further utilized to disclose genetic variants associated with phenotypic traits and define signatures of selection via GWAS and selective sweep, respectively.

Genotyping-by-sequencing and genomic selection applications in hexaploid triticale.

Triticale, a hybrid species between wheat and rye, is one of the newest additions to the plant kingdom with a very short history of improvement. It has very limited genomic resources because of its large and complex genome. Objectives of this study were to generate dense marker data, understand genetic diversity, population structure, linkage disequilibrium (LD), and estimate accuracies of commonly used genomic selection (GS) models on forage yield of triticale. Genotyping-by-sequencing (GBS), using PstI and MspI restriction enzymes for reducing genome complexity, was performed on a triticale diversity panel (n = 289). After filtering for biallelic loci with more than 70% genome coverage, and minor allele frequency (MAF) > 0.05, de novo variant calling identified 16,378 single nucleotide polymorphism (SNP) markers. Sequences of these variants were mapped to wheat and rye reference genomes to infer their homologous groups and chromosome positions. About 45% (7430), and 58% (9500) of the de novo identified SNPs were mapped to the wheat and rye reference genomes, respectively. Interestingly, 28.9% (2151) of the 7430 SNPs were mapped to the D genome of hexaploid wheat, indicating substantial substitution of the R genome with D genome in cultivated triticale. About 27% of marker pairs were in significant LD with an average r2 > 0.18 (P < 0.05). Genome-wide LD declined rapidly to r2 < 0.1 beyond 10 kb physical distance. The three sub-genomes (A, B, and R) showed comparable LD decay patterns. Genetic diversity and population structure analyses identified five distinct clusters. Genotype grouping did not follow prior winter vs spring-type classification. However, one of the clusters was largely dominated by winter triticale. GS accuracies were estimated for forage yield using three commonly used models with different training population sizes and marker densities. GS accuracy increased with increasing training population size while gain in accuracy tended to plateau with marker densities of 2000 SNPs or more. Average GS accuracy was about 0.52, indicating the potential of using GS in triticale forage yield improvement.

Genome-Wide Association Mapping of Late Blight Tolerance Trait in Potato (Solanum tuberosum L.).

Uncovering the genetic basis and optimizing the late blight tolerance trait in potatoes (Solanum tuberosum L.) are crucial for potato breeding. Late blight disease is one of the most significant diseases hindering potato production. The traits of late blight tolerance were evaluated for 284 potato cultivars to identify loci significantly associated with the late blight tolerance trait. Of all, 37 and 15 were the most tolerant to disease, and 107 and 30 were the most susceptible. A total of 22,489 high-quality single-nucleotide polymorphisms and indels were identified in 284 potato cultivars. All the potato cultivars were clustered into eight subgroups using population structure analysis and principal component analysis, which were consistent with the results of the phylogenetic tree analysis. The average genetic diversity for all 284 potato cultivars was 0.216, and the differentiation index of each subgroup was 0.025–0.149. Genome-wide linkage disequilibrium (LD) analysis demonstrated that the average LD was about 0.9 kb. A genome-wide association study using a mixed linear model identified 964 loci significantly associated with the late blight tolerance trait. Fourteen candidate genes for late blight tolerance traits were identified, including genes encoding late blight tolerance protein, chitinase 1, cytosolic nucleotide-binding site–leucine-rich repeat tolerance protein, protein kinase, ethylene-responsive transcription factor, and other potential plant tolerance-related proteins. This study provides novel insights into the genetic architecture of late blight tolerance traits and will be helpful for late blight tolerance in potato breeding.

Genome-wide association study reveals a quantitative trait locus and two candidate genes on Sus scrofa chromosome 5 affecting intramuscular fat content in Suhuai pigs

Intramuscular fat content (IFC) is an essential quantitative trait of meat, affecting multiple meat quality indicators. A certain amount of IFC could not only improve the sensory score of pork but also increase the flavour, tenderness, juiciness and shelf-life. To dissect the genetic determinants of IFC, two methods, including genome-wide efficient mixed-model analysis (GEMMA) and linkage disequilibrium adjusted kinships (LDAKs), were used to carry out genome-wide association studies for IFC in Suhuai pig population. A total of 14 and 18 significant single nucleotide polymorphisms (SNPs) were identified by GEMMA and LDAK, respectively. The results of these two methods were highly consistent and all 14 significant SNPs in GEMMA were detected by LDAK. Seven of the 18 SNPs reached the genome-wide significance level (P < 9.85E−07) while 11 cases reached the suggestive significance level (P < 1.77E−05). These significant SNPs were mainly distributed on Sus scrofa chromosome (SSC) 5, 3, and 7. Moreover, one locus resides in a 2.27 Mb (71.37–73.64 Mb) region on SSC5 harbouring 13 significant SNPs associated with IFC, and the lead SNP (rs81302978) also locates in this region. Linkage disequilibrium (LD) analysis showed that there were four pairs of complete LD (r2 = 1) among these 13 SNPs, and the remaining 9 SNPs with incomplete LD (r2 ≠ 1) were selected for subsequent analyses of IFC. Association analyses showed that 7 out of 9 SNPs were significantly associated with IFC (P < 0.05) in 330 Suhuai pigs, and the other 2 SNPs tended to reach a significant association level with IFC (P < 0.1). The phenotypic variance explained (PVE) range of these 9 SNPs was 0.92–3.55%. Meanwhile, the lead SNP was also significantly associated (rs81302978) with IFC (P < 0.05) in 378 commercial hybrid pigs (Pietrain × Duroc) × (Landrace × Yorkshire) (PDLY), and the PVE was 1.38%. Besides, two lipid metabolism-relevant candidate genes, the leucine rich repeat kinase 2 (LRRK2) and PDZ domain containing ring finger 4 (PDZRN4) were identified in the 2.27 Mb region on SSC5. In conclusion, our results may provide a set of markers useful for genetic improvement of IFC in pigs and will advance the genome selection process of IFC on pig breeding programmes.

A Novel Method for Mendelian Randomization Analyses With Pleiotropy and Linkage Disequilibrium in Genetic Variants From Individual Data.

Mendelian randomization makes use of genetic variants as instrumental variables to eliminate the influence induced by unknown confounders on causal estimation in epidemiology studies. However, with the soaring genetic variants identified in genome-wide association studies, the pleiotropy, and linkage disequilibrium in genetic variants are unavoidable and may produce severe bias in causal inference. In this study, by modeling the pleiotropic effect as a normally distributed random effect, we propose a novel mixed-effects regression model-based method PLDMR, pleiotropy and linkage disequilibrium adaptive Mendelian randomization, which takes linkage disequilibrium into account and also corrects for the pleiotropic effect in causal effect estimation and statistical inference. We conduct voluminous simulation studies to evaluate the performance of the proposed and existing methods. Simulation results illustrate the validity and advantage of the novel method, especially in the case of linkage disequilibrium and directional pleiotropic effects, compared with other methods. In addition, by applying this novel method to the data on Atherosclerosis Risk in Communications Study, we conclude that body mass index has a significant causal effect on and thus might be a potential risk factor of systolic blood pressure. The novel method is implemented in R and the corresponding R code is provided for free download.

Chromosome level reference of Atlantic halibut Hippoglossus hippoglossus provides insight into the evolution of sexual determination systems.

Changes in the genetic mechanisms that control sexual determination have occurred independently across the tree of life, and with exceptional frequency in teleost fishes. To investigate the genomic changes underlying the evolution of sexual determination, we sequenced a chromosome-level genome, multitissue transcriptomes, and reduced representation population data for the Atlantic halibut (Hippoglossus hippoglossus), which has an XY/XX sex determination mechanism and has recently diverged (0.9-3.8Ma) from the Pacific halibut (Hippoglossus stenolepis), which has a ZZ/ZW system. We used frequency and coverage-based population approaches to identify a putative sex-determining factor, GSDF. We characterized regions with elevated heterozygosity and linkage disequilibrium indicating suppression of recombination across a nascent sex chromosome. We detected testis-specific expression of GSDF, the sequence of which is highly conserved across flatfishes. Based on evidence from genome-wide association, coverage, linkage disequilibrium, testis and brain transcriptomes, and sequence conservation with other flatfishes, we propose a mechanism for the recent evolution of an XY sex-determination mechanism in Atlantic halibut. Changes to the ancestral sex-determining gene DMRT1 in regulating the downstream gene GSDF probably coincided with GSDF, or a proximal regulatory element of it, becoming the primary sex-determining factor. Our results suggest changes to a small number of elements can have drastic repercussions for the genomic substrate available to sex-specific evolutionary forces, providing insight into how certain elements repeatedly evolve to control sex across taxa. Our chromosome-level assembly, multitissue transcriptomes, and population genomic data provide a valuable resource and understanding of the evolution of sexual systems in fishes.

Drivers of linkage disequilibrium across a species' geographic range.

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species' range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Genome-wide association analysis of fiber fineness and yield in ramie (Boehmeria nivea) using SLAF-seq

Ramie (Boehmeria nivea L. Gaudich) is a major natural fiber crop cultivated in East Asia. The improvement of its fiber yield and fineness are important breeding goals. Fiber yield is a complex quantitative trait comprising ramet number, stem diameter, plant height, skin thickness, and fiber percentage. The fiber fineness is a crucial trait for ramie fiber quality. However, there are few association studies for fiber yield traits and fiber fineness on ramie, and lack high-density SNP maps in natural ramie population. Here, a panel of 112 core ramie germplasms were genotyped by 215,376 consistent single-nucleotide polymorphisms (SNPs) from specific-locus amplified-fragment sequencing (SLAF-seq), and used for genome-wide association study of fiber fineness and yield. Subsequently, the genetic diversity, linkage disequilibrium (LD), and population structure was conducted based on 215,376 SNPs. Population cluster analysis disclosed five subpopulations. Neighbor-joining (NJ) analysis revealed three major clusters. No obvious relationships were identified between them and their geographic origins. The genome-wide LD decayed to r2 = 0.1 was ~ 11.75 kbp in physical distance. One, seven, one, seven, and twenty-seven significant SNP marker associations were detected for fiber fineness (third season), stem diameter (third season), stem skin thickness (third season), fiber percentage (second season), and fiber percentage (third season), respectively. Two promising candidate genes, whole_GLEAN_10029622 and whole_GLEAN_10029638 resided in the significant trait-SNP association for fiber fineness (third season), which annotated as a cotton fiber-expressed protein and an Arabidopsis thaliana homebox protein ATH1, respectively and validated by qPCR. The identified loci or genes for fiber fineness and yield may provide the basis for future research on fiber fineness and yield and marker-assisted selection breeding for ramie.

Genome-Wide Association Study of Tan Spot Resistance in a Hexaploid Wheat Collection From Kazakhstan.

Tan spot, caused by Pyrenophora tritici-repentis, is a serious foliar disease of wheat in Kazakhstan with reported yield losses as high as 50% during epidemic years. Here, we report the evaluation of a collection of 191 hexaploid spring and winter wheat lines for tan spot resistance and its underlying genetic architecture using genome-wide association study (GWAS). Our wheat collection comprised candidate varieties from Kazakhstan, Russia, and CIMMYT. It was genotyped using the DArTseq technology and phenotyped for resistance to tan spot at seedling and adult plant stages in Kazakhstan. DArTseq SNPs revealed high genetic diversity (average polymorphic information content = 0.33) in the panel and genome-wide linkage disequilibrium decay at 22 Mb (threshold r2 = 0.1). Principal component analysis revealed a clear separation of Eurasian germplasm from CIMMYT and IWWIP lines. GWAS identified 34 marker-trait associations (MTA) for resistance to tan spot and the amount of phenotypic variation explained by these MTA ranged from 4% to 13.7%. Our results suggest the existence of novel valuable resistant alleles on chromosomes 3BS, and 5DL and 6AL for resistance to Race 1 and Race 5, respectively, in addition to known genes tsn1 and tsc2. On chromosome 6AL, a genomic region spanning 3 Mb was identified conferring resistance to both Race 1 and Race 5. Epistatic interaction of associated loci was revealed on chromosomes 1B, 5B, 7B, 5A, and 6A contributing to additional variation of 3.2–11.7%. Twenty-five lines with the best allele combinations of SNPs associated with resistance to both races have been identified as candidates for future variety release and breeding. The results of the present study will be further validated in other independent genetic backgrounds to be able to use markers in breeding.

Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil.

BackgroundPlasmodium vivax is a neglected human malaria parasite that causes significant morbidity in the Americas, the Middle East, Asia, and the Western Pacific. Population genomic approaches remain little explored to map local and regional transmission pathways of P. vivax across the main endemic sites in the Americas, where great progress has been made towards malaria elimination over the past decades.Methodology/Principal findingsWe analyze 38 patient-derived P. vivax genome sequences from Mâncio Lima (ML)–the Amazonian malaria hotspot next to the Brazil-Peru border—and 24 sequences from two other sites in Acre State, Brazil, a country that contributes 23% of malaria cases in the Americas. We show that the P. vivax population of ML is genetically diverse (π = 4.7 × 10−4), with a high polymorphism particularly in genes encoding proteins putatively involved in red blood cell invasion. Paradoxically, however, parasites display strong genome-wide linkage disequilibrium, being fragmented into discrete lineages that are remarkably stable across time and space, with only occasional recombination between them. Using identity-by-descent approaches, we identified a large cluster of closely related sequences that comprises 16 of 38 genomes sampled in ML over 26 months. Importantly, we found significant ancestry sharing between parasites at a large geographic distance, consistent with substantial gene flow between regional P. vivax populations.Conclusions/SignificanceWe have characterized the sustained expansion of highly inbred P. vivax lineages in a malaria hotspot that can seed regional transmission. Potential source populations in hotspots represent a priority target for malaria elimination in the Amazon.

Molecular diversity and genetic structure of wild rice accessions (Oryza rufipogon Griff.) in Guangdong Province, China, as revealed by SNP markers

Common wild rice (Oryza rufipogon Griff.) is a valuable germplasm resource for rice genetic improvement. However, it is endangered due to habitat loss. Assessment of its genetic diversity is essential for germplasm characterization, utilization, and conservation. In the present study, we developed a 79,422-single-nucleotide polymorphism genome-wide array by specific-locus amplified fragment sequencing (SLAF-seq) and used it to genotype 110 common wild rice accessions from different agroclimatic zones in Guangdong Province, China. The number of SNPs on each chromosome was consistent with the physical length of the chromosome, and the average marker density was approximately 4.83 kb/SNP. The genome-wide linkage disequilibrium decayed to its half maximum within 0.52 kb for O. rufipogon. The common wild rice accessions were clustered into four subgroups, with some overlap in principal component analysis and a neighbor-joining tree that was also confirmed by model-based admixture. Analysis of molecular variance revealed significant genetic differentiation among individuals, accounting for 90.64% of the total molecular variance; however, the differences between the subgroups (Enping, Lufeng, Suixi and Fogang) were not significant. The significant genetic variations in resources collected in different years were observed in Enping and Suixi subgroups, but there was no significant genetic variation in Fogang subgroup. The pairwise fixation index (FST) suggested that the populations from Enping and Lufeng were most differentiated while those from Enping and Fogang were least differentiated. Our results illustrated clear genetic relationships among 110 common wild rice accessions, showing a potential molecular marker-based strategy for preserving the genetic diversity and accelerating the systematic utilization of these important wild rice resources.

Molecular Markers Associated with Agro-Physiological Traits under Terminal Drought Conditions in Bread Wheat.

Terminal drought stress poses a big challenge to sustain wheat grain production in rain-fed environments. This study aimed to utilize the genetically diverse pre-breeding lines for identification of genomic regions associated with agro-physiological traits at terminal stage drought stress in wheat. A total of 339 pre-breeding lines panel derived from three-way crosses of ‘exotics × elite × elite’ lines were evaluated in field conditions at Obregon, Mexico for two years under well irrigated as well as drought stress environments. Drought stress was imposed at flowering by skipping the irrigations at pre and post anthesis stage. Results revealed that drought significantly reduced grain yield (Y), spike length (SL), number of grains spikes−1 (NGS) and thousand kernel weight (TKW), while kernel abortion (KA) was increased. Population structure analysis in this panel uncovered three sub-populations. Genome wide linkage disequilibrium (LD) decay was observed at 2.5 centimorgan (cM). The haplotypes-based genome wide association study (GWAS) identified significant associations of Y, SL, and TKW on three chromosomes; 4A (HB10.7), 2D (HB6.10) and 3B (HB8.12), respectively. Likewise, associations on chromosomes 6B (HB17.1) and 3A (HB7.11) were found for NGS while on chromosome 3A (HB7.12) for KA. The genomic analysis information generated in the study can be efficiently utilized to improve Y and/or related parameters under terminal stage drought stress through marker-assisted breeding.

Genome-wide linkage disequilibrium in South African abalone, Haliotis midae, and implications for understanding complex traits

The abalone, Haliotis midae, is the most import aquaculture species in South Africa. Nonetheless, genetic improvement of the species is hindered by late maturity and long generational intervals. The application of marker assisted- and genomic selection is one strategy to accelerate genetic gains for production traits. However, with limited understanding of linkage disequilibrium, it is unclear how such a molecular-driven selective breeding approach should be implemented. The aim of this study was therefore to characterise LD in two wild and two cultured populations of the South African abalone and to investigate how their respective population histories have shaped LD across the genome and on a locus-specific scale. One hundred and twelve mapped microsatellite markers were used to genotype the four population cohorts. Levels of genetic diversity were similarly high for all cohorts, although significant population structure was evident. The decay of LD with genetic distance was present in both cultured cohorts, but almost absent within the wild cohorts. Levels of LD were high in comparison with other aquaculture species, but could be explained by various population demographic events related to the over-exploitation of the species in the wild, domestication, and the reproductive strategy of abalone. A number of interesting patterns were observed surrounding candidate loci under selection, including possible hitch-hiking events and epistatic interactions between loci. In terms of the potential for association studies, the persistence of significant LD over relatively long distances within these cohorts suggested that fewer markers would be sufficient to ensure marker-QTL co-segregation, although the use of generationally discrete cohorts is advised.

GWAS to Identify Genetic Loci for Resistance to Yellow Rust in Wheat Pre-Breeding Lines Derived From Diverse Exotic Crosses.

Yellow rust (YR) or stripe rust, caused by Puccinia striformis f. sp tritici Eriks (Pst), is a major challenge to resistance breeding in wheat. A genome wide association study (GWAS) was performed using 22,415 single nucleotide polymorphism (SNP) markers and 591 haplotypes to identify genomic regions associated with resistance to YR in a subset panel of 419 pre-breeding lines (PBLs) developed at International Center for Maize and Wheat Improvement (CIMMYT). The 419 PBLs were derived from an initial set of 984 PBLs generated by a three-way crossing scheme (exotic/elite1//elite2) among 25 best elites and 244 exotics (synthetics, landraces) from CIMMYT’s germplasm bank. For the study, 419 PBLs were characterized with 22,415 high-quality DArTseq-SNPs and phenotyped for severity of YR disease at five locations in Mexico. A population structure was evident in the panel with three distinct subpopulations, and a genome-wide linkage disequilibrium (LD) decay of 2.5 cM was obtained. Across all five locations, 14 SNPs and 7 haplotype blocks were significantly (P < 0.001) associated with the disease severity explaining 6.0 to 14.1% and 7.9 to 19.9% of variation, respectively. Based on average LD decay of 2.5 cM, identified 14 SNP–trait associations were delimited to seven quantitative trait loci in total. Seven SNPs were part of the two haplotype blocks on chromosome 2A identified in haplotypes-based GWAS. In silico analysis of the identified SNPs showed hits with interesting candidate genes, which are related to pathogenic process or known to regulate induction of genes related to pathogenesis such as those coding for glunolactone oxidase, quinate O-hydroxycinnamoyl transferase, or two-component histidine kinase. The two-component histidine kinase, for example, acts as a sensor in the perception of phytohormones ethylene and cytokinin. Ethylene plays a very important role in regulation of multiple metabolic processes of plants, including induction of defense mechanisms mediated by jasmonate. The SNPs linked to the promising genes identified in the study can be used for marker-assisted selection.

P1542Identification of 13 novel susceptibility loci for early-onset myocardial infarction, hypertension, or chronic kidney disease in Japanese

Abstract Background Early-onset cardiovascular and renal diseases have a strong genetic component. Purpose To identify genetic variants that confer susceptibility to early-onset myocardial infarction (MI), hypertension, or chronic kidney disease (CKD) in Japanese. We have performed exome-wide association studies (EWASs) in subjects with early-onset forms of these diseases. Methods A total of 8093 individuals aged ≤65 years was enrolled in the study. The EWASs for MI, hypertension, and CKD were conducted with 6926 subjects (1152 cases, 5774 controls), 8080 subjects (3444 cases, 4636 controls), and 2556 subjects (1051 cases, 1505 controls), respectively. Genotyping of single nucleotide polymorphisms (SNPs) was performed with Illumina Human Exome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays. The relation of allele frequencies for 31,245, 31,276, or 31,514 SNPs that passed quality control to MI, hypertension, and CKD, respectively, was examined with Fisher's exact test. Bonferroni's correction for statistical significance of association was applied to compensate for multiple comparisons of genotypes with MI, hypertension, or CKD. Results The EWASs of allele frequencies revealed that 25, 11, and 11 SNPs were significantly associated with MI (P&lt;1.60 × 10–6), hypertension (P&lt;1.60 × 10–6), or CKD (P&lt;1.59 × 10–6), respectively. Multivariable logistic regression analysis with adjustment for covariates showed that all 25, 11, and 11 SNPs were significantly related to MI (P&lt;0.0005), hypertension (P&lt;0.0011), or CKD (P&lt;0.0011), respectively. After examination of results from previous genome-wide association studies and linkage disequilibrium of the identified SNPs, we newly identified 11 loci (TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP-GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, GOSR2, ZNF77), six loci (MOB3C-TMOD4, COL6A3, COL6A5, CXCL8-MARCH1, NFKBIL1-6p21.3-NCR3, PLCB2-VPS33B), and seven loci (MOB3C-TMOD4, COL6A3, COL6A5, ADGRL3-CXCL8-MARCH1, MUC17, PLCB2-VPS33B, ZNF77) that were significantly associated with MI, hypertension, or CKD, respectively. Furthermore, six genes (TMOD4, COL6A3, CXCL8, MARCH1, PLCB2, VPS33B) were significantly associated with MI, hypertension, and CKD; two genes (ADGRL3, ZNF77) with MI and CKD; and two genes (COL6A5, MOB3C) with hypertension and CKD. Network analysis showed that the 13, 10, or 11 genes associated with MI, hypertension, or CKD, respectively, in the present study had direct or indirect interactions with the corresponding sets of 50 genes previously shown to be associated with MI, hypertension, or CKD. Conclusion We newly identified 13 loci (MOB3C-TMOD4, COL6A3, ADGRL3-CXCL8-MARCH1, OR52E4, TCHP-GIT2, CCDC63, 12q24.1, OAS3, PLCB2-VPS33B, ZNF77, COL6A5, NFKBIL1-NCR3, MUC17) that confer susceptibility to early-onset MI, hypertension, or CKD. Determination of genotypes for the SNPs at these loci may prove informative for assessment of the genetic risk for MI, hypertension, or CKD in Japanese.

P708Identification of 26 novel loci that confer susceptibility to early-onset coronary artery disease in a Japanese population

Abstract Background Early-onset coronary artery disease (CAD) has a strong genetic component. Although genome-wide association studies have identified various genes and loci significantly associated with CAD mainly in European ancestry populations, genetic variants that contribute to susceptibility to this condition in Japanese individuals remain to be identified definitively. Purpose The purpose of the study was to identify genetic variants that confer susceptibility to early-onset CAD in Japanese. We have now performed exome-wide association studies (EWASs) in subjects with early-onset CAD and controls. Methods A total of 7256 individuals aged ≤65 years was enrolled in the study. The EWAS was conducted with 1482 subjects with CAD and 5774 controls. Genotyping of single nucleotide polymorphisms (SNPs) was performed with Illumina Human Exome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays. The relation of allele frequencies for 31,465 SNPs that passed quality control to CAD was examined with Fisher's exact test. To compensate for multiple comparisons of allele frequencies with CAD, we applied a false discovery rate (FDR) of &lt;0.05 for statistical significance of association. Results The relation of allele frequencies for 31,465 SNPs to CAD with the use of Fisher's exact test showed that 170 SNPs were significantly (FDR &lt;0.05) associated with CAD. Multivariable logistic regression analysis with adjustment for age, sex, and the prevalence of hypertension, diabetes mellitus, and dyslipidemia revealed that 162 SNPs were significantly (P&lt;0.05) related to CAD. A stepwise forward selection procedure was performed to examine the effects of genotypes for the 162 SNPs on CAD. The 54 SNPs were significant (P&lt;0.05) and independent [coefficient of determination (R2), 0.0008 to 0.0297] determinants of CAD. These SNPs together accounted for 15.5% of the cause of CAD. After examination of results from previous genome-wide association studies and linkage disequilibrium of the identified SNPs, we newly identified 21 genes (RNF2, YEATS2, USP45, ITGB8, TNS3, FAM170B-AS1, PRKG1, BTRC, MKI67, STIM1, OR52E4, KIAA1551, MON2, PLUT, LINC00354, TRPM1, ADAT1, KRT27, LIPE, GFY, EIF3L) and five chromosomal regions (2p13, 4q31.2, 5q12, 13q34, 20q13.2) that were significantly associated with CAD. Gene ontology analysis showed that various biological functions were predicted in the 18 genes identified in the present study. The network analysis revealed that the 18 genes had potential direct or indirect interactions with the 30 genes previously shown to be associated with CAD or with the 228 genes identified in previous genome-wide association studies of CAD. Conclusion We have newly identified 26 loci that confer susceptibility to CAD. Determination of genotypes for the SNPs at these loci may prove informative for assessment of the genetic risk for CAD in Japanese.

Coevolution of Sites under Immune Selection Shapes Epstein-Barr Virus Population Structure.

Epstein–Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host’s innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of worldwide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main subpopulations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV.