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.

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.

A comparison of random sequence reads versus 16S rDNA sequences for estimating the biodiversity of a metagenomic library

The construction of metagenomic libraries has permitted the study of microorganisms resistant to isolation and the analysis of 16S rDNA sequences has been used for over two decades to examine bacterial biodiversity. Here, we show that the analysis of random sequence reads (RSRs) instead of 16S is a suitable shortcut to estimate the biodiversity of a bacterial community from metagenomic libraries. We generated 10,010 RSRs from a metagenomic library of microorganisms found in human faecal samples. Then searched them using the program BLASTN against a prokaryotic sequence database to assign a taxon to each RSR. The results were compared with those obtained by screening and analysing the clones containing 16S rDNA sequences in the whole library. We found that the biodiversity observed by RSR analysis is consistent with that obtained by 16S rDNA. We also show that RSRs are suitable to compare the biodiversity between different metagenomic libraries. RSRs can thus provide a good estimate of the biodiversity of a metagenomic library and, as an alternative to 16S, this approach is both faster and cheaper.

Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

The mitochondrial targeting function of randomly generated peptide sequences correlates with predicted helical amphiphilicity.

A pool of oligonucleotides encoding a start methionine and nine random amino acids was inserted at the 5'-end of the gene for the yeast cytochrome oxidase subunit IV lacking its own mitochondrial targeting sequence. Approximately one-quarter of the randomly generated sequences targeted subunit IV to its correct intramitochondrial location in vivo. Sequence analysis of 89 randomly generated sequences showed that their efficiencies as mitochondrial targeting signals correlated with the potential to fold into an amphiphilic alpha-helix. Functional targeting sequences were enriched in arginine and isoleucine residues but contained few aspartate, glutamate, and proline residues. Nonfunctional sequences predicted to have significant helical amphiphilicity often had at least one acidic or multiple helix-breaking residues that would be expected to interfere with targeting functioning. These results support the hypothesis that the signal for targeting a protein into the mitochondrial matrix is usually a positively charged amphiphilic helix.

Prenatal exposure to traffic-related air pollution and ultrasound measures of fetal growth in the INMA Sabadell cohort

Background: Few studies have used longitudinal ultrasound measurements to assess the effect of traffic-related air pollution on fetal growth.Objective: We examined the relationship between exposure to nitrogen dioxide (NO2) and aromatic hydrocarbons [benzene, toluene, ethylbenzene, m/p-xylene, and o-xylene (BTEX)] on fetal growth assessed by 1,692 ultrasound measurements among 562 pregnant women from the Sabadell cohort of the Spanish INMA (Environment and Childhood) study.Methods: We used temporally adjusted land-use regression models to estimate exposures to NO2 and BTEX. We fitted mixed-effects models to estimate longitudinal growth curves for femur length (FL), head circumference (HC), abdominal circumference (AC), biparietal diameter (BPD), and estimated fetal weight (EFW). Unconditional and conditional SD scores were calculated at 12, 20, and 32 weeks of gestation. Sensitivity analyses were performed considering time–activity patterns during pregnancy.Results: Exposure to BTEX from early pregnancy was negatively associated with growth in BPD during weeks 20–32. None of the other fetal growth parameters were associated with exposure to air pollution during pregnancy. When considering only women who spent 2 hr/day in nonresidential outdoor locations, effect estimates were stronger and statistically significant for the association between NO2 and growth in HC during weeks 12–20 and growth in AC, BPD, and EFW during weeks 20–32.Conclusions: Our results lend some support to an effect of exposure to traffic-related air pollutants from early pregnancy on fetal growth during mid-pregnancy.

Estimation of the viscoelastic properties of vessel walls using a computational model and Doppler ultrasound

Human arteries affected by atherosclerosis are characterized by altered wall viscoelastic properties. The possibility of noninvasively assessing arterial viscoelasticity in vivo would significantly contribute to the early diagnosis and prevention of this disease. This paper presents a noniterative technique to estimate the viscoelastic parameters of a vascular wall Zener model. The approach requires the simultaneous measurement of flow variations and wall displacements, which can be provided by suitable ultrasound Doppler instruments. Viscoelastic parameters are estimated by fitting the theoretical constitutive equations to the experimental measurements using an ARMA parameter approach. The accuracy and sensitivity of the proposed method are tested using reference data generated by numerical simulations of arterial pulsation in which the physiological conditions and the viscoelastic parameters of the model can be suitably varied. The estimated values quantitatively agree with the reference values, showing that the only parameter affected by changing the physiological conditions is viscosity, whose relative error was about 27% even when a poor signal-to-noise ratio is simulated. Finally, the feasibility of the method is illustrated through three measurements made at different flow regimes on a cylindrical vessel phantom, yielding a parameter mean estimation error of 25%.

Ex-vivo ultrasound control of resection margins during partial nephrectomy

Introduction: Surgery represents the treatment of choice for localized renal cell neoplasia. Partial nephrectomy (PN) has widened its indications over the past two decades and has shown oncological results equivalent to radical nephrectomy for small tumors. The role of negative surgical margins has been widely debated. Intraoperative fresh frozen section analysis is shown to be unreliable, expensive, time-consuming and not well correlated to final pathology. The goal of the present study was to assess the feasibility of intraoperative ex-vivo ultrasound (US) control of resection margins and its correlation to margin status at definitive pathology in patients undergoing PN.Material and Methods: The study was carried out in our institution from February 2008 to March 2010. Patients undergoing PN for T1-T2 renal tumors were included. Ex vivo US was performed by one single senior radiologist. Considering its availability, not all consecutive eligible patients were included. PN was undertaken in a standardized technique applying the "minimal healthy tissue margin" technique. Once resected, the specimen was kept in a saline solution and ex-vivo US was performed to evaluate the whole tumor pseudocapsule.Results: Twelve patients (five women, age (mean}SD) 65}11 years) were included. Intraoperative ex-vivo US showed negative surgical margin in all cases. US duration ranged from 1 to 4 minutes, with a median time of 1 minute. Definitive histological analysis confirmed the presence of two angiomyolipoma, eight pT1a tumors, of which seven were clear cell carcinoma and one was a type II papillary tumor, one pT1b clear cell carcinoma and one pT2 chromophobe carcinoma (size 2.9}2.3 cm). Final pathology revealed R0 margins.Conclusion: Intraoperative ex-vivo US control of resection margins in patients undergoing PN is feasible, time-efficient and well correlated to definitive pathological examination with regards to margin status.

A trap system for the isolation of amino acid sequences inducing GPI anchor addition.

We designed a trap system to isolate different amino acid sequences which could target proteins to the cell surface via GPI anchor transfer. This selection procedure is based on the insertion of various sequences which regenerate a functional GPI anchor signal sequence and therefore provoke re-expression at the surface of a reporter molecule. Using this trap for cell surface targeting sequences, we could show the importance of the defined elements essential for GPI anchor addition. Such a system could be used for an exhaustive analysis of the carboxyl terminus structural requirements for GPI membrane anchoring.

Low estimates of dietary acid load are positively associated with bone ultrasound in women older than 75 years of age with a lifetime fracture

Dietary acid load from Western diets may be a risk factor for osteoporosis. It can be estimated by net endogenous acid production (NEAP). No data currently exists for NEAP estimates and bone indices in the very elderly (i.e. > or = 75 y). The aim of this study was to determine the association between NEAP estimates by using the potential renal acid load (PRAL) equation and quantitative bone ultrasound (QUS) measurements at the heel [broadband ultrasound attenuation (BUA)] in Caucasian women. We assessed NEAP and QUS in 401 very elderly Swiss ambulatory women. We evaluated dietary intake and NEAP estimates with a validated FFQ. QUS was measured using Achilles (Lunar). We identified 2 subgroups: 256 women (80.6 y +/- 3; BUA, 96.8 dB/MHz) with a fracture history and the remaining 145 (79.9 y SD 2.9; BUA, 101.7 dB/MHz) without. Women who reported having suffered a fracture had lower BUA (P < 0.001) than nonfractured women but did not differ in nutrient intakes and NEAP. Lower NEAP (P = 0.023) and higher potassium intake (P = 0.033) were correlated with higher BUA, which remained significant even after adjustment for age, BMI, and osteoporosis treatment. BUA was positively correlated with calcium (P = 0.016) and BMI (P < 0.001). Women who reported no fractures had no significant correlations between nutrient intake, NEAP, and BUA. Low nutritional acid load was correlated with higher BUA in very elderly women with a fracture history. Although relatively weak compared with age and BMI, this association was significant and may be an important additional risk factor that might be particularly relevant in frail patients with an already high fracture risk.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells.

BACKGROUND: Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian. RESULTS: We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types. CONCLUSIONS: The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Genometrics as an essential tool for the assembly of whole genome sequences: the example of the chromosome of Bifidobacterium longum NCC2705.

BACKGROUND: Analysis of the first reported complete genome sequence of Bifidobacterium longum NCC2705, an actinobacterium colonizing the gastrointestinal tract, uncovered its proteomic relatedness to Streptomyces coelicolor and Mycobacterium tuberculosis. However, a rapid scrutiny by genometric methods revealed a genome organization totally different from all so far sequenced high-GC Gram-positive chromosomes. RESULTS: Generally, the cumulative GC- and ORF orientation skew curves of prokaryotic genomes consist of two linear segments of opposite slope: the minimum and the maximum of the curves correspond to the origin and the terminus of chromosome replication, respectively. However, analyses of the B. longum NCC2705 chromosome yielded six, instead of two, linear segments, while its dnaA locus, usually associated with the origin of replication, was not located at the minimum of the curves. Furthermore, the coorientation of gene transcription with replication was very low. Comparison with closely related actinobacteria strongly suggested that the chromosome of B. longum was misassembled, and the identification of two pairs of relatively long homologous DNA sequences offers the possibility for an alternative genome assembly proposed here below. By genometric criteria, this configuration displays all of the characters common to bacteria, in particular to related high-GC Gram-positives. In addition, it is compatible with the partially sequenced genome of DJO10A B. longum strain. Recently, a corrected sequence of B. longum NCC2705, with a configuration similar to the one proposed here below, has been deposited in GenBank, confirming our predictions. CONCLUSION: Genometric analyses, in conjunction with standard bioinformatic tools and knowledge of bacterial chromosome architecture, represent fast and straightforward methods for the evaluation of chromosome assembly.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Improved smoking cessation in smokers given ultrasound photographs of their own atherosclerotic plaques.

BACKGROUND: We examined whether making smokers aware that they had developed peripheral atherosclerosis would improve smoking cessation. METHODS: Smokers selected from the general population were randomly allocated to undergo high-resolution B-mode ultrasonography of their carotid and femoral arteries. All smokers received quit-smoking counseling. Smokers with > or =1 atherosclerotic plaque were given two photographs of a plaque with a relevant explanation. Quit rates were assessed by telephone 6 months later. RESULTS: Seventy-nine smokers did not undergo ultrasonography (A). Among the 74 smokers submitted to ultrasonography, 20 had no plaque (B) and 54 had > or =1 plaque (C). Quit rates were, respectively, 6.3, 5.0, and 22.2% in groups A, B, and C. Quit rates were higher in smokers submitted to ultrasonography (B + C vs A; P = 0.031) and in those receiving photographs (C vs A + B; P = 0.003). Smoking cessation was independently associated with intervention C (OR = 6.2; 95% CI = 1.8-21) and a white-collar job but not with age or gender. CONCLUSIONS: Providing smokers with photographs demonstrating atherosclerosis on their own person was an effective adjunct to physician's advice to quit smoking. Since ultrasonography is used increasingly often in clinical practice for cardiovascular risk stratification, this can provide an additional opportunity and means to deter smokers from smoking.