Dominance of Sulfurospirillum in Metagenomes Associated with the Methane Ice Worm (Sirsoe methanicola).

Sirsoe methanicola, commonly known as the methane ice worm, is the only macrofaunal species known to inhabit the Gulf of Mexico methane hydrates. Little is known about this elusive marine polychaete that can colonize rich carbon and energy reserves. Metagenomic analysis of gut contents and worm fragments predicted diverse metabolic capabilities with the ability to utilize a range of nitrogen, sulfur, and organic carbon compounds through microbial taxa affiliated with Campylobacterales, Desulfobacterales, Enterobacterales, SAR324, Alphaproteobacteria, and Mycoplasmatales. Entomoplasmatales and Chitinivibrionales were additionally identified from extracted full-length 16S rRNA sequences, and read analysis identified 196 bacterial families. Overall, the microbial community appeared dominated by uncultured Sulfurospirillum, a taxon previously considered free-living rather than host-associated. Metagenome-assembled genomes (MAGs) classified as uncultured Sulfurospirillum predicted thiosulfate disproportionation and the reduction of tetrathionate, sulfate, sulfide/polysulfide, and nitrate. Microbial amino acid and vitamin B12 biosynthesis genes were identified in multiple MAGs, suggesting nutritional value to the host. Reads assigned to aerobic or anaerobic methanotrophic taxa were rare. IMPORTANCE Methane hydrates represent vast reserves of natural gas with roles in global carbon cycling and climate change. This study provided the first analysis of metagenomes associated with Sirsoe methanicola, the only polychaete species known to colonize methane hydrates. Previously unrecognized participation of Sulfurospirillum in a gut microbiome is provided, and the role of sulfur compound redox reactions within this community is highlighted. The comparative biology of S. methanicola is of general interest given research into the adverse effects of sulfide production in human gut microbiomes. In addition, taxonomic assignments are provided for nearly 200 bacterial families, expanding our knowledge of microbiomes in the deep sea.

Novel Dormancy Mechanism of Castration Resistance in Bone Metastatic Prostate Cancer Organoids.

Advanced prostate cancer (PCa) patients with bone metastases are treated with androgen pathway directed therapy (APDT). However, this treatment invariably fails and the cancer becomes castration resistant. To elucidate resistance mechanisms and to provide a more predictive pre-clinical research platform reflecting tumor heterogeneity, we established organoids from a patient-derived xenograft (PDX) model of bone metastatic prostate cancer, PCSD1. APDT-resistant PDX-derived organoids (PDOs) emerged when cultured without androgen or with the anti-androgen, enzalutamide. Transcriptomics revealed up-regulation of neurogenic and steroidogenic genes and down-regulation of DNA repair, cell cycle, circadian pathways and the severe acute respiratory syndrome (SARS)-CoV-2 host viral entry factors, ACE2 and TMPRSS2. Time course analysis of the cell cycle in live cells revealed that enzalutamide induced a gradual transition into a reversible dormant state as shown here for the first time at the single cell level in the context of multi-cellular, 3D living organoids using the Fucci2BL fluorescent live cell cycle tracker system. We show here a new mechanism of castration resistance in which enzalutamide induced dormancy and novel basal-luminal-like cells in bone metastatic prostate cancer organoids. These PDX organoids can be used to develop therapies targeting dormant APDT-resistant cells and host factors required for SARS-CoV-2 viral entry.

Top-Down Atmospheric Ionization Mass Spectrometry Microscopy Combined With Proteogenomics.

Mass spectrometry-based protein analysis has become an important methodology for proteogenomic mapping by providing evidence for the existence of proteins predicted at the genomic level. However, screening and identification of proteins directly on tissue samples, where histological information is preserved, remain challenging. Here we demonstrate that the ambient ionization source, nanospray desorption electrospray ionization (nanoDESI), interfaced with light microscopy allows for protein profiling directly on animal tissues at the microscopic scale. Peptide fragments for mass spectrometry analysis were obtained directly on ganglia of the medicinal leech (Hirudo medicinalis) without in-gel digestion. We found that a hypothetical protein, which is predicted by the leech genome, is highly expressed on the specialized neural cells that are uniquely found in adult sex segmental ganglia. Via this top-down analysis, a post-translational modification (PTM) of tyrosine sulfation to this neuropeptide was resolved. This three-in-one platform, including mass spectrometry, microscopy, and genome mining, provides an effective way for mappings of proteomes under the lens of a light microscope.

The WNT receptor FZD7 is required for maintenance of the pluripotent state in human embryonic stem cells.

WNT signaling is involved in maintaining stem cells in an undifferentiated state; however, it is often unclear which WNTs and WNT receptors are mediating these activities. Here we examined the role of the WNT receptor FZD7 in maintaining human embryonic stem cells (hESCs) in an undifferentiated and pluripotent state. FZD7 expression is significantly elevated in undifferentiated cells relative to differentiated cell populations, and interfering with its expression or function, either by short hairpin RNA-mediated knockdown or with a fragment antigen binding (Fab) molecule directed against FZD7, disrupts the pluripotent state of hESCs. The FZD7-specific Fab blocks signaling by Wnt3a protein by down-regulating FZD7 protein levels, suggesting that FZD7 transduces Wnt signals to activate Wnt/ß-catenin signaling. These results demonstrate that FZD7 encodes a regulator of the pluripotent state and that hESCs require endogenous WNT/ß-catenin signaling through FZD7 to maintain an undifferentiated phenotype.

Vespucci: a system for building annotated databases of nascent transcripts.

Global run-on sequencing (GRO-seq) is a recent addition to the series of high-throughput sequencing methods that enables new insights into transcriptional dynamics within a cell. However, GRO-sequencing presents new algorithmic challenges, as existing analysis platforms for ChIP-seq and RNA-seq do not address the unique problem of identifying transcriptional units de novo from short reads located all across the genome. Here, we present a novel algorithm for de novo transcript identification from GRO-sequencing data, along with a system that determines transcript regions, stores them in a relational database and associates them with known reference annotations. We use this method to analyze GRO-sequencing data from primary mouse macrophages and derive novel quantitative insights into the extent and characteristics of non-coding transcription in mammalian cells. In doing so, we demonstrate that Vespucci expands existing annotations for mRNAs and lincRNAs by defining the primary transcript beyond the polyadenylation site. In addition, Vespucci generates assemblies for un-annotated non-coding RNAs such as those transcribed from enhancer-like elements. Vespucci thereby provides a robust system for defining, storing and analyzing diverse classes of primary RNA transcripts that are of increasing biological interest.

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma.

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹8), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹°). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Hypothesis-independent pathway analysis implicates GABA and acetyl-CoA metabolism in primary open-angle glaucoma and normal-pressure glaucoma.

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. Using genome-wide association single-nucleotide polymorphism data from the Glaucoma Genes and Environment study and National Eye Institute Glaucoma Human Genetics Collaboration comprising 3,108 cases and 3,430 controls, we assessed biologic pathways as annotated in the KEGG database for association with risk of POAG. After correction for genic overlap among pathways, we found 4 pathways, butanoate metabolism (hsa00650), hematopoietic cell lineage (hsa04640), lysine degradation (hsa00310) and basal transcription factors (hsa03022) related to POAG with permuted p < 0.001. In addition, the human leukocyte antigen (HLA) gene family was significantly associated with POAG (p < 0.001). In the POAG subset with normal-pressure glaucoma (NPG), the butanoate metabolism pathway was also significantly associated (p < 0.001) as well as the MAPK and Hedgehog signaling pathways (hsa04010 and hsa04340), glycosaminoglycan biosynthesis-heparan sulfate pathway (hsa00534) and the phenylalanine, tyrosine and tryptophan biosynthesis pathway (hsa0400). The butanoate metabolism pathway overall, and specifically the aspects of the pathway that contribute to GABA and acetyl-CoA metabolism, was the only pathway significantly associated with both POAG and NPG. Collectively these results implicate GABA and acetyl-CoA metabolism in glaucoma pathogenesis, and suggest new potential therapeutic targets.

Association of CAV1/CAV2 genomic variants with primary open-angle glaucoma overall and by gender and pattern of visual field loss.

PURPOSE: The CAV1/CAV2 (caveolin 1 and caveolin 2) genomic region previously was associated with primary open-angle glaucoma (POAG), although replication among independent studies has been variable. The aim of this study was to assess the association between CAV1/CAV2 single nucleotide polymorphisms (SNPs) and POAG in a large case-control dataset and to explore associations by gender and pattern of visual field (VF) loss further. DESIGN: Case-control study. PARTICIPANTS: We analyzed 2 large POAG data sets: the Glaucoma Genes and Environment (GLAUGEN) study (976 cases, 1140 controls) and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium (2132 cases, 2290 controls). METHODS: We studied the association between 70 SNPs located within the CAV1/CAV2 genomic region in the GLAUGEN and NEIGHBOR studies, both genotyped on the Illumina Human 660WQuadv1C BeadChip array and imputed with the Markov Chain Haplotyping algorithm using the HapMap 3 reference panel. We used logistic regression models of POAG in the overall population and separated by gender, as well as by POAG subtypes defined by type of VF defect (peripheral or paracentral). Results from GLAUGEN and NEIGHBOR were meta-analyzed, and a Bonferroni-corrected significance level of 7.7 × 10(-4) was used to account for multiple comparisons. MAIN OUTCOME MEASURES: Overall POAG, overall POAG by gender, and POAG subtypes defined by pattern of early VF loss. RESULTS: We found significant associations between 10 CAV1/CAV2 SNPs and POAG (top SNP, rs4236601; pooled P = 2.61 × 10(-7)). Of these, 9 were significant only in women (top SNP, rs4236601; pooled P = 1.59 × 10(-5)). Five of the 10 CAV1/CAV2 SNPs were associated with POAG with early paracentral VF (top SNP, rs17588172; pooled P = 1.07 × 10(-4)), and none of the 10 were associated with POAG with peripheral VF loss only or POAG among men. CONCLUSIONS: CAV1/CAV2 SNPs were associated significantly with POAG overall, particularly among women. Furthermore, we found an association between CAV1/CAV2 SNPs and POAG with paracentral VF defects. These data support a role for caveolin 1, caveolin 2, or both in POAG and suggest that the caveolins particularly may affect POAG pathogenesis in women and in patients with early paracentral VF defects.

Illumina mate-paired DNA sequencing-library preparation using Cre-Lox recombination.

Standard Illumina mate-paired libraries are constructed from 3- to 5-kb DNA fragments by a blunt-end circularization. Sequencing reads that pass through the junction of the two joined ends of a 3-5-kb DNA fragment are not easy to identify and pose problems during mapping and de novo assembly. Longer read lengths increase the possibility that a read will cross the junction. To solve this problem, we developed a mate-paired protocol for use with Illumina sequencing technology that uses Cre-Lox recombination instead of blunt end circularization. In this method, a LoxP sequence is incorporated at the junction site. This sequence allows screening reads for junctions without using a reference genome. Junction reads can be trimmed or split at the junction. Moreover, the location of the LoxP sequence in the reads distinguishes mate-paired reads from spurious paired-end reads. We tested this new method by preparing and sequencing a mate-paired library with an insert size of 3 kb from Saccharomyces cerevisiae. We present an analysis of the library quality statistics and a new bio-informatics tool called DeLoxer that can be used to analyze an IlluminaCre-Lox mate-paired data set. We also demonstrate how the resulting data significantly improves a de novo assembly of the S. cerevisiae genome.

Paving the future: finding suitable ISMB venues.

The International Society for Computational Biology, ISCB, organizes the largest event in the field of computational biology and bioinformatics, namely the annual international conference on Intelligent Systems for Molecular Biology, the ISMB. This year at ISMB 2012 in Long Beach, ISCB celebrated the 20th anniversary of its flagship meeting. ISCB is a young, lean and efficient society that aspires to make a significant impact with only limited resources. Many constraints make the choice of venues for ISMB a tough challenge. Here, we describe those challenges and invite the contribution of ideas for solutions.

Estrogen pathway polymorphisms in relation to primary open angle glaucoma: an analysis accounting for gender from the United States.

PURPOSE: Circulating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender. METHODS: We included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP <22 mmHg (normal pressure glaucoma [NPG]) at diagnosis. We conducted these analyses for each gender separately and then jointly in men and women. RESULTS: Among women, the estrogen SNP pathway was associated with POAG overall (permuted p=0.006) and HPG (permuted p<0.001) but not NPG (permuted p=0.09). Interestingly, there was no relation between the estrogen SNP pathway and POAG when men were considered alone (permuted p>0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01). CONCLUSIONS: The estrogen SNP pathway was associated with POAG among women.

The medicinal leech genome encodes 21 innexin genes: different combinations are expressed by identified central neurons.

Gap junctional proteins are important components of signaling pathways required for the development and ongoing functions of all animal tissues, particularly the nervous system, where they function in the intracellular and extracellular exchange of small signaling factors and ions. In animals whose genomes have been sufficiently sequenced, large families of these proteins, connexins, pannexins, and innexins, have been found, with 25 innexins in the nematode Caenorhabditis elegans Starich et al. (Cell Commun Adhes 8: 311-314, 2001) and at least 37 connexins in the zebrafish Danio rerio Cruciani and Mikalsen (Biol Chem 388:253-264, 2009). Having recently sequenced the medicinal leech Hirudo verbana genome, we now report the presence of 21 innexin genes in this species, nine more than we had previously reported from the analysis of an EST-derived transcriptomic database Dykes and Macagno (Dev Genes Evol 216: 185-97, 2006); Macagno et al. (BMC Genomics 25:407, 2010). Gene structure analyses show that, depending on the leech innexin gene, they can contain from 0 to 6 introns, with closely related paralogs showing the same number of introns. Phylogenetic trees comparing Hirudo to another distantly related leech species, Helobdella robusta, shows a high degree of orthology, whereas comparison to other annelids shows a relatively low level. Comparisons with other Lophotrochozoans, Ecdyzozoans and with vertebrate pannexins suggest a low number (one to two) of ancestral innexin/pannexins at the protostome/deuterostome split. Whole-mount in situ hybridization for individual genes in early embryos shows that ∼50% of the expressed innexins are detectable in multiple tissues. Expression analyses using quantitative PCR show that ∼70% of the Hirudo innexins are expressed in the nervous system, with most of these detected in early development. Finally, quantitative PCR analysis of several identified adult neurons detects the presence of different combinations of innexin genes, a property that may underlie the participation of these neurons in different adult coupling circuits.

Biomedical impact of splicing mutations revealed through exome sequencing.

Splicing is a cellular mechanism, which dictates eukaryotic gene expression by removing the noncoding introns and ligating the coding exons in the form of a messenger RNA molecule. Alternative splicing (AS) adds a major level of complexity to this mechanism and thus to the regulation of gene expression. This widespread cellular phenomenon generates multiple messenger RNA isoforms from a single gene, by utilizing alternative splice sites and promoting different exon-intron inclusions and exclusions. AS greatly increases the coding potential of eukaryotic genomes and hence contributes to the diversity of eukaryotic proteomes. Mutations that lead to disruptions of either constitutive splicing or AS cause several diseases, among which are myotonic dystrophy and cystic fibrosis. Aberrant splicing is also well established in cancer states. Identification of rare novel mutations associated with splice-site recognition, and splicing regulation in general, could provide further insight into genetic mechanisms of rare diseases. Here, disease relevance of aberrant splicing is reviewed, and the new methodological approach of starting from disease phenotype, employing exome sequencing and identifying rare mutations affecting splicing regulation is described. Exome sequencing has emerged as a reliable method for finding sequence variations associated with various disease states. To date, genetic studies using exome sequencing to find disease-causing mutations have focused on the discovery of nonsynonymous single nucleotide polymorphisms that alter amino acids or introduce early stop codons, or on the use of exome sequencing as a means to genotype known single nucleotide polymorphisms. The involvement of splicing mutations in inherited diseases has received little attention and thus likely occurs more frequently than currently estimated. Studies of exome sequencing followed by molecular and bioinformatic analyses have great potential to reveal the high impact of splicing mutations underlying human disease.

Mild iron overload in an African American man with SLC40A1 D270V.

We report on a 46-year-old black man who resided in Alabama with normal transferrin saturation, mild hyperferritinemia, chronic hepatitis C, and 3+ iron in hepatocytes and Kupffer cells. Exome sequencing revealed heterozygosity for SLC40A1 D270V (exon 7, c.809A→T), a mutation previously reported only in 1 black patient with iron overload who resided in the Republic of South Africa. The present patient was also heterozygous for: heme transporter FLVCR1 novel allele P542S (exon 10, 1624C→T); FLVCR1 T544M (rs3207090); hemopexin (HPX) R371W (rs75307540); ferritin scavenger receptor (SCARA5) R471H (rs61737287); and transferrin receptor (TFRC) G420S (rs41295879). He had no HFE, TFR2,HJV, or HAMP mutations. D270V was not detected in 19 other African Americans with iron overload who resided in Alabama. The allele frequency of SLC40A1 D270V in 258 African American adults who participated in a health appraisal clinic was 0.0019 (95% confidence interval 0-0.0057). D270V could explain 'classical' ferroportin hemochromatosis phenotypes in some African Americans.

An RNA map predicting Nova-dependent splicing regulation.

Nova proteins are a neuron-specific alternative splicing factors. We have combined bioinformatics, biochemistry and genetics to derive an RNA map describing the rules by which Nova proteins regulate alternative splicing. This map revealed that the position of Nova binding sites (YCAY clusters) in a pre-messenger RNA determines the outcome of splicing. The map correctly predicted Nova's effect to inhibit or enhance exon inclusion, which led us to examine the relationship between the map and Nova's mechanism of action. Nova binding to an exonic YCAY cluster changed the protein complexes assembled on pre-mRNA, blocking U1 snRNP (small nuclear ribonucleoprotein) binding and exon inclusion, whereas Nova binding to an intronic YCAY cluster enhanced spliceosome assembly and exon inclusion. Assays of splicing intermediates of Nova-regulated transcripts in mouse brain revealed that Nova preferentially regulates removal of introns harbouring (or closest to) YCAY clusters. These results define a genome-wide map relating the position of a cis-acting element to its regulation by an RNA binding protein, namely that Nova binding to YCAY clusters results in a local and asymmetric action to regulate spliceosome assembly and alternative splicing in neurons.

Automated querying and identification of novel peptides using MALDI mass spectrometric imaging.

MSI is a molecular imaging technique that allows for the generation of topographic 2D maps for various endogenous and some exogenous molecules without prior specification of the molecule. In this paper, we start with the premise that a region of interest (ROI) is given to us based on preselected morphological criteria. Given an ROI, we develop a pipeline, first to determine mass values with distinct expression signatures, localized to the ROI, and second to identify the peptides corresponding to these mass values. To identify spatially differentiated masses, we implement a statistic that allows us to estimate, for each spectral peak, the probability that it is over- or under-expressed within the ROI versus outside. To identify peptides corresponding to these masses, we apply LC-MS/MS to fragment endogenous (nonprotease digested) peptides. A novel pipeline based on constructing sequence tags de novo from both original and decharged spectra and a subsequent database search is used to identify peptides. As the MSI signal and the identified peptide are only related by a single mass value, we isolate the corresponding transcript and perform a second validation via in situ hybridization of the transcript. We tested our approach, MSI-Query, on a number of ROIs in the medicinal leech, Hirudo medicinalis, including the central nervous system (CNS). The Hirudo CNS is capable of regenerating itself after injury, thus forming an important model system for neuropeptide identification. The pipeline helps identify a number of novel peptides. Specifically, we identify a gene that we name HmIF4, which is a member of the intermediate filament family involved in neural development and a second novel, uncharacterized peptide. A third peptide, derived from the histone H2B, is also identified, in agreement with the previously suggested role of histone H2B in axon targeting.

A Machine Learning Method to Identify Genetic Variants Potentially Associated With Alzheimer's Disease.

There is hope that genomic information will assist prediction, treatment, and understanding of Alzheimer's disease (AD). Here, using exome data from ∼10,000 individuals, we explore machine learning neural network (NN) methods to estimate the impact of SNPs (i.e., genetic variants) on AD risk. We develop an NN-based method (netSNP) that identifies hundreds of novel potentially protective or at-risk AD-associated SNPs (along with an effect measure); the majority with frequency under 0.01. For case individuals, the number of "protective" (or "at-risk") netSNP-identified SNPs in their genome correlates positively (or inversely) with their age of AD diagnosis and inversely (or positively) with autopsy neuropathology. The effect measure increases correlations. Simulations suggest our results are not due to genetic linkage, overfitting, or bias introduced by netSNP. These findings suggest that netSNP can identify SNPs associated with AD pathophysiology that may assist with the diagnosis and mechanistic understanding of the disease.

Dynamic transcriptome landscape in the song nucleus HVC between juvenile and adult zebra finches.

Male juvenile zebra finches learn to sing by imitating songs of adult males early in life. The development of the song control circuit and song learning and maturation are highly intertwined processes, involving gene expression, neurogenesis, circuit formation, synaptic modification, and sensory-motor learning. To better understand the genetic and genomic mechanisms underlying these events, we used RNA-Seq to examine genome-wide transcriptomes in the song control nucleus HVC of male juvenile (45 d) and adult (100 d) zebra finches. We report that gene groups related to axon guidance, RNA processing, lipid metabolism, and mitochondrial functions show enriched expression in juvenile HVC compared to the rest of the brain. As juveniles mature into adulthood, massive gene expression changes occur. Expression of genes related to amino acid metabolism, cell cycle, and mitochondrial function is reduced, accompanied by increased and enriched expression of genes with synaptic functions, including genes related to G-protein signaling, neurotransmitter receptors, transport of small molecules, and potassium channels. Unexpectedly, a group of genes with immune system functions is also developmentally regulated, suggesting potential roles in the development and functions of HVC. These data will serve as a rich resource for investigations into the development and function of a neural circuit that controls vocal behavior.

Construction of a medicinal leech transcriptome database and its application to the identification of leech homologs of neural and innate immune genes.

BACKGROUND: The medicinal leech, Hirudo medicinalis, is an important model system for the study of nervous system structure, function, development, regeneration and repair. It is also a unique species in being presently approved for use in medical procedures, such as clearing of pooled blood following certain surgical procedures. It is a current, and potentially also future, source of medically useful molecular factors, such as anticoagulants and antibacterial peptides, which may have evolved as a result of its parasitizing large mammals, including humans. Despite the broad focus of research on this system, little has been done at the genomic or transcriptomic levels and there is a paucity of openly available sequence data. To begin to address this problem, we constructed whole embryo and adult central nervous system (CNS) EST libraries and created a clustered sequence database of the Hirudo transcriptome that is available to the scientific community. RESULTS: A total of approximately 133,000 EST clones from two directionally-cloned cDNA libraries, one constructed from mRNA derived from whole embryos at several developmental stages and the other from adult CNS cords, were sequenced in one or both directions by three different groups: Genoscope (French National Sequencing Center), the University of Iowa Sequencing Facility and the DOE Joint Genome Institute. These were assembled using the phrap software package into 31,232 unique contigs and singletons, with an average length of 827 nt. The assembled transcripts were then translated in all six frames and compared to proteins in NCBI's non-redundant (NR) and to the Gene Ontology (GO) protein sequence databases, resulting in 15,565 matches to 11,236 proteins in NR and 13,935 matches to 8,073 proteins in GO. Searching the database for transcripts of genes homologous to those thought to be involved in the innate immune responses of vertebrates and other invertebrates yielded a set of nearly one hundred evolutionarily conserved sequences, representing all known pathways involved in these important functions. CONCLUSIONS: The sequences obtained for Hirudo transcripts represent the first major database of genes expressed in this important model system. Comparison of translated open reading frames (ORFs) with the other openly available leech datasets, the genome and transcriptome of Helobdella robusta, shows an average identity at the amino acid level of 58% in matched sequences. Interestingly, comparison with other available Lophotrochozoans shows similar high levels of amino acid identity, where sequences match, for example, 64% with Capitella capitata (a polychaete) and 56% with Aplysia californica (a mollusk), as well as 58% with Schistosoma mansoni (a platyhelminth). Phylogenetic comparisons of putative Hirudo innate immune response genes present within the Hirudo transcriptome database herein described show a strong resemblance to the corresponding mammalian genes, indicating that this important physiological response may have older origins than what has been previously proposed.