Target capture sequencing of SARS-CoV-2 genomes using the ONETest Coronaviruses Plus

We introduce a target capture next-generation sequencing methodology, the ONETest Coronaviruses Plus, to sequence the SARS-CoV-2 genome and select loci of other respiratory viruses. We applied the ONETest on 70 respiratory samples (collected in Florida, USA between May and July, 2020), in which SARS-CoV-2 had been detected by a PCR assay. For 48 of the samples, we also applied the ARTIC protocol. Of the 70 ONETest libraries, 45 (64%) had a (near-)complete sequence (>29,000 bases and >90% covered by >9 reads). Of the 48 ARTIC libraries, 25 (52%) had a (near-)complete sequence. In 19 out of 25 (76%) samples in which both the ONETest and ARTIC yielded (near-)complete sequences, the lineages assigned were identical. As a target capture approach, the ONETest is less prone to loss of sequence coverage than amplicon approaches, and thus can provide complete genomic information more often to track and monitor SARS-CoV-2 variants.

Iatrogenic endometriosis harbors somatic cancer-driver mutations

STUDY QUESTION

Does incisional endometriosis (IE) harbor somatic cancer-driver mutations?

SUMMARY ANSWER

We found that approximately one-quarter of IE cases harbor somatic-cancer mutations, which commonly affect components of the MAPK/RAS or PI3K-Akt-mTor signaling pathways.

WHAT IS KNOWN ALREADY

Despite the classification of endometriosis as a benign gynecological disease, it shares key features with cancers such as resistance to apoptosis and stimulation of angiogenesis and is well-established as the precursor of clear cell and endometrioid ovarian carcinomas. Our group has recently shown that deep infiltrating endometriosis (DE), a form of endometriosis that rarely undergoes malignant transformation, harbors recurrent somatic mutations.

STUDY DESIGN, SIZE, DURATION

In a retrospective study comparing iatrogenically induced and endogenously occurring forms of endometriosis unlikely to progress to cancer, we examined endometriosis specimens from 40 women with IE and 36 women with DE. Specimens were collected between 2004 and 2017 from five hospital sites in either Canada, Germany or the Netherlands. IE and DE cohorts were age-matched and all women presented with histologically typical endometriosis without known history of malignancy.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Archival tissue specimens containing endometriotic lesions were macrodissected and/or laser-capture microdissected to enrich endometriotic stroma and epithelium and a hypersensitive cancer hotspot sequencing panel was used to assess for presence of somatic mutations. Mutations were subsequently validated using droplet digital PCR. PTEN and ARID1A immunohistochemistry (IHC) were performed as surrogates for somatic events resulting in functional loss of respective proteins.

MAIN RESULTS AND THE ROLE OF CHANCE

Overall, we detected somatic cancer-driver events in 11 of 40 (27.5%) IE cases and 13 of 36 (36.1%) DE cases, including hotspot mutations in KRAS, ERBB2, PIK3CA and CTNNB1. Heterogeneous PTEN loss occurred at similar rates in IE and DE (7/40 vs 5/36, respectively), whereas ARID1A loss only occurred in a single case of DE. While rates of detectable somatic cancer-driver events between IE and DE are not statistically significant (P > 0.05), KRAS activating mutations were more prevalent in DE.

LIMITATIONS, REASONS FOR CAUTION

Detection of somatic cancer-driver events were limited to hotspots analyzed in our panel-based sequencing assay and loss of protein expression by IHC from archival tissue. Whole genome or exome sequencing, or epigenetic analysis may uncover additional somatic alterations. Moreover, because of the descriptive nature of this study, the functional roles of identified mutations within the context of endometriosis remain unclear and causality cannot be established.

WIDER IMPLICATIONS OF THE FINDINGS

The alterations we report may be important in driving the growth and survival of endometriosis in ectopic regions of the body. Given the frequency of mutation in surgically displaced endometrium (IE), examination of similar somatic events in eutopic endometrium, as well as clinically annotated cases of other forms of endometriosis, in particular endometriomas that are most commonly linked to malignancy, is warranted.

STUDY FUNDING/COMPETING INTEREST(S)

This study was funded by a Canadian Cancer Society Impact Grant [701603, PI Huntsman], Canadian Institutes of Health Research Transitional Open Operating Grant [MOP-142273, PI Yong], the Canadian Institutes of Health Research Foundation Grant [FDN-154290, PI Huntsman], the Canadian Institutes of Health Research Project Grant [PJT-156084, PIs Yong and Anglesio], and the Janet D. Cottrelle Foundation through the BC Cancer Foundation [PI Huntsman]. D.G. Huntsman is a co-founder and shareholder of Contextual Genomics Inc., a for profit company that provides clinical reporting to assist in cancer patient treatment. R. Aguirre-Hernandez, J. Khattra and L.M. Prentice have a patent MOLECULAR QUALITY ASSURANCE METHODS FOR USE IN SEQUENCING pending and are current (or former) employees of Contextual Genomics Inc. The remaining authors have no competing interests to declare.

TRIAL REGISTRATION NUMBER

Not applicable.

Divergent modes of clonal spread and intraperitoneal mixing in high-grade serous ovarian cancer

We performed phylogenetic analysis of high-grade serous ovarian cancers (68 samples from seven patients), identifying constituent clones and quantifying their relative abundances at multiple intraperitoneal sites. Through whole-genome and single-nucleus sequencing, we identified evolutionary features including mutation loss, convergence of the structural genome and temporal activation of mutational processes that patterned clonal progression. We then determined the precise clonal mixtures comprising each tumor sample. The majority of sites were clonally pure or composed of clones from a single phylogenetic clade. However, each patient contained at least one site composed of polyphyletic clones. Five patients exhibited monoclonal and unidirectional seeding from the ovary to intraperitoneal sites, and two patients demonstrated polyclonal spread and reseeding. Our findings indicate that at least two distinct modes of intraperitoneal spread operate in clonal dissemination and highlight the distribution of migratory potential over clonal populations comprising high-grade serous ovarian cancers.

Dynamics of genomic clones in breast cancer patient xenografts at single-cell resolution

Human cancers, including breast cancers, comprise clones differing in mutation content. Clones evolve dynamically in space and time following principles of Darwinian evolution, underpinning important emergent features such as drug resistance and metastasis. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours. However, the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours have not been systematically examined at single-cell resolution. Here we show, using deep-genome and single-cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all 15 cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (<5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single-cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies using patient-derived breast cancer xenoengraftment.

TITAN: inference of copy number architectures in clonal cell populations from tumor whole-genome sequence data

The evolution of cancer genomes within a single tumor creates mixed cell populations with divergent somatic mutational landscapes. Inference of tumor subpopulations has been disproportionately focused on the assessment of somatic point mutations, whereas computational methods targeting evolutionary dynamics of copy number alterations (CNA) and loss of heterozygosity (LOH) in whole-genome sequencing data remain underdeveloped. We present a novel probabilistic model, TITAN, to infer CNA and LOH events while accounting for mixtures of cell populations, thereby estimating the proportion of cells harboring each event. We evaluate TITAN on idealized mixtures, simulating clonal populations from whole-genome sequences taken from genomically heterogeneous ovarian tumor sites collected from the same patient. In addition, we show in 23 whole genomes of breast tumors that the inference of CNA and LOH using TITAN critically informs population structure and the nature of the evolving cancer genome. Finally, we experimentally validated subclonal predictions using fluorescence in situ hybridization (FISH) and single-cell sequencing from an ovarian cancer patient sample, thereby recapitulating the key modeling assumptions of TITAN.

PyClone: statistical inference of clonal population structure in cancer

We introduce PyClone, a statistical model for inference of clonal population structures in cancers. PyClone is a Bayesian clustering method for grouping sets of deeply sequenced somatic mutations into putative clonal clusters while estimating their cellular prevalences and accounting for allelic imbalances introduced by segmental copy-number changes and normal-cell contamination. Single-cell sequencing validation demonstrates PyClone's accuracy.

An integrated strategy to study muscle development and myofilament structure in Caenorhabditis elegans

A crucial step in the development of muscle cells in all metazoan animals is the assembly and anchorage of the sarcomere, the essential repeat unit responsible for muscle contraction. In Caenorhabditis elegans, many of the critical proteins involved in this process have been uncovered through mutational screens focusing on uncoordinated movement and embryonic arrest phenotypes. We propose that additional sarcomeric proteins exist for which there is a less severe, or entirely different, mutant phenotype produced in their absence. We have used Serial Analysis of Gene Expression (SAGE) to generate a comprehensive profile of late embryonic muscle gene expression. We generated two replicate long SAGE libraries for sorted embryonic muscle cells, identifying 7,974 protein-coding genes. A refined list of 3,577 genes expressed in muscle cells was compiled from the overlap between our SAGE data and available microarray data. Using the genes in our refined list, we have performed two separate RNA interference (RNAi) screens to identify novel genes that play a role in sarcomere assembly and/or maintenance in either embryonic or adult muscle. To identify muscle defects in embryos, we screened specifically for the Pat embryonic arrest phenotype. To visualize muscle defects in adult animals, we fed dsRNA to worms producing a GFP-tagged myosin protein, thus allowing us to analyze their myofilament organization under gene knockdown conditions using fluorescence microscopy. By eliminating or severely reducing the expression of 3,300 genes using RNAi, we identified 122 genes necessary for proper myofilament organization, 108 of which are genes without a previously characterized role in muscle. Many of the genes affecting sarcomere integrity have human homologs for which little or nothing is known.

Muscular diseases affect many people worldwide. While we have learned much about the sarcomere, the basic building block of muscle cells, there are still numerous questions that remain to be answered. We must learn more about proteins expressed in muscle and how they interact so that better treatments for myopathies can be developed. The nematode Caenorhabditis elegans is a valuable model organism for the study of muscle due to similarities between worm body wall muscle and vertebrate muscle, along with its semi-transparent cuticle that allows for visualization of muscle structures in live animals. We have used transcriptional profiling methods to identify the majority of genes that are expressed in the embryonic body wall muscle cells of C. elegans. To gain insight into possible functions performed by these genes and their corresponding proteins, we examined animals and muscle cells for abnormalities after the targeted inactivation of about 3,300 genes. We identified 122 genes necessary for proper myofilament organization, 108 of which had no previously characterized role in muscle. This approach proved to be a rapid and sensitive means to identify genes that affect muscle differentiation and sarcomere assembly.

Mutation of FOXL2 in granulosa-cell tumors of the ovary

BACKGROUND

Granulosa-cell tumors (GCTs) are the most common type of malignant ovarian sex cord-stromal tumor (SCST). The pathogenesis of these tumors is unknown. Moreover, their histopathological diagnosis can be challenging, and there is no curative treatment beyond surgery.

METHODS

We analyzed four adult-type GCTs using whole-transcriptome paired-end RNA sequencing. We identified putative GCT-specific mutations that were present in at least three of these samples but were absent from the transcriptomes of 11 epithelial ovarian tumors, published human genomes, and databases of single-nucleotide polymorphisms. We confirmed these variants by direct sequencing of complementary DNA and genomic DNA. We then analyzed additional tumors and matched normal genomic DNA, using a combination of direct sequencing, analyses of restriction-fragment-length polymorphisms, and TaqMan assays.

RESULTS

All four index GCTs had a missense point mutation, 402C-->G (C134W), in FOXL2, a gene encoding a transcription factor known to be critical for granulosa-cell development. The FOXL2 mutation was present in 86 of 89 additional adult-type GCTs (97%), in 3 of 14 thecomas (21%), and in 1 of 10 juvenile-type GCTs (10%). The mutation was absent in 49 SCSTs of other types and in 329 unrelated ovarian or breast tumors.

CONCLUSIONS

Whole-transcriptome sequencing of four GCTs identified a single, recurrent somatic mutation (402C-->G) in FOXL2 that was present in almost all morphologically identified adult-type GCTs. Mutant FOXL2 is a potential driver in the pathogenesis of adult-type GCTs.

LongSAGE profiling of nine human embryonic stem cell lines

To facilitate discovery of novel human embryonic stem cell (ESC) transcripts, we generated 2.5 million LongSAGE tags from 9 human ESC lines. Analysis of this data revealed that ESCs express proportionately more RNA binding proteins compared with terminally differentiated cells, and identified novel ESC transcripts, at least one of which may represent a marker of the pluripotent state.

A Modified Polymerase Chain Reaction-Long Serial Analysis of Gene Expression Protocol Identifies Novel Transcripts in Human CD34+Bone Marrow Cells

Transcriptome profiling offers a powerful approach to investigating developmental processes. Long serial analysis of gene expression (LongSAGE) is particularly attractive for this purpose because of its inherent quantitative features and independence of both hybridization variables and prior knowledge of transcript identity. Here, we describe the validation and initial application of a modified protocol for amplifying cDNA preparations from <10 ng of RNA (<10(3) cells) to allow representative LongSAGE libraries to be constructed from rare stem cell-enriched populations. Quantitative real-time polymerase chain reaction (Q-RT-PCR) analyses and comparison of tag frequencies in replicate LongSAGE libraries produced from amplified and nonamplified cDNA preparations demonstrated preservation of the relative levels of different transcripts originally present at widely varying levels. This PCR-LongSAGE protocol was then used to obtain a 200,000-tag library from the CD34+ subset of normal adult human bone marrow cells. Analysis of this library revealed many anticipated transcripts, as well as transcripts not previously known to be present in CD34+ hematopoietic cells. The latter included numerous novel tags that mapped to unique and conserved sites in the human genome but not previously identified as transcribed elements in human cells. Q-RT-PCR was used to demonstrate that 10 of these novel tags were expressed in cDNA pools and present in extracts of other sources of normal human CD34+ hematopoietic cells. These findings illustrate the power of LongSAGE to identify new transcripts in stem cell-enriched populations and indicate the potential of this approach to be extended to other sources of rare cells. Disclosure of potential conflicts of interest is found at the end of this article.

Identification of ciliary and ciliopathy genes in Caenorhabditis elegans through comparative genomics

Comparative genomic analysis of three nematode species identifies 93 genes that encode putative components of the ciliated neurons in C. elegans and are subject to the same regulatory control.

Background

The recent availability of genome sequences of multiple related Caenorhabditis species has made it possible to identify, using comparative genomics, similarly transcribed genes in Caenorhabditis elegans and its sister species. Taking this approach, we have identified numerous novel ciliary genes in C. elegans, some of which may be orthologs of unidentified human ciliopathy genes.

Results

By screening for genes possessing canonical X-box sequences in promoters of three Caenorhabditis species, namely C. elegans, C. briggsae and C. remanei, we identified 93 genes (including known X-box regulated genes) that encode putative components of ciliated neurons in C. elegans and are subject to the same regulatory control. For many of these genes, restricted anatomical expression in ciliated cells was confirmed, and control of transcription by the ciliogenic DAF-19 RFX transcription factor was demonstrated by comparative transcriptional profiling of different tissue types and of daf-19(+) and daf-19(-) animals. Finally, we demonstrate that the dye-filling defect of dyf-5(mn400) animals, which is indicative of compromised exposure of cilia to the environment, is caused by a nonsense mutation in the serine/threonine protein kinase gene M04C9.5.

Conclusion

Our comparative genomics-based predictions may be useful for identifying genes involved in human ciliopathies, including Bardet-Biedl Syndrome (BBS), since the C. elegans orthologs of known human BBS genes contain X-box motifs and are required for normal dye filling in C. elegans ciliated neurons.

Large-scale production of SAGE libraries from microdissected tissues, flow-sorted cells, and cell lines

We describe the details of a serial analysis of gene expression (SAGE) library construction and analysis platform that has enabled the generation of >298 high-quality SAGE libraries and >30 million SAGE tags primarily from sub-microgram amounts of total RNA purified from samples acquired by microdissection. Several RNA isolation methods were used to handle the diversity of samples processed, and various measures were applied to minimize ditag PCR carryover contamination. Modifications in the SAGE protocol resulted in improved cloning and DNA sequencing efficiencies. Bioinformatic measures to automatically assess DNA sequencing results were implemented to analyze the integrity of ditag structure, linker or cross-species ditag contamination, and yield of high-quality tags per sequence read. Our analysis of singleton tag errors resulted in a method for correcting such errors to statistically determine tag accuracy. From the libraries generated, we produced an essentially complete mapping of reliable 21-base-pair tags to the mouse reference genome sequence for a meta-library of approximately 5 million tags. Our analyses led us to reject the commonly held notion that duplicate ditags are artifacts. Rather than the usual practice of discarding such tags, we conclude that they should be retained to avoid introducing bias into the results and thereby maintain the quantitative nature of the data, which is a major theoretical advantage of SAGE as a tool for global transcriptional profiling.

The ELT-2 GATA-factor and the global regulation of transcription in the C. elegans intestine

A SAGE library was prepared from hand-dissected intestines from adult Caenorhabditis elegans, allowing the identification of >4000 intestinally-expressed genes; this gene inventory provides fundamental information for understanding intestine function, structure and development. Intestinally-expressed genes fall into two broad classes: widely-expressed "housekeeping" genes and genes that are either intestine-specific or significantly intestine-enriched. Within this latter class of genes, we identified a subset of highly-expressed highly-validated genes that are expressed either exclusively or primarily in the intestine. Over half of the encoded proteins are candidates for secretion into the intestinal lumen to hydrolyze the bacterial food (e.g. lysozymes, amoebapores, lipases and especially proteases). The promoters of this subset of intestine-specific/intestine-enriched genes were analyzed computationally, using both a word-counting method (RSAT oligo-analysis) and a method based on Gibbs sampling (MotifSampler). Both methods returned the same over-represented site, namely an extended GATA-related sequence of the general form AHTGATAARR, which agrees with experimentally determined cis-acting control sequences found in intestine genes over the past 20 years. All promoters in the subset contain such a site, compared to <5% for control promoters; moreover, our analysis suggests that the majority (perhaps all) of genes expressed exclusively or primarily in the worm intestine are likely to contain such a site in their promoters. There are three zinc-finger GATA-type factors that are candidates to bind this extended GATA site in the differentiating C. elegans intestine: ELT-2, ELT-4 and ELT-7. All evidence points to ELT-2 being the most important of the three. We show that worms in which both the elt-4 and the elt-7 genes have been deleted from the genome are essentially wildtype, demonstrating that ELT-2 provides all essential GATA-factor functions in the intestine. The SAGE analysis also identifies more than a hundred other transcription factors in the adult intestine but few show an RNAi-induced loss-of-function phenotype and none (other than ELT-2) show a phenotype primarily in the intestine. We thus propose a simple model in which the ELT-2 GATA factor directly participates in the transcription of all intestine-specific/intestine-enriched genes, from the early embryo through to the dying adult. Other intestinal transcription factors would thus modulate the action of ELT-2, depending on the worm's nutritional and physiological needs.

Transgene constructs in coho salmon (Oncorhynchus kisutch) are repeated in a head-to-tail fashion and can be integrated adjacent to horizontally-transmitted parasite DNA

Currently, little information is available regarding the molecular organization of integrated transgenes in genetically-engineered fish. We performed a detailed structural analysis of an inserted transgene in one strain (M77) of transgenic coho salmon (Oncorhynchus kisutch) containing a salmon growth hormone gene construct (OnMTGH1). Microinjected DNA was found to have inserted into a single site in the coho salmon genome, and was organized with four complete internal copies and two partial terminal copies of the OnMTGH1 construct. All construct copies were organized in a direct-tandem (head-to-tail) repeat fashion in strain M77 and five additional strains (one also possessed a second recombinant junction fragment). For strain M77, the junctions between the transgene insert and the insertion point within the wild-type genome were cloned from strain-specific cosmid libraries and sequenced, revealing that the transgene insertion was accompanied by a deletion of 587 bp of wild-type DNA as well as a small insertion (19 bp) of unknown DNA upstream and a 14 bp direct- tandem duplication of sequence downstream. Upstream and downstream wild-type DNA sequence contained several repetitive sequence elements based on Southern blot analysis and homology to repetitive sequences in GenBank. In the downstream flank, a pseudogene sequence was also identified which has high homology to the CA membrane protein gene from Schistosoma japonicum, a parasite closely related to Sanguinicola sp. parasites which infect salmonids. Whether the presence of an inserted transgene and the presence of potentially horizontally-transmitted DNA are indicative of a genomic region with a predisposition for insertion of foreign DNA requires further study. The information derived from this transgene structure provides information useful for comparison to other transgenic organisms and for determination of the mechanism of transgene integration in lower vertebrates.

A mouse atlas of gene expression: large-scale digital gene-expression profiles from precisely defined developing C57BL/6J mouse tissues and cells

We analyzed 8.55 million LongSAGE tags generated from 72 libraries. Each LongSAGE library was prepared from a different mouse tissue. Analysis of the data revealed extensive overlap with existing gene data sets and evidence for the existence of approximately 24,000 previously undescribed genomic loci. The visual cortex, pancreas, mammary gland, preimplantation embryo, and placenta contain the largest number of differentially expressed transcripts, 25% of which are previously undescribed loci.

Functional genomics of the cilium, a sensory organelle

Cilia and flagella play important roles in many physiological processes, including cell and fluid movement, sensory perception, and development. The biogenesis and maintenance of cilia depend on intraflagellar transport (IFT), a motility process that operates bidirectionally along the ciliary axoneme. Disruption in IFT and cilia function causes several human disorders, including polycystic kidneys, retinal dystrophy, neurosensory impairment, and Bardet-Biedl syndrome (BBS). To uncover new ciliary components, including IFT proteins, we compared C. elegans ciliated neuronal and nonciliated cells through serial analysis of gene expression (SAGE) and screened for genes potentially regulated by the ciliogenic transcription factor, DAF-19. Using these complementary approaches, we identified numerous candidate ciliary genes and confirmed the ciliated-cell-specific expression of 14 novel genes. One of these, C27H5.7a, encodes a ciliary protein that undergoes IFT. As with other IFT proteins, its ciliary localization and transport is disrupted by mutations in IFT and bbs genes. Furthermore, we demonstrate that the ciliary structural defect of C. elegans dyf-13(mn396) mutants is caused by a mutation in C27H5.7a. Together, our findings help define a ciliary transcriptome and suggest that DYF-13, an evolutionarily conserved protein, is a novel core IFT component required for cilia function.

Functional characterization of a catabolic plasmid from polychlorinated- biphenyl-degrading Rhodococcus sp. strain RHA1

Rhodococcus sp. strain RHA1, a potent polychlorinated-biphenyl (PCB)-degrading strain, contains three linear plasmids ranging in size from 330 to 1,100 kb. As part of a genome sequencing project, we report here the complete sequence and characterization of the smallest and least-well-characterized of the RHA1 plasmids, pRHL3. The plasmid is an actinomycete invertron, containing large terminal inverted repeats with a tightly associated protein and a predicted open reading frame (ORF) that is similar to that of a mycobacterial rep gene. The pRHL3 plasmid has 300 putative genes, almost 21% of which are predicted to have a catabolic function. Most of these are organized into three clusters. One of the catabolic clusters was predicted to include limonene degradation genes. Consistent with this prediction, RHA1 grew on limonene, carveol, or carvone as the sole carbon source. The plasmid carries three cytochrome P450-encoding (CYP) genes, a finding consistent with the high number of CYP genes found in other actinomycetes. Two of the CYP genes appear to belong to novel families; the third belongs to CYP family 116 but appears to belong to a novel class based on the predicted domain structure of its reductase. Analyses indicate that pRHL3 also contains four putative "genomic islands" (likely to have been acquired by horizontal transfer), insertion sequence elements, 19 transposase genes, and a duplication that spans two ORFs. One of the genomic islands appears to encode resistance to heavy metals. The plasmid does not appear to contain any housekeeping genes. However, each of the three catabolic clusters contains related genes that appear to be involved in glucose metabolism.

Systematic recovery and analysis of full-ORF human cDNA clones

The Mammalian Gene Collection (MGC) consortium (http://mgc.nci.nih.gov) seeks to establish publicly available collections of full-ORF cDNAs for several organisms of significance to biomedical research, including human. To date over 15,200 human cDNA clones containing full-length open reading frames (ORFs) have been identified via systematic expressed sequence tag (EST) analysis of a diverse set of cDNA libraries; however, further systematic EST analysis is no longer an efficient method for identifying new cDNAs. As part of our involvement in the MGC program, we have developed a scalable method for targeted recovery of cDNA clones to facilitate recovery of genes absent from the MGC collection. First, cDNA is synthesized from various RNAs, followed by polymerase chain reaction (PCR) amplification of transcripts in 96-well plates using gene-specific primer pairs flanking the ORFs. Amplicons are cloned into a sequencing vector, and full-length sequences are obtained. Sequences are processed and assembled using Phred and Phrap, and analyzed using Consed and a number of bioinformatics methods we have developed. Sequences are compared with the Reference Sequence (RefSeq) database, and validation of sequence discrepancies is attempted using other sequence databases including dbEST and dbSNP. Clones with identical sequence to RefSeq or containing only validated changes will become part of the MGC human gene collection. Clones containing novel splice variants or polymorphisms have also been identified. Our approach to clone recovery, applied at large scale, has the potential to recover many and possibly most of the genes absent from the MGC collection.

The Genome sequence of the SARS-associated coronavirus

We sequenced the 29,751-base genome of the severe acute respiratory syndrome (SARS)-associated coronavirus known as the Tor2 isolate. The genome sequence reveals that this coronavirus is only moderately related to other known coronaviruses, including two human coronaviruses, HCoV-OC43 and HCoV-229E. Phylogenetic analysis of the predicted viral proteins indicates that the virus does not closely resemble any of the three previously known groups of coronaviruses. The genome sequence will aid in the diagnosis of SARS virus infection in humans and potential animal hosts (using polymerase chain reaction and immunological tests), in the development of antivirals (including neutralizing antibodies), and in the identification of putative epitopes for vaccine development.

Recent advances in our knowledge of the Myxozoa

In the last few years two factors have helped to significantly advance our understanding of the Myxozoa. First, the phenomenal increase in fin fish aquaculture in the 1990s has lead to the increased importance of these parasites; in turn this has lead to intensified research efforts, which have increased knowledge of the development, diagnosis. and pathogenesis of myxozoans. The hallmark discovery in the 1980s that the life cycle of Myxobolus cerebralis requires development of an actinosporean stage in the oligochaete. Tubifex tubifex, led to the elucidation of the life cycles of several other myxozoans. Also, the life cycle and taxonomy of the enigmatic PKX myxozoan has been resolved: it is the alternate stage of the unusual myxozoan, Tetracapsula bryosalmonae, from bryozoans. The 18S rDNA gene of many species has been sequenced, and here we add 22 new sequences to the data set. Phylogenetic analyses using all these sequences indicate that: 1) the Myxozoa are closely related to Cnidaria (also supported by morphological data); 2) marine taxa at the genus level branch separately from genera that usually infect freshwater fishes; 3) taxa cluster more by development and tissue location than by spore morphology; 4) the tetracapsulids branched off early in myxozoan evolution, perhaps reflected by their having bryozoan, rather than annelid hosts; 5) the morphology of actinosporeans offers little information for determining their myxosporean counterparts (assuming that they exist); and 6) the marine actinosporeans from Australia appear to form a clade within the platysporinid myxosporeans. Ribosomal DNA sequences have also enabled development of diagnostic tests for myxozoans. PCR and in situ hybridisation tests based on rDNA sequences have been developed for Myxobolus cerebralis, Ceratomyxa shasta, Kudoa spp., and Tetracapsula bryosalmonae (PKX). Lectin-based and antibody tests have also been developed for certain myxozoans, such as PKX and C. shasta. We also review important diseases caused by myxozoans, which are emerging or re-emerging. Epizootics of whirling disease in wild rainbow trout (Oncorhynchus mykiss) have recently been reported throughout the Rocky Mountain states of the USA. With a dramatic increase in aquaculture of fishes using marine netpens, several marine myxozoans have been recognized or elevated in status as pathological agents. Kudoa thyrsites infections have caused severe post-harvest myoliquefaction in pen-reared Atlantic salmon (Salmo salar), and Ceratomyxa spp., Sphaerospora spp., and Myxidium leei cause disease in pen-reared sea bass (Dicentrarchus labrax) and sea bream species (family Sparidae) in Mediterranean countries.

Tetracapsula renicola n. sp. (Myxozoa : Saccosporidae); the PKX myxozoan--the cause of proliferative kidney disease of salmonid fishes

Proliferative kidney disease (PKD) of salmonid fishes is caused by the extrasporogonic stage of an enigmatic myxozoan, referred to as PKX. Sporogenesis occurs in the renal tubules, resulting in monosporous pseudoplasmodia. The spores are ovoid with indistinguishable valves and measure 12 microm in length and 7 microm in width. Two spherical polar capsules (2 microm diameter) with 4 coils occur at the anterior end of the spore. Prominent capsulogenic cell nuclei posterior to the polar capsules are evident in histological sections stained with hematoxylin and eosin. Regardless of the true nature of the valves (indistinguishable or absent), this myxozoan is morphologically distinct from all other described members of the phylum Myxozoa. Comparisons of small subunit rDNA sequences of PKX with other myxozoans demonstrated that it branches from all other members of the myxosporeans from fish examined thus far, including representatives of the phenotypically most closely related genera, Sphaerospora and Parvicapsula. Recent reports, based on rDNA comparisons, indicate that the alternate stage of PKX occurs in bryozoans, and that PKX clusters in a clade with Tetracapsula bryozoides. Our analyses and those of others, along with phenotypic observations, indicate that salmonids are the primary myxosporean hosts for PKX, that the cryptic spores of PKX in salmonids are the fully formed myxospores as they occur in the fish host, and that PKX represents distinct species that we previously place in the genus Tetracapsula in the family Saccosporidae. The latter 2 taxa were described based on stages from bryozoans, and the myxosporean stage in fish of the type species, T. bryozoides, has not been identified (if it exists). Thus, more complete resolution of the life cycle of both PKX and T. bryozoides, as well as more genetic data, are required to determine the precise relationship of these organisms.

A new Microsporidium sp. (microsporidia) from the musculature of the mountain whitefish Prosopium williamsoni from British Columbia: morphology and phylogeny

Few microsporidia have been reported from whitefish species (subfamily Coregoninae). For the most part, these microsporidia have been incompletely described. In a survey of parasites of mountain whitefish Prosopium williamsoni collected from Kootenay Lake, British Columbia, we encountered an unusual microsporidium infecting the endomysium of the skeletal musculature. Spores were uninucleate, ovoid to pyriform, and were 5.6 (5-7) microm x 3.2 (3-4) microm with 13-16 coils in the polar filament. We describe here this organism as a new species based on its site of development and its relationship among fish microsporidia based on small subunit ribosomal DNA sequence data, i.e., our analysis showed that it is not closely related to other microsporidia for which ribosomal DNA sequence is available thus far.

Taxonomy of Kudoa species (Myxosporea), using a small-subunit ribosomal DNA sequence

Myxosporeans in the genus Kudoa infect the flesh of many marine fishes and often cause unsightly lesions and softening of the flesh texture. We are particularly interested in K. thyrsites because it is associated with soft flesh in the Atlantic salmon (Salmo solar), an important commercial species in Canada. Sequences of the small-subunit (SSU) rDNA (about 1600 base pairs) were obtained from K. miniauriculata, K. amanuensis, and K. poniformis. We aligned these sequences with one obtained from coho salmon (Oncorhynchus kisutch) and designed "Kudoa general" primers (KUD1f and KUD2r). These primers, in combination with other general primers, were then used to obtain the SSU rDNA sequence of K. thyrsites from two host species, Atlantic salmon and tubesnout (Aulorhynchus flavidus), from British Columbia, Canada. Sequence comparisons of these isolates indicated that Kudoa species cluster by geographic location rather than by morphology of spores. The three species from the eastern Pacific were approximately 97% identical, whereas K. amamiensis (from Japan) was about 91% identical with these species. Sequence comparisons of K. thyrsites from Atlantic salmon and tubesnout revealed a difference of only 0.07% between these isolates. Comparison of SSU rDNA sequences from the four Kudoa species and Henneguya salminicolo analyzed in this study with those from other available myxosporean genera (Myxidium and Myxoholus) showed that taxonomic divisions at the order and suborder levels were consistent with classical views of the taxonomy of the Myxosporea. Using specific regions of the SSU rDNA, we also developed a sensitive and specific polymerase chain reaction test for detection of K. thyrsites.