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.

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.

Analysis of DNA polymorphism in ancient barley herbarium material: validation of the KASP SNP genotyping platform

The availability of crop specimens archived in herbaria and old seed collections represent valuable resources for the analysis of plant genetic diversity and crop domestication. The ability to extract ancient DNA (aDNA) from such samples has recently allowed molecular genetic investigations to be undertaken in ancient materials. While analyses of aDNA initially focused on the use of markers which occur in multiple copies such as the internal transcribed spacer region (ITS) within ribosomal DNA and those requiring amplification of short DNA regions of variable length such as simple sequence repeats (SSRs), emphasis is now moving towards the genotyping of single nucleotide polymorphisms (SNPs), traditionally undertaken in aDNA by Sanger sequencing. Here, using a panel of barley aDNA samples previously surveyed by Sanger sequencing for putative causative SNPs within the flowering-time gene PPD-H1, we assess the utility of the Kompetitive Allele Specific PCR (KASP) genotyping platform for aDNA analysis. We find KASP to out-perform Sanger sequencing in the genotyping of aDNA samples (78% versus 61% success, respectively), as well as being robust to contamination. The small template size (≥46 bp) and one-step, closed-tube amplification/genotyping process make this platform ideally suited to the genotypic analysis of aDNA, a process which is often hampered by template DNA degradation and sample cross-contamination. Such attributes, as well as its flexibility of use and relatively low cost, make KASP particularly relevant to the genetic analysis of aDNA samples. Furthermore, KASP provides a common platform for the genotyping and analysis of corresponding SNPs in ancient, landrace and modern plant materials. The extended haplotype analysis of PPD-H1 undertaken here (allelic variation at which is thought to be important for the spread of domestication and local adaptation) provides further resolution to the previously identified geographic cline of flowering-time allele distribution, illustrating how KASP can be used to aid genetic analyses of aDNA from plant species. We further demonstrate the utility of KASP by genotyping ten additional genetic markers diagnostic for morphological traits in barley, shedding light on the phenotypic traits, alleles and allele combinations present in these unviable ancient specimens, as well as their geographic distributions.

Land plants and DNA barcodes: short-term and long-term goals

Land plants have had the reputation of being problematic for DNA barcoding for two general reasons: (i) the standard DNA regions used in algae, animals and fungi have exceedingly low levels of variability and (ii) the typically used land plant plastid phylogenetic markers (e.g. rbcL, trnL-F, etc.) appear to have too little variation. However, no one has assessed how well current phylogenetic resources might work in the context of identification (versus phylogeny reconstruction). In this paper, we make such an assessment, particularly with two of the markers commonly sequenced in land plant phylogenetic studies, plastid rbcL and internal transcribed spacers of the large subunits of nuclear ribosomal DNA (ITS), and find that both of these DNA regions perform well even though the data currently available in GenBank/EBI were not produced to be used as barcodes and BLAST searches are not an ideal tool for this purpose. These results bode well for the use of even more variable regions of plastid DNA (such as, for example, psbA-trnH) as barcodes, once they have been widely sequenced. In the short term, efforts to bring land plant barcoding up to the standards being used now in other organisms should make swift progress. There are two categories of DNA barcode users, scientists in fields other than taxonomy and taxonomists. For the former, the use of mitochondrial and plastid DNA, the two most easily assessed genomes, is at least in the short term a useful tool that permits them to get on with their studies, which depend on knowing roughly which species or species groups they are dealing with, but these same DNA regions have important drawbacks for use in taxonomic studies (i.e. studies designed to elucidate species limits). For these purposes, DNA markers from uniparentally (usually maternally) inherited genomes can only provide half of the story required to improve taxonomic standards being used in DNA barcoding. In the long term, we will need to develop more sophisticated barcoding tools, which would be multiple, low-copy nuclear markers with sufficient genetic variability and PCR-reliability; these would permit the detection of hybrids and permit researchers to identify the 'genetic gaps' that are useful in assessing species limits.

Diversification of a large genus in a continental biodiversity hotspot: temporal and spatial origin of Muraltia (Polygalaceae) in the Cape of South Africa

Phylogenetic relationships in the largely South African genus Muraltia (Polygalaceae) are assessed based on DNA sequence data (nuclear ribosomal ITS, plastid atpB-rbcL spacer, trnL intron, and trnL-F spacer) for 73 of the 117 currently recognized species in the genus. The previously recognised subgenus Muraltia is monophyletic, but the South African endemic genus Nylandtia is embedded in Muraltia subgenus Psiloclada. Subgenus Muraltia is found to be sister to subgenus Psiloclada. Estimates show the beginning of diversification of the two subgenera in the early Miocene (Psiloclada, 19.3+/-3.4 Ma; Muraltia, 21.0+/-3.5 Ma) pre-dating the establishment of the Benguela current (intermittent in the middle to late Oligocene and markedly intensifying in the late Miocene), and summer-dry climate in the Cape region. However, the later increase in species numbers is contemporaneous with these climatic phenomena. Results of dispersal-vicariance analyses indicate that major clades in Muraltia diversified from the southwestern and northwestern Cape, where most of the species are found today.

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.

Development of a PCR-based assay for rapid and reliable identification of pathogenic Fusaria

Identification of Fusarium species has always been difficult due to confusing phenotypic classification systems. We have developed a fluorescent-based polymerase chain reaction assay that allows for rapid and reliable identification of five toxigenic and pathogenic Fusarium species. The species includes Fusarium avenaceum, F. culmorum, F. equiseti, F. oxysporum and F. sambucinum. The method is based on the PCR amplification of species-specific DNA fragments using fluorescent oligonucleotide primers, which were designed based on sequence divergence within the internal transcribed spacer region of nuclear ribosomal DNA. Besides providing an accurate, reliable, and quick diagnosis of these Fusaria, another advantage with this method is that it reduces the potential for exposure to carcinogenic chemicals as it substitutes the use of fluorescent dyes in place of ethidium, bromide. Apart from its multidisciplinary importance and usefulness, it also obviates the need for gel electrophoresis. (C) 2002 Published by Elsevier Science B.V. on behalf of the Federation of European Microbiological Societies.

Variation in 16S-23S rRNA intergenic spacer regions in Photobacterium damselae: a mosaic-like structure

Phenotypically, Photobacterium damselae subsp. piscicida and P. damselae subsp. damselae are easily distinguished. However, their 16S rRNA gene sequences are identical, and attempts to discriminate these two subspecies by molecular tools are hampered by their high level of DNA-DNA similarity. The 16S-23S rRNA internal transcribed spacers (ITS) were sequenced in two strains of Photobacterium damselae subsp. piscicida and two strains of P. damselae subsp. damselae to determine the level of molecular diversity in this DNA region. A total of 17 different ITS variants, ranging from 803 to 296 bp were found, some of which were subspecies or strain specific. The largest ITS contained four tRNA genes (tDNAs) coding for tRNA(Glu(UUC)), tRNA(LyS(UUU)), tRNA(Val(UAC)), and tRNA(Ala(GGC)). Five amplicons contained tRNA(Glu(UUC)) combined with two additional tRNA genes, including tRNA(Lys(UUU)), tRNA(Val(UAC)), or tRNA(Ala(UGC)). Five amplicons contained tRNA(Ile(GAU)) and tRNA(Ala(UGC)). Two amplicons contained tRNA(Glu(UUC)) and tRNA(Val(UGC)). Two different isoacceptor tRNA(Ala) genes (GGC and UGC anticodons) were found. The five smallest amplicons contained no tRNA genes. The tRNA-gene combinations tRNA(Glu(UUC)) -tRNA(Val(UAC)) -tRNA(Ala(UGC)) and tRNA(Glu(UUC)) -tRNA(Ala(UGC)) have not been previously reported in bacterial ITS regions. The number of copies of the ribosomal operon (rrn) in the P. damselae chromosome ranged from at least 9 to 12. For ITS variants coexisting in two strains of different subspecies or in strains of the same subspecies, nucleotide substitution percentages ranged from 0 to 2%. The main source of variation between ITS variants was due to different combinations of DNA sequence blocks, constituting a mosaic-like structure.

Molecular inventory of faecal microflora in patients with Crohn's disease

Intestinal microbial community is involved in the pathogenesis of Crohn's disease, but knowledge of its potential abnormalities has been limited by the impossibility to grow many dominant intestinal bacteria. Using sequence analysis of randomly cloned bacterial 16S ribosomal DNA, the dominant faecal species from four Crolin's disease patients and four controls were compared. Whereas marked inter-individual differences were observed in the faecal microflora of patients, three remained distantly related to controls on the basis of their operational taxonomic unit composition. Bacteroides vidgatus and closely related organisms represented the only molecular species shared by all patients and exhibited an unusually high rate of occurrence. Escherichia coli clones were isolated only in two patients with ileocolonic Crohn's disease. Moreover, numerous clones belonged to phylogenetic groups or species that are commonly not dominant in the faecal microflora of healthy subjects: Pectinatus, Sutterella, Verritcomicrobium, Fusobacterium, Clostridium disporicum, clostridium glycolicum, Clostridium ramosum, Clostridium innocuum and Clostridium perfringens. (C) 2004 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

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.

Phylogeny, function, and evolution of the cupins, a structurally conserved, functionally diverse superfamily of proteins

The cupin superfamily is a group of functionally diverse proteins that are found in all three kingdoms of life, Archaea, Eubacteria, and Eukaryota. These proteins have a characteristic signature domain comprising two histidine- containing motifs separated by an intermotif region of variable length. This domain consists of six beta strands within a conserved beta barrel structure. Most cupins, such as microbial phosphomannose isomerases (PMIs), AraC- type transcriptional regulators, and cereal oxalate oxidases (OXOs), contain only a single domain, whereas others, such as seed storage proteins and oxalate decarboxylases (OXDCs), are bi-cupins with two pairs of motifs. Although some cupins have known functions and have been characterized at the biochemical level, the majority are known only from gene cloning or sequencing projects. In this study, phylogenetic analyses were conducted on the conserved domain to investigate the evolution and structure/function relationships of cupins, with an emphasis on single- domain plant germin-like proteins (GLPs). An unrooted phylogeny of cupins from a wide spectrum of evolutionary lineages identified three main clusters, microbial PMIs, OXDCs, and plant GLPs. The sister group to the plant GLPs in the global analysis was then used to root a phylogeny of all available plant GLPs. The resulting phylogeny contained three main clades, classifying the GLPs into distinct subfamilies. It is suggested that these subfamilies correlate with functional categories, one of which contains the bifunctional barley germin that has both OXO and superoxide dismutase (SOD) activity. It is proposed that GLPs function primarily as SODs, enzymes that protect plants from the effects of oxidative stress. Closer inspection of the DNA sequence encoding the intermotif region in plant GLPs showed global conservation of thymine in the second codon position, a character associated with hydrophobic residues. Since many of these proteins are multimeric and enzymatically inactive in their monomeric state, this conservation of hydrophobicity is thought to be associated with the need to maintain the various monomer- monomer interactions. The type of structure-based predictive analysis presented in this paper is an important approach for understanding gene function and evolution in an era when genomes from a wide range of organisms are being sequenced at a rapid rate.

Inbreeding depression and founder diversity among captive and free-living populations of the endangered pink pigeon Columba mayeri

The endemic pink pigeon has recovered from less than 20 birds in the mid-1970s to 355 free-living individuals in 2003. A major concern for the species' recovery has been the potential genetic problem of inbreeding. Captive pink pigeons bred for reintroduction were managed to maximise founder representation and minimise inbreeding. In this paper, we quantify the effect of inbreeding on survival and reproductive parameters in captive and wild populations and quantify DNA sequence variation in the mitochondrial d-loop region for pink pigeon founders. Inbreeding affected egg fertility, squab, juvenile and adult survival, but effects were strongest in highly inbred birds (F≥0.25). Inbreeding depression was more apparent in free-living birds where even moderate levels of inbreeding affected survival, although highly inbred birds were equally compromised in both captive and wild populations. Mitochondrial DNA haplotypic diversity in pink pigeon founders is low, suggesting that background inbreeding is contributing to low fertility and depressed productivity in this species, as well as comparable survival of some groups of non-inbred and nominally inbred birds. Management of wild populations has boosted population growth and may be required long-term to offset the negative effects of inbreeding depression and enhance the species' survival.

Characterization and origin of infection of Rhizoctonia solani associated with Brassica oleracea crops in the UK

An extensive study was conducted to determine where in the production chain Rhizoctonia solani became associated with UK module-raised Brassica oleracea plants. In total, 2600 plants from 52 crops were sampled directly from propagators and repeat sampled from the field. Additional soil, compost and water samples were collected from propagation nurseries and screened using conventional agar isolation methods. No isolates of R. solani were recovered from any samples collected from propagation nurseries. Furthermore, nucleic acid preparations from samples of soil and compost from propagation nurseries gave negative results when tested for R. solani using real-time PCR. Conversely, R. solani was recovered from 116 of 1300 stem bases collected from field crops. All the data collected suggested R. solani became associated with B. oleracea in the field rather than during propagation. Parsimony and Bayesian phylogenetic studies of ribosomal DNA suggested the majority of further classified isolates belonged to anastomosis groups 2-1 (48/57) and AG-4HGII (8/57), groups known to be pathogenic on Brassica spp. in other countries. Many R. solani isolates were recovered from symptomless plant material and the possibilities for such an association are discussed.

The Fox genes of Branchiostoma floridae

The Fox genes are united by encoding a fork head domain, a deoxyribonucleic acid (DNA)-binding domain of the winged-helix type that marks these genes as encoding transcription factors. Vertebrate Fox genes are classified into 23 subclasses named from FoxA to FoxS. We have surveyed the genome of the amphioxus Branchiostoma floridae, identifying 32 distinct Fox genes representing 21 of these 23 subclasses. The missing subclasses, FoxR and FoxS, are specific to vertebrates, and in addition, B. floridae has one further group, FoxAB, that is not found in vertebrates. Hence, we conclude B. floridae has maintained a high level of Fox gene diversity. Expressed sequence tag and complementary DNA sequence data support the expression of 23 genes. Several linkages between B. floridae Fox genes were noted, including some that have evolved relatively recently via tandem duplication in the amphioxus lineage and others that are more ancient.