Comparative sequence analysis for Brassica oleracea with similar sequences in B. rapa and Arabidopsis thaliana

The power of retrotransposons in high-throughput genotyping and sequencing

The use of molecular markers has become an essential part of molecular genetics through their application in numerous fields, which includes identification of genes associated with targeted traits, operation of backcrossing programs, modern plant breeding, genetic characterization, and marker-assisted selection. Transposable elements are a core component of all eukaryotic genomes, making them suitable as molecular markers. Most of the large plant genomes consist primarily of transposable elements; variations in their abundance contribute to most of the variation in genome size. Retrotransposons are widely present throughout plant genomes, and replicative transposition enables them to insert into the genome without removing the original elements. Various applications of molecular markers have been developed that exploit the fact that these genetic elements are present everywhere and their ability to stably integrate into dispersed chromosomal localities that are polymorphic within a species. The ongoing development of molecular marker technologies is directly related to the deployment of high-throughput genotype sequencing platforms, and this research is of considerable significance. In this review, the practical application to molecular markers, which is a use of technology of interspersed repeats in the plant genome were examined using genomic sources from the past to the present. Prospects and possibilities are also presented.

Genome-wide characterization of DNA methylation, small RNA expression, and histone H3 lysine nine di-methylation in Brassica rapa L

Epigenetic gene regulation is crucial to plant life and can involve dynamic interactions between various histone modifications, DNA methylation, and small RNAs. Detailed analysis of epigenome information is anticipated to reveal how the DNA sequence of the genome is translated into the plant’s phenotype. The aim of this study was to map the DNA methylation state at the whole genome level and to clarify the relationship between DNA methylation and transcription, small RNA expression, and histone H3 lysine 9 di-methylation (H3K9me2) in Brassica rapa. We performed whole genome bisulfite sequencing, small RNA sequencing, and chromatin immunoprecipitation sequencing using H3K9me2 antibody in a Chinese cabbage inbred line, RJKB-T24, and examined the impact of epigenetic states on transcription. Cytosine methylation in DNA was analysed in different sequence contexts (CG, CHG, and CHH) (where H could be A, C, or T) and position (promoter, exon, intron, terminator, interspersed repeat regions), and the H3K9me2 and 24 nucleotide small interfering RNAs (24 nt-siRNA) were overlaid onto the B. rapa reference genome. The epigenome was compared with that of Arabidopsis thaliana and the relationship between the position of DNA methylation and gene expression, and the involvement of 24 nt siRNAs and H3K9me2 are discussed.

Genetic characterization and fine mapping for multi-inflorescence in Brassica napus L

A major QTL for multi-inflorescence was mapped to a 27.18-kb region on A05 in Brassica napus by integrating QTL mapping, microarray analysis and whole-genome sequencing. Multi-inflorescence is a desirable trait for the genetic improvement of rapeseed (Brassica napus L.). However, the genetic mechanism underlying the multi-inflorescence trait is not well understood. In the present study, a doubled haploid (DH) population derived from a cross between single- and multi-inflorescence lines was investigated for the penetrance of multi-inflorescence across 3 years and genotyped with 257 simple sequence repeat and sequence-related amplified polymorphism loci. A major quantitative trait locus (QTL) for penetrance of multi-inflorescence was mapped to a 9.31-Mb region on chromosome A05, explaining 45.81% of phenotypic variance on average. Subsequently, 13 single-inflorescence and 15 multi-inflorescence DH lines were genotyped with the Brassica microarray, and the QTL interval of multi-inflorescence was narrowed to a 0.74-Mb region with 37 successive single nucleotide polymorphisms between single- and multi-inflorescence groups. A 27.18-kb QTL interval was detected by screening 420 recessive F₂ individuals with genome-specific markers. These results will be valuable for gene cloning and molecular breeding of multi-inflorescence in rapeseed.

Glucosinolate biosynthesis in Eruca sativa

Glucosinolates (GSLs) are a highly important group of secondary metabolites in the Caparalles order, both due to their significance in plant-biome interactions and to their chemoprotective properties. This study identified genes involved in all steps of aliphatic and indolic GSL biosynthesis in Eruca sativa, a cultivated plant closely related to Arabidopsis thaliana with agronomic and nutritional value. The impact of nitrogen (N) and sulfur (S) availability on GSL biosynthetic pathways at a transcriptional level, and on the final GSL content of plant leaf and root tissues, was investigated. N and S supply had a significant and interactive effect on the GSL content of leaves, in a structure-specific and tissue-dependent manner; the metabolites levels were significantly correlated with the relative expression of the genes involved in their biosynthesis. A more complex effect was observed in roots, where aliphatic and indolic GSLs and related biosynthetic genes responded differently to the various nutritional treatments suggesting that nitrogen and sulfur availability are important factors that control plant GSL content at a transcriptional level. The biological activity of extracts derived from these plants grown under the specific nutritional schemes was examined. N and S availability were found to significantly affect the cytotoxicity of E. sativa extracts on human cancer cells, supporting the notion that carefully designed nutritional schemes can promote the accumulation of chemoprotective substances in edible plants.

Dietary Sulforaphane in Cancer Chemoprevention: The Role of Epigenetic Regulation and HDAC Inhibition

SIGNIFICANCE

Sulforaphane, produced by the hydrolytic conversion of glucoraphanin after ingestion of cruciferous vegetables, particularly broccoli and broccoli sprouts, has been extensively studied due to its apparent health-promoting properties in disease and limited toxicity in normal tissue. Recent Studies: Recent identification of a sub-population of tumor cells with stem cell-like self-renewal capacity that may be responsible for relapse, metastasis, and resistance, as a potential target of the dietary compound, may be an important aspect of sulforaphane chemoprevention. Evidence also suggests that sulforaphane may target the epigenetic alterations observed in specific cancers, reversing aberrant changes in gene transcription through mechanisms of histone deacetylase inhibition, global demethylation, and microRNA modulation.

CRITICAL ISSUES

In this review, we discuss the biochemical and biological properties of sulforaphane with a particular emphasis on the anticancer properties of the dietary compound. Sulforaphane possesses the capacity to intervene in multistage carcinogenesis through the modulation and/or regulation of important cellular mechanisms. The inhibition of phase I enzymes that are responsible for the activation of pro-carcinogens, and the induction of phase II enzymes that are critical in mutagen elimination are well-characterized chemopreventive properties. Furthermore, sulforaphane mediates a number of anticancer pathways, including the activation of apoptosis, induction of cell cycle arrest, and inhibition of NFκB.

FUTURE DIRECTIONS

Further characterization of the chemopreventive properties of sulforaphane and its capacity to be selectively toxic to malignant cells are warranted to potentially establish the clinical utility of the dietary compound as an anti-cancer compound alone, and in combination with clinically relevant therapeutic and management strategies.

High-density single nucleotide polymorphism (SNP) array mapping in Brassica oleracea: identification of QTL associated with carotenoid variation in broccoli florets

A high-resolution genetic linkage map of B. oleracea was developed from a B. napus SNP array. The work will facilitate genetic and evolutionary studies in Brassicaceae. A broccoli population, VI-158 × BNC, consisting of 150 F2:3 families was used to create a saturated Brassica oleracea (diploid: CC) linkage map using a recently developed rapeseed (Brassica napus) (tetraploid: AACC) Illumina Infinium single nucleotide polymorphism (SNP) array. The map consisted of 547 non-redundant SNP markers spanning 948.1 cM across nine chromosomes with an average interval size of 1.7 cM. As the SNPs are anchored to the genomic reference sequence of the rapid cycling B. oleracea TO1000, we were able to estimate that the map provides 96 % coverage of the diploid genome. Carotenoid analysis of 2 years data identified 3 QTLs on two chromosomes that are associated with up to half of the phenotypic variation associated with the accumulation of total or individual compounds. By searching the genome sequences of the two related diploid species (B. oleracea and B. rapa), we further identified putative carotenoid candidate genes in the region of these QTLs. This is the first description of the use of a B. napus SNP array to rapidly construct high-density genetic linkage maps of one of the constituent diploid species. The unambiguous nature of these markers with regard to genomic sequences provides evidence to the nature of genes underlying the QTL, and demonstrates the value and impact this resource will have on Brassica research.

Comparative genomics of Brassicaceae crops

The family Brassicaceae is one of the major groups of the plant kingdom and comprises diverse species of great economic, agronomic and scientific importance, including the model plant Arabidopsis. The sequencing of the Arabidopsis genome has revolutionized our knowledge in the field of plant biology and provides a foundation in genomics and comparative biology. Genomic resources have been utilized in Brassica for diversity analyses, construction of genetic maps and identification of agronomic traits. In Brassicaceae, comparative sequence analysis across the species has been utilized to understand genome structure, evolution and the detection of conserved genomic segments. In this review, we focus on the progress made in genetic resource development, genome sequencing and comparative mapping in Brassica and related species. The utilization of genomic resources and next-generation sequencing approaches in improvement of Brassica crops is also discussed.

Identification of suitable qPCR reference genes in leaves of Brassica oleracea under abiotic stresses

Real-time quantitative PCR is nowadays a standard method to study gene expression variations in various samples and experimental conditions. However, to interpret results accurately, data normalization with appropriate reference genes appears to be crucial. The present study describes the identification and the validation of suitable reference genes in Brassica oleracea leaves. Expression stability of eight candidates was tested following drought and cold abiotic stresses by using three different softwares (BestKeeper, NormFinder and geNorm). Four genes (BolC.TUB6, BolC.SAND1, BolC.UBQ2 and BolC.TBP1) emerged as the most stable across the tested conditions. Further gene expression analysis of a drought- and a cold-responsive gene (BolC.DREB2A and BolC.ELIP, respectively), confirmed the stability and the reliability of the identified reference genes when used for normalization in the leaves of B. oleracea. These four genes were finally tested upon a benzene exposure and all appeared to be useful reference genes along this toxicological condition. These results provide a good starting point for future studies involving gene expression measurement on leaves of B. oleracea exposed to environmental modifications.

High-throughput sequencing identification of novel and conserved miRNAs in the Brassica oleracea leaves

Background

Plant microRNAs are short (~21 nt) non-coding molecules that regulate gene expression by targeting the mRNA cleavage or protein translation inhibition. In this manner, they play many important roles in the cells of living organisms. One of the plant species in which the entire set of miRNAs has not been yet completely identified is Brassica oleracea var. capitata (cabbage). For this reason and for the economic and nutritional importance of this food crop, high-throughput small RNAs sequencing has been performed to discover the novel and conserved miRNAs in mature cabbage leaves.

Results

In this study, raw reads generated from three small RNA libraries were bioinformatically processed and further analyzed to select sequences homologous to known B. oleracea and other plant miRNAs. As a result of this analysis, 261 conserved miRNAs (belonging to 62 families) have been discovered. MIR169, MIR167 and MIR166 were the largest miRNA families, while the highest abundance molecules were miR167, miR166, miR168c and miR157a. Among the generated sequencing reads, miRNAs* were also found, such as the miR162c*, miR160a* and miR157a*. The unannotated tags were used in the prediction and evaluation of novel miRNAs, which resulted in the 26 potential miRNAs proposal. The expressions of 13 selected miRNAs were analyzed by northern blot hybridization. The target prediction and annotation for identified miRNAs were performed, according to which discovered molecules may target mRNAs encoding several potential proteins – e.g., transcription factors, polypeptides that regulate hormone stimuli and abiotic stress response, and molecules participating in transport and cell communication. Additionally, KEGG maps analysis suggested that the miRNAs in cabbage are involved in important processing pathways, including glycolysis, glycerolipid metabolism, flavonoid biosynthesis and oxidative phosphorylation.

Conclusions

Conclusively, for the first time, the large set of miRNAs was identified in mature cabbage leaves. Potential targets designation for these miRNAs may suggest their essential role in many plants primary biological processes. Presented study not only supplements the knowledge about B. oleracea miRNAs, but additionally it may be used in other research concerning the improvement of the cabbage cultivation.

Novel bioresources for studies of Brassica oleracea: identification of a kale MYB transcription factor responsible for glucosinolate production

Plants belonging to the Brassicaceae family exhibit species-specific profiles of glucosinolates (GSLs), a class of defence compounds against pathogens and insects. GSLs also exhibit various human health-promoting properties. Among them, glucoraphanin (aliphatic 4-methylsulphinylbutyl GSL) has attracted the most attention because it hydrolyses to form a potent anticancer compound. Increased interest in developing commercial varieties of Brassicaceae crops with desirable GSL profiles has led to attempts to identify genes that are potentially valuable for controlling GSL biosynthesis. However, little attention has been focused on genes of kale (Brassica oleracea var. acephala). In this study, we established full-length kale cDNA libraries containing 59 904 clones, which were used to generate an expressed sequence tag (EST) data set with 119 204 entries. The EST data set clarified genes related to the GSL biosynthesis pathway in kale. We specifically focused on BoMYB29, a homolog of Arabidopsis MYB29/PMG2/HAG3, not only to characterize its function but also to demonstrate its usability as a biological resource. BoMYB29 overexpression in wild-type Arabidopsis enhanced the expression of aliphatic GSL biosynthetic genes and the accumulation of aliphatic GSLs. When expressed in the myb28myb29 mutant, which exhibited no detectable aliphatic GSLs, BoMYB29 restored the expression of biosynthetic genes and aliphatic GSL accumulation. Interestingly, the ratio of methylsulphinyl GSL content, including glucoraphanin, to that of methylthio GSLs was greatly increased, indicating the suitability of BoMYB29 as a regulator for increasing methylsulphinyl GSL content. Our results indicate that these biological resources can facilitate further identification of genes useful for modifications of GSL profiles and accumulation in kale.

Promoter variation and transcript divergence in Brassicaceae lineages of FLOWERING LOCUS T

Brassica napus (AACC, 2n = 38), an oil crop of world-wide importance, originated from interspecific hybridization of B. rapa (AA, 2n = 20) and B. oleracea (CC, 2n = 18), and has six FLOWERING LOCUS T (FT) paralogues. Two located on the homeologous chromosomes A2 and C2 arose from a lineage distinct from four located on A7 and C6. A set of three conserved blocks A, B and C, which were found to be essential for FT activation by CONSTANS (CO) in Arabidopsis, was identified within the FT upstream region in B. napus and its progenitor diploids. However, on chromosome C2, insertion of a DNA transposable element (TE) and a retro-element in FT upstream blocks A and B contributed to significant structural divergence between the A and C genome orthologues. Phylogenetic analysis of upstream block A indicated the conserved evolutionary relationships of distinct FT genes within Brassicaceae. We conclude that the ancient At-α whole genome duplication contributed to distinct ancestral lineages for this key adaptive gene, which co-exist within the same genus. FT-A2 was found to be transcribed in all leaf samples from different developmental stages in both B. rapa and B. napus, whereas FT-C2 was not transcribed in either B. napus or B. oleracea. Silencing of FT-C2 appeared to result from TE insertion and consequent high levels of cytosine methylation in TE sequences within upstream block A. Interestingly, FT-A7/C6 paralogues were specifically silenced in winter type B. napus but abundantly expressed in spring type cultivars under vernalization-free conditions. Motif prediction indicated the presence of two CO protein binding sites within all Brassica block A and additional sites for FT activation in block C. We propose that the ancestral whole genome duplications have contributed to more complex mechanisms of floral regulation and niche adaptation in Brassica compared to Arabidopsis.

Molecular and functional characterization of broccoli EMBRYONIC FLOWER 2 genes

Polycomb group (PcG) proteins regulate major developmental processes in Arabidopsis. EMBRYONIC FLOWER 2 (EMF2), the VEFS domain-containing PcG gene, regulates diverse genetic pathways and is required for vegetative development and plant survival. Despite widespread EMF2-like sequences in plants, little is known about their function other than in Arabidopsis and rice. To study the role of EMF2 in broccoli (Brassica oleracea var. italica cv. Elegance) development, we identified two broccoli EMF2 (BoEMF2) genes with sequence homology to and a similar gene expression pattern to that in Arabidopsis (AtEMF2). Reducing their expression in broccoli resulted in aberrant phenotypes and gene expression patterns. BoEMF2 regulates genes involved in diverse developmental and stress programs similar to AtEMF2 in Arabidopsis. However, BoEMF2 differs from AtEMF2 in the regulation of flower organ identity, cell proliferation and elongation, and death-related genes, which may explain the distinct phenotypes. The expression of BoEMF2.1 in the Arabidopsis emf2 mutant (Rescued emf2) partially rescued the mutant phenotype and restored the gene expression pattern to that of the wild type. Many EMF2-mediated molecular and developmental functions are conserved in broccoli and Arabidopsis. Furthermore, the restored gene expression pattern in Rescued emf2 provides insights into the molecular basis of PcG-mediated growth and development.

BraSto, a Stowaway MITE from Brassica: recently active copies preferentially accumulate in the gene space

We characterized a Brassica miniature inverted repeat transposable element (MITE) from the Stowaway superfamily, designated BraSto (Bra ssica Sto waway). BraSto copy number was assessed using real-time quantitative PCR in the two diploid species B. rapa (genome A) and B. oleracea (genome C) and the corresponding allotetraploid species B. napus (genome AC). Phylogenetic relationships among a set of 131 BraSto copies were then analyzed. BraSto appears to have been only moderately amplified in the Brassica genome and was still active recently with marks of proliferation in both diploid Brassica species, which diverged 3.75 million years ago, but also in the allotetraploid species after reuniting of the two differentiated genomes. We characterized insertion sites for low-divergence BraSto copies among the gene space of the B. rapa genome using bioinformatics approaches. For BraSto copies localized nearby or within genes, we observed frequent associations of BraSto with putative promoters and regulatory regions of genes, but exclusion from coding regions. In addition, BraSto was significantly similar to several Brassica expressed sequence tags (ESTs), including stress-induced ESTs. We also demonstrated the enrichment of BraSto sequences in binding sites for transcription factors and other regulatory elements. Our results lead to the question of a role for BraSto in the regulation of gene expression: this putative role, if further confirmed experimentally, would help to obtain a new insight into the significance of MITEs in the functional plant genome.

Characterization of seedling proteomes and development of markers to distinguish the Brassica A and C genomes

The diploid species Brassica rapa (genome AA) and B. oleracea (genome CC) were compared by full-scale proteome analyses of seedling. A total of 28.2% of the proteins was common to both species, indicating the existence of a basal or ubiquitous proteome. However, a number of discriminating proteins (32.0%) and specific proteins (39.8%) of the Brassica A and C genomes, respectively, were identified, which could represent potentially species-specific functions. Based on these A or C genome-specific proteins, a number of PCR-based markers to distinguish B. rapa and B. oleracea species were also developed.

QTL analysis of leaf architecture

Leaf size and shape define the photosynthetic capability of a plant and have a significant impact on important agronomic traits, such as yield, quality, disease resistance and stress responses. Cultivated varieties of many plant species show remarkable variations in leaf morphology. Such variation usually exhibits a continuous phenotypic distribution and is controlled by the interaction of multiple genes known as quantitative trait loci (QTL). Here, we review several studies that evaluate natural variations in the leaf morphologies of crop species, as well as in the model plant Arabidopsis thaliana. The use of high-throughput, genome-wide approaches such as transcriptomics and metabolomics is helping to identify the nucleotide polymorphism(s) responsible for the phenotypic differences attributed to some of these QTL.

A Brassica rapa linkage map of EST-based SNP markers for identification of candidate genes controlling flowering time and leaf morphological traits

For identification of genes responsible for varietal differences in flowering time and leaf morphological traits, we constructed a linkage map of Brassica rapa DNA markers including 170 EST-based markers, 12 SSR markers, and 59 BAC sequence-based markers, of which 151 are single nucleotide polymorphism (SNP) markers. By BLASTN, 223 markers were shown to have homologous regions in Arabidopsis thaliana, and these homologous loci covered nearly the whole genome of A. thaliana. Synteny analysis between B. rapa and A. thaliana revealed 33 large syntenic regions. Three quantitative trait loci (QTLs) for flowering time were detected. BrFLC1 and BrFLC2 were linked to the QTLs for bolting time, budding time, and flowering time. Three SNPs in the promoter, which may be the cause of low expression of BrFLC2 in the early-flowering parental line, were identified. For leaf lobe depth and leaf hairiness, one major QTL corresponding to a syntenic region containing GIBBERELLIN 20 OXIDASE 3 and one major QTL containing BrGL1, respectively, were detected. Analysis of nucleotide sequences and expression of these genes suggested possible involvement of these genes in leaf morphological traits.

Development and evaluation of single-nucleotide polymorphism markers in allotetraploid rapeseed (Brassica napus L.)

Single-nucleotide polymorphisms (SNPs) and insertion-deletions (INDELs) are currently the important classes of genetic markers for major crop species. In this study, methods for developing SNP markers in rapeseed (Brassica napus L.) and their in silico mapping and use for genotyping are demonstrated. For the development of SNP and INDEL markers, 181 fragments from 121 different gene sequences spanning 86 kb were examined. A combination of different selection methods (genome-specific amplification, hetero-duplex analysis and sequence analysis) allowed the detection of 18 singular fragments that showed a total of 87 SNPs and 6 INDELs between 6 different rapeseed varieties. The average frequency of sequence polymorphism was estimated to be one SNP every 247 bp and one INDEL every 3,583 bp. Most SNPs and INDELs were found in non-coding regions. Polymorphism information content values for SNP markers ranged between 0.02 and 0.50 in a set of 86 varieties. Using comparative genetics data for B. napus and Arabidopsis thaliana, an allocation of SNP markers to linkage groups in rapeseed was achieved: a unique location was determined for seven gene sequences; two and three possible locations were found for six and four sequences, respectively. The results demonstrate the usefulness of existing genomic resources for SNP discovery in rapeseed.