Molecules and Morphologies: The Perpetuation of Pattern in the Ciliated Protozoa*†

SYNOPSIS. Three apparently conflicting generalizations concerning the relationships between molecules and cell structure may be derived from studies on cellular patterning in the ciliates. (A) Cells with identical genes and molecular composition may have different hereditary patterns. (B) Genes, through their prescribed molecular derivatives, constrain the modes of pattern permutation and define the states of greatest stability. (C) Cells with identical hereditary patterns may have entirely different genes and molecular compositions. These principles may be reconciled through the recognition that they are characteristically applicable over different time intervals. Hereditary differences within a clone and without related molecular differences (principle A) may persist for hundreds of cell generations, but they are resolved eventually within a constant environment (principle B) as the configuration of minimal free energy is approached. On an even longer time scale, molecular substitutions have occurred for many or most components of the cell, but these have been constrained by selective pressures on an ancient design (principle C) that disallow substitutions affecting certain form-function relations which have been elevated to an adaptive peak.     

Alloreactive T cells can distinguish between the same human class II MHC products on different B cell lines.

Certain allele-specific alloreactive T cell clones do not recognize the products expressed by some B cell lines that, according to typing methods other than sequencing, carry the allelic molecules recognized by these clones. In order to characterize the naturally occurring sequence polymorphisms putatively responsible for the differential allorecognition of these class II molecules, we have determined the third and/or second exon nucleotide sequences of HLA-DRB1, -DRB3/4/5, -DQB1, and -DQA1 genes from 35 representative lymphoblastoid cell lines. In some cases, the lack of recognition correlates with the presence of single amino acid substitutions in either the second or third hypervariable region (HVR) of the first domain of these molecules. In other cases, the differentially allorecognized class II molecules have identical second and/or first domain amino acid sequences. These findings indicate that a) class II MHC-alloreactive T cell clones can distinguish between molecules with identical amino acid sequences expressed by B cell lines established from unrelated individuals; b) allorecognition of class II molecules is sensitive to naturally occurring single amino acid substitutions in either the second HVR of class II molecules, which is unavailable to interact with TCR residues, or the third HVR. Our results also suggest that 1) in different B cell lines, identical class II molecules may present different endogenous peptides, which may behave as histocompatibility Ag; 2) the peptide-binding specificity of a class II molecule may be affected by amino acid substitutions in its second HVR (Ag-binding site); and 3) human class II allorecognition may be restricted by epitopes contributed by residues of their third HVR.     

The DNA-binding protein HU from mesophilic and thermophilic bacilli: gene cloning, overproduction and purification.

The major histone-like bacterial protein (HU)-encoding genes (hup) from five different Bacilli have been cloned, sequenced and overexpressed in Escherichia coli. The five Bacilli selected are closely related, but have different optimum growth temperatures: greater than 70 degrees C for Bacillus caldolyticus and B. caldotenax; 60-65 degrees C for B. stearothermophilus (Bst); 37 degrees C for B. subtilis and 30 degrees C for B. globigii. The deduced amino acid (aa) sequences from the three thermophiles are identical. Those from the two mesophiles are also identical and differ from those of the thermophiles at eleven aa positions. The mesophilic proteins have an extra two aa at the C terminus. Cells harbouring plasmids containing the hup genes can produce HU. An efficient purification scheme using cation-exchange chromatography and fast protein liquid chromatography is presented. This gives approx. 30-40 mg of more than 95% pure Bst HU per litre of E. coli culture.     

Evolution of duplicated alpha-tubulin genes in ciliates

Ciliates provide a powerful system to analyze the evolution of duplicated alpha-tubulin genes in the context of single-celled organisms. Genealogical analyses of ciliate alpha-tubulin sequences reveal five apparently recent gene duplications. Comparisons of paralogs in different ciliates implicate differing patterns of substitutions (e.g., ratios of replacement/synonymous nucleotides and radical/conservative amino acids) following duplication. Most substitutions between paralogs in Euplotes crassus, Halteria grandinella and Paramecium tetraurelia are synonymous. In contrast, alpha-tubulin paralogs within Stylonychia lemnae and Chilodonella uncinata are evolving at significantly different rates and have higher ratios of both replacement substitutions to synonymous substitutions and radical amino acid changes to conservative amino acid changes. Moreover, the amino acid substitutions in C. uncinata and S. lemnae paralogs are limited to short stretches that correspond to functionally important regions of the alpha-tubulin protein. The topology of ciliate alpha-tubulin genealogies are inconsistent with taxonomy based on morphology and other molecular markers, which may be due to taxonomic sampling, gene conversion, unequal rates of evolution, or asymmetric patterns of gene duplication and loss.     

Unequal rates of Y chromosome gene divergence during speciation of the family Ursidae

Evolution of the bear family Ursidae is well investigated in terms of morphological, paleontological, and genetic features. However, several phylogenetic ambiguities occur within the subfamily Ursinae (the family Ursidae excluding the giant panda and spectacled bear), which may correlate with behavioral traits of female philopatry and male-biased dispersal which form the basis of the observed matriarchal population structure in these species. In the process of bear evolution, we investigate the premise that such behavioral traits may be reflected in patterns of variation among genes with different modes of inheritance: matrilineal mitochondrial DNA (mtDNA), patrilineal Y chromosome, biparentally inherited autosomes, and the X chromosome. In the present study, we sequenced 3 Y-linked genes (3,453 bp) and 4 X-linked genes (4,960 bp) and reanalyzed previously published sequences from autosome genes (2,347 bp) in ursid species to investigate differences in evolutionary rates associated with patterns of inheritance. The results describe topological incongruence between sex-linked genes and autosome genes and between nuclear DNA and mtDNA. In more ancestral branches within the bear phylogeny, Y-linked genes evolved faster than autosome and X-linked genes, consistent with expectations based on male-driven evolution. However, this pattern changes among branches leading to each species within the lineage of Ursinae whereby the evolutionary rates of Y-linked genes have fewer than expected substitutions. This inconsistency between more recent nodes of the bear phylogeny with more ancestral nodes may reflect the influences of sex-biased dispersal as well as molecular evolutionary characteristics of the Y chromosome, and stochastic events in species natural history, and phylogeography unique to ursine bears.     

Divergence in expression between duplicated genes in Arabidopsis

New genes may arise through tandem duplication, dispersed small-scale duplication, and polyploidy, and patterns of divergence between duplicated genes may vary among these classes. We have examined patterns of gene expression and coding sequence divergence between duplicated genes in Arabidopsis thaliana. Due to the simultaneous origin of polyploidy-derived gene pairs, we can compare covariation in the rates of expression divergence and sequence divergence within this group. Among tandem and dispersed duplicates, much of the divergence in expression profile appears to occur at or shortly after duplication. Contrary to findings from other eukaryotic systems, there is little relationship between expression divergence and synonymous substitutions, whereas there is a strong positive relationship between expression divergence and nonsynonymous substitutions. Because this pattern is pronounced among the polyploidy-derived pairs, we infer that the strength of purifying selection acting on protein sequence and expression pattern is correlated. The polyploidy-derived pairs are somewhat atypical in that they have broader expression patterns and are expressed at higher levels, suggesting differences among polyploidy- and nonpolyploidy-derived duplicates in the types of genes that revert to single copy. Finally, within many of the duplicated pairs, 1 gene is expressed at a higher level across all assayed conditions, which suggests that the subfunctionalization model for duplicate gene preservation provides, at best, only a partial explanation for the patterns of expression divergence between duplicated genes.     

CG dinucleotide transitions in the factor IX gene account for about half of the point mutations in hemophilia B patients: a Seattle series

Summary Hemophilia B is due to multiple molecular defects in the factor IX gene. Over 80% of mutations are single base substitutions. By amplification and direct sequencing, 51 single base substitutions were found in the transcribed sequence of the factor IX genes of patients from 50 distinct families with hemophilia B. These include 30 mutations in 29 families not previously reported by us; of these, 12 are novel, i.e., not previously published in other series. Of the 51 substitutions in our overall series 23 (45%) occurred as C-to-T or G-to-A transitions at 11 sites within CG dinucleotides. It is estimated that CG transitions occur from one to two orders of magnitude more frequently than mutations in nucleotides that are not within a CG pair. More than one family had identical defects for 6 of the CG mutations. At 4 of these sites, most patients had different haplotypes compatible with distinct mutations. Non-CG-type mutations occurred thoughout the coding regions with only one mutation in more than one family. The latter included 7 families with a 397 Ile-to-Thr defect that all share a rare haplotype, suggesting a common ancestor.     

Effects of amino acid substitutions outside an antigenic site on protein binding to monoclonal antibodies of predetermined specificity obtained by peptide immunization: Demonstration with region 15–22 (antigenic site 1) of myoglobin

Monoclonal antibodies (mAbs) of predetermined specificity were prepared by immunizing with a free (i.e., not conjugated to any carrier) synthetic peptide representing region 15–22 (site 1) of sperm whale myoglobin (SpMb). The cross-reactions of Mb variants with three mAbs were studied in order to determine whether such interactions are influenced by substitutions outsde the site. Finback whale Mb, which has no substitutions within region 15–22, showed lower cross-reactivity and relative binding affinity than the reference antigen, SpMb. Bottle-nose Atlantic dolphin myoglobin (BdMb) and badger myoglobin (BgMb), although they have identical substitutions within region 15–22 (Ala-15 to Gly and Val-21 to Leu), showed very different binding properties. The cross-reaction of BdMb was quite comparable to that of SpMb, while that of BgMb was much lower. Since the two proteins have identical structures in regions 15–22, the differences in their cross-reactivities are readily attributed to the effects of substitutions outside this region. Another pair of myoglobins, horse myoglobins (HsMb) and chicken myoglobin (ChMb), also have two identical substitutions (Ala-15 to Gly and Val-21 to Ile) within region 15–22, but possessed different cross-reactivity. The results indicate that the reaction of mAbs, whose specificity is precisely known and predetermined by the immunizing free peptide, can be markedly affected by substitutions outside the indicated binding region on the protein.     

Sequence features and evolutionary mechanisms in the chicken avidin gene family.

The chicken avidin gene family comprises the avidin gene (avd) and several homologous avidin-related genes (avrs). The sequences of the avr genes are nearly identical to each other but exhibit nonrandomly distributed, frequently nonsynonymous nucleotide substitutions compared to avd. In this study, we determined the genetic distances and the phylogeny of the avd and avr genes and found differences between different exons and introns. Our results suggest the involvement of biased gene conversion in the evolution of the genes. Furthermore, one of the genes was identified as a putative fusion gene. The occurrence of both gene conversion and recombination supports the models suggesting a common initiation mechanism for conversion and crossing-over. The existence of avidin-related proteins (AVRs) is currently unknown, but the putative AVRs are expected to bind biotin similarly to avidin. However, the observed sequence differences may affect the stability and glycosylation patterns of the putative AVR proteins.     

The polymorphism and structural homology of storage polypeptides (hordein) coded by the Hor-2 locus in barley (Hordeum vulgare L.)

Two-dimensional mapping (isoelectric focusing followed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis) of the polypeptide components of “B” hordein fractions from eight barley varieties of widely different ancestry has been carried out. The relative positions of 47 different polypeptides were mapped, there being between 8 and 16 present in any one variety. The individual polypeptides differed in their distribution patterns; some were present in a number of varieties, while others were restricted to one or two. They also differed in their relative contributions to the total hordein fraction, both within and between varieties. The structural homology of the major polypeptides was compared by cleavage at methionine residues with cyanogen bromide and separation of the peptides on gradient gels. The polypeptides were classified into three groups which gave cleavage patterns with either two (class I), four (class II), or five (class III) low molecular weight bands. Class III polypeptides were found in all eight varieties, but in seven of the varieties class I or class II polypeptides were also present. With one exception, polypeptides migrating in the same position in different varieties gave identical or almost identical patterns. The three classes of polypeptides showed different distributions on the two-dimensional gels. Classes II and III polypeptides had a similar range of isoelectric points (pH 6.5–8.0), but all of the class II polypeptides were of slightly lower molecular weight. Class I polypeptides had a wider range of pI and molecular weight; the most alkaline and the lowest molecular weight polypeptides were in this group. The hordein fractions from a number of other barley varieties were compared with that of Julia. All had major polypeptides which migrated with ones present in Julia, but they differed in the relative amounts of these and in the absence of some polypeptides and the presence of others. B hordein is coded for by a single locus which has been suggested to be a complex multigenic family derived by duplication and divergence of a single gene. The data reported here provide support for this hypothesis and suggest that both mutations in the duplicated genes and recombination within the locus may have contributed to the polymorphism of the polypeptides.     

Evaluation of the flanking nucleotide sequences of sarcomeric hypertrophic cardiomyopathy substitution mutations

Hypertrophic cardiomyopathy (HCM) is a familial myocardial disease with a prevalence of 1 in 500. More than 400 causative mutations have been identified in 13 sarcomeric and myofilament related genes, 350 of these are substitution mutations within eight sarcomeric genes. Within a population, examples of recurring identical disease causing mutations that appear to have arisen independently have been noted as well as those that appear to have been inherited from a common ancestor. The large number of novel HCM mutations could suggest a mechanism of increased mutability within the sarcomeric genes. The objective of this study was to evaluate the most commonly reported HCM genes, beta myosin heavy chain (MYH7), myosin binding protein C, troponin I, troponin T, cardiac regulatory myosin light chain, cardiac essential myosin light chain, alpha tropomyosin and cardiac alpha-actin for sequence patterns surrounding the substitution mutations that may suggest a mechanism of increased mutability. The mutations as well as the 10 flanking nucleotides were evaluated for frequency of di-, tri- and tetranucleotides containing the mutation as well as for the presence of certain tri- and tetranculeotide motifs. The most common substitutions were guanine (G) to adenine (A) and cytosine (C) to thymidine (T). The CG dinucleotide had a significantly higher relative mutability than any other dinucleotide (p<0.05). The relative mutability of each possible trinucleotide and tetranucleotide sequence containing the mutation was calculated; none were at a statistically higher frequency than the others. The large number of G to A and C to T mutations as well as the relative mutability of CG may suggest that deamination of methylated CpG is an important mechanism for mutation development in at least some of these cardiac genes.     

Cell type-specific enhancement of hepatitis C virus internal ribosome entry site-directed translation due to 5' nontranslated region substitutions selected during passage of virus in lymphoblastoid cells

Low-level replication of hepatitis C virus (HCV) in cultured lymphoblastoid cells inoculated with H77 serum inoculum led to the appearance of new virus variants containing identical substitutions at three sites within the viral 5' nontranslated RNA (5'NTR): G(107)-->A, C(204)-->A, and G(243)-->A (N. Nakajima, M. Hijikata, H. Yoshikura, and Y. K. Shimizu, J. Virol. 70:3325-3329, 1996). These results suggest that virus with this 5'NTR sequence may have a greater capacity for replication in such cells, possibly due to more efficient cap-independent translation, since these nucleotide substitutions reside within the viral internal ribosome entry site (IRES). To test this hypothesis, we examined the translation of dicistronic RNAs containing upstream and downstream reporter sequences (Renilla and firefly luciferases, respectively) separated by IRES sequences containing different combinations of these substitutions. The activity of the IRES was assessed by determining the relative firefly and Renilla luciferase activities expressed in transfected cells. Compared with the IRES present in the dominant H77 quasispecies, an IRES containing all three nucleotide substitutions had significantly greater translational activity in three of five human lymphoblastoid cell lines (Raji, Bjab, and Molt4 but not Jurkat or HPBMa10-2 cells). In contrast, these substitutions did not enhance IRES activity in cell lines derived from monocytes or granulocytes (HL-60, KG-1, or THP-1) or hepatocytes (Huh-7) or in cell-free translation assays carried out with rabbit reticulocyte lysates. Each of the three substitutions was required for maximally increased translational activity in the lymphoblastoid cells. The 2- to 2.5-fold increase in translation observed with the modified IRES sequence may facilitate the replication of HCV, possibly accounting for differences in quasispecies variants recovered from liver tissue and peripheral blood mononuclear cells of the same patient.     

"Silent" sites in Drosophila genes are not neutral: evidence of selection among synonymous codons

The patterns of synonymous codon usage in 91 Drosophila melanogaster genes have been examined. Codon usage varies strikingly among genes. This variation is associated with differences in G+C content at silent sites, but (unlike the situation in mammalian genes) these differences are not correlated with variation in intron base composition and so are not easily explicable in terms of mutational biases. Instead, those genes with high G+C content at silent sites, resulting from a strong "preference" for a particular subset of the codons that are mostly C- ending, appear to be the more highly expressed genes. This suggests that G+C content is reduced in sequences where selective constraints are weaker, as indeed seen in a pseudogene. These and other data discussed are consistent with the effects of translational selection among synonymous codons, as seen in unicellular organisms. The existence of selective constraints on silent substitutions, which may vary in strength among genes, has implications for the use of silent molecular clocks.      

Identification of novel epigenetically modified genes in human melanoma via promoter methylation gene profiling

The inactivation of tumor-related genes through the aberrant methylation of promoter CpG islands is thought to contribute to tumor initiation and progression. We therefore investigated promoter methylation events involved in cutaneous melanoma by screening 30 genes of interest for evidence of promoter hypermethylation, examining 20 melanoma cell lines and 40 freshly procured melanoma samples. Utilizing quantitative methylation-specific PCR, we identified five genes (SOCS1, SOCS2, RAR-beta 2, TNFSF10C, and TNFSF10D) with hypermethylation frequencies ranging from 50% to 80% in melanoma cell lines as well as freshly procured tissue samples. Eighteen genes (LOX, RASSF1A, WFDC1, TM, APC, TFPI2, TNFSF10A, CDKN2A, MGMT, TIMP3, ASC, TPM1, IRF8, CIITA-PIV, CDH1, SYK, HOXB13, and DAPK1) were methylated at lower frequencies (2-30%). Two genes (CDKN1B and PTEN), previously reported as methylated in melanoma, and five other genes (RECK, IRF7, PAWR, TNFSF10B, and Rb) were not methylated in the samples screened here. Daughter melanoma cell lines showed identical methylation patterns when compared with original samples from which they were derived, as did synchronous metastatic lesions from the same patient. We identified four genes (TNFSF10C, TNFSF10D, LOX, and TPM1) that have never before been identified as hypermethylated in melanoma, with an overall methylation frequency of 60, 80, 50, and 10%, respectively, hypothesizing that these genes may play an important role in melanoma progression.     

Naturally occurring human parainfluenza type 3 viruses exhibit divergence in amino acid sequence of their fusion protein neutralization epitopes and cleavage sites.

Many human parainfluenza type 3 virus (PIV3) strains isolated from children with respiratory illness are resistant to neutralization by monoclonal antibodies (MAbs) which recognize epitopes in antigenic site A or B of the fusion (F) protein of the prototype 1957 PIV3 strain. The F protein genes of seven PIV3 clinical isolates were sequenced to determine whether their neutralization-resistant phenotypes were associated with specific differences in amino acids which are recognized by neutralizing MAbs. Several clinical strains which were resistant to neutralization by site A or B MAbs had amino acid differences at residues 398 or 73, respectively. These specific changes undoubtedly account for the neutralization-resistant phenotype of these isolates, since identical substitutions at residues 398 or 73 in MAb-selected escape mutants confer resistance to neutralization by site A or B MAbs. The existence of identical changes in naturally occurring and MAb-selected neutralization-resistant PIV3 strains raises the possibility that antigenically different strains may arise by immune selection during replication in partially immune children. Three of the seven clinical strains examined had differences in their F protein cleavage site sequence. Whereas the prototype PIV3 strain has the cleavage site sequence Arg-Thr-Lys-Arg, one clinical isolate had the sequence Arg-Thr-Arg-Arg and two isolates had the sequence Arg-Thr-Glu-Arg. The different cleavage site sequences of these viruses did not affect their level of replication in either continuous simian or bovine kidney cell monolayers (in the presence or absence of exogenous trypsin or plasmin) or in the upper or lower respiratory tract of rhesus monkeys. We conclude that two nonconsecutive basic residues within the F protein cleavage site are sufficient for efficient replication of human PIV3 in primates.     

Effects of amino acid substitutions outside an antigenic site on protein binding to monoclonal antibodies of predetermined specificity obtained by peptide immunization: demonstration with region 15-22 (antigenic site 1) of myoglobin.

Monoclonal antibodies (mAbs) of predetermined specificity were prepared by immunizing with a free (i.e., not conjugated to any carrier) synthetic peptide representing region 15–22 (site 1) of sperm whale myoglobin (SpMb). The cross-reactions of Mb variants with three mAbs were studied in order to determine whether such interactions are influenced by substitutions outsde the site. Finback whale Mb, which has no substitutions within region 15–22, showed lower cross-reactivity and relative binding affinity than the reference antigen, SpMb. Bottle-nose Atlantic dolphin myoglobin (BdMb) and badger myoglobin (BgMb), although they have identical substitutions within region 15–22 (Ala-15 to Gly and Val-21 to Leu), showed very different binding properties. The cross-reaction of BdMb was quite comparable to that of SpMb, while that of BgMb was much lower. Since the two proteins have identical structures in regions 15–22, the differences in their cross-reactivities are readily attributed to the effects of substitutions outside this region. Another pair of myoglobins, horse myoglobins (HsMb) and chicken myoglobin (ChMb), also have two identical substitutions (Ala-15 to Gly and Val-21 to Ile) within region 15–22, but possessed different cross-reactivity. The results indicate that the reaction of mAbs, whose specificity is precisely known and predetermined by the immunizing free peptide, can be markedly affected by substitutions outside the indicated binding region on the protein.     

CYP6B1 and CYP6B3 of the black swallowtail (Papilio polyxenes): adaptive evolution through subfunctionalization

Gene duplication provides essential material for functional divergence of proteins and hence allows organisms to adapt to changing environments. Following duplication events, redundant paralogs may undergo different evolutionary paths via processes known as nonfunctionalization, neofunctionalization, or subfunctionalization. Studies of adaptive evolution at the molecular level have progressed rapidly by computationally analyzing nucleotide substitution patterns but such studies are limited by the absence of information relating to alterations of function of the encoded enzymes. In this respect, evolution of the Papilio polyxenes cytochrome P450 monooxygenases (P450s) responsible for the adaptation of this insect to furanocoumarin-containing host plants provides an excellent model for elucidating the evolutionary fate of duplicated genes. Evidence from sequence and functional analysis in combination with molecular modeling indicates that the paralogous CYP6B1 and CYP6B3 genes in P. polyxenes have probably evolved via subfunctionalization after the duplication event by which they arose. Both enzymes have been under independent purifying selection as evidenced by the low dN/dS ratio in both the coding region and substrate recognition sites. Both enzymes have maintained their ability to metabolize linear and angular furanocoumarins albeit at different efficiencies. Comparisons of molecular models developed for the CYP6B3 and CYP6B1 proteins highlight differences in their binding modes that account for their different activities toward linear and angular furanocoumarins. That P. polyxenes maintains these 2 furanocoumarin-metabolizing loci with somewhat different activities and expression patterns provides this species with the potential to acquire P450s with novel functions while maintaining those most critical to its exclusive feeding on its current range of host plants.