Search

Stochastic mapping is a simulation-based method for probabilistically mapping substitution histories onto phylogenies according to continuous-time Markov models of evolution. This technique can be used to infer properties of the evolutionary process on the phylogeny and, unlike parsimony-based mapping, conditions on the observed data to randomly draw substitution mappings that do not necessarily require the minimum number of events on a tree. Most stochastic mapping applications simulate substitution mappings only to estimate the mean and/or variance of two commonly used mapping summaries: the number of particular types of substitutions (labeled substitution counts) and the time spent in a particular group of states (labeled dwelling times) on the tree. Fast, simulation-free algorithms for calculating the mean of stochastic mapping summaries exist. Importantly, these algorithms scale linearly in the number of tips/leaves of the phylogenetic tree. However, to our knowledge, no such algorithm exists for calculating higher-order moments of stochastic mapping summaries. We present one such simulation-free dynamic programming algorithm that calculates prior and posterior mapping variances and scales linearly in the number of phylogeny tips. Our procedure suggests a general framework that can be used to efficiently compute higher-order moments of stochastic mapping summaries without simulations. We demonstrate the usefulness of our algorithm by extending previously developed statistical tests for rate variation across sites and for detecting evolutionarily conserved regions in genomic sequences.

MicroRNA166 (miR166) is highly conserved and has diverse functions across plant species. The highbush blueberry (Vaccinium corymbosum) genome is thought to harbor 10 miRNA166 loci (Vco-miR166), but the extent of their evolutionary conservation or functional diversification remains unknown. In this study, we identified six additional Vco-miR166 loci based on conserved features of the miR166 family. Phylogenetic analyses showed that mature Vco-miR166s and their precursor cluster in several clades are evolutionary conserved with diverse species. The cis -regulatory elements in the Vco-miR166 promoters indicated functions related to different phytohormones and defense responses. We also identified putative targets of vco-miR166s, which targeted the same gene families, suggesting the functional conservation and diversification of Vco-miR166 family members. Furthermore, we examined the accumulation patterns of six mature Vco-miR166s in response to abiotic stresses by stem-loop reverse RT-qPCR, which revealed their upregulation under freezing, cold, and heat stress, while they were downregulated by drought compared to control growth conditions. However, Vco-miR166 members showed different expression patterns when exposed to salt stress. These results showed that conserved Vco-miR166 family members display functional diversification but also coordinately influence plant responses to abiotic stress. [ABSTRACT FROM AUTHOR]

Intrinsically disordered proteins (IDPs) effect biological function despite their sequence‐encoded lack of preference for stable three‐dimensional structure. Among their many functions, IDPs form membraneless cellular compartments through liquid–liquid phase separation (LLPS), also termed biomolecular condensation. The extent to which LLPS has been evolutionarily selected remains largely unknown, as the complexities of IDP evolution hamper progress. Unlike structured proteins, rapid sequence divergence typical of IDPs confounds inference of their biophysical or biological functions from comparative sequence analyses. Here, we leverage mitosis as a universal eukaryotic feature to interrogate condensate evolutionary history. We observe that evolution has conserved the ability for six homologs of the mitotic IDP BuGZ to undergo LLPS and to serve the same mitotic function, despite low sequence conservation. We also observe that cellular context may tune LLPS. The phylogenetic correlation of LLPS and mitotic function in one protein raises the possibility of an ancient evolutionary interplay between LLPS and biological function, dating back at least 1.6 billion years to the last common ancestor of plants and animals. [ABSTRACT FROM AUTHOR]

While the concept of "evolutionary conservation" has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or "derivedness". We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, "conservation" mainly considers genes or traits which species have in common, while "derivedness" additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, "transcriptomic derivedness index", for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits. [ABSTRACT FROM AUTHOR]

MicroRNA166 (miR166) is highly conserved and has diverse functions across plant species. The highbush blueberry (Vaccinium corymbosum) genome is thought to harbor 10 miRNA166 loci (Vco-miR166), but the extent of their evolutionary conservation or functional diversification remains unknown. In this study, we identified six additional Vco-miR166 loci based on conserved features of the miR166 family. Phylogenetic analyses showed that mature Vco-miR166s and their precursor cluster in several clades are evolutionary conserved with diverse species. The cis-regulatory elements in the Vco-miR166 promoters indicated functions related to different phytohormones and defense responses. We also identified putative targets of vco-miR166s, which targeted the same gene families, suggesting the functional conservation and diversification of Vco-miR166 family members. Furthermore, we examined the accumulation patterns of six mature Vco-miR166s in response to abiotic stresses by stem-loop reverse RT-qPCR, which revealed their upregulation under freezing, cold, and heat stress, while they were downregulated by drought compared to control growth conditions. However, Vco-miR166 members showed different expression patterns when exposed to salt stress. These results showed that conserved Vco-miR166 family members display functional diversification but also coordinately influence plant responses to abiotic stress.

Mutations that reduce the biosynthetic cost of ATP production or increase the gene translation efficiency (tAI) are favorable for rapid cell growth and proliferation and therefore likely to be observed in tumors. Whether the mutations in tumors optimize the trade-off between the ATP biosynthesis cost and gene translation efficiency by increasing the tAI/ATP ratio is currently unknown. We retrieved transcriptome data of normal and osteosarcoma tissue samples from humans and mice and identified tumor-specific mutations in each species by using stringent cutoffs and outgroup information. We compared the tAI/ATP values of genes before and after mutation. The tAI/ATP profile was found to be highly conserved in humans and mice, and also correlated with the essentiality of genes. Tumor-specific rather than shared mutations were found to lead to increased tAI/ATP values in both species. Thus, tumor-specific mutations were found to optimize the cost-efficiency trade-off by increasing the tAI/ATP ratio of genes in osteosarcoma. This may indicate an evolutionarily conserved mechanism that promotes tumorigenesis by facilitating rapid cell growth and proliferation. [ABSTRACT FROM AUTHOR]

Cytolytic T cell responses are predicted to be biased towards membrane proteins. The peptide-binding grooves of most alleles of histocompatibility complex class I (MHC-I) are relatively hydrophobic, therefore peptide fragments derived from human transmembrane helices (TMHs) are predicted to be presented more often as would be expected based on their abundance in the proteome. However, the physiological reason of why membrane proteins might be over-presented is unclear. In this study, we show that the predicted over-presentation of TMH-derived peptides is general, as it is predicted for bacteria and viruses and for both MHC-I and MHC-II, and confirmed by re-analysis of epitope databases. Moreover, we show that TMHs are evolutionarily more conserved, because single nucleotide polymorphisms (SNPs) are present relatively less frequently in TMH-coding chromosomal regions compared to regions coding for extracellular and cytoplasmic protein regions. Thus, our findings suggest that both cytolytic and helper T cells are more tuned to respond to membrane proteins, because these are evolutionary more conserved. We speculate that TMHs are less prone to mutations that enable pathogens to evade T cell responses. [ABSTRACT FROM AUTHOR]

Flowering is the key process for the sexual reproduction in seed plants. In gramineous crops, the process of flowering, which includes the actions of both glume opening and glume closing, is directly driven by the swelling and withering of lodicules due to the water flow into and out of lodicule cells. All these processes are considered to be controlled by aquaporins, which are the essential transmembrane proteins that facilitate the transport of water and other small molecules across the biological membranes. In the present study, the evolution of aquaporins and their contribution to flowering process in plants were investigated via an integration of genome-wide analysis and gene expression profiling. Across the barley genome, we found that HvTIP1;1 , HvTIP1;2 , HvTIP2;3 , and HvPIP2;1 were the predominant aquaporin genes in lodicules and significantly upregulated in responding to glume opening and closing, suggesting the importance of them in the flowering process of barley. Likewise, the putative homologs of the above four aquaporin genes were also abundantly expressed in lodicules of the other monocots like rice and maize and in petals of eudicots like cotton, tobacco, and tomato. Furthermore, all of them were mostly upregulated in responding to the process of floret opening, indicating a conserved function of these aquaporin proteins in plant flowering. The phylogenetic analysis based on the OneKP database revealed that the homologs of TIP1;1, TIP1;2, TIP2;3, and PIP2;1 were highly conserved during the evolution, especially in the angiosperm species, in line with their conserved function in controlling the flowering process. Taken together, it could be concluded that the highly evolutionary conservation of TIP1;1, TIP1;2, TIP2;3 and PIP2;1 plays important roles in the flowering process for both monocots and eudicots. [ABSTRACT FROM AUTHOR]

The histone variant H2A.Z is deposited into chromatin by specific machinery and is required for genome functions. Using a linker-mediated complex strategy combined with yeast genetic complementation, we demonstrate evolutionary conservation of H2A.Z together with its chromatin incorporation and functions. This approach is applicable to the evolutionary analyses of proteins that form complexes with interactors. [ABSTRACT FROM AUTHOR]

Thioredoxins, small disulphide-containing redox proteins, play an important role in the regulation of cellular thiol redox balance through their disulfide reductase activity. In this study, we have identified, cloned, purified and characterized thioredoxin 1 (Hv Trx1) from the Cnidarian Hydra vulgaris Ind-Pune. Bioinformatics analysis revealed that Hv Trx1 contains an evolutionarily conserved catalytic active site Cys-Gly-Pro-Cys and shows a closer phylogenetic relationship with vertebrate Trx1. Optimum pH and temperature for enzyme activity of purified Hv Trx1 was found to be pH 7.0 and 25°C, respectively. Enzyme activity decreased significantly at acidic or alkaline pH as well as at higher temperatures. Hv Trx1 was found to be expressed ubiquitously in whole mount in situ hybridization. Treatment of Hydra with hydrogen peroxide (H2O2), a highly reactive oxidizing agent, led to a significant increase in gene expression and enzyme activity of Trx1. Further experiments using PX12, an inhibitor of Trx1, indicated that Trx1 plays an important role in regeneration in Hydra. Finally, by using growth assay in Escherichia coli and wound healing assay in human colon cancer cells, we demonstrate that Hv Trx1 is functionally active in both prokaryotic and eukaryotic heterologous systems. [ABSTRACT FROM AUTHOR]

Due to the redundancy of the genetic code most amino acids are encoded by several 'synonymous' codons. These codons are used unevenly, and each organism demonstrates its own unique codon usage bias, where the 'preferred' codons are associated with tRNAs that are found in high concentrations. Therefore, for decades, the prevailing view had been that preferred and non-preferred codons are linked to high or slow translation rates, respectively. However, this simplified view is contrasted by the frequent failures of codon-optimization efforts and by evidence of non-preferred (i.e. 'slow') codons having specific roles important for efficient production of functional proteins. One such specific role of slower codons is the regulation of co-translational protein folding, a complex biophysical process that is very challenging to model or to measure. Here, we combined a genome-wide approach with experiments to investigate the role of slow codons in protein production and co-translational folding. We analysed homologous gene groups from divergent bacteria and identified positions of inter-species conservation of bias towards slow codons. We then generated mutants where the conserved slow codons are substituted with 'fast' ones, and experimentally studied the effects of these codon substitutions. Using cellular and biochemical approaches we find that at certain locations, slow-to-fast codon substitutions reduce protein expression, increase protein aggregation, and impair protein function. This report provides an approach for identifying functionally relevant regions with slower codons and demonstrates that such codons are important for protein expression and function. [ABSTRACT FROM AUTHOR]

Abstract It is now well admitted by ecologists that the conservation of biodiversity should imply preserving the evolutionary processes that will permit its adaptation to ongoing and future environmental changes. This is attested by the ever‐growing reference to the conservation of evolutionary potential in the scientific literature. The impression that one may have when reading papers is that conserving evolutionary potential can only be a good thing, whatever biological system is under scrutiny. However, different objectives, such as maintaining species richness versus ecosystem services, may express different, when not conflicting, underlying values attributed to biodiversity. For instance, biodiversity can be intrinsically valued, as worth it to be conserved per se, or it can be conserved as a means for human flourishing. Consequently, both the concept of evolutionary potential and the prescriptions derived from the commitment to conserve it remain problematic, due to a lack of explicit mention of the norms underlying different conservation visions. Here, we contend that those who advocate for the conservation of evolutionary potential should position their conception along four dimensions: what vehicles instantiate the evolutionary potential relevant to their normative commitment; what temporality is involved; how measurable evolutionary potential is, and what degree of human influence is tolerated. We need to address these dimensions if we are to determine why and when the maintenance of evolutionary potential is an appropriate target for the conservation of biodiversity.

Abstract When populations of a rare species are small, isolated and declining under climate change, some populations may become locally maladapted. Detecting this maladaptation may allow effective rapid conservation interventions, even if based on incomplete knowledge. Population maladaptation may be estimated by finding genome–environment associations (GEA) between allele frequencies and environmental variables across a local species range, and identifying populations whose allele frequencies do not fit with these trends. We can then design assisted gene flow strategies for maladapted populations, to adjust their allele frequencies, entailing lower levels of intervention than with undirected conservation action. Here, we investigate this strategy in Scottish populations of the montane plant dwarf birch (Betula nana). In genome‐wide restriction site‐associated single nucleotide polymorphism (SNP) data, we found 267 significant associations between SNP loci and environmental variables. We ranked populations by maladaptation estimated using allele frequency deviation from the general trends at these loci; this gave a different prioritization for conservation action than the Shapely Index, which seeks to preserve rare neutral variation. Populations estimated to be maladapted in their allele frequencies at loci associated with annual mean temperature were found to have reduced catkin production. Using an environmental niche modelling (ENM) approach, we found annual mean temperature (35%), and mean diurnal range (15%), to be important predictors of the dwarf birch distribution. Intriguingly, there was a significant correlation between the number of loci associated with each environmental variable in the GEA and the importance of that variable in the ENM. Together, these results suggest that the same environmental variables determine both adaptive genetic variation and species range in Scottish dwarf birch. We suggest an assisted gene flow strategy that aims to maximize the local adaptation of dwarf birch populations under climate change by matching allele frequencies to current and future environments.

Flowering is the key process for the sexual reproduction in seed plants. In gramineous crops, the process of flowering, which includes the actions of both glume opening and glume closing, is directly driven by the swelling and withering of lodicules due to the water flow into and out of lodicule cells. All these processes are considered to be controlled by aquaporins, which are the essential transmembrane proteins that facilitate the transport of water and other small molecules across the biological membranes. In the present study, the evolution of aquaporins and their contribution to flowering process in plants were investigated via an integration of genome-wide analysis and gene expression profiling. Across the barley genome, we found that HvTIP1;1, HvTIP1;2, HvTIP2;3, and HvPIP2;1 were the predominant aquaporin genes in lodicules and significantly upregulated in responding to glume opening and closing, suggesting the importance of them in the flowering process of barley. Likewise, the putative homologs of the above four aquaporin genes were also abundantly expressed in lodicules of the other monocots like rice and maize and in petals of eudicots like cotton, tobacco, and tomato. Furthermore, all of them were mostly upregulated in responding to the process of floret opening, indicating a conserved function of these aquaporin proteins in plant flowering. The phylogenetic analysis based on the OneKP database revealed that the homologs of TIP1;1, TIP1;2, TIP2;3, and PIP2;1 were highly conserved during the evolution, especially in the angiosperm species, in line with their conserved function in controlling the flowering process. Taken together, it could be concluded that the highly evolutionary conservation of TIP1;1, TIP1;2, TIP2;3 and PIP2;1 plays important roles in the flowering process for both monocots and eudicots.

The 'Nest' motif plays a functional role in protein owing to its ligand binding potential aided by geometric concavity. The presence of less favored left-handed conformation (L-state) in its structure makes this concavity possible and in shaping the native chemical environment amenable to stable binding interactions. To understand the persistent appearance of L-state torsion in the Nest motif, we analyzed 0.5μs Molecular Dynamics (MD) simulation trajectories of 35 six-residue peptides (out of a total of 50 large Nest sequences of ≥6 residues) identified in our previous study. Analysis of the MD trajectories of the individual peptides reveals initial L-state in 60% of the peptides persists for >40% of the trajectory. Further, Nests with different sequences appear to adopt a specific conformational state driven by the neighboring L-state residues. The sequences also possess short secondary structures and amino acid repeats, suggesting evolutionary conservation and the specific role of amino acids in locally predisposing the torsion angle to the L-state. These findings help us to understand how L-state conformation is an essential prerequisite in stabilizing the Nest motif and shed light on the sequence-structure-function paradigm in the rational design of peptides and peptidomimetics for therapeutics.Communicated by Ramaswamy H. Sarma.

In response to environmental stress, human cells have been shown to form reversible amyloid aggregates within the nucleus, termed amyloid bodies (A-bodies). These protective physiological structures share many of the biophysical characteristics associated with the pathological amyloids found in Alzheimer's and Parkinson's disease. Here, we show that A-bodies are evolutionarily conserved across the eukaryotic domain, with their detection in Drosophila melanogaster and Saccharomyces cerevisiae marking the first examples of these functional amyloids being induced outside of a cultured cell setting. The conditions triggering amyloidogenesis varied significantly among the species tested, with results indicating that A-body formation is a severe, but sublethal, stress response pathway that is tailored to the environmental norms of an organism. RNA-sequencing analyses demonstrate that the regulatory low-complexity long non-coding RNAs that drive A-body aggregation are both conserved and essential in human, mouse and chicken cells. Thus, the identification of these natural and reversible functional amyloids in a variety of evolutionarily diverse species highlights the physiological significance of this protein conformation, and will be informative in advancing our understanding of both functional and pathological amyloid aggregation events. [ABSTRACT FROM AUTHOR]

The mechanisms controlling membrane recognition by proteins with one hydrophobic stretch at their carboxyl terminus (tail anchor, TA) are poorly defined. The Escherichia coli TAs of ElaB and YqjD, which share sequential and structural similarity with the Saccharomyces cerevisiae TA of Fis1, were shown to localize to mitochondria. We show that YqjD and ElaB are directed by their TAs to bacterial cell poles. Fis1(TA) expressed in E. coli localizes like the endogenous TAs. The yeast and bacterial TAs are inserted in the E. coli inner membrane, and they all show affiliation to phosphatidic acid (PA), found in the membrane of the bacterial cell poles and of the yeast mitochondria. Our results suggest a mechanism for TA membrane recognition conserved from bacteria to mitochondria and raise the possibility that through their interaction with PA, and TAs play a role across prokaryotes and eukaryotes in controlling cell/organelle fate. [ABSTRACT FROM AUTHOR]

Mycobacterium tuberculosis harbours nine toxin-antitoxin (TA) systems of the MazEF family. MazEF TA modules are of immense importance due to the perceived role of the MazF toxin in M. tuberculosis persistence and disease. The MazE antitoxin has a disordered C-terminal domain that binds the toxin, MazF and neutralizes its endoribonuclease activity. However, the structure of most MazEF TA complexes remains unsolved till date, obscuring structural and functional information about the antitoxins. We present a facile method to identify toxin binding residues on the disordered antitoxin. Charged residue scanning mutagenesis was used to screen a yeast surface displayed MazE6 antitoxin library against its purified cognate partner, the MazF6 toxin. Binding residues were deciphered by probing the relative reduction in binding to the ligand by flow cytometry. We have used this to identify putative antitoxin interface residues and local structure attained by the antitoxin upon interaction in the MazEF6 TA system and the same methodology is readily applicable to other intrinsically disordered protein regions. [ABSTRACT FROM AUTHOR]

Nascent ribosomal 60S subunits undergo the last maturation steps in the cytoplasm. The last one involves removing the anti-association factor eIF6 from the 60S ribosomal surface by the joint action of the Elongation Factor-like 1 (EFL1) GTPase and the SBDS protein. Herein, we studied the evolutionary relationship of the EFL1 and EF-2 protein families and the functional conservation within EFL1 orthologues. Phylogenetic analysis demonstrated that the EFL1 proteins are exclusive of eukaryotes and share an evolutionary origin with the EF-2 and EF-G protein families. EFL1 proteins originated by gene duplication from the EF-2 proteins and specialized in ribosome maturation while the latter retained their function in translation. Some organisms have more than one EFL1 protein resulting from alternative splicing, while others are encoded in different genes originated by gene duplication. However, the function of these alternative EFL1 proteins is still unknown. We performed GTPase activity and complementation assays to study the functional conservation of EFL1 homologs alone and together with their SBDS counterparts. None of the orthologues or cross-species combinations could replace the function of the corresponding yeast EFL1•SBDS binomial. The complementation of SBDS interspecies chimeras indicates that domain 2 is vital for its function together with EFL1 and the 60S subunit. The results suggest a functional species-specificity and possible co-evolution between EFL1, SBDS, and the 60S ribosomal subunit. These findings set the basis for further studies directed to understand the molecular evolution of these proteins and their impact on ribosome biogenesis and disease. [ABSTRACT FROM AUTHOR]

Enhancers are often studied as noncoding regulatory elements that modulate the precise spatiotemporal expression of genes in a highly tissue-specific manner. This paradigm has been challenged by recent evidence of individual enhancers acting in multiple tissues or developmental contexts. However, the frequency of these enhancers with high degrees of "pleiotropy" out of all putative enhancers is not well understood. Consequently, it is unclear how the variation of enhancer pleiotropy corresponds to the variation in expression breadth of target genes. Here, we use multi-tissue chromatin maps from diverse human tissues to investigate the enhancer–gene interaction architecture while accounting for 1) the distribution of enhancer pleiotropy, 2) the variations of regulatory links from enhancers to target genes, and 3) the expression breadth of target genes. We show that most enhancers are tissue-specific and that highly pleiotropy enhancers account for <1% of all putative regulatory sequences in the human genome. Notably, several genomic features are indicative of increasing enhancer pleiotropy, including longer sequence length, greater number of links to genes, increasing abundance and diversity of encoded transcription factor motifs, and stronger evolutionary conservation. Intriguingly, the number of enhancers per gene remains remarkably consistent for all genes (∼14). However, enhancer pleiotropy does not directly translate to the expression breadth of target genes. We further present a series of Gaussian Mixture Models to represent this organization architecture. Consequently, we demonstrate that a modest trend of more pleiotropic enhancers targeting more broadly expressed genes can generate the observed diversity of expression breadths in the human genome. [ABSTRACT FROM AUTHOR]

Major depressive disorder (MDD) affects around 15% of the population at some stage in their lifetime. It can be gravely disabling and it is associated with increased risk of suicide. Genetics play an important role; however, there are additional environmental contributions to the pathogenesis. A number of possible risk genes that increase liability for developing symptoms of MDD have been identified in genome-wide association studies (GWAS). The goal of this study was to characterize the MDD risk genes with respect to the degree of evolutionary conservation in simpler model organisms such as Caenorhabditis elegans and zebrafish, the phenotypes associated with variation in these genes and the extent of network connectivity. The MDD risk genes showed higher conservation in C. elegans and zebrafish than genome-to-genome comparisons. In addition, there were recurring themes among the phenotypes associated with variation of these risk genes in C. elegans. The phenotype analysis revealed enrichment for essential genes with pleiotropic effects. Moreover, the MDD risk genes participated in more interactions with each other than did randomly-selected genes from similar-sized gene sets. Syntenic blocks of risk genes with common functional activities were also identified. By characterizing evolutionarily-conserved counterparts to the MDD risk genes, we have gained new insights into pathogenetic processes relevant to the emergence of depressive symptoms in man.

Interleukin-3 (IL-3) is a cytokine belonging to the family of common β (βc) and is involved in various biological systems. Its activity is mediated by the interaction with its receptor (IL-3R), a heterodimer composed of two distinct subunits: IL-3Rα and βc. IL-3 and its receptor, especially IL-3Rα, play a crucial role in pathologies like inflammatory diseases and therefore are interesting therapeutic targets. Here, we have performed an analysis of these proteins and their interaction based on structural and evolutionary information. We highlighted that IL-3 and IL-3Rα structural architectures are conserved across evolution and shared with other proteins belonging to the same βc family interleukin-5 (IL-5) and granulocyte–macrophage colony-stimulating factor (GM-CSF). The IL-3Rα/IL-3 interaction is mediated by a large interface in which most residues are surprisingly not conserved during evolution and across family members. In spite of this high variability, we suggested small regions constituted by few residues conserved during the evolution in both proteins that could be important for the binding affinity. [ABSTRACT FROM AUTHOR]

Major depressive disorder (MDD) affects around 15% of the population at some stage in their lifetime. It can be gravely disabling and it is associated with increased risk of suicide. Genetics play an important role; however, there are additional environmental contributions to the pathogenesis. A number of possible risk genes that increase liability for developing symptoms of MDD have been identified in genome-wide association studies (GWAS). The goal of this study was to characterize the MDD risk genes with respect to the degree of evolutionary conservation in simpler model organisms such as Caenorhabditis elegans and zebrafish, the phenotypes associated with variation in these genes and the extent of network connectivity. The MDD risk genes showed higher conservation in C. elegans and zebrafish than genome-to-genome comparisons. In addition, there were recurring themes among the phenotypes associated with variation of these risk genes in C. elegans. The phenotype analysis revealed enrichment for essential genes with pleiotropic effects. Moreover, the MDD risk genes participated in more interactions with each other than did randomly-selected genes from similar-sized gene sets. Syntenic blocks of risk genes with common functional activities were also identified. By characterizing evolutionarily-conserved counterparts to the MDD risk genes, we have gained new insights into pathogenetic processes relevant to the emergence of depressive symptoms in man. [ABSTRACT FROM AUTHOR]

The opioid receptor (OPR) family comprises the mu-, delta-, and kappa-opioid, and nociceptin receptors that belong to the superfamily of 7-transmembrane spanning G protein-coupled receptors (GPCRs). The mu-opioid receptor is the main target for clinically used opioid analgesics, and its biology has been extensively studied. The N-terminally truncated 6TM receptors isoform produced through alternative splicing of the OPRM1 gene displays unique signaling and analgesic properties, but it is unclear if other OPRs have the same ability. In this study, we have built a comprehensive map of alternative splicing events that produce 6TM receptor variants in all the OPRs and demonstrated their evolutionary conservation. We then obtained evidence for their translation through ribosomal footprint analysis. We discovered that N-terminally truncated 6TM GPCRs are rare in the human genome and OPRs are overrepresented in this group. Finally, we also observed a significant enrichment of 6TM GPCR genes among genes associated with pain, psychiatric disorders, and addiction. Understanding the biology of 6TM receptors and leveraging this knowledge for drug development should pave the way for novel therapies. [ABSTRACT FROM AUTHOR]

Background: Transcription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. It is critical to locate these TF-DNA interactions to understand transcriptional regulation. Efforts to predict bona fide TFBSs benefit from the availability of experimental data mapping DNA binding regions of TFs (chromatin immunoprecipitation followed by sequencing - ChIP-seq). Results: In this study, we processed ~ 10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~ 56 million TFBSs with experimental and computational support for direct TF-DNA interactions for 644 TFs in > 1000 cell lines and tissues. These TFBSs were used to predict > 197,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface (https://unibind.uio.no/), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions. Conclusions: UniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species. [ABSTRACT FROM AUTHOR]

Objectives: Bipolar disorder (BPD) is a highly heritable psychiatric disorder whose genetic complexity and pathogenetic mechanisms are still being unraveled. The main goal of this work was to characterize BPD risk‐gene candidates (identified by Nurnberger et al., JAMA Psychiatry 71:657, 2014, and Stahl et al., Nat. Genet. 51:793, 2019) with respect to their evolutionary conservation, associated phenotypes, and extent of gene‐gene interactions. Methods: Database searches and BLAST were used to identify homologous counterparts of human BPD risk genes in C. elegans, zebrafish, and Drosophila. Phenotypes associated with the C. elegans genes were annotated and searched. With GeneMANIA, we characterized and quantified gene‐gene interactions among members of the BPD gene set in comparison to randomly chosen gene sets of the same size. Results: BPD risk genes are highly conserved across species and are enriched for essential genes and genes associated with lethality and altered life span. They are significantly more interactive with each other in comparison to random genes. We identified syntenic blocks of risk genes, which provided potential insights into molecular pathways and co‐morbidities associated with BPD including coronary disease, obesity, and decreased life expectancy. Conclusions: BPD risk genes appear to be special in terms of their degree of conservation, interconnectedness, and pleiotropic effects that extend beyond a role in brain function. Key hub genes or pleiotropic regulatory components may represent attractive targets for future drug discovery. [ABSTRACT FROM AUTHOR]

The short-chain dehydrogenase/reductase (SDR) superfamily encompasses enzymes that play essential roles in the metabolism of steroid hormones and lipids. Despite an enigmatic function, recent genetic studies have linked the novel SDR 42 extended-1 (SDR42E1) gene to 25-hydroxyvitamin D levels. This study investigated the potential SDR42E1 functions and interactions with vitamin D using bioinformatics and molecular docking studies. Phylogenetic analysis unveiled that the nucleotide sequences of human SDR42E1 exhibit high evolutionary conservation across nematodes and fruit flies. Molecular docking analysis identified strong binding affinities between SDR42E1 and its orthologs with vitamin D 3 and essential precursors, 8-dehydrocholesterol, followed by 7-dehydrocholesterol and 25-hydroxyvitamin D. The hydrophobic interactions observed between the protein residues and vitamin D compounds supported the predicted transmembrane localization of SDR42E1. Our investigation provides valuable insights into the potential role of SDR42E1 in skin vitamin D biosynthesis throughout species. This provides the foundation for future research and development of targeted therapies for vitamin D deficiency and related health conditions. • SDR42E1 protein demonstrates a strong affinity for key vitamin D intermediates, notably 8-dehydrocholesterol. • The evolutionary conservation of hydrophobic interactions between vitamin D and SDR42E1 emphasizes its crucial role in regulating vitamin D skin synthesis. [ABSTRACT FROM AUTHOR]

Abstract Background Aquaporins (AQPs) are integral membrane proteins from a larger family of major intrinsic proteins (MIPs) and function in a huge variety of processes such as water transport, plant growth and stress response. The availability of the whole-genome data of different cotton species allows us to study systematic evolution and function of cotton AQPs on a genome-wide level. Results Here, a total of 53, 58, 113 and 111 AQP genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. A comprehensive analysis of cotton AQPs, involved in exon/intron structure, functional domains, phylogenetic relationships and gene duplications, divided these AQPs into five subfamilies (PIP, NIP, SIP, TIP and XIP). Comparative genome analysis among 30 species from algae to angiosperm as well as common tandem duplication events in 24 well-studied plants further revealed the evolutionary conservation of AQP family in the organism kingdom. Combining transcriptome analysis and Quantitative Real-time PCR (qRT-PCR) verification, most AQPs exhibited tissue-specific expression patterns both in G. raimondii and G. hirsutum. Meanwhile, a bias of time to peak expression of several AQPs was also detected after treating G. davidsonii and G. hirsutum with 200 mM NaCl. It is interesting that both PIP1;4 h/i/j and PIP2;2a/e showed the highly conserved tandem structure, but differentially contributed to tissue development and stress response in different cotton species. Conclusions These results demonstrated that cotton AQPs were structural conservation while experienced the functional differentiation during the process of evolution and domestication. This study will further broaden our insights into the evolution and functional elucidation of AQP gene family in cotton.

Abstract Background Structural homology modeling supported by bioinformatics analysis plays a key role in uncovering new molecular interactions within gene regulatory networks. Here, we have applied this powerful approach to analyze the molecular interactions orchestrating death receptor signaling networks. In particular, we focused on the molecular mechanisms of CD95-mediated NF-κB activation and the role of c-FLIP/NEMO interaction in the induction of this pathway. Results To this end, we have created the homology model of the c-FLIP/NEMO complex using the reported structure of the v-FLIP/NEMO complex, and rationally designed peptides targeting this complex. The designed peptides were based on the NEMO structure. Strikingly, the experimental in vitro validation demonstrated that the best inhibitory effects on CD95-mediated NF-κB activation are exhibited by the NEMO-derived peptides with the substitution D242Y of NEMO. Furthermore, we have assumed that the c-FLIP/NEMO complex is recruited to the DED filaments formed upon CD95 activation and validated this assumption in silico. Further insight into the function of c-FLIP/NEMO complex was provided by the analysis of evolutionary conservation of interacting regions which demonstrated that this interaction is common in distinct mammalian species. Conclusions Taken together, using a combination of bioinformatics and experimental approaches we obtained new insights into CD95-mediated NF-κB activation, providing manifold possibilities for targeting the death receptor network.

Membrane proteins are unique in that they interact with lipid bilayers, making them indispensable for transporting molecules and relaying signals between and across cells. Due to the significance of the protein's functions, mutations often have profound effects on the fitness of the host. This is apparent both from experimental studies, which implicated numerous missense variants in diseases, as well as from evolutionary signals that allow elucidating the physicochemical constraints that intermembrane and aqueous environments bring. In this review, we report on the current state of knowledge acquired on missense variants (referred to as to single amino acid variants) affecting membrane proteins as well as the insights that can be extrapolated from data already available. This includes an overview of the annotations for membrane protein variants that have been collated within databases dedicated to the topic, bioinformatics approaches that leverage evolutionary information in order to shed light on previously uncharacterized membrane protein structures or interaction interfaces, tools for predicting the effects of mutations tailored specifically towards the characteristics of membrane proteins as well as two clinically relevant case studies explaining the implications of mutated membrane proteins in cancer and cardiomyopathy. [ABSTRACT FROM AUTHOR]

The evolution of the mammalian cerebral cortex leading to humans involved a remarkable sophistication of developmental mechanisms. Specific adaptations of progenitor cell proliferation and neuronal migration mechanisms have been proposed to play major roles in this evolution of neocortical development. One of the central elements influencing neocortex development is the extracellular matrix (ECM). The ECM provides both a structural framework during tissue formation and to present signaling molecules to cells, which directly influences cell behavior and movement. Here we review recent advances in the understanding of the role of ECM molecules on progenitor cell proliferation and neuronal migration, and how these contribute to cerebral cortex expansion and folding. We discuss how transcriptomic studies in human, ferret and mouse identify components of ECM as being candidate key players in cortex expansion during development and evolution. Then we focus on recent functional studies showing that ECM components regulate cortical progenitor cell proliferation, neuron migration and the mechanical properties of the developing cortex. Finally, we discuss how these features differ between lissencephalic and gyrencephalic species, and how the molecular evolution of ECM components and their expression profiles may have been fundamental in the emergence and evolution of cortex folding across mammalian phylogeny.

While the concept of “evolutionary conservation” has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or “derivedness”. We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, “conservation” mainly considers genes or traits which species have in common, while “derivedness” additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, “transcriptomic derivedness index”, for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits.

Wenxiang diagrams illustrate protein helices as spirals on a plane and thus have the advantage over helical wheels of being planar graphs. Wenxiang 3.0 extends the original version by adding 3 major features: (1) individual amino acid residues can be colored according to their evolutionary conservation in comparative multiple sequence alignments using CONSURF encoding; (2) α, π, and 3/10 helices can be illustrated by overlaying arcs representative of the pitches of these helices; and, (3) the physico-chemical properties of amino acids residues in the protein sequence can be re-presented by colored geometric shapes. [ABSTRACT FROM AUTHOR]

During mate choice, three-spined stickleback (Gasterosteus aculeatus) females make use of male olfactory cues to determine the suitability of potential mates in a species- and habitat-dependent manner. This signal contains peptides reflecting the individual male's MHC profile and a so-called male validation factor (MVF) indicating the species identity of the sender. In a two-armed flow channel, we exposed gravid females to male odour. The odours from males that produced only the MVF but no MHC signal yet were used to examine whether the MVF of distant stickleback populations, one from Canada and one from Germany – genetically isolated since approximately 31–59 thousand years – has changed to such a degree that present-day females only react to the MVF produced by males from their population of origin. The fish used in the experiments were bred from in vitro-fertilized eggs in either country; those from Canada were shipped for analysis to the Plön laboratory to minimize experimental variability. Given the choice between the two MVFs, females did not significantly prefer the MVF of males of their own origin over the other. However, compared to water, both types of females preferred their own MVF significantly and the foreign MVF almost significantly, respectively, suggesting that they could smell MVFs from both origins. These results suggest that, despite long-term genetic isolation between Canadian and German populations, the three-spined stickleback male validation factor has remained constant and thus represents an evolutionarily conserved signal. [ABSTRACT FROM AUTHOR]

Aging is accompanied by the gradual deregulation of the transcriptome. However, whether age-dependent changes in the transcriptome are evolutionarily conserved or diverged remains largely unexplored. Here, we performed a meta-analysis examining the age-dependent changes in the transcriptome using publicly available datasets of 11 representative metazoans, ranging from Caenorhabditis elegans to humans. To identify the transcriptomic changes associated with aging, we analyzed various aspects of the transcriptome, including genome composition, RNA processing, and functional consequences. The use of introns and novel splice sites tended to increase with age, particularly in the brain. In addition, our analysis suggests that the age-dependent accumulation of premature termination codon-containing transcripts is a common feature of aging across multiple animal species. Using C. elegans as a test model, we showed that several splicing factors that are evolutionarily conserved and age-dependently downregulated were required to maintain a normal lifespan. Thus, aberrant RNA processing appears to be associated with aging and a short lifespan in various species., Competing Interests: Declaration of Competing Interests The authors have no potential conflicts of interest to disclose., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The hypothalamus is characterized by great neuronal diversity, with many neuropeptides and other neuromodulators being expressed within its multiple anatomical domains. The regulatory networks directing hypothalamic development have been studied in detail, but, for many neuron types, control of differentiation is still not understood. The highly conserved Brain-specific homeobox (Bsx) transcription factor has previously been described in regulating Agrp and Npy expression in the hypothalamic arcuate nucleus (ARC) in mice. While Bsx is expressed in many more subregions of both tuberal and mamillary hypothalamus, the functions therein are not known. Using genetic analyses in zebrafish, we show that most bsx expression domains are dependent on Nkx2.1 and Nkx2.4 homeodomain transcription factors, while a subset depends on Otp. We show that the anatomical pattern of the ventral forebrain appears normal in bsx mutants, but that Bsx is necessary for the expression of many neuropeptide encoding genes, including agrp, penka, vip, trh, npb, and nts, in distinct hypothalamic anatomical domains. We also found Bsx to be critical for normal expression of two Crh family members, crhb and uts1, as well as crhbp, in the hypothalamus and the telencephalic septal region. Furthermore, we demonstrate a crucial role for Bsx in serotonergic, histaminergic and nitrergic neuron development in the hypothalamus. We conclude that Bsx is critical for the terminal differentiation of multiple neuromodulatory cell types in the forebrain.

Knowledge of the detailed structure of macromolecular interactions is key to a better understanding and modulation of essential cellular functions and pathological situations. Great efforts are invested in the development of improved computational prediction methods, including binding site prediction and protein-protein docking. These tools should benefit from the inclusion of evolutionary information, since the pressure to maintain functional interactions leads to conservation signals on protein surfaces at interacting sites and coevolution between contacting positions. However, unveiling such constraints and finding the best way to integrate them into computational pipelines remains a challenging area of research. Here, we first introduce evolutionary properties of interface structures, focusing on recent work exploring evolutionary mechanisms at play in protein interfaces and characterizing the complexity of evolutionary signals, such as interface deep mutational scans. Then, we review binding site predictors and interface structure modeling methods that integrate conservation and coevolution as important ingredients to improve predictive capacity, ending with a section dedicated to the prediction of binding modes between a globular protein domain and a short motif located within an intrinsically disordered or flexible region. Throughout the review, we discuss practical applications and recent promising developments, in particular regarding the use of coevolutionary information for interface structural prediction and the integration of these coevolution signals with machine learning and deep learning algorithms. [ABSTRACT FROM AUTHOR]

• Zebrafish is becoming increasingly popular in behavioral neuroscience. • This review explores theoretical and practical reasons for it. • It emphasizes the need for the comparative approach. • It discusses questions about evolution and translational relevance. • It focusses on relational learning and memory. The zebrafish is becoming increasingly well utilized in several fields of biology, including behavioral neuroscience. This review, based upon an Outstanding Achievement Award lecture presented by the author at an IBNS conference, explores the potential reasons for the popularity of this species. First, some theoretical questions are discussed, including why a comparative approach using more than a single laboratory animal species (the mouse) may be important. Points about evolutionary conservation are considered. Discussion on why complexity may build over time, and what it means when a species is considered "primitive" vs "advanced", are also included. Arguments about how using a "primitive" species may work to the researcher's advantage are made in the context of translational relevance. Most discussions center around a complex functional aspect of the brain, relational learning and memory, and use empirical examples from the author's own laboratory. It is hoped that these examples will show how a simple vertebrate, like the zebrafish, may help the investigator address fundamental mechanistic questions about complex brain function and behavior. [ABSTRACT FROM AUTHOR]

It is now well admitted by ecologists that the conservation of biodiversity should imply preserving the evolutionary processes that will permit its adaptation to ongoing and future environmental changes. This is attested by the ever‐growing reference to the conservation of evolutionary potential in the scientific literature. The impression that one may have when reading papers is that conserving evolutionary potential can only be a good thing, whatever biological system is under scrutiny. However, different objectives, such as maintaining species richness versus ecosystem services, may express different, when not conflicting, underlying values attributed to biodiversity. For instance, biodiversity can be intrinsically valued, as worth it to be conserved per se, or it can be conserved as a means for human flourishing. Consequently, both the concept of evolutionary potential and the prescriptions derived from the commitment to conserve it remain problematic, due to a lack of explicit mention of the norms underlying different conservation visions. Here, we contend that those who advocate for the conservation of evolutionary potential should position their conception along four dimensions: what vehicles instantiate the evolutionary potential relevant to their normative commitment; what temporality is involved; how measurable evolutionary potential is, and what degree of human influence is tolerated. We need to address these dimensions if we are to determine why and when the maintenance of evolutionary potential is an appropriate target for the conservation of biodiversity. [ABSTRACT FROM AUTHOR]

The hypothalamus is characterized by great neuronal diversity, with many neuropeptides and other neuromodulators being expressed within its multiple anatomical domains. The regulatory networks directing hypothalamic development have been studied in detail, but, for many neuron types, control of differentiation is still not understood. The highly conserved Brain-specific homeobox (Bsx) transcription factor has previously been described in regulating Agrp and Npy expression in the hypothalamic arcuate nucleus (ARC) in mice. While Bsx is expressed in many more subregions of both tuberal and mamillary hypothalamus, the functions therein are not known. Using genetic analyses in zebrafish, we show that most bsx expression domains are dependent on Nkx2.1 and Nkx2.4 homeodomain transcription factors, while a subset depends on Otp. We show that the anatomical pattern of the ventral forebrain appears normal in bsx mutants, but that Bsx is necessary for the expression of many neuropeptide encoding genes, including agrp, penka , vip , trh , npb , and nts , in distinct hypothalamic anatomical domains. We also found Bsx to be critical for normal expression of two Crh family members, crhb and uts1 , as well as crhbp , in the hypothalamus and the telencephalic septal region. Furthermore, we demonstrate a crucial role for Bsx in serotonergic, histaminergic and nitrergic neuron development in the hypothalamus. We conclude that Bsx is critical for the terminal differentiation of multiple neuromodulatory cell types in the forebrain. [ABSTRACT FROM AUTHOR]

When populations of a rare species are small, isolated and declining under climate change, some populations may become locally maladapted. Detecting this maladaptation may allow effective rapid conservation interventions, even if based on incomplete knowledge. Population maladaptation may be estimated by finding genome–environment associations (GEA) between allele frequencies and environmental variables across a local species range, and identifying populations whose allele frequencies do not fit with these trends. We can then design assisted gene flow strategies for maladapted populations, to adjust their allele frequencies, entailing lower levels of intervention than with undirected conservation action. Here, we investigate this strategy in Scottish populations of the montane plant dwarf birch (Betula nana). In genome‐wide restriction site‐associated single nucleotide polymorphism (SNP) data, we found 267 significant associations between SNP loci and environmental variables. We ranked populations by maladaptation estimated using allele frequency deviation from the general trends at these loci; this gave a different prioritization for conservation action than the Shapely Index, which seeks to preserve rare neutral variation. Populations estimated to be maladapted in their allele frequencies at loci associated with annual mean temperature were found to have reduced catkin production. Using an environmental niche modelling (ENM) approach, we found annual mean temperature (35%), and mean diurnal range (15%), to be important predictors of the dwarf birch distribution. Intriguingly, there was a significant correlation between the number of loci associated with each environmental variable in the GEA and the importance of that variable in the ENM. Together, these results suggest that the same environmental variables determine both adaptive genetic variation and species range in Scottish dwarf birch. We suggest an assisted gene flow strategy that aims to maximize the local adaptation of dwarf birch populations under climate change by matching allele frequencies to current and future environments. [ABSTRACT FROM AUTHOR]

Patterns observed by examining the evolutionary relationships among proteins of common origin can reveal the structural and functional importance of specific residue positions. In particular, amino acids that are highly conserved (i.e., their positions evolve at a slower rate than other positions) are particularly likely to be of biological importance, for example, for ligand binding. ConSurf is a bioinformatics tool for accurately estimating the evolutionary rate of each position in a protein family. Here we introduce a new release of ConSurf‐DB, a database of precalculated ConSurf evolutionary conservation profiles for proteins of known structure. ConSurf‐DB provides high‐accuracy estimates of the evolutionary rates of the amino acids in each protein. A reliable estimate of a query protein's evolutionary rates depends on having a sufficiently large number of effective homologues (i.e., nonredundant yet sufficiently similar). With current sequence data, ConSurf‐DB covers 82% of the PDB proteins. It will be updated on a regular basis to ensure that coverage remains high—and that it might even increase. Much effort was dedicated to improving the user experience. The repository is available at https://consurfdb.tau.ac.il/. Broader audience: By comparing a protein to other proteins of similar origin, it is possible to determine the extent to which each amino acid position in the protein evolved slowly or rapidly. A protein's evolutionary profile can provide valuable insights: For example, amino acid positions that are highly conserved (i.e., evolved slowly) are particularly likely to be of structural and/or functional importance, for example, for ligand binding and catalysis. We introduce here a new and improved version of ConSurf‐DB, a continually updated database that provides precalculated evolutionary profiles of proteins with known structure. [ABSTRACT FROM AUTHOR]

Transcription factors bind sequence-specific sites in DNA to regulate gene transcription. Identifying transcription factor binding sites (TFBS’s) is an important step for understanding gene regulation. Although sophisticated in modeling TFBS’s and their combinatorial patterns, computational methods for TFBS detection and motif finding often make oversimplified homogeneous model assumptions for background sequences. Since nucleotide base composition varies across genomic regions, it is expected to be helpful for motif finding to incorporate the heterogeneity into background modeling. When sequences from multiple species are utilized, variation in evolutionary conservation violates the common assumption of an identical conservation level in multiple alignments. To handle both types of heterogeneity, we propose a generative model in which a segmented Markov chain is used to partition a multiple alignment into regions of homogeneous nucleotide base composition and a hidden Markov model (HMM) is employed to account for different conservation levels. Bayesian inference on the model is developed via Gibbs sampling with dynamic programming recursions. Simulation studies and empirical evidence from biological data sets reveal the dramatic effect of background modeling on motif finding, and demonstrate that the proposed approach is able to achieve substantial improvements over commonly used background models.

The myodural bridge (MDB) is confirmed that connecting the most of suboccipital muscles to the cervical dura mater through the posterior intervertebral spaces and widely exists in mammals and birds. In order to reveal whether the MDB is universally existing in amniota of vertebrates, we explored the existence and the morphological features of the MDB in the Trachemys scripta elegans. Twenty fresh red-eared slider specimens were observed by the gross anatomy dissection and histological analysis. In the results, three kind of muscles in the postoccipital region of the red-eared slider were found. The rectus capitis dorsum minor muscle originated from the posterior margin of the occiput (C0) and terminated at the spinous process of the atlas (C1). The transversospinales muscle was attached to the vertebral arch and the postzygapophysis of the atlas and extended to the spinous process of the axis (C2). The C2-C3 intertransversales muscle were extended from the postzygapophysis of C2 and the one of C3. The three muscles covered the dorsal interspaces among C0-C3, and meantime they were closely connected with dense connective tissues, which filled in these interspaces. Each of these thick dense connective tissue membranes sent off several short and strong fibrous bundles ventrally to merge with the cervical spinal dura mater. Furthermore the connective tissues connecting these muscles with cervical spinal dura mater directly were revealed under the microscopy and they consisted of parallel and intensive collagen fibers with orientation from dorsal to ventral. In conclusion, this study for the first time demonstrated the existence of the MDB in the testudines, in all of the dorsal atlantooccipital, atlantoaxial and C2-C3 intervertebral spaces. Based on our results and comparative anatomical evidences in recent year, it could be inferred that the MDB might be its highly conserved structure in the evolution of amniota. [ABSTRACT FROM AUTHOR]

Trypanosoma brucei, the agent of human sleeping sickness and ruminant nagana, is the most genetically tractable representative of the domain Excavata. It is evolutionarily very distant from humans, with a last common ancestor over 1 billion years ago. Frataxin, a highly conserved small protein involved in iron-sulfur cluster synthesis, is present in both organisms, and its deficiency is responsible for Friedreich's ataxia in humans. We have found that T. brucei growth-inhibition phenotype caused by down-regulated frataxin is rescued by means of human frataxin. The rescue is fully dependent on the human frataxin being imported into the trypanosome mitochondrion. Processing of the imported protein by mitochondrial processing peptidase can be blocked by mutations in the signal peptide, as in human cells. Although in human cells frataxin must be processed to execute its function, the same protein in the T. brucei mitochondrion is functional even in the absence of processing. Our results illuminate remarkable conservation of the mechanisms of mitochondrial protein import and processing.

Bone Morphogenetic Proteins (Bmps) are secreted growth factors that play crucial roles in animal development across the phylogenetic spectrum. Bmp signaling results in the phosphorylation and nuclear translocation of Smads, downstream signal transducers that bind DNA. In Drosophila, the zinc finger protein Schnurri (Shn) plays a key role in signaling by the Bmp2/Bmp4 homolog Decapentaplegic (Dpp), by forming a Shn/Smad repression complex on defined promoter elements in the brinker (brk) gene. Brk is a transcriptional repressor that downregulates Dpp target genes. Thus, brk inhibition by Shn results in the upregulation of Dpp-responsive genes. We present evidence that vertebrate Shn homologs can also mediate Bmp responsiveness through a mechanism similar to Drosophila Shn. We find that a Bmp response element (BRE) from the Xenopus Vent2 promoter drives Dpp-dependent expression in Drosophila. However, in sharp contrast to its activating role in vertebrates, the frog BRE mediates repression in Drosophila. Remarkably, despite these opposite transcriptional polarities, sequence changes that abolish cis-element activity in Drosophila also affect BRE function in Xenopus. These similar cis requirements reflect conservation of trans-acting factors, as human Shn1 (hShn1; HIVEP1) can interact with Smad1/Smad4 and assemble an hShn1/Smad complex on the BRE. Furthermore, both Shn and hShn1 activate the BRE in Xenopus embryos, and both repress brk and rescue embryonic patterning defects in shn mutants. Our results suggest that vertebrate Shn proteins function in Bmp signal transduction, and that Shn proteins recruit coactivators and co-repressors in a context-dependent manner, rather than acting as dedicated activators or repressors.