Comparative sequence analysis of the mouse and human Lgn1/SMA interval.

Human chromosome 5q11.2-q13.3 and its ortholog on mouse chromosome 13 contain candidate genes for an inherited human neurodegenerative disorder called spinal muscular atrophy (SMA) and for an inherited mouse susceptibility to infection with Legionella pneumophila (Lgn1). These homologous genomic regions also have unusual repetitive organizations that create practical difficulties in mapping and raise interesting issues about the evolutionary origin of the repeats. In an attempt to analyze this region in detail, and as a way to identify additional candidate genes for these diseases, we have determined the sequence of 179 kb of the mouse Lgn1/SMA interval. We have analyzed this sequence using BLAST searches and various exon prediction programs to identify potential genes. Since these methods can generate false-positive exon declarations, our alignments of the mouse sequence with available human orthologous sequence allowed us to discriminate rapidly among this collection of potential coding regions by indicating which regions were well conserved and were more likely to represent actual coding sequence. As a result of our analysis, we accurately mapped two additional genes in the SMA interval that can be tested for involvement in the pathogenesis of SMA. While no new Lgn1 candidates emerged, we have identified new genetic markers that exclude Smn as an Lgn1 candidate. In addition to providing important resources for studying SMA and Lgn1, our data provide further evidence of the value of sequencing the mouse genome as a means to help with the annotation of the human genomic sequence and vice versa.     

Evolutionary divergence of the APETALA1 and CAULIFLOWER proteins

Abstract APETALA1 (AP1) and CAULIFLOWER (CAL) are a pair of paralogous genes that were generated through the pre‐Brassicaceae whole‐genome duplication event. AP1 and CAL have both partially redundant and unique functions. Previous studies have shown that the K and C regions of their proteins are essential for the functional divergence. However, which differences in these regions are the major contributors and how the differences were accumulated remain unknown. In the present study, we compared the sequences of the two proteins and identified five gaps and 55 amino acid replacements between them. Investigation of genomic sequences further indicated that the differences in the proteins were caused by non‐synonymous substitutions and changes in exon–intron structures. Reconstruction of three‐dimensional structures revealed that the sequence divergence of AP1 and CAL has resulted in differences between the two in terms of the number, length, position and orientation of α‐helices, especially in the K and C regions. Comparisons of sequences and three‐dimensional structures of ancestral proteins with AP1 and CAL suggest that the ancestral AP1 protein experienced fewer changes, whereas the ancestral CAL protein accumulated more changes shortly after gene duplication, relative to their common ancestor. Thereafter, AP1‐like proteins experienced few mutations, whereas CAL‐like proteins were not conserved until the diversification of the Brassicaceae lineage I. This indicates that AP1‐ and CAL‐like proteins evolved asymmetrically after gene duplication. These findings provide new insights into the functional divergence of AP1 and CAL genes.     

Evolutionary conservation of gene structures of the Pax1/9 gene family

Based on amino acid sequence comparisons, Pax1 and Pax9 genes are considered to form a subgroup of vertebrate Pax genes. We show here that the gene structures of mouse Pax1, human PAX9 genes are similar to that of a single Pax1/9 related gene in Branchiostoma lanceolatum, AmphiPax1. This supports the hypothesis that Pax1 and Pax9 genes were derived from a single ancestral gene. A refined protein alignment of AmphiPax1, mouse Pax1 and human PAX9 proteins based on the determined exon boundaries indicates that sequence divergence at the C-termini may be related to the unique functions of the Pax1 and Pax9 genes in vertebrates. AmphiPax1 is expressed in adult amphioxus in the pharyngeal endoderm.     

Structures of mouse SOD1 and human/mouse SOD1 chimeras

Mutations in human copper-zinc superoxide dismutase (SOD1) cause an inherited form of amyotrophic lateral sclerosis (ALS). Inclusions enriched in pathogenic SOD1 accumulate in the spinal cords of transgenic mice expressing these proteins, but endogenous mouse SOD1 is not found as a component of these aggregates. In the accompanying paper, Karch and colleagues analyze aggregation propensities of human/mouse SOD1 chimeras in cell culture and identify two sequence elements in the human enzyme that seem to enhance its aggregation relative to the mouse enzyme. Here, we report the first structure of mouse SOD1 along with those of SOD1 chimeras in which residues 1-80 come from human SOD1 and residues 81-153 come from mouse SOD1 and vice versa. Taken together, the structural and cell-based data suggest a model in which residues Q42 and Q123 in mouse SOD1 modulate non-native SOD1-SOD1 intermolecular interactions at edge strands in the SOD1 Greek key β-barrel.     

DNA structures, repeat expansions and human hereditary disorders

Expansions of simple DNA repeats are responsible for more than two dozen hereditary disorders in humans, including fragile X syndrome, myotonic dystrophy, Huntington's disease, various spinocerebellar ataxias, Friedreich's ataxia and others. During the past decade, it became clear that unusual structural features of expandable repeats greatly contribute to their instability and could lead to their expansion. Furthermore, DNA replication, repair and recombination are implicated in the formation of repeat expansions, as shown in various experimental systems. The replication model of repeat expansion stipulates that unusual structures of expandable repeats stall replication fork progression, whereas extra repeats are added during replication fork restart. It also explains the bias toward repeat expansion or contraction that was observed in different organisms.     

Evolutionary significance of gene expression divergence

Recent large-scale studies of evolutionary changes in gene expression among mammalian species have led to the proposal that gene expression divergence may be neutral with respect to organismic fitness. Here, we employ a comparative analysis of mammalian gene sequence divergence and gene expression divergence to test the hypothesis that the evolution of gene expression is predominantly neutral. Two models of neutral gene expression evolution are considered: 1-purely neutral evolution (i.e., no selective constraint) of gene expression levels and patterns and 2-neutral evolution accompanied by selective constraint. With respect to purely neutral evolution, levels of change in gene expression between human-mouse orthologs are correlated with levels of gene sequence divergence that are determined largely by purifying selection. In contrast, evolutionary changes of tissue-specific gene expression profiles do not show such a correlation with sequence divergence. However, divergence of both gene expression levels and profiles are significantly lower for orthologous human-mouse gene pairs than for pairs of randomly chosen human and mouse genes. These data clearly point to the action of selective constraint on gene expression divergence and are inconsistent with the purely neutral model; however, there is likely to be a neutral component in evolution of gene expression, particularly, in tissues where the expression of a given gene is low and functionally irrelevant. The model of neutral evolution with selective constraint predicts a regular, clock-like accumulation of gene expression divergence. However, relative rate tests of the divergence among human-mouse-rat orthologous gene sets reveal clock-like evolution for gene sequence divergence, and to a lesser extent for gene expression level divergence, but not for the divergence of tissue-specific gene expression profiles. Taken together, these results indicate that gene expression divergence is subject to the effects of purifying selective constraint and suggest that it might also be substantially influenced by positive Darwinian selection.     

Evolutionary divergence of human chromosome 9 as revealed by the position of the ABL protooncogene in higher primates

Attempts to solve the fundamental questions regarding the descent of man are dogged by superstitions and unexamined orthodoxies. The origin of humans, established a decade ago based upon cytological analysis of ape chromosomes, continues to be called into question. Although molecular methods have provided a framework for tracing the paths of human evolution, conclusive evidence remains elusive. We have used a single ABL gene probe derived from human chromosome 9 to assess the direction of change in the equivalent ape chromosomes. This approach has resulted in a few surprises which again challenge the prevailing view of early primate evolution based solely on chromosome banding patterns. The ABL protooncogene is present on human chromosome 9 at band q34. Similar DNA sequences presumed to represent an ABL gene, are present on chromosome 11 in chimpanzee (Pan troglodytes) but at a different relative location, indicating that the mechanism of the origin of human chromosome 9 is far more complex than has previously been suggested. Nevertheless, in gorilla (Gorilla gorilla) and orangutan (Pongo pygmaeus), the equivalent to human chromosome band 9 q34 is apparently located on chromosome 13 at a putative telomeric position and no discernible differences could be established. Despite the presence of the ABL protooncogene on human equivalent ape chromosomes, molecular systematics will continue to generate enigmas in the evolutionary context until the entire genome is sequenced.     

Evolutionary divergence in the pan-Atlantic mangrove Avicennia germinans

Cuticular hydrocarbon composition was determined for 18 populations of Avicennia germinans (Avicenniaceae) collected from West Africa, Florida, the Pacific coast of Mexico and French Guyana. Variation in covariance structure among the hydrocarbons was evaluated from correlations among the traits for each of the populations and from a hierarchical common principal components analysis. Principal components ordinations of populations based on the 171 z * transformed correlation coefficients, onto which a minimal spanning tree was superimposed, suggested a network in which French Guyanan populations formed a nexus from which African populations differentiated biochemically from Floridan and Mexican groups. Common principal components analysis provided further support for the differentiation in covariance structure among these biochemical traits of African populations. Pacific coast Mexican populations appeared to have diverged less in covariation among traits from Atlantic American sources than did the northern (Floridan) and southern (French Guyanan) Atlantic American populations from one another. The complete isolation of Pacific coast sources after emergence of the Central American isthmus during the Pliocene provides a reference for comparison of levels of differentiation among regions that suggests an underlying pattern of divergence in this species. Partial Mantel matrix tests revealed significant geographic effects, with no significant climatic effects. Our results provide support for recognition of a major differentiation between Old and New World sources within the taxon and indicate the need for a full revision of this species.     

Evolutionary divergence plots of homologous proteins.

A simple and efficient method is described for analyzing quantitatively multiple protein sequence alignments and finding the most conserved blocks as well as the maxima of divergence within the set of aligned sequences. It consists of calculating the mean distance and the root-mean-square distance in each column of the multiple alignment, averaging the values in a window of defined length and plotting the results as a function of the position of the window. Due attention is paid to the presence of gaps in the columns. Several examples are provided, using the sequences of several cytochromes c, serine proteases, lysozymes and globins. Two distance matrices are compared, namely the matrix derived by Gribskov and Burgess from the Dayhoff matrix, and the Risler Structural Superposition Matrix. In each case, the divergence plots effectively point to the specific residues which are known to be essential for the catalytic activity of the proteins. In addition, the regions of maximum divergence are clearly delineated. Interestingly, they are generally observed in positions immediately flanking the most conserved blocks. The method should therefore be useful for delineating the peptide segments which will be good candidates for site-directed mutagenesis and for visualizing the evolutionary constraints along homologous polypeptide chains.     

Evolutionary divergence of giant tortoises in Galapagos

The giant tortoises in the Galapagos Archipelago diverge considerably in size, and in shape and other carapace characteristics. The saddleback morphotype is known only from insular faunas lacking large terrestrial predators (i.e. Galapagos and Mauritius) and in Galapagos is associated with xeric habitats where vertical feeding range and vertical reach in agonistic encounters are adaptive. The large domed morphotype is associated with relatively cool, mesic habitats where intraspecific competition for food and other resources may be less intense than in xeric habitats. Other external characteristics that differ between tortoise populations are also correlated with ecological variation. Tortoises have radiated into a mosaic of ecological conditions in the Galapagos but critical data are lacking on the role of genetic and environmental controls on phenotypic variation. Morphological divergence in tortoises is potentially a better indicator of present ecological conditions than of evolutionary relationships.     

Evolutionary divergence of chloroplast FAD synthetase proteins

Background  

Flavin adenine dinucleotide synthetases (FADSs) - a group of bifunctional enzymes that carry out the dual functions of riboflavin phosphorylation to produce flavin mononucleotide (FMN) and its subsequent adenylation to generate FAD in most prokaryotes - were studied in plants in terms of sequence, structure and evolutionary history.     

Evolutionary relationship of plant catalase genes inferred from exon-intron structures: isozyme divergence after the separation of monocots and dicots

In order to understand the molecular evolution of catalase genes in higher plants, we compared the exon-intron structures of 12 genomic sequences from six plant species. It was assumed that the putative single primordial catalase gene had seven introns, because only those catalase genes having this structure are found in the monocotyledonae and dicotyledonae classes. After the evolutionary divergence of monocots from dicots, consecutive duplication of the primordial gene followed by the differential loss of introns occurred in each class to form three (or possibly four in dicots) diverse isozyme genes. In monocots, three ancestral isozyme genes were formed before the divergence of ancestral rice and maize. One of the rice genes, CatA, has an entirely new short intron which was not found in any other plant catalase gene examined. We have investigated the existence of the intron in the CatA homolog in other rice species by polymerase chain reaction (PCR) analysis. One major PCR product was found with the genomic DNAs from O. sativa (indica and japonica types), O. rufipogon and O. glaberrima. DNAs from several accessions of O. longistaminata showed variation in both the number and size of the DNA fragments amplified. PCR analyses and sequencing of the PCR products revealed that there are several CatA homologs having different sequences in some accessions of O. longistaminata. We have extended our study to other species in the Poaceae. The results suggest that the gain of the intron, most likely by insertion of a retroposon, took place in the ancestral genome of rice after its evolutionary divergence from other ancestral cereals such as barley, wheat and oat. Received: 20 November 1997 / Accepted: 5 January 1998     

Global similarity and local divergence in human and mouse gene co-expression networks

Background  

A genome-wide comparative analysis of human and mouse gene expression patterns was performed in order to evaluate the evolutionary divergence of mammalian gene expression. Tissue-specific expression profiles were analyzed for 9,105 human-mouse orthologous gene pairs across 28 tissues. Expression profiles were resolved into species-specific coexpression networks, and the topological properties of the networks were compared between species.     

Crystal structures of mouse and human RP105/MD-1 complexes reveal unique dimer organization of the toll-like receptor family

The Toll-like receptor (TLR) 4/MD-2 heterodimer senses lipopolysaccharide (LPS). RP105 (radioprotective 105 kDa), a TLR-related molecule, is similar to TLR4 in that the extracellular leucine-rich repeats associate with MD-1, the MD-2-like molecule. MD-2 has a unique hydrophobic cavity that directly binds to lipid A, the active center of LPS. LPS-bound MD-2 opens the secondary interface with TLR4, leading to dimerization of TLR4/MD-2. MD-1 also has a hydrophobic cavity that accommodates lipid IVa, a precursor of lipid A, suggesting a role for the RP105/MD-1 heterodimer in sensing LPS or related microbial products. Little is known, however, about the structure of the RP105/MD-1 heterodimer or its oligomer. Here, we have determined the crystal structures of mouse and human RP105/MD-1 complexes at 1.9 and 2.8 Å resolutions, respectively. Both mouse and human RP105/MD-1 exhibit dimerization of the 1:1 RP105/MD-1 complex, demonstrating a novel organization. The “m”-shaped 2:2 RP105/MD-1 complex exhibits an inverse arrangement, with N-termini interacting in the middle. Thus, the dimerization interface of RP105/MD-1 is located on the opposite side of the complex, compared to the 2:2 TLR4/MD-2 complex. These results demonstrate that the 2:2 RP105/MD-1 complex is distinct from previously reported TLR dimers, including TLR4/MD-2, TLR1/TLR2, TLR2/TLR6, and TLR3, all of which facilitate homotypic or heterotypic interaction of the C-terminal cytoplasmic signaling domain.