Evolution and genetic population structure of prickly lettuce (Lactuca serriola) and its RGC2 resistance gene cluster

Genomic Signatures of Adaptation Across a Precipitation Gradient From Niche Centre to Niche Edge

Evaluating the potential for species to adapt to changing climates relies on understanding current patterns of adaptive variation and selection, which might vary in intensity across a species' niche, hence affecting our inference of where adaptation might be most important in the future. Here, we investigate the genetic basis of adaptation in Lactuca serriola along a steep precipitation gradient in Israel approaching the species' arid niche limit and use candidate loci to inform predictions of its past and future adaptive evolution. Environmental association analyses combined with generalised dissimilarity models revealed 108 candidate genes showing nonlinear shifts in allele frequencies across the gradient, with 66% of these genes under strong selection near the dry niche edge. We detected selection acting on genes with separate suites of biological functions, specifically related to phenology and responses to environmental stressors, including osmotic stress, at the dry niche edge, and related to biotic interactions and defence closer to the niche centre. The adaptive genetic composition of populations, as inferred through polygenic risk scores, point to intensified selection operating towards the dry niche edge. However, inference of past and future evolutionary change predicts larger adaptive shifts occurring in the mesic part of the range, which is most affected by climate change. Our study reveals that adaptive shifts in response to climate change can be heterogeneous across a species' range and not necessarily strongest near its niche edge.

Amplified Fragment Length Polymorphism: Applications and Recent Developments

AFLP or amplified fragment length polymorphism is a PCR-based molecular technique that uses selective amplification of a subset of digested DNA fragments from any source to generate and compare unique fingerprints of genomes. It is more efficient in terms of time, economy, reproducibility, informativeness, resolution, and sensitivity, compared to other popular DNA markers. Besides, it requires very small quantities of DNA and no prior genome information. This technique is widely used in plants for taxonomy, genetic diversity, phylogenetic analysis, construction of high-resolution genetic maps, and positional cloning of genes, to determine relatedness among cultivars and varietal identity, etc. The review encompasses in detail the various applications of AFLP in plants and the major advantages and disadvantages. The review also considers various modifications of this technique and novel developments in detection of polymorphism. A wet-lab protocol is also provided.

Molecular Biogeography of Prickly Lettuce (Lactuca serriola L.) Shows Traces of Recent Range Expansion

Prickly lettuce (Lactuca serriola L., Asteraceae), a wild relative of cultivated lettuce, is an autogamous species which greatly expanded throughout Western and Northern Europe during the last 2 centuries. Here, we present a large-scale biogeographic genetic analysis performed on a dataset represented by 2622 individuals from 110 wild European populations. Thirty-two maternally inherited chloroplast RFLP-markers and 10 nuclear microsatellite loci were used. Microsatellites revealed low genetic variation and high inbreeding coefficients within populations, as well as strong genetic differentiation between populations, which was in accordance with the autogamous breeding system. Analysis of molecular variance based clustering indicated the presence of 3 population clusters, which showed strong geographical patterns. One cluster occupied United Kingdom and part of Northern Europe, and characterized populations with a single predominant genotype. The second mostly combined populations from Northern Europe, while the third cluster grouped populations particularly from Southern Europe. Kriging of gene diversity for L. serriola corroborated northwards and westwards spread from Central (Eastern) Europe. Significant lower genetic diversity characterized the newly colonized parts of the range compared to the historical ones, confirming the importance of founder effects. Stronger pattern of isolation by distance was assessed in the newly colonized areas than in the historical areas (Mantel’s r = 0.20). In the newly colonized areas, populations at short geographic distances were genetically more similar than those in the historical areas. Our results corroborate the species’ recent and rapid northward and westward colonization from Eastern Europe, as well as a decrease of genetic diversity in recently established populations.

Distribution and haplotype diversity of WKS resistance genes in wild emmer wheat natural populations

The wheat stripe rust resistance gene Yr36 ( WKS1 ) with a unique kinase-START domain architecture is highly conserved in wild emmer wheat natural populations. Wild emmer wheat (Triticum dicoccoides) populations have developed various resistance strategies against the stripe rust pathogen Puccinia striiformis f. sp. tritici (Pst). The wild emmer gene, Yr36 (WKS1), which confers partial resistance to a broad spectrum of Pst races, is composed of a kinase and a START lipid-binding domain, a unique gene architecture found only in the Triticeae tribe. The analysis of 435 wild emmer accessions from a broad range of natural habitats revealed that WKS1 and its paralogue WKS2 are present only in the southern distribution range of wild emmer in the Fertile Crescent, supporting the idea that wheat domestication occurred in the northern populations. An analysis of full-length WKS1 sequence from 54 accessions identified 15 different haplotypes and very low-nucleotide diversity (π = 0.00019). The high level of WKS1 sequence conservation among wild emmer populations is in contrast to the high level of diversity previously observed in NB-LRR genes (e.g., Lr10 and Pm3). This phenomenon may reflect the different resistance mechanisms and different evolutionary pathways that shaped these genes, and may shed light on the evolution of genes that confer partial resistance to stripe rust. Only five WKS1 coding sequence haplotypes were revealed among all tested accessions, encoding four different putative WKS1 proteins (designated P0, P1, P2, and P3). Infection tests showed that P0, P1, and P3 haplotypes display a resistance response, while P2 displayed a susceptible response. These results show that the WKS1 proteins (P0, P1, and P3) can be useful to improve wheat resistance to stripe rust.

Genetic Control of Water and Nitrate Capture and Their Use Efficiency in Lettuce (Lactuca sativa L.)

Robustness in lettuce, defined as the ability to produce stable yields across a wide range of environments, may be associated with below-ground traits such as water and nitrate capture. In lettuce, research on the role of root traits in resource acquisition has been rather limited. Exploring genetic variation for such traits and shoot performance in lettuce across environments can contribute to breeding for robustness. A population of 142 lettuce cultivars was evaluated during two seasons (spring and summer) in two different locations under organic cropping conditions, and water and nitrate capture below-ground and accumulation in the shoots were assessed at two sampling dates. Resource capture in each soil layer was measured using a volumetric method based on fresh and dry weight difference in the soil for soil moisture, and using an ion-specific electrode for nitrate. We used these results to carry out an association mapping study based on 1170 single nucleotide polymorphism markers. We demonstrated that our indirect, high-throughput phenotyping methodology was reliable and capable of quantifying genetic variation in resource capture. QTLs for below-ground traits were not detected at early sampling. Significant marker-trait associations were detected across trials for below-ground and shoot traits, in number and position varying with trial, highlighting the importance of the growing environment on the expression of the traits measured. The difficulty of identifying general patterns in the expression of the QTLs for below-ground traits across different environments calls for a more in-depth analysis of the physiological mechanisms at root level allowing sustained shoot growth.

Genetic and biochemical evaluation of natural rubber from Eastern Washington prickly lettuce (Lactuca serriola L.)

Alternative sources of natural rubber are of importance due to economic, biological, and political threats that could diminish supplies of this resource. Prickly lettuce (Lactuca serriola L.) synthesizes long-chain natural rubber and was studied to determine underlying genetic and phenotypic characteristics of rubber biosynthesis. Genotypic and phenotypic analysis of an F2 segregating population using EST-SSR markers led to the discovery of genetic regions linked to natural rubber production. Interval mapping (IM) and multiple QTL mapping (MQM) identified several QTL in the mapping population that had significance based on LOD score thresholds. The discovered QTL and the corresponding local markers are genetic resources for understanding rubber biosynthesis in prickly lettuce and could be used in marker-assisted selection (MAS) breeding. Prickly lettuce is an excellent candidate for elucidating the rubber synthesis mechanism and has potential as a crop plant for rubber production.

Dynamic nucleotide-binding site and leucine-rich repeat-encoding genes in the grass family

The proper use of resistance genes (R genes) requires a comprehensive understanding of their genomics and evolution. We analyzed genes encoding nucleotide-binding sites and leucine-rich repeats in the genomes of rice (Oryza sativa), maize (Zea mays), sorghum (Sorghum bicolor), and Brachypodium distachyon. Frequent deletions and translocations of R genes generated prevalent presence/absence polymorphism between different accessions/species. The deletions were caused by unequal crossover, homologous repair, nonhomologous repair, or other unknown mechanisms. R gene loci identified from different genomes were mapped onto the chromosomes of rice cv Nipponbare using comparative genomics, resulting in an integrated map of 495 R loci. Sequence analysis of R genes from the partially sequenced genomes of an African rice cultivar and 10 wild accessions suggested that there are many additional R gene lineages in the AA genome of Oryza. The R genes with chimeric structures (termed type I R genes) are diverse in different rice accessions but only account for 5.8% of all R genes in the Nipponbare genome. In contrast, the vast majority of R genes in the rice genome are type II R genes, which are highly conserved in different accessions. Surprisingly, pseudogene-causing mutations in some type II lineages are often conserved, indicating that their conservations were not due to their functions. Functional R genes cloned from rice so far have more type II R genes than type I R genes, but type I R genes are predicted to contribute considerable diversity in wild species. Type I R genes tend to reduce the microsynteny of their flanking regions significantly more than type II R genes, and their flanking regions have slightly but significantly lower G/C content than those of type II R genes.

Two loci in sorghum with NB-LRR encoding genes confer resistance to Colletotrichum sublineolum

The aim of this work was to identify plant resistance genes to the sorghum anthracnose fungus Colletotrichum sublineolum. cDNA-AFLP transcript profiling on two contrasting sorghum genotypes inoculated with C. sublineolum generated about 3,000 informative fragments. In a final set of 126 sequenced genes, 15 were identified as biotic stress related. Seven of the plant-derived genes were selected for functional analysis using a Brome mosaic virus-based virus-induced gene silencing (VIGS) system followed by fungal inoculation and quantitative real-time PCR analysis. The candidate set comprised genes encoding resistance proteins (Cs1A, Cs2A), a lipid transfer protein (SbLTP1), a zinc finger-like transcription factor (SbZnTF1), a rice defensin-like homolog (SbDEFL1), a cell death related protein (SbCDL1), and an unknown gene harboring a casein kinase 2-like domain (SbCK2). Our results demonstrate that down-regulation of Cs1A, Cs2A, SbLTP1, SbZnF1 and SbCD1 via VIGS, significantly compromised the resistance response while milder effects were observed with SbDEFL1 and SbCK2. Expanded genome analysis revealed that Cs1A and Cs2A genes are located in two different loci on chromosome 9 closely linked with duplicated genes Cs1B and Cs2B, respectively. The nucleotide binding-leucine rich repeat (NB-LRR) encoding Cs gene sequence information is presently employed in regional breeding programs.

Genome diversity in wild grasses under environmental stress

Patterns of diversity distribution in the Isa defense locus in wild-barley populations suggest adaptive selection at this locus. The extent to which environmental selection may act at additional nuclear-encoded defense loci and within the whole chloroplast genome has now been examined by analyses in two grass species. Analysis of genetic diversity in wild barley (Hordeum spontaneum) defense genes revealed much greater variation in biotic stress-related genes than abiotic stress-related genes. Genetic diversity at the Isa defense locus in wild populations of weeping ricegrass [Microlaena stipoides (Labill.) R. Br.], a very distant wild-rice relative, was more diverse in samples from relatively hotter and drier environments, a phenomenon that reflects observations in wild barley populations. Whole-chloroplast genome sequences of bulked weeping ricegrass individuals sourced from contrasting environments showed higher levels of diversity in the drier environment in both coding and noncoding portions of the genome. Increased genetic diversity may be important in allowing plant populations to adapt to greater environmental variation in warmer and drier climatic conditions.

Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations

INTRODUCTION

Recombination is a key evolutionary factor enhancing diversity. However, the effect of recombination on diversity in inbreeding species is expected to be low. To estimate this effect, recombination and diversity patterns of Lr10 gene were studied in natural populations of the inbreeder species, wild emmer wheat (Triticum dicoccoides). Wild emmer wheat is the progenitor of most cultivated wheats and it harbors rich genetic resources for disease resistance. Lr10 is a leaf rust resistance gene encoding three domains: a coiled-coil, nucleotide-binding site, and leucine-rich repeat (CC-NBS-LRR).

RESULTS

Lr10 was sequenced from 58 accessions representing 12 diverse habitats in Israel. Diversity analysis revealed a high rate of synonymous and non-synonymous substitutions (d (S) = 0.029, d (N) = 0.018, respectively) in the NBS-LRR domains. Moreover, in contrast to other resistance genes, in Lr10 the CC domain was more diverse than the NBS-LRR domains (d (S) = 0.069 vs. 0.029, d (N) = 0.094 vs. 0.018) and was subjected to positive selection in some of the populations. Seventeen recombination events were detected between haplotypes, especially in the CC domain. Linkage disequilibrium (LD) analysis has shown a rapid decay from r (2) = 0.5 to r (2) = 0.1 within a 2-kb span.

CONCLUSION

These results suggest that recombination is a diversifying force for the R-gene, Lr10, in the selfing species T. dicoccoides. This is the first report of a short-range LD decay in wild emmer wheat.

EST-SSR development from 5 Lactuca species and their use in studying genetic diversity among L. serriola biotypes

Prickly lettuce (Lactuca serriola L.) is a problematic weed of Pacific Northwest and recently developed resistance to the auxinic herbicide 2,4-D. There are no publically available simple sequence repeat (SSR) markers to tag 2,4-D resistance genes in L. serriola. Therefore, a study was conducted to develop SSR markers from expressed sequence tags (ESTs) of 5 Lactuca species. A total of 15,970 SSRs were identified among 57,126 EST assemblies belonging to 5 Lactuca species. SSR-containing ESTs (SSR-ESTs) ranged from 6.23% to 7.87%, and SSR densities ranged from 1.28 to 2.51 kb(-1) among the ESTs of 5 Lactuca species. Trinucleotide repeats were the most abundant SSRs detected during the study. As a representative sample, 45 ESTs carrying class I SSRs (≥ 20 nucleotides) were selected for designing primers and were also searched against the dbEST entries for L. sativa and Helianthus annuus (≤ 10(-50); score ≥ 100). In silico analysis of 45 SSR-ESTs showed 82% conservation across species and 68% conservation across genera. Primer pairs synthesized for the above 45 EST-SSRs were used to study genetic diversity among a collection of 22 L. serriola biotypes. Comparison of the resultant dendrogram to that developed using phenotypic evaluation of the same subset of lines showed limited correspondence. Taken together, this study reported a collection of useful SSR markers for L. serriola, confirmed transferability of these markers within and across genera, and demonstrated their usefulness in studying genetic diversity.

The genomic architecture of disease resistance in lettuce

Genbank and The Compositae Genome Project database, containing over 42,000 lettuce unigenes from Lactuca sativa cv. Salinas and L. serriola accession UC96US23 were mined to identify 702 candidate genes involved in pathogen recognition (RGCs), resistance signal transduction, defense responses, and disease susceptibility. In addition, to identify sequences representing additional sub-families of nucleotide binding site (NBS)-leucine-rich repeat encoding genes; the major classes of resistance genes (R-genes), NBS-encoding sequences were amplified by PCR using degenerate oligonucleotides designed to NBS sub-families specific to the subclass Asteridae, which includes the Compositae family. These products were cloned and sequenced resulting in 18 novel NBS sequences from cv. Salinas and 15 novel NBS sequences from UC96US23. Using a variety of marker technologies, 294 of the 735 candidate disease resistance genes were mapped in our primary mapping population, which consisted of 119 F7 recombinant inbred lines derived from an interspecific cross between cv. Salinas and UC96US23. Using markers shared across multiple genetic maps, 36 resistance phenotypic loci, including two new loci for resistance to downy mildew and two quantitative trait loci for resistance to anthracnose were positioned onto the reference map to provide a global view of the genomic architecture of disease resistance in lettuce and to identify candidate genes for resistance phenotypes. The majority but not all of the resistance phenotypes were genetically associated with RGCs.

The genomic architecture of disease resistance in lettuce

Genbank and The Compositae Genome Project database, containing over 42,000 lettuce unigenes from Lactuca sativa cv. Salinas and L. serriola accession UC96US23 were mined to identify 702 candidate genes involved in pathogen recognition (RGCs), resistance signal transduction, defense responses, and disease susceptibility. In addition, to identify sequences representing additional sub-families of nucleotide binding site (NBS)-leucine-rich repeat encoding genes; the major classes of resistance genes (R-genes), NBS-encoding sequences were amplified by PCR using degenerate oligonucleotides designed to NBS sub-families specific to the subclass Asteridae, which includes the Compositae family. These products were cloned and sequenced resulting in 18 novel NBS sequences from cv. Salinas and 15 novel NBS sequences from UC96US23. Using a variety of marker technologies, 294 of the 735 candidate disease resistance genes were mapped in our primary mapping population, which consisted of 119 F7 recombinant inbred lines derived from an interspecific cross between cv. Salinas and UC96US23. Using markers shared across multiple genetic maps, 36 resistance phenotypic loci, including two new loci for resistance to downy mildew and two quantitative trait loci for resistance to anthracnose were positioned onto the reference map to provide a global view of the genomic architecture of disease resistance in lettuce and to identify candidate genes for resistance phenotypes. The majority but not all of the resistance phenotypes were genetically associated with RGCs.

Ecological genomics of natural plant populations: the Israeli perspective

The genomic era revolutionized evolutionary population biology. The ecological genomics of the wild progenitors of wheat and barley reviewed here was central in the research program of the Institute of Evolution, University of Haifa, since 1975 ( http://evolution.haifa.ac.il ). We explored the following questions: (1) How much of the genomic and phenomic diversity of wild progenitors of cultivars (wild emmer wheat, Triticum dicoccoides, the progenitor of most wheat, plus wild relatives of the Aegilops species; wild barley, Hordeum spontaneum, the progenitor of cultivated barley; wild oat, Avena sterilis, the progenitor of cultivated oats; and wild lettuce species, Lactuca, the progenitor and relatives of cultivated lettuce) are adaptive and processed by natural selection at both coding and noncoding genomic regions? (2) What is the origin and evolution of genomic adaptation and speciation processes and their regulation by mutation, recombination, and transposons under spatiotemporal variables and stressful macrogeographic and microgeographic environments? (3) How much genetic resources are harbored in the wild progenitors for crop improvement? We advanced ecological genetics into ecological genomics and analyzed (regionally across Israel and the entire Near East Fertile Crescent and locally at microsites, focusing on the "Evolution Canyon" model) hundreds of populations and thousands of genotypes for protein (allozyme) and deoxyribonucleic acid (DNA) (coding and noncoding) diversity, partly combined with phenotypic diversity. The environmental stresses analyzed included abiotic (climatic and microclimatic, edaphic) and biotic (pathogens, demographic) stresses. Recently, we introduced genetic maps, cloning, and transformation of candidate genes. Our results indicate abundant genotypic and phenotypic diversity in natural plant populations. The organization and evolution of molecular and organismal diversity in plant populations, at all genomic regions and geographical scales, are nonrandom and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection including diversifying, balancing, cyclical, and purifying selection regimes interacting with, but, ultimately, overriding the effects of mutation, migration, and stochasticity. The progenitors of cultivated plants harbor rich genetic resources and are the best hope for crop improvement by both classical and modern biotechnological methods. Future studies should focus on the interplay between structural and functional genome organization focusing on gene regulation.

Conditional neutrality at two adjacent NBS-LRR disease resistance loci in natural populations of Arabidopsis lyrata

We examined patterns of nucleotide diversity at a genomic region containing two linked candidate disease resistance (NBS-LRR) genes in seven populations of the outcrossing plant Arabidopsis lyrata. In comparison with two adjacent control genes and neutral reference genes across the genome, the NBS-LRR genes exhibited elevated nonsynonymous variation and a large number of major-effect polymorphisms causing early stop codons and/or frameshift mutations. In contrast, analysis of synonymous diversity provided no evidence that the region was subject to long-term balancing selection or recent selective sweeps in any of the seven populations surveyed. Also in contrast with earlier surveys of one of these R genes, there was no evidence that the resistance genes or the major-effect mutations were subject to elevated differentiation between populations. We suggest that conditional neutrality in the absence of the corresponding pathogen, rather than long-term balancing selection or local adaptation, may in some circumstances be a significant cause of elevated functional polymorphism at R genes. In contrast with the R genes, analysis of diversity and differentiation at the flanking FERONIA locus showed high population divergence, suggesting local adaptation on this locus controlling male-female signalling during fertilization.