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.

DNA barcodes for Mexican Cactaceae, plants under pressure from wild collecting

DNA barcodes could be a useful tool for plant conservation. Of particular importance is the ability to identify unknown plant material, such as from customs seizures of illegally collected specimens. Mexican cacti are an example of a threatened group, under pressure because of wild collection for the xeriscaping trade and private collectors. Mexican cacti also provide a taxonomically and geographically coherent group with which to test DNA barcodes. Here, we sample the matK barcode for 528 species of Cactaceae including approximately 75% of Mexican species and test the utility of the matK region for species-level identification. We find that the matK DNA barcode can be used to identify uniquely 77% of species sampled, and 79-87% of species of particular conservation importance. However, this is far below the desired rate of 95% and there are significant issues for PCR amplification because of the variability of primer sites. Additionally, we test the nuclear ITS regions for the cactus subfamily Opuntioideae and for the genus Ariocarpus (subfamily Cactoideae). We observed higher rates of variation for ITS (86% unique for Opuntioideae sampled) but a much lower PCR success, encountering significant intra-individual polymorphism in Ariocarpus precluding the use of this marker in this taxon. We conclude that the matK region should provide useful information as a DNA barcode for Cactaceae if the problems with primers can be addressed, but matK alone is not sufficiently variable to achieve species-level identification. Additional complementary regions should be investigated as ITS is shown to be unsuitable

Movement of newly assimilated 13C carbon in the grass Lolium perenne and its incorporation into rhizosphere microbial DNA

One of the key processes that drives rhizosphere microbial activity is the exudation of soluble organic carbon (C) by plant roots. We describe an experiment designed to determine the impact of defoliation on the partitioning and movement of C in grass (Lolium perenne L.), soil and grass-sterile sand microcosms, using a (13)CO(2) pulse-labelling method. The pulse-derived (13)C in the shoots declined over time, but that of the roots remained stable throughout the experiment. There were peaks in the atom% (13)C of rhizosphere CO(2) in the first few hours after labelling probably due to root respiration, and again at around 100 h. The second peak was only seen in the soil microcosms and not in those with sterilised sand as the growth medium, indicating possible microbial activity. Incorporation of the (13)C label into the microbial biomass increased at 100 h when incorporation into replicating cells, as indicated by the amounts of the label in the microbial DNA, started to increase. These results indicate that the rhizosphere environment is conducive to bacterial growth and replication. The results also show that defoliation had no impact on the pattern of movement of (13)C from plant roots into the microbial population in the rhizosphere.

Characterisation of an s-type low molecular weight glutenin subunit of wheat and its proline and glutamine-rich repetitive domain

The low molecular weight glutenin subunits (LMW-GS) are major components of the glutenin polymers which determine the elastomeric properties of wheat (Triticum aestivum L.) gluten and dough. They comprise a complex mixture of components and have proved to be difficult to purify for detailed characterisation. The mature LMW subunit proteins comprise two structural domains, with one domain consisting of repeated sequences based on short peptide motifs. DNA sequences encoding this domain and a whole subunit were expressed in Escherichia coli and the recombinant proteins purified. Detailed comparisons by spectroscopy (CD, FT-IR) and dynamic light scattering indicated that the repetitive and non-repetitive domains of the proteins formed different structures with the former having an extended conformation with an equilibrium between poly-L-proline II-like structure and type II’ b-turns, and the latter a more compact globular structure rich in a-helix. Although the structures of these two domains appear to form independently, dynamic light scattering of the whole subunit dissolved in trifluoroethanol(TFE) suggested that they interact, leading to a more compact conformation. These observations may have relevance to the role of the LMW-GS in gluten structure and functionality.

Phylogeographic analysis of the chloroplast DNA variation in wild common bean (Phaseolus vulgaris L.) in the Americas

The wild common bean (Phaseolus vulgaris) is widely but discontinuously distributed from northern Mexico to northern Argentina on both sides of the Isthmus of Panama. Little is known on how the species has reached its current disjunct distribution. In this research, chloroplast DNA polymorphisms in seven non-coding regions were used to study the history of migration of wild P. vulgaris between Mesoamerica and South America. A penalized likelihood analysis was applied to previously published Leguminosae ITS data to estimate divergence times between P. vulgaris and its sister taxa from Mesoamerica, and divergence times of populations within P. vulgaris. Fourteen chloroplast haplotypes were identified by PCR-RFLP and their geographical associations were studied by means of a Nested Clade Analysis and Mantel Tests. The results suggest that the haplotypes are not randomly distributed but occupy discrete parts of the geographic range of the species. The current distribution of haplotypes may be explained by isolation by distance and by at least two migration events between Mesoamerica and South America: one from Mesoamerica to South America and another one from northern South America to Mesoamerica. Age estimates place the divergence of P. vulgaris from its sister taxa from Mesoamerica at or before 1.3 Ma, and divergence of populations from Ecuador-northern Peru at or before 0.6 Ma. As these ages are taken as minimum divergence times, the influence of past events, such as the closure of the Isthmus of Panama and the final uplift of the Andes, on the migration history and population structure of this species cannot be disregarded.

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.

The isolation of a single spore isolate of Pasteuria penetrans and its pathogenicity on Meloidogyne javanica

We have obtained a single spore isolate of Pasteuria penetrans, derived by allowing a single spore to attach to a second-stage juvenile (J2) of the root-knot nematode Meloidogyne javanica. By analysing DNA sequences at three different loci we have obtained evidence that the isolate is, indeed, genetically pure. We compared the ability of the single spore isolate and the parent population from which it was selected to attach to and parasitise both the original population of M. javanica on which it was isolated and a single egg mass line derived from it. There was no difference in the attachment of spores of the single spore isolate to juveniles compared to the parental population, although there were higher numbers of both attaching to J2 of the single egg mass line compared to its parental population. Judging from the numbers of egg masses and Pasteuria-infected females, the single spore isolate was less pathogenic to the parental population of M. javanica than was the parental spore population.

DNA barcoding of a large genus, Aspalathus L. (Fabaceae)

The DNA barcode potential of three regions (the nuclear ribosomal ITS and the plastid psbA-trnH and trnT-trnL intergenic spacers) was investigated for the plant genus Aspalathus L. (Fabaceac: Crotalarieae). Aspalathus is a large genus (278 species) that revealed low levels of DNA variation in phylogenetic studies. In a 51-species dataset for the psbA-trnH and ITS regions, 45%, and 16% of sequences respectively were identical to the sequence of at least one other species, with two species undiscriminated even when the two regions were combined. In contrast, trnT-trnL, discriminated between all species in this dataset. In a larger ITS and trnT-trnL dataset. including a further 82 species. 7 species in five pairwise comparisons remained Undiscriminated when the two regions were combined. Four of the five pairs of species not discriminated by sequence data were readily distinguished using a combination of qualitative and quantitative morphological data. The difficulty of barcoding in this group is increased by the presence of intraspecific variation in all three regions studied. In the case of psbA-trnH, three intraspecific samples had a sequence identical to at least one other species. Overall, psbA-trnH. currently a candidate for plant barcoding, was the least discriminatory region in our study.

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.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

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.

High-resolution crystal structure of the intramolecular d(TpA) thymine-adenine photoadduct and its mechanistic implications

A high-resolution crystal structure is reported for d(TpA)*, the intramolecular thymine-adenine photoadduct that is produced by direct ultraviolet excitation of the dinucleoside monophosphate d(TpA). It confirms the presence of a central 1,3-diazacyclooctatriene ring linking the remnants of the T and A bases, as previously deduced from heteronuclear NMR measurements by Zhao et al. (The structure of d(TpA)*, the major photoproduct of thymidylyl-(3'-5')-deoxyadenosine. Nucleic Acids Res., 1996, 24, 1554-1560). Within the crystal, the d(TpA)* molecules exist as zwitterions with a protonated amidine fragment of the eight-membered ring neutralizing the charge of the internucleotide phosphate monoanion. The absolute configuration at the original thymine C5 and C6 atoms is determined as 5S,6R. This is consistent with d(TpA)* arising by valence isomerization of a precursor cyclobutane photoproduct with cis-syn stereochemistry that is generated by [2 + 2] photoaddition of the thymine 5,6-double bond across the C6 and C5 positions of adenine. This mode of photoaddition should be favoured by the stacked conformation of adjacent T and A bases in B-form DNA. It is probable that the primary photoreaction is mechanistically analogous to pyrimidine dimerization despite having a much lower quantum yield.

Conformational and hydration effects of site-selective sodium, calcium and strontium ion binding to the DNA Holliday junction structure d(TCGGTACCGA)(4)

The role of metal ions in determining the solution conformation of the Holliday junction is well established, but to date the picture of metal ion binding from structural studies of the four-way DNA junction is very incomplete. Here we present two refined structures of the Holliday junction formed by the sequence d(TCGGTACCGA) in the presence of Na+ and Ca2+, and separately with Sr2+ to resolutions of 1.85 Angstrom and 1.65 Angstrom, respectively. This sequence includes the ACC core found to promote spontaneous junction formation, but its structure has not previously been reported. Almost complete hydration spheres can be defined for each metal cation. The Na+ sites, the most convincing observation of such sites in junctions to date, are one on either face of the junction crossover region, and stabilise the ordered hydration inside the junction arms. The four Ca2+ sites in the same structure are at the CG/CG steps in the minor groove. The Sr2+ ions occupy the TC/AG, GG/CC, and TA/TA sites in the minor groove, giving ten positions forming two spines of ions, spiralling through the minor grooves within each arm of the stacked-X structure. The two structures were solved in the two different C2 lattices previously observed, with the Sr2+ derivative crystallising in the more highly symmetrical form with two-fold symmetry at its centre. Both structures show an opening of the minor groove face of the junction of 8.4degrees in the Ca2+ and Na+ containing structure, and 13.4degrees in the Sr2+ containing structure. The crossover angles at the junction are 39.3degrees and 43.3degrees, respectively. In addition to this, a relative shift in the base pair stack alignment of the arms of 2.3 Angstrom is observed for the Sr2+ containing structure only. Overall these results provide an insight into the so-far elusive stabilising ion structure for the DNA Holliday junction. (C) 2003 Elsevier Science Ltd. All rights reserved.