Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Orthologous analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources
PubMed2006

Identification and analysis of genes and pseudogenes within duplicated regions in the human and mouse genomes.

The identification and classification of genes and pseudogenes in duplicated regions still constitutes a challenge for standard automated genome annotation procedures. Using an integrated homology and orthology analysis independent of current gene annotation, we have identified 9,484 and 9,017 gene duplicates in human and mouse, respectively. On the basis of the integrity of their coding regions, we have classified them into functional and inactive duplicates, allowing us to define the first consistent and comprehensive collection of 1,811 human and 1,581 mouse unprocessed pseudogenes. Furthermore, of the total of 14,172 human and mouse duplicates predicted to be functional genes, as many as 420 are not included in current reference gene databases and therefore correspond to likely novel mammalian genes. Some of these correspond to partial duplicates with less than half of the length of the original source genes, yet they are conserved and syntenic among different mammalian lineages. The genes and unprocessed pseudogenes obtained here will enable further studies on the mechanisms involved in gene duplication as well as of the fate of duplicated genes.

Animals↗
PubMed2001

Automatic clustering of orthologs and in-paralogs from pairwise species comparisons.

Orthologs are genes in different species that originate from a single gene in the last common ancestor of these species. Such genes have often retained identical biological roles in the present-day organisms. It is hence important to identify orthologs for transferring functional information between genes in different organisms with a high degree of reliability. For example, orthologs of human proteins are often functionally characterized in model organisms. Unfortunately, orthology analysis between human and e.g. invertebrates is often complex because of large numbers of paralogs within protein families. Paralogs that predate the species split, which we call out-paralogs, can easily be confused with true orthologs. Paralogs that arose after the species split, which we call in-paralogs, however, are bona fide orthologs by definition. Orthologs and in-paralogs are typically detected with phylogenetic methods, but these are slow and difficult to automate. Automatic clustering methods based on two-way best genome-wide matches on the other hand, have so far not separated in-paralogs from out-paralogs effectively. We present a fully automatic method for finding orthologs and in-paralogs from two species. Ortholog clusters are seeded with a two-way best pairwise match, after which an algorithm for adding in-paralogs is applied. The method bypasses multiple alignments and phylogenetic trees, which can be slow and error-prone steps in classical ortholog detection. Still, it robustly detects complex orthologous relationships and assigns confidence values for both orthologs and in-paralogs. The program, called INPARANOID, was tested on all completely sequenced eukaryotic genomes. To assess the quality of INPARANOID results, ortholog clusters were generated from a dataset of worm and mammalian transmembrane proteins, and were compared to clusters derived by manual tree-based ortholog detection methods. This study led to the identification with a high degree of confidence of over a dozen novel worm-mammalian ortholog assignments that were previously undetected because of shortcomings of phylogenetic methods.A WWW server that allows searching for orthologs between human and several fully sequenced genomes is installed at http://www.cgb.ki.se/inparanoid/. This is the first comprehensive resource with orthologs of all fully sequenced eukaryotic genomes. Programs and tables of orthology assignments are available from the same location.

Algorithms↗
PubMed2026

Insights into glandular trichome biology from analysis of organ-specific gene expression programmes in cannabis, hop and tomato.

Glandular trichomes (GTs) are epidermal outgrowths in which diverse specialised (secondary) metabolites are synthesised and stored. Cannabis (Cannabis sativa L.) and its close relative hop (Humulus lupulus L.) have pharmaceutical and industrial significance due to the presence of these metabolites in their GTs. We examined the conservation or divergence of the specific transcriptional programmes underlying GT biology. To achieve this, we generated transcriptome atlases of trichomes, flower, leaf, stem and root for cannabis, hop and tomato. We found that 12.9, 10.1 and 16.8% of cannabis, hop and tomato genes, respectively, were expressed organ/tissue specifically across all organs/tissues. Transcription factors (TFs) on average accounted for 7.5% of the organ-specific transcriptome and likely regulate organ-specific functions. We also conducted weighted gene co-expression network analysis and gene regulatory network (GRN) analysis to identify key regulators of GT function across the species and validated our predictions by DNA affinity purification sequencing for a subset of the cannabis and tomato GT TFs. The GRNs specific to cannabis or hop GTs were enriched for TFs and target genes associated with specialised metabolism, reflecting their species-specific nature. Conversely, the shared GRN components (identified via orthology analysis) were involved in highly conserved processes, such as flavonoid biosynthesis, solute transport and metabolite storage. Together, these GRNs and the associated transcriptome atlases are valuable resources to improve our knowledge of GT function and organ-specific genome regulation.

Solanum lycopersicum↗
PubMed2006

Common genomic structure for the Lepidoptera cadherin-like genes.

A cadherin-like protein present in the midgut epithelial cells of Lepidoptera is associated with insect resistance to Bacillus thuringiensis Cry toxins. We describe for the first time the genes that encode the cadherin-like proteins in Ostrinia nubilalis, Helicoverpa armigera, and Bombyx mori, and analyze their organization. These genes encompass 19.6 kb, 20.0 kb, and 41.8 kb of genomic DNA, respectively, and despite the size heterogeneity, they are all composed of 35 exons that are linked by 34 introns. In contrast to the high variability noted for the sizes of the introns, the sizes of the coding exons were almost completely preserved among the three species, because the intronic sequences (except the first and last one) were inserted in homologous positions in the respective cDNA sequences. In all cases, the first intron (the longest one) was located in the 5'-untranslated region. These results point out a highly conserved structure that indicates that these genes are orthologous. Analysis of the gene sequences showed a common GC-rich region located in exon 31 and several tandem repeats that were specific for each gene. Partial sequences from transposable elements were found only in B. mori and accounted, in part, for the large size of the insect gene. The in silico analysis of the protein products of these genes showed a common structure and grouped these proteins into the protocadherin family.

3' Untranslated Regions↗
PubMed1999

Genome size as a mutation-selection-drift process.

A novel method for estimating neutral rates and patterns of DNA evolution in Drosophila takes advantage of the propensity of non-LTR retrotransposable elements to create nonfunctional, transpositionally inactive copies as a product of transposition. For many LINE elements, most copies present in a genome at any one time are nonfunctional "dead-on-arrival" (DOA) copies. Because these are off-shoots of active, transpositionally competent "master" lineages, in a gene tree of a LINE element from multiple samples from related species, the DOA lineages are expected to map to the terminal branches and the active lineages to the internal branches, the primary exceptions being when the sample includes DOA copies that are allelic or orthologous. Analysis of nucleotide substitutions and other changes along the terminal branches therefore allows estimation of the fixation process in the DOA copies, which are unconstrained with respect to protein coding; and under selective neutrality, the fixation process estimates the underlying mutational pattern. We have studied the retroelement Helena in Drosophila. An unexpectedly high rate of DNA loss was observed, yielding a half-life of unconstrained DNA sequences approximately 60-fold faster in Drosophila than in mammals. The high rate of DNA loss suggests a straightforward explanation of the seeming paradox that Drosophila has many fewer pseudogenes than found in mammalian species. Differential rates of deletion in different taxa might also contribute to the celebrated C-value paradox of why some closely related organisms can have very different DNA contents. New data presented here rule out the possibility that the transposition process itself is highly mutagenic, hence the observed linear relation between number of deletions and number of nucleotide substitutions is most easily explained by the hypothesis that both types of changes accumulate in unconstrained sequences over time.

Animals↗
PubMed1997

Phylogenetic footprinting of hypersensitive site 3 of the beta-globin locus control region.

Hypersensitive site 3 (HS3) of the beta-like globin locus control region has been implicated as an important regulator of the beta-like globin genes, but the trans factors that bind HS3 have only been partially characterized. Using a five-species alignment (human, galago, rabbit, goat, and mouse) that represents 370 million years of evolution, we have identified 24 phylogenetic footprints in the HS3 core and surrounding regions. Probes corresponding to the human sequence at each footprint have been used in binding studies to identify the nuclear factors that bind within and near these conserved sequence elements. Among the high-affinity interactions observed were several binding sites for proteins with repressor activity, including YY1, CCAAT displacement protein, and G1/G2 complexes (uncharacterized putative repressors) and several binding sites for the stage selector protein. To complement this analysis, orthologous galago sequences were also used to derive probes and the pattern of proteins binding to human and galago probes was compared. Binding interactions differing between these two species could be responsible for the different expression patterns shown by the two gamma genes (galago gamma is embryonic; human gamma is fetal). Alternatively, binding interactions that are conserved in the two species may be important in the regulation of common expression patterns (eg, repression of gamma in adult life).

Animals↗
PubMed1999

cDNA and genomic cloning of mouse aquaporin-2: functional analysis of an orthologous mutant causing nephrogenic diabetes insipidus.

As the first step in generating a transgenic mouse model of nephrogenic diabetes insipidus (NDI), we have analyzed the mouse aquaporin-2 (Aqp2) cDNA and gene and generated a mutated Aqp2 orthologous to NDI-causing human AQP2-T126M. Aqp2 cDNA was isolated from mouse kidney and encoded a 271-amino-acid protein with 90.4% identity to human AQP2. Expression in Xenopus oocytes indicated that Aqp2 encoded a mercurial-sensitive, water-selective channel. Northern blot analysis showed a single 1.7-kb Aqp2 transcript expressed only in kidney (medulla > cortex); transcript expression was increased approximately 20-fold in 48-h water-deprived mice. Immunoblot analysis revealed a 29-kDa glycoprotein in mouse kidney. Sequence comparison of the Aqp2 cDNA with a 5.5-kb mouse genomic DNA indicated three introns (lengths 2.4, 0.9, and 0.6 kb) separating four exons with boundaries at amino acids 120, 175, and 202. Genomic Southern blot analysis revealed a single-copy Aqp2 gene. The mutant Aqp2-T126M was water permeable when expressed in Xenopus oocytes, but was retained at the endoplasmic reticulum (ER) in transfected mammalian cells. The chemical chaperone glycerol produced a redistribution of Aqp2-T126M from ER to plasma membrane/endosomes. These results establish a basis for an Aqp2-T126M transgenic knock-in model of NDI.

Amino Acid Sequence↗
PubMed1996

Comparative analysis of 1196 orthologous mouse and human full-length mRNA and protein sequences.

A large set of mRNA and encoded protein sequences, from orthologous murine and human genes, was compiled to analyze statistical, biological, and evolutionary properties of coding and noncoding transcribed sequences. Protein sequence conservation varied between 36% and 100% identity, with an average value of 85%. The average degree of nucleotide sequence identity for the corresponding coding sequences was also approximately 85%, whereas 5' and 3' untranslated regions (UTRs) were less conserved, with aligned identities of 67% and 69%, respectively. For some mouse and human genes, nucleotide sequences are more highly conserved than the encoded protein sequences. A subset of 32 sequences, consisting of only mouse/human protein pairs for which the human sequence represents a positionally cloned disease gene, had properties very similar to the larger data set, suggesting that our data are representative of the genome as a whole. With respect to sequence conservation, two interesting outliers are the breast cancer (BRCAI) gene product and the testis-determining factor (SRY), both of which display among the lowest degrees of sequence identity. The occurrence of both introns and repetitive elements (e.g., Alu, Bl) in 5' and 3' UTRs was also studied. These results provide one benchmark for the "comparative genomics" of mice and humans, with practical implications for the cross-referencing of transcript maps. Also, they should prove useful in estimating the additional sampling diversity provided by mouse EST sequencing projects designed to complement the existing human cDNA collection.

Amino Acid Sequence↗
PubMed2003

A comparative analysis of an orthologous proteomic environment in the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe.

The sequential application of protein tagging, affinity purification, and mass spectrometry enables highly accurate charting of proteomic environments by the characterization of stable protein assemblies and the identification of subunits that are shared between two or more protein complexes, termed here "proteomic hyperlinks." We have charted the proteomic environments surrounding the histone methyltransferase, Set1, in both yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. Although the composition of these nonessential Set1 complexes is remarkably conserved, they differ with respect to their hyperlinks to their proteomic environments. We speculate that conservation of the core components of protein assemblies and variability of hyperlinks represents a general principle in the molecular organization of eukaryotic proteomes.

DNA-Binding Proteins↗
PubMed1999

Human and nematode orthologs--lessons from the analysis of 1800 human genes and the proteome of Caenorhabditis elegans.

Recently, we have defined and analyzed over 1800 orthologous human and rodent genes. Here we extend this work to compare human and Caenorhabditis elegans coding sequences. 1880 human proteins were compared with about 20000 predicted nematode proteins presumably comprising nearly the complete proteome of C. elegans. We found that 44% of human/rodent orthologs have convincing nematode counterparts. On average, the amino acid similarity and identity between aligned human and C. elegans orthologous gene products are 69.3% and 49.1% respectively, and the nucleotide identity is 49.8%. Detailed investigation of our results suggests that some nematode gene predictions are incorrect, leading to erroneous pairing with human genes (e.g. calcineurin and polymerase II elongation factor III). Furthermore, other proteins (i.e. homologs of human ribosomal proteins S20 and L41, thymosin) are missing entirely from the nematode proteome, suggesting that it may not be complete. These results underscore the fact that metazoan gene prediction is a very challenging task and that most computer-predicted nematode genes require supporting evidence of their existence from comparative genomics and/or laboratory investigation.

Animals↗
PubMed2002

Identification and characterization of a Brucella abortus ATP-binding cassette transporter homolog to Rhizobium meliloti ExsA and its role in virulence and protection in mice.

Brucella abortus is a facultative intracellular bacterial pathogen that causes abortion in domestic animals and undulant fever in humans. The mechanism of virulence of Brucella spp. is not fully understood yet. Furthermore, genes that allow Brucella to reach the intracellular niche and to interact with host cells need to be identified. Using the genomic survey sequence (GSS) approach, we identified the gene encoding an ATP-binding cassette (ABC) transporter of B. abortus strain S2308. The deduced amino acid sequence encoded by this gene exhibited 69 and 67% identity with the sequences of the ABC transporters encoded by the exsA genes of Rhizobium meliloti and Mesorhizobium loti, respectively. Additionally, B. abortus ExsA, like R. meliloti and M. loti ExsA, possesses ATP-binding motifs and the ABC signature domain features of a typical ABC transporter. Furthermore, ortholog group analysis placed B. abortus ExsA in ortholog group 6 of ABC transporters more likely to be involved in bacterial pathogenesis. In R. meliloti, ExsA is an exopolysaccharide transporter essential for alfalfa root nodule invasion and establishment of infection. To test the role of ExsA in Brucella pathogenesis, an exsA deletion mutant was constructed. Replacement of the wild-type exsA by recombination was demonstrated by Southern blot analysis of Brucella genomic DNA. Decreased survival in mice of the Brucella DeltaexsA mutant compared to the survival of parental strain S2308 demonstrated that ExsA is critical for full bacterial virulence. Additionally, the B. abortus exsA deletion mutant was used as a live vaccine. Challenge experiments revealed that the exsA mutant strain induced superior protective immunity in BALB/c mice compared to the protective immunity induced by strain S19 or RB51.

ATP-Binding Cassette Transporters↗
PubMed2003

Systematic variation of amino acid substitutions for stringent assessment of pairwise covariation.

During protein evolution, amino acids change due to a combination of functional constraints and genetic drift. Proteins frequently contain pairs of amino acids that appear to change together (covariation). Analysis of covariation from naturally occurring sets of orthologs cannot distinguish between residue pairs retained by functional requirements of the protein and those pairs existing due to changes along a common evolutionary path. Here, we have separated the two types of covariation by independently recombining every naturally occurring amino acid variant within a set of 15 subtilisin orthologs. Our analysis shows that in this family of subtilisin orthologs, almost all possible pairwise combinations of amino acids can coexist. This suggests that amino acid covariation found in the subtilisin orthologs is almost entirely due to common ancestral origin of the changes rather than functional constraints. We conclude that naturally occurring sequence diversity can be used to identify positions that can vary independently without destroying protein function.

Amino Acid Sequence↗
PubMed2005

PhyloGibbs: a Gibbs sampling motif finder that incorporates phylogeny.

A central problem in the bioinformatics of gene regulation is to find the binding sites for regulatory proteins. One of the most promising approaches toward identifying these short and fuzzy sequence patterns is the comparative analysis of orthologous intergenic regions of related species. This analysis is complicated by various factors. First, one needs to take the phylogenetic relationship between the species into account in order to distinguish conservation that is due to the occurrence of functional sites from spurious conservation that is due to evolutionary proximity. Second, one has to deal with the complexities of multiple alignments of orthologous intergenic regions, and one has to consider the possibility that functional sites may occur outside of conserved segments. Here we present a new motif sampling algorithm, PhyloGibbs, that runs on arbitrary collections of multiple local sequence alignments of orthologous sequences. The algorithm searches over all ways in which an arbitrary number of binding sites for an arbitrary number of transcription factors (TFs) can be assigned to the multiple sequence alignments. These binding site configurations are scored by a Bayesian probabilistic model that treats aligned sequences by a model for the evolution of binding sites and "background" intergenic DNA. This model takes the phylogenetic relationship between the species in the alignment explicitly into account. The algorithm uses simulated annealing and Monte Carlo Markov-chain sampling to rigorously assign posterior probabilities to all the binding sites that it reports. In tests on synthetic data and real data from five Saccharomyces species our algorithm performs significantly better than four other motif-finding algorithms, including algorithms that also take phylogeny into account. Our results also show that, in contrast to the other algorithms, PhyloGibbs can make realistic estimates of the reliability of its predictions. Our tests suggest that, running on the five-species multiple alignment of a single gene's upstream region, PhyloGibbs on average recovers over 50% of all binding sites in S. cerevisiae at a specificity of about 50%, and 33% of all binding sites at a specificity of about 85%. We also tested PhyloGibbs on collections of multiple alignments of intergenic regions that were recently annotated, based on ChIP-on-chip data, to contain binding sites for the same TF. We compared PhyloGibbs's results with the previous analysis of these data using six other motif-finding algorithms. For 16 of 21 TFs for which all other motif-finding methods failed to find a significant motif, PhyloGibbs did recover a motif that matches the literature consensus. In 11 cases where there was disagreement in the results we compiled lists of known target genes from the literature, and found that running PhyloGibbs on their regulatory regions yielded a binding motif matching the literature consensus in all but one of the cases. Interestingly, these literature gene lists had little overlap with the targets annotated based on the ChIP-on-chip data. The PhyloGibbs code can be downloaded from http://www.biozentrum.unibas.ch/~nimwegen/cgi-bin/phylogibbs.cgi or http://www.imsc.res.in/~rsidd/phylogibbs. The full set of predicted sites from our tests on yeast are available at http://www.swissregulon.unibas.ch.

Algorithms↗
PubMed2000

Isolation, characterization, expression and functional analysis of the zebrafish ortholog of MEN1.

Mutations in the MEN1 gene lead to an autosomal dominant disorder, multiple endocrine neoplasia type 1 (MEN1), which is characterized by tumors of the parathyroid, entero-pancreatic neuroendocrine, and pituitary tissues. The protein encoded by MEN1, 610-amino acid menin, resides primarily in the nucleus and binds to the transcription factor JunD, resulting in the repression of JunD-induced transcription. We report here a detailed characterization of the zebrafish men1 gene and its full-length (2551 nt) transcript, encoding a 617-amino acid protein with 67% identity and 80% similarity to human menin. Of the 81 missense mutations and in-frame deletions reported in MEN1 patients, 72 occur in residues that are identical in zebrafish, suggesting the importance of the conserved regions. The zebrafish men1 gene maps 61 cM from the top of linkage group 7 (LG7), a region that appears to show conserved synteny to the MEN1 loci at human 11q13. A 2.7-kb men1 message is detected at all stages of zebrafish development analyzed, from one-cell embryos to adult fish. Whole-mount in situ hybridization showed ubiquitous distribution of men1 message in zebrafish embryos at cleavage, blastula, gastrula, and early segmentation stages, with relatively abundant expression in blood cell progenitors (24 h post fertilization) and mesenchymal tissues (48 h post fertilization) at later stages. Zebrafish menin binds both human and mouse JunD, and represses JunD-induced transcription, indicating that the JunD-binding ability of menin is evolutionarily conserved.

Amino Acid Sequence↗
PubMed2002

Cloning and expression analysis of the chick ortholog of TBX22, the gene mutated in X-linked cleft palate and ankyloglossia.

T-box genes constitute a conserved gene family with important roles in many developmental processes. Several family members have been implicated in human congenital diseases. Recently, mutations in TBX22 were found to cause X-linked cleft palate (CPX and ankyloglossia), a semidominant X-linked disorder affecting formation of the secondary palate. Here, we have cloned the chick ortholog of human TBX22 and have analyzed its expression during embryogenesis. Expression is very prominent in the somites and in the myotome, and in the mandible and maxilla of the developing jaw. Other sites of expression include the limbs, the cranial mesenchyme and the eye. Hence, Tbx22 expression domains encompass the regions important for the development of the disease phenotype.

Amino Acid Sequence↗
PubMed1999

Mapping and expression analysis of the mouse ortholog of Xenopus Eomesodermin.

The T-box gene family has been conserved throughout metazoan evolution. The genes code for putative transcription factors which share a uniquely defining DNA binding domain, known as the T-box ([Bollag et al., 1994]). They are implicated in the control of diverse developmental processes by their highly specific expression patterns throughout gastrulation and organogenesis in mouse and other species ([Chapman et al., 1996]) ([Gibson-Brown et al., 1998]), and by mutations in T-box genes that have profound developmental effects ([Papaioannou, 1997]; [Chapman and Papaioannou, 1998]; [Papaioannou and Silver, 1998]). In this report, we describe the mapping and expression pattern of the mouse ortholog of a gene, Eomesodermin, first identified in Xenopus ([Ryan et al., 1996]). The mouse gene was previously reported ([Wattler et al., 1998]) under the name MmEomes. The gene maps to mouse chromosome 9 in a region syntenic with human chromosome 3p. Mouse eomesodermin is expressed in the trophoblast of the blastocyst and in its derivative, the chorionic ectoderm. At gastrulation, eomesodermin is expressed in the primitive streak and embryonic mesoderm as well, but this expression disappears prior to the end of gastrulation. Later, eomesodermin is expressed in the developing forebrain, in a pattern largely overlapping a closely related T-box gene, Tbr1 ([Bulfone et al., 1995]), and is also seen in a localized area of each limb.

Animals↗
PubMed1998

Evolutionary parameters of the transcribed mammalian genome: an analysis of 2,820 orthologous rodent and human sequences.

We have rigorously defined 2,820 orthologous mRNA and protein sequence pairs from rats, mice, and humans. Evolutionary rate analyses indicate that mammalian genes are evolving 17-30% more slowly than previous textbook values. Data are presented on the average properties of mRNA and protein sequences, on variations in sequence conservation in coding and noncoding regions, and on the absolute and relative frequencies of repetitive elements and splice sites in untranslated regions of mRNAs. Our data set contains 1,880 unique human/rodent sequence pairs that represent about 2-4% of all mammalian genes. Of the 1,880 human orthologs, 70% are present on a new gene map of the human genome, thus providing a valuable resource for cross-referencing human and rodent genomes. In addition to comparative mapping, these results have practical applications in the interpretation of noncoding sequence conservation between syntenic regions of human and mouse genomic sequence, and in the design and calibration of gene expression arrays.

Animals↗