A Bayesian network approach to feature selection in mass spectrometry data.

BACKGROUND: Time-of-flight mass spectrometry (TOF-MS) has the potential to provide non-invasive, high-throughput screening for cancers and other serious diseases via detection of protein biomarkers in blood or other accessible biologic samples. Unfortunately, this potential has largely been unrealized to date due to the high variability of measurements, uncertainties in the distribution of proteins in a given population, and the difficulty of extracting repeatable diagnostic markers using current statistical tools. With studies consisting of perhaps only dozens of samples, and possibly hundreds of variables, overfitting is a serious complication. To overcome these difficulties, we have developed a Bayesian inductive method which uses model-independent methods of discovering relationships between spectral features. This method appears to efficiently discover network models which not only identify connections between the disease and key features, but also organizes relationships between features--and furthermore creates a stable classifier that categorizes new data at predicted error rates. RESULTS: The method was applied to artificial data with known feature relationships and typical TOF-MS variability introduced, and was able to recover those relationships nearly perfectly. It was also applied to blood sera data from a 2004 leukemia study, and showed high stability of selected features under cross-validation. Verification of results using withheld data showed excellent predictive power. The method showed improvement over traditional techniques, and naturally incorporated measurement uncertainties. The relationships discovered between features allowed preliminary identification of a protein biomarker which was consistent with other cancer studies and later verified experimentally. CONCLUSIONS: This method appears to avoid overfitting in biologic data and produce stable feature sets in a network model. The network structure provides additional information about the relationships among features that is useful to guide further biochemical analysis. In addition, when used to classify new data, these feature sets are far more consistent than those produced by many traditional techniques.

Molecular pathology of prostate cancer.

This chapter includes discussion of the molecular pathology of tissue, blood, urine, and expressed prostatic secretions. Because we are unable to reliably image the disease in vivo, a 12 core method that oversamples the peripheral zone is widely used. This generates large numbers of cores that need to be carefully processed and sampled. In spite of the large number of tissue cores, the amount of tumor available for study is often quite limited. This is a particular challenge for research, as new biomarker assays will need to preserve tissue architecture intact for histopathology. Methods of processing and reporting pathology are discussed. With the exception of ductal variants, recognized subtypes of prostate cancer are largely confined to research applications, and most prostate cancers are acinar. Biomarker discovery in urine and expressed prostatic secretions would be useful since these are readily obtained and are proximate fluids. The well-known challenges of biomarker discovery in blood and urine are referenced and discussed. Mediators of carcinogenesis can serve as biomarkers as exemplified by mutations in PTEN and TMPRSS2:ERG fusion. The use of proteomics in biomarker discovery with an emphasis on imaging mass spectroscopy of tissues is discussed. Small RNAs are of great interest, however, their usefulness as biomarkers in clinical decision making remains the subject of ongoing research. The chapter concludes with an overview of blood biomarkers such as circulating nucleic acids and tumor cells and bound/free isoforms of prostate specific antigen (PSA).

Discrete serum protein signatures discriminate between human retrovirus-associated hematologic and neurologic disease.

The human T-cell leukemia virus type I (HTLV-I) is the causative agent for adult T-cell leukemia (ATL) and HTLV-I-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Approximately 5% of infected individuals will develop either disease and currently there are no diagnostic tools for early detection or accurate assessment of disease state. We have employed high-throughput expression profiling of serum proteins using mass spectrometry to identify protein expression patterns that can discern between disease states of HTLV-I-infected individuals. Our study group consisted of 42 ATL, 50 HAM/TSP, and 38 normal controls. Spectral peaks corresponding to peptide ions were generated from MS-TOF data. We applied Classification and Regression Tree analysis to build a decision algorithm, which achieved 77% correct classification rate across the three groups. A second cohort of 10 ATL, 10 HAM and 10 control samples was used to validate this result. Linear discriminate analysis was performed to verify and visualize class separation. Affinity and sizing chromatography coupled with tandem mass spectrometry was used to identify three peaks specifically overexpressed in ATL: an 11.7 kDa fragment of alpha trypsin inhibitor, and two contiguous fragments (19.9 and 11.9 kDa) of haproglobin-2. To the best of our knowledge, this is the first application of protein profiling to distinguish between two disease states resulting from a single infectious agent.

Proteomic approaches to biomarker discovery in prostate and bladder cancers.

Proteomic technologies, including high resolution two-dimensional electrophoresis (2-DE), antibody/protein arrays, and advances in mass spectrometry (MS), are providing the tools needed to discover and identify disease associated biomarkers. Although application of these technologies to search for potential diagnostic/prognostic biomarkers associated with prostate and bladder cancer have been somewhat limited to date, proteins either overexpressed or underexpressed have been detected in both these urological cancers. Recent advances in mass spectrometry, especially platforms that permit rapid "fingerprint" profiling of multiple biomarkers, and tandem mass spectrometers for protein identification, will most assuredly enhance the discovery, identification, and characterization of potential cancer associated biomarkers. Furthermore, application of laser capture microdissection microscopes has provided a rapid and reproducible approach to procure pure populations of cells. This technology coupled to 2-DE and MS has significantly aided the elucidation of the differential expression profiles between disease, benign and normal prostate and bladder cell populations. Finally, development and application of learning algorithms and bioinformatics to the data generated by these proteomic technologies will be essential in determining the clinical potential of a protein biomarker. The purpose of this review is to provide the reader with an overview of the application of these technologies in the search and identification of potential diagnostic/prognostic biomarkers for prostate and bladder cancers.

Properties of two EBV Mta nuclear export signal sequences.

The Epstein-Barr virus (EBV) Mta protein is a posttranscriptional regulator of EBV lytic gene expression that affects RNA splicing and transport. Mta mediates cytoplasmic accumulation of unspliced EBV replication gene transcripts and shuttles between the nucleus and cytoplasm. Mta contains a recognized leucine-rich, putative nuclear export signal (NES) between aa 227 and 236. Deletion of this signal sequence eliminated shuttling, while mutation of the core LXL motif in the putative NES diminished but did not abolish the ability of Mta to shuttle from donor to recipient cells in a heterokaryon assay. A double mutation of the LXL motif plus an upstream VTL motif eliminated shuttling, suggesting that Mta may have two NES motifs. In confirmation of this, transfer of either the sequence encoding the leucine-rich aa 227-236 motif or that encoding the adjacent hydrophobic aa 218-227 sequence to a GFP-NLS-pyruvate kinase reporter protein conferred the property of cytoplasmic accumulation onto the heterologous protein. Cytoplasmic accumulation of both the aa 225-237 and 218-227 containing reporters was minimal in the presence of the inhibitor leptomycin B, indicating that both motifs mediated Crm-1-dependent export. Mutations in the NES signal sequences abolished the ability of Mta to mediate cytoplasmic accumulation of BALF2 replication gene transcripts. This included mutation of the LXL motif which still showed cytoplasmic shuttling, suggesting that the NES mutations might have additional effects on Mta function. Wild-type Mta co-immunoprecipitated with the splicing factor SC35 and colocalized with SC35 in transfected cells, modifying endogenous SC35 distribution within the nucleus to give more intense, rounded spots. Interestingly, the NES mutant proteins appeared to have altered interactions with the splicing complex, binding more tightly to SC35 in co-immunoprecipitation assays. These observations suggest a linkage between the splicing and export functions of Mta.

Human T-cell leukemia virus type 1 Tax shuttles between functionally discrete subcellular targets.

Human T-cell leukemia virus type 1 (HTLV-1) Tax is a nuclear protein with striking pleiotropic functionality. We recently demonstrated that Tax localizes to a multicomponent nuclear speckled structure (Tax speckled structure [TSS]). Here, we examine these structures further and identify a partial overlap of TSS with transcription hot spots. We used a strategy of directed expression via fusion proteins to determine if these transcription sites are the subtargets within TSS required for Tax function. When fused to human immunodeficiency virus type 1 (HIV-1) Tat, the resulting Tat-Tax fusion protein displayed neither a Tat-like nor a Tax-like pattern but rather was targeted specifically to the transcription subsites. The Tat-Tax fusion was able to activate both the HIV-1 long terminal repeat (LTR) and the HTVL-1 LTR at the same level as the individual component; thus, targeting proteins to transcription hot spots was compatible with both Tax and Tat transcription function. In contrast, the fusion with HIV-1 Rev, Rev-Tax, resulted in a pattern of expression that was largely Rev-like (nucleolar and cytoplasmic). The reduced localization of Rev-Tax to transcription sites was reflected in a 10-fold drop in activation of the HTLV-1 LTR. However, there was no loss in the ability of Tax to activate via NF-kappaB. Thus, NF-kappaB-dependent Tax function does not require targeting of Tax to these transcription sites and suggests that activation via NF-kappaB is a cytoplasmic function. Selective mutation of the nuclear localization signal site in the Rev portion resulted in retargeting of Rev-Tax to TSS and subsequent restoration of transcription function, demonstrating that inappropriate localization preceded loss of function. Mutation of the nuclear export signal site in the Rev portion had no effect on transcription, although the relative amount of Rev-Tax in the cytoplasm was reduced. Finally, in explaining how Tax can occupy multiple subcellular sites, we show that Tax shuttles from the nucleus to the cytoplasm in a heterokaryon fusion assay. Thus, pleiotropic functionality by Tax is regulatable via shuttling between discrete subcellular compartments.

Analysis of potential phosphorylation sites in human T cell leukemia virus type 1 Tax.

The human T cell leukemia virus type 1 (HTLV-1) Tax is a phosphoprotein, however, the contribution of phosphorylation to Tax activity is unknown. Previous studies have shown that phosphorylation of Tax occurs on serine residue(s), within one tryptic fragment, in response to 4beta-phorbol-12beta-myristate-13alpha-acetate, in both mouse and human cells. Studies were conducted in multiple cell lines to identify the specific phosphorylated serines as a prelude to functional analysis. The phosphorylation pattern of Tax was found to be different in 293T and COS-7 cells in comparison with MT-4 and Px-1 cells. However, one tryptic fragment remained consistent in comigration analyses among all cell lines. Using selected Tax serine mutants a tryptic fragment containing a serine at residue 113 believed to be the site of phosphorylation of Tax did not comigrate with the common phosphorylated tryptic fragment. Analysis of selected Tax mutants for ability to trans-activate the cytomegalovirus promoter demonstrated mutation of serine 77 to alanine reduced trans-activation by 90% compared to wild-type Tax. However, examination of the phosphorylation pattern of the serine 77 mutant demonstrated that it is not the site of phosphorylation. These studies demonstrate the importance of using relevant cell lines to characterize the role of phosphorylation in protein function.

Mta has properties of an RNA export protein and increases cytoplasmic accumulation of Epstein-Barr virus replication gene mRNA.

The Epstein-Barr virus (EBV) Zta and Mta regulatory proteins were previously found to be required for efficient replication of oriLyt in cotransfection-replication assays, but the contribution of Mta to the replication process was unknown. We now demonstrate that Mta regulates replication gene expression. Using the polymerase processivity factor BMRF1 as an example, we found that in transfected cells, total BMRF1 mRNA levels were unaffected by Mta but that the amounts of cytoplasmic BMRF1 RNA and protein were greatly increased in the presence of Mta. Mta also increased cytoplasmic accumulation of the BALF2, BALF5, BSLF1, and BBLF4 replication gene mRNAs but did not affect cytoplasmic levels of BBLF2/3 mRNA. Thus, five of the six core replication genes require Mta for efficient accumulation of cytoplasmic RNA. The contribution of Mta to posttranscriptional RNA processing was examined. Examination of Mta localization in transfected cells by indirect immunofluorescence revealed that Mta colocalized with the splicing factor SC35. We also found that Mta has RNA binding activity. Glutathione S-transferase-Mta bound to BMRF1 and BMLF1 transcripts but not to a control cellular gene RNA. Mta contains a consensus leucine-rich nuclear export signal. Such signal sequences are characteristic of proteins that undergo nuclear export. Examination of Mta localization in a heterokaryon assay provided evidence that Mta shuttles between the nucleus and the cytoplasm. Our experiments indicate that Mta functions in RNA processing and transport and mediates cytoplasmic accumulation of a number of EBV early mRNAs.

The Epstein-Barr virus lytic transactivator Zta interacts with the helicase-primase replication proteins.

The Epstein-Barr virus transactivator Zta triggers lytic gene expression and is essential for replication of the lytic origin, oriLyt. Previous analysis indicated that the Zta activation domain contributed a replication-specific function. We now show that the Zta activation domain interacts with components of the EBV helicase-primase complex. The three helicase-primase proteins BBLF4 (helicase), BSLF1 (primase), and BBLF2/3 (primase-associated factor) were expressed fused to the Myc epitope. When expression plasmids for BBLF4 or BBLF2/3 plus BSLF1 (primase subcomplex) were separately transfected, the proteins localized to the cytoplasm. Interaction between Zta and the components of the helicase-primase complex was tested by examining the ability of Zta to alter the intracellular localization of these proteins. Cotransfection of Zta with Myc-BBLF4 resulted in nuclear translocation of Myc-BBLF4; similarly, cotransfection of Zta with the primase subcomplex led to nuclear translocation of the Myc-BSLF1 and Myc-BBLF2/3 proteins. This relocalization provides evidence for an interaction between Zta and the helicase and Zta and the primase subcomplex. An affinity assay using glutathione S-transferase-Zta fusion proteins demonstrated that Myc-BBLF4 and Myc-BBLF2/3 plus BSLF1 bound to the Zta activation domain (amino acids 1 to 133). In the nuclear relocalization assay, the amino-terminal 25 amino acids of Zta were required for efficient interaction with the primase subcomplex but not for interaction with BBLF4. Evidence for interaction between oriLyt bound Zta and the helicase-primase complex was obtained in a superactivation assay using an oriLyt-chloramphenicol acetyltransferase (CAT) reporter. Zta activated expression from a CAT reporter containing the complete oriLyt region and regulated by the oriLyt BHLF1 promoter. Cotransfection of the helicase-primase proteins, one of which was fused to a heterologous activation domain, led to Zta-dependent superactivation of CAT expression. This assay also provided evidence for an interaction between the single-stranded DNA binding protein, BALF2, and the Zta-tethered helicase-primase complex. The helicase-primase interaction is consistent with a role for Zta in stabilizing the formation of an origin-bound replication complex.

A human T-cell leukemia virus Tax variant incapable of activating NF-kappaB retains its immortalizing potential for primary T-lymphocytes.

The human T-cell leukemia virus type 1 (HTLV-1) transactivator (Tax) has been shown to interfere with regulated cellular proliferation. Many studies have focused on the ability of Tax to transform rodent fibroblasts; however, none has defined the molecular requirements for Tax transformation of human lymphoid cells. We show here that tax induces permanent growth of human primary T-lymphocytes by using a transformation/immortalization defective rhadinovirus vector. The cells phenotypically resemble HTLV-immortalized lymphocytes and contain episomally persisting recombinant rhadinoviral sequences, which stably express functional Tax protein. As Tax can activate major cellular signal transducing pathways including NF-kappaB and cAMP-responsive element binding protein (CREB), we asked for the relevance of these routes in the immortalization of human lymphocytes. By using Tax mutants that either activate exclusively CREB/activating transcription factor or are defective in activating this signaling pathway, we delineated that Tax can induce immortalization of primary human T-lymphocytes through a mechanism independent of NF-kappaB activation.

Modulation of Sp1 phosphorylation by human immunodeficiency virus type 1 Tat.

We previously reported (K. T. Jeang, R. Chun, N. H. Lin, A. Gatignol, C. G. Glabe, and H. Fan, J. Virol. 67: 6224-6233, 1993) that human immunodeficiency virus type 1 (HIV-1) Tat and Sp1 form a protein-protein complex. Here, we have characterized the physical interaction and a functional consequence of Tat-Sp1 contact. Using in vitro protein chromatography, we mapped the region in Tat that contacts Sp1 to amino acids 30 to 55. We found that in cell-free reactions, Tat augmented double-stranded DNA-dependent protein kinase (DNA-PK)-mediated Sp1 phosphorylation in a contact-dependent manner. Tat mutants that do not bind Sp1 failed to influence phosphorylation of the latter. In complementary experiments, we also found that Tat forms protein-protein contacts with DNA-PK. We confirmed that in HeLa and Jurkat cells, Tat expression indeed increased the intracellular amount of phosphorylated Sp1 in a manner consistent with the results of cell-free assays. Furthermore, using two phosphatase inhibitors and a kinase inhibitor, we demonstrated a modulation of reporter gene expression as a consequence of changes in Sp1 phosphorylation. Taken together, these findings suggest that activity at the HIV-1 promoter is influenced by phosphorylation of Sp1 which is affected by Tat and DNA-PK.

Divergent Subcellular Locations of HTLV-I Tax and Int-6: A Contrast between in vitro Protein-Protein Binding and Intracellular Protein Colocalization.

Protein-protein interactions define many important molecular and cellular processes in prokaryotic and eukaryotic biology. In trying to delineate the contact between two proteins, the yeast two-hybrid assay has emerged as a powerful technique. Complementing the yeast two-hybrid assay are in vitro techniques (e.g. GST-fusion-protein chromatography) that can also yield information on protein-protein associations. However, unambiguous functional significance to these interactions is best supported through a finding of colocalization of two proteins inside cells. In instances where two proteins interact in vitro but have divergent localizations within cells one needs to reconsider the biological importance of the former finding. Here, we present evidence for different subcellular locations of HTLV-I Tax and the Int-6 protein. We suggest a reexploration of the functional significance between Tax and Int-6 in cellular transformation.

HTLV-I Tax and Cytokeratin: Tax-Expressing Cells Show Morphological Changes in Keratin-Containing Cytoskeletal Networks.

Human T cell leukemia virus type I (HTLV-I) has been linked to the development of an aggressive lymphoproliferative disorder (adult T cell leukemia), a chronic neurodegenerative presentation (HTLV-I-associated myelopathy/tropical spastic paraparesis) and numerous less well-defined inflammatory conditions. The viral regulatory protein Tax has been implicated in cellular transformation events leading to the onset of adult T cell leukemia. Details on the stepwise processes through which Tax induces morphological changes in cells are poorly understood. We show here that Tax can bind to a class of intermediate filaments, the cytokeratins (Ker). Tax interacts with the 1B helical coil of keratin 8, a domain critical for higher-order intermediate filament matrix formation. Expression of Tax in epithelial cells visibly altered the structural pattern of the Ker network. In a T lymphocyte cell line, induction of Tax expression resulted in increased cellular adherence/invasion of Matrigel filters. We propose that one aspect of Tax function is the induction of morphological changes in cellular cytoskeletal structures. This finding for Tax-expressing cells might be one factor contributing directly to the pathogenesis of HTLV-I disease(s). Copyright 1997 S. Karger AG, Basel

Localization of human T-cell leukemia virus type 1 tax to subnuclear compartments that overlap with interchromatin speckles.

Tax, the virally encoded activator of the human T-cell leukemia virus type 1 long terminal repeats, regulates the expression of many cellular genes. This protein has been implicated in transformation events leading to the development of adult T-cell leukemia. Because subcellular localization contributes importantly to protein function, we determined the compartment(s) within the cell in which Tax is found. Using confocal microscopy, we found that Tax localizes to subnuclear domains which overlap with structures previously identified as interchromatin granules or spliceosomal speckles. These Tax speckled structures are coincident with a subset of nuclear transcriptional hot spots. Disruption of the Tax speckled structures by heat shock revealed the existence of different populations of Tax. One population of Tax is tightly associated with nuclear speckles. A second population exists outside of the speckles and is transcriptionally active for some promoters.

Human T-cell leukemia virus type I tax masks c-Myc function through a cAMP-dependent pathway.

Human T-cell leukemia virus type I Tax is a pleiotropic gene regulator that functions through CREB/ATF- and NF-kappaB-mediated pathways. In most contexts, Tax is a potent gene activator. Here, we describe an unexpected finding of Myc repression by Tax. In cells that overexpress human T-cell leukemia virus type I Tax, the detection of c-Myc protein in the nucleus by a monoclonal antibody was masked. Tax prevented immunological visualization of a Myc epitope contained within amino acids 45-104, resulting in interference with Myc function in transcription and in anchorage-independent cell growth. Tax did not affect steady-state protein levels since detection of c-Myc with other antibodies was unperturbed. Four observations suggest that this Tax-Myc interaction is mediated through CREB/ATF signal transduction. 1) Tax point mutants, selectively defective for activation of CREB/ATF but not NF-kappaB, failed to mask c-Myc; 2) masking of Myc was abolished when Tax-expressing cells were treated with protein kinase inhibitor H-9; 3) Tax-specific shielding of Myc is absent in cells (B1R) that are genetically defective for cAMP signaling; and 4) forskolin treatment of cells mimicked Tax in masking the Myc epitope. Considered collectively, these findings suggest a regulation of Myc function at the level of localized protein conformation.

HTLV-I and HTLV-II Tax: differences in induction of micronuclei in cells and transcriptional activation of viral LTRs.

Human T-cell leukemia virus (HTLV) types I and II are highly related viruses that differ in disease manifestations. HTLV-I has been linked unmistakably to adult T-cell leukemia-lymphoma. On the other hand, there is little evidence that prior infection with HTLV-II increases risk for lymphoproliferative disorders. Both viruses encode homologous transcriptional-activating proteins (respectively designated as Tax1 and Tax2) which have been suggested to be important mediators of viral pathogenesis. Previously, we reported that Tax1 is a potent inducer of micronuclei formation in cells. Here, we present evidence that Tax2 lacks micronuclei inductive ability. We contrast this phenotypic difference between Tax1 and Tax2 at the cellular level with their similarities at the molecular level in transcriptional activation.

The amino terminus of Tax is required for interaction with the cyclic AMP response element binding protein.

Tax of human T-cell leukemia virus type 1 was analyzed for interaction with the cyclic AMP response element binding protein (CREB) in vitro with and without Tax response element DNA. Mutations in the carboxy terminus of Tax (L296G and L320G) did not affect binding to CREB and led to supershifts. In contrast, mutants with changes in the amino-terminal cysteine-rich region lost the ability to bind to CREB. The S10A mutant protein bound moderately. Thus, the amino terminus of Tax is essential for Tax-CREB interaction.

Definition of a minimal activation domain in human T-cell leukemia virus type I Tax.

Fourteen mutants were used to delineate a minimal activation domain in the Tax protein of human T-cell leukemia virus type I. In an assay using a Gal4-Tax (GalTx) fusion protein and a responsive promoter containing Gal4 consensus binding sites, we found that activation was "squelched" by coexpression of wild-type Tax protein in trans. When Tax mutants were tested for squelching, many competed effectively against GalTx. However, those containing changes in amino acids 289 to 322 failed to inhibit activity. In particular, three mutants that were expressed stably, with changes at amino acids 289, 296, and 320 respectively, did not squelch GalTx activity. On the other hand, mutants with individual changes at amino acid 3, 9, 29, 41, 273, and 337 efficiently inhibited GalTx function. Three other mutants failed to be stably expressed. In separate experiments, when fused alone to the DNA-binding domain of Gal4, amino acids 289 to 322 of Tax conferred trans activation ability. This fusion protein was able to activate a core promoter. These findings suggest that amino acids 289 to 322 define a region that contacts an essential transcription factor and that this region is a modular activation domain.

Induction of micronuclei by HTLV-I Tax: a cellular assay for function.

Cellular chromosomal damage is ubiquitously seen in HTLV-I-transformed lymphocytes. It is also characteristic of cells that have been exposed to mutagens. A sensitive measurement for mutagen-induced DNA damage is the formation of micronuclei in treated cells. Because current evidence suggests that HTLV-I Tax is etiologically linked to transformation, we tested for its activity in inducing micronuclei. We show here that transfection into cells of a Tax-producing plasmid rapidly induced the formation of micronuclei. This effect cooperated with that of a mutagen (mitomycin C) and was correlated with the inherent trans-activation capacity of Tax. These findings suggest that a commonly used mutagen assay could be a quick biological test for putatively oncogenic proteins.

Mutational analysis of human T-cell leukemia virus type I Tax: regions necessary for function determined with 47 mutant proteins.

We have made 47 mutations that span the length of the human T-cell leukemia virus type I (HTLV-I) Tax open reading frame. Of the 47 mutations, 38 were substitutions of single amino acids, 5 were missense changes in two or more amino acids, and 4 were deletions. A subset of these mutations includes individual changes of all 26 naturally occurring serines to alanines. By assaying each mutant protein separately on the HTLV-I long terminal repeat (LTR) and the human immunodeficiency virus type 1 (HIV-1) LTR in parallel, we were able to identify regions of Tax selectively necessary for each promoter. A small region in the carboxyl terminus, amino acids 315 to 325, was found to be selectively important for activation of the HTLV-I LTR. Three changes at serine 113, serine 160, and serine 258 were found to specifically affect function on the HIV-1 LTR. Surprisingly, we found that the great preponderance of missense changes (32 of 42) in Tax did not affect function.