Functional Divergence Analysis Research Articles

Genome-Wide Identification, Evolution, and Expression Analysis of the MAPK Gene Family in Rosaceae Plants

Mitogen-activated protein kinases (MAPKs) are crucial regulators in coping with abiotic and biotic stresses, including drought, salinity, fungi, and pathogens. However, little is known about the characteristics, evolution process, and functional divergence of the MAPK gene family in Rosaceae plants. A total of 97 MAPK members were identified in six Rosaceae species, including 12 genes in Fragaria vesca, 22 genes in Malus domestica, 23 genes in Pyrus bretschneideri, 12 genes in Prunus mume, 14 genes in Prunus persica, and 14 genes in Rosa chinensis. All MAPK members of six Rosaceae plants were categorized into four clusters by the phylogenetic relationship analysis. Collinearity analysis discovered that both segmental duplication and tandem duplication contributed to the expansion of MAPK family genes in Rosaceae plants. And the analysis of motifs and gene structures indicated that the evolution of the MAPK gene family was highly conserved among phylogenetic clusters in Rosaceae species. In addition, the dN/dS rates of MAPK paralogous gene pairs were below one, suggesting the MAPK gene family in Rosaceae was driven by purifying selective pressure. Furthermore, functional divergence analysis discovered that 14 amino acid residues were detected as potentially key sites for functional divergence of MAPK family genes between different cluster pairs, specifically Type I functional divergence. The analysis of functional distance indicated that cluster C retained more of the original functional features, while cluster B exhibited functional specialization. Moreover, the expression profiles revealed that PmMAPK8, PmMAPK9, and PmMAPK10 were both highly expressed under drought stress and low temperature conditions. This study aims to comprehensively analyze the evolutionary process and functional analyses of the MAPK gene family in Rosaceae plants, which will lay the foundation for future studies into MAPK genes of Rosaceae in response to drought and cold stress.

Evolution and functional role prediction of the CYP6DE and CYP6DJ subfamilies in Dendroctonus (Curculionidae: Scolytinae) bark beetles.

Dendroctonus-bark beetles are natural components and key ecological agents of coniferous forests. They spend most of their lives under the bark, where they are exposed to highly toxic terpenes present in the oleoresin. Cytochrome P450 (CYP) is a multigene family involved in the detoxification of these compounds. It has been demonstrated that CYP6DE and CYP6DJ subfamilies hydroxylate monoterpenes, whose derivatives can act as pheromone synergist compounds or be pheromones themselves in these insects. Given the diversity and functional role of CYPs, we investigated whether these cytochromes have retained their function throughout the evolution of these insects. To test this hypothesis, we performed a Bayesian phylogenetic analysis to determine phylogenetic subgroups of cytochromes in these subfamilies. Subgroups were mapped and reconciled with the Dendroctonus phylogeny. Molecular docking analyses were performed with the cytochromes of each subgroup and enantiomers of α-pinene and β-pinene, (+)-3-carene, β-myrcene and R-(+)-limonene. In addition, functional divergence analysis was performed to identify critical amino acid sites that influence changes in catalytic site conformation and/or protein folding. Three and two phylogenetic subgroups were recovered for the CYP6DE and CYP6DJ subfamilies, respectively. Mapping and reconciliation analysis showed different gain and loss patterns for cytochromes of each subgroup. Functional predictions indicated that the cytochromes analyzed are able to hydroxylate all monoterpenes; however, they showed preferential affinities to different monoterpenes. Functional divergence analyses indicated that the CYP6DE subfamily has experimented type I and II divergence, whereas the CYP6DJ subfamily has evolved under strong functional constraints. Results suggest cytochromes of the CYP6DE subfamily evolve to reinforce their detoxifying capacity hydroxylating mainly α- and β-pinene to (+) and (-)-trans-verbenol, being the negative enantiomer used as a pheromone by several Dendroctonus species; whereas cytochromes of the CYP6DJ subfamily appear to retain their original function related to the detoxification of these compounds.

Functional mutation, splice, distribution, and divergence analysis of impactful genes associated with heart failure and other cardiovascular diseases

Cardiovascular disease (CVD) is caused by a multitude of complex and largely heritable conditions. Identifying key genes and understanding their susceptibility to CVD in the human genome can assist in early diagnosis and personalized treatment of the relevant patients. Heart failure (HF) is among those CVD phenotypes that has a high rate of mortality. In this study, we investigated genes primarily associated with HF and other CVDs. Achieving the goals of this study, we built a cohort of thirty-five consented patients, and sequenced their serum-based samples. We have generated and processed whole genome sequence (WGS) data, and performed functional mutation, splice, variant distribution, and divergence analysis to understand the relationships between each mutation type and its impact. Our variant and prevalence analysis found FLNA, CST3, LGALS3, and HBA1 linked to many enrichment pathways. Functional mutation analysis uncovered ACE, MME, LGALS3, NR3C2, PIK3C2A, CALD1, TEK, and TRPV1 to be notable and potentially significant genes. We discovered intron, 5ʹ Flank, 3ʹ UTR, and 3ʹ Flank mutations to be the most common among HF and other CVD genes. Missense mutations were less common among HF and other CVD genes but had more of a functional impact. We reported HBA1, FADD, NPPC, ADRB2, ADBR1, MYH6, and PLN to be consequential based on our divergence analysis.

Genetic basis underlying Lassa fever endemics in the Mano River region, West Africa

Lassa fever (LF), a haemorrhagic fever disease caused by Lassa virus (LASV), is a serious public health burden in West Africa. The Mano River region (Sierra Leone, Guinea, Liberia, and Côte d’Ivoire) has been an endemic focus of the disease over the past decades. Here, we deciphered the genetic basis underlying LF endemics in this region. Clade model and type I functional divergence analyses revealed that the major LASV group, Kenema sub-clade, which is currently circulating in the Eastern Province of Sierra Leone, has been affected by different selective pressure compared to isolates from the other areas with effects on the viral RNA-dependent RNA polymerase (L protein) and probably nucleoprotein (NP). Further, contingency analysis showed that, in the early endemic, the sub-clade has undergone adaptive diversification via acceleration of amino acid substitutions in L protein. These findings highlight the key viral factor and local adaptation regarding the endemicity of LF.

Functional Characterization of the Cystine-Rich-Receptor-like Kinases (CRKs) and Their Expression Response to Sclerotinia sclerotiorum and Abiotic Stresses in Brassica napus.

Cysteine-rich receptor-like kinases (CRKs) are transmembrane proteins that bind to the calcium ion to regulate stress-signaling and plant development-related pathways, as indicated by several pieces of evidence. However, the CRK gene family hasn't been inadequately examined in Brassica napus. In our study, 27 members of the CRK gene family were identified in Brassica napus, which are categorized into three phylogenetic groups and display synteny relationship to the Arabidopsis thaliana orthologs. All the CRK genes contain highly conserved N-terminal PKINASE domain; however, the distribution of motifs and gene structure were variable conserved. The functional divergence analysis between BnaCRK groups indicates a shift in evolutionary rate after duplication events, demonstrating that BnaCRKs might direct a specific function. RNA-Seq datasets and quantitative real-time PCR (qRT-PCR) exhibit the complex expression profile of the BnaCRKs in plant tissues under multiple stresses. Nevertheless, BnaA06CRK6-1 and BnaA08CRK8 from group B were perceived to play a predominant role in the Brassica napus stress signaling pathway in response to drought, salinity, and Sclerotinia sclerotiorum infection. Insights gained from this study improve our knowledge about the Brassica napus CRK gene family and provide a basis for enhancing the quality of rapeseed.

Genomic organization and gene evolution of two warm temperature acclimation proteins (Wap65s) of Micropterus salmoides and their responses to temperature and bacterial/viral infections

Warm temperature acclimation-related 65-kDa proteins (Wap65s) are fish plasma acute-phase glycoproteins homologous to hemopexin with high affinity and clearance for heme. The study characterized Mswap65-1 and Mswap65-2 genes in Micropterus salmoides. Structural analysis showed MsWap65s contained conserved heme-binding sites. MsWap65-1 had a chloride-binding site similar to hemopexin, while MsWap65-2 had an additional calcium-binding site. Phylogenetic and Ka/Ks analysis showed that fish Wap65s were evolutionarily conserved and underwent strong purifying selection. Functional divergence analysis indicated that fish Wap65-2 retained the putative function of ancestral Wap65, while Wap65-1 underwent neofunctional differentiation. QPCR showed Mswap65s were predominantly expressed in liver, but prolonged hyperthermy inhibited Mswap65-2 expression. Mswap65-2 expression was up-regulated in liver and spleen after Nocardia seriolae infection, while Mswap65-1 was down-regulated. MsWap65-2 may be associated with pathogenesis and play potential role in pathogen resistance. LMBV infection resulted in both significant downregulation of Mswap65s were both significantly down-regulated, with differences observed between sexes. We speculated the immune system might suppress expression after viral infection. Exogenous rMsWap65s were prepared, and injection of rMsWap65s alleviated phenylhydrazine-induced hemolysis and inhibited increases in heme, complement C3 and inflammatory symptoms. Our results contribute to an advanced understanding of the functions and mechanisms of MsWap65s in stress resistance.

Genome-Wide Characterization, Evolutionary Analysis of ARF Gene Family, and the Role of SaARF4 in Cd Accumulation of Sedum alfredii Hance.

Auxin response factors (ARFs) play important roles in plant development and environmental adaption. However, the function of ARFs in cadmium (Cd) accumulation are still unknown. Here, 23 SaARFs were detected in the genome of hyperaccumulating ecotype of Sedum alfredii Hance (HE), and they were not evenly distributed on the chromosomes. Their protein domains remained highly conservative. SaARFs in the phylogenetic tree can be divided into three groups. Genes in the group Ⅰ contained three introns at most. However, over ten introns were found in other two groups. Collinearity relationships were exhibited among ten SaARFs. The reasons for generating SaARFs may be segmental duplication and rearrangements. Collinearity analysis among different species revealed that more collinear genes of SaARFs can be found in the species with close relationships of HE. A total of eight elements in SaARFs promoters were related with abiotic stress. The qRT-PCR results indicated that four SaARFs can respond to Cd stress. Moreover, that there may be functional redundancy among six SaARFs. The adaptive selection and functional divergence analysis indicated that SaARF4 may undergo positive selection pressure and an adaptive-evolution process. Overexpressing SaARF4 effectively declined Cd accumulation. Eleven single nucleotide polymorphism (SNP) sites relevant to Cd accumulation can be detected in SaARF4. Among them, only one SNP site can alter the sequence of the SaARF4 protein, but the SaARF4 mutant of this site did not cause a significant difference in cadmium content, compared with wild-type plants. SaARFs may be involved in Cd-stress responses, and SaARF4 may be applied for decreasing Cd accumulation of plants.

Functional conservation and divergence in plant-specific GRF gene family revealed by sequences and expression analysis.

Unique to plants, growth regulatory factors (GRFs) play important roles in plant growth and reproduction. This study investigated the evolutionary and functional characteristics associated with plant growth. Using genome-wide analysis of 15 plant species, 173 members of the GRF family were identified and phylogenetically categorized into six groups. All members contained WRC and QLQ conserved domains, and the family’s expansion largely depended on segmental duplication. The promoter region of the GRF gene family mainly contained four types of cis-acting elements (light-responsive elements, development-related elements, hormone-responsive elements, and environmental stress-related elements) that are mainly related to gene expression levels. Functional divergence analysis revealed that changes in amino acid site evolution rate played a major role in the differentiation of the GRF gene family, with ten significant sites identified. Six significant sites were identified for positive selection. Moreover, the four groups of coevolutionary sites identified may play a key role in regulating the transcriptional activation of the GRF protein. Expression profiles revealed that GRF genes were generally highly expressed in young plant tissues and had tissue or organ expression specificity, demonstrating their functional conservation with distinct divergence. The results of these sequence and expression analyses are expected to provide molecular evolutionary and functional references for the plant GRF gene family.

Genome-Wide Prediction, Functional Divergence, and Characterization of Stress-Responsive BZR Transcription Factors in B. napus.

BRASSINAZOLE RESISTANT (BZR) are transcriptional factors that bind to the DNA of targeted genes to regulate several plant growth and physiological processes in response to abiotic and biotic stresses. However, information on such genes in Brassica napus is minimal. Furthermore, the new reference Brassica napus genome offers an excellent opportunity to systematically characterize this gene family in B. napus. In our study, 21 BnaBZR genes were distributed across 19 chromosomes of B. napus and clustered into four subgroups based on Arabidopsis thaliana orthologs. Functional divergence analysis among these groups evident the shifting of evolutionary rate after the duplication events. In terms of structural analysis, the BnaBZR genes within each subgroup are highly conserved but are distinctive within groups. Organ-specific expression analyses of BnaBZR genes using RNA-seq data and quantitative real-time polymerase chain reaction (qRT-PCR) revealed complex expression patterns in plant tissues during stress conditions. In which genes belonging to subgroups III and IV were identified to play central roles in plant tolerance to salt, drought, and Sclerotinia sclerotiorum stress. The insights from this study enrich our understanding of the B. napus BZR gene family and lay a foundation for future research in improving rape seed environmental adaptability.

Molecular epidemiology of yellow fever virus in Africa: A perspective of the phylogeographic split between East/Central African and West African lineages

Yellow fever (YF) is a major public-health problem in Africa. Yellow fever virus (YFV), the etiological agent responsible for the disease, exhibits clear delineation of phylogeography between East/Central Africa and West Africa. In order to decipher the genetic nature of the YFV epidemic between these areas, we performed a genome-wide study on its African isolates using the McDonald-Kreitman (MK) test in combination with the type II functional divergence analysis. The results showed that adaptive genetic diversifications have occurred on viral nonstructural protein 1 (NS1) and NS5, which are essential for viral genome replication and immune antagonism, with the East/Central African–West African epidemic split. On both proteins, a number of amino acid replacements have been favored by functional divergence. These findings could help to bridge the gap between the phylogeographic delineation and niche adaptation underlying the YFV-epidemic across Africa and shed light on viral determinants of this process.

Genome-wide identification and evolutionary analysis of MLO gene family in Rosaceae plants

Mutants lacking wild-type MLO (Mildew resistance Locus O) proteins show broad-spectrum resistance to the powdery mildew fungus, and dysregulated cell death control, with spontaneous cell death in response to developmental or abiotic stimuli. In order to understand the evolution and divergence patterns of the MLO gene family in Rosaceae plants, we analysed systematically genome-wide data from Fragaria vesca, Prunus persica, Prunus mume, Malus domestica, Pyrus bretschneideri and Rubus occidentalis based on bioinformatics methods. Using three phylogenetic methods (the neighbour-joining, maximum likelihood, and Bayesian methods), we identified 117 MLO genes from 6 Rosaceae species. The results of all three phylogenetic analysis methods supported that these genes were divided into six clades. Conserved motif analysis found that only motif 2 was present in all MLO proteins and had 3 nearly invariant amino acid residues. The findings indicated that motif 2 might be shared by the MLO gene family. The structural features of these genes showed large variations in sequence length among different species, although the lengths and the numbers of exons exhibited high degrees of similarity. Selective pressure analysis showed extremely significant differences in all 6 clades, with 2, 1, and 1 site(s) under significant positive selection detected in clades III, IV, and VI, respectively. These positive selection sites were important driving forces for the promotion of the functional differentiation of the MLO genes. Functional divergence analysis showed that the significantly divergent sites were located within the domains of the MLO genes. Functional distance analysis showed that the clade V had more conservative functions and might have retained more original functions during the evolutionary process. However, clade I may have undergone extensive altered functional constraints as a specialised functional role. Moreover, the most original function of the MLO genes in Rosaceae could be related to the evolution of their resistance to powdery mildew, which then gradually evolved into functions such as the regulation of flower development, the control of root morphology, and seed evolution due to the different evolutionary rates after gene duplication. These results provide a theoretical basis for further studies of the molecular evolutionary patterns of the plant MLO gene family.

Evolution and functional divergence of the ERBB receptor family.

Aim: The ERBB gene family plays an important rolein cell proliferation and differentiation, and aberrant activations could result in tumorigenesis, which makes this gene family an attractive drug target in the area of precision oncology. Materials & methods: Functional divergence analysis and conservation analysis were performed using ClustalW, MEGA7 and DIVERGE3 software. Results: One hundred and forty fivefunctional divergence residues sites, 94 totally conserved sites and averagely 345 conserved sites of individual gene member were obtained. Some have been reported to play rolein drug binding, tumorigenesis and drug resistance. Conclusion: Functional divergence residues with high posterior probabilities and conserved residues may possess certain functions, and aberrant alterations may confer drug resistance or contribute to tumorigenesis.

Functional divergence and adaptive selection of KNOX gene family in plants.

KNOTTED-like homeodomain (KNOX) genes are transcriptional regulators that play an important role in morphogenesis. In the present study, a comparative analysis was performed to investigate the molecular evolution of the characteristics of the KNOX gene family in 10 different plant species. We identified 129 KNOX gene family members, which were categorized into two subfamilies based on multiple sequence alignment and phylogenetic tree reconstruction. Several segmental duplication pairs were found, indicating that different species share a common expansion model. Functional divergence analysis identified the 15 and 52 amino acid sites with significant changes in evolutionary rates and amino acid physicochemical properties as functional divergence sites. Additional selection analysis showed that 14 amino acid sites underwent positive selection during evolution, and two groups of co-evolutionary amino acid sites were identified by Coevolution Analysis using Protein Sequences software. These sites could play critical roles in the molecular evolution of the KNOX gene family in these species. In addition, the expression profiles of KNOX duplicated genes demonstrated functional divergence. Taken together, these results provide novel insights into the structural and functional evolution of the KNOX gene family.

New insights into the evolution and expression dynamics of invertase gene family in Solanum lycopersicum

Invertases catalyze the irreversible hydrolysis of sucrose into glucose and fructose and thus play key roles in carbon metabolism and plant development. To gain insights into their evolutional and functional relationships, we conducted genome-wide analyses of invertase genes in tomato and other species, focusing on their evolution and expression dynamics. The analyses unexpectedly identified in the tomato genome 5 pseudo invertsase sequences and 5 non-functional cell wall invertases (CWINs) lacking the critical β-fructosidase motif or other amino acids required for hydrolyzing sucrose. Based on their phylogeny relationship and exon–intron structure, we speculated that the invertase gene family could arose from different ancestral genes. The acid invertase gene family, comprised of CWIN and vacuolar invertase (VIN), expanded through segmental and tandem duplication. Analysis of functional divergence suggests site-specific shifted evolutionary rate (Type-I) have played an important role in evolutionary novelties after acid invertase gene duplication in plants. Finally, paralogs within each of the CWIN, VIN and CIN subfamilies exhibited diverse expression responses to the same set of stress treatments including salt and temperature stresses, probably reflecting functional adaptability of the invertase genes during evolution.

Genome-wide analysis of wheat DNA-binding with one finger (Dof) transcription factor genes: evolutionary characteristics and diverse abiotic stress responses

BackgroundDNA binding with one finger (Dof) transcription factors play important roles in plant growth and abiotic stress responses. Although genome-wide identification and analysis of the DOF transcription factor family has been reported in other species, no relevant studies have emerged in wheat. The aim of this study was to investigate the evolutionary and functional characteristics associated with plant growth and abiotic stress responses by genome-wide analysis of the wheat Dof transcription factor gene family.ResultsUsing the recently released wheat genome database (IWGSC RefSeq v1.0), we identified 96 wheat Dof gene family members, which were phylogenetically clustered into five distinct subfamilies. Gene duplication analysis revealed a broad and heterogeneous distribution of TaDofs on the chromosome groups 1 to 7, and obvious tandem duplication genes were present on chromosomes 2 and 3.Members of the same gene subfamily had similar exon-intron structures, while members of different subfamilies had obvious differences. Functional divergence analysis indicated that type-II functional divergence played a major role in the differentiation of the TaDof gene family. Positive selection analysis revealed that the Dof gene family experienced different degrees of positive selection pressure during the process of evolution, and five significant positive selection sites (30A, 31 T, 33A, 102G and 104S) were identified. Additionally, nine groups of coevolving amino acid sites, which may play a key role in maintaining the structural and functional stability of Dof proteins, were identified. The results from the RNA-seq data and qRT-PCR analysis revealed that TaDof genes exhibited obvious expression preference or specificity in different organs and developmental stages, as well as in diverse abiotic stress responses. Most TaDof genes were significantly upregulated by heat, PEG and heavy metal stresses.ConclusionsThe genome-wide analysis and identification of wheat DOF transcription factor family and the discovery of important amino acid sites are expected to provide new insights into the structure, evolution and function of the plant Dof gene family.

Apolipoprotein D in Lepidoptera: Evolution and functional divergence

Apolipoprotein D (ApoD) plays important roles in response to injury, cell differentiation, lifespan extension, and increasing stress resistance. However, the evolutionary mechanism of ApoD in insects remains largely unelucidated. We conducted a comprehensive study of the molecular evolution and functional divergence of ApoD in insects. A type I functional divergence analysis revealed significant differences among insect ApoD homologs, suggesting that they underwent functional divergence. We demonstrated that lepidopteran insects have three genes that are close homologs to ApoD and show divergences in sequence, expression pattern, and protein-protein interaction. Furthermore, positive selection was detected in lepidopteran ApoD2, and positively selected sites were located around the pocket and loop domains, which might result in conformational changes and affect binding properties. Moreover, we showed that the three ApoDs in Bombyx mori were significantly regulated by environmental stress. Thus, this work illustrates the dialectical relationship between genetic diversity and functional conservation of ApoD and highlights its unique functions in the stress response of Lepidoptera.

Evolutionary divergence of function and expression of laccase genes in plants

Laccases (LACs) are versatile enzymes that catalyze oxidation of a wide range of substrates, thereby functioning in regulation of plant developmental processes and stress responses. However, with a few exceptions, the function of most LACs remains unclear in plants. In this study, we newly identified 4, 12, 22, 26, 27, 28 and 49 LAC genes for Physcomitrella patens, Amborella trichopoda, Zeamays, Ricinus communis, Vitis vinifera, Triticum aestivum and Glycine max, on the basis of exhaustive homologous sequence searches. In these plants, LACs differ greatly in sequence length and physical properties, such as molecular weight and theoretical isoelectric point (pI), but majority of them contain a signal peptide at their N-terminus. The originality of LACs could be traced back to as early as the emergence of moss. Plant LACs are clearly divided into seven distinct classes, where six ancient LACs should be present prior to the divergence of gymnosperms and angiosperms. Functional divergence analysis reveal that functional differentiation should occur among different groups of LACs because of altered selective constraints working on some critical amino acid sites (CAASs) within conserved laccase domains during evolution. Soybean and maize LACs have significantly different exon frequency (6.08 vs 4.82), and they are unevenly distributed and tend to form gene clusters on some chromosomes. Further analysis shows that the expansion of LAC gene family would be due toextensive tandem and chromosomal segmental duplications in the two plant species. Interestingly, *81.6% and 36.4% of soybean and maize LACs are potential targets of miRNAs, such as miR397a/b, miR408d, or miR528a/b etc. Both soybean and maize LACs are tissue specifically and developmental-specifically expressed, and are in response to different external abiotic and biotic stressors. These results suggest a diversity of functions of plant LAC genes, which will broaden our understanding and lay solid foundation for further investigating their biological functions in plants.

Comparative functional genomics analysis of cytochrome P450 gene superfamily in wheat and maize

BackgroundThe cytochrome P450s (CYP450s) as the largest enzyme family of plant metabolism participate in various physiological processes, whereas no study has demonstrated interest in comprehensive comparison of the genes in wheat and maize. Genome-wide survey, characterization and comparison of wheat and maize CYP450 gene superfamily are useful for genetic manipulation of the Gramineae crops.ResultsIn total, 1285 and 263 full-length CYP450s were identified in wheat and maize, respectively. According to standard nomenclature, wheat CYP450s (TaCYP450s) were categorized into 45 families, while maize CYP450s (ZmCYP450s) into 43 families. A comprehensive analysis of wheat and maize CYP450s, involved in functional domains, conserved motifs, phylogeny, gene structures, chromosome locations and duplicated events was performed. The result showed that each family/subfamily in both species exhibited characteristic features, suggesting their phylogenetic relationship and the potential divergence in their functions. Functional divergence analysis at the amino acid level of representative clans CYP51, CYP74 and CYP97 in wheat, maize and rice identified some critical amino acid sites that are responsible for functional divergence of a gene family. Expression profiles of Ta-, ZmCYP450s were investigated using RNA-seq data, which contribute to infer the potential functions of the genes during development and stress responses. We found in both species CYP450s had preferential expression in specific tissues, and many tissue-specific genes were identified. Under water-deficit condition, 82 and 39 significantly differentially expressed CYP450s were respectively detected in wheat and maize. These genes may have some roles in protecting plants against drought damage. Thereinto, fourteen CYP450s were selected to validate their expression level through qRT-PCR. To further elucidating molecular mechanisms of CYP450 action, gene co-expression network was constructed. In total, 477 TaCYP450s were distributed in 22 co-expression modules, and some co-expressed genes that likely take part in the same biochemical pathway were identified. For instance, the expression of TaCYP74A98_4D was highly correlated with TaLOX9, TaLOX36, TaLOX39, TaLOX44 and TaOPR8, and all of them may be involved in jasmonate (JA) biosynthesis. TaCYP73A201_3A showed coexpression with TaPAL1.25, TaCCoAOMT1.2, TaCOMT.1, TaCCR1.6 and TaLAC5, which probably act in the wheat stem and/or root lignin synthesis pathway.ConclusionOur study first established systematic information about evolutionary relationship, expression pattern and function characterization of CYP450s in wheat and maize.

Genome-wide phylogenetic and structural analysis reveals the molecular evolution of the ABA receptor gene family.

The plant hormone abscisic acid (ABA) plays a crucial role during the plant life cycle as well as in adaptive responses to environmental stresses. The core regulatory components of ABA signaling in plants are the pyrabactin resistance1/PYR1-like/regulatory component of ABA receptor family (PYLs), which comprise the largest plant hormone receptor family known. They act as negative regulators of members of the protein phosphatase type 2C family. Due to the biological importance of PYLs, many researchers have focused on their genetic redundancy and consequent functional divergence. However, little is understood of their evolution and its impact on the generation of regulatory diversity. In this study, we identify positive selection and functional divergence in PYLs through phylogenetic reconstruction, gene structure and expression pattern analysis, positive selection analysis, functional divergence analysis, and structure comparison. We found the correlation of desensitization of PYLs under specific modifications in the molecular recognition domain with functional diversification. Hence, an interesting antagonistic co-evolutionary mechanism is proposed for the functional diversification of ABA receptor family proteins. We believe a compensatory evolutionary pathway may have occurred.

Unraveling the evolutionary origin of ELR motif using fish CXC chemokine CXCL8

Chemokines are chemotactic proteins involved in host defense through the migration of immune-regulatory cells to the site of infection. Interleukin-8 (CXCL8/IL8) is the most studied “ELR-CXC chemokine/neutrophil activating chemokine (NAC) that regulate neutrophil trafficking during infections and inflammation by binding to its cognate G-protein coupled receptors CXCR1/CXCR2. The “ELR” motif of NAC chemokines is essential for the CXCR1/CXCR2 receptor activation. In order to understand the evolutionary origin of “ELR” motif in the CXC chemokines, a thorough evolutionary study of CXCL8 gene from various fishes and primates was performed. Phylogenetic analysis revealed that the CXCL8 gene can be classified into four distinct lineages (CXCL8-L1a, CXCL8-L1b, CXCL8-L2, and CXCL8-L3), where CXCL8-L1a is the fastest evolving lineage and CXCL8-L3 is the slowest. Selection analysis suggested that The “ELR/DLR” motif containing branches (gadoid and coelacanth) are positively selected. The probable evolutionary trend of “ELR” motif suggested that this motif in ancestor CXCL8 is evolved from the GGR of Lamprey (Agnatha), followed by duplication giving rise to two main motifs in CXCL8 “NXH” in L3 lineage and “ELR/DLR” in L1a/L1b lineages. Although, structural analysis suggested that the overall topology of the CXCL8 proteins is similar, differences do exist at the individual structural elements among the members of different lineages. Functional distance analysis suggested that the CXCL8-L3 lineage is more distant compared to the CXCL8-L1a and L1b lineages from the inferred ancestor. Functional divergence analysis between different lineages suggested that most of the selected residues are important for receptor or glycosaminoglycan binding. Such a functional diversification can be attributed to the novel set of functions adopted by CXCL8 in various species.