Differential evolution of members of the rhomboid gene family with conservative and divergent patterns

Comprehensive Review on Plant Cytochrome P450 Evolution: Copy Number, Diversity, and Motif Analysis From Chlorophyta to Dicotyledoneae

Cytochrome P450 enzymes (CYPs) are widely distributed among various plant groups and constitute approximately 1% of the total number of protein-coding genes. Extensive studies suggest that CYPs are involved in nearly all molecular processes that occur in plants. Over the past two decades, the identification of CYP genes has expanded rapidly, with more than 40,000 CYP genes and 819 CYP families being discovered. Copy number variation is a significant evolutionary characteristic of gene families, yet a systematic characterization of the copy evolution patterns in plant CYP gene families has been lacking, resulting in confusion and challenges in understanding CYP functions. To address these concerns, this review provides comprehensive statistics and analyses of the copy number and diversity of almost all plant CYP gene families, focusing on CYP evolution from Chlorophyta to Dicotyledoneae. Additionally, we examined the subfamily characteristics of certain CYP families with restricted copy changes and identified several CYP subfamilies that play pivotal roles in this event. Furthermore, we analyzed the structural conservation of CYPs across different taxa and compiled a comprehensive database to support plant CYP studies. Our analysis revealed differences in the six core conserved motifs of plant CYP proteins among various clans and plant taxa, while demonstrating similar conservation patterns for the ERR (glutamic acid-arginine-arginine) triad motifs. These findings will significantly facilitate the understanding of plant CYP gene evolution and metabolic diversity and serve as a valuable reference for researchers studying CYP enzymes.

Genome-wide identification, gene expression and haplotype analysis of the rhomboid-like gene family in wheat (Triticum aestivum L.)

The rhomboid-like (RBL) gene encodes serine protease, which plays an important role in the response to cell development and diverse stresses. However, genome-wide identification, expression profiles, and haplotype analysis of the RBL family genes have not been performed in wheat (Triticum aestivum L.). This study investigated the phylogeny and diversity of the RBL family genes in the wheat genome through various approaches, including gene structure analysis, evolutionary relationship analysis, promoter cis-acting element analysis, expression pattern analysis, and haplotype analysis. The 41 TaRBL genes were identified and divided into five subfamilies in the wheat genome. RBL family genes were expanded through segmented duplication and purification selection. The cis-element analysis revealed their involvement in various stress responses and plant development. The results of RNA-seq and quantitative real-time-PCR showed that TaRBL genes displayed higher expression levels in developing spike/grain and were differentially regulated under polyethylene glycol, NaCl, and abscisic acid treatments, indicating their roles in grain development and abiotic stress response. A kompetitive allele-specific PCR molecular marker was developed to confirm the single nucleotide polymorphism of TaRBL14a gene in 263 wheat accessions. We found that the elite haplotype TaRBL14a-Hap2 showed a significantly higher 1000-grain weight than TaRBL14a-Hap11 in at least three environments, and the TaRBL14a-Hap2 was positively selected in wheat breeding. The findings will provide a good insight into the evolutionary and functional characteristics of the TaRBL genes family in wheat and lay the foundation for future exploration of the regulatory mechanisms of TaRBL genes in plant growth and development, as well as their response to abiotic stresses.

Physico-chemical features and functional relevance of tomato rhomboid proteases

In plants, regulated intramembrane proteolysis (RIP) is crucial for proper growth, development, and stress management. Rhomboid proteases (RPs) residing in the membrane play a vital role in orchestrating RIP. Although RPs can be found in most sequenced genomes, tomato rhomboids (SlRPs) have not yet been studied. Using alternative and comprehensive strategies, we found ten SlRPs encoded in the tomato genome. These SlRPs possess signature motifs and transmembrane domains, showing structural similarity to other members of the RP family. Also, SlRPs are genetically related to other known RPs of the Solanaceae family. Seven of the SlRPs retain serine-histidine catalytic dyads, making them proteolytically active, while three iRhoms lack the dyad and other structural motifs. Although SlRPs could have functional redundancy, their distribution and expression pattern indicate tissue specificity and responsiveness to specific external stimuli. The presence of development and stress-response-related cis-elements in the promoters of SlRPs supports this view. Furthermore, our strategically designed substrate-reporter assay shows that SlRPs have proteolytic activity similar to that of known RPs. This study provides a detailed understanding of all SlRPs and their physico-chemical features, shedding light on their involvement in physiological processes.

Genome-wide identification and expression pattern profiling of the ATP-binding cassette gene family in tea plant (Camelliasinensis)

The ATP-binding cassette (ABC) gene family is one of the largest and oldest protein families, consisting of ATP-driven transporters facilitating substrate transportation across cell membranes. However, little is known about the evolution and biological function of the ABC gene family in tea plants. In this study, we performed a genome-wide identification and expression analysis of genes encoding ABC transporter proteins in Camellia sinensis. Our analysis of 170 ABC genes revealed that CsABCs were unevenly distributed across 15 chromosomes, with an amino acid length ranging from 188 to 2489 aa, molecular weight ranging from 20.29 to 277.34 kDa, and an isoelectric point ranging from 4.89 to 10.63. Phylogenetic analysis showed that CsABCs were divided into eight subfamilies, among which the ABCG subfamily was the most abundant. Furthermore, the subcellular localization of CsABCs indicated that they were present in various organelles. Collinearity analysis between the tea plant and Arabidopsis thaliana genomes revealed that the CsABC genes were homologous to the AtABC genes. Large gene fragment duplication analysis identified ten gene pairs as tandem repeats, and interaction network analysis demonstrated that CsABCs interacted with various types of target genes, with protein interactions also occurring within the family. Tissue expression analysis indicated that CsABCs were highly expressed in roots, stems, and leaves and were easily induced by drought and cold stress. Moreover, qRT-PCR analysis of the relative expression level of the gene under drought and cold stress correlated with the sequencing results. Identifying ABC genes in tea plants lays a foundation for the classification and functional analysis of ABC family genes, which can facilitate molecular breeding and the development of new tea varieties.

Evolution and functional diversification of catalase genes in the green lineage

Background

Catalases (CATs) break down hydrogen peroxide into water and oxygen to prevent cellular oxidative damage, and play key roles in the development, biotic and abiotic stresses of plants. However, the evolutionary relationships of the plant CAT gene family have not been systematically reported.

Results

Here, we conducted genome-wide comparative, phylogenetic, and structural analyses of CAT orthologs from 29 out of 31 representative green lineage species to characterize the evolution and functional diversity of CATs. We found that CAT genes in land plants were derived from core chlorophytes and detected a lineage-specific loss of CAT genes in Fabaceae, suggesting that the CAT genes in this group possess divergent functions. All CAT genes were split into three major groups (group α, β1, and β2) based on the phylogeny. CAT genes were transferred from bacteria to core chlorophytes and charophytes by lateral gene transfer, and this led to the independent evolution of two types of CAT genes: α and β types. Ten common motifs were detected in both α and β groups, and β CAT genes had five unique motifs, respectively. The findings of our study are inconsistent with two previous hypotheses proposing that (i) new CAT genes are acquired through intron loss and that (ii) the Cys-343 residue is highly conserved in plants. We found that new CAT genes in most higher plants were produced through intron acquisition and that the Cys-343 residue was only present in monocots, Brassicaceae and Pp_CatX7 in P. patens, which indicates the functional specificity of the CATs in these three lineages. Finally, our finding that CAT genes show high overall sequence identity but that individual CAT genes showed developmental stage and organ-specific expression patterns suggests that CAT genes have functionally diverged independently.

Conclusions

Overall, our analyses of the CAT gene family provide new insights into their evolution and functional diversification in green lineage species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08621-6.

Phylogeny and evolution of plant Phytochrome Interacting Factors (PIFs) gene family and functional analyses of PIFs in Brachypodium distachyon

Plant PIFs have been characterized, WGDs contributed to the expansion of class II PIFs; BdPIFs localized in the nucleus; BdPIF4/5C most likely response to high temperature and light stress. Phytochrome interacting factors (PIFs) belong to a small subset of basic helix-loop-helix (bHLH) transcription factors (TFs). As cellular signaling hubs, PIFs integrate multiple external and internal signals to orchestrate the regulation of the transcriptional network, thereby actuating the pleiotropic aspects of downstream morphogenesis. Nevertheless, the origin, phylogeny and function of plant PIFs are not well understood. To elucidate their evolution history and biological function, the comprehensive genomic analysis of the PIF genes was conducted using 40 land plant genomes plus additionally four alga lineages and also performed their gene organizations, sequence features and expression patterns in different subfamilies. In this study, phylogenetic analysis displayed that 246 PIF gene members retrieved from all embryophytes could be divided into three main clades, which were further felled into five distinct classes (Class I-V). The duplications of Class II PIFs were associated specially with whole genome duplication (WGD) events during the plant evolution process. Sequence analysis showed that PIF proteins had a conserved APB motif, and its crucial amino acid residues were relatively high proportion in the average abundance. As expected, subcellular localization analysis revealed that all BdPIF proteins were localized to the nucleus. Especially, BdPIF4/5C showed the highest expression level at high temperature, and the most significant hypocotyl elongation phenotype of overexpression of BdPIFs in Arabidopsis, which was consistent with the function and phenotype of AtPIF4. In brief, our findings provide a novel perspective on the origin and evolutionary history of plant PIFs, and lays a foundation for further investigation on its functions in plant growth and development.

In Silico Predictions of Ecological Plasticity Mediated by Protein Family Expansions in Early-Diverging Fungi

Early-diverging fungi (EDF) are ubiquitous and versatile. Their diversity is reflected in their genome sizes and complexity. For instance, multiple protein families have been reported to expand or disappear either in particular genomes or even whole lineages. The most commonly mentioned are CAZymes (carbohydrate-active enzymes), peptidases and transporters that serve multiple biological roles connected to, e.g., metabolism and nutrients intake. In order to study the link between ecology and its genomic underpinnings in a more comprehensive manner, we carried out a systematic in silico survey of protein family expansions and losses among EDF with diverse lifestyles. We found that 86 protein families are represented differently according to EDF ecological features (assessed by median count differences). Among these there are 19 families of proteases, 43 CAZymes and 24 transporters. Some of these protein families have been recognized before as serine and metallopeptidases, cellulases and other nutrition-related enzymes. Other clearly pronounced differences refer to cell wall remodelling and glycosylation. We hypothesize that these protein families altogether define the preliminary fungal adaptasome. However, our findings need experimental validation. Many of the protein families have never been characterized in fungi and are discussed in the light of fungal ecology for the first time.

Comprehensive Analysis of Evolutionary Characterization and Expression for Monosaccharide Transporter Family Genes in Nelumbo nucifera

Sugar transporters, an important class of transporters for sugar function, regulate many processes associated with growth, maturation, and senescence processes in plants. In this study, a total of 35 NuMSTs were identified in the Nelumbo nucifera genome and grouped by conserved domains and phylogenetic analysis. Additionally, we identified 316 MST genes in 10 other representative plants and performed a comparative analysis with Nelumbo nucifera genes, including evolutionary trajectory, gene duplication, and expression pattern. A large number of analyses across plants and algae indicated that the MST family could have originated from STP and Glct, expanding to form STP and SFP by dispersed duplication. Finally, a quantitative real-time polymerase chain reaction and cis-element analysis showed that some of them may be regulated by plant hormones (e.g., abscisic acid), biotic stress factors, and abiotic factors (e.g., drought, excessive cold, and light). We found that under the four abiotic stress conditions, only NuSTP5 expression was upregulated, generating a stress response, and ARBE and LTR were present in NuSTP5. In summary, our findings are significant for understanding and exploring the molecular evolution and mechanisms of NuMSTs in plants.

InpactorDB: A Classified Lineage-Level Plant LTR Retrotransposon Reference Library for Free-Alignment Methods Based on Machine Learning

Long terminal repeat (LTR) retrotransposons are mobile elements that constitute the major fraction of most plant genomes. The identification and annotation of these elements via bioinformatics approaches represent a major challenge in the era of massive plant genome sequencing. In addition to their involvement in genome size variation, LTR retrotransposons are also associated with the function and structure of different chromosomal regions and can alter the function of coding regions, among others. Several sequence databases of plant LTR retrotransposons are available for public access, such as PGSB and RepetDB, or restricted access such as Repbase. Although these databases are useful to identify LTR-RTs in new genomes by similarity, the elements of these databases are not fully classified to the lineage (also called family) level. Here, we present InpactorDB, a semi-curated dataset composed of 130,439 elements from 195 plant genomes (belonging to 108 plant species) classified to the lineage level. This dataset has been used to train two deep neural networks (i.e., one fully connected and one convolutional) for the rapid classification of these elements. In lineage-level classification approaches, we obtain up to 98% performance, indicated by the F1-score, precision and recall scores.

Comparative analysis of plant MKK gene family reveals novel expansion mechanism of the members and sheds new light on functional conservation

Background

Mitogen-activated protein kinase (MAPK) cascades play critical functions in almost every aspect of plant growth and development, which regulates many physiological and biochemical processes. As a middle nodal point of the MAPK cascades, although evolutionary analysis of MKK from individual plant families had some reports, their evolutionary history in entire plants is still not clear.

Results

To better understand the evolution and function of plant MKKs, we performed systematical molecular evolutionary analysis of the MAPKK gene family and also surveyed their gene organizations, sequence features and expression patterns in different subfamilies. Phylogenetic analysis showed that plant MAPKK fall into five different groups (Group A–E). Majority orthology groups seemed to be a single or low-copy genes in all plant species analyzed in Group B, C and D, whereas group A MKKs undergo several duplication events, generating multiple gene copies. Further analysis showed that these duplication events were on account of whole genome duplications (WGDs) in plants and the duplicate genes maybe have undergone functional divergence. We also found that group E MKKs had mutation with one change of serine or theronine might lead to inactivity originated through the ancient tandem duplicates in monocots. Moreover, we also identified MKK3 integrated NTF2 domain that might have gradually lost the cytoplasmic-nuclear trafficking activity, which suggests that they may involve with the gene function more and more sophistication in the evolutionary process. Moreover, expression analyses indicated that plant MKK genes play probable roles in UV-B signaling.

Conclusion

In general, ancient gene and genome duplications are significantly conducive to the expansion of the plant MKK gene family. Our study reveals two distinct evolutionary patterns for plant MKK proteins and sheds new light on the functional evolution of this gene family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4793-8) contains supplementary material, which is available to authorized users.

Comprehensive analyses of the BES1 gene family in Brassica napus and examination of their evolutionary pattern in representative species

BACKGROUND

The BES1 gene family, an important class of plant-specific transcription factors, play key roles in the BR signal pathway in plants, regulating various development processes. Until now, there has been no comprehensive analysis of the BES1 gene family in Brassica napus, and a cross-genome exploration of their origin, copy number changes, and functional innovation in plants was also not available.

RESULTS

We identified 28 BES1 genes in B. napus from its two subgenomes (AA and CC). We found that 71.43% of them were duplicated in the tetraploidization, and their gene expression showed a prominent subgenome bias in the roots. Additionally, we identified 104 BES1 genes in another 18 representative angiosperms and performed a comparative analysis with B. napus, including evolutionary trajectory, gene duplication, positive selection, and expression pattern. Exploiting the available genome datasets, we performed a large-scale analysis across plants and algae suggested that the BES1 gene family could have originated from group F, expanding to form other groups (A to E) by duplicating or alternatively deleting some domains. We detected an additional domain containing M4 to M8 in exclusively groups F1 and F2. We found evidence that whole-genome duplication (WGD) contributed the most to the expansion of this gene family among examined dicots, while dispersed duplication contributed the most to expansion in certain monocots. Moreover, we inferred that positive selection might have occurred on major phylogenetic nodes during the evolution of plants.

CONCLUSIONS

Grossly, a cross-genome comparative analysis of the BES1 genes in B. napus and other species sheds light on understanding its copy number expansion, natural selection, and functional innovation.

Genome-Wide Identification, Evolution, and Expression Analysis of the ATP-Binding Cassette Transporter Gene Family in Brassica rapa

ATP-binding cassette (ABC) proteins can act as transporters of different substrates across biological membranes by hydrolyzing ATP. However, little information is available about ABC transporters in Brassica rapa, an important leafy vegetable. In the present study, we carried out genome-wide identification, characterization and molecular evolution analyses of ABC gene family in B. rapa and 9 other plant species. A total of 179 B. rapa ABC genes (BraABCs) were identified. Among them, 173 BraABCs were identified on 10 chromosomes. Based on phylogenetic analysis and domain organization, the BraABC family could be grouped into eight subfamilies. BraABCs in the same subfamily showed similar motif composition and exon-intron organization. Common and unique cis-elements involved in the transcriptional regulation were also identified in the promoter regions of BraABCs. Tissue-expression analysis of BraABCs demonstrated their diverse spatiotemporal expression profiles. Influences of the whole genome triplication (WGT) on the evolution of BraABCs were studied in detail. BraABCs were preferentially retained compared with their neighboring genes during diploidization after WGT. Synteny analysis identified 76 pairs of syntenic BraABC paralogs among the three subgenomes of B. rapa, and 10 paralog pairs underwent positive selection with ω (= Ka/Ks) ratios greater than 1. Analyses of the expression patterns of syntenic BraABC paralogs pairs across five tissues and under stress treatments revealed their functional conservation, sub-functionalization, neo-functionalization and pseudogenization during evolution. Our study presents a comprehensive overview of the ABC gene family in B. rapa and will be helpful for the further functional study of BraABCs in plant growth, development, and stress responses.

Substrates and physiological functions of secretase rhomboid proteases

Rhomboids are conserved intramembrane serine proteases with widespread functions. They were the earliest discovered members of the wider rhomboid-like superfamily of proteases and pseudoproteases. The secretase class of rhomboid proteases, distributed through the secretory pathway, are the most numerous in eukaryotes, but our knowledge of them is limited. Here we aim to summarise all that has been published on secretase rhomboids in a concise encyclopaedia of the enzymes, their substrates, and their biological roles. We also discuss emerging themes of how these important enzymes are regulated.

Developmental roles of Rhomboid proteases

Rhomboid proteins have emerged as one of the most tantalizing and diverse families of proteases. Gene duplication events and structural alterations have sculpted the varied roles of this protein family, maintaining a conserved structural core throughout the bacterial, plant and animal kingdoms. Unresolved questions pop up at many junctions. This review will focus on a distinct class of Rhomboid proteins that plays an essential role in development. It will outline the diverse mechanisms by which these proteins are regulated, and the implications on the biological processes they control. While most of the review will deal with Rhomboids in Drosophila, a system that has been studied in the greatest detail, it will also explore parallels and differences in the function of Rhomboids in the flour beetle T. casteneum and the worm C. elegans.

Holophytochrome-Interacting Proteins in Physcomitrella: Putative Actors in Phytochrome Cytoplasmic Signaling

Phytochromes are the principle photoreceptors in light-regulated plant development, primarily acting via translocation of the light-activated photoreceptor into the nucleus and subsequent gene regulation. However, several independent lines of evidence indicate unambiguously that an additional cytoplasmic signaling mechanism must exist. Directional responses in filament tip cells of the moss Physcomitrella patens are steered by phy4 which has been shown to interact physically with the blue light receptor phototropin at the plasma membrane. This complex might perceive and transduce vectorial information leading to cytoskeleton reorganization and finally a directional growth response. We developed yeast two-hybrid procedures using photochemically functional, full-length phy4 as bait in Physcomitrella cDNA library screens and growth assays under different light conditions, revealing Pfr-dependent interactions possibly associated with phytochrome cytoplasmic signaling. Candidate proteins were then expressed in planta with fluorescent protein tags to determine their intracellular localization in darkness and red light. Of 14 candidates, 12 were confirmed to interact with phy4 in planta using bimolecular fluorescence complementation. We also used database information to study their expression patterns relative to those of phy4. We discuss the likely functional characteristics of these holophytochrome-interacting proteins (HIP's) and their possible roles in signaling.

Large-Scale Analyses of Angiosperm Nucleotide-Binding Site-Leucine-Rich Repeat Genes Reveal Three Anciently Diverged Classes with Distinct Evolutionary Patterns

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.

Analysis of Arabidopsis floral transcriptome: detection of new florally expressed genes and expansion of Brassicaceae-specific gene families

The flower is essential for sexual reproduction of flowering plants and has been extensively studied. However, it is still not clear how many genes are expressed in the flower. Here, we performed RNA-seq analysis as a highly sensitive approach to investigate the Arabidopsis floral transcriptome at three developmental stages. We provide evidence that at least 23, 961 genes are active in the Arabidopsis flower, including 8512 genes that have not been reported as florally expressed previously. We compared gene expression at different stages and found that many genes encoding transcription factors are preferentially expressed in early flower development. Other genes with expression at distinct developmental stages included DUF577 in meiotic cells and DUF220, DUF1216, and Oleosin in stage 12 flowers. DUF1216 and DUF577 are Brassicaceae specific, and together with other families experienced expansion within the Brassicaceae lineage, suggesting novel/greater roles in Brassicaceae floral development than other plants. The large dataset from this study can serve as a resource for expression analysis of genes involved in flower development in Arabidopsis and for comparison with other species. Together, this work provides clues regarding molecular networks underlying flower development.