Simultaneous cloning of multiple nuclear genes by pooling PCR products of variable size: a cost‐effective method of improving efficiency in large‐scale genetic analyses

I present a simple approach to overcome the high cost and low efficiency of cloning polymerase chain reaction (PCR) products for individuals in wide‐scale population genetic analyses. The methodology reduces the number of cloning reactions per individual by engineering a suite of genetic markers that differ in size and pooling these PCR products prior to cloning. Alleles from each gene are then recovered by screening transformed bacterial colonies and identifying the inserts corresponding to each gene based on size. I demonstrate the utility of this technique by presenting the results I obtained from cloning four nuclear genes in 118 individuals from three species of sea urchins (Strongylocentrotus purpuratus, S. droebachiensis and S. pallidus). Of the 472 different PCR products I cloned, I recovered at least one allele for 432 of them (91.5%) by screening between 16 and 32 bacterial colonies for each individual. There existed a bias with respect to recovery efficiency: the two largest fragments (1130–800 bp) were recovered 100% of the time, while the two smaller fragments (580–650 bp) were recovered in 85.6% and 81.4% of the experiments, respectively. I discuss the promise of this application for wide‐scale genetic analyses.     

Sex identification of Eurasian otter (<Emphasis Type="Italic">Lutra lutra</Emphasis>) non-invasive DNA samples using ZFX/ZFY sequences

We developed a new PCR/RFLP system targeted to amplify and cut a segment of the ZFX/ZFY gene in non-invasive otter (Lutra lutra) samples. This assay produced one sex-specific fragment in females (XX genotypes) and two fragments in males (XY genotypes), and is intrinsically more reliable than alternative systems (e.g., SRY genes) that are based on the amplification of a single Y-linked fragment. The new primer pair correctly identified the sex of 23 DNA extracted from sexed otter tissues. This procedure was used to sex unknown DNA extracted from otter scats. A multi-tube approach led to identify the sex of 72/91 (79%) samples, using a minimum of two PCR replicates. This procedure is currently used in non-invasive genetic monitoring of Italian otter populations.     

Evaluation of cleaved amplified polymorphic sequence markers for <Emphasis Type="Italic">Chamaecyparis obtusa</Emphasis> based on expressed sequence tag information from <Emphasis Type="Italic">Cryptomeria japonica</Emphasis>

We have developed and evaluated sequence-tagged site (STS) primers based on expressed sequence-tag information derived from sugi (Cryptomeria japonica) for use in hinoki (Chamaecyparis obtusa), a species that belongs to a different family (although it appears to be fairly closely related to sugi). Of the 417 C. japonica STS primer pairs we screened, 120 (~30%) were transferable and provided specific PCR amplification products from 16 C. obtusa plus trees. We used haploid megagametophytes to investigate the homology of 80 STS fragments between C. obtusa and C. japonica and to identify orthologous loci. Nearly 90% of the fragments showed high (>70%) degrees of similarity between the species, and 35 STSs indicated homology to entries with the same putative function in a public DNA database. Of the 120 STS fragments amplified, 72 showed restriction fragment length polymorphisms; in addition, the CC2430 primers detected amplicon length polymorphism. We assessed the inheritance pattern of 27 cleaved amplified polymorphic sequence markers, using 20 individuals from the segregation population. All the markers analyzed were consistent with the marker inheritance patterns obtained from the screening panel, and no markers (except CC2716) showed significant (P<0.01) deviation from the expected segregation ratio. In total, 136 polymorphic markers were developed using C. japonica-based STS primers without any sequence modification. In addition, the applicability of STS-based markers developed in one species to other species was found to closely reflect the evolutionary distance between the species, which is roughly concordant with the difference between their rbcL sequences. We plan to use these markers for genetic studies in C. obtusa. Most of the markers should also provide reliable anchor loci for comparative mapping studies of the C. obtusa and C. japonica genomes.     

Overall Informativity,OI, in DNA Polymorphisms Revealed by inter-AluPCR: Detection of Genomic Rearrangements

We studied two systems of multilocus markers revealed by PCR using primers directing amplification betweenAlurepeats in a tail-to-tail orientation. Genomic polymorphisms were detected as the presence or absence of the electrophoretic bands representing DNA fragments of a given length. A total of 104 such fragments segregating as Mendelian markers in a panel of eight CEPH families were analyzed by two-point linkage analysis. Fifty-one of these fragments were localized with respect to CEPH markers; they represented 33 loci, 7 of which were multiallelic. Locus-specific oligonucleotides were developed and used as hybridization probes to identify the mapped loci within a complex pattern of inter-AluPCR products. A great proportion of inter-AluPCR polymorphisms represented length variants within amplified DNA segments, while others were presumably due to mutations within the priming sites. To describe the expected number of informative loci per typing experiment we introduced a parameter called overall informativity (OI), which provides a single measure of the multiplex ratio and the informativity of markers contributing to a multilocus system (OIof a single locus is equivalent to its heterozygosity and cannot exceed 0.5 for a biallelic codominant marker). HighOIvalues (5.8 and 11.5) of the two presented systems of inter-AluPCR markers of random chromosomal distribution render them suitable for mapping genomic rearrangements such as genomic deletions in tumoral tissues. This was illustrated by the detection of loss of heterozygosity in the 9q22–qter region in sporadic colon cancer.     

Amplification and detection of polymorphic sequence-tagged sites in<Emphasis Type="Italic">Lathyrus sativus</Emphasis>

A simple procedure was developed to convertLathyrus sativus defence-related expressed sequence tags (ESTs) into mappable genetic markers by using PCR. Twenty-nine STS primer pairs were generated on the basis of sequence information from anL. sativus cDNA library. These primers were used to screen for polymorphisms between 2L. sativus accessions, ATC 80878 and ATC 80407, resistant and susceptible, respectively, toMycosphaerella pinodes infection. All 29 primer pairs amplified PCR products in both accessions, 11 of which amplified multiple RAPD-like products. The remaining 18 primer pairs amplified single monomorphic products. Following cloning, sequencing, and database searches, 17 of 18 PCR products were confirmed to have amplified the targeted genome region. Ten of these 17 STS primer pairs revealed polymorphisms between ATC 80878 and ATC 80407 when PCR products were digested with a range of restriction endonucleases. These results suggest that the STS-based PCR analysis will be useful for generating informative molecular markers inL. sativus for future genome mapping experiments.     

Tagging and mapping of a rice gall midge resistance gene, <Emphasis Type="Italic">Gm8</Emphasis>, and development of SCARs for use in marker-aided selection and gene pyramiding

Using amplified fragment length polymorphisms (AFLPs) and random amplified polymorphic DNAs (RAPDs), we have tagged and mapped Gm8, a gene conferring resistance to the rice gall midge (Orseolia oryzae), a major insect pest of rice, onto rice chromosome 8. Using AFLPs, two fragments, AR257 and AS168, were identified that were linked to the resistant and susceptible phenotypes, respectively. Another resistant phenotype-specific marker, AP19₅₈₇, was also identified using RAPDs. SCAR primers based on the sequence of the fragments AR257 and AS168 failed to reveal polymorphism between the resistant and the susceptible parents. However, PCR using primers based on the regions flanking AR257 revealed polymorphism that was phenotype-specific. In contrast, PCR carried out using primers flanking the susceptible phenotype-associated fragment AS168 produced a monomorphic fragment. Restriction digestion of these monomorphic fragments revealed polymorphism between the susceptible and resistant parents. Nucleotide BLAST searches revealed that the three fragments show strong homology to rice PAC and BAC clones that formed a contig representing the short arm of chromosome 8. PCR amplification using the above-mentioned primers on a larger population, derived from a cross between two indica rice varieties, Jhitpiti (resistant parent) and TN1 (susceptible parent), showed that there is a tight linkage between the markers and the Gm8 locus. These markers, therefore, have potential for use in marker-aided selection and pyramiding of Gm8 along with other previously tagged gall midge resistance genes [Gm2, Gm4(t), and Gm7].The nucleotide sequence data reported here will appear in the EMBL, GenBank and DDBJ nucleotide sequence databases under the accession numbers AY545920–AY545923     

Discordance in the distribution of markers of different inheritance systems (nDNA,mtDNA, and Chromosomes) in the superspecies complex <Emphasis Type="Italic">Mus musculus</Emphasis> as a result of extensive hybridization in primorye

The genetic structure of eight Mus musculus L. populations in Primorskii krai was studied with the use of taxon-specific markers of different inheritance systems: nDNA (RAPD), mtDNA (D-loop), and chromosomes. The results obtained demonstrate that although the compared nuclear marker characteristics (nDNA and chromosomes) have the same basis they are not linke with each other and, moreover, are often mutually inconsistent. Discordance in the inheritance of the marker characteristics in most of the animals studied is a result of extensive hybridization involving two to four house mouse subspecies. To identify taxonspecific nuclear markers revealed by RAPD, some RAPD PCR products were cloned, and their localization on chromosomes was determined. It was found that some fragments similar in size consist of two different comigrating sequences that are localized on different chromosomes and belong to different subspecies. All sequenced anonymous markers are localized in protein-coding genes. The functions of genes containing the marker sequences have been established. Differences in the taxon-specific RAPD fragments are associated with changes in the structure of important functional genes, and this can be considered as a significant genetic marker.     

Novel DNA markers for taxonomic characterization and identification of <Emphasis Type="Italic">Fusarium</Emphasis> species

Ten Fusarium sporotrichioides strains from different geographic regions were analyzed by RAPD in order to detect DNA loci potentially suitable as new markers for taxonomic characterization and identification of toxigenic Fusarium fungi. Three monomorphic fragments were selected from PCR amplificates obtained with one of the standard RAPD primers and sequenced. Analysis of the sequences enabled the development of specific SCAR markers for identification of Fusarium fungi at the level of species groups characterized by similar profiles of produced mycotoxins.     

Development of a microsatellite marker set applicable to genome-wide screening of cynomolgus monkeys (<Emphasis Type="Italic">Macaca fascicularis</Emphasis>)

To develop a microsatellite marker set applicable to genome-wide screening of cynomolgus monkeys (Macaca fascicularis), 148 microsatellite markers were selected from the human genome database. The polymorphisms and inheritance of PCR products were determined by screening twenty unrelated monkeys and by analysis of three families, respectively. As a result, 106 primers (72%) gave PCR products of the size expected for humans and rhesus monkeys. Among these products, polymorphism and single-gene inheritance in cynomolgus monkeys was observed for 66 markers (62%). The average number of alleles at the 66 polymorphic loci was 5.86 (range 2–10), and average heterozygosity was 0.63 (range 0.10–0.88). This is the first report of microsatellite markers for cynomolgus monkeys. Chromosomal mapping of these markers is now in progress.     

Novel molecular markers differentiate Oncorhynchus mykiss (rainbow trout and steelhead) and the O. clarki (cutthroat trout) subspecies

A suite of 26 PCR‐based markers was developed that differentiates rainbow (Oncorhynchus mykiss) and coastal cutthroat trout (O. clarki clarki). The markers also differentiated rainbow from other cutthroat trout subspecies (O. clarki), and several of the markers differentiated between cutthroat trout subspecies. This system has numerous positive attributes, including: nonlethal sampling, high species‐specificity and products that are easily identified and scored using agarose gel electrophoresis. The methodology described for developing the markers can be applied to virtually any system in which numerous markers are desired for identifying or differentiating species or subspecies.     

Fine mapping and targeted SNP survey using rice-wheat gene colinearity in the region of the <Emphasis Type="Italic">Bo1</Emphasis> boron toxicity tolerance locus of bread wheat

Toxicity due to high levels of soil boron (B) represents a significant limitation to cereal production in some regions, and the Bo1 gene provides a major source of B toxicity tolerance in bread wheat (Triticum aestivum L.). A novel approach was used to develop primers to amplify and sequence gene fragments specifically from the Bo1 region of the hexaploid wheat genome. Single-nucleotide polymorphisms (SNPs) identified were then used to generate markers close to Bo1 on the distal end of chromosome 7BL. In the 16 gene fragments totaling 19.6 kb, SNPs were observed between the two cultivars Cranbrook and Halberd at a low frequency (one every 613 bp). Furthermore, SNPs were distributed unevenly, being limited to only two genes. In contrast, RFLP provided a much greater number of genetic markers, with every tested gene identifying polymorphism. Bo1 previously known only as a QTL was located as a discrete Mendelian locus. In total, 28 new RFLP, PCR and SSR markers were added to the existing map. The 1.8 cM Bo1 interval of wheat corresponds to a 227 kb section of rice chromosome 6L encoding 21 predicted proteins with no homology to any known B transporters. The co-dominant PCR marker AWW5L7 co-segregated with Bo1 and was highly predictive of B tolerance status within a set of 94 Australian bread wheat cultivars and breeding lines. The markers and rice colinearity described here represent tools that will assist B tolerance breeding and the positional cloning of Bo1. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

High-resolution genetic mapping of the pepper resistancelocus Bs3 governing recognition of the Xanthomonascampestris pv vesicatora AvrBs3 protein

The pepper (Capsicum annuum) Bs3 gene confers resistance to Xanthomonas campestris pv vesicatoria strains expressing the avirulence protein AvrBs3. Using amplified fragment length polymorphism (AFLP) and bulked DNA templates from resistant and susceptible plants we identified markers linked to Bs3 and defined a 2.1-cM interval containing the target gene. Bs3-linked AFLP fragments were cloned and conformity of isolated PCR products with the desired markers was determined by hybridisation to membrane-bound AFLP reactions. AFLP markers flanking the target gene were converted into locus-specific PCR-based markers. These markers were employed for the analysis of 790 plants segregating for Bs3, resulting in a linkage map with a genetic resolution of 0.13 cM. Mapping of Bs3-linked markers in tomato placed them to a syntenic interval on tomato chromosome 2. Received: 15 October 1999 / Accepted: 29 November 1999     

Characterisation and marker development for low molecular weight glutenin genes from<Emphasis Type="Italic"> Glu-A3</Emphasis> alleles of bread wheat (<Emphasis Type="Italic">Triticum aestivum.</Emphasis> L)

PCR was used to amplify low-molecular-weight (LMW) glutenin genes from the Glu-A3 loci of hexaploid wheat cultivars containing different Glu-A3 alleles. The complete coding sequence of one LMW glutenin gene was obtained for each of the seven alleles Glu-A3a to Glu-A3g. Chromosome assignment of PCR products using Chinese Spring nulli-tetrasomic lines confirmed the amplified products were from chromosome 1A. All sequences were classified as LMW-i-type genes based on the presence of an N-terminal isoleucine residue and eight cysteine residues located within the C-terminal domain of the predicted, mature amino acid sequence. All genes contained a single uninterrupted open reading frame, including the sequence from the Glu-A3e allele, for which no protein product has been identified. Comparison of LMW glutenin gene sequences obtained from different alleles showed a wide range of sequence identity between the genes, with between 1 and 37 single nucleotide polymorphisms and between one and five insertion/deletion events between genes from different alleles. Allele-specific PCR markers were designed based on the DNA polymorphisms identified between the LMW glutenin genes, and these markers were validated against a panel of cultivars containing different Glu-A3 alleles. This collection of markers represents a valuable resource for use in marker-assisted breeding to select for specific alleles of this important quality-determining locus in bread wheat.Communicated by P. Langridge     

PCR-CTPP: a rapid and reliable genotyping technique based on ZFX/ZFY alleles for sex identification of tiger (<Emphasis Type="Italic">Panthera tigris</Emphasis>) and four other endangered felids

An inexpensive, time-saving and reliable method, polymerase chain reaction with confronting two-pair primers (PCR-CTPP), was developed for sex identification in tiger (Panthera tigris) based on zinc finger alleles (ZFX/ZFY). A site of “C/G” transversion representing fixed differences that discriminated between ZFX and ZFY exons among felids was identified for primers designing. This primer set was successfully tested on samples including blood, shed hairs, dried skin, and stool which contained potential contamination caused by prey DNA. Cross species tests shown that this primer set was also useful for sex identification in four other endangered felids.     

Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias

Studies of insect assemblages are suited to the simultaneous DNA‐based identification of multiple taxa known as metabarcoding. To obtain accurate estimates of diversity, metabarcoding markers ideally possess appropriate taxonomic coverage to avoid PCR‐amplification bias, as well as sufficient sequence divergence to resolve species. We used in silico PCR to compare the taxonomic coverage and resolution of newly designed insect metabarcodes (targeting 16S) with that of existing markers [16S and cytochrome oxidase c subunit I (COI)] and then compared their efficiency in vitro. Existing metabarcoding primers amplified in silico <75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided >90% coverage. Furthermore, metabarcodes targeting COI appeared to introduce taxonomic PCR‐amplification bias, typically amplifying a greater percentage of Lepidoptera and Diptera species, while failing to amplify certain orders in silico. To test whether bias predicted in silico was observed in vitro, we created an artificial DNA blend containing equal amounts of DNA from 14 species, representing 11 insect orders and one arachnid. We PCR‐amplified the blend using five primer sets, targeting either COI or 16S, with high‐throughput amplicon sequencing yielding more than 6 million reads. In vitro results typically corresponded to in silico PCR predictions, with newly designed 16S primers detecting 11 insect taxa present, thus providing equivalent or better taxonomic coverage than COI metabarcodes. Our results demonstrate that in silico PCR is a useful tool for predicting taxonomic bias in mixed template PCR and that researchers should be wary of potential bias when selecting metabarcoding markers.     

Sex identification of northern ungulates using low quality and quantity DNA

We evaluated PCR primer sets to determine the most effective technique for identifying sex of northern ungulates. We sought markers that required only a single pair of primers to amplify both X- and Y-linked alleles; that amplified X- and Y-linked products that were easily distinguishable using agarose gel electrophoresis; and that produced short amplicons amenable to amplification using DNA of poor quality and low quantity, as is often found in non-invasively collected samples such as feces. Primer pairs KY1/KY2 and SE47/SE48, which amplify X- and Y-specific alleles of the amelogenin gene, met our criteria and were tested for moose (Alces alces), mountain goat (Oreamnos americanus), Sitka black-tailed deer (Odocoileus hemionus sitkensis), and caribou (Rangifer tarandus). KY primers amplified shorter PCR products than did SE primers; moreover, SE primers inconsistently amplified certain Y-chromosome products, creating potential for misidentification of sex. DNA fragments amplified using KY primers were sequenced for each species, allowing us to characterize a 45-bp deletion for Y-linked alleles (136-bp product) relative to X-linked alleles (181-bp product) in all species and a 9-bp deletion in the X-linked allele of moose relative to other species. This is the first sex-determination technique using PCR reported for several ungulate species of Alaska. Although other protocols exist for cervids and bovids, this is the first report of markers meeting the aforementioned criteria for Odocoileus, the most abundant and intensively managed genus of large mammals in North America.     

PCR-based markers for the powdery mildew resistance gene<Emphasis Type="Italic"> Pm4a</Emphasis> in wheat

Gene tagging is the basis of marker-assisted selection and map-based cloning. To develop PCR-based markers for Pm4a, a dominant powdery mildew resistance gene of wheat, we surveyed 46 group 2 microsatellite markers between Pm4a near-isogenic line (NIL) CI 14124 and the recurrent parent Chancellor (Cc). One of the markers, gwm356, detected polymorphism and was used for genotyping an F₂ population of 85 plants derived from CI 14124 × Cc. Linkage mapping indicated that Xgwm356 was linked to Pm4a at a distance of 4.8 cM. To identify more PCR-based markers for Pm4a, we also converted the restriction fragment length polymorphism marker BCD1231 linked to it into a sequence-tagged site (STS) marker. The STS primer designed based on the end sequences of BCD1231 amplified an approximately 1.6-kb monomorphic band in both parents. Following digestion of the products with the four-cutter enzymes HaeIII and MspI, it was discovered that the band from CI 14124 consisted of at least two products, one of which had a digestion pattern different from the band from Cc. In the F₂ population, the cleaved polymorphism revealed by the STS marker between the parents co-segregated with powdery mildew resistance. To design Pm4a-specific PCR markers, the 1.6-kb band from Cc and the fragment associated with Pm4a in CI 14124 were sequenced and compared. Based on these sequences a new PCR marker was identified, which detected a 470-bp product only in the Pm4a-containing lines. These PCR-based markers provide a cost-saving option for marker-assisted selection of Pm4a.Communicated by F. Salamini     

Detection of section-specific random amplified polymorphic DNA (RAPD) markers in Lilium

Random amplified polymorphic DNA (RAPD) markers were utilized for the identification of Lilium species and inter-specific hybrids. The optimum annealing temperature of the polymerase chain reaction (PCR) for the RAPD assay in Lilium was 54 °C, which is relatively higher than the temperature used for other genera reported by previous researchers. Among 76 primers used to amplify genomic DNA by PCR, 18 primers (24%) generated polymorphic DNA fragments in Lilium species and hybrids. Cultivars were also identified by RAPD markers. Some amplified fragments were unique to species of each section and to hybrids derived from these species; that is, they were the section-specific DNA markers. Sections, Sinomartagon, Leucolirion b, Leucolirion a and Archelirion could be identified by 6 section-specific markers amplified with five primers. Seven inter-section hybrids showed the section-specific bands of both parental sections, indicating that these markers would be useful for identifying the parental sections of inter-section hybrids.     

Polymorphic sequence‐characterized codominant loci in the chestnut blight fungus,Cryphonectria parasitica

Studies on the population biology of the chestnut blight fungus, Cryphonectria parasitica, have previously been carried out with dominant restriction fragment length polymorphism (RFLP) fingerprinting markers. In this study, we describe the development of 11 codominant markers from randomly amplified polymorphic DNAs (RAPDs). RAPD fragments were cloned and sequenced, and polymerase chain reaction (PCR) primers were designed flanking insertions/deletions. Primers labelled with fluorescent dyes were combined in multiplex reactions to assay five or six loci simultaneously in a capillary sequencing system. These codominant markers have the potential to complement RFLP methods for studying C. parasitica.