Genome-wide identification and expression analysis of chitinase gene family in Brassica rapa reveals its role in clubroot resistance

Genome-Wide Mining of Chitinase Diversity in the Marine Diatom Thalassiosira weissflogii and Functional Characterization of a Novel GH19 Enzyme

Chitin represents a globally abundant marine polymer with significant ecological and biotechnological value. β-chitin is an important carbon fixation product of diatoms and has a greater range of applications than α- and γ-chitin. However, there has been a paucity of research on the characterization of chitin-related enzymes from β-chitin producers. In this study, we performed a genome-wide identification of 38 putative chitinase genes in Thalassiosira weissflogii, a key producer of β-chitin. Through comprehensive analyses of phylogenetic relationships, conserved motifs, structural domains, and subcellular localization predictions, we revealed that T. weissflogii possesses evolutionarily distinct GH18 and GH19 chitinase families exhibiting unique motif and domain configurations. Subcellular localization predictions showed that most TwChis were presumed to be located in the chloroplast, with a few being present in the nucleus and extracellular. The enzymatic activity of TwChi2, a GH19 chitinase, showed that TwChi2 was a member of exochitinase (EC 3.2.1.201) with strong thermal stability (40 °C) and broad substrate adaptability of hydrolyzing bipolymer, 1% and 5% colloidal chitin, α-chitin and β-chitin. Altogether, we analyzed the chitinase gene family and characterized a highly active exochitinase from T. weissflogii, which can catalyze the degradation of both chitin polymers and chitin oligosaccharides. The relevant results lay a foundation for the internal regulation mechanism of chitin metabolism in diatoms and provide a candidate enzyme for the green industrial preparation of high-value chitin oligosaccharides.

Plant microbiota feedbacks through dose-responsive expression of general non-self response genes

The ability of plants to perceive and react to biotic and abiotic stresses is critical for their health. We recently identified a core set of genes consistently induced by members of the leaf microbiota, termed general non-self response (GNSR) genes. Here we show that GNSR components conversely impact leaf microbiota composition. Specific strains that benefited from this altered assembly triggered strong plant responses, suggesting that the GNSR is a dynamic system that modulates colonization by certain strains. Examination of the GNSR to live and inactivated bacteria revealed that bacterial abundance, cellular composition and exposure time collectively determine the extent of the host response. We link the GNSR to pattern-triggered immunity, as diverse microbe- or danger-associated molecular patterns cause dynamic GNSR gene expression. Our findings suggest that the GNSR is the result of a dose-responsive perception and signalling system that feeds back to the leaf microbiota and contributes to the intricate balance of plant-microbiome interactions.

Genome-wide analysis of the NYN domain gene family in Brassica napus and its function role in plant growth and development

The NYN domain gene family consists of genes that encode ribonucleases that are characterized by a newly identified NYN domain. Members of the family were widely distributed in all life kingdoms and play a crucial role in various RNA regulation processes, although the wide genome overview of the NYN domain gene family is not yet available in any species. Rapeseed (Brassica napus L.), a polyploid model species, is an important oilseed crop. Here, the phylogenetic analysis of these BnaNYNs revealed five distinct groups strongly supported by gene structure, conserved domains, and conserved motifs. The survey of the expansion of the gene family showed that the birth of BnaNYNs is explained by various duplication events. Furthermore, tissue-specific expression analysis, protein-protein interaction prediction, and cis-element prediction suggested a role for BnaNYNs in plant growth and development. Interestingly, the data showed that three tandem duplicated BnaNYNs (TDBs) exhibited distinct expression patterns from those other BnaNYNs and had a high similarity in protein sequence level. Furthermore, the analysis of one of these TDBs, BnaNYN57, showed that overexpression of BnaNYN57 in Arabidopsis thaliana and B. napus accelerated plant growth and significantly increased silique length, while RNA interference resulted in the opposite growth pattern. It suggesting a key role for the TDBs in processes related to plant growth and development.

Systematic analysis and functional characterization of the chitinase gene family in Fagopyrum tataricum under salt stress

BACKGROUND

Chitinases (CHIs) are glycosidases that degrade chitin, playing critical roles in plant responses to both abiotic and biotic stress. Despite their importance, the CHI family's systematic analysis and evolutionary pattern in F. tataricum (Tartary buckwheat) yet to be explored.

RESULTS

This study analyzed their phylogenetic relationships, conserved motifs, gene structures, syntenic relationships, physiological functions, and biochemical properties. This research identified 26 FtCHIs and examined their expression patterns under different salt stress conditions and across various tissues. Differential expression analysis revealed a significant upregulation of multiple FtCHIs in response to salt stress, which RT-qPCR further validated. Additionally, subcellular localization experiments demonstrated that Ft_chitinaseIV-2 is localized in vacuoles. The results of transient·transformation showed that·overexpression of Ft_chitinaseIV-2 could·enhance the salt tolerance of plants.

CONCLUSIONS

The findings provide new insights into the role of CHIs in stress tolerance and lay the groundwork for future research on the functional characterization of FtCHIs. Understanding the molecular mechanisms of CHI-mediated stress responses could contribute to developing stress-resistant crops.

Characterization of the chitinase gene family in Saccharum reveals the disease resistance mechanism of ScChiVII1

KEY MESSAGE

A chitinase gene ScChiVII1 which is involved in defense against pathogen stress was characterized in sugarcane. Chitinases, a subclass of pathogenesis-related proteins, catalyze chitin hydrolysis and play a key role in plant defense against chitin-containing pathogens. However, there is little research on disease resistance analysis of chitinase genes in sugarcane, and the systematic identification of their gene families has not been reported. In this study, 85 SsChi and 23 ShChi genes, which were divided into 6 groups, were identified from the wild sugarcane species Saccharum spontaneum and Saccharum hybrid cultivar R570, respectively. Transcriptome analysis and real-time quantitative PCR revealed that SsChi genes responded to smut pathogen stress. The chitinase crude extracted from the leaves of transgenic Nicotiana benthamiana plants overexpressing ScChiVII1 (a homologous gene of SsChi22a) inhibited the hyphal growth of Fusarium solani var. coeruleum and Sporisorium scitamineum. Notably, the chitinase and catalase activities and the jasmonic acid content in the leaves of ScChiVII1 transgenic N. benthamiana increased after inoculation with F solani var. coeruleum, but the salicylic acid, hydrogen peroxide, and malondialdehyde contents decreased. Comprehensive RNA sequencing of leaves before (0 day) and after inoculation (2 days) revealed that ScChiVII1 transgenic tobacco enhanced plant disease resistance by activating transcription factors and disease resistance-related signaling pathways, and modulating the expression of genes involved in the hypersensitive response and ethylene synthesis pathways. Taken together, this study provides comprehensive information on the chitinase gene family and offers potential genetic resources for disease resistance breeding in sugarcane.

Transcriptome Analysis of Chinese Cabbage Infected with Plasmodiophora Brassicae in the Primary Stage

Clubroot disease caused by the infection of Plasmodiophora brassicae is widespread in China, and significantly reduces the yield of Chinese cabbage (Brassica rapa L. ssp. pekinensis). However, the resistance mechanism of Chinese cabbage against clubroot disease is still unclear. It is important to exploit the key genes that response to early infection of P. brassicae. In this study, it was found that zoospores were firstly invaded hair roots on the 8th day after inoculating with 1 × 107 spores/mL P. brassicae. Transcriptome analysis found that the early interaction between Chinese cabbage and P. brassicae caused the significant expression change of some defense genes, such as NBS-LRRs and pathogenesis-related genes, etc. The above results were verified by quantitative reverse-transcription polymerase chain reaction (qRT-PCR). Otherwise, peroxidase (POD) salicylic acid (SA) and jasmonic acid (JA) were also found to be important signal molecules in the resistance to clubroot disease in Chinese cabbage. This study provides important clues for understanding the resistance mechanism of Chinese cabbage against clubroot disease.

Bioinformatics and functional analysis of EDS1 genes in Brassica napus in response to Plasmodiophora brassicae infection

Enhanced Disease Susceptibility 1 (EDS1) is a key regulator of plant-pathogen-associated molecular pattern-triggered immunity (PTI) and effector-triggered immunity (ETI) responses. In the Brassica napus genome, we identified six novel EDS1 genes, among which four were responsive to clubroot infection, a major rapeseed disease resistant to chemical control. Developing resistant cultivars is a potent and economically viable strategy to control clubroot infection. Bioinformatics analysis revealed conserved domains and structural uniformity in Bna-EDS1 homologs. Bna-EDS1 promoters harbored elements associated with diverse phytohormones and stress responses, highlighting their crucial roles in plant defense. A functional analysis was performed with Bna-EDS1 overexpression and RNAi transgenic lines. Bna-EDS1 overexpression boosted resistance to clubroot and upregulated defense-associated genes (PR1, PR2, ICS1, and CBP60), while Bna-EDS1 RNAi increased plant susceptibility, indicating suppression of the defense signaling pathway downstream of NBS-LRRs. RNA-Seq analysis identified key transcripts associated with clubroot resistance, including phenylpropanoid biosynthesis. Activation of SA regulator NPR1, defense signaling markers PR1 and PR2, and upregulation of MYC-TFs suggested that EDS1-mediated clubroot resistance potentially involves the SA pathway. Our findings underscore the pivotal role of Bna-EDS1-dependent mechanisms in resistance of B. napus to clubroot disease, and provide valuable insights for fortifying resistance against Plasmodiophora brassicae infection in rapeseed.

Comparative Analysis of Transcriptomes Reveals Pathways and Verifies Candidate Genes for Clubroot Resistance in Brassica oleracea

Clubroot, a soil-borne disease caused by Plasmodiophora brassicae, is one of the most destructive diseases of Brassica oleracea all over the world. However, the mechanism of clubroot resistance remains unclear. In this research, transcriptome sequencing was conducted on root samples from both resistant (R) and susceptible (S) B. oleracea plants infected by P. brassicae. Then the comparative analysis was carried out between the R and S samples at different time points during the infection stages to reveal clubroot resistance related pathways and candidate genes. Compared with 0 days after inoculation, a total of 4991 differential expressed genes were detected from the S pool, while only 2133 were found from the R pool. Gene function enrichment analysis found that the effector-triggered immunity played a major role in the R pool, while the pathogen-associated molecular pattern triggered immune response was stronger in the S pool. Simultaneously, candidate genes were identified through weighted gene co-expression network analysis, with Bol010786 (CNGC13) and Bol017921 (SD2-5) showing potential for conferring resistance to clubroot. The findings of this research provide valuable insights into the molecular mechanisms underlying clubroot resistance and present new avenues for further research aimed at enhancing the clubroot resistance of B. oleracea through breeding.

Genome-Wide Identification, Expression, and Protein Analysis of CKX and IPT Gene Families in Radish (Raphanus sativus L.) Reveal Their Involvement in Clubroot Resistance

Cytokinins (CKs) are a group of phytohormones that are involved in plant growth, development, and disease resistance. The isopentenyl transferase (IPT) and cytokinin oxidase/dehydrogenase (CKX) families comprise key enzymes controlling CK biosynthesis and degradation. However, an integrated analysis of these two gene families in radish has not yet been explored. In this study, 13 RsIPT and 12 RsCKX genes were identified and characterized, most of which had four copies in Brassica napus and two copies in radish and other diploid Brassica species. Promoter analysis indicated that the genes contained at least one phytohormone or defense and stress responsiveness cis-acting element. RsIPTs and RsCKXs were expanded through segmental duplication. Moreover, strong purifying selection drove the evolution of the two gene families. The expression of the RsIPT and RsCKX genes distinctly showed diversity in different tissues and developmental stages of the root. Expression profiling showed that RsCKX1-1/1-2/1-3 was significantly upregulated in club-resistant materials during primary infection, suggesting their vital function in clubroot resistance. The interaction network of CKX proteins with similar 3D structures also reflected the important role of RsCKX genes in disease resistance. This study provides a foundation for further functional study on the IPT and CKX genes for clubroot resistance improvement in Raphanus.

Characterization of soybean chitinase genes induced by rhizobacteria involved in the defense against Fusarium oxysporum

Rhizobacteria are capable of inducing defense responses via the expression of pathogenesis-related proteins (PR-proteins) such as chitinases, and many studies have validated the functions of plant chitinases in defense responses. Soybean (Glycine max) is an economically important crop worldwide, but the functional validation of soybean chitinase in defense responses remains limited. In this study, genome-wide characterization of soybean chitinases was conducted, and the defense contribution of three chitinases (GmChi01, GmChi02, or GmChi16) was validated in Arabidopsis transgenic lines against the soil-borne pathogen Fusarium oxysporum. Compared to the Arabidopsis Col-0 and empty vector controls, the transgenic lines with GmChi02 or GmChi16 exhibited fewer chlorosis symptoms and wilting. While GmChi02 and GmChi16 enhanced defense to F. oxysporum, GmChi02 was the only one significantly induced by Burkholderia ambifaria. The observation indicated that plant chitinases may be induced by different rhizobacteria for defense responses. The survey of 37 soybean chitinase gene expressions in response to six rhizobacteria observed diverse inducibility, where only 10 genes were significantly upregulated by at least one rhizobacterium and 9 genes did not respond to any of the rhizobacteria. Motif analysis on soybean promoters further identified not only consensus but also rhizobacterium-specific transcription factor-binding sites for the inducible chitinase genes. Collectively, these results confirmed the involvement of GmChi02 and GmChi16 in defense enhancement and highlighted the diverse inducibility of 37 soybean chitinases encountering F. oxysporum and six rhizobacteria.

Genome-Wide Identification and Expression Analysis of Chitinase Genes in Watermelon under Abiotic Stimuli and Fusarium oxysporum Infection

Chitinases, which catalyze the hydrolysis of chitin, the primary components of fungal cell walls, play key roles in defense responses, symbiotic associations, plant growth, and stress tolerance. In this study, 23 chitinase genes were identified in watermelon (Citrullus lanatus [Thunb.]) and classified into five classes through homology search and phylogenetic analysis. The genes with similar exon-intron structures and conserved domains were clustered into the same class. The putative cis-elements involved in the responses to phytohormone, stress, and plant development were identified in their promoter regions. A tissue-specific expression analysis showed that the ClChi genes were primarily expressed in the roots (52.17%), leaves (26.09%), and flowers (34.78%). Moreover, qRT-PCR results indicate that ClChis play multifaceted roles in the interaction between plant/environment. More ClChi members were induced by Race 2 of Fusarium oxysporum f. sp. niveum, and eight genes were expressed at higher levels on the seventh day after inoculation with Races 1 and 2, suggesting that these genes play a key role in the resistance of watermelon to Fusarium wilt. Collectively, these results improve knowledge of the chitinase gene family in watermelon species and help to elucidate the roles played by chitinases in the responses of watermelon to various stresses.

Functional verification and screening of protein interacting with the slPHB3

slPHB3 was cloned from Salix linearistipularis, the amino acid sequence blast and phylogenetic tree analysis showed that slPHB3 has the most similarity with PHB3 from Populus trichocarpa using DNAMAN software and MEGA7 software. RT-qPCR results confirmed that the expression of slPHB3 was induced obviously under stress treatments. The growth of recombinant yeast cells was better than that of the control group under the stress treatment, indicating that slPHB3 may be involved in the stress response of yeast cells. The transgenic tobacco was treated with different concentrations of NaCl, NaHCO3 and H2O2, fresh weigh of overexpression tobacco were heavier than wild-types. The results showed that transgenic tobacco was more tolerant to salt and oxidation than wild-type tobacco. Expression of important genes including NHX1 and SOS1 in salt stress response pathways are steadily higher in overexpression tobacco than that in wild-types. We identified 17 proteins interacting with slPHB3 by yeast two-hybrid technique, most of these proteins were relation to the stresses. The salt tolerance of slPHB3 expressing yeast and slPHB3 overexpressing plants were better than that of the control. Ten stress-related proteins may interact with slPHB3, which preliminarily indicated that slPHB3 had a certain response relationship with salt stress. The study of slPHB3 under abiotic stress can improve our understanding of PHB3 gene function.

Application of the NanoString nCounter System as an Alternative Method to Investigate Molecular Mechanisms Involved in Host Plant Responses to Plasmodiophora brassicae

Clubroot, caused by the soilborne pathogen Plasmodiophora brassicae, is an important disease of canola (Brassica napus) and other crucifers. The recent application of RNA sequencing (RNA-seq) technologies to study P. brassicae−host interactions has generated large amounts of gene expression data, improving knowledge of the molecular mechanisms of pathogenesis and host resistance. Quantitative PCR (qPCR) analysis has been widely applied to examine the expression of a limited number of genes and to validate the results of RNA-seq studies, but may not be ideal for analyzing larger suites of target genes or increased sample numbers. Moreover, the need for intermediate steps such as cDNA synthesis may introduce variability that could affect the accuracy of the data generated by qPCR. Here, we report the validation of gene expression data from a previous RNA-seq study of clubroot using the NanoString nCounter System, which achieves efficient gene expression quantification in a fast and simple manner. We first confirm the robustness of the NanoString system by comparing the results with those generated by qPCR and RNA-seq and then discuss the importance of some candidate genes for resistance or susceptibility to P. brassicae in the host. The results show that the expression of genes measured using NanoString have a high correlation with the values obtained using the other two technologies, with R > 0.90 and p < 0.01, and the same expression patterns for most genes. The three methods (qPCR, RNA-seq, and NanoString) were also compared in terms of laboratory procedures, time, and cost. We propose that the NanoString nCounter System is a robust, sensitive, highly reproducible, and simple technology for gene expression analysis. NanoString could become a common alternative to qPCR to validate RNA-seq data or to create panels of genes for use as markers of resistance/susceptibility when plants are challenged with different P. brassicae pathotypes.

Genetic and molecular analysis reveals that two major loci and their interaction confer clubroot resistance in canola introgressed from rutabaga

Clubroot disease caused by Plasmodiophora brassicae is one of the serious threats to canola (Brassica napus L. subsp. napus) production. The evolution of new pathotypes rendering available resistances ineffective compel the introgression of new resistance into canola and extend our understanding of the genetic and molecular basis of the resistance. In this paper, we report the genetic and molecular basis of clubroot resistance in canola, introgressed from a rutabaga (B. napus L. subsp. rapifera Metzg. 'Polycross'), by using a doubled-haploid (DH) mapping population. Whole-genome resequencing (WGRS)-based bulked segregant analysis followed by genetic mapping and expression analysis of the genes in resistant and susceptible DH lines at 7 and 14 d after inoculation were carried out. Following this approach, two major quantitative trait loci (QTL) located at 14.41-15.44 Mb of A03 and at 9.96-11.09 Mb of A08 chromosomes and their interaction was observed to confer resistance to pathotypes 3H, 3A, and 3D. Analysis of the genes from the two QTL regions suggested that decreased expression of sugar transporter genes (BnaA03g29290D and BnaA03g29310D) may play an important role in resistance conferred by the A03 QTL, while increased expression of the toll/interleukin-1 receptor (TIR)-nucleotide binding (NB)-leucine rich repeat (LRR) (TNL) genes (BnaA08g10100D, BnaA08g09220D, and BnaA08g10540D) could be the major determinant of the resistance conferred by the A08 QTL. Single-nucleotide polymorphism (SNP) allele-specific polymerase chain reaction (PCR)-based markers, which could be detected by agarose gel electrophoresis, were also developed from the two QTL regions for use in breeding including pyramiding of multiple clubroot resistance genes.

Genome-Wide Identification and Expression Analyses of the Chitinase Gene Family in Response to White Mold and Drought Stress in Soybean (Glycine max)

Chitinases are enzymes catalyzing the hydrolysis of chitin that are present on the cell wall of fungal pathogens. Here, we identified and characterized the chitinase gene family in cultivated soybean (Glycine max L.) across the whole genome. A total of 38 chitinase genes were identified in the whole genome of soybean. Phylogenetic analysis of these chitinases classified them into five separate clusters, I–V. From a broader view, the I–V classes of chitinases are basically divided into two mega-groups (X and Y), and these two big groups have evolved independently. In addition, the chitinases were unevenly and randomly distributed in 17 of the total 20 chromosomes of soybean, and the majority of these chitinase genes contained few introns (≤2). Synteny and duplication analysis showed the major role of tandem duplication in the expansion of the chitinase gene family in soybean. Promoter analysis identified multiple cis-regulatory elements involved in the biotic and abiotic stress response in the upstream regions (1.5 kb) of chitinase genes. Furthermore, qRT-PCR analysis showed that pathogenic and drought stress treatment significantly induces the up-regulation of chitinase genes belonging to specific classes at different time intervals, which further verifies their function in the plant stress response. Hence, both in silico and qRT-PCR analysis revealed the important role of the chitinases in multiple plant defense responses. However, there is a need for extensive research efforts to elucidate the detailed function of chitinase in various plant stresses. In conclusion, our investigation is a detailed and systematic report of whole genome characterization of the chitinase family in soybean.

What Can We Learn from -Omics Approaches to Understand Clubroot Disease?

Clubroot is one of the most economically significant diseases worldwide. As a result, many investigations focus on both curing the disease and in-depth molecular studies. Although the first transcriptome dataset for the clubroot disease describing the clubroot disease was published in 2006, many different pathogen-host plant combinations have only recently been investigated and published. Articles presenting -omics data and the clubroot pathogen Plasmodiophora brassicae as well as different host plants were analyzed to summarize the findings in the richness of these datasets. Although genome data for the protist have only recently become available, many effector candidates have been identified, but their functional characterization is incomplete. A better understanding of the life cycle is clearly required to comprehend its function. While only a few proteome studies and metabolome analyses were performed, the majority of studies used microarrays and RNAseq approaches to study transcriptomes. Metabolites, comprising chemical groups like hormones were generally studied in a more targeted manner. Furthermore, functional approaches based on such datasets have been carried out employing mutants, transgenic lines, or ecotypes/cultivars of either Arabidopsis thaliana or other economically important host plants of the Brassica family. This has led to new discoveries of potential genes involved in disease development or in (partial) resistance or tolerance to P. brassicae. The overall contribution of individual experimental setups to a larger picture will be discussed in this review.

Characterization of the Chitinase Gene Family in Mulberry (Morus notabilis) and MnChi18 Involved in Resistance to Botrytis cinerea

Chitinase is a hydrolase that uses chitin as a substrate. It plays an important role in plant resistance to fungal pathogens by degrading chitin. Here, we conducted bioinformatics analysis and transcriptome data analysis of the mulberry (Morus notabilis) chitinase gene family to determine its role in the resistance to Botrytis cinerea. A total of 26 chitinase genes were identified, belonging to the GH18 and GH19 families. Among them, six chitinase genes were differentially expressed under the infection of B. cinerea. MnChi18, which significantly responded to B. cinerea, was heterologously expressed in Arabidopsis (Arabidopsis thaliana). The resistance of MnChi18 transgenic Arabidopsis to B. cinerea was significantly enhanced, and after inoculation with B. cinerea, the activity of catalase (CAT) increased and the content of malondialdehyde (MDA) decreased. This shows that overexpression of MnChi18 can protect cells from damage. In addition, our study also indicated that MnChi18 may be involved in B. cinerea resistance through other resistance-related genes. This study provides an important basis for further understanding the function of mulberry chitinase.

A Novel Target (Oxidation Resistant 2) in Arabidopsis thaliana to Reduce Clubroot Disease Symptoms via the Salicylic Acid Pathway without Growth Penalties

The clubroot disease (Plasmodiophora brassicae) is one of the most damaging diseases worldwide among brassica crops. Its control often relies on resistant cultivars, since the manipulation of the disease hormones, such as salicylic acid (SA) alters plant growth negatively. Alternatively, the SA pathway can be increased by the addition of beneficial microorganisms for biocontrol. However, this potential has not been exhaustively used. In this study, a recently characterized protein Oxidation Resistant 2 (OXR2) from Arabidopsis thaliana is shown to increase the constitutive pathway of SA defense without decreasing plant growth. Plants overexpressing AtOXR2 (OXR2-OE) show strongly reduced clubroot symptoms with improved plant growth performance, in comparison to wild type plants during the course of infection. Consequently, oxr2 mutants are more susceptible to clubroot disease. P. brassicae itself was reduced in these galls as determined by quantitative real-time PCR. Furthermore, we provide evidence for the transcriptional downregulation of the gene encoding a SA-methyltransferase from the pathogen in OXR2-OE plants that could contribute to the phenotype.

Pathogenesis-Related Genes of PR1, PR2, PR4, and PR5 Families Are Involved in the Response to Fusarium Infection in Garlic (Allium sativum L.)

Plants of the genus Allium developed a diversity of defense mechanisms against pathogenic fungi of the genus Fusarium, including transcriptional activation of pathogenesis-related (PR) genes. However, the information on the regulation of PR factors in garlic (Allium sativum L.) is limited. In the present study, we identified AsPR genes putatively encoding PR1, PR2, PR4, and PR5 proteins in A. sativum cv. Ershuizao, which may be involved in the defense against Fusarium infection. The promoters of the AsPR1-5 genes contained jasmonic acid-, salicylic acid-, gibberellin-, abscisic acid-, auxin-, ethylene-, and stress-responsive elements associated with the response to plant parasites. The expression of AsPR1c, d, g, k, AsPR2b, AsPR5a, c (in roots), and AsPR4a(c), b, and AsPR2c (in stems and cloves) significantly differed between garlic cultivars resistant and susceptible to Fusarium rot, suggesting that it could define the PR protein-mediated protection against Fusarium infection in garlic. Our results provide insights into the role of PR factors in A. sativum and may be useful for breeding programs to increase the resistance of Allium crops to Fusarium infections.

A phylogenetic study of the members of the MAPK and MEK families across Viridiplantae

Protein phosphorylation is regulated by the activity of enzymes generically known as kinases. One of those kinases is Mitogen-Activated Protein Kinases (MAPK), which operate through a phosphorylation cascade conformed by members from three related protein kinase families namely MAPK kinase kinase (MEKK), MAPK kinase (MEK), and MAPK; these three acts hierarchically. Establishing the evolution of these proteins in the plant kingdom is an interesting but complicated task because the current MAPK, MAPKK, and MAPKKK subfamilies arose from duplications and subsequent sub-functionalization during the early stage of the emergence of Viridiplantae. Here, an in silico genomic analysis was performed on 18 different plant species, which resulted in the identification of 96 genes not previously annotated as components of the MAPK (70) and MEK (26) families. Interestingly, a deeper analysis of the sequences encoded by such genes revealed the existence of putative domains not previously described as signatures of MAPK and MEK kinases. Additionally, our analysis also suggests the presence of conserved activation motifs besides the canonical TEY and TDY domains, which characterize the MAPK family.

Genome-Wide Identification and Expression of Chitinase Class I Genes in Garlic (Allium sativum L.) Cultivars Resistant and Susceptible to Fusarium proliferatum

Vegetables of the Allium genus are prone to infection by Fusarium fungi. Chitinases of the GH19 family are pathogenesis-related proteins inhibiting fungal growth through the hydrolysis of cell wall chitin; however, the information on garlic (Allium sativum L.) chitinases is limited. In the present study, we identified seven class I chitinase genes, AsCHI1–7, in the A. sativum cv. Ershuizao genome, which may have a conserved function in the garlic defense against Fusarium attack. The AsCHI1–7 promoters contained jasmonic acid-, salicylic acid-, gibberellins-, abscisic acid-, auxin-, ethylene-, and stress-responsive elements associated with defense against pathogens. The expression of AsCHI2, AsCHI3, and AsCHI7 genes was constitutive in Fusarium-resistant and -susceptible garlic cultivars and was mostly induced at the early stage of F. proliferatum infection. In roots, AsCHI2 and AsCHI3 mRNA levels were increased in the susceptible and decreased in the resistant cultivar, whereas in cloves, AsCHI7 and AsCHI5 expression was decreased in the susceptible but increased in the resistant plants, suggesting that these genes are involved in the garlic response to Fusarium proliferatum attack. Our results provide insights into the role of chitinases in garlic and may be useful for breeding programs to increase the resistance of Allium crops to Fusarium infections.

Genome-wide identification and expression analysis of the Brassica oleracea L. chitin-binding genes and response to pathogens infections

Chitinase family genes were involved in the response of Brassica oleracea to Fusarium wilt, powdery mildew, black spot and downy mildew. Abstract Chitinase, a category of pathogenesis-related proteins, is believed to play an important role in defending against external stress in plants. However, a comprehensive analysis of the chitin-binding gene family has not been reported to date in cabbage (Brassica oleracea L.), especially regarding the roles that chitinases play in response to various diseases. In this study, a total of 20 chitinase genes were identified using a genome-wide search method. Phylogenetic analysis was employed to classify these genes into two groups. The genes were distributed unevenly across six chromosomes in cabbage, and all of them contained few introns (≤ 2). The results of collinear analysis showed that the cabbage genome contained 1-5 copies of each chitinase gene (excluding Bol035470) identified in Arabidopsis. The heatmap of the chitinase gene family showed that these genes were expressed in various tissues and organs. Two genes (Bol023322 and Bol041024) were relatively highly expressed in all of the investigated tissues under normal conditions, exhibiting the expression characteristics of housekeeping genes. In addition, under four different stresses, namely, Fusarium wilt, powdery mildew, black spot and downy mildew, we detected 9, 5, 8 and 8 genes with different expression levels in different treatments, respectively. Our results may help to elucidate the roles played by chitinases in the responses of host plants to various diseases.

Genome-wide identification of chitinase genes in Thalassiosira pseudonana and analysis of their expression under abiotic stresses

BACKGROUND

The nitrogen-containing polysaccharide chitin is the second most abundant biopolymer on earth and is found in the cell walls of diatoms, where it serves as a scaffold for biosilica deposition. Diatom chitin is an important source of carbon and nitrogen in the marine environment, but surprisingly little is known about basic chitinase metabolism in diatoms.

RESULTS

Here, we identify and fully characterize 24 chitinase genes from the model centric diatom Thalassiosira pseudonana. We demonstrate that their expression is broadly upregulated under abiotic stresses, despite the fact that chitinase activity itself remains unchanged, and we discuss several explanations for this result. We also examine the potential transcriptional complexity of the intron-rich T. pseudonana chitinase genes and provide evidence for two separate tandem duplication events during their evolution.

CONCLUSIONS

Given the many applications of chitin and chitin derivatives in suture production, wound healing, drug delivery, and other processes, new insight into diatom chitin metabolism has both theoretical and practical value.

Genome-wide identification and expression profiling of chitinase genes in tea (Camellia sinensis (L.) O. Kuntze) under biotic stress conditions

Chitinases are a diverse group of enzymes having the ability to degrade chitin. Chitin is the second most abundant polysaccharide on earth, predominantly found in insect exoskeletons and fungal cell walls. In this study, we performed a genome-wide search for chitinase genes and identified a total of 49 chitinases in tea. These genes were categorized into 5 classes, where an expansion of class V chitinases has been observed in comparison to other plant species. Extensive loss of introns in 46% of the GH18 chitinases indicates that an evolutionary pressure is acting upon these genes to lose introns for rapid gene expression. The promoter upstream regions in 65% of the predicted chitinases contain methyl-jasmonate, salicylic acid and defense responsive cis-acting elements, which may further illustrate the possible role of chitinases in tea plant's defense against various pests and pathogens. Differential expression analysis revealed that transcripts of two GH19 chitinases TEA028279 and TEA019397 got upregulated during three different fungal infections in tea. While GH19 chitinase TEA031377 showed an increase in transcript abundance in the two insect infested tea tissues. Semi-quantitative RT-PCR analysis revealed that five GH19 chitinases viz. TEA018892, TEA031484, TEA28279, TEA033470 and TEA031277 showed significant increase in expression in the tea plants challenged with a biotrophic pathogen Exobasidium vexans. The study endeavours in highlighting biotic stress responsive defensive role of chitinase genes in tea.

Genome-Wide Identification and Evolution of Receptor-Like Kinases (RLKs) and Receptor like Proteins (RLPs) in Brassica juncea

Brassica juncea, an allotetraploid species, is an important germplasm resource for canola improvement, due to its many beneficial agronomic traits, such as heat and drought tolerance and blackleg resistance. Receptor-like kinase (RLK) and receptor-like protein (RLP) genes are two types of resistance gene analogues (RGA) that play important roles in plant innate immunity, stress response and various development processes. In this study, genome wide analysis of RLKs and RLPs is performed in B. juncea. In total, 493 RLKs (LysM-RLKs and LRR-RLKs) and 228 RLPs (LysM-RLPs and LRR-RLPs) are identified in the genome of B. juncea, using RGAugury. Only 13.54% RLKs and 11.79% RLPs are observed to be grouped within gene clusters. The majority of RLKs (90.17%) and RLPs (52.83%) are identified as duplicates, indicating that gene duplications significantly contribute to the expansion of RLK and RLP families. Comparative analysis between B. juncea and its progenitor species, B. rapa and B. nigra, indicate that 83.62% RLKs and 41.98% RLPs are conserved in B. juncea, and RLPs are likely to have a faster evolution than RLKs. This study provides a valuable resource for the identification and characterisation of candidate RLK and RLP genes.

Comparative Transcriptome Analysis of Rutabaga (Brassica napus) Cultivars Indicates Activation of Salicylic Acid and Ethylene-Mediated Defenses in Response to Plasmodiophora brassicae

Clubroot, caused by Plasmodiophora brassicae Woronin, is an important soilborne disease of Brassica napus L. and other crucifers. To improve understanding of the mechanisms of resistance and pathogenesis in the clubroot pathosystem, the rutabaga (B. napus subsp. rapifera Metzg) cultivars ‘Wilhelmsburger’ (resistant) and ‘Laurentian’ (susceptible) were inoculated with P. brassicae pathotype 3A and their transcriptomes were analyzed at 7, 14, and 21 days after inoculation (dai) by RNA sequencing (RNA-seq). Thousands of transcripts with significant changes in expression were identified in each host at each time-point in inoculated vs. non-inoculated plants. Molecular responses at 7 and 14 dai supported clear differences in the clubroot response mechanisms of the two genotypes. Both the resistant and the susceptible cultivars activated receptor-like protein (RLP) genes, resistance (R) genes, and genes involved in salicylic acid (SA) signaling as clubroot defense mechanisms. In addition, genes related to calcium signaling and genes encoding leucine-rich repeat (LRR) receptor kinases, the respiratory burst oxidase homolog (RBOH) protein, and transcription factors such as WRKYs, ethylene responsive factors, and basic leucine zippers (bZIPs), appeared to be upregulated in ‘Wilhelmsburger’ to restrict P. brassicae development. Some of these genes are essential components of molecular defenses, including ethylene (ET) signaling and the oxidative burst. Our study highlights the importance of activation of genes associated with SA- and ET-mediated responses in the resistant cultivar. A set of candidate genes showing contrasting patterns of expression between the resistant and susceptible cultivars was identified and includes potential targets for further study and validation through approaches such as gene editing.

Transcriptomic Analysis Reveals Candidate Genes Responsive to Sclerotinia scleroterum and Cloning of the Ss-Inducible Chitinase Genes in Morus laevigata

Sclerotinia sclerotiorum (Ss) is a devastating fungal pathogen that causes Sclerotinia stem rot in rapeseed (Brassica napus), and is also detrimental to mulberry and many other crops. A wild mulberry germplasm, Morus laevigata, showed high resistance to Ss, but the molecular basis for the resistance is largely unknown. Here, the transcriptome response characteristics of M. laevigata to Ss infection were revealed by RNA-seq. A total of 833 differentially expressed genes (DEGs) were detected after the Ss inoculation in the leaf of M. laevigata. After the GO terms and KEGG pathways enrichment analyses, 42 resistance-related genes were selected as core candidates from the upregulated DEGs. Their expression patterns were detected in the roots, stems, leaves, flowers, and fruits of M. laevigata. Most of them (30/42) were specifically or mainly expressed in flowers, which was consistent with the fact that Ss mainly infects plants through floral organs, and indicated that Ss-resistance genes could be induced by pathogen inoculation on ectopic organs. After the Ss inoculation, these candidate genes were also induced in the two susceptible varieties of mulberry, but the responses of most of them were much slower with lower extents. Based on the expression patterns and functional annotation of the 42 candidate genes, we cloned the full-length gDNA and cDNA sequences of the Ss-inducible chitinase gene set (MlChi family). Phylogenetic tree construction, protein interaction network prediction, and gene expression analysis revealed their special roles in response to Ss infection. In prokaryotic expression, their protein products were all in the form of an inclusion body. Our results will help in the understanding of the molecular basis of Ss-resistance in M. laevigata, and the isolated MlChi genes are candidates for the improvement in plant Ss-resistance via biotechnology.

Evaluation of Brassica oleracea accessions for resistance to Plasmodiophora brassicae and identification of genomic regions associated with resistance

Clubroot disease caused by Plasmodiophora brassicae is a challenge to Brassica crop production. Breakdown of resistance controlled by major genes of the Brassica A genome has been reported. Therefore, identification of resistance in the Brassica C genome is needed to broaden the genetic base of resistance in Brassica napus canola. In this study, we evaluated 135 Brassica oleracea accessions, belonging to eight variants of this species to identify resistant accessions as well as to identify the genomic regions associated with resistance to two recently evolved P. brassicae pathotypes, F3-14 (3A) and F-359-13 (5X L-G2). Resistance to these pathotypes was observed more frequently in var. acephala (kale) followed by var. capitata (cabbage); few accessions also carried resistance to both pathotypes. Association mapping using single nucleotide polymorphism (SNP) markers developed through genotyping by sequencing technique identified 10 quantitative trait loci (QTL) from six C-genome chromosomes to be associated with resistance to these pathotypes; among these, two QTL associated with resistance to 3A and one QTL associated with resistance to 5X L-G2 carried ≥3 SNP markers. The 10 QTL identified in this study individually accounted for 8%–18% of the total phenotypic variance. Thus, the results from this study can be used in molecular breeding of Brassica crops for resistance to this disease.

Effects of enhanced UV-B radiation on the interaction between rice and Magnaporthe oryzae in Yuanyang terrace

Enhanced ultraviolet-B (UV-B) radiation affected the growth of rice and Magnaporthe oryzae, and changed the interactions between them. Increased UV-B radiation (5.0 kJ m-2 d-1) on rice leaves in a Yuanyang terrace was conducted before, during, and after infection of the leaves with Magnaporthe oryzae. The relationship between rice blast and UV-B radiation on the disease resistance of rice and the pathogenicity of M. oryzae was studied, and the effects of enhanced UV-B radiation on the interactions between rice and M. oryzae were analysed. The results indicated the following: (1) enhanced UV-B radiation significantly reduced the rice blast disease index, but as infection progressed, the inhibitory effect of UV-B radiation on the disease was weakened. (2) UV-B radiation treatment before infection with M. oryzae (UV-B + M.) significantly increased the activity of the enzymes related to disease resistance (phenylalanine ammonia lyase, lipoxygenase, chitinase, and β-1,3-glucanase), and the plants exhibited light-induced resistance. (3) Exposure to UV-B radiation after M. oryzae infection (M. + UV-B) did not induce disease course-related protein (PR) activity, but the content of soluble sugar increased. The osmotic stress caused by pathogenic fungi infection was alleviated by active accumulation of soluble sugar; due to this lack of nutrients, it was difficult for the rice blast fungus to grow. (4) Enhanced UV-B radiation significantly inhibited the production of conidia by M. oryzae, and the expression of the pathogenic genes Chitinase, MGP1, MAGB, and CPKA was significantly downregulated. The pathogenicity of M. oryzae was reduced by UV-B radiation. The resistance of rice leaves was weakened by simultaneous exposure to UV-B radiation and M. oryzae (UV-B/M.). Hence, UV-B radiation can weaken the infectivity of M. oryzae, improve the resistance of traditional rice, and contain the disease.

Genome-wide identification and gene expression analysis of SOS family genes in tuber mustard (Brassica juncea var. tumida)

The Salt Overly Sensitive (SOS) pathway in Arabidopsis thaliana plays important roles in maintaining appropriate ion homeostasis in the cytoplasm and regulating plant tolerance to salinity. However, little is known about the details regarding SOS family genes in the tuber mustard crop (Brassica juncea var. tumida). Here, 12 BjSOS family genes were identified in the B. juncea var. tumida genome including two homologous genes of SOS1, one and three homologs of SOS2 and SOS3, two homologs of SOS4, two homologs of SOS5 and two homologs of SOS6, respectively. The results of conserved motif analysis showed that these SOS homologs contained similar protein structures. By analyzing the cis-elements in the promoters of those BjSOS genes, several hormone- and stress-related cis-elements were found. The results of gene expression analysis showed that the homologous genes were induced by abiotic stress and pathogen. These findings indicate that BjSOS genes play crucial roles in the plant response to biotic and abiotic stresses. This study provides valuable information for further investigations of BjSOS genes in tuber mustard.

Comprehensive Analysis of the Chitinase Gene Family in Cucumber (Cucumis sativus L.): From Gene Identification and Evolution to Expression in Response to Fusarium oxysporum

Chitinases, a subgroup of pathogenesis-related proteins, are responsible for catalyzing the hydrolysis of chitin. Accumulating reports indicate that chitinases play a key role in plant defense against chitin-containing pathogens and are therefore good targets for defense response studies. Here, we undertook an integrated bioinformatic and expression analysis of the cucumber chitinases gene family to identify its role in defense against Fusarium oxysporum f. sp. cucumerinum. A total of 28 putative chitinase genes were identified in the cucumber genome and classified into five classes based on their conserved catalytic and binding domains. The expansion of the chitinase gene family was due mainly to tandem duplication events. The expression pattern of chitinase genes was organ-specific and 14 genes were differentially expressed in response to F. oxysporum challenge of fusarium wilt-susceptible and resistant lines. Furthermore, a class I chitinase, CsChi23, was constitutively expressed at high levels in the resistant line and may play a crucial role in building a basal defense and activating a rapid immune response against F. oxysporum. Whole-genome re-sequencing of both lines provided clues for the diverse expression patterns observed. Collectively, these results provide useful genetic resource and offer insights into the role of chitinases in cucumber-F. oxysporum interaction.

Genome-Wide Identification and Expression Analysis of the Metacaspase Gene Family in Gossypium Species

Metacaspases (MCs) are cysteine proteases that are important for programmed cell death (PCD) in plants. In this study, we identified 89 MC genes in the genomes of four Gossypium species (Gossypium raimondii, Gossypium barbadense, Gossypium hirsutum, and Gossypium arboreum), and classified them as type-I or type-II genes. All of the type-I and type-II MC genes contain a sequence encoding the peptidase C14 domain. During developmentally regulated PCD, type-II MC genes may play an important role related to fiber elongation, while type-I genes may affect the thickening of the secondary wall. Additionally, 13 genes were observed to be differentially expressed between two cotton lines with differing fiber strengths, and four genes (GhMC02, GhMC04, GhMC07, and GhMC08) were predominantly expressed in cotton fibers at 5–30 days post-anthesis (DPA). During environmentally induced PCD, the expression levels of four genes were affected in the root, stem, and leaf tissues within 6 h of an abiotic stress treatment. In general, the MC gene family affects the development of cotton fibers, including fiber elongation and fiber thickening while four prominent fiber- expressed genes were identified. The effects of the abiotic stress and hormone treatments imply that the cotton MC gene family may be important for fiber development. The data presented herein may form the foundation for future investigations of the MC gene family in Gossypium species.

Genome-Wide Identification and Expression Analyses of the Chitinases under Cold and Osmotic stress in Ammopiptanthus nanus

Chitinase is a kind of hydrolase with chitin as a substrate and is proposed to play an essential role in plant defense system by functioning against fungal pathogens through degrading chitin. Recent studies indicated chitinase is also involved in abiotic stress response in plants, helping plants to survive in stressful environments. A. nanus, a rare evergreen broad-leaved shrub distrusted in deserts in Central Asia, exhibits a high level of tolerance to drought and low temperature stresses. To identify the chitinase gene involved in drought and low temperature responses in A. nanus, we performed genome-wide identification, classification, sequence alignment, and spatio-temporal gene expression analysis of the chitinases in A. nanus under osmotic and low temperature stress. A total of 32 chitinase genes belonging to glycosyl hydrolase 18 (GH18) and GH19 families were identified from A. nanus. Class III chitinases appear to be amplified quantitatively in A. nanus, and their genes carry less introns, indicating their involvement in stress response in A. nanus. The expression level of the majority of chitinases varied in leaves, stems, and roots, and regulated under environmental stress. Some chitinases, such as EVM0022783, EVM0020238, and EVM0003645, are strongly induced by low temperature and osmotic stress, and the MYC/ICE1 (inducer of CBF expression 1) binding sites in promoter regions may mediate the induction of these chitinases under stress. These chitinases might play key roles in the tolerance to these abiotic stress in A. nanus and have potential for biotechnological applications. This study provided important data for understanding the biological functions of chitinases in A. nanus.

Genome-Wide Identification and Gene Expression Analysis of ABA Receptor Family Genes in Brassica juncea var. tumida

Abscisic acid (ABA) plays important roles in multiple physiological processes, such as plant response to stresses and plant development. The ABA receptors pyrabactin resistance (PYR)/ PYR1-like (PYL)/regulatory components of ABA receptor (RCAR) play a crucial role in ABA perception and signaling. However, little is known about the details regarding PYL family genes in Brassica juncea var. tumida. Here, 25 PYL family genes were identified in B. juncea var. tumida genome, including BjuPYL3, BjuPYL4s, BjuPYL5s, BjuPYL6s, BjuPYL7s, BjuPYL8s, BjuPYL10s, BjuPYL11s, and BjuPYL13. The results of phylogenic analysis and gene structure showed that the PYL family genes performed similar gene characteristics. By analyzing cis-elements in the promoters of those BjuPYLs, several hormone and stress related cis-elements were found. The results of gene expression analysis showed that the ABA receptor homologous genes were induced by abiotic and biotic stress. The tissue-specific gene expression patterns of BjuPYLs also suggested those genes might regulate the stem swelling during plant growth. These findings indicate that BjuPYLs are involved in plant response to stresses and organ development. This study provides valuable information for further functional investigations of PYL family genes in B. juncea var. tumida.

Genome-wide identification and characterization of Chitinase gene family in Brassica juncea and Camelina sativa in response to Alternaria brassicae

Chitinases belong to the group of Pathogenesis-related (PR) proteins that provides protection against fungal pathogens. This study presents the, genome-wide identification and characterization of chitinase gene family in two important oilseed crops B. juncea and C. sativa belonging to family Brassicaceae. We have identified 47 and 79 chitinase genes in the genomes of B. juncea and C. sativa, respectively. Phylogenetic analysis of chitinases in both the species revealed four distinct sub-groups, representing different classes of chitinases (I-V). Microscopic and biochemical study reveals the role of reactive oxygen species (ROS) scavenging enzymes in disease resistance of B. juncea and C. sativa. Furthermore, qRT-PCR analysis showed that expression of chitinases in both B. juncea and C. sativa was significantly induced after Alternaria brassicae infection. However, the fold change in chitinase gene expression was considerably higher in C. sativa compared to B. juncea, which further proves their role in C. sativa disease resistance to A. brassicae. This study provides comprehensive analysis on chitinase gene family in B. juncea and C. sativa and in future may serve as a potential candidate for improving disease resistance in B. juncea through transgenic approach.

Genome-Wide Analysis of Auxin Receptor Family Genes in Brassica juncea var. tumida

Transport inhibitor response 1/auxin signaling f-box proteins (TIR1/AFBs) play important roles in the process of plant growth and development as auxin receptors. To date, no information has been available about the characteristics of the TIR1/AFB gene family in Brassica juncea var. tumida. In this study, 18 TIR1/AFB genes were identified and could be clustered into six groups. The genes are located in 11 of 18 chromosomes in the genome of B. juncea var. tumida, and similar gene structures are found for each of those genes. Several cis-elements related to plant response to phytohormones, biotic stresses, and abiotic stresses are found in the promoter of BjuTIR1/AFB genes. The results of qPCR analysis show that most genes have differential patterns of expression among six tissues, with the expression levels of some of the genes repressed by salt stress treatment. Some of the genes are also responsive to pathogen Plasmodiophora brassicae treatment. This study provides valuable information for further studies as to the role of BjuTIR1/AFB genes in the regulation of plant growth, development, and response to abiotic stress.

Genome-Wide Identification of Cyclophilin Gene Family in Cotton and Expression Analysis of the Fibre Development in Gossypium barbadense

Cyclophilins (CYPs) are a member of the immunophilin superfamily (in addition to FKBPs and parvulins) and play a significant role in peptidyl-prolyl cis-trans isomerase (PPIase) activity. Previous studies have shown that CYPs have important functions in plants, but no genome-wide analysis of the cotton CYP gene family has been reported, and the specific biological function of this gene is still elusive. Based on the release of the cotton genome sequence, we identified 75, 78, 40 and 38 CYP gene sequences from G. barbadense, G. hirsutum, G. arboreum, and G. raimondii, respectively; 221 CYP genes were unequally located on chromosomes. Phylogenetic analysis showed that 231 CYP genes clustered into three major groups and eight subgroups. Collinearity analysis showed that segmental duplications played a significant role in the expansion of CYP members in cotton. There were light-responsiveness, abiotic-stress and hormone-response elements upstream of most of the CYPs. In addition, the motif composition analysis revealed that 49 cyclophilin proteins had extra domains, including TPR (tetratricopeptide repeat), coiled coil, U-box, RRM (RNA recognition motif), WD40 (RNA recognition motif) and zinc finger domains, along with the cyclophilin-like domain (CLD). The expression patterns based on qRT-PCR showed that six CYP expression levels showed greater differences between Xinhai21 (long fibres, G. barbadense) and Ashmon (short fibres, G. barbadense) at 10 and 20 days postanthesis (DPA). These results signified that CYP genes are involved in the elongation stage of cotton fibre development. This study provides a valuable resource for further investigations of CYP gene functions and molecular mechanisms in cotton.

Mining of Brassica-Specific Genes (BSGs) and Their Induction in Different Developmental Stages and under Plasmodiophora brassicae Stress in Brassica rapa

Orphan genes, also called lineage-specific genes (LSGs), are important for responses to biotic and abiotic stresses, and are associated with lineage-specific structures and biological functions. To date, there have been no studies investigating gene number, gene features, or gene expression patterns of orphan genes in Brassica rapa. In this study, 1540 Brassica-specific genes (BSGs) and 1824 Cruciferae-specific genes (CSGs) were identified based on the genome of Brassica rapa. The genic features analysis indicated that BSGs and CSGs possessed a lower percentage of multi-exon genes, higher GC content, and shorter gene length than evolutionary-conserved genes (ECGs). In addition, five types of BSGs were obtained and 145 out of 529 real A subgenome-specific BSGs were verified by PCR in 51 species. In silico and semi-qPCR, gene expression analysis of BSGs suggested that BSGs are expressed in various tissue and can be induced by Plasmodiophora brassicae. Moreover, an A/C subgenome-specific BSG, BSGs1, was specifically expressed during the heading stage, indicating that the gene might be associated with leafy head formation. Our results provide valuable biological information for studying the molecular function of BSGs for Brassica-specific phenotypes and biotic stress in B. rapa.