Assessing probability of ancestry using simple sequence repeat profiles: applications to maize hybrids and inbreds

Determination of parentage is fundamental to the study of biology and to applications such as the identification of pedigrees. Limitations to studies of parentage have stemmed from the use of an insufficient number of hypervariable loci and mismatches of alleles that can be caused by mutation or by laboratory error and that can generate false exclusions. Furthermore, most studies of parentage have been limited to comparisons of small numbers of specific parent-progeny triplets thereby precluding large-scale surveys of candidates where there may be no prior knowledge of parentage. We present an algorithm that can determine probability of parentage in circumstances where there is no prior knowledge of pedigree and that is robust in the face of missing data or mistyped data. We present data from 54 maize hybrids and 586 maize inbreds that were profiled using 195 SSR loci including simulations of additional levels of missing and mistyped data to demonstrate the utility and flexibility of this algorithm.     

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct "Seq-to-SSR" approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus, and Clark's Nutcracker, Nucifraga columbiana), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample--a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.     

Evaluation of Polymorphism at Microsatellite Loci of Spring Durum Wheat (Triticum durum Desf.) Varieties and the Use of SSR-Based Analysis in Phylogenetic Studies

Polymorphism at 28 SSR loci was analyzed and described in 45 cultivars of spring durum wheat created in the former USSR and Russia during the last 80 years. Each cultivar was shown to have a unique allele combination. This allows SSR markers to be used to identify durum wheat varieties. Meanwhile, these markers can hardly be used to detect phylogenetic relationships among varieties and to specify their pedigrees, because genetic distances calculated on the basis of these markers do not correlate with the distance calculated by coefficient of parentage.     

Evaluation of polymorphism at microsatellite loci of spring durum wheat (Triticum durum Desf.) varieties and the use of SSR-based analysis in phylogenetic studies

Polymorphism at 28 SSR loci was analyzed and described in 45 cultivars of spring durum wheat created in the former USSR and Russia during the last 80 years. Each cultivar was shown to have a unique allele combination. This allows SSR markers to be used to identify durum wheat varieties. Meanwhile, these markers can hardly be used to detect phylogenetic relationships among varieties and to specify their pedigrees, because genetic distances calculated on the basis of these markers do not correlate with the distance calculated by coefficient of parentage.     

Statistical confidence in parentage analysis with incomplete sampling: how many loci and offspring are needed?

We have recently presented models to estimate parentage in breeding systems with multiple mating and incomplete sampling of the candidate parents. Here we provide formulas to calculate the statistical confidence and the optimal trade-off between the number of loci and offspring. These calculations allow an understanding of the statistical significance of the parentage estimates as well as the appropriate sampling regime required to obtain a desired level of confidence. We show that the trade-off generally depends on the parentage of the putative parents. When parentage is low, sampling effort should concentrate on increasing the number of loci. Otherwise, there are similar benefits from increasing the number of loci or offspring. We demonstrate these methods using genetic data from a nest of the bluegill sunfish (Lepomis macrochirus).     

辣椒优良自交系间遗传差异的分子分析

作物自交系间遗传差异的分析与评价是杂种优势育种和杂交育种成功的基础。鲜食尖椒类(Capsicum an-nuumvar.longum)品种是我国辣椒生产的主要品种类型之一。针对我国鲜食尖椒的育种目标,以国内外10份尖椒优良自交系为材料,利用相关序列扩增多态性(SRAP)和简单序列重复(SSR)标记技术对其进行了遗传差异分析。结果显示:SRAP技术具有较高的位点和多态性检测能力,平均每次检测的位点数和多态性位点数分别为34个和10个,是SSR的10倍和5倍;辣椒自交系间基于SRAP标记的遗传距离和基于SSR标记的遗传距离之间的相关程度较低(r=0.144);基于SRAP标记和SSR标记联合数据计算的遗传距离,10个尖椒自交系被分为3大类,这种分类结果与辣椒杂种优势育种实践相一致。本研究结果表明,SRAP具有较高的遗传分析效力;基于不同分子标记遗传分析结果的差异与标记间共享位点的多少有关;10个尖椒自交系的分类结果可用于指导育种实践。     

Selection of high heterozygosity popcorn varieties in Brazil based on SSR markers

We analyzed genetic structure and diversity among eight populations of popcorn, using SSR loci as genetic markers. Our objectives were to select SSR loci that could be used to estimate genetic diversity within popcorn populations, and to analyze the genetic structure of promising populations with high levels of heterozygosity that could be used in breeding programs. Fifty-seven alleles (3.7 alleles per locus) were detected; the highest effective number of alleles (4.21) and the highest gene diversity (0.763) were found for the Umc2226 locus. A very high level of population differentiation was found (F(ST) = 0.3664), with F(ST) for each locus ranging from 0.1029 (Umc1664) to 0.6010 (Umc2350). This analysis allowed us to identify SSR loci with high levels of heterozygosity and heterozygous varieties, which could be selected for production of inbred lines and for developing new cultivars.     

Using molecular sizes of simple sequence repeats vs. discrete binned data in assessing probability of ancestry: application to maize hybrids

Most inferential methods for profiling genotypes based upon the use of DNA fragments use molecular-size data transcribed into discrete bins, which are intervals of DNA fragment sizes. Categorizing into bins is labor intensive with inevitable arbitrariness that may vary between laboratories. We describe and evaluate an algorithm for determining probabilities of parentage based on raw molecular-size data without establishing bins. We determine the standard deviation of DNA fragment size and assess the association of standard deviation with fragment size. We consider a pool of potential ancestors for an index line that is a hybrid with unknown pedigree. We evaluate the identification of inbred parents of maize hybrids with simple sequence repeat data in the form of actual molecular sizes received from two laboratories. We find the standard deviation to be essentially constant over the molecular weight. We compare these results with those of parallel analyses based on these same data that had been transcribed into discrete bins by the respective laboratories. The conclusions were quite similar in the two cases, with excellent performance using either binned or molecular-size data. We demonstrate the algorithm's utility and robustness through simulations of levels of missing and misscored molecular-size data.     

APIS: An auto‐adaptive parentage inference software that tolerates missing parents

In the context of parentage assignment using genomic markers, key issues are genotyping errors and an absence of parent genotypes because of sampling, traceability or genotyping problems. Most likelihood‐based parentage assignment software programs require a priori estimates of genotyping errors and the proportion of missing parents to set up meaningful assignment decision rules. We present here the R package APIS, which can assign offspring to their parents without any prior information other than the offspring and parental genotypes, and a user‐defined, acceptable error rate among assigned offspring. Assignment decision rules use the distributions of average Mendelian transmission probabilities, which enable estimates of the proportion of offspring with missing parental genotypes. APIS has been compared to other software (CERVUS, VITASSIGN), on a real European seabass (Dicentrarchus labrax) single nucleotide polymorphism data set. The type I error rate (false positives) was lower with APIS than with other software, especially when parental genotypes were missing, but the true positive rate was also lower, except when the theoretical exclusion power reached 0.99999. In general, APIS provided assignments that satisfied the user‐set acceptable error rate of 1% or 5%, even when tested on simulated data with high genotyping error rates (1% or 3%) and up to 50% missing sires. Because it uses the observed distribution of Mendelian transmission probabilities, APIS is best suited to assigning parentage when numerous offspring (>200) are genotyped. We have demonstrated that APIS is an easy‐to‐use and reliable software for parentage assignment, even when up to 50% of sires are missing.     

Development of SSR markers and construction of a linkage map in jute

Jute is an important natural fibre crop, which is only second to cotton in its importance at the global level. It is mostly grown in Indian subcontinent and has been recently used for the development of genomics resources.We recently initiated a programme to develop simple sequence repeat markers and reported a set of 2469 SSR that were developed using four SSR-enriched libraries (Mir et al. 2009). In this communication, we report an additional set of 607 novel SSR in 393 SSR containing sequences. However, primers could be designed for only 417 potentially useful SSR. Polymorphism survey was carried out for 374 primer pairs using two parental genotypes (JRO 524 and PPO4) of a mapping population developed for fibre fineness; only 66 SSR were polymorphic. Owing to a low level of polymorphism between the parental genotypes and a high degree of segregation distortion in recombinant inbred lines, genotypic data of only 53 polymorphic SSR on the mapping population consisting of 120 RIL could be used for the construction of a linkage map; 36 SSR loci were mapped on six linkage groups that covered a total genetic distance of 784.3 cM. Hopefully, this map will be enriched with more SSR loci in future and will prove useful for identification of quantitative trait loci/genes for molecular breeding involving improvement of fibre fineness and other related traits in jute.     

Mining and validation of pyrosequenced simple sequence repeats (SSRs) from American cranberry (<Emphasis Type="Italic">Vaccinium macrocarpon</Emphasis> Ait.)

The American cranberry (Vaccinium macrocarpon Ait.) is a major commercial fruit crop in North America, but limited genetic resources have been developed for the species. Furthermore, the paucity of codominant DNA markers has hampered the advance of genetic research in cranberry and the Ericaceae family in general. Therefore, we used Roche 454 sequencing technology to perform low-coverage whole genome shotgun sequencing of the cranberry cultivar ‘HyRed’. After de novo assembly, the obtained sequence covered 266.3 Mb of the estimated 540–590 Mb in cranberry genome. A total of 107,244 SSR loci were detected with an overall density across the genome of 403 SSR/Mb. The AG repeat was the most frequent motif in cranberry accounting for 35% of all SSRs and together with AAG and AAAT accounted for 46% of all loci discovered. To validate the SSR loci, we designed 96 primer-pairs using contig sequence data containing perfect SSR repeats, and studied the genetic diversity of 25 cranberry genotypes. We identified 48 polymorphic SSR loci with 2–15 alleles per locus for a total of 323 alleles in the 25 cranberry genotypes. Genetic clustering by principal coordinates and genetic structure analyzes confirmed the heterogeneous nature of cranberries. The parentage composition of several hybrid cultivars was evident from the structure analyzes. Whole genome shotgun 454 sequencing was a cost-effective and efficient way to identify numerous SSR repeats in the cranberry sequence for marker development.     

利用SSR标记鉴定西瓜杂交种纯度的研究

以2个西瓜杂交品种(系)的种子黑公子和04-17及其亲本为材料,用SSR标记技术研究杂种与其双亲之间的扩增谱带多态性,以甄别真假杂种.结果发现,所试验的52对SSR引物中有13对引物分别在2个西瓜杂交种和其双亲之间存在扩增条带的多态性,表现为:多数SSR引物对自交系的扩增只出现1条带,但部分引物在某些自交系中扩增出2条带,杂交种条带均为父母本的互补型,很适合做杂交种纯度鉴定.用引物CMCT134b对黑公子和引物CMGA165对04-17进行了各100粒单种子SSR鉴定,所测纯度分别为96%和100%,与田间纯度95.6%和99.7%非常接近,表明SSR标记技术在西瓜杂交种子纯度室内快速检测中的应用前景.     

Analysis of genetic diversity and population structure within Florida coconut (<Emphasis Type="Italic">Cocos nucifera</Emphasis> L.) germplasm using microsatellite DNA,with special emphasis on the Fiji Dwarf cultivar

Using 15 simple sequence repeat (SSR) microsatellite DNA loci, we analyzed genetic variation within Cocos nucifera germplasm collections at two locations in south Florida, representing eight cultivars. The loci were also used in a parentage analysis of progeny of the 'Fiji Dwarf' variety at both locations. A total of 67 alleles were detected, with eight the highest number at any one locus. These loci identified 83 of the 110 individual palms. Gene diversity of the 15 loci ranged from 0.778 to 0.223, with a mean of 0.574. 'Fiji Dwarf', 'Malayan Dwarf', 'Green Ni?o' and 'Red Spicata' cultivars resolve as distinct clusters in a neighbor joining tree using modified Rogers distance, while the tall varieties form two aggregates. The highest gene diversity was found in the tall cultivars (H = 0.583 cumulatively), and the lowest in the 'Malayan Dwarf' (H = 0.202). After the tall coconuts, the 'Fiji Dwarf' was most genetically diverse (H = 0.436), and had the largest number of unique alleles. Genetic identity is highest among the 'Malayan Dwarf' phenotypes, and between the tall varieties. The 'Red Malayan Dwarf' is genetically distinct from the 'Green' and 'Yellow Malayan Dwarf' phenotypes, which cannot be distinguished with the SSR loci used. Off-type 'Malayan Dwarf' phenotypes (putative hybrids with talls) can be identified genotypically. Parentage analyses of 30 'Fiji Dwarf' progeny propagated from five adults surrounded by other cultivars estimate that only 20% of the progeny were out-crossed to the other varieties, while 40-46% were possible selfs. This suggests that a seed-production orchard of the variety maintained at reasonable distance from other varieties, will likely yield only 'Fiji Dwarf' genotypes. Our data are discussed in the context of hypotheses of coconut dissemination around the world.     

A bayesian framework for parentage analysis: the value of genetic and other biological data

We develop fractional allocation models and confidence statistics for parentage analysis in mating systems. The models can be used, for example, to estimate the paternities of candidate males when the genetic mother is known or to calculate the parentage of candidate parent pairs when neither is known. The models do not require two implicit assumptions made by previous models, assumptions that are potentially erroneous. First, we provide formulas to calculate the expected parentage, as opposed to using a maximum likelihood algorithm to calculate the most likely parentage. The expected parentage is superior as it does not assume a symmetrical probability distribution of parentage and therefore, unlike the most likely parentage, will be unbiased. Second, we provide a mathematical framework for incorporating additional biological data to estimate the prior probability distribution of parentage. This additional biological data might include behavioral observations during mating or morphological measurements known to correlate with parentage. The value of multiple sources of information is increased accuracy of the estimates. We show that when the prior probability of parentage is known, and the expected parentage is calculated, fractional allocation provides unbiased estimates of the variance in reproductive success, thereby correcting a problem that has previously plagued parentage analyses. We also develop formulas to calculate the confidence interval in the parentage estimates, thus enabling the assessment of precision. These confidence statistics have not previously been available for fractional models. We demonstrate our models with several biological examples based on data from two fish species that we study, coho salmon (Oncorhychus kisutch) and bluegill sunfish (Lepomis macrochirus). In coho, multiple males compete to fertilize a single female's eggs. We show how behavioral observations taken during spawning can be combined with genetic data to provide an accurate calculation of each male's paternity. In bluegill, multiple males and multiple females may mate in a single nest. For a nest, we calculate the fertilization success and the 95% confidence interval of each candidate parent pair.     

Isolation and characterization of simple sequence repeat loci in the gray tree frog, Hyla chrysoscelis.

A gray tree frog (Hyla chrysoscelis) genomic library was constructed and characterized with regard to the incidence and complexity of simple sequence repeat (SSR) loci. The partial genomic library, containing approximately 10,000 clones with an average-sized insert of 350 bp, was screened with six SSR repeat oligonucleotides (AC, AG, ACG, AGC, AAC, and AAG). Screening identified 31 unique positive clones containing 41 SSR loci. Sequences of tandemly arrayed dinucleotide repeats were more common (36 of 41) than trinucleotide repeats. Twenty-six loci were identified using the AC dinucleotide probe, while 7 loci were identified using the AG dinucleotide probe. An additional 3 AT dinucleotide loci were serendipitously identified. The AT repeats generally comprised the longest dinucleotide repeat loci. The SSR repeat loci reported here should provide potent markers for identity, parentage, and short-lineage determinations in large-scale experiments using gray tree frogs.     

Genetic mapping of agronomic traits in false flax (Camelina sativa subsp. sativa).

The crucifer oilseed plant false flax (Camelina sativa subsp. sativa) possesses numerous valuable agronomic attributes that make it attractive as an alternative spring-sown crop for tight crop rotations. The oil of false flax is particularly rich in polyunsaturated C18-fatty acids, making it a valuable renewable feedstock for the oleochemical industry. Because of the minimal interest in the crop throughout the 20th century, breeding efforts have been limited. In this study, a genetic map for C. sativa was constructed, using amplified fragment length polymorphism (AFLP) markers, in a population of recombinant inbred lines that were developed, through single-seed descent, from a cross between 'Lindo' and 'Licalla', 2 phenotypically distinct parental varieties. Three Brassica simple sequence repeat (SSR) markers were also integrated into the map, and 1 of these shows linkage to oil-content loci in both C. sativa and Brassica napus. Fifty-five other SSR primer combinations showed monomorphic amplification products, indicating partial genome homoeology with the Brassica species. Using data from field trials with different fertilization treatments (0 and 80 kg N/ha) at multiple locations over 3 years, the map was used to localize quantitative trait loci (QTLs) for seed yield, oil content, 1000-seed mass, and plant height. Some yield QTLs were found only with the N0 treatment, and might represent loci contributing to the competitiveness of false flax in low-nutrient soils. The results represent a starting point for future marker-assisted breeding.     

Rapid and efficient resolution of parentage by amplification of short tandem repeats.

Short tandem repeat (STR) loci are highly informative polymorphic loci that are gaining popularity for identity testing. We have conducted parentage testing by using nine STR loci on 50 paternity trios that had been previously tested using VNTR loci. These nine unlinked STR loci are amplified in three multiplex reactions and, when examined for genetic informativeness, provide a combined average power of exclusion of 99.73% (Caucasian data). The informative value of the selected loci is based on extensive STR typing of four racial/ethnic populations. In 37 of the 50 cases, paternity could not be excluded by any of the loci. In the remaining 13 cases, paternity was excluded by at least two of the STR markers. The probability of paternity calculated for the alleged father of each matching trio was > 99% in 36 of the 37 inclusion cases. All data agreed with the results reported using VNTR loci and conventional Southern technology. Our studies validate the use of DNA typing with STR loci for parentage testing, thus providing an accurate, highly sensitive, and rapid assay.     

Mapping quantitative trait loci with dominant and missing markers in various crosses from two inbred lines

Dominant phenotype of a genetic marker provides incomplete information about the marker genotype of an individual. A consequence of using this incomplete information for mapping quantitative trait loci (QTL) is that the inference of the genotype of a putative QTL flanked by a marker with dominant phenotype will depend on the genotype or phenotype of the next marker. This dependence can be extended further until a marker genotype is fully observed. A general algorithm is derived to calculate the probability distribution of the genotype of a putative QTL at a given genomic position, conditional on all observed marker phenotypes in the region with dominant and missing marker information for an individual. The algorithm is implemented for various populations stemming from two inbred lines in the context of mapping QTL. Simulation results show that if only a proportion of markers contain missing or dominant phenotypes, QTL mapping can be almost as efficient as if there were no missing information in the data. The efficiency of the analysis, however, may decrease substantially when a very large proportion of markers contain missing or dominant phenotypes and a genetic map has to be reconstructed first on the same data as well. So it is important to combine dominant markers with codominant markers in a QTL mapping study. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analysis of genetic diversity and relationships among waxy maize inbred lines in Korea using SSR markers

Information regarding the genetic diversity and genetic relationships among elite inbred lines is necessary to improve new cultivars in maize breeding programs. In this study, genetic diversity and genetic relationships were investigated among 84 waxy maize inbred lines using 50 SSR markers. A total of 269 alleles were identified at all the loci with an average of 5.38 and a range between 2 and 13 alleles per locus. The gene diversity values varied from 0.383 to 0.923 with an average of 0.641. The cluster tree generated using the described SSR markers recognized two major groups at 32% genetic similarity. Group I included 33 inbred lines while group II included 51 inbred lines. The clustering patterns of most of the waxy maize inbred lines did not clearly agree with their source, pedigree or geographic location. The average GS among all inbred lines was 35.7 ± 10.8. Analysis of waxy maize inbred lines collected from Korea and China at 50 SSR loci revealed higher values of average number of alleles (4.9) and gene diversity (0.638) in Korean inbred lines as compared to Chinese inbred lines (3.5 and 0.563, respectively). The information obtained from the present studies would be very useful for maize breeding programs in Korea.     

DNA-identification of rice varieties (Oryza sativa L.) of Ukrainian selection

Genetic diversity of 11 rice varieties of the Institute of rice UAAS was assessed using a set of 25 SSR loci. Based on analyses of 12 polymorphic microsatillite loci the procedures for composing genetic formulas of varieties and their identification have been elaborated. UPGMA-cluster-analysis based on genetic distance coefficients clearly separated all the varieties into two groups, and showed that the Ukrainian rice varieties are closely related. Although the genetic diversity was low, SSRs proved to be an efficient tool in assessing the genetic diversity of rice genotypes.