Genome-Wide Association Mapping Reveals the Genetic Control Underlying Branch Angle in Rapeseed (Brassica napus L.)

Homoeologous exchanges contribute to branch angle variations in rapeseed: Insights from transcriptome, QTL-seq and gene functional analysis

Branch angle (BA) is a critical morphological trait that significantly influences planting density, light interception and ultimately yield in plants. Despite its importance, the regulatory mechanism governing BA in rapeseed remains poorly understood. In this study, we generated 109 transcriptome data sets for 37 rapeseed accessions with divergent BA phenotypes. Relative to adaxial branch segments, abaxial segments accumulated higher levels of auxin and exhibited lower expression of six TCP1 homologues and one GA20ox3. A co-expression network analysis identified two modules highly correlated with BA. The modules contained homologues to known BA control genes, such as FUL, YUCCA6, TCP1 and SGR3. Notably, a homoeologous exchange (HE), occurring at the telomeres of A09, was prevalent in large BA accessions, while an A02-C02 HE was common in small BA accessions. In their corresponding regions, these HEs explained the formation of hub gene hotspots in the two modules. QTL-seq analysis confirmed that the presence of a large A07-C06 HE (~8.1 Mb) was also associated with a small BA phenotype, and BnaA07.WRKY40.b within it was predicted as candidate gene. Overexpressing BnaA07.WRKY40.b in rapeseed increased BA by up to 20°, while RNAi- and CRISPR-mediated mutants (BnaA07.WRKY40.b and BnaC06.WRKY40.b) exhibited decreased BA by up to 11.4°. BnaA07.WRKY40.b was exclusively localized to the nucleus and exhibited strong expression correlations with many genes related to gravitropism and plant architecture. Taken together, our study highlights the influence of HEs on rapeseed plant architecture and confirms the role of WRKY40 homologues as novel regulators of BA.

Genome-wide association study and genomic selection of flax powdery mildew in Xinjiang Province

Flax powdery mildew (PM), caused by Oidium lini, is a globally distributed fungal disease of flax, and seriously impairs its yield and quality. To data, only three resistance genes and a few putative quantitative trait loci (QTL) have been reported for flax PM resistance. To dissect the resistance mechanism against PM and identify resistant genetic regions, based on four years of phenotypic datasets (2017, 2019 to 2021), a genome-wide association study (GWAS) was performed on 200 flax core accessions using 674,074 SNPs and 7 models. A total of 434 unique quantitative trait nucleotides (QTNs) associated with 331 QTL were detected. Sixty-four loci shared in at least two datasets were found to be significant in haplotype analyses, and 20 of these sites were shared by multiple models. Simultaneously, a large-effect locus (qDI 11.2) was detected repeatedly, which was present in the mapping study of flax pasmo resistance loci. Oil flax had more QTL with positive-effect or favorable alleles (PQTL) and showed higher PM resistance than fiber flax, indicating that effects of these QTL were mainly additive. Furthermore, an excellent resistant variety C120 was identified and can be used to promote planting. Based on 331 QTLs identified through GWAS and the statistical model GBLUP, a genomic selection (GS) model related to flax PM resistance was constructed, and the prediction accuracy rate was 0.96. Our results provide valuable insights into the genetic basis of resistance and contribute to the advancement of breeding programs.

Root system architecture change in response to waterlogging stress in a 448 global collection of rapeseeds (Brassica napus L.)

MAIN CONCLUSIONS

A novel image-based screening method for precisely identifying genotypic variations in rapeseed RSA under waterlogging stress was developed. Five key root traits were confirmed as good indicators of waterlogging and might be employed in breeding, particularly when using the MFVW approach. Waterlogging is a vital environmental factor that has detrimental effects on the growth and development of rapeseed (Brassica napus L.). Plant roots suffer from hypoxia under waterlogging, which ultimately confers yield penalty. Therefore, it is crucially important to understand the genetic variation of root system architecture (RSA) in response to waterlogging stress to guide the selection of new tolerant cultivars with favorable roots. This research was conducted to investigate RSA traits using image-based screening techniques to better understand how RSA changes over time during waterlogging at the seedling stage. First, we performed a t-test by comparing the relative root trait value between four tolerant and four sensitive accessions. The most important root characteristics associated with waterlogging tolerance at 12 h are total root length (TRL), total root surface area (TRSA), total root volume (TRV), total number of tips (TNT), and total number of forks (TNF). The root structures of 448 rapeseed accessions with or without waterlogging showed notable genetic diversity, and all traits were generally restrained under waterlogging conditions, except for the total root average diameter. Additionally, according to the evaluation and integration analysis of 448 accessions, we identified that five traits, TRL, TRSA, TRV, TNT, and TNF, were the most reliable traits for screening waterlogging-tolerant accessions. Using analysis of the membership function value (MFVW) and D-value of the five selected traits, 25 extremely waterlogging-tolerant materials were screened out. Waterlogging significantly reduced RSA, inhibiting root growth compared to the control. Additionally, waterlogging increased lipid peroxidation, accompanied by a decrease in the activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT). This study effectively improves our understanding of the response of RSA to waterlogging. The image-based screening method developed in this study provides a new scientific guidance for quickly examining the basic RSA changes and precisely predicting waterlogging-tolerant rapeseed germplasms, thus expanding the genetic diversity of waterlogging-tolerant rapeseed germplasm available for breeding.

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.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

IGT/LAZY genes are differentially influenced by light and required for light-induced change to organ angle

BACKGROUND

Plants adjust their growth orientations primarily in response to light and gravity signals. Considering that the gravity vector is fixed and the angle of light incidence is constantly changing, plants must somehow integrate these signals to establish organ orientation, commonly referred to as gravitropic set-point angle (GSA). The IGT gene family contains known regulators of GSA, including the gene clades LAZY, DEEPER ROOTING (DRO), and TILLER ANGLE CONTROL (TAC).

RESULTS

Here, we investigated the influence of light on different aspects of GSA phenotypes in LAZY and DRO mutants, as well as the influence of known light signaling pathways on IGT gene expression. Phenotypic analysis revealed that LAZY and DRO genes are collectively required for changes in the angle of shoot branch tip and root growth in response to light. Single lazy1 mutant branch tips turn upward in the absence of light and in low light, similar to wild-type, and mimic triple and quadruple IGT mutants in constant light and high-light conditions, while triple and quadruple IGT/LAZY mutants show little to no response to changing light regimes. Further, the expression of IGT/LAZY genes is differentially influenced by daylength, circadian clock, and light signaling.

CONCLUSIONS

Collectively, the data show that differential expression of LAZY and DRO genes are required to enable plants to alter organ angles in response to light-mediated signals.

Peach LAZY1 and DRO1 protein-protein interactions and co-expression with PRAF/RLD family support conserved gravity-related protein interactions across plants

IGT/LAZY proteins play a central role in determining gravitropic set point angle and orientation of lateral organs across plant species. Recent work in model systems has demonstrated that interactions between IGT/LAZY proteins and BREVIS RADIX (BRX)-domain containing proteins, such as PH, RCC1, AND FYVE/RCC1-LIKE DOMAIN (PRAF/RLD), and BREVIS RADIX LIKE (BRXL) family members, are mechanistically important for setting gravitropic set point angle. Here, we identified peach PRAF/RLD proteins as interactors of the peach IGT/LAZY proteins PpeLAZY1 and DEEPER ROOTING 1 (PpeDRO1) from a yeast-two-hybrid screen. We also show that the BRX domains of these interacting proteins have high sequence similarity with PRAF/RLD and BRX family proteins from rice and Arabidopsis. Further, PpeLAZY1 and the peach PRAF/RLD interactors are all expressed at relatively high levels in leaf, meristem, and shoot tip tissues. Together, this evidence supports the importance and conservation of IGT/LAZY-BRX-domain interactions, which underlie setting gravitropic set point angle across angiosperms.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3-45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

Genome-wide association study reveals genomic loci influencing agronomic traits in Ethiopian sorghum (Sorghum bicolor (L.) Moench) landraces

Uncovering the genetic basis of agronomic traits in sorghum landraces that have adapted to various agro-climatic conditions would contribute to sorghum improvement efforts around the world. To identify quantitative trait nucleotides (QTNs) associated with nine agronomic traits in a panel of 304 sorghum accessions collected from diverse environments across Ethiopia (considered to be the center of origin and diversity), multi-locus genome-wide association studies (ML-GWAS) were performed using 79,754 high quality single nucleotide polymorphism (SNP) markers. Association analyses using six ML-GWAS models identified a set of 338 significantly (LOD ≥ 3)-associated QTNs for nine agronomic traits of sorghum accessions evaluated in two environments (E1 and E2) and their combined dataset (Em). Of these, 121 reliable QTNs, including 13 for flowering time (DF), 13 for plant height (PH), 9 for tiller number (TN), 15 for panicle weight (PWT), 30 for grain yield per panicle (GYP), 12 for structural panicle mass (SPM), 13 for hundred seed weight (HSW), 6 for grain number per panicle (GNP), and 10 for panicle exertion (PE) were consistently detected by at least three ML-GWAS methods and/or in two different environments. Notably, Ethylene responsive transcription factor gene AP2/ERF, known for regulation of plant growth, and the sorghum Terminal flower1/TF1 gene, which functions in the control of floral architecture, were identified as strong candidate genes associated with PH and HSW, respectively. This study provides an entry point for further validation studies to elucidate complex mechanisms controlling important agronomic traits in sorghum.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01381-5.

Branch angle and leaflet shape are associated with canopy coverage in soybean

Early canopy coverage is a desirable trait that is a major determinant of yield in soybean (Glycine max). Variation in traits comprising shoot architecture can influence canopy coverage, canopy light interception, canopy-level photosynthesis, and source-sink partitioning efficiency. However, little is known about the extent of phenotypic diversity of shoot architecture traits and their genetic control in soybean. Thus, we sought to understand the contribution of shoot architecture traits to canopy coverage and to determine the genetic control of these traits. We examined the natural variation for shoot architecture traits in a set of 399 diverse maturity group I soybean (SoyMGI) accessions to identify relationships between traits, and to identify loci that are associated with canopy coverage and shoot architecture traits. Canopy coverage was correlated with branch angle, number of branches, plant height, and leaf shape. Using previously collected 50K single nucleotide polymorphism data, we identified quantitative trait locus (QTL) associated with branch angle, number of branches, branch density, leaflet shape, days to flowering, maturity, plant height, number of nodes, and stem termination. In many cases, QTL intervals overlapped with previously described genes or QTL. We also found QTL associated with branch angle and leaflet shape located on chromosomes 19 and 4, respectively, and these QTL overlapped with QTL associated with canopy coverage, suggesting the importance of branch angle and leaflet shape in determining canopy coverage. Our results highlight the role individual architecture traits play in canopy coverage and contribute information on their genetic control that could help facilitate future efforts in their genetic manipulation.

The Role of Italy in the Use of Advanced Plant Genomic Techniques on Fruit Trees: State of the Art and Future Perspectives

Climate change is deeply impacting the food chain production, lowering quality and yield. In this context, the international scientific community has dedicated many efforts to enhancing resilience and sustainability in agriculture. Italy is among the main European producers of several fruit trees; therefore, national research centers and universities undertook several initiatives to maintain the specificity of the 'Made in Italy' label. Despite their importance, fruit crops are suffering from difficulties associated with the conventional breeding approaches, especially in terms of financial commitment, land resources availability, and long generation times. The 'new genomic techniques' (NGTs), renamed in Italy as 'technologies for assisted evolution' (TEAs), reduce the time required to obtain genetically improved cultivars while precisely targeting specific DNA sequences. This review aims to illustrate the role of the Italian scientific community in the use of NGTs, with a specific focus on Citrus, grapevine, apple, pear, chestnut, strawberry, peach, and kiwifruit. For each crop, the key genes and traits on which the scientific community is working, as well as the technological improvements and advancements on the regeneration of local varieties, are presented. Lastly, a focus is placed on the legal aspects in the European and in Italian contexts.

BSA‑seq and genetic mapping identified candidate genes for branching habit in peanut

The candidate gene AhLBA1 controlling lateral branch angel of peanut was fine-mapped to a 136.65-kb physical region on chromosome 15 using the BSA-seq and QTL mapping. Lateral branch angel (LBA) is an important plant architecture trait of peanut, which plays key role in lodging, peg soil penetration and pod yield. However, there are few reports of fine mapping and quantitative trait loci (QTLs)/cloned genes for LBA in peanut. In this project, a mapping population was constructed using a spreading variety Tifrunner and the erect variety Fuhuasheng. Through bulked segregant analysis sequencing (BSA-seq), a major gene related to LBA, named as AhLBA1, was preliminarily mapped at the region of Chr.15: 150-160 Mb. Then, using traditional QTL approach, AhLBA1 was narrowed to a 1.12 cM region, corresponding to a 136.65-kb physical interval of the reference genome. Of the nine genes housed in this region, three of them were involved in hormone metabolism and regulation, including one "F-box protein" and two "2-oxoglutarate (2OG) and Fe(II)-dependent oxygenase (2OG oxygenase)" encoding genes. In addition, we found that the level of some classes of cytokinin (CK), auxin and ethylene showed significant differences between spreading and erect peanuts at the junction of main stem and lateral branch. These findings will aid further elucidation of the genetic mechanism of LBA in peanut and facilitating marker-assisted selection (MAS) in the future breeding program.

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.

The Citrus sinensis TILLER ANGLE CONTROL 1 (CsTAC1) gene regulates tree architecture in sweet oranges by modulating the endogenous hormone content

Citrus is a major fruit crop cultivated on a global scale. Citrus trees are long lived perennials with a large canopy. Understanding the genetic control of tree architecture could provide tools for breeding and selection of citrus cultivars suitable for high density planting with improved light exposure. Tree architecture is modulated by the TILLER ANGLE CONTROL 1 (TAC1) gene which plays an important role in the regulation of the shoot angle. Herein, we used CRISPR/Cas9 technology to knockout the CsTAC1 gene for the biochemical and molecular analysis of its function. Nine transgenic lines were obtained, and five edited plants were confirmed based on T7EI mismatch detection assay and Sanger sequencing. The transgenic citrus lines exhibited pleiotropic phenotypes, including differences in branch angle and stem growth. Additionally, silencing CsTAC1 led to enhanced CsLAZY1 transcript levels in the tested lines. Analysis of the phytohormonal profile revealed that TAC1-edited plants exhibited lower auxin contents and increased cytokinin levels in the leaves compared to the wild-type plants. The GA₇ gibberellin level was enhanced in most of the edited lines. Collectively, TAC1 affects branch angle in association with hormone signals in citrus.

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.

Genetic dissection of the natural variation of ovule number per ovary in oilseed rape germplasm (Brassica napus L.)

Oilseed rape is one of the world's largest oil and industrial crops, providing humans with various products, such as vegetable oil and biofuel. Ovules are the direct precursors of seeds, and ovule number per ovary (ONPO) largely determines seed number per fruit that affects both yield and fitness of seed crops. The ONPO shows wide variation in oilseed rape, whereas the underlying genes and mechanisms are poorly known. The present study performed the genetic, physiological and transcriptomic analyses of ovule number per ovary using an association panel and the extreme lines. The ONPO of 327 accessions planted in four environments showed a large variation from 19.2 to 43.8, indicating a great potential for the further genetic improvement of ovule number. The genome-wide association study (GWAS) identified a total of 43 significant SNP markers. Further, these SNPs were integrated into 18 association loci, which were distributed on chromosomes A01, A03, A06, A07, A09, C01, C03, C06, C07, and C09, explaining 4.3-11.5% of the phenotypic variance. The ONPO decreased as their appearance order on the inflorescence and was associated with the level of several types of endogenous phytohormones but not related to leaf area and photosynthetic rate. Comparative transcriptomic analysis identified a total of 4,449 DEGs enriched in 30 classes, including DNA, RNA, protein, signaling, transport, development, cell wall, lipid metabolism, and secondary metabolism. Nearly half of DEGs were involved in the known pathways in regulating ovule number, of which 12 were homologous to know ovule number regulating genes, indicating a strong link between the identified DEGs and ovule number. A total of 73 DEGs were located within the genomic regions of association loci, of which six were identified as candidates based on functional annotation. These results provide useful information for the further genetic improvement of ovule and seed number in oilseed rape.

GWAS and bulked segregant analysis reveal the Loci controlling growth habit-related traits in cultivated Peanut (Arachis hypogaea L.)

Background

Peanut (Arachis hypogaea L.) is a grain legume crop that originated from South America and is now grown around the world. Peanut growth habit affects the variety’s adaptability, planting patterns, mechanized harvesting, disease resistance, and yield. The objective of this study was to map the quantitative trait locus (QTL) associated with peanut growth habit-related traits by combining the genome-wide association analysis (GWAS) and bulked segregant analysis sequencing (BSA-seq) methods.

Results

GWAS was performed with 17,223 single nucleotide polymorphisms (SNPs) in 103 accessions of the U.S. mini core collection genotyped using an Affymetrix version 2.0 SNP array. With a total of 12,342 high-quality polymorphic SNPs, the 90 suggestive and significant SNPs associated with lateral branch angle (LBA), main stem height (MSH), lateral branch height (LBL), extent radius (ER), and the index of plant type (IOPT) were identified. These SNPs were distributed among 15 chromosomes. A total of 597 associated candidate genes may have important roles in biological processes, hormone signaling, growth, and development. BSA-seq coupled with specific length amplified fragment sequencing (SLAF-seq) method was used to find the association with LBA, an important trait of the peanut growth habit. A 4.08 Mb genomic region on B05 was associated with LBA. Based on the linkage disequilibrium (LD) decay distance, we narrowed down and confirmed the region within the 160 kb region (144,193,467–144,513,467) on B05. Four candidate genes in this region were involved in plant growth. The expression levels of Araip.E64SW detected by qRT-PCR showed significant difference between ‘Jihua 5’ and ‘M130’.

Conclusions

In this study, the SNP (AX-147,251,085 and AX-144,353,467) associated with LBA by GWAS was overlapped with the results in BSA-seq through combined analysis of GWAS and BSA-seq. Based on LD decay distance, the genome range related to LBA on B05 was shortened to 144,193,467–144,513,467. Three candidate genes related to F-box family proteins (Araip.E64SW, Araip.YG1LK, and Araip.JJ6RA) and one candidate gene related to PPP family proteins (Araip.YU281) may be involved in plant growth and development in this genome region. The expression analysis revealed that Araip.E64SW was involved in peanut growth habits. These candidate genes will provide molecular targets in marker-assisted selection for peanut growth habits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08640-3.

Construction of a high-density genetic map using specific-length amplified fragment markers and identification of QTLs for branching angle in poplar

Branching angle is a critical factor that determines the morphological establishment and is a typical quantitative trait controlled by multiple genes. In this study, we used SLAF-seq to construct a high-density genetic map, to investigate the genetic architecture of branching angle in poplar (Populus leucopyramidalis). A total of 240,672 SLAF tags were obtained, including 103,691 polymorphic SLAF tags. After filtering, 53,407 polymorphic markers were sorted into eight segregation types, and 11,162 of them were used to construct the genetic map. 8447 were on the female parent map, 8532 were on the male parent map, and 11,162 were on the integrated map. The marker coverage was 4820.84 and 5044.80 cM for the female and male maps, and 3142.61 cM for the integrated map. The average intervals between two adjacent mapped markers were 0.55, 0.59, and 0.28 cM for the three maps, respectively. Two quantitative trait loci (QTLs) were detected. Seven markers that exceeded the threshold in these two regions were considered as being associated with branching angle and the phenotypic variance explained by each of these marker was 10.64-11.66%. After functional annotation, we identified 15 candidate genes and analyzed the expression of candidate genes in narrow and wide crown progenies by qRT-PCR. These results show that the combination of QTL and SLAF-seq will contribute to future breeding plans in poplar breeding, especially in narrow crown poplar breeding.

Developing Chenopodium ficifolium as a potential B genome diploid model system for genetic characterization and improvement of allotetraploid quinoa (Chenopodium quinoa)

BACKGROUND

Quinoa (Chenopodium quinoa) is a high-value grain known for its excellent nutritional balance. It is an allotetraploid species (AABB, 2n = 4x = 36) formed by the hybridization between AA and BB genome diploid (2n = 2x = 18) species. This study reports genetic studies in Chenopodium ficifolium as a potential B genome diploid model system to simplify the genetic studies of quinoa including gene identification and marker-assisted breeding.

RESULTS

Portsmouth, New Hampshire and Quebec City, Quebec accessions of C. ficifolium were used to develop an F2 population segregating for agronomically relevant traits including flowering time, plant height, the number of branches, branch angle, and internode length. Marker-trait associations were identified for the FLOWERING LOCUS T-LIKE 1 (FTL1) marker gene, where the alternate alleles (A1/A2) were segregating among the F2 generation plants in association with flowering time, plant height, and the number of branches. There was a strong correlation of the flowering time trait with both plant height and the number of branches. Thus, a possible multifaceted functional role for FTL1 may be considered. The parental Portsmouth and Quebec City accessions were homozygous for the alternate FTL1 alleles, which were found to be substantially diverged. SNPs were identified in the FTL1 coding sequence that could have some functional significance in relation to the observed trait variation.

CONCLUSION

These results draw further attention to the possible functional roles of the FTL1 locus in Chenopodium and justify continued exploration of C. ficifolium as a potential diploid model system for the genetic study of quinoa. We expect our findings to aid in quinoa breeding as well as to any studies related to the Chenopodium genus.

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.

Proteomic analysis of a clavata-like phenotype mutant in Brassica napus

Rapeseed is one of important oil crops in China. Better understanding of the regulation network of main agronomic traits of rapeseed could improve the yielding of rapeseed. In this study, we obtained an influrescence mutant that showed a fusion phenotype, similar with the Arabidopsis clavata-like phenotype, so we named the mutant as Bnclavata-like (Bnclv-like). Phenotype analysis illustrated that abnormal development of the inflorescence meristem (IM) led to the fused-inflorescence phenotype. At the stage of protein abundance, major regulators in metabolic processes, ROS metabolism, and cytoskeleton formation were seen to be altered in this mutant. These results not only revealed the relationship between biological processes and inflorescence meristem development, but also suggest bioengineering strategies for the improved breeding and production of Brassica napus.

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.

Genomic Prediction Accuracy of Seven Breeding Selection Traits Improved by QTL Identification in Flax

Molecular markers are one of the major factors affecting genomic prediction accuracy and the cost of genomic selection (GS). Previous studies have indicated that the use of quantitative trait loci (QTL) as markers in GS significantly increases prediction accuracy compared with genome-wide random single nucleotide polymorphism (SNP) markers. To optimize the selection of QTL markers in GS, a set of 260 lines from bi-parental populations with 17,277 genome-wide SNPs were used to evaluate the prediction accuracy for seed yield (YLD), days to maturity (DTM), iodine value (IOD), protein (PRO), oil (OIL), linoleic acid (LIO), and linolenic acid (LIN) contents. These seven traits were phenotyped over four years at two locations. Identification of quantitative trait nucleotides (QTNs) for the seven traits was performed using three types of statistical models for genome-wide association study: two SNP-based single-locus (SS), seven SNP-based multi-locus (SM), and one haplotype-block-based multi-locus (BM) models. The identified QTNs were then grouped into QTL based on haplotype blocks. For all seven traits, 133, 355, and 1208 unique QTL were identified by SS, SM, and BM, respectively. A total of 1420 unique QTL were obtained by SS+SM+BM, ranging from 254 (OIL, LIO) to 361 (YLD) for individual traits, whereas a total of 427 unique QTL were achieved by SS+SM, ranging from 56 (YLD) to 128 (LIO). SS models alone did not identify sufficient QTL for GS. The highest prediction accuracies were obtained using single-trait QTL identified by SS+SM+BM for OIL (0.929 ± 0.016), PRO (0.893 ± 0.023), YLD (0.892 ± 0.030), and DTM (0.730 ± 0.062), and by SS+SM for LIN (0.837 ± 0.053), LIO (0.835 ± 0.049), and IOD (0.835 ± 0.041). In terms of the number of QTL markers and prediction accuracy, SS+SM outperformed other models or combinations thereof. The use of all SNPs or QTL of all seven traits significantly reduced the prediction accuracy of traits. The results further validated that QTL outperformed high-density genome-wide random markers, and demonstrated that the combined use of single and multi-locus models can effectively identify a comprehensive set of QTL that improve prediction accuracy, but further studies on detection and removal of redundant or false-positive QTL to maximize prediction accuracy and minimize the number of QTL markers in GS are warranted.

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.

De novo transcriptome characterisation of two auxin-related genes associated with plant growth habit in Astragalus adsurgens Pall

Astragalus adsurgens Pall., a perennial legume native to China, is commonly used as a forage crop. And it has great value for sustainable development of grasslands in arid and semi arid regions. However, to date, little is known regarding the A. adsurgens genome, and no studies have determined whether it would be possible to improve the germplasm of A. adsurgens through genetic modification. In this study, we used an RNA-seq protocol to generate a de novo transcriptome including 151,516 unigenes of A. adsurgens. We compared the transcriptomes of A. adsurgens having different growth habits (prostrate/erect) and identified 14,133 single nucleotide polymorphism sites (SNP) in 8,139 unigenes. Differential expression gene (DEG) analysis suggested that 10,982 unigenes were up-regulated in the prostrate plant relative to the erect plant, while 10,607 unigenes were down-regulated. Of the 21,589 DEG, Unigene72782_All (LAX4) and CL12494.Contig3_All (TIR1), an auxin transporter gene and an auxin transport inhibitor gene, respectively, were predicted to influence the growth habit of A. adsurgens, which were verified by qRT-PCR in these phenotypes. These results suggest that auxin transport was more active in the prostrate plant than in the erect plant, resulting in asymmetric distribution of auxin that affects the growth habit of A. adsurgens. Overall, this study may provide a basis for future research on key genes in A. adsurgens and may deepen our understanding of the molecular mechanisms regulating plant growth habit.

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.

Comparing Single-SNP, Multi-SNP, and Haplotype-Based Approaches in Association Studies for Major Traits in Barley

The multiple single nucleotide polymorphism (multi-SNP) and haplotype-based approaches that jointly consider multiple markers unveiled a larger number of associations, some of which were shared with the single-SNP approach. A larger overlap of quantitative trait loci (QTLs) between the single-SNP and haplotype-based approaches was obtained than with the multi-SNP approach. Despite a limited overlap between the QTLs detected by these approaches, each uncovered QTLs reported previously, suggesting that each approach is capable of uncovering a different subset of QTLs. We demonstrated the efficiency of an integrated genome-wide association study (GWAS) procedure, combining single-locus and multilocus approaches to improve the capacity and reliability of association analysis to detect key QTLs. The efficiency of barley breeding programs may be improved by the practical use of QTLs identified in this study. Genome-wide association studies (GWAS) have been widely used to identify quantitative trait loci (QTLs) underlying complex agronomic traits. The conventional GWAS model is based on a single-locus model, which may prove inaccurate if a trait is controlled by multiple loci, which is the case for most agronomic traits in barley (Hordeum vulgare L.). Additionally, an individual single nucleotide polymorphism (SNP) will prove incapable of capturing underlying allelic diversity. A multilocus model could potentially represent a better alternative for QTL identification. This study aimed to explore different GWAS approaches (single-SNP, multi-SNP, and haplotype-based) to establish SNP-trait associations and to potentially describe the complex genetic architecture of seven key traits in spring barley. The multi-SNP and haplotype-based approaches unveiled a larger number of significant associations, some of which were shared with the single-SNP approach. Globally, the multi-SNP approach explained more of the phenotypic variance (cumulative R² ) and provided the best fit with the genetic model [Bayesian information criterion (BIC)]. Compared with the multi-SNP approach, the single-SNP and haplotype-based approaches were relatively similar in terms of cumulative R² and BIC, with an improvement with the haplotype-based approach. Despite limited overlap between detected QTLs, each approach discovered QTLs that had been validated previously, suggesting that each approach can uncover a different subset of QTLs. An integrated GWAS procedure, considering single-locus and multilocus GWAS approaches jointly, may improve the capacity of association studies to detect key QTLs and to provide a more complete picture of the genetic architecture of complex traits in barley.

Identification of genomewide single-nucleotide polymorphisms associated with presummer, summer and autumn bolls in upland cotton

Presummer, summer, and autumn bolls (PSB, SB and AB, respectively) in cotton are related to both maturity and yield. Therefore, studying their genetic basis is important for breeding purposes. In this study, we developed an association analysis panel consisting of 169 upland cotton accessions. The panel was phenotyped for PSB, SB and AB across four environments and genotyped using a Cotton SNP80K array. Single-nucleotide polymorphisms (SNPs) associated with these three traits were identified by a genomewide association study. A total of 53,848 high-quality SNPs were screened, and 91 significant trait-associated SNPs were detected. Of the 91 SNPs 33 were associated with PSB, 21 with SB and 37 with AB. Three SNPs for PSB (TM10410, TM13158 and TM21762) and five for AB (TM13730, TM13733, TM13834, TM29666 and TM43214) were repeatedly detected in two environments or by two methods. These eight SNPs exhibited high phenotypic variation of more than 10%, thus allowing their use formarker-assisted selection. The candidate genes for target traits were also identified. These findings provide a theoretical basis for the improvement of early maturity and yield in cotton breeding programmes.

Genome-Wide Association Studies for Pasmo Resistance in Flax (Linum usitatissimum L.)

Pasmo is one of the most widespread diseases threatening flax production. To identify genetic regions associated with pasmo resistance (PR), a genome-wide association study was performed on 370 accessions from the flax core collection. Evaluation of pasmo severity was performed in the field from 2012 to 2016 in Morden, MB, Canada. Genotyping-by-sequencing has identified 258,873 single nucleotide polymorphisms (SNPs) distributed on all 15 flax chromosomes. Marker-trait associations were identified using ten different statistical models. A total of 692 unique quantitative trait nucleotides (QTNs) associated with 500 putative quantitative trait loci (QTL) were detected from six phenotypic PR datasets (five individual years and average across years). Different QTNs were identified with various statistical models and from individual PR datasets, indicative of the complementation between analytical methods and/or genotype × environment interactions of the QTL effects. The single-locus models tended to identify large-effect QTNs while the multi-loci models were able to detect QTNs with smaller effects. Among the putative QTL, 67 had large effects (3-23%), were stable across all datasets and explained 32-64% of the total variation for PR in the various datasets. Forty-five of these QTL spanned 85 resistance gene analogs including a large toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (TNL) type gene cluster on chromosome 8. The number of QTL with positive-effect or favorite alleles (NPQTL) in accessions was significantly correlated with PR (R 2 = 0.55), suggesting that these QTL effects are mainly additive. NPQTL was also significantly associated with morphotype (R 2 = 0.52) and major QTL with positive effect alleles were present in the fiber type accessions. The 67 large effect QTL are suited for marker-assisted selection and the 500 QTL for effective genomic prediction in PR molecular breeding.

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.

TILLER ANGLE CONTROL 1 modulates plant architecture in response to photosynthetic signals

Light serves as an important environmental cue in regulating plant architecture. Previous work had demonstrated that both photoreceptor-mediated signaling and photosynthesis play a role in determining the orientation of plant organs. TILLER ANGLE CONTROL 1 (TAC1) was recently shown to function in setting the orientation of lateral branches in diverse plant species, but the degree to which it plays a role in light-mediated phenotypes is unknown. Here, we demonstrated that TAC1 expression was light dependent, as expression was lost under continuous dark or far-red growth conditions, but did not drop to these low levels during a diurnal time course. Loss of TAC1 in the dark was gradual, and experiments with photoreceptor mutants indicated this was not dependent upon red/far-red or blue light signaling, but partially required the signaling integrator CONSTITUTIVE PHOTOMORPHOGENESIS 1 (COP1). Overexpression of TAC1 partially prevented the narrowing of branch angles in the dark or under far-red light. Treatment with the carotenoid biosynthesis inhibitor norflurazon or the PSII inhibitor 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU) led to loss of TAC1 expression similar to dark or far-red conditions, but expression increased in response to the PSI inhibitor paraquat. Treatment of adult plants with norflurazon resulted in upward growth angle of branch tips. Our results indicate that TAC1 plays an important role in modulating plant architecture in response to photosynthetic signals.

Alteration of TAC1 expression in Prunus species leads to pleiotropic shoot phenotypes

Prunus persica (peach) trees carrying the "Pillar" or "Broomy" trait (br) have vertically oriented branches caused by loss-of-function mutations in a gene called TILLER ANGLE CONTROL 1 (TAC1). TAC1 encodes a protein in the IGT gene family that includes LAZY1 and DEEPER ROOTING 1 (DRO1), which regulate lateral branch and root orientations, respectively. Here we found that some of the native TAC1 alleles in the hexaploid plum species Prunus domestica, which has a naturally more upright stature, contained a variable length trinucleotide repeat within the same exon 3 region previously found to be disrupted in pillar peach trees. RNAi silencing of TAC1 in plum resulted in trees with severely vertical branch orientations similar to those in pillar peaches but with an even narrower profile. In contrast, PpeTAC1 overexpression in plum led to trees with wider branch angles and more horizontal branch orientations. Pillar peach trees and transgenic plum lines exhibited pleiotropic phenotypes, including differences in trunk and branch diameter, stem growth, and twisting branch phenotypes. Expression profiling of pillar peach trees revealed differential expression of numerous genes associated with biotic and abiotic stress, hormone responses, plastids, reactive oxygen, secondary, and cell wall metabolism. Collectively, the data provide important clues for understanding TAC1 function and show that alteration of TAC1 expression may have broad applicability to agricultural and ornamental tree industries.

Association Mapping Analysis of Fatty Acid Content in Different Ecotypic Rapeseed Using mrMLM

Brassica napus L. is a widely cultivated oil crop and provides important resources of edible vegetable oil, and its quality is determined by fatty acid composition and content. To explain the genetic basis and identify more minor loci for fatty acid content, the multi-locus random-SNP-effect mixed linear model (mrMLM) was used to identify genomic regions associated with fatty acid content in a genetically diverse population of 435 rapeseed accessions, including 77 winter-type, 55 spring-type, and 303 semi-winter-type accessions grown in different environments. A total of 149 quantitative trait nucleotides (QTNs) were found to be associated with fatty acid content and composition, including 34 QTNs that overlapped with the previously reported loci, and 115 novel QTNs. Of these, 35 novel QTNs, located on chromosome A01, A02, A03, A05, A06, A09, A10, and C02, respectively, were repeatedly detected across different environments. Subsequently, we annotated 95 putative candidate genes by BlastP analysis using sequences from Arabidopsis thaliana homologs of the identified regions. The candidate genes included 34 environmentally-insensitive genes (e.g., CER4, DGK2, KCS17, KCS18, MYB4, and TT16) and 61 environment-sensitive genes (e.g., FAB1, FAD6, FAD7, KCR1, KCS9, KCS12, and TT1) as well as genes invloved in the fatty acid biosynthesis. Among these, BnaA08g08280D and BnaC03g60080D differed in genomic sequence between the high- and low-oleic acid lines, and might thus be the novel alleles regulating oleic acid content. Furthermore, RT-qPCR analysis of these genes showed differential expression levels during seed development. Our results highlight the practical and scientific value of mrMLM or QTN detection and the accuracy of linking specific QTNs to fatty acid content, and suggest a useful strategy to improve the fatty acid content of B. napus seeds by molecular marker-assisted breeding.

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

A high-density SNP map was constructed and several novel QTL for branch angle across six environments in Brassica napus were identified. Branch angle is a major determinant for the ideotype of a plant, while the mechanisms underlying this trait in Brassica napus remain elusive. Herein, we developed one doubled haploid population from a cross involving one Capsella bursa-pastoris derived B. napus intertribal introgression line with the compressed branches and wooden stems, and constructed a high-density SNP map covering the genetic distance of 2242.14 cM, with an average marker interval of 0.73 cM. After phenotypic measurements across six environments, the inclusive composite interval mapping algorithm was conducted to analyze the QTL associated with branch angle. In single-environment analysis, a total of 17 QTL were detected and mainly distributed on chromosomes A01, A03, A09 and C03. Of these, three major QTL, qBA.A03-2, qBA.C03-3 and qBA.C03-4 were steadily expressed, each explaining more than 10% of the phenotypic variation in at least two environments. Compared with other results on rapeseed branch angle, these major QTL were newly detected. In QTL by environment interactions (QEI) mapping, 10 QTL were identified, and the QTL average effect and QEI effect were estimated. Of these, 7 QTL were detected in both single-environment analysis and QEI mapping. Based on the physical positions of SNPs and the functional annotation of the Arabidopsis thaliana genome, 27 genes within the QTL regions were selected as candidate genes, including early auxin-responsive genes, small auxin-up RNA, auxin/indoleacetic acid and gretchenhagen-3. These results may pave the way for deciphering the genetic control of branch angle in B. napus.

Single-Locus and Multi-Locus Genome-Wide Association Studies in the Genetic Dissection of Fiber Quality Traits in Upland Cotton (Gossypium hirsutum L.)

A major breeding target in Upland cotton (Gossypium hirsutum L.) is to improve the fiber quality. To address this issue, 169 diverse accessions, genotyped by 53,848 high-quality single-nucleotide polymorphisms (SNPs) and phenotyped in four environments, were used to conduct genome-wide association studies (GWASs) for fiber quality traits using three single-locus and three multi-locus models. As a result, 342 quantitative trait nucleotides (QTNs) controlling fiber quality traits were detected. Of the 342 QTNs, 84 were simultaneously detected in at least two environments or by at least two models, which include 29 for fiber length, 22 for fiber strength, 11 for fiber micronaire, 12 for fiber uniformity, and 10 for fiber elongation. Meanwhile, nine QTNs with 10% greater sizes (R²) were simultaneously detected in at least two environments and between single- and multi-locus models, which include TM80185 (D13) for fiber length, TM1386 (A1) and TM14462 (A6) for fiber strength, TM18616 (A7), TM54735 (D3), and TM79518 (D12) for fiber micronaire, TM77489 (D12) and TM81448 (D13) for fiber uniformity, and TM47772 (D1) for fiber elongation. This indicates the possibility of marker-assisted selection in future breeding programs. Among 455 genes within the linkage disequilibrium regions of the nine QTNs, 113 are potential candidate genes and four are promising candidate genes. These findings reveal the genetic control underlying fiber quality traits and provide insights into possible genetic improvements in Upland cotton fiber quality.