Intron loss and gain during evolution of the catalase gene family in angiosperms

Comprehensive Analysis of the Catalase (CAT) Gene Family and Expression Patterns in Rubber Tree (Hevea brasiliensis) under Various Abiotic Stresses and Multiple Hormone Treatments

Catalase (CAT) is one of the key enzymes involved in antioxidant defense systems and mainly scavenges H2O2 and plays a vital role in plant growth, development, and various adverse stresses. To date, a systematic study of the CAT gene family in rubber tree has not been reported. In this study, five HbCAT gene family members were identified from the rubber tree genome, and these were mainly clustered into two subfamilies. Gene structure and motif analysis showed that exon-intron and motif patterns were conserved across different plant species. Sequence analysis revealed that HbCAT proteins contain one active catalytic site, one heme-ligand signature sequence, three conserved amino acid residues (His, Tyr, and Asn), and one peroxisome-targeting signal 1 (PTS1) sequence. Fragment duplication is a selection pressure for the evolution of the HbCAT family based on Ka/Ks values. Analysis of cis-acting elements in the promoters indicated that HbCAT gene expression might be regulated by abscisic acid (ABA), salicylic acid (SA), and MYB transcription factors; furthermore, these genes might be involved in plant growth, development, and abiotic stress responses. A tissue-specific expression analysis showed that HbCATs gradually increased with leaf development and were highly expressed in mature leaves. Gene expression profiling exhibited the differential expression of the HbCATs under cold, heat, drought, and NaCl stresses. Our results provide comprehensive information about the HbCAT gene family, laying the foundation for further research on its function in rubber tree.

Genome-Wide Identification and Expression Analysis of the Sucrose Synthase Gene Family in Sweet Potato and Its Two Diploid Relatives

Sucrose synthases (SUS; EC 2.4.1.13) encoded by a small multigene family are the central system of sucrose metabolism and have important implications for carbon allocation and energy conservation in nonphotosynthetic cells of plants. Though the SUS family genes (SUSs) have been identified in several plants, they have not been explored in sweet potato. In this research, nine, seven and seven SUSs were identified in the cultivated sweet potato (Ipomoea batatas, 2n = 6x = 90) as well as its two diploid wild relatives I. trifida (2n = 2x = 30) and I. triloba (2n = 2x = 30), respectively, and divided into three subgroups according to their phylogenetic relationships. Their protein physicochemical properties, chromosomal localization, phylogenetic relationship, gene structure, promoter cis-elements, protein interaction network and expression patterns were systematically analyzed. The results indicated that the SUS gene family underwent segmental and tandem duplications during its evolution. The SUSs were highly expressed in sink organs. The IbSUSs especially IbSUS2, IbSUS5 and IbSUS7 might play vital roles in storage root development and starch biosynthesis. The SUSs could also respond to drought and salt stress responses and take part in hormone crosstalk. This work provides new insights for further understanding the functions of SUSs and candidate genes for improving yield, starch content, and abiotic stress tolerance in sweet potatoes.

Genome-wide survey of catalase genes in Brassica rapa, Brassica oleracea, and Brassica napus: identification, characterization, molecular evolution, and expression profiling of BnCATs in response to salt and cadmium stress

Catalase (CAT, EC 1.11.1.6), one of the most important antioxidant enzymes, can control excess levels of H₂O₂ produced under oxidative stress in plants. In this study, 16, 8, and 7 CAT genes in the genome of Brassica napus, B. rapa, and B. oleracea were identified, respectively. Phylogenetic studies showed that CATs could be divided into two main groups, each containing specific monocotyledon and dicotyledon subgroups. Motifs, gene structure, and intron phase of CATs in B. napus, Brassica rapa, and Brassica oleracea are highly conserved. Analysis of codon usage bias showed the mutation pressure and natural selection of the codon usage of CATs. Segmental duplication and polyploid were major factors in the expansion of this gene family in B. napus, and genes have experienced  negative selection during evolution. Existence of hormones and stress-responsive cis-elements and identifying miRNA molecules affecting CATs showed that these genes are complexly regulated at the transcriptional and posttranscriptional levels. Based on RNA-seq data, CATs are divided into two groups; the first group has moderate and specific expression in flowers, leaves, stems, and roots, while the second group shows expression in most tissues. qRT-PCR analysis showed that the expression of these genes is dynamic and has a specific expression consistent with other CAT genes in response to salinity and cadmium (Cd) stresses. These results provide information for further investigation of the function of CAT genes in response to stresses and the development of tolerant  plants.

Genome-wide analysis of Catalase gene family reveal insights into abiotic stress response mechanism in Brassica juncea and B. rapa

Environmental stresses affect the yield and productivity of Brassica crops. Catalases are important antioxidant enzymes involved in reducing excess hydrogen peroxide produced by environmental stresses. In the present study, nine and seven CAT family members in two oilseed Brassica species (B. juncea and B. rapa) were identified with complete characterization based on gene and protein structure. Phylogenetic classification categorized CAT proteins into three classes and differentiated the monocot and dicot-specific CAT proteins. Further, the gene and protein characterizations revealed a high degree of conservation across the CAT family members. Differences were observed in the CAT-HEME binding affinity in CAT1, CAT2, and CAT3 isozymes, which could suggest their differential enzyme activities in different conditions. Furthermore, protein-protein interaction with other antioxidant proteins suggested their coordinated role in ROS scavenging mechanisms. Notably, the differential gene expression of BjuCATs and BraCATs and CAT enzyme activities suggested their crucial roles in major abiotic stresses faced by Brassica species. Promoter analysis in BjuCATs and BraCATs suggested the presence of abiotic-stress responsive cis-regulatory elements. Gene regulatory network analysis suggested miRNA and TF mediated stress response in BjuCATs and BraCATs. CAT family screening and characterization in Brassica sp. has established a basic ground for further functional validation in abiotic and heavy-metal stresses which can help in developing stress tolerant crops.

Functional characterization of TBL genes revealed the role of GhTBL7 and GhTBL58 in cotton fiber elongation

TBL (Trichome Birefringence Like) gene family members are involved in trichome initiation and xylan acetylation in several plant species. In our research, we identified 102 TBLs from G. hirsutum. The phylogenetic tree classified TBL genes into five groups. Collinearity analysis of TBL genes indicated 136 paralogous gene pairs in G. hirsutum. Gene duplication indicated that WGD or segmental duplication contributed to the GhTBL gene family expansion. Promoter cis-elements of GhTBLs were related to growth and development, seed-specific regulation, light, and stress responses. GhTBL genes (GhTBL7, GhTBL15, GhTBL21, GhTBL25, GhTBL45, GhTBL54, GhTBL67, GhTBL72, and GhTBL77) exhibited upregulated response under exposure to cold, heat, NaCl, and PEG. GhTBL genes exhibited high expression during fiber development stages. Two GhTBL genes (GhTBL7 and GhTBL58) showed differential expression at 10 DPA fiber, as 10 DPA is a fast fiber elongation stage and fiber elongation is a very important stage of cotton fiber development. Subcellular localization of GhTBL7 and GhTBL58 revealed that these genes reside inside the cell membrane. Promoter GUS activity of GhTBL7 and GhTBL58 exhibited deep staining in roots. To further validate the role of these genes in cotton fiber elongation, we silenced these genes and observed a significant reduction in the fiber length at 10 DPA. In conclusion, the functional study of cell membrane-associated genes (GhTBL7 and GhTBL58) showed deep staining in root tissues and potential function during cotton fiber elongation at 10 DPA fiber.

Comprehensive Genome-Wide Identification and Expression Profiling of Eceriferum (CER) Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant surfaces are covered with cuticle wax and are the first barrier between a plant and environmental stresses. Eceriferum (CER) is an important gene family involved in wax biosynthesis and stress resistance. In this study, for the first time, 34 CER genes were identified in the passion fruit (Passiflora edulis) genome, and PeCER proteins varied in physicochemical properties. A phylogenetic tree was constructed and divided into seven clades to identify the evolutionary relationship with other plant species. Gene structure analyses revealed that conserved motifs ranged from 1 to 24, and that exons ranged from 1 to 29. The cis-element analysis provides insight into possible roles of PeCER genes in plant growth, development and stress responses. The syntenic analysis revealed that segmental (six gene pairs) and tandem (six gene pairs) gene duplication played an important role in the expansion of PeCER genes and underwent a strong purifying selection. In addition, 12 putative ped-miRNAs were identified to be targeting 16 PeCER genes, and PeCER6 was the most targeted by four miRNAs including ped-miR157a-5p, ped-miR164b-5p, ped-miR319b, and ped-miR319l. Potential transcription factors (TFs) such as ERF, AP2, MYB, and bZIP were predicted and visualized in a TF regulatory network interacting with PeCER genes. GO and KEGG annotation analysis revealed that PeCER genes were highly related to fatty acid, cutin, and wax biosynthesis, plant-pathogen interactions, and stress response pathways. The hypothesis that most PeCER proteins were predicted to localize to the plasma membrane was validated by transient expression assays of PeCER32 protein in onion epidermal cells. qRT-PCR expression results showed that most of the PeCER genes including PeCER1, PeCER11, PeCER15, PeCER17, and PeCER32 were upregulated under drought and Fusarium kyushuense stress conditions compared to controls. These findings provide a foundation for further studies on functions of PeCER genes to further facilitate the genetic modification of passion fruit wax biosynthesis and stress resistance.

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.

Catalase (CAT) Gene Family in Wheat (Triticum aestivum L.): Evolution, Expression Pattern and Function Analysis

Catalases (CATs) are present in almost all living organisms and play important roles in plant development and response to various stresses. However, there is relatively little information on CAT genes in wheat and related Triticeae species. A few studies on CAT family genes in wheat have been reported. In this study, ten CAT proteins (TaCATs) were identified in wheat and classified into three groups based on their phylogenetic features and sequence analysis. The analysis of the structure and motif composition of the TaCAT proteins suggested that a segmental duplication event occurred in the TaCAT gene family. Collinearity relationship analysis among different species showed that there were three orthologous CAT genes in rice and in maize. By analyzing the cis-elements in the promoter regions, we speculated that TaCAT genes expression might be regulated by light, oxygen deficit, methyl jasmonate and abscisic acid, and by transcription factors such as MYB. A Gene Ontology (GO)-based analysis showed that TaCAT proteins may be related to the response to various stresses, are cytoplasm localized, and may function as antioxidant enzymes. RT-qPCR and transcriptome data analyses exhibited distinct expression patterns of TaCAT genes in different tissues and in response to various treatments. In this study, a comprehensive analysis of wheat CAT genes was performed, enriching our knowledge of CAT genes and providing a foundation for further functional analyses of this gene family in wheat.

Genome-Wide Identification and Expression Profiling of KCS Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions

Plant and fruit surfaces are covered with cuticle wax and provide a protective barrier against biotic and abiotic stresses. Cuticle wax consists of very-long-chain fatty acids (VLCFAs) and their derivatives. β-Ketoacyl-CoA synthase (KCS) is a key enzyme in the synthesis of VLCFAs and provides a precursor for the synthesis of cuticle wax, but the KCS gene family was yet to be reported in the passion fruit (Passiflora edulis). In this study, thirty-two KCS genes were identified in the passion fruit genome and phylogenetically grouped as KCS1-like, FAE1-like, FDH-like, and CER6-like. Furthermore, thirty-one PeKCS genes were positioned on seven chromosomes, while one PeKCS was localized to the unassembled genomic scaffold. The cis-element analysis provides insight into the possible role of PeKCS genes in phytohormones and stress responses. Syntenic analysis revealed that gene duplication played a crucial role in the expansion of the PeKCS gene family and underwent a strong purifying selection. All PeKCS proteins shared similar 3D structures, and a protein-protein interaction network was predicted with known Arabidopsis proteins. There were twenty putative ped-miRNAs which were also predicted that belong to nine families targeting thirteen PeKCS genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation results were highly associated with fatty acid synthase and elongase activity, lipid metabolism, stress responses, and plant-pathogen interaction. The highly enriched transcription factors (TFs) including ERF, MYB, Dof, C2H2, TCP, LBD, NAC, and bHLH were predicted in PeKCS genes. qRT-PCR expression analysis revealed that most PeKCS genes were highly upregulated in leaves including PeKCS2, PeKCS4, PeKCS8, PeKCS13, and PeKCS9 but not in stem and roots tissues under drought stress conditions compared with controls. Notably, most PeKCS genes were upregulated at 9th dpi under Fusarium kyushuense biotic stress condition compared to controls. This study provides a basis for further understanding the functions of KCS genes, improving wax and VLCFA biosynthesis, and improvement of passion fruit resistance.

Identification and Analysis of GhEXO Gene Family Indicated That GhEXO7_At Promotes Plant Growth and Development Through Brassinosteroid Signaling in Cotton (Gossypium hirsutum L.)

Brassinosteroids (BRs), an efficient plant endogenous hormone, significantly promotes plant nutrient growth adapting to biological and abiotic adversities. BRs mainly promote plant cell elongation by regulating gene expression patterns. EXORDIUM (EXO) genes have been characterized as the indicators of BR response genes. Cotton, an ancient crop, is of great economic value and its fibers can be made into all kinds of fabrics. However, EXO gene family genes have not been full identified in cotton. 175 EXO genes were identified in nine plant species, of which 39 GhEXO genes in Gossypium hirsutum in our study. A phylogenetic analysis grouped all of the proteins encoded by the EXO genes into five major clades. Sequence identification of conserved amino acid residues among monocotyledonous and dicotyledonous species showed a high level of conservation across the N and C terminal regions. Only 25% the GhEXO genes contain introns besides conserved gene structure and protein motifs distribution. The 39 GhEXO genes were unevenly distributed on the 18 At and Dt sub-genome chromosomes. Most of the GhEXO genes were derived from gene duplication events, while only three genes showed evidence of tandem duplication. Homologous locus relationships showed that 15 GhEXO genes are located on collinear blocks and that all orthologous/paralogous gene pairs had Ka > Ks values, indicating purifying selection pressure. The GhEXO genes showed ubiquitous expression in all eight tested cotton tissues and following exposure to three phytohormones, IAA, GA, and BL. Furthermore, GhEXO7_At was mainly expressed in response to BL treatment, and was predominantly expressed in the fibers. GhEXO7_At was found to be a plasma membrane protein, and its ectopic expression in Arabidopsis mediated BR-regulated plant growth and development with altered expression of DWF4, CPD, KCS1, and EXP5. Additionally, the functions of GhEXO7_At were confirmed by virus-induced gene silencing (VIGS) in cotton. This study will provide important genetic resources for future cotton breeding programs.

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.

Molecular characterization revealed the role of catalases under abiotic and arsenic stress in bread wheat (Triticum aestivum L.)

Catalases are crucial antioxidant enzymes that reduce the excessive level of H₂O₂ caused by various environmental stresses and metal toxicity and hence protect the plant cells. In this study, a total of ten TaCAT genes, forming three homeologous groups, were identified in the genome of bread wheat (Triticum aestivum L.) and named as per the wheat gene symbolization guidelines. The identified catalases were characterized for various structural and physicochemical features. The proximal active-site (F(D/A)RERIPERVVHAKGASA) and heme-ligand (R(I/V)F(S/A)Y(A/S)DTQ) signature motifs, catalytic residues and peroxisomal targeting peptides were found conserved. Phylogenetic analysis clustered TaCATs into three classes, which showed conserved functional specialization based on their tissue specific expression. Modulated spatio-temporal expression of various TaCAT genes and alteration in total catalase enzyme activity during heat, drought, salt and arsenic (As^III and As^V) treatment suggested their roles in abiotic stress response and arsenic tolerance. Molecular cloning and overexpression of TaCAT3-B gene in Escherichia coli showed tolerance against heat, drought, salt and varied concentrations of As^III and As^V treatments. The results further confirmed their role in stress tolerance and recommended that these genes can be used in future stress management strategies for the development of abiotic and arsenic stress resistant transgenic crops.

Molecular cloning and expression analysis of sucrose phosphate synthase genes in cassava (Manihot esculenta Crantz)

Sucrose phosphate synthase (SPS), a key rate-limiting enzyme in the sucrose biosynthesis pathway in plants, is encoded by a multi-gene family. Until recently, the identification and characterization of the SPS gene family have been performed for dozens of plant species; however, few studies have involved a comprehensive analysis of the SPS family members in tropical crops, such as cassava (Manihot esculenta Crantz). In the current study, five SPS genes (MeSPS1, MeSPS2, MeSPS3, MeSPS4, and MeSPS5) were isolated from cassava, and their sequence characteristics were comprehensively characterized. These MeSPS genes were found distributed on five chromosomes (Chr2, Chr14, Chr15, Chr16, and Chr18). Phylogenetic analysis showed that the MeSPS protein sequences were clustered into three families, together with other SPS sequences from both dicot and monocot species (families A, B, and C). The spatio-temporal expression pattern analysis of MeSPS genes showed a tissue-specific and partially overlapping expression pattern, with the genes mainly expressed in source tissues during cassava growth and development. Correlation analysis revealed that the expression of MeSPS genes correlated positively with root starch content, indicating that the expression of MeSPS genes might accelerate the rate of starch accumulation in the roots of cassava plants.

Programmed cell death (PCD) control in plants: New insights from the Arabidopsis thaliana deathosome

Programmed cell death (PCD) is a genetically controlled process that leads to cell suicide in both eukaryotic and prokaryotic organisms. In plants PCD occurs during development, defence response and when exposed to adverse conditions. PCD acts controlling the number of cells by eliminating damaged, old, or unnecessary cells to maintain cellular homeostasis. Unlike in animals, the knowledge about PCD in plants is limited. The molecular network that controls plant PCD is poorly understood. Here we present a review of the current mechanisms involved with the genetic control of PCD in plants. We also present an updated version of the AtLSD1 deathosome, which was previously proposed as a network controlling HR-mediated cell death in Arabidopsis thaliana. Finally, we discuss the unclear points and open questions related to the AtLSD1 deathosome.

Genome-wide identification and analysis of the sucrose synthase gene family in cassava (Manihot esculenta Crantz)

Sucrose synthase (SUS), a key enzyme of the sucrose metabolism pathway, is encoded by a multi-gene family in plants. To date, dozens of SUS gene families have been characterized in various plant genomes. However, only a few studies have performed comprehensive analyses in tropical crops like cassava (Manihot esculenta Crantz). In the present study, seven non-redundant members of the SUS gene family (MeSUS1-7) were identified and characterized from the cassava genome. The MeSUS genes were distributed on five chromosomes (Chr1, Chr2, Chr3, Chr14, and Chr16) and the encoded proteins could be classified into three major groups with other SUS proteins from both dicot and monocot species (SUS I, SUS II, and SUS III). The spatio-temporal expression profiles of MeSUS genes showed a developmental stage-dependent, partially overlapping pattern, mainly expressed in the source and sink tissues. Cold and drought treatments significantly induced the expressions of MeSUS2, MeSUS4, MeSUS6, and MeSUS7 and the activities of the encoded enzymes, indicating that these genes may play crucial roles in resistance against abiotic stresses. These results provide new insights into the multifaceted role of the SUS gene family members in various physiological processes, especially sucrose transport and starch accumulation in cassava roots.

Plant catalases as NO and H2S targets

Catalase is a powerful antioxidant metalloenzyme located in peroxisomes which also plays a central role in signaling processes under physiological and adverse situations. Whereas animals contain a single catalase gene, in plants this enzyme is encoded by a multigene family providing multiple isoenzymes whose number varies depending on the species, and their expression is regulated according to their tissue/organ distribution and the environmental conditions. This enzyme can be modulated by reactive oxygen and nitrogen species (ROS/RNS) as well as by hydrogen sulfide (H2S). Catalase is the major protein undergoing Tyr-nitration [post-translational modification (PTM) promoted by RNS] during fruit ripening, but the enzyme from diverse sources is also susceptible to undergo other activity-modifying PTMs. Data on S-nitrosation and persulfidation of catalase from different plant origins are given and compared here with results from obese children where S-nitrosation of catalase occurs. The cysteine residues prone to be S-nitrosated in catalase from plants and from bovine liver have been identified. These evidences assign to peroxisomes a crucial statement in the signaling crossroads among relevant molecules (NO and H2S), since catalase is allocated in these organelles. This review depicts a scenario where the regulation of catalase through PTMs, especially S-nitrosation and persulfidation, is highlighted.

Genome-wide expression analysis suggests glutaredoxin genes response to various stresses in cotton

Oxidative stress reflects an imbalance between the systemic manifestation of reactive oxygen species (ROS) and a biological system's ability to readily detoxify the reactive intermediates or to repair the resulting damage. Glutaredoxins (GRXs) are ubiquitous oxidoreductase enzymes involved in diverse cellular processes and play a key role in oxidative stress responsive mechanisms. This study was aimed to explore the structure-function relationship and to provide a framework for functional validation and biochemical characterization of various GRX members. In this study, our analysis revealed the presence of 127 genes encoding GRX proteins in G. hirsutum. A total of 758 genes from two typical monocot and nine dicot species were naturally divided into four classes based on phylogenetic analysis. The classification was supported with organization of conserved protein motifs and sequence logos comparison between cotton, rice and Arabidopsis. Cotton GRX gene family has underwent strong purifying selection with limited functional divergence. A good collinearity was observed in the synteny analysis of four Gossypium species. Majority of cotton GRXs were influenced by various phytohormones and abiotic stress conditions during expression analysis, suggesting an important role of GRX proteins in response to oxidative stress. Cis-regulatory elements, gene enrichments and co-expression network analysis also support their predicted role against various abiotic stresses. Whole genome and segmental duplication were determined to be the two major impetuses for the expansion of gene numbers during the evolution. The identification of GRX genes showing differential expression in specific tissues or in response to environmental stimuli provides a new avenue for in-depth characterization of selected genes of importance. This study will further broaden our insights into the evolution and functional elucidation of GRX gene family in cotton.

Genome-wide identification of GhAAI genes reveals that GhAAI66 triggers a phase transition to induce early flowering

Plants undergo a phase transition from vegetative to reproductive development that triggers floral induction. Genes containing an AAI (α-amylase inhibitor) domain form a large gene family, but there have been no comprehensive analyses of this gene family in any plant species. Here, we identified 336 AAI genes from nine plant species including122 AAI genes in cotton (Gossypium hirsutum). The AAI gene family has evolutionarily conserved amino acid residues throughout the plant kingdom. Phylogenetic analysis classified AAI genes into five major clades with significant polyploidization and showing effects of genome duplication. Our study identified 42 paralogous and 216 orthologous gene pairs resulting from segmental and whole-genome duplication, respectively, demonstrating significant contributions of gene duplication to expansion of the cotton AAI gene family. Further, GhAAI66 was preferentially expressed in flower tissue and as responses to phytohormone treatments. Ectopic expression of GhAAI66 in Arabidopsis and silencing in cotton revealed that GhAAI66 triggers a phase transition to induce early flowering. Further, GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of RNA sequencing data and qRT-PCR (quantitative reverse transcription-PCR) analysis indicated that GhAAI66 integrates multiple flower signaling pathways including gibberellin, jasmonic acid, and floral integrators to trigger an early flowering cascade in Arabidopsis. Therefore, characterization of the AAI family provides invaluable insights for improving cotton breeding.

The evolution of the plastid phosphate translocator family

The plastid phosphate translocators evolved in algae but diversified into several groups, which adopted different physiological functions by extensive gene duplications and losses in Streptophyta. The plastid phosphate translocators (pPT) are a family of transporters involved in the exchange of metabolites and inorganic phosphate between stroma and cytosol. Based on their substrate specificities, they were divided into four subfamilies named TPT, PPT, GPT and XPT. To analyse the occurrence of these transporters in different algae and land plant species, we identified 652 pPT genes in 101 sequenced genomes for phylogenetic analysis. The first three subfamilies are found in all species and evolved before the split of red and green algae while the XPTs were derived from the duplication of a GPT gene at the base of Streptophyta. The analysis of the intron-exon structures of the pPTs corroborated these findings. While the number and positions of introns are conserved within each subfamily, they differ between the subfamilies suggesting an insertion of the introns shortly after the three subfamilies evolved. During angiosperm evolution, the subfamilies further split into different groups (TPT1-2, PPT1-3, GPT1-6). Angiosperm species differ significantly in the total number of pPTs, with many species having only a few, while several plants, especially crops, have a higher number, pointing to the importance of these transporters for improved source-sink strength and yield. The differences in the number of pPTs can be explained by several small-scale gene duplications and losses in plant families or single species, but also by whole genome duplications, for example, in grasses. This work could be the basis for a comprehensive analysis of the molecular and physiological functions of this important family of transporters.

Plant peroxisomes at the crossroad of NO and H2 O2 metabolism

Plant peroxisomes are subcellular compartments involved in many biochemical pathways during the life cycle of a plant but also in the mechanism of response against adverse environmental conditions. These organelles have an active nitro-oxidative metabolism under physiological conditions but this could be exacerbated under stress situations. Furthermore, peroxisomes have the capacity to proliferate and also undergo biochemical adaptations depending on the surrounding cellular status. An important characteristic of peroxisomes is that they have a dynamic metabolism of reactive nitrogen and oxygen species (RNS and ROS) which generates two key molecules, nitric oxide (NO) and hydrogen peroxide (H₂ O₂ ). These molecules can exert signaling functions by means of post-translational modifications that affect the functionality of target molecules like proteins, peptides or fatty acids. This review provides an overview of the endogenous metabolism of ROS and RNS in peroxisomes with special emphasis on polyamine and uric acid metabolism as well as the possibility that these organelles could be a source of signal molecules involved in the functional interconnection with other subcellular compartments.

Genome-Wide Identification and Characterization of the PERK Gene Family in Gossypium hirsutum Reveals Gene Duplication and Functional Divergence

Proline-rich extensin-like receptor kinases (PERKs) are an important class of receptor kinases in plants. Receptor kinases comprise large gene families in many plant species, including the 15 PERK genes in Arabidopsis. At present, there is no comprehensive published study of PERK genes in G. hirsutum. Our study identified 33 PERK genes in G. hirsutum. Phylogenetic analysis of conserved PERK protein sequences from 15 plant species grouped them into four well defined clades. The GhPERK gene family is an evolutionarily advanced gene family that lost its introns over time. Several cis-elements were identified in the promoter regions of the GhPERK genes that are important in regulating growth, development, light responses and the response to several stresses. In addition, we found evidence for gene loss or addition through segmental or whole genome duplication in cotton. Gene duplication and synteny analysis identified 149 orthologous/paralogous gene pairs. Ka/Ks values show that most GhPERK genes experienced strong purifying selection during the rapid evolution of the gene family. GhPERK genes showed high expression levels in leaves and during ovule development. Furthermore, the expression of GhPERK genes can be regulated by abiotic stresses and phytohormone treatments. Additionally, PERK genes could be involved in several molecular, biological and physiological processes that might be the result of functional divergence.

Cyanobacterial Catalase Activity Prevents Oxidative Stress Induced by Pseudomonas fluorescens DUS1-27 from Inhibiting Brassica napus L. (canola) Growth

Plant growth-promoting bacteria (PGPB) inhabit the rhizosphere of plants and are capable of enhancing plant growth through a number of mechanisms. A strain of Pseudomonas fluorescens DUS1-27 was identified as a potential PGPB candidate based on its ability to increase the growth of Brassica napus L. (canola) over that of uninoculated control plants in a soil-based system. The same P. fluorescens isolate was found to reduce plant growth in a hydroponic growth system, with plants showing the symptoms of a microbe-associated molecular pattern (MAMP) response to the bacteria. The amperometric quantification of H2O2, fluorescence-based total peroxidase assays, and quantification of catalase gene expression levels using qRT-PCR revealed that oxidative stress reduced plant growth in the hydroponic system. The addition of the cyanobacterium Nostoc punctiforme (known to have high catalase activity levels) in the hydroponic system as a co-inoculant reduced oxidative stress (49.7% decrease in H2O2 concentrations) triggered by the addition of P. fluorescens DUS1-27, thereby enabling plants to grow larger than uninoculated control plants. These results show the advantage of inoculating with multiple bacteria to promote plant growth and, for the first time, demonstrate that N. punctiforme beneficially assists plants under oxidative stress through its catalase activity in planta.

Genome-Wide Identification of PIFs in Grapes (Vitis vinifera L.) and Their Transcriptional Analysis under Lighting/Shading Conditions

Phytochrome-interacting factors (PIFs), as the basic helix⁻loop⁻helix (bHLH) transcription factors, are the primary signaling partners for phytochromes (PHY) that play a key role in PHY-mediated light signal transduction. At present, there are few studies on PIFs in fruit trees. In order to clarify the status of PIFs in grapevines, we identified members of the grape PIFs family and conducted phylogenetic and expression analysis. We identified PIF1, PIF3, PIF4, and PIF7 in PIFs families of the grapevine (Vitis vinifera L.), which were distributed on four different chromosomes with similar gene structures. Except for the closer relationship with PIF1 of citrus, PIFs of grape were distant from the other fruit species such as apple, pear, peach, and strawberry. The VvPIFs (except VvPIF4) were located in the syntenic block with those from Arabidopsisthaliana, Solanum lycopersicum, or Citrus sinensis. In addition to PIF1, all PIFs in grapevines have conserved active PHYB binding (APB) sequences. VvPIF1 has a conserved PIF1-specific active PHYA binding (APA) sequence, while amino acid mutations occurred in the specific APA sequence in VvPIF3. Interestingly, two specific motifs were found in the PIF4 amino acid sequence. The photoreceptor-related elements in the VvPIFs promoter region were the most abundant. PIF1, LONG HYPOCOTYL 5 (HY5) and PIF3, PIF4, GIBBERELLIC ACID INSENSITIVE 1 (GAI1) may interact with each other and participate together in light signal transduction. The relative expression levels of the VvPIFs showed diverse patterns in the various organs at different developmental stages, of which PIF4 was most highly expressed. Prior to maturation, the expression of PIF4 and PIF7 in the skin of the different cultivars increased, while the expression of all PIFs in the flesh decreased. The transcription level of PIFs in grape leaves was sensitive to changes in lighting and shading. Shading treatment was beneficial for enhancing the transcription level of VvPIFs, but the effect on VvPIF3 and VvPIF4 was time-controlled. We concluded that PIFs in grapevines are both conservative and species-specific. The identification and analysis of grape PIFs could provide a theoretical foundation for the further construction of grape light regulation networks.

Comprehensive analysis and transcript profiling of Arabidopsis thaliana and Oryza sativa catalase gene family suggests their specific roles in development and stress responses

Stress induces the generation of Reactive Oxygen Species (ROS) that ultimately hampers the growth, development, and productivity of the plant. As an antioxidant enzyme, catalase converts hydrogen peroxide to water and keeps ROS level down to protect cells from stress-induced apoptosis. Here, a genome-wide analysis of catalase gene family has been performed in two model plants- Arabidopsis thaliana and Oryza sativa. Both Arabidopsis and rice has a small family of three and four genes, respectively; that code for seven proteins each. Detailed analysis of these members in terms of their structure, duplication, chromosomal position and proteins subcellular localization, as well as expression profiling under various developmental and environmental cues, was performed. Catalase proteins were mostly found to be localized in the cytoplasm, followed by peroxisome and mitochondria. Phylogenetically plant catalases showed strong divergence from their non-plant counterparts. Expression profiling revealed that AtCAT3 and OsCATA are the constitutively expressive member; while AtCAT2, OsCATA, and OsCATC are the stress-responsive members. Moreover, an altered level of total rice catalase enzyme activity and H₂O₂ level was observed under various abiotic stress conditions. This indicates the stress-responsive transcriptome as well as proteome alteration of catalase in the plant.

Genetic variation of bread wheat accessions in response to the cereal cyst nematode, Heterodera filipjevi

Bread wheat, Triticum aestivum, produces large edible grains and is consumed by 75% of the world’s populations. Cereal cyst nematodes have a global distribution and cause significant economic yield losses in many countries. Therefore, there is an urgent need to identify new resistance sources. In this study, the genetic diversity of 43 wheat accessions in response to cereal cyst nematode, Heterodera filipjevi Isfahan pathotype, was assessed using a simple sequence repeat (SSR) marker. Seven primers were used, out of which five primers showed polymorphisms. Alleles per primer varied from one to three per locus (mean 2.85). The highest and lowest polymorphic information content of 0.81 and 0.44 (mean 0.66) were related to Xgwm 3012DL and Xgwm147, respectively. Genetic similarity was 29-88% between accessions. SSR analysis divided the accessions into five main groups. Resistant cultivars ‘Bam’ and ‘Behrang’ possessed both Cre1 and Cre8 resistant genes. The Cre3 and Cat genes were partially sequenced in five cultivars of different responses to H. filipjevi. The nucleotide sequences were compared to Cre3 and Cat homologues, indicating 93-100% and 86-92% homology, respectively. The MEGA program showed highest similarity of Cre3 and Cat genes amplified with the resistance gene analogues (RGA14) in the wheat and Cat3-A1 gene in ‘Carnamah’. This research showed that SRR markers could efficiently verify genetic diversity between wheat accessions, and the known resistance genes (Cre genes) against the cereal cyst nematodes could not control the H. filipjevi Isfahan pathotype populations, except the Cre1 gene.

Parallel loss of introns in the ABCB1 gene in angiosperms

BACKGROUND

The presence of non-coding introns is a characteristic feature of most eukaryotic genes. While the size of the introns, number of introns per gene and the number of intron-containing genes can vary greatly between sequenced eukaryotic genomes, the structure of a gene with reference to intron presence and positions is typically conserved in closely related species. Unexpectedly, the ABCB1 (ATP-Binding Cassette Subfamily B Member 1) gene which encodes a P-glycoprotein and underlies dwarfing traits in maize (br2), sorghum (dw3) and pearl millet (d2) displayed considerable variation in intron composition.

RESULTS

An analysis of the ABCB1 gene structure in 80 angiosperms revealed that the number of introns ranged from one to nine. All introns in ABCB1 underwent either a one-time loss (single loss in one lineage/species) or multiple independent losses (parallel loss in two or more lineages/species) with the majority of losses occurring within the grass family. In contrast, the structure of the closest homolog to ABCB1, ABCB19, remained constant in the majority of angiosperms analyzed. Using known phylogenetic relationships within the grasses, we determined the ancestral branch-points where the losses occurred. Intron 7, the longest intron, was lost in only a single species, Mimulus guttatus, following duplication of ABCB1. Semiquantitative PCR showed that the M. guttatus ABCB1 gene copy without intron 7 had significantly lower transcript levels than the gene copy with intron 7. We further demonstrated that intron 7 carried two motifs that were highly conserved across the monocot-dicot divide.

CONCLUSIONS

The ABCB1 gene structure is highly dynamic, while the structure of ABCB19 remained largely conserved through evolution. Precise removal of introns, preferential removal of smaller introns and presence of at least 2 bp of microhomology flanking most introns indicated that intron loss may have predominantly occurred through non-homologous end-joining (NHEJ) repair of double strand breaks. Lack of microhomology in the exon upstream of lost phase I introns was likely due to release of the selective constraint on the penultimate base (3rd base in codon) of the terminal codon by the splicing machinery. In addition to size, the presence of regulatory motifs will make introns recalcitrant to loss.

Analysis of aquaporins from the euryhaline barnacle Balanus improvisus reveals differential expression in response to changes in salinity

Barnacles are sessile macro-invertebrates, found along rocky shores in coastal areas worldwide. The euryhaline bay barnacle Balanus improvisus (Darwin, 1854) (= Amphibalanus improvisus) can tolerate a wide range of salinities, but the molecular mechanisms underlying the osmoregulatory capacity of this truly brackish species are not well understood. Aquaporins are pore-forming integral membrane proteins that facilitate transport of water, small solutes and ions through cellular membranes, and that have been shown to be important for osmoregulation in many organisms. The knowledge of the function of aquaporins in crustaceans is, however, limited and nothing is known about them in barnacles. We here present the repertoire of aquaporins from a thecostracan crustacean, the barnacle B. improvisus, based on genome and transcriptome sequencing. Our analyses reveal that B. improvisus contains eight genes for aquaporins. Phylogenetic analysis showed that they represented members of the classical water aquaporins (Aqp1, Aqp2), the aquaglyceroporins (Glp1, Glp2), the unorthodox aquaporin (Aqp12) and the arthropod-specific big brain aquaporin (Bib). Interestingly, we also found two big brain-like proteins (BibL1 and BibL2) constituting a new group of aquaporins not yet described in arthropods. In addition, we found that the two water-specific aquaporins were expressed as C-terminal splice variants. Heterologous expression of some of the aquaporins followed by functional characterization showed that Aqp1 transported water and Glp2 water and glycerol, agreeing with the predictions of substrate specificity based on 3D modeling and phylogeny. To investigate a possible role for the B. improvisus aquaporins in osmoregulation, mRNA expression changes in adult barnacles were analysed after long-term acclimation to different salinities. The most pronounced expression difference was seen for AQP1 with a substantial (>100-fold) decrease in the mantle tissue in low salinity (3 PSU) compared to high salinity (33 PSU). Our study provides a base for future mechanistic studies on the role of aquaporins in osmoregulation.

Genome-Wide Analysis of the Sucrose Synthase Gene Family in Grape (Vitis vinifera): Structure, Evolution, and Expression Profiles

Sucrose synthase (SS) is widely considered as the key enzyme involved in the plant sugar metabolism that is critical to plant growth and development, especially quality of the fruit. The members of SS gene family have been identified and characterized in multiple plant genomes. However, detailed information about this gene family is lacking in grapevine (Vitis vinifera L.). In this study, we performed a systematic analysis of the grape (V. vinifera) genome and reported that there are five SS genes (VvSS1-5) in the grape genome. Comparison of the structures of grape SS genes showed high structural conservation of grape SS genes, resulting from the selection pressures during the evolutionary process. The segmental duplication of grape SS genes contributed to this gene family expansion. The syntenic analyses between grape and soybean (Glycine max) demonstrated that these genes located in corresponding syntenic blocks arose before the divergence of grape and soybean. Phylogenetic analysis revealed distinct evolutionary paths for the grape SS genes. VvSS1/VvSS5, VvSS2/VvSS3 and VvSS4 originated from three ancient SS genes, which were generated by duplication events before the split of monocots and eudicots. Bioinformatics analysis of publicly available microarray data, which was validated by quantitative real-time reverse transcription PCR (qRT-PCR), revealed distinct temporal and spatial expression patterns of VvSS genes in various tissues, organs and developmental stages, as well as in response to biotic and abiotic stresses. Taken together, our results will be beneficial for further investigations into the functions of SS gene in the processes of grape resistance to environmental stresses.

Phylogenetic Analysis, Lineage-Specific Expansion and Functional Divergence of seed dormancy 4-Like Genes in Plants

The rice gene seed dormancy 4 (OsSdr4) functions in seed dormancy and is a major factor associated with pre-harvest sprouting (PHS). Although previous studies of this protein family were reported for rice and other species, knowledge of the evolution of genes homologous to OsSdr4 in plants remains inadequate. Fifty four Sdr4-like (hereafter designated Sdr4L) genes were identified in nine plant lineages including 36 species. Phylogenetic analysis placed these genes in eight subfamilies (I-VIII). Genes from the same lineage clustered together, supported by analysis of conserved motifs and exon-intron patterns. Segmental duplications were present in both dicot and monocot clusters, while tandemly duplicated genes occurred only in monocot clusters indicating that both tandem and segmental duplications contributed to expansion of the grass I and II subfamilies. Estimation of the approximate ages of the duplication events indicated that ancestral Sdr4 genes evolved from a common angiosperm ancestor, about 160 million years ago (MYA). Moreover, diversification of Sdr4L genes in mono and dicot plants was mainly associated with genome-wide duplication and speciation events. Functional divergence was observed in all subfamily pairs, except IV/VIIIa. Further analysis indicated that functional constraints between subfamily pairs I/II, I/VIIIb, II/VI, II/VIIIb, II/IV, and VI/VIIIb were statistically significant. Site and branch-site model analyses of positive selection suggested that these genes were under strong adaptive selection pressure. Critical amino acids detected for both functional divergence and positive selection were mostly located in the loops, pointing to functional importance of these regions in this protein family. In addition, differential expression studies by transcriptome atlas of 11 Sdr4L genes showed that the duplicated genes may have undergone divergence in expression between plant species. Our findings showed that Sdr4L genes are functionally divergent and positively selected. These may contribute to further functional analysis and molecular evolution of Sdr4L gene families in land plants.

Identification of a member of the catalase multigene family on wheat chromosome 7A associated with flour b* colour and biological significance of allelic variation

Carotenoids (especially lutein) are known to be the pigment source for flour b* colour in bread wheat. Flour b* colour variation is controlled by a quantitative trait locus (QTL) on wheat chromosome 7AL and one gene from the carotenoid pathway, phytoene synthase, was functionally associated with the QTL on 7AL in some, but not all, wheat genotypes. A SNP marker within a sequence similar to catalase (Cat3-A1snp) derived from full-length (FL) cDNA (AK332460), however, was consistently associated with the QTL on 7AL and implicated in regulating hydrogen peroxide (H2O2) to control carotenoid accumulation affecting flour b* colour. The number of catalase genes on chromosome 7AL was investigated in this study to identify which gene may be implicated in flour b* variation and two were identified through interrogation of the draft wheat genome survey sequence consisting of five exons and a further two members having eight exons identified through comparative analysis with the single catalase gene on rice chromosome 6, PCR amplification and sequencing. It was evident that the catalase genes on chromosome 7A had duplicated and diverged during evolution relative to its counterpart on rice chromosome 6. The detection of transcripts in seeds, the co-location with Cat3-A1snp marker and maximised alignment of FL-cDNA (AK332460) with cognate genomic sequence indicated that TaCat3-A1 was the member of the catalase gene family associated with flour b* colour variation. Re-sequencing identified three alleles from three wheat varieties, TaCat3-A1a, TaCat3-A1b and TaCat3-A1c, and their predicted protein identified differences in peroxisomal targeting signal tri-peptide domain in the carboxyl terminal end providing new insights into their potential role in regulating cellular H2O2 that contribute to flour b* colour variation.

Analysis of the sucrose synthase gene family in tobacco: structure, phylogeny, and expression patterns

MAIN CONCLUSION

Provide an evolutionary and an empirical molecular genetic foundation of the Sus gene family in tobacco and will be beneficial for further investigations of Sus gene functions Sucrose synthase (Sus) has been well characterized as the key enzyme participating in sucrose metabolism, and the gene family encoding different Sus isozymes has been cloned and characterized in several plant species. However, scant information about this gene family is available to date in tobacco. Here, we identified 14, 6, and 7 Sus genes in the genomes of Nicotiana tabacum, N. sylvestris and N. tomentosiformis, respectively. These tobacco Sus family members shared high levels of similarity in their nucleotide and amino acid sequences. Phylogenetic analysis revealed distinct evolutionary paths for the tobacco Sus genes. Sus1-4, Sus5, and Sus6-7 originated from three Sus precursors, respectively, which were generated by duplication before the split of monocots and eudicots. There were two additional duplications, before and after the differentiation of the Solanaceae, which separately gave rise to Sus3/4 and Sus1/2. Gene exon/intron structure analysis showed that the tobacco Sus genes contain varying numbers of conserved introns, resulting from intron loss under different selection pressures during the course of evolution. The expression patterns of the NtSus genes differed from each other in various tobacco tissues. Transcripts of Ntab0259170 and Ntab0259180 were detected in leaves at all tested developmental stages, suggesting that these two genes play a predominant role in sucrose metabolism during leaf development. Expression of Ntab0288750 and Ntab0234340 were conspicuously induced by low temperature and virus treatment, indicating that these two isozymes are important in meeting the increased glycolytic demand that occurs during abiotic stress. Our results provide an evolutionary and an empirical molecular genetic foundation of the Sus gene family in tobacco, and will be beneficial for further investigations of Sus gene functions in the processes of tobacco leaf development and tobacco resistance to environmental stresses.

Reannotation and extended community resources for the genome of the non-seed plant Physcomitrella patens provide insights into the evolution of plant gene structures and functions

Background

The moss Physcomitrella patens as a model species provides an important reference for early-diverging lineages of plants and the release of the genome in 2008 opened the doors to genome-wide studies. The usability of a reference genome greatly depends on the quality of the annotation and the availability of centralized community resources. Therefore, in the light of accumulating evidence for missing genes, fragmentary gene structures, false annotations and a low rate of functional annotations on the original release, we decided to improve the moss genome annotation.

Results

Here, we report the complete moss genome re-annotation (designated V1.6) incorporating the increased transcript availability from a multitude of developmental stages and tissue types. We demonstrate the utility of the improved P. patens genome annotation for comparative genomics and new extensions to the cosmoss.org resource as a central repository for this plant “flagship” genome. The structural annotation of 32,275 protein-coding genes results in 8387 additional loci including 1456 loci with known protein domains or homologs in Plantae. This is the first release to include information on transcript isoforms, suggesting alternative splicing events for at least 10.8% of the loci. Furthermore, this release now also provides information on non-protein-coding loci. Functional annotations were improved regarding quality and coverage, resulting in 58% annotated loci (previously: 41%) that comprise also 7200 additional loci with GO annotations. Access and manual curation of the functional and structural genome annotation is provided via the http://www.cosmoss.org model organism database.

Conclusions

Comparative analysis of gene structure evolution along the green plant lineage provides novel insights, such as a comparatively high number of loci with 5’-UTR introns in the moss. Comparative analysis of functional annotations reveals expansions of moss house-keeping and metabolic genes and further possibly adaptive, lineage-specific expansions and gains including at least 13% orphan genes.

Genome-wide analysis of the Sus gene family in cotton

Sucrose synthase (Sus) is a key enzyme in plant sucrose metabolism. In cotton, Sus (EC 2.4.1.13) is the main enzyme that degrades sucrose imported into cotton fibers from the phloem of the seed coat. This study demonstrated that the genomes of Gossypium arboreum L., G. raimondii Ulbr., and G. hirsutum L., contained 8, 8, and 15 Sus genes, respectively. Their structural organizations, phylogenetic relationships, and expression profiles were characterized. Comparisons of genomic and coding sequences identified multiple introns, the number and positions of which were highly conserved between diploid and allotetraploid cotton species. Most of the phylogenetic clades contained sequences from all three species, suggesting that the Sus genes of tetraploid G. hirsutum derived from those of its diploid ancestors. One Sus group (Sus I) underwent expansion during cotton evolution. Expression analyses indicated that most Sus genes were differentially expressed in various tissues and had development-dependent expression profiles in cotton fiber cells. Members of the same orthologous group had very similar expression patterns in all three species. These results provide new insights into the evolution of the cotton Sus gene family, and insight into its members' physiological functions during fiber growth and development.

Analyses of the sucrose synthase gene family in cotton: structure, phylogeny and expression patterns

BACKGROUND

In plants, sucrose synthase (Sus) is widely considered as a key enzyme involved in sucrose metabolism. Several paralogous genes encoding different isozymes of Sus have been identified and characterized in multiple plant genomes, while limited information of Sus genes is available to date for cotton.

RESULTS

Here, we report the molecular cloning, structural organization, phylogenetic evolution and expression profiles of seven Sus genes (GaSus1 to 7) identified from diploid fiber cotton (Gossypium arboreum). Comparisons between cDNA and genomic sequences revealed that the cotton GaSus genes were interrupted by multiple introns. Comparative screening of introns in homologous genes demonstrated that the number and position of Sus introns are highly conserved among Sus genes in cotton and other more distantly related plant species. Phylogenetic analysis showed that GaSus1, GaSus2, GaSus3, GaSus4 and GaSus5 could be clustered together into a dicot Sus group, while GaSus6 and GaSus7 were separated evenly into other two groups, with members from both dicot and monocot species. Expression profiles analyses of the seven Sus genes indicated that except GaSus2, of which the transcripts was undetectable in all tissues examined, and GaSus7, which was only expressed in stem and petal, the other five paralogues were differentially expressed in a wide ranges of tissues, and showed development-dependent expression profiles in cotton fiber cells.

CONCLUSIONS

This is a comprehensive study of the Sus gene family in cotton plant. The results presented in this work provide new insights into the evolutionary conservation and sub-functional divergence of the cotton Sus gene family in response to cotton fiber growth and development.

Deciphering the genomic structure, function and evolution of carotenogenesis related phytoene synthases in grasses

BACKGROUND

Carotenoids are isoprenoid pigments, essential for photosynthesis and photoprotection in plants. The enzyme phytoene synthase (PSY) plays an essential role in mediating condensation of two geranylgeranyl diphosphate molecules, the first committed step in carotenogenesis. PSY are nuclear enzymes encoded by a small gene family consisting of three paralogous genes (PSY1-3) that have been widely characterized in rice, maize and sorghum.

RESULTS

In wheat, for which yellow pigment content is extremely important for flour colour, only PSY1 has been extensively studied because of its association with QTLs reported for yellow pigment whereas PSY2 has been partially characterized. Here, we report the isolation of bread wheat PSY3 genes from a Renan BAC library using Brachypodium as a model genome for the Triticeae to develop Conserved Orthologous Set markers prior to gene cloning and sequencing. Wheat PSY3 homoeologous genes were sequenced and annotated, unravelling their novel structure associated with intron-loss events and consequent exonic fusions. A wheat PSY3 promoter region was also investigated for the presence of cis-acting elements involved in the response to abscisic acid (ABA), since carotenoids also play an important role as precursors of signalling molecules devoted to plant development and biotic/abiotic stress responses. Expression of wheat PSYs in leaves and roots was investigated during ABA treatment to confirm the up-regulation of PSY3 during abiotic stress.

CONCLUSIONS

We investigated the structural and functional determinisms of PSY genes in wheat. More generally, among eudicots and monocots, the PSY gene family was found to be associated with differences in gene copy numbers, allowing us to propose an evolutionary model for the entire PSY gene family in Grasses.

Phylogenetic distribution of intron positions in alpha-amylase genes of bilateria suggests numerous gains and losses

Most eukaryotes have at least some genes interrupted by introns. While it is well accepted that introns were already present at moderate density in the last eukaryote common ancestor, the conspicuous diversity of intron density among genomes suggests a complex evolutionary history, with marked differences between phyla. The question of the rates of intron gains and loss in the course of evolution and factors influencing them remains controversial. We have investigated a single gene family, alpha-amylase, in 55 species covering a variety of animal phyla. Comparison of intron positions across phyla suggests a complex history, with a likely ancestral intronless gene undergoing frequent intron loss and gain, leading to extant intron/exon structures that are highly variable, even among species from the same phylum. Because introns are known to play no regulatory role in this gene and there is no alternative splicing, the structural differences may be interpreted more easily: intron positions, sizes, losses or gains may be more likely related to factors linked to splicing mechanisms and requirements, and to recognition of introns and exons, or to more extrinsic factors, such as life cycle and population size. We have shown that intron losses outnumbered gains in recent periods, but that "resets" of intron positions occurred at the origin of several phyla, including vertebrates. Rates of gain and loss appear to be positively correlated. No phase preference was found. We also found evidence for parallel gains and for intron sliding. Presence of introns at given positions was correlated to a strong protosplice consensus sequence AG/G, which was much weaker in the absence of intron. In contrast, recent intron insertions were not associated with a specific sequence. In animal Amy genes, population size and generation time seem to have played only minor roles in shaping gene structures.

Size polymorphism in alleles of the myoglobin gene from biomphalaria mollusks

Introns are common among all eukaryotes, while only a limited number of introns are found in prokaryotes. Globin and globin-like proteins are widely distributed in nature, being found even in prokaryotes and a wide range of patterns of intron-exon have been reported in several eukaryotic globin genes. Globin genes in invertebrates show considerable variation in the positions of introns; globins can be found without introns, with only one intron or with three introns in different positions. In this work we analyzed the introns in the myoglobin gene from Biomphalaria glabrata, B. straminea and B. tenagophila. In the Biomphalaria genus, the myoglobin gene has three introns; these were amplified by PCR and analyzed by PCR-RFLP. Results showed that the size (number or nucleotides) and the nucleotide sequence of the coding gene of the myoglobin are variable in the three species. We observed the presence of size polymorphisms in intron 2 and 3; this characterizes a homozygous/heterozygous profile and it indicates the existence of two alleles which are different in size in each species of Biomphalaria. This polymorphism could be explored for specific identification of Biomphalaria individuals.

Genome-wide analysis of the chalcone synthase superfamily genes of Physcomitrella patens

Enzymes of the chalcone synthase (CHS) superfamily catalyze the production of a variety of secondary metabolites in bacteria, fungi and plants. Some of these metabolites have played important roles during the early evolution of land plants by providing protection from various environmental assaults including UV irradiation. The genome of the moss, Physcomitrella patens, contains at least 17 putative CHS superfamily genes. Three of these genes (PpCHS2b, PpCHS3 and PpCHS5) exist in multiple copies and all have corresponding ESTs. PpCHS11 and probably also PpCHS9 encode non-CHS enzymes, while PpCHS10 appears to be an ortholog of plant genes encoding anther-specific CHS-like enzymes. It was inferred from the genomic locations of genes comprising it that the moss CHS superfamily expanded through tandem and segmental duplication events. Inferred exon-intron architectures and results from phylogenetic analysis of representative CHS superfamily genes of P. patens and other plants showed that intron gain and loss occurred several times during evolution of this gene superfamily. A high proportion of P. patens CHS genes (7 of 14 genes for which the full sequence is known and probably 3 additional genes) are intronless, prompting speculation that CHS gene duplication via retrotransposition has occurred at least twice in the moss lineage. Analyses of sequence similarities, catalytic motifs and EST data indicated that a surprisingly large number (as many as 13) of the moss CHS superfamily genes probably encode active CHS. EST distribution data and different light responsiveness observed with selected genes provide evidence for their differential regulation. Observed diversity within the moss CHS superfamily and amenability to gene manipulation make Physcomitrella a highly suitable model system for studying expansion and functional diversification of the plant CHS superfamily of genes.

Multiple gains of spliceosomal introns in a superfamily of vertebrate protease inhibitor genes

Background

Intron gains reportedly are very rare during evolution of vertebrates, and the mechanisms underlying their creation are largely unknown. Previous investigations have shown that, during metazoan radiation, the exon-intron patterns of serpin superfamily genes were subject to massive changes, in contrast to many other genes.

Results

Here we investigated intron dynamics in the serpin superfamily in lineages pre- and postdating the split of vertebrates. Multiple intron gains were detected in a group of ray-finned fishes, once the canonical groups of vertebrate serpins had been established. In two genes, co-occurrence of non-standard introns was observed, implying that intron gains in vertebrates may even happen concomitantly or in a rapidly consecutive manner. DNA breakage/repair processes associated with genome compaction are introduced as a novel factor potentially favoring intron gain, since all non-canonical introns were found in a lineage of ray-finned fishes that experienced genomic downsizing.

Conclusion

Multiple intron acquisitions were identified in serpin genes of a lineage of ray-finned fishes, but not in any other vertebrates, suggesting that insertion rates for introns may be episodically increased. The co-occurrence of non-standard introns within the same gene discloses the possibility that introns may be gained simultaneously. The sequences flanking the intron insertion points correspond to the proto-splice site consensus sequence MAG↑N, previously proposed to serve as intron insertion site. The association of intron gains in the serpin superfamily with a group of fishes that underwent genome compaction may indicate that DNA breakage/repair processes might foster intron birth.

Endogenous mechanisms for the origins of spliceosomal introns

Over 30 years since their discovery, the origin of spliceosomal introns remains uncertain. One nearly universally accepted hypothesis maintains that spliceosomal introns originated from self-splicing group-II introns that invaded the uninterrupted genes of the last eukaryotic common ancestor (LECA) and proliferated by "insertion" events. Although this is a possible explanation for the original presence of introns and splicing machinery, the emphasis on a high number of insertion events in the genome of the LECA neglects a considerable body of empirical evidence showing that spliceosomal introns can simply arise from coding or, more generally, nonintronic sequences within genes. After presenting a concise overview of some of the most common hypotheses and mechanisms for intron origin, we propose two further hypotheses that are broadly based on central cellular processes: 1) internal gene duplication and 2) the response to aberrant and fortuitously spliced transcripts. These two nonmutually exclusive hypotheses provide a powerful way to explain the establishment of spliceosomal introns in eukaryotes without invoking an exogenous source.

Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Background

The 12-oxo-phytodienoic acid reductases (OPRs) are enzymes that catalyze the reduction of double-bonds in α, β-unsaturated aldehydes or ketones and are part of the octadecanoid pathway that converts linolenic acid to jasmonic acid. In plants, OPRs belong to the old yellow enzyme family and form multigene families. Although discoveries about this family in Arabidopsis and other species have been reported in some studies, the evolution and function of multiple OPRs in plants are not clearly understood.

Results

A comparative genomic analysis was performed to investigate the phylogenetic relationship, structural evolution and functional divergence among OPR paralogues in plants. In total, 74 OPR genes were identified from 11 species representing the 6 major green plant lineages: green algae, mosses, lycophytes, gymnosperms, monocots and dicots. Phylogenetic analysis showed that seven well-conserved subfamilies exist in plants. All OPR genes from green algae were clustered into a single subfamily, while those from land plants fell into six other subfamilies, suggesting that the events leading to the expansion of the OPR family occurred in land plants. Further analysis revealed that lineage-specific expansion, especially by tandem duplication, contributed to the current OPR subfamilies in land plants after divergence from aquatic plants. Interestingly, exon/intron structure analysis showed that the gene structures of OPR paralogues exhibits diversity in intron number and length, while the intron positions and phase were highly conserved across different lineage species. These observations together with the phylogenetic tree revealed that successive single intron loss, as well as indels within introns, occurred during the process of structural evolution of OPR paralogues. Functional divergence analysis revealed that altered functional constraints have occurred at specific amino acid positions after diversification of the paralogues. Most notably, significant functional divergence was also found in all pairs, except for the II/IV, II/V and V/VI pairs. Strikingly, analysis of the site-specific profiles established by posterior probability revealed that the positive-selection sites and/or critical amino acid residues for functional divergence are mainly distributed in α-helices and substrate binding loop (SBL), indicating the functional importance of these regions for this protein family.

Conclusion

This study highlights the molecular evolution of the OPR gene family in all plant lineages and indicates critical amino acid residues likely relevant for the distinct functional properties of the paralogues. Further experimental verification of these findings may provide valuable information on the OPRs' biochemical and physiological functions.

AtMKK1 mediates ABA-induced CAT1 expression and H2O2 production via AtMPK6-coupled signaling in Arabidopsis

Catalase controls cellular H(2)O(2) and plays important roles in the adaptation of plants to various stresses, but little is known about the signaling events that lead to the expression of CAT1 and the production of H(2)O(2). Here we report the dependence of CAT1 expression and H(2)O(2) production on a mitogen-activated protein kinase (MAPK) cascade. CAT1 transcript was induced in an ABA-dependent way and the induction was abolished in the T-DNA insertion mutant mkk1 (SALK_015914), while AtMKK1 overexpression significantly enhanced the ABA-induced CAT1 expression and H(2)O(2) production. AtMPK6, another component in the MAPK cascade, was also involved: mpk6 mutant blocked and overexpressing AtMPK6 enhanced the ABA-dependent expression of CAT1 and H(2)O(2) production. The activity of AtMPK6 was increased by ABA in an AtMKK1-dependent manner. These data clearly suggest an ABA-dependent signaling pathway connecting CAT1 expression through a phosphorelay including AtMKK1 and AtMPK6. In further support of this view, mkk1 mutant reduced both the sensitivity to ABA during germination and the drought tolerance of seedlings, whereas the AtMKK1 overexpression line showed the opposite responses when compared with the wild type. The data suggest AtMKK1-AtMPK6 to be a key module in an ABA-dependent signaling cascade causing H(2)O(2) production and stress responses.

Origins of human malaria: rare genomic changes and full mitochondrial genomes confirm the relationship of Plasmodium falciparum to other mammalian parasites but complicate the origins of Plasmodium vivax

Despite substantial work, the phylogeny of malaria parasites remains debated. The matter is complicated by concerns about patterns of evolution in potentially strongly selected genes as well as the extreme AT bias of some Plasmodium genomes. Particularly contentious has been the position of the most virulent human parasite Plasmodium falciparum, whether grouped with avian parasites or within a larger clade of mammalian parasites. Here, we study 3 classes of rare genomic changes, as well as the sequences of mitochondrial ribosomal RNA (rRNA) genes. We report 3 lines of support for a clade of mammalian parasites: 1) we find no instances of spliceosomal intron loss in a hypothetical ancestor of P. falciparum and the avian parasite Plasmodium gallinaceum, suggesting against a close relationship between those species; 2) we find 4 genomic mitochondrial indels supporting a mammalian clade, but none grouping P. falciparum with avian parasites; and 3) slowly evolving mitochondrial rRNA sequences support a mammalian parasite clade with 100% posterior probability. We further report a large deletion in the mitochondrial large subunit rRNA gene, which suggests a subclade including both African and Asian parasites within the clade of closely related primate malarias. This contrasts with previous studies that provided strong support for separate Asian and African clades, and reduces certainty about the historical and geographic origins of Plasmodium vivax. Finally, we find a lack of synapomorphic gene losses, suggesting a low rate of ancestral gene loss in Plasmodium.

Identification and characterization of the Arabidopsis gene encoding the tetrapyrrole biosynthesis enzyme uroporphyrinogen III synthase

UROS (uroporphyrinogen III synthase; EC 4.2.1.75) is the enzyme responsible for the formation of uroporphyrinogen III, the precursor of all cellular tetrapyrroles including haem, chlorophyll and bilins. Although UROS genes have been cloned from many organisms, the level of sequence conservation between them is low, making sequence similarity searches difficult. As an alternative approach to identify the UROS gene from plants, we used functional complementation, since this does not require conservation of primary sequence. A mutant of Saccharomyces cerevisiae was constructed in which the HEM4 gene encoding UROS was deleted. This mutant was transformed with an Arabidopsis thaliana cDNA library in a yeast expression vector and two colonies were obtained that could grow in the absence of haem. The rescuing plasmids encoded an ORF (open reading frame) of 321 amino acids which, when subcloned into an Escherichia coli expression vector, was able to complement an E. coli hemD mutant defective in UROS. Final proof that the ORF encoded UROS came from the fact that the recombinant protein expressed with an N-terminal histidine-tag was found to have UROS activity. Comparison of the sequence of AtUROS (A. thaliana UROS) with the human enzyme found that the seven invariant residues previously identified were conserved, including three shown to be important for enzyme activity. Furthermore, a structure-based homology search of the protein database with AtUROS identified the human crystal structure. AtUROS has an N-terminal extension compared with orthologues from other organisms, suggesting that this might act as a targeting sequence. The precursor protein of 34 kDa translated in vitro was imported into isolated chloroplasts and processed to the mature size of 29 kDa. Confocal microscopy of plant cells transiently expressing a fusion protein of AtUROS with GFP (green fluorescent protein) confirmed that AtUROS was targeted exclusively to chloroplasts in vivo.

The necrotroph Botrytis cinerea induces a non-host type II resistance mechanism in Pinus pinaster suspension-cultured cells

Models of non-host resistance have failed to account for the pathogenicity of necrotrophic agents. During the interaction of Pinus pinaster (maritime pine) with the non-host necrotrophic pathogen Botrytis cinerea, the generation and scavenging of reactive oxygen species (ROS) and the induction of the hypersensitive response (HR) were analyzed. Elicitation of maritime pine suspended cells with B. cinerea spores resulted in the biphasic induction of ROS. The phase I oxidative burst was dependent on calcium influx, while the phase II oxidative burst also depended on NADPH oxidase, protein kinase activity, and de novo transcription and protein synthesis. A decline was observed in catalase (CAT) and superoxide dismutase (SOD) activity, together with the down-regulation of Fe-Sod1, chlCu, Zn-Sod1 and csApx1, suggesting a coordinated response towards a decrease in the ROS-scavenging capacity of maritime pine cells during challenge. Following the second oxidative burst, programmed cell death events characteristic of the HR were observed. The results suggest the ROS-mediated and cell-breach-independent activation of Type II non-host resistance during the P. pinaster-B. cinerea interaction.