Journal of Molecular Evolution最新文献

The voltage-sensing domain (VSD) is a module capable of responding to changes in the membrane potential through conformational changes and facilitating electromechanical coupling to open a pore gate, activate proton permeation pathways, or promote enzymatic activity in some membrane-anchored phosphatases. To carry out these functions, this module acts cooperatively through conformational changes. The VSD is formed by four transmembrane segments (S1-S4) but the S4 segment is critical since it carries positively charged residues, mainly Arg or Lys, which require an aqueous environment for its proper function. The discovery of this module in voltage-gated ion channels (VGICs), proton channels (Hv1), and voltage sensor-containing phosphatases (VSPs) has expanded our understanding of the principle of modularity in the voltage-sensing mechanism of these proteins. Here, by sequence comparison and the evaluation of the relationship between sequence composition, intrinsic flexibility, and structural analysis in 14 selected representatives of these three major protein groups, we report five interesting differences in the folding patterns of the VSD both in prokaryotes and eukaryotes. Our main findings indicate that this module is highly conserved throughout the evolutionary scale, however: (1) segments S1 to S3 in eukaryotes are significantly more hydrophobic than those present in prokaryotes; (2) the S4 segment has retained its hydrophilic character; (3) in eukaryotes the extramembranous linkers are significantly larger and more flexible in comparison with those present in prokaryotes; (4) the sensors present in the kHv1 proton channel and the ciVSP phosphatase, both of eukaryotic origin, exhibit relationships of flexibility and folding patterns very close to the typical ones found in prokaryotic voltage sensors; and (5) archaeal channels KvAP and MVP have flexibility profiles which are clearly contrasting in the S3-S4 region, which could explain their divergent activation mechanisms. Finally, to elucidate the obscure origins of this module, we show further evidence for a possible connection between voltage sensors and TolQ proteins.

For several decades, it has been known that a substantial number of genes within human DNA exhibit overlap; however, the biological and evolutionary significance of these overlaps remain poorly understood. This study focused on investigating specific instances of overlap where the overlapping DNA region encompasses the coding DNA sequences (CDSs) of protein-coding genes. The results revealed that proteins encoded by overlapping CDSs exhibit greater disorder than those from nonoverlapping CDSs. Additionally, these DNA regions were identified as GC-rich. This could be partially attributed to the absence of stop codons from two distinct reading frames rather than one. Furthermore, these regions were found to harbour fewer single-nucleotide polymorphism (SNP) sites, possibly due to constraints arising from the overlapping state where mutations could affect two genes simultaneously.While elucidating these properties, the NR1D1-THRA gene pair emerged as an exceptional case with highly structured proteins and a distinctly conserved sequence across eutherian mammals. Both NR1D1 and THRA are nuclear receptors lacking a ligand-binding domain at their C-terminus, which is the region where these gene pairs overlap. The NR1D1 gene is involved in the regulation of circadian rhythm, while the THRA gene encodes a thyroid hormone receptor, and both play crucial roles in various physiological processes. This study suggests that, in addition to their well-established functions, the specifically overlapping CDS regions of these genes may encode protein segments with additional, yet undiscovered, biological roles.

Folds are the architecture and topology of a protein domain. Categories of folds are very few compared to the astronomical number of sequences. Eukaryotes have more protein folds than Archaea and Bacteria. These folds are of two types: shared with Archaea and/or Bacteria on one hand and specific to eukaryotic clades on the other hand. The first kind of folds is inherited from the first endosymbiosis and confirms the mixed origin of eukaryotes. In a dataset of 1073 folds whose presence or absence has been evidenced among 210 species equally distributed in the three super-kingdoms, we have identified 28 eukaryotic folds unambiguously inherited from Bacteria and 40 eukaryotic folds unambiguously inherited from Archaea. Compared to previous studies, the repartition of informational function is higher than expected for folds originated from Bacteria and as high as expected for folds inherited from Archaea. The second type of folds is specifically eukaryotic and associated with an increase of new folds within eukaryotes distributed in particular clades. Reconstructed ancestral states coupled with dating of each node on the tree of life provided fold appearance rates. The rate is on average twice higher within Eukaryota than within Bacteria or Archaea. The highest rates are found in the origins of eukaryotes, holozoans, metazoans, metazoans stricto sensu, and vertebrates: the roots of these clades correspond to bursts of fold evolution. We could correlate the functions of some of the fold synapomorphies within eukaryotes with significant evolutionary events. Among them, we find evidence for the rise of multicellularity, adaptive immune system, or virus folds which could be linked to an ecological shift made by tetrapods.

Results from phylogenetic analyses that study the evolution of species according to their biological characteristics are frequently structured as phylogenetic trees. One of the most widely used methods for reconstructing them is the distance-based method known as the neighbor-joining (NJ) algorithm. It is known that the NJ algorithm can produce different phylogenetic trees depending on the order of the taxa in the input matrix of evolutionary distances, because the method only yields bifurcating branches or dichotomies. According to this, results and conclusions published in articles that only calculate one of the possible dichotomic phylogenetic trees are somehow biased. We have generalized the formulas used in the NJ algorithm to cope with Multifurcating branches or polytomies, and we have called this new variant of the method the multifurcating neighbor-joining (MFNJ) algorithm. Instead of the dichotomic phylogenetic trees reconstructed by the NJ algorithm, the MFNJ algorithm produces polytomic phylogenetic trees. The main advantage of using the MFNJ algorithm is that only one phylogenetic tree can be obtained, which makes the experimental section of any study completely reproducible and unbiased to external issues such as the input order of taxa.

Investigations of the molecular mechanisms behind detection of short, and particularly ultraviolet, wavelengths in arthropods have relied heavily on studies from insects due to the relative ease of heterologous expression of modified opsin proteins in model organisms like Drosophila. However, species outside of the Insecta can provide information on mechanisms for spectral tuning as well as the evolutionary history of pancrustacean visual pigments. Here we investigate the basis of spectral tuning in malacostracan short wavelength sensitive (SWS) opsins using phylogenetic comparative methods. Tuning sites that may be responsible for the difference between ultraviolet (UV) and violet visual pigment absorbance in the Malacostraca are identified, and the idea that an amino acid polymorphism at a single site is responsible for this shift is shown to be unlikely. Instead, we suggest that this change in absorbance is accomplished through multiple amino acid substitutions. On the basis of our findings, we conducted further surveys to identify spectral tuning mechanisms in the order Stomatopoda where duplication of UV opsins has occurred. Ancestral state reconstructions of stomatopod opsins from two main clades provide insight into the amino acid changes that lead to differing absorption by the visual pigments they form, and likely contribute the basis for the wide array of UV spectral sensitivities found in this order.

The genetic basis underlying adaptive physiological mechanisms has been extensively explored in mammals after colonizing the seas. However, independent lineages of aquatic mammals exhibit complex patterns of secondary colonization in freshwater environments. This change in habitat represents new osmotic challenges, and additional changes in key systems, such as the osmoregulatory system, are expected. Here, we studied the selective regime on coding and regulatory regions of 20 genes related to the osmoregulation system in strict aquatic mammals from independent evolutionary lineages, cetaceans, and sirenians, with representatives in marine and freshwater aquatic environments. We identified positive selection signals in genes encoding the protein vasopressin (AVP) in mammalian lineages with secondary colonization in the fluvial environment and in aquaporins for lineages inhabiting the marine and fluvial environments. A greater number of sites with positive selection signals were found for the dolphin species compared to the Amazonian manatee. Only the AQP5 and AVP genes showed selection signals in more than one independent lineage of these mammals. Furthermore, the vasopressin gene tree indicates greater similarity in river dolphin sequences despite the independence of their lineages based on the species tree. Patterns of distribution and enrichment of Transcription Factors in the promoter regions of target genes were analyzed and appear to be phylogenetically conserved among sister species. We found accelerated evolution signs in genes ACE, AQP1, AQP5, AQP7, AVP, NPP4, and NPR1 for the fluvial mammals. Together, these results allow a greater understanding of the molecular bases of the evolution of genes responsible for osmotic control in aquatic mammals.

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Venomous marine gastropods of the family Conidae are among the most diversified predators in marine realm-in large due to their complex venoms. Besides being a valuable source of bioactive neuropeptides conotoxins, cone-snails venoms are an excellent model for molecular evolution studies, addressing origin of key innovations. However, these studies are handicapped by scarce current knowledge on the tissues involved in venom production, as it is generally assumed the sole prerogative of the venom gland (VG). The role of other secretory glands that are present in all Conus species (salivary gland, SG) or only in some species (accessory salivary gland, ASG) remains poorly understood. Here, for the first time, we carry out a detailed analysis of the VG, SG, and ASG transcriptomes in the vermivorous Conus virgo. We detect multiple transcripts clusters in both the SG and ASG, whose annotations imply venom-related functions. Despite the subsets of transcripts highly-expressed in the VG, SG, and ASG being very distinct, SG expresses an L-, and ASG-Cerm08-, and MEFRR- superfamily conotoxins, all previously considered specific for VG. We corroborate our results with the analysis of published SG and VG transcriptomes from unrelated fish-hunting C. geographus, and C. striatus, possibly fish-hunting C. rolani, and worm-hunting Conus quercinus. In spite of low expression levels of conotoxins, some other specific clusters of putative venom-related peptides are present and may be highly expressed in the SG of these species. Further functional studies are necessary to determine the role that these peptides play in envenomation. In the meantime, our results show importance of routine multi-tissue sampling both for accurate interpretation of tissue-specific venom composition in cone-snails, and for better understanding origin and evolution of venom peptides genes.

Hypoxia-inducible factor 1 (HIF-1) is a crucial transcriptional factor that can restore oxygen balance in the body by regulating multiple vital activities. Two HIF-1α copies were retained in cyprinid fish after experiencing a teleost-specific genome duplication. How the "divergent collaboration" of HIF-1αA and HIF-1αB proceeds in regulating mitophagy and apoptosis under hypoxic stress in cells of cyprinid fish remains unclear. In this study, zebrafish HIF-1αA/B expression plasmids were constructed and transfected into the epithelioma papulosum cyprini cells and were subjected to hypoxic stress. HIF-1αA induced apoptosis through promoting ROS generation and mitochondrial depolarization when cells were subjected to oxygen deficiency. Conversely, HIF-1αB was primarily responsible for mitophagy induction, prompting ATP production to mitigate apoptosis. HIF-1αA did not induce mitophagy in the mitochondria and lysosomes co-localization assay but it was involved in the regulation of different mitophagy pathways. Over-expression of HIF-1αA increased the expression of bnip3, fundc1, Beclin1, and foxo3, suggesting it has a dual role in mitochondrial autophagy and cell death. Each duplicated copy also experienced functional divergence and target shifting in the regulation of complexes in the mitochondrial electron transport chain (ETC). Our findings shed light on the post-subfunctionalization function of HIF-1αA and HIF-1αB in zebrafish to fine-tune regulation of mitophagy and apoptosis following hypoxia exposure.

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject.