Clinical utility of next-generation sequencing in the diagnosis of hereditary haemolytic anaemias.

Hereditary haemolytic anaemias are genetically and phenotypically heterogeneous disorders characterized by increased red cell destruction, with consequences ranging from innocuous to severe life-threatening anaemia. Diagnostic laboratories endeavour to assist clinicians reach the exact patient diagnosis by using tests principally based on morphological and biochemical techniques. However, these routine studies may be inconclusive, particularly in newborn infants and when transfusions have recently been administered. Large numbers and size of the potentially involved genes also impose a practical challenge for molecular diagnosis using routine sequencing approaches. To overcome these diagnostic shortcomings, we have utilized next-generation sequencing to provide a high-throughput, highly sensitive assay. We developed a panel interrogating 28 genes encoding cytoskeletal proteins and enzymes with sequencing coverage of the coding regions, splice site junctions, deep intronic and regulatory regions. We then evaluated 19 samples, including infants with unexplained extreme hyperbilirubinaemia and patients with transfusion-dependent haemolytic anaemia. Where possible, inheritance patterns of pathogenic mutations were determined by sequencing of immediate relatives. We conclude that this next-generation sequencing panel could be a cost-effective approach to molecular diagnosis of hereditary haemolytic anaemia, especially when the family history is uninformative or when routine laboratory testing fails to identify the causative haemolytic process.

Prophenoloxidases 1 and 2 from the spruce budworm, Choristoneura fumiferana: molecular cloning and assessment of transcriptional regulation by a polydnavirus.

Immune challenge in arthropods is frequently accompanied by melanization of the hemolymph, a reaction triggered by the activation of prophenoloxidase (PPO). Because their immature stages are spent inside the hemocoel of insect larvae, endoparasitoids have evolved strategies to escape or counter melanin formation. Very little molecular information is available on these endoparasitoid counterstrategies. We have sought to shed light on the inhibition of melanization in the spruce budworm, Choristoneura fumiferana, by the parasitic wasp Tranosema rostrale, by cloning two host PPO homologs and studying their transcriptional regulation after parasitization. The two polypeptides are encoded by transcripts of approximately 3.3 kb (for CfPPO1) and 3.0 kb (for CfPPO2) and possess structural features typical of other insect PPOs. While there appears to be a single CfPPO2 gene in the C. fumiferana genome, we detected three CfPPO1 mRNA variants displaying insertions/deletions in the 3' untranslated region, suggesting that there may be more than one CfPPO1 gene copy. Both CfPPO1 and CfPPO2 were expressed at high levels in C. fumiferana 6th instars, and parasitization by T. rostrale had no apparent impact on the level of their transcripts. Injection of a large dose (0.5 female-equivalent) of polydnavirus-laden calyx fluid extracted from T. rostrale, which is known to inhibit melanization in C. fumiferana, only caused a transient decrease in CfPPO1 and CfPPO2 transcript accumulation at 2-3 d post injection. It thus appears that transcriptional downregulation of C. fumiferana PPO by T. rostrale plays a minor role in the inhibition of hemolymph melanization in this host-parasitoid system.

Differential analysis for high density tiling microarray data.

BACKGROUND: High density oligonucleotide tiling arrays are an effective and powerful platform for conducting unbiased genome-wide studies. The ab initio probe selection method employed in tiling arrays is unbiased, and thus ensures consistent sampling across coding and non-coding regions of the genome. These arrays are being increasingly used to study the associated processes of transcription, transcription factor binding, chromatin structure and their association. Studies of differential expression and/or regulation provide critical insight into the mechanics of transcription and regulation that occurs during the developmental program of a cell. The time-course experiment, which comprises an in-vivo system and the proposed analyses, is used to determine if annotated and un-annotated portions of genome manifest coordinated differential response to the induced developmental program. RESULTS: We have proposed a novel approach, based on a piece-wise function - to analyze genome-wide differential response. This enables segmentation of the response based on protein-coding and non-coding regions; for genes the methodology also partitions differential response with a 5' versus 3' versus intra-genic bias. CONCLUSION: The algorithm built upon the framework of Significance Analysis of Microarrays, uses a generalized logic to define regions/patterns of coordinated differential change. By not adhering to the gene-centric paradigm, discordant differential expression patterns between exons and introns have been identified at a FDR of less than 12 percent. A co-localization of differential binding between RNA Polymerase II and tetra-acetylated histone has been quantified at a p-value < 0.003; it is most significant at the 5' end of genes, at a p-value < 10-13. The prototype R code has been made available as supplementary material [see Additional file 1].

Tiling array-driven elucidation of transcriptional structures based on maximum-likelihood and Markov models.

Tiling arrays of high-density oligonucleotide probes spanning the entire genome are powerful tools for the discovery of new genes. However, it is difficult to determine the structure of the spliced product of a structurally unknown gene from noisy array signals only. Here we introduce a statistical method that estimates the precise splicing points and the exon/intron structure of a structurally unknown gene by maximizing the odds or the ratio of posterior probabilities of the structure under the observation of array signal intensities and nucleic acid sequences. Our method more accurately predicted the gene structures than the simple threshold-based method, and more correctly estimated the expression values of structurally unknown genes than the window-based method. It was observed that the Markov model contributed to the precision of splice points, and that the statistical significance of expression (P-value) represented the reliability of the estimated gene structure and expression value well. We have implemented the method as a program ARTADE (ARabidopsis Tiling Array-based Detection of Exons) and applied it to the Arabidopsis thaliana whole-genome array data analysis. The database of the predicted results and the ARTADE program are available at http://omicspace.riken.jp/ARTADE/.

Annotation of cis-regulatory elements by identification, subclassification, and functional assessment of multispecies conserved sequences.

An important step toward improving the annotation of the human genome is to identify cis-acting regulatory elements from primary DNA sequence. One approach is to compare sequences from multiple, divergent species. This approach distinguishes multispecies conserved sequences (MCS) in noncoding regions from more rapidly evolving neutral DNA. Here, we have analyzed a region of approximately 238kb containing the human alpha globin cluster that was sequenced and/or annotated across the syntenic region in 22 species spanning 500 million years of evolution. Using a variety of bioinformatic approaches and correlating the results with many aspects of chromosome structure and function in this region, we were able to identify and evaluate the importance of 24 individual MCSs. This approach sensitively and accurately identified previously characterized regulatory elements but also discovered unidentified promoters, exons, splicing, and transcriptional regulatory elements. Together, these studies demonstrate an integrated approach by which to identify, subclassify, and predict the potential importance of MCSs.

AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints.

We present a WWW server for AUGUSTUS, a software for gene prediction in eukaryotic genomic sequences that is based on a generalized hidden Markov model, a probabilistic model of a sequence and its gene structure. The web server allows the user to impose constraints on the predicted gene structure. A constraint can specify the position of a splice site, a translation initiation site or a stop codon. Furthermore, it is possible to specify the position of known exons and intervals that are known to be exonic or intronic sequence. The number of constraints is arbitrary and constraints can be combined in order to pin down larger parts of the predicted gene structure. The result then is the most likely gene structure that complies with all given user constraints, if such a gene structure exists. The specification of constraints is useful when part of the gene structure is known, e.g. by expressed sequence tag or protein sequence alignments, or if the user wants to change the default prediction. The web interface and the downloadable stand-alone program are available free of charge at http://augustus.gobics.de/submission.

SLAM web server for comparative gene finding and alignment.

SLAM is a program that simultaneously aligns and annotates pairs of homologous sequences. The SLAM web server integrates SLAM with repeat masking tools and the AVID alignment program to allow for rapid alignment and gene prediction in user submitted sequences. Along with annotations and alignments for the submitted sequences, users obtain a list of predicted conserved non-coding sequences (and their associated alignments). The web site also links to whole genome annotations of the human, mouse and rat genomes produced with the SLAM program. The server can be accessed at http://bio.math.berkeley.edu/slam.