Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae

Genome assembly and multiomic analyses reveal insights into flower and bark colors of Lagerstroemia excelsa

Lagerstroemia excelsa is a unique plant species from China, holds a significant aesthetic and economic value, and plays a crucial role in landscape architecture and horticulture. Thus far, there is little genetic and genomic information available about this species, which limits its use in development of new cultivars. In this study, a high-quality genome map of L. excelsa was obtained via whole-genome sequencing. Results showed that its genome size is about 330.4 Mb and a scaffold mapping rate is approximately 97.20%, resulting in 24 pseudochromosomes. L. excelsa might have undergone a recent whole-genome triplication event and diverged from the pomegranate about 32.3 million years ago (MYA). Subsequently, the divergence time between L. indica and L. excelsa was around 5.9 MYA. The transcriptomic and metabolomic analyses of L. excelsa and L. indica indicated that the chalcone synthase pathway may play a key role in regulating flower color differentiation between the two species. Additionally, a transcription factor LeMYB103 may be involved in regulating anthocyanin synthesis by interacting with LeMYB66, resulting in the accumulation of anthocyanins in the stem bark. This study is the first step toward genomic analysis of L. excelsa, which may provide a foundation for further molecular investigation of this species and offer valuable insights into the molecular mechanisms underlying the flower and stem bark colors in L. excelsa, two important ornamental traits in Lagerstroemia breeding.

Different orthology inference algorithms generate similar predicted orthogroups among Brassicaceae species

Premise

Orthology inference is crucial for comparative genomics, and multiple algorithms have been developed to identify putative orthologs for downstream analyses. Despite the abundance of proposed solutions, including publicly available benchmarks, it is difficult to assess which tool is most suitable for plant species, which commonly have complex genomic histories.

Methods

We explored the performance of four orthology inference algorithms-OrthoFinder, SonicParanoid, Broccoli, and OrthNet-on eight Brassicaceae genomes in two groups: one group comprising only diploids and another set comprising the diploids, two mesopolyploids, and one recent hexaploid genome.

Results

The composition of the orthogroups reflected the species' ploidy and genomic histories, with the diploid set having a higher proportion of identical orthogroups. While the diploid + higher ploidy set had a lower proportion of orthogroups with identical compositions, the average degree of similarity between the orthogroups was not different from the diploid set.

Discussion

Three algorithms-OrthoFinder, SonicParanoid, and Broccoli-are helpful for initial orthology predictions. Results produced using OrthNet were generally outliers but could still provide detailed information about gene colinearity. With our Brassicaceae dataset, slight discrepancies were found across the orthology inference algorithms, necessitating additional analyses such as tree inference to fine-tune results.

Heat shock responsive genes in Brassicaceae: genome-wide identification, phylogeny, and evolutionary associations within and between genera

Heat stress affects the growth and development of Brassicaceae crops. Plant breeders aim to mitigate the effects of heat stress by selecting for heat stress tolerance, but the genes responsible for heat stress in Brassicaceae remain largely unknown. During heat stress, heat shock proteins (HSPs) function as molecular chaperones to aid in protein folding, and heat shock transcription factors (HSFs) serve as transcriptional regulators of HSP expression. We identified 5002 heat shock related genes, including HSPs and HSFs, across 32 genomes in Brassicaceae. Among these, 3347 genes were duplicated, with segmented duplication primarily contributing to their expansion. We identified 466 physical gene clusters, including 240 homogenous clusters and 226 heterogeneous clusters, shedding light on the organization of heat shock related genes. Notably, 37 genes were co-located with published thermotolerance quantitative trait loci, which supports their functional role in conferring heat stress tolerance. This study provides a comprehensive resource for the identification of functional Brassicaceae heat shock related genes, elucidates their clustering and duplication patterns and establishes the genomic foundation for future heat tolerance research. We hypothesise that genetic variants in HSP and HSF genes in certain species have potential for improving heat stress tolerance in Brassicaceae crops.

Expanding the Triangle of U: Comparative analysis of the Hirschfeldia incana genome provides insights into chromosomal evolution, phylogenomics and high photosynthesis-related traits

BACKGROUND AND AIMS

The Brassiceae tribe encompasses many economically important crops and exhibits high intraspecific and interspecific phenotypic variation. After a shared whole-genome triplication (WGT) event (Br-α, ~15.9 million years ago), differential lineage diversification and genomic changes contributed to an array of divergence in morphology, biochemistry, and physiology underlying photosynthesis-related traits. Here, the C3 species Hirschfeldia incana is studied as it displays high photosynthetic rates under high-light conditions. Our aim was to elucidate the evolution that gave rise to the genome of H. incana and its high-photosynthesis traits.

METHODS

We reconstructed a chromosome-level genome assembly for H. incana (Nijmegen, v2.0) using nanopore and chromosome conformation capture (Hi-C) technologies, with 409Mb in size and an N50 of 52Mb (a 10× improvement over the previously published scaffold-level v1.0 assembly). The updated assembly and annotation was subsequently employed to investigate the WGT history of H. incana in a comparative phylogenomic framework from the Brassiceae ancestral genomic blocks and related diploidized crops.

KEY RESULTS

Hirschfeldia incana (x=7) shares extensive genome collinearity with Raphanus sativus (x=9). These two species share some commonalities with Brassica rapa and B. oleracea (A genome, x=10 and C genome, x=9, respectively) and other similarities with B. nigra (B genome, x=8). Phylogenetic analysis revealed that H. incana and R. sativus form a monophyletic clade in between the Brassica A/C and B genomes. We postulate that H. incana and R. sativus genomes are results of hybridization or introgression of the Brassica A/C and B genome types. Our results might explain the discrepancy observed in published studies regarding phylogenetic placement of H. incana and R. sativus in relation to the "Triangle of U" species. Expression analysis of WGT retained gene copies revealed sub-genome expression divergence, likely due to neo- or sub-functionalization. Finally, we highlighted genes associated with physio-biochemical-anatomical adaptive changes observed in H. incana which likely facilitate its high-photosynthesis traits under high light.

CONCLUSIONS

The improved H. incana genome assembly, annotation and results presented in this work will be a valuable resource for future research to unravel the genetic basis of its ability to maintain a high photosynthetic efficiency in high-light conditions and thereby improve photosynthesis for enhanced agricultural production.

Molecular evolution analysis of MYB5 in Brassicaceae with specific focus on seed coat color of Brassica napus

BACKGROUND

MYB transcription factors are splay a vital role in plant biology, with previous research highlighting the significant impact of the R2R3-MYB-like transcription factor MYB5 on seed mucilage biosynthesis, trichome branching, and seed coat development. However, there is a dearth of studies investigating its role in the regulation of proanthocyanidin (PA) biosynthesis.

RESULTS

In this study, a total of 51 MYB5 homologous genes were identified across 31 species belonging to the Brassicaceae family, with particular emphasis on Brassica napus for subsequent investigation. Through phylogenetic analysis, these genes were categorized into four distinct subclasses. Protein sequence similarity and identity analysis demonstrated a high degree of conservation of MYB5 among species within the Brassicaceae family. Additionally, the examination of selection pressure revealed that MYB5 predominantly underwent purifying selection during its evolutionary history, as indicated by the Ka/Ks values of all MYB5 homologous gene pairs being less than one. Notably, we observed a higher rate of non-synonymous mutations in orthologous genes compared to paralogous genes, and the Ka/Ks value displayed a stronger correlation with Ka. In B. napus, an examination of expression patterns in five tissues revealed that MYB5 exhibited particularly high expression in the black seed coat. The findings from the WGCNA demonstrated a robust correlation between MYB5 and BAN(ANR) associated with PA biosynthesis in the black seed coat, providing further evidence of their close association and co-expression. Furthermore, the results obtained from of the analysis of protein interaction networks offer supplementary support for the proposition that MYB5 possesses the capability to interact with transcriptional regulatory proteins, specifically TT8 and TT2, alongside catalytic enzymes implicated in the synthesis of PAs, thereby making a contribution to the biosynthesis of PAs. These findings imply a plausible and significant correlation between the nuique expression pattern of MYB5 and the pigmentation of rapeseed coats. Nevertheless, additional research endeavors are imperative to authenticate and substantiate these findings.

CONCLUSIONS

This study offers valuable insights into the genetic evolution of Brassicaceae plants, thereby serving as a significant reference for the genetic enhancement of Brassicaceae seed coat color.

Integrative Physiological and Transcriptome Analysis Reveals the Mechanism of Cd Tolerance in Sinapis alba

Recently, pollution caused by the heavy metal Cd has seriously affected the environment and agricultural crops. While Sinapis alba is known for its edible and medicinal value, its tolerance to Cd and molecular response mechanism remain unknown. This study aimed to analyze the tolerance of S. alba to Cd and investigate its molecular response mechanism through transcriptomic and physiological indicators. To achieve this, S. alba seedlings were treated with different concentrations of CdCl2 (0.25 mmol/L, 0.5 mmol/L, and 1.0 mmol/L) for three days. Based on seedling performance, S. alba exhibited some tolerance to a low concentration of Cd stress (0.25 mmol/L CdCl2) and a strong Cd accumulation ability in its roots. The activities and contents of several antioxidant enzymes generally exhibited an increase under the treatment of 0.25 mmol/L CdCl2 but decreased under the treatment of higher CdCl2 concentrations. In particular, the proline (Pro) content was extremely elevated under the 0.25 and 0.5 mmol/L CdCl2 treatments but sharply declined under the 1.0 mmol/L CdCl2 treatment, suggesting that Pro is involved in the tolerance of S. alba to low concentration of Cd stress. In addition, RNA sequencing was utilized to analyze the gene expression profiles of S. alba exposed to Cd (under the treatment of 0.25 mmol/L CdCl2). The results indicate that roots were more susceptible to disturbance from Cd stress, as evidenced by the detection of 542 differentially expressed genes (DEGs) in roots compared to only 37 DEGs in leaves. GO and KEGG analyses found that the DEGs induced by Cd stress were primarily enriched in metabolic pathways, plant hormone signal transduction, and the biosynthesis of secondary metabolites. The key pathway hub genes were mainly associated with intracellular ion transport and cell wall synthesis. These findings suggest that S. alba is tolerant to a degree of Cd stress, but is also susceptible to the toxic effects of Cd. Furthermore, these results provide a theoretical basis for understanding Cd tolerance in S. alba.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Evolutionary Dynamics of FLC-like MADS-Box Genes in Brassicaceae

MADS-box genes encode transcription factors that play important roles in the development and evolution of plants. There are more than a dozen clades of MADS-box genes in angiosperms, of which those with functions in the specification of floral organ identity are especially well-known. From what has been elucidated in the model plant Arabidopsis thaliana, the clade of FLC-like MADS-box genes, comprising FLC-like genes sensu strictu and MAF-like genes, are somewhat special among the MADS-box genes of plants since FLC-like genes, especially MAF-like genes, show unusual evolutionary dynamics, in that they generate clusters of tandemly duplicated genes. Here, we make use of the latest genomic data of Brassicaceae to study this remarkable feature of the FLC-like genes in a phylogenetic context. We have identified all FLC-like genes in the genomes of 29 species of Brassicaceae and reconstructed the phylogeny of these genes employing a Maximum Likelihood method. In addition, we conducted selection analyses using PAML. Our results reveal that there are three major clades of FLC-like genes in Brassicaceae that all evolve under purifying selection but with remarkably different strengths. We confirm that the tandem arrangement of MAF-like genes in the genomes of Brassicaceae resulted in a high rate of duplications and losses. Interestingly, MAF-like genes also seem to be prone to transposition. Considering the role of FLC-like genes sensu lato (s.l.) in the timing of floral transition, we hypothesize that this rapid evolution of the MAF-like genes was a main contributor to the successful adaptation of Brassicaceae to different environments.

Sequencing and Analysis of Complete Chloroplast Genomes Provide Insight into the Evolution and Phylogeny of Chinese Kale (Brassica oleracea var. alboglabra)

Chinese kale is a widely cultivated plant in the genus Brassica in the family Brassicaceae. The origin of Brassica has been studied extensively, but the origin of Chinese kale remains unclear. In contrast to Brassica oleracea, which originated in the Mediterranean region, Chinese kale originated in southern China. The chloroplast genome is often used for phylogenetic analysis because of its high conservatism. Fifteen pairs of universal primers were used to amplify the chloroplast genomes of white-flower Chinese kale (Brassica oleracea var. alboglabra cv. Sijicutiao (SJCT)) and yellow-flower Chinese kale (Brassica oleracea var. alboglabra cv. Fuzhouhuanghua (FZHH)) via PCR. The lengths of the chloroplast genomes were 153,365 bp (SJCT) and 153,420 bp (FZHH) and both contained 87 protein-coding genes and eight rRNA genes. There were 36 tRNA genes in SJCT and 35 tRNA genes in FZHH. The chloroplast genomes of both Chinese kale varieties, along with eight other Brassicaceae, were analyzed. Simple sequence repeats, long repeats, and variable regions of DNA barcodes were identified. An analysis of inverted repeat boundaries, relative synonymous codon usage, and synteny revealed high similarity among the ten species, albeit the slight differences that were observed. The Ka/Ks ratios and phylogenetic analysis suggest that Chinese kale is a variant of B. oleracea. The phylogenetic tree shows that both Chinese kale varieties and B. oleracea var. oleracea were clustered in a single group. The results of this study suggest that white and yellow flower Chinese kale comprise a monophyletic group and that their differences in flower color arose late in the process of artificial cultivation. Our results also provide data that will aid future research on genetics, evolution, and germplasm resources of Brassicaceae.