Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Genetic diversity and population structure assessment of Western Canadian barley cooperative trials

Studying the population structure and genetic diversity of historical datasets is a proposed use for association analysis. This is particularly important when the dataset contains traits that are time-consuming or costly to measure. A set of 96 elite barley genotypes, developed from eight breeding programs of the Western Canadian Cooperative Trials were used in the current study. Genetic diversity, allelic variation, and linkage disequilibrium (LD) were investigated using 5063 high-quality SNP markers via the Illumina 9K Barley Infinium iSelect SNP assay. The distribution of SNPs markers across the barley genome ranged from 449 markers on chromosome 1H to 1111 markers on chromosome 5H. The average polymorphism information content (PIC) per locus was 0.275 and ranged from 0.094 to 0.375. Bayesian clustering in STRUCTURE and principal coordinate analysis revealed that the populations are differentiated primarily due to the different breeding program origins and ear-row type into five subpopulations. Analysis of molecular variance based on PhiPT values suggested that high values of genetic diversity were observed within populations and accounted for 90% of the total variance. Subpopulation 5 exhibited the most diversity with the highest values of the diversity indices, which represent the breeding program gene pool of AFC, AAFRD, AU, and BARI. With increasing genetic distance, the LD values, expressed as r², declined to below the critical r^{2 =} 0.18 after 3.91 cM, and the same pattern was observed on each chromosome. Our results identified an important pattern of genetic diversity among the Canadian barley panel that was proposed to be representative of target breeding programs and may have important implications for association mapping in the future. This highlight, that efforts to identify novel variability underlying this diversity may present practical breeding opportunities to develop new barley genotypes.

The Population Structure of a Globe Artichoke Worldwide Collection, as Revealed by Molecular and Phenotypic Analyzes

The knowledge of the organization of the domesticated gene pool of crop species is an essential requirement to understand crop evolution, to rationalize conservation programs, and to support practical decisions in plant breeding. Here, we integrate simple sequence repeat (SSR) analysis and phenotypic characterization to investigate a globe artichoke collection that comprises most of the varieties cultivated worldwide. We show that the cultivated gene pool of globe artichoke includes five distinct genetic groups associated with the major phenotypic typologies: Catanesi (which based on our analysis corresponds to Violetti di Provenza), Spinosi, Violetti di Toscana, Romaneschi, and Macau. We observed that 17 and 11% of the molecular and phenotypic variance, respectively, is between these groups, while within groups, strong linkage disequilibrium and heterozygote excess are evident. The divergence between groups for quantitative traits correlates with the average broad-sense heritability within the groups. The phenotypic divergence between groups for both qualitative and quantitative traits is strongly and positively correlated with SSR divergence (FST) between groups. All this implies a low population size and strong bottleneck effects, and indicates a long history of clonal propagation and selection during the evolution of the domesticated gene pool of globe artichoke. Moreover, the comparison between molecular and phenotypic population structures suggests that harvest time, plant architecture (i.e., plant height, stem length), leaf spininess, head morphology (i.e., head shape, bract shape, spininess) together with the number of heads per plant were the main targets of selection during the evolution of the cultivated germplasm. We emphasize our findings in light of the potential exploitation of this collection for association mapping studies.

The influence of breeding history, origin and growth type on population structure of barley as revealed by SSR markers

Natural and mass selection during domestication and cultivation favored particular traits of interest in barley. In the present study, population structure, and genetic relationships among 144 accessions of barley landraces and breeding materials from various countries were studied using a set of 77 and 72 EST-SSR and gSSR markers, respectively distributed on seven chromosomes of barley. In total, 262 and 429 alleles were amplified in 77 EST-SSRs and 72 gSSR loci, respectively. Out of which, 185 private/group-specific alleles were identified in the landraces compared with 14 in "cultivar and advanced breeding lines", indicating the possibility to introgress favorite alleles from landraces into breeding materials. Comparative analysis of genetic variation among breeding materials, Iranian landraces, and exotic landraces revealed higher genetic diversity in Iranian landraces compared with others. A total of 37, 15, and 14 private/group-specific alleles were identified in Iranian landraces, exotic landraces, and breeding materials, respectively. The most likely groups for 144 barley genotypes were three as inferred using model- and distance-based clustering as well as principal coordinate analysis which assigned the landraces and breeding materials into separate groups. The distribution of alleles was found to be correlated with population structure, domestication history and eco-geographical factors. The high allelic richness in the studied set of barley genotype provides insights into the available diversity and allows the construction of core groups based on maximizing allelic diversity for use in barley breeding programs.

Environmental Association Identifies Candidates for Tolerance to Low Temperature and Drought

Barley (Hordeum vulgare ssp. vulgare) is cultivated from the equator to the Arctic Circle. The wild progenitor species, Hordeum vulgare ssp. spontaneum, occupies a relatively narrow latitudinal range (∼30 - 40° N) primarily at low elevation (< 1,500 m). Adaptation to the range of cultivation has occurred over ∼8,000 years. The genetic basis of adaptation is amenable to study through environmental association. An advantage of environmental association in a well-characterized crop is that many loci that contribute to climatic adaptation and abiotic stress tolerance have already been identified. This provides the opportunity to determine if environmental association approaches effectively identify these loci of large effect. Using published genotyping from 7,864 SNPs in 803 barley landraces, we examined allele frequency differentiation across multiple partitions of the data and mixed model associations relative to bioclimatic variables. Using newly generated resequencing data from a subset of these landraces, we tested for linkage disequilibrium (LD) between SNPs queried in genotyping and SNPs in neighboring loci. Six loci previously reported to contribute to adaptive differences in flowering time and abiotic stress in barley and six loci previously identified in other plant species were identified in our analyses. In many cases, patterns of LD are consistent with the causative variant occurring in the immediate vicinity of the queried SNP. The identification of barley orthologs to well-characterized genes may provide a new understanding of the nature of adaptive variation and could permit a more targeted use of potentially adaptive variants in barley breeding and germplasm improvement.

A tiered approach to genome-wide association analysis for the adherence of hulls to the caryopsis of barley seeds reveals footprints of selection

BACKGROUND

Seeds of domesticated barley are grouped into two distinct types, which differ in morphology. Caryopses covered by adaxial (palea) and abaxial (lemma) hulls that tightly adhere to the pericarp at maturity give rise to hulled seeds whereas caryopses without adhering hulls give rise to naked seeds. The naked caryopsis character is an essential trait regarding the end use of barley.

RESULTS

To uncover the genetic basis of the trait, a genome-wide association study (GWAS) has been performed in a panel comprising 222 2-rowed and 303 6-rowed spring barley landrace accessions. In addition to the well-described Nud locus on chromosome 7H, three novel loci showed strong associations with the trait: the first locus on 2H was specifically detected in 6-rowed accessions, the second locus on 3H was found in 2-rowed accessions from Eurasia and the third locus on 6H was revealed in 6-rowed accessions from Ethiopia. PCR analysis of naked accessions also confirmed the absence of a 17 kb region harboring the Nud gene on chromosome 7H for all but one naked accession. The latter was characterized by a slightly variant phenotype of the caryopsis.

CONCLUSION

Our findings provide evidence of the pervasiveness of the 17 kb deletion in spring barley from different geographic regions and at the same time reveal genomic footprints of selection in naked barley, which follow both geographic and morphological patterns.

Spatial genetic structure in wild cardoon, the ancestor of cultivated globe artichoke: Limited gene flow, fragmentation and population history

Nuclear and chloroplast markers and phenotypic characters were integrated to analyse the population genetic structure of wild cardoon, Cynara cardunculus var. sylvestris, the ancestor of cultivated globe artichoke, Cynara cardunculus var. scolymus on the island of Sardinia, Italy. The spatial scale ranged from a few metres to ∼200km. Wild cardoon appears to be genetically fragmented, with significant genetic divergence at various scales, indicating that gene flow is insufficient to counterbalance the effects of genetic drift or founder effects. Divergence between populations was higher for chloroplast (40%) than for nuclear markers (15%), suggesting that gene flow via seed was lower than via pollen. Two main genetic groups were detected; these correlated with differences in flowering time, capitula size, glossiness, and anthocyanin pigmentation. A complex population structure of wild cardoon emerged over small spatial scales, likely resulting from the interplay between gene dispersal, colonisation history and selective forces. Indeed, Sardinia appears to be a 'hybrid zone' of different gene pools. The island has unique diverse germplasm that has originated from hybridisation among different gene pools. The sampling of seeds from a few plants but from many sites is suggested as the best strategy to harvest the genetic diversity of wild cardoon.

Genetic variation, population structure and linkage disequilibrium in Switchgrass with ISSR, SCoT and EST-SSR markers

Background

To evaluate genetic variation, population structure, and the extent of linkage disequilibrium (LD), 134 switchgrass (Panicum virgatum L.) samples were analyzed with 51 markers, including 16 ISSRs, 20 SCoTs, and 15 EST-SSRs.

Results

In this study, a high level of genetic variation was observed in the switchgrass samples and they had an average Nei’s gene diversity index (H) of 0.311. A total of 793 bands were obtained, of which 708 (89.28 %) were polymorphic. Using a parameter marker index (MI), the efficiency of the three types of markers (ISSR, SCoT, and EST-SSR) in the study were compared and we found that SCoT had a higher marker efficiency than the other two markers. The 134 switchgrass samples could be divided into two sub-populations based on STRUCTURE, UPGMA clustering, and principal coordinate analyses (PCA), and upland and lowland ecotypes could be separated by UPGMA clustering and PCA analyses. Linkage disequilibrium analysis revealed an average r² of 0.035 across all 51 markers, indicating a trend of higher LD in sub-population 2 than that in sub-population 1 (P < 0.01).

Conclusions

The population structure revealed in this study will guide the design of future association studies using these switchgrass samples.

Co-evolution in a landrace meta-population: two closely related pathogens interacting with the same host can lead to different adaptive outcomes

We examined the local adaptation patterns in a system comprising several interconnected heterogeneous plant populations from which populations of two phylogenetically closely related pathogens were also sampled. The host is Hordeum vulgare (cultivated barley); the pathogens are Pyrenophora teres f. teres (net form) and Pyrenophora teres f. maculata (spot form), the causal agents of barley net blotch. We integrated two approaches, the comparison between the population structures of the host and the pathogens, and a cross-inoculation test. We demonstrated that two closely related pathogens with very similar niche specialisation and life-styles can give rise to different co-evolutionary outcomes on the same host. Indeed, we detected local adaptation for the net form of the pathogen but not for the spot form. We also provided evidence that an a-priori well-known resistance quantitative-trait-locus on barley chromosome 6H is involved in the co-evolutionary ‘arms race’ between the plant and the net-form pathogen. Moreover, data suggested latitudinal clines of host resistance and that different ecological conditions can result in differential selective pressures at different sites. Our data are of interest for on-farm conservation of plant genetic resources, as also in establishing efficient breeding programs and strategies for deployment of resistance genes of P. teres.

Genetic diversity and population structure in a legacy collection of spring barley landraces adapted to a wide range of climates

Global environmental change and increasing human population emphasize the urgent need for higher yielding and better adapted crop plants. One strategy to achieve this aim is to exploit the wealth of so called landraces of crop species, representing diverse traditional domesticated populations of locally adapted genotypes. In this study, we investigated a comprehensive set of 1485 spring barley landraces (Lrc1485) adapted to a wide range of climates, which were selected from one of the largest genebanks worldwide. The landraces originated from 5° to 62.5° N and 16° to 71° E. The whole collection was genotyped using 42 SSR markers to assess the genetic diversity and population structure. With an average allelic richness of 5.74 and 372 alleles, Lrc1485 harbours considerably more genetic diversity than the most polymorphic current GWAS panel for barley. Ten major clusters defined most of the population structure based on geographical origin, row type of the ear and caryopsis type – and were assigned to specific climate zones. The legacy core reference set Lrc648 established in this study will provide a long-lasting resource and a very valuable tool for the scientific community. Lrc648 is best suited for multi-environmental field testing to identify candidate genes underlying quantitative traits but also for allele mining approaches.

Farmers without borders-genetic structuring in century old barley (Hordeum vulgare)

The geographic distribution of genetic diversity can reveal the evolutionary history of a species. For crop plants, phylogeographic patterns also indicate how seed has been exchanged and spread in agrarian communities. Such patterns are, however, easily blurred by the intense seed trade, plant improvement and even genebank conservation during the twentieth century, and discerning fine-scale phylogeographic patterns is thus particularly challenging. Using historical crop specimens, these problems are circumvented and we show here how high-throughput genotyping of historical nineteenth century crop specimens can reveal detailed geographic population structure. Thirty-one historical and nine extant accessions of North European landrace barley (Hordeum vulgare L.), in total 231 individuals, were genotyped on a 384 single nucleotide polymorphism assay. The historical material shows constant high levels of within-accession diversity, whereas the extant accessions show more varying levels of diversity and a higher degree of total genotype sharing. Structure, discriminant analysis of principal components and principal component analysis cluster the accessions in latitudinal groups across country borders in Finland, Norway and Sweden. FST statistics indicate strong differentiation between accessions from southern Fennoscandia and accessions from central or northern Fennoscandia, and less differentiation between central and northern accessions. These findings are discussed in the context of contrasting historical records on intense within-country south to north seed movement. Our results suggest that although seeds were traded long distances, long-term cultivation has instead been of locally available, possibly better adapted, genotypes.

Population structure of barley landrace populations and gene-flow with modern varieties

Landraces are heterogeneous plant varieties that are reproduced by farmers as populations that are subject to both artificial and natural selection. Landraces are distinguished by farmers due to their specific traits, and different farmers often grow different populations of the same landrace. We used simple sequence repeats (SSRs) to analyse 12 barley landrace populations from Sardinia from two collections spanning 10 years. We analysed the population structure, and compared the population diversity of the landraces that were collected at field level (population). We used a representative pool of barley varieties for diversity comparisons and to analyse the effects of gene flow from modern varieties. We found that the Sardinian landraces are a distinct gene pool from those of both two-row and six-row barley varieties. There is also a low, but significant, mean level and population-dependent level of introgression from the modern varieties into the Sardinian landraces. Moreover, we show that the Sardinian landraces have the same level of gene diversity as the representative sample of modern commercial varieties grown in Italy in the last decades, even within population level. Thus, these populations represent crucial sources of germplasm that will be useful for crop improvement and for population genomics studies and association mapping, to identify genes, loci and genome regions responsible for adaptive variations. Our data also suggest that landraces are a source of valuable germplasm for sustainable agriculture in the context of future climate change, and that in-situ conservation strategies based on farmer use can preserve the genetic identity of landraces while allowing adaptation to local environments.

Twentieth-century changes in the genetic composition of Swedish field pea metapopulations

Landrace crops are formed by local adaptation, genetic drift and gene flow through seed exchange. In reverse, the study of genetic structure between landrace populations can reveal the effects of these forces over time. We present here the analysis of genetic diversity in 40 Swedish field pea (Pisum sativum L.) populations, either available as historical seed samples from the late nineteenth century or as extant gene bank accessions assembled in the late twentieth century. The historical material shows constant high levels of within-population diversity, whereas the extant accessions show varying, and overall lower, levels of within-population diversity. Structure and principal component analysis cluster most accessions, both extant and historical, in groups after geographical origin. County-wise analyses of the accessions show that the genetic diversity of the historical accessions is largely overlapping. In contrast, most extant accessions show signs of genetic drift. They harbor a subset of the alleles found in the historical accessions and are more differentiated from each other. These results reflect how, historically present metapopulations have been preserved during the twentieth century, although as genetically isolated populations.

DNA bending-induced phase transition of encapsidated genome in phage λ

The DNA structure in phage capsids is determined by DNA-DNA interactions and bending energy. The effects of repulsive interactions on DNA interaxial distance were previously investigated, but not the effect of DNA bending on its structure in viral capsids. By varying packaged DNA length and through addition of spermine ions, we transform the interaction energy from net repulsive to net attractive. This allowed us to isolate the effect of bending on the resulting DNA structure. We used single particle cryo-electron microscopy reconstruction analysis to determine the interstrand spacing of double-stranded DNA encapsidated in phage λ capsids. The data reveal that stress and packing defects, both resulting from DNA bending in the capsid, are able to induce a long-range phase transition in the encapsidated DNA genome from a hexagonal to a cholesteric packing structure. This structural observation suggests significant changes in genome fluidity as a result of a phase transition affecting the rates of viral DNA ejection and packaging.