Predicting protein secondary structure with a nearest-neighbor algorithm.

We have developed a new method for protein secondary structure prediction that achieves accuracies as high as 71.0%, the highest value yet reported. The main component of our method is a nearest-neighbor algorithm that uses a more sophisticated treatment of the feature space than standard nearest-neighbor methods. It calculates distance tables that allow it to produce real-valued distances between amino acid residues, and attaches weights to the instances to further modify the the structure of feature space. The algorithm, which is closely related to the memory-based reasoning method of Zhang et al., is simple and easy to train, and has also been applied with excellent results to the problem of identifying DNA promoter sequences.     

Genome re-annotation: a wiki solution?

The annotation of most genomes becomes outdated over time, owing in part to our ever-improving knowledge of genomes and in part to improvements in bioinformatics software. Unfortunately, annotation is rarely if ever updated and resources to support routine reannotation are scarce. Wiki software, which would allow many scientists to edit each genome's annotation, offers one possible solution.     

Quake: quality-aware detection and correction of sequencing errors

We introduce Quake, a program to detect and correct errors in DNA sequencing reads. Using a maximum likelihood approach incorporating quality values and nucleotide specific miscall rates, Quake achieves the highest accuracy on realistically simulated reads. We further demonstrate substantial improvements in de novo assembly and SNP detection after using Quake. Quake can be used for any size project, including more than one billion human reads, and is freely available as open source software from .     

<Emphasis Type="Italic">α</Emphasis>-Amylase inhibitor-1 gene from <Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> expressed in <Emphasis Type="Italic">Coffea arabica</Emphasis>plants inhibits α-amylases from the coffee berry borer pest

Background  

Coffee is an important crop and is crucial to the economy of many developing countries, generating around US70 billion per year. There are 115 species in the < i > Coffea < /i > genus, but only two, < i > C. arabica < /i > and < i > C. canephora < /i > , are commercially cultivated. Coffee plants are attacked by many pathogens and insect-pests, which affect not only the production of coffee but also its grain quality, reducing the commercial value of the product. The main insect-pest, the coffee berry borer ( < i > Hypotheneumus hampei < /i > ), is responsible for worldwide annual losses of around US70 billion per year. There are 115 species in the Coffea genus, but only two, C. arabica and C. canephora, are commercially cultivated. Coffee plants are attacked by many pathogens and insect-pests, which affect not only the production of coffee but also its grain quality, reducing the commercial value of the product. The main insect-pest, the coffee berry borer (Hypotheneumus hampei), is responsible for worldwide annual losses of around US500 million. The coffee berry borer exclusively damages the coffee berries, and it is mainly controlled by organochlorine insecticides that are both toxic and carcinogenic. Unfortunately, natural resistance in the genus Coffea to H. hampei has not been documented. To overcome these problems, biotechnological strategies can be used to introduce an α-amylase inhibitor gene (α-AI1), which confers resistance against the coffee berry borer insect-pest, into C. arabica plants.     

Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote

The ciliate Tetrahymena thermophila is a model organism for molecular and cellular biology. Like other ciliates, this species has separate germline and soma functions that are embodied by distinct nuclei within a single cell. The germline-like micronucleus (MIC) has its genome held in reserve for sexual reproduction. The soma-like macronucleus (MAC), which possesses a genome processed from that of the MIC, is the center of gene expression and does not directly contribute DNA to sexual progeny. We report here the shotgun sequencing, assembly, and analysis of the MAC genome of T. thermophila, which is approximately 104 Mb in length and composed of approximately 225 chromosomes. Overall, the gene set is robust, with more than 27,000 predicted protein-coding genes, 15,000 of which have strong matches to genes in other organisms. The functional diversity encoded by these genes is substantial and reflects the complexity of processes required for a free-living, predatory, single-celled organism. This is highlighted by the abundance of lineage-specific duplications of genes with predicted roles in sensing and responding to environmental conditions (e.g., kinases), using diverse resources (e.g., proteases and transporters), and generating structural complexity (e.g., kinesins and dyneins). In contrast to the other lineages of alveolates (apicomplexans and dinoflagellates), no compelling evidence could be found for plastid-derived genes in the genome. UGA, the only T. thermophila stop codon, is used in some genes to encode selenocysteine, thus making this organism the first known with the potential to translate all 64 codons in nuclear genes into amino acids. We present genomic evidence supporting the hypothesis that the excision of DNA from the MIC to generate the MAC specifically targets foreign DNA as a form of genome self-defense. The combination of the genome sequence, the functional diversity encoded therein, and the presence of some pathways missing from other model organisms makes T. thermophila an ideal model for functional genomic studies to address biological, biomedical, and biotechnological questions of fundamental importance.     

JIGSAW,GeneZilla, and GlimmerHMM: puzzling out the features of human genes in the ENCODE regions

Background

Predicting complete protein-coding genes in human DNA remains a significant challenge. Though a number of promising approaches have been investigated, an ideal suite of tools has yet to emerge that can provide near perfect levels of sensitivity and specificity at the level of whole genes. As an incremental step in this direction, it is hoped that controlled gene finding experiments in the ENCODE regions will provide a more accurate view of the relative benefits of different strategies for modeling and predicting gene structures.

Results

Here we describe our general-purpose eukaryotic gene finding pipeline and its major components, as well as the methodological adaptations that we found necessary in accommodating human DNA in our pipeline, noting that a similar level of effort may be necessary by ourselves and others with similar pipelines whenever a new class of genomes is presented to the community for analysis. We also describe a number of controlled experiments involving the differential inclusion of various types of evidence and feature states into our models and the resulting impact these variations have had on predictive accuracy.

Conclusion

While in the case of the non-comparative gene finders we found that adding model states to represent specific biological features did little to enhance predictive accuracy, for our evidence-based 'combiner' program the incorporation of additional evidence tracks tended to produce significant gains in accuracy for most evidence types, suggesting that improved modeling efforts at the hidden Markov model level are of relatively little value. We relate these findings to our current plans for future research.

     

Antizyme is a target of sex-lethal in the Drosophila germline and appears to act downstream of hedgehog to regulate sex-lethal and cyclin B

The sex determination master switch, Sex-lethal, has been shown to regulate the mitosis of early germ cells in Drosophila melanogaster. Sex-lethal is an RNA binding protein that regulates splicing and translation of specific targets in the soma, but the germline targets are unknown. In an experiment aimed at identifying targets of Sex-lethal in early germ cells, the RNA encoded by gutfeeling, the Drosophila homolog of Ornithine Decarboxylase Antizyme, was isolated. gutfeeling interacts genetically with Sex-lethal. It is not only a target of Sex-lethal, but also appears to regulate the nuclear entry and overall levels of Sex-lethal in early germ cells. This regulation of Sex-lethal by gutfeeling appears to occur downstream of the Hedgehog signal. We also show that Hedgehog, Gutfeeling, and Sex-lethal function to regulate Cyclin B, providing a link between Sex-lethal and mitosis.     

Genome sequence of Chlamydophila caviae (Chlamydia psittaci GPIC): examining the role of niche-specific genes in the evolution of the Chlamydiaceae

The genome of Chlamydophila caviae (formerly Chlamydia psittaci, GPIC isolate) (1 173 390 nt with a plasmid of 7966 nt) was determined, representing the fourth species with a complete genome sequence from the Chlamydiaceae family of obligate intracellular bacterial pathogens. Of 1009 annotated genes, 798 were conserved in all three other completed Chlamydiaceae genomes. The C.caviae genome contains 68 genes that lack orthologs in any other completed chlamydial genomes, including tryptophan and thiamine biosynthesis determinants and a ribose-phosphate pyrophosphokinase, the product of the prsA gene. Notable amongst these was a novel member of the virulence-associated invasin/intimin family (IIF) of Gram-negative bacteria. Intriguingly, two authentic frameshift mutations in the ORF indicate that this gene is not functional. Many of the unique genes are found in the replication termination region (RTR or plasticity zone), an area of frequent symmetrical inversion events around the replication terminus shown to be a hotspot for genome variation in previous genome sequencing studies. In C.caviae, the RTR includes several loci of particular interest including a large toxin gene and evidence of ancestral insertion(s) of a bacteriophage. This toxin gene, not present in Chlamydia pneumoniae, is a member of the YopT effector family of type III-secreted cysteine proteases. One gene cluster (guaBA-add) in the RTR is much more similar to orthologs in Chlamydia muridarum than those in the phylogenetically closest species C.pneumoniae, suggesting the possibility of horizontal transfer of genes between the rodent-associated Chlamydiae. With most genes observed in the other chlamydial genomes represented, C.caviae provides a good model for the Chlamydiaceae and a point of comparison against the human atherosclerosis-associated C.pneumoniae. This crucial addition to the set of completed Chlamydiaceae genome sequences is enabling dissection of the roles played by niche-specific genes in these important bacterial pathogens.