Domestication reshaped the genetic basis of inbreeding depression in a maize landrace compared to its wild relative, teosinte

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. The extent to which the genetic load of mutations contributing to inbreeding depression is due to large-effect mutations versus variants with very small individual effects is unknown and may be affected by population history. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits, congruent with the greater segregating genetic load in the maize population that we predicted from sequence data. Parental breeding values were highly consistent between outcross and selfed offspring, indicating that additive effects determine most of the genetic value even in the presence of strong inbreeding depression. We developed a novel linkage scan to identify quantitative trait loci (QTL) representing large-effect rare variants carried by only a single parent, which were more important in teosinte than maize. Teosinte also carried more putative juvenile-acting lethal variants identified by segregation distortion. These results suggest a mixture of mostly polygenic, small-effect partially recessive effects in linkage disequilibrium underlying inbreeding depression, with an additional contribution from rare larger-effect variants that was more important in teosinte but depleted in maize following the domestication bottleneck. Purging associated with the maize domestication bottleneck may have selected against some large effect variants, but polygenic load is harder to purge and overall segregating mutational burden increased in maize compared to teosinte.

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. Mating of close relatives increases the probability that an individual inherits two non-functioning mutations at the same gene, resulting in lower fitness of such matings. We do not know the extent to which inbreeding depression is due to mutations with large-effects versus small-effect polygenic variants. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits and we found that this was consistent with higher predicted ‘genetic load’ in maize based solely on the evolutionary conservation of the sequence variants observed in the population. We also mapped genome positions associated with inbreeding depression, identifying more and larger-effect genetic variants in teosinte than maize. These results suggest that during domestication, some of the rare large-effect variants in teosinte were bred out, but many genetic variants of small effects on inbreeding depression increased in frequency maize.

Genomic prediction across years in a maize doubled haploid breeding program to accelerate early-stage testcross testing

Genomic selection with a multiple-year training population dataset could accelerate early-stage testcross testing by skipping the first-stage yield testing, which significantly saves the time and cost of early-stage testcross testing. With the development of doubled haploid (DH) technology, the main task for a maize breeder is to estimate the breeding values of thousands of DH lines annually. In early-stage testcross testing, genomic selection (GS) offers the opportunity of replacing expensive multiple-environment phenotyping and phenotypic selection with lower-cost genotyping and genomic estimated breeding value (GEBV)-based selection. In the present study, a total of 1528 maize DH lines, phenotyped in multiple-environment trials in three consecutive years and genotyped with a low-cost per-sample genotyping platform of rAmpSeq, were used to explore how to implement GS to accelerate early-stage testcross testing. Results showed that the average prediction accuracy estimated from the cross-validation schemes was above 0.60 across all the scenarios. The average prediction accuracies estimated from the independent validation schemes ranged from 0.23 to 0.32 across all the scenarios, when the one-year datasets were used as training population (TRN) to predict the other year data as testing population (TST). The average prediction accuracies increased to a range from 0.31 to 0.42 across all the scenarios, when the two-years datasets were used as TRN. The prediction accuracies increased to a range from 0.50 to 0.56, when the TRN consisted of two-years of breeding data and 50% of third year's data converted from TST to TRN. This information showed that GS with a multiple-year TRN set offers the opportunity to accelerate early-stage testcross testing by skipping the first-stage yield testing, which significantly saves the time and cost of early-stage testcross testing.

Genotyping by sequencing for SNP marker development in onion

Onion (Allium cepa) is not highly tractable for development of molecular markers due to its large (16 gigabases per 1C) nuclear genome. Single nucleotide polymorphisms (SNPs) are useful for genetic characterization and marker-aided selection of onion because of codominance and common occurrence in elite germplasm. We completed genotyping by sequencing (GBS) to identify SNPs in onion using 46 F₂ plants, parents of the F₂ plants (Ailsa Craig 43 and Brigham Yellow Globe 15-23), two doubled haploid (DH) lines (DH2107 and DH2110), and plants from 94 accessions in the USDA National Plant Germplasm System (NPGS). SNPs were called using the TASSEL 3.0 Universal Network Enabled Analysis (UNEAK) bioinformatics pipeline. Sequences from the F₂ and DH plants were used to construct a pseudo-reference genome against which genotypes from all accessions were scored. Quality filters were used to identify a set of 284 high quality SNPs, which were placed onto an existing genetic map for the F₂ family. Accessions showed a moderate level of diversity (mean H_e = 0.341) and evidence of inbreeding (mean F = 0.592). GBS is promising for SNP discovery in onion, although lack of a reference genome required extensive custom scripts for bioinformatics analyses to identify high quality markers.

High-Density Genetic Linkage Mapping of Lepidium Based on Genotyping-by-Sequencing SNPs and Segregating Contig Tag Haplotypes

Lepidium campestre has been targeted for domestication as future oilseed and catch crop. Three hundred eighty plants comprising genotypes of L. campestre, Lepidium heterophyllum, and their interspecific F2 mapping population were genotyped using genotyping by sequencing (GBS), and the generated polymorphic markers were used for the construction of high-density genetic linkage map. TASSEL-GBS, a reference genome-based pipeline, was used for this analysis using a draft L. campestre whole genome sequence. The analysis resulted in 120,438 biallelic single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) above 0.01. The construction of genetic linkage map was conducted using MSTMap based on phased SNPs segregating in 1:2:1 ratio for the F2 individuals, followed by genetic mapping of segregating contig tag haplotypes as dominant markers against the linkage map. The final linkage map consisted of eight linkage groups (LGs) containing 2,330 SNP markers and spanned 881 Kosambi cM. Contigs (10,302) were genetically mapped to the eight LGs, which were assembled into pseudomolecules that covered a total of ∼120.6 Mbp. The final size of the pseudomolecules ranged from 9.4 Mbp (LG-4) to 20.4 Mpb (LG-7). The following major correspondence between the eight Lepidium LGs (LG-1 to LG-8) and the five Arabidopsis thaliana (At) chromosomes (Atx-1-Atx-5) was revealed through comparative genomics analysis: LG-1&2_Atx-1, LG-3_Atx-2&3, LG-4_Atx-2, LG-5_Atx-2&Atx-3, LG-6_Atx-4&5, LG-7_Atx-4, and LG-8_Atx-5. This analysis revealed that at least 66% of the sequences of the LGs showed high collinearity with At chromosomes. The sequence identity between the corresponding regions of the LGs and At chromosomes ranged from 80.6% (LG-6) to 86.4% (LG-8) with overall mean of 82.9%. The map positions on Lepidium LGs of the homologs of 24 genes that regulate various traits in A. thaliana were also identified. The eight LGs revealed in this study confirm the previously reported (1) haploid chromosome number of eight in L. campestre and L. heterophyllum and (2) chromosomal fusion, translocation, and inversion events during the evolution of n = 8 karyotype in ancestral species shared by Lepidium and Arabidopsis to n = 5 karyotype in A. thaliana. This study generated highly useful genomic tools and resources for Lepidium that can be used to accelerate its domestication.

Haplotyping the Vitis collinear core genome with rhAmpSeq improves marker transferability in a diverse genus

Transferable DNA markers are essential for breeding and genetics. Grapevine (Vitis) breeders utilize disease resistance alleles from congeneric species ~20 million years divergent, but existing Vitis marker platforms have cross-species transfer rates as low as 2%. Here, we apply a marker strategy targeting the inferred Vitis core genome. Incorporating seven linked-read de novo assemblies and three existing assemblies, the Vitis collinear core genome is estimated to converge at 39.8 Mb (8.67% of the genome). Adding shotgun genome sequences from 40 accessions enables identification of conserved core PCR primer binding sites flanking polymorphic haplotypes with high information content. From these target regions, we develop 2,000 rhAmpSeq markers as a PCR multiplex and validate the panel in four biparental populations spanning the diversity of the Vitis genus, showing transferability increases to 91.9%. This marker development strategy should be widely applicable for genetic studies in many taxa, particularly those ~20 million years divergent.

Identification of genetic variants associated with maize flowering time using an extremely large multi-genetic background population

Flowering time is one of the major adaptive traits in domestication of maize and an important selection criterion in breeding. To detect more maize flowering time variants we evaluated flowering time traits using an extremely large multi- genetic background population that contained more than 8000 lines under multiple Sino-United States environments. The population included two nested association mapping (NAM) panels and a natural association panel. Nearly 1 million single-nucleotide polymorphisms (SNPs) were used in the analyses. Through the parallel linkage analysis of the two NAM panels, both common and unique flowering time regions were detected. Genome wide, a total of 90 flowering time regions were identified. One-third of these regions were connected to traits associated with the environmental sensitivity of maize flowering time. The genome-wide association study of the three panels identified nearly 1000 flowering time-associated SNPs, mainly distributed around 220 candidate genes (within a distance of 1 Mb). Interestingly, two types of regions were significantly enriched for these associated SNPs - one was the candidate gene regions and the other was the approximately 5 kb regions away from the candidate genes. Moreover, the associated SNPs exhibited high accuracy for predicting flowering time.

Recombination in diverse maize is stable, predictable, and associated with genetic load

Among the fundamental evolutionary forces, recombination arguably has the largest impact on the practical work of plant breeders. Varying over 1,000-fold across the maize genome, the local meiotic recombination rate limits the resolving power of quantitative trait mapping and the precision of favorable allele introgression. The consequences of low recombination also theoretically extend to the species-wide scale by decreasing the power of selection relative to genetic drift, and thereby hindering the purging of deleterious mutations. In this study, we used genotyping-by-sequencing (GBS) to identify 136,000 recombination breakpoints at high resolution within US and Chinese maize nested association mapping populations. We find that the pattern of cross-overs is highly predictable on the broad scale, following the distribution of gene density and CpG methylation. Several large inversions also suppress recombination in distinct regions of several families. We also identify recombination hotspots ranging in size from 1 kb to 30 kb. We find these hotspots to be historically stable and, compared with similar regions with low recombination, to have strongly differentiated patterns of DNA methylation and GC content. We also provide evidence for the historical action of GC-biased gene conversion in recombination hotspots. Finally, using genomic evolutionary rate profiling (GERP) to identify putative deleterious polymorphisms, we find evidence for reduced genetic load in hotspot regions, a phenomenon that may have considerable practical importance for breeding programs worldwide.

Grapevine powdery mildew resistance and susceptibility loci identified on a high-resolution SNP map

Improved efficacy and durability of powdery mildew resistance can be enhanced via knowledge of the genetics of resistance and susceptibility coupled with the development of high-resolution maps to facilitate the stacking of multiple resistance genes and other desirable traits. We studied the inheritance of powdery mildew (Erysiphe necator) resistance and susceptibility of wild Vitis rupestris B38 and cultivated V. vinifera 'Chardonnay', finding evidence for quantitative variation. Molecular markers were identified using genotyping-by-sequencing, resulting in 16,833 single nucleotide polymorphisms (SNPs) based on alignment to the V. vinifera 'PN40024' reference genome sequence. With an average density of 36 SNPs/Mbp and uniform coverage of the genome, this 17K set was used to identify 11 SNPs on chromosome 7 associated with a resistance locus from V. rupestris B38 and ten SNPs on chromosome 9 associated with a locus for susceptibility from 'Chardonnay' using single marker association and linkage disequilibrium analysis. Linkage maps for V. rupestris B38 (1,146 SNPs) and 'Chardonnay' (1,215 SNPs) were constructed and used to corroborate the 'Chardonnay' locus named Sen1 (Susceptibility to Erysiphe necator 1), providing the first insight into the genetics of susceptibility to powdery mildew from V. vinifera. The identification of markers associated with a susceptibility locus in a V. vinifera background can be used for negative selection among breeding progenies. This work improves our understanding of the nature of powdery mildew resistance in V. rupestris B38 and 'Chardonnay', while applying next-generation sequencing tools to advance grapevine genomics and breeding.

Comprehensive genotyping of the USA national maize inbred seed bank

BACKGROUND

Genotyping by sequencing, a new low-cost, high-throughput sequencing technology was used to genotype 2,815 maize inbred accessions, preserved mostly at the National Plant Germplasm System in the USA. The collection includes inbred lines from breeding programs all over the world.

RESULTS

The method produced 681,257 single-nucleotide polymorphism (SNP) markers distributed across the entire genome, with the ability to detect rare alleles at high confidence levels. More than half of the SNPs in the collection are rare. Although most rare alleles have been incorporated into public temperate breeding programs, only a modest amount of the available diversity is present in the commercial germplasm. Analysis of genetic distances shows population stratification, including a small number of large clusters centered on key lines. Nevertheless, an average fixation index of 0.06 indicates moderate differentiation between the three major maize subpopulations. Linkage disequilibrium (LD) decays very rapidly, but the extent of LD is highly dependent on the particular group of germplasm and region of the genome. The utility of these data for performing genome-wide association studies was tested with two simply inherited traits and one complex trait. We identified trait associations at SNPs very close to known candidate genes for kernel color, sweet corn, and flowering time; however, results suggest that more SNPs are needed to better explore the genetic architecture of complex traits.

CONCLUSIONS

The genotypic information described here allows this publicly available panel to be exploited by researchers facing the challenges of sustainable agriculture through better knowledge of the nature of genetic diversity.

Population genomic and genome-wide association studies of agroclimatic traits in sorghum

Accelerating crop improvement in sorghum, a staple food for people in semiarid regions across the developing world, is key to ensuring global food security in the context of climate change. To facilitate gene discovery and molecular breeding in sorghum, we have characterized ~265,000 single nucleotide polymorphisms (SNPs) in 971 worldwide accessions that have adapted to diverse agroclimatic conditions. Using this genome-wide SNP map, we have characterized population structure with respect to geographic origin and morphological type and identified patterns of ancient crop diffusion to diverse agroclimatic regions across Africa and Asia. To better understand the genomic patterns of diversification in sorghum, we quantified variation in nucleotide diversity, linkage disequilibrium, and recombination rates across the genome. Analyzing nucleotide diversity in landraces, we find evidence of selective sweeps around starch metabolism genes, whereas in landrace-derived introgression lines, we find introgressions around known height and maturity loci. To identify additional loci underlying variation in major agroclimatic traits, we performed genome-wide association studies (GWAS) on plant height components and inflorescence architecture. GWAS maps several classical loci for plant height, candidate genes for inflorescence architecture. Finally, we trace the independent spread of multiple haplotypes carrying alleles for short stature or long inflorescence branches. This genome-wide map of SNP variation in sorghum provides a basis for crop improvement through marker-assisted breeding and genomic selection.

Comparative population genomics of maize domestication and improvement

Domestication and plant breeding are ongoing 10,000-year-old evolutionary experiments that have radically altered wild species to meet human needs. Maize has undergone a particularly striking transformation. Researchers have sought for decades to identify the genes underlying maize evolution, but these efforts have been limited in scope. Here, we report a comprehensive assessment of the evolution of modern maize based on the genome-wide resequencing of 75 wild, landrace and improved maize lines. We find evidence of recovery of diversity after domestication, likely introgression from wild relatives, and evidence for stronger selection during domestication than improvement. We identify a number of genes with stronger signals of selection than those previously shown to underlie major morphological changes. Finally, through transcriptome-wide analysis of gene expression, we find evidence both consistent with removal of cis-acting variation during maize domestication and improvement and suggestive of modern breeding having increased dominance in expression while targeting highly expressed genes.

Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.

A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species

Advances in next generation technologies have driven the costs of DNA sequencing down to the point that genotyping-by-sequencing (GBS) is now feasible for high diversity, large genome species. Here, we report a procedure for constructing GBS libraries based on reducing genome complexity with restriction enzymes (REs). This approach is simple, quick, extremely specific, highly reproducible, and may reach important regions of the genome that are inaccessible to sequence capture approaches. By using methylation-sensitive REs, repetitive regions of genomes can be avoided and lower copy regions targeted with two to three fold higher efficiency. This tremendously simplifies computationally challenging alignment problems in species with high levels of genetic diversity. The GBS procedure is demonstrated with maize (IBM) and barley (Oregon Wolfe Barley) recombinant inbred populations where roughly 200,000 and 25,000 sequence tags were mapped, respectively. An advantage in species like barley that lack a complete genome sequence is that a reference map need only be developed around the restriction sites, and this can be done in the process of sample genotyping. In such cases, the consensus of the read clusters across the sequence tagged sites becomes the reference. Alternatively, for kinship analyses in the absence of a reference genome, the sequence tags can simply be treated as dominant markers. Future application of GBS to breeding, conservation, and global species and population surveys may allow plant breeders to conduct genomic selection on a novel germplasm or species without first having to develop any prior molecular tools, or conservation biologists to determine population structure without prior knowledge of the genome or diversity in the species.

Genetic properties of the maize nested association mapping population

Maize genetic diversity has been used to understand the molecular basis of phenotypic variation and to improve agricultural efficiency and sustainability. We crossed 25 diverse inbred maize lines to the B73 reference line, capturing a total of 136,000 recombination events. Variation for recombination frequencies was observed among families, influenced by local (cis) genetic variation. We identified evidence for numerous minor single-locus effects but little two-locus linkage disequilibrium or segregation distortion, which indicated a limited role for genes with large effects and epistatic interactions on fitness. We observed excess residual heterozygosity in pericentromeric regions, which suggested that selection in inbred lines has been less efficient in these regions because of reduced recombination frequency. This implies that pericentromeric regions may contribute disproportionally to heterosis.

The genetic architecture of maize flowering time

Flowering time is a complex trait that controls adaptation of plants to their local environment in the outcrossing species Zea mays (maize). We dissected variation for flowering time with a set of 5000 recombinant inbred lines (maize Nested Association Mapping population, NAM). Nearly a million plants were assayed in eight environments but showed no evidence for any single large-effect quantitative trait loci (QTLs). Instead, we identified evidence for numerous small-effect QTLs shared among families; however, allelic effects differ across founder lines. We identified no individual QTLs at which allelic effects are determined by geographic origin or large effects for epistasis or environmental interactions. Thus, a simple additive model accurately predicts flowering time for maize, in contrast to the genetic architecture observed in the selfing plant species rice and Arabidopsis.

Population structure and genetic diversity of New World maize races assessed by DNA microsatellites

Because of the economic importance of maize and its scientific importance as a model system for studies of domestication, its evolutionary history is of general interest. We analyzed the population genetic structure of maize races by genotyping 964 individual plants, representing almost the entire set of ∼350 races native to the Americas, with 96 microsatellites. Using Bayesian clustering, we detected four main clusters consisting of highland Mexican, northern United States (US), tropical lowland, and Andean races. Phylogenetic analysis indicated that the southwestern US was an intermediary stepping stone between Mexico and the northern US. Furthermore, southeastern US races appear to be of mixed northern flint and tropical lowland ancestry, while lowland middle South American races are of mixed Andean and tropical lowland ancestry. Several cases of post-Columbian movement of races were detected, most notably from the US to South America. Of the four main clusters, the highest genetic diversity occurs in highland Mexican races, while diversity is lowest in the Andes and northern US. Isolation by distance appears to be the main factor underlying the historical diversification of maize. We identify highland Mexico and the Andes as potential sources of genetic diversity underrepresented among elite lines used in maize breeding programs.

Panzea: an update on new content and features

Panzea (http://www.panzea.org), the public web site of the project ‘Molecular and Functional Diversity in the Maize Genome’, has expanded over the past two years in data content, display tools and informational sections. The most significant data content expansions occurred for the single nucleotide polymorphism (SNP), sequencing, isozyme and phenotypic data types. We have enhanced our existing web display tools and have launched a number of new tools for data display and analysis. For example, we have implemented one that allows users to find polymorphisms between two accessions, a geographic map tool to visualize the geographic distribution of SNPs, simple sequence repeats (SSRs) and isozyme alleles and a graphical view of the placement of Panzea markers and genes/loci on genetic and physical maps. One goal of the informatics component of our project has been to generate code that can be used by other groups. We have enhanced our existing code base and have made our new tools available. Finally, we have also made available new informational sections as part of our educational and outreach efforts.

Genetic Structure Inside a Declining Red Oak Community in Old-Growth Forest

Problems with oak regeneration have been documented in the last 50 years at numerous sites in the Midwestern United States. We applied nuclear microsatellites to examine the demographic and fine-scale spatial genetic structure of red oaks in two old-growth stands in Indiana. Oaks in one stand have declined in numbers over the past several decades whereas oaks in the other, smaller stand have increased. Large amounts of genetic variation were maintained within stands, and there was slight but significant differentiation among stands. There was significant but weak isolation by distance genetic structure within the large stand, likely reflecting family structure. No significant differences exist in allele frequencies or in levels of genetic diversity between cohorts that remain well represented within each stand, even between medium-sized adults and those antedating European settlement of the area. However, a virtual absence of smaller size classes in the forest interior of the large stand represents the early stages of a genetic bottleneck in what had been the core habitat of this stand. Whether future generations of this old-growth stand will retain the present genetic character depends on the oaks regenerating at the forest margins, absent any major changes in disturbance regimes. Similar demographic and genetic dynamics are likely occurring in a large number of remnant oak forests across the Midwest.

Effects of native forest regeneration practices on genetic diversity in Eucalyptus consideniana

Impacts of forest harvesting and regeneration practices on genetic diversity in the Australian native forest species Eucalyptus consideniana Maiden (yertchuk) were examined using 29 Mendelian DNA markers (18 RFLPs and 11 microsatellites). Two replicate logging coupes were studied from each of the two most commonly employed silvicultural treatments: clear felling with aerial re-sowing and the seed tree system. For each coupe, genetic diversity measures were compared between a sample of the sapling regeneration and a corresponding control sample from bordering unharvested trees. When calculations were performed over all 29 loci, significant reductions of allelic richness (AR), effective number of alleles (AE) and/or expected heterozygosity (HE) were detected on one or both of seed tree coupes, but on neither of the clear falls. When calculations were performed over the 11 microsatellites alone, all three measures, AR, AE and HE, were significantly reduced on both of the seed replicates but on neither of the two clear falls. In contrast, when the RFLPs were examined separately, there were no significant reductions of diversity on either of the two seed tree coupes or on the two clear falls. These results suggest that genetic erosion is more likely under the seed tree system than under clear-felling with aerial re-sowing and that there is greater statistical power to detect it with microsatellites than with RFLPs. A Monte Carlo simulation to test the statistical significance of the number of apparently lost or gained alleles showed that significant losses of alleles above specified threshold frequencies occurred only in the two seed tree replicates. Three of the four control and regeneration population pairs were significantly differentiated, as indicated by exact tests or by pairwise FST estimates. Comparisons of CONTML dendrograms, constructed for the regeneration populations only versus the control populations only, indicated that genetic drift was significantly promoted under forest management. No significant decreases in observed heterozygosity, or increases in the panmictic index (f), were observed in any of the comparisons suggesting that inbreeding was not promoted by a single rotation of forest management.

Prospects for inferring pairwise relationships with single nucleotide polymorphisms

An extraordinarily large number of single nucleotide polymorphisms (SNPs) are now available in humans as well as in other model organisms. Technological advancements may soon make it feasible to assay hundreds of SNPs in virtually any organism of interest. One potential application of SNPs is the determination of pairwise genetic relationships in populations without known pedigrees. Although microsatellites are currently the marker of choice for this purpose, the number of independently segregating microsatellite markers that can be feasibly assayed is limited. Thus, it can be difficult to distinguish reliably some classes of relationship (e.g. full-sibs from half-sibs) with microsatellite data alone. We assess, via Monte Carlo computer simulation, the potential for using a large panel of independently segregating SNPs to infer genetic relationships, following the analytical approach of Blouin et al. (1996). We have explored a 'best case scenario' in which 100 independently segregating SNPs are available. For discrimination among single-generation relationships or for the identification of parent-offspring pairs, it appears that such a panel of moderately polymorphic SNPs (minor allele frequency of 0.20) will provide discrimination power equivalent to only 16-20 independently segregating microsatellites. Although newly available analytical methods that can account for tight genetic linkage between markers will, in theory, allow improved estimation of relationships using thousands of SNPs in highly dense genomic scans, in practice such studies will only be feasible in a handful of model organisms. Given the comparable amount of effort required for the development of both types of markers, it seems that microsatellites will remain the marker of choice for relationship estimation in nonmodel organisms, at least for the foreseeable future.

Dinucleotide microsatellites from Eucalyptus sieberi: inheritance, diversity, and improved scoring of single-base differences

Eight dinucleotide microsatellites were developed in Eucalyptus sieberi L. Johnson (silvertop ash), a member of the subgenus Eucalyptus. Transfer of six of these to the subgenus Symphyomyrtus and their Mendelian inheritance are demonstrated using a full-sib cross in Eucalyptus nitens. Genetic diversity parameters are presented for the eight loci based on a sample of 100 old-growth E. sieberi trees from a single natural stand. One locus, Es266, had an atypically high fixation index, and significantly deviated from Hardy-Weinberg equilibrium genotypic proportions, indicating the likely presence of null alleles. Two of the loci, Es076 and Es140, had many alleles that differed in size by only a single base pair, possibly because of short poly(A) or poly(T) stretches in their flanking regions. These two loci were by far the most polymorphic, but were difficult to score reliably on a capillary DNA sequencer. Reliability of scoring of these two one-base microsatellite loci was markedly improved by the incorporation of internal reference alleles into each sample analysed.Key words: SSRs, single base pair alleles, null alleles, internal reference alleles.

Dinucleotide microsatellites from Eucalyptus sieberi: inheritance, diversity, and improved scoring of single-base differences

Eight dinucleotide microsatellites were developed in Eucalyptus sieberi L. Johnson (silvertop ash), a member of the subgenus Eucalyptus. Transfer of six of these to the subgenus Symphyomyrtus and their Mendelian inheritance are demonstrated using a full-sib cross in Eucalyptus nitens. Genetic diversity parameters are presented for the eight loci based on a sample of 100 old-growth E. sieberi trees from a single natural stand. One locus, Es266, had an atypically high fixation index, and significantly deviated from Hardy-Weinberg equilibrium genotypic proportions, indicating the likely presence of null alleles. Two of the loci, Es076 and Es140, had many alleles that differed in size by only a single base pair, possibly because of short poly(A) or poly(T) stretches in their flanking regions. These two loci were by far the most polymorphic, but were difficult to score reliably on a capillary DNA sequencer. Reliability of scoring of these two one-base microsatellite loci was markedly improved by the incorporation of internal reference alleles into each sample analysed.

Single tree genetic linkage mapping in conifers using haploid DNA from megagametophytes

We report a unique application of the Random Amplified Polymorphic DNA (RAPD) technique for genetic linkage mapping of a single spruce tree using haploid DNA from megagametophyte tissue of individual seeds. Sixty-one segregating loci were analysed for reproducibility, inheritance and linkage. Forty-seven of the 61 markers were distributed into 12 linkage groups and covered 873.8 cM. The 14 markers not associated with any linkage group should be assigned to linkage groups as more markers are added to the map. This new approach quickly provides molecular genetic markers that are simple to evaluate for constructing genetic linkage maps and for other related genetic studies in forest tree species.

RNA editing in the mitochondria of a conifer

A 712-base portion of the mitochondrial gene coxI and the corresponding portion of the coxI transcript were amplified by PCR and by RT-PCR, respectively, from the gymnosperm western red cedar. Sequence comparison of amplified coxI DNA and mRNA revealed 26 C-to-U differences that are best explained by RNA editing of the type known to occur in angiosperms. This finding suggests that mitochondrial RNA-editing of the C-to-U type arose before the divergence of gymnosperms and angiosperms and can be considered a feature common to all higher plants.

Segregation of random amplified DNA markers in F1 progeny of conifers

The recently developed approach to deriving genetic markers via amplification of random DNA segments with single primers of arbitrary nucleotide sequence was tested for its utility in genetic linkage mapping studies with conifers. Reaction conditions were optimized to reproducibly yield clean and specific amplification products. Template DNA from several genotypes of Douglas-fir (Pseudotsuga menziesii) and white spruce (Picea glauca) were tested against eight ten-base oligonucleotide primers. Most of the tested primer/parent tree combinations yielded polymorphic PCR products ("RAPD" markers). Selected primers were then used in PCR reactions with template DNA isolated from offspring in Douglas-fir and black spruce diallel crosses among the same parental lines. The diallel study confirmed the appropriate inheritance of RAPD markers in the F1 generation. The value of these dominant RAPD markers for genetic linkage mapping in trees was established from both theoretical and applied perspectives.