Metodi computazionali per l'annotazione di genomi e proteomi

Il progresso tecnologico nel campo della biologia molecolare, pone la comunità scientifica di fronte all’esigenza di dare un’interpretazione all’enormità di sequenze biologiche che a mano a mano vanno a costituire le banche dati, siano esse proteine o acidi nucleici. In questo contesto la bioinformatica gioca un ruolo di primaria importanza. Un nuovo livello di possibilità conoscitive è stato introdotto con le tecnologie di Next Generation Sequencing (NGS), per mezzo delle quali è possibile ottenere interi genomi o trascrittomi in poco tempo e con bassi costi. Tra le applicazioni del NGS più rilevanti ci sono senza dubbio quelle oncologiche che prevedono la caratterizzazione genomica di tessuti tumorali e lo sviluppo di nuovi approcci diagnostici e terapeutici per il trattamento del cancro. Con l’analisi NGS è possibile individuare il set completo di variazioni che esistono nel genoma tumorale come varianti a singolo nucleotide, riarrangiamenti cromosomici, inserzioni e delezioni. Va però sottolineato che le variazioni trovate nei geni vanno in ultima battuta osservate dal punto di vista degli effetti a livello delle proteine in quanto esse sono le responsabili più dirette dei fenotipi alterati riscontrabili nella cellula tumorale. L’expertise bioinformatica va quindi collocata sia a livello dell’analisi del dato prodotto per mezzo di NGS ma anche nelle fasi successive ove è necessario effettuare l’annotazione dei geni contenuti nel genoma sequenziato e delle relative strutture proteiche che da esso sono espresse, o, come nel caso dello studio mutazionale, la valutazione dell’effetto della variazione genomica. È in questo contesto che si colloca il lavoro presentato: da un lato lo sviluppo di metodologie computazionali per l’annotazione di sequenze proteiche e dall’altro la messa a punto di una pipeline di analisi di dati prodotti con tecnologie NGS in applicazioni oncologiche avente come scopo finale quello della individuazione e caratterizzazione delle mutazioni genetiche tumorali a livello proteico.

In vivo konditionale Depletion von latentem murinen Cytomegalovirus

Die Lunge stellt einen Hauptort der CMV-Latenz dar. Die akute CMV-Infektion wird durch infiltrierende antivirale CD8 T-Zellen terminiert. Das virale Genom verbleibt jedoch im Lungengewebe in einem nicht replikativen Zustand, der Latenz, erhalten. Es konnte bereits gezeigt werden, dass während der Latenz die Major Immediate Early- (MIE) Gene ie1- und ie2 sporadisch transkribiert werden. Bisher konnte diese beginnende Reaktivierung latenter CMV-Genome nur in einer Momentaufnahme gezeigt werden (Kurz et al., 1999; Grzimek et al., 2001; Simon et al., 2005; zur Übersicht: Reddehase et al., 2008). Die sporadische Expression der MIE-Gene führt jedoch zur Präsentation eines antigenen IE1-Peptids und somit zur Stimulation antiviraler IE1-Peptid-spezifischer CD8 T-Zellen, die durch ihre Effektorfunktion die beginnende Reaktivierung wieder beenden. Dies führte uns zu der Hypothese, dass MIE-Genexpression über einen Zeitraum betrachtet (period prevalence) häufiger stattfindet als es in einer Momentaufnahme (point prevalence) beobachtet werden kann.rnrnUm die Häufigkeit der MIE-Genexpression in der Dynamik in einem definierten Zeitraum zu erfassen, sollte eine Methode entwickelt werden, welche es erstmals ermöglicht, selektiv und konditional transkriptionell aktive Zellen sowohl während der akuten Infektion als auch während der Latenz auszulöschen. Dazu wurde mit Hilfe der Zwei-Schritt BAC-Mutagenese ein rekombinantes death-tagged Virus hergestellt, welches das Gen für den Diphtherie Toxin Rezeptor (DTR) unter Kontrolle des ie2-Promotors (P2) enthält. Ist der P2 transkriptionell aktiv, wird der DTR an der Zelloberfläche präsentiert und die Zelle wird suszeptibel für den Liganden Diphtherie Toxin (DT). Durch Gabe von DT werden somit alle Zellen ausgelöscht, in denen virale Genome transkriptionell aktiv sind. Mit zunehmender Dauer der DT-Behandlung sollte also die Menge an latenten viralen Genomen abnehmen.rnrnIn Western Blot-Analysen konnte das DTR-Protein bereits 2h nach der Infektion nachgewiesen werden. Die Präsentation des DTR an der Zelloberfläche wurde indirekt durch dessen Funktionalität bewiesen. Das rekombinante Virus konnte in Fibroblasten in Gegenwart von DT nicht mehr replizieren. In akut infizierten Tieren konnte die virale DNA-Menge durch eine einmalige intravenöse (i.v.) DT-Gabe signifikant reduziert werden. Verstärkt wurde dieser Effekt durch eine repetitive i.v. DT-Gabe. Auch während der Latenz gelang es, die Zahl der latenten viralen Genome durch repetitive i.v. und anschließende intraperitoneale (i.p.) DT-Gabe zu reduzieren, wobei wir abhängig von der Dauer der DT-Gabe eine Reduktion um 60\% erreichen konnten. Korrespondierend zu der Reduktion der DNA-Menge sank auch die Reaktivierungshäufigkeit des rekombinanten Virus in Lungenexplantatkulturen. rnrnrnUm die Reaktivierungshäufigkeit während der Latenz berechnen zu können, wurde durch eine Grenzverdünnungsanalyse die Anzahl an latenten viralen Genomen pro Zelle bestimmt. Dabei ergab sich eine Kopienzahl von 9 (6 bis 13). Ausgehend von diesen Ergebnissen lässt sich berechnen, dass, bezogen auf die gesamte Lunge, in dem getesteten Zeitraum von 184h durch die DT-Behandlung 1.000 bis 2.500 Genome pro Stunde ausgelöscht wurden. Dies entspricht einer Auslöschung von 110 bis 280 MIE-Gen-exprimierenden Lungenzellen pro Stunde. Damit konnte in dieser Arbeit erstmals die Latenz-assoziierte Genexpression in ihrer Dynamik dargestellt werden.rn

Charakterisierung der Latenz des murinen Cytomegalovirus in vivo

Die primäre, produktive Cytomegalovirus (CMV)-Infektion wird im immunkompetenten Patienten effizient durch antivirale CD8+ T-Zellen kontrolliert. Das virale Genom besitzt jedoch die Fähigkeit, in einem nicht replikativen, Latenz genannten Zustand, in gewissen Zelltypen zu persistieren, ohne dass infektiöse Nachkommenviren produziert werden. Die molekularen Mechanismen, welche der Etablierung und Aufrechterhaltung der Latenz zugrundeliegen, sind noch weitestgehend unbekannt. Es gibt Hinweise darauf, dass zelluläre Verteidigungsmechanismen die Zirkularisierung und Chromatinisierung viraler Genome hervorrufen und dadurch die virale Genexpression größtenteils verhindert wird (Marks & Spector, 1984; Reeves et al., 2006).rnAllerdings liegen die Genome nicht in einem komplett inaktiven Zustand vor. Vielmehr konnte für das murine CMV (mCMV) bereits die sporadische Transkription der Gene ie1 und ie2 während der Latenz nachgewiesen werden (Kurz et al., 1999; Grzimek et al., 2001).rnIn der vorliegenden Arbeit wurde zum ersten Mal eine umfassende in vivo Latenz-Analyse zur Charakterisierung der viralen Transkription in einer Kinetik anhand der alle drei kinetischen Klassen repräsentierenden Transkripte IE1, IE3, E1, m164, M105 und M86 vorgenommen.rnNach Latenz-Etablierung, verifiziert durch Abwesenheit von infektiösem Virus, konnten alle getesteten Transkripte in der Lunge quantifiziert werden. Interessanterweise war die transkriptionelle Aktivität zu keinem Analyse-Zeitpunkt mit der klassischen IE-E-L-Kinetik der produktiven Infektion kompatibel. Stattdessen lag eine stochastische Transkript-Expression vor, deren Aktivität mit voranschreitender Zeit immer weiter abnahm.rnWährend der Latenz exprimierte Transkripte, die für antigene Peptide kodieren, können infizierte Zellen für das Immunsystem sichtbar machen, was zu einer fortwährenden Restimulation des memory T-Zell-pools führen würde. Durch zeitgleiche Analyse der Transkript-Expression, sowie der Frequenzen Epitop-spezifischer CD8+ T-Zellen während der Latenz (IE1, m164, M105), wurde eine möglicher Zusammenhang zwischen der transkriptionellen Aktivität und der Expansion des memory T-Zell-pools untersucht. Die weitere Charakterisierung von Subpopulationen der Epitop-spezifischen CD8+ T-Zellen identifizierte die SLECs (short-lived-effector cells; CD127low CD62Llow KLRG1high) als die dominante Population in Lunge und Milz während der mCMV-Latenz.rnIn einem weiteren Teil der Arbeit sollte untersucht werden, ob IE-Genexpression zur Etablierung von Latenz notwendig ist. Mit Hilfe der Rekombinanten mCMV-Δie2-DTR, die die Gensequenz des Diphtherietoxin-Rezeptors (DTR) anstelle des Gens ie2 trägt, konnten infizierte, DTR exprimierende Zellen durch eine DT-Applikation konditional depletiert werden.rnIm latent infizierbaren Zelltyp der Leber, den LSECs (liver sinusoidal endothelial cells) wurde die virale Load durch 90-stündige DT–Applikation nach mCMV-Δie2-DTR Infektion auf das Level latent infizierter LSECs reduziert. Diese Daten sprechen für die Hypothese eines von Beginn an inaktiven Genoms, das keine IE-Genexpression zur Latenz-Etablierung benötigt. Zusätzlich stellt dieser Ansatz ein neues Tier-Modell zur Latenz-Etablierung dar. Verringerte Wartezeiten bis zur vollständigen Latenz-Etablierung, im Vergleich zum bisherigen Knochenmarktransplantations-Modell, könnten anfallende Tierhaltungskosten erheblich reduzieren und das Voranschreiten der Forschung beschleunigen.

Phylogenomische Untersuchungen zur Aufklärung des Metazoa-Stammbaums

Zentrales Thema der Arbeit war die Aufklärung von Verwandtschaftsverhältnissen im „Tree of Life“ der vielzelligen Tiere (Metazoa) unter Einsatz großer DNA-Sequenzdatensätze und phylogenomischer Methoden. Zur Untersuchung der internen Phylogenie der Syndermata (= meist freilebende Rädertiere („Rotifera“) + endoparasitische Kratzwürmer (Acanthocephala)) sowie ihrer Position im Metazoen-Stammbaum wurden insgesamt sieben neue mitochondriale (mt) Genome sowie neue Transkriptom-Sequenzdaten von sieben verschiedenen Syndermata-Spezies generiert und/oder analysiert. Die Stammbaumrekonstruktionen auf Grundlage dieser sowie orthologer Sequenzen anderer Spezies in Form von phylogenomischen Datensätzen mit bis zu 82.000 Aminosäurepositionen ergaben folgende Aussagen zur Evolution: (i) Innerhalb der Acanthocephala bilden monophyletische Palaeacanthocephala das Schwestertaxon zu den Eoacanthocephala. Die Archiacanthocephala sind Schwestertaxon zu allen vorgenannten. (ii) Innerhalb der Syndermata bilden die epizoisch lebenden Seisonidea das Schwestertaxon zu den endoparasitischen Acanthocephala (= Pararotatoria), die Bdelloidea sind das Schwestertaxon zu den Pararotatoria (= Hemirotifera) und die Monogononta das Schwestertaxon zu den Hemirotifera. Die klassischen Eurotatoria (= Bdelloidea + Monogononta) sind demnach paraphyletisch. (iii) Innerhalb der Metazoa bilden die Syndermata gemeinsam mit den Gnathostomulida die Gnathifera. Diese sind die Schwestergruppe zu allen anderen Spiralia-Taxa, welche sich in Rouphozoa (= Platyhelminthes + Gastrotricha) sowie die Lophotrochozoa aufspalten. Die Platyzoa (= Gnathifera + Platyhelminthes + Gastrotricha) sind demnach paraphyletisch. Diese phylogenetischen Hypothesen wurden im Hinblick auf ihre Implikationen für die Evolution morphologischer und ökologischer Merkmale interpretiert. Demnach sind während der Evolution dieser Tiergruppen mehrfach sekundäre Verlustereignisse von komplexen morphologischen Merkmalen aufgetreten (laterale sensorische Organe innerhalb der Acanthocephala und das Räderorgan (Corona) innerhalb der Syndermata), was die Verwendung dieser Merkmale im Sinne einer klassisch-morphologischen Phylogenetik kritisch erscheinen lässt. Der Endoparasitismus der Acanthocephala hat sich wahrscheinlich über ein epizoisches Zwischenstadium, wie man es heute noch bei den Seisonidea findet, entwickelt. Der letzte gemeinsame Vorfahre der Spiralia war vermutlich klein und unsegmentiert und besaß keine echte Leibeshöhle (Coelom). Demnach hätten sich Segmentierung und Coelome innerhalb der Metazoa mehrfach unabhängig voneinander (konvergent) entwickelt. Die Arbeit beinhaltete folgende weitere, zum Teil methodische Aspekte: (i) die Analyse der Architektur der mt Genome der Monogononta bestätigte die aberrante Organisation in zwei Subgenomen für die Brachionidae. (ii) Eine Prüfung der Tauglichkeit ribosomaler Proteine für molekular-phylogenetische Arbeiten ergab das Vorhandensein widersprüchlicher phylogenetischer Signale in diesen speziellen Proteinsequenzen. (iii) Es konnte nachgewiesen werden, dass systematische Fehler wie „long-branch attraction“ bei der Positionierung der Syndermata im Stammbaum der Metazoa eine große Rolle spielen und adressiert werden müssen.

RNA interference in mammalian cell systems

In the last decade, few areas of biology have been transformed as thoroughly as RNA molecular biology. Without any doubt, one of the most significant advances has been the discovery of small (20-30 nucleotide) noncoding RNAs that regulate genes and genomes. The effects of small RNAs on gene expression and control are generally inhibitory, and the corresponding regulatory mechanisms are therefore collectively subsumed under the heading of RNA silencing and/or RNA interference. Two primary categories of these small RNAs - short interfering RNAs (siRNAs) and microRNAs (miRNAs) - act in both somatic and germline lineages of eukaryotic species to regulate endogenous genes and to defend the genome from invasive nucleic acids. Recent advances have revealed unexpected diversity in their biogenesis pathways and the regulatory mechanisms that they access. Our understanding of siRNA and miRNA-based regulation has direct implications for fundamental biology as well as disease aetiology and treatment as it is discussed in this review on 'new techniques in molecular biology'.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

Comparative genomics of metabolic networks of free-living and parasitic eukaryotes

Background Obligate endoparasites often lack particular metabolic pathways as compared to free-living organisms. This phenomenon comprises anabolic as well as catabolic reactions. Presumably, the corresponding enzymes were lost in adaptation to parasitism. Here we compare the predicted core metabolic graphs of obligate endoparasites and non-parasites (free living organisms and facultative parasites) in order to analyze how the parasites' metabolic networks shrunk in the course of evolution. Results Core metabolic graphs comprising biochemical reactions present in the presumed ancestor of parasites and non-parasites were reconstructed from the Kyoto Encyclopedia of Genes and Genomes. While the parasites' networks had fewer nodes (metabolites) and edges (reactions), other parameters such as average connectivity, network diameter and number of isolated edges were similar in parasites and non-parasites. The parasites' networks contained a higher percentage of ATP-consuming reactions and a lower percentage of NAD-requiring reactions. Control networks, shrunk to the size of the parasites' by random deletion of edges, were scale-free but exhibited smaller diameters and more isolated edges. Conclusions The parasites' networks were smaller than those of the non-parasites regarding number of nodes or edges, but not regarding network diameters. Network integrity but not scale-freeness has acted as a selective principle during the evolutionary reduction of parasite metabolism. ATP-requiring reactions in particular have been retained in the parasites' core metabolism while NADH- or NADPH-requiring reactions were lost preferentially.

Glycomic analysis of human mast cells, eosinophils and basophils

In allergic diseases such as asthma, eosinophils, basophils and mast cells, through release of preformed and newly generated mediators, granule proteins and cytokines, are recognized as key effector cells. While their surface protein phenotypes, mediator release profiles, ontogeny, cell trafficking and genomes have been generally explored and compared, there has yet to be any thorough analysis and comparison of their glycomes. Such studies are critical to understand the contribution of carbohydrates to the induction and regulation of allergic inflammatory responses and are now possible using improved technologies for detecting and characterizing cell-derived glycans. We thus report here the application of high-sensitivity mass spectrometric-based glycomics methodologies to the analysis of N-linked glycans derived from isolated populations of human mast cells, eosinophils and basophils. The samples were subjected to matrix-assisted laser desorption ionization (MALDI) time-of-flight (TOF) screening analyses and MALDI-TOF/TOF sequencing studies. Results reveal substantive quantities of terminal N-acetylglucosamine containing structures in both the eosinophil and the basophil samples, whereas mast cells display greater relative quantities of sialylated terminal epitopes. For the first time, we characterize the cell surface glycan structures of principal allergic effector cells, which by interaction with glycan-binding proteins (e.g. lectins) have the possibility to dictate cellular functions, and might thus have important implications for the pathogenesis of inflammatory and allergic diseases.

qPCR-based mitochondrial DNA quantification: Influence of template DNA fragmentation on accuracy

Real-time PCR (qPCR) is the method of choice for quantification of mitochondrial DNA (mtDNA) by relative comparison of a nuclear to a mitochondrial locus. Quantitative abnormal mtDNA content is indicative of mitochondrial disorders and mostly confines in a tissue-specific manner. Thus handling of degradation-prone bioptic material is inevitable. We established a serial qPCR assay based on increasing amplicon size to measure degradation status of any DNA sample. Using this approach we can exclude erroneous mtDNA quantification due to degraded samples (e.g. long post-exicision time, autolytic processus, freeze-thaw cycles) and ensure abnormal DNA content measurements (e.g. depletion) in non-degraded patient material. By preparation of degraded DNA under controlled conditions using sonification and DNaseI digestion we show that erroneous quantification is due to the different preservation qualities of the nuclear and the mitochondrial genome. This disparate degradation of the two genomes results in over- or underestimation of mtDNA copy number in degraded samples. Moreover, as analysis of defined archival tissue would allow to precise the molecular pathomechanism of mitochondrial disorders presenting with abnormal mtDNA content, we compared fresh frozen (FF) with formalin-fixed paraffin-embedded (FFPE) skeletal muscle tissue of the same sample. By extrapolation of measured decay constants for nuclear DNA (λnDNA) and mtDNA (λmtDNA) we present an approach to possibly correct measurements in degraded samples in the future. To our knowledge this is the first time different degradation impact of the two genomes is demonstrated and which evaluates systematically the impact of DNA degradation on quantification of mtDNA copy number.

Mitochondrial tRNA import and its consequences for mitochondrial translation

The mitochondrial genomes of most eukaryotes lack a variable number of tRNA genes. This lack is compensated for by import of a small fraction of the corresponding cytosolic tRNAs. There are two broad mechanisms for the import of tRNAs into mitochondria. In the first one, the tRNA is coimported together with a mitochondrial precursor protein along the protein import pathway. It applies to the yeast tRNA(Lys) and has been elucidated in great detail. In the second more vaguely defined mechanism, which is mainly found in plants and protozoa, tRNAs are directly imported independent of cytosolic factors. However, results in plants indicate that direct import of tRNAs may nevertheless require some components of the protein import machinery. All imported tRNAs in all systems are of the eukaryotic type but need to be functionally integrated into the mitochondrial translation system of bacterial descent. For some tRNAs, this is not trivial and requires unique evolutionary adaptations.

Prevalence of genes encoding extracellular proteases in Staphylococcus aureus - important targets triggering immune response in vivo

Proteases of Staphylococcus aureus have long been considered to function as important virulence factors, although direct evidence of the role of particular enzymes remains incomplete and elusive. Here, we sought to provide a collective view of the prevalence of extracellular protease genes in genomes of commensal and pathogenic strains of S. aureus and their expression in the course of human and mouse infection. Data on V8 protease, staphopains A and B, aureolysin, and the recently described and poorly characterized group of six Spl proteases are provided. A phylogenetically diverse collection of 167 clinical isolates was analyzed, resulting in the comprehensive genetic survey of the prevalence of protease-encoding genes. No correlation between identified gene patterns with specific infections was established. Humoral response against the proteases of interest was examined in the sera derived from human patients and from a model mouse infection. The analysis suggests that at least some, if not all, tested proteases are expressed and secreted during the course of infection. Overall, the results presented in this study support the hypothesis that the secretory proteases as a group may contribute to the virulence of S. aureus.

Utilization of Probabilistic Models in Short Read Assembly from Second-Generation Sequencing

With the advent of cheaper and faster DNA sequencing technologies, assembly methods have greatly changed. Instead of outputting reads that are thousands of base pairs long, new sequencers parallelize the task by producing read lengths between 35 and 400 base pairs. Reconstructing an organism’s genome from these millions of reads is a computationally expensive task. Our algorithm solves this problem by organizing and indexing the reads using n-grams, which are short, fixed-length DNA sequences of length n. These n-grams are used to efficiently locate putative read joins, thereby eliminating the need to perform an exhaustive search over all possible read pairs. Our goal was develop a novel n-gram method for the assembly of genomes from next-generation sequencers. Specifically, a probabilistic, iterative approach was utilized to determine the most likely reads to join through development of a new metric that models the probability of any two arbitrary reads being joined together. Tests were run using simulated short read data based on randomly created genomes ranging in lengths from 10,000 to 100,000 nucleotides with 16 to 20x coverage. We were able to successfully re-assemble entire genomes up to 100,000 nucleotides in length.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Genome-Wide Survey of Pseudogenes in 80 Fully Re-sequenced Arabidopsis thaliana Accessions

Pseudogenes (Ψs), including processed and non-processed Ψs, are ubiquitous genetic elements derived from originally functional genes in all studied genomes within the three kingdoms of life. However, systematic surveys of non-processed Ψs utilizing genomic information from multiple samples within a species are still rare. Here a systematic comparative analysis was conducted of Ψs within 80 fully re-sequenced Arabidopsis thaliana accessions, and 7546 genes, representing ~28% of the genomic annotated open reading frames (ORFs), were found with disruptive mutations in at least one accession. The distribution of these Ψs on chromosomes showed a significantly negative correlation between Ψs/ORFs and their local gene densities, suggesting a higher proportion of Ψs in gene desert regions, e.g. near centromeres. On the other hand, compared with the non-Ψ loci, even the intact coding sequences (CDSs) in the Ψ loci were found to have shorter CDS length, fewer exon number and lower GC content. In addition, a significant functional bias against the null hypothesis was detected in the Ψs mainly involved in responses to environmental stimuli and biotic stress as reported, suggesting that they are likely important for adaptive evolution to rapidly changing environments by pseudogenization to accumulate successive mutations.

The Legal Protection of Fundamental Human Values in Central and Eastern Europe and Modern Biotechnology: Towards European Harmonization

The protection of the fundamental human values (life, bodily integrity, human dignity, privacy) becomes imperative with the rapid progress in modern biotechnology, which can result in major alterations in the genetic make-up of organisms. It has become possible to insert human genes into pigs so that their internal organs coated in human proteins are more suitable for transplantation into humans (xenotransplantation), and micro-organisms that cam make insulin have been created, thus changing the genetic make-up of humans. At the end of the 1980s, the Central and Eastern European (CEE) countries either initiated new legislation or started to amend existing laws in this area (clinical testing of drugs, experiments on man, prenatal genetic diagnosis, legal protection of the embryo/foetus, etc.). The analysis here indicates that the CEE countries have not sufficiently adjusted their regulations to the findings of modern biotechnology, either because of the relatively short period they have had to do so, or because there are no definite answers to the questions which modern biotechnology has raised (ethical aspects of xenotransplantation, or of the use of live-aborted embryonic or foetal tissue in neuro-transplantation, etc.). In order to harmonise the existing regulations in CEE countries with respect to the EU and supranational contexts, two critical issues should be taken into consideration. The first is the necessity for CEE countries to recognise the place of humans within the achievements of modern biotechnology (a broader affirmation of the principle of autonomy, an explicit ban on the violation of the genetic identity of either born or unborn life, etc.). The second concerns the definition of the status of different biotechnological procedures and their permissibility (gene therapy, therapeutic genomes, xenotransplantation, etc.). The road towards such answers may be more easily identified once all CEE countries become members of the Council of Europe and express their wish to join the EU, which in turn presupposes taking over the entire body of EU legislation.