A Numerical Model Supports the Evolutionary Advantage of Recombination Plasticity in Shifting Environments.

AbstractNumerous empirical studies have witnessed an increase in meiotic recombination rate in response to physiological stress imposed by unfavorable environmental conditions. Thus, inherited plasticity in recombination rate is hypothesized to be evolutionarily advantageous in changing environments. Previous theoretical models proceeded from the assumption that organisms increase their recombination rate when the environment becomes more stressful and demonstrated the evolutionary advantage of such a form of plasticity. Here, we numerically explore a complementary scenario-when the plastic increase in recombination rate is triggered by the environmental shifts. Specifically, we assume increased recombination in individuals developing in a different environment than their parents and, optionally, also in offspring of such individuals. We show that such shift-inducible recombination is always superior when the optimal constant recombination implies an intermediate rate. Moreover, under certain conditions, plastic recombination may also appear beneficial when the optimal constant recombination is either zero or free. The advantage of plastic recombination was better predicted by the range of the population's mean fitness over the period of environmental fluctuations, compared with the geometric mean fitness. These results hold for both panmixia and partial selfing, with faster dynamics of recombination modifier alleles under selfing. We think that recombination plasticity can be acquired under the control of environmentally responsive mechanisms, such as chromatin epigenetics remodeling.

Modeling the evolution of recombination plasticity: A prospective review.

Meiotic recombination is one of the main sources of genetic variation, a fundamental factor in the evolutionary adaptation of sexual eukaryotes. Yet, the role of variation in recombination rate and other recombination features remains underexplored. In this review, we focus on the sensitivity of recombination rates to different extrinsic and intrinsic factors. We briefly present the empirical evidence for recombination plasticity in response to environmental perturbations and/or poor genetic background and discuss theoretical models developed to explain how such plasticity could have evolved and how it can affect important population characteristics. We highlight a gap between the evidence, which comes mostly from experiments with diploids, and theory, which typically assumes haploid selection. Finally, we formulate open questions whose solving would help to outline conditions favoring recombination plasticity. This will contribute to answering the long-standing question of why sexual recombination exists despite its costs, since plastic recombination may be evolutionary advantageous even in selection regimes rejecting any non-zero constant recombination.

Genetic architecture of rust resistance in a wheat (Triticum turgidum) diversity panel.

Introduction: Wheat rust diseases are widespread and affect all wheat growing areas around the globe. Breeding strategies focus on incorporating genetic disease resistance. However, pathogens can quickly evolve and overcome the resistance genes deployed in commercial cultivars, creating a constant need for identifying new sources of resistance. Methods: We have assembled a diverse tetraploid wheat panel comprised of 447 accessions of three Triticum turgidum subspecies and performed a genome-wide association study (GWAS) for resistance to wheat stem, stripe, and leaf rusts. The panel was genotyped with the 90K Wheat iSelect single nucleotide polymorphism (SNP) array and subsequent filtering resulted in a set of 6,410 non-redundant SNP markers with known physical positions. Results: Population structure and phylogenetic analyses revealed that the diversity panel could be divided into three subpopulations based on phylogenetic/geographic relatedness. Marker-trait associations (MTAs) were detected for two stem rust, two stripe rust and one leaf rust resistance loci. Of them, three MTAs coincide with the known rust resistance genes Sr13, Yr15 and Yr67, while the other two may harbor undescribed resistance genes. Discussion: The tetraploid wheat diversity panel, developed and characterized herein, captures wide geographic origins, genetic diversity, and evolutionary history since domestication making it a useful community resource for mapping of other agronomically important traits and for conducting evolutionary studies.

The genomics of linkage drag in inbred lines of sunflower.

Crop wild relatives represent valuable sources of alleles for crop improvement, including adaptation to climate change and emerging diseases. However, introgressions from wild relatives might have deleterious effects on desirable traits, including yield, due to linkage drag. Here, we analyzed the genomic and phenotypic impacts of wild introgressions in inbred lines of cultivated sunflower to estimate the impacts of linkage drag. First, we generated reference sequences for seven cultivated and one wild sunflower genotype, as well as improved assemblies for two additional cultivars. Next, relying on previously generated sequences from wild donor species, we identified introgressions in the cultivated reference sequences, as well as the sequence and structural variants they contain. We then used a ridge-regression best linear unbiased prediction (BLUP) model to test the effects of the introgressions on phenotypic traits in the cultivated sunflower association mapping population. We found that introgression has introduced substantial sequence and structural variation into the cultivated sunflower gene pool, including >3,000 new genes. While introgressions reduced genetic load at protein-coding sequences, they mostly had negative impacts on yield and quality traits. Introgressions found at high frequency in the cultivated gene pool had larger effects than low-frequency introgressions, suggesting that the former likely were targeted by artificial selection. Also, introgressions from more distantly related species were more likely to be maladaptive than those from the wild progenitor of cultivated sunflower. Thus, breeding efforts should focus, as far as possible, on closely related and fully compatible wild relatives.

Evolutionary Conservation and Transcriptome Analyses Attribute Perenniality and Flowering to Day-Length Responsive Genes in Bulbous Barley (Hordeum bulbosum).

Rapid population growth and dramatic climatic turnovers are challenging global crop production. These challenges are spurring plant breeders to enhance adaptation and sustainability of major crops. One intriguing approach is to turn annual systems into perennial ones, yet long-term classical breeding efforts to induce perenniality have achieved limited success. Here, we report the results of our investigation of the genetic basis of bulb formation in the nonmodel organism Hordeum bulbosum, a perennial species closely related to barley. To identify candidate genes that regulate bulb formation in H. bulbosum, we applied two complementary approaches. First, we explored the evolutionary conservation of expressed genes among annual Poaceae species. Next, we assembled a reference transcriptome for H. bulbosum and conducted a differential expression (DE) analysis before and after stimulating bulb initiation. Low conservation was identified in genes related to perenniality in H. bulbosum compared with other species, including bulb development and sugar accumulation genes. We also inspected these genes using a DE analysis, which enabled identification of additional genes responsible for bulb initiation and flowering regulation. We propose a molecular model for the regulation of bulb formation involving storage organ development and starch biosynthesis genes. The high conservation observed along a major part of the pathway between H. bulbosum and barley suggests a potential for the application of biotechnological techniques to accelerate breeding toward perenniality in barley.

Emmer Wheat Eco-Geographic and Genomic Congruence Shapes Phenotypic Performance under Mediterranean Climate.

Emmer wheat (Triticum turgidum ssp. dicoccum) is one of the world's oldest domesticated crops, and it harbors a potentially rich reservoir of agronomic and nutritional quality trait variations. The growing global demand for plant-based health-food niche markets has promoted new commercial interest in ancient grains, including Emmer wheat. Although T. dicoccum can also perform well under harsh environments, its cultivation along the Mediterranean agro-ecosystems is sparse. Here, we analyze a unique tetraploid wheat collection (n = 121) representing a wide geographic range of Emmer accessions, using 9897 DArTseq markers and on-field phenotypic characterization to quantify the extent of diversity among populations and the interactions between eco-geographic, genetic, and phenotypic attributes. Population genomic inferences based on the DArTseq data indicated that the collection could be split into four distinguished clusters in accordance with their eco-geographic origin although significant phenotypic variation was observed within clusters. Superior early vegetative vigor, shorter plant height, and early phenology were observed among emmer wheat accessions from Ethiopia compared to accessions from northern regions. This adaptive advantage highlights the potential of emmer wheat as an exotic germplasm for wheat improvement through breeding. The direct integration of such germplasm into conventional or organic farming agro-systems under the Mediterranean basin climate is also discussed.

Are we there yet? Driving the road to evolutionary graph-pangenomics.

With increase in the number of sequenced genomes, it is now recognized that graph-based pangenomes can provide a comprehensive platform to study diversity in a population or species, from point mutations to large chromosomal rearrangements. By incorporating concepts from graph theory, a graph-pangenome can be studied directly to identify genomic regions and genes that underlie important evolutionary processes and traits. Here, I discuss how basic concepts in graph theory can be implemented to address questions in evolutionary genomics and guide future breeding efforts. Despite its compelling versatility, a graph-pangenome assembly is still challenging especially in species with large complex genomes. As technology is rapidly improving, the graph-pangenome is expected to become a central platform in genomics studies and applications. Thus, development of tools and methods that exploit the graph structure are urged to pave the route to evolutionary graph-pangenomics.

Wild sunflower goes viral: Citizen science and comparative genomics allow tracking the origin and establishment of invasive sunflower in the Levant.

Globalization and intensified volume of trade and transport around the world are accelerating the rate of biological invasions. It is therefore increasingly important to understand the processes through which invasive species colonize new habitats, often to the detriment of native flora. The initial steps of an invasion are particularly critical, as the introduced species relies on limited genetic diversity to adapt to a new environment. However, our understanding of this critical stage of the invasion is currently limited. We used a citizen science approach and social media to survey the distribution of invasive sunflower in Israel. We then sampled and sequenced a representative collection and compared it with available genomic data sets of North American wild sunflower, landraces and cultivars. We show that invasive wild sunflower is rapidly establishing throughout Israel, probably from a single, recent introduction from Texas, while maintaining high genetic diversity through ongoing gene flow. Since its introduction, invasive sunflower has spread quickly to most regions, and differentiation was detected despite extensive gene flow between clusters. Our findings suggest that rapid spread followed by continuous gene flow between diverging populations can serve as an efficient mechanism for maintaining sufficient genetic diversity at the early stages of invasion, promoting rapid adaptation and establishment in the new territory.

Biological and Transcriptomic Characterization of Pre-Haustorial Resistance to Sunflower Broomrape (Orobanche cumana W.) in Sunflowers (Helianthus annuus).

Infestations with sunflower broomrape (Orobanche cumana Wallr.), an obligatory root parasite, constitute a major limitation to sunflower production in many regions around the world. Breeding for resistance is the most effective approach to reduce sunflower broomrape infestation, yet resistance mechanisms are often broken by new races of the pathogen. Elucidating the mechanisms controlling resistance to broomrape at the molecular level is, thus, a desirable way to obtain long-lasting resistance. In this study, we investigated broomrape resistance in a confectionery sunflower cultivar with a robust and long-lasting resistance to sunflower broomrape. Visual screening and histological examination of sunflower roots revealed that penetration of the broomrape haustorium into the sunflower roots was blocked at the cortex, indicating a pre-haustorial mechanism of resistance. A comparative RNA sequencing between broomrape-resistant and -susceptible accessions allowed the identification of genes that were significantly differentially expressed upon broomrape infestation. Among these genes were ß-1,3-endoglucanase, ß-glucanase, and ethylene-responsive transcription factor 4 (ERF4). These genes were previously reported to be pathogenesis-related in other plant species. This transcriptomic investigation, together with the histological examinations, led us to conclude that the resistance mechanism involves the identification of the broomrape and the consequent formation of a physical barrier that prevents the establishment of the broomrape into the sunflower roots.

Fitness dependence preserves selection for recombination across diverse mixed mating strategies.

Meiotic recombination and the factors affecting its rate and fate in nature have inspired many studies in theoretical evolutionary biology. Classical theoretical models have inferred that recombination can be favored under a rather restricted parameter range. Thus, the ubiquity of recombination in nature remains an open question. However, these models assumed constant recombination with an equal rate across all individuals within the population, whereas empirical evidence suggests that recombination may display certain sensitivity to ecological stressors and/or genotype fitness. Models assuming condition-dependent recombination show that such a strategy can often be favored over constant recombination. Moreover, in our recent model with panmictic populations subjected to purifying selection, fitness-dependent recombination was quite often favored even when any constant recombination was rejected. By using numerical modeling, we test whether such a 'recombination-rescuing potential' of fitness dependence holds also beyond panmixia, given the recognized effect of mating strategy on the evolution of recombination. We show that deviations from panmixia generally increase the recombination-rescuing potential of fitness dependence, with the strongest effect under intermediate selfing or high clonality. We find that under partial clonality, the evolutionary advantage of fitness-dependent recombination is determined mostly by selection against heterozygotes and additive-by-additive epistasis, while under partial selfing, additive-by-dominance epistasis is also a driver.

Demographic and ecogeographic factors limit wild grapevine spread at the southern edge of its distribution range.

The spatial distribution of plants is constrained by demographic and ecogeographic factors that determine the range and abundance of the species. Wild grapevine (Vitis vinifera ssp. sylvestris) is distributed from Switzerland in the north to Israel in the south. However, little is known about the ecogeographic constraints of this species and its genetic and phenotypic characteristics, especially at the southern edge of its distribution range in the Levant region. In this study, we explore the population structure of southern Levantine wild grapevines and the correlation between demographic and ecogeographic characteristics. Based on our genetic analysis, the wild grapevine populations in this region can be divided into two major subgroups in accordance with a multivariate spatial and ecogeographical clustering model. The identified subpopulations also differ in morphological traits, mainly leaf hairiness which may imply adaptation to environmental stress. The findings suggest that the Upper Jordan River population was spread to the Sea of Galilee area and that a third smaller subpopulation at the south of the Golan Heights may represent a distinguished gene pool or a recent establishment of a new population. A spatial distribution model indicated that distance to water sources, Normalized difference vegetation index, and precipitation are the main environmental factors constraining V. v. sylvestris distribution at its southern distribution range. These factors in addition to limited gene flow between populations prevent further spread of wild grapevines southwards to semi-arid regions.

Tapping Diversity From the Wild: From Sampling to Implementation.

The diversity observed among crop wild relatives (CWRs) and their ability to flourish in unfavorable and harsh environments have drawn the attention of plant scientists and breeders for many decades. However, it is also recognized that the benefit gained from using CWRs in breeding is a potential rose between thorns of detrimental genetic variation that is linked to the trait of interest. Despite the increased interest in CWRs, little attention was given so far to the statistical, analytical, and technical considerations that should guide the sampling design, the germplasm characterization, and later its implementation in breeding. Here, we review the entire process of sampling and identifying beneficial genetic variation in CWRs and the challenge of using it in breeding. The ability to detect beneficial genetic variation in CWRs is strongly affected by the sampling design which should be adjusted to the spatial and temporal variation of the target species, the trait of interest, and the analytical approach used. Moreover, linkage disequilibrium is a key factor that constrains the resolution of searching for beneficial alleles along the genome, and later, the ability to deplete linked deleterious genetic variation as a consequence of genetic drag. We also discuss how technological advances in genomics, phenomics, biotechnology, and data science can improve the ability to identify beneficial genetic variation in CWRs and to exploit it in strive for higher-yielding and sustainable crops.

Genetic and phenotypic analyses indicate that resistance to flooding stress is uncoupled from performance in cultivated sunflower.

Given the rising risk of extreme weather caused by climate change, enhancement of abiotic stress resistance in crops is increasingly urgent. But will the development of stress-resistant cultivars come at the cost of yield under ideal conditions? We hypothesize that this need not be inevitable, because resistance alleles with minimal pleiotropic costs may evade artificial selection and be retained in crop germplasm. Genome-wide association (GWA) analyses for variation in plant performance and flooding response were conducted in cultivated sunflower, a globally important oilseed. We observed broad variation in flooding responses among genotypes. Flooding resistance was not strongly correlated with performance in control conditions, suggesting no inherent trade-offs. Consistent with this finding, we identified a subset of loci conferring flooding resistance, but lacking antagonistic effects on growth. Genetic diversity loss at candidate genes underlying these loci was significantly less than for other resistance genes during cultivated sunflower evolution. Despite bottlenecks associated with domestication and improvement, low-cost resistance alleles remain within the cultivated sunflower gene pool. Thus, development of cultivars that are both flooding-tolerant and highly productive should be straightforward. Results further indicate that estimates of pleiotropic costs from GWA analyses explain, in part, patterns of diversity loss in crop genomes.

Genomic sequence and copy number evolution during hybrid crop development in sunflowers.

Hybrid crops, an important part of modern agriculture, rely on the development of male and female heterotic gene pools. In sunflowers, heterotic gene pools were developed through the use of crop-wild relatives to produce cytoplasmic male sterile female and branching, fertility restoring male lines. Here, we use genomic data from a diversity panel of male, female, and open-pollinated lines to explore the genetic changes brought during modern improvement. We find the male lines have diverged most from their open-pollinated progenitors and that genetic differentiation is concentrated in chromosomes, 8, 10 and 13, due to introgressions from wild relatives. Ancestral variation from open-pollinated varieties almost universally evolved in parallel for both male and female lines suggesting little or no selection for heterotic overdominance. Furthermore, we show that gene content differs between the male and female lines and that differentiation in gene content is concentrated in high FST regions. This means that the introgressions that brought branching and fertility restoration to the male lines, brought with them different gene content from the ancestral haplotypes, including the removal of some genes. Although we find no evidence that gene complementation genomewide is responsible for heterosis between male and female lines, several of the genes that are largely absent in either the male or female lines are associated with pathogen defense, suggesting complementation may be functionally relevant for crop breeders.

Sunflower pan-genome analysis shows that hybridization altered gene content and disease resistance.

Domesticated plants and animals often display dramatic responses to selection, but the origins of the genetic diversity underlying these responses remain poorly understood. Despite domestication and improvement bottlenecks, the cultivated sunflower remains highly variable genetically, possibly due to hybridization with wild relatives. To characterize genetic diversity in the sunflower and to quantify contributions from wild relatives, we sequenced 287 cultivated lines, 17 Native American landraces and 189 wild accessions representing 11 compatible wild species. Cultivar sequences failing to map to the sunflower reference were assembled de novo for each genotype to determine the gene repertoire, or 'pan-genome', of the cultivated sunflower. Assembled genes were then compared to the wild species to estimate origins. Results indicate that the cultivated sunflower pan-genome comprises 61,205 genes, of which 27% vary across genotypes. Approximately 10% of the cultivated sunflower pan-genome is derived through introgression from wild sunflower species, and 1.5% of genes originated solely through introgression. Gene ontology functional analyses further indicate that genes associated with biotic resistance are over-represented among introgressed regions, an observation consistent with breeding records. Analyses of allelic variation associated with downy mildew resistance provide an example in which such introgressions have contributed to resistance to a globally challenging disease.

The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution.

The domesticated sunflower, Helianthus annuus L., is a global oil crop that has promise for climate change adaptation, because it can maintain stable yields across a wide variety of environmental conditions, including drought. Even greater resilience is achievable through the mining of resistance alleles from compatible wild sunflower relatives, including numerous extremophile species. Here we report a high-quality reference for the sunflower genome (3.6 gigabases), together with extensive transcriptomic data from vegetative and floral organs. The genome mostly consists of highly similar, related sequences and required single-molecule real-time sequencing technologies for successful assembly. Genome analyses enabled the reconstruction of the evolutionary history of the Asterids, further establishing the existence of a whole-genome triplication at the base of the Asterids II clade and a sunflower-specific whole-genome duplication around 29 million years ago. An integrative approach combining quantitative genetics, expression and diversity data permitted development of comprehensive gene networks for two major breeding traits, flowering time and oil metabolism, and revealed new candidate genes in these networks. We found that the genomic architecture of flowering time has been shaped by the most recent whole-genome duplication, which suggests that ancient paralogues can remain in the same regulatory networks for dozens of millions of years. This genome represents a cornerstone for future research programs aiming to exploit genetic diversity to improve biotic and abiotic stress resistance and oil production, while also considering agricultural constraints and human nutritional needs.

Genomic analysis of 6,000-year-old cultivated grain illuminates the domestication history of barley.

The cereal grass barley was domesticated about 10,000 years before the present in the Fertile Crescent and became a founder crop of Neolithic agriculture. Here we report the genome sequences of five 6,000-year-old barley grains excavated at a cave in the Judean Desert close to the Dead Sea. Comparison to whole-exome sequence data from a diversity panel of present-day barley accessions showed the close affinity of ancient samples to extant landraces from the Southern Levant and Egypt, consistent with a proposed origin of domesticated barley in the Upper Jordan Valley. Our findings suggest that barley landraces grown in present-day Israel have not experienced major lineage turnover over the past six millennia, although there is evidence for gene flow between cultivated and sympatric wild populations. We demonstrate the usefulness of ancient genomes from desiccated archaeobotanical remains in informing research into the origin, early domestication and subsequent migration of crop species.

Hybridization and extinction.

Hybridization may drive rare taxa to extinction through genetic swamping, where the rare form is replaced by hybrids, or by demographic swamping, where population growth rates are reduced due to the wasteful production of maladaptive hybrids. Conversely, hybridization may rescue the viability of small, inbred populations. Understanding the factors that contribute to destructive versus constructive outcomes of hybridization is key to managing conservation concerns. Here, we survey the literature for studies of hybridization and extinction to identify the ecological, evolutionary, and genetic factors that critically affect extinction risk through hybridization. We find that while extinction risk is highly situation dependent, genetic swamping is much more frequent than demographic swamping. In addition, human involvement is associated with increased risk and high reproductive isolation with reduced risk. Although climate change is predicted to increase the risk of hybridization-induced extinction, we find little empirical support for this prediction. Similarly, theoretical and experimental studies imply that genetic rescue through hybridization may be equally or more probable than demographic swamping, but our literature survey failed to support this claim. We conclude that halting the introduction of hybridization-prone exotics and restoring mature and diverse habitats that are resistant to hybrid establishment should be management priorities.

Why Close a Bacterial Genome? The Plasmid of Alteromonas Macleodii HOT1A3 is a Vector for Inter-Specific Transfer of a Flexible Genomic Island.

Genome sequencing is rapidly becoming a staple technique in environmental and clinical microbiology, yet computational challenges still remain, leading to many draft genomes which are typically fragmented into many contigs. We sequenced and completely assembled the genome of a marine heterotrophic bacterium, Alteromonas macleodii HOT1A3, and compared its full genome to several draft genomes obtained using different reference-based and de novo methods. In general, the de novo assemblies clearly outperformed the reference-based or hybrid ones, covering >99% of the genes and representing essentially all of the gene functions. However, only the fully closed genome (∼4.5 Mbp) allowed us to identify the presence of a large, 148 kbp plasmid, pAM1A3. While HOT1A3 belongs to A. macleodii, typically found in surface waters ("surface ecotype"), this plasmid consists of an almost complete flexible genomic island (fGI), containing many genes involved in metal resistance previously identified in the genomes of Alteromonas mediterranea ("deep ecotype"). Indeed, similar to A. mediterranea, A. macleodii HOT1A3 grows at concentrations of zinc, mercury, and copper that are inhibitory for other A. macleodii strains. The presence of a plasmid encoding almost an entire fGI suggests that wholesale genomic exchange between heterotrophic marine bacteria belonging to related but ecologically different populations is not uncommon.

Transcriptome profiling of wheat glumes in wild emmer, hulled landraces and modern cultivars.

BACKGROUND: Wheat domestication is considered as one of the most important events in the development of human civilization. Wheat spikelets have undergone significant changes during evolution under domestication, resulting in soft glumes and larger kernels that are released easily upon threshing. Our main goal was to explore changes in transcriptome expression in glumes that accompanied wheat evolution under domestication. METHODS: A total of six tetraploid wheat accessions were selected for transcriptome profiling based on their rachis brittleness and glumes toughness. RNA pools from glumes of the central spikelet at heading time were used to construct cDNA libraries for sequencing. The trimmed reads from each library were separately aligned to the reference sub-genomes A and B, which were extracted from wheat survey sequence. Differentially expression analysis and functional annotation were performed between wild and domesticated wheat, to identity candidate genes associated with evolution under domestication. Selected candidate genes were validated using real time PCR. RESULTS: Transcriptome profiles of wild emmer wheat, wheat landraces, and wheat cultivars were compared using next generation sequencing (RNA-seq). We have found a total of 194,893 transcripts, of which 73,150 were shared between wild, landraces, and cultivars. From 781 differentially expressed genes (DEGs), 336 were down-regulated and 445 were up-regulated in the domesticated compared to wild wheat genotypes. Gene Ontology (GO) annotation assigned 293 DEGs (37.5 %) to GO term groups, of which 134 (17.1 %) were down-regulated and 159 (20.4 %) up-regulated in the domesticated wheat. Some of the down-regulated DEGs in domesticated wheat are related to the biosynthetic pathways that eventually define the mechanical strength of the glumes, such as cell wall, lignin, pectin and wax biosynthesis. The reduction in gene expression of such genes, may explain the softness of the glumes in the domesticated forms. In addition, we have identified genes involved in nutrient remobilization that may affect grain size and other agronomic traits evolved under domestication. CONCLUSIONS: The comparison of RNA-seq profiles between glumes of wheat groups differing in glumes toughness and rachis brittleness revealed a few DEGs that may be involved in glumes toughness and nutrient remobilization. These genes may be involved in processes of wheat improvement under domestication.