Novel and major QTL for branch angle detected by using DH population from an exotic introgression in rapeseed (Brassica napus L.)

Identification of a major QTL and candidate genes analysis for branch angle in rapeseed (Brassica napus L.) using QTL-seq and RNA-seq

Introduction

Branching angle is an essential trait in determining the planting density of rapeseed (Brassica napus L.) and hence the yield per unit area. However, the mechanism of branching angle formation in rapeseed is not well understood.

Methods

In this study, two rapeseed germplasm with extreme branching angles were used to construct an F2 segregating population; then bulked segregant analysis sequencing (BSA-seq) and quantitative trait loci (QTL) mapping were utilized to localize branching anglerelated loci and combined with transcriptome sequencing (RNA-seq) and quantitative real-time PCR (qPCR) for candidate gene mining.

Results and discussion

A branching angle-associated quantitative trait loci (QTL) was mapped on chromosome C3 (C3: 1.54-2.65 Mb) by combining BSA-seq as well as traditional QTL mapping. A total of 54 genes had SNP/Indel variants within the QTL interval were identified. Further, RNA-seq of the two parents revealed that 12 of the 54 genes were differentially expressed between the two parents. Finally, we further validated the differentially expressed genes using qPCR and found that six of them presented consistent differential expression in all small branching angle samples and large branching angles, and thus were considered as candidate genes related to branching angles in rapeseed. Our results introduce new candidate genes for the regulation of branching angle formation in rapeseed, and provide an important reference for the subsequent exploration of its formation mechanism.

3D genome structural variations play important roles in regulating seed oil content of Brassica napus

Dissecting the complex regulatory mechanism of seed oil content (SOC) is one of the main research goals in Brassica napus. Increasing evidence suggests that genome architecture is linked to multiple biological functions. However, the effect of genome architecture on SOC regulation remains unclear. Here, we used high-throughput chromatin conformation capture to characterize differences in the three-dimensional (3D) landscape of genome architecture of seeds from two B. napus lines, N53-2 (with high SOC) and Ken-C8 (with low SOC). Bioinformatics analysis demonstrated that differentially accessible regions and differentially expressed genes between N53-2 and Ken-C8 were preferentially enriched in regions with quantitative trait loci (QTLs)/associated genomic regions (AGRs) for SOC. A multi-omics analysis demonstrated that expression of SOC-related genes was tightly correlated with genome structural variations in QTLs/AGRs of B. napus. The candidate gene BnaA09g48250D, which showed structural variation in a QTL/AGR on chrA09, was identified by fine-mapping of a KN double-haploid population derived from hybridization of N53-2 and Ken-C8. Overexpression and knockout of BnaA09g48250D led to significant increases and decreases in SOC, respectively, in the transgenic lines. Taken together, our results reveal the 3D genome architecture of B. napus seeds and the roles of genome structural variations in SOC regulation, enriching our understanding of the molecular mechanisms of SOC regulation from the perspective of spatial chromatin structure.

High-Throughput Association Mapping in Brassica napus L.: Methods and Applications

Oil seed rape (Braasica napus L.) is ranked second among oil seed crops cultivated globally for edible oil for human, and seed cake for animal consumption. Recent genetic and genomics advancements highlighted the diversity that exists within B. napus, which is largely discovered using the most promising genetic markers called single nucleotide polymorphism (SNP). Their calling rate is also enhanced to ~100 folds after the continuous advancements in the next generation sequencing (NGS) technologies. As the high throughput of NGS resulted in multi-Giga bases data, the detailed quality control (QC) prior to downstream analyses is a pre-requisite. It mainly involved the removal of false positives, missing proportions, filtering of low-quality SNPs, and adjustments of minor-allele frequency and heterozygosity. After marker-trait association, for conformation of target SNPs, validations of SNPs can be performed using various methods, especially allele-specific PCR assay-based methods have been utilized for SNP genotyping of genes targeting agronomic traits and somaclonal variations occurred during transgenic studies. In the present study, the authors mainly argue on the genotypic progress, and pipelines/methods that are being used for detection, calling, filtering, and validation of SNPs. Also, insight is provided into the application of SNPs in linkage and association mapping, including QTL mapping and genome-wide association studies targeting mainly developmental traits related to the root system and plant architecture, flowering time, silique, and oil quality. Briefly, the present study provides the recent information and recommendations on the SNP genotyping methods and its applications, which can be useful for marker-assisted breeding in B. napus and other crops.

Quantitative Trait Loci for Genotype and Genotype by Environment Interaction Effects for Seed Yield Plasticity to Terminal Water-Deficit Conditions in Canola (Brassica napus L.)

Canola plants suffer severe crop yield and oil content reductions when exposed to water-deficit conditions, especially during the reproductive stages of plant development. There is a pressing need to develop canola cultivars that can perform better under increased water-deficit conditions with changing weather patterns. In this study, we analysed genetic determinants for the main effects of quantitative trait loci (QTL), (Q), and the interaction effects of QTL and Environment (QE) underlying seed yield and related traits utilising 223 doubled haploid (DH) lines of canola in well-watered and water-deficit conditions under a rainout shelter. Moderate water-deficit at the pre-flowering stage reduced the seed yield to 40.8%. Multi-environmental QTL analysis revealed 23 genomic regions associated with days to flower (DTF), plant height (PH) and seed yield (SY) under well-watered and water-deficit conditions. Three seed yield QTL for main effects were identified on chromosomes A09, C03, and C09, while two were related to QE interactions on A02 and C09. Two QTL regions were co-localised to similar genomic regions for flowering time and seed yield (A09) and the second for plant height and chlorophyll content. The A09 QTL was co-located with a previously mapped QTL for carbon isotope discrimination (Δ¹³C) that showed a positive relationship with seed yield in the same population. Opposite allelic effects for plasticity in seed yield were identified due to QE interactions in response to water stress on chromosomes A02 and C09. Our results showed that QTL's allelic effects for DTF, PH, and SY and their correlation with Δ¹³C are stable across environments (field conditions, previous study) and contrasting water regimes (this study). The QTL and DH lines that showed high yield under well-watered and water-deficit conditions could be used to manipulate water-use efficiency for breeding improved canola cultivars.

Genome-wide association and genomic prediction for resistance to southern corn rust in DH and testcross populations

Southern corn rust (SCR), caused by Puccinia polysora Underw, is a destructive disease that can severely reduce grain yield in maize (Zea mays L.). Owing to P. polysora being multi-racial, it is very important to explore more resistance genes and develop more efficient selection approaches in maize breeding programs. Here, four Doubled Haploid (DH) populations with 384 accessions originated from selected parents and their 903 testcross hybrids were used to perform genome-wide association (GWAS). Three GWAS processes included the additive model in the DH panel, additive and dominant models in the hybrid panel. As a result, five loci were detected on chromosomes 1, 7, 8, 8, and 10, with P-values ranging from 4.83×10^-7 to 2.46×10^-41. In all association analyses, a highly significant locus on chromosome 10 was detected, which was tight chained with the known SCR resistance gene RPPC and RPPK. Genomic prediction (GP), has been proven to be effective in plant breeding. In our study, several models were performed to explore predictive ability in hybrid populations for SCR resistance, including extended GBLUP with different genetic matrices, maker based prediction models, and mixed models with QTL as fixed factors. For GBLUP models, the prediction accuracies ranged from 0.56-0.60. Compared with traditional prediction only with additive effect, prediction ability was significantly improved by adding additive-by-additive effect (P-value< 0.05). For maker based models, the accuracy of BayesA and BayesB was 0.65, 8% higher than other models (i.e., RRBLUP, BRR, BL, BayesC). Finally, by adding QTL into the mixed linear prediction model, the accuracy can be further improved to 0.67, especially for the G_A model, the prediction performance can be increased by 11.67%. The prediction accuracy of the BayesB model can be further improved significantly by adding QTL information (P-value< 0.05). This study will provide important valuable information for understanding the genetic architecture and the application of GP for SCR in maize breeding.

Genetic dissection of branch architecture in oilseed rape (Brassica napus L.) germplasm

Branch architecture is an important factor influencing rapeseed planting density, mechanized harvest, and yield. However, its related genes and regulatory mechanisms remain largely unknown. In this study, branch angle (BA) and branch dispersion degree (BD) were used to evaluate branch architecture. Branch angle exhibited a dynamic change from an increase in the early stage to a gradual decrease until reaching a stable state. Cytological analysis showed that BA variation was mainly due to xylem size differences in the vascular bundle of the branch junction. The phenotypic analysis of 327 natural accessions revealed that BA in six environments ranged from 24.3° to 67.9°, and that BD in three environments varied from 4.20 cm to 21.4 cm, respectively. A total of 115 significant loci were detected through association mapping in three models (MLM, mrMLM, and FarmCPU), which explained 0.53%-19.4% of the phenotypic variations. Of them, 10 loci were repeatedly detected in different environments and models, one of which qBAD.A03-2 was verified as a stable QTL using a secondary segregation population. Totally, 1066 differentially expressed genes (DEGs) were identified between branch adaxial- and abaxial- sides from four extremely large or small BA/BD accessions through RNA sequencing. These DEGs were significantly enriched in the pathways related to auxin biosynthesis and transport as well as cell extension such as indole alkaloid biosynthesis, other glycan degradation, and fatty acid elongation. Four known candidate genes BnaA02g16500D (PIN1), BnaA03g10430D (PIN2), BnaC03g06250D (LAZY1), and BnaC06g20640D (ARF17) were identified by both GWAS and RNA-seq, all of which were involved in regulating the asymmetric distribution of auxins. Our identified association loci and candidate genes provide a theoretical basis for further study of gene cloning and genetic improvement of branch architecture.

Genome-wide association study identifies candidate genes and favorable haplotypes for seed yield in Brassica napus

Oilseed rape (Brassica napus L.) is one of the most essential oil crops. Genetic improvement of seed yield (SY) is a major aim of B. napus breeding. Several studies have been reported on the genetic mechanisms of SY of B. napus. Here, a genome-wide association study (GWAS) of SY was conducted using a panel of 403 natural accessions of B. napus, with more than five million high-quality single-nucleotide polymorphisms (SNPs). A total of 1773 significant SNPs were detected associated with SY, and 783 significant SNPs were co-located with previously reported QTLs. The lead SNPs chrA01__8920351 and chrA02__4555979 were jointly detected in Trial 2_2 and Trial 2_mean value, and in Trial 1_2 and Trial 1_mean value, respectively. Subsequently, two candidate genes of BnaA01g17200D and BnaA02g08680D were identified through combining transcriptome, candidate gene association analysis, and haplotype analysis. BnaA09g10430D detected through lead SNP chrA09__5160639 was associated with SY of B. napus. Our results provide valuable information for studying the genetic control of seed yield in B. napus and valuable genes, haplotypes, and cultivars resources for the breeding of high seed yield B. napus cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01332-6.

Genome-wide association study and transcriptome analysis reveal key genes controlling fruit branch angle in cotton

Fruit branch angle (FBA), a pivotal component of cotton plant architecture, is vital for field and mechanical harvesting. However, the molecular mechanism of FBA formation is poorly understood in cotton. To uncover the genetic basis for FBA formation in cotton, we performed a genome-wide association study (GWAS) of 163 cotton accessions with re-sequencing data. A total of 55 SNPs and 18 candidate genes were significantly associated with FBA trait. By combining GWAS and transcriptome analysis, four genes underlying FBA were identified. An FBA-associated candidate gene Ghi_A09G08736, which is homologous to SAUR46 in Arabidopsis thaliana, was detected in our study. In addition, transcriptomic evidence was provided to show that gravity and light were implicated in the FBA formation. This study provides new insights into the genetic architecture of FBA that informs architecture breeding in cotton.

Further insight into decreases in seed glucosinolate content based on QTL mapping and RNA-seq in Brassica napus L

The QTL hotspots determining seed glucosinolate content instead of only four HAG1 loci and elucidation of a potential regulatory model for rapeseed SGC variation. Glucosinolates (GSLs) are amino acid-derived, sulfur-rich secondary metabolites that function as biopesticides and flavor compounds, but the high seed glucosinolate content (SGC) reduces seed quality for rapeseed meal. To dissect the genetic mechanism and further reduce SGC in rapeseed, QTL mapping was performed using an updated high-density genetic map based on a doubled haploid (DH) population derived from two parents that showed significant differences in SGC. In 15 environments, a total of 162 significant QTLs were identified for SGC and then integrated into 59 consensus QTLs, of which 32 were novel QTLs. Four QTL hotspot regions (QTL-HRs) for SGC variation were discovered on chromosomes A09, C02, C07 and C09, including seven major QTLs that have previously been reported and four novel major QTLs in addition to HAG1 loci. SGC was largely determined by superimposition of advantage allele in the four QTL-HRs. Important candidate genes directly related to GSL pathways were identified underlying the four QTL-HRs, including BnaC09.MYB28, BnaA09.APK1, BnaC09.SUR1 and BnaC02.GTR2a. Related differentially expressed candidates identified in the minor but environment stable QTLs indicated that sulfur assimilation plays an important rather than dominant role in SGC variation. A potential regulatory model for rapeseed SGC variation constructed by combining candidate GSL gene identification and differentially expressed gene analysis based on RNA-seq contributed to a better understanding of the GSL accumulation mechanism. This study provides insights to further understand the genetic regulatory mechanism of GSLs, as well as the potential loci and a new route to further diminish the SGC in rapeseed.

Identification of a Major QTL and Validation of Related Genes for Tiller Angle in Rice Based on QTL Analysis

An ideal plant architecture is an important condition to achieve high crop yields. The tiller angle is an important and complex polygenic trait of rice (Oryza sativa L.) plant architecture. Therefore, the discovery and identification of tiller angle-related genes can aid in the improvement of crop architecture and yield. In the present study, 222 SSR markers were used to establish a high-density genetic map of rice doubled haploid population, and a total of 8 quantitative trait loci (QTLs) were detected based on the phenotypic data of the tiller angle and tiller crown width over 2 years. Among them, four QTLs (qTA9, qCW9, qTA9-1, and qCW9-1) were overlapped at marker interval RM6235-RM24288 on chromosome 9 with a large effect value regarded as a stable major QTL. The selected promising related genes were further identified by relative gene expression analysis, which gives us a basis for the future cloning of these genes. Finally, OsSAURq9, which belongs to the SMALL AUXIN UP RNA (SAUR), an auxin-responsive protein family, was selected as a target gene. Overall, this work will help broaden our knowledge of the genetic control of tiller angle and tiller crown width, and this study provides both a good theoretical basis and a new genetic resource for the breeding of ideal-type rice.

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.

Identification of quantitative trait loci (QTLs) and candidate genes of seed Iron and zinc content in soybean [Glycine max (L.) Merr.]

BACKGROUND

Deciphering the hereditary mechanism of seed iron (Fe) and zinc (Zn) content in soybean is important and sustainable to address the "hidden hunger" that presently affects approximately 2 billion people worldwide. Therefore, in order to detect genomic regions related to soybean seed Fe and Zn content, a recombinant inbred line (RIL) population with 248 lines was assessed in four environments to detect Quantitative Trait Loci (QTLs) related to soybean seed Fe and Zn content.

RESULT

Wide variation was found in seed Fe and Zn content in four environments, and genotype, environment, and genotype × environment interactions had significant influences on both the seed Fe and Zn content. A positive correlation was observed between seed Fe content and seed Zn content, and broad-sense heritability (H2) of seed Fe and Zn content were 0.73 and 0.75, respectively. In this study, five QTLs for seed Fe content were detected with 4.57 - 32.71% of phenotypic variation explained (PVE) and logarithm of odds (LOD) scores ranging from 3.60 to 33.79. Five QTLs controlling the seed Zn content were detected, and they individually explained 3.35 to 26.48% of the phenotypic variation, with LOD scores ranging from 3.64 to 20.4. Meanwhile, 409,541 high-quality single-nucleotide variants (SNVs) and 85,102 InDels (except intergenic regions) between two bi-parental lines were identified by whole genome resequencing. A total of 12 candidate genes were reported in one major QTL for seed Fe content and two major QTLs for seed Zn content, with the help of RNA-Seq analysis, gene ontology (GO) enrichment, gene annotation, and bi-parental whole genome sequencing (WGS) data.

CONCLUSIONS

Limited studies were performed about microelement of soybean, so these results may play an important role in the biofortification of Fe and Zn and accelerate the development of marker-assisted selection (MAS) for breeding soybeans fortified with iron and zinc.

QTL Identification for Stem Fiber, Strength and Rot Resistance in a DH Population from an Alien Introgression of Brassica napus

Stem fiber, stem strength and stem-rot resistance are important agronomic traits in Brassica napus. To understand the molecular mechanism that controls the stem-related traits, we investigated the stem lignin (ADL), cellulose (Cel), hemicellulose (Hem) content, S/G monolignol ratio (SG), stem breaking force (BF), breaking strength (F) and Sclerotinia sclerotiorum resistance (SSR). Each trait was significantly positively or negatively correlated with more than three of the other six traits. QTL mapping for ADL, Cel, Hem, SG, BF, F and SSR were performed using a doubled haploid population derived from an intertribal B. napus introgression line 'Y689' crossed with B. napus cv. 'Westar'. A total of 67 additive QTL were identified and integrated into 55 consensus QTL by meta-analysis. Among the 55 consensus QTL, 23 (41.8%) QTL were co-located and were integrated into 11 unique QTL. The QTL by environment (Q × E) interactions were analyzed and 22 combined QTL were identified. In addition, candidate genes within the QTL intervals were proposed based on the known function of Arabidopsis orthologs. These results provided valuable information for improving lodging resistance, S. sclerotiorum resistance and mechanized harvesting of B. napus.

The Bnapus50K array: a quick and versatile genotyping tool for Brassica napus genomic breeding and research

Rapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for B. napus genomic breeding and research.

QTL Mapping and Diurnal Transcriptome Analysis Identify Candidate Genes Regulating Brassica napus Flowering Time

Timely flowering is important for seed formation and maximization of rapeseed (Brassica napus) yield. Here, we performed flowering-time quantitative trait loci (QTL) mapping using a double haploid (DH) population grown in three environments to study the genetic architecture. Brassica 60 K Illumina Infinium™ single nucleotide polymorphism (SNP) array and simple sequence repeat (SSR) markers were used for genotyping of the DH population, and a high-density genetic linkage map was constructed. QTL analysis of flowering time from the three environments revealed five consensus QTLs, including two major QTLs. A major QTL located on chromosome A03 was detected specifically in the semi-winter rapeseed growing region, and the one on chromosome C08 was detected in all environments. Ribonucleic acid sequencing (RNA-seq) was performed on the parents’ leaves at seven time-points in a day to determine differentially expressed genes (DEGs). The biological processes and pathways with significant enrichment of DEGs were obtained. The DEGs in the QTL intervals were analyzed, and four flowering time-related candidate genes were found. These results lay a foundation for the genetic regulation of rapeseed flowering time and create a rapeseed gene expression library for seven time-points in a day.

Identification and fine mapping of qSB.A09, a major QTL that controls shoot branching in Brassica rapa ssp. chinensis Makino

QTL mapping plus bulked segregant analysis revealed a major QTL for shoot branching in non-heading Chinese cabbage. The candidate gene was then identified using sequence alignment and expression analysis. Shoot branching is a complex quantitative trait that contributes to plant architecture and ultimately yield. Although many studies have examined branching in grain crops, the genetic control of shoot branching in vegetable crops such as Brassica rapa L. ssp. chinensis remains poorly understood. In this study, we used bulked segregant analysis (BSA) of an F₂ population to detect a major quantitative trait locus (QTL) for shoot branching, designated shoot branching 9 (qSB.A09) on the long arm of chromosome A09 in Brassica rapa L. ssp. chinensis. In addition, traditional QTL mapping of the F₂ population revealed six QTLs in different regions. Of these, the mapping region on chromosome A09 was consistent with the results of BSA-seq analysis, as well as being stable over the 2-year study period, explaining 19.37% and 22.18% of the phenotypic variation across multiple genetic backgrounds. Using extreme recombinants, qSB.A09 was further delimited to a 127-kb genomic region harboring 28 annotated genes. We subsequently identified the GRAS transcription factor gene Bra007056 as a potential candidate gene; Bra007056 is an ortholog of MONOCULM 1 (MOC1), the key gene that controls tillering in rice. Quantitative RT-PCR further revealed that expression of Bra007056 was positively correlated with the shoot branching phenotype. Furthermore, an insertion/deletion marker specific to Bra007056 co-segregated with the shoot branching trait in the F₂ populations. Overall, these results provide the basis for elucidating the molecular mechanism of shoot branching in Brassica rapa ssp. chinensis Makino.

Construction of a High-Density Genetic Map and Identification of Loci Related to Hollow Stem Trait in Broccoli (Brassic oleracea L. italica)

A high-quality genetic map is important for mapping of compound traits. In this study, a genetic map was constructed based on the reference genome TO1000 after specific locus amplified fragment (SLAF) sequencing in a double-haploid segregation population of broccoli, and loci controlling hollow stem trait were identified in the genetic map. The genetic map contains 4,787 SLAF markers, with a mean marker distance of 0.22 cM and the mean sequencing depths of 91.14-fold in the maternal line, 88.97-fold in the paternal line and 17.11-fold in each DH progeny. A locus controlling the hollow stem trait, QHS.C09-2, which could explain 14.1% of the phenotypic variation, was steadily detected on the linkage group nine in the indicated data of 3 years' trials and BLUE analysis. The genetic map could lay an important foundation for mapping of compound traits, and mapping of hollow stem trait would be basis to clone the genes related to hollow stems in broccoli.

An APETALA1 ortholog affects plant architecture and seed yield component in oilseed rape (Brassica napus L.)

BACKGROUND

Increasing the productivity of rapeseed as one of the widely cultivated oil crops in the world is of upmost importance. As flowering time and plant architecture play a key role in the regulation of rapeseed yield, understanding the genetic mechanism underlying these traits can boost the rapeseed breeding. Meristem identity genes are known to have pleiotropic effects on plant architecture and seed yield in various crops. To understand the function of one of the meristem identity genes, APETALA1 (AP1) in rapeseed, we performed phenotypic analysis of TILLING mutants under greenhouse conditions. Three stop codon mutant families carrying a mutation in Bna.AP1.A02 paralog were analyzed for different plant architecture and seed yield-related traits.

RESULTS

It was evident that stop codon mutation in the K domain of Bna.AP1.A02 paralog caused significant changes in flower morphology as well as plant architecture related traits like plant height, branch height, and branch number. Furthermore, yield-related traits like seed yield per plant and number of seeds per plants were also significantly altered in the same mutant family. Apart from phenotypic changes, stop codon mutation in K domain of Bna.AP1.A02 paralog also altered the expression of putative downstream target genes like Bna.TFL1 and Bna.FUL in shoot apical meristem (SAM) of rapeseed. Mutant plants carrying stop codon mutations in the COOH domain of Bna.AP1.A02 paralog did not have a significant effect on plant architecture, yield-related traits or the expression of the downstream targets.

CONCLUSIONS

We found that Bna.AP1.A02 paralog has pleiotropic effect on plant architecture and yield-related traits in rapeseed. The allele we found in the current study with a beneficial effect on seed yield can be incorporated into rapeseed breeding pool to develop new varieties.

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.