Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

Background

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F₂ intercross population and a Chinese Sutai half-sib population.

Results

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F₂ animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

Conclusions

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0089-5) contains supplementary material, which is available to authorized users.     

Evaluation of European Schizophrenia GWAS Loci in Asian Populations via Comprehensive Meta-Analyses

Schizophrenia is a severe and highly heritable neuropsychiatric disorder. Recent genetic analyses including genome-wide association studies (GWAS) have implicated multiple genome-wide significant variants for schizophrenia among European populations. However, many of these risk variants were not largely validated in other populations of different ancestry such as Asians. To validate whether these European GWAS significant loci are associated with schizophrenia in Asian populations, we conducted a systematic literature search and meta-analyses on 19 single nucleotide polymorphisms (SNPs) in Asian populations by combining all available case-control and family-based samples, including up to 30,000 individuals. We employed classical fixed (or random) effects inverse variance weighted methods to calculate summary odds ratios (ORs) and 95 % confidence intervals (CIs). Among the 19 GWAS loci, we replicated the risk associations of nine markers (e.g., SNPs at VRK2, ITIH3/4, NDST3, NOTCH4) surpassing significance level (two-tailed P?<?0.05), and three additional SNPs in MIR137 and ZNF804A also showed trend associations (one-tailed P?<?0.05). These risk associations are in the same directions of allelic effects between Asian replication samples and initial European GWAS findings, and the successful replications of these GWAS loci in a different ethnic group provide stronger evidence for their clinical associations with schizophrenia. Further studies, focusing on the molecular mechanisms of these GWAS significant loci, will become increasingly important for understanding of the pathogenesis to schizophrenia.     

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 analysis of multi-environmental spring wheat trials identifies genomic regions for locus-specific trade-offs for grain weight and grain number

Key message

GWAS on multi-environment data identified genomic regions associated with trade-offs for grain weight and grain number.

Abstract

Grain yield (GY) can be dissected into its components thousand grain weight (TGW) and grain number (GN), but little has been achieved in assessing the trade-off between them in spring wheat. In the present study, the Wheat Association Mapping Initiative (WAMI) panel of 287 elite spring bread wheat lines was phenotyped for GY, GN, and TGW in ten environments across different wheat growing regions in Mexico, South Asia, and North Africa. The panel genotyped with the 90 K Illumina Infinitum SNP array resulted in 26,814 SNPs for genome-wide association study (GWAS). Statistical analysis of the multi-environmental data for GY, GN, and TGW observed repeatability estimates of 0.76, 0.62, and 0.95, respectively. GWAS on BLUPs of combined environment analysis identified 38 loci associated with the traits. Among them four loci—6A (85 cM), 5A (98 cM), 3B (99 cM), and 2B (96 cM)—were associated with multiple traits. The study identified two loci that showed positive association between GY and TGW, with allelic substitution effects of 4% (GY) and 1.7% (TGW) for 6A locus and 0.2% (GY) and 7.2% (TGW) for 2B locus. The locus in chromosome 6A (79–85 cM) harbored a gene TaGW2-6A. We also identified that a combination of markers associated with GY, TGW, and GN together explained higher variation for GY (32%), than the markers associated with GY alone (27%). The marker-trait associations from the present study can be used for marker-assisted selection (MAS) and to discover the underlying genes for these traits in spring wheat.

     

Nucleotide polymorphisms of candidate genes of adaptive significance in the ural populations of <Emphasis Type="Italic">Larix sibirica</Emphasis> Ledeb.

The objectives of conservation and sustainable forest management require in depth study of genomes of woody plants and definition of their intraspecific genetic diversity. In recent years, an approach was developed based on the study of “candidate genes” that can potentially be involved in the formation of adaptive traits. In this study, we investigated nucleotide polymorphism of several adaptive candidate genes in the populations of Siberian larch (Larix sibirica Ledeb.) in the Urals. Representatives of this genus are among the most valuable and widely distributed forest tree species in Russia. From ten selected gene loci in the genome of L. sibirica, we isolated and investigated three loci, one of which (ABA-WDS) was sequenced in L. sibirica for the first time. The total length of the analyzed sequence in each individual amounted to 2865 bp. The length of locus alignment was from 360 bp to 1395 bp. In total, we identified 200 polymorphic positions. The most conservative is locus 4CL1-363, and the most polymorphic is locus sSPcDFD040B03103-274. The studied populations of L. sibirica are characterized by a high level of nucleotide polymorphism in comparison with other species and genuses (Picea, Pinus, Pseudotsuga, Abies) conifers plants (Hd = 0.896; π = 0.007; θ_W = 0.015). The most selectively neutral polymorphism (D _T =–0.997) was attributed to locus 4CL1-363, and polymorphism with high probability of adaptability (D _T =–1.807) was determined for the ABA-WDS locus. We identified 54 SNP markers, only five of which were nonsynonymous (9.26%) replacements. The average frequency of SNPs in the three studied loci of L. sibirica was one SNP in 53 bp. We detected unique SNP markers for eight populations, which could potentially be used to identify populations. Populations that are characterized by the highest number of unique SNP markers can be recommended for selection in order to preserve the gene pool of the species.     

Interactions between <Emphasis Type="Italic">Glu-1</Emphasis> and <Emphasis Type="Italic">Glu-3</Emphasis> loci and associations of selected molecular markers with quality traits in winter wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.) DH lines

The quality of wheat depends on a large complex of genes and environmental factors. The objective of this study was to identify quantitative trait loci controlling technological quality traits and their stability across environments, and to assess the impact of interaction between alleles at loci Glu-1 and Glu-3 on grain quality. DH lines were evaluated in field experiments over a period of 4 years, and genotyped using simple sequence repeat markers. Lines were analysed for grain yield (GY), thousand grain weight (TGW), protein content (PC), starch content (SC), wet gluten content (WG), Zeleny sedimentation value (ZS), alveograph parameter W (APW), hectolitre weight (HW), and grain hardness (GH). A number of QTLs for these traits were identified in all chromosome groups. The Glu-D1 locus influenced TGW, PC, SC, WG, ZS, APW, GH, while locus Glu-B1 affected only PC, ZS, and WG. Most important marker-trait associations were found on chromosomes 1D and 5D. Significant effects of interaction between Glu-1 and Glu-3 loci on technological properties were recorded, and in all types of this interaction positive effects of Glu-D1 locus on grain quality were observed, whereas effects of Glu-B1 locus depended on alleles at Glu-3 loci. Effects of Glu-A3 and Glu-D3 loci per se were not significant, while their interaction with alleles present at other loci encoding HMW and LMW were important. These results indicate that selection of wheat genotypes with predicted good bread-making properties should be based on the allelic composition both in Glu-1 and Glu-3 loci, and confirm the predominant effect of Glu-D1d allele on technological properties of wheat grains.     

The effect of QTL-rich region polymorphisms identified by targeted DNA-seq on pig production traits

The aim of the present study was to analyse the effect of PLCD4, PECR, FN1 and PNKD mutations on pig productive traits and tested the usefulness of targeted enrichment DNA sequencing method as tool for preselection of genetic markers. The potential genetic markers for pig productive traits were identified by using targeted enrichment DNA sequencing of chromosome 15 region that is QTL-rich. The selected mutations were genotyped by using HRM, Sanger sequencing and PCR-ACRS methods. The association study was performed by using GLM model in SAS program and included over 500 pigs of 5 populations maintained in Poland. The variation (C/T) of PLCD4 gene affected feed conversion, intramuscular fat and water exudation. The PNKD mutations were associated with texture parameters measured after cooking. In turn, the variation rs792423408 (C/T) in the FN1 gene influenced toughness measured in semimembranosus muscle and growth traits that was observed particularly in Duroc breed. Summarizing, the investigated gene variants delivered valuable information that could be used during developing SNP microarray for genomic estimated breeding value procedure in pigs. Moreover, the study showed that the TEDNA-seq method could be used to preselect the molecular markers associated with pig traits. However, the further association study that included large number animal populations is necessary.     

Structural organization of fatty acid desaturase loci in linseed lines with contrasting linolenic acid contents

Flax (Linum usitatissimum L.), the richest crop source of omega-3 fatty acids (FAs), is a diploid plant with an estimated genome size of ~370 Mb and is well suited for studying genomic organization of agronomically important traits. In this study, 12 bacterial artificial chromosome clones harbouring the six FA desaturase loci sad1, sad2, fad2a, fad2b, fad3a and fad3b from the conventional variety CDC Bethune and the high linolenic acid line M5791 were sequenced, analysed and compared to determine the structural organization of these loci and to gain insights into the genetic mechanisms underlying FA composition in flax. With one gene every 3.2–4.6 kb, the desaturase loci have a higher gene density than the genome’s average of one gene per 7.8–8.2 kb. The gene order and orientation across the two genotypes were generally conserved with the exception of the sad1 locus that was predicted to have additional genes in CDC Bethune. High sequence conservation in both genic and intergenic regions of the sad and fad2b loci contrasted with the significant level of variation of the fad2a and fad3 loci, with SNPs being the most frequently observed mutation type. The fad2a locus had 297 SNPs and 36 indels over ~95 kb contrasting with the fad2b locus that had a mere seven SNPs and four indels in ~110 kb. Annotation of the gene-rich loci revealed other genes of known role in lipid or carbohydrate metabolic/catabolic pathways. The organization of the fad2b locus was particularly complex with seven copies of the fad2b gene in both genotypes. The presence of Gypsy, Copia, MITE, Mutator, hAT and other novel repeat elements at the desaturase loci was similar to that of the whole genome. This structural genomic analysis provided some insights into the genomic organization and composition of the main desaturase loci of linseed and of their complex evolution through both tandem and whole genome duplications.     

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.     

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.     

Genetic structure of the populations of <Emphasis Type="Italic">Dactylorhiza ochroleuca</Emphasis> and <Emphasis Type="Italic">D. incarnata</Emphasis> (Orchidaceae) in the area of their joint growth in Russia and Belarus

We carried out an allozyme analysis to investigate polymorphism and genetic structure of the populations of D. incarnata and D. ochroleuca in regions of their joint growth in Russia and Belarus. We found that D. ochroleuca individuals in the populations of the Urals and Siberia, which are distant fragments from the main range of the species, do not differ significantly from individuals within the main part of the area (Belarus) on the basis of the allelic composition of eight gene loci. We revealed that D. ochroleuca and D. incarnata are differentiated by different alleles of the GDH locus. Thus, we established a genetic marker suitable to distinguish these closely related taxa. In addition to the GDH locus, D. ochroleuca and D. incarnata in the places of their joint growth, differ in the allelic structure of the PGI and NADHD loci. D. incarnata from the Urals and Siberia were polymorphic for both loci, and individuals from Belarus were polymorphic for one locus (PGI). In contrast, all D. ochroleuca individuals growing in sympatric populations with polymorphic D. incarnata were homozygous for the same alleles. Thus, comparison of the genetic structure of D. ochroleuca and D. incarnata points to the existence of a genetic isolation and a functioning isolation mechanism even under conditions of their joint growth. We found that the GDH locus in D. incarnata is polymorphic only in populations which grow together with D. ochroleuca, with exception a few examples. Thus, we conclude that variability of the GDH locus in D. incarnata is associated with hybridization with D. ochroleuca.     

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.     

Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc?×?(Landrace × Yorkshire) commercial population

Background

Understanding the genetic mechanisms that underlie meat quality traits is essential to improve pork quality. To date, most quantitative trait loci (QTL) analyses have been performed on F₂ crosses between outbred pig strains and have led to the identification of numerous QTL. However, because linkage disequilibrium is high in such crosses, QTL mapping precision is unsatisfactory and only a few QTL have been found to segregate within outbred strains, which limits their use to improve animal performance. To detect QTL in outbred pig populations of Chinese and Western origins, we performed genome-wide association studies (GWAS) for meat quality traits in Chinese purebred Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) (DLY) commercial population.

Methods

Three hundred and thirty six Chinese Erhualian and 610 DLY pigs were genotyped using the Illumina PorcineSNP60K Beadchip and evaluated for 20 meat quality traits. After quality control, 35 985 and 56 216 single nucleotide polymorphisms (SNPs) were available for the Chinese Erhualian and DLY datasets, respectively, and were used to perform two separate GWAS. We also performed a meta-analysis that combined P-values and effects of 29 516 SNPs that were common to Erhualian, DLY, F₂ and Sutai pig populations.

Results

We detected 28 and nine suggestive SNPs that surpassed the significance level for meat quality in Erhualian and DLY pigs, respectively. Among these SNPs, ss131261254 on pig chromosome 4 (SSC4) was the most significant (P = 7.97E-09) and was associated with drip loss in Erhualian pigs. Our results suggested that at least two QTL on SSC12 and on SSC15 may have pleiotropic effects on several related traits. All the QTL that were detected by GWAS were population-specific, including 12 novel regions. However, the meta-analysis revealed seven novel QTL for meat characteristics, which suggests the existence of common underlying variants that may differ in frequency across populations. These QTL regions contain several relevant candidate genes.

Conclusions

These findings provide valuable insights into the molecular basis of convergent evolution of meat quality traits in Chinese and Western breeds that show divergent phenotypes. They may contribute to genetic improvement of purebreds for crossbred performance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0120-x) contains supplementary material, which is available to authorized users.     

De novo development and characterization of polymorphic microsatellite markers in a schilbid catfish, <Emphasis Type="Italic">Silonia silondia</Emphasis> (Hamilton, 1822) and their validation for population genetic studies

The stock characterization of wild populations of Silonia silondia is important for its scientific management. At present, the information on genetic parameters of S. silondia is very limited. The species-specific microsatellite markers were developed in current study. The validated markers were used to genotype individuals from four distant rivers. To develop de novo microsatellite loci, an enriched genomic library was constructed for S. silondia using affinity–capture approach. The markers were validated for utility in population genetics. A total number of 76 individuals from four natural riverine populations were used to generate data for population analysis. The screening of isolated repeat sequences yielded eleven novel polymorphic microsatellite loci. The microsatellite loci exhibited high level of polymorphism, with 6–24 alleles per locus and the PIC value ranged from 0.604 to 0.927. The observed (Ho) and expected (He) heterozygosities ranged from 0.081 to 0.84 and 0.66 to 0.938, respectively. The AMOVA analysis indicated significant genetic differentiation among riverine populations (overall F_ST = 0.075; P < 0.0001) with maximum variation (92.5 %) within populations. Cross-priming assessment revealed successful amplification (35–38 %) of heterologous loci in four related species viz. Clupisoma garua, C. taakree, Ailia coila and Eutropiichthys vacha. The results demonstrated that these de novo polymorphic microsatellite loci are promising for population genetic variation and diversity studies in S. silondia. Cross-priming results indicated that these primers can help to get polymorphic microsatellite loci in the related catfish species of family Schilbidae.     

Genome-wide association study for growth and fatness traits in Chinese Sujiang pigs

Growth and fatness traits are complex and economically important traits in the pig industry. The molecular basis underlying porcine growth and fatness traits remains largely unknown. To uncover genetic loci and candidate genes for these traits, we explored the GeneSeek GGP Porcine 80K SNP chip to perform a GWAS for seven growth and fatness traits in 365 individuals from the Sujiang pig, a recently developed breed in China. We identified two, 17, one and 11 SNPs surpassing the suggestively significant threshold (P < 1.86 × 10⁻⁵) for body weight, chest circumference, chest width and backfat thickness respectively. Of these SNPs, 20 represent novel genetic loci, and five and four SNPs were respectively associated with chest circumference and backfat thickness at a genome-wide significant threshold (P < 9.31 × 10⁻⁷). Eight SNPs had a pleiotropic effect on both chest circumference and backfat thickness. The most remarkable locus resided in a region between 72.95 and 76.27 Mb on pig chromosome 4, harboring a number of previously reported quantitative trait loci related to backfat deposition. In addition to two reported genes (PLAG1 and TAS2R38), we identified four genes including GABRB3, ZNF106, XKR4 and MGAM as novel candidates for body weight and backfat thickness at the mapped loci. Our findings provide insights into the genetic architecture of porcine growth and fatness traits and potential markers for selective breeding of Chinese Sujiang pigs.     

Genome-wide association and validation of key loci for yield-related traits in wheat founder parent Xiaoyan 6

Xiaoyan 6, one of the most important founder parents in wheat, possesses many superior agronomic traits and has played a crucial role in Chinese wheat breeding programs. In this study, a panel of 66 elite wheat accessions derived from Xiaoyan 6 was planted in four growing seasons; genome-wide association study (GWAS) was performed for six yield-related traits using the wheat 90K genotyping assay. A total of 803 significant marker-trait associations (MTAs) that explained up to 35.0% of the phenotypic variation were detected. Of these, the locus QTkw-5B which contains 19 MTAs for thousand kernel weight (TKW) was consistently detected in three growing seasons and confirmed in a recombinant inbred line (RIL) population by developing simple sequence repeats (SSR) and kompetitive allele-specific PCR (KASP) markers. The locus QPh-3A containing eight repetitive MTAs for plant height (PH) was consistently identified in all the four growing seasons and validated in a RIL population by developing SSR markers. The transmission of Xiaoyan 6 allele indicated that the favorite allele of QPh-3A was strongly selected in breeding programs. Comparing with previous studies, QTkw-5B and QPh-3A should be novel QTL. The locus QFss-2D for fertile spikelet number per spike (FSS) was identified and then validated in three bi-parental populations. This locus controlled various spike-related traits and may be a key spike polymorphic locus. This study could provide insight into dissecting yield-related traits in the breeding population and reliable molecular markers that might be valuable for marker-assisted selection in wheat high-yield breeding programs.     

Genome-wide association mapping for seedling and field resistance to <Emphasis Type="Italic">Puccinia striiformis</Emphasis> f. sp. <Emphasis Type="Italic">tritici</Emphasis> in elite durum wheat

Key message

Genome-wide association analysis in tetraploid wheat revealed novel and diverse loci for seedling and field resistance to stripe rust in elite spring durum wheat accessions from worldwide.

Abstract

Improving resistance to stripe rust, caused by Puccinia striiformis f. sp. tritici, is a major objective for wheat breeding. To identify effective stripe rust resistance loci, a genome-wide association study (GWAS) was conducted using 232 elite durum wheat (Triticum turgidum ssp. durum) lines from worldwide breeding programs. Genotyping with the 90 K iSelect wheat single nucleotide polymorphism (SNP) array resulted in 11,635 markers distributed across the genome. Response to stripe rust infection at the seedling stage revealed resistant and susceptible accessions present in rather balanced frequencies for the six tested races, with a higher frequency of susceptible responses to United States races as compared to Italian races (61.1 vs. 43.1% of susceptible accessions). Resistance at the seedling stage only partially explained adult plant resistance, which was found to be more frequent with 67.7% of accessions resistant across six nurseries in the United States. GWAS identified 82 loci associated with seedling stripe rust resistance, five of which were significant at the false discovery rate adjusted P value <0.1 and 11 loci were detected for the field response at the adult plant stages in at least two environments. Notably, Yrdurum-1BS.1 showed the largest effect for both seedling and field resistance, and is therefore considered as a major locus for resistance in tetraploid wheat. Our GWAS study is the first of its kind for stripe rust resistance in tetraploid wheat and provides an overview of resistance in elite germplasm and reports new loci that can be used in breeding resistant cultivars.

     

Genomewide association study for economic traits in the large yellow croaker with different numbers of extreme phenotypes

A traditional genomewide association study (GWAS) detects genotype–phenotype associations by the vast number of genotyped individuals. This method requires large-scale samples and considerable sequencing costs. Extreme phenotypic sampling proposes make GWAS more cost-efficient and are applied more widely. With extreme phenotypic sampling, we performed a GWAS for n-3 highly unsaturated fatty acids (HUFA) and eviscerated weight (EW) traits in the large yellow croaker population. Of the 32,249 and 29,748 detected SNPs for the two traits, three candidate regions were found in each trait. Three candidate regions associated with HUFA were known near genes on chromosomes 4 and 11, and three candidate regions were on chromosome 6, and 15 for the EW trait. By combing through our GWAS results and the biological functional analysis of the genes, we suggest that the FABP, DGAT, ATP8B1, FAF2 and CERS2 genes,  as well as the IGF2, BORA, CYP1A1, GRTP1 and HOX genes are promising candidate genes for n-3 HUFA and EW, respectively, in the large yellow croaker. Moreover, compared with the different numbers of the extreme phenotypic sampling, we conclude that 60% of the extreme phenotypic subsample can obtain a similar result as GWAS with whole phenotypes. Thus, extreme phenotypic sampling could save 40% of the cost for genotyping and DNA extraction without loss of the candidate regions and functional genes. Our study may provide a basis for further genomic breeding and a reference for others who want to perform GWAS with extreme phenotypes.