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.

Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves.

Powdery mildew is an important disease of wheat caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici. This pathogen invades exclusively epidermal cells after penetrating directly through the cell wall. Because powdery mildew colonizes exclusively epidermal cells, it is of importance not only to identify genes which are activated, but also to monitor tissue specificity of gene activation. Acquired resistance of wheat to powdery mildew can be induced by a previous inoculation with the non-host pathogen B. graminis f. sp. hordei, the causal agent of barley powdery mildew. The establishment of the resistant state is accompanied by the activation of genes. Here we report the tissue-specific cDNA-AFLP analysis and cloning of transcripts accumulating 6 and 24 h after the resistance-inducing inoculation with B. graminis f. sp. hordei. A total of 25,000 fragments estimated to represent about 17,000 transcripts were displayed. Out of these, 141 transcripts, were found to accumulate after Bgh inoculation using microarray hybridization analysis. Forty-four accumulated predominantly in the epidermis whereas 76 transcripts accumulated mostly in mesophyll tissue.

The loss of GLUT2 expression in the pancreatic beta-cells of diabetic db/db mice is associated with an impaired DNA-binding activity of islet-specific trans-acting factors.

GLUT2 expression is reduced in the pancreatic beta-cells of several diabetic animals. The transcriptional control of the gene in beta-cells involves at least two islet-specific DNA-binding proteins, GTIIa and PDX-1, which also transactivates the insulin, somatostatin and glucokinase genes. In this report, we assessed the DNA-binding activities of GTIIa and PDX-1 to their respective cis-elements of the GLUT2 promoter using nuclear extracts prepared from pancreatic islets of 12 week old db/db diabetic mice. We show that the decreased GLUT2 mRNA expression correlates with a decrease of the GTIIa DNA-binding activity, whereas the PDX-1 binding activity is increased. In these diabetic animals, insulin mRNA expression remains normal. The adjunction of dexamethasone to isolated pancreatic islets, a treatment previously shown to decrease PDX-1 expression in the insulin-secreting HIT-T15 cells, has no effect on the GTIIa and PDX-1 DNA-binding activities. These data suggest that the decreased activity of GTIIa, in contrast to PDX-1, may be a major initial step in the development of the beta-cell dysfunction in this model of diabetes.

FANCM regulates DNA chain elongation and is stabilized by S-phase checkpoint signalling.

FANCM binds and remodels replication fork structures in vitro. We report that in vivo, FANCM controls DNA chain elongation in an ATPase-dependent manner. In the presence of replication inhibitors that do not damage DNA, FANCM counteracts fork movement, possibly by remodelling fork structures. Conversely, through damaged DNA, FANCM promotes replication and recovers stalled forks. Hence, the impact of FANCM on fork progression depends on the underlying hindrance. We further report that signalling through the checkpoint effector kinase Chk1 prevents FANCM from degradation by the proteasome after exposure to DNA damage. FANCM also acts in a feedback loop to stabilize Chk1. We propose that FANCM is a ringmaster in the response to replication stress by physically altering replication fork structures and by providing a tight link to S-phase checkpoint signalling.

Reassessing the role of mitochondrial DNA mutations in autism spectrum disorder

Background: There is increasing evidence that impairment of mitochondrial energy metabolism plays an important role in the pathophysiology of autism spectrum disorders (ASD; OMIM number: 209850). A significant proportion of ASD cases display biochemical alterations suggestive of mitochondrial dysfunction and several studies have reported that mutations in the mitochondrial DNA (mtDNA) molecule could be involved in the disease phenotype. Methods: We analysed a cohort of 148 patients with idiopathic ASD for a number of mutations proposed in the literature as pathogenic in ASD. We also carried out a case control association study for the most common European haplogroups (hgs) and their diagnostic single nucleotide polymorphisms (SNPs) by comparing cases with 753 healthy and ethnically matched controls.Results: We did not find statistical support for an association between mtDNA mutations or polymorphisms and ASD.Conclusions: Our results are compatible with the idea that mtDNA mutations are not a relevant cause of ASD and the frequent observation of concomitant mitochondrial dysfunction and ASD could be due to nuclear factors influencing mitochondrion functions or to a more complex interplay between the nucleus and the mitochondrion/mtDNA.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

Structure and expression profile of the Arabidopsis PHO1 gene family indicates a broad role in inorganic phosphate homeostasis.

PHO1 has been recently identified as a protein involved in the loading of inorganic phosphate into the xylem of roots in Arabidopsis. The genome of Arabidopsis contains 11 members of the PHO1 gene family. The cDNAs of all PHO1 homologs have been cloned and sequenced. All proteins have the same topology and harbor a SPX tripartite domain in the N-terminal hydrophilic portion and an EXS domain in the C-terminal hydrophobic portion. The SPX and EXS domains have been identified in yeast (Saccharomyces cerevisiae) proteins involved in either phosphate transport or sensing or in sorting proteins to endomembranes. The Arabidopsis genome contains additional proteins of unknown function containing either a SPX or an EXS domain. Phylogenetic analysis indicated that the PHO1 family is subdivided into at least three clusters. Reverse transcription-PCR revealed a broad pattern of expression in leaves, roots, stems, and flowers for most genes, although two genes are expressed exclusively in flowers. Analysis of the activity of the promoter of all PHO1 homologs using promoter-beta-glucuronidase fusions revealed a predominant expression in the vascular tissues of roots, leaves, stems, or flowers. beta-Glucuronidase expression is also detected for several promoters in nonvascular tissue, including hydathodes, trichomes, root tip, root cortical/epidermal cells, and pollen grains. The expression pattern of PHO1 homologs indicates a likely role of the PHO1 proteins not only in the transfer of phosphate to the vascular cylinder of various tissues but also in the acquisition of phosphate into cells, such as pollen or root epidermal/cortical cells.

Study of common database feeding with results coming from different analytical methods in the framework of illicit drugs chemical profiling

Analytical results harmonisation is investigated in this study to provide an alternative to the restrictive approach of analytical methods harmonisation which is recommended nowadays for making possible the exchange of information and then for supporting the fight against illicit drugs trafficking. Indeed, the main goal of this study is to demonstrate that a common database can be fed by a range of different analytical methods, whatever the differences in levels of analytical parameters between these latter ones. For this purpose, a methodology making possible the estimation and even the optimisation of results similarity coming from different analytical methods was then developed. In particular, the possibility to introduce chemical profiles obtained with Fast GC-FID in a GC-MS database is studied in this paper. By the use of the methodology, the similarity of results coming from different analytical methods can be objectively assessed and the utility in practice of database sharing by these methods can be evaluated, depending on profiling purposes (evidential vs. operational perspective tool). This methodology can be regarded as a relevant approach for database feeding by different analytical methods and puts in doubt the necessity to analyse all illicit drugs seizures in one single laboratory or to implement analytical methods harmonisation in each participating laboratory.

Forensic identification of urine samples: a comparison between nuclear and mitochondrial DNA markers.

Urine samples from 20 male volunteers of European Caucasian origin were stored at 4 degrees C over a 4-month period in order to compare the identification potential of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) markers. The amount of nDNA recovered from urines dramatically declined over time. Consequently, nDNA likelihood ratios (LRs) greater than 1,000 were obtained for 100, 70 and 55% of the urines analysed after 6, 60 and 120 days, respectively. For the mtDNA, HVI and HVII sequences were obtained for all samples tested, whatever the period considered. Nevertheless, the highest mtDNA LR of 435 was relatively low compared to its nDNA equivalent. Indeed, LRs obtained with only three nDNA loci could easily exceed this value and are quite easier to obtain. Overall, the joint use of nDNA and mtDNA markers enabled the 20 urine samples to be identified, even after the 4-month period.

Viral transport of DNA damage that mimics a stalled replication fork.

Adeno-associated virus type 2 (AAV2) infection incites cells to arrest with 4N DNA content or die if the p53 pathway is defective. This arrest depends on AAV2 DNA, which is single stranded with inverted terminal repeats that serve as primers during viral DNA replication. Here, we show that AAV2 DNA triggers damage signaling that resembles the response to an aberrant cellular DNA replication fork. UV treatment of AAV2 enhances the G2 arrest by generating intrastrand DNA cross-links which persist in infected cells, disrupting viral DNA replication and maintaining the viral DNA in the single-stranded form. In cells, such DNA accumulates into nuclear foci with a signaling apparatus that involves DNA polymerase delta, ATR, TopBP1, RPA, and the Rad9/Rad1/Hus1 complex but not ATM or NBS1. Focus formation and damage signaling strictly depend on ATR and Chk1 functions. Activation of the Chk1 effector kinase leads to the virus-induced G2 arrest. AAV2 provides a novel way to study the cellular response to abnormal DNA replication without damaging cellular DNA. By using the AAV2 system, we show that in human cells activation of phosphorylation of Chk1 depends on TopBP1 and that it is a prerequisite for the appearance of DNA damage foci.

EV02: a Phase I trial to compare the safety and immunogenicity of HIV DNA-C prime-NYVAC-C boost to NYVAC-C alone.

The aim of this randomised controlled trial was to see if the addition of 4 mg/ml DNA-C priming given by the intramuscular route at weeks 0 and 4 to NYVAC-C at weeks 20 and 24, safely increased the proportion of participants with HIV-specific T-cell responses measured by the interferon (IFN)-gamma ELISpot assay at weeks 26 and/or 28 compared to NYVAC-C alone. Although 2 individuals discontinued after the first DNA-C due to adverse events (1 vaso-vagal; 1 transient, asymptomatic elevation in alanine transaminase), the vaccines were well tolerated. Three others failed to complete the regimen (1 changed her mind; 2 lost to follow-up). Of the 35 that completed the regimen 90% (18/20) in the DNA-C group had ELISpot responses compared to 33% (5/15) that received NYVAC-C alone (p=0.001). Responses were to envelope in the majority (21/23). Of the 9 individuals with responses to envelope and other peptides, 8 were in the DNA-C group. These promising results suggest that DNA-C was an effective priming agent, that merits further investigation.

Nonshivering thermogenesis capacity associated to mitochondrial DNA haplotypes and gender in the greater white-toothed shrew, Crocidura russula.

A selection gradient was recently suggested as one possible cause for a clinal distribution of mitochondrial DNA (mtDNA) haplotypes along an altitudinal transect in the greater white-toothed shrew, Crocidura russula (Ehinger et al. 2002). One mtDNA haplotype (H1) rare in lowland, became widespread when approaching the altitudinal margin of the distribution. As H1 differs from the main lowland haplotype by several nonsynonymous mutations (including on ATP6), and as mitochondria play a crucial role in metabolism and thermogenesis, distribution patterns might stem from differences in the thermogenic capacity of different mtDNA haplotypes. In order to test this hypothesis, we measured the nonshivering thermogenesis (NST) associated with different mtDNA haplotypes. Sixty-two shrews, half of which had the H1 haplotype, were acclimated in November at semioutdoor conditions and measured for NST throughout winter. Our results showed the crucial role of NST for winter survival in C. russula. The individuals that survived winter displayed a higher significant increase in NST during acclimation, associated with a significant gain in body mass, presumably from brown fat accumulation. The NST capacity (ratio of NST to basal metabolic rate) was exceptionally high for such a small species. NST was significantly affected by a gender x haplotype interaction after winter-acclimation: females bearing the H1 haplotype displayed a better thermogenesis at the onset of the breeding season, while the reverse was true for males. Altogether, our results suggest a sexually antagonistic cyto-nuclear selection on thermogenesis.

Tissue biomarker development in a multicentre trial context: a feasibility study on the PETACC3 stage II and III colon cancer adjuvant treatment trial.

PURPOSE: We evaluated the feasibility of biomarker development in the context of multicenter clinical trials. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded (FFPE) tissue samples were collected from a prospective adjuvant colon cancer trial (PETACC3). DNA was isolated from tumor as well as normal tissue and used for analysis of microsatellite instability, KRAS and BRAF genotyping, UGT1A1 genotyping, and loss of heterozygosity of 18 q loci. Immunohistochemistry was used to test expression of TERT, SMAD4, p53, and TYMS. Messenger RNA was retrieved and tested for use in expression profiling experiments. RESULTS: Of the 3,278 patients entered in the study, FFPE blocks were obtained from 1,564 patients coming from 368 different centers in 31 countries. In over 95% of the samples, genomic DNA tests yielded a reliable result. Of the immmunohistochemical tests, p53 and SMAD4 staining did best with reliable results in over 85% of the cases. TERT was the most problematic test with 46% of failures, mostly due to insufficient tissue processing quality. Good quality mRNA was obtained, usable in expression profiling experiments. CONCLUSIONS: Prospective clinical trials can be used as framework for biomarker development using routinely processed FFPE tissues. Our results support the notion that as a rule, translational studies based on FFPE should be included in prospective clinical trials.

Traducción comentada de “Genetically Modified Organisms (GMOs): Transgenic Crops and Recombinant DNA Technology”

The purpose of this paper is to provide a translation into Spanish of a review articleabout genetically modified organisms (GMOs) entitled “Genetically ModifiedOrganisms (GMOs): Transgenic Crops and Recombinant DNA Technology” publishedby the well-known scientific journal Nature. In a world where English has become thelingua franca when it comes to transferring scientific knowledge and information, itmust be taken into account that not everyone—from scientist to the general public—hasa good enough command of English so that they can feel comfortable enough reading inthis language. Translators are consequently needed resulting from a great demand oftranslation activity into, for example, Spanish. This is the reason why the proposedSpanish translation is followed by a detailed analysis emphasizing the difficulties andproblems that characterize scientific—and also general—translation (i.e. terminology,syntax, semantics, pragmatics, and ideology), for which different approaches as how tosolve them are provided. On the basis of the analysis, it can be concluded thatexperience will be of much help to scientific translators, given that specificterminological knowledge and style requirements must always be born in mind whentranslating in this field. Moreover, this paper is intended to serve as a guide forTranslation students specializing in the field of science and the expectation is to helpthem make the right decisions when it comes to translating. However, it is clear that itcan only be thought of as an introduction that should be completed with further researchand documentation tasks in order to offer a complete reference tool: the ultimatehandbook of scientific translation.