Rank-based genome-wide analysis reveals the association of Ryanodine receptor-2 gene variants with childhood asthma among human populations

Background

The standard approach to determine unique or shared genetic factors across populations is to identify risk alleles in one population and investigate replication in others. However, since populations differ in DNA sequence information, allele frequencies, effect sizes, and linkage disequilibrium patterns, SNP association using a uniform stringent threshold on p values may not be reproducible across populations. Here, we developed rank-based methods to investigate shared or population-specific loci and pathways for childhood asthma across individuals of diverse ancestry. We performed genome-wide association studies on 859,790 SNPs genotyped in 527 affected offspring trios of European, African, and Hispanic ancestry using publically available asthma database in the Genotypes and Phenotypes database.

Results

Rank-based analyses showed that there are shared genetic factors for asthma across populations, more at the gene and pathway levels than at the SNP level. Although the top 1,000 SNPs were not shared, 11 genes (RYR2, PDE4D, CSMD1, CDH13, ROBO2, RBFOX1, PTPRD, NPAS3, PDE1C, SEMA5A, and CTNNA2) mapped by these SNPs were shared across populations. Ryanodine receptor 2 (RYR2, a statin response-related gene) showed the strongest association in European (p value?=?2.55?×?10^?7) and was replicated in African (2.57?×?10^?4) and Hispanic (1.18?×?10^?3) Americans. Imputation analyses based on the 1000 Genomes Project uncovered additional RYR2 variants associated with asthma. Network and functional ontology analyses revealed that RYR2 is an integral part of dermatological or allergic disorder biological networks, specifically in the functional classes involving inflammatory, eosinophilic, and respiratory diseases.

Conclusion

Our rank-based genome-wide analysis revealed for the first time an association of RYR2 variants with asthma and replicated previously discovered PDE4D asthma gene across human populations. The replication of top-ranked asthma genes across populations suggests that such loci are less likely to be false positives and could indicate true associations. Variants that are associated with asthma across populations could be used to identify individuals who are at high risk for asthma regardless of genetic ancestry.

     

Identification and validation of genomic regions influencing kernel zinc and iron in maize

Key message

Genome-wide association study (GWAS) on 923 maize lines and validation in bi-parental populations identified significant genomic regions for kernel-Zinc and-Iron in maize.

Abstract

Bio-fortification of maize with elevated Zinc (Zn) and Iron (Fe) holds considerable promise for alleviating under-nutrition among the world’s poor. Bio-fortification through molecular breeding could be an economical strategy for developing nutritious maize, and hence in this study, we adopted GWAS to identify markers associated with high kernel-Zn and Fe in maize and subsequently validated marker-trait associations in independent bi-parental populations. For GWAS, we evaluated a diverse maize association mapping panel of 923 inbred lines across three environments and detected trait associations using high-density Single nucleotide polymorphism (SNPs) obtained through genotyping-by-sequencing. Phenotyping trials of the GWAS panel showed high heritability and moderate correlation between kernel-Zn and Fe concentrations. GWAS revealed a total of 46 SNPs (Zn-20 and Fe-26) significantly associated (P?≤?5.03?×?10^?05) with kernel-Zn and Fe concentrations with some of these associated SNPs located within previously reported QTL intervals for these traits. Three double-haploid (DH) populations were developed using lines identified from the panel that were contrasting for these micronutrients. The DH populations were phenotyped at two environments and were used for validating significant SNPs (P?≤?1?×?10^?03) based on single marker QTL analysis. Based on this analysis, 11 (Zn) and 11 (Fe) SNPs were found to have significant effect on the trait variance (P?≤?0.01, R²?≥?0.05) in at least one bi-parental population. These findings are being pursued in the kernel-Zn and Fe breeding program, and could hold great value in functional analysis and possible cloning of high-value genes for these traits in maize.

     

A genome-wide association study on growth traits in orange-spotted grouper (<Emphasis Type="Italic">Epinephelus coioides</Emphasis>) with RAD-seq genotyping

The orange-spotted grouper, Epinephelus coioides, is one of the most popular fish in China and Southeast Asian countries because of its important economic value. However, molecular mechanism underlying the growth of orange-spotted grouper has never been fully understood. Herein, we performed a genome-wide association study (GWAS) on a natural population of 198 individuals aiming to screen the whole genome of orange-spotted grouper for identification of growth-related loci by restrictionsite associated DNA sequencing. In this research, 261,366 single nucleotide polymorphisms (SNPs) were developed, in which 110 SNPs were identified to be correlated with growth and 20 SNPs were further confirmed to be associated with both body weight and total length. From these identified SNPs, we annotated a total of 34 genes, including adgrb2, csnkza1, cers5, col22a1, creb5, dnd1, dzank1, dnai1, npy2r, fat3, lrrk2, lrp5, map3k9, and so on. Among these candidate genes, npy2r (neuropeptide Y receptor Y2) was reported to play a critical role in growth of the orange-spotted grouper. In addition, population structure, principal component analysis, kinship matrix and linkage disequilibrium were examined to verify the accuracy and reliability of our GWAS results. Our data will also provide a valuable genetic resource for further marker-assisted selection program to improve growth quality in groupers.     

Identification of genetic loci associated with abdominal visceral adiposity in Korean populations

Abdominal obesity is characterized by accumulation of subcutaneous and visceral fat in the abdomen and has been reported to be largely responsible for many metabolic and vascular diseases. Although substantial effort has been dedicated to identification of genetic factors associated with abdominal obesity, as measured by the waist-hip ratio and waist circumference, only a few studies have explored associations with visceral fat accumulation in the abdomen. Furthermore, genetic studies of abdominal visceral adiposity conducted in Asian ethnic groups are rare. To gain insight into the genetic basis for visceral adiposity in Asian subjects, we conducted genome-wide association analysis for a pool of 1594 Korean subjects. Abdominal visceral fat area was estimated by computed tomography. After adjustment for age, linear association analysis identified three loci showing suggestive evidence of association (P?<?5?×?10^?6) in ASIC2, SLC35F3, and 5q14.2. Stratification by sex revealed one female-specific locus (rs17104731) located near LINC01519 with a genome-wide significant association for visceral adiposity (P?=?4.66?×?10^?8). Since visceral fat has been suggested to influence metabolic traits, we analyzed associations of the loci identified in this study with metabolic indicators, such as glucose, insulin, and lipid levels, and markers of kidney function. A locus (rs6699737) in SLC35F3 showed a nominal association (P?<?5?×?10^?2) with alanine transaminase, aspartate transaminase, and fasting plasma insulin. In addition, the linear association test using genetic risk score demonstrated that visceral adiposity loci detected in this study had a cumulative effect on abdominal visceral fat area, waist-hip-ratio, total cholesterol, and low density lipoprotein cholesterol. In summary, this study reports new loci associated with visceral adiposity and provides evidence supporting involvement of these loci in several metabolic traits in Korean populations.     

Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

Association study of genetic markers of schizophrenia and its cognitive endophenotypes

A replicative analysis of associations of 15 SNPs located in the regions of 11 genes (TCF4, VRK2, NOTCH4, ZNF804A, AGBL1, RELN, ZFP64P1, KCNB2, CSMD1, CPVL, NRIP1) and three intergenic regions (SLCO6A1/LINCOO491, LOC105376248/LOC105376249, SPA17/NRGN) with schizophrenia was conducted in the Russian population of the Siberian region. These SNPs were previously identified in genome-wide association studies (GWAS) of schizophrenia and cognitive abnormalities. The present study confirmed associations of KCNB2 rs2247572, CSMD1 rs2616984, and intergenic rs12807809 located in SPA17/NRGN with schizophrenia. It was established that the frequency of the CSMD1 rs2616984 G/G genotype was higher in patients compared to the control group (OR = 1.73; CI: 1.14–2.62; р = 0.0337). The frequencies of the KCNB2 rs2247572 TT genotype (OR = 0.41; CI: 0.20–0.87; р = 0.0485) and intergenic rs12807809 CT genotype located in SPA17/NRGN (OR = 0.70; CI: 0.53–0.94; р = 0.0464) were significantly decreased in patients compared to the control group.     

Multivariate eQTL mapping uncovers functional variation on the X-chromosome associated with complex disease traits

Very few studies have investigated the associations between genetic polymorphisms and gene expression on the X-chromosome. This is a major bottleneck when conducting functional follow-up studies of trait-associated variants, as those identified in genome-wide association studies (GWAS). We used a multivariate approach to test the association between individual single nucleotide polymorphisms (SNPs) and exon expression levels measured in 356 Epstein–Barr virus-transformed lymphoblastoid cell lines (LCLs) from the Geuvadis RNA sequencing project to identify SNPs associated with variation in gene expression on the X-chromosome, which we refer to as eSNPs. At an FDR of 5 %, we discovered 548 independent [linkage disequilibrium (LD) r ² < 0.1] eSNPs on the X-chromosome. Of these, 35 were in LD (r ² > 0.8) with previously published disease- or trait-associated variants identified through GWAS. One of the strongest eSNPs identified was rs35975601, which was associated with F8A1 expression (p value = 3 × 10^?20) and was in LD with a type 1 diabetes risk variant. Additionally, we identified a number of genes for which eSNPs were in LD with multiple diseases or traits, including DNASE1L1 which was mapped to bilirubin levels, type 1 diabetes and schizophrenia. Our results also indicate that multivariate exon-level analysis provides a more powerful approach than univariate gene-level analysis, particularly when SNPs influence the expression of different exons with different magnitude and/or direction of effect. The associations identified in our study may provide new insights into the molecular process by which gene expression may contribute to trait variation or disease risk in humans.     

Genome-wide association study for soybean cyst nematode resistance in Chinese elite soybean cultivars

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is a highly recalcitrant endoparasite of soybean roots, causing more yield loss than any other pest. To identify quantitative trait loci (QTL) controlling resistance to SCN (HG type 2.5.7, race 1), a genome-wide association study (GWAS) was performed. The association panel, consisting of 120 Chinese soybean cultivars, was genotyped with 7189 single nucleotide polymorphism (SNPs). A total of 6204 SNPs with minor allele frequency >0.05 were used to estimate linkage disequilibrium (LD) and population structure. The mean level of LD measured by r ² declined very rapidly to half its maximum value (0.51) at 220 kb. The overall population structure was approximately coincident with geographic origin. The GWAS results identified 13 SNPs in 7 different genomic regions significantly associated with SCN resistance. Of these, three SNPs were localized in previously mapped QTL intervals, including rhg1 and Rhg4. The GWAS results also detected 10 SNPs in 5 different genomic regions associated with SCN resistance. The identified loci explained an average of 95.5% of the phenotypic variance. The proportion of phenotypic variance was due to additive genetic variance of the validated SNPs. The present study identified multiple new loci and refined chromosomal regions of known loci associated with SCN resistance. The loci and trait-associated SNPs identified in this study can be used for developing soybean cultivars with durable resistance against SCN.     

Genome-wide association study identified genetic variations and candidate genes for plant architecture component traits in Chinese upland cotton

Key message

Thirty significant associations between 22 SNPs and five plant architecture component traits in Chinese upland cotton were identified via GWAS. Four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits. A candidate gene, Gh_D03G0922, might be responsible for plant height in upland cotton.

Abstract

A compact plant architecture is increasingly required for mechanized harvesting processes in China. Therefore, cotton plant architecture is an important trait, and its components, such as plant height, fruit branch length and fruit branch angle, affect the suitability of a cultivar for mechanized harvesting. To determine the genetic basis of cotton plant architecture, a genome-wide association study (GWAS) was performed using a panel composed of 355 accessions and 93,250 single nucleotide polymorphisms (SNPs) identified using the specific-locus amplified fragment sequencing method. Thirty significant associations between 22 SNPs and five plant architecture component traits were identified via GWAS. Most importantly, four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits, and these SNPs were harbored in one linkage disequilibrium block. Furthermore, 21 candidate genes for plant architecture were predicted in a 0.95-Mb region including the four peak SNPs. One of these genes (Gh_D03G0922) was near the significant SNP D03_31584163 (8.40 kb), and its Arabidopsis homologs contain MADS-box domains that might be involved in plant growth and development. qRT-PCR showed that the expression of Gh_D03G0922 was upregulated in the apical buds and young leaves of the short and compact cotton varieties, and virus-induced gene silencing (VIGS) proved that the silenced plants exhibited increased PH. These results indicate that Gh_D03G0922 is likely the candidate gene for PH in cotton. The genetic variations and candidate genes identified in this study lay a foundation for cultivating moderately short and compact varieties in future Chinese cotton-breeding programs.

     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.     

Potential association between <Emphasis Type="Italic">ITPKC</Emphasis> genetic variations and Hirschsprung disease

Hirschsprung disease (HSCR) is a congenital and complex disorder characterized by intestinal obstruction due to the absence of enteric neurons along variable lengths of the hindgut. Our recent genome-wide association study (GWAS) has revealed regional associations with HSCR at several loci of inositol-trisphosphate 3-kinase C (ITPKC). For fine mapping, we additionally selected and genotyped a total of 12 single nucleotide polymorphisms (SNPs) of ITPKC in 187 HSCR patients and 283 unaffected controls, and performed a further combined imputation analysis based on genotype data from this second stage of fine mapping and our previous GWAS stage, totaling 902 subjects (187 HSCR cases and 715 controls). As a result, several SNPs (minimum P?=?0.004) and a haplotype (P?=?0.02) were found to be significantly associated with HSCR. In further in silico analyses to ascertain the potential functions of the significant variants, the change from the common allele to the rare allele of the highly conserved nonsynonymous rs76785336 showed a difference in mRNA folding structure. In the case of intronic SNPs, rs2607420 with a high consensus value was predicted to be a new splice site. Although this study has limitations (such as lack of functional evaluations, small number of cases, and further need of replication in other cohorts), our findings suggest that genetic variants of ITPKC may have a potential association with HSCR susceptibility and/or developmental diseases related to enteric nervous system development.     

Mouse genome-wide association study identifies polymorphisms on chromosomes 4, 11, and 15 for age-related cardiac fibrosis

Dystrophic cardiac calcinosis (DCC), also called epicardial and myocardial fibrosis and mineralization, has been detected in mice of a number of laboratory inbred strains, most commonly C3H/HeJ and DBA/2J. In previous mouse breeding studies between these DCC susceptible and the DCC-resistant strain C57BL/6J, 4 genetic loci harboring genes involved in DCC inheritance were identified and subsequently termed Dyscalc loci 1 through 4. Here, we report susceptibility to cardiac fibrosis, a sub-phenotype of DCC, at 12 and 20 months of age and close to natural death in a survey of 28 inbred mouse strains. Eight strains showed cardiac fibrosis with highest frequency and severity in the moribund mice. Using genotype and phenotype information of the 28 investigated strains, we performed genome-wide association studies (GWAS) and identified the most significant associations on chromosome (Chr) 15 at 72 million base pairs (Mb) (P < 10^?13) and Chr 4 at 122 Mb (P < 10^?11) and 134 Mb (P < 10^?7). At the Chr 15 locus, Col22a1 and Kcnk9 were identified. Both have been reported to be morphologically and functionally important in the heart muscle. The strongest Chr 4 associations were located approximately 6 Mb away from the Dyscalc 2 quantitative trait locus peak within the boundaries of the Extl1 gene and in close proximity to the Trim63 and Cap1 genes. In addition, a single-nucleotide polymorphism association was found on chromosome 11. This study provides evidence for more than the previously reported 4 genetic loci determining cardiac fibrosis and DCC. The study also highlights the power of GWAS in the mouse for dissecting complex genetic traits.     

Genetic variants influencing lipid levels and risk of dyslipidemia in Chinese population

Recently, several human genetic and genomewide association studies (GWAS) have discovered many genetic loci that are associated with the concentration of the blood lipids. To confirm the reported loci in Chinese population, we conducted a cross-section study to analyse the association of 25 reported SNPs, genotyped by the ABI SNaPshot method, with the blood levels of total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TG) in 1900 individuals by multivariate analysis. Logistic regression was applied to assess the association of the genetic loci with the risk of different types of dyslipidemia. Our study has convincingly identified that 12 of 25 studied SNPs were strongly associated with one or more blood lipid parameters (TC, LDL, HDL and TG). Among the 12 associated SNPs, 10 significantly influence the risk of one or more types of dyslipidemia. We firstly found four SNPs (rs12654264 in HMGCR; rs2479409 in PCSK9; rs16996148 in CILP2, PBX4; rs4420638 in APOE-C1-C4-C2) robustly and independently associate with four types of dyslipidemia (MHL, mixed hyperlipidemia; IHTC, isolated hypercholesterolemia; ILH, isolated low HDL-C; IHTG, isolated hypertriglyceridemia). Our results suggest that genetic susceptibility is different on the same candidate locus for the different populations. Meanwhile, most of the reported genetic variants strongly influence one or more plasma lipid levels and the risk of dyslipidemia in Chinese population.     

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^?5) and 88 SNPs had suggestive (P < 10^?4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.     

A genome-wide association study identifies genomic loci associated with backfat thickness,carcass weight,and body weight in two commercial pig populations

Growth and fatness traits are economically important in the pig industry. To dissect the genetic architecture of these traits in commercial pigs, we conducted a genome-wide association study (GWAS) for carcass weight, backfat thickness, and body weight in two commercial populations: Duroc × (Landrace × Yorkshire) (DLY) and Duroc populations. To enhance the detection power, three GWAS approaches including single-trait GWAS, multi-trait GWAS and meta-analysis were used in this study. A total of 13 suggestive loci were identified on nine chromosomes. The most significant locus was detected at 272.05 Mb on SSC1, and it was associated with backfat thickness at the first rib in the DLY population. Three genes at the identified loci (TBC1D1, BAAT and PHLPP1) were highlighted as functionally plausible candidate genes for pig growth and fatness traits. Genome-wide significant locus was not evidenced in this study, indicating that large populations are required to identify QTL with minor effects on growth and fatness traits in commercial pig populations, in which intensively artificial selections have been imposed on these traits and small genetical variances usually retain in these traits.     

Genome-wide association study of heading and flowering dates and construction of its prediction equation in Chinese common wheat

Heading date is one of the most important traits in wheat breeding as it affects adaptation and yield potential. A genome-wide association study (GWAS) using the 90 K iSelect SNP genotyping assay indicated that a total of 306 loci were significantly associated with heading and flowering dates in 13 environments in Chinese common wheat from the Yellow and Huai wheat region. Of these, 105 loci were significantly correlated with both heading and flowering dates and were found in clusters on chromosomes 2, 5, 6, and 7. Based on differences in distribution of the vernalization and photoperiod genes among chromosomes, arms, or block regions, 13 novel, environmentally stable genetic loci were associated with heading and flowering dates, including RAC875_c41145_189 on 1DS, RAC875_c50422_299 on 2BL, and RAC875_c48703_148 on 2DS, that accounted for more than 20% phenotypic variance explained (PVE) of the heading/flowering date in at least four environments. GWAS and t test of a combination of SNPs and vernalization and photoperiod alleles indicated that the Vrn-B1, Vrn-D1, and Ppd-D1 genes significantly affect heading and flowering dates in Chinese common wheat. Based on the association of heading and flowering dates with the vernalization and photoperiod alleles at seven loci and three significant SNPs, optimal linear regression equations were established, which show that of the seven loci, the Ppd-D1 gene plays the most important role in modulating heading and flowering dates in Chinese wheat, followed by Vrn-B1 and Vrn-D1. Additionally, three novel genetic loci (RAC875_c41145_189, Excalibur_c60164_137, and RAC875_c50422_299) also show important effect on heading and flowering dates. Therefore, Ppd-D1, Vrn-B1, Vrn-D1, and the novel genetic loci should be further investigated in terms of improving heading and flowering dates in Chinese wheat. Further quantitative analysis of an F₁₀ recombinant inbred lines population identified a major QTL that controls heading and flowering dates within the Ppd-D1 locus with PVEs of 28.4% and 34.0%, respectively; this QTL was also significantly associated with spike length, peduncle length, fertile spikelets number, cold resistance, and tiller number.     

Association mapping utilizing diverse barley lines reveals net form net blotch seedling resistance/susceptibility loci

Key message

A diverse collection of barley lines was phenotyped with three North American Pyrenophora teres f. teres isolates and association analyses detected 78 significant marker-trait associations at 16 genomic loci.

Abstract

Pyrenophora teres f. teres is a necrotrophic fungal pathogen and the causal agent of the economically important foliar disease net form net blotch (NFNB) of barley. The deployment of effective and durable resistance against P. teres f. teres has been hindered by the complexity of quantitative resistance and susceptibility. Several bi-parental mapping populations have been used to identify QTL associated with NFNB disease on all seven barley chromosomes. Here, we report the first genome-wide association study (GWAS) to detect marker-trait associations for resistance or susceptibility to P. teres f. teres. Geographically diverse barley genotypes from a world barley core collection (957) were genotyped with the Illumina barley iSelect chip and phenotyped with three P. teres f. teres isolates collected in two geographical regions of the USA (15A, 6A and LDNH04Ptt19). The best of nine regression models tested were identified for each isolate and used for association analysis resulting in the identification of 78 significant marker-trait associations (MTA; ?log10p value?>3.0). The MTA identified corresponded to 16 unique genomic loci as determined by analysis of local linkage disequilibrium between markers that did not meet a correlation threshold of R ²?≥?0.1, indicating that the markers represented distinct loci. Five loci identified represent novel QTL and were designated QRptts-3HL, QRptts-4HS, QRptts-5HL.1, QRptts-5HL.2, and QRptts-7HL.1. In addition, 55 of the barley lines examined exhibited a high level of resistance to all three isolates and the SNP markers identified will provide useful genetic resources for barley breeding programs.

     

Genomewide association analysis for awn length linked to the seed shattering gene <Emphasis Type="Italic">qSH1</Emphasis> in rice

Awn is one of the most important domesticated traits in rice (Oryza sativa). Understanding the genetic basis of awn length is important for grain harvest and production, because long awn length is disadvantageous for both grain harvest and milling. We investigated the awn length of 529 rice cultivars and performed a Genomewide association studies (GWAS) in the indica and japonica subpopulations, and the whole population. In total, we found 17 loci associated with awn length. Of these loci, seven were linked to previously reported quantitative trait loci, and one was linked to the awn gene An-1. Nine novel loci were repeatedly identified in different environments. One of the nine associations was identified in both the whole and japonica populations. Special interest was the detection of the most significant association SNP, sf0136352825, which was less than 95 kb from the seed shattering gene qSH1. These results may provide potentially favourable haplotypes for molecular breeding in rice.     

Evidence of multiple introductions and genetic admixture of the Asian brush-clawed shore crab <Emphasis Type="Italic">Hemigrapsus takanoi</Emphasis> (Decapoda: Brachyura: Varunidae) along the Northern European coast

The Asian brush-clawed shore crab Hemigrapsus takanoi is a non-indigenous species along the Northern European coast. Although the history of range expansion of European H. takanoi has been well-documented, little is known about the genetic compositions of either the introduced European populations or the native Asian ones. We therefore collected H. takanoi broadly from their native Asian sites and introduced European ranges, and genotyped them by sequencing the mitochondrial 16S RNA gene and by analyzing nuclear microsatellite loci. Our results revealed that the H. takanoi Bay of Seine (France) populations consisted of a genetic admixture between populations in Japan and those in the Yellow Sea region. These French populations should be carefully monitored in the future, since the genetic admixture of multiple source populations may accelerate range expansion in non-indigenous organisms. Our results also suggested that shipping lines from East Asia were more probable vectors than historical juvenile oyster transportations from Japan for the foundation of present European H. takanoi populations. Interestingly, gene flow between populations in Japan and those in the Yellow Sea region (i.e., domestic invasion) was not observed despite the higher potential for artificial translocations via shipping lines in the native Asian range compared with those from Asia to Europe. The lack of domestic invasions implied that intra-specific priority effects of the resident H. takanoi populations played an important role in preventing the successful colonization of artificially-transferred individuals.     

Association of the <Emphasis Type="Italic">rs1870634</Emphasis> Variant in Long Intergenic Non-protein Coding RNA 841 with Coronary Artery Disease: A GWAS-Replication Study in an Iranian Population

Recent genome-wide association studies (GWAS) identified a list of single-nucleotide polymorphisms (SNPs) associated with coronary artery disease (CAD). Replication of GWAS findings in different population corroborated the observed association in the parent GWAS. In this study, we aimed to replicate the association of rs1870634, a GWAS identified SNP, to CAD in an Iranian population. The study population consisted of 267 subjects undergoing coronary angiography coronary angiography including 155 CAD patients and 112 non-CAD age- and gender-matched controls. The genotype determination of rs1870634 SNP performed using high-resolution melting analysis (HRM) technique. Our results revealed that the GG genotype frequency was significantly higher in CAD patients compared with controls (P?=?0.03). The results of binary logistic regression suggested that this genotype was significantly associated with CAD risk adjustment for age, BMI, sex, TC, and LDL-C lipid levels (OR of 2.78, 95% CI (1.10–7.01), P?=?0.03). Moreover, our results showed that the GG+TG genotypes were 2.52 times more likely to develop CAD (95% CI 1.05–6.03) than TT genotype carriers after adjusting for age, sex, and lipid profiles (P?=?0.037). These data showed that the GG genotype could be associated with increased risk of CAD in a sample of Iranian population.