Genome-wide DNA methylation profiling by modified reduced representation bisulfite sequencing in Brassica rapa suggests that epigenetic modifications play a key role in polyploid genome evolution

Brassica rapa includes some of the most important vegetables worldwide as well as oilseed crops. The complete annotated genome sequence confirmed its paleohexaploid origins and provides opportunities for exploring the detailed process of polyploid genome evolution. We generated a genome-wide DNA methylation profile for B. rapa using a modified reduced representation bisulfite sequencing (RRBS) method. This sampling represented 2.24% of all CG loci (2.5 x 105), 2.16% CHG (2.7 x 105) and 1.68% CHH loci (1.05 x 105) (where H = A, T or C). Our sampling of DNA methylation in B. rapa indicated that 52.4% of CG sites were present as 5mCG, with 31.8% of CHG and 8.3% of CHH. It was found that genic regions of single copy genes had significantly higher methylation compared to those of two or three copy genes. Differences in degree of genic DNA methylation were observed in a hierarchical relationship corresponding to the relative age of the three ancestral subgenomes, primarily accounted by single-copy genes. RNA-seq analysis revealed that overall the level of transcription was negatively correlated with mean gene methylation content and depended on copy number or associated with the different subgenomes. These results provide new insights into the role epigenetic variation plays in polyploid genome evolution, and suggest an alternative mechanism for duplicate gene loss

Universal endogenous gene controls for bisulphite conversion in analysis of plant DNA methylation

Accurate analysis of DNA methylation by bisulphite sequencing depends on the complete conversion of all cytosines into uracil. Until now there has been no standard or universal gene identified as an endogenous control to monitor the conversion frequency in plants. Here, we report the development of PCR based assays for one nuclear gene IND (INDEHISCENT) and two mitochondrial genes, NAD (NICOTINAMIDE ADENINE DINUCLEOTIDE) and ATP1 (ATPase SUBUNIT 1). We demonstrated their efficacy as bisulphite conversion controls in Brassica and other plant taxa. The target regions amplified by four primer pairs were found to be consistently free from DNA methylation. Primer pairs for IND.a and NAD were effective within Brassica species, whereas two primer pairs for ATP1 provided reliable controls across a representative range of dicot and monocot angiosperm species. These primer sets may therefore be adopted as controls in plant methylation analysis for a wide range of studies

Genetic Properties of a Nested Association Mapping Population Constructed With Semi-Winter and Spring Oilseed Rapes

Nested association mapping (NAM) populations have been widely applied to dissect the genetic basis of complex quantitative traits in a variety of crops. In this study, we developed a Brassica napus NAM (BN-NAM) population consisting of 15 recombination inbred line (RIL) families with 2,425 immortal genotypes. Fifteen high-density genetic linkage maps were constructed by genotyping by sequencing (GBS) based on all RIL families, with further integration into a joint linkage map (JLM) having 30,209 unique markers in common with multiple linkage maps. Furthermore, an ultra-density whole-genome variation map was constructed by projecting 4,444,309 high-quality variants onto the JLM. The NAM population captured a total of 88,542 recombination events (REs). The uneven distribution of recombination rate along chromosomes is positively correlated with the densities of genes and markers, but negatively correlated with the density of transposable elements and linkage disequilibrium (LD). Analyses of population structure and principal components revealed that the BN-NAM population could be divided into three groups with weak stratification. The LD decay distance across genome varied between 170 and 2,400 Kb, with LD decay more rapid in the A than in the C sub-genome. The pericentromeric regions contained large LD blocks, especially in the C sub-genome. This NAM population provides a valuable resource for dissecting the genetic basis of important traits in rapeseed, especially in semi-winter oilseed rape

Development and Validation of an Effective CRISPR/Cas9 Vector for Efficiently Isolating Positive Transformants and Transgene-Free Mutants in a Wide Range of Plant Species

The CRISPR/Cas9 technique is a highly valuable tool in creating new materials for both basic and applied researches. Previously, we succeeded in effectively generating mutations in Brassica napus using an available CRISPR/Cas9 vector pKSE401, while isolation of Cas9-free mutants is laborious and inefficient. Here, we inserted a fluorescence tag (sGFP) driven by the constitutive 35S promoter into pKSE401 to facilitate a visual screen of mutants. This modified vector was named pKSE401G and tested in several dicot plant species, including Arabidopsis, B. napus, Fragaria vesca (strawberry), and Glycine max (soybean). Consequently, GFP-positive plants were readily identified through fluorescence screening in all of these species. Among these GFP-positive plants, the average mutation frequency ranged from 20.4 to 52.5% in Arabidopsis and B. napus with stable transformation, and was 90.0% in strawberry and 75.0% in soybean with transient transformation, indicating that the editing efficiency resembles that of the original vector. Moreover, transgene-free mutants were sufficiently identified in Arabidopsis in the T2 generation and B. napus in the T1 generation based on the absence of GFP fluorescence, and these mutants were stably transmissible to next generation without newly induced mutations. Collectively, pKSE401G provides us an effective tool to readily identify positive primary transformants and transgene-free mutants in later generations in a wide range of dicot plant species

Construction of an integrated genetic linkage map for the A genome of Brassica napus using SSR markers derived from sequenced BACs in B. rapa

Background The Multinational Brassica rapa Genome Sequencing Project (BrGSP) has developed valuable genomic resources, including BAC libraries, BAC-end sequences, genetic and physical maps, and seed BAC sequences for Brassica rapa. An integrated linkage map between the amphidiploid B. napus and diploid B. rapa will facilitate the rapid transfer of these valuable resources from B. rapa to B. napus (Oilseed rape, Canola). Results In this study, we identified over 23,000 simple sequence repeats (SSRs) from 536 sequenced BACs. 890 SSR markers (designated as BrGMS) were developed and used for the construction of an integrated linkage map for the A genome in B. rapa and B. napus. Two hundred and nineteen BrGMS markers were integrated to an existing B. napus linkage map (BnaNZDH). Among these mapped BrGMS markers, 168 were only distributed on the A genome linkage groups (LGs), 18 distrubuted both on the A and C genome LGs, and 33 only distributed on the C genome LGs. Most of the A genome LGs in B. napus were collinear with the homoeologous LGs in B. rapa, although minor inversions or rearrangements occurred on A2 and A9. The mapping of these BAC-specific SSR markers enabled assignment of 161 sequenced B. rapa BACs, as well as the associated BAC contigs to the A genome LGs of B. napus. Conclusion The genetic mapping of SSR markers derived from sequenced BACs in B. rapa enabled direct links to be established between the B. napus linkage map and a B. rapa physical map, and thus the assignment of B. rapa BACs and the associated BAC contigs to the B. napus linkage map. This integrated genetic linkage map will facilitate exploitation of the B. rapa annotated genomic resources for gene tagging and map-based cloning in B. napus, and for comparative analysis of the A genome within Brassica species

Identification of Brassica napus lysophosphatidylcholine acyltransferase genes through yeast functional screening

Acyl-CoA:lysophosphatidylcholine acyltransferase (LPCAT), which acylates lysophosphatidylcholine (LPC) to produce phosphatidylcholine (PC), is a key enzyme in the Lands cycle. There is evidence that acyl exchange involving LPCAT is a prevailing metabolic process during triacylglycerol (TAG) synthesis in seeds. In this study, by complementing the yeast lca1\u394 mutant deficient in LPCAT activity with an Arabidopsis seedling cDNA library, it was found that the previously reported lysophospholipid acyltransferases (LPLATs), At1g12640 and At1g63050, were the only two acyltransferase genes that restored hyposensitivity of the lca1\u394 mutant to lyso-platelet-activating factor (lyso-PAF). A developing seed cDNA library from Brassica napus L. cv Hero was constructed to further explore the heterologous yeast complementation approach. Three B. napusLPCAT homologs were identified, of which BnLPCAT1-1 and BnLPCAT1-2 are orthologous to ArabidopsisAtLPLAT1 (At1g12640) while BnLPCAT2 is an ortholog of AtLPLAT2 (At1g63050). The proteins encoded by BnLPCAT1-1 and BnLPCAT2 were chosen for further study. Enzymatic assays demonstrated that both proteins exhibited a substrate preference for LPCs and unsaturated fatty acyl-CoAs. In addition to the enzymatic properties of plant lysophosphatidylcholine acyltransferases uncovered in this study, this report describes a useful technique that facilitates subsequent analyses into the role of LPCATs in PC turnover and seed oil biosynthesis.Peer reviewed: YesNRC publication: Ye

CCS52A2/FZR1, a cell cycle regulator, is an essential factor for shoot apical meristem maintenance in <it>Arabidopsis thaliana</it>

Abstract Background Cell division and cell fate decisions regulate organ formation and function in plant growth and development. It is still unclear how specific meristematic regulatory networks operate with the cell cycle machinery to translate stem cell identity and maintenance into cellular behavior. In this study, we address these questions by analysis of a shoot apex defective mutant, namely xcm9. Results Phenotypic analysis of the xcm9 mutant reveals concomitant premature termination of floral shoots with frequent bifurcation of the shoot apices, stems, and flowers. Microscopic observations show irregular cell organization in shoot apical meristems of xcm9. Positional cloning revealed that xcm9 is a loss of function allele of the CCS52A2/FZR1 gene, which has previously been implicated in root development. Expression analysis demonstrated that CCS52A2 maintains a higher transcriptional expression level in actively dividing tissue. Genetic studies indicated that the CCS52A2 gene functions together with WUSCHEL (WUS) and CLAVATA3 (CLV3) in regulating the development of the shoot meristem, and also contributes to this regulation together with the chromatin remodeling pathway. In addition, fewer xcm9 cells express CYCLIN B1:1, showing that cell cycle progression is disrupted in the mutant. Conclusion We propose that the CCS52A2 gene is a mediator that functions together with meristematic genes to regulate meristem organization, and cross-functions with chromatin regulators in cell cycle progression during shoot apical meristem development.</p