Detection and isolation of nuclear haplotypes by PCR-SSCP

SSCP (single-strand conformational polymorphism) is used widely in the field of human biomedicine, but its potential as a population genetics tool for the recovery of nuclear gene genealogies remains to be realized. We describe and illustrate a use for SSCP in the physical isolation of nuclear haplotypes that circumvents several difficulties associated with more conventional cloning procedures. The DNA sequence can be determined directly from the isolated haplotypes and used for phylogenetic inference. SSCP provides a convenient first step toward generating nuclear genealogies for population studies.     

The challenge for genetic epidemiologists: how to analyze large numbers of SNPs in relation to complex diseases

Background

Molecular genetic approaches have much to offer population biology. Despite recent advances, convenient techniques to develop and screen highly-resolving markers can be limiting for some applications and taxa. We describe an improved PCR-based, cloning-free, nuclear marker development procedure, in which single-stranded conformation polymorphism (SSCP) plays a central role. Sequence-variable alleles at putative nuclear loci are simultaneously identified and isolated from diploid tissues. Based on a multiple allele alignment, locus-specific primers are designed in conserved regions, minimizing 'null' alleles. Using two undescribed endemic Australian Collembola as exemplars, we outline a comprehensive approach to generating and validating suites of codominant, sequence-yielding nuclear loci for previously unstudied invertebrates.

Results

Six markers per species were developed without any baseline genetic information. After evaluating the characteristics of each new locus via SSCP pre-screening, population samples were genotyped on the basis of either DNA sequence, restriction site, or insertion/deletion variation, depending on which assay was deemed most appropriate. Polymorphism was generally high (mean of nine alleles per locus), and the markers were capable of resolving population structuring over very fine spatial scales (<100 km). SSCP coupled with targeted DNA sequencing was used to obtain genotypic, genic and genealogical information from six loci (three per species). Phylogeographic analysis identified introns as being most informative.

Conclusion

The comprehensive approach presented here feasibly overcomes technical hurdles of (i) developing suitably polymorphic nuclear loci for non-model organisms, (ii) physically isolating nuclear allele haplotypes from diploid tissues without cloning, and (iii) genotyping population samples on the basis of nuclear DNA sequence variation.     

The single strand conformational analysis of cattle and human single nucleotide polymorphisms may be biased towards specific sequence motifs that minimize local secondary structure of single strand DNA

Single strand conformational polymorphisms (SSCP) resulting from point mutations were found to be associated preferentially with two DNA sequence motifs. These motifs are (1) three or more of the same base but in which the polymorphism is not due to length variation and (2) a region of polypurine or polypyrimidine bases. These motifs were identified after SSCP alleles from cattle were sequenced. The sequence difference and flanking sequence for each single nucleotide polymorphism are shown. The motifs were also found in SSCP from humans chosen at random from the literature, in which the alleles had been sequenced. There is a low level of complementarity of adjacent bases in these motifs and they should represent regions of low secondary structure in the single stranded DNA. Regions of high secondary structure, such as palindromes, were found in the same sample to have allelic variation that was not detected by SSC analysis. These results give a rule of thumb for selecting the particular part of a DNA fragment to be selected for testing for polymorphisms, but this rule clashes with rules used to design primers to amplify sequences using the PCR, namely, minimise hydrogen bonding within and between primers and reduce self-complementarity.     

Direct comparison of selected methods for genetic categorisation of Cryptosporidium parvum and Cryptosporidium hominis species

A study was undertaken to compare the performance of five different molecular methods (available in four different laboratories) for the identification of Cryptosporidium parvum and Cryptosporidium hominis and the detection of genetic variation within each of these species. The same panel of oocyst DNA samples derived from faeces (n=54; coded blindly) was sent for analysis by: (i) DNA sequence analysis of a fragment of the HSP70 gene; (ii) DNA sequence analysis and the ssrRNA gene in laboratory 1; (iii) single-strand conformation polymorphism analysis of part of the ssrRNA; (iv) SSCP analysis of the second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA region in laboratory 2; (v) 60 kDa glycoprotein (gp60) gene sequencing with prior species determination using PCR with restriction fragment length polymorphism analysis of the ssrRNA gene in laboratory 3; and (vi) multilocus genotyping at three microsatellite markers in laboratory 4. For detecting variation within C. parvum and C. hominis, SSCP analysis of ITS-2 was considered to have superior utility and determined 'subgenotypes' in samples containing DNA from both species. SSCP was also most cost effective in terms of time, cost and consumables. Sequence analysis of gp60 and microsatellite markers ML1, ML2 and 'gp15' provided good comparators for the SSCP of ITS-2. However, applicability of these methods to other Cryptosporidium species or genotypes and to environmental samples needs to be evaluated. This trial provided, for the first time, a direct comparison of multiple methods for the genetic characterisation of C. parvum and C. hominis samples. A protocol has been established for the international distribution of samples for the characterisation of Cryptosporidium. This can be applied in further evaluation of molecular methods by investigation of a larger number of unrelated samples to establish sensitivity, typability, reproducibility and discriminatory power based on internationally accepted methods for evaluation of microbial typing schemes.     

Minisatellite isoalleles can be distinguished by single-stranded conformational polymorphism analysis in agarose gels.

Minisatellite isoallelism, i.e. the occurrence of minisatellite alleles with different internal sequence composition but indistinguishable length, is a common limitation of minisatellite allele length analysis. Internal sequence variation can be used to distinguish such isoalleles, provided that detailed sequence knowledge of its basis is available. We now show that minisatellite isoalleles can also be simply resolved by single-stranded conformational polymorphisms (SSCP) arising during agarose gel electrophoresis. SSCP on agarose gels can be used to distinguish minisatellite isoalleles either after PCR amplification, or by standard Southern blot analysis of genomic DNA.     

Electrophoretic analysis of genetic variability within Cryptosporidium parvum from imported and autochthonous cases of human cryptosporidiosis in the United Kingdom

Cryptosporidium parvum oocyst DNA samples (n = 184) from humans with cryptosporidiosis contracted during foreign travel or during outbreaks in the United Kingdom were characterized genetically and categorized by single-strand conformation polymorphism (SSCP)-based analysis of the small-subunit gene (pSSU) (approximately 300 bp) and second internal transcribed spacer (pITS-2) (approximately 230 bp) of nuclear ribosomal DNA. The two recognized genotypes (types 1 and 2) of C. parvum could be readily differentiated by a distinct electrophoretic shift in the pSSU SSCP profile, associated with a nucleotide difference of approximately 1.3 to 1.7%. Of the 102 samples from cases contracted during foreign travel, 88 (86.3%) were identified as C. parvum type 1 and 14 (13.7%) were identified as type 2. For outbreak samples, unequivocal differentiation between type 1 (n = 20; one child nursery outbreak) and type 2 (n = 62; two waterborne outbreaks) was also achieved. Nucleotide variation in pITS-2 (both within and among samples representing each genotype) was substantially greater (10 to 13 different profiles for each genotype, relating to sequence differences of approximately 1 to 42%) than that in pSSU. SSCP analysis of pITS-2 for all samples revealed that some profiles had a broad geographical distribution whereas others were restricted to particular locations, suggesting a link between some subgenotypes and the geographical origin or source. Comparative denaturing polyacrylamide gel electrophoretic analysis revealed the same genotypic identification and a similar subgenotypic classification of samples as SSCP analysis. The findings of this study, particularly the detection of intragenotypic variation by SSCP, should have significant diagnostic implications for investigating transmission patterns and the monitoring of outbreaks.     

Polymorphic molecular markers from anonymous nuclear DNA for genetic analysis of populations

A simple method is presented for developing polymorphic, anonymous DNA markers suitable for population genetic studies. Anonymous DNA fragments are screened for sequence variability using a common mutation detection technique (single strand conformation polymorphism analysis; SSCP) and locus-specific PCR primers are designed for polymorphic DNA fragments. Detection of the markers by SSCP analysis coupled with sequence analysis of SSCP variants allows rapid screening while retaining information about the genealogical relationship among alleles. Variability detected for six markers was assessed in rainbow trout Oncorhynchus mykiss and was compared with variability detected by similar analysis of intron loci. Between three and 12 distinct alleles were observed at each marker locus, and average within-population heterozygosity ranged from 0.12 to 0.44. Advantages and limitations of the methodology for population genetic analysis are discussed.     

Optimization of single-strand conformation polymorphism and sequence analysis of the mitochondrial control region in Pagellus bogaraveo (Sparidae, Teleostei): rationalized tools in fish population biology

We report the isolation and sequencing of 400–550 base pairs (bp) of the mitochondrial DNA (mtDNA) control region of eight species of Sparidae (Perciformes, Teleostei). This sequence information allowed us to design specific primers to one of these species (Pagellus bogaraveo). The new set of primers was used to test a rationalized approach to study the mtDNA nucleotide variability at the intraspecific level. The single-strand conformation polymorphism (SSCP) technique was applied to detect sequence variation in two non-overlapping fragments of the control region of 32 individuals of P. bogaraveo. To assess the sensitivity of the method, the nucleotide sequence of the analysed region was determined for all the specimens. The results showed that, for one of the two fragments, SSCP analysis was able to detect 100% of the underlying genetic variability. In sharp contrast, nucleotide variation of the second DNA fragment was completely unresolved by SSCP under different experimental conditions. This suggests that the resolution power of SSCP is crucially dependent on the nature of the fragment subjected to the analysis; therefore, a preliminary test of the sensitivity of the method should be performed on each specific DNA fragment before starting a large-scale survey. A rationalized approach, combining the SSCP technique and a simplified sequencing procedure, is proposed for studying intraspecific polymorphism at the mtDNA control region in fish.     

Linkage disequilibria and the site frequency spectra in the su(s) and su(w(a)) regions of the Drosophila melanogaster X chromosome

Over the last decade, surveys of DNA sequence variation in natural populations of several Drosophila species and other taxa have established that polymorphism is reduced in genomic regions characterized by low rates of crossing over per physical length. Parallel studies have also established that divergence between species is not reduced in these same genomic regions, thus eliminating explanations that rely on a correlation between the rates of mutation and crossing over. Several theoretical models (directional hitchhiking, background selection, and random environment) have been proposed as population genetic explanations. In this study samples from an African population (n = 50) and a European population (n = 51) were surveyed at the su(s) (1955 bp) and su(w(a)) (3213 bp) loci for DNA sequence polymorphism, utilizing a stratified SSCP/DNA sequencing protocol. These loci are located near the telomere of the X chromosome, in a region of reduced crossing over per physical length, and exhibit a significant reduction in DNA sequence polymorphism. Unlike most previously surveyed, these loci reveal substantial skews toward rare site frequencies, consistent with the predictions of directional hitchhiking and random environment models and inconsistent with the general predictions of the background selection model (or neutral theory). No evidence for excess geographic differentiation at these loci is observed. Although linkage disequilibrium is observed between closely linked sites within these loci, many recombination events in the genealogy of the sampled alleles can be inferred and the genomic scale of linkage disequilibrium, measured in base pairs between sites, is the same as that observed for loci in regions of normal crossing over. We conclude that gene conversion must be high in these regions of low crossing over.     

Extensive intraspecific polymorphism detected by SSCP at the nuclear C-mos gene in the endemic Iberian lizard Lacerta schreiberi

C-mos is a highly conserved intronless gene that has proved useful in the analysis of ancient phylogenetic relationships within vertebrates. We selected the Iberian endemic Schreiber's green lizard (Lacerta schreiberi) that persisted in allopatric refugia since the late Pliocene to investigate the utility of the C-mos nuclear gene for intraspecific phylogeographic studies. Our combination of DNA sequencing with the high resolving power of single-strand conformational polymorphism (SSCP) effectively discriminated four common alleles showing strong population structuring (F(ST) = 0.46). In addition, reconstruction of allele phylogenetic relationships further improved our understanding of C-mos spatial patterns of variation and allowed a comparison with previously described mitochondrial DNA data. Finally, limited sequencing of an extended C-mos fragment in six additional Lacerta species showed extensive polymorphism, to our knowledge representing a rare example of variation in a highly conserved nuclear gene.     

Invasion genetics of the Mediterranean fruit fly: variation in multiple nuclear introns

Biological invasions generally start from low initial population sizes, leading to reduced genetic variation in nuclear and especially mitochondrial DNA. Consequently, genetic approaches for the study of invasion history and population structure are difficult. An extreme example is the Mediterranean fruit fly, Ceratitis capitata (Medfly), for which successive invasions during this century have resulted in a loss of 60% of ancestral genetic variation in isozymes and 75% of variation in mitochondrial DNA. Using Medflies as an example, we present a new approach to invasion genetics that measures DNA sequence variation within introns from multiple nuclear loci. These loci are so variable that even relatively recently founded Medfly populations within California and Hawaii retain ample genetic diversity. Invading populations have only lost 35% of the ancestral genetic variation. Intron variation will allow high-resolution genetic characterization of invading populations in both natural and managed systems, although non-equilibrium methods of analysis may be necessary if the genetic diversity represents sorting ancestral polymorphism.     

Mutational analysis of SRY: nonsense and missense mutations in XY sex reversal

Summary XY females (n=17) were analysed for mutations in SRY (sex-determining region Y gene), a gene that has recently been equated with the testis determining factor (TDF). SRY sequences were amplified by the polymerase chain reaction (PCR) and analysed by both the single strand conformational polymorphism assay (SSCP) and DNA sequencing. The DNA from two individuals gave altered SSCP patterns; only these two individuals showed any DNA sequence variation. In both cases, a single base change was found, one altering a tryptophan codon to a stop codon, the other causing a glycine to arginine amino acid substitution. These substitutions lie in the high mobility group (HMG)-related box of the SRY protein, a potential DNA-binding domain. The corresponding regions of DNA from the father of one individual and the paternal uncle of the other, were sequenced and found to be normal. Thus, in both cases, sex reversal is associated with de novo mutations in SRY. Combining this data with two previously published reports, a total of 40 XY females have now been analysed for mutations in SRY. The number of de novo mutations in SRY is now doubled to four, adding further strength to the argument that SRY is TDF.     

Development of primer sets for PCR amplification of the PgiC gene in ferns

Polymerase chain reaction (PCR)-based nuclear DNA markers were developed for fern species. We first determined the partial nucleotide sequence of cDNA of the pgiC gene encoding cytosolic phosphoglucose isomerase from Dryopteris caudipinna, and then PCR primers for exon-primed, intron-crossing (EPIC) amplifications were designed. The EPIC primers are universally applicable to the most derived indusiate fern families such as Dryopteridaceae, Thelypteridaceae, and Woodsiaceae. The PCR products of primers 14F/16R containing two introns are moderate in size (534 bp-ca.1000 bp) and are possibly of value in phylogenetic reconstruction at specific and generic levels. Codominant nuclear DNA markers applicable to the estimation of mating systems and other population genetic studies were also developed by a combination of single-strand conformation polymorphism (SSCP) and EPIC amplification using primers 14F/15R and 15F/16R. In order to provide a case study using these markers, allelic variation of PCR products using 15F/16R was examined in populations of Arachniodes standishii (Dryopteridaceae).     

A survey of chromosomal and nucleotide sequence variation in Drosophila miranda

There have recently been several studies of the evolution of Y chromosome degeneration and dosage compensation using the neo-sex chromosomes of Drosophila miranda as a model system. To understand these evolutionary processes more fully, it is necessary to document the general pattern of genetic variation in this species. Here we report a survey of chromosomal variation, as well as polymorphism and divergence data, for 12 nuclear genes of D. miranda. These genes exhibit varying levels of DNA sequence polymorphism. Compared to its well-studied sibling species D. pseudoobscura, D. miranda has much less nucleotide sequence variation, and the effective population size of this species is inferred to be several-fold lower. Nevertheless, it harbors a few inversion polymorphisms, one of which involves the neo-X chromosome. There is no convincing evidence for a recent population expansion in D. miranda, in contrast to D. pseudoobscura. The pattern of population subdivision previously observed for the X-linked gene period is not seen for the other loci, suggesting that there is no general population subdivision in D. miranda. However, data on an additional region of period confirm population subdivision for this gene, suggesting that local selection is operating at or near period to promote differentiation between populations.     

Conflicting patterns of mitochondrial and nuclear DNA diversity in Phylloscopus warblers

Molecular variation is often used to infer the demographic history of species, but sometimes the complexity of species history can make such inference difficult. The willow warbler, Phylloscopus trochilus, shows substantially less geographical variation than the chiffchaff, Phylloscopus collybita, both in morphology and in mitochondrial DNA (mtDNA) divergence. We therefore predicted that the willow warbler should harbour less nuclear DNA diversity than the chiffchaff. We analysed sequence data obtained from multiple samples of willow warblers and chiffchaffs for the mtDNA cytochrome b gene and four nuclear genes. We confirmed that the mtDNA diversity among willow warblers is low (pi = 0.0021). Sequence data from three nuclear genes (CHD-Z, AFLP-WW1 and MC1R) not linked to the mitochondria demonstrated unexpectedly high nucleotide diversity (pi values of 0.0172, 0.0141 and 0.0038) in the willow warbler, on average higher than the nucleotide diversity for the chiffchaff (pi values of 0.0025, 0.0017 and 0.0139). In willow warblers, Tajima's D analyses showed that the mtDNA diversity, but not the nuclear DNA diversity, has been reduced relative to the neutral expectation of molecular evolution, suggesting the action of a selective sweep affecting the maternally inherited genes. The large nuclear diversity seen within willow warblers is not compatible with processes of neutral evolution occurring in a population with a constant population size, unless the long-term effective population size has been very large (N(e) > 10(6)). We suggest that the contrasting patterns of genetic diversity in the willow warbler may reflect a more complex evolutionary history, possibly including historical demographic fluctuations or historical male-biased introgression of nuclear genes from a differentiated population of Phylloscopus warblers.     

Polymorphism of the IGHA gene in sheep

Genetic variation in immunoglobulin A, the most abundant immunoglobulin in mammalian cells, has not been reported in ruminants. In this study, variation in the immunoglobulin heavy alpha chain constant gene (IGHA) of sheep was investigated by amplification of a fragment that included the hinge coding sequence, followed by single-strand conformational polymorphism (SSCP) analysis and DNA sequencing. Three novel sequences, each characterized by unique SSCP banding patterns, were identified. One or two sequences were detected in individual sheep and all the sequences identified shared high homology to the published ovine and bovine IGHA sequences, suggesting that these sequences represent allelic variants of the IGHA gene in sheep. Sequence alignment showed that these sequences differed mainly in the 3′ end of exon 1 and in the coding sequence of the hinge region. There was either a deletion or an insertion of two codons in the hinge coding region in these allelic variants. Codon usage in the hinge coding region was quite different from that in the non-hinge coding regions of the gene, suggesting different evolution of the IGHA hinge sequence. Three novel amino acid sequences of ovine IGHA were also predicted, and variation in these sequences might not only affect antigen recognition but also susceptibility to cleavage by bacterial or parasitic proteases. Nucleotide sequence data reported in this paper have been submitted to the NCBI GenBank nucleotide sequence database and have been assigned the accession nos. AY956424–AY956426.     

Within-population variation in cytoplasmic genes affects female life span and aging in Drosophila melanogaster

It has been suggested that mitochondrial DNA (mtDNA) may play an important role in aging. Yet, few empirical studies have tested this hypothesis, partly because the degree of sequence polymorphism in mtDNA is assumed to be low. However, low sequence variation may not necessarily translate into low phenotypic variation. Here, we report an experiment that tests whether there is within-population variation in cytoplasmic genes for female longevity and senescence. To achieve this, we randomly selected 25 "mitochondrial founders" from a single, panmictic population of Drosophila melanogaster and used these founders to generate distinct "mt" lines in which we controlled for the nuclear background by successive backcrossing. Potential confounding effects of cytoplasmically transmitted bacteria were eliminated by tetracycline treatment. The mt lines were then assayed for differences in longevity, Gompertz intercept (frailty), and demographic rate of change in mortality with age (rate-of-senescence) in females. We found significant cytoplasmic effects on all three variables. This provides evidence that genetic variation in cytoplasmic genes, presumably mtDNA, contributes to variation in female mortality and aging.     

家马核基因组中线粒体核内插入序列分析

在各种真核生物核基因组中,存在一些由线粒体基因组转移进入核基因组中的DNA片段,这些被认为是分子化石的片段叫做线粒体核内插入序列(Numt)。由于Numt与真实的线粒体序列高度相似,因此它的存在必然会成为PCR扩增线粒体DNA的不利因素。利用已经公布的家马(Equus caballus)基因组序列(2007年9月公布,GenBank登录号为NC_009144-NC_009175)对家马Numt进行了深入分析,共发现200个可能的Numt,长度范围为29到3727bp,其中有10个的长度大于800bp。分析结果显示由于不存在线粒体控制区域的疑似Numt,因此对基于此区域的群体遗传学研究不会产生影响。本研究还发现在家马进化过程中,第1号和27号染色体更倾向于接受线粒体序列的转移。以上结果将为今后马科动物的研究提供重要的参考信息,有助于避免在线粒体DNA研究中由于Numt污染的存在而得出错误的实验结果。     

Genetic variation of geminiviruses: comparison between sexual and asexual host plant populations

One of the most promising hypotheses for the evolution of sex is that sexual reproduction is advantageous because it increases the rate of adaptive evolution in response to parasites. To investigate this advantage of sex, we compared genetic variation of geminiviruses infecting sexual and asexual populations of Eupatorium (Asteraceae). The infection frequency was 37.5% in the sexual population and 87.8% in the asexual population. The lower infection frequency in the sexual population might be the result of higher genetic diversity of host plants. If geminiviruses have diverged to counter defence systems of genetically variable hosts, genetic diversity of viruses is expected to be higher in sexual host populations than in asexual host populations. To test this expectation, we used single-strand conformation polymorphism (SSCP) analysis to examine genetic diversity of the geminiviruses in a DNA region containing the open-reading frame (ORF) C4 gene, which is known to function as a host range determinant. As predicted, higher genetic diversity of viruses was observed in the sexual population: three SSCP types were found in the asexual population while six types were found in the sexual population. Sequencing of the polymerase chain reaction (PCR) products revealed further genetic diversity. Phylogenetic analysis of the sequences showed that the SSCP types belonged to four different clades. Several SSCP types from the same clade were found in the sexual population, whereas the asexual population included only one SSCP type from each clade. Amino acid replacements of ORF C4 are suggested to be accelerated in the sexual population. This evidence supports the hypothesis that sexual reproduction is advantageous as a defence against epidemic disease.     

Nuclear gene sequences and DNA variation of Cryptomeria japonica samples from the postglacial period

Genomic DNA was extracted from heartwood blocks of six Cryptomeria japonica individuals that had been buried (in an area now covered by rice fields) for about 3600 years. Attempts were made to determine the sequences of five nuclear genes following polymerase chain reaction amplification, using previously obtained C. japonica expressed sequence tag (EST) information. We detected 15 nucleotide substitutions and four insertion/deletions (indels) in a partial GapC gene sequence among 13 individuals of the buried and an extant population, which allowed us to estimate the extent of DNA variation within the buried populations, and the level of genetic differentiation between the buried population and the extant population growing in a neighbouring area. For the entire haplotypes of the GapC region, pi and theta nucleotide diversity estimates were 0.0063 and 0.0010, respectively, when both populations were included, while corresponding figures for the buried population alone were 0.0009 and 0.0017. Estimates of DNA divergence statistics (dXY = 0.0062, dA = 0.0005, FST = 0.0832 and KST = 0.0935) suggest that differentiation between the two populations was not great. However, permutation tests gave FST and KST values rejecting the null hypothesis (that populations were not differentiated) at the 5% and 1% probability levels, respectively. The significant genetic differentiation between the two populations was mainly caused by differences in haplotype diversity. The significant level of haplotype diversity in the extant population compared to the buried population might be the result of gene flow from neighbouring artificial forests. Alternatively, it is possible that we failed to detect all the DNA variation in the buried population because of clonal growth in the buried population.