Toward Guidelines for Harvest Intensities and Regeneration Targets with Minimal Impact Upon Retained Genetic Diversity in Central Hardwood Tree Species

There is an urgent need for a coordinated and systematic approach to the in situ conservation of the genetic resources of commercially important forest tree species in the Central Hardwoods. Effective in situ management of genetic resources would benefit from clear guidelines for how many adult trees can be harvested with minimal impact on allelic diversity. We are constructing a computer model for this purpose, and present preliminary results based upon replicate harvests of a virtual forest stand consisting of 200 adult trees. Our model explores how much regeneration is needed so that there is no more than a 10 percent risk of retaining less than 90 percent of the original allelic diversity. In the absence of regeneration, up to 55 percent of the adult trees can be harvested without exceeding the 10 percent risk level. At higher harvest intensities, locally-derived regeneration is needed to replace the alleles removed from the adult population. When all 200 adult trees are harvested, the 10 percent risk level is not exceeded if there are at least 116 regenerants, provided that these are derived from pre-harvest random mating among the adults. In the presence of substantial pollen flow from a genetically differentiated outside pollen source (e.g., 10-20 percent pollen flow), the minimum amount of regeneration needed is reduced. This indicates that outside pollen can be more efficient, relative to pollen from within the stand, at replacing alleles lost from the adult population

Different histories but similar genetic diversity and structure for black walnut in Indiana and Missouri

—Missouri and Indiana have markedly different histories of glaciation and recolonization by forest trees. These states also differ in land use patterns and degree of anthropogenic landscape change such as forest fragmentation. To determine the overall effects of these and other demographic differences on the levels of genetic diversity and structure in black walnut (Juglans nigra L.) more than 550 total black walnut trees from nine populations in Indiana and 10 in Missouri were sampled and analyzed using 12 nuclear microsatellite loci. Although genetic diversity parameters such as allelic richness and expected heterozygosity were high overall, they varied little among populations and their mean values for the two states were not significantly different. Pairwise genetic distance values between all population pairs ranged from 0.012-0.159, but no significant pattern of isolation by distance was detected. The estimate of the degree of genetic differentiation between states (FPT = 0.0009) was very small and not significant, indicating that differences between states explained an inconsequential portion of the total variance. The observed low levels of local and regional genetic structure indicate that high levels of pollen flow have buffered black walnut from the genetic consequences of founder effects and genetic drift in both geologic and recent time scales

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, smalleffect 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

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

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

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 F2 plants, parents of the F2 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 F2 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 F2 family. Accessions showed a moderate level of diversity (mean He = 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.The accepted manuscript in pdf format is listed with the files at the bottom of this page. The presentation of the authors' names and (or) special characters in the title of the manuscript may differ slightly between what is listed on this page and what is listed in the pdf file of the accepted manuscript; that in the pdf file of the accepted manuscript is what was submitted by the author

TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline.

Genotyping by sequencing (GBS) is a next generation sequencing based method that takes advantage of reduced representation to enable high throughput genotyping of large numbers of individuals at a large number of SNP markers. The relatively straightforward, robust, and cost-effective GBS protocol is currently being applied in numerous species by a large number of researchers. Herein we describe a bioinformatics pipeline, TASSEL-GBS, designed for the efficient processing of raw GBS sequence data into SNP genotypes. The TASSEL-GBS pipeline successfully fulfills the following key design criteria: (1) Ability to run on the modest computing resources that are typically available to small breeding or ecological research programs, including desktop or laptop machines with only 8-16 GB of RAM, (2) Scalability from small to extremely large studies, where hundreds of thousands or even millions of SNPs can be scored in up to 100,000 individuals (e.g., for large breeding programs or genetic surveys), and (3) Applicability in an accelerated breeding context, requiring rapid turnover from tissue collection to genotypes. Although a reference genome is required, the pipeline can also be run with an unfinished "pseudo-reference" consisting of numerous contigs. We describe the TASSEL-GBS pipeline in detail and benchmark it based upon a large scale, species wide analysis in maize (Zea mays), where the average error rate was reduced to 0.0042 through application of population genetic-based SNP filters. Overall, the GBS assay and the TASSEL-GBS pipeline provide robust tools for studying genomic diversity

Independent Molecular Basis of Convergent Highland Adaptation in Maize

Convergent evolution is the independent evolution of similar traits in different species or lineages of the same species; this often is a result of adaptation to similar environments, a process referred to as convergent adaptation. We investigate here the molecular basis of convergent adaptation in maize to highland climates in Mesoamerica and South America, using genome-wide SNP data. Taking advantage of archaeological data on the arrival of maize to the highlands, we infer demographic models for both populations, identifying evidence of a strong bottleneck and rapid expansion in South America. We use these models to then identify loci showing an excess of differentiation as a means of identifying putative targets of natural selection and compare our results to expectations from recently developed theory on convergent adaptation. Consistent with predictions across a wide parameter space, we see limited evidence for convergent evolution at the nucleotide level in spite of strong similarities in overall phenotypes. Instead, we show that selection appears to have predominantly acted on standing genetic variation and that introgression from wild teosinte populations appears to have played a role in highland adaptation in Mexican maize