Genome-Wide Identification and Analysis of CC-NBS-LRR Family in Response to Downy Mildew and Black Rot in Chinese Cabbage

The nucleotide-binding site–leucine-rich repeat (NBS–LRR) gene family is the largest group of plant disease resistance (R) genes widespread in response to viruses, bacteria, and fungi usually involved in effector triggered immunity (ETI). Forty members of the Chinese cabbage CC type NBS–LRR family were investigated in this study. Gene and protein characteristics, such as distributed locations on chromosomes and gene structures, were explored through comprehensive analysis. CC–NBS–LRR proteins were classified according to their conserved domains, and the phylogenetic relationships of CC–NBS–LRR proteins in Brassica rapa, Arabidopsis thaliana, and Oryza sativa were compared. Moreover, the roles of BrCC–NBS–LRR genes involved in pathogenesis-related defense were studied and analyzed. First, the expression profiles of BrCC–NBS–LRR genes were detected by inoculating with downy mildew and black rot pathogens. Second, sensitive and resistant Chinese cabbage inbred lines were screened by downy mildew and black rot. Finally, the differential expression levels of BrCC–NBS–LRR genes were monitored at 0, 1, 3, 6, 12 and 24 h for short and 0, 3, 5, 7, 10 and 14 days for long inoculation time. Our study provides information on BrCC–NBS–LRR genes for the investigation of the functions and mechanisms of CC-NBS-LRR genes in Chinese cabbage.     

Genome-Wide Identification,Classification and Expression Analysis of m6A Gene Family in Solanum lycopersicum

Advanced knowledge of messenger RNA (mRNA) N⁶-methyladenosine (m⁶A) and DNA N⁶-methyldeoxyadenosine (6 mA) redefine our understanding of these epigenetic modifications. Both m⁶A and 6mA carry important information for gene regulation, and the corresponding catalytic enzymes sometimes belong to the same gene family and need to be distinguished. However, a comprehensive analysis of the m⁶A gene family in tomato remains obscure. Here, 24 putative m⁶A genes and their family genes in tomato were identified and renamed according to BLASTP and phylogenetic analysis. Chromosomal location, synteny, phylogenetic, and structural analyses were performed, unravelling distinct evolutionary relationships between the MT-A70, ALKBH, and YTH protein families, respectively. Most of the 24 genes had extensive tissue expression, and 9 genes could be clustered in a similar expression trend. Besides, SlYTH1 and SlYTH3A showed a different expression pattern in leaf and fruit development. Additionally, qPCR data revealed the expression variation under multiple abiotic stresses, and LC-MS/MS determination exhibited that the cold stress decreased the level of N⁶ 2′-O dimethyladenosine (m⁶Am). Notably, the orthologs of newly identified single-strand DNA (ssDNA) 6mA writer–eraser–reader also existed in the tomato genome. Our study provides comprehensive information on m⁶A components and their family proteins in tomato and will facilitate further functional analysis of the tomato N⁶-methyladenosine modification genes.     

Genome-Wide Identification and Analysis of the Polycomb Group Family in Medicago truncatula

Polycomb group (PcG) proteins, which are important epigenetic regulators, play essential roles in the regulatory networks involved in plant growth, development, and environmental stress responses. Currently, as far as we know, no comprehensive and systematic study has been carried out on the PcG family in Medicago truncatula. In the present study, we identified 64 PcG genes with distinct gene structures from the M. truncatula genome. All of the PcG genes were distributed unevenly over eight chromosomes, of which 26 genes underwent gene duplication. The prediction of protein interaction network indicated that 34 M. truncatula PcG proteins exhibited protein–protein interactions, and MtMSI1;4 and MtVRN2 had the largest number of protein–protein interactions. Based on phylogenetic analysis, we divided 375 PcG proteins from 27 species into three groups and nine subgroups. Group I and Group III were composed of five components from the PRC1 complex, and Group II was composed of four components from the PRC2 complex. Additionally, we found that seven PcG proteins in M. truncatula were closely related to the corresponding proteins of Cicer arietinum. Syntenic analysis revealed that PcG proteins had evolved more conservatively in dicots than in monocots. M. truncatula had the most collinearity relationships with Glycine max (36 genes), while collinearity with three monocots was rare (eight genes). The analysis of various types of expression data suggested that PcG genes were involved in the regulation and response process of M. truncatula in multiple developmental stages, in different tissues, and for various environmental stimuli. Meanwhile, many differentially expressed genes (DEGs) were identified in the RNA-seq data, which had potential research value in further studies on gene function verification. These findings provide novel and detailed information on the M. truncatula PcG family, and in the future it would be helpful to carry out related research on the PcG family in other legumes.     

Genome-Wide Analysis of Type-III Polyketide Synthases in Wheat and Possible Roles in Wheat Sheath-Blight Resistance

The enzymes in the chalcone synthase family, also known as type-III polyketide synthases (PKSs), play important roles in the biosynthesis of various plant secondary metabolites and plant adaptation to environmental stresses. There have been few detailed reports regarding the gene and tissue expression profiles of the PKS (TaPKS) family members in wheat (Triticum aestivum L.). In this study, 81 candidate TaPKS genes were identified in the wheat genome, which were designated as TaPKS1–81. Phylogenetic analysis divided the TaPKS genes into two groups. TaPKS gene family expansion mainly occurred via tandem duplication and fragment duplication. In addition, we analyzed the physical and chemical properties, gene structures, and cis-acting elements of TaPKS gene family members. RNA-seq analysis showed that the expression of TaPKS genes was tissue-specific, and their expression levels differed before and after infection with Rhizoctonia cerealis. The expression levels of four TaPKS genes were also analyzed via qRT-PCR after treatment with methyl jasmonate, salicylic acid, abscisic acid, and ethylene. In the present study, we systematically identified and analyzed TaPKS gene family members in wheat, and our findings may facilitate the cloning of candidate genes associated with resistance to sheath blight in wheat.     

Dynamic Expression,Differential Regulation and Functional Diversity of the CNGC Family Genes in Cotton

Cyclic nucleotide-gated channels (CNGCs) constitute a family of non-selective cation channels that are primarily permeable to Ca²⁺ and activated by the direct binding of cyclic nucleotides (i.e., cAMP and cGMP) to mediate cellular signaling, both in animals and plants. Until now, our understanding of CNGCs in cotton (Gossypium spp.) remains poorly addressed. In the present study, we have identified 40, 41, 20, 20, and 20 CNGC genes in G. hirsutum, G. barbadense, G. herbaceum, G. arboreum, and G. raimondii, respectively, and demonstrated characteristics of the phylogenetic relationships, gene structures, chromosomal localization, gene duplication, and synteny. Further investigation of CNGC genes in G. hirsutum, named GhCNGC1-40, indicated that they are not only extensively expressed in various tissues and at different developmental stages, but also display diverse expression patterns in response to hormones (abscisic acid, salicylic acid, methyl jasmonate, ethylene), abiotic (salt stress) and biotic (Verticillium dahlia infection) stimuli, which conform with a variety of cis-acting regulatory elements residing in the promoter regions; moreover, a set of GhCNGCs are responsive to cAMP signaling during cotton fiber development. Protein–protein interactions supported the functional aspects of GhCNGCs in plant growth, development, and stress responses. Accordingly, the silencing of the homoeologous gene pair GhCNGC1&18 and GhCNGC12&31 impaired plant growth and development; however, GhCNGC1&18-silenced plants enhanced Verticillium wilt resistance and salt tolerance, whereas GhCNGC12&31-silenced plants had opposite effects. Together, these results unveiled the dynamic expression, differential regulation, and functional diversity of the CNGC family genes in cotton. The present work has laid the foundation for further studies and the utilization of CNGCs in cotton genetic improvement.     

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1–cullin–FBX) belonging to the E3 ligase family, which is involved in the ubiquitin–proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.     

Catalase (CAT) Gene Family in Rapeseed (Brassica napus L.): Genome-Wide Analysis,Identification, and Expression Pattern in Response to Multiple Hormones and Abiotic Stress Conditions

Catalase (CAT) is an antioxidant enzyme expressed by the CAT gene family and exists in almost all aerobic organisms. Environmental stresses induce the generation of reactive oxygen species (ROS) that eventually hinder plant growth and development. The CAT enzyme translates the hydrogen peroxide (H₂O₂) to water (H₂O) and reduce the ROS levels to shelter the cells’ death. So far, the CAT gene family has not been reported in rapeseed (Brassica napus L.). Therefore, a genome-wide comprehensive analysis was conducted to classify the CAT genes in the rapeseed genome. The current study identified 14 BnCAT genes in the rapeseed genome. Based on phylogenetic and synteny analysis, the BnCATs belong to four groups (Groups I–IV). A gene structure and conserved motif analysis showed that Group I, Group II, and Group IV possess almost the same intron/exon pattern, and an equal number of motifs, while Group III contains diverse structures and contain 15 motifs. By analyzing the cis-elements in the promoters, we identified five hormone-correlated responsive elements and four stress-related responsive elements. Further, six putative bna-miRNAs were also identified, targeting three genes (BnCAT4, BnCAT6, and BnCAT8). Gene ontology (GO) enrichment analysis showed that the BnCAT genes were largely related to cellular organelles, ROS response, stimulus response, stress response, and antioxidant enzymes. Almost 10 BnCAT genes showed higher expression levels in different tissues, i.e., root, leaf, stem, and silique. The expression analysis showed that BnCAT1–BnCAT3 and BnCAT11–BnCAT13 were significantly upregulated by cold, salinity, abscisic acid (ABA), and gibberellic acid (GA) treatment, but not by drought and methyl jasmonate (MeJA). Notably, most of the genes were upregulated by waterlogging stress, except BnCAT6, BnCAT9, and BnCAT10. Our results opened new windows for future investigations and provided insights into the CAT family genes in rapeseed.     

Sequence Analysis and Functional Verification of the Effects of Three Key Structural Genes,PdTHC2’GT,PdCHS and PdCHI,on the Isosalipurposide Synthesis Pathway in Paeonia delavayi var. lutea

Isosalipurposide (ISP) is the most important yellow pigment in tree peony. In ISP biosynthesis, CHS catalyzes 1-molecule coumaroyl-CoA and 3-molecule malonyl-CoA to form 2′,4′,6′,4-tetrahyroxychalcone (THC), and THC generates a stable ISP in the vacuole under the action of chalcone2′-glucosyltransferases (THC2′GT). In tree peony, the details of the THC2’GT gene have not yet been reported. In this study, the candidate THC2’GT gene (PdTHC2’GT) in Paeonia delavayi var. lutea was screened. At the same time, we selected the upstream CHS gene (PdCHS) and the competitive CHI gene (PdCHI) to study the biosynthesis pathway of ISP. We successfully cloned three genes and sequenced them; subcellular localization showed that the three genes were located in the nucleus and cytoplasm. The overexpression of PdTHC2’GT in tobacco caused the accumulation of ISP in tobacco petals, which indicated that PdTHC2’GT was the key structural gene in the synthesis of ISP. After the overexpression of PdCHS and PdCHI in tobacco, the accumulation of anthocyanins in tobacco petals increased to different degrees, showing the role of PdCHS and PdCHI in anthocyanin accumulation. The analysis of NtCHS and NtCHI of transgenic tobacco lines by qRT-PCR showed that the THC2’GT gene could increase the expression of CHS. THC2’GT and CHI were found to be competitive; hence, the overexpression of THC2’GT could lead to a decrease in CHI expression. The CHS gene and CHI gene could increase the expression of each other. In conclusion, we verified the key structural gene PdTHC2’GT and studied the operation of the genes in its upstream and competitive pathway, providing a new perspective for the biosynthesis of ISP and new candidate genes for the directional breeding of tree peony.     

Molecular Identification of the G-Protein Genes and Their Expression Profiles in Response to Nitrogen Deprivation in Brassica napus

Heterotrimeric guanine nucleotide binding protein (G-protein) consisting of Gα, Gβ, and Gγ subunits is one of the key signal transducers in plants. Recent studies indicated that G-protein has been proposed as an important mediator of nitrogen responses in rice, wheat, and Arabidopsis. However, little is known about these G-proteins in Brassica napus (B. napus), except for three identified G-proteins, BnGA1, BnGB1, and BnGG2. Therefore, the aim of the present study is to characterize the members of the G-protein gene family in allotetraploid B. napus and to analyze their expression profiles in response to nitrogen deprivation. In total, 21 G-protein family members were identified in B. napus, encoding two Gα, six Gβ, and 13 Gγ. Sequence and phylogenetic analyses showed that although genome-wide triploid events increased the number of genes encoding Gα, Gβ, and Gγ subunits, the gene structure and protein properties of the genes encoding each G-protein subunit were extremely conserved. Collinearity analysis showed that most G-protein genes in B. napus had syntenic relationships with G-protein members of Arabidopsis, Brassica rape (B. rapa), and Brassica oleracea (B. oleracea). Expression profile analysis indicated that Gα and C-type Gγ genes (except BnGG10 and BnGG12 were highly expressed in flower and ovule) were barely expressed in most organs, whereas most Gβ and A-type Gγ genes tended to be highly expressed in most organs. G-protein genes also showed various expression patterns in response to nitrogen-deficient conditions. Under nitrogen deficiency, Gα and five C-type Gγ genes were upregulated initially in roots, while in leaves, Gα was downregulated initially and five C-type Gγ genes were highly expressed in different times. These results provide a complex genetic dissection of G-protein genes in B. napus, and insight into the biological functions of G-protein genes in response to nitrogen deficiency.