Utilization of Probabilistic Models in Short Read Assembly from Second-Generation Sequencing

With the advent of cheaper and faster DNA sequencing technologies, assembly methods have greatly changed. Instead of outputting reads that are thousands of base pairs long, new sequencers parallelize the task by producing read lengths between 35 and 400 base pairs. Reconstructing an organism’s genome from these millions of reads is a computationally expensive task. Our algorithm solves this problem by organizing and indexing the reads using n-grams, which are short, fixed-length DNA sequences of length n. These n-grams are used to efficiently locate putative read joins, thereby eliminating the need to perform an exhaustive search over all possible read pairs. Our goal was develop a novel n-gram method for the assembly of genomes from next-generation sequencers. Specifically, a probabilistic, iterative approach was utilized to determine the most likely reads to join through development of a new metric that models the probability of any two arbitrary reads being joined together. Tests were run using simulated short read data based on randomly created genomes ranging in lengths from 10,000 to 100,000 nucleotides with 16 to 20x coverage. We were able to successfully re-assemble entire genomes up to 100,000 nucleotides in length.

Application of Discrete Fourier Inter-Coefficient Difference for Assessing Genetic Sequence Similarity

Digital signal processing (DSP) techniques for biological sequence analysis continue to grow in popularity due to the inherent digital nature of these sequences. DSP methods have demonstrated early success for detection of coding regions in a gene. Recently, these methods are being used to establish DNA gene similarity. We present the inter-coefficient difference (ICD) transformation, a novel extension of the discrete Fourier transformation, which can be applied to any DNA sequence. The ICD method is a mathematical, alignment-free DNA comparison method that generates a genetic signature for any DNA sequence that is used to generate relative measures of similarity among DNA sequences. We demonstrate our method on a set of insulin genes obtained from an evolutionarily wide range of species, and on a set of avian influenza viral sequences, which represents a set of highly similar sequences. We compare phylogenetic trees generated using our technique against trees generated using traditional alignment techniques for similarity and demonstrate that the ICD method produces a highly accurate tree without requiring an alignment prior to establishing sequence similarity.

Short communication: duplication in the 5'-flanking region of the beta-lactoglobulin gene is linked to the BLG A allele

beta-Lactoglobulin (beta-LG) is the major whey protein in the milk of cows and other ruminants. It is well established that the predominant genetic variants beta-LG A and B are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants. However, the genetic basis for the differentially expressed BLG A and B alleles is still elusive. In this study additional genetic variation further upstream in the 5'-flanking region of the BLG gene was identified, including 6 single nucleotide substitutions, a single nucleotide deletion, and a 7-bp duplication. Comparison of DNA sequences showed that the investigated 5'-flanking region is highly conserved between ruminants, and the duplication g.-1885_-1879dupCTCTCGC and the substitution g.-1888A>G are only found in the BLG A and D alleles in cattle. The cytosine at position g.-1957 and the thymines at positions g.-2008 and g.-2049 are only found in BLG B alleles of cattle. It is suggested that the described genetic variability contributes to the differential allelic expression of the BLG gene.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardilpulcher species complex endemic to Lake Tanganyika

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past. (C) 2007 Elsevier Inc. All rights reserved.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Phylogeography of silver pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression, and riverine barriers

The role of Pleistocene glacial cycles in forming the contemporary genetic structure of organisms has been well studied in China with a particular focus on the Tibetan Plateau. However, China has a complex topography and diversity of local climates, and how glacial cycles may have shaped the subtropical and tropical biota of the region remains mostly unaddressed. To investigate the factors that affected the phylogeography and population history of a widely distributed and nondeciduous forest species, we analysed morphological characters, mitochondrial DNA sequences and nuclear microsatellite loci in the Silver Pheasant (Lophura nycthemera). In a pattern generally consistent with phenotypic clusters, but not nominal subspecies, deeply divergent mitochondrial lineages restricted to different geographic regions were detected. Coalescent simulations indicated that the time of main divergence events corresponded to major glacial periods in the Pleistocene and gene flow was only partially lowered by drainage barriers between some populations. Intraspecific cytonuclear discordance was revealed in mitochondrial lineages from Hainan Island and the Sichuan Basin with evidence of nuclear gene flow from neighbouring populations into the latter. Unexpectedly, hybridization was revealed in Yingjiang between the Silver Pheasant and Kalij Pheasant (Lophura leucomelanos) with wide genetic introgression at both the mtDNA and nuclear levels. Our results highlight a novel phylogeographic pattern in a subtropical area generated from the combined effects of climate oscillation, partial drainage barriers and interspecific hybridization. Cytonuclear discordance combined with morphological differentiation implies that complex historical factors shaped the divergence process in this biodiversity hot spot area of southern China.

Nano-mechanical transduction of polymer micro-cantilevers to detect bio-molecular interactions

Using variothermal polymer micro-injection molding, disposable arrays of eight polymer micro-cantilevers each 500 μm long, 100 μm wide and 25 μm thick were fabricated. The present study took advantage of an easy flow grade polypropylene. After gold coating for optical read-out and asymmetrical sensitization, the arrays were introduced into the Cantisens(®) Research system to perform mechanical and functional testing. We demonstrate that polypropylene cantilevers can be used as biosensors for medical purposes in the same manner as the established silicon ones to detect single-stranded DNA sequences and metal ions in real-time. A differential signal of 7 nm was detected for the hybridization of 1 μM complementary DNA sequences. For 100 nM copper ions the differential signal was found to be (36 ± 5) nm. Nano-mechanical sensing of medically relevant, nanometer-size species is essential for fast and efficient diagnosis.

Mutagenicity of herpes simplex virus type 1 in a shuttle vector

The shuttle vector plasmid pZ189 was used to find the kinds of mutations that are induced by herpes simplex virus type-1 (HSV-1). In cells infected by HSV-1 the frequency of mutation in supF gene, the mutagenesis marker, was increased over background by from two- to seven-fold, reaching 0.14-0.45%. No increase was induced by infection by vaccinia virus under the same conditions. Mutagenesis was an early event, showing a four-fold increase in mutation frequency at only two hours after infection, and peaking at a seven-fold increase at four hours after infection. DNA sequencing and gel electrophoresis analysis were performed on 105 HSV-1 induced mutants and 65 spontaneous mutants and provided the following information: (1) A change in plasmid size was seen in 54% of HSV-1 related mutants, compared with only 37% of spontaneous mutants. (2) Among point mutations, the predominant type was G:C to A:T transition, which accounted for 51% of point mutations in mutants isolated from cells infected with HSV-1, and 32% of point mutations in spontaneous mutants. (3) Deletions of DNA were seen in HSV-1 related mutants at a frequency of 40%, compared with 29% in spontaneous mutants. The HSV-1 related deletions were about half the length of spontaneous mutants and three contained short filler sequences. (4) Fifteen (15%) of HSV-1 induced mutants revealed the altered restriction patterns on agarose gel electrophoresis analysis and were due either to rearrangements of plasmid DNA, and/or to insertion of sequences derived from chromosomal DNA (seven plasmids). No insertions of DNA from HSV-1 were detected. Among spontaneous mutants, only 5 (7.7%) were rearrangements and none had inserted chromosomal DNA. (5) DNA sequence analysis of seven plasmids with inserted chromosomal DNA revealed that four cases had repetitive DNA sequences integrated and the other three were unidentified sequences from the GenBank database. Three repetitive DNA included $\alpha$ satellite, Alu and KpnI family sequences. The other sequence was identified as tRNA-like component. The observed mutations have implications for the mechanism of malignant transformation of cells by HSV-1. ^

Analysis of a human placental lactogen gene promoter

Human placental lactogen (hPL) and human growth hormone (hGH) comprise a multigene family that share $>$90% nucleic acid sequence homology including 500 bp of 5$\sp\prime$ flanking sequence. Despite these similarities, hGH is produced in the anterior pituitary while hPL is expressed in the placenta. For most genes studied to date, regulation of expression occurs by alterations at the level of transcriptional initiation. Nuclear proteins bind specific DNA sequences in the promoter to regulate gene expression. In this study, the hPL$\sb3$ promoter was analyzed for DNA sequences that contribute to its expression. The interaction between the hPL$\sb3$ promoter and nuclear proteins was examined using nuclear extracts from placental and non-placental cells.^ To identify regulatory elements in the promoter of the hPL$\sb3$ gene, 5$\sp\prime$ deletion mutants were constructed by cleaving 1200 bp of upstream sequence with various restriction enzymes. These DNA fragments were ligated 5$\sp\prime$ to a promoterless bacterial gene chloramphenicol acetyltransferase (CAT) and transfected into JEG-3 cells, a human placental choriocarcinoma cell line. The level of CAT activity reflects the ability of the promoter mutants to activate transcription. Deletion of the sequence between $-$142 bp and $-$129 bp, relative to the start of transcription, resulted in an 8-fold decrease in CAT activity. Nuclear proteins from JEG-3, HeLa, and HepG2 (human liver cells), formed specific binding complexes with this region of the hPL$\sb3$ promoter, as shown by gel mobility shift assay. The $-$142 bp to $-$129 bp region contains a sequence similar to that of a variant binding site for the transcription factor Sp1. Sp1-like proteins were identified by DNA binding assay, in the nuclear extracts of the three cell lines. A series of G nucleotides in the hPL$\sb3$ promoter regulatory region were identified by methylation interference assay to interact with the DNA-binding proteins and the pattern obtained is similar to that for other Sp1 binding sites that have been studied. This suggests that hPL$\sb3$ may be transcriptionally regulated by Sp1 or a Sp1-like transacting factor. ^

Molecular characterization of a gene required for entry into S phase of the mitotic cell cycle in the yeast {\it Saccharomyces cerevisiae\/}

A strain of Saccaromyces cerevisiae (SC3B) with a temperature sensitive defect in the synthesis of DNA has been isolated. This defect is due to a single recessive mutation in a gene named INS1 required for the initiation of S phase. Arrested cells carrying the ins1$\sp{ts}$ allele are defective in the completion of G1 to S phase transition events including SPB duplication or separation, initiation of DNA synthesis, normal control of budding, and bud neck stability. The mutation and a gene which complements the mutation were mapped to chromosome IV. The complementing gene was proved to be the wild type allele of the temperature sensitive mutation by genetic linkage of an integrated clone. A very low abundance 4.2 kb RNA message was observed in the strain SC3B which increased greatly in this strain transformed with a multiple copy plasmid carrying the complementing clone. The wild type gene was sequenced and found to encode a 1268 amino acid protein of with a molecular weight of 142,655 Daltons. Computer assisted searches for similar DNA sequences revealed no significant homology matches. However, searches for protein sequence homology revealed a protein (the DIS3 gene product of S. pombe) with a similar sequence over a 534 amino acid stretch to the predicted INS1 gene product. A later search revealed a near identical sequence for a gene (SRK1) also isolated from S. cerevisiae. ^

The role of the amino terminus of Bacillus subtilis sigmaK in transcription initiation

The sigma (σ) subunit of eubacterial RNA polymerase (RNAP) is required for specific recognition of promoter DNA sequences and transcription initiation. Regulation of bacterial gene expression can be achieved by modulating a factor activity. The Bacillus subtilis sporulation a σ factor, σ K, controls gene expression of the late sporulation regulon. σ K is synthesized as an inactive precursor protein, pro-σ K, with a 20 amino acid pro sequence. Proteolytic processing of the pro sequence produces the active form, σK, which is able to bind to the core subunits of RNAP to direct gene expression. Thus, the pro sequence renders σK inactive in vivo. After processing, the amino terminus of σK consists of region 1.2, which is conserved among various σ factors. To understand the role of the amino terminus of σK, namely the pro sequence and region 1.2, mutagenesis of both regions was pursued. NH 2-terminal truncations of pro-σK were constructed to address how the pro sequence silences σK activity. The work described here shows that the pro sequence inhibits the ability of σ K to associate with the core subunits and that a deletion of only six amino acids of the pro sequence is sufficient to activate pro-σ K for DNA binding and transcription initiation to levels similar to σ K. Additionally, site directed mutagenesis was used to obtain single amino acid substitutions in region 1.2 to address the role of region 1.2 in σ K transcriptional activity. Two mutations were isolated, converting a lysine (K) to an alanine (A) at position three, and an asparagine (N) to a tyrosine (Y) at position five, both of which alter the efficiency of transcription initiation by RNAP containing the mutant σKs. Surprisingly, σ KK3A increased transcript production when compared to wild type. This increase is due to improvement in DNA affinity and increased stability of RNAP-DNA promoter open complexes. σKN5Y showed a decrease in transcription activity that is related to defects in the ability of RNAP to make the transition from the closed to open RNAP-DNA complex. Results of both the pro sequence and region 1.2 analyses indicate that the amino terminus of σK is important for transcription activity and this work adds to the increasing body of evidence that the amino termini of many σ factors modulate transcription initiation by RNAP. ^

PHENOTYPIC ASSOCIATION OF GENE AMPLIFICATION WITH MULTIPLE DRUG RESISTANCE IN VINCRISTINE RESISTANT CHINESE HAMSTER OVARY CELLS

Resistance of tumors to pharmacologic agents poses a significant problem in the treatment of human malignancies. This study overviews the scope of clinical resistance and focuses upon current research attempts toward investigation of the phenomenon of multidrug resistance (MDR).^ The objective of this investigation was to determine whether gene amplification had a role in the development of the MDR phenotype in Chinese hamster ovary cells (CHO) primarily selected for resistance to vincristine (VCR). A DNA fragment, previously shown to be amplified in two independently derived Chinese hamster cell lines exhibiting the MDR phenotype, was also amplified in VCR hamster lines. Sequences flanking this fragment were shown to contain coding information for a 4.3 kb transcript overproduced in VCR cells. These sequences were not enriched in double minute DNA preparations isolated from VCR cells. There was an approximately forty-fold increase in both the level of gene amplification and transcript overproduction in the VCR cell lines, independent of the level of primary resistance. This DNA amplification and overproduction of the 4.3 kb transcript was also demonstrated in CHO cells independently selected for resistance to Adriamycin and vinblastine.^ All the DNA sequences of two hamster cDNA clones containing 785 and 932 base pair inserts showed direct homology to the published mouse mdr sequences (about 90%). This sequence conservation held for only portions of the gene when the human mdr1 sequences were compared with those from either the mouse or hamster.^ Somatic cell hybrids, constructed between VCR CHO cells and sensitive murine cells, were used to determine whether there was a functional relationship between the chromosome bearing the amplified sequences and the MDR phenotype. Concordant segregation between vincristine resistance, the MDR phenotype, the presence of MDR-associated amplified sequences, overexpression of the mRNA encoded by these sequences, overexpression of the mRNA encoded by these sequences, and CHO chromosome Z1 was consistent with the hypothesis that there is an amplified gene on chromosome Z1 of the VCR CHO cells which is responsible for MDR in these cells. ^

Mycobacterium ulcerans persistence at a village water source of Buruli ulcer patients

Buruli ulcer (BU), a neglected tropical disease of the skin and subcutaneous tissue, is caused by Mycobacterium ulcerans and is the third most common mycobacterial disease after tuberculosis and leprosy. While there is a strong association of the occurrence of the disease with stagnant or slow flowing water bodies, the exact mode of transmission of BU is not clear. M. ulcerans has emerged from the environmental fish pathogen M. marinum by acquisition of a virulence plasmid encoding the enzymes required for the production of the cytotoxic macrolide toxin mycolactone, which is a key factor in the pathogenesis of BU. Comparative genomic studies have further shown extensive pseudogene formation and downsizing of the M. ulcerans genome, indicative for an adaptation to a more stable ecological niche. This has raised the question whether this pathogen is still present in water-associated environmental reservoirs. Here we show persistence of M. ulcerans specific DNA sequences over a period of more than two years at a water contact location of BU patients in an endemic village of Cameroon. At defined positions in a shallow water hole used by the villagers for washing and bathing, detritus remained consistently positive for M. ulcerans DNA. The observed mean real-time PCR Ct difference of 1.45 between the insertion sequences IS2606 and IS2404 indicated that lineage 3 M. ulcerans, which cause human disease, persisted in this environment after successful treatment of all local patients. Underwater decaying organic matter may therefore represent a reservoir of M. ulcerans for direct infection of skin lesions or vector-associated transmission.

Hierarchical self-assembly of nucleotide-appended oligopyrenotides into defined supramolecular objects

Supramolecular DNA assembly blends DNA building blocks with synthetic organic and inorganic molecules giving structural and functional advantages both to the initial self-assembly process and to the final construct. Synthetic molecules can bring a number of additional interactions into DNA nanotechnology. Incorporating extended aromatic molecules as connectors of DNA strands allows folding of these strands through π-π stacking (DNA “foldamers”). In previous work it was shown that short oligopyrenotides (phosphodiester-linked pyrene oligomers) behave as staircase-like foldamers, which cooperatively self-assemble into two-dimensional supramolecular polymers in aqueous medium. Herein, we demonstrate that a 10-mer DNA-sequence modified with 7 pyrene units (see illustration) forms dimensionally-defined supramolecular polymers under thermodynamic conditions in water. We present the self-assembly behavior, morphological studies, and the spectroscopic properties of the investigated DNA-sequences (illustrative AFM picture shown below).

Identification of Clostridium chauvoei in cultures and clinical material from blackleg using PCR.

An identification system for Clostridium chauvoei, using PCR amplification of the 16S rRNA gene (rrs) with specific oligonucleotide primers and subsequent restriction digestion of the amplification product is described. The specific oligonucleotide primers were designed based on the rrs gene sequences of C. chauvoei by comparing it to the DNA sequences of the rrs genes of its most closely related species Clostridium septicum and Clostridium carnis. A subsequent restriction digestion of the 960 bp amplification product was used in order to unambiguously identify C. chauvoei. The developed identification system was evaluated on clinical material during a recent outbreak of blackleg in cattle. Thereby, C. chauvoei was identified as the etiologic agent of the outbreak either directly from clinical samples of muscle, liver, spleen and kidney or from primary cultures made with this material. A comparison of the newly developed method with standard diagnostic tools for C. chauvoei showed that it has advantages over the immunofluorescence and is, therefore, a useful option to it. Moreover, the assay is a valuable tool for the phylogenetic identification of C. chauvoei which can assist to substitute the fastidious traditional identification methods and replace laboratory animal testing currently used.