Combining DNA pooling with selective recombinant genotyping for increased efficiency in fine mapping

One of the key steps in positional cloning and marker-aided selection is to identify marker(s) tightly linked to the target gene (i.e., fine mapping). Selective genotyping such as selective recombinant genotyping (SRG) is commonly used in fine mapping for cost-saving. To further decrease genotyping effort and rapidly screen for tightly linked markers, we propose here a combined DNA pooling and SRG strategy. A two-stage pooled genotyping can be used for identifying recombinants between a pair of flanking markers more efficiently, and a joint use of bulked DNA analysis and two-stage pooling can also save cost for genotyping recombinants. The combined DNA pooling and SRG strategy can further be extended to fine mapping for polygenic traits. The numerical results based on hypothetical scenarios and an illustrative application to fine mapping of a mutant gene, called xl(t), in rice suggest that the proposed strategy can remarkably reduce genotyping amount compared with the conventional SRG.     

Sequence-specific amplification polymorphisms (SSAPs): a multi-locus approach for analyzing transposon insertions

Multiple copies of transposable elements, inserted at random around the host genome, can be used as molecular markers. Sequence-specific amplification polymorphisms (SSAPs) amplify the region between a PCR primer site near the end of an element and an adjacent restriction site in the flanking genomic DNA. Each amplified insertion is revealed as a band on a sequencing gel or capillary electrophoresis, and a genomic DNA sample produces a characteristic fingerprint of bands. Here, we explain the SSAP marker development method, which can be performed for any species, with recommendations for experimental parameters in several plant species. It takes about 1-2 weeks to complete the whole SSAP marker procedure.     

Characterization and genetic mapping of a short, highly repeated, interspersed DNA sequence from rice (Oryza sativa L.).

Summary A short, highly repeated, interspersed DNA sequence from rice was characterized using a combination of techniques and genetically mapped to rice chromosomes by restriction fragment length polymorphism (RFLP) analysis. A consensus sequence (GGC)_n, where n varies from 13–16, for the repeated sequence family was deduced from sequence analysis. Southern blot analysis, restriction mapping of repeat element-containing genomic clones, and DNA sequence analysis indicated that the repeated sequence is interspersed in the rice genome, and is heterogeneous and divergent. About 200000 copies are present in the rice genome. Single copy sequences flanking the repeat element were used as RFLP markers to map individual repeat elements. Eleven such repeat elements were mapped to seven different chromosomes. The strategy for characterization of highly dispersed repeated DNA and its uses in genetic mapping, DNA fingerprinting, and evolutionary studies are discussed.     

Bayesian method for gene detection and mapping, using a case and control design and DNA pooling

Association mapping studies aim to determine the genetic basis of a trait. A common experimental design uses a sample of unrelated individuals classified into 2 groups, for example cases and controls. If the trait has a complex genetic basis, consisting of many quantitative trait loci (QTLs), each group needs to be large. Each group must be genotyped at marker loci covering the region of interest; for dense coverage of a large candidate region, or a whole-genome scan, the number of markers will be very large. The total amount of genotyping required for such a study is formidable. A laboratory effort efficient technique called DNA pooling could reduce the amount of genotyping required, but the data generated are less informative and require novel methods for efficient analysis. In this paper, a Bayesian statistical analysis of the classic model of McPeek and Strahs is proposed. In contrast to previous work on this model, I assume that data are collected using DNA pooling, so individual genotypes are not directly observed, and also account for experimental errors. A complete analysis can be performed using analytical integration, a propagation algorithm for a hidden Markov model, and quadrature. The method developed here is both statistically and computationally efficient. It allows simultaneous detection and mapping of a QTL, in a large-scale association mapping study, using data from pooled DNA. The method is shown to perform well on data sets simulated under a realistic coalescent-with-recombination model, and is shown to outperform classical single-point methods. The method is illustrated on data consisting of 27 markers in an 880-kb region around the CYP2D6 gene.     

Old can be new again: HAPPY whole genome sequencing, mapping and assembly

During the last three decades, both genome mapping and sequencing methods have advanced significantly to provide a foundation for scientists to understand genome structures and functions in many species. Generally speaking, genome mapping relies on genome sequencing to provide basic materials, such as DNA probes and markers for their localizations, thus constructing the maps. On the other hand, genome sequencing often requires a high-resolution map as a skeleton for whole genome assembly. However, both genome mapping and sequencing have never come together in one pipeline. After reviewing mapping and next-generation sequencing methods, we would like to share our thoughts with the genome community on how to combine the HAPPY mapping technique with the new-generation sequencing, thus integrating two systems into one pipeline, called HAPPY pipeline. The pipeline starts with preparation of a HAPPY panel, followed by multiple displacement amplification for producing a relatively large quantity of DNA. Instead of conventional marker genotyping, the amplified panel DNA samples are subject to new-generation sequencing with barcode method, which allows us to determine the presence/absence of a sequence contig as a traditional marker in the HAPPY panel. Statistical analysis will then be performed to infer how close or how far away from each other these contigs are within a genome and order the whole genome sequence assembly as well. We believe that such a universal approach will play an important role in genome sequencing, mapping, and assembly of many species; thus advancing genome science and its applications in biomedicine and agriculture.     

An entropy-based measure for QTL mapping using extreme samples of population

Quantitative trait locus (QTL) mapping can be accomplished through the method of selective genotyping, which is based on the differences of frequencies between an upper sample and a lower sample in population. However, amplifying the differences in marker allele frequencies in extreme samples may increase the probability for QTL mapping. Shannon entropy, which is a nonlinear function of allele frequencies, can be used to amplify the differences in marker allele frequencies. In this paper, we present a novel measure for linkage disequilibrium (LD) between a marker and single QTL, that is based on the comparison of the entropy and conditional entropy in a marker in extreme samples of population. This measure of LD between the marker and the trait locus can be used when the marker allele frequencies are known in the extreme samples of a population. We investigate the mapping performance in both analytic and simulation scenarios of a single QTL linked to a single marker. Our results show that the measure has very reasonable performance. In addition, a simulation study is performed on the basis of the haplotype frequencies of 10 SNPs of angiotensin-I converting enzyme (ACE) genes.     

Alumorphs--human DNA polymorphisms detected by polymerase chain reaction using Alu-specific primers

The simultaneous analysis of multiple loci could substantially increase the efficiency of mapping studies. Toward this goal, we used the polymerase chain reaction to amplify multiple DNA fragments originating from dispersed genomic segments that are flanked by Alu repeats. Analysis of different human DNA samples revealed numerous amplification products distinguishable by size, some of which vary between individuals. A family study demonstrated that these polymorphic fragments are inherited in a Mendelian fashion. Because of the ubiquitous distribution of Alu repeats, these markers, called "alumorphs," could be useful for linkage mapping of the human genome. A major advantage of alumorphs is that no prior knowledge of DNA sequence of marker loci is required. This approach may find general application for any genome where interspersed repetitive sequences are found.     

High-throughput RFLP genotyping method for large genomes based on a chemiluminescent detection system

Genetic linkage mapping based on RFLPs is a valuable genomics tool for studying organisms with no genome sequence information. However, the generally used Southern hybridization method based on the radioisotope³²P is not ideal for genotyping large mapping populations. We have overcome limitations of the alternative chemiluminescent detection system and developed a high-throughput RFLP genotyping method suitable for large-scale mapping studies of large genomes. Important elements in our process are PCR labeling of probes, complete removal of post-PCR unincorporated nucleotides via column-based purification methods, use of a 1:4 DIG-[11]-dUTP:dTTP ratio, and using a rocker instead of an orbital shaker during hybridization and post-hybridization processing of membranes. Using this method, we mapped the large genome of the homosporous fern speciesCeratopteris richardii by genotyping a mapping population of 513 doubled haploid line (DHL) progeny of a cross between two completely homozygous parental lines. Our genotyping method can robustly detect sub-picogram quantities of DNA fragments from a large number of samples and can be applied to linkage mapping studies of other organisms with large genomes.     

The Use of a Stringent Selection System Allows the Identification of DNA Elements that Augment Gene Expression

The use of high stringency selection systems often results in the induction of very few recombinant mammalian cell lines, which limits the ability to isolate a cell line with favorable characteristics. The employment of for instance STAR elements in DNA constructs elevates the induced number of colonies and also the protein expression levels in these colonies. Here, we describe a method to systematically identify genomic DNA elements that are able to induce many stably transfected mammalian cell lines. We isolated genomic DNA fragments upstream from the human Rb1 and p73 gene loci and cloned them around an expression cassette that contains a very stringent selection marker. Due to the stringency of the selection marker, hardly any colony survives without flanking DNA elements. We tested fourteen ~3500 bp DNA stretches from the Rb1 and p73 loci. Only two ~3500 bp long DNA fragments, called Rb1E and Rb1F, induced many colonies in the context of the stringent selection system and these colonies displayed high protein expression levels. Functional analysis showed that the Rb1 DNA fragments contained no enhancer, promoter, or STAR activity. Our data show the potential of a methodology to identify novel gene expression augmenting DNA elements in an unbiased manner.     

An alternative inverse PCR (IPCR) method to amplify DNA sequences flanking Tn5 transposon insertions

We have developed an alternative method to amplify DNA sequences flanking Tn5 transposon insertions. This method relies on the identical sequences of inverted terminal repeats, located at the 5' and 3' ends of Tn5, to determine the location and orientation of a transposon insertion within a restriction endonuclease fragment. From this information, PCR primers can be designed to selectively amplify by inverse PCR the DNA flanking one side of the transposon. This method avoids the problem of amplifying or cloning long sequences flanking Tn5. To demonstrate the applicability of this method, we generated Tn5 transposon mutants of Pseudomonas abietaniphila BKME-9 which no longer grew on dehydroabietic acid (DhA). The flanking sequence of one of the mutant (strain BKME-941) which accumulated 7-oxoDhA, was amplified.     

A new family of repetitive, retroposon-like sequences in the genome of the rainbow trout

We have identified a new family of interspersed, moderately repetitive DNA elements, termed the RSg-1 family, in the genome of the rainbow trout. Two of the elements examined here are situated upstream of sequences which code for trout nuclear proteins; a protamine gene (p101) and the clustered histone H4 gene. Sequence comparison of various RSg-1 elements indicated a high degree of nucleotide sequence homology between different members of the family. These repetitive elements exhibit well defined 3' ends which contain poly(A) segments preceded by the consensus polyadenylation signal AATAAA. Sequences flanking the 3' end of the poly(A) tract also conform to a consensus sequence. A similar sequence is also found flanking the 5' terminus of the element in the protamine clone p101, and thus may represent a target-site duplication generated upon insertion of the element into the genome. These characteristics, together with the heterogeneous nature of the 5' ends of the elements, are reminiscent of processed pseudogenes and retroposons such as the mammalian L1 family of interspersed repetitive elements.     

The efficacy of detecting variants with small effects on the Affymetrix 6.0 platform using pooled DNA

Genome-wide genotyping of a cohort using pools rather than individual samples has long been proposed as a cost-saving alternative for performing genome-wide association (GWA) studies. However, successful disease gene mapping using pooled genotyping has thus far been limited to detecting common variants with large effect sizes, which tend not to exist for many complex common diseases or traits. Therefore, for DNA pooling to be a viable strategy for conducting GWA studies, it is important to determine whether commonly used genome-wide SNP array platforms such as the Affymetrix 6.0 array can reliably detect common variants of small effect sizes using pooled DNA. Taking obesity and age at menarche as examples of human complex traits, we assessed the feasibility of genome-wide genotyping of pooled DNA as a single-stage design for phenotype association. By individually genotyping the top associations identified by pooling, we obtained a 14- to 16-fold enrichment of SNPs nominally associated with the phenotype, but we likely missed the top true associations. In addition, we assessed whether genotyping pooled DNA can serve as an inexpensive screen as the second stage of a multi-stage design with a large number of samples by comparing the most cost-effective 3-stage designs with 80% power to detect common variants with genotypic relative risk of 1.1, with and without pooling. Given the current state of the specific technology we employed and the associated genotyping costs, we showed through simulation that a design involving pooling would be 1.07 times more expensive than a design without pooling. Thus, while a significant amount of information exists within the data from pooled DNA, our analysis does not support genotyping pooled DNA as a means to efficiently identify common variants contributing small effects to phenotypes of interest. While our conclusions were based on the specific technology and study design we employed, the approach presented here will be useful for evaluating the utility of other or future genome-wide genotyping platforms in pooled DNA studies.     

Efficiency of selective genotyping for genetic analysis of complex traits and potential applications in crop improvement

Selective genotyping of individuals from the two tails of the phenotypic distribution of a population provides a cost efficient alternative to analysis of the entire population for genetic mapping. Past applications of this approach have been confounded by the small size of entire and tail populations, and insufficient marker density, which result in a high probability of false positives in the detection of quantitative trait loci (QTL). We studied the effect of these factors on the power of QTL detection by simulation of mapping experiments using population sizes of up to 3,000 individuals and tail population sizes of various proportions, and marker densities up to one marker per centiMorgan using complex genetic models including QTL linkage and epistasis. The results indicate that QTL mapping based on selective genotyping is more powerful than simple interval mapping but less powerful than inclusive composite interval mapping. Selective genotyping can be used, along with pooled DNA analysis, to replace genotyping the entire population, for mapping QTL with relatively small effects, as well as linked and interacting QTL. Using diverse germplasm including all available genetics and breeding materials, it is theoretically possible to develop an “all-in-one plate” approach where one 384-well plate could be designed to map almost all agronomic traits of importance in a crop species. Selective genotyping can also be used for genomewide association mapping where it can be integrated with selective phenotyping approaches. We also propose a breeding-to-genetics approach, which starts with identification of extreme phenotypes from segregating populations generated from multiple parental lines and is followed by rapid discovery of individual genes and combinations of gene effects together with simultaneous manipulation in breeding programs.     

Standardization and conversion of marker polymorphism measures

Large scale gene mapping efforts in domestic animals have generated and mapped a large number of genetic markers that are useful for mapping quantitative trait and disease loci and for DNA diagnostic purposes such as parentage testing. Marker polymorphism is an important criterion for selecting genetic markers in planning experiment for mapping quantitative trait loci or for DNA diagnostic purposes. Current formulations of marker polymorphism measures are functions of marker allele frequencies. In this study, two measures of marker polymorphism that are available from gene mapping studies and do not require allele frequencies were proposed and analyzed: the observed polymorphic information content (PIC) and the observed family information content (FIC). The observed FIC was more stable than the observed PIC because the observed FIC is unaffected by the variation in the frequency of heterozygous parents. However, both FIC and PIC are dependent on the gene mapping design. The effective number of alleles is recommended as a tool to standardize marker polymorphism measures so that polymorphism of different markers can be compared on an equal basis, and to obtain a new polymorphism measure (such an exclusion probability) from an existing measure (such as FIC). The usage of the effective number of alleles to standardize FIC, PIC and exclusion probabilities is illustrated using genetic markers in a published linkage map.     

Development and initial characterization of a HAPPY panel for mapping the X. tropicalis genome

HAPPY mapping was designed to pursue the analysis of approximately random HAPloid DNA breakage samples using the PolYmerase chain reaction for mapping genomes. In the present study, we improved the method and integrated two other molecular techniques into the process: whole genome amplification and the Sequenom SNP (single nucleotide polymorphism) genotyping assay in order to facilitate whole genome mapping of X. tropicalis. The former technique amplified enough DNA materials to genotype a large number of markers, while the latter allowed for relatively high throughput marker genotyping with multiplex assays on the HAPPY lines. A total of 58 X. tropicalis genes were genotyped on an initial panel of 383 HAPPY lines, which contributed to formation of a working panel of 146 lines. Further genotyping of 29 markers on the working panel led to construction of a HAPPY map for the X. tropicalis genome. We believe that our improved HAPPY method described in the present study has paved the way for the community to map different genomes with a simple, but powerful approach.     

Calculation of multipoint likelihoods using flanking marker data: a simulation study

The calculation of multipoint likelihoods is computationally challenging, with the exact calculation of multipoint probabilities only possible on small pedigrees with many markers or large pedigrees with few markers. This paper explores the utility of calculating multipoint likelihoods using data on markers flanking a hypothesized position of the trait locus. The calculation of such likelihoods is often feasible, even on large pedigrees with missing data and complex structures. Performance characteristics of the flanking marker procedure are assessed through the calculation of multipoint heterogeneity LOD scores on data simulated for Genetic Analysis Workshop 14 (GAW14). Analysis is restricted to data on the Aipotu population on chromosomes 1, 3, and 4, where chromosomes 1 and 3 are known to contain disease loci. The flanking marker procedure performs well, even when missing data and genotyping errors are introduced.     

Recent amplification of the kangaroo endogenous retrovirus, KERV, limited to the centromere

Mammalian retrotransposons, transposable elements that are processed through an RNA intermediate, are categorized as short interspersed elements (SINEs), long interspersed elements (LINEs), and long terminal repeat (LTR) retroelements, which include endogenous retroviruses. The ability of transposable elements to autonomously amplify led to their initial characterization as selfish or junk DNA; however, it is now known that they may acquire specific cellular functions in a genome and are implicated in host defense mechanisms as well as in genome evolution. Interactions between classes of transposable elements may exert a markedly different and potentially more significant effect on a genome than interactions between members of a single class of transposable elements. We examined the genomic structure and evolution of the kangaroo endogenous retrovirus (KERV) in the marsupial genus Macropus. The complete proviral structure of the kangaroo endogenous retrovirus, phylogenetic relationship among relative retroviruses, and expression of this virus in both Macropus rufogriseus and M. eugenii are presented for the first time. In addition, we show the relative copy number and distribution of the kangaroo endogenous retrovirus in the Macropus genus. Our data indicate that amplification of the kangaroo endogenous retrovirus occurred in a lineage-specific fashion, is restricted to the centromeres, and is not correlated with LINE depletion. Finally, analysis of KERV long terminal repeat sequences using massively parallel sequencing indicates that the recent amplification in M. rufogriseus is likely due to duplications and concerted evolution rather than a high number of independent insertion events.     

Standardization and conversion of marker polymorphism measures

Abstract

Large scale gene mapping efforts in domestic animals have generated and mapped a large number of genetic markers that are useful for mapping quantitative trait and disease loci and for DNA diagnostic purposes such as parentage testing. Marker polymorphism is an important criterion for selecting genetic markers in planning experiment for mapping quantitative trait loci or for DNA diagnostic purposes. Current formulations of marker polymorphism measures are functions of marker allele frequencies. In this study, two measures of marker polymorphism that are available from gene mapping studies and do not require allele frequencies were proposed and analyzed: the observed polymorphic information content (PIC) and the observed family information content (FIC). The observed FIC was more stable than the observed PIC because the observed FIC is unaffected by the variation in the frequency of heterozygous parents. However, both FIC and PIC are dependent on the gene mapping design. The effective number of alleles is recommended as a tool to standardize marker polymorphism measures so that polymorphism of different markers can be compared on an qual basis, and to obtain a new polymorphism measure (such an exclusion probability) from an existing measure (such as FIC). The usage of the effective number of alleles to standardize FIC, PIC and exclusion probabilities is illustrated using genetic markers in a published linkage map.     

Locus-specific chromatin profiling of evolutionarily young transposable elements

Despite a vast expansion in the availability of epigenomic data, our knowledge of the chromatin landscape at interspersed repeats remains highly limited by difficulties in mapping short-read sequencing data to these regions. In particular, little is known about the locus-specific regulation of evolutionarily young transposable elements (TEs), which have been implicated in genome stability, gene regulation and innate immunity in a variety of developmental and disease contexts. Here we propose an approach for generating locus-specific protein–DNA binding profiles at interspersed repeats, which leverages information on the spatial proximity between repetitive and non-repetitive genomic regions. We demonstrate that the combination of HiChIP and a newly developed mapping tool (PAtChER) yields accurate protein enrichment profiles at individual repetitive loci. Using this approach, we reveal previously unappreciated variation in the epigenetic profiles of young TE loci in mouse and human cells. Insights gained using our method will be invaluable for dissecting the molecular determinants of TE regulation and their impact on the genome.     

Characterization of equine microsatellites and microsatellite-linked repetitive elements (eMLREs) by efficient cloning and genotyping methods.

We performed efficient cloning and genotyping methods for isolation of a large number of polymorphic microsatellites. The methods contain the time-efficient cloning method of constructing microsatellite-enriched libraries and the economic genotyping method of fluorescent labeling of PCR products. Eighty novel equine microsatellites cloned were efficiently isolated from the enrichment library and analyzed for genotype polymorphism. Of these, 72 microsatellites were analyzed with a good resolution. The average heterozygosity of all loci was 0.52, and the number of alleles ranged from one to 9 with an average of 4.5 alleles. The other eight loci showed multiple bands of PCR products, suggesting the occurrence of microsatellites in a repetitive element, in which the number of microsatellite repeats varies among different members of the repetitive element. We found five homologous groups at flanking regions in comparison with the flanking regions of microsatellites from DNA databases. One of them showed homology to equine repetitive element-2. In the other four homologous groups, the two groups were named equine microsatellite-linked repetitive element-1 (eMLRE-1) and equine microsatellite-linked repetitive element-2 (eMLRE-2) as novel equine repetitive elements identified from equine genome. These data should help the analysis of equine DNA sequences and the design of equine genome markers.