Association Mapping of Flowering Time QTLs and Insight into Their Contributions to Rapeseed Growth Habits

Discovery of common loci and candidate genes for controlling salt-alkali tolerance and yield-related traits in Brassica napus L

KEY MESSAGE

Common loci and candidate genes for controlling salt-alkali tolerance and yield-related traits were identified in Brassica napus combining QTL mapping with transcriptome under salt and alkaline stresses. The yield of rapeseed (Brassica napus L.) is determined by multiple yield-related traits, which are susceptible to environmental factors. Many yield-related quantitative trait loci (QTLs) have been reported in Brassica napus; however, no studies have been conducted to investigate both salt-alkali tolerance and yield-related traits simultaneously. Here, specific-locus amplified fragment sequencing (SLAF-seq) technologies were utilized to map the QTLs for salt-alkali tolerance and yield-related traits. A total of 65 QTLs were identified, including 30 QTLs for salt-alkali tolerance traits and 35 QTLs for yield-related traits, accounting for 7.61-27.84% of the total phenotypic variations. Among these QTLs, 18 unique QTLs controlling two to four traits were identified by meta-analysis. Six novel and unique QTLs were detected for salt-alkali tolerance traits. By comparing these unique QTLs for salt-alkali tolerance traits with those previously reported QTLs for yield-related traits, seven co-localized chromosomal regions were identified on A09 and A10. Combining QTL mapping with transcriptome of two parents under salt and alkaline stresses, thirteen genes were identified as the candidates controlling both salt-alkali tolerance and yield. These findings provide useful information for future breeding of high-yield cultivars resistant to alkaline and salt stresses.

Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica napus

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Genome-Wide Association Mapping Unravels the Genetic Control of Seed Vigor under Low-Temperature Conditions in Rapeseed (Brassica napus L.)

Low temperature inhibits rapid germination and successful seedling establishment of rapeseed (Brassica napus L.), leading to significant productivity losses. Little is known about the genetic diversity for seed vigor under low-temperature conditions in rapeseed, which motivated our investigation of 13 seed germination- and emergence-related traits under normal and low-temperature conditions for 442 diverse rapeseed accessions. The stress tolerance index was calculated for each trait based on performance under non-stress and low-temperature stress conditions. Principal component analysis of the low-temperature stress tolerance indices identified five principal components that captured 100% of the seedling response to low temperature. A genome-wide association study using ~8 million SNP (single-nucleotide polymorphism) markers identified from genome resequencing was undertaken to uncover the genetic basis of seed vigor related traits in rapeseed. We detected 22 quantitative trait loci (QTLs) significantly associated with stress tolerance indices regarding seed vigor under low-temperature stress. Scrutiny of the genes in these QTL regions identified 62 candidate genes related to specific stress tolerance indices of seed vigor, and the majority were involved in DNA repair, RNA translation, mitochondrial activation and energy generation, ubiquitination and degradation of protein reserve, antioxidant system, and plant hormone and signal transduction. The high effect variation and haplotype-based effect of these candidate genes were evaluated, and high priority could be given to the candidate genes BnaA03g40290D, BnaA06g07530D, BnaA09g06240D, BnaA09g06250D, and BnaC02g10720D in further study. These findings should be useful for marker-assisted breeding and genomic selection of rapeseed to increase seed vigor under low-temperature stress.

Genome-Wide Association Study for Grain Micronutrient Concentrations in Wheat Advanced Lines Derived From Wild Emmer

Wheat is one of the important staple crops as the resources of both food and micronutrient for most people of the world. However, the levels of micronutrients (especially Fe and Zn) in common wheat are inherently low. Biofortification is an effective way to increase the micronutrient concentration of wheat. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, AABB, 2n = 4x = 28) is an important germplasm resource for wheat micronutrients improvement. In the present study, a genome-wide association study (GWAS) was performed to characterize grain iron, zinc, and manganese concentration (GFeC, GZnC, and GMnC) in 161 advanced lines derived from wild emmer. Using both the general linear model and mixed linear model, we identified 14 high-confidence significant marker-trait associations (MTAs) that were associated with GFeC, GZnC, and GMnC of which nine MTAs were novel. Six MTAs distributed on chromosomes 3B, 4A, 4B, 5A, and 7B were significantly associated with GFeC. Three MTAs on 1A and 2A were significantly associated with GZnC and five MTAs on 1B were significantly associated with GMnC. These MTAs show no negative effects on thousand kernel weight (TKW), implying the potential value for simultaneous improvement of micronutrient concentrations and TKW in breeding. Meanwhile, the GFeC, GZnC and GMnC are positively correlated, suggesting that these traits could be simultaneously improved. Genotypes containing high-confidence MTAs and 61 top genotypes with a higher concentration of grain micronutrients were recommended for wheat biofortification breeding. A total of 38 candidate genes related to micronutrient concentrations were identified. These candidates can be classified into four main groups: enzymes, transporter proteins, MYB transcription factor, and plant defense responses proteins. The MTAs and associated candidate genes provide essential information for wheat biofortification breeding through marker-assisted selection (MAS).

Genome-Wide Association Mapping for Heat Stress Responsive Traits in Field Pea

Environmental stress hampers pea productivity. To understand the genetic basis of heat resistance, a genome-wide association study (GWAS) was conducted on six stress responsive traits of physiological and agronomic importance in pea, with an objective to identify the genetic loci associated with these traits. One hundred and thirty-five genetically diverse pea accessions from major pea growing areas of the world were phenotyped in field trials across five environments, under generally ambient (control) and heat stress conditions. Statistical analysis of phenotype indicated significant effects of genotype (G), environment (E), and G × E interaction for all traits. A total of 16,877 known high-quality SNPs were used for association analysis to determine marker-trait associations (MTA). We identified 32 MTAs that were consistent in at least three environments for association with the traits of stress resistance: six for chlorophyll concentration measured by a soil plant analysis development meter; two each for photochemical reflectance index and canopy temperature; seven for reproductive stem length; six for internode length; and nine for pod number. Forty-eight candidate genes were identified within 15 kb distance of these markers. The identified markers and candidate genes have potential for marker-assisted selection towards the development of heat resistant pea cultivars.

Genome-wide association study (GWAS) reveals genetic loci of lead (Pb) tolerance during seedling establishment in rapeseed (Brassica napus L.)

BACKGROUND

Lead (Pb) pollution in soil has become one of the major environmental threats to plant growth and human health. Safe utilization of Pb contaminated soil by phytoremediation require Pb-tolerant rapeseed (Brassica napus L.) accessions. However, breeding of new B. napus cultivars tolerance to Pb stress has been restricted by limited knowledge on molecular mechanisms involved in Pb tolerance. This work was carried out to identify genetic loci related to Pb tolerance during seedling establishment in rapeseed.

RESULTS

Pb tolerance, which was assessed by quantifying radicle length (RL) under 0 or 100 mg/L Pb stress condition, shown an extensive variation in 472 worldwide-collected rapeseed accessions. Based on the criterion of relative RL > 80%, six Pb-tolerant genotypes were selected. Four quantitative trait loci (QTLs) associated with Pb tolerance were identified by Genome-wide association study. The expression level of nine promising candidate genes, including GSTUs, BCATs, UBP13, TBR and HIPP01, located in these four QTL regions, were significantly higher or induced by Pb in Pb-tolerant accessions in comparison to Pb-sensitive accessions.

CONCLUSION

To our knowledge, this is the first study on Pb-tolerant germplasms and genomic loci in B. napus. The findings can provide valuable genetic resources for the breeding of Pb-tolerant B. napus cultivars and understanding of Pb tolerance mechanism in Brassica species.

Mapping of quantitative trait loci for traits linked to fusarium head blight in barley

Fusarium head blight (FHB) is a devastating disease occurring in small grain cereals worldwide. The disease results in the reduction of grain yield, and mycotoxins accumulated in grain are also harmful to both humans and animals. It has been reported that response to pathogen infection may be associated with the morphological and developmental traits of the host plant, e.g. earliness and plant height. Despite many studies, effective markers for selection of barley genotypes with increased resistance to FHB have not been developed. In the present study, we investigated 100 recombinant inbred lines (RIL) of spring barley. Plants were examined in field conditions (three locations) in a completely randomized design with three replications. Barley genotypes were artificially infected with spores of Fusarium culmorum before heading. Apart from the main phenotypic traits (plant height, spike characteristic, grain yield), infected kernels were visually scored and the content of deoxynivalenol (DON) mycotoxin was investigated. A set of 70 Quantitative Trait Loci (QTLs) were detected through phenotyping of the mapping population in field conditions and genotyping using a barley Ilumina 9K iSelect platform. Six loci were detected for the FHB index on chromosomes 2H, 3H, 5H, and 7H. A region on the short arm of chromosome 2H was detected in which many QTLs associated with FHB- and yield-related traits were found. This study confirms that agromorphological traits are tightly related to FHB and should be taken into consideration when breeding barley plants for FHB resistance.

Strong temporal dynamics of QTL action on plant growth progression revealed through high-throughput phenotyping in canola

A major challenge of plant biology is to unravel the genetic basis of complex traits. We took advantage of recent technical advances in high-throughput phenotyping in conjunction with genome-wide association studies to elucidate genotype-phenotype relationships at high temporal resolution. A diverse Brassica napus population from a commercial breeding programme was analysed by automated non-invasive phenotyping. Time-resolved data for early growth-related traits, including estimated biovolume, projected leaf area, early plant height and colour uniformity, were established and complemented by fresh and dry weight biomass. Genome-wide SNP array data provided the framework for genome-wide association analyses. Using time point data and relative growth rates, multiple robust main effect marker-trait associations for biomass and related traits were detected. Candidate genes involved in meristem development, cell wall modification and transcriptional regulation were detected. Our results demonstrate that early plant growth is a highly complex trait governed by several medium and many small effect loci, most of which act only during short phases. These observations highlight the importance of taking the temporal patterns of QTL/allele actions into account and emphasize the need for detailed time-resolved analyses to effectively unravel the complex and stage-specific contributions of genes affecting growth processes that operate at different developmental phases.

Integrating GWAS, linkage mapping and gene expression analyses reveals the genetic control of growth period traits in rapeseed (Brassica napus L.)

BACKGROUND

Brassica napus is one of the most important oilseed crops, and also an important biofuel plant due to its low air pollution and renewability. Growth period are important traits that affect yield and are crucial for its adaptation to different environments in B. napus.

RESULTS

To elucidate the genetic basis of growth period traits, genome-wide association analysis (GWAS) and linkage mapping were employed to detect the quantitative trait loci (QTL) for days to initial flowering (DIF), days to final flowering (DFF), flowering period (FP), maturity time (MT), and whole growth period (GP). A total of 146 SNPs were identified by association mapping, and 83 QTLs were identified by linkage mapping using the RIL population. Among these QTLs, 19 were pleiotropic SNPs related to multiple traits, and six (q18DFF.A03-2, q18MT.A03-2, q17DFF.A05-1, q18FP.C04, q17DIF.C05 and q17GP.C09) were consistently detected using both mapping methods. Additionally, we performed RNA sequencing to analyze the differential expression of gene (DEG) transcripts between early- and late-flowering lines selected from the RIL population, and the DEGs were integrated with association mapping and linkage analysis to confirm their roles in the growth period. Consequently, 12 candidate genes associated with growth period traits were identified in B. napus. Among these genes, seven have polymorphic sites in the coding sequence and the upstream 2-kb sequence based on the resequencing data. The haplotype BnaSOC1.A05-Haplb and BnaLNK2.C06-Hapla showed more favorable phenotypic traits.

CONCLUSIONS

The candidate genes identified in this study will contribute to our genetic understanding of growth period traits and can be used as targets for target mutations or marker-assisted breeding for rapeseed adapted to different environments.

Identification of SNP loci and candidate genes related to four important fatty acid composition in Brassica napus using genome wide association study

Rapeseed oil (canola, Brassica napus L.) is an important healthy vegetable oil throughout the world, the nutritional and economical value of which largely depends on its seed fatty acid composition. In this study, based on 201,187 SNP markers developed from the SLAF-seq (specific locus amplified fragment sequencing), a genome wide association study of four important fatty acid content traits (erucic acid, oleic acid, linoleic acid and linolenic acid) in a panel of 300 inbred lines of rapeseed in two environments (JXAU and JXRIS) was carried out. A total of 148 SNP loci significantly associated with these traits were detected by MLM model analysis respectively, and 30 SNP loci on A08 and C03 chromosomes were detected in three traits of erucic acid, oleic acid and linoleic acid contents simultaneously. Furthermore, 108 highly favorable alleles for increasing oleic acid and linoleic acid content, also for decreasing erucic acid content simultaneously were observed. By a basic local alignment search tool (BLAST) search with in a distance of 100 Kb around these significantly SNP-trait associations, we identified 20 orthologs of the functional candidate genes related to fatty acid biosynthesis, including the known vital fatty acid biosynthesis genes of BnaA.FAE1 and BnaC. FAE1 on the A08 and C03 chromosomes, and other potential candidate genes involving in the fatty acid biosynthesis pathway, such as the orthologs genes of FAD2, LACS09, KCS17, CER4, TT16 and ACBP5. This study lays a basis for uncovering the genetic variations and the improvement of fatty acid composition in B. napus.

Identification and application of markers closely linked to the restorer gene (Rfm) in rapeseed (Brassica napus L.)

The Mutsu-Isuzu cytoplasmic male sterility (MI CMS) system is one of the three-line hybrid systems used in China. As we know, the hybrid system is tightly associated with the yield variation in F1 heterosis, while the restorer gene for the MI CMS (Rfm) has not been finely mapped for further application in marker-assisted selection (MAS). In this study, the sets of near-isogenic lines (NILs) of Rfm in two different genetic backgrounds were hybridized with the genome-wide 60 K single-nucleotide polymorphism (SNP) chip of Brassica for screening the possible associated genomic region of Rfm. Through screening genotypes with SNP loci and sequencing the candidate loci, one 2.5 Mb physical region (covering three scaffolds) on chrA09 was identified as the candidate for the Rfm region. Then, the SSR markers for the target scaffolds were used to detect the recombination in an F2 population and narrowed the Rfm gene within the genetic distance of 0.52 cM, equivalent to a 350 kb physical segment. Moreover, the markers were tested to improve new elite restoration lines and to assess the percentage of hybrid seeds. Our results could potentially accelerate the map-based cloning of the Rfm gene to benefit rapeseed breeding.

Genome-Wide Association Study Reveals Novel Genomic Regions Associated With High Grain Protein Content in Wheat Lines Derived From Wild Emmer Wheat

Grain protein content (GPC) and yield are of two important traits in wheat, but their negative correlation has hampered their simultaneous improvement in conventional breeding. Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is an important genetic resource for wheat quality improvement. In this study, we report a genome-wide association study (GWAS) using 13116 DArT-seq markers to characterize GPC in 161 wheat lines derived from wild emmer. Using a general linear model, we identified 141 markers that were significantly associated with GPC, and grouped into 48 QTL regions. Using both general linear model and mixed linear model, we identified four significant markers that were grouped into two novel QTL regions on chromosomes 2BS (QGpc.cd1-2B.1) and 7BL (QGpc.cd1-7B.2). The two QTLs have no negative effects on thousand kernel weight (TKW) and should be useful for simultaneous improvement of GPC and TKW in wheat breeding. Searches of public databases revealed 61 putative candidate/flanking genes related to GPC. The putative proteins of interest were grouped in four main categories: enzymes, kinase proteins, metal transport-related proteins, and disease resistance proteins. The linked markers and associated candidate genes provide essential information for cloning genes related to high GPC and performing marker-assisted breeding in wheat.

Translating Flowering Time From Arabidopsis thaliana to Brassicaceae and Asteraceae Crop Species

Flowering and seed set are essential for plant species to survive, hence plants need to adapt to highly variable environments to flower in the most favorable conditions. Endogenous cues such as plant age and hormones coordinate with the environmental cues like temperature and day length to determine optimal time for the transition from vegetative to reproductive growth. In a breeding context, controlling flowering time would help to speed up the production of new hybrids and produce high yield throughout the year. The flowering time genetic network is extensively studied in the plant model species Arabidopsis thaliana, however this knowledge is still limited in most crops. This article reviews evidence of conservation and divergence of flowering time regulation in A. thaliana with its related crop species in the Brassicaceae and with more distant vegetable crops within the Asteraceae family. Despite the overall conservation of most flowering time pathways in these families, many genes controlling this trait remain elusive, and the function of most Arabidopsis homologs in these crops are yet to be determined. However, the knowledge gathered so far in both model and crop species can be already exploited in vegetable crop breeding for flowering time control.

Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar

High-quality and dense genetic maps were constructed, and leaf shape variation was dissected by QTL mapping in poplar. Species in the genus Populus, also known as poplars, are important woody species and considered model plants for perennial trees. High-density genetic maps are valuable genomic resources for population genetics. Here, we generated a high-quality and dense genetic map for an F1 poplar population using high-throughput NGS-based genotyping. A total of 92,097 high-quality SNP markers were developed by stringent filtering and identification. In total, 889 and 1650 SNPs formed the female and male genetic maps, respectively. To test the application of the genetic maps, QTL mapping of leaf shape was conducted for this F1 population. A total of nine parameters were scored for leaf shape variation in three different environments. Combining genetic maps and measurements of the nine leaf shape parameters, we mapped a total of 42 significant QTLs. The highest LOD score of all QTLs was 9.2, and that QTL explained the most (15.13%) trait variation. A total of nine QTLs could be detected in at least two environments, and they were located in two genomic regions. Within these two QTL regions, some candidate genes for regulating leaf shape were predicted through functional annotation. The successful mapping of leaf shape QTLs demonstrated the utility of our genetic maps. According to the performance of this study, we were able to provide high-quality and dense genetic maps and dissect the leaf shape variation in poplar.

A comprehensive and precise set of intervarietal substitution lines to identify candidate genes and quantitative trait loci in oilseed rape (Brassica napus L.)

A set of intervarietal substitution lines were developed in rapeseed by recurrent backcrossing and marker-assisted selection and employed for mapping both qualitative and quantitative traits. Intervarietal substitution lines (ISLs) may be assembled into advanced secondary mapping populations that have remarkable potential for resolving trait loci and mapping candidate genes. To facilitate the identification of important genes in oilseed rape (canola, Brassica napus), we developed 89 ISLs using an elite cultivar 'Zhongyou 821' (ZY821) as the recipient and a re-synthesized line 'No.2127' as the donor. In the whole process of ISLs development, the target chromosome segments were selected based on the genotypes of 300 microsatellite markers evenly distributed across the genome. Eighty-nine ISLs fixed at BC₅F₄ were genotyped by sequencing using double digestion to survey the lengths of target substitution segments from the donor parent and the background segments from the recurrent parent. The total length of the substituted chromosome segments was 3030.27 Mb, representing 3.56 × of the Darmor-bzh reference genome sequence (version 4.1). Gene mapping was conducted for two qualitative traits, flower colour and seed-coat colour, and nine quantitative traits including yield- and quality-related traits, with 19 QTLs identified for the latter. Overlapping substitution segments were identified for flower colour and seed-coat colour loci, as well as for QTLs consistently detected in 2 or 3 years. These results demonstrate the value of these ISLs for locus resolution and subsequent cloning, targeted mutation or editing of genes controlling important traits in oilseed rape.

Sequence variation and functional analysis of a FRIGIDA orthologue (BnaA3.FRI) in Brassica napus

BACKGROUND

Allelic variation at the FRIGIDA (FRI) locus is a major contributor to natural variation of flowering time and vernalization requirement in Arabidopsis thaliana. Dominant FRI inhibits flowering by activating the expression of the MADS box transcriptional repressor FLOWERING LOCUS C (FLC), which represses flowering prior to vernalization. Four FRI orthologues had been identified in the domesticated amphidiploid Brassica napus. Linkage and association studies had revealed that one of the FRI orthologues, BnaA3.FRI, contributes to flowering time variation and crop type differentiation.

RESULTS

Sequence analyses indicated that three out of the four BnaFRI paralogues, BnaA3.FRI, BnaA10.FRI and BnaC3.FRI, contained a large number of polymorphic sites. Haplotype analysis in a panel of 174 B. napus accessions using PCR markers showed that all the three paralogues had a biased distribution of haplotypes in winter type oilseed rape (P < 0.01). Association analysis indicated that only BnaA3.FRI contributes to flowering time variation in B. napus. In addition, transgenic functional complementation demonstrated that mutations in the coding sequence of BnaA3.FRI lead to weak alleles, and subsequently to flowering time variation.

CONCLUSION

This study for the first time provides a molecular basis for flowering time control by BnaA3.FRI in B. napus, and will facilitate predictive oilseed rape breeding to select varieties with favorable flowering time and better adaption to latitude and seasonal shifts due to changing climate.

QTL Alignment for Seed Yield and Yield Related Traits in Brassica napus

Worldwide consumption of oil is increasing with the growing population in need for edible oil and the expansion of industry using biofuels. Then, demand for high yielding varieties of oil crops is always increasing. Brassica napus (rapeseed) is one of the most important oil crop in the world, therefore, increasing rapeseed yield through breeding is inevitable in order to cater for the high demand of vegetable oil and high-quality protein for live stocks. Quantitative trait loci (QTL) analysis is a powerful tool to identify important loci and which is also valuable for molecular marker assisted breeding. Seed-yield (SY) is a complex trait that is controlled by multiple loci and is affected directly by seed weight, seeds per silique and silique number. Some yield-related traits, such as plant height, biomass yield, flowering time, and so on, also affect the SY indirectly. This study reports the assembly of QTLs identified for seed-yield and yield-related traits in rapeseed, in one unique map. A total of 972 QTLs for seed-yield and yield-related were aligned into the physical map of B. napus Darmor-bzh and 92 regions where 198 QTLs overlapped, could be discovered on 16 chromosomes. Also, 147 potential candidate genes were discovered in 65 regions where 131 QTLs overlapped, and might affect nine different traits. At the end, interaction network of candidate genes was studied, and showed nine genes that could highly interact with the other genes, and might have more influence on them. The present results would be helpful to develop molecular markers for yield associated traits and could be used for breeding improvement in B. napus.

Major Co-localized QTL for Plant Height, Branch Initiation Height, Stem Diameter, and Flowering Time in an Alien Introgression Derived Brassica napus DH Population

Plant height (PH), branch initiation height (BIH), and stem diameter (SD) are three stem-related traits that play crucial roles in plant architecture and lodging resistance. Herein, we show one doubled haploid (DH) population obtained from a cross between Y689 (one Capsella bursa-pastoris derived Brassica napus intertribal introgression) and Westar (B. napus cultivar) that these traits were significantly positively correlated with one another and with flowering time (FT). Based on a high-density SNP map, a total of 102 additive quantitative trait loci (QTL) were identified across six environments. Seventy-two consensus QTL and 49 unique QTL were identified using a two-round strategy of QTL meta-analysis. Notably, a total of 19 major QTL, including 11 novel ones, were detected for these traits, which comprised two QTL clusters on chromosomes A02 and A07. Conditional QTL mapping was performed to preliminarily evaluate the genetic basis (pleiotropy or tight linkage) of the co-localized QTL. In addition, QTL by environment interactions (QEI) mapping was performed to verify the additive QTL and estimate the QEI effect. In the genomic regions of all major QTL, orthologs of the genes involved in phytohormone biosynthesis, phytohormone signaling, flower development, and cell differentiation in Arabidopsis were proposed as candidate genes. Of these, BnaA02g02560, an ortholog of Arabidopsis GASA4, was suggested as a candidate gene for PH, SD, and FT; and BnaA02g08490, an ortholog of Arabidopsis GNL, was associated with PH, BIH and FT. These results provide useful information for further genetic studies on stem-related traits and plant growth adaptation.

An Integration of Genome-Wide Association Study and Gene Co-expression Network Analysis Identifies Candidate Genes of Stem Lodging-Related Traits in Brassica napus

Lodging is a persistent problem which severely reduce yield and impair seed quality in rapeseed (Brassica napus L.). Enhancing stem strength (SS) has proven to be an effective approach to decrease lodging risk. In the present study, four interrelated stem lodging-related traits, including stem breaking resistance (SBR), stem diameter (SD), SS, and lodging coefficient (LC), were investigated among 472 rapeseed accessions. A genome-wide association study (GWAS) using Brassica 60K SNP array for stem lodging-related traits identified 67 significantly associated quantitative trait loci (QTLs) and 71 candidate genes. In parallel, a gene co-expression network based on transcriptome sequencing was constructed. The module associated with cellulose biosynthesis was highlighted. By integrating GWAS and gene co-expression network analysis, some promising candidate genes, such as ESKIMO1 (ESK1, BnaC08g26920D), CELLULOSE SYNTHASE 6 (CESA6, BnaA09g06990D), and FRAGILE FIBER 8 (FRA8, BnaC04g39510D), were prioritized for further research. These findings revealed the genetic basis underlying stem lodging and provided worthwhile QTLs and genes information for genetic improvement of stem lodging resistance in B. napus.

Genome-wide association study reveals genetic loci and candidate genes for average daily gain in Duroc pigs

Objective

Average daily gain (ADG) is an important target trait of pig breeding programs. We aimed to identify single nucleotide polymorphisms (SNPs) and genomic regions that are associated with ADG in the Duroc pig population.

Methods

We performed a genome-wide association study involving 390 Duroc boars and by using the PorcineSNP60K Beadchip and two linear models.

Results

After quality control, we detected 3,5971 SNPs, which included seven SNPs that are significantly associated with the ADG of pigs. We identified six quantitative trait loci (QTL) regions for ADG. These QTLs included four previously reported QTLs on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, and SSC13, as well as two novel QTLs on SSC6 and SSC16. In addition, we selected six candidate genes (general transcription factor 3C polypeptide 5, high mobility group AT-hook 2, nicotinamide phosphoribosyltransferase, oligodendrocyte transcription factor 1, pleckstrin homology and RhoGEF domain containing G4B, and ENSSSCG00000031548) associated with ADG on the basis of their physiological roles and positional information. These candidate genes are involved in skeletal muscle cell differentiation, diet-induced obesity, and nervous system development.

Conclusion

This study contributes to the identification of the casual mutation that underlies QTLs associated with ADG and to future pig breeding programs based on marker-assisted selection. Further studies are needed to elucidate the role of the identified candidate genes in the physiological processes involved in ADG regulation.

Genomic Prediction of Genotypic Effects with Epistasis and Environment Interactions for Yield-Related Traits of Rapeseed (Brassica napus L.)

Oilseed rape (Brassica napus) is an economically important oil crop, yet the genetic architecture of its complex traits remain largely unknown. Here, genome-wide association study was conducted for eight yield-related traits to dissect the genetic architecture of additive, dominance, epistasis, and their environment interaction. Additionally, the optimal genotype combination and the breeding value of superior line, superior hybrid and existing best line in mapping population were predicted for each trait in two environments based on the predicted genotypic effects. As a result, 17 quantitative trait SNPs (QTSs) were identified significantly for target traits with total heritability varied from 58.47 to 87.98%, most of which were contributed by dominance, epistasis, and environment-specific effects. The results indicated that non-additive effects were large contributions to heritability and epistasis, and also noted that environment interactions were important variants for oilseed breeding. Our study facilitates the understanding of genetic basis of rapeseed yield trait, helps to accelerate rapeseed breading, and also offers a roadmap for precision plant breeding via marker-assisted selection.

Joint genome-wide association and transcriptome sequencing reveals a complex polygenic network underlying hypocotyl elongation in rapeseed (Brassica napus L.)

Hypocotyl elongation is considered an important typical seedling trait contributing directly to an increase in and stabilization of the yield in Brassica napus, but its molecular genetic mechanism is poorly understood. In the present study, hypocotyl lengths of 210 lines were measured in an illuminated culture room. A genome-wide association study (GWAS) was performed with 23,435 single nucleotide polymorphisms (SNPs) for hypocotyl length. Three lines with long hypocotyl length and three lines with short hypocotyl length from one doubled haploid line (DH) population were used for transcriptome sequencing. A GWAS followed by transcriptome analysis identified 29 differentially expressed genes associated with significant SNPs in B. napus. These genes regulate hypocotyl elongation by mediating flowering morphogenesis, circadian clock, hormone biosynthesis, or important metabolic signaling pathways. Among these genes, BnaC07g46770D negatively regulates hypocotyl elongation directly, as well as flowering time. Our results indicate that a joint GWAS and transcriptome analysis has significant potential for identifying the genes responsible for hypocotyl elongation; The extension of hypocotyl is a complex biological process regulated by a polygenic network.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

Genome-Wide Association Study Reveals Candidate Genes for Control of Plant Height, Branch Initiation Height and Branch Number in Rapeseed (Brassica napus L.)

Plant architecture is crucial for rapeseed yield and is determined by plant height (PH), branch initiation height (BIH), branch number (BN) and leaf and inflorescence morphology. In this study, we measured three major factors (PH, BIH, and BN) in a panel of 333 rapeseed accessions across 4 years. A genome-wide association study (GWAS) was performed via Q + K model and the panel was genotyped using the 60 k Brassica Infinium SNP array. We identified seven loci for PH, four for BIH, and five for BN. Subsequently, by determining linkage disequilibrium (LD) decay associated with 38 significant SNPs, we gained 31, 15, and 17 candidate genes for these traits, respectively. We also showed that PH is significantly correlated with BIH, while no other correlation was revealed. Notably, a GA signaling gene (BnRGA) and a flowering gene (BnFT) located on chromosome A02 were identified as the most likely candidate genes associated with PH regulation. Furthermore, a meristem initiation gene (BnLOF2) and a NAC domain transcriptional factor (BnCUC3) that may be associated with BN were identified on the chromosome A07. This study reveals novel insight into the genetic control of plant architecture and may facilitate marker-based breeding for rapeseed.