Haplotype structure at seven barley genes: relevance to gene pool bottlenecks, phylogeny of ear type and site of barley domestication

Archaeological remains indicate that the origin of western agriculture occurred in a brief period about 10,500 years ago in a region of the Middle East known as the Fertile Crescent, where the wild progenitors of several key agricultural cereal species are endemic. Domestication entailed the appearance of agronomic traits such as seed size and threshability. For a representative sample of 20 domesticated barley (Hordeum vulgare) lines, including 13 two-rowed and 7 six-rowed varieties, we determined the haplotypes at seven loci-Adh2, Adh3, Amy1, Dhn9, GAPDH, PEPC and WAXY encompassing 5,616 bases per line-and compared them to the haplotypes at the same loci for 25 wild forms (Hordeum spontaneum) collected within and outside the Fertile Crescent. In comparisons of wild versus domesticated barley, the number of haplotypes (70 vs. 17), average nucleotide diversity, pi, (0.0077 vs. 0.0028), and Watterson's theta at silent sites (0.0104 vs. 0.0028) was reduced in domesticated lines. Two loci, Amy1 and PEPC, were monomorphic in domesticated lines; Amy1 and GAPDH produced significant values of Tajima's D. At GAPDH, pi was slightly higher in domesticated than wild forms, due to divergent high-frequency haplotypes; for the remaining six loci, 87% of nucleotide diversity has been lost in the domesticated forms. Bottlenecks acting on neutrally evolving loci either during the domestication process, during subsequent breeding, or both, are sufficient to account for reduced diversity and the results of Tajima's test, without the need to evoke selection at these loci. Phylogenetic networks data uncover distinct wild and domesticated barley genotypes and suggest that barley may have been domesticated in the Jordan valley. Because, based on AFLP data, the domesticated Turkish cultivars had a genetic basis as large as that present in large germplasm collections, all comparisons provided in this paper are of general value more than being restricted to the Turkish barley germplasm.     

Genetic Variation within and Between Domesticated Chinook Salmon, Oncorhynchus tshawytscha, Strains and their Progenitor Populations

Domesticated chinook salmon strains in British Columbia (BC), Canada are believed to have originated primarily from populations of the Big Qualicum (BQ) River and Robertson Creek (RC) on Vancouver Island in the early 1980s. The number of parental fish that gave rise to the domesticated strains and their subsequent breeding history during approximately five ensuing generations of domestication were not documented. Genetic variation at 13 microsatellite loci was examined in samples from two domesticated strains and the two progenitor populations to determine the genetic relationships among them. The domesticated strains had lower allelic diversity and tended to have lower levels of expected heterozygosity than did the BQ and RC progenitor populations. Only three alleles over all 13 loci were detected in the domesticated strains that were not present in the BQ and RC samples, whereas the progenitor strains possessed over 25 (BQ) and 43 (RC) private alleles. Genetic distance and F_ST values also indicated a closer relationship of the domesticated strains with the BQ than the RC population. One domesticated strain had a significant excess of heterozygosity compared with that expected under conditions of mutation-drift equilibrium, indicative of a recent genetic bottleneck. Genetic differentiation between the domesticated strains was as great as that distinguishing them from the progenitor populations, indicating that the genetic base of domesticated chinook salmon could be increased by hybridization. The existence of genetically distinct domesticated strains of chinook salmon in coastal BC generates the need for an evaluation of potential genetic interactions between domesticated escapees and natural spawning populations.     

Characterization of 20 highly variable tetranucleotide microsatellite loci for bull trout (Salvelinus confluentus) and cross‐amplification in other Salvelinus species

We describe the isolation and development of 20 polymorphic tetranucleotide microsatellite loci for bull trout (Salvelinus confluentus). In a sample of 76 fish, we observed between four and 23 alleles per locus. These loci show greater levels of genetic diversity than other loci previously examined in bull trout. Seventeen of these loci were also polymorphic in at least one of the following species: Salvelinus alpinus, Salvelinus malma, Salvelinus namaycush and Salvelinus fontinalis. These loci will aid in our understanding of the molecular ecology of bull trout and other char.     

A Comparative Study of the Biodesulfurization Efficiency of Acidithiobacillus ferrooxidans LY01 Cells Domesticated with Ferrous Iron and Pyrite

The high-sulfur coal desulfurization process completed by A. ferrooxidans LY01 cells domesticated with either ferrous iron [Fe(II)] or pyrite (FeS₂) was investigated in this article. The desulfurization rate for 13 d was as high as 67.8% for the LY01 cells domesticated with pyrite but was only 45.6% for the LY01 cells domesticated with Fe(II). Bacterial adsorption experiments indicated that the bacterial adsorption quantity onto the pyrite particles was similar to the desulfurization efficiency. FTIR analysis showed that chemical composition of the two cell types was similar, but the LY01 cells domesticated with pyrite had higher levels of hydrophobic aromatic R-O groups than cells domesticated with Fe(II). The amount of extracellular polymeric substances (EPS) from the pyrite-domesticated LY01 cells was 1820 μg C/10¹⁰ cells, which was five times more than the amount of EPS in the Fe(II)-domesticated cells; the EPS readily bound Fe(III) with a maximum binding capacity of 0.21 mg Fe(III) per mg C EPS. Strains of pyrite-domesticated LY01 with a high amount of Fe(III) in their EPS possess greater oxidation activity than Fe(II)-domesticated strains with fewer Fe(III). These experiments showed the importance of the substrate-specific differences in the oxidative activity of A. ferrooxidans LY01. In addition, this study provides theoretical guidance for the future optimization of the biodesulfurization process.     

Patterns of allozyme variation in cultivated and wild Sorghum bicolor

Summary Patterns of allozyme variation were surveyed in collections of cultivated and wild sorghum from Africa, the Middle East, and Asia. Data for 30 isozyme loci from a total of 2067 plants representing 429 accessions were analyzed. Regional levels of genetic diversity in the cultivars are greater in northern and central Africa compared to southern Africa, the Middle East, or Asia. The spatial distribution of individual alleles at the most variable loci was studied by plotting allele frequencies on geographic maps covering the distribution of sorghum. Generally, many of the alleles with frequencies below 0.25 are localized in specific portions of the range and are commonly present in more than one race in that region. Several alleles occur in both wild and cultivated sorghum of one region and are absent from sorghum elsewhere, suggesting local introgression between the wild and cultivated forms. Although the same most common allele was found in the wild and cultivated gene pools at 29 of the 30 loci, phenetic analyses separated the majority of wild collections from the cultivars, indicating that the two gene pools are distinct. Wild sorghum from northeast and central Africa exhibits greater genetic similarities to the cultivars compared to wild sorghum of northwest or southern Africa. This is consistent with the theory that wild sorghum of northeast-central Africa is ancestral to domesticated sorghum. Wild sorghums of race arundinaceum of northwest Africa and race virgatum from Egypt are shown to be genetically distinct from both other forms of wild sorghum and from the cultivars. Suggestions for genetic conservation are presented in light of these data.     

Identification of genetic determinants of IGF‐1 levels and longevity among mouse inbred strains

The IGF‐1 signaling pathway plays an important role in regulating longevity. To identify the genetic loci and genes that regulate plasma IGF‐1 levels, we intercrossed MRL/MpJ and SM/J, inbred mouse strains that differ in IGF‐1 levels. Quantitative trait loci (QTL) analysis of IGF‐1 levels of these F2 mice detected four QTL on chromosomes (Chrs) 9 (48 Mb), 10 (86 Mb), 15 (18 Mb), and 17 (85 Mb). Haplotype association mapping of IGF‐1 levels in 28 domesticated inbred strains identified three suggestive loci in females on Chrs 2 (13 Mb), 10 (88 Mb), and 17 (28 Mb) and in four males on Chrs 1 (159 Mb), 3 (52 and 58 Mb), and 16 (74 Mb). Except for the QTL on Chr 9 and 16, all loci co‐localized with IGF‐1 QTL previously identified in other mouse crosses. The most significant locus was the QTL on Chr 10, which contains the Igf1 gene and which had a LOD score of 31.8. Haplotype analysis among 28 domesticated inbred strains revealed a major QTL on Chr 10 overlapping with the QTL identified in the F2 mice. This locus showed three major haplotypes; strains with haplotype 1 had significantly lower plasma IGF‐1 and extended longevity (P < 0.05) than strains with haplotype 2 or 3. Bioinformatic analysis, combined with sequencing and expression studies, showed that Igf1 is the most likely QTL gene, but that other genes may also play a role in this strong QTL.     

Hitchhiking: A Comparison of Linkage and Partial Selfing

Genetic hitchhiking occurs when alleles at unselected loci are changed in frequency because of an association with alleles at a selected locus. This association may be mediated either by linkage or partial selfing (inbreeding) and can affect the gene frequency and gametic disequilibrium at the neutral loci. Hitchhiking from partial selfing (unlinked loci) occurs more quickly than linkage hitchhiking and generally has a greater effect. In addition, partial-selfing hitchhiking can cause increases or changes in sign in gametic disequilibrium between neutral loci. The effects of the two types of hitchhiking with different levels of dominance, zygotic frequencies and number of selected loci are also examined. The general conditions for linkage and partial-selfing hitchhiking are outlined and the implications of hitchhiking are discussed for marker or electrophoretic loci.     

Asymmetry of gene flow and differential geographical structure of molecular diversity in wild and domesticated common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.) from Mesoamerica

Using amplified fragment length polymorphisms (AFLPs), we analyzed the genetic structure of wild and domesticated common bean (Phaseolus vulgaris L.) from Mesoamerica at different geographical levels to test the hypothesis of asymmetric gene flow and investigate the origin of weedy populations. We showed both by phenetic and admixture population analyses that gene flow is about three- to four-fold higher from domesticated to wild populations than in the reverse direction. This result, combined with other work, points to a displacement of genetic diversity in wild populations due to gene flow from the domesticated populations. The weedy populations appear to be genetically intermediate between domesticated and wild populations, suggesting that they originated by hybridization between wild and domesticated types rather than by escape from cultivation. In addition, the domesticated bean races were genetically similar confirming a single domestication event for the Mesoamerican gene pool. Finally, the genetic diversity of the domesticated bean population showed a lower level of geographic structure in comparison to that of the wild populations.     

Cloning and characterization of microsatellite loci in channel catfish, Ictalurus punctatus

Short tandem repeat (microsatellite) loci were cloned from the channel catfish, Ictalurus punctatus , genome for use as molecular markers for genetic improvement of this important agricultural species. Plasmid clones containing catfish genomic DNA inserts were identified, by hybridization with tandem repeat DNA probes, and sequenced using automated laser fluorescence. A feral population of catfish displayed levels of heterozygosity greater than 0·7 for 13 of 22 loci and heterozygosity greater than 0·5 for 20 of 22 loci. Allelic polymorphism ranged from three to 17 alleles per locus in the feral population. Populations of domestic, farm-raised catfish and a research strain displayed levels of heterozygosity similar to the feral population. Non-invasive tissue sampling provided abundant material for the polymerase chain reaction-based genotype assay. The microsatellite loci will be useful in the molecular characterization and genetic improvement of channel catfish populations.     

Functional Traits Differ between Cereal Crop Progenitors and Other Wild Grasses Gathered in the Neolithic Fertile Crescent

The reasons why some plant species were selected as crops and others were abandoned during the Neolithic emergence of agriculture are poorly understood. We tested the hypothesis that the traits of Fertile Crescent crop progenitors were advantageous in the fertile, disturbed habitats surrounding early settlements and in cultivated fields. We screened functional traits related to competition and disturbance in a group of grass species that were increasingly exploited by early plant gatherers, and that were later domesticated (crop progenitors); and in a set of grass species for which there is archaeological evidence of gathering, but which were never domesticated (wild species). We hypothesised that crop progenitors would have greater seed mass, growth rate, height and yield than wild species, as these traits are indicative of greater competitive ability, and that crop progenitors would be more resilient to defoliation. Our results show that crop progenitors have larger seed mass than wild species, germinate faster and have greater seedling size. Increased seed size is weakly but positively correlated with a higher growth rate, which is primarily driven by greater biomass assimilation per unit leaf area. Crop progenitors also tend to have a taller stature, greater grain yield and higher resilience to defoliation. Collectively, the data are consistent with the hypothesis that adaptations to competition and disturbance gave crop progenitors a selective advantage in the areas surrounding early human settlements and in cultivated environments, leading to their adoption as crops through processes of unconscious selection.     

Genetic structure and diversity of natural and domesticated populations of Citrus medica L. in the Eastern Himalayan region of Northeast India

Citron (Citrus medica L.) is a medicinally important species of citrus native to India and occurs in natural forests and home gardens in the foothills of the eastern Himalayan region of northeast India. The wild populations of citron in the region have undergone rapid decline due to natural and anthropogenic disturbances and most of the remaining individuals of citron are found in fragmented natural forests and home gardens in the region. In order to assess the genetic structure and diversity of citron in wild and domesticated populations, we analyzed 219 individuals of C. medica collected from four wild and eight domesticated populations using microsatellite markers. The genetic analysis based on five polymorphic microsatellite loci revealed an average of 13.40 allele per locus. The mean observed and expected heterozygosity values ranged between 0.220–0.540 and 0.438–0.733 respectively among the wild and domesticated populations. Domesticated populations showed close genetic relationships as compared to wild populations and pairwise Nei's genetic distance ranged from 0.062 to 2.091 among wild and domesticated populations. Analysis of molecular variance (AMOVA) showed higher genetic diversity among‐ than within populations. The analysis of population structure revealed five groups. Mixed ancestry of few individuals of different populations revealed exchange of genetic materials among farmers in the region. Citron populations in the region show high genetic variation. The knowledge gained through this study is invaluable for devising genetically sound strategies for conservation of citron genetic resources in the region.     

Microsatellite loci in Bactris gasipaes (Arecaceae): their isolation and characterization

We constructed a microsatellite‐enriched genomic library for Bactris gasipaes, an economically important, domesticated palm. We developed 18 polymorphic microsatellite markers, and found an average of seven alleles per locus in a sample of 14 individuals selected from a germplasm bank. Cross‐species amplification was evaluated in six other Bactris species. The loci detected will permit population studies and germplasm characterization, and can be used for genetic analyses in related species.     

The Rosy Region of Drosophila Melanogaster and Drosophila Simulans. I. Contrasting Levels of Naturally Occurring DNA Restriction Map Variation and Divergence

A 40-kb region around the rosy and snake loci was analyzed for restriction map variation among 60 lines of Drosophila melanogaster and 30 lines of Drosophila simulans collected together at a single locality in Raleigh, North Carolina. DNA sequence variation in D. simulans was estimated to be 6.3 times greater than in D. melanogaster (heterozygosities per nucleotide of 1.9% vs. 0.3%). This result stands in marked contrast to results of studies of phenotypic variation including proteins (allozymes), morphology and chromosome arrangements which are generally less variable and less geographically differentiated in D. simulans. Intraspecific polymorphism is not distributed uniformly over the 40-kb region. The level of heterozygosity per nucleotide varies more than 12-fold across the region in D. simulans, being highest over the hsc2 gene. Similar, though less extreme, variation in heterozygosity is also observed in D. melanogaster. Average interspecific divergence (corrected for intraspecific polymorphism) averaged 3.8%. The pattern of interspecific divergence over the 40-kb region shows some disparities with the spatial distribution of intraspecific variation, but is generally consistent with selective neutrality predictions: the most polymorphic regions within species are generally the most divergent between species. Sequence-length polymorphism is observed for D. melanogaster to be at levels comparable to other gene regions in this species. In contrast, no sequence length variation was observed among D. simulans chromosomes (limit of resolution approximately 100 bp). These data indicate that transposable elements play at best a minor role in the generation of naturally occurring genetic variation in D. simulans compared to D. melanogaster. We hypothesize that differences in species effective population size are the major determinant of the contrasting levels and patterns of DNA sequence and insertion/deletion variation that we report here and the patterns of allozyme and morphological variation and differentiation reported by other workers for these two species.     

Structure of genetic diversity in the two major gene pools of common bean (Phaseolus vulgaris L., Fabaceae)

Domesticated materials with well-known wild relatives provide an experimental system to reveal how human selection during cultivation affects genetic composition and adaptation to novel environments. In this paper, our goal was to elucidate how two geographically distinct domestication events modified the structure and level of genetic diversity in common bean. Specifically, we analyzed the genome-wide genetic composition at 26, mostly unlinked microsatellite loci in 349 accessions of wild and domesticated common bean from the Andean and Mesoamerican gene pools. Using a model-based approach, implemented in the software STRUCTURE, we identified nine wild or domesticated populations in common bean, including four of Andean and four of Mesoamerican origins. The ninth population was the putative wild ancestor of the species, which was classified as a Mesoamerican population. A neighbor-joining analysis and a principal coordinate analysis confirmed genetic relationships among accessions and populations observed with the STRUCTURE analysis. Geographic and genetic distances in wild populations were congruent with the exception of a few putative hybrids identified in this study, suggesting a predominant effect of isolation by distance. Domesticated common bean populations possessed lower genetic diversity, higher F _ST, and generally higher linkage disequilibrium (LD) than wild populations in both gene pools; their geographic distributions were less correlated with genetic distance, probably reflecting seed-based gene flow after domestication. The LD was reduced when analyzed in separate Andean and Mesoamerican germplasm samples. The Andean domesticated race Nueva Granada had the highest F _ST value and widest geographic distribution compared to other domesticated races, suggesting a very recent origin or a selection event, presumably associated with a determinate growth habit, which predominates in this race. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evaluating the utility of microsatellites for investigations of autopolyploid taxa

Autopolyploid taxa present numerous challenges for population genetic analyses due to difficulties determining allele dosage. Dosage ambiguity hinders accurate assessment of allele frequencies, multilocus genotypes (MLGTs), as well as levels and patterns of clonality. The pervasiveness of polyploidy in the evolutionary history of plant taxa makes this a recurring problem. Whereas diploidization of loci may occur over time, duplication of at least some loci is still frequently evident. Fortunately, with high-quality allozyme gels, it is possible to accurately infer allele dosage and, thus, determine exact MLGTs. However, accurately assessing dosage of microsatellite peaks is nearly impossible when studying wild populations with a large number of alleles per locus. Even if precise knowledge of genotypes is not required, for comparable numbers of alleles per locus and loci, the number of "phenotypes" is always lower with microsatellites than allozymes due to the inability to assess allele dosage. Microsatellite loci typically have more alleles per locus relative to allozymes although fewer loci are generally employed. Here, we present a mathematical model for comparing the relative utility of simple sequence repeat (SSR) versus allozyme markers to discriminate MLGTs. For example, the average plant allozyme study (2.6 alleles per locus, 10 polymorphic loci) has better discriminating power than SSR markers with 10 alleles at each of 3 loci, 9 alleles at 4 loci, 6 alleles at 5 loci, 5 alleles at 6 loci, and 4 alleles at 8 loci, demonstrating the value of assessing the relative discriminating power of these markers.     

Dominant vs. codominant genetic markers in the estimation of male mating success

Although a major component of fitness, male reproductive success is generally extremely difficult to estimate. As a result, genetic methods and maximum likelihood models have been developed to estimate male parentage, but all are limited in practice by the degree of genetic variation observable. Scoring individuals phenotypically at a large number of random loci exhibiting dominance (e.g. RAPD markers) may provide a means of detecting sufficient genetic variation. Dominance, however, represents a loss of information and therefore greater variation in the estimate of paternity. A mixture model describing mating in a population is presented to quantify the trade-off between marker types when estimates of male fertility are sought. A sample size 1.5-2.0 times greater is required for dominant markers under some conditions to obtain the same confidence in fertility estimates as for codominant markers, although with large sample sizes the fertility estimates are similar for either marker type. Since the number of dominant DN A markers is not limited in the same manner as is the number of codominant protein markers, one's confidence in the estimates can be increased above that possible from proteins by surveying additional loci. However, for a fixed sample size a trade-off exists between the number of progeny assayed per female and the number of loci surveyed. In many cases more progeny per female provide better estimates of fertility than more loci.     

The impact of microsatellite electromorph size homoplasy on multilocus population structure estimates in a tropical tree (Corythophora alta) and an anadromous fish (Morone saxatilis)

Microsatellite allelic states are determined by electrophoretic sizing of polymerase chain reaction fragments to define electromorphs. Numerous studies have documented that identical microsatellite electromorphs are potentially heterogeneous at the DNA sequence level, a phenomenon called electromorph size homoplasy. Few studies have examined the impact of electromorph size homoplasy on estimates of population genetic parameters. We investigated the frequency of microsatellite electromorph size homoplasy for 12 loci in the tropical tree Corythophora alta and 11 loci in the anadromous fish Morone saxatilis by sequencing 14-23 homozygotes per locus sampled from multiple populations for a total of 453 sequences. Sequencing revealed no homoplasy for M. saxatilis loci. Seven C. alta loci exhibited homoplasy, including single and compound repeat motifs both with and without interruptions. Between 12.5 and 42.9% of electromorphs sampled per locus showed size homoplasy. Two methods of correction for homoplasy in C. alta generally produced little or no change in single-locus estimates of RST, except for two loci in which some additional differentiation among populations was revealed. Twelve-locus estimates of RST (including the seven loci corrected for homoplasy) were slightly greater than estimates from uncorrected data, although the 95% confidence intervals overlapped. The frequency of methodological errors such as clerical mistakes or sample mislabelling per genotype scored was estimated at 5.4 and 7.3% for C. alta and M. saxatilis, respectively. Simulations showed that the increase in RST produced by homoplasy correction was only slightly larger than variation in RST estimates expected to be caused by methodological errors.     

Evolutionarily advanced ant farmers rear polyploid fungal crops

Innovative evolutionary developments are often related to gene or genome duplications. The crop fungi of attine fungus‐growing ants are suspected to have enhanced genetic variation reminiscent of polyploidy, but this has never been quantified with cytological data and genetic markers. We estimated the number of nuclei per fungal cell for 42 symbionts reared by 14 species of Panamanian fungus‐growing ants. This showed that domesticated symbionts of higher attine ants are polykaryotic with 7–17 nuclei per cell, whereas nonspecialized crops of lower attines are dikaryotic similar to most free‐living basidiomycete fungi. We then investigated how putative higher genetic diversity is distributed across polykaryotic mycelia, using microsatellite loci and evaluating models assuming that all nuclei are either heterogeneously haploid or homogeneously polyploid. Genetic variation in the polykaryotic symbionts of the basal higher attine genera Trachymyrmex and Sericomyrmex was only slightly enhanced, but the evolutionarily derived crop fungi of Atta and Acromyrmex leaf‐cutting ants had much higher genetic variation. Our opposite ploidy models indicated that the symbionts of Trachymyrmex and Sericomyrmex are likely to be lowly and facultatively polyploid (just over two haplotypes on average), whereas Atta and Acromyrmex symbionts are highly and obligatorily polyploid (ca. 5–7 haplotypes on average). This stepwise transition appears analogous to ploidy variation in plants and fungi domesticated by humans and in fungi domesticated by termites and plants, where gene or genome duplications were typically associated with selection for higher productivity, but allopolyploid chimerism was incompatible with sexual reproduction.     

Accounting for epistasis in linkage analysis of general pedigrees

Complex traits are generally believed to be influenced by multiple loci. Identification of loci involved in complex traits is more difficult for interacting than for additive loci. Here we describe an extension of the program lm_twoqtl in the package MORGAN to handle two quantitative trait loci (QTLs) with gene-gene interaction. We investigate whether parametric linkage analysis that accounts for such epistasis improves prospects for linkage detection and accuracy of localization of QTLs. Through use of simulated data we show that analysis that accounts for epistasis provides higher lod scores and better localization than does analysis without epistasis. In addition, we demonstrate that the difference between lod scores in the presence vs. absence of use of an interaction model in analysis is greater in extended than in nuclear pedigrees.