WormBase: a modern Model Organism Information Resource.

WormBase (https://wormbase.org/) is a mature Model Organism Information Resource supporting researchers using the nematode Caenorhabditis elegans as a model system for studies across a broad range of basic biological processes. Toward this mission, WormBase efforts are arranged in three primary facets: curation, user interface and architecture. In this update, we describe progress in each of these three areas. In particular, we discuss the status of literature curation and recently added data, detail new features of the web interface and options for users wishing to conduct data mining workflows, and discuss our efforts to build a robust and scalable architecture by leveraging commercial cloud offerings. We conclude with a description of WormBase's role as a founding member of the nascent Alliance of Genome Resources.

WormBase 2017: molting into a new stage.

WormBase (http://www.wormbase.org) is an important knowledge resource for biomedical researchers worldwide. To accommodate the ever increasing amount and complexity of research data, WormBase continues to advance its practices on data acquisition, curation and retrieval to most effectively deliver comprehensive knowledge about Caenorhabditis elegans, and genomic information about other nematodes and parasitic flatworms. Recent notable enhancements include user-directed submission of data, such as micropublication; genomic data curation and presentation, including additional genomes and JBrowse, respectively; new query tools, such as SimpleMine, Gene Enrichment Analysis; new data displays, such as the Person Lineage browser and the Summary of Ontology-based Annotations. Anticipating more rapid data growth ahead, WormBase continues the process of migrating to a cutting-edge database technology to achieve better stability, scalability, reproducibility and a faster response time. To better serve the broader research community, WormBase, with five other Model Organism Databases and The Gene Ontology project, have begun to collaborate formally as the Alliance of Genome Resources.

Correction to: Use of multicolor fluorescence in situ hybridization to detect deletions in clinical tissue sections.

Figure 2 is incorrect in the original version of this article. The correct figure 2 is provided below.

Use of multicolor fluorescence in situ hybridization to detect deletions in clinical tissue sections.

A variety of laboratory methods are available for the detection of deletions of tumor suppressor genes and losses of their proteins. The clinical utility of fluorescence in situ hybridization (FISH) for the identification of deletions of tumor suppressor genes has previously been limited by difficulties in the interpretation of FISH signal patterns. The first deletion FISH assays using formalin-fixed paraffin-embedded tissue sections had to deal with a significant background level of signal losses affecting nuclei that are truncated by the cutting process of slide preparation. Recently, more efficient probe designs, incorporating probes adjacent to the tumor suppressor gene of interest, have increased the accuracy of FISH deletion assays so that true chromosomal deletions can be readily distinguished from the false signal losses caused by sectioning artifacts. This mini-review discusses the importance of recurrent tumor suppressor gene deletions in human cancer and reviews the common FISH methods being used to detect the genomic losses encountered in clinical specimens. The use of new probe designs to recognize truncation artifacts is illustrated with a four-color PTEN FISH set optimized for prostate cancer tissue sections. Data are presented to show that when section thickness is reduced, the frequency of signal truncation losses is increased. We also provide some general guidelines that will help pathologists and cytogeneticists run routine deletion FISH assays and recognize sectioning artifacts. Finally, we summarize how recently developed sequence-based approaches are being used to identify recurrent deletions using small DNA samples from tumors.

WormBase 2016: expanding to enable helminth genomic research.

WormBase (www.wormbase.org) is a central repository for research data on the biology, genetics and genomics of Caenorhabditis elegans and other nematodes. The project has evolved from its original remit to collect and integrate all data for a single species, and now extends to numerous nematodes, ranging from evolutionary comparators of C. elegans to parasitic species that threaten plant, animal and human health. Research activity using C. elegans as a model system is as vibrant as ever, and we have created new tools for community curation in response to the ever-increasing volume and complexity of data. To better allow users to navigate their way through these data, we have made a number of improvements to our main website, including new tools for browsing genomic features and ontology annotations. Finally, we have developed a new portal for parasitic worm genomes. WormBase ParaSite (parasite.wormbase.org) contains all publicly available nematode and platyhelminth annotated genome sequences, and is designed specifically to support helminth genomic research.

Exploring high dimensional data with Butterfly: a novel classification algorithm based on discrete dynamical systems.

MOTIVATION: We introduce a novel method for visualizing high dimensional data via a discrete dynamical system. This method provides a 2D representation of the relationship between subjects according to a set of variables without geometric projections, transformed axes or principal components. The algorithm exploits a memory-type mechanism inherent in a certain class of discrete dynamical systems collectively referred to as the chaos game that are closely related to iterative function systems. The goal of the algorithm was to create a human readable representation of high dimensional patient data that was capable of detecting unrevealed subclusters of patients from within anticipated classifications. This provides a mechanism to further pursue a more personalized exploration of pathology when used with medical data. For clustering and classification protocols, the dynamical system portion of the algorithm is designed to come after some feature selection filter and before some model evaluation (e.g. clustering accuracy) protocol. In the version given here, a univariate features selection step is performed (in practice more complex feature selection methods are used), a discrete dynamical system is driven by this reduced set of variables (which results in a set of 2D cluster models), these models are evaluated for their accuracy (according to a user-defined binary classification) and finally a visual representation of the top classification models are returned. Thus, in addition to the visualization component, this methodology can be used for both supervised and unsupervised machine learning as the top performing models are returned in the protocol we describe here. RESULTS: Butterfly, the algorithm we introduce and provide working code for, uses a discrete dynamical system to classify high dimensional data and provide a 2D representation of the relationship between subjects. We report results on three datasets (two in the article; one in the appendix) including a public lung cancer dataset that comes along with the included Butterfly R package. In the included R script, a univariate feature selection method is used for the dimension reduction step, but in the future we wish to use a more powerful multivariate feature reduction method based on neural networks (Kriesel, 2007). AVAILABILITY AND IMPLEMENTATION: A script written in R (designed to run on R studio) accompanies this article that implements this algorithm and is available at http://butterflygeraci.codeplex.com/. For details on the R package or for help installing the software refer to the accompanying document, Supporting Material and Appendix.

WormBase 2024: status and transitioning to Alliance infrastructure.

WormBase has been the major repository and knowledgebase of information about the genome and genetics of Caenorhabditis elegans and other nematodes of experimental interest for over 2 decades. We have 3 goals: to keep current with the fast-paced C. elegans research, to provide better integration with other resources, and to be sustainable. Here, we discuss the current state of WormBase as well as progress and plans for moving core WormBase infrastructure to the Alliance of Genome Resources (the Alliance). As an Alliance member, WormBase will continue to interact with the C. elegans community, develop new features as needed, and curate key information from the literature and large-scale projects.

MicroRNA signature helps distinguish early from late biochemical failure in prostate cancer.

PURPOSE: Prostate-specific antigen testing has led to overtreatment of prostate cancer (PCa). Only a small subset of PCa patients will have an aggressive disease that requires intensive therapy, and there is currently no biomarker to predict disease aggressiveness at the time of surgery. MicroRNAs (miRNAs) are reported to be involved in PCa pathogenesis. METHODS: This study involved 105 participants. For the discovery phase, prostatectomy samples were dichotomized to high-risk (n = 27, biochemical failure <36 months after prostatectomy) and low-risk groups (n = 14, ≥ 36 months without biochemical failure). Expression of 754 mature miRNAs was compared between the 2 groups. Linear regression models were built to accurately predict biochemical failure risk. miRNA mimics were transfected into PCa model cell lines to test effects on proliferation and to deduce responding signaling pathways. RESULTS: We identified 25 differentially expressed miRNAs between the biochemical failure risk groups. Based on the expression of 2-3 miRNAs, 3 logistic regression models were developed, each with a high positive predictive value. Candidate miRNAs and the best-performing model were also verified on an independent PCa set. miRNA-152, featured in the models, was further investigated by using cell line models and was shown to affect cell proliferation. Predicted interaction between miR-152 and (mRNA)ERBB3 (erythroblastic leukemia viral oncogene homolog 3) was experimentally validated in vitro. CONCLUSIONS: miRNAs can help to predict biochemical failure risk at the time of prostatectomy.

Identification of the IGF1/PI3K/NF κB/ERK gene signalling networks associated with chemotherapy resistance and treatment response in high-grade serous epithelial ovarian cancer.

BACKGROUND: Resistance to platinum-based chemotherapy remains a major impediment in the treatment of serous epithelial ovarian cancer. The objective of this study was to use gene expression profiling to delineate major deregulated pathways and biomarkers associated with the development of intrinsic chemotherapy resistance upon exposure to standard first-line therapy for ovarian cancer. METHODS: The study cohort comprised 28 patients divided into two groups based on their varying sensitivity to first-line chemotherapy using progression free survival (PFS) as a surrogate of response. All 28 patients had advanced stage, high-grade serous ovarian cancer, and were treated with standard platinum-based chemotherapy. Twelve patient tumours demonstrating relative resistance to platinum chemotherapy corresponding to shorter PFS (< eight months) were compared to sixteen tumours from platinum-sensitive patients (PFS > eighteen months). Whole transcriptome profiling was performed using an Affymetrix high-resolution microarray platform to permit global comparisons of gene expression profiles between tumours from the resistant group and the sensitive group. RESULTS: Microarray data analysis revealed a set of 204 discriminating genes possessing expression levels which could influence differential chemotherapy response between the two groups. Robust statistical testing was then performed which eliminated a dependence on the normalization algorithm employed, producing a restricted list of differentially regulated genes, and which found IGF1 to be the most strongly differentially expressed gene. Pathway analysis, based on the list of 204 genes, revealed enrichment in genes primarily involved in the IGF1/PI3K/NF κB/ERK gene signalling networks. CONCLUSIONS: This study has identified pathway specific prognostic biomarkers possibly underlying a differential chemotherapy response in patients undergoing standard platinum-based treatment of serous epithelial ovarian cancer. In addition, our results provide a pathway context for further experimental validations, and the findings are a significant step towards future therapeutic interventions.

PTEN genomic deletions that characterize aggressive prostate cancer originate close to segmental duplications.

Deletion of PTEN at 10q23.3 occurs in ∼40% of human prostate cancers and is associated with aggressive metastatic potential, poor prognosis, and androgen-independence. This high frequency of recurrent PTEN deletions in prostate cancer suggests there may be unusual genomic features close to this locus that facilitate DNA alteration at 10q23.3. To explore possible mechanisms for deletions in the PTEN region, a meta-analysis of 311 published human genome array datasets was conducted and determined that the minimal prostate cancer-associated deletion at 10q23.3 corresponds to ∼2.06 MB region flanked by BMPR1A and FAS. On a separate cohort comprising an additional 330 tumors, four-color fluorescence in situ hybridization analysis using probes for BMPR1A, FAS, cen(10), and PTEN showed that 132 of 330 (40%) tumors had PTEN loss, 50 (15%) of which were homozygous losses (comprising in total 100 deletion events). Breakpoints between PTEN and BMPR1A or FAS were subsequently mapped in 100 homozygous and 82 hemizygous PTEN losses, revealing that 125/182 PTEN microdeletions occurred within the 940 kB interval between BMPR1A and PTEN. Furthermore, this breakpoint interval coincides with a repeat-rich region of 414 kB containing the SD17 and SD18 segmental duplications, which contain at least 13 homologous inverted repeat sequences. Together, these data suggest that a strong selective growth advantage for loss of PTEN and upregulation of PI3K/AKT, combined with the close proximity of PTEN to a large unstable segment of repeated DNA comprising SD17 and SD18, can lead to recurrent microdeletions of the PTEN gene in prostate cancer. © 2011 Wiley Periodicals, Inc.

Transcriptomic and metabolomic analysis of Yukon Thellungiella plants grown in cabinets and their natural habitat show phenotypic plasticity.

BACKGROUND: Thellungiella salsuginea is an important model plant due to its natural tolerance to abiotic stresses including salt, cold, and water deficits. Microarray and metabolite profiling have shown that Thellungiella undergoes stress-responsive changes in transcript and organic solute abundance when grown under controlled environmental conditions. However, few reports assess the capacity of plants to display stress-responsive traits in natural habitats where concurrent stresses are the norm. RESULTS: To determine whether stress-responsive changes observed in cabinet-grown plants are recapitulated in the field, we analyzed leaf transcript and metabolic profiles of Thellungiella growing in its native Yukon habitat during two years of contrasting meteorological conditions. We found 673 genes showing differential expression between field and unstressed, chamber-grown plants. There were comparatively few overlaps between genes expressed under field and cabinet treatment-specific conditions. Only 20 of 99 drought-responsive genes were expressed both in the field during a year of low precipitation and in plants subjected to drought treatments in cabinets. There was also a general pattern of lower abundance among metabolites found in field plants relative to control or stress-treated plants in growth cabinets. Nutrient availability may explain some of the observed differences. For example, proline accumulated to high levels in cold and salt-stressed cabinet-grown plants but proline content was, by comparison, negligible in plants at a saline Yukon field site. We show that proline accumulated in a stress-responsive manner in Thellungiella plants salinized in growth cabinets and in salt-stressed seedlings when nitrogen was provided at 1.0 mM. In seedlings grown on 0.1 mM nitrogen medium, the proline content was low while carbohydrates increased. The relatively higher content of sugar-like compounds in field plants and seedlings on low nitrogen media suggests that Thellungiella shows metabolic plasticity in response to environmental stress and that resource availability can influence the expression of stress tolerance traits under field conditions. CONCLUSION: Comparisons between Thellungiella plants responding to stress in cabinets and in their natural habitats showed differences but also overlap between transcript and metabolite profiles. The traits in common offer potential targets for improving crops that must respond appropriately to multiple, concurrent stresses.

EMT transcription factors snail and slug directly contribute to cisplatin resistance in ovarian cancer.

BACKGROUND: The epithelial to mesenchymal transition (EMT) is a molecular process through which an epithelial cell undergoes transdifferentiation into a mesenchymal phenotype. The role of EMT in embryogenesis is well-characterized and increasing evidence suggests that elements of the transition may be important in other processes, including metastasis and drug resistance in various different cancers. METHODS: Agilent 4 × 44 K whole human genome arrays and selected reaction monitoring mass spectrometry were used to investigate mRNA and protein expression in A2780 cisplatin sensitive and resistant cell lines. Invasion and migration were assessed using Boyden chamber assays. Gene knockdown of snail and slug was done using targeted siRNA. Clinical relevance of the EMT pathway was assessed in a cohort of primary ovarian tumours using data from Affymetrix GeneChip Human Genome U133 plus 2.0 arrays. RESULTS: Morphological and phenotypic hallmarks of EMT were identified in the chemoresistant cells. Subsequent gene expression profiling revealed upregulation of EMT-related transcription factors including snail, slug, twist2 and zeb2. Proteomic analysis demonstrated up regulation of Snail and Slug as well as the mesenchymal marker Vimentin, and down regulation of E-cadherin, an epithelial marker. By reducing expression of snail and slug, the mesenchymal phenotype was largely reversed and cells were resensitized to cisplatin. Finally, gene expression data from primary tumours mirrored the finding that an EMT-like pathway is activated in resistant tumours relative to sensitive tumours, suggesting that the involvement of this transition may not be limited to in vitro drug effects. CONCLUSIONS: This work strongly suggests that genes associated with EMT may play a significant role in cisplatin resistance in ovarian cancer, therefore potentially leading to the development of predictive biomarkers of drug response or novel therapeutic strategies for overcoming drug resistance.

Isolation of a Human Betaretrovirus from Patients with Primary Biliary Cholangitis.

A human betaretrovirus (HBRV) has been linked with the autoimmune liver disease, primary biliary cholangitis (PBC), and various cancers, including breast cancer and lymphoma. HBRV is closely related to the mouse mammary tumor virus, and represents the only exogenous betaretrovirus characterized in humans to date. Evidence of infection in patients with PBC has been demonstrated through the identification of proviral integration sites in lymphoid tissue, the major reservoir of infection, as well as biliary epithelium, which is the site of the disease process. Accordingly, we tested the hypothesis that patients with PBC harbor a transmissible betaretrovirus by co-cultivation of PBC patients' lymph node homogenates with the HS578T breast cancer line. Because of the low level of HBRV replication, betaretrovirus producing cells were subcloned to optimize viral isolation and production. Evidence of infection was provided by electron microscopy, RT-PCR, in situ hybridization, cloning of the HBRV proviral genome and demonstration of more than 3400 integration sites. Further evidence of viral transmissibility was demonstrated by infection of biliary epithelial cells. While HBRV did not show a preference for integration proximal to specific genomic features, analyses of common insertion sites revealed evidence of integration proximal to cancer associated genes. These studies demonstrate the isolation of HBRV with features similar to mouse mammary tumor virus and confirm that patients with PBC display evidence of a transmissible viral infection.

WormBase in 2022-data, processes, and tools for analyzing Caenorhabditis elegans.

WormBase (www.wormbase.org) is the central repository for the genetics and genomics of the nematode Caenorhabditis elegans. We provide the research community with data and tools to facilitate the use of C. elegans and related nematodes as model organisms for studying human health, development, and many aspects of fundamental biology. Throughout our 22-year history, we have continued to evolve to reflect progress and innovation in the science and technologies involved in the study of C. elegans. We strive to incorporate new data types and richer data sets, and to provide integrated displays and services that avail the knowledge generated by the published nematode genetics literature. Here, we provide a broad overview of the current state of WormBase in terms of data type, curation workflows, analysis, and tools, including exciting new advances for analysis of single-cell data, text mining and visualization, and the new community collaboration forum. Concurrently, we continue the integration and harmonization of infrastructure, processes, and tools with the Alliance of Genome Resources, of which WormBase is a founding member.

Design and verification of a pangenome microarray oligonucleotide probe set for Dehalococcoides spp.

Dehalococcoides spp. are an industrially relevant group of Chloroflexi bacteria capable of reductively dechlorinating contaminants in groundwater environments. Existing Dehalococcoides genomes revealed a high level of sequence identity within this group, including 98 to 100% 16S rRNA sequence identity between strains with diverse substrate specificities. Common molecular techniques for identification of microbial populations are often not applicable for distinguishing Dehalococcoides strains. Here we describe an oligonucleotide microarray probe set designed based on clustered Dehalococcoides genes from five different sources (strain DET195, CBDB1, BAV1, and VS genomes and the KB-1 metagenome). This "pangenome" probe set provides coverage of core Dehalococcoides genes as well as strain-specific genes while optimizing the potential for hybridization to closely related, previously unknown Dehalococcoides strains. The pangenome probe set was compared to probe sets designed independently for each of the five Dehalococcoides strains. The pangenome probe set demonstrated better predictability and higher detection of Dehalococcoides genes than strain-specific probe sets on nontarget strains with <99% average nucleotide identity. An in silico analysis of the expected probe hybridization against the recently released Dehalococcoides strain GT genome and additional KB-1 metagenome sequence data indicated that the pangenome probe set performs more robustly than the combined strain-specific probe sets in the detection of genes not included in the original design. The pangenome probe set represents a highly specific, universal tool for the detection and characterization of Dehalococcoides from contaminated sites. It has the potential to become a common platform for Dehalococcoides-focused research, allowing meaningful comparisons between microarray experiments regardless of the strain examined.

FOXQ1 is Differentially Expressed Across Breast Cancer Subtypes with Low Expression Associated with Poor Overall Survival.

PURPOSE: Forkhead box Q1 (FOXQ1) has been shown to contribute to the development and progression of cancers, including ovarian and breast cancer (BC). However, research exploring FOXQ1 expression, copy number variation (CNV), and prognostic value across different BC subtypes is limited. Our purpose was to evaluate FOXQ1 mRNA expression, CNV, and prognostic value across BC subtypes. MATERIALS AND METHODS: We determined FOXQ1 expression and CNV in BC patient tumors using RT-qPCR and qPCR, respectively. We also analyzed FOXQ1 expression and CNV in BC cell lines in the CCLE database using K-means clustering. The prognostic value of FOXQ1 expression in the TCGA-BRCA database was assessed using univariate and multivariate Cox's regression analysis as well as using the online tools OncoLnc, GEPIA, and UALCAN. RESULTS: Our analyses reveal that FOXQ1 mRNA is differentially expressed between different subtypes of BC and is significantly decreased in luminal BC and HER2 patients when compared to normal breast tissue samples. Furthermore, analysis of BC cell lines showed that FOXQ1 mRNA expression was independent of CNV. Moreover, patients with low FOXQ1 mRNA expression had significantly poorer overall survival compared to those with high FOXQ1 mRNA expression. Finally, low FOXQ1 expression had a critical impact on the prognostic values of BC patients and was an independent predictor of overall survival when it was adjusted for BC subtypes and to two other FOX genes, FOXF2 and FOXM1. CONCLUSION: Our study reveals for the first time that FOXQ1 is differentially expressed across BC subtypes and that low expression of FOXQ1 is indicative of poor prognosis in patients with BC.

The accuracy of several multiple sequence alignment programs for proteins.

BACKGROUND: There have been many algorithms and software programs implemented for the inference of multiple sequence alignments of protein and DNA sequences. The "true" alignment is usually unknown due to the incomplete knowledge of the evolutionary history of the sequences, making it difficult to gauge the relative accuracy of the programs. RESULTS: We tested nine of the most often used protein alignment programs and compared their results using sequences generated with the simulation software Simprot which creates known alignments under realistic and controlled evolutionary scenarios. We have simulated more than 30,000 alignment sets using various evolutionary histories in order to define strengths and weaknesses of each program tested. We found that alignment accuracy is extremely dependent on the number of insertions and deletions in the sequences, and that indel size has a weaker effect. We also considered benchmark alignments from the latest version of BAliBASE and the results relative to BAliBASE- and Simprot-generated data sets were consistent in most cases. CONCLUSION: Our results indicate that employing Simprot's simulated sequences allows the creation of a more flexible and broader range of alignment classes than the usual methods for alignment accuracy assessment. Simprot also allows for a quick and efficient analysis of a wider range of possible evolutionary histories that might not be present in currently available alignment sets. Among the nine programs tested, the iterative approach available in Mafft (L-INS-i) and ProbCons were consistently the most accurate, with Mafft being the faster of the two.

Interphase FISH analysis of PTEN in histologic sections shows genomic deletions in 68% of primary prostate cancer and 23% of high-grade prostatic intra-epithelial neoplasias.

Prostate cancer (CaP) is characterized by the accumulation of both genetic and epigenetic alterations that transform premalignant lesions to invasive carcinoma. However, the molecular events underlying this critical transition are poorly understood. One of the important genes that might play a role in CaP development is the PTEN gene. At the present time, there has been no systematic analysis of the incidence of genomic PTEN deletion by fluorescence in situ hybridization (FISH) in CaP and associated preneoplastic histologic lesions. This study assesses the frequency of PTEN deletion by interphase FISH analysis in CaP and prostatic intra-epithelial neoplasia (PIN). Dual-color FISH was performed using DNA probes for bands 10q23.3 (PTEN locus) and chromosome 10 centromere using 35 radical prostatectomy specimens. PTEN deletions were not found in 3/3 of stroma, 6/6 samples of benign glandular epithelium, and 12/12 samples of low-grade PIN. However, PTEN deletions were found in 3/13 (23%) of high-grade PIN and 24/35 (68%) of CaP. Concordance was observed between PTEN deletion status and the overall cellular PTEN protein expression levels, as assessed by immunohistochemistry. The high frequency of PTEN deletion observed in CaP versus precursor lesions implicates a pivotal role for PTEN haploinsufficiency in the transition from preneoplastic PIN to CaP. Moreover, this observation is an important consideration for novel therapeutic trials in CaP in which biologic efficacy is influenced by the activity level of PTEN. These findings draw attention to the usefulness of this relatively simple FISH assay for future applications in clinical laboratories.

Metaphase and array comparative genomic hybridization: unique copy number changes and gene amplification of medulloblastomas in South America.

Tumors of the central nervous system are the second most frequent malignancy of childhood, accounting for the majority of cancer-related deaths in this age group. Among these tumors, medulloblastomas (MB) remain in need of further genomic characterization toward understanding of pathogenesis and outcome predictors. Eight pediatric embryonal brain tumors were analyzed: five MB (one being desmoplastic), one PNET, one medulloepithelioma, and one ependymoblastoma. Analyses identified genomic imbalances, including the gain of 16p and the nonsyntenic coamplification of MYCN and TERT loci. More detailed FISH analysis showed that coamplification of MYCN and TERT in one of the MBs manifested as dispersed nuclear speckling, consistent with the presence of double minute chromosomes. There was considerable cell-to-cell copy number heterogeneity present, but it was clear that both genes were amplified concordantly. The amplification of oncogenes seems to play an important role in the pathogenesis of MB, and the association between MYCN and TERT amplifications and poor prognosis has not been well recognized. The uncharacteristic pattern of genomic imbalances detected in MB tumors may be a reflection of the characteristics of these tumors occurring in South America.

SIMPROT: using an empirically determined indel distribution in simulations of protein evolution.

BACKGROUND: General protein evolution models help determine the baseline expectations for the evolution of sequences, and they have been extensively useful in sequence analysis and for the computer simulation of artificial sequence data sets. RESULTS: We have developed a new method of simulating protein sequence evolution, including insertion and deletion (indel) events in addition to amino-acid substitutions. The simulation generates both the simulated sequence family and a true sequence alignment that captures the evolutionary relationships between amino acids from different sequences. Our statistical model for indel evolution is based on the empirical indel distribution determined by Qian and Goldstein. We have parameterized this distribution so that it applies to sequences diverged by varying evolutionary times and generalized it to provide flexibility in simulation conditions. Our method uses a Monte-Carlo simulation strategy, and has been implemented in a C++ program named Simprot. CONCLUSION: Simprot will be useful for testing methods of analysis of protein sequence families particularly alignment methods, phylogenetic tree building, detection of recombination and horizontal gene transfer, and homology detection, where knowing the true course of sequence evolution is essential.