The airway microbiome and disease.

Although traditionally thought to be sterile, accumulating evidence now supports the concept that our airways harbor a microbiome. Thus far, studies have focused upon characterizing the bacterial constituents of the airway microbiome in both healthy and diseased lungs, but what perhaps provides the greatest impetus for the exploration of the airway microbiome is that different bacterial phyla appear to dominate diseased as compared with healthy lungs. As yet, there is very limited evidence supporting a functional role for the airway microbiome, but continued research in this direction is likely to provide such evidence, particularly considering the progress that has been made in understanding host-microbe mutualism in the intestinal tract. In this review, we highlight the major advances that have been made discovering and describing the airway microbiome, discuss the experimental evidence that supports a functional role for the microbiome in health and disease, and propose how this emerging field is going to impact clinical practice.

Antigen Receptor Signaling to NF-kappa B via CARMA1, BCL10, and MALT1

The signaling pathway controlling antigen receptor-induced regulation of the transcription factor NF-kappa B plays a key role in lymphocyte activation and development and the generation of lymphomas. Work of the past decade has led to dramatic progress in the identification and characterization of new players in the pathway. Moreover, novel enzymatic activities relevant for this pathway have been discovered, which represent interesting drug targets for immuno-suppression or lymphoma treatment. Here, we summarize these findings and give an outlook on interesting open issues that need to be addressed in the future.

Biochemical characterization of recombinant human telomerase

Résumé Les télomères sont les structures ADN-protéines des extrémités des chromosomes des eucaryotes. L'ADN télomérique est constitué de courtes séquences répétitives. L'intégrité des télomères est essentielle pour protéger les extrémités des chromosomes contre les systèmes de dégradations et pour les distinguer des cassures de l'ADN double brin. Parce que la machinerie de la réplication de l'ADN n'est pas capable de répliquer l'extrémité des chromosomes, les télomères raccourcissent au fur et à mesure des cycles de réplication. Dès que les télomères atteignent une longueur critique, leur structure protectrice est perdue. Cela induit un signal de dommage de l'ADN et l'arrêt du cycle cellulaire. Pour contrebalancer le raccourcissement des télomères, les cellules qui s'auto régénèrent, dont les cellules de la moelle osseuse, les lymphocytes activés et 80-90% des cellules cancéreuses, expriment la télomérase. C'est une ribonucléoprotéine qui a la capacité de synthétiser des séquences télomériques par transcription inverse d'une courte séquence contenue dans sa propre sous-unité ARN avec laquelle elle est associée. La télomérase humaine est une enzyme processive au niveau de l'addition des nucléotides et aussi des répétitions télomériques. La télomérase de levure et la télomérase humaine sont toutes deux dimériques et il a été montré que la télomérase humaine recombinante contient deux ARN qui coopèrent pour fonctionner ainsi que deux sous-unités catalytiques. Cependant, il n'a pas encore été montré quel est le rôle de la dimérisation dans l'activité de la télomérase. Afin d'élucider ce rôle, nous avons exprimé, reconstitué et purifié la télomérase humaine dimérique recombinante. Et pour étudier l'effet d'ARN mutants sur l'activité de la télomérase, nous avons développé une méthode pour reconstituer et enrichir en hétérodimères de télomérase. Les hétérodimères contiennent une sous-unité ARN sauvage et une sous-unité ARN mutée au niveau de la séquence de la matrice. Sur l'ARN muté nous avons introduit une étiquette aptamer ARN-S1 puis nous avons purifié la télomérase via l'etiquette Si. Nous avons montré que la dimérisation est essentielle pour l'activité de la télomérase. Nos données indiquent que chaque télomérase du dimère allonge leur substrat, l'ADN télomérique, indépendamment l'une de l'autre à chaque cycle d'élongation mais que l'addition itérative de répétitions télomériques nécessite une coopération entre les deux télomérases du dimère. Nous proposons donc un modèle dans lequel les deux télomérases du dimères se lient et allongent deux substrats télomères et que pendant l'élongation processive les deux enzymes subissent un changement de conformation de manière coordonnée, ce changement va permettre le repositionnement des substrats pour d'autres cycles d'additions de répétitions télomériques. Dyskeratosis congenita est une maladie mortelle due majoritairement au disfonctionnement de la moelle osseuse. Dans la forme autosomale de la maladie, l'ARN de la télomérase contient des mutations. En utilisant notre système de reconstitution, nous avons montré que ces ARN mutés, qui ont perdu leur activité enzymatique dans le cas d'un homodimère de mutants, sont dominant négatifs quand ils sont présents dans les hétérodimères sauvage/mutant. Cet effet trans-dominant négatif pourrait contribuer à la progression de la maladie. Abstract Telomeres are protein-DNA structures at the ends of linear eukaryotic chromosomes. The telomeric DNA consists of tandemly repeated sequences. Telomeric integrity is essential to protect chromosomal ends from nucleolytic degradation and to prevent their recognition as DNA double strand breaks. Due to the inability of the conventional DNA replication machinery to replicate terminal DNA stretches, telomeres shorten with continuous rounds of DNA replication. As soon as telomeres reach a critical length, their protective structure is lost and the deprotected telomeres will induce a DNA damage response leading to cell cycle arrest. To counteract telomere shortening, self-renewing cells, including bone marrow cells, activated lymphocytes and 80-90% of cancer cells express the cellular reverse transcriptase telomerase, which has the capacity to synthesize telomeric repeats by reverse transcription of a short template sequence encoded by its stably associated RNA subunit. Human telomerase is a processive enzyme for nucleotide as well as repeat addition. Both yeast and human telomerase are dimeric enzymes and recombinant human telomerase has been shown to contain two functionally cooperating RNAs and most probably also two protein subunits. However, it has remained unclear how dimerization may contribute to telomerase activity. To study the role of dimerization, we expressed, reconstituted and purified recombinant human telomerase. We also developed a new method to reconstitute and enrich for telomerase heterodimers containing wild-type (wt) and mutant telomerase RNA subunits. To this end we introduced an S1-RNA-aptamer tag into telomerase RNA and purified telomerase reconstituted with a mixture of untagged and tagged RNA via the S1-tag. Using this experimental system, we introduced template mutations in the tagged RNA subunit and examined the effect of mutant RNAs on wt telomerase activity in wt/mutant heterodimers. We obtained evidence that dimerization is essential for telomerase activity. Our data indicate that the two subunits elongate telomere substrates independently of each other during single rounds of elongation, but that iterative addition of telomeric repeats requires cooperation between the two subunits. We suggest a model, in which dimeric telomerases bind and elongate two telomere substrates and that the two subunits undergo coordinated conformational changes during processive elongation that enable repositioning the substrates for subsequent rounds of repeat addition. Dyskeratosis congenita is a multisystemic disease with bone marrow failure as the major cause of death. The autosomal form of this disease was found to harbor mutations in the telomerase RNA. Using our reconstitution system, we tested whether mutant dyskeratosis telomerase RNAs behaved in a dominant negative manner. We observed that dyskeratosis telomerase RNA mutants, which lacked enzymatic activity were dominant negative, when present in wt/ mutant heterodimers. The transdominant negative effect of these mutants may contribute to disease progression.

Significant genetic and phenotypic changes arising from clonal growth of a single spore of an arbuscular mycorrhizal fungus over multiple generations.

Arbuscular mycorrhizal fungi (AMF) are highly successful plant symbionts. They reproduce clonally producing multinucleate spores. It has been suggested that some AMF harbor genetically different nuclei. However, recent advances in sequencing the Glomus irregulare genome have indicated very low within-fungus polymorphism. We tested the null hypothesis that, with no genetic differences among nuclei, no significant genetic or phenotypic variation would occur among clonal single spore lines generated from one initial AMF spore. Furthermore, no additional variation would be expected in the following generations of single spore lines. Genetic diversity contained in one initial spore repeatedly gave rise to genetically different variants of the fungus with novel phenotypes. The genetic changes represented quantitative changes in allele frequencies, most probably as a result of changes in the frequency of genetic variation partitioned on different nuclei. The genetic and phenotypic variation is remarkable, given that it arose repeatedly from one clonal individual. Our results highlight the dynamic nature of AMF genetics. Even though within-fungus genetic variation is low, some is probably partitioned among nuclei and potentially causes changes in the phenotype. Our results are important for understanding AMF genetics, as well as for researchers and biotechnologists hoping to use AMF genetic diversity for the improvement of AMF inoculum.

Interactions between feather-degrading bacteria and an avian host : zebra finches (Taeniopygia guttata)

Abstract: Birds harbor a variety of bacteria on their plumage, some of which can degrade feathers in vitro. Whether these keratinolytic bacteria are active on live birds and can effect feather degradation on birds is debatable. The effect of such bacteria on the body condition and behavior of birds, is unknown. Using a community of feather-degrading bacteria (EB), we investigate the interaction between the activity and load of such bacteria, on the morphology, body condition, and behavior of zebra finches (Taeniopygia guttata). In Chapter 2, we find that the elevated loads of such microbes lead to a reduction in the expression of morphological traits, such as male bill color (a sexually selected trait) and uropygial gland volume, without reducing body mass, or evoking a cellular immune response. We also suggest the presence of a carotenoid based defense response in hosts, to such elevated loads of microbes and document a sex-based difference in the source of carotenoids used for such a response. In Chapter 3, we investigated the effect of EB loads on male mate choice of zebra finches, wherein male choice of females with elevated and un-altered bacterial loads, varied with male size. We found that larger males preferred females with higher bacterial load and smaller males preferred females with lower bacterial load. Chapter 4 demonstrates that the presence of melanin in feathers reduces the growth and activity of the community of feather-degrading bacteria (EB) and that the EB community can effect feather degradation in humid conditions, without broth. Additional results also demonstrate that the EB community consists of bacteria that can attach themselves to feathers on live birds and those that can live freely on avian plumage. Finally, chapter 5 demonstrates that the self-maintenance, social and sexual behaviors of birds are implicated in the infection and horizontal transmission of bacteria. It also suggests a linked oral - faecal - genital mode of transmission of pathogens in birds. These results demonstrate that differential loads of normal flora of vertebrate hosts can effect changes in their morphology and behavior. They also shed light on the role of feather-degrading bacteria in the evolution of melanin polymorphism in birds and suggest that bacteria can be active on live birds. This thesis also highlights the importance of social and, sexual behaviors of birds, in epidemiology. Résumé: Les Oiseaux ont dans leur plumage diverses bactéries dont certaines dégradent les plumes in vitro, néanmoins. Il n'est pas clair, au vu de précédentes études, si ces bactéries kératinolytiques sont actives sur des oiseaux vivants, et si celles-ci dégradent effectivement le plumage de leur hôte, L'effet de ces bactéries sur la condition corporelle ainsi que le comportement des oiseaux n'est pas connu. A l'aide d'une communauté de bactéries dégradant les plumes (EB), non pathogènes, nous examinons les interactions entre l'activité et la charge bactérienne sur la morphologie, la condition corporelle et le comportement du diamant mandarins (Taeniopygia guttata). Dans le chapitre 2, nous montrons qu'une charge élevée de ces microbes mène à une réduction de l'expression de certains traits morphologiques, tels que la couleur du bec chez le mâle (un trait soumis à sélection sexuelle), ainsi que le volume de la glande uropygienne, sans qu'il y ait une réduction de la masse corporelle, ni déclenchement d'une réponse immune cellulaire. Nos données suggèrent la présence d'une défense chez l'hôte à des charges élevées de bactéries basée sur la présence de caroténoïdes. Nous montrons, de plus une différence liée au sexe dans la source des caroténoïdes utilisé pour cette réponse. Dans le chapitre 3 nous examinons l'influence de la charge bactérienne EB sur le choix des mâles chez le diamant mandarins. Des femelles avec une charge bactérienne normale et augmentée sont choisies par les mâles et ce choix varie avec la taille des mâles. Nous avons mis en évidence que les grands mâles préfèrent les femelles avec une charge bactérienne plus élevée. Les petits mâles préfèrent les femelles avec une charge bactérienne réduite. Le chapitre 4 démontre que la présence de mélanine dans les plumes réduit la croissance et l'activité de la communauté de bactéries dégradant le plumage (EB), et que cette communauté EB peut dégrader les plumes dans des conditions humides, sans milieu de culture liquide. De plus nous montrons que cette communauté consiste en des bactéries qui peuvent s'attacher sur les plumes d'oiseaux vivants ainsi que des bactéries libres. Pour finir nous montrons dans le chapitre 5 que la maintenance corporelle, l'interaction sociale et le comportement sexuel de ces oiseaux sont impliqués dans l'infection et la transmission horizontale de ces bactéries. Nos données suggèrent une transmission orale-fécale-génitale des pathogènes chez les oiseaux. Ces résultats montrent que des charges différentes de la flore bactérienne habituelle et non pathogène de vertébrés peuvent affecter leur morphologie et leur comportement. Ils éclaircissent également le rôle des bactéries dégradant les plumes dans l'évolution du polymorphisme mélanique chez les oiseaux et suggèrent que ces bactéries peuvent être actives sur des oiseaux vivants. Cette thèse souligne également l'importance du comportement social et sexuel des oiseaux dans l'épidémiologie.

Multifocal epithelial tumors and field cancerization: stroma as a primary determinant.

It is increasingly evident that cancer results from altered organ homeostasis rather than from deregulated control of single cells or groups of cells. This applies especially to epithelial cancer, the most common form of human solid tumors and a major cause of cancer lethality. In the vast majority of cases, in situ epithelial cancer lesions do not progress into malignancy, even if they harbor many of the genetic changes found in invasive and metastatic tumors. While changes in tumor stroma are frequently viewed as secondary to changes in the epithelium, recent evidence indicates that they can play a primary role in both cancer progression and initiation. These processes may explain the phenomenon of field cancerization, i.e., the occurrence of multifocal and recurrent epithelial tumors that are preceded by and associated with widespread changes of surrounding tissue or organ "fields."

Levels of PSP toxins in bivalves exposed to natural blooms of Alexandrium minutum in Catalan harbours

The development, validation, comparison and evaluation of analytical methods for marine toxins rely on the availability of toxic material. Within the project JACUMAR PSP, our interest is mainly focused on autochthonous bivalve species with the toxic profile of Alexandrium minutum, since this is the principal species involved regionally in PSP outbreaks. Mussels and oysters were exposed during few days in the harbor of Vilanova i la Geltrú, to blooms reaching a maximum A. minutum concentration of 200,000 cells L-1 in 2008, and 40,000 and 800,000 cells L-1, in 2009. Mussels, oysters and clams were exposed to one bloom of 22,000 cells L-1 in the harbor of Cambrils in 2009. In all situations higher toxic levels analyzed by HPLC-FD with postcolumn oxidation were observed in mussels (i.e. 1,200-2,500 μg eq. STX kg-1) than in oysters (i.e. 60-800 μg eq. STX kg-1) exposed to the same bloom. Blooms with higher concentrations of A. minutum did not correspond to higher levels of PSP toxins in bivalves. These differences may be explained by differences in A. minutum population dynamics, toxin production or in the physiological state or behaviour of shellfish. These results confirm that mussels concentrate more PSP toxins from A. minutum than oysters and clams.

On the cellular site of two-pore channel TPC1 action in the Poaceae.

The slow vacuolar (SV) channel has been characterized in different dicots by patch-clamp recordings. This channel represents the major cation conductance of the largest organelle in most plant cells. Studies with the tpc1-2 mutant of the model dicot plant Arabidopsis thaliana identified the SV channel as the product of the TPC1 gene. By contrast, research on rice and wheat TPC1 suggested that the monocot gene encodes a plasma membrane calcium-permeable channel. To explore the site of action of grass TPC1 channels, we expressed OsTPC1 from rice (Oryza sativa) and TaTPC1 from wheat (Triticum aestivum) in the background of the Arabidopsis tpc1-2 mutant. Cross-species tpc1 complementation and patch-clamping of vacuoles using Arabidopsis and rice tpc1 null mutants documented that both monocot TPC1 genes were capable of rescuing the SV channel deficit. Vacuoles from wild-type rice but not the tpc1 loss-of-function mutant harbor SV channels exhibiting the hallmark properties of dicot TPC1/SV channels. When expressed in human embryonic kidney (HEK293) cells OsTPC1 was targeted to Lysotracker-Red-positive organelles. The finding that the rice TPC1, just like those from the model plant Arabidopsis and even animal cells, is localized and active in lyso-vacuolar membranes associates this cation channel species with endomembrane function.

Carotid Dose Spring in Definitive Irradiation of T1N0 Squamous Cell Laryngeal Carcinoma Using Helical Tomotherapy

Purpose/Objective(s): To implement a carotid dose sparing protocol using helical Tomotherapy in T1N0 squamous cell laryngeal carcinoma.Materials/Methods: Between July and August 2010, 7 men with stage T1N0 laryngeal carcinoma were included in this study. Age ranged from 47 - 74 years. Staging included endoscopic examination, CT-scan and MRI when indicated. Planned irradiation dose was 70 Gy in 35 fractions over 7 weeks. A simple treatment planning algorithm for carotid sparing was used: maximum point dose to the carotids 35 Gy, to the spinal cord 30 Gy, and 100% PTV volume to be covered with 95% of the prescribed dose. Carotid volume of interest extended to 1 cm above and below of the PTV. Doses to the carotid arteries, to the critical organs, and to the planned target volume (PTV) with our standard laryngeal irradiation protocol was compared. Daily megavoltage scans were obtained before each fraction. When necessary, the Planned Adaptive software (TomoTherapy Inc., Madison, WI) was used to evaluatethe need for a re-planning, which has never been indicated. Dose data were extracted using the VelocityAI software (Atlanta, GA), and data normalization and dose-volume histogram (DVH) interpolation were realized using the Igor Pro software (Portland, OR).Results:A significant (p\0.05) carotid dose sparing compared to our standard protocol with an average maximum point dose of 38.3 Gy (standard deviation [SD] 4.05 Gy), average mean dose of 18.59 Gy (SD 0.83 Gy) was achieved. In all patients, 95% of the carotid volume received less than 28.4 Gy (SD 0.98 Gy). The average maximum point dose to the spinal cord was 25.8 Gy (SD 3.24 Gy). PTV was fully covered with more than 95% of the prescribed dose for all patients with an average maximum point dose of 74.1 Gy and the absolute maximum dose in a single patient of 75.2 Gy. To date, the clinical outcomes have been excellent. Three patients (42%) developed stage 1 mucositis that was conservatively managed, and all the patients presented a mild to moderate dysphonia. All adverse effects resolved spontaneously in the month following the end of treatment. Early local control rate is 100% considering a 4 - 5 months post treatment follow-up.Conclusions: Helical Tomotherapy allows a clinically significant decrease of carotid irradiation dose compared to standard irradiation protocols with an acceptable spinal cord dose tradeoff. Moreover, this technique allows the PTV to be homogenously covered with a curative irradiation dose. Daily control imaging brings added security margins especially when working with high dose gradients. Further investigations and follow-up are underway to better evaluate the late clinical outcomes especially the local control rate, late laryngeal and vascular toxicity, and expected potential impact on cerebrovascular events.

Gènere i identitats en la construcció d'una comunitat colonial a Empúries (segles VI-V a.C.)

El projecte de recerca ha analitzat la complexitat ètnica, social i de gènere de l’assentament colonial d’Empúries a partir d’un exhaustiu estudi contextual dels materials ceràmics d’ús quotidià provinents de diferents sectors excavats als anys vuitanta a la Neàpolis empuritana (N-7000, N-5000 i N-1000). Aquests sectors datats els segles V i IV a.C., corresponen a un moment en el qual s’està construint una identitat colonial diferenciada en l’espai emporità. La comparació dels aixovars domèstics usats quotidianament pels habitants d’aquest tres sectors situats en punts distants de la Neàpolis mostra que estem davant d’un assentament colonial molt heterogeni, on conflueixen materialitats, tradicions tecnològiques i pràctiques quotidianes iberes, gregues i híbrides. Els estudis de materials confirmen que no va existir una segregació espacial entre poblacions d’origen grec i iber en aquest assentament i suggereixen la cohabitació de gent –tant d’homes com de dones- d’origen grec i iber a les diferents àrees estudiades d’aquest assentament portuari. En el registre material ceràmic de totes les zones estudiades s’observa un predomini dels productes ceràmics de tradició grega (colonials o àtics) en els serveis de taula, però no així en els estris utilitzats a la cuina, que són majoritàriament de tradició ibera. També els materials ceràmics relacionats amb l’emmagatzematge domèstic i amb el transport i el comerç (àmfores) són predominantment ibers. Aquest patró suggereix que els ibers que van cohabitar amb els grecs al port emporità no van assumir en tots els casos posicions subalternes, sinó que tingueren un rol rellevant en la gestió econòmica i comercial d’aquest espai portuari.

Comparison of splice sites in mammals and chicken

We have carried out an initial analysis of the dynamics of the recent evolution of the splice-sites sequences on a large collection of human, rodent (mouse and rat), and chicken introns. Our results indicate that the sequences of splice sites are largely homogeneous within tetrapoda. We have also found that orthologous splice signals between human and rodents and within rodents are more conserved than unrelated splice sites, but the additional conservation can be explained mostly by background intron conservation. In contrast, additional conservation over background is detectable in orthologous mammalian and chicken splice sites. Our results also indicate that the U2 and U12 intron classes seem to have evolved independently since the split of mammals and birds; we have not been able to find a convincing case of interconversion between these two classes in our collections of orthologous introns. Similarly, we have not found a single case of switching between AT-AC and GT-AG subtypes within U12 introns, suggesting that this event has been a rare occurrence in recent evolutionary times. Switching between GT-AG and the noncanonical GC-AG U2 subtypes, on the contrary, does not appear to be unusual; in particular, T to C mutations appear to be relatively well tolerated in GT-AG introns with very strong donor sites.

SGP-1: prediction and validation of homologous genes based on sequence alignments

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of SGP-1 depends little on species-specific properties such as codon usage or the nucleotide distribution. SGP-1 may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.