Rate of gene sequence evolution and species diversification in flowering plants: a re-evaluation

Barraclough and co-workers (in a paper published in 1996) observed that there was a significant positive correlation between the rate of evolution of the rbcL chloroplast gene within families of flowering plants and the number of species in those families. We tested three additional data sets of our own (based on both plastid and nuclear genes) and used methods designed specifically for the comparison of sister families (based on random speciation and extinction). We show that, over all sister groups, the correlation between the rate of gene evolution and an increased diversity is not always present. Despite tending towards a positive association, the observation of individual probabilities presents a U-shaped distribution of association (i.e. it can be either significantly positive or negative). We discuss the influence of both phylogenetic sampling and applied taxonomies on the results.

Advances in Sequence Analysis: Theory, Methods, Applications

This book gives a general view of sequence analysis, the statistical study of successions of states or events. It includes innovative contributions on life course studies, transitions into and out of employment, contemporaneous and historical careers, and political trajectories. The approach presented in this book is now central to the life-course perspective and the study of social processes more generally. This volume promotes the dialogue between approaches to sequence analysis that developed separately, within traditions contrasted in space and disciplines. It includes the latest developments in sequential concepts, coding, atypical datasets and time patterns, optimal matching and alternative algorithms, survey optimization, and visualization. Field studies include original sequential material related to parenting in 19th-century Belgium, higher education and work in Finland and Italy, family formation before and after German reunification, French Jews persecuted in occupied France, long-term trends in electoral participation, and regime democratization. Overall the book reassesses the classical uses of sequences and it promotes new ways of collecting, formatting, representing and processing them. The introduction provides basic sequential concepts and tools, as well as a history of the method. Chapters are presented in a way that is both accessible to the beginner and informative to the expert.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes.

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Neuroprotection in cerebral ischemia with XG-102, a new peptide inhibitor of JNK

Abstract Stroke or cerebrovascular accident, whose great majority is of ischemic nature, is the third leading cause of mortality and long lasting disability in industrialised countries. Resulting from the loss of blood supply to the brain depriving cerebral tissues of oxygen and glucose, it induces irreversible neuronal damages. Despite the large amount of research carried out into the causes and pathogenic features of cerebral ischemia the progress toward effective treatments has been poor. Apart the clot-busting drug tissue-type plasminogen activator (tPA) as effective therapy for acute stroke (reperfusion by thrombolysis) but limited to a low percentage of patients, there are currently no other approved medical treatments. The need for new therapy strategies is therefore imperative. Neuronal death in cerebral ischemia is among others due to excitotoxic mechanisms very early after stroke onset. One of the main involved molecular pathways leading to excitotoxic cell death is the c-Jun NH2-terminal kinase (JNK) pathway. Several studies have already shown the efficacy of a neuroprotective agent of a new type, a dextrogyre peptide synthesized in the retro inverso form (XG102, formerly D-JNKI1), which is protease-resistant and cell-penetrating and that selectively and strongly blocks the access of JNK to many of its targets. A powerful protection was observed with this compound in several models of ischemia (Borsello et al. 2003;Hirt et al. 2004). This chimeric compound, made up of a 10 amino acid TAT transporter sequence followed by a 20 amino acids JNK binding domain (JBD) sequence from JNK inhibitor protein (JIP) molecule, induced both a major reduction in lesion size and improved functional outcome. Moreover it presents a wide therapeutic window. XG-102 has proved its powerful efficacy in an occlusion model of middle cerebral artery in mice with intracérebroventricular (i.c.v.) injection but in order to be able to consider the development of this drug for human ischemic stroke it was therefore necessary to determine the feasibility of its systemic administration. The studies being the subject of this thesis made it possible to show a successful neuroprotection with XG-102 administered systemically after transient mouse middle cerebral artery occlusion (MCAo). Moreover our data. provided information about the feasibility to combine XG-102 with tPA without detrimental action on cell survival. By combining the benefits from a reperfusion treatment with the effects of a neuroprotective compound, it would represent the advantage of bringing better chances to protect the cerebral tissue. Résumé L'attaque cérébrale ou accident vasculaire cérébral, dont la grande majorité est de nature ischémique, constitue la troisième cause de mortalité et d'infirmité dans les pays industrialisés. Résultant de la perte d'approvisionnement de sang au cerveau privant les tissus cérébraux d'oxygène et de glucose, elle induit des dommages neuronaux irréversibles. En dépit du nombre élevé de recherches effectuées pour caractériser les mécanismes pathogènes de l'ischémie. cérébrale, les progrès vers des traitements efficaces restent pauvres. Excepté l'activateur tissulaire du plasminogène (tPA) dont le rôle est de désagréger les caillots sanguins et employé comme thérapie efficace contre l'attaque cérébrale aiguë (reperfusion par thrombolyse) mais limité à un faible pourcentage de patients, il n'y a actuellement aucun autre traitement médical approuvé. Le besoin de nouvelles stratégies thérapeutiques est par conséquent impératif. La mort neuronale dans l'ischémie cérébrale est entre autres due à des mécanismes excitotoxiques survenant rapidement après le début de l'attaque cérébrale. Une des principales voies moléculaires impliquée conduisant à la mort excitotoxique des cellules est la voie de la c-Jun NH2terminal kinase (JNK). Plusieurs études ont déjà montré l'efficacité d'un agent neuroprotecteur d'un nouveau type, un peptide dextrogyre synthétisé sous la forme retro inverso (XG-102, précédemment D-JNKI1) résistant aux protéases, capable de pénétrer dans les cellules et de bloquer sélectivement et fortement l'accès de JNK à plusieurs de ses cibles. Une puissante protection a été observée avec ce composé dans plusieurs modèles d'ischémie (Borsello et al. 2003;Hirt et al. 2004). Ce composé chimérique, construit à partir d'une séquence TAT de 10 acides aminés suivie par une séquence de 20 acides aminés d'un domaine liant JNK (JBD) issu de la molécule JNK protéine inhibitrice. (JIP), induit à la fois une réduction importante de la taille de lésion et un comportement fonctionnel amélioré. De plus il présente une fenêtre thérapeutique étendue. XG-102 a prouvé sa puissante efficacité dans un modèle d'occlusion de l'artère cérébrale moyenne chez la souris avec injection intracerebroventriculaire (i.c.v.) mais afin de pouvoir envisager le développement de ce composé pour l'attaque cérébrale chez l'homme, il était donc nécessaire de déterminer la faisabilité de son administration systémique. Les études faisant l'objet de cette thèse ont permis de montrer une neuroprotection importante avec XG-102 administré de façon systémique après l'occlusion transitoire de l'artère cérébrale moyenne chez la souris (MCAo). De plus nos données ont fourni des informations quant à la faisabilité de combiner XG-102 et tPA, démontrant une protection efficace par XG-102 malgré l'action nuisible du tPA sur la survie des cellules. En combinant les bénéfices de la reperfusion avec les effets d'un composé neurooprotecteur, cela représenterait l'avantage d'apporter des meilleures chances de protéger le tissu cérébral.

Usefulness of double locus sequence typing (DLST) for regional and international epidemiological surveillance of methicilin-resistant Staphylococcus aureus.

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of nosocomial infections worldwide. To differentiate reliably among S. aureus isolates, we recently developed double locus sequence typing (DLST) based on the analysis of partial sequences of clfB and spa genes. In the present study, we evaluated the usefulness of DLST for epidemiological investigations of MRSA by routinely typing 1242 strains isolated in Western Switzerland. Additionally, particular local and international collections were typed by pulsed field gel electrophoresis (PFGE) and DLST to check the compatibility of DLST with the results obtained by PFGE, and for international comparisons. Using DLST, we identified the major MRSA clones of Western Switzerland, and demonstrated the close relationship between local and international clones. The congruence of 88% between the major PFGE and DLST clones indicated that our results obtained by DLST were compatible with earlier results obtained by PFGE. DLST could thus easily be incorporated in a routine surveillance procedure. In addition, the unambiguous definition of DLST types makes this method more suitable than PFGE for long-term epidemiological surveillance. Finally, the comparison of the results obtained by DLST, multilocus sequence typing, PFGE, Staphylococcal cassette chromosome mec typing and the detection of Panton-Valentine leukocidin genes indicated that no typing scheme should be used on its own. It is only the combination of data from different methods that gives the best chance of describing precisely the epidemiology and phylogeny of MRSA.

Complex organization of CTF/NF-I, C/EBP, and HNF3 binding sites within the promoter of the liver-specific vitellogenin gene.

Vitellogenin genes are expressed specifically in the liver of female oviparous vertebrates under the strict control of estrogen. To explain this tissue-specific expression, we performed a detailed analysis of the Xenopus laevis vitellogenin gene B1 promoter by DNase I footprinting and gel mobility-shift assays. We characterized five binding sites for the ubiquitous factor CTF/NF-I. Two of these sites are close to the TATA-box, whereas the others are located on both sides of the estrogen responsive unit formed by two imperfect estrogen response elements. Moreover two liver-enriched factors, C/EBP and HNF3, were found to interact with multiple closely spaced proximal promoter elements in the first 100 base pairs upstream of the TATA-box. To confirm the physiological significance of this in vitro analysis, in vivo DNase I footprinting experiments were carried out using the ligation-mediated polymerase chain reaction technique. The various cis-elements characterized in vitro as binding sites for known transcription factors and more particularly for liver-enriched transcription factors are efficiently recognized in vivo as well, suggesting that they play an important role in the control of the liver-specific vitellogenin gene B1 expression.

The transcriptional repressor REST determines the cell-specific expression of the human MAPK8IP1 gene encoding IB1 (JIP-1).

Islet-brain 1 (IB1) is the human and rat homologue of JIP-1, a scaffold protein interacting with the c-Jun amino-terminal kinase (JNK). IB1 expression is mostly restricted to the endocrine pancreas and to the central nervous system. Herein, we explored the transcriptional mechanism responsible for this preferential islet and neuronal expression of IB1. A 731-bp fragment of the 5' regulatory region of the human MAPK8IP1 gene was isolated from a human BAC library and cloned upstream of a luciferase reporter gene. This construct drove high transcriptional activity in both insulin-secreting and neuron-like cells but not in unrelated cell lines. Sequence analysis of this promoter region revealed the presence of a neuron-restrictive silencer element (NRSE) known to bind repressor zinc finger protein REST. This factor is not expressed in insulin-secreting and neuron-like cells. By mobility shift assay, we confirmed that REST binds to the NRSE present in the IB1 promoter. Once transiently transfected in beta-cell lines, the expression vector encoding REST repressed IB1 transcriptional activity. The introduction of a mutated NRSE in the 5' regulating region of the IB1 gene abolished the repression activity driven by REST in insulin-secreting beta cells and relieved the low transcriptional activity of IB1 observed in unrelated cells. Moreover, transfection in non-beta and nonneuronal cell lines of an expression vector encoding REST lacking its transcriptional repression domain relieved IB1 promoter activity. Last, the REST-mediated repression of IB1 could be abolished by trichostatin A, indicating that deacetylase activity is required to allow REST repression. Taken together, these data establish a critical role for REST in the control of the tissue-specific expression of the human IB1 gene.

Three-dimensional sequence stratigraphy and subtle stratigraphic traps associated with systems tracts: West Cameron region, offshore Louisiana, Gulf of Mexico

Three-dimensional sequence stratigraphy is a potent exploration and development tool for the discovery of subtle stratigraphic traps. Reservoir morphology, heterogeneity and subtle stratigraphic trapping mechanisms can be better understood through systematic horizontal identification of sedimentary facies of systems tracts provided by three-dimensional attribute maps used as an important complement to the sequential analysis on the two-dimensional seismic lines and the well log data. On new prospects as well as on already-producing fields, the additional input of sequential analysis on three-dimensional data enables the identification, location and precise delimitation of new potentially productive zones. The first part of this paper presents four typical horizontal seismic facies assigned to the successive systems tracts of a third- or fourth-order sequence deposited in inner to outer neritic conditions on a elastic shelf. The construction of this synthetic representative sequence is based on the observed reproducibility of the horizontal seismic facies response to cyclic eustatic events on more than 35 sequences registered in the Gulf coast Plio-Pleistocene and Late Miocene, offshore Louisiana in the West Cameron region of the Gulf of Mexico. The second part shows how three-dimensional sequence stratigraphy can contribute in localizing and understanding sedimentary facies associated with productive zones. A case study in the early Middle Miocene Cibicides opima sands shows multiple stacked gas accumulations in the top slope fan, prograding wedge and basal transgressive systems tract of the third-order sequence between SB15.5 and SB 13.8 Ma.

Decoding sequence learning from single-trial intracranial EEG in humans.

We propose and validate a multivariate classification algorithm for characterizing changes in human intracranial electroencephalographic data (iEEG) after learning motor sequences. The algorithm is based on a Hidden Markov Model (HMM) that captures spatio-temporal properties of the iEEG at the level of single trials. Continuous intracranial iEEG was acquired during two sessions (one before and one after a night of sleep) in two patients with depth electrodes implanted in several brain areas. They performed a visuomotor sequence (serial reaction time task, SRTT) using the fingers of their non-dominant hand. Our results show that the decoding algorithm correctly classified single iEEG trials from the trained sequence as belonging to either the initial training phase (day 1, before sleep) or a later consolidated phase (day 2, after sleep), whereas it failed to do so for trials belonging to a control condition (pseudo-random sequence). Accurate single-trial classification was achieved by taking advantage of the distributed pattern of neural activity. However, across all the contacts the hippocampus contributed most significantly to the classification accuracy for both patients, and one fronto-striatal contact for one patient. Together, these human intracranial findings demonstrate that a multivariate decoding approach can detect learning-related changes at the level of single-trial iEEG. Because it allows an unbiased identification of brain sites contributing to a behavioral effect (or experimental condition) at the level of single subject, this approach could be usefully applied to assess the neural correlates of other complex cognitive functions in patients implanted with multiple electrodes.

The gene for the ligand binding chain of the human interferon gamma receptor.

In order to characterize the gene encoding the ligand binding (1(st); alpha) chain of the human IFN-gamma receptor, two overlapping cosmid clones were analyzed. The gene spans over 25 kilobases (kb) of the genomic DNA and has seven exons. The extracellular domain is encoded by exons 1 to 5 and by part of exon 6. The transmembrane region is also encoded by exon 6. Exon 7 encodes the intracellular domain and the 3' untranslated portion. The gene was located on chromosome 6q23.1, as determined by in situ hybridization. The 4 kb region upstream (5') of the gene was sequenced and analyzed for promoter activity. No consensus-matching TATA or CAAT boxes in the 5' region were found. Potential binding sites for Sp1, AP-1, AP-2, and CREB nuclear factors were identified. Compatible with the presence of the Sp1/AP-2 sites and the lack of TATA box, S1-nuclease mapping experiments showed multiple transcription initiation sites. Promoter activity of the 5' flanking region was analyzed with two different reporter genes: the Escherichia coli chloramphenicol acetyltransferase and human growth hormone. The smallest 5' region of the gene that still had full promoter activity was 692 base pairs in length. In addition, we found sequences belonging to the oldest family of Alu repeats, 2 - 3 kb upstream of the gene, which could be useful for genetic studies.

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.

Islands of euchromatin-like sequence and expressed polymorphic sequences within the short arm of human chromosome 21

The goals of the human genome project did not include sequencing of the heterochromatic regions. We describe here an initial sequence of 1.1 Mb of the short arm of human chromosome 21 (HSA21p), estimated to be 10% of 21p. This region contains extensive euchromatic-like sequence and includes on average one transcript every 100 kb. These transcripts show multiple inter- and intrachromosomal copies, and extensive copy number and sequence variability. The sequencing of the "heterochromatic" regions of the human genome is likely to reveal many additional functional elements and provide important evolutionary information.

Prominent use of distal 5’ transcription start sites and discovery of a large number of additional exons in ENCODE regions

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

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.