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.

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.

Molecular evidence for adaptive evolution of drought tolerance in wild cereals.

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel.

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

Genome-Wide Association Study in Bread Wheat Identifies Genomic Regions Associated with Grain Yield and Quality under Contrasting Water Availability.

The objectives of this study were to identify genetic loci in the bread wheat genome that would influence yield stability and quality under water stress, and to identify accessions that can be recommended for cultivation in dry and hot regions. We performed a genome-wide association study (GWAS) using a panel of 232 wheat accessions spanning diverse ecogeographic regions. Plants were evaluated in the Israeli Northern Negev, under two environments: water-limited (D; 250 mm) and well-watered (W; 450 mm) conditions; they were genotyped with ~71,500 SNPs derived from exome capture sequencing. Of the 14 phenotypic traits evaluated, 12 had significantly lower values under D compared to W conditions, while the values for two traits were higher under D. High heritability (H2 = 0.5-0.9) was observed for grain yield, spike weight, number of grains per spike, peduncle length, and plant height. Days to heading and grain yield could be partitioned based on accession origins. GWAS identified 154 marker-trait associations (MTAs) for yield and quality-related traits, 82 under D and 72 under W, and identified potential candidate genes. We identified 24 accessions showing high and/or stable yields under D conditions that can be recommended for cultivation in regions under the threat of global climate change.

Variation in phosphorus and sulfur content shapes the genetic architecture and phenotypic associations within the wheat grain ionome.

Dissection of the genetic basis of wheat ionome is crucial for understanding the physiological and biochemical processes underlying mineral accumulation in seeds, as well as for efficient crop breeding. Most of the elements essential for plants are metals stored in seeds as chelate complexes with phytic acid or sulfur-containing compounds. We assume that the involvement of phosphorus and sulfur in metal chelation is the reason for strong phenotypic correlations within ionome. Adjustment of element concentrations for the effect of variation in phosphorus and sulfur seed content resulted in drastic change of phenotypic correlations between the elements. The genetic architecture of wheat grain ionome was characterized by quantitative trait loci (QTL) analysis using a cross between durum and wild emmer wheat. QTL analysis of the adjusted traits and two-trait analysis of the initial traits paired with either P or S considerably improved QTL detection power and accuracy, resulting in the identification of 105 QTLs and 617 QTL effects for 11 elements. Candidate gene search revealed some potential functional associations between QTLs and corresponding genes within their intervals. Thus, we have shown that accounting for variation in P and S is crucial for understanding of the physiological and genetic regulation of mineral composition of wheat grain ionome and can be implemented for other plants.

Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster.

Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).

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.

Genomic Architecture of Phenotypic Plasticity in Response to Water Stress in Tetraploid Wheat.

Phenotypic plasticity is one of the main mechanisms of adaptation to abiotic stresses via changes in critical developmental stages. Altering flowering phenology is a key evolutionary strategy of plant adaptation to abiotic stresses, to achieve the maximum possible reproduction. The current study is the first to apply the linear regression residuals as drought plasticity scores while considering the variation in flowering phenology and traits under non-stress conditions. We characterized the genomic architecture of 17 complex traits and their drought plasticity scores for quantitative trait loci (QTL) mapping, using a mapping population derived from a cross between durum wheat (Triticum turgidum ssp. durum) and wild emmer wheat (T. turgidum ssp. dicoccoides). We identified 79 QTLs affected observed traits and their plasticity scores, of which 33 reflected plasticity in response to water stress and exhibited epistatic interactions and/or pleiotropy between the observed and plasticity traits. Vrn-B3 (TaTF1) residing within an interval of a major drought-escape QTL was proposed as a candidate gene. The favorable alleles for most of the plasticity QTLs were contributed by wild emmer wheat, demonstrating its high potential for wheat improvement. Our study presents a new approach for the quantification of plant adaptation to various stresses and provides new insights into the genetic basis of wheat complex traits under water-deficit stress.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments.

KEY MESSAGE: Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F6 RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

The transcriptome landscape of the carcinogenic treatment response in the blind mole rat: insights into cancer resistance mechanisms.

BACKGROUND: Spalax, the blind mole rat, developed an extraordinary cancer resistance during 40 million years of evolution in a subterranean, hypoxic, thus DNA damaging, habitat. In 50 years of Spalax research, no spontaneous cancer development has been observed. The mechanisms underlying this resistance are still not clarified. We investigated the genetic difference between Spalax and mice that might enable the Spalax relative resistance to cancer development. We compared Spalax and mice responses to a treatment with the carcinogen 3-Methylcholantrene, as a model to assess Spalax' cancer-resistance. RESULTS: We compared RNA-Seq data of untreated Spalax to Spalax with a tumor and identified a high number of differentially expressed genes. We filtered these genes by their expression in tolerant Spalax that resisted the 3MCA, and in mice, and found 25 genes with a consistent expression pattern in the samples susceptible to cancer among species. Contrasting the expressed genes in Spalax with benign granulomas to those in Spalax with malignant fibrosarcomas elucidated significant differences in several pathways, mainly related to the extracellular matrix and the immune system. We found a central cluster of ECM genes that differ greatly between conditions. Further analysis of these genes revealed potential microRNA targets. We also found higher levels of gene expression of some DNA repair pathways in Spalax than in other murines, like the majority of Fanconi Anemia pathway. CONCLUSION: The comparison of the treated with the untreated tissue revealed a regulatory complex that might give an answer how Spalax is able to restrict the tumor growth. By remodeling the extracellular matrix, the possible growth is limited, and the proliferation of cancer cells was potentially prevented. We hypothesize that this regulatory cluster plays a major role in the cancer resistance of Spalax. Furthermore, we identified 25 additional candidate genes that showed a distinct expression pattern in untreated or tolerant Spalax compared to animals that developed a developed either a benign or malignant tumor. While further study is necessary, we believe that these genes may serve as candidate markers in cancer detection.

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity.

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

Features of the organization of bread wheat chromosome 5BS based on physical mapping.

BACKGROUND: The IWGSC strategy for construction of the reference sequence of the bread wheat genome is based on first obtaining physical maps of the individual chromosomes. Our aim is to develop and use the physical map for analysis of the organization of the short arm of wheat chromosome 5B (5BS) which bears a number of agronomically important genes, including genes conferring resistance to fungal diseases. RESULTS: A physical map of the 5BS arm (290 Mbp) was constructed using restriction fingerprinting and LTC software for contig assembly of 43,776 BAC clones. The resulting physical map covered ~ 99% of the 5BS chromosome arm (111 scaffolds, N50 = 3.078 Mb). SSR, ISBP and zipper markers were employed for anchoring the BAC clones, and from these 722 novel markers were developed based on previously obtained data from partial sequencing of 5BS. The markers were mapped using a set of Chinese Spring (CS) deletion lines, and F2 and RICL populations from a cross of CS and CS-5B dicoccoides. Three approaches have been used for anchoring BAC contigs on the 5BS chromosome, including clone-by-clone screening of BACs, GenomeZipper analysis, and comparison of BAC-fingerprints with in silico fingerprinting of 5B pseudomolecules of T. dicoccoides. These approaches allowed us to reach a high level of BAC contig anchoring: 96% of 5BS BAC contigs were located on 5BS. An interesting pattern was revealed in the distribution of contigs along the chromosome. Short contigs (200-999 kb) containing markers for the regions interrupted by tandem repeats, were mainly localized to the 5BS subtelomeric block; whereas the distribution of larger 1000-3500 kb contigs along the chromosome better correlated with the distribution of the regions syntenic to rice, Brachypodium, and sorghum, as detected by the Zipper approach. CONCLUSION: The high fingerprinting quality, LTC software and large number of BAC clones selected by the informative markers in screening of the 43,776 clones allowed us to significantly increase the BAC scaffold length when compared with the published physical maps for other wheat chromosomes. The genetic and bioinformatics resources developed in this study provide new possibilities for exploring chromosome organization and for breeding applications.

Multiple Avirulence Loci and Allele-Specific Effector Recognition Control the Pm3 Race-Specific Resistance of Wheat to Powdery Mildew.

In cereals, several mildew resistance genes occur as large allelic series; for example, in wheat (Triticum aestivum and Triticum turgidum), 17 functional Pm3 alleles confer agronomically important race-specific resistance to powdery mildew (Blumeria graminis). The molecular basis of race specificity has been characterized in wheat, but little is known about the corresponding avirulence genes in powdery mildew. Here, we dissected the genetics of avirulence for six Pm3 alleles and found that three major Avr loci affect avirulence, with a common locus_1 involved in all AvrPm3-Pm3 interactions. We cloned the effector gene AvrPm3(a2/f2) from locus_2, which is recognized by the Pm3a and Pm3f alleles. Induction of a Pm3 allele-dependent hypersensitive response in transient assays in Nicotiana benthamiana and in wheat demonstrated specificity. Gene expression analysis of Bcg1 (encoded by locus_1) and AvrPm3 (a2/f2) revealed significant differences between isolates, indicating that in addition to protein polymorphisms, expression levels play a role in avirulence. We propose a model for race specificity involving three components: an allele-specific avirulence effector, a resistance gene allele, and a pathogen-encoded suppressor of avirulence. Thus, whereas a genetically simple allelic series controls specificity in the plant host, recognition on the pathogen side is more complex, allowing flexible evolutionary responses and adaptation to resistance genes.

Divergence of Drosophila melanogaster repeatomes in response to a sharp microclimate contrast in Evolution Canyon, Israel.

Repeat sequences, especially mobile elements, make up large portions of most eukaryotic genomes and provide enormous, albeit commonly underappreciated, evolutionary potential. We analyzed repeatomes of Drosophila melanogaster that have been diverging in response to a microclimate contrast in Evolution Canyon (Mount Carmel, Israel), a natural evolutionary laboratory with two abutting slopes at an average distance of only 200 m, which pose a constant ecological challenge to their local biotas. Flies inhabiting the colder and more humid north-facing slope carried about 6% more transposable elements than those from the hot and dry south-facing slope, in parallel to a suite of other genetic and phenotypic differences between the two populations. Nearly 50% of all mobile element insertions were slope unique, with many of them disrupting coding sequences of genes critical for cognition, olfaction, and thermotolerance, consistent with the observed patterns of thermotolerance differences and assortative mating.

Adaptive patterns in the p53 protein sequence of the hypoxia- and cancer-tolerant blind mole rat Spalax.

BACKGROUND: The subterranean blind mole rat, Spalax (genus Nannospalax) endures extreme hypoxic conditions and fluctuations in oxygen levels that threaten DNA integrity. Nevertheless, Spalax is long-lived, does not develop spontaneous cancer, and exhibits an outstanding resistance to carcinogenesis in vivo, as well as anti-cancer capabilities in vitro. We hypothesized that adaptations to similar extreme environmental conditions involve common mechanisms for overcoming stress-induced DNA damage. Therefore, we aimed to identify shared features among species that are adapted to hypoxic stress in the sequence of the tumor-suppressor protein p53, a master regulator of the DNA-damage response (DDR). RESULTS: We found that the sequences of p53 transactivation subdomain 2 (TAD2) and tetramerization and regulatory domains (TD and RD) are more similar among hypoxia-tolerant species than expected from phylogeny. Specific positions in these domains composed patterns that are more frequent in hypoxia-tolerant species and have proven to be good predictors of species' classification into stress-related categories. Some of these positions, which are known to be involved in the interactions between p53 and critical DDR proteins, were identified as positively selected. By 3D modeling of p53 interactions with the coactivator p300 and the DNA repair protein RPA70, we demonstrated that, compared to humans, these substitutions potentially reduce the binding of these proteins to Spalax p53. CONCLUSIONS: We conclude that extreme hypoxic conditions may have led to convergent evolutionary adaptations of the DDR via TAD2 and TD/RD domains of p53.

Rapid genomic changes in Drosophila melanogaster adapting to desiccation stress in an experimental evolution system.

BACKGROUND: Experimental evolution studies, coupled with whole genome resequencing and advances in bioinformatics, have become a powerful tool for exploring how populations respond to selection at the genome-wide level, complementary to genome-wide association studies (GWASs) and linkage mapping experiments as strategies to connect genotype and phenotype. In this experiment, we analyzed genomes of Drosophila melanogaster from lines evolving under long-term directional selection for increased desiccation resistance in comparison with control (no-selection) lines. RESULTS: We demonstrate that adaptive responses to desiccation stress have exerted extensive footprints on the genomes, manifested through a high degree of fixation of alleles in surrounding neighborhoods of eroded heterozygosity. These patterns were highly convergent across replicates, consistent with signatures of 'soft' selective sweeps, where multiple alleles present as standing genetic variation become beneficial and sweep through the replicate populations at the same time. Albeit much less frequent, we also observed line-unique sweep regions with zero or near-zero heterozygosity, consistent with classic, or 'hard', sweeps, where novel rather than pre-existing adaptive mutations may have been driven to fixation. Genes responsible for cuticle and protein deubiquitination seemed to be central to these selective sweeps. High divergence within coding sequences between selected and control lines was also reflected by significant results of the McDonald-Kreitman and Ka/Ks tests, showing that as many as 347 genes may have been under positive selection. CONCLUSIONS: Desiccation stress, a common challenge to many organisms inhabiting dry environments, proves to be a very potent selecting factor having a big impact on genome diversity.

Genome differentiation of Drosophila melanogaster from a microclimate contrast in Evolution Canyon, Israel.

The opposite slopes of "Evolution Canyon" in Israel have served as a natural model system of adaptation to a microclimate contrast. Long-term studies of Drosophila melanogaster populations inhabiting the canyon have exhibited significant interslope divergence in thermal and drought stress resistance, candidate genes, mobile elements, habitat choice, mating discrimination, and wing-shape variation, all despite close physical proximity of the contrasting habitats, as well as substantial interslope migration. To examine patterns of genetic differentiation at the genome-wide level, we used high coverage sequencing of the flies' genomes. A total of 572 genes were significantly different in allele frequency between the slopes, 106 out of which were associated with 74 significantly overrepresented gene ontology (GO) terms, particularly so with response to stimulus and developmental and reproductive processes, thus corroborating previous observations of interslope divergence in stress response, life history, and mating functions. There were at least 37 chromosomal "islands" of interslope divergence and low sequence polymorphism, plausible signatures of selective sweeps, more abundant in flies derived from one (north-facing) of the slopes. Positive correlation between local recombination rate and the level of nucleotide polymorphism was also found.

Experimental evolution of recombination and crossover interference in Drosophila caused by directional selection for stress-related traits.

BACKGROUND: Population genetics predicts that tight linkage between new and/or pre-existing beneficial and deleterious alleles should decrease the efficiency of natural selection in finite populations. By decoupling beneficial and deleterious alleles and facilitating the combination of beneficial alleles, recombination accelerates the formation of high-fitness genotypes. This may impose indirect selection for increased recombination. Despite the progress in theoretical understanding, interplay between recombination and selection remains a controversial issue in evolutionary biology. Even less satisfactory is the situation with crossover interference, which is a deviation of double-crossover frequency in a pair of adjacent intervals from the product of recombination rates in the two intervals expected on the assumption of crossover independence. Here, we report substantial changes in recombination and interference in three long-term directional selection experiments with Drosophila melanogaster: for desiccation (~50 generations), hypoxia, and hyperoxia tolerance (>200 generations each). RESULTS: For all three experiments, we found a high interval-specific increase of recombination frequencies in selection lines (up to 40-50% per interval) compared to the control lines. We also discovered a profound effect of selection on interference as expressed by an increased frequency of double crossovers in selection lines. Our results show that changes in interference are not necessarily coupled with increased recombination. CONCLUSIONS: Our results support the theoretical predictions that adaptation to a new environment can promote evolution toward higher recombination. Moreover, this is the first evidence of selection for different recombination-unrelated traits potentially leading, not only to evolution toward increased crossover rates, but also to changes in crossover interference, one of the fundamental features of recombination.

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.