COLD1 confers chilling tolerance in rice.

Rice is sensitive to cold and can be grown only in certain climate zones. Human selection of japonica rice has extended its growth zone to regions with lower temperature, while the molecular basis of this adaptation remains unknown. Here, we identify the quantitative trait locus COLD1 that confers chilling tolerance in japonica rice. Overexpression of COLD1(jap) significantly enhances chilling tolerance, whereas rice lines with deficiency or downregulation of COLD1(jap) are sensitive to cold. COLD1 encodes a regulator of G-protein signaling that localizes on plasma membrane and endoplasmic reticulum (ER). It interacts with the G-protein α subunit to activate the Ca(2+) channel for sensing low temperature and to accelerate G-protein GTPase activity. We further identify that a SNP in COLD1, SNP2, originated from Chinese Oryza rufipogon, is responsible for the ability of COLD(jap/ind) to confer chilling tolerance, supporting the importance of COLD1 in plant adaptation.

GW9 determines grain size and floral organ identity in rice.

Grain weight is an important determinant of grain yield. However, the underlying regulatory mechanisms for grain size remain to be fully elucidated. Here, we identify a rice mutant grain weight 9 (gw9), which exhibits larger and heavier grains due to excessive cell proliferation and expansion in spikelet hull. GW9 encodes a nucleus-localized protein containing both C2H2 zinc finger (C2H2-ZnF) and VRN2-EMF2-FIS2-SUZ12 (VEFS) domains, serving as a negative regulator of grain size and weight. Interestingly, the non-frameshift mutations in C2H2-ZnF domain result in increased plant height and larger grain size, whereas frameshift mutations in both C2H2-ZnF and VEFS domains lead to dwarf and malformed spikelet. These observations indicated the dual functions of GW9 in regulating grain size and floral organ identity through the C2H2-ZnF and VEFS domains, respectively. Further investigation revealed the interaction between GW9 and the E3 ubiquitin ligase protein GW2, with GW9 being the target of ubiquitination by GW2. Genetic analyses suggest that GW9 and GW2 function in a coordinated pathway controlling grain size and weight. Our findings provide a novel insight into the functional role of GW9 in the regulation of grain size and weight, offering potential molecular strategies for improving rice yield.

Formin protein DRT1 affects gross morphology and chloroplast relocation in rice.

Plant height and tiller number are two major factors determining plant architecture and yield. However, in rice (Oryza sativa), the regulatory mechanism of plant architecture remains to be elucidated. Here, we reported a recessive rice mutant presenting dwarf and reduced tillering phenotypes (drt1). Map-based cloning revealed that the phenotypes are caused by a single point mutation in DRT1, which encodes the Class I formin protein O. sativa formin homolog 13 (OsFH13), binds with F-actin, and promotes actin polymerization for microfilament organization. DRT1 protein localized on the plasma membrane (PM) and chloroplast (CP) outer envelope. DRT1 interacted with rice phototropin 2 (OsPHOT2), and the interaction was interrupted in drt1. Upon blue light stimulus, PM localized DRT1 and OsPHOT2 were translocated onto the CP membrane. Moreover, deficiency of DRT1 reduced OsPHOT2 internalization and OsPHOT2-mediated CP relocation. Our study suggests that rice formin protein DRT1/OsFH13 is necessary for plant morphology and CP relocation by modulating the actin-associated cytoskeleton network.

The LARGE2-APO1/APO2 regulatory module controls panicle size and grain number in rice.

Panicle size and grain number are important agronomic traits and influence grain yield in rice (Oryza sativa), but the molecular and genetic mechanisms underlying panicle size and grain number control remain largely unknown in crops. Here we report that LARGE2 encodes a HECT-domain E3 ubiquitin ligase OsUPL2 and regulates panicle size and grain number in rice. The loss of function large2 mutants produce large panicles with increased grain number, wide grains and leaves, and thick culms. LARGE2 regulates panicle size and grain number by repressing meristematic activity. LARGE2 is highly expressed in young panicles and grains. Biochemical analyses show that LARGE2 physically associates with ABERRANT PANICLE ORGANIZATION1 (APO1) and APO2, two positive regulators of panicle size and grain number, and modulates their stabilities. Genetic analyses support that LARGE2 functions with APO1 and APO2 in a common pathway to regulate panicle size and grain number. These findings reveal a novel genetic and molecular mechanism of the LARGE2-APO1/APO2 module-mediated control of panicle size and grain number in rice, suggesting that this module is a promising target for improving panicle size and grain number in crops.

Carotenoid isomerase regulates rice tillering and grain productivity by its biosynthesis pathway.

Carotenoid isomerase activity and carotenoid content maintain the appropriate tiller number, photosynthesis, and grain yield. Interactions between the strigolactone and abscisic acid pathways regulates tiller formation.

NLG1, encoding a mitochondrial membrane protein, controls leaf and grain development in rice.

BACKGROUND: Mitochondrion is the key respiratory organ and participate in multiple anabolism and catabolism pathways in eukaryote. However, the underlying mechanism of how mitochondrial membrane proteins regulate leaf and grain development remains to be further elucidated. RESULTS: Here, a mitochondria-defective mutant narrow leaf and slender grain 1 (nlg1) was identified from an EMS-treated mutant population, which exhibits narrow leaves and slender grains. Moreover, nlg1 also presents abnormal mitochondria structure and was sensitive to the inhibitors of mitochondrial electron transport chain. Map-based cloning and transgenic functional confirmation revealed that NLG1 encodes a mitochondrial import inner membrane translocase containing a subunit Tim21. GUS staining assay and RT-qPCR suggested that NLG1 was mainly expressed in leaves and panicles. The expression level of respiratory function and auxin response related genes were significantly down-regulated in nlg1, which may be responsible for the declination of ATP production and auxin content. CONCLUSIONS: These results suggested that NLG1 plays an important role in the regulation of leaf and grain size development by maintaining mitochondrial homeostasis. Our finding provides a novel insight into the effects of mitochondria development on leaf and grain morphogenesis in rice.

SEMI-ROLLED LEAF 10 stabilizes catalase isozyme B to regulate leaf morphology and thermotolerance in rice (Oryza sativa L.).

Plant architecture and stress tolerance play important roles in rice breeding. Specific leaf morphologies and ideal plant architecture can effectively improve both abiotic stress resistance and rice grain yield. However, the mechanism by which plants simultaneously regulate leaf morphogenesis and stress resistance remains elusive. Here, we report that SRL10, which encodes a double-stranded RNA-binding protein, regulates leaf morphology and thermotolerance in rice through alteration of microRNA biogenesis. The srl10 mutant had a semi-rolled leaf phenotype and elevated sensitivity to high temperature. SRL10 directly interacted with catalase isozyme B (CATB), and the two proteins mutually increased one other's stability to enhance hydrogen peroxide (H2 O2 ) scavenging, thereby contributing to thermotolerance. The natural Hap3 (AGC) type of SRL10 allele was found to be present in the majority of aus rice accessions, and was identified as a thermotolerant allele under high temperature stress in both the field and the growth chamber. Moreover, the seed-setting rate was 3.19 times higher and grain yield per plant was 1.68 times higher in near-isogenic line (NIL) carrying Hap3 allele compared to plants carrying Hap1 allele under heat stress. Collectively, these results reveal a new locus of interest and define a novel SRL10-CATB based regulatory mechanism for developing cultivars with high temperature tolerance and stable yield. Furthermore, our findings provide a theoretical basis for simultaneous breeding for plant architecture and stress resistance.

Control of Grain Size and Weight by the GSK2-LARGE1/OML4 Pathway in Rice.

Regulation of grain size is crucial for improving crop yield and is also a basic aspect in developmental biology. However, the genetic and molecular mechanisms underlying grain size control in crops remain largely unknown despite their central importance. Here, we report that the MEI2-LIKE PROTEIN4 (OML4) encoded by the LARGE1 gene is phosphorylated by GLYCOGEN SYNTHASE KINASE2 (GSK2) and negatively controls grain size and weight in rice (Oryza sativa). Loss of function of OML4 leads to large and heavy grains, while overexpression of OML4 causes small and light grains. OML4 regulates grain size by restricting cell expansion in the spikelet hull. OML4 is expressed in developing panicles and grains, and the GFP-OML4 fusion protein is localized in the nuclei. Biochemical analyses show that the GSK2 physically interacts with OML4 and phosphorylates it, thereby possibly influencing the stability of OML4. Genetic analyses support that GSK2 and OML4 act, at least in part, in a common pathway to control grain size in rice. These results reveal the genetic and molecular mechanism of a GSK2-OML4 regulatory module in grain size control, suggesting that this pathway is a suitable target for improving seed size and weight in crops.

Karrikin Signaling Acts Parallel to and Additively with Strigolactone Signaling to Regulate Rice Mesocotyl Elongation in Darkness.

Seedling emergence in monocots depends mainly on mesocotyl elongation, requiring coordination between developmental signals and environmental stimuli. Strigolactones (SLs) and karrikins are butenolide compounds that regulate various developmental processes; both are able to negatively regulate rice (Oryza sativa) mesocotyl elongation in the dark. Here, we report that a karrikin signaling complex, DWARF14-LIKE (D14L)-DWARF3 (D3)-O. sativa SUPPRESSOR OF MAX2 1 (OsSMAX1) mediates the regulation of rice mesocotyl elongation in the dark. We demonstrate that D14L recognizes the karrikin signal and recruits the SCFD3 ubiquitin ligase for the ubiquitination and degradation of OsSMAX1, mirroring the SL-induced and D14- and D3-dependent ubiquitination and degradation of D53. Overexpression of OsSMAX1 promoted mesocotyl elongation in the dark, whereas knockout of OsSMAX1 suppressed the elongated-mesocotyl phenotypes of d14l and d3 OsSMAX1 localizes to the nucleus and interacts with TOPLESS-RELATED PROTEINs, regulating downstream gene expression. Moreover, we showed that the GR24 enantiomers GR245DS and GR24 ent-5DS specifically inhibit mesocotyl elongation and regulate downstream gene expression in a D14- and D14L-dependent manner, respectively. Our work revealed that karrikin and SL signaling play parallel and additive roles in modulating downstream gene expression and negatively regulating mesocotyl elongation in the dark.

WHITE AND LESION-MIMIC LEAF1, encoding a lumazine synthase, affects reactive oxygen species balance and chloroplast development in rice.

The riboflavin derivatives flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) are essential cofactors for enzymes in multiple cellular processes. Characterizing mutants with impaired riboflavin metabolism can help clarify the role of riboflavin in plant development. Here, we characterized a rice (Oryza sativa) white and lesion-mimic (wll1) mutant, which displays a lesion-mimic phenotype with white leaves, chlorophyll loss, chloroplast defects, excess reactive oxygen species (ROS) accumulation, decreased photosystem protein levels, changes in expression of chloroplast development and photosynthesis genes, and cell death. Map-based cloning and complementation test revealed that WLL1 encodes lumazine synthase, which participates in riboflavin biosynthesis. Indeed, the wll1 mutant showed riboflavin deficiency, and application of FAD rescued the wll1 phenotype. In addition, transcriptome analysis showed that cytokinin metabolism was significantly affected in wll1 mutant, which had increased cytokinin and δ-aminolevulinic acid contents. Furthermore, WLL1 and riboflavin synthase (RS) formed a complex, and the rs mutant had a similar phenotype to the wll1 mutant. Taken together, our findings revealed that WLL1 and RS play pivotal roles in riboflavin biosynthesis, which is necessary for ROS balance and chloroplast development in rice.

Disruption of EARLY LESION LEAF 1, encoding a cytochrome P450 monooxygenase, induces ROS accumulation and cell death in rice.

Lesion-mimic mutants (LMMs) provide a valuable tool to reveal the molecular mechanisms determining programmed cell death (PCD) in plants. Despite intensive research, the mechanisms behind PCD and the formation of lesions in various LMMs still remain to be elucidated. Here, we identified a rice (Oryza sativa) LMM, early lesion leaf 1 (ell1), cloned the causal gene by map-based cloning, and verified this by complementation. ELL1 encodes a cytochrome P450 monooxygenase, and the ELL1 protein was located in the endoplasmic reticulum. The ell1 mutant exhibited decreased chlorophyll contents, serious chloroplast degradation, upregulated expression of chloroplast degradation-related genes, and attenuated photosynthetic protein activity, indicating that ELL1 is involved in chloroplast development. RNA sequencing analysis showed that genes related to oxygen binding were differentially expressed in ell1 and wild-type plants; histochemistry and paraffin sectioning results indicated that hydrogen peroxide (H2 O2 ) and callose accumulated in the ell1 leaves, and the cell structure around the lesions was severely damaged, which indicated that reactive oxygen species (ROS) accumulated and cell death occurred in the mutant. TUNEL staining and comet experiments revealed that severe DNA degradation and abnormal PCD occurred in the ell1 mutants, which implied that excessive ROS accumulation may induce DNA damage and ROS-mediated cell death in the mutant. Additionally, lesion initiation in the ell1 mutant was light dependent and temperature sensitive. Our findings revealed that ELL1 affects chloroplast development or function, and that loss of ELL1 function induces ROS accumulation and lesion formation in rice.

UDP-N-acetylglucosamine pyrophosphorylase enhances rice survival at high temperature.

High-temperature stress inhibits normal cellular processes and results in abnormal growth and development in plants. However, the mechanisms by which rice (Oryza sativa) copes with high temperature are not yet fully understood. In this study, we identified a rice high temperature enhanced lesion spots 1 (hes1) mutant, which displayed larger and more dense necrotic spots under high temperature conditions. HES1 encoded a UDP-N-acetylglucosamine pyrophosphorylase, which had UGPase enzymatic activity. RNA sequencing analysis showed that photosystem-related genes were differentially expressed in the hes1 mutant at different temperatures, indicating that HES1 plays essential roles in maintaining chloroplast function. HES1 expression was induced under high temperature conditions. Furthermore, loss-of-function of HES1 affected heat shock factor expression and its mutation exhibited greater vulnerability to high temperature. Several experiments revealed that higher accumulation of reactive oxygen species occurred in the hes1 mutant at high temperature. Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) and comet experiments indicated that the hes1 underwent more severe DNA damage at high temperature. The determination of chlorophyll content and chloroplast ultrastructure showed that more severe photosystem defects occurred in the hes1 mutant under high temperature conditions. This study reveals that HES1 plays a key role in adaptation to high-temperature stress in rice.

The kinesin-13 protein BR HYPERSENSITIVE 1 is a negative brassinosteroid signaling component regulating rice growth and development.

Phytohormones performed critical roles in regulating plant architecture and thus determine grain yield in rice. However, the roles of brassinosteroids (BRs) compared to other phytohormones in shaping rice architecture are less studied. In this study, we report that BR hypersensitive1 (BHS1) plays a negative role in BR signaling and regulate rice architecture. BHS1 encodes the kinesin-13a protein and regulates grain length. We found that bhs1 was hypersensitive to BR, while BHS1-overexpression was less sensitive to BR compare to WT. BHS1 was down-regulated at RNA and protein level upon exogenous BR treatment, and proteasome inhibitor MG132 delayed the BHS1 degradation, indicating that both the transcriptional and posttranscriptional regulation machineries are involved in BHS1-mediated regulation of plant growth and development. Furthermore, we found that the BR-induced degradation of BHS1 was attenuated in Osbri1 and Osbak1 mutants, but not in Osbzr1 and Oslic mutants. Together, these results suggest that BHS1 is a novel component which is involved in negative regulation of the BR signaling downstream player of BRI1.

Integrated Multi-Omics Perspective to Strengthen the Understanding of Salt Tolerance in Rice.

Salt stress is one of the major constraints to rice cultivation worldwide. Thus, the development of salt-tolerant rice cultivars becomes a hotspot of current rice breeding. Achieving this goal depends in part on understanding how rice responds to salt stress and uncovering the molecular mechanism underlying this trait. Over the past decade, great efforts have been made to understand the mechanism of salt tolerance in rice through genomics, transcriptomics, proteomics, metabolomics, and epigenetics. However, there are few reviews on this aspect. Therefore, we review the research progress of omics related to salt tolerance in rice and discuss how these advances will promote the innovations of salt-tolerant rice breeding. In the future, we expect that the integration of multi-omics salt tolerance data can accelerate the solution of the response mechanism of rice to salt stress, and lay a molecular foundation for precise breeding of salt tolerance.

Genome-Wide Association Study Reveals Novel QTLs and Candidate Genes for Grain Number in Rice.

Grain number per panicle (GNPP), determined mainly by panicle branching, is vital for rice yield. The dissection of the genetic basis underlying GNPP could help to improve rice yield. However, genetic resources, including quantitative trait loci (QTL) or genes for breeders to enhance rice GNPP, are still limited. Here, we conducted the genome-wide association study (GWAS) on the GNPP, primary branch number (PBN), and secondary branch number (SBN) of 468 rice accessions. We detected a total of 18 QTLs, including six for GNPP, six for PBN, and six for SBN, in the whole panel and the indica and japonica subpanels of 468 accessions. More importantly, qPSG1 was a common QTL for GNPP, PBN, and SBN and was demonstrated by chromosome segment substitution lines (CSSLs). Considering gene annotation, expression, and haplotype analysis, seven novel and strong GNPP-related candidate genes were mined from qPSG1. Our results provide clues to elucidate the molecular regulatory network of GNPP. The identified QTLs and candidate genes will contribute to the improvement of GNPP and rice yield via molecular marker-assisted selection (MAS) breeding and genetic engineering techniques.

Primary leaf-type ferredoxin 1 participates in photosynthetic electron transport and carbon assimilation in rice.

Ferredoxins (Fds) play a crucial role in photosynthesis by regulating the distribution of electrons to downstream enzymes. Multiple Fd genes have been annotated in the Oryza sativa L. (rice) genome; however, their specific functions are not well understood. Here, we report the functional characterization of rice Fd1. Sequence alignment, phylogenetic analysis of seven rice Fd proteins and quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that rice Fd1 is a primary leaf-type Fd. Electron transfer assays involving NADP+ and cytochrome c indicated that Fd1 can donate electrons from photosystem I (PSI) to ferredoxin-NADP+ reductase. Loss-of-function fd1 mutants showed chlorosis and seedling lethality at the three-leaf stage. The deficiency of Fd1 impaired photosynthetic electron transport, which affected carbon assimilation. Exogenous glucose treatment partially restored the mutant phenotype, suggesting that Fd1 plays an important role in photosynthetic electron transport in rice. In addition, the transcript levels of Fd-dependent genes were affected in fd1 mutants, and the trend was similar to that observed in fdc2 plants. Together, these results suggest that OsFd1 is the primary Fd in photosynthetic electron transport and carbon assimilation in rice.

iTRAQ-based proteomic analysis provides insights into the molecular mechanisms of rice formyl tetrahydrofolate deformylase in salt response.

MAIN CONCLUSION: A new molecular mechanism of tetrahydrofolate deformylase involved in the salt response presumably affects mitochondrial and chloroplast function by regulating energy metabolism and accumulation of reactive oxygen species. High salinity severely restrains plant growth and development, consequently leading to a reduction in grain yield. It is therefore critical to identify the components involved in plant salt resistance. In our previous study, we identified a rice leaf early-senescence mutant hpa1, which encodes a formyl tetrahydrofolate deformylase (Xiong et al. in Sci China Life Sci 64(5):720-738, 2021). Here, we report that HPA1 also plays a role in the salt response. To explore the molecular mechanism of HPA1 in salt resistance, we attempted to identify the differentially expressed proteins between wild type and hpa1 mutant for salinity treatment using an iTRAQ-based comparative protein quantification approach. A total of 4598 proteins were identified, of which 279 were significantly altered, including 177 up- and 102 down-regulated proteins. A functional analysis suggested that the 279 differentially expressed proteins are involved mainly in the regulation of oxidative phosphorylation, phenylpropanoid biosynthesis, photosynthesis, posttranslational modifications, protein turnover and energy metabolism. Moreover, a deficiency in HPA1 impaired chlorophyll metabolism and photosynthesis in chloroplasts and affected the electron flow of the electron transport chain in mitochondria. These changes led to abnormal energy metabolism and accumulation of reactive oxygen species, which may affect the permeability and integrity of cell membranes, leading to cell death. In addition, the results were verified by transcriptional or physiological experiments. Our results provide an insight into a new molecular mechanism of the tetrahydrofolate cycle protein formyl tetrahydrofolate deformylase, which is involved in the salt response, presumably by affecting mitochondrial and chloroplast function regulating energy metabolism and accumulation of reactive oxygen species under salt stress.

Cytokinin oxidase/dehydrogenase OsCKX11 coordinates source and sink relationship in rice by simultaneous regulation of leaf senescence and grain number.

The flag leaf and grain belong to the source and sink, respectively, of cereals, and both have a bearing on final yield. Premature leaf senescence significantly reduces the photosynthetic rate and severely lowers crop yield. Cytokinins play important roles in leaf senescence and determine grain number. Here, we characterized the roles of the rice (Oryza sativa L.) cytokinin oxidase/dehydrogenase OsCKX11 in delaying leaf senescence, increasing grain number, and coordinately regulating source and sink. OsCKX11 was predominantly expressed in the roots, leaves, and panicles and was strongly induced by abscisic acid and leaf senescence. Recombinant OsCKX11 protein catalysed the degradation of various types of cytokinins but showed preference for trans-zeatin and cis-zeatin. Cytokinin levels were significantly increased in the flag leaves of osckx11 mutant compared to those of the wild type (WT). In the osckx11 mutant, the ABA-biosynthesizing genes were down-regulated and the ABA-degrading genes were up-regulated, thereby reducing the ABA levels relative to the WT. Thus, OsCKX11 functions antagonistically between cytokinins and ABA in leaf senescence. Moreover, osckx11 presented with significantly increased branch, tiller, and grain number compared with the WT. Collectively, our findings reveal that OsCKX11 simultaneously regulates photosynthesis and grain number, which may provide new insights into leaf senescence and crop molecular breeding.

A ß-ketoacyl carrier protein reductase confers heat tolerance via the regulation of fatty acid biosynthesis and stress signaling in rice.

Heat stress is a major environmental threat affecting crop growth and productivity. However, the molecular mechanisms associated with plant responses to heat stress are poorly understood. Here, we identified a heat stress-sensitive mutant, hts1, in rice. HTS1 encodes a thylakoid membrane-localized ß-ketoacyl carrier protein reductase (KAR) involved in de novo fatty acid biosynthesis. Phylogenetic and bioinformatic analysis showed that HTS1 probably originated from streptophyte algae and is evolutionarily conserved in land plants. Thermostable HTS1 is predominantly expressed in green tissues and strongly induced by heat stress, but is less responsive to salinity, cold and drought treatments. An amino acid substitution at A254T in HTS1 causes a significant decrease in KAR enzymatic activity and, consequently, impairs fatty acid synthesis and lipid metabolism in the hts1 mutant, especially under heat stress. Compared to the wild-type, the hts1 mutant exhibited heat-induced higher H2 O2 accumulation, a larger Ca2+ influx to mesophyll cells, and more damage to membranes and chloroplasts. Also, disrupted heat stress signaling in the hts1 mutant depresses the transcriptional activation of HsfA2s and the downstream target genes. We suggest that HTS1 is critical for underpinning membrane stability, chloroplast integrity and stress signaling for heat tolerance in rice.

Combinatorial Evolution of a Terpene Synthase Gene Cluster Explains Terpene Variations in Oryza.

Terpenes are specialized metabolites ubiquitously produced by plants via the action of terpene synthases (TPSs). There are enormous variations in the types and amounts of terpenes produced by individual species. To understand the mechanisms responsible for such vast diversity, here we investigated the origin and evolution of a cluster of tandemly arrayed TPS genes in Oryza In the Oryza species analyzed, TPS genes occur as a three-TPS cluster, a two-TPS cluster, and a single TPS gene in five, one, and one species, respectively. Phylogenetic analysis revealed the origins of the two-TPS and three-TPS clusters and the role of species-specific losses of TPS genes. Within the three-TPS clusters, one orthologous group exhibited conserved catalytic activities. The other two groups, both of which contained pseudogenes and/or nonfunctional genes, exhibited distinct profiles of terpene products. Sequence and structural analyses combined with functional validation identified several amino acids in the active site that are critical for catalytic activity divergence of the three orthologous groups. In the five Oryza species containing the three-TPS cluster, their functional TPS genes showed both conserved and species-specific expression patterns in insect-damaged and untreated plants. Emission patterns of volatile terpenes from each species were largely consistent with the expression of their respective TPS genes and the catalytic activities of the encoded enzymes. This study indicates the importance of combinatorial evolution of TPS genes in determining terpene variations among individual species, which includes gene duplication, retention/loss/degradation of duplicated genes, varying selection pressure, retention/divergence in catalytic activities, and divergence in expression regulation.