Transcriptional regulation of transcription factor genes WRI1 and LAFL during Brassica napus seed development.

The WRINKLED1 (WRI1) and LAFL [LEAFY COTYLEDON1 (LEC1), ABSCISIC ACID INSENSITIVE3 (ABI3), FUSCA3 (FUS3), and LEC2] transcription factors play essential roles in governing seed development and oil biosynthesis. To gain a comprehensive understanding of the transcriptional regulation of WRI1 and LAFL, we conducted genome-wide association studies for the expression profiles of WRI1 and LAFL in developing seeds at 20 and 40 days after flowering (DAF) using 302 rapeseed (Brassica napus) accessions. We identified a total of 237 expression quantitative trait nucleotides (eQTNs) and 51 expression QTN-by-environment interactions (eQEIs) associated with WRI1 and LAFL. Around these eQTNs and eQEIs, we pinpointed 41 and 8 candidate genes with known transcriptional regulations or protein interactions with their expression traits, respectively. Based on RNA-seq and ATAC-seq data, we employed the XGBoost and Basenji models which predicted 15 candidate genes potentially regulating the expression of WRI1 and LAFL. We further validated the predictions via tissue expression profile, haplotype analysis, and expression correlation analysis, and verified the transcriptional activation activity of BnaC03.MYB56 (R2R3-MYB transcription factor 56) on the expression of BnaA09.LEC1 by dual-luciferase reporter and yeast one-hybrid assays. BnaA10.AGL15 (AGAMOUS-LIKE 15), BnaC04.VAL1 (VIVIPAROUS1/ABSCISIC ACID INSENSITIVE3-LIKE 1), BnaC03.MYB56, and BnaA10.MYB56 were co-expressed with WRI1 and LAFL at 20 DAF in M35, a key module for seed development and oil biosynthesis. We further validated the positive regulation of MYB56 on seed oil accumulation using Arabidopsis (Arabidopsis thaliana) mutants. This study not only delivers a framework for future eQEI identification but also offers insights into the developmental regulation of seed oil accumulation.

Identification of QTNs, QTN-by-environment interactions, and their candidate genes for salt tolerance related traits in soybean.

BACKGROUND: Salt stress significantly reduces soybean yield. To improve salt tolerance in soybean, it is important to mine the genes associated with salt tolerance traits. RESULTS: Salt tolerance traits of 286 soybean accessions were measured four times between 2009 and 2015. The results were associated with 740,754 single nucleotide polymorphisms (SNPs) to identify quantitative trait nucleotides (QTNs) and QTN-by-environment interactions (QEIs) using three-variance-component multi-locus random-SNP-effect mixed linear model (3VmrMLM). As a result, eight salt tolerance genes (GmCHX1, GsPRX9, Gm5PTase8, GmWRKY, GmCHX20a, GmNHX1, GmSK1, and GmLEA2-1) near 179 significant and 79 suggested QTNs and two salt tolerance genes (GmWRKY49 and GmSK1) near 45 significant and 14 suggested QEIs were associated with salt tolerance index traits in previous studies. Six candidate genes and three gene-by-environment interactions (GEIs) were predicted to be associated with these index traits. Analysis of four salt tolerance related traits under control and salt treatments revealed six genes associated with salt tolerance (GmHDA13, GmPHO1, GmERF5, GmNAC06, GmbZIP132, and GmHsp90s) around 166 QEIs were verified in previous studies. Five candidate GEIs were confirmed to be associated with salt stress by at least one haplotype analysis. The elite molecular modules of seven candidate genes with selection signs were extracted from wild soybean, and these genes could be applied to soybean molecular breeding. Two of these genes, Glyma06g04840 and Glyma07g18150, were confirmed by qRT-PCR and are expected to be key players in responding to salt stress. CONCLUSIONS: Around the QTNs and QEIs identified in this study, 16 known genes, 6 candidate genes, and 8 candidate GEIs were found to be associated with soybean salt tolerance, of which Glyma07g18150 was further confirmed by qRT-PCR.

Novel resistance loci for quantitative resistance to Septoria tritici blotch in Asian wheat (Triticum aestivum) via genome-wide association study.

BACKGROUND: Septoria tritici blotch (STB) disease causes yield losses of up to 50 per cent in susceptible wheat cultivars and can reduce wheat production. In this study, genomic architecture for adult-plant STB resistance in a Septoria Association Mapping Panel (SAMP) having 181 accessions and genomic regions governing STB resistance in a South Asian wheat panel were looked for. RESULTS: Field experiments during the period from 2019 to 2021 revealed those certain accessions, namely BGD52 (CHIR7/ANB//CHIR1), BGD54 (CHIR7/ANB//CHIR1), IND92 (WH 1218), IND8 (DBW 168), and IND75 (PBW 800), exhibited a high level of resistance. Genetic analysis revealed the presence of 21 stable quantitative trait nucleotides (QTNs) associated with resistance to STB (Septoria tritici blotch) on all wheat chromosomes, except for 2D, 3A, 3D, 4A, 4D, 5D, 6B, 6D, and 7A. These QTNs were predominantly located in chromosome regions previously identified as associated with STB resistance. Three Quantitative Trait Loci (QTNs) were found to have significant phenotypic effects in field evaluations. These QTNs are Q.STB.5A.1, Q.STB.5B.1, and Q.STB.5B.3. Furthermore, it is possible that the QTNs located on chromosomes 1A (Q.STB.1A.1), 2A (Q.STB_DH.2A.1, Q.STB.2A.3), 2B (Q.STB.2B.4), 5A (Q.STB.5A.1, Q.STB.5A.2), and 7B (Q.STB.7B.2) could potentially be new genetic regions associated with resistance. CONCLUSION: Our findings demonstrate the importance of Asian bread wheat as a source of STB resistance alleles and novel stable QTNs for wheat breeding programs aiming to develop long-lasting and wide-ranging resistance to Zymoseptoria tritici in wheat cultivars.

New Insights into the Genetic Basis of Lysine Accumulation in Rice Revealed by Multi-Model GWAS.

Lysine is an essential amino acid that cannot be synthesized in humans. Rice is a global staple food for humans but has a rather low lysine content. Identification of the quantitative trait nucleotides (QTNs) and genes underlying lysine content is crucial to increase lysine accumulation. In this study, five grain and three leaf lysine content datasets and 4,630,367 single nucleotide polymorphisms (SNPs) of 387 rice accessions were used to perform a genome-wide association study (GWAS) by ten statistical models. A total of 248 and 71 common QTNs associated with grain/leaf lysine content were identified. The accuracy of genomic selection/prediction RR-BLUP models was up to 0.85, and the significant correlation between the number of favorable alleles per accession and lysine content was up to 0.71, which validated the reliability and additive effects of these QTNs. Several key genes were uncovered for fine-tuning lysine accumulation. Additionally, 20 and 30 QTN-by-environment interactions (QEIs) were detected in grains/leaves. The QEI-sf0111954416 candidate gene LOC_Os01g21380 putatively accounted for gene-by-environment interaction was identified in grains. These findings suggested the application of multi-model GWAS facilitates a better understanding of lysine accumulation in rice. The identified QTNs and genes hold the potential for lysine-rich rice with a normal phenotype.

A Systemic Investigation of Genetic Architecture and Gene Resources Controlling Kernel Size-Related Traits in Maize.

Grain yield is the most critical and complex quantitative trait in maize. Kernel length (KL), kernel width (KW), kernel thickness (KT) and hundred-kernel weight (HKW) associated with kernel size are essential components of yield-related traits in maize. With the extensive use of quantitative trait locus (QTL) mapping and genome-wide association study (GWAS) analyses, thousands of QTLs and quantitative trait nucleotides (QTNs) have been discovered for controlling these traits. However, only some of them have been cloned and successfully utilized in breeding programs. In this study, we exhaustively collected reported genes, QTLs and QTNs associated with the four traits, performed cluster identification of QTLs and QTNs, then combined QTL and QTN clusters to detect consensus hotspot regions. In total, 31 hotspots were identified for kernel size-related traits. Their candidate genes were predicted to be related to well-known pathways regulating the kernel developmental process. The identified hotspots can be further explored for fine mapping and candidate gene validation. Finally, we provided a strategy for high yield and quality maize. This study will not only facilitate causal genes cloning, but also guide the breeding practice for maize.

Powerful Identification of Large Quantitative Trait Loci Using Genome-Wide R/glmnet-Based Regression.

R/glmnet has been successfully applied to jointly mapped multiple quantitative trait loci for linkage analysis, along with statistical inference for quantitative trait loci candidates with nonzero genetic effects using R/lm for normally distributed traits, R/glm for discrete traits, and R/coxph for survival times. In this study, we extended R/glmnet to a genome-wide association study by means of parallel computation. A multi-locus genome-wide association study for high-throughput single-nucleotide polymorphisms was implemented in the "Multi-Runking" software written within the R workspace. This software can better detect common and large quantitative trait nucleotides and more accurately estimate than genome-wide mixed model analysis for one single-nucleotide polymorphism at a time and linear mixed models-least absolute shrinkage and selection operator. Its applicability and utility were demonstrated by multi-locus genome-wide association studies for the simulated and real traits distributed normally, binary traits, and survival times.

Genetic dissection of flowering time in flax (Linum usitatissimum L.) through single- and multi-locus genome-wide association studies.

In a rapidly changing climate, flowering time (FL) adaptation is important to maximize seed yield in flax (Linum usitatissimum L.). However, our understanding of the genetic mechanism underlying FL in this multipurpose crop remains limited. With the aim of dissecting the genetic architecture of FL in flax, a genome-wide association study (GWAS) was performed on 200 accessions of the flax core collection evaluated in four environments. Two single-locus and six multi-locus models were applied using 70,935 curated single nucleotide polymorphism (SNP) markers. A total of 40 quantitative trait nucleotides (QTNs) associated with 27 quantitative trait loci (QTL) were identified in at least two environments. The number of QTL with positive-effect alleles in accessions was significantly correlated with FL (r = 0.77 to 0.82), indicating principally additive gene actions. Nine QTL were significant in at least three of the four environments accounting for 3.06-14.71% of FL variation. These stable QTL spanned regions that harbored 27 Arabidopsis thaliana and Oryza sativa FL-related orthologous genes including FLOWERING LOCUS T (Lus10013532), FLOWERING LOCUS D (Lus10028817), transcriptional regulator SUPERMAN (Lus10021215), and gibberellin 2-beta-dioxygenase 2 (Lus10037816). In silico gene expression analysis of the 27 FL candidate gene orthologous suggested that they might play roles in the transition from vegetative to reproductive phase, flower development and fertilization. Our results provide new insights into the QTL architecture of flowering time in flax, identify potential candidate genes for further studies, and demonstrate the effectiveness of combining different GWAS models for the genetic dissection of complex traits.

Mapping QTL/QTN and mining candidate genes for plant height and its response to planting densities in soybean [Glycine max (L.) Merr.] through a FW-RIL population.

Plant height (PH) determines the morphology and seed yield of soybean, so it is an important breeding target, which is controlled by multiple genes and affected by plant density. In this research, it was used about a four-way recombinant inbred lines (FW-RIL) with 144 families constructed by double cross (Kenfeng 14 × Kenfeng 15) × (Heinong 48 × Kenfeng 19) as experimental materials, with the purpose to map QTL/QTN associated with PH under densities of 2.2×105 plant/ha (D1) and 3×105 plant/ha (D2) in five environments. The results showed that response of PH to densities varied in accordance to genotypes among environments. A total of 26 QTLs and 13 QTNs were identified specifically in D1; 20 QTLs and 21 QTNs were identified specifically in D2. Nine QTLs and one QTN were discovered commonly in two densities. Fifteen QTLs and 9 QTNs were repeatedly detected by multiple statistical methods, densities, or environments, which could be considered stable. Eighteen QTLs were detected, as well as 7 QTNs underlying responses of PH to density increment. Six QTNs, co-located in the interval of QTL, were detected in more than two environments or methods with a longer genome length over 3000 kb. Based on gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, five genes were predicted as candidates, which were likely to be involved in growth and development of PH. The results will help elucidate the genetic basis and improve molecular assistant selection of PH. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01209-0.

Short communication: Replication of genome-wide association studies for milk production traits in Chinese Holstein by an efficient rotated linear mixed model.

Milk is regarded as an important nutrient for humans, and Chinese Holstein cows provide high-quality milk for billions of Chinese people. Therefore, detecting quantitative trait nucleotides (QTN) or candidate genes for milk production traits in Chinese Holstein is important. In this study, we performed genome-wide association studies (GWAS) in a Chinese Holstein population of 6,675 cows and 71,633 SNP using deregressed proofs (DRP) as phenotypes to replicate our previous study in a population of 1,815 cows and 39,163 SNP using estimated breeding values (EBV) as phenotypes. The associations between 3 milk production traits-milk yield (MY), fat percentage (FP), and protein percentage (PP)-and the SNP were determined by using an efficient rotated linear mixed model, which benefits from linear transformations of genomic estimated values and Eigen decomposition of the genomic relationship matrix algorithm. In total, we detected 94 SNP that were significantly associated with one or more milk production traits, including 7 SNP for MY, 76 for FP, and 36 for PP; 87% of these SNP were distributed across Bos taurus autosomes 14 and 20. In total, 83 SNP were found to be located within the reported quantitative trait loci (QTL) regions, and one novel segment (between 1.41 and 1.49 Mb) on chromosome 14 was significantly associated with FP, which could be an important candidate QTL region. In addition, the detected intervals were narrowed down from the reported regions harboring causal variants. The top significant SNP for the 3 traits was ARS-BFGL-NGS-4939, which is located within the DGAT1 gene. Five detected genes (CYHR1, FOXH1, OPLAH, PLEC, VPS28) have effects on all 3 traits. Our study provides a suite of QTN, candidate genes, and a novel QTL associated with milk production traits, and thus forms a solid basis for genomic selection and molecular breeding for milk production traits in Chinese Holstein.

QTL mapping of fungicide sensitivity reveals novel genes and pleiotropy with melanization in the pathogen Zymoseptoria tritici.

A major problem associated with the intensification of agriculture is the emergence of fungicide resistance. Azoles are ergosterol biosynthesis inhibitors that have been widely used in agriculture and medicine since the 1970s, leading to emergence of increasingly resistant fungal populations. The known genetic mechanisms underlying lower azole sensitivity include mutations affecting the CYP51 gene that encodes the target protein, but in many cases azole resistance is a more complex trait with an unknown genetic basis. We used quantitative trait locus (QTL) mapping to identify genes affecting azole sensitivity in two crosses of Zymoseptoria tritici, the most damaging wheat pathogen in Europe. Restriction site associated DNA sequencing (RADseq) was used to genotype 263 (cross 1) and 261 (cross 2) progeny at ∼ 8500 single nucleotide polymorphisms (SNP) and construct two dense linkage maps. Azole sensitivity was assessed using high-throughput digital image analysis of colonies growing on Petri dishes with or without the fungicide propiconazole. We identified three QTLs for azole sensitivity, including two that contained novel fungicide sensitivity genes. One of these two QTLs contained only 16 candidate genes, among which four most likely candidates were identified. The third QTL contained ERG6, encoding another protein involved in ergosterol biosynthesis. Known genes in QTLs affecting colony growth included CYP51 and PKS1, a gene affecting melanization in Z. tritici. PKS1 showed compelling evidence for pleiotropy, with a rare segregating allele that increased melanization while decreasing growth rate and propiconazole sensitivity. This study resolved the genetic architecture of an important agricultural trait and led to identification of novel genes that are likely to affect azole sensitivity in Z. tritici. It also provided insight into fitness costs associated with lowered azole sensitivity and suggests a novel fungicide mixture strategy.

Short communication: Validation of 4 candidate causative trait variants in 2 cattle breeds using targeted sequence imputation.

Most association studies for pinpointing trait-associated variants are performed within breed. The availability of sequence data from key ancestors of several cattle breeds now enables immediate assessment of the frequency of trait-associated variants in populations different from the mapping population and their imputation into large validation populations. The objective of this study was to validate the effects of 4 putatively causative variants on milk production traits, male fertility, and stature in German Fleckvieh and Holstein-Friesian animals using targeted sequence imputation. We used whole-genome sequence data of 456 animals to impute 4 missense mutations in DGAT1, GHR, PRLR, and PROP1 into 10,363 Fleckvieh and 8,812 Holstein animals. The accuracy of the imputed genotypes exceeded 95% for all variants. Association testing with imputed variants revealed consistent antagonistic effects of the DGAT1 p.A232K and GHR p.F279Y variants on milk yield and protein and fat contents, respectively, in both breeds. The allele frequency of both polymorphisms has changed considerably in the past 20 yr, indicating that they were targets of recent selection for milk production traits. The PRLR p.S18N variant was associated with yield traits in Fleckvieh but not in Holstein, suggesting that it may be in linkage disequilibrium with a mutation affecting yield traits rather than being causal. The reported effects of the PROP1 p.H173R variant on milk production, male fertility, and stature could not be confirmed. Our results demonstrate that population-wide imputation of candidate causal variants from sequence data is feasible, enabling their rapid validation in large independent populations.

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis.

Introduction: An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security. Methods: In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season. Results: The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation. Discussion: The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies.

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

Single nucleotide polymorphism profile for quantitative trait nucleotide in populations with small effective size and its impact on mapping and genomic predictions.

Increasing SNP density by incorporating sequence information only marginally increases prediction accuracies of breeding values in livestock. To find out why, we used statistical models and simulations to investigate the shape of distribution of estimated SNP effects (a profile) around quantitative trait nucleotides (QTNs) in populations with a small effective population size (Ne). A QTN profile created by averaging SNP effects around each QTN was similar to the shape of expected pairwise linkage disequilibrium (PLD) based on Ne and genetic distance between SNP, with a distinct peak for the QTN. Populations with smaller Ne showed lower but wider QTN profiles. However, adding more genotyped individuals with phenotypes dragged the profile closer to the QTN. The QTN profile was higher and narrower for populations with larger compared to smaller Ne. Assuming the PLD curve for the QTN profile, 80% of the additive genetic variance explained by each QTN was contained in ± 1/Ne Morgan interval around the QTN, corresponding to 2 Mb in cattle and 5 Mb in pigs and chickens. With such large intervals, identifying QTN is difficult even if all of them are in the data and the assumed genetic architecture is simplistic. Additional complexity in QTN detection arises from confounding of QTN profiles with signals due to relationships, overlapping profiles with closely spaced QTN, and spurious signals. However, small Ne allows for accurate predictions with large data even without QTN identification because QTNs are accounted for by QTN profiles if SNP density is sufficient to saturate the segments.

Association mapping and candidate gene identification for yield traits in European hazelnut (Corylus avellana L.).

European hazelnut (Corylus avellana L.) is an important nut crop due to its nutritional benefits, culinary uses, and economic value. Türkiye is the leading producer of hazelnut, followed by Italy and the United States. Quantitative trait locus studies offer promising opportunities for breeders and geneticists to identify genomic regions controlling desirable traits in hazelnut. A genome-wide association analysis was conducted with 5,567 single nucleotide polymorphisms on a Turkish core set of 86 hazelnut accessions, revealing 189 quantitative trait nucleotides (QTNs) associated with 22 of 31 traits (p < 2.9E-07). These QTNs were associated with plant and leaf, phenological, reproductive, nut, and kernel traits. Based on the close physical distance of QTNs associated with the same trait, we identified 23 quantitative trait loci. Furthermore, we identified 23 loci of multiple QTs comprising chromosome locations associated with more than one trait at the same position or in close proximity. A total of 159 candidate genes were identified for 189 QTNs, with 122 of them containing significant conserved protein domains. Some candidate matches to known proteins/domains were highly significant, suggesting that they have similar functions as their matches. This comprehensive study provides valuable insights for the development of breeding strategies and the improvement of hazelnut and enhances the understanding of the genetic architecture of complex traits by proposing candidate genes and potential functions.

Identification of QTNs, QTN-by-environment interactions, and their candidate genes for grain size traits in main crop and ratoon rice.

Although grain size is an important quantitative trait affecting rice yield and quality, there are few studies on gene-by-environment interactions (GEIs) in genome-wide association studies, especially, in main crop (MC) and ratoon rice (RR). To address these issues, the phenotypes for grain width (GW), grain length (GL), and thousand grain weight (TGW) of 159 accessions of MC and RR in two environments were used to associate with 2,017,495 SNPs for detecting quantitative trait nucleotides (QTNs) and QTN-by-environment interactions (QEIs) using 3VmrMLM. As a result, 64, 71, 67, 72, 63, and 56 QTNs, and 0, 1, 2, 2, 2, and 1 QEIs were found to be significantly associated with GW in MC (GW-MC), GL-MC, TGW-MC, GW-RR, GL-RR, and TGW-RR, respectively. 3, 4, 7, 2, 2, and 4 genes were found to be truly associated with the above traits, respectively, while 2 genes around the above QEIs were found to be truly associated with GL-RR, and one of the two known genes was differentially expressed under two soil moisture conditions. 10, 7, 1, 8, 4, and 3 candidate genes were found by differential expression and GO annotation analysis to be around the QTNs for the above traits, respectively, in which 6, 3, 1, 2, 0, and 2 candidate genes were found to be significant in haplotype analysis. The gene Os03g0737000 around one QEI for GL-MC was annotated as salt stress related gene and found to be differentially expressed in two cultivars with different grain sizes. Among all the candidate genes around the QTNs in this study, four were key, in which two were reported to be truly associated with seed development, and two (Os02g0626100 for GL-MC and Os02g0538000 for GW-MC) were new. Moreover, 1, 2, and 1 known genes, along with 8 additional candidate genes and 2 candidate GEIs, were found to be around QTNs and QEIs for GW, GL, and TGW, respectively in MC and RR joint analysis, in which 3 additional candidate genes were key and new. Our results provided a solid foundation for genetic improvement and molecular breeding in MC and RR.

Genome-Wide Association Studies Using 3VmrMLM Model Provide New Insights into Branched-Chain Amino Acid Contents in Rice Grains.

Rice (Oryza sativa L.) is a globally important food source providing carbohydrates, amino acids, and dietary fiber for humans and livestock. The branched-chain amino acid (BCAA) level is a complex trait related to the nutrient quality of rice. However, the genetic mechanism underlying the BCAA (valine, leucine, and isoleucine) accumulation in rice grains remains largely unclear. In this study, the grain BCAA contents and 239,055 SNPs of a diverse panel containing 422 rice accessions were adopted to perform a genome-wide association study (GWAS) using a recently proposed 3VmrMLM model. A total of 357 BCAA-content-associated main-effect quantitative trait nucleotides (QTNs) were identified from 15 datasets (12 BCAA content datasets and 3 BLUP datasets of BCAA). Furthermore, the allelic variation of two novel candidate genes, LOC_Os01g52530 and LOC_Os06g15420, responsible for the isoleucine (Ile) content alteration were identified. To reveal the genetic basis of the potential interactions between the gene and environmental factor, 53 QTN-by-environment interactions (QEIs) were detected using the 3VmrMLM model. The LOC_Os03g24460, LOC_Os01g55590, and LOC_Os12g31820 were considered as the candidate genes potentially contributing to the valine (Val), leucine (Leu), and isoleucine (Ile) accumulations, respectively. Additionally, 10 QTN-by-QTN interactions (QQIs) were detected using the 3VmrMLM model, which were putative gene-by-gene interactions related to the Leu and Ile contents. Taken together, these findings suggest that the implementation of the 3VmrMLM model in a GWAS may provide new insights into the deeper understanding of BCAA accumulation in rice grains. The identified QTNs/QEIs/QQIs serve as potential targets for the genetic improvement of rice with high BCAA levels.

Identification of QTNs, QTN-by-environment interactions for plant height and ear height in maize multi-environment GWAS.

Plant height (PH) and ear height (EH) are important traits associated with biomass, lodging resistance, and grain yield in maize. There were strong effects of genotype x environment interaction (GEI) on plant height and ear height of maize. In this study, 203 maize inbred lines were grown at five locations across China's Spring and Summer corn belts, and plant height (PH) and ear height (EH) phenotype data were collected and grouped using GGE biplot. Five locations fell into two distinct groups (or mega environments) that coincide with two corn ecological zones called Summer Corn Belt and Spring Corn Belt. In total, 73,174 SNPs collected using GBS sequencing platform were used as genotype data and a recently released multi-environment GWAS software package IIIVmrMLM was employed to identify QTNs and QTN x environment (corn belt) interaction (QEIs); 12 and 11 statistically significant QEIs for PH and EH were detected respectively and their phenotypic effects were further partitioned into Add*E and Dom*E components. There were 28 and 25 corn-belt-specific QTNs for PH and EH identified, respectively. The result shows that there are a large number of genetic loci underlying the PH and EH GEIs and IIIVmrMLM is a powerful tool in discovering QTNs that have significant QTN-by-Environment interaction. PH and EH candidate genes were annotated based on transcriptomic analysis and haplotype analysis. EH related-QEI S10_135 (Zm00001d025947, saur76, small auxin up RNA76) and PH related-QEI S4_4 (Zm00001d049692, mads32, encoding MADS-transcription factor 32), and corn-belt specific QTNs including S10_4 (Zm00001d023333, sdg127, set domain gene127) and S7_1 (Zm00001d018614, GLR3.4, and glutamate receptor 3.4 or Zm00001d018616, DDRGK domain-containing protein) were reported, and the relationship among GEIs, QEIs and phenotypic plasticity and their biological and breeding implications were discussed.

Association mapping reveals novel genes and genomic regions controlling grain size architecture in mini core accessions of Indian National Genebank wheat germplasm collection.

Wheat (Triticum aestivum L.) is a staple food crop for the global human population, and thus wheat breeders are consistently working to enhance its yield worldwide. In this study, we utilized a sub-set of Indian wheat mini core germplasm to underpin the genetic architecture for seed shape-associated traits. The wheat mini core subset (125 accessions) was genotyped using 35K SNP array and evaluated for grain shape traits such as grain length (GL), grain width (GW), grain length, width ratio (GLWR), and thousand grain weight (TGW) across the seven different environments (E1, E2, E3, E4, E5, E5, E6, and E7). Marker-trait associations were determined using a multi-locus random-SNP-effect Mixed Linear Model (mrMLM) program. A total of 160 non-redundant quantitative trait nucleotides (QTNs) were identified for four grain shape traits using two or more GWAS models. Among these 160 QTNs, 27, 36, 38, and 35 QTNs were associated for GL, GW, GLWR, and TGW respectively while 24 QTNs were associated with more than one trait. Of these 160 QTNs, 73 were detected in two or more environments and were considered reliable QTLs for the respective traits. A total of 135 associated QTNs were annotated and located within the genes, including ABC transporter, Cytochrome450, Thioredoxin_M-type, and hypothetical proteins. Furthermore, the expression pattern of annotated QTNs demonstrated that only 122 were differentially expressed, suggesting these could potentially be related to seed development. The genomic regions/candidate genes for grain size traits identified in the present study represent valuable genomic resources that can potentially be utilized in the markers-assisted breeding programs to develop high-yielding varieties.

Multi-locus genome-wide association study of fusarium head blight in relation to days to anthesis and plant height in a spring wheat association panel.

Fusarium head blight (FHB) is a highly destructive fungal disease of wheat to which host resistance is quantitatively inherited and largely influenced by the environment. Resistance to FHB has been associated with taller height and later maturity; however, a further understanding of these relationships is needed. An association mapping panel (AMP) of 192 predominantly Canadian spring wheat was genotyped with the wheat 90K single-nucleotide polymorphism (SNP) array. The AMP was assessed for FHB incidence (INC), severity (SEV) and index (IND), days to anthesis (DTA), and plant height (PLHT) between 2015 and 2017 at three Canadian FHB-inoculated nurseries. Seven multi-environment trial (MET) datasets were deployed in a genome-wide association study (GWAS) using a single-locus mixed linear model (MLM) and a multi-locus random SNP-effect mixed linear model (mrMLM). MLM detected four quantitative trait nucleotides (QTNs) for INC on chromosomes 2D and 3D and for SEV and IND on chromosome 3B. Further, mrMLM identified 291 QTNs: 50 (INC), 72 (SEV), 90 (IND), 41 (DTA), and 38 (PLHT). At two or more environments, 17 QTNs for FHB, DTA, and PLHT were detected. Of these 17, 12 QTNs were pleiotropic for FHB traits, DTA, and PLHT on chromosomes 1A, 1D, 2D, 3B, 5A, 6B, 7A, and 7B; two QTNs for DTA were detected on chromosomes 1B and 7A; and three PLHT QTNs were located on chromosomes 4B and 6B. The 1B DTA QTN and the three pleiotropic QTNs on chromosomes 1A, 3B, and 6B are potentially identical to corresponding quantitative trait loci (QTLs) in durum wheat. Further, the 3B pleiotropic QTN for FHB INC, SEV, and IND co-locates with TraesCS3B02G024900 within the Fhb1 region on chromosome 3B and is ~3 Mb from a cloned Fhb1 candidate gene TaHRC. While the PLHT QTN on chromosome 6B is putatively novel, the 1B DTA QTN co-locates with a disease resistance protein located ~10 Mb from a Flowering Locus T1-like gene TaFT3-B1, and the 7A DTA QTN is ~5 Mb away from a maturity QTL QMat.dms-7A.3 of another study. GWAS and QTN candidate genes enabled the characterization of FHB resistance in relation to DTA and PLHT. This approach should eventually generate additional and reliable trait-specific markers for breeding selection, in addition to providing useful information for FHB trait discovery.