Genome-scale targeted mutagenesis in Brassica napus using a pooled CRISPR library.

The recently constructed mutant libraries of diploid crops by the CRISPR-Cas9 system have provided abundant resources for functional genomics and crop breeding. However, because of the genome complexity, it is a big challenge to accomplish large-scale targeted mutagenesis in polyploid plants. Here, we demonstrate the feasibility of using a pooled CRISPR library to achieve genome-scale targeted editing in an allotetraploid crop of Brassica napus A total of 18,414 sgRNAs were designed to target 10,480 genes of interest, and afterward, 1104 regenerated transgenic plants harboring 1088 sgRNAs were obtained. Editing interrogation results revealed that 93 of the 178 genes were identified as mutated, thus representing an editing efficiency of 52.2%. Furthermore, we have discovered that Cas9-mediated DNA cleavages tend to occur at all the target sites guided by the same individual sgRNA, a novel finding in polyploid plants. Finally, we show the strong capability of reverse genetic screening for various traits with the postgenotyped plants. Several genes, which might dominate the fatty acid profile and seed oil content and have yet to be reported, were unveiled from the forward genetic studies. Our research provides valuable resources for functional genomics, elite crop breeding, and a good reference for high-throughput targeted mutagenesis in other polyploid plants.

Genetic dissection of Brassica napus seed vigor after aging.

KEY MESSAGE: Stable and novel QTLs that affect seed vigor under different storage durations were discovered, and BnaOLE4, located in the interval of cqSW-C2-3, increased seed vigor after aging. Seed vigor is an important trait in crop breeding; however, the underlying molecular regulatory mechanisms governing this trait in rapeseed remain largely unknown. In the present study, vigor-related traits were analyzed in seeds from a doubled haploid (DH) rapeseed (Brassica napus) population grown in 2 different environments using seeds stored for 7, 5, and 3 years under natural storage conditions. A total of 229 quantitative trait loci (QTLs) were identified and were found to explain 3.78%-17.22% of the phenotypic variance for seed vigor-related traits after aging. We further demonstrated that seed vigor-related traits were positively correlated with oil content (OC) but negatively correlated with unsaturated fatty acids (FAs). Some pleiotropic QTLs that collectively regulate OC, FAs, and seed vigor, such as uq.A8, uq.A3-2, uq.A9-2, and uq.C3-1, were identified. The transcriptomic results from extreme pools of DH lines with distinct seed vigor phenotypes during accelerated aging revealed that various biological pathways and metabolic processes (such as glutathione metabolism and reactive oxygen species) were involved in seed vigor. Through integration of QTL analysis and RNA-Seq, a regulatory network for the control of seed vigor was constructed. Importantly, a candidate (BnaOLE4) from cqSW-C2-3 was selected for functional analysis, and transgenic lines overexpressing BnaOLE4 showed increased seed vigor after artificial aging. Collectively, these results provide novel information on QTL and potential candidate genes for molecular breeding for improved seed storability.

Class A lysophosphatidic acid acyltransferase 2 from Camelina sativa promotes very long-chain fatty acids accumulation in phospholipid and triacylglycerol.

Very long-chain fatty acids (VLCFAs) are important industrial raw materials and can be produced by genetically modified oil plants. For a long time, class A lysophosphatidic acid acyltransferase (LPAT) was considered unable to promote the accumulation of VLCFA in oil crops. The bottlenecks that the transgenic high VLCFA lines have an oil content penalty and the low amount of VLCFA in phosphatidylcholine remains intractable. In the present study, a class A LPAT2 from Camelina sativa (CsaLPAT2) promoting VLCFAs accumulation in phospholipid was found. Overexpression of CsaLPAT2 alone in Arabidopsis seeds significantly increased the VLCFA content in triacylglycerol, including C20:0, C20:2, C20:3, C22:0, and C22:1. The proportion of phosphatidic acid molecules containing VLCFAs in transgenic seeds reached up to 45%, which was 2.8-fold greater than that in wild type. The proportion of phosphatidylcholine and diacylglycerol molecules containing VLCFAs also increased significantly. Seed size in CsaLPAT2 transgenic lines showed a slight increase without an oil content penalty. The total phospholipid content in the seed of the CsaLPAT2 transgenic line was significantly increased. Furthermore, the function of class A LPAT in promoting the accumulation of VLCFAs is conserved in the representative oil crops of Brassicaceae, such as C. sativa, Arabidopsis thaliana, Brassica napus, Brassica rapa, and Brassica oleracea. The findings of this study provide a promising gene resource for the production of VLCFAs.

Characterization of Oil Body and Starch Granule Dynamics in Developing Seeds of Brassica napus.

Brassica napus is the most important oilseed crop in the world, and the lipid was stored in the oil body (OB) in the form of triacylglycerol. At present, most of studies on the relationship between oil body morphology and seed oil content in B. napus was focused on mature seeds. In the present study, the OBs in different developing seeds of B. napus with relatively high oil content (HOC) of about 50% and low oil content (LOC) of about 39% were analyzed. It was revealed that the size of OBs was first increased and then decreased in both materials. And in late seed developmental stages, the average OB size of rapeseed with HOC was higher than that of LOC, while it was reversed in the early seed developmental stages. No significant difference was observed on starch granule (SG) size in HOC and LOC rapeseed. Further results indicated that the expression of genes that involved in malonyl-CoA metabolism, fatty acid carbon chain extension, lipid metabolism, and starch synthesis in the rapeseed with HOC was higher than that of rapeseed with LOC. These results give some new insight for understanding the dynamics of OBs and SGs in embryos of B. napus.

Ultra-high α-linolenic acid accumulating developmental defective embryo was rescued by lysophosphatidic acid acyltransferase 2.

For decades, genetic engineering approaches to produce unusual fatty acids (UFAs) in crops has reached a bottleneck, including reduced seed oil production and seed vigor. Currently, plant models in the field of research are primarily used to investigate defects in oil production and seedling development, while the role of UFAs in embryonic developmental defects remains unknown. In this study, we developed a transgenic Arabidopsis plant model, in which the embryo exhibits severely wrinkled appearance owing to α-linolenic acid (ALA) accumulation. RNA-sequencing analysis in the defective embryo suggested that brassinosteroid synthesis, FA synthesis and photosynthesis were inhibited, while FA degradation, endoplasmic reticulum stress and oxidative stress were activated. Lipidomics analysis showed that ultra-accumulated ALA is released from phosphatidylcholine as a free FA in cells, inducing severe endoplasmic reticulum and oxidative stress. Furthermore, we identified that overexpression of lysophosphatidic acid acyltransferase 2 rescued the defective phenotype. In the rescue line, the pool capacity of the Kennedy pathway was increased, and the esterification of ALA indirectly to triacylglycerol was enhanced to avoid stress. This study provides a plant model that aids in understanding the molecular mechanism of embryonic developmental defects and generates strategies to produce higher levels of UFAs.

In Silico Analysis of Fatty Acid Desaturases Structures in Camelina sativa, and Functional Evaluation of Csafad7 and Csafad8 on Seed Oil Formation and Seed Morphology.

Fatty acid desaturases add a second bond into a single bond of carbon atoms in fatty acid chains, resulting in an unsaturated bond between the two carbons. They are classified into soluble and membrane-bound desaturases, according to their structure, subcellular location, and function. The orthologous genes in Camelina sativa were identified and analyzed, and a total of 62 desaturase genes were identified. It was revealed that they had the common fatty acid desaturase domain, which has evolved separately, and the proteins of the same family also originated from the same ancestry. A mix of conserved, gained, or lost intron structure was obvious. Besides, conserved histidine motifs were found in each family, and transmembrane domains were exclusively revealed in the membrane-bound desaturases. The expression profile analysis of C. sativa desaturases revealed an increase in young leaves, seeds, and flowers. C. sativa ω3-fatty acid desaturases CsaFAD7 and CsaDAF8 were cloned and the subcellular localization analysis showed their location in the chloroplast. They were transferred into Arabidopsis thaliana to obtain transgenic lines. It was revealed that the ω3-fatty acid desaturase could increase the C18:3 level at the expense of C18:2, but decreases in oil content and seed weight, and wrinkled phenotypes were observed in transgenic CsaFAD7 lines, while no significant change was observed in transgenic CsaFAD8 lines in comparison to the wild-type. These findings gave insights into the characteristics of desaturase genes, which could provide an excellent basis for further investigation for C. sativa improvement, and overexpression of ω3-fatty acid desaturases in seeds could be useful in genetic engineering strategies, which are aimed at modifying the fatty acid composition of seed oil.

Genome-wide identification and functional analysis of oleosin genes in Brassica napus L.

BACKGROUND: Rapeseed is the third largest oil seed crop in the world. The seeds of this plant store lipids in oil bodies, and oleosin is the most important structural protein in oil bodies. However, the function of oleosin in oil crops has received little attention. RESULTS: In the present study, 48 oleosin sequences from the Brassica napus genome were identified and divided into four lineages (T, U, SH, SL). Synteny analysis revealed that most of the oleosin genes were conserved, and all of these genes experienced purifying selection during evolution. Three and four important oleosin genes from Arabidopsis and B. napus, respectively, were cloned and analyzed for function in Arabidopsis. Overexpression of these oleosin genes in Arabidopsis increased the seed oil content slightly, except for BnaOLE3. Further analysis revealed that the average oil body size of the transgenic seeds was slightly larger than that of the wild type (WT), except for BnaOLE1. The fatty acid profiles showed that the linoleic acid content (13.3% at most) increased and the peanut acid content (11% at most) decreased in the transgenic lines. In addition, the seed size and thousand-seed weight (TSW) also increased in the transgenic lines, which could lead to increased total lipid production. CONCLUSION: We identified oleosin genes in the B. napus genome, and overexpression of oleosin in Arabidopsis seeds increased the seed weight and linoleic acid content (13.3% at most).

Drought-responsive genes, late embryogenesis abundant group3 (LEA3) and vicinal oxygen chelate, function in lipid accumulation in Brassica napus and Arabidopsis mainly via enhancing photosynthetic efficiency and reducing ROS.

Drought is an abiotic stress that affects plant growth, and lipids are the main economic factor in the agricultural production of oil crops. However, the molecular mechanisms of drought response function in lipid metabolism remain little known. In this study, overexpression (OE) of different copies of the drought response genes LEA3 and VOC enhanced both drought tolerance and oil content in Brassica napus and Arabidopsis. Meanwhile, seed size, membrane stability and seed weight were also improved in OE lines. In contrast, oil content and drought tolerance were decreased in the AtLEA3 mutant (atlea3) and AtVOC-RNAi of Arabidopsis and in both BnLEA-RNAi and BnVOC-RNAi B. napus RNAi lines. Hybrids between two lines with increased or reduced expression (LEA3-OE with VOC-OE, atlea3 with AtVOC-RNAi) showed corresponding stronger trends in drought tolerance and lipid metabolism. Comparative transcriptomic analysis revealed the mechanisms of drought response gene function in lipid accumulation and drought tolerance. Gene networks involved in fatty acid (FA) synthesis and FA degradation were up- and down-regulated in OE lines, respectively. Key genes in the photosynthetic system and reactive oxygen species (ROS) metabolism were up-regulated in OE lines and down-regulated in atlea3 and AtVOC-RNAi lines, including LACS9, LIPASE1, PSAN, LOX2 and SOD1. Further analysis of photosynthetic and ROS enzymatic activities confirmed that the drought response genes LEA3 and VOC altered lipid accumulation mainly via enhancing photosynthetic efficiency and reducing ROS. The present study provides a novel way to improve lipid accumulation in plants, especially in oil production crops.

Genome-wide identification and Phylogenic analysis of kelch motif containing ACBP in Brassica napus.

BACKGROUND: Acyl-coA binding proteins (ACBPs) bind long chain acyl-CoA esters with very high affinity. Their possible involvement in fatty acid transportation from the plastid to the endoplasmic reticulum, prior to the formation of triacylglycerol has been suggested. Four classes of ACBPs were identified in Arabidopsis thaliana: the small ACBPs, the large ACBPs, the ankyrin repeats containing ACBPs and the kelch motif containing ACBPs. They differed in structure and in size, and showed multiple important functions. In the present study, Brassica napus ACBPs were identified and characterized. RESULTS: Eight copies of kelch motif ACBPs were cloned, it showed that B. napus ACBPs shared high amino acid sequence identity with A. thaliana, Brassica rapa and Brassica oleracea. Furthermore, phylogeny based on domain structure and comparison map showed the relationship and the evolution of ACBPs within Brassicaceae family: ACBPs evolved into four separate classes with different structure. Chromosome locations comparison showed conserved syntenic blocks. CONCLUSIONS: ACBPs were highly conserved in Brassicaceae. They evolved from a common ancestor, but domain duplication and rearrangement might separate them into four distinct classes, with different structure and functions. Otherwise, B. napus inherited kelch motif ACBPs from ancestor conserving chromosomal location, emphasizing preserved synteny block region. This study provided a first insight for exploring ACBPs in B. napus, which supplies a valuable tool for crop improvement in agriculture.

New insights into the genetic networks affecting seed fatty acid concentrations in Brassica napus.

BACKGROUND: Rapeseed (B. napus, AACC, 2n = 38) is one of the most important oil seed crops in the world, it is also one of the most common oil for production of biodiesel. Its oil is a mixture of various fatty acids and dissection of the genetic network for fatty acids biosynthesis is of great importance for improving seed quality. RESULTS: The genetic basis of fatty acid biosynthesis in B. napus was investigated via quantitative trail locus (QTL) analysis using a doubled haploid (DH) population with 202 lines. A total of 72 individual QTLs and a large number pairs of epistatic interactions associated with the content of 10 different fatty acids were detected. A total of 234 homologous genes of Arabidopsis thaliana that are involved in fatty acid metabolism were found within the confidence intervals (CIs) of 47 QTLs. Among them, 47 and 15 genes homologous to those of B. rapa and B. oleracea were detected, respectively. After the QTL mapping, the epistatic and the candidate gene interaction analysis, a potential regulatory pathway controlling fatty acid biosynthesis in B. napus was constructed, including 50 enzymes encoded genes and five regulatory factors (LEC1, LEC2, FUS3, WRI1 and ABI3). Subsequently, the interaction between these five regulatory factors and the genes involved in fatty acid metabolism were analyzed. CONCLUSIONS: In this study, a potential regulatory pathway controlling the fatty acid was constructed by QTL analysis and in silico mapping analysis. These results enriched our knowledge of QTLs for fatty acids metabolism and provided a new clue for genetic engineering fatty acids composition in B. napus.

Proteomic and comparative genomic analysis of two Brassica napus lines differing in oil content.

Ultrastructural observations, combined with proteomic and comparative genomic analyses, were applied to interpret the differences in protein composition and oil-body characteristics of mature seed of two Brassica napus lines with high and low oil contents of 55.19% and 36.49%, respectively. The results showed that oil bodies were arranged much closer in the high than in the low oil content line, and differences in cell size and thickness of cell walls were also observed. There were 119 and 32 differentially expressed proteins (DEPs) of total and oil-body proteins identified. The 119 DEPs of total protein were mainly involved in the oil-related, dehydration-related, storage and defense/disease, and some of these may be related to oil formation. The DEPs involved with dehydration-related were both detected in total and oil-body proteins for high and low oil lines and may be correlated with the number and size of oil bodies in the different lines. Some genes that corresponded to DEPs were confirmed by quantitative trait loci (QTL) mapping analysis for oil content. The results revealed that some candidate genes deduced from DEPs were located in the confidence intervals of QTL for oil content. Finally, the function of one gene that coded storage protein was verified by using a collection of Arabidopsis lines that can conditionally express the full length cDNA from developing seeds of B. napus.

New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus.

BACKGROUND: Increasing seed oil content is the most important breeding goal in Brassica napus, and phenotyping is crucial to dissect its genetic basis in crops. To date, QTL mapping for oil content has been based on whole seeds, and the lipid distribution is far from uniform in different tissues of seeds in B. napus. In this case, the phenotype based on whole seeds was unable to sufficiently reveal the complex genetic characteristics of seed oil content. RESULTS: Here, the three-dimensional (3D) distribution of lipid was determined for B. napus seeds by magnetic resonance imaging (MRI) and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues, the outer cotyledon (OC), inner cotyledon (IC), radicle (R) and seed coat (SC), which explained up to 13.76% of the phenotypic variation. Notably, 14 tissue-specific QTLs were reported for the first time, 7 of which were novel. Moreover, haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Furthermore, tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in the IC, OC and R than in the SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in the QTLs for OC and IC. CONCLUSIONS: The present study provides further insight into the genetic basis of seed oil content at the tissue-specific level.

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

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

Lysophosphatidic acid acyltransferase 2 and 5 commonly, but differently, promote seed oil accumulation in Brassica napus.

BACKGROUND: Increasing seed oil content (SOC) of Brassica napus has become one of the main plant breeding goals over the past decades. Lysophosphatidic acid acyltransferase (LPAT) performs an important molecular function by regulating the production of phosphatidic acid (PA), a key intermediate in the synthesis of membrane and storage lipids. However, the mechanism underlying the effect of LPAT on the SOC of B. napus remains unclear. RESULTS: In the present study, significant elevation of SOC was achieved by overexpressing BnLPAT2 and BnLPAT5 in B. napus. RNAi and CRISPR-Cas9 were also successfully used to knock down and knock out these two genes in B. napus where SOC significantly decreased. Meanwhile, we found an accumulation of lipid droplets and oil bodies in seeds of BnLPAT2 and BnLPAT5 overexpression lines, whereas an increase of sugar and protein in Bnlpat2 and Bnlpat5 mutant seeds. Sequential transcriptome analysis was further performed on the developing seeds of the BnLPAT2 and BnLPAT5 overexpression, knockdown, and knockout rapeseed lines. Most differentially expressed genes (DEGs) that were expressed in the middle and late stages of seed development were enriched in photosynthesis and lipid metabolism, respectively. The DEGs involved in fatty acid and lipid biosynthesis were active in the overexpression lines but were relatively inactive in the knockdown and knockout lines. Further analysis revealed that the biological pathways related to fatty acid/lipid anabolism and carbohydrate metabolism were specifically enriched in the BnLPAT2 overexpression lines. CONCLUSIONS: BnLPAT2 and BnLPAT5 are essential for seed oil accumulation. BnLPAT2 preferentially promoted diacylglycerol synthesis to increase SOC, whereas BnLPAT5 tended to boost PA synthesis for membrane lipid generation. Taken together, BnLPAT2 and BnLPAT5 can jointly but differently promote seed oil accumulation in B. napus. This study provides new insights into the potential mechanisms governing the promotion of SOC by BnLPAT2 and BnLPAT5 in the seeds of B. napus.

Genetic and Biochemical Investigation of Seed Fatty Acid Accumulation in Arabidopsis.

As a vegetable oil, consisting principally of triacylglycerols, is the major storage form of photosynthetically-fixed carbon in oilseeds which are of significant agricultural and industrial value. Photosynthesis in chlorophyll-containing green seeds, along with photosynthesis in leaves and other green organs, generates ATP and reductant (NADPH and NADH) needed for seed fatty acid production. However, contribution of seed photosynthesis to fatty acid accumulation in seeds have not been well-defined. Here, we report the contribution of seed-photosynthesis to fatty acid production by probing segregating green (photosynthetically-competent) and non-green or yellow (photosynthetically-non-competent) seeds in siliques of an Arabidopsis chlorophyll synthase mutant. Using this mutant, we found that yellow seeds lacking photosynthetic capacity reached 80% of amounts of oil in green seeds at maturity. Combining this with studies using shaded siliques, we determined that seed-photosynthesis accounts for 20% and silique and leaf/stem photosynthesis each account for ~40% of the ATP and reductant for seed oil production. Transmission electron microscopy (TEM) and pyridine nucleotides and ATP analyses revealed that seed photosynthesis provides ATP and reductant for oil production mostly during early development, as evidenced by delayed oil accumulation in non-green seeds. Transcriptomic analyses suggests that the oxidative pentose phosphate pathway could be the source of carbon, energy and reductants required for fatty acid synthesis beyond the early stages of seed development.

Δ6-Desaturase from Mortierella alpina: cDNA cloning, expression, and phylogenetic analysis.

The open reading frame of the Δ(6)-desaturase gene was isolated from Mortierella alpina W15 and the gene was cloned into a pPIC3.5K vector. The vector was transformed into Pichia pastoris GS115 and expression was induced with methanol. The Δ(6)-desaturase expressed in P. pastoris GS115 catalyzed the conversion of linoleic acid to γ-linolenic acid but not the conversion of α-linolenic acid to octadecatetraenoic acid. The results indicate that the Δ(6)-desaturase gene from M. alpina W15 has substrate specificity in different organisms. Phylogenetic analysis revealed that Δ(6)-desaturase genes can be divided into four monophyletic groups. This work paves the way for further study of the functions of Δ(6)-desaturase in fatty acid metabolism and its three-dimensional structure.

Author Correction: Integration of metabolome and transcriptome reveals flavonoid accumulation in the intergeneric hybrid between Brassica rapa and Raphanus sativus.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A convenient, rapid and efficient method for establishing transgenic lines of Brassica napus.

BACKGROUND: Brassica napus is an important oilseed crop that offers a considerable amount of biomass for global vegetable oil production. The establishment of an efficient genetic transformation system with a convenient transgenic-positive screening method is of great importance for gene functional analysis and molecular breeding. However, to our knowledge, there are few of the aforementioned systems available for efficient application in B. napus. RESULTS: Based on the well-established genetic transformation system in B. napus, five vectors carrying the red fluorescence protein encoding gene from Discosoma sp. (DsRed) were constructed and integrated into rapeseed via Agrobacterium-mediated hypocotyl transformation. An average of 59.1% tissues were marked with red fluorescence by the visual screening method in tissue culture medium, 96.1% of which, on average, were amplified with the objective genes from eight different rapeseed varieties. In addition, the final transgenic-positive efficiency of the rooted plantlets reached up to 90.7% from red fluorescence marked tissues, which was much higher than that in previous reports. Additionally, visual screening could be applicable to seedlings via integration of DsRed, including seed coats, roots, hypocotyls and cotyledons during seed germination. These results indicate that the highly efficient genetic transformation system combined with the transgenic-positive visual screening method helps to conveniently and efficiently obtain transgenic-positive rapeseed plantlets. CONCLUSION: A rapid, convenient and highly efficient method was developed to obtain transgenic plants, which can help to obtain the largest proportion of transgene-positive regenerated plantlets, thereby avoiding a long period of plant regeneration. The results of this study will benefit gene functional studies especially in high-throughput molecular biology research.

Effective editing for lysophosphatidic acid acyltransferase 2/5 in allotetraploid rapeseed (Brassica napus L.) using CRISPR-Cas9 system.

BACKGROUND: Brassica napus is one of the most important oilseed crops, and can supply considerable amounts of edible oil as well as provide raw materials for the production of biodiesel in the biotechnology industry. Lysophosphatidic acid acyltransferase (LPAT), a key enzyme in the Kennedy pathway, catalyses fatty acid chains into 3-phosphoglycerate and promotes further production of oil in the form of triacylglycerol. However, because B. napus is an allotetraploid with two subgenomes, the precise genes which involved in oil production remain unclear due to the intractability of efficiently knocking out all copies with high genetic redundancy. Therefore, a robust gene editing technology is necessary for gene function analysis. RESULTS: An efficient gene editing technology was developed for the allotetraploid plant B. napus using the CRISPR-Cas9 system. Previous studies showed poor results in either on-target or off-target activity in B. napus. In the present study, four single-gRNAs and two multi-gRNAs were deliberately designed from the conserved coding regions of BnLPAT2 which has seven homologous genes, and BnLPAT5, which has four homologous genes. The mutation frequency was found to range from 17 to 68%, while no mutation was observed in the putative off-target sites. The seeds of the Bnlpat2/Bnlpat5 mutant were wizened and showed enlarged oil bodies, disrupted distribution of protein bodies and increased accumulation of starch in mature seeds. The oil content decreased, with an average decrease of 32% for Bnlpat2 lines and 29% for Bnlpat5 lines in single-gRNA knockout lines, and a decline of 24% for Bnlpat2 mutant lines (i.e., g123) and 39% for Bnlpat2/Bnlpat5 double mutant lines (i.e., g134) in multi-gRNA knockout lines. CONCLUSIONS: Seven BnLPAT2 homologous genes and four BnLPAT5 homologous genes were cleaved completely using the CRISPR-Cas9 system, which indicated that it is effective for editing all homologous genes in allotetraploid rapeseed, despite the relatively low sequence identities of both gene families. The size of the oil bodies increased significantly while the oil content decreased, confirming that BnLPAT2 and BnLPAT5 play a role in oil biosynthesis. The present study lays a foundation for further oil production improvement in oilseed crop species.