Moderate Salinity Stress Increases the Seedling Biomass in Oilseed Rape (Brassica napus L.).

Oilseed rape (Brassica napus L.), an important oil crop of the world, suffers various abiotic stresses including salinity stress during the growth stage. While most of the previous studies paid attention to the adverse effects of high salinity stress on plant growth and development, as well as their underlying physiological and molecular mechanisms, less attention was paid to the effects of moderate or low salinity stress. In this study, we first tested the effects of different concentrations of NaCl solution on the seedling growth performance of two oilseed rape varieties (CH336, a semi-winter type, and Bruttor, a spring type) in pot cultures. We found that moderate salt concentrations (25 and 50 mmol L-1 NaCl) can stimulate seedling growth by a significant increase (10~20%, compared to controls) in both above- and underground biomasses, as estimated at the early flowering stage. We then performed RNA-seq analyses of shoot apical meristems (SAMs) from six-leaf-aged seedlings under control (CK), low (LS, 25 mmol L-1), and high (HS, 180 mmol L-1) salinity treatments in the two varieties. The GO and KEGG enrichment analyses of differentially expressed genes (DEGs) demonstrated that such a stimulating effect on seedling growth by low salinity stress may be caused by a more efficient capacity for photosynthesis as compensation, accompanied by a reduced energy loss for the biosynthesis of secondary metabolites and redirecting of energy to biomass formation. Our study provides a new perspective on the cultivation of oilseed rape in saline regions and new insights into the molecular mechanisms of salt tolerance in Brassica crops. The candidate genes identified in this study can serve as targets for molecular breeding selection and genetic engineering toward enhancing salt tolerance in B. napus.

Structural variations and environmental specificities of flowering time-related genes in Brassica napus.

KEY MESSAGE: We found that the flowering time order of accessions in a genetic population considerably varied across environments, and homolog copies of essential flowering time genes played different roles in different locations. Flowering time plays a critical role in determining the life cycle length, yield, and quality of a crop. However, the allelic polymorphism of flowering time-related genes (FTRGs) in Brassica napus, an important oil crop, remains unclear. Here, we provide high-resolution graphics of FTRGs in B. napus on a pangenome-wide scale based on single nucleotide polymorphism (SNP) and structural variation (SV) analyses. A total of 1337 FTRGs in B. napus were identified by aligning their coding sequences with Arabidopsis orthologs. Overall, 46.07% of FTRGs were core genes and 53.93% were variable genes. Moreover, 1.94%, 0.74%, and 4.49% FTRGs had significant presence-frequency differences (PFDs) between the spring and semi-winter, spring and winter, and winter and semi-winter ecotypes, respectively. SNPs and SVs across 1626 accessions of 39 FTRGs underlying numerous published qualitative trait loci were analyzed. Additionally, to identify FTRGs specific to an eco-condition, genome-wide association studies (GWASs) based on SNP, presence/absence variation (PAV), and SV were performed after growing and observing the flowering time order (FTO) of plants in a collection of 292 accessions at three locations in two successive years. It was discovered that the FTO of plants in a genetic population changed a lot across various environments, and homolog copies of some key FTRGs played different roles in different locations. This study revealed the molecular basis of the genotype-by-environment (G × E) effect on flowering and recommended a pool of candidate genes specific to locations for breeding selection.

Genome-wide-association study and transcriptome analysis reveal the genetic basis controlling the formation of leaf wax in Brassica napus.

Cuticular wax protects plants from various biotic and abiotic stresses. However, the genetic network of wax biosynthesis and the environmental factors influencing leaf wax production in rapeseed (Brassica napus) remains unclear. Here, we demonstrated the role of leaf wax in the resistance to Sclerotinia infection in rapeseed. We found that leaves grown under high light intensity had higher expression of genes involved in wax biosynthesis, and produced more wax on the leaf surface, compared with those grown under low light conditions. Genome-wide association study (GWAS) identified 89 single nucleotide polymorphisms significantly associated with leaf wax coverage. A cross-analysis between GWAS and differentially expressed genes (DEGs) in the leaf epidermis of the accessions with contrasting differences in wax content revealed 17 candidate genes that control this variation in rapeseed. Selective sweep analysis combined with DEG analysis unveiled 510 candidate genes with significant selective signatures. From the candidate genes, we selected BnaA02.LOX4, a putative lipoxygenase, and BnaCnn.CER1, BnaA02.CER3, BnaC02.CER3, and BnaA01.CER4 (ECERIFERUM1-4) that were putatively responsible for wax biosynthesis, to analyse the allelic forms and haplotypes corresponding to high or low leaf wax coverage. These data enrich our knowledge about wax formation, and provide a gene pool for breeding an ideal leaf wax content in rapeseed.

Evolutionary Analysis of the YABBY Gene Family in Brassicaceae.

The YABBY gene family is one of the plant transcription factors present in all seed plants. The family members were extensively studied in various plants and shown to play important roles in plant growth and development, such as the polarity establishment in lateral organs, the formation and development of leaves and flowers, and the response to internal plant hormone and external environmental stress signals. In this study, a total of 364 YABBY genes were identified from 37 Brassicaceae genomes, of which 15 were incomplete due to sequence gaps, and nine were imperfect (missing C2C2 zinc-finger or YABBY domain) due to sequence mutations. Phylogenetic analyses resolved these YABBY genes into six compact clades except for a YAB3-like gene identified in Aethionema arabicum. Seventeen Brassicaceae species each contained a complete set of six basic YABBY genes (i.e., 1 FIL, 1 YAB2, 1 YAB3, 1 YAB5, 1 INO and 1 CRC), while 20 others each contained a variable number of YABBY genes (5-25) caused mainly by whole-genome duplication/triplication followed by gene losses, and occasionally by tandem duplications. The fate of duplicate YABBY genes changed considerably according to plant species, as well as to YABBY gene type. These YABBY genes were shown to be syntenically conserved across most of the Brassicaceae species, but their functions might be considerably diverged between species, as well as between paralogous copies, as demonstrated by the promoter and expression analysis of YABBY genes in two Brassica species (B. rapa and B. oleracea). Our study provides valuable insights for understanding the evolutionary story of YABBY genes in Brassicaceae and for further functional characterization of each YABBY gene across the Brassicaceae species.

Genome wide survey, evolution and expression analysis of PHD finger genes reveal their diverse roles during the development and abiotic stress responses in Brassica rapa L.

BACKGROUND: Plant homeodomain (PHD) finger proteins are widely present in all eukaryotes and play important roles in chromatin remodeling and transcriptional regulation. The PHD finger can specifically bind a number of histone modifications as an "epigenome reader", and mediate the activation or repression of underlying genes. Many PHD finger genes have been characterized in animals, but only few studies were conducted on plant PHD finger genes to this day. Brassica rapa (AA, 2n = 20) is an economically important vegetal, oilseed and fodder crop, and also a good model crop for functional and evolutionary studies of important gene families among Brassica species due to its close relationship to Arabidopsis thaliana. RESULTS: We identified a total of 145 putative PHD finger proteins containing 233 PHD domains from the current version of B. rapa genome database. Gene ontology analysis showed that 67.7% of them were predicted to be located in nucleus, and 91.3% were predicted to be involved in protein binding activity. Phylogenetic, gene structure, and additional domain analyses clustered them into different groups and subgroups, reflecting their diverse functional roles during plant growth and development. Chromosomal location analysis showed that they were unevenly distributed on the 10 B. rapa chromosomes. Expression analysis from RNA-Seq data showed that 55.7% of them were constitutively expressed in all the tested tissues or organs with relatively higher expression levels reflecting their important housekeeping roles in plant growth and development, while several other members were identified as preferentially expressed in specific tissues or organs. Expression analysis of a subset of 18 B. rapa PHD finger genes under drought and salt stresses showed that all these tested members were responsive to the two abiotic stress treatments. CONCLUSIONS: Our results reveal that the PHD finger genes play diverse roles in plant growth and development, and can serve as a source of candidate genes for genetic engineering and improvement of Brassica crops against abiotic stresses. This study provides valuable information and lays the foundation for further functional determination of PHD finger genes across the Brassica species.

Comprehensive genomic survey, structural classification and expression analysis of C2H2 zinc finger protein gene family in Brassica rapa L.

C2H2 zinc finger protein (ZFP) genes have been extensively studied in many organisms and can function as transcription factors and be involved in many biological processes including plant growth and development and stress responses. In the current study, a comprehensive genomics analysis of the C2H2-ZFP genes in B. rapa was performed. A total of 301 B. rapa putative C2H2-ZFP (BrC2H2-ZFP) genes were identified from the available Brassica genome databases, and further characterized through analysis of conserved amino acid residues in C2H2-ZF domains and their organization, subcellular localization, phylogeny, additional domain, chromosomal location, synteny relationship, Ka/Ks ratio, and expression pattern. We also analyzed the expression patterns of eight B. rapa C2H2-ZFP genes under salt and drought stress conditions by using qRT-PCR technique. Our results showed that about one-third of these B. rapa C2H2-ZFP genes were originated from segmental duplication caused by the WGT around 13 to 17 MYA, one-third of them were highly and consecutively expressed in all tested tissues, and 92% of them were located in nucleus by prediction supporting then their functional roles as transcription factors, of which some may play important roles in plant growth and development. The Ka/Ks ratios of 264 orthologous C2H2-ZFP gene pairs between A. thaliana and B. rapa were all, except two, inferior to 1 (varied from 0.0116 to 1.4919, with an average value of 0.3082), implying that these genes had mainly experienced purifying selection during species evolution. The estimated divergence times of the same set of gene pairs ranged from 6.23 to 38.60 MY, with an average value of 18.29 MY, indicating that these gene members have undergone different selective pressures resulting in different evolutionary rates during species evolution. In addition, a few of these B. rapa C2H2-ZFPs were shown to be involved in stress responses in a similar way as their orthologs in A. thaliana. Comparison between A. thaliana and B. rapa orthologous C2H2-ZFP genes showed that the majority of these C2H2-ZFP gene members encodes proteins with conserved subcellular localization and functional domains between the two species but differed in their expression patterns in five tissues or organs. Thus, our study provides valuable information for further functional determination of each C2H2-ZFP gene across the Brassica species, and may help to select the appropriate gene targets for further in-depth studies, and genetic engineering and improvement of Brassica crops.

Comprehensive Genomic Survey, Characterization and Expression Analysis of the HECT Gene Family in Brassica rapa L. and Brassica oleracea L.

The HECT-domain protein family is one of the most important classes of E3 ligases. While the roles of this family in human diseases have been intensively studied, the information for plant HECTs is limited. In the present study, we performed the identification of HECT genes in Brassica rapa and Brassica oleracea, followed by analysis of phylogeny, gene structure, additional domains, putative cis-regulatory elements, chromosomal location, synteny, and expression. Ten and 13 HECT genes were respectively identified in B. rapa and B. oleracea and then resolved into seven groups along with their Arabidopsis orthologs by phylogenetic analysis. This classification is well supported by analyses of gene structure, motif composition within the HECT domain and additional protein domains. Ka/Ks ratio analysis showed that these HECT genes primarily underwent purifying selection with varied selection pressures resulting in different rates of evolution. RNA-Seq data analysis showed that the overwhelming majority of them were constitutively expressed in all tested tissues. qRT-PCR based expression analysis of the 10 B. rapa HECT genes under salt and drought stress conditions showed that all of them were responsive to the two stress treatments, which was consistent with their promoter sequence analysis revealing the presence of an important number of phytohormone-responsive and stress-related cis-regulatory elements. Our study provides useful information and lays the foundation for further functional determination of each HECT gene across the Brassica species.

Genome-wide survey, characterization, and expression analysis of RING finger protein genes in Brassica oleracea and their syntenic comparison to Brassica rapa and Arabidopsis thaliana.

The ubiquitin-mediated post-translational regulatory pathway regulates a broad range of cell functions in all eukaryotes. It requires the involvement of a large number of E3 ligases, of which more than one third belongs to the RING protein family as in Arabidopsis thaliana. In this study, a total of 756 RING domains in 734 predicted proteins were identified in Brassica oleracea. Their encoding genes were characterized by RING domain type, additional domain, and expression pattern, and mapped on the nine chromosomes of B. oleracea. Comparison of these results with B. rapa and A. thaliana revealed some common as well as species-specific features. Our results showed that the differential gene loss following the whole genome triplication has largely contributed to the RING protein gene number variation among these species, although other factors such as tandem duplication, RING domain loss, or modification had also contributed to this variation. Analysis of RNA-seq data showed that these RING protein genes were functionally diversified and involved in all the stages of plant growth and development, and that the triplicated members were also diverged in expression with one member often more dominantly expressed over the two others in the majority of cases. Our study lays the foundation for further functional determination of each RING protein gene among species of the genus Brassica.

Molecular insights into the origin of the brown rust resistance gene Bru1 among Saccharum species.

KEY MESSAGE: Analysis of 387 sugarcane clones using Bru 1 diagnostic markers revealed two possible sources of Bru 1 in Chinese cultivars: one from Saccharum spontaneum and another from Saccharum robustum of New Guinea. Sugarcane brown rust (SBR) is an important fungal disease in many sugarcane production areas around the world, and can cause considerable yield losses in susceptible sugarcane cultivars. One major SBR resistance gene, named Bru1, initially identified from cultivar R570, was shown to be a major SBR resistance source in most of the sugarcane producing areas of the world. In this study, by using the two Bru1-associated markers, R12H16 and 9O20-F4, we surveyed the presence of Bru1 in a Chinese sugarcane germplasm collection of 387 clones, consisting of 228 hybrid cultivars bred by different Chinese sugarcane breeding establishments, 54 exotic hybrid cultivars introduced from other countries and 105 clones of sugarcane ancestral species. The Bru1-bearing haplotype was detected in 43.4% of Chinese sugarcane cultivars, 20.4% of exotic hybrid cultivars, and only 3.8% of ancestral species. Among the 33 Chinese cultivars for which phenotypes of resistance to SBR were available, Bru1 was present in 69.2% (18/26) of the resistant clones. Analyses of the allelic sequence variations of R12H16 and 9O20-F4 suggested two possible sources of Bru1 in Chinese cultivars: one from S. spontaneum and another from S. robustum of New Guinea. In addition, we developed an improved Bru1 diagnostic marker, 9O20-F4-HaeIII, which can eliminate all the false results of 9O20-F4-RsaI observed among S. spontaneum, as well as a new dominant Bru1 diagnostic marker, R12E03-2, from the BAC ShCIR12E03. Our results provide valuable information for further efforts of breeding SBR-resistant varieties, searching new SBR resistance sources and cloning of Bru1 in sugarcane.

Genome-wide identification, evolution and expression analysis of RING finger protein genes in Brassica rapa.

More and more RING finger genes were found to be implicated in various important biological processes. In the present study, a total of 731 RING domains in 715 predicted proteins were identified in Brassica rapa genome (AA, 2n = 20), which were further divided into eight types: RING-H2 (371), RING-HCa (215), RING-HCb (47), RING-v (44), RING-C2 (38), RING-D (10), RING-S/T (5) and RING-G (1). The 715 RING finger proteins were further classified into 51 groups according to the presence of additional domains. 700 RING finger protein genes were mapped to the 10 chromosomes of B. rapa with a range of 47 to 111 genes for each chromosome. 667 RING finger protein genes were expressed in at least one of the six tissues examined, indicating their involvement in various physiological and developmental processes in B. rapa. Hierarchical clustering analysis of RNA-seq data divided them into seven major groups, one of which includes 231 members preferentially expressed in leaf, and constitutes then a panel of gene candidates for studying the genetic and molecular mechanisms of leafy head traits in Brassica crops. Our results lay the foundation for further studies on the classification, evolution and putative functions of RING finger protein genes in Brassica species.

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome.

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Biogeographical variation and population genetic structure of Sporisorium scitamineum in Mainland China: insights from ISSR and SP-SRAP markers.

A total of 100 Sporisorium scitamineum isolates were investigated by inter simple sequence repeat (ISSR) and single primer-sequence related amplified polymorphism (SP-SRAP) markers. These isolates were clearly assorted into three distinct clusters regardless of method used: either cluster analysis or by principal component analysis (PCA) of the ISSR, SP-SRAP, or ISSR + SP-SRAP data set. The total gene diversity (H t) and gene diversity between subpopulations (H s) were estimated to be 0.34 to 0.38 and 0.22 to 0.29, respectively, by analyzing separately the ISSR and SP-SRAP data sets, and to be 0.26-0.36 by analyzing ISSR + SP-SRAP data set. The gene diversity attributable to differentiation among populations (G st) was estimated to be 0.35 and 0.22, and the gene flow (Nm) was 0.94 and 1.78, respectively, when analyzing separately ISSR and SP-SRAP data set, and was 0.27 and 1.33, respectively, when analyzing ISSR + SP-SRAP data set. Our study showed that there is considerable genetic variation in the analyzed 100 isolates, and the environmental heterogeneity has played an important role for this observed high degree of variation. The genetic differentiation of sugarcane smut fungus depends to a large extent on the heterogeneity of their habitats and is the result of long-term adaptations of pathogens to their ecological environments.

Genetic analysis of diversity within a Chinese local sugarcane germplasm based on start codon targeted polymorphism.

In-depth information on sugarcane germplasm is the basis for its conservation and utilization. Data on sugarcane molecular markers are limited for the Chinese sugarcane germplasm collections. In the present study, 20 start codon targeted (SCoT) marker primers were designed to assess the genetic diversity among 107 sugarcane accessions within a local sugarcane germplasm collection. These primers amplified 176 DNA fragments, of which 163 were polymorphic (92.85%). Polymorphic information content (PIC) values ranged from 0.783 to 0.907 with a mean of 0.861. Unweighted pair group method of arithmetic averages (UPGMA) cluster analysis of the SCoT marker data divided the 107 sugarcane accessions into six clusters at 0.674 genetic similarity coefficient level. Relatively abundant genetic diversity was observed among ROC22, ROC16, and ROC10, which occupied about 80% of the total sugarcane acreage in China, indicating their potential breeding value on Mainland China. Principal component analysis (PCA) partitioned the 107 sugarcane accessions into two major groups, the Domestic Group and the Foreign Introduction Group. Each group was further divided based on institutions, where the sugarcane accessions were originally developed. The knowledge of genetic diversity among the local sugarcane germplasm provided foundation data for managing sugarcane germplasm, including construction of a core collection and regional variety distribution and subrogation.

Prevalence of gene expression additivity in genetically stable wheat allohexaploids.

The reprogramming of gene expression appears as the major trend in synthetic and natural allopolyploids where expression of an important proportion of genes was shown to deviate from that of the parents or the average of the parents. In this study, we analyzed gene expression changes in previously reported, highly stable synthetic wheat allohexaploids that combine the D genome of Aegilops tauschii and the AB genome extracted from the natural hexaploid wheat Triticum aestivum. A comprehensive genome-wide analysis of transcriptional changes using the Affymetrix GeneChip Wheat Genome Array was conducted. Prevalence of gene expression additivity was observed where expression does not deviate from the average of the parents for 99.3% of 34,820 expressed transcripts. Moreover, nearly similar expression was observed (for 99.5% of genes) when comparing these synthetic and natural wheat allohexaploids. Such near-complete additivity has never been reported for other allopolyploids and, more interestingly, for other synthetic wheat allohexaploids that differ from the ones studied here by having the natural tetraploid Triticum turgidum as the AB genome progenitor. Our study gave insights into the dynamics of additive gene expression in the highly stable wheat allohexaploids.

A dominant point mutation in a RINGv E3 ubiquitin ligase homoeologous gene leads to cleistogamy in Brassica napus.

In the allopolyploid Brassica napus, we obtained a petal-closed flower mutation by ethyl methanesulfonate mutagenesis. Here, we report cloning and characterization of the Bn-CLG1A (CLG for cleistogamy) gene and the Bn-clg1A-1D mutant allele responsible for the cleistogamy phenotype. Bn-CLG1A encodes a RINGv E3 ubiquitin ligase that is highly conserved across eukaryotes. In the Bn-clg1A-1D mutant allele, a C-to-T transition converts a Pro at position 325 to a Leu (P325L), causing a dominant mutation leading to cleistogamy. B. napus and Arabidopsis thaliana plants transformed with a Bn-clg1A-1D allele show cleistogamous flowers, and characterization of these flowers suggests that the Bn-clg1A-1D mutation causes a pronounced negative regulation of cutin biosynthesis or loading and affects elongation or differentiation of petal and sepal cells. This results in an inhibition or a delay of petal development, leading to folded petals. A homoeologous gene (Bn-CLG1C), which shows 99.5% amino acid identity and is also constitutively and equally expressed to the wild-type Bn-CLG1A gene, was also identified. We showed that P325L is not a loss-of-function mutation and did not affect expression of Bn-clg1A-1D or Bn-CLG1C. Our findings suggest that P325L is a gain-of-function semidominant mutation, which led to either hyper- or neofunctionalization of a redundant homoeologous gene.