Ultraconserved elements in insect genomes: A highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to Occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.

Characterization of highly conserved microsatellite loci in Araucaria cunninghamii and related species

Ten microsatellite loci are described in Araucaria cunninghamii, the first reported in the Araucariaceae. Eight were tested in sections Eutacta and Bunya, which diverged more than 200 MYA, and to the sister genus Agathis. Specific amplification products within the expected size range were obtained for six to eight loci in section Eutacta (depending on species), five loci in section Bunya and three. loci in Agathis. Two of the loci (CRCAc1 and CRCAc2, both GA repeats) produced specific amplification products in all taxa, with orthology confirmed by sequence analysis. The repeats were perfect in all taxa. The flanking sequences were extremely conserved, with sequence divergence of 0% to 2.0% within Araucaria species and 2.9% to 7.5% between Araucaria and Agathis. These microsatellites represent some of the most conserved microsatellite loci reported in plants. This may be due to a low evolutionary rate in Araucariaceae genome or the loci may be closely associated with highly conserved, unreported genes.

Identification of two evolutionarily conserved and functional regulatory elements in intron 2 of the human BRCA1 gene

Cross-species comparative genomics is a powerful strategy for identifying functional regulatory elements within noncoding DNA. In this paper, comparative analysis of human and mouse intronic sequences in the breast cancer susceptibility gene (BRCA1) revealed two evolutionarily conserved noncoding sequences (CNS) in intron 2, 5 kb downstream of the core BRCA1 promoter. The functionality of these elements was examined using homologous-recombination-based mutagenesis of reporter gene-tagged cosmids incorporating these regions and flanking sequences from the BRCA1 locus. This showed that CNS-1 and CNS-2 have differential transcriptional regulatory activity in epithelial cell lines. Mutation of CNS-1 significantly reduced reporter gene expression to 30% of control levels. Conversely mutation of CNS-2 increased expression to 200% of control levels. Regulation is at the level of transcription and shows promoter specificity. Both elements also specifically bind nuclear proteins in vitro. These studies demonstrate that the combination of comparative genomics and functional analysis is a successful strategy to identify novel regulatory elements and provide the first direct evidence that conserved noncoding sequences in BRCA1 regulate gene expression. (c) 2005 Elsevier Inc. All rights reserved.

The common tetratricopeptide repeat acceptor site for steroid receptor-associated immunophilins and Hop is located in the dimerization domain of hsp90

Structurally related tetratricopeptide repeat motifs in steroid receptor-associated immunophilins and the STI1 homolog, Hop, mediate the interaction with a common cellular target, hsp90, We have identified the binding domain in hsp90 for cyclophilin 40 (CyP40) using a two-hybrid system screen of a mouse cDNA library. All isolated clones encoded the intact carboxyl terminus of hsp90 and overlapped with a common region corresponding to amino acids 558-724 of murine hsp84, The interaction was confirmed in vitro with bacterially expressed CyP40 and deletion mutants of hsp90 beta and was delineated further to a 124-residue COOH-terminal segment of hsp90, Deletion of the conserved MEEVD sequence at the extreme carboxyl terminus of hsp90 precludes interaction with CyP40, signifying an important role for this motif in hsp90 function. We show that CyP40 and Hop display similar interaction profiles with hsp90 truncation mutants and present evidence for the direct competition of Hop and FK506-binding protein 52 with CyP40 for binding to the hsp90 COOH-terminal region. Our results are consistent with a common tetratricopeptide repeat interaction site for Hop and steroid receptor associated immunophilins within a discrete COOH-terminal domain of hsp90. This region of hsp90 mediates ATP-independent chaperone activity, overlaps the hsp90 dimerization domain, and includes structural elements important for steroid receptor interaction.

Cloning and functional analysis of the Sry-related HMG box gene, Sox18

The Sox gene family (Sry like HMG box gene) is characterised by a conserved DNA sequence encoding a domain of approximately 80 amino acids which is responsible for sequence specific DNA binding. We initially published the identification and partial cDNA sequence of murine Sox18, a new member of this gene family, isolated from a cardiac cDNA library. This sequence allowed us to classify Sox18 into the F sub-group of Sox proteins, along with Sox7 and Sox17. Recently, we demonstrated that mutations in the Sox18 activation domain underlie cardiovascular and hair follicle defects in the mouse mutation, ragged (Ra) (Pennisi et al., 2000. Mutations in Sox18 underlie cardiovascular and hair follicle defecs in ragged mice. Nat. Genet. 24, 434-437). Ra homozygotes lack vibrissae and coat hairs, have generalised oedema and an accumulation of chyle in the peritoneum. Here we have investigated the genomic sequences encoding Sox18. Screening of a mouse genomic phage library identified four overlapping clones, we sequenced a 3.25 kb XbaI fragment that defined the entire coding region and approximately 1.5 kb of 5' flanking sequences. This identified (i) an additional 91 amino acids upstream of the previously designated methionine start codon in the original cDNA, and (ii);ln intron encoded within the HMG box/DNA binding domain in exactly the same position as that found in the Sox5, -13 and -17 genes. The Sox18 gene encodes a protein of 468 aa. We present evidence that suggests HAF-2, the human HMG-box activating factor-2 protein, is the orthologue of murine Sox18. HAF-2 has been implicated in the regulation of the Human IgH enhancer in a B cell context. Random mutagenesis coupled with GAL4 hybrid analysis in the activation domain between amino acids 252 and 346, of Sox18, implicated the phosphorylation motif, SARS, and the region between amino acid residues 313 and 346 as critical components of Sox18 mediated transactivation. Finally, we examined the expression of Sox18 in multiple adult mouse tissues using RT-PCR. Low-moderate expression was observed in spleen, stomach, kidney, intestine, skeletal muscle and heart. Very abundant expression was detected in lung tissue. (C) 2001 Elsevier Science B.V. All rights reserved.

Conserved structural and sequence elements implicated in the processing of gene-encoded circular proteins

The cyclotides are the largest family of naturally occurring circular proteins. The mechanism by which the termini of these gene-encoded proteins are linked seamlessly with a peptide bond to form a circular backbone is unknown. Here we report cyclotide-encoding cDNA sequences from the plant Viola odorata and compare them with those from an evolutionarily distinct species, Oldenlandia affinis. Individual members of this multigene family encode one to three mature cyclotide domains. These domains are preceded by N-terminal repeat regions (NTRs) that are conserved within a plant species but not between species. We have structurally characterized peptides corresponding to these NTRs and show that, despite them having no sequence homology, they form a structurally conserved alpha-helical motif. This structural conservation suggests a vital role for the NTR in the in vivo folding, processing, or detoxification of cyclotide domains from the precursor protein.

Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function

The mammalian transcriptome contains many nonprotein-coding RNAs (ncRNAs), but most of these are of unclear significance and lack strong sequence conservation, prompting suggestions that they might be non-functional. However, certain long functional ncRNAs such as Air and Xist are also poorly conserved. In this article, we systematically analyzed the conservation of several groups of functional ncRNAs, including miRNAs, snoRNAs and longer ncRNAs whose function has been either documented or confidently predicted. As expected, miRNAs and snoRNAs were highly conserved. By contrast, the longer functional non-micro, non-sno ncRNAs were much less conserved with many displaying rapid sequence evolution. Our findings suggest that longer ncRNAs are under the influence of different evolutionary constraints and that the lack of conservation displayed by the thousands of candidate ncRNAs does not necessarily signify an absence of function.

Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA

Two small RNAs regulate the timing of Caenorhabditis elegans development(1,2). Transition from the first to the second larval stage fates requires the 22-nucleotide lin-4 RNA(1,3,4), and transition from late larval to adult cell fates requires the 21-nucleotide let-7 RNA 2. The lin-4 and let-7 RNA genes are not homologous to each other, but are each complementary to sequences in the 3' untranslated regions of a set of protein-coding target genes that are normally negatively regulated by the RNAs1,2,5,6. Here we have detected let-7 RNAs of similar to 21 nucleotides in samples from a wide range of animal species, including vertebrate, ascidian, hemichordate, mollusc, annelid and arthropod, but not in RNAs from several cnidarian and poriferan species, Saccharomyces cerevisiae, Escherichia coli or Arabidopsis. We did not detect lin-4 RNA in these species. We found that let-7 temporal regulation is also conserved: let-7 RNA expression is first detected at late larval stages in C. elegans and Drosophila, at 48 hours after fertilization in zebrafish, and in adult stages of annelids and molluscs. The let-7 regulatory RNA may control late temporal transitions during development across animal phylogeny.

The uncharged surface features surrounding the active site of Escherichia coli DsbA are conserved and are implicated in peptide binding

DsbA is a protein-folding catalyst from the periplasm of Escherichia coli that interacts with newly translocated polypeptide substrate and catalyzes the formation of disulfide bonds in these secreted proteins. The precise nature of the interaction between DsbA and unfolded substrate is not known. Here, we give a detailed analysis of the DsbA crystal structure, now refined to 1.7 Angstrom, and present a proposal for its interaction with peptide. The crystal structure of DsbA implies flexibility between the thioredoxin and helical domains that may be an important feature for the disulfide transfer reaction. A hinge point for domain motion is identified-the typo IV beta-turn Phe 63-Met 64-Gly 65-Gly 66, which connects the two domains. Three unique features on the active site surface of the DsbA molecule-a groove, hydrophobic pocket, and hydrophobic patch-form an extensive uncharged surface surrounding the active-sits disulfide. Residues that contribute to these surface features are shown to be generally conserved in eight DsbA homologues. Furthermore, the residues immediately surrounding the active-site disulfide are uncharged in all nine DsbA proteins. A model for DsbA-peptide interaction has been derived from the structure of a human thioredoxin:peptide complex. This shows that peptide could interact with DsbA in a manner similar to that with thioredoxin. The active-site disulfide and all three surrounding uncharged surface features of DsbA could, in principle, participate in the binding or stabilization of peptide.

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% AST). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.

The sequence of a 200 kb portion of a Plasmodium vivax chromosome reveals a high degree of conservation with Plasmodium falciparum chromosome 3

Within a 199 866 base pair (bp) portion of a Plasmodium vivax chromosome we identified a conserved linkage group consisting of at least 41 genes homologous to Plasmodium falciparum genes located on chromosome 3. There were no P. vivax homologues of the P. falciparum cytoadherence-linked asexual genes clag 3.2, clag 3.1 and a var C pseudogene found on the P. vivax chromosome. Within the conserved linkage group, the gene order and structure are identical to those of P. falciparum chromosome 3. This conserved linkage group may extend to as many as 190 genes. The subtelomeric regions are different in size and the P. vivax segment contains genes for which no P. falciparum homologues have been identified to date. The size difference of at least 900 kb between the homologous P. vivax chromosome and P. falciparum chromosome 3 is presumably due to a translocation. There is substantial sequence divergence with a much higher guanine + cytosine (G + C) content in the DNA and a preference for amino acids using GC-rich codons in the deduced proteins of P. vivax. This structural conservation of homologous genes and their products combined with sequence divergence at the nucleotide level makes the P. vivax genome a powerful tool for comparative analyses of Plasmodium genomes. (C) 2001 Elsevier Science B.V. All rights reserved.

Ex vivo analysis of T-cell responses to Epstein-Barr virus-encoded oncogene latent membrane protein 1 reveals highly conserved epitope sequences in virus isolates from diverse geographic regions

Epstein-Barr virus (EBV)-encoded oncogene latent membrane protein (LMP) 1, which is consistently expressed in multiple EBV-associated malignancies, has been proposed as a potential target antigen for any future vaccine designed to control these malignancies. However, the high degree of genetic variation in the LMP1 sequence has been considered a major impediment for its use as a potential immunotherapeutic target for the treatment of EBV-associated malignancies. In the present study, we have employed a highly efficient strategy, based on ex vivo functional assays, to conduct an extensive sequence-wide analysis of LMP1-specific T-cell responses in a large panel of healthy virus carriers of diverse ethnic origin and nasopharyngeal carcinoma patients. By comparing the frequencies of T cells specific for overlapping peptides spanning LMP1, we mapped a number of novel HLA class I- and class II-restricted LMP1 T-cell epitopes, including an epitope with dual HLA class I restriction. More importantly, extensive sequence analysis of LMP1 revealed that the majority of the T-cell epitopes were highly conserved in EBV isolates from Caucasian, Papua New Guinean, African, and Southeast Asian populations, while unique geographically constrained genetic variation was observed within one HLA A2 supertype-restricted epitope. These findings indicate that conserved LMP1 epitopes should be considered in designing epitope-based immunotherapeutic strategies against EBV-associated malignancies in different ethnic populations.

Etiology of ovarian failure in blepharophimosis ptosis epicanthus inversus syndrome: FOXL2 is a conserved, early-acting gene in vertebrate ovarian development

Blepharophimosis ptosis epicanthus inversus syndrome (BPES) is a human disorder caused by mutations in the forkhead transcription factor gene FOXL2 and is characterized by facial dysmorphology combined in some cases with ovarian failure. To better understand the role of FOXL2 in the etiology of ovarian failure in BPES, we examined its expression in embryonic ovaries of mice, chickens, and red-eared slider turtles, representatives of three phylogenetically distant vertebrate groups that have different mechanisms of sex determination. Expression of Foxl2 was detected in early ovaries of all three species around the time of sex determination and was associated with both somatic and germ cell populations in mice. Expression was sexually dimorphic in all cases. Sequence analysis of turtle and chicken FoxL2 orthologues indicated an unusually high degree of structural conservation during evolution. FoxL2 was found to be autosomal in chickens, and therefore unlikely to represent the dominant ovarian-determining gene that has been postulated to exist as a possible explanation for female heterogamety in birds. Our observations suggest that BPES may result from early abnormalities in regulating the development of the fetal ovary, rather than premature degeneration of the postnatal or adult ovary. Further, our results suggest that FOXL2 is a highly conserved early regulator of vertebrate ovarian development.

Functional recycling of C2 domains throughout evolution: A comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level,the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various subfamilies. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces' our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics. (C) 2003 Elsevier Ltd. All rights reserved.