Prediction of total genetic value using genome-wide dense marker maps

Recent advances in molecular genetic techniques will make dense marker maps available and genotyping many individuals for these markers feasible. Here we attempted to estimate the effects of approximately 50,000 marker haplotypes simultaneously from a limited number of phenotypic records. A genome of 1000 cM was simulated with a marker spacing of 1 cM. The markers surrounding every 1-cM region were combined into marker haplotypes. Due to finite population size N(e) = 100, the marker haplotypes were in linkage disequilibrium with the QTL located between the markers. Using least squares, all haplotype effects could not be estimated simultaneously. When only the biggest effects were included, they were overestimated and the accuracy of predicting genetic values of the offspring of the recorded animals was only 0.32. Best linear unbiased prediction of haplotype effects assumed equal variances associated to each 1-cM chromosomal segment, which yielded an accuracy of 0.73, although this assumption was far from true. Bayesian methods that assumed a prior distribution of the variance associated with each chromosome segment increased this accuracy to 0.85, even when the prior was not correct. It was concluded that selection on genetic values predicted from markers could substantially increase the rate of genetic gain in animals and plants, especially if combined with reproductive techniques to shorten the generation interval.     

Constructing dense genetic linkage maps

This paper describes a novel combination of techniques for the construction of dense genetic linkage maps. The construction of such maps is hampered by the occurrence of even small proportions of typing errors. Simulated annealing is used to obtain the best map according to the optimality criterion: the likelihood or the total number of recombination events. Spatial sampling of markers is used to obtain a framework map. The construction of a framework map is essential if the steps used for simulated annealing are required to be simple. For missing-data imputation the Gibbs sampler is used. Map construction using simulated annealing and missing-data imputation are used in an iterative way. In order to obtain some measure of precision of the genetic linkage map obtained, the Metropolis-Hastings algorithm is used to obtain posterior intervals for the positions of markers. The process of map construction is embedded in a framework of pre-mapping and post-mapping diagnostics. The techniques described are illustrated using a practical application. Received: 1 June 2000 / Accepted: 21 September 2000     

Natural selection affects multiple aspects of genetic variation at putatively neutral sites across the human genome

A major question in evolutionary biology is how natural selection has shaped patterns of genetic variation across the human genome. Previous work has documented a reduction in genetic diversity in regions of the genome with low recombination rates. However, it is unclear whether other summaries of genetic variation, like allele frequencies, are also correlated with recombination rate and whether these correlations can be explained solely by negative selection against deleterious mutations or whether positive selection acting on favorable alleles is also required. Here we attempt to address these questions by analyzing three different genome-wide resequencing datasets from European individuals. We document several significant correlations between different genomic features. In particular, we find that average minor allele frequency and diversity are reduced in regions of low recombination and that human diversity, human-chimp divergence, and average minor allele frequency are reduced near genes. Population genetic simulations show that either positive natural selection acting on favorable mutations or negative natural selection acting against deleterious mutations can explain these correlations. However, models with strong positive selection on nonsynonymous mutations and little negative selection predict a stronger negative correlation between neutral diversity and nonsynonymous divergence than observed in the actual data, supporting the importance of negative, rather than positive, selection throughout the genome. Further, we show that the widespread presence of weakly deleterious alleles, rather than a small number of strongly positively selected mutations, is responsible for the correlation between neutral genetic diversity and recombination rate. This work suggests that natural selection has affected multiple aspects of linked neutral variation throughout the human genome and that positive selection is not required to explain these observations.     

Finding associations in dense genetic maps: a genetic algorithm approach

Large-scale association studies hold promise for discovering the genetic basis of common human disease. These studies will consist of a large number of individuals, as well as large number of genetic markers, such as single nucleotide polymorphisms (SNPs). The potential size of the data and the resulting model space require the development of efficient methodology to unravel associations between phenotypes and SNPs in dense genetic maps. Our approach uses a genetic algorithm (GA) to construct logic trees consisting of Boolean expressions involving strings or blocks of SNPs. These blocks or nodes of the logic trees consist of SNPs in high linkage disequilibrium (LD), that is, SNPs that are highly correlated with each other due to evolutionary processes. At each generation of our GA, a population of logic tree models is modified using selection, cross-over and mutation moves. Logic trees are selected for the next generation using a fitness function based on the marginal likelihood in a Bayesian regression frame-work. Mutation and cross-over moves use LD measures to pro pose changes to the trees, and facilitate the movement through the model space. We demonstrate our method and the flexibility of logic tree structure with variable nodal lengths on simulated data from a coalescent model, as well as data from a candidate gene study of quantitative genetic variation.     

Estimating quantitative genetic parameters using sibships reconstructed from marker data

Previous techniques for estimating quantitative genetic parameters, such as heritability in populations where exact relationships are unknown but are instead inferred from marker genotypes, have used data from individuals on a pairwise level only. At this level, families are weighted according to the number of pairs within which each family appears, hence by size rather than information content, and information from multiple relationships is lost. Estimates of parameters are therefore not the most efficient achievable. Here, Markov chain Monte Carlo techniques have been used to partition the population into complete sibships, including, if known, prior knowledge of the distribution of family sizes. These pedigrees have then been used with restricted maximum likelihood under an animal model to estimate quantitative genetic parameters. Simulations to compare the properties of parameter estimates with those of existing techniques indicate that the use of sibship reconstruction is superior to earlier methods, having lower mean square errors and showing nonsignificant downward bias. In addition, sibship reconstruction allows the estimation of population allele frequencies that account for the relationships within the sample, so prior knowledge of allele frequencies need not be assumed. Extensions to these techniques allow reconstruction of half sibships when some or all of the maternal genotypes are known.     

Social heterosis and the maintenance of genetic diversity at the genome level

Social heterosis is when individuals in groups or neighbourhoods receive a mutualistic benefit from across‐individual genetic diversity. Although it can be a viable evolutionary mechanism to maintain allelic diversity at a given locus, its efficacy at maintaining genome‐wide diversity is in question when multiple loci are being simultaneously selected. Therefore, we modelled social heterosis in a population of haploid genomes of two‐ or three‐linked loci. With such linkages, social heterosis decreases gametic diversity, but maintains allelic diversity. Genomes tend to survive as complimentary pairs, with alternate alleles at each locus (e.g. the pair AbC and aBc). The outcomes of selection appear similar to fitness epistasis but are novel in the sense that phenotypic interactions occur across rather than within individuals. The model’s results strongly suggest that strong linkage across gene loci actually increases the probability that social heterosis maintains significant genetic diversity at the level of the genome.     

Preliminary interspecific genetic maps of the populus genome constructed from RAPD markers.

We have constructed RAPD-based linkage maps for an interspecific cross between two species of the genus Populus (P. adenopoda and P. alba), based on a double pseudo-test-cross strategy. Of a total of 360 polymorphic fragments scored, 290 showed a test-cross configuration, corresponding to DNA polymorphisms heterozygous in one parent and null in the other. In the female parent, P. adenopoda, 82 markers were grouped in 19 different linkage groups (553 cM), whereas in the male parent P. alba, 197 markers established a much more complete framework map with an observed genome length of 2300 cM covering 87% of the total P. alba genome. The larger number of test-cross markers detected for the P. alba parent than for the P. adenopoda parent might be due to a higher level of heterozygosity in the former than in the latter. In this study, we detected only a small percentage (2%) of the intercross dominant markers heterozygous in both parents and segregating 3:1 in the progeny. The further focus in this mapping study should be on the identification of more intercross markers, to align the two parent-specific maps into a consensus map for mapping important genes causing species differentiation during long evolutionary divergences.     

Analysis of genetic diversity in Ganoderma population with a novel molecular marker SRAP

Genetic marker technology designed to detect naturally occurring polymorphisms at the DNA level had become an invaluable and revolutionizing tool for both applied and basic studies of fungi. To eliminate the confusion on the taxonomy of Ganoderma strains, in this study, a collection of 31 accessions representative of morphotypes and some unclassified types was used for analyzing molecular diversity using a novel molecular marker sequence-related amplified polymorphism (SRAP). This collection included commercial cultivars and wild varieties that represented the great diversification of types from different countries and regions. The experimental results showed that 50 out of 95 combinations of primers turned out to be polymorphic, and 85 polymorphism bands were obtained using six combinations. Based on the appearances of markers, the genetic similarity coefficients were calculated, and genetic variations were observed (0∼1) among the 31 different Ganoderma strains. The group of Ganoderma lucidum showed significant differences from the group of Ganoderma sinense. Moreover, G. lucidum in China was also different from G. lucidum in Yugoslavia. At the same time, cluster analysis successfully categorized these 31 Ganoderma strains into five groups. These results revealed the genetic diversity of Ganoderma strains and their correlation with geographic environments. It also suggested SRAP marker could be used in the taxonomic analysis of fungi. To our knowledge, this is the first application of SRAP marker on the systematics of Ganoderma strains within basidiomycetes.     

Castor bean organelle genome sequencing and worldwide genetic diversity analysis

Castor bean is an important oil-producing plant in the Euphorbiaceae family. Its high-quality oil contains up to 90% of the unusual fatty acid ricinoleate, which has many industrial and medical applications. Castor bean seeds also contain ricin, a highly toxic Type 2 ribosome-inactivating protein, which has gained relevance in recent years due to biosafety concerns. In order to gain knowledge on global genetic diversity in castor bean and to ultimately help the development of breeding and forensic tools, we carried out an extensive chloroplast sequence diversity analysis. Taking advantage of the recently published genome sequence of castor bean, we assembled the chloroplast and mitochondrion genomes extracting selected reads from the available whole genome shotgun reads. Using the chloroplast reference genome we used the methylation filtration technique to readily obtain draft genome sequences of 7 geographically and genetically diverse castor bean accessions. These sequence data were used to identify single nucleotide polymorphism markers and phylogenetic analysis resulted in the identification of two major clades that were not apparent in previous population genetic studies using genetic markers derived from nuclear DNA. Two distinct sub-clades could be defined within each major clade and large-scale genotyping of castor bean populations worldwide confirmed previously observed low levels of genetic diversity and showed a broad geographic distribution of each sub-clade.     

Estimating optimum number of marker loci for genetic analyses in Cynodon accessions

This study was carried out to determine genome sampling size of Cynodon accessions collected from Turkey for estimating genetic relatioships. One hundred and eighty-two Cynodon accessions collected in Turkey from an area south of the Taurus Mountains along the Mediterranean cost and ten known genotypes were genotyped using sequence related amplified polymorphism (SRAP) and the samples include diploids, triploids, tetraploids, pentaploids, and hexaploids. In this study, 182-markers were available from a previous work. We created 7 different data file, which include 5, 10, 20, 40, 60, 80, and 100 markers randomly selected from a 182-marker data file. First, similarity matrices was produced based on Dice’ and, alternatively, simple matching’ similarity coefficients for each data files. Then, a two-way mantel test was used to estimate Mantel correlation coefficients between each pair of similarity matrices. Mantel correlation coefficients between 5 and 10, 10 and 20, 20 and 40, 40 and 60, 60 and 80, and 80 and 100 marker-based simple matching’ similarity matrices were 0.718, 0.802, 0.826, 0.929, 0.977, and 0.975, respectively. Correlation estimates for the matrices based on Dice’s procedure for the same data sets were 0.863, 0.916, 0.930, 0.966, 0.965, and 0.984. The lower number of markers (10 or 20) was required with Dice’s procedure to achieve the higher correlation values (r > 0.9), and about random genome sampling size of 10 or 20 markers appeared to be sufficent in Cynodon accessions with a series of ploidy. This approach may apply to other plant species with ploidy series.     

Estimating genetic diversity of Arabidopsis thaliana ecotypes with amplified fragment length polymorphisms (AFLP)

The extensive natural variation of Arabidopsis thaliana ecotypes is being increasingly exploited as a source of variants of genes which control (agronomically) important traits. We have subjected 19 different Arabidopsis thaliana ecotypes to an analysis using the anplified fragment length polymorphism (AFLP) technique in order to estimate their genetic diversity. The genetic diversity was estimated applying the method of Nei and Li (1979) and a modified version of it and using 471 informative polymorphisms. The data obtained revealed that within this small set of ecotypes a group of three ecotypes and a further single ecotype exhibit considerable genetic diversity in comparison to the others. These ecotypes clustered at positions significantly separated from the bulk of the ecotypes in the generated similarity plots. The analysis demonstrated the usefulness of the AFLP method for determinating intraspecies genetic diversity as exemplified with Arabidopsis thaliana ecotypes. Results are discussed and compared with data obtained with other methods. Received: 18 June 1999 / Accepted: 28 July 1999     

Cloning of a polymorphic canine genetic marker which maps to human chromosome 9

We describe the cloning of a novel canine polymorphic genetic marker which maps to human chromosome 9. The sequence is 2092 bp, 59% GC rich, and contains three GC boxes. Chemilumin-escent probing of zooblots showed evolutionary conservation. Dogs have three Bam HI alleles: 2.3 kb, 2.1 kb and 1.7 kb. Allele frequencies in 17 unrelated dogs representing 13 breeds are presented. Polymorphism for the 1.7-kb allele in beagles is common. The 2.1-kb allele is probably the ancestral allele since it is the most common and is also noted in the Cape hunting dog. Interestingly, in more than 50 dogs tested to date, the 2.3-kb allele has been found only in miniature and giant schnauzers. This points to a common origin for these two breeds.     

EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome

In a panel of seven genotypes, 437 expressed sequence tag (EST)-derived DNA fragments were sequenced. Single nucleotide polymorphisms (SNPs) that were polymorphic between the parents of three mapping populations were mapped by heteroduplex analysis and a genome-wide consensus map comprising 216 EST-derived SNPs and 4 InDel (insertion/deletion) markers was constructed. The average frequency of SNPs amounted to 1/130 bp and 1/107.8 bp for a set of randomly selected and a set of mapped ESTs, respectively. The calculated nucleotide diversities (π) ranged from 0 to 40.0 × 10⁻³ (average 3.1 × 10⁻³) and 0.52 × 10⁻³ to 39.51 × 10^–3 (average 4.37 × 10⁻³) for random and mapped ESTs, respectively. The polymorphism information content value for mapped SNPs ranged from 0.24 to 0.50 with an average of 0.34. As expected, combination of SNPs present in an amplicon (haplotype) exhibited a higher information content ranging from 0.24 to 0.85 with an average of 0.50. Cleaved amplified polymorphic sequence assays (including InDels) were designed for a total of 87 (39.5%) SNP markers. The high abundance of SNPs in the barley genome provides avenues for the systematic development of saturated genetic maps and their integration with physical maps. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Both R. Kota and R.K. Varshney contributed equally to this work.     

Toward integration of comparative genetic, physical, diversity, and cytomolecular maps for grasses and grains, using the sorghum genome as a foundation

The small genome of sorghum (Sorghum bicolor L. Moench.) provides an important template for study of closely related large-genome crops such as maize (Zea mays) and sugarcane (Saccharum spp.), and is a logical complement to distantly related rice (Oryza sativa) as a "grass genome model." Using a high-density RFLP map as a framework, a robust physical map of sorghum is being assembled by integrating hybridization and fingerprint data with comparative data from related taxa such as rice and using new methods to resolve genomic duplications into locus-specific groups. By taking advantage of allelic variation revealed by heterologous probes, the positions of corresponding loci on the wheat (Triticum aestivum), rice, maize, sugarcane, and Arabidopsis genomes are being interpolated on the sorghum physical map. Bacterial artificial chromosomes for the small genome of rice are shown to close several gaps in the sorghum contigs; the emerging rice physical map and assembled sequence will further accelerate progress. An important motivation for developing genomic tools is to relate molecular level variation to phenotypic diversity. "Diversity maps," which depict the levels and patterns of variation in different gene pools, shed light on relationships of allelic diversity with chromosome organization, and suggest possible locations of genomic regions that are under selection due to major gene effects (some of which may be revealed by quantitative trait locus mapping). Both physical maps and diversity maps suggest interesting features that may be integrally related to the chromosomal context of DNA-progress in cytology promises to provide a means to elucidate such relationships. We seek to provide a detailed picture of the structure, function, and evolution of the genome of sorghum and its relatives, together with molecular tools such as locus-specific sequence-tagged site DNA markers and bacterial artificial chromosome contigs that will have enduring value for many aspects of genome analysis.     

Development of a high-resolution molecular marker for tracking Rhizosolenia setigera genetic diversity

Journal of Applied Phycology - Rhizosolenia setigera is cosmopolitan diatom with broad distribution in many ocean regions with diverse environmental features, suggesting that R. setigera may have...     

基于RAPD标记的南苍术居群遗传多样性分析

采用RAPD标记技术对分布于江苏小九华山、小汤山和湖山,安徽金寨和芜湖以及湖北保康和英山的7个南苍术〔Atractylodes lancea(Thunb.)DC.〕野生居群的28个单株基因组总DNA进行PCR扩增,在此基础上分析居群的遗传多样性及遗传分化,并采用聚类分析法对居群的遗传关系进行分析。结果表明:用18条RAPD引物共扩增出193条带,其中多态性条带111条,多态性条带百分率(PPB)为57.51%;平均每条引物扩增出10.72条带,其中多态性条带6.17条。从省级水平看,安徽居群的PPB、有效等位基因数(Ne)、Nei’s基因多样性指数(H)和Shannon信息指数(I)均最低,而湖北居群的Ne、H和I均最高,但江苏居群的PPB最高;从居群水平看,湖北保康居群的PPB、Ne、H和I均最高,而安徽金寨居群均最低。7个居群的基因分化系数和基因流分别为0.206 5和1.921 5,说明7个居群总遗传变异的20.65%存在于居群间、79.35%存在于居群内。7个居群间的遗传距离为0.150 7～0.252 1,其中,安徽金寨和芜湖居群间最小(0.150 7),江苏湖山和安徽芜湖居群间最大(0.252 1)。基于遗传距离的聚类分析结果表明:7个居群可分为2组,湖北保康居群单独成组,其他6个居群聚为另一组;来自同一居群的单株均聚在一起。研究结果提示:南苍术居群间的遗传多样性较低,居群间无明显的遗传分化。     

Comparative analysis of the Oenococcus oeni pan genome reveals genetic diversity in industrially-relevant pathways

ABSTRACT: BACKGROUND: Oenococcus oeni, a member of the lactic acid bacteria, is one of a limited number of microorganisms that not only survive, but actively proliferate in wine. It is is also unusual as, unlike the majority of bacteria present in wine, it is beneficial to wine quality rather than causing spoilage. These benefits are realised primarily through catalysing malolactic fermentation, but also through imparting other positive sensory properties. However, many of these industrially-important secondary attributes have been shown to be strain-dependent and their genetic basis it yet to be determined. RESULTS: In order to investigate the scale and scope of genetic variation in O. oeni, we have performed whole-genome sequencing on eleven strains of this bacterium, bringing the total number of strains for which genome sequences are available to fourteen. While any single strain of O. oeni was shown to contain around 1800 protein-coding genes, in-depth comparative annotation based on genomic synteny and protein orthology identified over 2800 orthologous open reading frames that comprise the pan genome of this species, and less than 1200 genes that make up the conserved genomic core present in all of the strains. The expansion of the pan genome relative to the coding potential of individual strains was shown to be due to the varied presence and location of multiple distinct bacteriophage sequences and also in various metabolic functions with potential impacts on the industrial performance of this species, including cell wall exopolysaccharide biosynthesis, sugar transport and utilisation and amino acid biosynthesis. CONCLUSIONS: By providing a large cohort of sequenced strains, this study provides a broad insight into the genetic variation present within O. oeni. This data is vital to understanding and harnessing the phenotypic variation present in this economically-important species.     

Isochore chromosome maps of the human genome

The human genome is a mosaic of isochores, which are long DNA segments (300 kbp) relatively homogeneous in G+C. Human isochores were first identified by density-gradient ultracentrifugation of bulk DNA, and differ in important features, e.g. genes are found predominantly in the GC-richest isochores. Here, we use a reliable segmentation method to partition the longest contigs in the human genome draft sequence into long homogeneous genome regions (LHGRs), thereby revealing the isochore structure of the human genome. The advantages of the isochore maps presented here are: (1) sequence heterogeneities at different scales are shown in the same plot; (2) pair-wise compositional differences between adjacent regions are all statistically significant; (3) isochore boundaries are accurately defined to single base pair resolution; and (4) both gradual and abrupt isochore boundaries are simultaneously revealed. Taking advantage of the wide sample of genome sequence analyzed, we investigate the correspondence between LHGRs and true human isochores revealed through DNA centrifugation. LHGRs show many of the typical isochore features, mainly size distribution, G+C range, and proportions of the isochore classes. The relative density of genes, Alu and long interspersed nuclear element repeats and the different types of single nucleotide polymorphisms on LHGRs also coincide with expectations in true isochores. Potential applications of isochore maps range from the improvement of gene-finding algorithms to the prediction of linkage disequilibrium levels in association studies between marker genes and complex traits. The coordinates for the LHGRs identified in all the contigs longer than 2 Mb in the human genome sequence are available at the online resource on isochore mapping: http://bioinfo2.ugr.es/isochores.     

A functional genetic marker for apple red skin coloration across different environments

The red skin color desired by most apple consumers is not easy to achieve in warm climates, as the expression of MYB10, which regulates red pigmentation in apple, is influenced negatively by high temperatures. We describe the development and validation of a genetic marker for red skin coloration that effectively predicts color in a warm summer environment in Spain, as well as more temperate climates in New Zealand and Italy. Following the determination of a major-effect quantitative trait locus (QTL) controlling red skin coloration on linkage group (LG)9, using four segregating populations grown in New Zealand, and screened using the IRSC apple 8-K single-nucleotide polymorphism (SNP) array, the most significant SNP marker (ss475879531) was transformed into a marker suitable for use in a real-time PCR assay. This marker was validated using five apple seedling populations growing in a warm summer environment in Spain, demonstrating that the marker system efficiently predicts red skin coloration and can be used for marker assisted selection, even under conditions considered adverse for skin color development.