Distinct functional constraints partition sequence conservation in a cis-regulatory element

Different functional constraints contribute to different evolutionary rates across genomes. To understand why some sequences evolve faster than others in a single cis-regulatory locus, we investigated function and evolutionary dynamics of the promoter of the Caenorhabditis elegans unc-47 gene. We found that this promoter consists of two distinct domains. The proximal promoter is conserved and is largely sufficient to direct appropriate spatial expression. The distal promoter displays little if any conservation between several closely related nematodes. Despite this divergence, sequences from all species confer robustness of expression, arguing that this function does not require substantial sequence conservation. We showed that even unrelated sequences have the ability to promote robust expression. A prominent feature shared by all of these robustness-promoting sequences is an AT-enriched nucleotide composition consistent with nucleosome depletion. Because general sequence composition can be maintained despite sequence turnover, our results explain how different functional constraints can lead to vastly disparate rates of sequence divergence within a promoter.     

Identification of an orthologous mammalian cytokeratin gene. High degree of intron sequence conservation during evolution of human cytokeratin 10

Among the human acidic (type I) cytokeratins, components 10 and 11 are especially interesting, as they are under various kinds of expression control. They are synthesized in the suprabasal cell layers of certain stratified epithelia, notably epidermis, in an endogenous differentiation program; they are expressed in certain epithelial tumours but not in others; they can appear de novo in certain pathological situations such as in squamous metaplasias; and their expression in vivo and in vitro is under positive influence of extracellular calcium concentrations and is reduced in the presence of vitamin A or other retinoids. To provide a basis for studies of the various regulatory elements, we have isolated the human gene encoding cytokeratin 10, using a cDNA probe derived from the corresponding bovine gene, and have sequenced the mRNA coding region as well as adjacent regions approximately 1500 bases 5' upstream and 1000 bases 3' downstream. The eight exons encode a polypeptide 59,535 Mr, i.e. somewhat larger than the corresponding bovine and murine proteins. The deduced amino acid sequences display a high degree of homology, which is not restricted to the exons and the 5' and 3' adjacent regions but, surprisingly, is also evident in the seven introns, some of which contain extended sequence elements with 70% identical nucleotides and more, i.e. similar to the homology in the adjacent exons. This exceptionally high level of conservation of intron sequences is discussed in relation to the recently accumulating evidence of the occurrence of intron sequences important in the regulation of the expression of members of other multigene families during development.     

Remarkable sequence conservation of the last intron in the PKD1 gene

The last intron of the PKD1 gene (intron 45) was found to have exceptionally high sequence conservation across four mammalian species: human, mouse, rat, and dog. This conservation did not extend to the comparable intron in pufferfish. Pairwise comparisons for intron 45 showed 91% identity (human vs. dog) to 100% identity (mouse vs. rat) for an average for all four species of 94% identity. In contrast, introns 43 and 44 of the PKD1 gene had average pairwise identities of 57% and 54%, and exons 43, 44, and 45 and the coding region of exon 46 had average pairwise identities of 80%, 84%, 82%, and 80%. Intron 45 is 90 to 95 bp in length, with the major region of sequence divergence being in a central 4-bp to 9-bp variable region. RNA secondary structure analysis of intron 45 predicts a branching stem-loop structure in which the central variable region lies in one loop and the putative branch point sequence lies in another loop, suggesting that the intron adopts a specific stem-loop structure that may be important for its removal. Although intron 45 appears to conform to the class of small, G-triplet-containing introns that are spliced by a mechanism utilizing intron definition, its high sequence conservation may be a reflection of constraints imposed by a unique mechanism that coordinates splicing of this last PKD1 intron with polyadenylation.     

engrailed sequences from four centipede orders: strong sequence conservation, duplications and phylogeny

We cloned and sequenced parts of the homologues of the engrailed gene from eight species of centipedes (Chilopoda), an arthropod group with very diverse patterns of body segmentation. We found very high sequence conservation and two independent instances of gene duplication (in Lithobius forficatus and Geophilus carpophagus). Gene phylogeny based on available engrailed sequences agrees with expected topology, but for two minor exceptions, both probably due to long branch effects.     

Comparative phosphoproteomics reveals evolutionary and functional conservation of phosphorylation across eukaryotes

Background  

Reversible phosphorylation of proteins is involved in a wide range of processes, ranging from signaling cascades to regulation of protein complex assembly. Little is known about the structure and evolution of phosphorylation networks. Recent high-throughput phosphoproteomics studies have resulted in the rapid accumulation of phosphopeptide datasets for many model organisms. Here, we exploit these novel data for the comparative analysis of phosphorylation events between different species of eukaryotes.     

Human dihydrofolate reductase gene organization. Extensive conservation of the G + C-rich 5' non-coding sequence and strong intron size divergence from homologous mammalian genes

The complete human dihydrofolate reductase (DHFR) gene has been cloned from four recombinant lambda libraries constructed with the DNA from a methotrexate-resistant human cell line with amplified DHFR genes. The detailed organization of the gene has been determined by restriction mapping of the cloned fragments and DNA sequencing of all the protein coding regions and adjacent intron segments, and shown to correspond to that of the native human DHFR gene. The gene spans a length of approximately 29 X 10(3) bases from the ATG initiator codon to the end of the 3' untranslated region, and contains five introns that interrupt the protein coding sequence. The number and positions of introns are identical to those found in the mouse gene. By contrast, the size of the homologous introns (with the exception of the first one) varies greatly, up to several fold, in the genes from man, mouse and Chinese hamster; the intron sequences also exhibit a great divergence, except in the junction regions. A striking sequence homology, extending over several hundred nucleotides, exists between the human and mouse gene 5' non-coding regions. These regions are characterized by an unusually high G + C content, 72% and 66% in the human and mouse genes, respectively, which is maintained in the first coding segment and first intron, and is in sharp contrast to the relatively low G + C content (approximately 40%) of the remainder of the gene.     

Prediction of functional sites by analysis of sequence and structure conservation

We present a method for prediction of functional sites in a set of aligned protein sequences. The method selects sites which are both well conserved and clustered together in space, as inferred from the 3D structures of proteins included in the alignment. We tested the method using 86 alignments from the NCBI CDD database, where the sites of experimentally determined ligand and/or macromolecular interactions are annotated. In agreement with earlier investigations, we found that functional site predictions are most successful when overall background sequence conservation is low, such that sites under evolutionary constraint become apparent. In addition, we found that averaging of conservation values across spatially clustered sites improves predictions under certain conditions: that is, when overall conservation is relatively high and when the site in question involves a large macromolecular binding interface. Under these conditions it is better to look for clusters of conserved sites than to look for particular conserved sites.     

Mouse asialoglycoprotein receptor cDNA sequence: conservation of receptor genes during mammalian evolution

The asialoglycoprotein receptor internalizes galactose-terminated glycoproteins into mammalian hepatocytes for degradation in lysosomes. We report the cloning and sequencing of one murine asialoglycoprotein receptor cDNA which exhibits homology with rat and human receptor forms. Conserved regions may correlate with functional domains. The membrane-bound M (mouse) HL polypeptide does not contain a cleavable N-terminal signal sequence and is probably anchored to the membrane via an internal insertion sequence.     

Molecular features and functional constraints in the evolution of the mammalian X chromosome

Recent advances in genomic sequencing of multiple organisms have fostered significant advances in our understanding of the evolution of the sex chromosomes. The integration of this newly available sequence information with functional data has facilitated a considerable refinement of our conceptual framework of the forces driving this evolution. Here we address multiple functional constraints that were encountered in the evolution of the X chromosome and the impact that this evolutionary history has had on its modern behavior.     

The tooth of perfection: functional and spatial constraints on mammalian tooth shape

This paper addresses the question of how close mammalian teeth are to ideal functional forms. An 'ideal' form is a morphology predicted to be the best functional shape according to information of the relationships between shape and function. Deviations from an ideal form are likely to indicate the presence of developmental or genetic constraints on form. Model tools were constructed to conform to functional principles from engineering and dental studies. The final model shapes are very similar to several mammalian tooth forms (carnassial teeth and tribosphenic-like cusps), suggesting that these tooth forms very closely approach ideal functional forms. Further evidence that these tooth forms are close to ideal comes from the conservation over 140 million years, the independent derivation and/or the occurrence over a size range of several orders of magnitude of these basic tooth forms. One of the main functional shapes derived here is the 'protoconoid', a fundamental design for double-bladed tools that fits a large number of functional parameters. This shape occurs in tooth forms such as tribosphenic, dilambdodont and zalambdodont. This study extends our understanding of constraints on tooth shape in terms of geometry (how space influences tooth shape) and function (how teeth divide food).  © 2003 The Linnean Society of London . Biological Journal of the Linnean Society , 2003, 78 , 173–191.     

The mammalian Rab family of small GTPases: definition of family and subfamily sequence motifs suggests a mechanism for functional specificity in the Ras superfamily

The Rab/Ypt/Sec4 family forms the largest branch of the Ras superfamily of GTPases, acting as essential regulators of vesicular transport pathways. We used the large amount of information in the databases to analyse the mammalian Rab family. We defined Rab-conserved sequences that we designate Rab family (RabF) motifs using the conserved PM and G motifs as "landmarks". The Rab-specific regions were used to identify new Rab proteins in the databases and suggest rules for nomenclature. Surprisingly, we find that RabF regions cluster in and around switch I and switch II regions, i.e. the regions that change conformation upon GDP or GTP binding. This finding suggests that specificity of Rab-effector interaction cannot be conferred solely through the switch regions as is usually inferred. Instead, we propose a model whereby an effector binds to RabF (switch) regions to discriminate between nucleotide-bound states and simultaneously to other regions that confer specificity to the interaction, possibly Rab subfamily (RabSF) specific regions that we also define here. We discuss structural and functional data that support this model and its general applicability to the Ras superfamily of proteins.     

Structural and functional conservation of Mycobacterium tuberculosis GroEL paralogs suggests that GroEL1 Is a chaperonin

GroEL is a group I chaperonin that facilitates protein folding and prevents protein aggregation in the bacterial cytosol. Mycobacteria are unusual in encoding two or more copies of GroEL in their genome. While GroEL2 is essential for viability and likely functions as the general housekeeping chaperonin, GroEL1 is dispensable, but its structure and function remain unclear.Here, we present the 2.2-Å resolution crystal structure of a 23-kDa fragment of Mycobacterium tuberculosis GroEL1 consisting of an extended apical domain. Our X-ray structure of the GroEL1 apical domain closely resembles those of Escherichia coli GroEL and M. tuberculosis GroEL2, thus highlighting the remarkable structural conservation of bacterial chaperonins. Notably, in our structure, the proposed substrate-binding site of GroEL1 interacts with the N-terminal region of a symmetry-related neighboring GroEL1 molecule. The latter is consistent with the known GroEL apical domain function in substrate binding and is supported by results obtained from using peptide array technology. Taken together, these data show that the apical domains of M. tuberculosis GroEL paralogs are conserved in three-dimensional structure, suggesting that GroEL1, like GroEL2, is a chaperonin.     

功能型分子标记(ISAP)的开发及评价

分子标记在图谱构建, QTL分析, 基因定位以及标记辅助育种中起着越来越重要的作用。研究者都期望一个分子标记位点代表一个特定的基因, 甚至与某种性状联系起来, 这样, 通过对某个分子标记的筛选即能对性状进行筛选, 此即功能型分子标记。而目前广泛使用的基于PCR基础的分子标记如RAPD、SSR、AFLP等或是扩增非编码区域, 或是随机在基因组中扩增, 得到的位点一般与目标性状基因距离较远,这使得分子标记在应用上与其目标有一定的偏差。研究建立了一种基于基因中内含子序列的功能型分子标记, 试图使标记位点与基因序列联系起来以达到其功能型的目的。它利用内含子剪接位置的保守一致序列作为其引物的核心序列,其上下游引物均为18 bp, 上下游引物间通过组合配对的方式作为扩增的引物对, 对真核生物的基因序列进行扩增。为了验证ISAP的功能性, 研究设计了17条引物(9条上游引物, 8条下游引物, 共计72个引物组合)对棉花F2群体进行扩增并构建遗传连锁图谱, 其中67个显示了多态性, 共得到212个位点。我们用此212个位点连同164个SRAP位点构建了一张包含276个位点的遗传连锁图谱, ISAP标记在整个连锁群中分布比较均匀,部分区域呈现标记高饱和现象, 可能为编码序列富集区。另外对20个片段进行测序的结果表明, 85％的序列显示了与已公布EST序列的同源性, 说明扩增是跨越了外显子进行的, 得到的序列与表达序列紧密连锁。结果显示, ISAP标记是简单, 可靠, 具有较高多态性, 并且扩增基因区域的一种功能型分子标记。同时, 还使用ISAP标记对其他植物进行了扩增, 取得了良好的效果。     

Phylogeny, sequence conservation, and functional complementation of the SBDS protein family

The Shwachman-Bodian-Diamond syndrome (SBDS) protein family occurs widely in nature, although its function has not been determined. Comprehensive database searches revealed SBDS homologues from 159 species, including examples from all sequenced archaeal and eukaryotic genomes and all eukaryotic kingdoms. Sequence alignment with ClustalX and MUSCLE algorithms led to the identification of conserved residues that occurred predominantly in the amino-terminal FYSH domain where they appeared to contribute to protein folding or stability. Only SBDS residue Gly91 was invariant in all species. Four distantly related protists were found to have two divergent SBDS genes in their genomes. In each case, phylogenetic analyses and the identification of shared sequence features suggested that one gene was derived from lateral gene transfer. We also identified a shared C-terminal zinc finger domain fusion in flowering plants and chromalveolates that may shed light on the function of the protein family and the evolutionary histories of these kingdoms. To assess the extent of SBDS functional conservation, we carried out complementation studies of SBDS homologues and interspecies chimeras in Saccharomyces cerevisiae. We determined that the FYSH domain was widely interchangeable among eukaryotes, while domain 2 imparted species specificity to protein function. Domain 3 was largely dispensable for function in our yeast complementation assay. Overall, the phylogeny of SBDS was shared with a group of proteins that were markedly enriched for RNA metabolism and/or ribosome-associated functions. These findings link Shwachman-Diamond syndrome to other bone marrow failure syndromes with defects in nucleolus-associated processes, including Diamond-Blackfan anemia, cartilage-hair hypoplasia, and dyskeratosis congenita.     

PEDF and the serpins: phylogeny, sequence conservation, and functional domains

Pigment epithelium derived factor (PEDF) is non-inhibitory serpin with neurotrophic and antiangiogenic functions. In this study, we have assembled PEDF sequences for 9 additional species by data base mining and performed cross-species alignment for 14 PEDF sequences to identify conserved structural domains. We found evolutionary conservation of a leader sequence, a single C-terminal glycosylation site, collagen-binding residues, and four specific conserved PEDF peptides. The C-terminus, 384--415 and an N-terminal region 78--95, show close homology with many other serpins, and there is strong conservation of 39 of 51 consensus key residues involved in serpin structure and function. Two peptide regions, 40--67 and 277--301, are unique to PEDF but conserved in all species. Conserved residues at the N-terminus, helix d (hD), and helix A (hA) of PEDF form a structure similar to the heparin-binding groove of other serpins. We identified a motif in PEDF that is homologous to the nuclear localization signals of other proteins. A bitopographical localization of PEDF was confirmed by immunocytochemistry and Western blots. Our results suggest that secretion is required for PEDF's activity, that PEDF can migrate to the nucleus, and that PEDF has structural and functional features more common with inhibitory serpins.     

High conservation of sequences involved in cystic fibrosis mutations in five mammalian species

Several mutations have been identified in the first nucleocide binding fold (NBF) of the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) gene. We have analyzed the DNA sequences of exons 10 and 11 in five different mammalian species, marmoset, mouse, cow, pig, and sheep; the amino acid conservation studied for nine disease mutations; and two “benign” mutations. For exon 10,87% homology at the DNA level and 93.5% at the amino acid level were found for these species. For exon 11, the lowest homology (70%), as found in mouse and the highest in marmoset (93%), whereas the amino acid sequence conservation ranged from 82.5 to 100%. All codons involved in CF mutations are highly conserved throughout evolution.