Homoplasy corrected estimation of genetic similarity from AFLP bands, and the effect of the number of bands on the precision of estimation

AFLP is a DNA fingerprinting technique, resulting in binary band presence–absence patterns, called profiles, with known or unknown band positions. We model AFLP as a sampling procedure of fragments, with lengths sampled from a distribution. Bands represent fragments of specific lengths. We focus on estimation of pairwise genetic similarity, defined as average fraction of common fragments, by AFLP. Usual estimators are Dice (D) or Jaccard coefficients. D overestimates genetic similarity, since identical bands in profile pairs may correspond to different fragments (homoplasy). Another complicating factor is the occurrence of different fragments of equal length within a profile, appearing as a single band, which we call collision. The bias of D increases with larger numbers of bands, and lower genetic similarity. We propose two homoplasy- and collision-corrected estimators of genetic similarity. The first is a modification of D, replacing band counts by estimated fragment counts. The second is a maximum likelihood estimator, only applicable if band positions are available. Properties of the estimators are studied by simulation. Standard errors and confidence intervals for the first are obtained by bootstrapping, and for the second by likelihood theory. The estimators are nearly unbiased, and have for most practical cases smaller standard error than D. The likelihood-based estimator generally gives the highest precision. The relationship between fragment counts and precision is studied using simulation. The usual range of band counts (50–100) appears nearly optimal. The methodology is illustrated using data from a phylogenetic study on lettuce

Diversity and linkage disequilibrium analysis within a selected set of cultivated tomatoes

Within the Dutch genomics initiative the “Centre for Biosystems Genomics” (CBSG) a major research effort is directed at the identification and unraveling of processes and mechanisms affecting fruit quality in tomato. The basis of this fruit quality program was a diverse set of 94 cultivated tomato cultivars, representing a wide spectrum of phenotypes for quality related traits. This paper describes a diversity study performed on these cultivars, using information of 882 AFLP markers, of which 304 markers had a known map position. The AFLP markers were scored as much as possible in a co-dominant fashion. We investigated genome distribution and coverage for the mapped markers and conclude that it proved difficult to arrive at a dense and uniformly distributed coverage of the genome with markers. Mapped markers and unmapped markers were used to investigate population structure. A clear substructure was observed which seemed to coincide with a grouping based on fruit size. Finally, we studied amount and decay of linkage disequilibrium (LD) along the chromosomes. LD was observed over considerable (genetic) distances. We discuss the feasibility of marker-trait association studies and conclude that the amount of genetic variation in our set of cultivars is limited, but that there exists scope for association studies

Codominant scoring of AFLP in association panels

A study on the codominant scoring of AFLP markers in association panels without prior knowledge on genotype probabilities is described. Bands are scored codominantly by fitting normal mixture models to band intensities, illustrating and optimizing existing methodology, which employs the EM-algorithm. We study features that improve the performance of the algorithm, and the unmixing in general, like parameter initialization, restrictions on parameters, data transformation, and outlier removal. Parameter restrictions include equal component variances, equal or nearly equal distances between component means, and mixing probabilities according to Hardy–Weinberg Equilibrium. Histogram visualization of band intensities with superimposed normal densities, and optional classification scores and other grouping information, assists further in the codominant scoring. We find empirical evidence favoring the square root transformation of the band intensity, as was found in segregating populations. Our approach provides posterior genotype probabilities for marker loci. These probabilities can form the basis for association mapping and are more useful than the standard scoring categories A, H, B, C, D. They can also be used to calculate predictors for additive and dominance effects. Diagnostics for data quality of AFLP markers are described: preference for three-component mixture model, good separation between component means, and lack of singletons for the component with highest mean. Software has been developed in R, containing the models for normal mixtures with facilitating features, and visualizations. The methods are applied to an association panel in tomato, comprising 1,175 polymorphic markers on 94 tomato hybrids, as part of a larger study within the Dutch Centre for BioSystems Genomics

A hybrid BAC physical map of potato: a framework for sequencing a heterozygous genome

<p>Abstract</p> <p>Background</p> <p>Potato is the world's third most important food crop, yet cultivar improvement and genomic research in general remain difficult because of the heterozygous and tetraploid nature of its genome. The development of physical map resources that can facilitate genomic analyses in potato has so far been very limited. Here we present the methods of construction and the general statistics of the first two genome-wide BAC physical maps of potato, which were made from the heterozygous diploid clone RH89-039-16 (RH).</p> <p>Results</p> <p>First, a gel electrophoresis-based physical map was made by AFLP fingerprinting of 64478 BAC clones, which were aligned into 4150 contigs with an estimated total length of 1361 Mb. Screening of BAC pools, followed by the KeyMaps <it>in silico </it>anchoring procedure, identified 1725 AFLP markers in the physical map, and 1252 BAC contigs were anchored the ultradense potato genetic map. A second, sequence-tag-based physical map was constructed from 65919 whole genome profiling (WGP) BAC fingerprints and these were aligned into 3601 BAC contigs spanning 1396 Mb. The 39733 BAC clones that overlap between both physical maps provided anchors to 1127 contigs in the WGP physical map, and reduced the number of contigs to around 2800 in each map separately. Both physical maps were 1.64 times longer than the 850 Mb potato genome. Genome heterozygosity and incomplete merging of BAC contigs are two factors that can explain this map inflation. The contig information of both physical maps was united in a single table that describes hybrid potato physical map.</p> <p>Conclusions</p> <p>The AFLP physical map has already been used by the Potato Genome Sequencing Consortium for sequencing 10% of the heterozygous genome of clone RH on a BAC-by-BAC basis. By layering a new WGP physical map on top of the AFLP physical map, a genetically anchored genome-wide framework of 322434 sequence tags has been created. This reference framework can be used for anchoring and ordering of genomic sequences of clone RH (and other potato genotypes), and opens the possibility to finish sequencing of the RH genome in a more efficient way via high throughput next generation approaches.</p

Unexpected Fine-Scale Population Structure in a Broadcast-Spawning Antarctic Marine Mollusc

Several recent empirical studies have challenged the prevailing dogma that broadcast-spawning species exhibit little or no population genetic structure by documenting genetic discontinuities associated with large-scale oceanographic features. However, relatively few studies have explored patterns of genetic differentiation over fine spatial scales. Consequently, we used a hierarchical sampling design to investigate the basis of a weak but significant genetic difference previously reported between Antarctic limpets (Nacella concinna) sampled from Adelaide and Galindez Islands near the base of the Antarctic Peninsula. Three sites within Ryder Bay, Adelaide Island (Rothera Point, Leonie and Anchorage Islands) were each sub-sampled three times, yielding a total of 405 samples that were genotyped at 155 informative Amplified Fragment Length Polymorphisms (AFLPs). Contrary to our initial expectations, limpets from Anchorage Island were found to be subtly, but significantly distinct from those sampled from the other sites. This suggests that local processes may play an important role in generating fine-scale population structure even in species with excellent dispersal capabilities, and highlights the importance of sampling at multiple spatial scales in population genetic surveys

Uncoupling of sexual reproduction from homologous recombination in homozygous Oenothera species

Salient features of the first meiotic division are independent segregation of chromosomes and homologous recombination (HR). In non-sexually reproducing, homozygous species studied to date HR is absent. In this study, we constructed the first linkage maps of homozygous, bivalent-forming Oenothera species and provide evidence that HR was exclusively confined to the chromosome ends of all linkage groups in our population. Co-segregation of complementary DNA-based markers with the major group of AFLP markers indicates that HR has only a minor role in generating genetic diversity of this taxon despite its efficient adaptation capability. Uneven chromosome condensation during meiosis in Oenothera may account for restriction of HR. The use of plants with ancient chromosomal arm arrangement demonstrates that limitation of HR occurred before and independent from species hybridizations and reciprocal translocations of chromosome arms—a phenomenon, which is widespread in the genus. We propose that consecutive loss of HR favored the evolution of reciprocal translocations, beneficial superlinkage groups and ultimately permanent translocation heterozygosity

Population structure and linkage disequilibrium unravelled in tetraploid potato

Association mapping is considered to be an important alternative strategy for the identification of quantitative trait loci (QTL) as compared to traditional QTL mapping. A necessary prerequisite for association analysis to succeed is detailed information regarding hidden population structure and the extent of linkage disequilibrium. A collection of 430 tetraploid potato cultivars, comprising two association panels, has been analysed with 41 AFLP® and 53 SSR primer combinations yielding 3364 AFLP fragments and 653 microsatellite alleles, respectively. Polymorphism information content values and detected number of alleles for the SSRs studied illustrate that commercial potato germplasm seems to be equally diverse as Latin American landrace material. Genome-wide linkage disequilibrium (LD)—reported for the first time for tetraploid potato—was observed up to approximately 5 cM using r2 higher than 0.1 as a criterion for significant LD. Within-group LD, however, stretched on average twice as far when compared to overall LD. A Bayesian approach, a distance-based hierarchical clustering approach as well as principal coordinate analysis were adopted to enquire into population structure. Groups differing in year of market release and market segment (starch, processing industry and fresh consumption) were repeatedly detected. The observation of LD up to 5 cM is promising because the required marker density is not likely to disable the possibilities for association mapping research in tetraploid potato. Population structure appeared to be weak, but strong enough to demand careful modelling of genetic relationships in subsequent marker-trait association analyses. There seems to be a good chance that linkage-based marker-trait associations can be identified at moderate marker densities

Geographic genetic structure of Iberian columbines (gen. Aquilegia)

Southern European columbines (genus Aquilegia)are involved in active processes of diversification, and the Iberian Peninsula offers a privileged observatory to witness the process. Studies on Iberian columbines have provided significant advances on species diversification,but we still lack a complete perspective of the genetic diversification in the Iberian scenario. This work explores how genetic diversity of the genus Aquilegia is geographically structured across the Iberian Peninsula. We used Bayesian clustering methods, principal coordinates analyses, and NJ phenograms to assess the genetic relationships among 285 individuals from 62 locations and detect the main lineages. Genetic diversity of Iberian columbines consists of five geographically structured lineages, corresponding to different Iberian taxa. Differentiation among lineages shows particularly complex admixture patterns at Northeast and highly homogeneous toward Northwest and Southeast. This geographic genetic structure suggests the existence of incomplete lineage sorting and interspecific hybridization as could be expected in recent processes of diversification under the influence of quaternary postglacial migrations. This scenario is consistent with what is proposed by the most recent studies on European and Iberian columbines, which point to geographic isolation and divergent selection by habitat specialization as the main diversification drivers of the Iberian Aquilegia complex