Genome-wide association study dissects the genetic control of plant height and branch number in response to low-phosphorus stress in Brassica napus

Background and AimsOilseed rape (Brassica napus) is one of the most important oil crops worldwide. Phosphorus (P) deficiency severely decreases the plant height and branch number of B. napus. However, the genetic bases controlling plant height and branch number in B. napus under P deficiency remain largely unknown. This study aims to mine candidate genes for plant height and branch number by genome-wide association study (GWAS) and determine low-P-tolerance haplotypes.MethodsAn association panel of B. napus was grown in the field with a low P supply (P, 0 kg ha⁻¹) and a sufficient P supply (P, 40 kg ha⁻¹) across 2 years and plant height and branch number were investigated. More than five million single-nucleotide polymorphisms (SNPs) were used to conduct GWAS of plant height and branch number at two contrasting P supplies.Key ResultsA total of 2127 SNPs were strongly associated (P < 6·25 × 10⁻⁰⁷) with plant height and branch number at two P supplies. There was significant correlation between phenotypic variation and the number of favourable alleles of associated loci on chromosomes A10 (chrA10_821671) and C08 (chrC08_27999846), which will contribute to breeding improvement by aggregating these SNPs. BnaA10g09290D and BnaC08g26640D were identified to be associated with chrA10_821671 and chrC08_27999846, respectively. Candidate gene association analysis and haplotype analysis showed that the inbred lines carrying ATT at BnaA10g09290Hap1 and AAT at BnaC08g26640Hap1 had greater plant height than lines carrying other haplotype alleles at low P supply.ConclusionOur results demonstrate the power of GWAS in identifying genes of interest in B. napus and provided insights into the genetic basis of plant height and branch number at low P supply in B. napus. Candidate genes and favourable haplotypes may facilitate marker-based breeding efforts aimed at improving P use efficiency in B. napus.     

Genetic analysis of tolerance to low-phosphorus stress in maize using restriction fragment length polymorphisms

Summary An understanding of the genetic nature underlying tolerance to low-phosphorus (low-P) stress could aid in the efficient development of tolerant plant strains. The objective of this study was to identify the number of loci in a maize (Zea mays L.) population segregating for tolerance to low-P stress, their approximate location, and the magnitude of their effect.Seventy-seven restriction fragment length polymorphisms (RFLPs) were identified and scored in a maize F₂ population derived from a cross between line NY821 and line H99. The F₂ individuals were self-pollinated to produce F₃ families. Ninety F₃ families were grown in a sand-alumina system, which simulated diffusion-limited, low-P soil conditions. The F₃ families were evaluated for vegetative growth in a controlled-environment experiment. To identify quantitative trait loci (QTLs) underlying tolerance to low-P stress, the mean phenotypic performances of the F₃ families were contrasted based on genotypic classification at each of 77 RFLP marker loci.Six RFLP marker loci were significantly associated with performance under low-P stress (P<0.01). One marker locus accounted for 25% of the total phenotypic variation. Additive gene action was predominant for all of the QTLs identified. Significant marker loci were located on four separate chromosomes representing five unlinked genomic regions. Two marker loci were associated with an additive by additive epistatic interaction. A multiple regression model including three marker loci and the significant epistatic interaction accounted for 46% of the total phenotypic variation. Heterozygosity per se was not predictive of phenotypic performance.     

Meta-analysis combined with syntenic metaQTL mining dissects candidate loci for maize yield

Yield potential and stability improvement with the goal of ensuring global food security is an important priority. Yield has a quantitative nature and is controlled by quantitative trait loci (QTL) and environmental factors. An increasingly large number of maize yield QTL have been identified, and how to integrate and re-analyze them is challenging. To this end, we tried to combine QTL meta-analysis with homology-based cloning techniques to dissect candidate loci/genes for maize yield. We first collected maize yield-related QTL from public resources. Then, 351 collected QTL were iteratively projected and meta-analyzed to obtain metaQTL (MQTL). A total of 54 MQTL were identified and tended to cluster in the maize genome. Seven MQTL containing ten maize orthologs of rice yield genes were dissected and temporarily termed syntenic MQTL. Maize orthologs of three functionally-characterized rice yield genes, GIF1, WFP/IPA1, and DEP1, were specially selected to undergo phylogenetic, proliferation, and selective pattern analysis. The results showed that maize orthologs were closely related to rice yield genes and subjected to mixed selective pressures, including positive selection during selective sweeps. The power of the combined techniques reported here was primarily validated not only by the congruency of MQTL and recently reported maize yield QTL but also by mined syntenic MQTL containing the well-characterized Miniature1 (Mn1) gene for maize kernel size and weight determination. Maize MQTL, especially syntenic MQTL regions, could serve not only for QTL fine-mapping and cloning but also for the marker-assisted selection breeding program. The maize yield candidate loci/genes presented here also deserve further investigation and will provide clues to the molecular bases of grain yield. Additionally, the combined technique described here will find its way into further quantitative trait research.     

Genome-wide association mapping identifies markers associated with cane yield components and sucrose traits in the Louisiana sugarcane core collection

Sugarcane is an economically important crop for both food and biofuel industries. Marker-assisted breeding in sugarcane is becoming a reality with the recent development and deployment of markers linked with disease resistance genes. Large linkage disequilibrium in sugarcane makes genome-wide association studies (GWAS) a better alternative to biparental mapping to identify markers associated with agronomic traits. GWAS was conducted on a Louisiana core collection to identify marker-trait associations (MTA) for 11 cane yield and sucrose traits using single nucleotide polymorphism (SNP) and insertion-deletion (Indel) markers. Significant (P < .05) MTAs were identified for all traits where the top ranked markers explained up to 15% of the total phenotypic variation. High correlations (0.732 to 0.999) were observed between sucrose traits and 56 markers were found consistent across multiple traits. These markers following validation in more diverse populations could be used in marker-assisted selection of clones in sugarcane breeding program in Louisiana and elsewhere.     

Meta-analysis and candidate gene mining of low-phosphorus tolerance in maize

Plants with tolerance to low-phosphorus （P） can grow better under Iow-P conditions, and understanding of genetic mechanisms of Iow-P tolerance can not only facilitate identifying relevant genes but also help to develop Iow-P tolerant cultivars. QTL meta-analysis was conducted after a comprehensive review of the reports on O, TL mapping for Iow-P tolerance-related traits in maize. Meta-analysis pro- duced 23 consensus QTL （cQTL）, 17 of which located in similar chromosome regions to those previously reported to influence root traits. Meanwhile, candidate gene mining yielded 215 genes, 22 of which located in the cQTL regions. These 22 genes are homologous to 14 functionally characterized genes that were found to participate in plant Iow-P tolerance, including genes encoding miR399s, Pi transporters and purple acid phosphatases. Four cQTL loci （cQTL2-L cQTLS-3, cQTL6-2, and cQTL10-2） may play important roles for low-P tolerance because each contains more original QTL and has better consistency across previous reports.     

Genome-wide association mapping of yield components and drought tolerance-related traits in cotton

Drought causes serious yield losses in cotton production throughout the world. Association mapping allows identification and localization of the genes controlling drought-related traits which will be helpful in cotton breeding. In the present study, genetic diversity analysis and association mapping of yield and drought traits were performed on a panel of 99 upland cotton genotypes using 177 SSR (simple sequence repeat) markers. Yield parameters and drought tolerance-related traits were evaluated for two seasons under two watering regimes: water-stressed and well-watered. The traits included seed cotton yield (SCY), lint yield (LY), lint percentage (LP), water-use efficiency (WUE), yield potential (YP), yield reduction (YR), yield index (YI), drought sensitivity index (DSI), stress tolerance index (STI), harmonic mean (HM), and geometric mean productivity (GMP). The genotypes with the least change in seed cotton yield under drought stress were Zeta 2, Delcerro, Nazilli 87, and DAK 66/3 which were also the most water-use efficient cultivars. The average genetic diversity of the panel was 0.38. The linkage disequilibrium decayed relatively rapidly at 20–30 cM (r²?≥?0.5). We identified 30 different SSR markers associated with the traits. Fifteen and 23 SSR markers were linked to the traits under well-watered and water-stress conditions, respectively. To our knowledge, most of these quantitative yield and drought tolerance-associated loci were newly identified. The genetic diversity and association mapping results should facilitate the development of drought-tolerant cotton lines with high yield in molecular breeding programs.     

Genome-wide association mapping of yield and yield components of spring wheat under contrasting moisture regimes

Key message

A stable QTL that may be used in marker-assisted selection in wheat breeding programs was detected for yield, yield components and drought tolerance-related traits in spring wheat association mapping panel.

Abstract

Genome-wide association mapping has become a widespread method of quantitative trait locus (QTL) identification for many crop plants including wheat (Triticum aestivum L.). Its benefit over traditional bi-parental mapping approaches depends on the extent of linkage disequilibrium in the mapping population. The objectives of this study were to determine linkage disequilibrium decay rate and population structure in a spring wheat association mapping panel (n = 285–294) and to identify markers associated with yield and yield components, morphological, phenological, and drought tolerance-related traits. The study was conducted under fully irrigated and rain-fed conditions at Greeley, CO, USA and Melkassa, Ethiopia in 2010 and 2011 (five total environments). Genotypic data were generated using diversity array technology markers. Linkage disequilibrium decay rate extended over a longer genetic distance for the D genome (6.8 cM) than for the A and B genomes (1.7 and 2.0 cM, respectively). Seven subpopulations were identified with population structure analysis. A stable QTL was detected for grain yield on chromosome 2DS both under irrigated and rain-fed conditions. A multi-trait region significant for yield and yield components was found on chromosome 5B. Grain yield QTL on chromosome 1BS co-localized with harvest index QTL. Vegetation indices shared QTL with harvest index on chromosome 1AL and 5A. After validation in relevant genetic backgrounds and environments, QTL detected in this study for yield, yield components and drought tolerance-related traits may be used in marker-assisted selection in wheat breeding programs.     

Genome-wide association study identifies genetic loci associated with iron deficiency

The existence of multiple inherited disorders of iron metabolism in man, rodents and other vertebrates suggests genetic contributions to iron deficiency. To identify new genomic locations associated with iron deficiency, a genome-wide association study (GWAS) was performed using DNA collected from white men aged≥25 y and women≥50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF)≤12 μg/L (cases) and iron replete controls (SF>100 μg/L in men, SF>50 μg/L in women). Regression analysis was used to examine the association between case-control status (336 cases, 343 controls) and quantitative serum iron measures and 331,060 single nucleotide polymorphism (SNP) genotypes, with replication analyses performed in a sample of 71 cases and 161 controls from a population of white male and female veterans screened at a US Veterans Affairs (VA) medical center. Five SNPs identified in the GWAS met genome-wide statistical significance for association with at least one iron measure, rs2698530 on chr. 2p14; rs3811647 on chr. 3q22, a known SNP in the transferrin (TF) gene region; rs1800562 on chr. 6p22, the C282Y mutation in the HFE gene; rs7787204 on chr. 7p21; and rs987710 on chr. 22q11 (GWAS observed P<1.51×10(-7) for all). An association between total iron binding capacity and SNP rs3811647 in the TF gene (GWAS observed P=7.0×10(-9), corrected P=0.012) was replicated within the VA samples (observed P=0.012). Associations with the C282Y mutation in the HFE gene also were replicated. The joint analysis of the HEIRS and VA samples revealed strong associations between rs2698530 on chr. 2p14 and iron status outcomes. These results confirm a previously-described TF polymorphism and implicate one potential new locus as a target for gene identification.     

Genome-wide association study identifies variants associated with hair length in Brangus cattle

Thermal stress limits beef cattle production and a shorter hair coat is a key thermoregulative adaptation that allows cattle to lose heat more efficiently. The objective of this study was to identify genetic variants associated with the length of the undercoat and topcoat of cattle utilizing 1456 Brangus heifers genotyped with the Bovine GGP F250 array. Seven SNPs in the PCCA gene were significantly associated with undercoat length. PCCA belongs to the biotin transport and metabolism pathway. Biotin deficiency has been reported to cause hair loss. Four SNPs in an 110 kb including a missense mutation in the PRLR gene were significantly associated with topcoat length. Whereas the association of this polymorphism with hair length is novel, the SLICK mutation in PRLR has previously been demonstrated to significantly impact hair length in cattle. These newly detected genetic variants may contribute to a shorter hair coat and more thermotolerant animals.     

Genome-wide association study of the level of blood components in Pekin ducks

Blood components are considered to reflect nutrient metabolism and immune activity in both humans and animals. In this study, we measured 12 blood components in Pekin ducks and performed genome-wide association analysis to identify the QTLs (quantitative trait locus) using a genotyping-by-sequencing strategy. A total of 54 QTLs were identified for blood components. One genome-wide significant QTL for alkaline phosphatase was identified within the intron-region of the OTOG gene (P = 1.31E-07). Moreover, 21 genome-wide significant SNPs for the level of serum cholinesterase were identified on six different scaffolds. In addition, for serum calcium, one genome-wide significant QTL was identified in the upstream region of gene RAB11B. These results provide new markers for functional studies in Pekin ducks, and several candidate genes were identified, which may provide additional insights into specific mechanisms for blood metabolism in ducks and their potential application for duck breeding programs.     

Genome-wide association study identifies four loci associated with eruption of permanent teeth

The sequence and timing of permanent tooth eruption is thought to be highly heritable and can have important implications for the risk of malocclusion, crowding, and periodontal disease. We conducted a genome-wide association study of number of permanent teeth erupted between age 6 and 14 years, analyzed as age-adjusted standard deviation score averaged over multiple time points, based on childhood records for 5,104 women from the Danish National Birth Cohort. Four loci showed association at P<5×10(-8) and were replicated in four independent study groups from the United States and Denmark with a total of 3,762 individuals; all combined P-values were below 10(-11). Two loci agreed with previous findings in primary tooth eruption and were also known to influence height and breast cancer, respectively. The two other loci pointed to genomic regions without any previous significant genome-wide association study results. The intronic SNP rs7924176 in ADK could be linked to gene expression in monocytes. The combined effect of the four genetic variants was most pronounced between age 10 and 12 years, where children with 6 to 8 delayed tooth eruption alleles had on average 3.5 (95% confidence interval: 2.9-4.1) fewer permanent teeth than children with 0 or 1 of these alleles.     

Genome-wide association study for candidate genes controlling seed yield and its components in rapeseed (Brassica napus subsp. napus)

Genetic improvement of seed yield per plant (SY) is one of the major objectives in Brassica napus breeding programme. SY, being a complex quantitative trait is directly and indirectly influenced by yield-component traits such as siliqua length (SL), number of seeds per siliqua (NSS), and thousand seed weight (TSW). Therefore, concurrent improvement in SL, NSS and TSW can lead to higher SY in B. napus. This study was conducted to identify significant SNPs and putative candidate genes governing SY and its component traits (SL, NSS, TSW). All these traits were evaluated in a diverse set of 200 genotypes representing diversity from wide geographical locations. Of these, a set of 125 genotypes were chosen based on pedigree diversity and multi-location trait variation for genotyping by sequencing (GBS). Best linear unbiased predictors (BLUPs) of all the traits were used for genome-wide association study (GWAS) with 85,126 SNPs obtained from GBS. A total of 16, 18, 27 and 18 SNPs were found to be significantly associated for SL, NSS, TSW and SY respectively. Based on linkage disequilibrium decay analysis, 150 kb genomic region flanking the SNP was used for the identification of underlying candidate genes for each test trait. Important candidate genes involved in phytohormone signaling (WAT1, OSR1, ARR8, CKX1, REM7, REM9, BG1) and seed storage proteins (Cruciferin) were found to have significant influence on seed weight and yield. Genes involved in sexual reproduction and fertilization (PERK7, PERK13, PRK3, GATA15, NFD6) were found to determine the number of seeds per siliqua. Several genes found in this study namely ATS3A, CKX1, SPL2, SPL6, SPL9, WAT1 showed pleiotropic effect with yield component traits. Significant SNPs and putative candidate genes identified for SL, NSS, TSW and SY could be used in marker-assisted breeding for improvement of crop yield in B. napus. Genotypes identified with high SL, NSS, TSW and SY could serve as donors in crop improvement programs in B. napus.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-01060-9.     

Genome-wide association study identifies novel loci associated with circulating phospho- and sphingolipid concentrations

Phospho- and sphingolipids are crucial cellular and intracellular compounds. These lipids are required for active transport, a number of enzymatic processes, membrane formation, and cell signalling. Disruption of their metabolism leads to several diseases, with diverse neurological, psychiatric, and metabolic consequences. A large number of phospholipid and sphingolipid species can be detected and measured in human plasma. We conducted a meta-analysis of five European family-based genome-wide association studies (N = 4034) on plasma levels of 24 sphingomyelins (SPM), 9 ceramides (CER), 57 phosphatidylcholines (PC), 20 lysophosphatidylcholines (LPC), 27 phosphatidylethanolamines (PE), and 16 PE-based plasmalogens (PLPE), as well as their proportions in each major class. This effort yielded 25 genome-wide significant loci for phospholipids (smallest P-value = 9.88×10⁻²⁰⁴) and 10 loci for sphingolipids (smallest P-value = 3.10×10⁻⁵⁷). After a correction for multiple comparisons (P-value<2.2×10⁻⁹), we observed four novel loci significantly associated with phospholipids (PAQR9, AGPAT1, PKD2L1, PDXDC1) and two with sphingolipids (PLD2 and APOE) explaining up to 3.1% of the variance. Further analysis of the top findings with respect to within class molar proportions uncovered three additional loci for phospholipids (PNLIPRP2, PCDH20, and ABDH3) suggesting their involvement in either fatty acid elongation/saturation processes or fatty acid specific turnover mechanisms. Among those, 14 loci (KCNH7, AGPAT1, PNLIPRP2, SYT9, FADS1-2-3, DLG2, APOA1, ELOVL2, CDK17, LIPC, PDXDC1, PLD2, LASS4, and APOE) mapped into the glycerophospholipid and 12 loci (ILKAP, ITGA9, AGPAT1, FADS1-2-3, APOA1, PCDH20, LIPC, PDXDC1, SGPP1, APOE, LASS4, and PLD2) to the sphingolipid pathways. In large meta-analyses, associations between FADS1-2-3 and carotid intima media thickness, AGPAT1 and type 2 diabetes, and APOA1 and coronary artery disease were observed. In conclusion, our study identified nine novel phospho- and sphingolipid loci, substantially increasing our knowledge of the genetic basis for these traits.     

Genome-wide association study (GWAS) of salt tolerance in worldwide soybean germplasm lines

Salt is a severe abiotic stress causing soybean yield loss in saline soils and irrigated fields. Marker-assisted selection (MAS) is a powerful genomic tool for improving the efficiency of breeding salt-tolerant soybean varieties. The objectives of this study were to uncover novel single-nucleotide polymorphisms (SNPs) and quantitative trait loci (QTLs) associated with salt tolerance and to confirm the previously identified genomic regions and SNPs for salt tolerance. A total of 283 diverse soybean plant introductions (PIs) were screened for salt tolerance in the greenhouse based on leaf chloride concentrations and leaf chlorophyll concentrations after 12–18 days of 120-mM NaCl treatment. A total of 33,009 SNPs across 283 genotypes from the Illumina Infinium SoySNP50K BeadChip database were employed in the association analysis with leaf chloride concentrations and leaf chlorophyll concentrations. Genome-wide association mapping showed that 45 SNPs representing nine genomic regions on chromosomes (Chr.) 2, 3, 7, 8, 10, 13, 14, 16, and 20 were significantly associated with both leaf chloride concentrations and leaf chlorophyll concentrations in 2014, 2015, and combined years. A total of 31 SNPs on Chr. 3 were mapped at or near the previously reported major salt tolerance QTL. The significant SNP on Chr. 2 was also in proximity to the previously reported SNP for salt tolerance. The other significant SNPs represent seven putative novel QTLs for salt tolerance. The significant SNP markers on Chr. 2, 3, 14, 16, and 20, which were identified in both general linear model and mixed linear model, were highly recommended for MAS in breeding salt-tolerant soybean varieties.     

Genome-wide association study reveals novel loci associated with fat-deposition and meat-quality traits in Pekin ducks

Meat-quality traits play an essential role in meat poultry production. To determine the genetic mechanisms of meat quality in Pekin ducks, we performed a large-scale GWAS to identify quantitative trait loci affecting meat quality in Pekin ducks. We measured 10 traits in 542 Pekin ducks and genotyped each duck using genotyping-by-sequencing. The genetic parameters (genomic heritability, genetic correlation) for 10 meat-quality related traits were evaluated. Based on the large genotype–phenotype dataset, we performed GWASs for all of these traits. A total of 33 significant QTL (P < 3.03 × 10⁻⁵) across 13 chromosomes were identified by loci-based analysis. Some newly identified candidate genes were discovered for fat-deposition and meat-quality traits, including PAG1 for body weight and eviscerated weight, INTU and NUP35 for abdominal fat weight and ratio, NUP3 and ARHGDIB for skin fat weight and ratio, GOLGA5 for breast muscle toughness and breast tenderness, and CTDSPL and PKP1 for breast muscle thickness. The current study is the first systematic report regarding duck meat quality.     

Genome-wide association study to identify genetic loci associated with gastrointestinal nematode resistance in Katahdin sheep

Resistance to gastrointestinal nematodes has previously been shown to be a moderately heritable trait in some breeds of sheep, but the mechanisms of resistance are not well understood. Selection for resistance currently relies upon faecal egg counts (FEC), blood packed cell volumes and FAMACHA visual indicator scores of anaemia. Identifying genomic markers associated with disease resistance would potentially improve the selection process and provide a more reliable means of classifying and understanding the biology behind resistant and susceptible sheep. A GWAS was conducted to identify possible genetic loci associated with resistance to Haemonchus contortus in Katahdin sheep. Forty animals were selected from the top and bottom 10% of estimated breeding values for FEC from a total pool of 641 sires and ram lambs. Samples were genotyped using Applied Biosystems™ Axiom™ Ovine Genotyping Array (50K) consisting of 51 572 SNPs. Following quality control, 46 268 SNPs were included in subsequent analyses. Analyses were conducted using a linear regression model in plink v1.90 and a single-locus mixed model in snp and variation suite . Genome-wide significance was determined by a Bonferroni correction for multiple testing. Using linear regression, loci on chromosomes 2, 3, 16, 23 and 24 were significantly associated at the genome level with FEC estimated breeding values, and we identified a region on chromosome 2 that was significant using both statistical analyses. We suggest a potential role for the gene DIS3L2 for gastrointestinal nematode resistance in Katahdin sheep, although further research is needed to validate these findings.     

Genome-wide association study identified a narrow chromosome 1 region associated with chicken growth traits

Chicken growth traits are important economic traits in broilers. A large number of studies are available on finding genetic factors affecting chicken growth. However, most of these studies identified chromosome regions containing putative quantitative trait loci and finding causal mutations is still a challenge. In this genome-wide association study (GWAS), we identified a narrow 1.5 Mb region (173.5-175 Mb) of chicken (Gallus gallus) chromosome (GGA) 1 to be strongly associated with chicken growth using 47,678 SNPs and 489 F2 chickens. The growth traits included aggregate body weight (BW) at 0-90 d of age measured weekly, biweekly average daily gains (ADG) derived from weekly body weight, and breast muscle weight (BMW), leg muscle weight (LMW) and wing weight (WW) at 90 d of age. Five SNPs in the 1.5 Mb KPNA3-FOXO1A region at GGA1 had the highest significant effects for all growth traits in this study, including a SNP at 8.9 Kb upstream of FOXO1A for BW at 22-48 d and 70 d, a SNP at 1.9 Kb downstream of FOXO1A for WW, a SNP at 20.9 Kb downstream of ENSGALG00000022732 for ADG at 29-42 d, a SNP in INTS6 for BW at 90 d, and a SNP in KPNA3 for BMW and LMW. The 1.5 Mb KPNA3-FOXO1A region contained two microRNA genes that could bind to messenger ribonucleic acid (mRNA) of IGF1, FOXO1A and KPNA3. It was further indicated that the 1.5 Mb GGA1 region had the strongest effects on chicken growth during 22-42 d.