Genetics for clinicians: from candidate genes to whole genome scans (technological advances).

Human genetics has progressed at an unprecedented pace during the past 10 years. DNA microarrays currently allow screening of the entire human genome with high level of coverage and we are now entering the era of high-throughput sequencing. These remarkable technical advances are influencing the way medical research is conducted and have boosted our understanding of the structure of the human genome as well as of disease biology. In this context, it is crucial for clinicians to understand the main concepts and limitations of modern genetics. This review will describe key concepts in genetics, including the different types of genetic markers in the human genome, review current methods to detect DNA variation, describe major online public databases in genetics, explain key concepts in statistical genetics and finally present commonly used study designs in clinical and epidemiological research. This review will therefore concentrate on human genetic variation analysis.

Long-distance wound signalling in arabidopsis

The leaves of all plants use elaborate and inducible defence systems to protect themselves. A wide variety of such defences are known and they include defence chemicals such as alkaloids, phenolics and terpenes, physical structures ranging from fibre cells to silica deposits, and a wide variety of defence proteins many of which target digestive processes in herbivores. It has long been known that the defence responses of plants under attack by insects are not restricted to the site of attack. Instead, if a leaf is damaged, defence can be triggered in other parts of the plant body, for example in distal leaves or even in roots and flowers. This raises the question of what are the organ-to-organ signals that coordinate this process. Several hypotheses have been proposed. These include the long-distance transfer of chemical signals through the plant vasculature, hydraulic signals that may transit through the xylem, and electrical signals that would move through living tissues such as the phloem. Much evidence for each of these scenarios has been published. In this thesis we took advantage of the fact that many plant defence responses are regulated by a signal transduction pathway based on a molecule called jasmonic acid. We used this molecule, one of its derivatives (jasmonoyl-isoleucine), and some of the genes it regulates as markers. Using these we investigated the possible role of the electrical signals in the leaf- to-leaf activation of the jasmonate pathway. We found that feeding insects stimulate easily detected electrical activity in the leaves of Arabidopsis thaliana and we used non-invasive surface electrodes to record this activity. This approach showed that jasmonate pathway activity and the electrical activity provoked by mechanical wounding occurred within identical spatial boundaries. Measurements of the apparent speed of surface potentials agreed well with previous velocity estimates for the speed of leaf-to-leaf signals that activate the jasmonate pathway. Using this knowledge we were able to investigate the effects of current injection into Arabidopsis leaves. This resulted in the strong expression of many jasmonate-regulated genes. All these results showed that electrical activity and the activation of jasmonate signalling were highly correlated. In order to test for possible causal links between the two processes, we conducted a small-scale reverse genetic screen on a series of T-DNA insertion mutants in ion channel genes and in other genes encoding proteins such as proton pumps. This screen, which was based on surface potential measurements, revealed that mutations in genes related to ionotropic glutamate receptors in animals had impaired electrical activity after wounding. Combining mutation of two of these glutamate-receptor-like genes in a double mutant reduced the response of leaves to current injection. When a leaf of this double mutant was wounded it failed to transmit a long-distance signal to a distal leaf. This result distinguished the double mutant from the wild-type plant and provides the first genetic evidence that electrical signalling is necessary to coordinate defence responses between organs in plants. - Les feuilles des plantes disposent de systèmes de défense inductibles très élaborés. Un grand nombre de ces systèmes de défenses sont connus et sont basés sur des composés chimiques comme les alcaloïdes, les composés phénoliques ou les terpènes, des systèmes physiques allant de la production de cellules fibreuses aux cristaux de silice ainsi qu'un grand nombre de protéines de défense ciblant le processus digestif des herbivores. Il est connu dépuis longtemps que la réponse défensive de la plante face à l'attaque pas un insecte n'est pas seulement localisée au niveau de la zone d'attaque. A la place, si une feuille est attaquée, les systèmes de défense peuvent être activés ailleurs dans la plante, comme par exemple dans d'autres feuilles, les racines ou même les fleurs. Ces observations soulèvent la question de la nature des signaux d'organes à organes qui régulent ces systèmes. Plusieurs hypothèses ont été formulées; une ou plusieures molécules pourraient être véhiculées dans la plante grâce au système vasculaire, un signal hydraulique transmis au travers du xylème ou encore des signaux électriques transmis par les cellules comme dans le phloème par exemple. De nombreuses études ont été publiées sur ces différentes hypothèses. Dans ce travail de thèse, nous avons choisi d'utiliser à notre avantage le fait que de nombreuses réponses de défense de la plante sont régulées par une même voie de signalisation utilisant l'acide jasmonique. Nous avons utilisé comme marqueurs cette molécule, un de ses dérivés (le jasmonoyl-isoleucine) ainsi que certains des gènes que l'acide jasmonique régule. Nous avons alors testé l'implication de la transmission de signaux électriques dans l'activation de la voie du jasmonate de feuille à feuille. Nous avons découvert que les insectes qui se nourrissent de feuilles d'Arabidopsis thaliana activent un signal électrique que nous avons pu mesurer grâce à une technique non invasive d'électrodes de surface. Les enregistrements ont montré que la génération de signaux électriques et l'activation de la voie du jasmonate avaient lieu aux mêmes endroits. La mesure de la vitesse de déplacement des impulsions électriques correspond aux estimations faites concernant l'activation de la voie du jasmonate. Grâce à cela, nous avons pu tester l'effet d'injection de courant électrique dans les feuilles d'Arabidopsis. La conséquence a été une forte expression de nombreux gènes de la voie du jasmonate, suggérant une forte corrélation entre l'activité électrique et l'activation de la voie du jasmonate. Afin de tester le lien de cause entre ces deux phénomènes, nous avons entrepris un criblage génétique sur une série de mutants d'insertion à l'ADN-T dans des gènes de canaux ioniques et d'autres gènes d'intérêt comme les gènes des pompes à protons. Ce criblage, basé sur la mesure de potentiels de surface, a permis de montrer que plusieurs mutations de gènes liés aux récepteurs au glutamate ionotropique présentent une baisse drastique de leurs activités électriques après une blessure mécanique des feuilles par rapport au type sauvage. Par la combinaison de deux mutations de ces récepteurs au glutamate en un double mutant, on obtient une réponse à la stimulation électrique encore plus faible. Quand une feuille du double mutant est blessée, elle est incapable de transmettre un signal à longue distance vers une feuille éloignée. Ce résultat permet de distinguer le double mutant de la plante sauvage et amène la première preuve génétique que l'activité électrique est nécessaire pour coordonner les réponses de défense entre les organes chez les plantes.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.

Evaluation de la sarcopénie de la personne agée par la bio-impédance spectroscopique [Multispectral bioimpedance and sarcopenia in the elderly].

The definition of sarcopenia includes both a loss of muscle strength and a decline in functional quality in addition to the loss of muscle protein mass. Multispectral bioimpendance allows bedside assessment of muscle mass. Using this new tool, we performed a pilot study to look for a possible correlation between muscle mass and various tests of muscle strength (grip strength, key-pitch, tip-pinch) and with functional tests (walk speed on 10 meters and Tinetti test). Our study demonstrates a good correlation between muscle mass determined by spectroscopic bioimpendance and muscle strength assessment, but no correlation with functional tests.

Sox genes in C. elegans: sox-2 role in postembryonic development

Report for the scientific sojourn carried out at the Columbia University, United States, from 2010 to 2012. Expression of SoxB genes correlates with the commitment of cells to a neural fate; however, the relevance of SoxB proteins in early vertebrate neurogenesis has been difficult to prove genetically due to embryonic lethality and presumed redundant functions. The nematode C. Elegants has only 5 sox genes: sox-2 and sox-3 form the SoxB group while sem-2, sox-4 and egl-13 belong to other Sox groups. Our results show that sox-2 and sem-2 are the sox genes expressed earliest and in a broader manner during embryogenesis, being expressed in several neuronal progenitors. sox-3, sox-4 and egl-13 are expressed in few cells during late embryogenesis, when most neurons are already born. Both sox-2 and sem-2 null mutants are early larval lethal but do not show neuronal specification defects during embryonic development as indicated by quantification of a panneuronal reporter. Potential redundancy or compensatory mechanisms between different sox genes have been ruled out, strongly suggesting that sox genes are not required for specification of embryonically-derived neurons. However, at the first larval stage there are still several blast cells that will give rise to different postembryonic lineages, which generate several neurons amongst other cell types. nterestingly, sox-2 is expressed in many of these progenitor cells. Using mosaic analysis we have so far identified neurons derived from two different postembryonic lineages which fail to be generated in C. elegans sox-2 mutants. These results support the idea that postembryonic progenitor competence is compromised in the absence of sox-2.

Nanomotors: new strategies of (bio)sensing based on motion

Nanomotors are nanoscale devices capable of converting energy into movement and forces. Among them, self-propelled nanomotors offer considerable promise for developing new and novel bioanalytical and biosensing strategies based on the direct isolation of target biomolecules or changes in their movement in the presence of target analytes. The mainachievements of this project consists on the development of receptor-functionalized nanomotors that offer direct and rapid target detection, isolation and transport from raw biological samples without preparatory and washing steps. For example, microtube engines functionalized with aptamer, antibody, lectin and enzymes receptors were used for the direct isolation of analytes of biomedical interest, including proteins and whole cells, among others. A target protein was also isolated from a complex sample by using an antigen-functionalized microengine navigating into the reservoirs of a lab-on-a-chip device. The new nanomotorbased target biomarkers detection strategy not only offers highly sensitive, rapid, simple and low cost alternative for the isolation and transport of target molecules, but also represents a new dimension of analytical information based on motion. The recognition events can be easily visualized by optical microscope (without any sophisticated analytical instrument) to reveal the target presence and concentration. The use of artificial nanomachines has shown not only to be useful for (bio)recognition and (bio)transport but also for detection of environmental contamination and remediation. In this context, micromotors modified with superhydrophobic layer demonstrated that effectively interacted, captured, transported and removed oil droplets from oil contaminated samples. Finally, a unique micromotor-based strategy for water-quality testing, that mimics live-fish water-quality testing, based on changes in the propulsion behavior of artificial biocatalytic microswimmers in the presence of aquatic pollutants was also developed. The attractive features of the new micromachine-based target isolation and signal transduction protocols developed in this project offer numerous potential applications in biomedical diagnostics, environmental monitoring, and forensic analysis.

Algorithmic aspects of bio-inspired string operations

We present building blocks for algorithms for the efficient reduction of square factor, i.e. direct repetitions in strings. So the basic problem is this: given a string, compute all strings that can be obtained by reducing factors of the form zz to z. Two types of algorithms are treated: an offline algorithm is one that can compute a data structure on the given string in advance before the actual search for the square begins; in contrast, online algorithms receive all input only at the time when a request is made. For offline algorithms we treat the following problem: Let u and w be two strings such that w is obtained from u by reducing a square factor zz to only z. If we further are given the suffix table of u, how can we derive the suffix table for w without computing it from scratch? As the suffix table plays a key role in online algorithms for the detection of squares in a string, this derivation can make the iterated reduction of squares more efficient. On the other hand, we also show how a suffix array, used for the offline detection of squares, can be adapted to the new string resulting from the deletion of a square. Because the deletion is a very local change, this adaption is more eficient than the computation of the new suffix array from scratch.

La vénerie royale / Robert de Salnove ; précédée d'une notice bio-bibliographique sur l'auteur, par le commandant G. de Marolles...

Collection : Les maîtres de la vénerie ; 4

Repression of flagellar genes in exponential phase by CsgD and CpxR, two crucial modulators of Escherichia coli biofilm formation.

Escherichia coli adapts its lifestyle to the variations of environmental growth conditions, swapping between swimming motility or biofilm formation. The stationary-phase sigma factor RpoS is an important regulator of this switch, since it stimulates adhesion and represses flagellar biosynthesis. By measuring the dynamics of gene expression, we show that RpoS inhibits the transcription of the flagellar sigma factor, FliA, in exponential growth phase. RpoS also partially controls the expression of CsgD and CpxR, two transcription factors important for bacterial adhesion. We demonstrate that these two regulators repress the transcription of fliA, flgM, and tar and that this regulation is dependent on the growth medium. CsgD binds to the flgM and fliA promoters around their -10 promoter element, strongly suggesting direct repression. We show that CsgD and CpxR also affect the expression of other known modulators of cell motility. We propose an updated structure of the regulatory network controlling the choice between adhesion and motility.

Biogeochemical evolution of a marine shore tailings deposit during bioremediation

Summary : Mining activities produce enormous amounts of waste material known as tailings which are composed of fine to medium size particles. These tailings often contain sulfides, which oxidation can lead to acid and metal contamination of water; therefore they need to be remediated. In this work a tailings bioremediation approach was investigated by an interdisciplinary study including geochemistry, mineralogy and microbiology. The aim of the work was to study the effect of the implementation of wetland above oxidizing tailings on the hydrogeology and the biogeochemical element cycles, and to assess the system evolution over time. To reach these goals, biogeochemical processes occurring in a marine shore tailings deposit were investigated. The studied tailings deposit is located at the Bahìa de Ite, Pacific Ocean, southern Peru, where between 1940 and 1996 the tailings were discharged from the two porphyry copper mines Cuajone and Toquepala. After the end of deposition, a remediation approach was initiated in 1997 with a wetland implementation above the oxidizing tailings. Around 90% of the tailings deposits (total 16 km2) were thus remediated, except the central delta area and some areas close to the shoreline. The multi-stable isotope study showed that the tailings were saturated with fresh water in spite of the marine setting, due to the high hydraulic gradient resulting from the wetland implementation. Submarine groundwater discharge (SGD) was the major source of SO4 2-, C1-, Na+, Fe2+, and Mn2+ input into the tailings at the original shelf-seawater interface. The geochemical study (aquatic geochemistry and X-Ray diffraction (XRD) and sequential extractions from the solid fraction) showed that iron and sulfur oxidation were the main processes in the non-remediated tailings, which showed a top a low-pH oxidation zone with strong accumulation of efflorescent salts at the surface due to capillary upward transport of heavy metals (Fe, Cu, Zn, Mn, Cd, Co, and Ni) in the arid climate. The study showed also that the implementation of the wetland resulted in very low concentrations of heavy metals in solution (mainly under the detection limit) due to the near neutral pH and more reducing conditions (100-150 mV). The heavy metals, which were taken from solution, precipitated as hydroxides and sulfides or were bound to organic matter. The bacterial community composition analysis by Terminal Restriction Fragment Length Polymorphism (T-RFLP) and cloning and sequencing of 16S rRNA genes combined with a detailed statistical analysis revealed a high correlation between the bacterial distribution and the geochemical variables. Acidophilic autotrophic oxidizing bacteria were dominating the oxidizing tailings, whereas neutrophilic and heterotrophic reducing bacteria were driving the biogeochemical processes in the remediated tailings below the wetland. At the subsurface of the remediated tailings, an iron cycling was highlighted with oxidation and reduction processes due to micro-aerophilic niches provided by the plant rhizosphere in this overall reducing environment. The in situ bioremediation experiment showed that the main parameter to take into account for the effectiveness was the water table and chemistry which controls the system. The constructed remediation cells were more efficient and rapid in metal removal when saturation conditions were available. This study showed that the bioremediation by wetland implementation could be an effective and rapid treatment for some sulfidic mine tailings deposits. However, the water saturation of the tailings has to be managed on a long-term basis in order to guarantee stability. Résumé : L'activité minière produit d'énormes quantités de déchets géologiques connus sous le nom de « tailings » composées de particules de taille fine à moyenne. Ces déchets contiennent souvent des sulfures dont l'oxydation conduit à la formation d'effluents acides contaminés en métaux, d'où la nécessité d'effectuer une remédiation des sites de stockage concernés. Le but de ce travail est dans un premier temps d'étudier l'effet de la bio-remédiation d'un dépôt de tailings oxydés sur l'hydrogéologie du système et les cycles biogéochimiques des éléments et en second lieu, d'évaluer l'évolution du processus de remédiation dans le temps. Le site étudié dans ce travail est situé dans la Bahía de Ite, au sud du Pérou, au bord de l'Océan Pacifique. Les déchets miniers en question sont déposés dans un environnement marin. De 1940 à 1996, les déchets de deux mines de porphyre cuprifère - Cuajone et Toquepala - ont été acheminés sur le site via la rivière Locumba. En 1997, une première remédiation a été initiée avec la construction d'une zone humide sur les tailings. Depuis, environ 90% de la surface du dépôt (16 km2) a été traité, les parties restantes étant la zone centrale du delta du Locumba et certaines zones proches de la plage. Malgré la proximité de l'océan, les études isotopiques menées dans le cadre de ce travail ont montré que les tailings étaient saturés en eau douce. Cette saturation est due à la pression hydraulique résultant de la mise en place des zones humides. Un écoulement d'eau souterrain sous-marin a été à détecté à l'interface entre les résidus et l'ancien fond marin. En raison de la géologie locale, il constitue une source d'entrée de SO4 2-, Cl-, Na+, FeZ+, et Mn2+ dans le système. L'analyse de la géochimie aquatique, la Diffraction aux Rayons X (XRD) et l'extraction séquentielle ont montré que l'oxydation du fer et .des sulfures est le principal processus se produisant dans les déchets non remédiés. Ceci a entraîné le développement d'une zone d'oxydation à pH bas induisant une forte accumulation des sels efflorescents, conséquence de la migration capillaire des métaux lourds (Fe, Cu, Zn, Mn, Cd, Co et Ni) de la solution vers la surface dans ce climat aride. Cette étude a montré également que la construction de la zone humide a eu comme résultats une précipitation des métaux dans des phases minérales en raison du pH neutre et des conditions réductrices (100-150mV). Les métaux lourds ont précipité sous la forme d'hydroxydes et de sulfures ou sont adsorbés à la matière organique. L'analyse de la composition de la communauté bactérienne à l'aide la technique T-RFLP (Terminal Restriction Fragment Length Polymorphism) et par le clonage/séquençage des gènes de l'ARNr 16S a été combinée à une statistique détaillée. Cette dernière a révélé une forte corrélation entre la distribution de bactéries spécifiques et la géochimie : Les bactéries autotrophes acidophiles dominent dans les déchets oxydés non remédiés, tandis que des bactéries hétérotrophes neutrophiles ont mené les processus microbiens dans les déchets remédiés sous la zone humide. Sous la surface de la zone humide, nos analyses ont également mis en évidence un cycle du fer par des processus d'oxydoréduction rendus possibles par la présence de niches micro-aérées par la rhizosphère dans cet environnement réducteur. L'expérience de bio-remédiation in situ a montré que les paramètres clés qui contrôlent l'efficacité du traitement sont le niveau de la nappe aquifère et la chimie de l'eau. Les cellules de remédiation se sont montrées plus efficaces et plus rapides lorsque le système a pu être saturé en eau. Finalement, cette étude a montré que la bio-remédiation de déchets miniers par la construction de zones humides est un moyen de traitement efficace, rapide et peu coûteux. Cependant, la saturation en eau du système doit être gérée sur le long terme afin de garantir la stabilité de l'ensemble du système.

The CNGCb and CNGCd genes from Physcomitrella patens moss encode for thermosensory calcium channels responding to fluidity changes in the plasma membrane.

Land plants need precise thermosensors to timely establish molecular defenses in anticipation of upcoming noxious heat waves. The plasma membrane-embedded cyclic nucleotide-gated Ca(2+) channels (CNGCs) can translate mild variations of membrane fluidity into an effective heat shock response, leading to the accumulation of heat shock proteins (HSP) that prevent heat damages in labile proteins and membranes. Here, we deleted by targeted mutagenesis the CNGCd gene in two Physcomitrella patens transgenic moss lines containing either the heat-inducible HSP-GUS reporter cassette or the constitutive UBI-Aequorin cassette. The stable CNGCd knockout mutation caused a hyper-thermosensitive moss phenotype, in which the heat-induced entry of apoplastic Ca(2+) and the cytosolic accumulation of GUS were triggered at lower temperatures than in wild type. The combined effects of an artificial membrane fluidizer and elevated temperatures suggested that the gene products of CNGCd and CNGCb are paralogous subunits of Ca(2+)channels acting as a sensitive proteolipid thermocouple. Depending on the rate of temperature increase, the duration and intensity of the heat priming preconditions, terrestrial plants may thus acquire an array of HSP-based thermotolerance mechanisms against upcoming, otherwise lethal, extreme heat waves.

Bioinformática: interfaz web para mostrar las relaciones de expresión entre grupos de genes

El projecte que es presenta en aquesta memòria està composat per un preproces i una aplicació web. L’objectiu del treball realizat és mostrar les relacions d’expresió entre grups de gens. A més de dissenyar l’aplicació web en : http://revolutionresearch.uab.es que permetrà l’estudi dels resultats obtinguts.

Network-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanoma.

Cancer genomes frequently contain somatic copy number alterations (SCNA) that can significantly perturb the expression level of affected genes and thus disrupt pathways controlling normal growth. In melanoma, many studies have focussed on the copy number and gene expression levels of the BRAF, PTEN and MITF genes, but little has been done to identify new genes using these parameters at the genome-wide scale. Using karyotyping, SNP and CGH arrays, and RNA-seq, we have identified SCNA affecting gene expression ('SCNA-genes') in seven human metastatic melanoma cell lines. We showed that the combination of these techniques is useful to identify candidate genes potentially involved in tumorigenesis. Since few of these alterations were recurrent across our samples, we used a protein network-guided approach to determine whether any pathways were enriched in SCNA-genes in one or more samples. From this unbiased genome-wide analysis, we identified 28 significantly enriched pathway modules. Comparison with two large, independent melanoma SCNA datasets showed less than 10% overlap at the individual gene level, but network-guided analysis revealed 66% shared pathways, including all but three of the pathways identified in our data. Frequently altered pathways included WNT, cadherin signalling, angiogenesis and melanogenesis. Additionally, our results emphasize the potential of the EPHA3 and FRS2 gene products, involved in angiogenesis and migration, as possible therapeutic targets in melanoma. Our study demonstrates the utility of network-guided approaches, for both large and small datasets, to identify pathways recurrently perturbed in cancer.

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.

Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.