A comparative genomics strategy for targeted discovery of single-nucleotide polymorphisms and conserved-noncoding sequences in orphan crops

Completed genome sequences provide templates for the design of genome analysis tools in orphan species lacking sequence information. To demonstrate this principle, we designed 384 PCR primer pairs to conserved exonic regions flanking introns, using Sorghum/Pennisetum expressed sequence tag alignments to the Oryza genome. Conserved-intron scanning primers (CISPs) amplified single-copy loci at 37% to 80% success rates in taxa that sample much of the approximately 50-million years of Poaceae divergence. While the conserved nature of exons fostered cross-taxon amplification, the lesser evolutionary constraints on introns enhanced single-nucleotide polymorphism detection. For example, in eight rice (Oryza sativa) genotypes, polymorphism averaged 12.1 per kb in introns but only 3.6 per kb in exons. Curiously, among 124 CISPs evaluated across Oryza, Sorghum, Pennisetum, Cynodon, Eragrostis, Zea, Triticum, and Hordeum, 23 (18.5%) seemed to be subject to rigid intron size constraints that were independent of per-nucleotide DNA sequence variation. Furthermore, we identified 487 conserved-noncoding sequence motifs in 129 CISP loci. A large CISP set (6,062 primer pairs, amplifying introns from 1,676 genes) designed using an automated pipeline showed generally higher abundance in recombinogenic than in nonrecombinogenic regions of the rice genome, thus providing relatively even distribution along genetic maps. CISPs are an effective means to explore poorly characterized genomes for both DNA polymorphism and noncoding sequence conservation on a genome-wide or candidate gene basis, and also provide anchor points for comparative genomics across a diverse range of species.     

General-use polymerase chain reaction primers for amplification and direct sequencing of enolase, a single-copy nuclear gene, from different animal phyla

In contrast to mitochondrial DNA, remarkably few general-use primer sets are available for single-copy nuclear genes across animal phyla. Here, we present a primer set that yields a c. 364-bp coding fragment of the metabolic gene enolase, which includes an intron in some taxa. In species where introns are absent or have few insertions/deletions, the amplified fragment can be sequenced directly for phylogenetic or population analysis. Between-species variation in the coding region occurs widely at third codon positions, even between closely related taxa, making the fragment useful for species-level systematics. In low gene-flow species, the primers may also be of use for population genetics, as intraspecific polymorphisms occur at several silent positions in the taxa examined.     

Characterisation of type-I thionin loci from the A, B, D and R genomes of wheat and rye

 DNA sequences encoding type-I thionins were isolated from Triticum aestivum L. cv ‘Chinese Spring’ using PCR with consensus primers. Blunt-end cloning, sequencing and PCR-based chromosome assignment of these fragments uncovered the three orthologous sequences corresponding to the single-copy genes at the Pur-1 loci on each of the group-1 chromosomes. Comparison with two previously published cDNA sequences revealed the presence of two introns that contain most of the polymorphic nucleotide sites. The observed orthologous DNA sequence variation among Pur-1 loci, encoded by each of the A, B and D genomes, enabled us to establish interlocus relationships and to construct locus-specific primer sets. Analogously, the Pur-R1 sequence from rye was isolated, and a locus-specific primer pair was constructed as well. Hence, four locus-specific primer sets are now available as molecular markers for the homoeologous 1AL, 1BL, 1DL and 1RL chromosome arms. Amplification from several diploid and tetraploid wheat species showed that the primers can be used as molecular tools for studying wheat phylogeny. Received: 30 January 1997 / Accepted: 23 June 1997     

Multiple nuclear-gene phylogenies: application to pinnipeds and comparison with a mitochondrial DNA gene phylogeny

Phylogenetic analyses of closely related species should use information from multiple, independent genes with relatively high rates of sequence evolution. To investigate species for which there are few prior sequence data for single-copy nuclear (scnDNA) genes, primers for gene amplification can be designed to highly conserved regions of exons in order to amplify both coding (exons) and noncoding (introns) sequences. We have explored this approach in a phylogenetic analysis of six species of pinnipeds that, together with terrestrial carnivore outgroups, encompass divergence times < or = 40-50 Mya. We sequenced one intron from each of the aldolase A (ALD-A), aldolase C (ALD-C), and histone H2AF genes; one exon from the major-histocompatibility-complex DQA gene; a H2AF processed pseudogene (psi H2AF); and, for comparison with the nuclear genes, the 5' portion of the mitochondrial DNA (mtDNA) control region. The pinniped psi H2AF genes were found to be of limited use because they were paralogous with the gene in the outgroup. The rate of silent substitution in scnDNA (primarily introns) was 5-10-fold lower than that for mtDNA control region I, and scnDNA sequence divergence increased linearly with time < or = 40-50 Mya. Alleles at three polymorphic scnDNA loci (ALD-A, H2AF, and DQA) in the southern elephant seal were paraphyletic with respect to the allele from the closely related northern elephant seal, while the more numerous mtDNA alleles were monophyletic. This we attribute to the consequences of a higher mutation rate rather than to a lower effective population size of mtDNA compared with scnDNA. Within the short (i.e., < 500-bp) sequences of individual scnDNA sequences, phylogenetically informative variation was insufficient to obtain robust phylogenies. However, the combined scnDNA sequences produced a well-supported phylogeny congruent with that derived from mtDNA. This analysis illustrates the high resolution of mtDNA sequences compared with a similar length of scnDNA sequence, but it also demonstrates the utility of combining information from multiple short scnDNA sequences obtained using broadly applicable primers.      

Rapid development of multiple nuclear loci for phylogenetic analysis using genomic resources: an example from squamate reptiles

Recently, as genome-scale data have become available for more organisms, the development of phylogenetic markers from nuclear protein-coding loci (NPCL) has become more tractable. However, new methods are needed to efficiently sort the large number of genes from genomic databases into more limited sets appropriate for particular phylogenetic questions, while avoiding introns and paralogs. Here we describe a general methodology for identifying candidate single-copy NPCL from genomic databases. Our method uses information from reference genomes to identify genes with relatively large continuous protein-coding regions (i.e., 700bp). BLAST comparisons are used to help avoid genes with paralogous copies or close relatives (i.e., gene families) that might confound phylogenetic analyses. Exon boundary information is used to identify appropriately spaced potential priming sites. Using this method, we have developed over 25 novel NPCL, which span a variety of desirable evolutionary rates for phylogenetic analyses. Although targeted for higher-level phylogenetics of squamate reptiles, many of these loci appear to be useful across and within other vertebrate clades (e.g., amphibians), and some are relatively rapidly evolving and may be useful for closely-related species (e.g., within genera). This general method can be used whenever large-scale genomic data are available for an appropriate reference species (not necessarily within the focal clade). The method is also well suited for the development of intron regions for lower-level phylogenetic and phylogeographic studies. We provide an online database of alignments and suggested primers for approximately 85 NPCL that should be useful across vertebrates.     

Endochitinase gene-based phylogenetic analysis of Trichoderma

DNA sequences of the single-copy gene coding for the 42 kDa endochitinase enzyme (EC 3.2.1.14) were used for phylogenetic analysis in Trichoderma. A set of 12 primers was developed and the entire gene was sequenced for 16 strains, and nucleotide and deduced amino acid sequences were compared to data from GenBank for additional Trichoderma strains. Analysis of the sequences revealed parsimony informative variation from 2.4 to 43.6% depending on the part of the gene (exons/introns) and the taxonomic level considered. Results are discussed in comparison to previous data from ITS-1 and ITS-2 rDNA sequencing and suggest the 42 kDa endochitinase gene as a potential molecular marker for reconstructing phylogenetic relationships in the genus Trichoderma at species level.     

Development of PCR primer sets for intron 1 of the low-copy gene LEAFY in Davalliaceae

? Premise of the study: Primers were designed for amplifying intron 1 of the single-copy nuclear LEAFY gene for species of Davalliaceae. ? Methods and Results: New primer sets were designed and successfully amplified for intron 1 of the LEAFY gene in 13 species representing the five genera of Davalliaceae. The orthology of these sequences was further confirmed by phylogenetic analyses. Site variation in LEAFY intron 1 sequences across genera of the Davalliaceae and among accessions of the Humata repens complex were 18% and 8%, respectively. Such variation was greater than that for the cpDNA atpB-rbcL intergenic spacer region across the same taxa and accessions. ? Conclusions: Using our newly designed primers, intron 1 of the LEAFY gene could be amplified for all species tested. In addition, this single-copy, biparentally inherited, and quickly evolving region showed considerable potential for addressing infraspecific-level questions.     

A smörgåsbord of markers for avian ecology and evolution

The polymerase chain reaction has been a boon to the study of molecular ecology and population genetics of birds. But the nagging truth is that for many bird species, the number of polymerase chain reaction (PCR) primer pairs that one can pick off the shelf and expect to amplify their target loci with ease is frustratingly small. Now, studying DNA sequence variation in natural populations of birds just got a whole lot easier. This issue of Molecular Ecology reports a large-scale bioinformatics search for exonic sequences conserved between the chicken and zebra finch genomes and flanking polymorphic introns that has generated a staggering 242 PCR primer pairs that readily amplify their single-copy target loci in five avian species spanning ~100 million years of avian evolution ( Backström et al . 2008 ). As proof of principle, these primers have also been used to survey the genomic landscape in over 110 kb of intronic sequence in the collared flycatcher, a model species in ecology and evolution. These resources pave the way for easy multilocus study of evolving populations and lineages of birds, and bring the goal of quickly turning nonmodel species in to ecological genomic models tantalizingly close.     

Utility of Low-Copy Nuclear Gene Sequences in Plant Phylogenetics

Low-copy nuclear genes in plants are a rich source of phylogenetic information. They hold a great potential to improve the robustness of phylogenetic reconstruction at all taxonomic levels, especially where universal markers such as cpDNA and nrDNA are unable to generate strong phylogenetic hypotheses. Low-copy nuclear genes, however, remain underused in plant phylogenetic studies due to practical and theoretical complications in unraveling the evolutionary dynamics of nuclear gene families. The lack of the universal markers or universal PCR primers of low-copy nuclear genes has also hampered their phylogenetic utility. It has recently become clear that low-copy nuclear genes are particularly helpful in resolving close interspecific relationships and in reconstructing allopolyploidization in plants. Gene markers that are widely, if not universally, useful have begun to emerge. Although utilizing low-copy nuclear genes usually requires extra lab work such as designing PCR primers, PCR-cloning, and/or Southern blotting, rapid accumulation of gene sequences in the databases and advances in cloning techniques have continued to make such studies more feasible. With the growing number of theoretical studies devoted to the gene tree and species tree problem, a solid foundation for reconstructing complex plant phylogenies based on multiple gene trees began to build. It is also realized increasingly that fast evolving introns of the low-copy nuclear genes will provide much needed phylogenetic information around the species boundary and allow us to address fundamental questions concerning processes of plant speciation. Phylogenetic and molecular evolutionary analyses of developmentally important genes will add a new dimension to systematic and evolutionary studies of plant diversity.     

Utility of low-copy nuclear gene sequences in plant phylogenetics

Low-copy nuclear genes in plants are a rich source of phylogenetic information. They hold a great potential to improve the robustness of phylogenetic reconstruction at all taxonomic levels, especially where universal markers such as cpDNA and nrDNA are unable to generate strong phylogenetic hypotheses. Low-copy nuclear genes, however, remain underused in plant phylogenetic studies due to practical and theoretical complications in unraveling the evolutionary dynamics of nuclear gene families. The lack of the universal markers or universal PCR primers of low-copy nuclear genes has also hampered their phylogenetic utility. It has recently become clear that low-copy nuclear genes are particularly helpful in resolving close interspecific relationships and in reconstructing allopolyploidization in plants. Gene markers that are widely, if not universally, useful have begun to emerge. Although utilizing low-copy nuclear genes usually requires extra lab work such as designing PCR primers, PCR-cloning, and/or Southern blotting, rapid accumulation of gene sequences in the databases and advances in cloning techniques have continued to make such studies more feasible. With the growing number of theoretical studies devoted to the gene tree and species tree problem, a solid foundation for reconstructing complex plant phylogenies based on multiple gene trees began to build. It is also realized increasingly that fast evolving introns of the low-copy nuclear genes will provide much needed phylogenetic information around the species boundary and allow us to address fundamental questions concerning processes of plant speciation. Phylogenetic and molecular evolutionary analyses of developmentally important genes will add a new dimension to systematic and evolutionary studies of plant diversity.     

Targeting optimal introns for phylogenetic analyses in non-model taxa: experimental results in Asian pitvipers

Nuclear introns are increasingly used as phylogenetic markers. Here, we present a multidisciplinary approach towards optimal locus selection and amplification using Asian pitvipers as an example of a non‐model taxon, and raise the profile of length variant heterozygotes (LVHs) in intron loci. Taxon‐specific primers were identified using a bioinformatic approach, and also designed from existing exon primed, intron crossing (EPIC) primer amplifications. Eleven further universal EPIC primer pairs were assayed using a range of PCR optimization strategies. Taxon‐specific primers yielded the most consistent amplifications, but assaying a large number of universal EPIC primers yielded another appropriate locus for phylogenetic purposes. Modified Taq DNA polymerases such as JumpStart?Taq either significantly improved the specificity and yield of EPIC PCR amplifications (of low copy number nuclear targets), or resulted in amplifications that were not significantly worse than those derived from a generic Taq DNA polymerase. Finally, LVHs were detected in all loci that were sequenced suggesting that they are relatively common in introns. This study provides an efficient and cost effective template for the successful identification of intron markers for molecular systematics which is universally applicable to other non‐model taxon groups. © The Willi Hennig Society 2005.     

Discovery of paralogous nuclear gene sequences coding for the second-largest subunit of RNA polymerase II (RPB2) and their phylogenetic utility in gentianales of the asterids

Paralogous sequences of the RPB2 gene are demonstrated in the angiosperm order Gentianales. Two different copies were found by using different PCR primer pairs targeting a region that corresponds to exons 22-24 in the Arabidopsis RPB2 gene. One of the copies (RPB2-d) lacks introns in this region, whereas the other has introns at locations corresponding to those of green plants previously investigated. When analyzed with other available RPB2 sequences from this region, all 28 RPB2-d sequences obtained from the Gentianales and the four sequences from the Lamiales form a monophyletic group, together with a previously published tomato cDNA sequence. The substitution patterns, relative rates of change, and nucleotide compositions of the two paralogous RPB2 exon regions are similar, and none of them shows any signs of being a pseudogene. Although multiple copies of similar, paralogous sequences can confound phylogenetic interpretations, the lack of introns in RPB2-d make a priori homology assessment easy. The phylogenetic utility of RPB2-d within the Gentianales is evaluated in comparison with the chloroplast genes ndhF and rbcL. The hierarchical information in the RPB2-d region sequenced is more incongruent with that of the plastid genes than the plastid genes are with each other as determined by incongruence length difference tests. In contrast to the plastid genes, parsimony-informative third codon positions of RPB2 have a significantly higher rate of change than first and second positions. Topologically, the trees from the three genes are similar, and the differences are usually only weakly supported. In terms of support, RPB2 gives the highest jackknife support per sequenced nucleotide, whereas ndhF gives the highest Bremer support per sequenced nucleotide. The RPB2-d locus has the potential to be a valuable nuclear marker for determination of phylogenetic relationships within the euasterid I group of plants.     

PCR‐based isolation of multigene families: lessons from the avian MHC class IIB

The amount of sequence data available today highly facilitates the access to genes from many gene families. Primers amplifying the desired genes over a range of species are readily obtained by aligning conserved gene regions, and laborious gene isolation procedures can often be replaced by quicker PCR‐based approaches. However, in the case of multigene families, PCR‐based approaches bear the often ignored risk of incomplete isolation of family members. This problem is most prominent in gene families with highly variable and thus unpredictable number of gene copies among species, such as in the major histocompatibility complex (MHC). In this study, we (i) report new primers for the isolation of the MHC class IIB (MHCIIB) gene family in birds and (ii) share our experience with isolating MHCIIB genes from an unprecedented number of avian species from all over the avian phylogeny. We report important and usually underappreciated problems encountered during PCR‐based multigene family isolation and provide a collection of measures to help significantly improving the chance of successfully isolating complete multigene families using PCR‐based approaches.     

鲸类c - mos 基因的序列变异性及在系统发生分析中应用的初步研究

通过猪等的c-mos 基因保守区序列设计了用于扩增鲸类c-mos 基因的引物。应用此引物扩增并测定了齿鲸类5 个科12 个种546 bp 的c - mos 基因编码区序列。结果表明鲸类的c-mos 基因遗传变异水平较低。在系统发生重建中, 同科的物种聚为单系, 但不能很好地解决科下亚科间的关系, 提示c-mos 基因仅适于鲸类科级以上阶元的系统发生研究。     

Universal primers for fluorescent labelling of PCR fragments--an efficient and cost-effective approach to genotyping by fluorescence

Directly labelling locus‐specific primers for microsatellite analysis is expensive and a common limitation to small‐budget molecular ecology projects. More cost‐effective end‐labelling of PCR products can be achieved through a three primer PCR approach, involving a fluorescently labelled universal primer in combination with modified locus‐specific primers with 5′ universal primer sequence tails. This technique has been widely used but has been limited largely due to a lack of available universal primers suitable for co‐amplifying large numbers of size overlapping loci and without requiring locus‐specific PCR conditions to be modified. In this study, we report a suite of four high‐performance universal primers that can be employed in a three primer PCR approach for efficient and cost‐effective fluorescent end‐labelling of PCR fragments. Amplification efficiency is maximized owing to high universal primer Tm values (approximately 60+ °C) that enhance primer versatility and enable higher annealing temperatures to be employed compared with commonly used universal primers such as M13. We demonstrate that these universal primers can be combined with multiple fluorophores to co‐amplify multiple loci efficiently via multiplex PCR. This method provides a level of multiplexing and PCR efficiency similar to microsatellite fluorescent detection assays using directly labelled primers while dramatically reducing project costs. Primer performance is tested using several alternative PCR strategies that involve both single and multiple fluorophores in single and multiplex PCR across a wide range of taxa.     

Exon-primed, intron-crossing (EPIC) loci for five nuclear genes in deep-sea protobranch bivalves: primer design, PCR protocols and locus utility

We describe PCR primers and amplification protocols developed to obtain introns from conserved nuclear genes in deep-sea protobranch bivalves. Because almost no sequence data for protobranchs are publically available, mollusk and other protostome sequences from GenBank were used to design degenerate primers, making these loci potentially useful in other invertebrate taxa. Amplification and sequencing success varied across the test group of 30 species, and we present five loci spanning this range of outcomes. Intron presence in the targeted regions also varied across genes and species, often within single genera; for instance, the calmodulin and β-tubulin loci contained introns with high frequency, whereas the triose phosphate isomerase locus never contained an intron. In introns for which we were able to obtain preliminary estimates of polymorphism levels in single species, polymorphism was greater than traditional mitochondrial loci. These markers will greatly increase the ability to assess population structure in the ecologically important protobranchs, and may prove useful in other taxa as well.     

The D4 set: primers that target highly variable intron loops in plant chloroplast genomes

Chloroplast group II introns offer high-quality, rapidly evolving single-copy loci for comparative sequence analysis. These introns feature diagnostic secondary structures with loops that are among the least evolutionarily constrained sequence in plastomes. We exploited these structures to develop universal primers that amplify and sequence the large Domain IV (D4) loop in several angiosperm introns. With a single sequence read, we recover 300-600 nucleotides of highly variable sequence across angiosperms, with rates of change that are equal to or higher than many of the best known intergenic spacers in plant chloroplast genomes.     

EST derived PCR-based markers for functional gene homologues in cotton.

We investigated the utility of the Gossypium arboreum EST sequences in the GenBank database for developing PCR-based markers targeting known-function genes in cultivated tetraploid cottons, G. hirsutum and G. barbadense. Four hundred sixty-five randomly selected ESTs from this library were subjected to BLASTn search against all GenBank databases, of which putative function was assigned to 93 ESTs based on high nucleotide homology to previously studied genes. PCR primers were synthesized for 89 of the known-function ESTs. A total of 57 primer pairs amplified G. arboreum genomic DNA, but only 39 amplified in G. hirsutum and G. barbadense, suggesting that sequence divergence may be a factor causing non-amplification for some sites. DNA sequence analysis showed that most primer pairs were targeting the expected homologous loci. While the amplified products that were of larger size than the corresponding EST sequences contain introns, the primer pairs with a smaller amplicon than predicted from the flanking EST sequences did not amplify the expected orthologous gene sequences. Among the 39 primer pairs that amplified tetraploid cotton DNA, 3 detected amplicon size polymorphisms and 10 detected polymorphisms after digestion with one of six restriction enzymes. Ten of the polymorphic loci were subsequently mapped to an anchor RFLP map. Digestion of PCR-amplified sequences offers one means by which cotton genes can be mapped to their chromosomal locations more quickly and economically than by RFLP analysis.     

Ketosynthase Domain Probes Identify Two Subclasses of Fungal Polyketide Synthase Genes

Analysis of fungal polyketide synthase gene sequences suggested that these might be divided into two subclasses, designated WA-type and MSAS-type. Two pairs of degenerate PCR primers (LC1 and LC2c, LC3 and LC5c) were designed for the amplification of ketosynthase domain fragments from fungal PKS genes in each of these subclasses. Both primer pairs were shown to amplify one or more PCR products from the genomes of a range of ascomycetous Deuteromycetes and Southern blot analysis confirmed that the products obtained with each pair of primers emanated from distinct genomic loci. PCR products obtained from Penicillium patulum and Aspergillus parasiticus with the LC1/2c primer pair and from Phoma sp. C2932 with both primer pairs were cloned and sequenced; the deduced protein sequences were highly homologous to the ketosynthase domains of other fungal PKS genes. Genes from which LC1/2c fragments were amplified (WA-type) were shown by a phylogenetic analysis to be closely related to fungal PKS genes involved in pigment and aflatoxin biosynthetic pathways, whereas the gene from which the LC3/5c fragment was amplified (MSAS-type) was shown to be closely related to genes encoding 6-methylsalicylic acid synthase (MSAS). The phylogenetic tree strongly supported the division of fungal PKS genes into two subclasses. The LC-series primers may be useful molecular tools to facilitate the cloning of novel fungal polyketide synthase genes.