An update of DIVERGE software for functional divergence analysis of protein family

Associating protein sequence positions with the modulation of quantitative phenotypes

Although protein sequences encode the information for folding and function, understanding their link is not an easy task. Unluckily, the prediction of how specific amino acids contribute to these features is still considerably impaired. Here, we developed a simple algorithm that finds positions in a protein sequence with potential to modulate the studied quantitative phenotypes. From a few hundred protein sequences, we perform multiple sequence alignments, obtain the per-position pairwise differences for both the sequence and the observed phenotypes, and calculate the correlation between these last two quantities. We tested our methodology with four cases: archaeal Adenylate Kinases and the organisms optimal growth temperatures, microbial rhodopsins and their maximal absorption wavelengths, mammalian myoglobins and their muscular concentration, and inhibition of HIV protease clinical isolates by two different molecules. We found from 3 to 10 positions tightly associated with those phenotypes, depending on the studied case. We showed that these correlations appear using individual positions but an improvement is achieved when the most correlated positions are jointly analyzed. Noteworthy, we performed phenotype predictions using a simple linear model that links per-position divergences and differences in the observed phenotypes. Predictions are comparable to the state-of-art methodologies which, in most of the cases, are far more complex. All of the calculations are obtained at a very low information cost since the only input needed is a multiple sequence alignment of protein sequences with their associated quantitative phenotypes. The diversity of the explored systems makes our work a valuable tool to find sequence determinants of biological activity modulation and to predict various functional features for uncharacterized members of a protein family.

Haematococcus lacustris Carotenogensis: A Historical Event of Primary to Secondary Adaptations to Earth's Oxygenation

(1) Background: Oxygen has exerted a great effect in shaping the environment and driving biological diversity in Earth's history. Green lineage has evolved primary and secondary carotenoid biosynthetic systems to adapt to Earth's oxygenation, e.g., Haematococcus lacustris, which accumulates the highest amount of secondary astaxanthin under stresses. The two systems are controlled by lycopene ε-cyclase (LCYE) and β-cyclase (LCYB), which leave an important trace in Earth's oxygenation. (2) Objectives: This work intends to disclose the underlying molecular evolutionary mechanism of Earth's oxygenation in shaping green algal carotenogensis with a special focus on lycopene cyclases. (3) Methods: The two kinds of cyclases were analyzed by site-directed mutagenesis, phylogeny, divergence time and functional divergence. (4) Results: Green lineage LCYEs appeared at ~1.5 Ga after the first significant appearance and accumulation of atmospheric oxygen, the so-called Great Oxygenation Event (GOE), from which LCYBs diverged by gene duplication. Bacterial β-bicyclases evolved from β-monocyclase. Enhanced catalytic activity accompanied evolutionary transformation from ε-/β-monocyclase to β-bicyclase. Strong positive selection occurred in green lineage LCYEs after the GOE and in algal LCYBs during the second oxidation, the Neoproterozoic Oxygenation Event (NOE). Positively selected sites in the catalytic cavities of the enzymes controlled the mono-/bicyclase activity, respectively. Carotenoid profiling revealed that oxidative adaptation has been wildly preserved in evolution. (5) Conclusions: the functionalization of the two enzymes is a result of primary to secondary adaptations to Earth's oxygenation.

Genome-wide analysis of plant specific YABBY transcription factor gene family in carrot (Dacus carota) and its comparison with Arabidopsis

YABBY gene family is a plant-specific transcription factor with DNA binding domain involved in various functions i.e. regulation of style, length of flowers, and polarity development of lateral organs in flowering plants. Computational methods were utilized to identify members of the YABBY gene family, with Carrot (Daucus carota) 's genome as a foundational reference. The structure of genes, location of the chromosomes, protein motifs and phylogenetic investigation, syntony and transcriptomic analysis, and miRNA targets were analyzed to unmask the hidden structural and functional characteristics YABBY gene family in Carrots. In the following research, it has been concluded that 11 specific YABBY genes irregularly dispersed on all 9 chromosomes and proteins assembled into five subgroups i.e. AtINO, AtCRC, AtYAB5, AtAFO, and AtYAB2, which were created on the well-known classification of Arabidopsis. The wide ranges of YABBY genes in carrots were dispersed due to segmental duplication, which was detected as prevalent when equated to tandem duplication. Transcriptomic analysis showed that one of the DcYABBY genes was highly expressed during anthocyanin pigmentation in carrot taproots. The cis-regulatory elements (CREs) analysis unveiled elements that particularly respond to light, cell cycle regulation, drought induce ability, ABA hormone, seed, and meristem expression. Furthermore, a relative study among Carrot and Arabidopsis genes of the YABBY family indicated 5 sub-families sharing common characteristics. The comprehensive evaluation of YABBY genes in the genome provides a direction for the cloning and understanding of their functional properties in carrots. Our investigations revealed genome-wide distribution and role of YABBY genes in the carrots with best-fit comparison to Arabidopsis thaliana.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

Rapid evolution of genes with anti-cancer functions during the origins of large bodies and cancer resistance in elephants

Elephants have emerged as a model system to study the evolution of body size and cancer resistance because, despite their immense size, they have a very low prevalence of cancer. Previous studies have found that duplication of tumor suppressors at least partly contributes to the evolution of anti-cancer cellular phenotypes in elephants. Still, many other mechanisms must have contributed to their augmented cancer resistance. Here, we use a suite of codon-based maximum-likelihood methods and a dataset of 13,310 protein-coding gene alignments from 261 Eutherian mammals to identify positively selected and rapidly evolving elephant genes. We found 496 genes (3.73% of alignments tested) with statistically significant evidence for positive selection and 660 genes (4.96% of alignments tested) that likely evolved rapidly in elephants. Positively selected and rapidly evolving genes are statistically enriched in gene ontology terms and biological pathways related to regulated cell death mechanisms, DNA damage repair, cell cycle regulation, epidermal growth factor receptor (EGFR) signaling, and immune functions, particularly neutrophil granules and degranulation. All of these biological factors are plausibly related to the evolution of cancer resistance. Thus, these positively selected and rapidly evolving genes are promising candidates for genes contributing to elephant-specific traits, including the evolution of molecular and cellular characteristics that enhance cancer resistance.

Genome-wide in-silico analysis of ethylene biosynthesis gene family in Musa acuminata L. and their response under nutrient stress

Ethylene is a gaseous phytohormone involved in plants' growth and developmental processes, including seed germination, root initiation, fruit ripening, flower and leaf senescence, abscission, and stress responses. Ethylene biosynthesis (EB) gene analysis in response to nitrogen (N) and potassium (K) stress has not yet been conducted in Musa acuminata (banana) roots. The genome mining of banana (Musa acuminata L.) revealed 14 putative 1-aminocyclopropane-1-carboxylate synthase (ACS), 10 1-aminocyclopropane-1-carboxylate oxidase (ACO), and 3 Ethylene overproducer 1 (ETO1) genes. ACS, ACO, and ETO1 proteins possessed amino acid residues ranging from 422-684, 636-2670, and 893-969, respectively, with molecular weight (Mw) ranging from 4.93-7.55 kD, 10.1-8.3 kD and 10.1-10.78 kD. The number of introns present in ACS, ACO, and ETO1 gene sequences ranges from 0-14, 1-6, and 0-6, respectively. The cis-regulatory element analysis revealed the presence of light-responsive, abscisic acid, seed regulation, auxin-responsive, gibberellin element, endosperm-specific, anoxic inducibility, low-temperature responsiveness, salicylic acid responsiveness, meristem-specific and stress-responsive elements. Comprehensive phylogenetic analyses ACS, ACO, and ETO1 genes of Banana with Arabidopsis thaliana revealed several orthologs and paralogs assisting in understanding the putative functions of these genes. The expression profile of Musa acuminata genes in root under normal and low levels of nitrogen and potassium shows that MaACS14 and MaACO6 expressed highly at normal nitrogen supply. MaACS1 expression was significantly upregulated at low potassium levels, whereas, MaACO6 gene expression was significantly downregulated. The functional divergence and site-specific selective pressures on specific gene sequences of banana have been investigated. The bioinformatics-based genome-wide assessment of the family of banana attempted in the present study could be a significant step for deciphering novel ACS, ACO, and ETO1 genes based on genome-wide expression profiling.

Parallel evolution of picobirnaviruses from distinct ancestral origins

IMPORTANCE

Picobirnaviruses (PBVs) are highly heterogeneous viruses encoding a capsid and RdRp. Detected in a wide variety of animals with and without disease, their association with gastrointestinal and respiratory infections, and consequently their public health importance, has rightly been questioned. Determining the "true" host of Picobirnavirus lies at the center of this debate, as evidence exists for them having both vertebrate and prokaryotic origins. Using integrated and time-stamped phylogenetic approaches, we show they are contemporaneous viruses descending from two different ancestors: avian Reovirus and fungal Partitivirus. The fungal PBV-R2 species emerged with a single segment (RdRp) until it acquired a capsid from vertebrate PBV-R1 and PBV-R3 species. Protein and RNA folding analyses revealed how the former came to resemble the latter over time. Thus, parallel evolution from disparate hosts has driven the adaptation and genetic diversification of the Picobirnaviridae family.

HDAC1: An Essential and Conserved Member of the Diverse Zn2+-Dependent HDAC Family Driven by Divergent Selection Pressure

Zn2+-dependent histone deacetylases (HDACs) are enzymes that regulate gene expression by removing acetyl groups from histone proteins. These enzymes are essential in all living systems, playing key roles in cancer treatment and as potential pesticide targets. Previous phylogenetic analyses of HDAC in certain species have been published. However, their classification and evolutionary origins across biological kingdoms remain unclear, which limits our understanding of them. In this study, we collected the HDAC sequences from 1451 organisms and performed analyses. The HDACs are found to diverge into three classes and seven subclasses under divergent selection pressure. Most subclasses show species specificity, indicating that HDACs have evolved with high plasticity and diversification to adapt to different environmental conditions in different species. In contrast, HDAC1 and HDAC3, belonging to the oldest class, are conserved and crucial in major kingdoms of life, especially HDAC1. These findings lay the groundwork for the future application of HDACs.

Slc11 Synapomorphy: A Conserved 3D Framework Articulating Carrier Conformation Switch

Transmembrane carriers of the Slc11 family catalyze proton (H+)-dependent uptake of divalent metal ions (Me2+) such as manganese and iron-vital elements coveted during infection. The Slc11 mechanism of high-affinity Me2+ cell import is selective and conserved between prokaryotic (MntH) and eukaryotic (Nramp) homologs, though processes coupling the use of the proton motive force to Me2+ uptake evolved repeatedly. Adding bacterial piracy of Nramp genes spread in distinct environmental niches suggests selective gain of function that may benefit opportunistic pathogens. To better understand Slc11 evolution, Alphafold (AF2)/Colabfold (CF) 3D predictions for bacterial sequences from sister clades of eukaryotic descent (MCb and MCg) were compared using both native and mutant templates. AF2/CF model an array of native MCb intermediates spanning the transition from outwardly open (OO) to inwardly open (IO) carriers. In silico mutagenesis targeting (i) a set of (evolutionarily coupled) sites that may define Slc11 function (putative synapomorphy) and (ii) residues from networked communities evolving during MCb transition indicates that Slc11 synapomorphy primarily instructs a Me2+-selective conformation switch which unlocks carrier inner gate and contributes to Me2+ binding site occlusion and outer gate locking. Inner gate opening apparently proceeds from interaction between transmembrane helix (h) h5, h8 and h1a. MCg1 xenologs revealed marked differences in carrier shape and plasticity, owing partly to an altered intramolecular H+ network. Yet, targeting Slc11 synapomorphy also converted MCg1 IO models to an OO state, apparently mobilizing the same residues to control gates. But MCg1 response to mutagenesis differed, with extensive divergence within this clade correlating with MCb-like modeling properties. Notably, MCg1 divergent epistasis marks the emergence of the genus Bordetella-Achromobacter. Slc11 synapomorphy localizes to the 3D areas that deviate least among MCb and MCg1 models (either IO or OO) implying that it constitutes a 3D network of residues articulating a Me2+-selective carrier conformation switch which is maintained in fast-evolving clades at the cost of divergent epistatic interactions impacting carrier shape and dynamics.

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.

Identification of key amino acid residues in AtUMAMIT29 for transport of glucosinolates

Glucosinolates are key defense compounds of plants in Brassicales order, and their accumulation in seeds is essential for the protection of the next generation. Recently, members of the Usually Multiple Amino acids Move In and Out Transporter (UMAMIT) family were shown to be essential for facilitating transport of seed-bound glucosinolates from site of synthesis within the reproductive organ to seeds. Here, we set out to identify amino acid residues responsible for glucosinolate transport activity of the main seed glucosinolate exporter UMAMIT29 in Arabidopsis thaliana. Based on a predicted model of UMAMIT29, we propose that the substrate transporting cavity consists of 51 residues, of which four are highly conserved residues across all the analyzed homologs of UMAMIT29. A comparison of the putative substrate binding site of homologs within the brassicaceous-specific, glucosinolate-transporting clade with the non-brassicaceous-specific, non-glucosinolate-transporting UMAMIT32 clade identified 11 differentially conserved sites. When each of the 11 residues of UMAMIT29 was individually mutated into the corresponding residue in UMAMIT32, five mutant variants (UMAMIT29#V27F, UMAMIT29#M86V, UMAMIT29#L109V, UMAMIT29#Q263S, and UMAMIT29#T267Y) reduced glucosinolate transport activity over 75% compared to wild-type UMAMIT29. This suggests that these residues are key for UMAMIT29-mediated glucosinolate transport activity and thus potential targets for blocking the transport of glucosinolates to the seeds.

Genetic diversity of Toll-like receptor 9 in swamp eels (Monopterus albus)

The swamp eel, Monopterus albus, is an important aquaculture species in Asia (mainly China) whose production has seriously suffered from infectious diseases. In spite of the critical requirement for aquaculture practices, to date there is scant information on its immune defence. Here, the genetic characteristics of Toll-like receptor 9 (TLR9), which plays crucial roles in the initiation of host defence against microbial invasion, were analysed. It exhibits a striking lack of genetic variation resulting from a recent demographic bottleneck. A comparison with the homologue of M. javanensis revealed that replacement but not silent differences have nonrandomly accumulated in the coding sequences at the early stage following their split from a common ancestor. Furthermore, the replacements relevant to the type II functional divergence have mainly occurred in structural motifs mediating ligand recognition and receptor homodimerization. These results provide hints to understand the diversity-based strategy of TLR9 in the arms race against pathogens. Furthermore, the findings reported here give credence to the importance of basic immunology knowledge, especially for the key elements, in genetic engineering and breeding for disease resistance in the eel and other fishes.

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.

Large-scale analyses of angiosperm Flowering Locus T genes reveal duplication and functional divergence in monocots

FLOWERING LOCUS T (FT) are well-known key genes for initiating flowering in plants. Delineating the evolutionary history and functional diversity of FT genes is important for understanding the diversification of flowering time and how plants adapt to the changing surroundings. We performed a comprehensive phylogenetic analysis of FT genes in 47 sequenced flowering plants and the 1,000 Plant Transcriptomes (1KP) database with a focus on monocots, especially cereals. We revealed the evolutionary history of FT genes. The FT genes in monocots can be divided into three clades (I, II, and III), whereas only one monophyletic group was detected in early angiosperms, magnoliids, and eudicots. Multiple rounds of whole-genome duplications (WGD) events followed by gene retention contributed to the expansion and variation of FT genes in monocots. Amino acid sites in the clade II and III genes were preferentially under high positive selection, and some sites located in vital domain regions are known to change functions when mutated. Clade II and clade III genes exhibited high variability in important regions and functional divergence compared with clade I genes; thus, clade I is more conserved than clade II and III. Genes in clade I displayed higher expression levels in studied organs and tissues than the clade II and III genes. The co-expression modules showed that some of the FT genes might have experienced neofunctionalization and subfunctionalization, such as the acquisition of environmental resistance. Overall, FT genes in monocots might form three clades by the ancient gene duplication, and each clade was subsequently subjected to different selection pressures and amino acid substitutions, which eventually led to different expression patterns and functional diversification. Our study provides a global picture of FT genes' evolution in monocots, paving a road for investigating FT genes' function in future.

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.

Birth-and-death evolution of ribonuclease 9 genes in Cetartiodactyla

RNase9 plays a reproductive function and has been recognized as an important member of the ribonuclease (RNase) A superfamily, a gene family that is widely used as a model for molecular evolutionary studies. Here, we identified 178 RNase9 genes from 95 Cetartiodactyla species that represent all four lineages and 21 families of this clade. Unexpectedly, RNase9 experienced an evolutionary scenario of "birth and death" in Ruminantia, and expression analyses showed that duplicated RNase9A and RNase9B genes are expressed in reproductive tissues (epididymis, vas deferens or prostate). This expression pattern combined with the estimate that these genes duplicated during the middle Eocene, a time when Ruminantia become a successful lineage, suggests that the RNase9 gene duplication might have been advantageous for promoting sperm motility and male fertility as an adaptation to climate seasonality changes of this period. In contrast, all RNase9 genes were lost in the Cetacean lineage, which might be associated with their high levels of prostatic lesions and lower reproductive rates as adaptations to a fully aquatic environment and a balance to the demands of ocean resources. This study reveals a complex and intriguing evolutionary history and functional divergence for RNase9 in Cetartiodactyla, providing new insights into the evolution of the RNaseA superfamily and molecular mechanisms for organismal adaptations to the environment.

Nramp: Deprive and conquer?

Solute carriers 11 (Slc11) evolved from bacterial permease (MntH) to eukaryotic antibacterial defense (Nramp) while continuously mediating proton (H⁺)-dependent manganese (Mn²⁺) import. Also, Nramp horizontal gene transfer (HGT) toward bacteria led to mntH polyphyly. Prior demonstration that evolutionary rate-shifts distinguishing Slc11 from outgroup carriers dictate catalytic specificity suggested that resolving Slc11 family tree may provide a function-aware phylogenetic framework. Hence, MntH C (MC) subgroups resulted from HGTs of prototype Nramp (pNs) parologs while archetype Nramp (aNs) correlated with phagocytosis. PHI-Blast based taxonomic profiling confirmed MntH B phylogroup is confined to anaerobic bacteria vs. MntH A (MA)’s broad distribution; suggested niche-related spread of MC subgroups; established that MA-variant MH, which carries ‘eukaryotic signature’ marks, predominates in archaea. Slc11 phylogeny shows MH is sister to Nramp. Site-specific analysis of Slc11 charge network known to interact with the protonmotive force demonstrates sequential rate-shifts that recapitulate Slc11 evolution. 3D mapping of similarly coevolved sites across Slc11 hydrophobic core revealed successive targeting of discrete areas. The data imply that pN HGT could advantage recipient bacteria for H⁺-dependent Mn²⁺ acquisition and Alphafold 3D models suggest conformational divergence among MC subgroups. It is proposed that Slc11 originated as a bacterial stress resistance function allowing Mn²⁺-dependent persistence in conditions adverse for growth, and that archaeal MH could contribute to eukaryogenesis as a Mn²⁺ sequestering defense perhaps favoring intracellular growth-competent bacteria.

Genome-wide identification of wheat ABC1K gene family and functional dissection of TaABC1K3 and TaABC1K6 involved in drought tolerance

Activity of BC1 complex kinase (ABC1K) serves as an atypical kinase family involved in plant stress resistance. This study identified 44 ABC1K genes in the wheat genome, which contained three clades (I-III). TaABC1K genes generally had similar structural features, but differences were present in motif and exon compositions from different clade members. More type II functional divergence sites were detected between clade I and clade III and no positive selection site were found in TaABC1K family. The three-dimensional structure prediction by Alphafold2 showed that TaABC1K proteins had more α-helixes with a relatively even distribution, and different clade members had differences in the content of secondary structures. The cis-acting element analysis showed that TaABC1K genes contained abundant cis-acting elements related to plant hormones and environmental stress response in the promoter region, and generally displayed a significantly upregulated expression under drought stress. In particular, both TaABC1K3 and TaABC1K6 genes from clade I was highly induced by drought stress, and their overexpression in yeast and Arabidopsis enhanced drought tolerance by suppressing active oxygen burst and reducing photosynthesis impairment. Meanwhile, TaABC1K3 and TaABC1K6 could, respectively, complement the function of Arabidopsis abc1k3 and abc1k6 mutants and reduce photosynthesis damage caused by drought stress.

Biochemical Evolution of a Potent Target of Mosquito Larvicide, 3-Hydroxykynurenine Transaminase

A specific mosquito enzyme, 3-hydroxykynurenine transaminase (HKT), is involved in the processing of toxic metabolic intermediates of the tryptophan metabolic pathway. The HKT enzymatic product, xanthurenic acid, is required for Plasmodium spp. development in the mosquito vectors. Therefore, an inhibitor of HKT may not only be a mosquitocide but also a malaria-transmission blocker. In this work, we present a study investigating the evolution of HKT, which is a lineage-specific duplication of an alanine glyoxylate aminotransferases (AGT) in mosquitoes. Synteny analyses, together with the phylogenetic history of the AGT family, suggests that HKT and the mosquito AGTs are paralogous that were formed via a duplication event in their common ancestor. Furthermore, 41 amino acid sites with significant evidence of positive selection were identified, which could be responsible for biochemical and functional evolution and the stability of conformational stabilization. To get a deeper understanding of the evolution of ligands’ capacity and the ligand-binding mechanism of HKT, the sequence and the 3D homology model of the common ancestor of HKT and AGT in mosquitoes, ancestral mosquito AGT (AncMosqAGT), were inferred and built. The homology model along with 3-hydroxykynurenine, kynurenine, and alanine were used in docking experiments to predict the binding capacity and ligand-binding mode of the new substrates related to toxic metabolites detoxification. Our study provides evidence for the dramatic biochemical evolution of the key detoxifying enzyme and provides potential sites that could hinder the detoxification function, which may be used in mosquito larvicide and design.

Evolution of osmosensing OSCA1 Ca2+ channel family coincident with plant transition from water to land

Water is crucial to plant growth, development, and environmental adaptation. Water stress triggers cytosolic Ca²⁺ ([Ca²⁺ ]_i ) increases, and the osmosensor OSCA1 (REDUCED-HYPEROSMOLALITY-INDUCED-[Ca²⁺ ]_i -INCREASE 1), a member of the OSCA family, perceives the initial water stress and governs its downstream responses. OSCA homologs exist in eukaryotes and largely radiate in higher plants. However, it is enigmatic whether the OSCA family is crucial for plant evolution from aqueous to terrestrial environments and for the subsequent adaptation on land. Here, we carried out the first phylogenetic and molecular evolutionary analyses of the OSCA family. The family originated and diversified during the early evolution of protists, and three more lineages were established (a) in plants, (b) in fungi, and (c) in a complex clade of several major eukaryotic lineages. The chlorophyte algal cluster is directly basal to streptophyte-specific Clades 1-3, consistent with plant transition from water to land. The Clades 1-3 present different gene expansion pattern and together with previous functional analysis of OSCAs reveal that they probably have evolved diverse functions in respond to various mechanical stresses during the independent evolution of land plant clades. Moreover, variable selection pressures on different land plant lineages were explored. OSCAs in early land plants (mosses and lycophytes) were under decelerated evolution, whereas OSCAs in seed plants showed accelerated evolution. Together, we hypothesize OSCAs have evolved to sense water stress in the ancestor of euphyllophytes, which occupies typical leaves, typical roots, and phloem tissues, all of which require osmosensors to maintain water balance and food conduction through plant bodies.

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.

Inter-paralog amino acid inversion events in large phylogenies of duplicated proteins

Connecting protein sequence to function is becoming increasingly relevant since high-throughput sequencing studies accumulate large amounts of genomic data. In order to go beyond the existing database annotation, it is fundamental to understand the mechanisms underlying functional inheritance and divergence. If the homology relationship between proteins is known, can we determine whether the function diverged? In this work, we analyze different possibilities of protein sequence evolution after gene duplication and identify “inter-paralog inversions”, i.e., sites where the relationship between the ancestry and the functional signal is decoupled. The amino acids in these sites are masked from being recognized by other prediction tools. Still, they play a role in functional divergence and could indicate a shift in protein function. We develop a method to specifically recognize inter-paralog amino acid inversions in a phylogeny and test it on real and simulated datasets. In a dataset built from the Epidermal Growth Factor Receptor (EGFR) sequences found in 88 fish species, we identify 19 amino acid sites that went through inversion after gene duplication, mostly located at the ligand-binding extracellular domain. Our work uncovers an outcome of protein duplications with direct implications in protein functional annotation and sequence evolution. The developed method is optimized to work with large protein datasets and can be readily included in a targeted protein analysis pipeline.

Proteins are critical components of living systems because they facilitate most biological processes like protein synthesis, DNA replication, chemical catalysis, etc. Proteins are encoded in their genes. During evolution, genes accumulate mutations that get translated at the protein level. These mutations can be “neutral” if they do not affect the protein function immediately and directly; otherwise, mutations can be functional if they directly modify protein function. An event that provides an opportunity to study protein function is gene duplication namely, when two copies of a gene encoding the same protein appear. One copy of the protein often retains the same function while the other is free to diverge and specialize to a different function. This work sheds light on an alternative outcome of gene duplication that might be critical to discern between neutral and functional mutations. By looking at 88 fish genomes, we found proteins in which the evolution of their sequences does not follow the expected pattern of divergence after gene duplication. In this case, the protein sequence of a subgroup of species diverges in the copy expected to retain its function, while the sequence is retained in the expectedly divergent one. We called this event “inter-paralog amino acid inversion”. Our data shows that this “inversion” event is correlated to function, and its detection has to be considered for assigning protein functions correctly.

Genome-Wide Analysis of Serine Hydroxymethyltransferase Genes in Triticeae Species Reveals That TaSHMT3A-1 Regulates Fusarium Head Blight Resistance in Wheat

Serine hydroxymethyltransferase (SHMT) plays a pivotal role in cellular one-carbon, photorespiration pathways and it influences the resistance to biotic and abiotic stresses. However, the function of SHMT proteins in wheat remains largely unexplored. In the present study, SHMT genes in five Triticeae species, Oryza sativa, and four dicotyledon species were identified based on whole genome information. The origin history of the target gene was traced by micro-collinearity analysis. Gene expression patterns of TaSHMTs in different tissues, various biotic stresses, exogenous hormones, and two biotic stresses were determined by Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The function of the selected TaSHMT3A-1 was studied by barley stripe mosaic virus-induced gene silencing in common wheat Bainong207. A total of 64 SHMT members were identified and further classified into two main classes based on the structure of SHMT proteins. The gene structure and motif composition analyses revealed that SHMTs kept relatively conserved within the same subclasses. Interestingly, there was a gene, TdSHMT7B-1, on chromosome 7B of Triticum dicoccoides, but there was no SHMT gene on chromosome 7 of other analyzed Triticeae species; TdSHMT7B-1 had fewer exons and conserved motifs than the genes in the same subclass, suggesting that the gene of TdSHMT7B-1 has a notable evolutionary progress. The micro-collinearity relationship showed that no homologs of TaSHMT3A-1 and its two neighboring genes were found in the collinearity region of Triticum urartu, and there were 27 genes inserted into the collinearity region of T. urartu. Furthermore, qRT-PCR results showed that TaSHMT3A-1 was responsive to abiotic stresses (NaCl and cold), abscisic acid, methyl jasmonate, and hydrogen peroxide. Significantly, upon Fusarium graminearum infection, the expression of TaSHMT3A-1 was highly upregulated in resistant cultivar Sumai3. More importantly, silencing of TaSHMT3A-1 compromises Fusarium head blight resistance in common wheat Bainong207. Our new findings suggest that the TaSHMT3A-1 gene in wheat plays an important role in resistance to Fusarium head blight. This provides a valuable reference for further study on the function of this gene family.

Evolutionary Analysis of Dipeptidyl Peptidase I

Human dipeptidyl peptidase I (DPPI) belongs to the family of papain-like cysteine peptidases. Its distinctive features are the unique exclusion domain which enables the eponymous activity and homotetramerization of DPPI, and its dependence on chloride ions for enzymatic activity. The oligomeric state of DPPI is unique in this family of predominantly monomeric peptidases. However, a distant DPPI ortholog from Plasmodium falciparum has been shown to be monomeric, indicating that the oligomeric state of DPPI varies between lineages. The aim of this work was to study the evolution of DPPI, with particular attention to the structural features that determine its characteristic enzymatic activity and preferences, and to reconstruct the evolution of its oligomerization. We analyzed fifty-seven selected sequences of DPPI and confirmed its presence in three lineages, namely, Amorphea (including animals and Amoebozoa), Alveolates and the metamonad Giardia. The amino acid residues that bind the chloride ion are highly conserved in all species, indicating that the dependence on chloride ions for activity is an evolutionarily conserved feature of DPPI. The number of N-glycosylation sites is significantly increased in animals, particularly vertebrates. Analysis of homology models and subunit contacts suggests that oligomerization is likely restricted to DPPIs in the Amorphea group.

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.

The Hair Cell α9α10 Nicotinic Acetylcholine Receptor: Odd Cousin in an Old Family

Nicotinic acetylcholine receptors (nAChRs) are a subfamily of pentameric ligand-gated ion channels with members identified in most eumetazoan clades. In vertebrates, they are divided into three subgroups, according to their main tissue of expression: neuronal, muscle and hair cell nAChRs. Each receptor subtype is composed of different subunits, encoded by paralogous genes. The latest to be identified are the α9 and α10 subunits, expressed in the mechanosensory hair cells of the inner ear and the lateral line, where they mediate efferent modulation. α9α10 nAChRs are the most divergent amongst all nicotinic receptors, showing marked differences in their degree of sequence conservation, their expression pattern, their subunit co-assembly rules and, most importantly, their functional properties. Here, we review recent advances in the understanding of the structure and evolution of nAChRs. We discuss the functional consequences of sequence divergence and conservation, with special emphasis on the hair cell α9α10 receptor, a seemingly distant cousin of neuronal and muscle nicotinic receptors. Finally, we highlight potential links between the evolution of the octavolateral system and the extreme divergence of vertebrate α9α10 receptors.

Genome-Wide Analysis of Gene Families of Pattern Recognition Receptors in Fig Wasps (Hymenoptera, Chalcidoidea)

Pattern recognition receptors (PRRs) play important roles in detecting pathogens and initiating the innate immune response. Different evolutionary histories of pollinators and non-pollinators may result in different immune recognition systems. A previous study had reported that there were significant differences in peptidoglycan recognition proteins (PGRPs) between pollinators and non-pollinators in gene number and lineage of specific genes. In this study, based on the genomic data of 12 fig wasp species, with seven pollinators and five non-pollinators, we investigated the evolution patterns of PRRs, such as Gram-negative bacteria-binding proteins (GNBPs), C-type lectins (CTLs), scavenger receptors class B (SCRBs), fibrinogen-related proteins (FREPs), galectins, and thioester-containing proteins (TEPs). Our results showed that pollinators had no GNBP, but non-pollinators all had two gene members, which were clustered into two different clades in the phylogenetic tree, with each clade having specific domain and motif characteristics. The analysis of CTL and SCRB gene families also showed that there were lineage-specific genes and specific expansion in non-pollinators. Our results showed that there were significant differences in immune recognition between pollinators and non-pollinators, and we concluded that they had undergone flexible adaptive evolution in different environments. Our study can provide more molecular evidence for future functional studies on the immune system of fig wasps.

TwinCons: Conservation score for uncovering deep sequence similarity and divergence

We have developed the program TwinCons, to detect noisy signals of deep ancestry of proteins or nucleic acids. As input, the program uses a composite alignment containing pre-defined groups, and mathematically determines a 'cost' of transforming one group to the other at each position of the alignment. The output distinguishes conserved, variable and signature positions. A signature is conserved within groups but differs between groups. The method automatically detects continuous characteristic stretches (segments) within alignments. TwinCons provides a convenient representation of conserved, variable and signature positions as a single score, enabling the structural mapping and visualization of these characteristics. Structure is more conserved than sequence. TwinCons highlights alternative sequences of conserved structures. Using TwinCons, we detected highly similar segments between proteins from the translation and transcription systems. TwinCons detects conserved residues within regions of high functional importance for the ribosomal RNA (rRNA) and demonstrates that signatures are not confined to specific regions but are distributed across the rRNA structure. The ability to evaluate both nucleic acid and protein alignments allows TwinCons to be used in combined sequence and structural analysis of signatures and conservation in rRNA and in ribosomal proteins (rProteins). TwinCons detects a strong sequence conservation signal between bacterial and archaeal rProteins related by circular permutation. This conserved sequence is structurally colocalized with conserved rRNA, indicated by TwinCons scores of rRNA alignments of bacterial and archaeal groups. This combined analysis revealed deep co-evolution of rRNA and rProtein buried within the deepest branching points in the tree of life.

Genome and transcriptome analysis to understand the role diversification of cytochrome P450 gene under excess nitrogen treatment

BACKGROUND

Panax notoginseng (Burk.) F. H. Chen (P. notoginseng) is a medicinal plant. Cytochrome P450 (CYP450) monooxygenase superfamily is involved in the synthesis of a variety of plant hormones. Studies have shown that CYP450 is involved in the synthesis of saponins, which are the main medicinal component of P. notoginseng. To date, the P. notoginseng CYP450 family has not been systematically studied, and its gene functions remain unclear.

RESULTS

In this study, a total of 188 PnCYP genes were identified, these genes were divided into 41 subfamilies and clustered into 9 clans. Moreover, we identified 40 paralogous pairs, of which only two had Ka/Ks ratio greater than 1, demonstrating that most PnCYPs underwent purification selection during evolution. In chromosome mapping and gene replication analysis, 8 tandem duplication and 11 segmental duplication events demonstrated that PnCYP genes were continuously replicating during their evolution. Gene ontology (GO) analysis annotated the functions of 188 PnCYPs into 21 functional subclasses, suggesting the functional diversity of these gene families. Functional divergence analyzed the members of the three primitive branches of CYP51, CYP74 and CYP97 at the amino acid level, and found some critical amino acid sites. The expression pattern of PnCYP450 related to nitrogen treatment was studied using transcriptome sequencing data, 10 genes were significantly up-regulated and 37 genes were significantly down-regulated. Combined with transcriptome sequencing analysis, five potential functional genes were screened. Quantitative real-time PCR (qRT-PCR) indicated that these five genes were responded to methyl jasmonate (MEJA) and abscisic acid (ABA) treatment.

CONCLUSIONS

These results provide a valuable basis for comprehending the classification and biological functions of PnCYPs, and offer clues to study their biological functions in response to nitrogen treatment.

An Evolutionary Perspective of the Lipocalin Protein Family

The protein family of Lipocalins is ubiquitously present throughout the tree of life, with the exception of the phylum Archaea. Phylogenetic relationships of chordate Lipocalins have been proposed in the past based on protein sequence similarities, but their highly divergent primary structures and a shortage of experimental annotations in genome projects have precluded a well-supported hypothesis for their evolution. In this work we propose a novel topology for the phylogenetic tree of chordate Lipocalins, inferred from multiple amino acid sequence alignments. Sixteen jawed vertebrates with fair coverage by genomic sequencing were compared. The selected species span an evolutionary range of ∼400 million years, allowing for a balanced representation of all major vertebrate clades. A consensus phylogenetic tree is proposed following a comparison of sequence-based maximum-likelihood trees and protein structure dendrograms. This new phylogeny suggests an APOD-like common ancestor in early chordates, which gave rise, via whole-genome or tandem duplications, to the six Lipocalins currently present in fish (APOD, RBP4, PTGDS, AMBP, C8G, and APOM). Further gene duplications of APOM and PTGDS resulted in the altogether 15 Lipocalins found in contemporary mammals. Insights into the functional impact of relevant amino acid residues in early diverging Lipocalins are also discussed. These results should foster the experimental exploration of novel functions alongside the identification of new members of the Lipocalin family.

Genome-wide identification and characterization of the Lateral Organ Boundaries Domain (LBD) gene family in polyploid wheat and related species

BACKGROUND

Wheat (Triticum aestivum) originated from three different diploid ancestral grass species and experienced two rounds of polyploidization. Exploring how certain wheat gene subfamilies have expanded during the evolutionary process is of great importance. The Lateral Organ Boundaries Domain (LBD) gene family encodes plant-specific transcription factors that share a highly conserved LOB domain and are prime candidates for this, as they are involved in plant growth, development, secondary metabolism and stress in various species.

METHODS

Using a genome-wide analysis of high-quality polyploid wheat and related species genome sequences, a total of 228 LBD members from five Triticeae species were identified, and phylogenetic relationship analysis of LBD members classified them into two main classes (classes I and II) and seven subgroups (classes I a-e, II a and II b).

RESULTS

The gene structure and motif composition analyses revealed that genes that had a closer phylogenetic relationship in the same subgroup also had a similar gene structure. Macrocollinearity and microcollinearity analyses of Triticeae species suggested that some LBD genes from wheat produced gene pairs across subgenomes of chromosomes 4A and 5A and that the complex evolutionary history of TaLBD4B-9 homologs was a combined result of chromosome translocation, polyploidization, gene loss and duplication events. Public RNA-seq data were used to analyze the expression patterns of wheat LBD genes in various tissues, different developmental stages and following abiotic and biotic stresses. Furthermore, qRT-PCR results suggested that some TaLBDs in class II responded to powdery mildew, regulated reproductive growth and were involved in embryo sac development in common wheat.

Important innate differences in determining symbiotic responsiveness in host and non-hosts of arbuscular mycorrhiza

Genetic components that regulate arbuscular mycorrhizal (AM) interactions in hosts and non-hosts are not completely known. Comparative transcriptomic analysis was combined with phylogenetic studies to identify the factors that distinguish AM host from non-host. Mycorrhized host, non-mycorrhized host and non-host cultivars of tomato (Solanum lycopersicum) were subjected to RNA seq analysis. The top 10 differentially expressed genes were subjected to extensive in silico phylogenetic analysis along with 10 more candidate genes that have been previously reported for AM-plant interactions. Seven distantly related hosts and four non-hosts were selected to identify structural differences in selected gene/protein candidates. The screened genes/proteins were subjected to MEME, CODEML and DIVERGE analysis to identify evolutionary patterns that differentiate hosts from non-hosts. Based on the results, candidate genes were categorized as highly influenced (SYMRK and CCaMK), moderately influenced and minimally influenced by evolutionary constraints. We propose that the amino acid and nucleotide changes specific to non-hosts are likely to correspond to aberrations in functionality towards AM symbiosis. This study paves way for future research aimed at understanding innate differences in genetic make-up of AM hosts and non-hosts, in addition to the theory of gene losses from the "AM-symbiotic toolkit".

Evolution and functional divergence of the ERBB receptor family

Aim: The ERBB gene family plays an important role in 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 five functional divergence residues sites, 94 totally conserved sites and averagely 345 conserved sites of individual gene member were obtained. Some have been reported to play role in 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.

Evolution of kaiA, a key circadian gene of cyanobacteria

The circadian system of cyanobacteria is built upon a central oscillator consisting of three genes, kaiA, kaiB, and kaiC. The KaiA protein plays a key role in phosphorylation/dephosphorylation cycles of KaiC, which occur over the 24-h period. We conducted a comprehensive evolutionary analysis of the kaiA genes across cyanobacteria. The results show that, in contrast to the previous reports, kaiA has an ancient origin and is as old as cyanobacteria. The kaiA homologs are present in nearly all analyzed cyanobacteria, except Gloeobacter, and have varying domain architecture. Some Prochlorococcales, which were previously reported to lack the kaiA gene, possess a drastically truncated homolog. The existence of the diverse kaiA homologs suggests significant variation of the circadian mechanism, which was described for the model cyanobacterium, Synechococcus elongatus PCC7942. The major structural modifications in the kaiA genes (duplications, acquisition and loss of domains) have apparently been induced by global environmental changes in the different geological periods.

Gene duplication and subsequent functional diversification of maltase in fig wasp (Chalcidoidea, Hymenoptera)

Maltase can catalyze the hydrolysis of α-1,4-glucosidic linkages and release α-d-glucoses that are used as a source of energy by insects. Maltase has been extensively studied in Lepidoptera and Diptera, while the characterization and evolutionary history of maltase are largely unknown in Hymenoptera. Here, we undertook a bioinformatics study and identified 105 maltase genes in 12 fig wasp species. Together with the maltase genes of Nasonia vitripennis and Apis mellifera, phylogenetic analysis showed that all the maltase genes were clustered into three clades. Clade I and III included maltase genes from all the fig wasp species, while clade II contained the maltase genes from non-pollinating fig wasps (NPFWs) only. Interestingly, the maltase genes located in clade II were intronless. Fig pollinators and NPFWs had lineage-specific gene expansion in clade I and II respectively, which were mainly derived from tandem duplications. The three clades displayed distinct gene structures. Furthermore, maltase showed significant functional divergence among the three clades and the critical amino acid sites were detected. These sites could be responsible for the ligand-binding preference and hydrolytic specificity. Overall, our results demonstrated that maltase might contribute to the discrepancy of life histories and feeding regimes between fig pollinators and NPFWs.

Distinct Evolutionary Trajectories of Neuronal and Hair Cell Nicotinic Acetylcholine Receptors

The expansion and pruning of ion channel families has played a crucial role in the evolution of nervous systems. Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels with distinct roles in synaptic transmission at the neuromuscular junction, the central and peripheral nervous system, and the inner ear. Remarkably, the complement of nAChR subunits has been highly conserved along vertebrate phylogeny. To ask whether the different subtypes of receptors underwent different evolutionary trajectories, we performed a comprehensive analysis of vertebrate nAChRs coding sequences, mouse single-cell expression patterns, and comparative functional properties of receptors from three representative tetrapod species. We found significant differences between hair cell and neuronal receptors that were most likely shaped by the differences in coexpression patterns and coassembly rules of component subunits. Thus, neuronal nAChRs showed high degree of coding sequence conservation, coupled to greater coexpression variance and conservation of functional properties across tetrapod clades. In contrast, hair cell α9α10 nAChRs exhibited greater sequence divergence, narrow coexpression pattern, and great variability of functional properties across species. These results point to differential substrates for random change within the family of gene paralogs that relate to the segregated roles of nAChRs in synaptic transmission.

Comparative Genomics within and across Bilaterians Illuminates the Evolutionary History of ALK and LTK Proto-Oncogene Origination and Diversification

Comparative genomic analyses have enormous potential for identifying key genes central to human health phenotypes, including those that promote cancers. In particular, the successful development of novel therapeutics using model species requires phylogenetic analyses to determine molecular homology. Accordingly, we investigate the evolutionary histories of anaplastic lymphoma kinase (ALK)—which can underlie tumorigenesis in neuroblastoma, nonsmall cell lung cancer, and anaplastic large-cell lymphoma—its close relative leukocyte tyrosine kinase (LTK) and their candidate ligands. Homology of ligands identified in model organisms to those functioning in humans remains unclear. Therefore, we searched for homologs of the human genes across metazoan genomes, finding that the candidate ligands Jeb and Hen-1 were restricted to nonvertebrate species. In contrast, the ligand augmentor (AUG) was only identified in vertebrates. We found two ALK-like and four AUG-like protein-coding genes in lamprey. Of these six genes, only one ALK-like and two AUG-like genes exhibited early embryonic expression that parallels model mammal systems. Two copies of AUG are present in nearly all jawed vertebrates. Our phylogenetic analysis strongly supports the presence of previously unrecognized functional convergences of ALK and LTK between actinopterygians and sarcopterygians—despite contemporaneous, highly conserved synteny of ALK and LTK. These findings provide critical guidance regarding the propriety of fish and mammal models with regard to model organism-based investigation of these medically important genes. In sum, our results provide the phylogenetic context necessary for effective investigations of the functional roles and biology of these critically important receptors.

A Phylogenetic Rate Parameter Indicates Different Sequence Divergence Patterns in Orthologs and Paralogs

Heterotachy-the change in sequence evolutionary rate over time-is a common feature of protein molecular evolution. Decades of studies have shed light on the conditions under which heterotachy occurs, and there is evidence that site-specific evolutionary rate shifts are correlated with changes in protein function. Here, we present a large-scale, computational analysis using thousands of protein sequence alignments from animal and plant proteomes, representing genes related either by orthology (speciation events) or paralogy (gene duplication), to compare sequence divergence patterns in orthologous vs. paralogous sequence alignments. We use sequence-based phylogenetic analyses to infer overall sequence divergence (tree length/number of sequences) and to fit site-specific rates to a discrete gamma distribution with a shape parameter α. This inference method is applied to real protein sequence alignments, as well as alignments simulated under various models of protein sequence evolution. Our simulations indicate that sequence divergence and the α parameter are positively correlated when sequences evolve with heterotachy, meaning that inferred site rate distributions appear more uniform as sequences diverge. Divergence and α are also positively correlated in both orthologous and paralogous genes, but the average increase in α (as a function of divergence) is significantly higher in paralogous protein alignments than in orthologous alignments. This result is consistent with the widely held view that recently duplicated proteins initially evolve under relaxed selective pressure, promoting functional divergence by accumulation of amino acid replacements, and hence experience more evolutionary rate fluctuations than orthologous proteins. We discuss these findings in the context of the ortholog conjecture, a long-standing assumption in molecular evolution, which posits that protein sequences related by orthology tend to be more functionally conserved than paralogous proteins.

Genome-wide identification of ATP binding cassette (ABC) transporter and heavy metal associated (HMA) gene families in flax (Linum usitatissimum L.)

Background

The recent release of the reference genome sequence assembly of flax, a self-pollinated crop with 15 chromosome pairs, into chromosome-scale pseudomolecules enables the characterization of gene families. The ABC transporter and HMA gene families are important in the control of cadmium (Cd) accumulation in crops. To date, the genome-wide analysis of these two gene families has been successfully conducted in some plant species, but no systematic evolutionary analysis is available for the flax genome.

Results

Here we describe the ABC transporter and HMA gene families in flax to provide a comprehensive overview of its evolution and some support towards the functional annotation of its members. The 198 ABC transporter and 12 HMA genes identified in the flax genome were classified into eight ABC transporter and four HMA subfamilies based on their phylogenetic analysis and domains’ composition. Nine of these genes, i.e., LuABCC9, LuABCC10, LuABCG58, LuABCG59, LuABCG71, LuABCG72, LuABCG73, LuHMA3, and LuHMA4, were orthologous with the Cd associated genes in Arabidopsis, rice and maize. Ten motifs were identified from all ABC transporter and HMA genes. Also, several motifs were conserved among genes of similar length, but each subfamily each had their own motif structures. Both the ABC transporter and HMA gene families were highly conserved among subfamilies of flax and with those of Arabidopsis. While four types of gene duplication were observed at different frequencies, whole-genome or segmental duplications were the most frequent with 162 genes, followed by 29 dispersed, 14 tandem and 4 proximal duplications, suggesting that segmental duplications contributed the most to the expansion of both gene families in flax. The rates of non-synonymous to synonymous (Ka/Ks) mutations of paired duplicated genes were for the most part lower than one, indicative of a predominant purifying selection. Only five pairs of genes clearly exhibited positive selection with a Ka/Ks ratio greater than one. Gene ontology analyses suggested that most flax ABC transporter and HMA genes had a role in ATP binding, transport, catalytic activity, ATPase activity, and metal ion binding. The RNA-Seq analysis of eight different organs demonstrated diversified expression profiling patterns of the genes and revealed their functional or sub-functional conservation and neo-functionalization.

Conclusion

Characterization of the ABC transporter and HMA gene families will help in the functional analysis of candidate genes in flax and other crop species.

Evolutionary history of histamine receptors: Early vertebrate origin and expansion of the H3-H4 subtypes

Histamine receptors belonging to the superfamily of G protein-coupled receptors (GPCRs) mediate the diverse biological effects of biogenic histamine. They are classified into four phylogenetically distinct subtypes H₁-H₄, each with a different binding affinity for histamine and divergent downstream signaling pathways. Here we present the evolutionary history of the histamine receptors using a phylogenetic approach complemented with comparative genomics analyses of the sequences, gene structures, and synteny of gene neighborhoods. The data indicate the earliest emergence of histamine-mediated GPCR signaling by a H₂ in a prebilaterian ancestor. The analyses support a revised classification of the vertebrate H₃-H₄ receptor subtypes. We demonstrate the presence of the H₄ across vertebrates, contradicting the currently held notion that H₄ is restricted to mammals. These non-mammalian vertebrate H₄ orthologs have been mistaken for H₃. We also identify the presence of a new H₃ subtype (H₃B), distinct from the canonical H₃ (H₃A), and propose that the H₃A, H₃B, and H₄ likely emerged from a H₃ progenitor through the 1R/2R whole genome duplications in an ancestor of the vertebrates. It is apparent that the ability of the H₁, H₂, and H_3-4 to bind histamine was acquired convergently. We identified genomic signatures suggesting that the H₁ and H₃-H₄ shared a last common ancestor with the muscarinic receptor in a bilaterian predecessor whereas, the H₂ and the α-adrenoreceptor shared a progenitor in a prebilaterian ancestor. Furthermore, site-specific analysis of the vertebrate subtypes revealed potential residues that may account for the functional divergence between them.

Characteristics and Expression Analyses of Trehalose-6-Phosphate Synthase Family in Prunus mume Reveal Genes Involved in Trehalose Biosynthesis and Drought Response

Trehalose and its key synthase (trehalose-6-phosphate synthase, TPS) can improve the drought tolerance of plants. However, little is known about the roles of trehalose and the TPS family in Prunus mume response to drought. In our study, we discovered that the trehalose content in leaf, root, and stem tissues significantly increased in P. mume in response to drought. Therefore, the characteristics and functions of the TPS family are worth investigating in P. mume. We identified nine TPS family members in P. mume, which were divided into two sub-families and characterized by gene structure, promoter elements, protein conserved domains, and protein motifs. We found that the Hydrolase_3 domain and several motifs were highly conserved in Group II instead of Group I. The distinctions between the two groups may result from selective constraints, which we estimated by the d_N/d_S (ω) ratio. The ω values of all the PmTPS family gene pairs were evaluated as less than 1, indicating that purity selection facilitated their divergence. A phylogenetic tree was constructed using 92 TPSs from 10 Rosaceae species, which were further divided into five clusters. Based on evolutionary analyses, the five clusters of TPS family proteins mainly underwent varied purity selection. The expression patterns of PmTPSs under drought suggested that the TPS family played an important role in the drought tolerance of P. mume. Combining the expression patterns of PmTPSs and the trehalose content changes in leaf, stem, and root tissues under normal conditions and drought stress, we found that the PmTPS2 and PmTPS6 mainly function in the trehalose biosynthesis in P. mume. Our findings not only provide valuable information about the functions of trehalose and TPSs in the drought response of P. mume, but they also contribute to the future drought breeding of P. mume.

Identification of biochemically neutral positions in liver pyruvate kinase

Understanding how each residue position contributes to protein function has been a long-standing goal in protein science. Substitution studies have historically focused on conserved protein positions. However, substitutions of nonconserved positions can also modify function. Indeed, we recently identified nonconserved positions that have large substitution effects in human liver pyruvate kinase (hLPYK), including altered allosteric coupling. To facilitate a comparison of which characteristics determine when a nonconserved position does vs does not contribute to function, the goal of the current work was to identify neutral positions in hLPYK. However, existing hLPYK data showed that three features commonly associated with neutral positions-high sequence entropy, high surface exposure, and alanine scanning-lacked the sensitivity needed to guide experimental studies. We used multiple evolutionary patterns identified in a sequence alignment of the PYK family to identify which positions were least patterned, reasoning that these were most likely to be neutral. Nine positions were tested with a total of 117 amino acid substitutions. Although exploring all potential functions is not feasible for any protein, five parameters associated with substrate/effector affinities and allosteric coupling were measured for hLPYK variants. For each position, the aggregate functional outcomes of all variants were used to quantify a "neutrality" score. Three positions showed perfect neutral scores for all five parameters. Furthermore, the nine positions showed larger neutral scores than 17 positions located near allosteric binding sites. Thus, our strategy successfully enriched the dataset for positions with neutral and modest substitutions.

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.

Rheostat positions: A new classification of protein positions relevant to pharmacogenomics

To achieve the full potential of pharmacogenomics, one must accurately predict the functional outcomes that arise from amino acid substitutions in proteins. Classically, researchers have focused on understanding the consequences of individual substitutions. However, literature surveys have shown that most substitutions were created at evolutionarily conserved positions. Awareness of this bias leads to a shift in perspective, from considering the outcomes of individual substitutions to understanding the roles of individual protein positions. Conserved positions tend to act as “toggle” switches, with most substitutions abolishing function. However, nonconserved positions have been found equally capable of affecting protein function. Indeed, many nonconserved positions act like functional dimmer switches (“rheostat” positions): this is revealed when multiple substitutions are made at a single position. Each substitution has a different functional outcome; the set of substitutions spans a range of outcomes. Finally, some nonconserved positions appear neutral, capable of accommodating all amino acid types without modifying function. This paper reviews the currently-known properties of rheostat positions, with examples shown for pyruvate kinase, organic anion transporting polypeptide 1B1, the beta-lactamase inhibitory protein, and angiotensin-converting enzyme 2. Outcomes observed for rheostat positions have implications for the rational design of drug analogs and allosteric drugs. Furthermore, this new framework—comprising three types of protein positions—provides a new approach to interpreting disease and population-based databases of amino acid changes. In conclusion, although a full understanding of substitution outcomes at rheostat positions poses a challenge, utilization of this new frame of reference will further advance the application of pharmacogenomics.

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.