Median network analysis of defectively sequenced entire mitochondrial genomes from early and contemporary disease studies

Sequence analysis of the mitochondrial genome has become a routine method in the study of mitochondrial diseases. Quite often, the sequencing efforts in the search of pathogenic or disease-associated mutations are affected by technical and interpretive pr

BarraCUDA – a Fast Sequence Mapping Software using Graphics Processing Units

High-throughput DNA sequencing (HTS) instruments today are capable of generating millions of sequencing reads in a short period of time, and this represents a serious challenge to current bioinformatics pipeline in processing such an enormous amount of data in a fast and economical fashion. Modern graphics cards are powerful processing units that consist of hundreds of scalar processors in parallel in order to handle the rendering of high-definition graphics in real-time. It is this computational capability that we propose to harness in order to accelerate some of the time-consuming steps in analyzing data generated by the HTS instruments. We have developed BarraCUDA, a novel sequence mapping software that utilizes the parallelism of NVIDIA CUDA graphics cards to map sequencing reads to a particular location on a reference genome. While delivering a similar mapping fidelity as other mainstream programs , BarraCUDA is a magnitude faster in mapping throughput compared to its CPU counterparts. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the mapping throughput. BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the mapping of millions of sequencing reads generated by HTS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology. BarraCUDA is currently available at http://seqbarracuda.sf.net

BarraCUDA - a fast short read sequence aligner using graphics processing units.

BACKGROUND: With the maturation of next-generation DNA sequencing (NGS) technologies, the throughput of DNA sequencing reads has soared to over 600 gigabases from a single instrument run. General purpose computing on graphics processing units (GPGPU), extracts the computing power from hundreds of parallel stream processors within graphics processing cores and provides a cost-effective and energy efficient alternative to traditional high-performance computing (HPC) clusters. In this article, we describe the implementation of BarraCUDA, a GPGPU sequence alignment software that is based on BWA, to accelerate the alignment of sequencing reads generated by these instruments to a reference DNA sequence. FINDINGS: Using the NVIDIA Compute Unified Device Architecture (CUDA) software development environment, we ported the most computational-intensive alignment component of BWA to GPU to take advantage of the massive parallelism. As a result, BarraCUDA offers a magnitude of performance boost in alignment throughput when compared to a CPU core while delivering the same level of alignment fidelity. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the alignment throughput. CONCLUSIONS: BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the alignment of millions of sequencing reads generated by NGS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology.BarraCUDA is currently available from http://seqbarracuda.sf.net.

Microsatellite marker isolation and cultured strain identification in Carassius auratus gibelio

Microsatellites have become the preferred molecular markers for strain selection and genetic breeding in fish. In this study a total of 105 microsatellites were isolated and identified in gibel carp (Carassius auratus gibelio) by microsatellite sequence searches in GenBank and other databases and by screening and sequencing of positive clones from the genomic library enriched for AG and GATA repeats. Moreover, nineteen microsatellites were randomly selected to design locus-specific primer pairs, and these were successfully used to identify and discriminate different cultured strains of gibel carp including strains A, D, L, and F. Three different types of microsatellite pattern were distinguished by the number and length of fragments amplified from the 19 primer pairs, and some microsatellite primer pairs were found to produce different microsatellite patterns among strains and strain-specific fragments. In addition, some duplicated alleles were also detected in two microsatellite patterns. Therefore, the current study provides direct molecular markers to discriminate among different cultured strains for selective breeding and aquaculture practice of gibel carp.

Molecular evolution of connective tissue growth factor in Cyprinidae (Teleostei : Cypriniformes)

Connective tissue growth factor (CTGF) plays an important role in regulation of cell growth, differentiation, apoptosis and individual development in animals. The study of sequences variation and molecular evolution of CTGF gene across various species of the cyprinid could be helpful for understanding of speciation and gene divergence in this kind of fish. In this study, 19 novel sequences of CTGF gene were obtained from the representative species of the family Cyprinidae using PCR amplification, cloning and sequencing. Phylogenetic relationships of Cyprinidae were reconstructed by neighbor-joining (NJ), maximum parsimony (MP), maximum likelihood (ML) and Bayesian method. Oryzias latipes from the family Cyprinodontidae was assigned to be the outgroup taxon. Leuciscini and Barbini were clustered into the monophyletic lineages, respectively, with the high nodal supports. The estimation of the ratio of non-synonymous to synonymous substitution (dN/dS) for the various branches indicated that there stood the different evolution rates between the Leuciscini and the Barbini. With the ratio of dN/dS of the Leuciscini being lower than that of the Barbini, species within the Barbini were demonstrated to be subjected to the relatively less selection pressure and under the relaxable evolution background. A 6 by indel (insertion/deletion) was found at the 5' end of CTGF gene of Cyprinidae, and this 6 by deletion only appeared in the Leuciscini, which is a typical characteristic of the Leuciscini and provides evidence for the monophylogeny of the Leuciscini. For the amino acid sequences of CTGF protein, the most variations and indels were distributed in the signal region and IGFBP region of this protein, implying that these variations were correlated with the regulation of the CTGF gene expression and protein activity. (c) 2007 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Phylogenetic relationships of Cyprinidae (Teleostei : Cypriniformes) inferred from the partial S6K1 gene sequences and implication of indel sites in intron 1

The family Cyprinidae is widely distributed in East Asia, and has the important phylogenetic significance in the fish evolution. In this study, the 5' end partial sequences (containing exon 1, exon 2 and indel 1) of S6K1 gene were obtained from 30 representative species in Cyprinidae and outgroup using PCR amplification and sequencing. The phylogenetic relationships of Cyprinidae were reconstructed with neighbor joining (NJ), maximum parsimony (MP), maximum likelihood (ML), and Bayesian methods. Myxocyprinus asiaticus (Catostomidae) was assigned to the outgroup taxon. Similar phylogenetic relationships within the family Cyprinidae were achieved with the four analyses. Leuciscini and Barbini were monophyletic lineages respectively with the high nodal supports. Leuciscini comprises Hypophthalmichthyinae, Xenocyprinae, Cultrinae, Gobioninae, Acheilognathinae and East Asian species of Leuciscinae and Danioninae. Monophyly of East Asian clade was supported with high nodal support. Barbini comprises Schizothoracinae, Barbinae, Cyprininae and Labeoninae. The monophyletic lineage consisting of Danio rerio, D. myersi, and Rasbora trilineata was basal in the tree. In addition, the large fragment indels in intron 1 were analyzed to improve the understanding of Cyprinidae relationships. The results showed that the large fragment indels were correlated with the relations among species. Some conserved regions in intron 1 were thought to be involved in the functional regulation. However, no correlation was found between sequence variations and species characteristic size.

Site-directed gene integration in transgenic zebrafish mediated by cre recombinase using a combination of mutant Lox sites

With current gene-transfer techniques in fish, insertion of DNA into the genome occurs randomly and in many instances at multiple sites. Associated position effects, copy number differences, and multiple gene interactions make gene expression experiments difficult to interpret and fish phenotype less predictable. To meet different fish engineering needs, we describe here a gene targeting model in zebrafish. At first, four target zebrafish lines, each harboring a single genomic lox71 target site, were generated by zebrafish transgenesis. The zygotes of transgenic zebrafish lines were coinjected with capped Cre mRNA and a knockin vector pZklox66RFP. Site-specific integration event happened from one target zebrafish line. In this line two integrant zebrafish were obtained from more than 80,000 targeted embryos (integrating efficiency about 10(-4) to 10(-5)) and confirmed to have a sole copy of the integrating DNA at the target genome site. Genomic polymerase chain reaction analysis and DNA sequencing verified the correct gene target events where lox71 and lox66 have accurately recombined into double mutant lox72 and wild-type loxP. Each integrant zebrafish chosen for analysis harbored the transgene rfp at the designated egfp concatenates. Although the Cre-mediated recombination is site specific, it is dependent on a randomly placed target site. That is, a genomic target cannot be preselected for integration based solely on its sequence. Conclusively, an rfp reporter gene was successfully inserted into the egfp target locus of zebrafish genome by Cre-lox-mediated recombination. This site-directed knockin system using the lox71/lox66 combination should be a promising gene-targeting platform serving various purposes in fish genetic engineering.

Characterization of complete genome sequence of the spring viremia of carp virus isolated from common carp (Cyprinus carpio) in China

The complete genome of spring viraemia of carp virus (SVCV) strain A-1 isolated from cultured common carp (Cyprinus carpio) in China was sequenced and characterized. Reverse transcription-polymerase chain reaction (RT-PCR) derived clones were constructed and the DNA was sequenced. It showed that the entire genome of SVCV A-1 consists of 11,100 nucleotide base pairs, the predicted size of the viral RNA of rhabdoviruses. However, the additional insertions in bp 4633-4676 and bp 4684-4724 of SVCV A-1 were different from the other two published SVCV complete genomes. Five open reading frames (ORFs) of SVCV A-1 were identified and further confirmed by RT-PCR and DNA sequencing of their respective RT-PCR products. The 5 structural proteins encoded by the viral RNA were ordered 3'-N-P-M-G-L-5'. This is the first report of a complete genome sequence of SVCV isolated from cultured carp in China. Phylogenetic analysis indicates that SVCV A-1 is closely related to the members of the genus Vesiculovirus, family Rhabdoviridae.

Construction and characterization of a Lipotes vexillifer genomic DNA BAC library

We constructed a genomic DNA library for Lipotes vexillifer (L. vexillifer), the Baiji or Yangtze River dolphin, one of the most endangered mammals in the world. The library consists of 149,000 BAC clones, with an average insert size of 83 kb, representing approximately 3.4 haploid genome equivalents. PCR amplification of four known L. vexillifer genes yielded two to four positive clones each. To demonstrate the utility of this library, we isolated and sequenced the L. vexillifer alpha lactalbumin gene, which is a gene specific to mammals and one which has been widely used as molecular tool in phylogenetic analysis. We also end-sequenced 20 randomly selected clones, resulting in the identification of at least five new L. vexilliter genes, five SSR loci, and one SINE locus. These results suggest that this library is a valuable resource for candidate gene cloning, physical mapping, and genome sequencing of this important and threatened species.

Isolation and characterization of a rhabdovirus from co-infection of two viruses in mandarin fish

Co-infection of two viruses has been observed in mandarin fish (Siniperca chuatsi), but the two viruses have not been characterized. In this study, a rhabdovirus has been isolated from the co-infected two viruses extracted from the diseased mandarin fish, and its morphological structure and partial biochemical and biophysical characteristics have been observed and analyzed. The isolated rhabdovirus has a typical bullet shape, and is therefore called S. chttatsi rhabdovirus (SCRV). And, the isolated rhabdovirus produced a higher titer (10(8.5) TCID50 ml(-1)) than did the co-infecting viruses (10(6.5) TCID50 ml(-1)). Subsequently, the viral genome RNA was extracted, and used as template to clone the complete nucleoprotein (N) gene by RT-PCR amplification. Cloning and sequencing of the SCRV N protein revealed 42%-31% amino acid identities to that of trout rhabdovirus 903/87 and the rhabdoviruses in genus Vesiculovirus. SDS-PAGE separation of the isolated SCRV and other two rhabdoviruses also revealed obvious polypeptide profile difference. Moreover, the anti-SCRV N protein antibody was prepared, and the anti-SCRV N protein antibody only could recognize the SCRV N protein, whereas no antigenicity was detected in other two rhabdoviruses. The data suggested that the SCRV should be a rhabdovirus member related to the genus Vesiculovirus in the Rhabdoviridae. (c) 2006 Elsevier B.V. All rights reserved.

Molecular cloning of growth hormone receptor (GHR) from common carp (Cyprinus carpio L.) and identification of its two forms of mRNA transcripts

The cDNA of growth hormone receptor (GHR) was cloned from the liver of 2-year common carp (Cyprinus carpio L.) by reverse transcription-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA end (RACE). Its open reading frame (ORF) of 1806 nucleotides is translated into a putative peptide of 602 amino acids, including an extracellular ligand-binding domain of 244 amino acids (aa), a single transmembrane domain of 24 aa and an intracellular signal-transduction domain of 334 aa. Sequence analysis indicated that common carp GHR is highly homologous to goldfish (Carassius auratus) GHR at both gene and protein levels. Using a pair of gene-specific primers, a GHR fragment was amplified from the cDNA of 2-year common carp, a 224 bp product was identified in liver and a 321 bp product in other tissues. The sequencing of the products and the partial genomic DNA indicated that the difference in product size was the result of a 97 bp intron that alternatively spliced. In addition, the 321 bp fragment could be amplified from all the tissues of 4-month common carp including liver, demonstrating the occurrence of the alternative splicing of this intron during the development of common carp. Moreover, a semi-quantitative RT-PCR was performed to analyze the expression level of GHR in tissues of 2-year common carp and 4-month common carp. The result revealed that in the tissues of gill, thymus and brain, the expression level of GHR in 2-year common carp was significantly tower than that of 4-month common carp.

Genomic organization, gene duplication, and expression analysis of interleukin-1 beta in channel catfish (Ictalurus punctatus)

Interleukin-1 beta (IL-1 beta) is one of the pivotal early response pro-inflammatory cytokines that enables organisms to respond to infection and induces a cascade of reactions leading to inflammation. In spite of its importance and two decades of studies in the mammalian species, genes encoding IL-1 beta were not identified from non-mammalian species until recently. Recent research, particularly with genomic approaches, has led to sequencing of IL-1 beta from many species. Clinical studies also Suggested IL-1 beta as an immunoreagulatory molecule potentially useful for enhancing vaccination. However, no IL-1 beta genes have been identified from channel catfish, the primary aquaculture species from the United States. In this study, we identified two distinct cDNAs encoding catfish IL-1 beta. Their encoding genes were identified, sequenced, and characterized. The catfish IL-1 beta genes were assigned to bacterial artificial chromosome (BAC) clones. Genomic studies indicated that the IL-1 beta genes were tandemly duplicated on the same chromosome. Phylogenetic analysis of various IL-1 beta genes indicated the possibility of recent species-specific gene duplications in channel catfish, and perhaps also in swine and carp. Expression analysis indicated that both IL-1 beta genes were expressed, but exhibited distinct expression profiles in various catfish tissues, and after bacterial infection with Edwardsiella ictaluri. (c) 2005 Elsevier Ltd. All rights reserved.

Identification of one intron loss and phylogenetic evolution of Dfak gene in the Drosophila melanogaster species group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5-10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.

Multiplex detection of mutations in clinical isolates of rifampin-resistant Mycobacterium tuberculosis by short oligonucleotide Ligation assay on DNA chips

A new approach, short-oligonucleotide-ligation assay on DNA chip (SOLAC), is developed to detect mutations in rifampin-resistant Mycobacterium tuberculosis. The method needs only four common probes to detect 15 mutational variants of the rpoB gene within 12 h. Fifty-five rifampin-resistant M. tuberculosis isolates were analyzed, resulting in 87.3% accuracy and 83.6% concordance relative to DNA sequencing.

Genetic structure of Chinese sucker population Myxocyprinus asiaticus in the Yangtze River based on mitochondrial DNA marker

The sequencing analysis of the mitochondrial DNA control region (mtCR DNA) was performed to assess the genetic divergence and population structure of the Chinese sucker Myxocyprinus asiaticus (Cypriniformes Catostomidae) using four sample lots from natural populations of the Yangtze River. The mtCR DNA sequences of approximately 920 base pairs were obtained. A total of 223 nucleotide positions were polymorphic, and these defined 39 haplotypes. Of the 39 haplotypes, 37 (90%) were not shared, and among the populations as a whole there was little sharing of haplotypes. The average haplotype diversity (0.958) and the average nucleotide diversity (0.052) indicated a higher level of genetic diversity of Chinese sucker through the river. Analysis of molecular variation (AMOVA) of data revealed significant partitioning of variance (P<0.001) among populations (60.29%), and within populations (39.71%). The topology according to the neighbor joining and maximum parsimony methods showed mosaic composition of the 39 haplotypes, suggesting that the populations wore not completely divergent. The pairwise F statistic values, however, indicated that the population structuring existed to some extent among the geographic populations. There was a positive relationship between the aquatic distance and the genetic distance (Fst) among the populations (P<0.05). Based on our data, it is suggested that genetic drift, gene flow, and stochastic events are the possible factors influencing the population structure and genetic variation.