SeqDoC: rapid SNP and mutation detection by direct comparison of DNA sequence chromatograms

Background: This paper describes SeqDoC, a simple, web-based tool to carry out direct comparison of ABI sequence chromatograms. This allows the rapid identification of single nucleotide polymorphisms (SNPs) and point mutations without the need to install or learn more complicated analysis software. Results: SeqDoC produces a subtracted trace showing differences between a reference and test chromatogram, and is optimised to emphasise those characteristic of single base changes. It automatically aligns sequences, and produces straightforward graphical output. The use of direct comparison of the sequence chromatograms means that artefacts introduced by automatic base-calling software are avoided. Homozygous and heterozygous substitutions and insertion/deletion events are all readily identified. SeqDoC successfully highlights nucleotide changes missed by the Staden package 'tracediff' program. Conclusion: SeqDoC is ideal for small-scale SNP identification, for identification of changes in random mutagenesis screens, and for verification of PCR amplification fidelity. Differences are highlighted, not interpreted, allowing the investigator to make the ultimate decision on the nature of the change.

Discovery of genes associated with fruit ripening in Carica papaya using expressed sequence tags

To identify genes involved in papaya fruit ripening, a total of 1171 expressed sequence tags (ESTs) were generated from randomly selected clones of two independent fruit cDNA libraries derived from yellow and red-fleshed fruit varieties. The most abundant sequences encoded: chitinase, 1-aminocyclopropane- 1-carboxylic acid (ACC) oxidase, catalase and methionine synthase, respectively. DNA sequence comparisons identified ESTs with significant similarity to genes associated with fruit softening, aroma and colour biosynthesis. Putative cell wall hydrolases, cell membrane hydrolases, and ethylene synthesis and regulation sequences were identified with predicted roles in fruit softening. Expressed papaya genes associated with fruit aroma included isoprenoid biosynthesis and shikimic acid pathway genes and proteins associated with acyl lipid catabolism. Putative fruit colour genes were identified due to their similarity with carotenoid and chlorophyll biosynthesis genes from other plant species. © 2005 Elsevier Ireland Ltd. All rights reserved.

Identification of the dominant translation start site in the attB1 sequence of the pET-DEST42 Gateway vector

Gateway technology is a powerful system for converting a single entry vector into a wide variety of expression vectors. We expressed recombinant influenza matrix protein M1 (FMP), a potent antigen for cytotoxic T cells, using the Gateway vector pET-DEST42 containing the FMP cDNA, and purified the expressed FMP as a single 32 kDa recombinant protein. N-terminal and internal protein sequencing, however, showed that the recombinant FMP contained an extra 10 amino acids fused to the N-terminal of native FMP. Further investigation of the DNA sequence adjacent to the 5'-FMP cDNA indicated that the TTG in the attB1 site (30bp upstream of the ATG in the 5'-FMP cDNA) behaved as a dominant translation start site, resulting in a 10 amino acid extension of the recombinant FMP. Thus, it is possible that recombinant proteins produced by this Gateway vector contain unexpected vector-derived peptides, which may affect experimental outcomes. (c) 2006 Elsevier Inc. All rights reserved.

Polymorphisms in the 5 '-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 50 segment, spanning common DNA sequence variations, T-261G, A-161T, and -182INS/DEL-181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype -261G_-182INS-181_A-161 enhanced transcriptional activity 2.3-fold compared with the haplotype T-261_-182INS-181_A-161. Conversely, -161T reversed this, and the net effect when -261G and -161T were in the same haplotype (-261G_-182INS-181_-161T) was equivalent to the major haplotype (T-261_-182INS-181_A-161). Electrophoretic mobility shift experiments showed that -261G and -161T modify the binding of transcription factors (TFs): -261G generates a new AP2 binding site, while alleles A-161 and -161T exhibit different binding characteristics to AP1. T-261G and A-161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.

Highly organized structure in the non-coding region of the psbA minicircle from clade C Symbiodinium

The chloroplast genes of dinoflagellates are distributed among small, circular dsDNA molecules termed minicircles. In this paper, we describe the structure of the non-coding region of the psbA minicircle from Symbiodinium. DNA sequence was obtained from five Symbiodinium strains obtained from four different coral host species (Goniopora tenuidens, Heliofungia actiniformis, Leptastrea purpurea and Pocillopora damicornis), which had previously been determined to be closely related using LSU rDNA region D1/D2 sequence analysis. Eight distinct sequence blocks, consisting of four conserved cores interspersed with two metastable regions and flanked by two variable regions, occurred at similar positions in all strains. Inverted repeats (IRs) occurred in tandem or 'twin' formation within two of the four cores. The metastable regions also consisted of twin IRs and had modular behaviour, being either fully present or completely absent in the different strains. These twin IRs are similar in sequence to double-hairpin elements (DHEs) found in the mitochondrial genomes of some fungi, and may be mobile elements or may serve a functional role in recombination or replication. Within the central unit (consisting of the cores plus the metastable regions), all IRs contained perfect sequence inverses, implying they are highly evolved. IRs were also present outside the central unit but these were imperfect and possessed by individual strains only. A central adenine-rich sequence most closely resembled one in the centre of the non-coding part of Amphidinium operculatum minicircles, and is a potential origin of replication. Sequence polymorphism was extremely high in the variable regions, suggesting that these regions may be useful for distinguishing strains that cannot be differentiated using molecular markers currently available for Symbiodinium.

Targeting a complex transcriptome: The construction of the mouse full-length cDNA encyclopedia

We report the construction of the mouse full-length cDNA encyclopedia, the most extensive view of a complex transcriptome, on the basis of preparing and sequencing 246 libraries. Before cloning, cDNAs were enriched in full-length by Cap-Trapper, and in most cases, aggressively subtracted/normalized. We have produced 1,442,236 successful 3'-end sequences clustered into 171,144 groups, from which 60,770 clones were fully sequenced cDNAs annotated in the FANTOM-2 annotation. We have also produced 547,149 5' end reads, which clustered into 124,258 groups. Altogether, these cDNAs were further grouped in 70,000 transcriptional units (TU), which represent the best coverage of a transcriptome so far. By monitoring the extent of normalization/subtraction, we define the tentative equivalent coverage (TEC), which was estimated to be equivalent to >12,000,000 ESTs derived from standard libraries. High coverage explains discrepancies between the very large. numbers of clusters (and TUs) of this project, which also include non-protein-coding RNAs, and the lower gene number estimation of genome annotations. Altogether, S'-end clusters identify regions that are potential promoters for 8637 known genes and S'-end clusters suggest the presence of almost 63,000 transcriptional starting points. An estimate of the frequency of polyadenylation signals suggests that at least half of the singletons in the EST set represent real mRNAs. Clones accounting for about half of the predicted TUs await further sequencing. The continued high-discovery rate suggests that the task of transcriptome discovery is not yet complete.

Development and evaluation of an automated annotation pipeline and cDNA annotation system

Manual curation has long been held to be the gold standard for functional annotation of DNA sequence. Our experience with the annotation of more than 20,000 full-length cDNA sequences revealed problems with this approach, including inaccurate and inconsistent assignment of gene names, as well as many good assignments that were difficult to reproduce using only computational methods. For the FANTOM2 annotation of more than 60,000 cDNA clones, we developed a number of methods and tools to circumvent some of these problems, including an automated annotation pipeline that provides high-quality preliminary annotation for each sequence by introducing an uninformative filter that eliminates uninformative annotations, controlled vocabularies to accurately reflect both the functional assignments and the evidence supporting them, and a highly refined, Web-based manual annotation tool that allows users to view a wide array of sequence analyses and to assign gene names and putative functions using a consistent nomenclature. The ultimate utility of our approach is reflected in the low rate of reassignment of automated assignments by manual curation. Based on these results, we propose a new standard for large-scale annotation, in which the initial automated annotations are manually investigated and then computational methods are iteratively modified and improved based on the results of manual curation.

Identification and cloning of the gene encoding BmpC: an outer-membrane lipoprotein associated with Brachyspira pilosicoli membrane vesicles

The intestinal spirochaete Brachyspira pilosicoli causes colitis in a wide variety of host species. Little is known about the structure or protein constituents of the B. pilosicoli outer membrane (OM). To identify surface-exposed proteins in this species, membrane vesicles were isolated from B. pilosicoli strain 95-1000 cells by osmotic lysis in dH(2)O followed by isopycnic centrifugation in sucrose density gradients. The membrane vesicles were separated into a high-density fraction (HDMV; p = 1.18 g CM-3) and a low-density fraction (LDMV; rho=1.12 g cm(-3)). Both fractions were free of flagella and soluble protein contamination. LDMV contained predominantly OM markers (lipo-oligosaccharide and a 29 kDa B. pilosicoli OM protein) and was used as a source of antigens to produce mAbs. Five B. pilosicoli-specific mAbs reacting with proteins with molecular masses of 23, 24, 35, 61 and 79 kDa were characterized. The 23 kDa protein was only partially soluble in Triton X-114, whereas the 24 and 35 kDa proteins were enriched in the detergent phase, implying that they were integral membrane proteins or lipoproteins. All three proteins were localized to the B. pilosicoli OM by immunogold labelling using specific mAbs. The gene encoding the abundant, surface-exposed 23 kDa protein was identified by screening a B. pilosicoli 95-1000 genome library with the mAb and was expressed in Escherichia coli. Sequence analysis showed that it encoded a unique lipoprotein, designated BmpC. Recombinant BmpC partitioned predominantly in the OM fraction of E. coli strain SOLR. The mAb to BmpC was used to screen a collection of 13 genetically heterogeneous strains of B. pilosicoli isolated from five different host species. Interestingly, only strain 95-1000 was reactive with the mAb, indicating that either the surface-exposed epitope on BmpC is variable between strains or that the protein is restricted in its distribution within B. pilosicoli.

Characterization of Streptococcus bovis from the rumen of the dromedary camel and Rusa deer

Aims: Isolation and characterization of Streptococcus bovis from the dromedary camel and Rusa deer. Methods and Results: Bacteria were isolated from the rumen contents of four camels and two deer fed lucerne hay by culturing on the semi-selective medium MRS agar. Based on Gram morphology and RFLP analysis seven isolates, MPR1, MPR2, MPR3, MPR4, MPR5, RD09 and RD11 were selected and putatively identified as Streptococcus. The identity of these isolates was later confirmed by comparative DNA sequence analysis of the 16S rRNA gene with the homologous sequence from S. bovis strains, JB1, C14b1, NCFB2476, SbR1, SbR7 and Sb5, from cattle and sheep, and the Streptococcus equinus strain NCD01037T. The percentage similarity amongst all strains was >99%, confirming the identification of the camel isolates as S. bovis. The strains were further characterized by their ability to utilize a range of carbohydrates, the production of volatile fatty acids (VFA) and lactate and the determination of the doubling time in basal medium 10 supplemented with glucose. All the isolates produced L-lactate as a major fermentation end product, while four of five camel isolates produced VFA. The range of carbohydrates utilized by all the strains tested, including those from cattle and sheep were identical, except that all camel isolates and the deer isolate RD11 were additionally able to utilize arabinose. Conclusions: Streptococcus bovis was successfully isolated from the rumen of camels and deer, and shown by molecular and biochemical characterization to be almost identical to S. bovis isolates from cattle and sheep. Significance and Impact of the Study: Streptococcus bovis is considered a key lactic acid producing bacterium from the gastrointestinal tract of ruminants, and has been implicated as a causative agent of lactic acidosis. This study is the first report of the isolation and characterization of S. bovis from the dromedary camel and Rusa deer, and suggests a major contributive role of this bacterium to fermentative acidosis.

The genetic diversity of lactic acid producing bacteria in the equine gastrointestinal tract

Seventy-two lactic acid producing bacterial isolates (excluding streptococci) were cultured from the gastrointestinal tract of six horses. Two of the horses were orally dosed with raftilose to induce lactic acidosis and laminitis while the remaining four were maintained on a roughage diet. Near complete 16S rDNA was amplified by PCR from the genomic DNA of each isolate. Following RFLP analysis with the restriction enzymes MboI, HhaI and HinfI, the PCR products from the IS isolates that produced L- and/or D-lactate were subsequently cloned and sequenced. DNA sequence analysis indicated that the majority of the isolates were closely related to species within the genus Lactobacillus, including Lactobacillus salivarius, Lactobacillus mucosae and Lactobacillus delbrueckii. Four isolates were closely related to Mitsuokella jalaludinii. Lactic acid producing bacteria (LAB) from the equine gastrointestinal tract was dominated by representatives from the genus Lactobacillus, but also included D-lactate-producing bacteria closely related to M. jalaludinii. Identification and characterization of LAB from the equine gastrointestinal tract should contribute to our understanding and management of fermentative acidosis, ulceration of the stomach and laminitis. (c) 2005 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Inherent size constraints on prokaryote gene networks due to "accelerating" growth

Networks exhibiting accelerating growth have total link numbers growing faster than linearly with network size and either reach a limit or exhibit graduated transitions from nonstationary-to-stationary statistics and from random to scale-free to regular statistics as the network size grows. However, if for any reason the network cannot tolerate such gross structural changes then accelerating networks are constrained to have sizes below some critical value. This is of interest as the regulatory gene networks of single-celled prokaryotes are characterized by an accelerating quadratic growth and are size constrained to be less than about 10,000 genes encoded in DNA sequence of less than about 10 megabases. This paper presents a probabilistic accelerating network model for prokaryotic gene regulation which closely matches observed statistics by employing two classes of network nodes (regulatory and non-regulatory) and directed links whose inbound heads are exponentially distributed over all nodes and whose outbound tails are preferentially attached to regulatory nodes and described by a scale-free distribution. This model explains the observed quadratic growth in regulator number with gene number and predicts an upper prokaryote size limit closely approximating the observed value. (c) 2005 Elsevier GmbH. All rights reserved.

Molecular characterization and cancer risk associated with BRCA1 and BRCA2 splice site variants identified in multiple-case breast cancer families

Genetic screening of women from multiple-case breast cancer families and other research-based endeavors have identified an extensive collection of germline variations of BRCA1 and BRCA2 that can be classified as deleterious and have clinical relevance. For some variants, such as those in the conserved intronic splice site regions which are highly likely to alter splicing, it is not possible to classify them based on the identified DNA sequence variation alone. We studied 11 multiple-case breast cancer families carrying seven distinct splice site region genetic alterations in BRCA1 or BRCA2 (BRCA1, c.IVS6-2delA, c.IVS9-2A>C, c.IVS4-1G>T, c.IVS20+1G>A and BRCA2, c.IVS17-1G>C, c.IVS20+1G>A, c.IVS7-1G>A) and applied SpliceSiteFinder to predict possible changes in efficiency of splice donor and acceptor sites, characterized the transcripts, and estimated the average age-specific cumulative risk (penetrance) using a modified segregation analysis. SpliceSiteFinder predicted and we identified transcipts that illustrated that all variants caused exon skipping, and all but two led to frameshifts. The risks of breast cancer to age 70 yrs, averaged over all variants, over BRCA1 variants alone, and over BRCA2 variants alone, were 73% (95% confidence interval 47-93), 64% (95%CI 28-96) and 79% (95%CI 48-98) respectively (all P

hnRNP A2, a potential ssDNA/RNA molecular adapter at the telomere

The heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multi-tasking protein that acts in the cytoplasm and nucleus. We have explored the possibility that this protein is associated with telomeres and participates in their maintenance. Rat brain hnRNP A2 was shown to have two nucleic acid binding sites. In the presence of heparin one site binds single-stranded oligodeoxyribonucleotides irrespective of sequence but not the corresponding oligoribonucleotides. Both the hnRNP A2-binding cis-acting element for the cytoplasmic RNA trafficking element, A2RE, and the ssDNA telomere repeat match a consensus sequence for binding to a second sequence-specific site identified by mutational analysis. hnRNP A2 protected the telomeric repeat sequence, but not the complementary sequence, against DNase digestion: the glycine-rich domain was found to be necessary, but not sufficient, for protection. The N-terminal RRM (RNA recognition motif) and tandem RRMs of hnRNP A2 also bind the single-stranded, template-containing segment of telomerase RNA. hnRNP A2 colocalizes with telomeric chromatin in the subset of PML bodies that are a hallmark of ALT cells, reinforcing the evidence for hnRNPs having a role in telomere maintenance. Our results support a model in which hnRNP A2 acts as a molecular adapter between single-stranded telomeric repeats, or telomerase RNA, and another segment of ssDNA.

Epigenetics of lung cancer

Epigenetics is the study of heritable changes in gene expression that occur without changes in DNA sequence. It has a role in determining when and where a gene is expressed during development. Perhaps the most well known epigenetic mechanism is DNA methylation whereby cytosines at position 5 in CpG dinucleotides are methylated. Histone modification is another form of epigenetic control, which is quite complex and diverse. Histones and DNA make up the nucleosome which is the structural unit of chromatin which are involved in packaging DNA. Apart from the crucial role epigenetics plays in embryonic development, transcription, chromatin structure, X chromosome inactivation and genomic imprinting, its role in an increasing number of human diseases is more and more recognized. These diseases include cancer, and lung cancer in particular has been increasingly studied for the potential biological role of epigenetic changes with the promise of better and novel diagnostic and therapeutic tools.

Characterization of Kudoa monodactyli n. sp (Myxosporea : Multivalvulida) from the muscle of Monodactylus argenteus (Teleostei : Monodactylidae) from Moreton Bay, Queensland, Australia

Kudoa monodactyli n. sp. is described from the somatic musculature of Monodactylus argenteus from several localities in southern Queensland, Australia. This is the first record of a myxozoan parasite from the family Monodactylidae. The spores typically have five polar capsules, making this species similar to the four other five-valved Kudoa species (K. neurophila, K. muscularis, K. shulmani, K. cutanea) that have been described to date. However, morphometric measurements particularly of spore length and width make the species from M. argenteus distinct from the other species. Comparison of the small subunit ribosomal DNA sequence of this species with its congeners for which sequence data are available, provides further evidence of novelty. Kudoa monodactyli n. sp. displays 38 (of 1,554) nucleotide differences compared with rDNA sequence of Kudoa neurophila, which on phylogenetic analysis places these species in clades exclusive of each other. Phylogenetic analyses also provide evidence that the number of valves per spore in this genus is an imperfect indicator of relatedness.