Genome-wide Association Mapping Research Articles

Glial expression of a steroidogenic enzyme underlies natural variation in hitchhiking behavior

Phoresy is an interspecies interaction that facilitates spatial dispersal by attaching to a more mobile species. Hitchhiking species have evolved specific traits for physical contact and successful phoresy, but the regulatory mechanisms involved in such traits and their evolution are largely unexplored. The nematode Caenorhabditis elegans displays a hitchhiking behavior known as nictation during its stress-induced developmental stage. Dauer-specific nictation behavior has an important role in natural C. elegans populations, which experience boom-and-bust population dynamics. In this study, we investigated the nictation behavior of 137 wild C. elegans strains sampled throughout the world. We identified species-wide natural variation in nictation and performed a genome-wide association mapping. We show that the variants in the promoter of nta-1, encoding a putative steroidogenic enzyme, underlie differences in nictation. This difference is due to the changes in nta-1 expression in glial cells, which implies that glial steroid metabolism regulates phoretic behavior. Population genetic analysis and geographic distribution patterns suggest that balancing selection maintained two nta-1 haplotypes that existed in ancestral C. elegans populations. Our findings contribute to further understanding of the molecular mechanism of species interaction and the maintenance of genetic diversity within natural populations.

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

IntroductionEffective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms.MethodsWe investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables.ResultsSnap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Determination of genomic regions associated with early storage root formation and bulking in cassava.

Early cassava storage root formation and bulking is a medium of escape that farmers and processors tend to adopt in cases of abiotic and biotic stresses like drought, flood, and destruction by domestic animals. In this study, 220 cassava genotypes from the International Institute of Tropical Agriculture (IITA), National Root Crops Research Institute (NRCRI), International Center for Tropical Agriculture (CIAT), local farmers (from farmer's field), and NextGen project were evaluated in three locations (Umudike, Benue, and Ikenne). The trials were laid out using a split plot in a randomized incomplete block design (alpha lattice) with two replications in 2 years. The storage roots for each plant genotype were sampled or harvested at 3, 6, 9, and 12 month after planting (MAP). All data collected were analyzed using the R-statistical package. The result showed moderate to high heritability among the traits, and there were significant differences (p< 0.05) among the performances of the genotypes. The genome-wide association mapping using the BLINK model detected 45 single-nucleotide polymorphism (SNP) markers significantly associated with the four early storage root bulking and formation traits on Chromosomes 1, 2, 3, 4, 5, 6, 8, 9, 10, 13, 14, 17, and 18.A total of 199 putative candidate genes were found to be directly linked to early storage root bulking and formation. The functions of these candidate genes were further characterized to regulate i) phytohormone biosynthesis, ii) cellular growth and development, and iii) biosynthesis of secondary metabolites for accumulation of starch and defense. Genome-wide association study (GWAS) also revealed the presence of four pleiotropic SNPs, which control starch content, dry matter content, dry yield, and bulking and formation index. The information on the GWAS could be used to develop improved cassava cultivars by breeders. Five genotypes (W940006, NR090146, TMS982123, TMS13F1060P0014, and NR010161) were selected as the best early storage root bulking and formation genotypes across the plant age. These selected cultivars should be used as sources of early storage root bulking and formation in future breeding programs.

Reshaping the Primary Cell Wall: Dual Effects on Plant Resistance to Ralstonia solanacearum and Heat Stress Response.

Temperature elevation drastically affects plant defense responses to Ralstonia solanacearum and inhibits the major source of resistance in Arabidopsis thaliana, mediated by the receptor pair RRS1-R/RPS4. In this study, we refined a previous Genome-Wide Association (GWA) mapping analysis by using a local score approach and detected the primary cell wall CESA3 gene as a major gene involved in plant response to R. solanacearum at both 27°C and elevated temperature, 30°C. We functionally validated CESA3 as a susceptibility gene involved in resistance to R. solanacearum at both 27°C and 30°C through a reverse genetic approach. We provide evidence that the cesa3mre1 mutant enhances resistance to bacterial disease and that resistance is associated with an alteration of root cell morphology conserved at elevated temperature. However, even by forcing the entry of the bacterium to bypass the primary cell wall barrier, the cesa3mre1 mutant still showed enhanced resistance to R. solanacearum with delayed onset of bacterial wilt symptoms. We demonstrated that the cesa3mre1 mutant had constitutive expression of the defense-related gene VSP1 which is up-regulated at elevated temperature and that during infection its expression level is maintained higher than in the wild-type Col-0. In conclusion, this study reveals that alteration of the primary cell wall by mutating the cellulose synthase subunit CESA3 contributes to enhanced resistance to R. solanacearum, remaining effective under heat stress. We expect that these results will help to identify robust genetic sources of resistance to R. solanacearum in the context of global warming.

Natural genetic variation underlying the negative effect of elevated CO2 on ionome composition in Arabidopsis thaliana.

The elevation of atmospheric CO2 leads to a decline in plant mineral content, which might pose a significant threat to food security in coming decades. Although few genes have been identified for the negative effect of elevated CO2 on plant mineral composition, several studies suggest the existence of genetic factors. Here, we performed a large-scale study to explore genetic diversity of plant ionome responses to elevated CO2, using six hundred Arabidopsis thaliana accessions, representing geographical distributions ranging from worldwide to regional and local environments. We show that growth under elevated CO2 leads to a global decrease of ionome content, whatever the geographic distribution of the population. We observed a high range of genetic diversity, ranging from the most negative effect to resilience or even to a benefit in response to elevated CO2. Using genome-wide association mapping, we identified a large set of genes associated with this response, and we demonstrated that the function of one of these genes is involved in the negative effect of elevated CO2 on plant mineral composition. This resource will contribute to understand the mechanisms underlying the effect of elevated CO2 on plant mineral nutrition, and could help towards the development of crops adapted to a high-CO2 world.

Natural genetic variation underlying the negative effect of elevated CO2 on ionome composition in Arabidopsis thaliana

The elevation of atmospheric CO2 leads to a decline in plant mineral content, which might pose a significant threat to food security in coming decades. Although few genes have been identified for the negative effect of elevated CO2 on plant mineral composition, several studies suggest the existence of genetic factors. Here, we performed a large-scale study to explore genetic diversity of plant ionome responses to elevated CO2, using six hundred Arabidopsis thaliana accessions, representing geographical distributions ranging from worldwide to regional and local environments. We show that growth under elevated CO2 leads to a global decrease of ionome content, whatever the geographic distribution of the population. We observed a high range of genetic diversity, ranging from the most negative effect to resilience or even to a benefit in response to elevated CO2. Using genome-wide association mapping, we identified a large set of genes associated with this response, and we demonstrated that the function of one of these genes is involved in the negative effect of elevated CO2 on plant mineral composition. This resource will contribute to understand the mechanisms underlying the effect of elevated CO2 on plant mineral nutrition, and could help towards the development of crops adapted to a high-CO2 world.

CgSSR marker based trait associations for Sitotroga cerealella (Olivier) resistance in stored paddy rice

Rice is stored for a longer time and many primary and secondary insect pests attack grains, causing both quantitative and qualitative losses. Physicochemical characteristics of grain influence stored grain insect susceptibility and resistance. Identification of gene/Quantitative Trait Loci (QTLs) and introgression in the cultivars would provide the desired solution for varietal resistance against S. cerealella. A group of 80 rice varieties was used for genetic analysis for the identification of candidate genes responsible for Sitotroga cerealella resistance. Heat maps derived from values in the kinship matrix were used to group the varieties into three subpopulations. The largest subpopulation has 34 varieties, while the lowest had 19. The first three principal components (PCs) explained 16, 9, and 6 percent of total variation, respectively, as indicated by the Principal Component Analysis (PCA). The PCA results complemented the heat map results in defining population structure. Genome-wide association mapping can detect marker-trait associations/QTLs in a heterogeneous population of breeding lines, despite population structure and unequal family relatedness. As shown by the Genomic Association and Prediction Integrated Tool (GAPIT) heat map, there are three subpopulations in the association mapping panel. Two cgSSR markers, related to gene SDG725 and FLO2, were linked to the attributes and responsible for conferring resistance against S. cerealella infestation. This is the first report of association and that would help in the development of cultivar with S. cerealella resistance.

Evolutionary patterns of DNA base composition at polymorphic sites highlight the role of the environment in shaping barley and rice genomes.

Insights into changes in genome base composition underlying crop domestication can be gained by using comparative genomics. With this approach, previous studies have reported that crop genomes during domestication accumulate more nucleotides adenine (A) and thymine (T) (termed as [AT]-increase) across polymorphic sites. However, the potential influence of the environment or its factors, for example, solar ultraviolet (UV) radiation and temperature, on the [AT]-increase has not been well elucidated. Here, we investigated the [AT]-increase in barley (Hordeum vulgare L.) and rice (Oryza sativa L.) and the association with natural environments, where accessions are distributed. With 12,798,376 and 2,861,535 single-nucleotide polymorphisms from 368 barley and 1375 rice accessions, respectively, we discovered that [AT] increases from wild accessions to improved cultivars, and genomic regions with larger [AT]-increase tend to have higher UV-related motif frequencies, suggesting solar UV radiation as a potential factor in driving genome variation. To link [AT] change with the geographic distribution, we gathered georeferenced accessions and examined their local environments. Interestingly, negative correlations between [AT] and environmental factors were observed (r=-0.39 ∼ -0.75) and modern accessions with higher [AT] values, as compared with wild relatives, are from the environments with lower solar UV radiation or lower temperature. With [AT] and environmental factors as phenotypes, genome-wide association mapping identified three candidate genes that have the potential to contribute to [AT] variation under the effect of environmental conditions. Our findings provide genomic and environmental insights into evolutionary pattern of DNA base composition and underlying mechanisms.

Genome-Wide Association Study of Rice Diversity Panel Reveals New QTLs for Tolerance to Water Deficit Under the Egyptian Conditions

Drought has a significant impact on rice yield by restricting the crop's ability to grow and develop. Producing rice cultivars adapted to water deficit conditions is still the main interest of rice breeders and geneticists. To address this challenge, a set of 413 highly diverse rice populations were evaluated under normal and water deficit conditions for two growing seasons of 2021 and 2022. High genetic variation was found among genotypes for all studied traits. The heritability estimates ranged from 0.82 (panicle length) to 0.95 (plant height). Sterility percentage (SET%) was the most trait affected by water deficit in two growing seasons. 22 Rice genotypes were classified as drought tolerant in both years. Genome-wide association mapping was performed for all traits in the two growing seasons under both conditions using a total of 700,000 SNPs. The GWAS results revealed important and major SNPs associated with all traits. 26 Significant SNPs with stable allele effects were found to be associated with yield traits under water deficit conditions in both years. The results of this study provided rice genotypes that can be adapted under water deficit conditions and important stable SNP markers that can be used for marker-assisted selection after validation in different genetic backgrounds.

No Common Candidate Genes for Resistance to Fusarium graminearum, F. proliferatum, F. sporotrichioides, and F. subglutanins in Soybean (Glycine max L.) Accessions from Maturity Groups 0 and I: Findings from Genome-Wide Association Mapping.

Seedling diseases and root rot, caused by species of Fusarium, can limit soybean (Glycine max L.) production in the United States. Currently, there are few commercially available cultivars resistant to Fusarium. This study was conducted to assess the resistance of soybean maturity group (MG) accessions from 0 and I to Fusarium proliferatum, F. sporotrichioides, and F. subglutinans, as well as to identify common quantitative trait loci (QTL) for resistance to these pathogens, in addition to F. graminearum, using a genome-wide association study (GWAS). A total of 155, 91, and 48 accessions from the USDA soybean germplasm collection from maturity groups 0 and I were screened with a single isolate each of F. proliferatum, F. sporotrichioides, and F. subglutinans, respectively, using the inoculum layer inoculation method in the greenhouse. The disease severity was assessed 21 days post-inoculation and analyzed using non-parametric statistics to determine the relative treatment effects (RTE). Eleven and seven accessions showed significantly lower RTEs when inoculated with F. proliferatum and F. subglutinans, respectively, compared to the susceptible cultivar 'Williams 82'. One accession was significantly less susceptible to both F. proliferatum and F. subglutinans. The GWAS conducted with 41,985 single-nucleotide markers identified one QTL associated with resistance to both F. proliferatum and F. sporotrichioides, as well as another QTL for resistance to both F. subglutinans and F. graminearum. However, no common QTLs were identified for the four pathogens. The USDA accessions and QTLs identified in this study can be utilized to selectively breed resistance to multiple species of Fusarium.

Leveraging the potential of big genomic and phenotypic data for genome-wide association mapping in wheat

Genome-wide association mapping studies (GWAS) based on Big Data are a potential approach to improve marker-assisted selection in plant breeding. The number of available phenotypic and genomic data sets in which medium-sized populations of several hundred individuals have been studied is rapidly increasing. Combining these data and using them in GWAS could increase both the power of QTL discovery and the accuracy of estimation of underlying genetic effects, but is hindered by data heterogeneity and lack of interoperability. In this study, we used genomic and phenotypic data sets, focusing on Central European winter wheat populations evaluated for heading date. We explored strategies for integrating these data and subsequently the resulting potential for GWAS. Establishing interoperability between data sets was greatly aided by some overlapping genotypes and a linear relationship between the different phenotyping protocols, resulting in high quality integrated phenotypic data. In this context, genomic prediction proved to be a suitable tool to study relevance of interactions between genotypes and experimental series, which was low in our case. Contrary to expectations, fewer associations between markers and traits were found in the larger combined data than in the individual experimental series. However, the predictive power based on the marker-trait associations of the integrated data set was higher across data sets. Therefore, the results show that the integration of medium-sized to Big Data is an approach to increase the power to detect QTL in GWAS. The results encourage further efforts to standardize and share data in the plant breeding community.

Genome-wide association study of plant color in Sorghum bicolor.

Sorghum plant color is the leaf sheath/leaf color and is associated with seed color, tannin and phenol content, head blight disease incidence, and phytoalexin production. In this study, we evaluated plant color of the sorghum mini core collection by scoring leaf sheath/leaf color at maturity as tan, red, or purple across three testing environments and performed genome-wide association mapping (GWAS) with 6,094,317 SNPs markers. Eight loci, one each on chromosomes 1, 2, 4, and 6 and two on chromosomes 5 and 9, were mapped. All loci contained one to three candidate genes. In qPC5-1, Sobic.005G165632 and Sobic.005G165700 were located in the same linkage disequilibrium (LD) block. In qPC6, Sobic.006G149650 and Sobic.006G149700 were located in the different LD block. The single peak in qPC6 covered one gene, Sobic.006G149700, which was a senescence regulator. We found a loose correlation between the degree of linkage and tissue/organ expression of the underlying genes possibly related to the plant color phenotype. Allele analysis indicated that none of the linked SNPs can differentiate between red and purple accessions whereas all linked SNPs can differentiate tan from red/purple accessions. The candidate genes and SNP markers may facilitate the elucidation of plant color development as well as molecular plant breeding.

Genome-wide association mapping in rice MAGIC indica panel detects QTLs and genes for broad-spectrum resistance breeding against African bacterial blight

Rice (Oryza sativa) bacterial blight (BB) is a direct threat to production in Africa. The extent of damage on the continent is not known. We conducted a genome-wide association study (GWAS) to seek out Quantitative trait loci (QTL) that confer the broad-spectrum resistance (BSR) to rice Xanthomonas oryzae pv. oryzae in the rice Multiparent advanced generation intercross (MAGIC) population. The diverse genetic composition and high levels of recombination, are potential resources for the identification of QTL for BSR. Dealing with the ever-evolving Xoo requires identification of novel R gene and resistance quantitative trait loci (QTL). We conducted a genome-wide association study (GWAS) using 250 rice MAGIC indica panel to evaluate for resistance to three Xoo strains from Cote D'Ivoire, namely XooCIkor, XooCIgag and XooCIcii2 which are prevalent highly virulent pathotypes. Using the 3 K RG 4.8mio filtered SNP Dataset, a total of 3 QTL associated with BB resistance on chromosome 9, 6 and 5 were identified. The results of the GWAS identified novel loci that gives useful insight for more investigation and also validated known genes underlying resistance to virulent African strains of BB. Genetic donors identified in this will help plant breeders in improving and broadening resistance to BB.

Genome-wide association studies and transcriptomics reveal mechanisms explaining the diversity of wheat root responses to nutrient availability.

Nutrient availability profoundly influences plant root system architecture, which critically determines crop productivity. While Arabidopsis has provided important insights into the genetic responses to nutrient deficiency, translating this knowledge to crops, particularly wheat, remains a subject of inquiry. Here, examining a diverse wheat population under varying nitrogen (N), phosphorus (P), potassium (K), and iron (Fe) levels, we uncover a spectrum of root responses, spanning from growth inhibition to stimulation, highlighting genotype-specific strategies. Furthermore, we reveal a nuanced interplay between macronutrient deficiency (N, P, and K) and Fe availability, emphasizing the central role of Fe in modulating root architecture. Through genome-wide association mapping, we identify 11 quantitative trait loci underlying root traits under varying nutrient availabilities, including homologous genes previously validated in Arabidopsis, supporting our findings. In addition, utilizing transcriptomics, ROS imaging, and antioxidant treatment, we uncover that wheat root growth inhibition by nutrient deficiency is attributed to ROS accumulation, akin to the role of ROS in governing Arabidopsis root responses to nutrient deficiency. Therefore, our study reveals the conservation of molecular and physiological mechanisms between Arabidopsis and wheat to adjust root growth to nutrient availability, paving the way for targeted crop improvement strategies aimed at increasing nutrient use efficiency.

Genome-wide identification of the key kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L.).

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (G. hirsutum L.). Results showed that 159 kinesin genes, including 15 genes of the kinesin-13 gene subfamily, were identified in upland cotton; of which 157 kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 kinesin genes in upland cotton, including 10 kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 kinesin genes were significantly associated with three fiber traits, among which a kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

Genetic variation within a stick insect species associated with community-level traits.

Phenotypic variation within species can affect the ecological dynamics of populations and communities. Characterizing the genetic variation underlying such effects can help parse the roles of genetic evolution and plasticity in 'eco-evolutionary dynamics' and inform how genetic variation may shape patterns of evolution. Here we employ genome-wide association (GWA) methods in Timema cristinae stick insects and their co-occurring arthropod communities to identify genetic variation associated with community-level traits. Previous studies have shown that maladaptation (i.e., imperfect crypsis) of T. cristinae can reduce the abundance and species richness of other arthropods due to an increase in bird predation. Whether genetic variation that is independent from crypsis has similar effects is unknown and was tested here using genome-wide genotyping-by-sequencing data of stick insects, arthropod community information, and GWA mapping with Bayesian sparse linear mixed models. We find associations between genetic variation in stick insects and arthropod community traits. However, these associations disappeared when host-plant traits are accounted for. We thus use path analysis to disentangle interrelationships among stick-insect genetic variation, host-plant traits and community traits. This revealed that host-plant size has large effects on arthropod communities, while genetic variation in stick insects has a smaller, but still significant effect. Our findings demonstrate that: (1) genetic variation in a species can be associated with community-level traits, but that (2) interrelationships among multiple factors may need to be analyzed to disentangle whether such associations represent causal relationships. This work helps to build a framework for genomic studies of eco-evolutionary dynamics.

Genome-wide association mapping and genomic prediction of stalk rot in two mid-altitude tropical maize populations

Maize stalk rot reduces grain yield and quality. Information about the genetics of resistance to maize stalk rot could help breeders design effective breeding strategies for the trait. Genomic prediction may be a more effective breeding strategy for stalk-rot resistance than marker-assisted selection. We performed a genome-wide association study (GWAS) and genomic prediction of resistance in testcross hybrids of 677 inbred lines from the Tuxpeño and non-Tuxpeño heterotic pools grown in three environments and genotyped with 200,681 single-nucleotide polymorphisms (SNPs). Eighteen SNPs associated with stalk rot shared genomic regions with gene families previously associated with plant biotic and abiotic responses. More favorable SNP haplotypes traced to tropical than to temperate progenitors of the inbred lines. Incorporating genotype-by-environment (G × E) interaction increased genomic prediction accuracy.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic β-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.

Genome-wide association study for resistances to yellow rust, powdery mildew, and Septoria tritici blotch in cultivated emmer

Emmer is a progenitor of bread wheat and evolved in the Levant together with the yellow rust (YR), powdery mildew (PM) fungi, and a precursor of Zymoseptoria tritici causing Septoria tritici blotch (STB). We performed a genome-wide association mapping for the three disease resistances with 143 cultivated emmer accessions in multi-environmental trials. Significant (P < 0.001) genotypic variation was found with high heritabilities for the resistances to the two biotrophs and a moderate heritability for STB resistance. For YR, PM, and STB severity nine, three, and seven marker-trait associations, respectively, were detected that were significant across all environments. Most of them were of low to moderate effect, but for PM resistance a potentially new major gene was found on chromosome 7AS. Genomic prediction abilities were high throughout for all three resistances (≥ 0.8) and decreased only slightly for YR and PM resistances when the prediction was done for the second year with the first year as training set (≥ 0.7). For STB resistance prediction ability was much lower in this scenario (0.4). Despite this, genomic selection should be advantageous given the large number of small QTLs responsible for quantitative disease resistances. A challenge for the future is to combine these multiple disease resistances with better lodging tolerance and higher grain yield.

Genome-wide association mapping for seedling and adult resistance to powdery mildew in barley

Key messageTwo new major QTL were identified for powdery mildew resistance. We confirmed that the QTL on 7HS contributed mainly to the adult-plant resistance, while another one on chromosome arm 1HS made a significant contribution to the seedling resistance.Powdery mildew (PM), caused by Blumeria hordei, can occur at all post emergent stages of barley and constantly threatens crop production. To identify more genes for effective resistance to powdery mildew for use in breeding programs, 696 barley accessions collected from different regions of the world were evaluated for PM resistance at seedling and adult growth stages in three different states of Australia. These barley accessions were genotyped using DArTSeq with over 18,000 markers for a genome-wide association study (GWAS). Using the FarmCPU model, 54 markers showed significant associations with PM resistance scored at the seedling and adult-plant stages in different states of Australia. Another 40 markers showed tentative associations (LOD > 4.0) with resistance. These markers are distributed across all seven barley chromosomes. Most of them were grouped into eleven QTL regions, coinciding with the locations of most of the reported resistance genes. Two major MTAs were identified on chromosome arms 3HS and 5HL, with one on 3HS contributing to adult plant resistance and the one on 5HL to both seedling and adult plant resistance. An MTA on 7HS contributed mainly to the adult-plant resistance, while another one on chromosome arm 1HS made a significant contribution to the seedling resistance.