Next Article in Journal
The Role of the NLRP3 Inflammasome and Programmed Cell Death in Acute Liver Injury
Previous Article in Journal
The Search of Association of HLA Class I and Class II Alleles with COVID-19 Mortality in the Russian Cohort
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 4
10.3390/ijms24043069
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fine Mapping of a Pleiotropic Locus (BnUD1) Responsible for the Up-Curling Leaves and Downward-Pointing Siliques in Brassica napus

by
Mao Yang
†,
Jun Chen
†,
Yuqing Chang
,
Shubei Wan
,
Zisu Zhao
,
Fei Ni
and
Rongzhan Guan
*
State Key Laboratory for Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing 210095, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(4), 3069; https://doi.org/10.3390/ijms24043069
Submission received: 22 December 2022 / Revised: 16 January 2023 / Accepted: 17 January 2023 / Published: 4 February 2023
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Leaves and siliques are important organs associated with dry matter biosynthesis and vegetable oil accumulation in plants. We identified and characterized a novel locus controlling leaf and silique development using the Brassica napus mutant Bnud1, which has downward-pointing siliques and up-curling leaves. The inheritance analysis showed that the up-curling leaf and downward-pointing silique traits are controlled by one dominant locus (BnUD1) in populations derived from NJAU5773 and Zhongshuang 11. The BnUD1 locus was initially mapped to a 3.99 Mb interval on the A05 chromosome with a BC6F2 population by a bulked segregant analysis-sequencing approach. To more precisely map BnUD1, 103 InDel primer pairs uniformly covering the mapping interval and the BC5F3 and BC6F2 populations consisting of 1042 individuals were used to narrow the mapping interval to a 54.84 kb region. The mapping interval included 11 annotated genes. The bioinformatic analysis and gene sequencing data suggested that BnaA05G0157900ZS and BnaA05G0158100ZS may be responsible for the mutant traits. Protein sequence analyses showed that the mutations in the candidate gene BnaA05G0157900ZS altered the encoded PME in the trans-membrane region (G45A), the PMEI domain (G122S), and the pectinesterase domain (G394D). In addition, a 573 bp insertion was detected in the pectinesterase domain of the BnaA05G0157900ZS gene in the Bnud1 mutant. Other primary experiments indicated that the locus responsible for the downward-pointing siliques and up-curling leaves negatively affected the plant height and 1000-seed weight, but it significantly increased the seeds per silique and positively affected photosynthetic efficiency to some extent. Furthermore, plants carrying the BnUD1 locus were compact, implying they may be useful for increasing B. napus planting density. The findings of this study provide an important foundation for future research on the genetic mechanism regulating the dicotyledonous plant growth status, and the Bnud1 plants can be used directly in breeding.

1. Introduction

Leaf morphology affects dicotyledonous plant photosynthetic activities and the accumulation of dry matter. The slight up-curling of leaves may increase the overall light energy use efficiency and improve plant yields. Leaf shape formation is a complex developmental process involving leaf primordium formation, polarity establishment, and cell differentiation. The curled leaf phenotype is caused by mutations to genes related to leaf development [1,2], including genes encoding the transcription factors Class III HOMEODOMAIN LEUCINE-ZIPPER (HD-Zip III) [3,4,5], KANADI (KAN) [6,7,8], WUSCHEL RELATED HOMEOBOX (WOX) [9], and TB1-CYC-PCF (TCP) [10]. These transcription factors modulate leaf polarity establishment and cause leaves to curl upward by regulating the expression of leaf asymmetry development-related genes, including ASYMMETRIC LEAVES 1 (AS1) and AS2 [11]. Plant hormone biosynthesis and signal transduction is another major factor influencing leaf shapes. For example, auxin response factors (e.g., ARF3) can induce leaf curling by affecting the mutual antagonism between HD-ZIP III and KAN transcription factors [12]. The UCU1 gene, which encodes the SHAGGY/GSK3 protein, is involved in auxin and Brassinosteroid (BR) signal transduction. A mutation to this gene can lead to the formation of leaves that curl downward as well as dwarfism [13]. Both miRNA165 and miRNA166 affect leaf development by regulating HD-ZIP III gene expression [4,14,15], whereas miRNA160 causes leaves to curl by regulating the expression of the auxin-responsive genes ARF10 and ARF17 [16,17,18]. In addition, a few non-coding RNAs reportedly contribute to the establishment of adaxial-abaxial polarity, which may lead to the development of curled leaves [19]. The leaf curling-related genes have also been associated with cutin and cuticular wax production and cell wall synthesis.
Siliques are specific to plants in the family Brassicaceae, which consists of numerous important species, including Arabidopsis thaliana, Brassica napus, and Brassica rapa. Siliques are organs that store seeds, which contain oil and protein. Silique traits, such as length, diameter, angle (upright, downward pointing, or oblique), seed number, stalk, and fruit shaft, directly affect the seed yield [20,21]. Accordingly, the genes controlling silique growth and development must be identified and functionally annotated. Previous studies on silique traits focused on the seed number per silique, silique length, silique number per plant, 1000-seed weight of siliques, and silique dehiscence resistance [22,23,24,25,26]. However, the genetic basis for the silique angle remains unknown. An earlier study on the A. thaliana brevipedicellus mutant revealed that a mutation to the AtBP gene encoding a homeodomain protein (KNAT1) is responsible for the formation of a downward-pointing pedicel and flower [27]. The AtBP gene is expressed in the peripheral zone of the shoot apical meristem and in the cortical cell layers of the inflorescence stem (peduncle) and pedicel [28]. A mutation to the BP gene alters its asymmetrical expression patterns in the pedicel abaxial and adaxial sides, thereby severely affecting pedicel cell differentiation, elongation, and growth, leading to downward-pointing pedicels [27].
Cell and organ shapes in plants are dependent on the chemical structure and mechanical properties of the extracellular matrix. The pectin matrix, which is the main load-bearing plant cell wall component, is constantly re-modeled to enable plant morphological development. This re-modeling is regulated by several loosening and stiffening agents, including pectin methylesterase (PME) and calcium ions. Specifically, PMEs are ubiquitous enzymes that can catalyze the demethylesterification of homogalacturonan (HG) [29]. The demethylesterified HG can either form Ca2+ bonds, which promote the formation of the so-called ‘egg-box’ model structure that solidifies the cell wall, or it serves as the substrate of other pectin-degrading enzymes that loosen the cell wall [30,31]. The continuous and precise spatiotemporal regulation of the activities of these agents is necessary for proper morphogenesis at the cell and tissue levels [32,33].
The down-regulated expression of PME genes may decrease the extent of the pectin demethylation in the cell wall and restrict the binding between Ca2+ and PMEs, leading to severe decreases in the cell wall texture and rigidity, which may substantially alter the plant organ phenotype. However, relatively few PME genes have been functionally characterized. In A. thaliana, a PME may also contribute to the downward-pointing flower and silique phenotype [34,35]. The overexpression of pectin methylesterase inhibitor (PMEI)-encoding genes and the down-regulated expression of PME genes in wild-type A. thaliana result in plants with curled leaves, a convoluted shoot, and misshapen siliques [34].
Oilseed rape (Brassica napus L.) is one of the most important oil crops worldwide because it is a source of high-quality edible oil that is consumed by humans, while also serving as protein-rich feed for animals and a raw material for industrial processes. As the global population continues to increase, there is a growing demand for higher yielding oilseed rape cultivars. The leaf and silique morphological traits can directly affect photosynthesis and rapeseed yields. Hence, studies on rapeseed mutants with abnormal leaf and silique traits may be relevant for the breeding of dicotyledonous plants with enhanced characteristics and future research on silique developmental biology.
Leaf morphological mutations in B. napus (e.g., upward and downward curling and wrinkling) are related to photosynthetic efficiency and affect the yield. Previous research on leaf curling revealed the associated loci and genes, including a dominant locus (BnDWF/DCL1) [36], semi-dominant genes (sca and ds-4) [37,38], and the genes (BnUC1, BnUC2, and BnUC3) [39,40,41] responsible for the up-curling of B. napus leaves. Analyses of gene functions demonstrated that abnormalities in some leaf flatness-related traits are caused by the deficient expression of genes encoding proteins mediating phytohormone (e.g., auxin and brassinolide) biosynthesis or signaling [13,42,43,44]. Defects in genes directly affecting phytohormone regulatory pathways usually result in severe morphological changes that lead to abnormal plant development, with potential implications for breeding. However, identifying loci associated with leaf morphological variations may help developmental biologists and breeders attempting to develop new varieties with enhanced photosynthetic activities.
In the present study, we analyzed a B. napus cv. NJAU5773 mutant (Bnud1) that had downward-pointing siliques and up-curled leaves using a genetic and molecular approach. We fine mapped a dominant locus (BnUD1), identified candidate genes for the mutated pleiotropic traits on the basis of sequencing and gene expression experiments, and explored the effects of BnUD1 on agronomic traits. Our findings may serve as the foundation for the breeding of B. napus lines with potentially improved morphological characteristics and may be useful for elucidating the genetic mechanism underlying leaf and silique development.

2. Results

2.1. Performance of the Bnud1 Mutant

Brassica napus NJAU5773 (named Bnud1 in this study) has up-curling leaves at the seedling stage, which is in contrast to the normal flat leaves of the canola variety Zhongshuang 11 (ZS11) (Figure 1a). Additionally, NJAU5773 has downward-pointing flowers that develop into downward-pointing siliques (Figure 1b,c). Plants with up-curling leaves and downward-pointing siliques exhibit semi-dwarfism in a segregating population derived from NJAU5773. To examine the Bnud1 mutant, near isogenic lines (ZS11-UD1 and ZS11) were developed via marker-assisted selection.

2.2. Inheritance of the Up-Curling Leaf and Downward-Pointing Silique Traits

In the current study, NJAU5773 was reciprocally crossed with canola variety ZS11, which has a sequenced genome (http://brassicadb.org/brad/, accessed on 15 January 2023), to generate F1 plants and progeny populations for genetic analyses. Investigations involving segregating populations, which included F2, and successively backcrossed populations indicated that the segregation ratio (i.e., plants with up-curling leaves and downward-pointing siliques vs. wild-type plants) was consistent with the Mendel inheritance ratio (as suggested by the Chi-square test) (Table 1). Accordingly, the up-curling leaf and downward-pointing silique traits appear to be controlled by a dominant locus (i.e., BnUD1).

2.3. Mapping of the BnUD1 Locus

To map the BnUD1 locus, a BC6F2 population consisting of 308 plants was derived from the cross between the two parents (NJAU5773 and ZS11). Thirty plants with up-curling leaves and downward-pointing siliques and 30 plants with flat leaves and upright siliques were selected from the BC6F2 population to scan the BnUD1 locus via a bulked segregant analysis (BSA) approach. The two pooled samples were sequenced, with 30.06× sequence coverage for plants with the BnUD1 locus and 21.65× sequence coverage for plants lacking the BnUD1 locus. The clean reads were aligned to the B. napus cv. ZS11 reference genome (ZS11-v20200127; http://cbi.hzau.edu.cn/cgi-bin/rape/download_ext, accessed on 15 January 2023). This helped to identify the single nucleotide polymorphisms (SNPs) between the two pools. A total of 444,519 SNPs and 79,892 small insertion/deletions (InDels) were detected (Table S1). The ΔSNP index was calculated and the following two segments of the A05 chromosome were considered as candidate BnUD1-harboring regions: 7,094,487–11,086,188 bp (approximately 3.99 Mb) and 37,261,266–37,670,607 bp (approximately 0.41 Mb) (Figure S1; Table S2).
To more precisely map the BnUD1 locus, 49 primer pairs for InDel markers were designed to uniformly cover the preliminary mapping interval on the basis of the BSA sequencing results. Fifteen markers targeting the smaller candidate region on the A05 chromosome (37,261,266–37,670,607 bp) were polymorphic (i.e., ID-62, ID-63, ID-64, ID-65, ID-66, ID-67, ID-68, ID-70, ID-73, ID-74, ID-77, ID-78, ID-79, ID-81, and ID82) (Table S3), whereas 12 polymorphic InDel markers were obtained for the larger candidate region on the A05 chromosome (7,094,487–11,086,188 bp; i.e., InDel-2, InDel-3, InDel-5, InDel-6, InDel-7, InDel-8, InDel-9, InDel-11, InDel-12, InDel-13, InDel-15, and InDel-16) (Figure 2a; Table S3). These polymorphic markers were used to genotype 1042 individuals in the BC5F3 and BC6F2 populations. The mapping using the JoinMap 4.1 software indicated that the BnUD1 locus was located in a mapping interval between InDel-13 and InDel-15 (1.63 Mb) on the A05 chromosome (Figure 2a), which eliminated one of the candidate regions (37,261,266–37,670,607 bp). To fine map the BnUD1 locus, 29 InDel marker primers were designed to target the narrowed mapping interval. Of these markers, 10 were polymorphic (InDel-47, InDel-50, InDel-55, InDel-56, InDel-57, InDel-65, InDel-67, InDel-70, InDel-72, and InDel-73). The analysis of the segregating population using these polymorphic markers generated a 466.58 kb mapping interval (Figure 2b), which was used to develop four InDel markers (InDel-75, InDel-78, InDel-80, and InDel-81). With these four markers, the BnUD1 locus was mapped to a 54.84 kb interval between InDel-78 and InDel-15 (Figure 2c and Figure 3). The other markers developed according to the BSA sequencing results were unable to further narrow the mapping interval.

2.4. Gene Cloning

The homologous segment sequences in the fine mapping interval were downloaded from the B. napus cv. ZS11 genome database. The genes in the interval were annotated on the basis of the A. thaliana genome to identify the candidate genes associated with the up-curling leaf and downward-pointing silique traits. The interval harbored 11 annotated genes (Table 2). To analyze the candidate genes in the BnUD1 locus, we cloned the 11 genes from the two parents of the mapping populations (NJAU5773 and ZS11). The subsequent alignment of these sequences indicated that BnaA05G0157900ZS, BnaA05G0158100ZS, and BnaA05G0158300ZS differed between the two parents (Figure S2), whereas the other examined genes within the mapping interval were identical in the two parents. The comparison of the BnaA05G0157900ZS sequences revealed that the encoded protein in NJAU5773 has three amino acid (AA) substitutions, of which the Gly-to-Ala mutation (G45A) occurred in the trans-membrane region, the Gly-to-Ser mutation (G122S) was detected in the PMEI domain, and the Gly-to-Asp mutation (G394D) was present in the pectinesterase domain (Figure S2a). In addition, a 573 bp segment was inserted into the pectinesterase domain of BnaA05G0157900ZS gene in the Bnud1 mutant. Next, we developed an InDel marker (BnA05ID1) specific for this insertion (Table S3). The marker co-segregated with the up-curling leaf and downward-pointing silique phenotypes in the segregating populations (Figure 4). These sequence alterations may affect leaf and silique trait formation in the Bnud1 mutant.
The protein encoded by BnaA05G0158100ZS in the Bnud1 mutant included one substitution (Val-to-Asp) at AA position 146, which is in the RING domain of a RING/U-box superfamily protein (Figure S2b). This substitution may lead to a change in function. The examination of BnaA05G0158300ZS in the Bnud1 mutant detected three AA substitutions, but they were not located in the IENR2 domain encoded specifically by this gene (Figure S2c). The results of the sequence analyses suggested that BnaA05G0157900ZS and BnaA05G0158100ZS are the most likely candidate genes in the BnUD1 locus.

2.5. Candidate Gene Analysis

To determine whether the gene mutations are associated with the mutant traits, we further examined the gene functions by conducting a bioinformatic analysis as well as gene expression experiments using a pair of isogenic lines with the ZS11 genetic background. The BnaA05G0157900ZS gene is homologous to AT1G53840, which encodes PECTIN METHYLESTERASE 1 (PME1). The PMEs can catalyze the specific demethylesterification of HG to form load-bearing Ca2+ crosslinks that affect the texture and rigidity of the cell wall [29,30,31]. Thus, they are reportedly involved in various biological processes, including cell wall expansion, seed germination, and hypocotyl elongation [45,46,47]. The overexpression of PMEI-encoding genes in wild-type A. thaliana plants leads to the production of curled leaves, a convoluted shoot, and downward-pointing siliques [34,35]. In the present study, the four sequence changes in the Bnud1 mutant BnaA05G0157900ZS gene were located in the gene-coding region or intron region. On the basis of the previously reported potential function of the protein encoded by BnaA05G0157900ZS, these sequence mutations are probably responsible for the observed mutated traits.
The mutation in the Bnud1 mutant BnaA05G0158100ZS gene was located in the sequence encoding the RING/U-box domain. This gene is homologous to AT1G53820, which encodes the biological stress-responsive RING/U-box superfamily member ARABIDOPSIS TOXICOS EN LEVADURA 60 (AL60). This gene has not been linked to substantial alterations to leaf and silique morphological characteristics [48,49,50,51,52]. The BnaA05G0158300ZS gene in Bnud1 has three mutations that are not within the sequence encoding the conserved IENR2 domains. Moreover, BnaA05G0158300ZS is homologous to AT1G53800, which was annotated as a gene involved in sarcomeric titin assembly during cardiac myofibrillogenesis in animals [53,54]. However, AT1G53800 homologs in plants have not been investigated.
The expression levels of the three mutated genes in the mapping interval were determined by performing a quantitative real-time polymerase chain reaction (qRT-PCR) analysis of leaf samples from the four plants with the Bnud1 mutated traits and the four plants with the wild-type traits in the BC6F2 population. Both BnaA05G0157900ZS and BnaA05G0158100ZS were expressed at significantly lower levels in the plants with the Bnud1 mutated traits than in the plants with the wild-type traits (Figure 5). In contrast, the BnaA05G0158300ZS expression level did not differ between the two plant types. This implied that BnaA05G0158300ZS is not associated with the formation of mutated traits (Figure 5). These results suggested that BnaA05G0157900ZS and BnaA05G0158100ZS may be responsible for the up-curling leaf and downward-pointing silique traits of the Bnud1 mutant.

2.6. Agronomic Traits

To evaluate the effects of the BnUD1 locus on plant agronomic traits, 30 plants with the Bnud1 mutated traits and 30 plants with the wild-type traits were randomly sampled from the BC6F2 population derived from the cross between NJAU5773 and the recurrent parent ZS11. The values for some of the agronomic traits, including plant height, branch height, main inflorescence length, number of first effective branches, and 1000-seed weight, were significantly lower for the plants with up-curling leaves and downward-pointing siliques than for the plants with flat leaves and upright siliques (Table 3). The other agronomic traits, including stem diameter, siliques of the main inflorescence, total siliques per plant, and silique length, did not differ between the plants with and without the BnUD1 locus (Table 3). Thus, the BnUD1 locus appeared to negatively affect the plant stature, resulting in a compact architecture. Additionally, the number of seeds per silique was significantly higher for the plants with the BnUD1 locus than for the plants without the BnUD1 locus, reflecting the positive effects of the BnUD1 locus on the seed yield.

2.7. Determination of the Chlorophyll Content and Photosynthetic Efficiency

Compared with the plants with the wild-type traits at the seedling stage, the leaf chlorophyll (Chl) a, Chl b, and total Chl contents as well as the Chl a/b ratio were significantly higher for the plants with the Bnud1 mutated traits selected from the BC6F2 population derived from the cross between the NJAU5773 and the recurrent parent ZS11 (Table 4). This result indicated that the Bnud1 mutated traits were associated with increases in the leaf Chl content. The leaf net photosynthetic rate, stomatal conductance, and concentration of intercellular CO2 were significantly higher in the plants with the BnUD1 locus than in the plants without the BnUD1 locus. However, the leaf transpiration rate did not differ between the plants with the BnUD1 locus and the plants without the BnUD1 locus (Table 5). Hence, the BnUD1 locus may lead to increased photosynthetic efficiency.

3. Discussion

Leaves and siliques are important photosynthesis-related organs that influence the agronomic value of crops. Up-curling leaves and downward-pointing siliques are typically the result of the heteromorphic development of plant organs. However, the genetic mechanism underlying the formation of up-curling leaves is generally unrelated to the genetic basis of silique formation.
There has recently been some progress in the research on the genes related to leaf curling in A. thaliana and rice [1,6,7,9,55]. For example, the MYB transcription factor-encoding PHANTASTICA (PHAN) gene was first reported to be associated with leaf adaxial-abaxial polarity in Antirrhinum majus [56,57]. In A. thaliana, AS1 is homologous to PHAN and regulates leaf polarity by forming a protein complex with the plant-specific lateral organ boundaries (LOB) family protein AS2 [11,55]. Several transcription factors (e.g., HD-ZIP III, KAN, and YABBY) can bind to the AS1 gene promotor to regulate the establishment of leaf abaxial polarity in A. thaliana, resulting in leaf-rolling behavior [4,58,59,60,61,62]. Auxin response factors, such as ARF3, ARF4, and ARF2, can form a complex with KAN to modulate leaf abaxial polarity [12,63,64]. MicroRNAs, including miR165 and miR166 targeting HD-ZIP III [62,65] and miR390 triggering the production of phasi-RNAs from TAS3 trans-acting short interfering RNA transcripts, also influence leaf-rolling behaviors by repressing ARF activities [66,67].
Silique morphological development does not generally depend on leaf developmental genes. Previous research on silique formation focused on silique yield-associated traits, such as seed number and weight within siliques and silique length. Some genes regulating silique development have been identified, including the silique length-related gene REPLUMLESS (RPL) [68], silique shattering-related genes SHATTERPROOF1 (SHP1) and SHP2 [69], INDEHISCENT (IND) [70], ALCATRAZ (ALC) [71], and silique polarity-associated KNOX genes [27]. To date, how silique angles form remains unknown. In A. thaliana, the KNOX gene BREVIPEDICELLUS (BP), which is expressed downstream of the auxin and AS1 pathways, may encode a protein that regulates silique polarity signals by inhibiting auxin transport, thereby affecting silique and pedicel morphological development [27,72].
Mutations that alter the ability of PMEs to catalyze the demethylesterification of the cell wall HG may modify cell wall components, resulting in up-curling leaves and silique and stem morphological abnormalities [29,34,35]. In this study, a PME gene contributing to the up-curling leaf and downward-pointing silique traits was identified (Figure 2 and Figure S2a). More specifically, the marker analysis revealed that mutations in the PME gene BnaA05G0157900ZS may help to explain the formation of downward-pointing siliques and up-curling leaves (Figure 3 and Figure 4). Another candidate gene (BnaA05G0158100ZS) for the mutated traits encodes a RING/U-box-containing protein associated with biological stress responses (Table 2). Although this gene may be part of a complex genetic regulatory system, whether it is directly associated with the mutated traits is unclear.
The mutant NJAU5773, which was originally identified in a breeding population and obtained via consecutive generations of selfing, has three mutated genes in the mapping interval (Figure S2). When the two candidate genes for the BnUD1 locus-associated traits were used as queries to screen for the corresponding gene sequences in the B. napus pan-genome database (http://cbi.hzau.edu.cn/bnapus/), the mutated candidate genes were undetectable in other genomes. The segment harboring the BnUD1 locus is located in the A sub-genome of B. napus, similar to the homologous segment in the A genome of B. rapa. However, the candidate genes differed from the corresponding homologous genes. Therefore, the candidate gene mutations are unique to NJAU5773 and may be the product of natural evolutionary events.
Moderately up-curled leaves and downward-pointing siliques can theoretically increase light transmittance and the light saturation point, thereby increasing the overall photosynthetic efficiency. The increased planting density associated with appropriately up-curled leaves and downward-pointing siliques may also positively affect the harvest index [73,74]. Thus, the up-curling leaf trait should be explored in more detail. To date, three loci related to the up-curling leaf trait have been identified in B. napus (i.e., BnUC1, BnUC2, and BnUC3) [39,40,41]. In present work, we found that the Bnud1 mutant had an increased photosynthetic efficiency at seedling stage, and relatively small plant architecture (Figure 1; Table 5). The BnUD1 locus is expected to be applied to increase the efficiency of leaf photosynthetic activities and decrease plant height (Table 3 and Table 5), but the underlying mechanisms will need to be characterized in future investigations. Unlike the other loci mediating the up-curling of leaves, the BnUD1 locus also controls the formation of downward-pointing siliques and increases in the number of seeds per silique (Table 3). However, the reason for the increase in the number of seeds per silique is unclear. We speculate that the influx of carbon is enhanced in the downward-pointing siliques. Alternatively, the increase in seed production may be the result of increased pollination efficiency.

4. Materials and Methods

4.1. Plant Materials

The double-low B. napus (oilseed rape) line NJAU5773 with up-curling leaves (before the budding stage) and downward-pointing siliques obtained from our germplasm and canola variety ZS11 provided by Nanjing Agricultural University were used as the parents to produce the F1 population. The F1 individuals were selfed to generate F2 mapping populations and backcrossed with the recurrent parent ZS11 (female) to construct the mapping populations. The selfed and backcrossed populations were examined to calculate the segregation ratio of plants with up-curling leaves and downward-pointing siliques to plants with flat leaves and upright siliques. The BC5F3 and BC6F2 populations were used for the preliminary and fine mapping of the BnUD1 locus. The plants with up-curling leaves and downward-pointing siliques and the plants with flat leaves and upright siliques in the BC6F2 population were used for the BSA, qRT-PCR, and analyses of the Chl content, photosynthetic efficiency, and agronomic traits.
All materials were grown on the research farm of Nanjing Agricultural University (Nanjing, China). Plants were cultivated in 2.5 m rows, with 15 plants per row and 0.4 m between rows.

4.2. Genetic Analysis

To determine the number of genes controlling the up-curling leaf and downward-pointing silique phenotype of the Bnud1 mutant, all generations, including the F1, F2, and BC1 populations and the derived lines (BC1-BC6, BC5F3, and BC6F2), were grown in the field. Leaf and silique morphology were assessed at the seedling and maturity stages, respectively. Chi-square tests were performed using the segregation data in each population to analyze the genetic regulation of the up-curling leaf and downward-pointing silique traits.

4.3. Bulked Segregant Analysis

To map the BnUD1 locus, a BC6F2 population was developed from the cross between the Bnud1 mutant with up-curling leaves and downward-pointing siliques and wild-type plants with flat leaves and upright siliques. Genomic DNA was extracted from young leaves using cetyl-trimethylammonium bromide (CTAB). The BC6F2 family population was identified according to a BSA. Equal amounts of DNA from 30 plants with up-curling leaves and downward-pointing siliques and 30 wild-type plants from the BC6F2 population were pooled to form the Bnud1 mutant trait bulk (UDB) and the wild-type trait bulk (WTB), respectively. The polymorphisms between the bulks (UDB and WTB) were screened using InDel markers from a previous study [75].
To efficiently develop linked markers in the target region in B. napus, a BSA sequencing (BSA-seq) experiment was performed. Genomic DNA from the different bulks was subjected to a whole-genome sequencing analysis. Short-insert (350–450 bp) sequencing libraries were constructed from approximately 2 μg parental genomic DNA using the TruSeq® DNA Sample Preparation Kit (Illumina, San Diego, CA, USA). The quantified libraries were sequenced on the HiSeq 3000 platform (Illumina) to produce 150 bp paired-end reads. The InDels and SNPs were called as previously described [76,77]. The de novo assembled ZS11 genome was used as the reference for calculating the SNP index of UDB and WTB. The ΔSNP index was calculated by subtracting the SNP index for UDB from the SNP index for WTB.

4.4. Mapping of the BnUD1 Locus

The BSA-seq results for the mapping interval were used to identify 5345 SNPs/1275 InDels and 721 SNPs/199 InDels covering 3.99 Mb and 0.41 Mb intervals between UDB and WTB (Figure S1; Table S2). Next, 103 differential InDel sites (at positions 37,261,266–37,670,607 bp and 7,094,487–11,086,188 bp of A05) were used to develop molecular markers that could narrow the mapping interval. The InDel markers were used to analyze the BC5F3 and BC6F2 populations to fine map the BnUD1 locus (Table S3). Finally, a fine linkage map for the locus associated with the up-curling leaf and downward-pointing silique traits was constructed using the JoinMap 4.1 software and the polymorphic InDel markers [78]. The primers used are listed in Table S3.
The PCR conditions for the molecular marker experiments were as follows: 94 °C for 5 min; 35 cycles of 94 °C for 30 s, annealing temperature of each InDel marker for 30 s, and 72 °C for 30 s; 72 °C for 10 min.

4.5. Identification of Genes in the Mapping Interval and Comparative Sequencing

The sequences of the fine mapping interval on the B. napus cv. ZS11 A05 chromosome were downloaded from the Brassicaceae database (http://brassicadb.org/brad/, accessed on 15 January 2023) and the B. napus pan-genome information online resource (http://cbi.hzau.edu.cn/bnapus/, accessed on 15 January 2023) to identify the genes in the mapping interval. The genes detected in the mapping interval were annotated on the basis of B. napus cv. ZS11 annotated genes.
Extracted DNA was digested using RNase I (Takara, Dalian, China) to remove RNA. The RNA-free DNA samples were genotyped by genome sequencing. Total RNA was extracted from the leaves at the seedling stage using the RNAprep Pure Plant Kit (BioTeke, Beijing, China). First-strand cDNA was synthesized from the RNA using a reverse transcription kit (Takara, Tokyo, Japan). All 11 genes identified in the mapped interval were cloned from the two parents with gene-specific primers designed using the Primer Premier 5.0 software [79] (Tables S4 and S5). The PCR amplifications were performed as previously described [80]. The amplified fragments were inserted into the pEASY-Blunt Cloning Kit vector (TransGen, Beijing, China) and sequenced. The resulting sequences were aligned using Clustal X1.83 software [81].

4.6. Quantitative Real-Time PCR Analysis

A qRT-PCR analysis was conducted to compare the expression levels of the three mutated genes in the mapping interval between the plants with the Bnud1 mutated traits and the plants with the wild-type traits in the BC6F2 population. Leaves were collected from four plants at the five-leaf stage and then used to prepare the cDNA template for the qRT-PCR analysis, which was completed with gene-specific primers (Table S6). The gene expression levels were normalized against the expression of an BnActin gene (i.e., housekeeping gene) (Table S6). The qRT-PCR was performed using the SYBR Green Real-time PCR Master mix and the CFX96-2 PCR system (Bio-Rad, Hercules, CA, USA) and the relative expression levels were analyzed as described previously [82]. Relative expression levels were calculated according to the 2−ΔΔCt method, with the actin gene serving as the internal control. Four biological replicates were used.

4.7. Agronomic Trait Analysis

To investigate the effects of the BnUD1 locus on plant agronomic traits, 30 plants with the Bnud1 mutated traits and 30 plants with the wild-type traits were randomly selected from the BC6F2 population. The examined agronomic traits included plant height, branch height, main inflorescence length, stem diameter, number of first effective branches, number of siliques on the main inflorescence, total number of siliques per plant, silique length, seeds per silique, and 1000-seed weight. The mean values for all agronomic traits were compared between the plants with Bnud1 mutated traits and the plants with wild-type traits by t-tests.

4.8. Determination of the Chlorophyll Content and Photosynthetic Efficiency

Fifteen homozygous plants with the Bnud1 mutated traits and 15 plants with the wild-type traits were randomly selected from the BC6F2 population at the seedling stage to measure the Chl contents. Specifically, Chl was extracted from 0.2 g fresh leaves using 50 mL 80% acetone, after which the Chl content was determined using the Alpha-1500 spectrophotometer (LASPEC, Shanghai, China). The leaf Chl a, Chl b, and total Chl contents were measured as previously described [83,84].
Six plants with the BnUD1 locus and six plants without the BnUD1 locus were randomly selected from the BC6F2 population at the seedling stage for an analysis of photosynthetic efficiency. The photosynthetic characteristics of the plants were determined using the Li-Cor 6400 portable photosynthesis system (Li-Cor Inc., Lincoln, NE, USA), with the built-in light source set at 1000 μmol photons m−2 s−1 at 23 °C as previously described. All measurements were completed between 9:00 a.m. and 11:00 a.m. [84].

5. Conclusions

A new mutant B. napus plant (NJAU5773) with up-curling leaves and downward-pointing siliques was identified in our B. napus germplasm. Inheritance studies showed that the up-curling leaf and downward-pointing silique traits were controlled by one dominant locus, which was mapped to a 54.84 kb interval on the BnA05 chromosome using BSA-seq and InDel markers. Both BnaA05G0157900ZS and BnaA05G0158100ZS were identified as candidate genes on the basis of sequencing and gene expression analyses. The examination of individual plants revealed the BnUD1 locus had positive effects on photosynthetic efficiency. In conclusion, the findings of this study provide a theoretical foundation for elucidating the mechanism mediating the up-curling leaf and downward-pointing silique traits, which may be relevant for B. napus breeding programs.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24043069/s1.

Author Contributions

R.G. conceived and designed the study. M.Y. and J.C. performed many of the experiments and wrote the manuscript. R.G. advised on the experiments and modified the manuscript. Y.C., S.W., Z.Z. and F.N., took part in the DNA extraction, marker experiments, and sequence analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported financially by the National Natural Science Foundation of China (32171974), and the Fundamental Research Funds for the Central Universities (KYZZ2022003). The funders provided the financial support to the research, but had no role in the design of the study, analysis, interpretations of data, or in writing the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fei, D.; Guan, C.; Jiao, Y. Molecular mechanisms of leaf morphogenesis. Mol. Plant 2018, 11, 1117–1134. [Google Scholar]
  2. Conklin, P.A.; Josh, S.; Shujie, L.; Scanlon, M.J. On the mechanisms of development in monocot and eudicot leaves. New Phytol. 2019, 221, 706–724. [Google Scholar] [CrossRef] [PubMed]
  3. Li, Y.Y.; Shen, A.; Xiong, W.; Sun, Q.L.; Luo, Q.; Song, T.; Li, Z.L.; Luan, W.J. Overexpression of OsHox32 results in pleiotropic effects on plant type architecture and leaf development in Rice. Rice 2016, 9, 46. [Google Scholar] [CrossRef]
  4. Rong, F.; Chen, F.; Huang, L.; Zhang, J.; Zhang, C.; Hou, D.; Cheng, Z.; Weng, Y.; Chen, P.; Li, Y. A mutation in class III homeodomain-leucine zipper (HD-ZIP III) transcription factor results in curly leaf (cul) in cucumber (Cucumis sativus L.). Theor. Appl. Genet. 2018, 132, 113–123. [Google Scholar] [CrossRef]
  5. Merelo, P.; Ram, H.; Caggiano, M.P.; Ohno, C.; Heisler, M.G. Regulation of MIR165/166 by class II and class III homeodomain leucine zipper proteins establishes leaf polarity. Proc. Natl. Acad. Sci. USA 2016, 113, 11973–11978. [Google Scholar] [CrossRef]
  6. Zhang, G.H.; Xu, Q.; Zhu, X.D.; Qian, Q.; Xue, H.W. SHALLOT-LIKE1 is a KANADI transcription factor that modulates Rice leaf rolling by regulating leaf Abaxial cell development. Plant Cell 2009, 21, 719–735. [Google Scholar] [CrossRef]
  7. Huang, T.; Harrar, Y.; Lin, C.; Reinhart, B.; Newell, N.R.; Talaverarauh, F.; Hokin, S.A.; Barton, M.K.; Kerstetter, R.A. Arabidopsis KANADI1 acts as a transcriptional repressor by interacting with a specific cis-element and regulates auxin biosynthesis, transport, and signaling in opposition to HD-ZIPIII factors. Plant Cell 2014, 26, 246. [Google Scholar] [CrossRef]
  8. Merelo, P.; Paredes, E.B.; Heisler, M.G.; Wenkel, S. The shady side of leaf development: The role of the REVOLUTA/KANADI1 module in leaf patterning and auxin-mediated growth promotion. Curr. Opin. Plant Biol. 2017, 35, 111–116. [Google Scholar] [CrossRef]
  9. Cho, S.H.; Yoo, S.C.; Zhang, H.; Pandeya, D.; Koh, H.J.; Hwang, J.Y.; Kim, G.T.; Paek, N.C. The rice narrow leaf2 and narrow leaf3 loci encode WUSCHEL-related homeobox 3A (OsWOX3A) and function in leaf, spikelet, tiller and lateral root development. New Phytol. 2013, 198, 1071–1084. [Google Scholar] [CrossRef] [PubMed]
  10. Koyama, T.; Sato, F.; Ohme-Takagi, M. Roles of miR319 and TCP transcription factors in leaf development. Plant Physiol. 2017, 175, 874–885. [Google Scholar] [CrossRef]
  11. Fu, Y.; Xu, L.; Xu, B.; Yang, L.; Ling, Q.; Wang, H.; Huang, H. Genetic interactions between leaf polarity-controlling genes and ASYMMETRIC LEAVES1 and 2 in Arabidopsis leaf patterning. Plant Cell Physiol. 2007, 48, 724–735. [Google Scholar] [CrossRef] [Green Version]
  12. Pekker, I.; Alvarez, J.P.; Eshed, Y. Auxin response factors mediate Arabidopsis organ asymmetry via modulation of KANADI activity. Plant Cell 2005, 17, 2899–2910. [Google Scholar] [CrossRef] [PubMed]
  13. Pérez-Pérez, J.M.; Ponce, M.R.; Micol, J.L. The UCU1 Arabidopsis Gene Encodes a SHAGGY/GSK3-like Kinase Required for Cell Expansion along the Proximodistal Axis. Dev. Biol. 2002, 242, 161–173. [Google Scholar] [CrossRef]
  14. Prigge, M.J.; Otsuga, D.; Alonso, J.M.; Ecker, J.R.; Drews, G.N.; Clark, S.E. Class III homeodomain-leucine zipper gene family members have overlapping, antagonistic, and distinct roles in Arabidopsis development. Plant Cell 2005, 17, 61–76. [Google Scholar] [CrossRef] [PubMed]
  15. Tatematsu, K.; Toyokura, K.; Miyashima, S.; Nakajima, K.; Okada, K. A molecular mechanism that confines the activity pattern of miR165 in Arabidopsis leaf primordia. Plant J. 2015, 82, 596–608. [Google Scholar] [CrossRef]
  16. Mallory, A.C.; Vaucheret, H. Functions of microRNAs and related small RNAs in plants. Nat. Genet. 2006, 38, 31–36. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, P.P.; Montgomery, T.A.; Fahlgren, N.; Kasschau, K.D.; Nonogaki, H.; Carrington, J.C. Repression of AUXIN RESPONSE FACTOR10 by microRNA160 is critical for seed germination and post-germination stages. Plant J. 2007, 52, 133–146. [Google Scholar] [CrossRef]
  18. Shen, X.X.; He, J.Q.; Ping, Y.K.; Guo, J.X.; Hou, N.; Cao, F.G.; Li, X.W.; Geng, D.L.; Wang, S.C.; Chen, P.X.; et al. The positive feedback regulatory loop of miR160-Auxin Response Factor 17-HYPONASTIC LEAVES 1 mediates drought tolerance in apple trees. Plant Physiol. 2022, 188, 1686–1708. [Google Scholar] [CrossRef]
  19. Nagasaki, H.; Itoh, J.-I.; Hayashi, K.; Hibara, K.-I.; Satoh-Nagasawa, N.; Nosaka, M.; Mukouhata, M.; Ashikari, M.; Kitano, H.; Matsuoka, M.; et al. The small interfering RNA production pathway is required for shoot meristem initiation in rice. Proc. Natl. Acad. Sci. USA 2007, 104, 14867–14871. [Google Scholar] [CrossRef]
  20. Chen, X.J.; Qi, C.K.; Pu, H.M.; Zhang, J.F.; Gao, J.Q.; Fu, S.Z. Evaluation of lodging resistance in rapeseed (Brassica napus L.) and relationship between plant architecture and lodging resistance. Chin. J. Oil Crop Sci. 2007, 29, 54–57. [Google Scholar]
  21. Li, N.; Shi, J.; Wang, X.; Liu, G.; Wang, H. A combined linkage and regional association mapping validation and fine mapping of two major pleiotropic QTLs for seed weight and silique length in rapeseed (Brassica napus L.). BMC Plant Biol. 2014, 14, 114. [Google Scholar] [CrossRef]
  22. Khan, F.; Ali, S.; Shakeel, A.; Saeed, A.; Abbas, G. Correlation analysis of some quantitative characters in Brassica napus L. J. Agric. Res. 2006, 44, 7–14. [Google Scholar]
  23. Wang, X.; Chen, L.; Wang, A.; Wang, H.; Tian, J.; Zhao, X.; Chao, H.; Zhao, Y.; Zhao, W.; Xiang, J.; et al. Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus. BMC Plant Biol. 2016, 16, 71. [Google Scholar] [CrossRef] [Green Version]
  24. Sadat, H.A.; Nematzadeh, G.A.; Jelodar, N.B.; Chapi, O.G. Genetic evaluation of yield and yield components at advanced generations in rapeseed (Brassica napus L.). Afr. J. Agric. Res. 2010, 5, 1958–1964. [Google Scholar]
  25. Samizadeh, H.; Samadi, B.Y.; Behamta, M.; Taleii, A.; Stringam, G. Study of pod length trait in doubled haploid Brassica napus population by molecular markers. J. Agric. Sci. 2007, 9, 129–136. [Google Scholar]
  26. Zhang, L.; Yang, G.; Liu, P.; Hong, D.; Li, S.; He, Q. Genetic and correlation analysis of silique-traits in Brassica napus L. by quantitative trait locus mapping. Theor. Appl. Genet. 2011, 122, 21–31. [Google Scholar] [CrossRef]
  27. Venglat, S.P.; Dumonceaux, T.; Rozwadowski, K.; Parnell, L.; Babic, V.; Keller, W.; Martienssen, R.; Selvaraj, G.; Datla, R. The homeobox gene BREVIPEDICELLUS is a key regulator of inflorescence architecture in Arabidopsis. Proc. Natl. Acad. Sci. USA 2002, 99, 4730–4735. [Google Scholar] [CrossRef] [PubMed]
  28. Lincoln, C.; Long, J.; Yamaguchi, J.; Serikawa, K.; Hake, S. A knotted 1-like homeobox gene in Arabidopsis is expressed in the vegetative meristem and dramatically alters leaf morphology when overexpressed in transgenic plants. Plant Cell 1994, 6, 1859–1876. [Google Scholar]
  29. Nari, J.; Noat, G.; Ricard, J. Pectin methylesterase, metal ions and plant cell-wall extension. Hydrolysis of pectin by plant cell-wall pectin methylesterase. Biochem. J. 1991, 279, 343–350. [Google Scholar] [CrossRef] [PubMed]
  30. Pelloux, J.; Rustérucci, C.; Mellerowicz, E.J. New insights into pectin methylesterase structure and function. Trends Plant Sci. 2007, 12, 267–277. [Google Scholar] [CrossRef] [PubMed]
  31. Röckel, N.; Wolf, S.; Kost, B.; Rausch, T.; Greiner, S. Elaborate spatial patterning of cell-wall PME and PMEI at the pollen tube tip involves PMEI endocytosis, and reflects the distribution of esterified and de-esterified pectins. Plant J. 2008, 53, 133–143. [Google Scholar] [CrossRef] [PubMed]
  32. Chebli, Y.; Geitmann, A. Cellular growth in plants requires regulation of cell wall biochemistry. Curr. Opin. Cell Biol. 2017, 44, 28–35. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, F.; Chen, W.; Traas, J. Mechanical signaling in plant morphogenesis. Curr. Opin. Genet. Dev. 2018, 51, 26–30. [Google Scholar] [CrossRef] [PubMed]
  34. Wolf, S.; Mravec, J.; Greiner, S.; Mouille, G.; Höfte, H. Plant cell wall homeostasis is mediated by brassinosteroid feedback signaling. Curr. Biol. 2012, 22, 1732–1737. [Google Scholar] [CrossRef]
  35. Müller, K.; Levesque-Tremblay, G.; Fernandes, A.; Wormit, A.; Bartels, S.; Usadel, B.; Kermode, A. Overexpression of a pectin methylesterase inhibitor in Arabidopsis thaliana leads to altered growth morphology of the stem and defective organ separation. Plant Signal. Behav. 2013, 8, e26464. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, Y.M.; Chen, W.J.; Pu, C.; Wan, S.B.; Mao, Y.; Wang, M.; Guan, R.Z. Mapping amajor QTL responsible for dwarf architecture in Brassica napus using a single-nucleotide polymorphism marker approach. BMC Plant Biol. 2016, 16, 178. [Google Scholar] [CrossRef] [Green Version]
  37. Li, H.; Li, J.; Song, J.; Zhao, B.; Guo, C.; Wang, B.; Zhang, Q.; Wang, J.; King, G.J.; Liu, K. An auxin signaling gene BnaA3.IAA7 contributes to improved plant architecture and yield heterosis in rapeseed. New Phytol. 2019, 222, 837–851. [Google Scholar] [CrossRef]
  38. Zhao, B.; Wang, B.; Li, Z.; Guo, T.; Zhao, J.; Guan, Z.; Liu, K. Identification and characterization of a new dwarf locus DS-4 encoding an aux/IAA7 protein in Brassica napus. Theor. Appl. Genet. 2019, 132, 1435–1449. [Google Scholar] [CrossRef]
  39. Yang, M.; Huang, C.W.; Wang, M.M.; Fan, H.; Wan, S.B.; Wang, Y.; He, J.; Guan, R.Z. Fine mapping of an up-curling leaf locus (BnUC1) in Brassica napus. BMC Plant Biol. 2019, 19, 324. [Google Scholar] [CrossRef]
  40. Huang, C.; Yang, M.; Shao, D.; Wang, Y.; Wan, S.; He, J.; Meng, Z.; Guan, R. Fine mapping of the BnUC2 locus related to leaf up-curling and plant semi-dwarfing in Brassica napus. BMC Genom. 2020, 21, 530. [Google Scholar] [CrossRef]
  41. Wan, S.; Qin, Z.; Jiang, X.; Yang, M.; Chen, W.; Wang, Y.; Ni, F.; Guan, Y.; Guan, R.Z. Identification and Fine Mapping of a Locus Related to Leaf Up-Curling Trait (Bnuc3) in Brassica napus. Int. J. Mol. Sci. 2021, 22, 11693. [Google Scholar] [CrossRef]
  42. Mai, Y.X.; Wang, L.; Yang, H.Q. A gain-of-function mutation in IAA7/AXR2 confers late flowering under short-day light in Arabidopsis. J. Integr. Plant Biol. 2011, 53, 480–592. [Google Scholar] [CrossRef] [PubMed]
  43. Liu, S.; Hu, Q.; Sha, L.; Li, Q.; Yang, X.; Wang, X.; Wang, S. Expression of wild-type PtrIAA14.1, a poplar aux/IAA gene causes morphological changes in Arabidopsis. Front. Plant Sci. 2015, 6, 388. [Google Scholar] [CrossRef]
  44. Hou, Y.; Li, H.; Zhai, L.; Xie, X.; Li, X.; Bian, S. Identification and functional characterization of the Aux/IAA gene VcIAA27 in blueberry. Plant Signal. Behav. 2020, 15, 1700327. [Google Scholar] [CrossRef] [PubMed]
  45. Wen, F.; Zhu, Y.; Hawes, M.C. Effect of pectin methylesterase gene expression on pea root development. Plant Cell 1999, 11, 1129–1140. [Google Scholar] [CrossRef]
  46. Hasunuma, T.; Fukusaki, E.I.; Kobayashi, A. Expression of fungal pectin methylesterase in transgenic tobacco leads to alteration in cell wall metabolism and a dwarf phenotype. J. Biotechnol. 2004, 111, 241–251. [Google Scholar] [CrossRef]
  47. Prasanna, V.; Prabha, T.N.; Tharanathan, R.N. Fruit ripening phenomena—An overview. Crit. Rev. Food Sci. Nutr. 2007, 47, 1–19. [Google Scholar] [CrossRef]
  48. Yang, R.; Wang, T.; Shi, W.; Li, S.; Liu, Z.; Wang, J.; Yang, Y. E3 ubiquitin ligase ATL61 acts as a positive regulator in abscisic acid mediated drought response in Arabidopsis. Biochem. Biophys. Res. Commun. 2020, 528, 292–298. [Google Scholar] [CrossRef] [PubMed]
  49. Cui, L.H.; Min, H.J.; Yu, S.G.; Byun, M.Y.; Oh, T.R.; Lee, A.; Yang, H.W.; Kim, W.T. OsATL38 mediates mono-ubiquitination of the 14-3-3 protein OsGF14d and negatively regulates the cold stress response in rice. J. Exp. Bot. 2022, 73, 307–323. [Google Scholar] [CrossRef] [PubMed]
  50. Kong, F.; Ramonell, K. Arabidopsis Toxicos en Levadura 12 Modulates Salt Stress and ABA Responses in Arabidopsis thaliana. Int. J. Mol. Sci. 2022, 23, 7290. [Google Scholar] [CrossRef]
  51. Aoyama, S.; Terada, S.; Sanagi, M.; Hasegawa, Y.; Lu, Y.; Morita, Y.; Chiba, Y.; Sato, T.; Yamaguchi, J. Membrane-localized ubiquitin ligase ATL15 functions in sugar-responsive growth regulation in Arabidopsis. Biochem. Biophys. Res. Commun. 2017, 491, 33–39. [Google Scholar] [CrossRef]
  52. Ramaiah, M.; Jain, A.; Yugandhar, P.; Raghothama, K.G. ATL8, a RING E3 ligase, modulates root growth and phosphate homeostasis in Arabidopsis. Plant Physiol. Biochem. 2022, 179, 90–99. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, S.M.; Lo, M.C.; Shang, C.; Kao, S.C.; Tseng, Y.Z. Role of M-line proteins in sarcomeric titin assembly during cardiac myofibrillogenesis. J. Cell. Biochem. 1998, 71, 82–95. [Google Scholar] [CrossRef]
  54. Qadota, H.; Matsunaga, Y.; Nguyen, K.Q.; Mattheyses, A.; Hall, D.H.; Benian, G.M. High-resolution imaging of muscle attachment structures in Caenorhabditis elegans. Cytoskeleton 2017, 74, 426–442. [Google Scholar] [CrossRef]
  55. Iwakawa, H.; Takahashi, H.; Machida, Y.; Machida, C. Roles of ASYMMETRIC LEAVES2 (AS2) and Nucleolar Proteins in the Adaxial–Abaxial Polarity Specification at the Perinucleolar Region in Arabidopsis. Int. J. Mol. Sci. 2020, 21, 7314. [Google Scholar] [CrossRef] [PubMed]
  56. Waites, R.; Hudson, A. Phantastica: A gene required for dorsoventrality of leaves in Antirrhinum majus. Development 1995, 121, 2143–2154. [Google Scholar] [CrossRef]
  57. Waites, R.; Selvadurai, H.R.; Oliver, I.R.; Hudson, A. The PHANTASTICA Gene Encodes a MYB Transcription Factor Involved in Growth and Dorsoventrality of Lateral Organs in Antirrhinum. Cell 1998, 93, 779–789. [Google Scholar] [CrossRef]
  58. McConnell, J.R.; Emery, J.; Eshed, Y.; Bao, N.; Bowman, J.; Barton, M.K. Role of PHABULOSA and PHAVOLUTA in deter-mining radial patterning in shoots. Nature 2001, 411, 709–713. [Google Scholar] [CrossRef]
  59. Itoh, J.I.; Hibara, K.-I.; Sato, Y.; Nagato, Y. Developmental Role and Auxin Responsiveness of Class III Homeodomain Leucine Zipper Gene Family Members in Rice. Plant Physiol. 2008, 147, 1960–1975. [Google Scholar] [CrossRef] [PubMed]
  60. Manuela, D.; Xu, M. Patterning a Leaf by Establishing Polarities. Front. Plant Sci. 2020, 11, 568730. [Google Scholar] [CrossRef]
  61. Wu, G.; Lin, W.-C.; Huang, T.; Poethig, R.S.; Springer, P.S.; Kerstetter, R.A. KANADI1 regulates adaxial-abaxial polarity in Arabidopsis by directly repressing the transcription of ASYMMETRIC LEAVES2. Proc. Natl. Acad. Sci. USA 2008, 105, 16392–16397. [Google Scholar] [CrossRef]
  62. Emery, J.F.; Floyd, S.K.; Alvarez, J.; Eshed, Y.; Hawker, N.P.; Izhaki, A.; Baum, S.F.; Bowman, J.L. Radial patterning of Ara-bidopsis shoots by class IIIHD-ZIP and KANADI genes. Curr. Biol. 2003, 13, 1768–1774. [Google Scholar] [CrossRef]
  63. Guan, C.; Wu, B.; Yu, T.; Wang, Q.; Krogan, N.T.; Liu, X.; Jiao, Y. Spatial Auxin Signaling Controls Leaf Flattening in Arabidopsis. Curr. Biol. 2017, 27, 2940–2950. [Google Scholar] [CrossRef]
  64. Kelley, D.; Arreola, A.; Gallagher, T.L.; Gasser, C.S. ETTIN (ARF3) physically interacts with KANADI proteins to form a functional complex essential for integument development and polarity determination in Arabidopsis. Development 2012, 139, 1105–1109. [Google Scholar] [CrossRef]
  65. Du, F.; Gong, W.; Boscá, S.; Tucker, M.; Laux, T. Dose-dependent AGO1-mediated inhibition of the miRNA165/166 pathway modulates stem cell maintenance in arabidopsis shoot apical meristem. Plant Commun. 2019, 1, 100002. [Google Scholar] [CrossRef] [PubMed]
  66. Kirolinko, C.; Hobecker, K.; Wen, J.; Mysore, K.S.; Niebel, A.; Blanco, F.A.; Zanetti, M.E. Auxin Response Factor 2 (ARF2), ARF3,and ARF4 Mediate Both Lateral Root and Nitrogen Fixing Nodule Development in Medicago truncatula. Front. Plant Sci. 2021, 12, 659061. [Google Scholar] [CrossRef]
  67. Yifhar, T.; Pekker, I.; Peled, D.; Friedlander, G.; Pistunov, A.; Sabban, M.; Wachsman, G.; Alvarez, J.P.; Amsellem, Z.; Eshed, Y. Failure of the tomato trans-acting short interfering RNA program to regulate AUXIN RESPONSE FACTOR3 and ARF4 un-derlies the wiry leaf syndrome. Plant Cell 2012, 24, 3575–3589. [Google Scholar] [CrossRef] [PubMed]
  68. Roeder, A.H.; Ferrándiz, C.; Yanofsky, M.F. The role of the REPLUMLESS homeodomain protein in patterning the Arabidopsis fruit. Curr. Biol. 2003, 13, 1630–1635. [Google Scholar] [CrossRef]
  69. Liljegren, S.J.; Ditta, G.S.; Eshed, Y.; Savidge, B.; Bowman, J.L.; Yanofsky, M.F. SHATTERPROOF MADS-box genes control seed dispersal in Arabidopsis. Nature 2000, 404, 766–770. [Google Scholar] [CrossRef]
  70. Liljegren, S.J.; Roeder, A.H.; Kempin, S.A.; Gremski, K.; Østergaard, L.; Guimil, S.; Reyes, D.K.; Yanofsky, M.F. Control of fruit patterning in Arabidopsis by INDEHISCENT. Cell 2004, 116, 843–853. [Google Scholar] [CrossRef]
  71. Rajani, S.; Sundaresan, V. The Arabidopsis myc/bHLH gene ALCATRAZ enables cell separation in fruit dehiscence. Curr. Biol. 2001, 11, 1914–1922. [Google Scholar] [CrossRef] [PubMed]
  72. Qi, J.; Wu, B.; Feng, S.; Lü, S.; Guan, C.; Zhang, X.; Qiu, D.; Hu, Y.; Zhou, Y.; Li, C.; et al. Mechanical regulation of organ asymmetry in leaves. Nat. Plants. 2017, 3, 724–3733. [Google Scholar] [CrossRef] [PubMed]
  73. Duncan, W.G. Leaf Angles, Leaf Area, and Canopy Photosynthesis 1. Crop Sci. 1971, 11, 482–485. [Google Scholar] [CrossRef]
  74. Liu, X.; Li, M.; Liu, K.; Tang, D.; Sun, M.; Li, Y.; Shen, Y.; Du, G.; Cheng, Z. Semi-Rolled Leaf2modulates rice leaf rolling by regulating abaxial side cell differentiation. J. Exp. Bot. 2016, 67, 2139–2150. [Google Scholar] [CrossRef] [PubMed]
  75. Mahmood, S.; Li, Z.; Yue, X.; Wang, B.; Chen, J.; Liu, K.D. Development of INDELs markers in oilseed rape (Brassica napus L.) using re-sequencing data. Mol. Breed. 2016, 36, 79. [Google Scholar] [CrossRef]
  76. Takagi, H.; Abe, A.; Yoshida, K.; Kosugi, S.; Natsume, S.; Mitsuoka, C.; Uemura, A.; Utsushi, H.; Tamiru, M.; Takuno, S.; et al. QTL-seq: Rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 2013, 74, 174–183. [Google Scholar] [CrossRef]
  77. Wang, B.; Wu, Z.K.; Li, Z.; Zhang, Q.H.; Hu, J.L.; Xiao, Y.J.; Cai, D.F.; Wu, J.S.; King, G.; Li, H.T.; et al. Dissection of the genetic architecture of three seed-quality traits and consequences for breeding in Brassica napus. Plant Biotechnol. J. 2018, 16, 1336–1348. [Google Scholar] [CrossRef]
  78. Ooijen, J.V. JoinMap 4, Software for the Calculation of Genetic Linkage Maps in Experimental Populations; Kyazma B.V.: Wageningen, The Netherlands, 2006. [Google Scholar]
  79. Singh, V.K.; Mangalam, A.K.; Dwivedi, S.; Naik, S. Primer premier: Program for design of degenerate primers from a protein sequence. Biotechniques 1998, 24, 318–319. [Google Scholar] [CrossRef]
  80. Wang, Y.; Wan, S.; Fan, H.; Yang, M.; Li, W.; Guan, R. A sulfotransferase gene BnSOT-like1 has a minor genetic effect on seed glucosinolate content in Brassica napus. Crop J. 2020, 8, 855–865. [Google Scholar] [CrossRef]
  81. Thompson, J.D.; Gibson, T.J.; Frédéric, P.; Franois, J.; Higgins, D.G. The CLUSTAL_X Windows Interface. Flexible Strategies for Multiple Sequence Alignment Aided by Quality Analysis Tools. Nucleic Acids Res. 1997, 24, 4876–4882. [Google Scholar] [CrossRef]
  82. Huggett, J.F.; Foy, C.A.; Vladimir, B.; Kerry, E.; Garson, J.A.; Ross, H.; Jan, H.; Mikael, K.; Mueller, R.D.; Tania, N. The digital MIQE guidelines: Minimum information for publication of quantitative digital PCR experiments. Clin. Chem. 2013, 59, 892–902. [Google Scholar] [CrossRef] [PubMed]
  83. Arnon, D.I. Copper enzymes in isolated chloroplasts: Polyphenoloxidase in Beta vulgaris. Plant Physiol. 1949, 24, 1–15. [Google Scholar] [CrossRef] [PubMed]
  84. Zhang, R.; Liu, H.H.; Zhao, H.X.; Hu, S. Comparison of two protein extraction methods for proteomic analysis of chlorophyll-deficient mutants in Brassica juncea L. Prog. Biochem. Biophys. 2010, 37, 1025–1032. [Google Scholar] [CrossRef]
Ijms 24 03069 g001 550
Figure 1. Performance of the parents NJAU5773 and ZS11 for constructing mapping population. (a) shows the leaves of the parent NJAU5773 with up-curling leaves (right) and the parent ZS11 with flat leaves (left) at the seeding stage. (b) shows the flowers of the parent NJAU5773 with downward-pointing flowers (right) and the parent ZS11 with normal flowers (left) at the flower stage. (c) shows the siliques of the parent NJAU5773 with downward-pointing silique (right) and the parent ZS11 with upright silique (left) at the maturity stage.
Figure 1. Performance of the parents NJAU5773 and ZS11 for constructing mapping population. (a) shows the leaves of the parent NJAU5773 with up-curling leaves (right) and the parent ZS11 with flat leaves (left) at the seeding stage. (b) shows the flowers of the parent NJAU5773 with downward-pointing flowers (right) and the parent ZS11 with normal flowers (left) at the flower stage. (c) shows the siliques of the parent NJAU5773 with downward-pointing silique (right) and the parent ZS11 with upright silique (left) at the maturity stage.
Ijms 24 03069 g001
Ijms 24 03069 g002 550
Figure 2. Mapping the leaf up-curling and silique downward-pointing locus BnUD1. (a) indicates that the BnUD1 locus was preliminary mapped in an interval of 1.63 Mb (shaded in green) between InDel markers InDel-13 and Indel-15 (black lines) on A05 chromosome. (b) shows that the BnUD1 locus was blue-colored in the physical interval of 466.58 kb (shaded in green) with facilitation of the developed InDel markers. (c) shows that the BnUD1 locus was fine mapped in the physical interval of 54.84 kb (shaded in green). (d) shows that BnUD1 was in a 0.135 cM region between InDel markers Indel-78 and Indel-15 (red lines), shaded in green. The BnUD1 locus is highlighted in red.
Figure 2. Mapping the leaf up-curling and silique downward-pointing locus BnUD1. (a) indicates that the BnUD1 locus was preliminary mapped in an interval of 1.63 Mb (shaded in green) between InDel markers InDel-13 and Indel-15 (black lines) on A05 chromosome. (b) shows that the BnUD1 locus was blue-colored in the physical interval of 466.58 kb (shaded in green) with facilitation of the developed InDel markers. (c) shows that the BnUD1 locus was fine mapped in the physical interval of 54.84 kb (shaded in green). (d) shows that BnUD1 was in a 0.135 cM region between InDel markers Indel-78 and Indel-15 (red lines), shaded in green. The BnUD1 locus is highlighted in red.
Ijms 24 03069 g002
Ijms 24 03069 g003 550
Figure 3. Experimental results of polymorphic markers Indel-78 and Indel-15. Marker scan with progeny BC5F3 and BC6F2 populations derived from the parents ZS11 and NJAU5773 was conducted. P1 and P2 indicates PCR products from the parents ZS11 and NJAU5773 plants, respectively. (a) Shows partial results for InDel marker Indel-78 in which 1–3 and 6 denote the PCR products from the heterozygous plants in the progeny populations, and 4 and 5 denotes products from the homozygous plants with flatten leaves and upright siliques and plants with up-curling leaves and downward-pointing siliques, respectively. (b) Shows partial results for InDel marker Indel-15 in which 7 and 9–11 denote the PCR products from the heterozygous plants in the progeny populations, and 8 and 12 denotes products from the homozygous plants with flatten leaves and upright siliques and plants with up-curling leaves and downward-pointing siliques, respectively.
Figure 3. Experimental results of polymorphic markers Indel-78 and Indel-15. Marker scan with progeny BC5F3 and BC6F2 populations derived from the parents ZS11 and NJAU5773 was conducted. P1 and P2 indicates PCR products from the parents ZS11 and NJAU5773 plants, respectively. (a) Shows partial results for InDel marker Indel-78 in which 1–3 and 6 denote the PCR products from the heterozygous plants in the progeny populations, and 4 and 5 denotes products from the homozygous plants with flatten leaves and upright siliques and plants with up-curling leaves and downward-pointing siliques, respectively. (b) Shows partial results for InDel marker Indel-15 in which 7 and 9–11 denote the PCR products from the heterozygous plants in the progeny populations, and 8 and 12 denotes products from the homozygous plants with flatten leaves and upright siliques and plants with up-curling leaves and downward-pointing siliques, respectively.
Ijms 24 03069 g003
Ijms 24 03069 g004 550
Figure 4. Partial molecular marker experimental results of the InDel marker BnA05ID1. ‘M’ denotes 5000 marker. The number 1, 3–5, 7–9, 11, 16–17, and 19 denote the PCR products from homozygous plants with up-curling leaf and downward-pointing silique. The number 6, 10, 12, 14, 18, 22, and 23 denote the PCR products from heterozygous plants with up-curling leaf and downward-pointing silique. The number 2, 13, 15, 20, and 21 denote the PCR products from homozygous plants with flatten leaf and upright silique.
Figure 4. Partial molecular marker experimental results of the InDel marker BnA05ID1. ‘M’ denotes 5000 marker. The number 1, 3–5, 7–9, 11, 16–17, and 19 denote the PCR products from homozygous plants with up-curling leaf and downward-pointing silique. The number 6, 10, 12, 14, 18, 22, and 23 denote the PCR products from heterozygous plants with up-curling leaf and downward-pointing silique. The number 2, 13, 15, 20, and 21 denote the PCR products from homozygous plants with flatten leaf and upright silique.
Ijms 24 03069 g004
Ijms 24 03069 g005 550
Figure 5. Differential expression of 3 mutated genes in the mapping interval in leaves of the advanced backcross population. The Actin of B. napus was used as the reference gene for normalization in qRT-PCR experiments. The relative mRNA level of genes in plants with Bnud1 mutated traits (grey) was calculated in reference to gene expression level of plants with wild-type traits (black). Values shown are means ± SD (n = 4). * Denotes significant at probability level of 0.05.
Figure 5. Differential expression of 3 mutated genes in the mapping interval in leaves of the advanced backcross population. The Actin of B. napus was used as the reference gene for normalization in qRT-PCR experiments. The relative mRNA level of genes in plants with Bnud1 mutated traits (grey) was calculated in reference to gene expression level of plants with wild-type traits (black). Values shown are means ± SD (n = 4). * Denotes significant at probability level of 0.05.
Ijms 24 03069 g005
Table
Table 1. Inheritance of the leaf up-curling and silique downward-pointing traits in populations derived from the two parents in B. napus.
Table 1. Inheritance of the leaf up-curling and silique downward-pointing traits in populations derived from the two parents in B. napus.
PopulationUDWTTotalExpectation χ 2
F130030
RF130030
F2100421423:11.43
BC158701281:11.13
BC259731321:11.48
BC355651201:10.83
BC463531161:10.86
BC574601341:11.46
BC653651181:11.22
BC5F35272077343:13.56
BC6F2220883083:11.91
UD denotes plants with up-curling leaf and downward-pointing silique plants; WT denotes that with flatten leaves and upright siliques.
Table
Table 2. Genes in the mapped interval and the Brassica within the interval on the A05 chromosome.
Table 2. Genes in the mapped interval and the Brassica within the interval on the A05 chromosome.
Gene in B. napus. cv ZS11Gene in B. napus. cv DarmorHomologue in A. thalianaGene Function
BnaA05G0157300ZSBnaA05g14090DAT2G20619Plant thionin family protein
BnaA05G0157400ZS Unknown
BnaA05G0157500ZS AT4G18570Tetratricopeptide repeat (TPR)-like superfamily protein
BnaA05G0157600ZS AT4G18570Tetratricopeptide repeat (TPR)-like superfamily protein
BnaA05G0157700ZSBnaA05g14100DAT1G53860Remorin family protein
BnaA05G0157800ZSBnaA05g14110DAT1G5385020S proteasome alpha subunit E1
BnaA05G0157900ZSBnaA05g14120DAT1G53840Pectin methylesterase 1
BnaA05G0158000ZSBnaA05g14130DAT1G53830Pectin methylesterase 2
BnaA05G0158100ZSBnaA05g14140DAT1G53820RING/U-box superfamily protein
BnaA05G0158200ZSBnaA05g14150DAT5G38830Cysteinyl-tRNA synthetase, class Ia family protein
BnaA05G0158300ZSBnaA05g14160DAT1G53800Muscle M-line assembly protein
Table
Table 3. Agronomic traits comparison between plants with BnUD1 locus and that without BnUD1 locus in the BC6F2 population.
Table 3. Agronomic traits comparison between plants with BnUD1 locus and that without BnUD1 locus in the BC6F2 population.
TraitPlants without BnUD1 LocusPlants with BnUD1 Locus
Plant height (cm)180.76 ± 1.95144.61 ± 6.49 *
Branch height (cm)56.86 ± 5.0333.70 ± 4.55 *
Main inflorescence length (cm)73.94 ± 2.4466.50 ± 8.41 *
Stem diameter (mm)24.37 ± 1.8523.10 ± 0.66
Number of first effective branch 8.80 ± 0.846.88 ± 0.99 *
Siliques of main inflorescence72.63 ± 6.7871.60 ± 4.45
Total siliques per plant346.20 ± 18.83363.13 ± 29.86
Silique length9.23 ± 0.639.37 ± 0.78
Seeds per siliques26.94 ± 1.9430.67 ± 2.54 *
1000-seed weight (g)5.42 ± 0.184.13 ± 0.15 *
* Indicates significant at the 0.05 probability level by t-test. Data are shown as mean ± SD (n = 30 for each sample).
Table
Table 4. Leaf chlorophyll contents in leaves of plants with BnUD1 locus and that without BnUD1 locus.
Table 4. Leaf chlorophyll contents in leaves of plants with BnUD1 locus and that without BnUD1 locus.
GenotypeChl a (mg/g)Chl b (mg/g)Total Chl a/b Ratio
Plants without BnUD1 locus1.59 ± 0.330.96 ± 0.322.55 ± 0.621.75 ± 0.37
Plants with BnUD1 locus3.39 ± 0.36 *1.80 ± 0.23 *5.19 ± 0.57 *1.89 ± 0.27 *
* Indicates significant at the 0.05 probability level. Mean ± standard deviation (SD) under sample size. (n = 15 for each sample).
Table
Table 5. Photosynthetic indicators in leaves of plants with BnUD1 locus and that without BnUD1 locus.
Table 5. Photosynthetic indicators in leaves of plants with BnUD1 locus and that without BnUD1 locus.
GenotypeNPR µmol CO2 m−2 s−1SC mol H2O m−2 s−1ICC µmol CO2 mol−1TR mmol H2O m−2 s−1
Plants without BnUD1 locus8.96 ± 0.480.24 ± 0.03368.86 ± 4.822.46 ± 0.31
Plants with BnUD1 locus11.89 ± 0.76 *0.42 ± 0.02 *425.50 ± 6.37 *2.89 ± 0.37
Data are presented as means ± SD. * Indicates significant at 0.05 probability level. NPR, SC, ICC, and TR denotes net photosynthetic rate, stomatal conductance, intercellular CO2 concentration and transpiration rate, respectively. (n = 6 for each sample).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, M.; Chen, J.; Chang, Y.; Wan, S.; Zhao, Z.; Ni, F.; Guan, R. Fine Mapping of a Pleiotropic Locus (BnUD1) Responsible for the Up-Curling Leaves and Downward-Pointing Siliques in Brassica napus. Int. J. Mol. Sci. 2023, 24, 3069. https://doi.org/10.3390/ijms24043069

AMA Style

Yang M, Chen J, Chang Y, Wan S, Zhao Z, Ni F, Guan R. Fine Mapping of a Pleiotropic Locus (BnUD1) Responsible for the Up-Curling Leaves and Downward-Pointing Siliques in Brassica napus. International Journal of Molecular Sciences. 2023; 24(4):3069. https://doi.org/10.3390/ijms24043069

Chicago/Turabian Style

Yang, Mao, Jun Chen, Yuqing Chang, Shubei Wan, Zisu Zhao, Fei Ni, and Rongzhan Guan. 2023. "Fine Mapping of a Pleiotropic Locus (BnUD1) Responsible for the Up-Curling Leaves and Downward-Pointing Siliques in Brassica napus" International Journal of Molecular Sciences 24, no. 4: 3069. https://doi.org/10.3390/ijms24043069

APA Style

Yang, M., Chen, J., Chang, Y., Wan, S., Zhao, Z., Ni, F., & Guan, R. (2023). Fine Mapping of a Pleiotropic Locus (BnUD1) Responsible for the Up-Curling Leaves and Downward-Pointing Siliques in Brassica napus. International Journal of Molecular Sciences, 24(4), 3069. https://doi.org/10.3390/ijms24043069

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop