Detecting recombination in 4-taxa DNA sequence alignments with Bayesian hidden Markov models and Markov chain Monte Carlo

This article presents a statistical method for detecting recombination in DNA sequence alignments, which is based on combining two probabilistic graphical models: (1) a taxon graph (phylogenetic tree) representing the relationship between the taxa, and (2) a site graph (hidden Markov model) representing interactions between different sites in the DNA sequence alignments. We adopt a Bayesian approach and sample the parameters of the model from the posterior distribution with Markov chain Monte Carlo, using a Metropolis-Hastings and Gibbs-within-Gibbs scheme. The proposed method is tested on various synthetic and real-world DNA sequence alignments, and we compare its performance with the established detection methods RECPARS, PLATO, and TOPAL, as well as with two alternative parameter estimation schemes.

Structural characterization of a new crystal form of the four-way Holliday junction formed by the DNA sequence d(CCGGTACCGG)2: sequence versus lattice?

DNA-strand exchange is a vital step in the recombination process, of which a key intermediate is the four-way DNA Holliday junction formed transiently in most living organisms. Here, the single-crystal structure at a resolution of 2.35 Å of such a DNA junction formed by d(CCGGTACCGG)2, which has crystallized in a more highly symmetrical packing mode to that previously observed for the same sequence, is presented. In this case, the structure is isomorphous to the mismatch sequence d(CCGGGACCGG)2, which reveals the roles of both lattice and DNA sequence in determining the junction geometry. The helices cross at the larger angle of 43.0° (the previously observed angle for this sequence was 41.4°) as a right-handed X. No metal cations were observed; the crystals were grown in the presence of only group I counter-cations.

Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

C. colchicum - variation measured by DNA sequence

Cyclamen colchicum has a mixed history in the hands of botanists. This paper examines the genetic identity of a group of wild Cyclamen populations from the Caucasus to discover whther they are Cyclamen colchicum, C. purpurascens or a mixture of the two. The collections supplemented by material collected at the type locality for C. colchicum, proved to be a single but variable genetic group of C. colchicum that was distinct from C. purpurascens.

Isolation and characterization of an AINTEGUMENTA-like gene in different coconut (Cocos nucifera L.) varieties from Sri Lanka

Development of an efficient tissue culture protocol in coconut is hampered by numerous technical constraints. Thus a greater understanding of the fundamental aspects of embryogenesis is essential. The role of AINTEGUMENTA-like genes in embryogenesis has been elucidated not only in model plants but also in economically important crops. A coconut gene, CnANT, that encodes two APETALA2 (AP2) domains and a conserved linker region similar to those of the BABY BOOM transcription factor was cloned, characterized, and its tissue specific expression was examined. The full-length cDNA of 1,780 bp contains a 1,425-bp open reading frame that encodes a putative peptide of 474 amino acids. The genomic DNA sequence includes 2,317 bp and consists of nine exons interrupted by eight introns. The exon/intron organization of CnANT is similar to that of homologous genes in other plant species. Analysis of differential tissue expression by real-time polymerase chain reaction indicated that CnANT is expressed more highly in in vitro grown tissues than in other vegetative tissues. Sequence comparison of the genomic sequence of CnANT in different coconut varieties revealed one single nucleotide polymorphism and one indel in the first exon and first intron, respectively, which differentiate the Tall group of trees from Dwarfs. The indel sequence, which can be considered a simple sequence repeats marker, was successfully used to distinguish the Tall and Dwarf groups as well as to develop a marker system, which may be of value in the identification of parental varieties that are used in coconut breeding programs in Sri Lanka.

Sequence analysis and distribution of an IS3-like insertion element isolated from Salmonella enteritidis

The nucleotide sequence of a 3 kb region immediately upstream of the sef operon operon of Salmonella enteritidis was determined. A 1230 base pair insertion sequence which shared sequence identity (> 75%) with members of the IS3 family was revealed. This element, designated IS1230, had almost identical (90% identity) terminal inverted repeats to Escherichia coli IS3 but unlike other IS3-like sequences lacked the two characteristic open reading frames which encode the putative transposase. S. enteritidis possessed only one copy of this insertion sequence although Southern hybridisation analysis of restriction digests of genomic DNA revealed another fragment located in a region different from the sef operon which hybridised weakly which suggested the presence of an IS1230 homologue. The distribution of IS1230 and IS1230-like elements was shown to be widespread amongst salmonellas and the patterns of restriction fragments which hybridised differed significantly between Salmonella serotypes and it is suggested that IS1230 has potential for development as a differential diagnostic tool.

Sequence analysis and distribution in Salmonella enterica serovars of IS3-like elements

The genome of Salmonella enterica serovar Enteritidis was shown to possess three IS3-like insertion elements, designated IS1230A, B and C, and each was cloned and their respective deoxynucleotide sequences determined. Mutations in elements IS1230A and B resulted in frameshifts in the open reading frames that encoded a putative transposase to be inactive. IS1230C was truncated at nucleotide 774 relative to IS1230B and therefore did not possess the 3' terminal inverted repeat. The three IS1230 derivatives were closely related to each other based on nucleotide sequence similarity. IS1230A was located adjacent to the sef operon encoding SEF14 fimbriae located at minute 97 of the genome of S. Enteritidis. IS1230B was located adjacent to the umuDC operon at minute 42.5 on the genome, itself located near to one terminus of an 815-kb genome inversion of S. Enteritidis relative to S. Typhimurium. IS1230C was located next to attB, the bacteriophage P22 attachment site, and proB, encoding gamma-glutamyl phosphate reductase. A truncated 3' remnant of IS1230, designated IS1230T, was identified in a clinical isolate of S. Typhimurium DT193 strain 2391. This element was located next to attB adjacent to which were bacteriophage P22-like sequences. Southern hybridisation of total genomic DNA from eighteen phage types of S. Enteritidis and eighteen definitive types of S. Typhimurium showed similar, if not identical, restriction fragment profiles in the respective serovars when probed with IS1230A.

Effect of corn silage from an herbicide-tolerant genetically modified variety on milk production and absence of transgenic DNA in milk

Data from 60 multiparous Holstein cows were used in a 12-wk continuous design feeding trial. Cows were allocated to 1 of 4 experimental treatments (T1 to T4). In T1 and T2, the total mixed ration (TMR) contained either corn silage from the genetically modified (GM) variety Chardon Liberty Link, which is tolerant to the herbicide glufosinate ammonium, or its near isogenic nonGM counterpart, whereas the TMR used in T3 and T4 contained corn silage from the commercially available nonGM varieties Fabius and Antares, respectively. The objectives of the study were to determine if the inserted gene produced a marked effect on chemical composition, nutritive value, feed intake, and milk production, and to determine if transgenic DNA and the protein expressed by the inserted gene could be detected in bovine milk. The nutritive value, fermentation characteristics, mineral content, and amino acid composition of all 4 silages were similar. There were no significant treatment effects on milk yield, milk composition, and yield of milk constituents, and the dry matter (DM) intake of the GM variety was not significantly different from the 2 commercial varieties. However, although the DM intake noted for the nonGM near-isogenic variety was similar to the commercial varieties, it was significantly lower when compared with the GM variety. Polymerase chain reaction analyses of milk samples collected at wk 1, 6, and 12 of the study showed that none of the 90 milk samples tested positive, above a detection limit of 2.5 ng of total genomic DNA/mL of milk, for either tDNA (event T25) or the single-copy endogenous Zea mays gene, alcohol dehydrogenase. Using ELISA assays, the protein expressed by the T25 gene was not detected in milk.

Effect of corn silage from an herbicide-tolerant genetically modified variety on milk production and absence of transgenic DNA in milk

Data from 60 multiparous Holstein cows were used in a 12-wk continuous design feeding trial. Cows were allocated to 1 of 4 experimental treatments (T1 to T4). In T1 and T2, the total mixed ration (TMR) contained either corn silage from the genetically modified (GM) variety Chardon Liberty Link, which is tolerant to the herbicide glufosinate ammonium, or its near isogenic nonGM counterpart, whereas the TMR used in T3 and T4 contained corn silage from the commercially available nonGM varieties Fabius and Antares, respectively. The objectives of the study were to determine if the inserted gene produced a marked effect on chemical composition, nutritive value, feed intake, and milk production, and to determine if transgenic DNA and the protein expressed by the inserted gene could be detected in bovine milk. The nutritive value, fermentation characteristics, mineral content, and amino acid composition of all 4 silages were similar. There were no significant treatment effects on milk yield, milk composition, and yield of milk constituents, and the dry matter (DM) intake of the GM variety was not significantly different from the 2 commercial varieties. However, although the DM intake noted for the nonGM near-isogenic variety was similar to the commercial varieties, it was significantly lower when compared with the GM variety. Polymerase chain reaction analyses of milk samples collected at wk 1, 6, and 12 of the study showed that none of the 90 milk samples tested positive, above a detection limit of 2.5 ng of total genomic DNA/mL of milk, for either tDNA (event T25) or the single-copy endogenous Zea mays gene, alcohol dehydrogenase. Using ELISA assays, the protein expressed by the T25 gene was not detected in milk.

Phylogeny of Pelargonium (Geraniaceae) based on DNA sequences from three genomes

Phylogenetic hypotheses for the largely South African genus Pelargonium L'Hér. (Geraniaceae) were derived based on DNA sequence data from nuclear, chloroplast and mitochondrial encoded regions. The datasets were unequally represented and comprised cpDNA trnL-F sequences for 152 taxa, nrDNA ITS sequences for 55 taxa, and mtDNA nad1 b/c exons for 51 taxa. Phylogenetic hypotheses derived from the separate three datasets were overall congruent. A single hypothesis synthesising the information in the three datasets was constructed following a total evidence approach and implementing dataset specific stepmatrices in order to correct for substitution biases. Pelargonium was found to consist of five main clades, some with contrasting evolutionary patterns with respect to biogeographic distributions, dispersal capacity, pollination biology and karyological diversification. The five main clades are structured in two (subgeneric) clades that correlate with chromosome size. One of these clades includes a "winter rainfall clade" containing more than 70% of all currently described Pelargonium species, and all restricted to the South African Cape winter rainfall region. Apart from (woody) shrubs and small herbaceous rosette subshrubs, this clade comprises a large "xerophytic" clade including geophytes, stem and leaf succulents, harbouring in total almost half of the genus. This clade is considered to be the result of in situ proliferation, possibly in response to late-Miocene and Pliocene aridification events. Nested within it is a radiation comprising c. 80 species from the geophytic Pelargonium section Hoarea, all characterised by the possession of (a series of) tunicate tubers.

Clostridium bolteae sp. nov., isolated from human sources

Seven obligately anaerobic, gram-positive, rod-shaped, spore-forming organisms isolated from human sources were characterized using phenotypic and molecular taxonomic methods. Comparative 16S rRNA gene sequencing showed that the strains were genetically highly related to each other (displaying >99% sequence similarity) and represent a previously unknown sub-line within the Clostridium coccoides rRNA group of organisms. Strains of the unidentified bacterium used carbohydrate as fermentable substrates, producing acetic acid and lactic acid as the major products of glucose metabolism. The closest described species to the novel bacterium corresponded to Clostridium clostridioforme, although a 16S rRNA sequence divergence of 3% demonstrated they represent different species. Genomic DNA-DNA pairing studies confirmed the separateness of the unknown species and Clostridium clostridioforme. Based on phenotypic and phylogenetic evidence, it is therefore proposed that the unknown bacterium, be classified as Clostridium bolteae sp. nov. The type strain of Clostridium bolteae is WAL 16351(T) (= ATCC(T) = BAA-613(T), CCUG(T) = 46953(T)).

Identification and characterisation of five novel miniature inverted-repeat transposable elements (MITEs) in amphioxus (Branchiostoma floridae)

As the sister group to vertebrates, amphioxus is consistently used as a model of genome evolution for understanding the invertebrate/vertebrate transition. The amphioxus genome has not undergone massive duplications like those in the vertebrates or disruptive rearrangements like in the genome of Ciona, a urochordate, making it an ideal evolutionary model. Transposable elements have been linked to many genomic evolutionary changes including increased genome size, modified gene expression, massive gene rearrangements, and possibly intron evolution. Despite their importance in genome evolution, few previous examples of transposable elements have been identified in amphioxus. We report five novel Miniature Inverted-repeat Transposable Elements (MITEs) identified by an analysis of amphioxus DNA sequence, which we have named LanceleTn-1, LanceleTn-2, LanceleTn-3a, LanceleTn-3b and LanceleTn-4. Several of the LanceleTn elements were identified in the amphioxus ParaHox cluster, and we suggest these have had important implications for the evolution of this highly conserved gene cluster. The estimated high copy numbers of these elements implies that MITEs are probably the most abundant type of mobile element in amphioxus, and are thus likely to have been of fundamental importance in shaping the evolution of the amphioxus genome.

Sequence heterogeneity of an mpb70 gene analogue in Mycobacterium kansasii

The protein antigen MPB70 is a major component of culture supernatants of Mycobacterium bovis and Is an active ingredient of bovine PPD used for skin-testing cattle for tuberculosis. we have shown that Mycobacterium kansasii possesses a similar gene that cross-reacts in a PCR test for M. bovis. Single strand conformational polymorphism analysis, and the DNA sequence of the PCR product, shows differences between M. kansasii strains, supporting the suggestion that M. kansasii is not a homogeneous species.

A maize bacterial artificial chromosome (BAC) library from the European flint inbred line F2

We report here the construction and characterisation of a BAC library from the maize flint inbred line F2, widely used in European maize breeding programs. The library contains 86,858 clones with an average insert size of approximately 90 kb, giving approximately 3.2-times genome coverage. High-efficiency BAC cloning was achieved through the use of a single size selection for the high-molecular-weight genomic DNA, and co-transformation of the ligation with yeast tRNA to optimise transformation efficiency. Characterisation of the library showed that less than 0.5% of the clones contained no inserts, while 5.52% of clones consisted of chloroplast DNA. The library was gridded onto 29 nylon filters in a double-spotted 8 × 8 array, and screened by hybridisation with a number of single-copy and gene-family probes. A 3-dimensional DNA pooling scheme was used to allow rapid PCR screening of the library based on primer pairs from simple sequence repeat (SSR) and expressed sequence tag (EST) markers. Positive clones were obtained in all hybridisation and PCR screens carried out so far. Six BAC clones, which hybridised to a portion of the cloned Rp1-D rust resistance gene, were further characterised and found to form contigs covering most of this complex resistance locus.