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.

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.

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.

The evolutionary advantage of fitness-dependent recombination in diploids: A deterministic mutation-selection balance model

Recombination's omnipresence in nature is one of the most intriguing problems in evolutionary biology. The question of why recombination exhibits certain general features is no less interesting than that of why it exists at all. One such feature is recombination's fitness dependence (FD). The so far developed population genetics models have focused on the evolution of FD recombination mainly in haploids, although the empirical evidence for this phenomenon comes mostly from diploids. Using numerical analysis of modifier models for infinite panmictic populations, we show here that FD recombination can be evolutionarily advantageous in diploids subjected to purifying selection. We ascribe this advantage to the differential rate of disruption of lower- versus higher-fitness genotypes, which can be manifested in selected systems with at least three loci. We also show that if the modifier is linked to such selected system, it can additionally benefit from modifying this linkage in a fitness-dependent manner. The revealed evolutionary advantage of FD recombination appeared robust to crossover interference within the selected system, either positive or negative. Remarkably, FD recombination was often favored in situations where any constant nonzero recombination was evolutionarily disfavored, implying a relaxation of the rather strict constraints on major parameters (e.g., selection intensity and epistasis) required for the evolutionary advantage of nonzero recombination formulated by classical models.

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.

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

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 F₆ 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.

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

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 F₆ 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.

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.

Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family

Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.

What drives the evolution of condition-dependent recombination in diploids? Some insights from simulation modelling

While the evolutionary advantages of non-zero recombination rates have prompted diverse theoretical explanations, the evolution of essential recombination features remains underexplored. We focused on one such feature, the condition dependence of recombination, viewed as the variation in within-generation sensitivity of recombination to external (environment) and/or internal (genotype) conditions. Limited empirical evidence for its existence comes mainly from diploids, whereas theoretical models show that it only easily evolves in haploids. The evolution of condition-dependent recombination can be explained by its advantage for the selected system (indirect effect), or by benefits to modifier alleles, ensuring this strategy regardless of effects on the selected system (direct effect). We considered infinite panmictic populations of diploids exposed to a cyclical two-state environment. Each organism had three selected loci. Examining allele dynamics at a fourth, selectively neutral recombination modifier locus, we frequently observed that a modifier allele conferring condition-dependent recombination between the selected loci displaced the allele conferring the optimal constant recombination rate. Our simulations also confirm the results of theoretical studies showing that condition-dependent recombination cannot evolve in diploids on the basis of direct fitness-dependent effects alone. Therefore, the evolution of condition-dependent recombination in diploids can be driven by indirect effects alone, i.e. by modifier effects on the selected system.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Regulation of gene expression and RNA editing in Drosophila adapting to divergent microclimates

Determining the mechanisms by which a species adapts to its environment is a key endeavor in the study of evolution. In particular, relatively little is known about how transcriptional processes are fine-tuned to adjust to different environmental conditions. Here we study Drosophila melanogaster from 'Evolution Canyon' in Israel, which consists of two opposing slopes with divergent microclimates. We identify several hundred differentially expressed genes and dozens of differentially edited sites between flies from each slope, correlate these changes with genetic differences, and use CRISPR mutagenesis to validate that an intronic SNP in prominin regulates its editing levels. We also demonstrate that while temperature affects editing levels at more sites than genetic differences, genetically regulated sites tend to be less affected by temperature. This work shows the extent to which gene expression and RNA editing differ between flies from different microclimates, and provides insights into the regulation responsible for these differences.

A High-Density Genetic Map of Wild Emmer Wheat from the Karaca Dağ Region Provides New Evidence on the Structure and Evolution of Wheat Chromosomes

Wild emmer (Triticum turgidum ssp. dicoccoides) is a progenitor of all cultivated wheat grown today. It has been hypothesized that emmer was domesticated in the Karaca Dağ region in southeastern Turkey. A total of 445 recombinant inbred lines of T. turgidum ssp. durum cv. 'Langdon' x wild emmer accession PI 428082 from this region was developed and genotyped with the Illumina 90K single nucleotide polymorphism Infinium assay. A genetic map comprising 2,650 segregating markers was constructed. The order of the segregating markers and an additional 8,264 co-segregating markers in the Aegilops tauschii reference genome sequence was used to compare synteny of the tetraploid wheat with the Brachypodium distachyon, rice, and sorghum. These comparisons revealed the presence of 15 structural chromosome rearrangements, in addition to the already known 4A-5A-7B rearrangements. The most common type was an intra-chromosomal translocation in which the translocated segment was short and was translocated only a short distance along the chromosome. A large reciprocal translocation, one small non-reciprocal translocation, and three large and one small paracentric inversions were also discovered. The use of inversions for a phylogeny reconstruction in the Triticum-Aegilops alliance was illustrated. The genetic map was inconsistent with the current model of evolution of the rearranged chromosomes 4A-5A-7B. Genetic diversity in the rearranged chromosome 4A showed that the rearrangements might have been contemporary with wild emmer speciation. A selective sweep was found in the centromeric region of chromosome 4A in Karaca Dağ wild emmer but not in 4A of T. aestivum. The absence of diversity from a large portion of chromosome 4A of wild emmer, believed to be ancestral to all domesticated wheat, is puzzling.

CYCLICAL ENVIRONMENTAL CHANGES AS A FACTOR MAINTAINING GENETIC POLYMORPHISM. 2. DIPLOID SELECTION FOR AN ADDITIVE TRAIT

The subject of this paper is polymorphism maintenance due to stabilizing selection with a moving optimum. It was shown that in case of two-locus additive control of the selected trait, global polymorphism is possible only when the geometric mean fitnesses of double homozygotes averaged over the period are lower than that of the single heterozygotes and of the double heterozygote (with a multiplier [1 - r]^p , which depends on recombination rate r and period length p). But local stability of polymorphism cannot be excluded even if geometric mean fitnesses of all double homozygotes are higher than that of all heterozygotes. We proved, that for logarithmically convex fitness functions, cyclical changes of the optimum cannot help in polymorphism maintenance in case of additive control of the selected trait by two equal loci. However, within the same class of fitness functions, nonequal gene action and/or dominance effect for one or both loci may lead to local polymorphism stability with large enough polymorphism attracting domain. The higher the intensity of selection and closer the linkage between selected loci the larger is this domain. Note that even simple cyclical selection could result in two forms of polymorphic limiting behavior: (a) usually expected forced cycle with a period equal to that of environmental changes; and (b) "supercycles," nondumping auto-oscillations with a period comprising of hundreds of forced oscillation periods.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0743-8) contains supplementary material, which is available to authorized users.

Genome-Wide Differentiation of Various Melon Horticultural Groups for Use in GWAS for Fruit Firmness and Construction of a High Resolution Genetic Map

Melon (Cucumis melo L.) is a phenotypically diverse eudicot diploid (2n = 2x = 24) has climacteric and non-climacteric morphotypes and show wide variation for fruit firmness, an important trait for transportation and shelf life. We generated 13,789 SNP markers using genotyping-by-sequencing (GBS) and anchored them to chromosomes to understand genome-wide fixation indices (Fst) between various melon morphotypes and genomewide linkage disequilibrium (LD) decay. The F_ST between accessions of cantalupensis and inodorus was 0.23. The F_ST between cantalupensis and various agrestis accessions was in a range of 0.19-0.53 and between inodorus and agrestis accessions was in a range of 0.21-0.59 indicating sporadic to wide ranging introgression. The EM (Expectation Maximization) algorithm was used for estimation of 1436 haplotypes. Average genome-wide LD decay for the melon genome was noted to be 9.27 Kb. In the current research, we focused on the genome-wide divergence underlying diverse melon horticultural groups. A high-resolution genetic map with 7153 loci was constructed. Genome-wide segregation distortion and recombination rate across various chromosomes were characterized. Melon has climacteric and non-climacteric morphotypes and wide variation for fruit firmness, a very important trait for transportation and shelf life. Various levels of QTLs were identified with high to moderate stringency and linked to fruit firmness using both genome-wide association study (GWAS) and biparental mapping. Gene annotation revealed some of the SNPs are located in β-D-xylosidase, glyoxysomal malate synthase, chloroplastic anthranilate phosphoribosyltransferase, and histidine kinase, the genes that were previously characterized for fruit ripening and softening in other crops.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1996-0) contains supplementary material, which is available to authorized users.

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.

RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum

BACKGROUND

The evolutionary basis of reproductive success in different environments is of major interest in the study of plant adaptation. Since the reproductive stage is particularly sensitive to drought, genes affecting reproductive success during this stage are key players in the evolution of adaptive mechanisms. We used an ecological genomics approach to investigate the reproductive response of drought-tolerant and sensitive wild barley accessions originating from different habitats in the Levant.

RESULTS

We sequenced mRNA extracted from spikelets at the flowering stage in drought-treated and control plants. The barley genome was used for a reference-guided assembly and differential expression analysis. Our approach enabled to detect biological processes affecting grain production under drought stress. We detected novel candidate genes and differentially expressed alleles associated with drought tolerance. Drought associated genes were shown to be more conserved than non-associated genes, and drought-tolerance genes were found to evolve more rapidly than other drought associated genes.

CONCLUSIONS

We show that reproductive success under drought stress is not a habitat-specific trait but a shared physiological adaptation that appeared to evolve recently in the evolutionary history of wild barley. Exploring the genomic basis of reproductive success under stress in crop wild progenitors is expected to have considerable ecological and economical applications.

Evaluation of marker-assisted selection for the stripe rust resistance gene Yr15, introgressed from wild emmer wheat

Stripe rust disease is caused by the fungus Puccinia striiformis f. sp. tritici and severely threatens wheat worldwide, repeatedly breaking resistance conferred by resistance genes and evolving more aggressive strains. Wild emmer wheat, Triticum dicoccoides, is an important source for novel stripe rust resistance (Yr) genes. Yr15, a major gene located on chromosome 1BS of T. dicoccoides, was previously reported to confer resistance to a broad spectrum of stripe rust isolates, at both seedling and adult plant stages. Introgressions of Yr15 into cultivated T. aestivum bread wheat and T. durum pasta wheat that began in the 1980s are widely used. In the present study, we aimed to validate SSR markers from the Yr15 region as efficient tools for marker-assisted selection (MAS) for introgression of Yr15 into wheat and to compare the outcome of gene introgression by MAS and by conventional phenotypic selection. Our findings establish the validity of MAS for introgression of Yr15 into wheat. We show that the size of the introgressed segment, defined by flanking markers, varies for both phenotypic selection and MAS. The genetic distance of the MAS marker from Yr15 and the number of backcross steps were the main factors affecting the length of the introgressed donor segments. Markers Xbarc8 and Xgwm493, which are the nearest flanking markers studied, were consistent and polymorphic in all 34 introgressions reported here and are therefore the most recommended markers for the introgression of Yr15 into wheat cultivars. Introgression directed by markers, rather than by phenotype, will facilitate simultaneous selection for multiple stripe rust resistant genes and will help to avoid escapees during the selection process.

Two-stage genome-wide search for epistasis with implementation to Recombinant Inbred Lines (RIL) populations

OBJECTIVE AND METHODS

This paper proposes an inegrative two-stage genome-wide search for pairwise epistasis on expression quantitative trait loci (eQTL). The traits are clustered into multi-trait complexes that account for correlations between them that may result from common epistasis effects. The search is done by first screening for epistatic regions and then using dense markers within the identified regions, resulting in substantial reduction in the number of tests for epistasis. The FDR is controlled using a hierarchical procedure that accounts for the search structure. Each combination of trait and marker-pair is tested using a model that accounts for both statistical and functional interpretations of epistasis and considers orthogonal effects, such that their contributions to heritability can be estimated individually. We examine the impact of using multi-trait complexes rather than single traits, and of using a hierarchical search for epistasis rather than skipping the initial screen for epistatic regions. We apply the proposed algorithm on Arabidopsis transcription data.

PRINCIPAL FINDINGS

Both epistasis detection power and heritability contributed by epistasis increased when using multi-trait complexes rather than single traits. Epistatic effects common to the eQTLs included in the complexes have higher chance of being identified by analysis of multi-trait complexes, particularly when epistatic effects on individual traits are small. Compared to direct testing for all potential epistatic effects, the hierarchical search was substantially more powerful in detecting epistasis, while controlling the FDR at the desired level. Association in functional roles within genomic regions was observed, supporting an initial screen for epistatic QTLs.

Threshing efficiency as an incentive for rapid domestication of emmer wheat

BACKGROUND AND AIMS

The harvesting method of wild and cultivated cereals has long been recognized as an important factor in the emergence of domesticated non-shattering ear genotypes. This study aimed to quantify the effects of spike brittleness and threshability on threshing time and efficiency in emmer wheat, and to evaluate the implications of post-harvest processes on domestication of cereals in the Near East.

METHODS

A diverse collection of tetraploid wheat genotypes, consisting of Triticum turgidum ssp. dicoccoides - the wild progenitor of domesticated wheat - traditional landraces, modern cultivars (T. turgidum ssp. durum) and 150 recombinant (wild × modern) inbred lines, was used in replicated controlled threshing experiments to quantify the effects of spike brittleness and threshability on threshing time and efficiency.

KEY RESULTS

The transition from a brittle hulled wild phenotype to non-brittle hulled phenotype (landraces) was associated with an approx. 30 % reduction in threshing time, whereas the transition from the latter to non-brittle free-threshing cultivars was associated with an approx. 85 % reduction in threshing time. Similar trends were obtained with groups of recombinant inbred lines showing extreme phenotypes of brittleness and threshability.

CONCLUSIONS

In tetraploid wheat, both non-brittle spike and free-threshing are labour-saving traits that increase the efficiency of post-harvest processing, which could have been an incentive for rapid domestication of the Near Eastern cereals, thus refuting the recently proposed hypothesis regarding extra labour associated with the domesticated phenotype (non-brittle spike) and its presumed role in extending the domestication episode time frame.

SNP discovery and chromosome anchoring provide the first physically-anchored hexaploid oat map and reveal synteny with model species

A physically anchored consensus map is foundational to modern genomics research; however, construction of such a map in oat (Avena sativa L., 2n = 6x = 42) has been hindered by the size and complexity of the genome, the scarcity of robust molecular markers, and the lack of aneuploid stocks. Resources developed in this study include a modified SNP discovery method for complex genomes, a diverse set of oat SNP markers, and a novel chromosome-deficient SNP anchoring strategy. These resources were applied to build the first complete, physically-anchored consensus map of hexaploid oat. Approximately 11,000 high-confidence in silico SNPs were discovered based on nine million inter-varietal sequence reads of genomic and cDNA origin. GoldenGate genotyping of 3,072 SNP assays yielded 1,311 robust markers, of which 985 were mapped in 390 recombinant-inbred lines from six bi-parental mapping populations ranging in size from 49 to 97 progeny. The consensus map included 985 SNPs and 68 previously-published markers, resolving 21 linkage groups with a total map distance of 1,838.8 cM. Consensus linkage groups were assigned to 21 chromosomes using SNP deletion analysis of chromosome-deficient monosomic hybrid stocks. Alignments with sequenced genomes of rice and Brachypodium provide evidence for extensive conservation of genomic regions, and renewed encouragement for orthology-based genomic discovery in this important hexaploid species. These results also provide a framework for high-resolution genetic analysis in oat, and a model for marker development and map construction in other species with complex genomes and limited resources.

New Diversity Arrays Technology (DArT) markers for tetraploid oat (Avena magna Murphy et Terrell) provide the first complete oat linkage map and markers linked to domestication genes from hexaploid A. sativa L

Nutritional benefits of cultivated oat (Avena sativa L., 2n = 6x = 42, AACCDD) are well recognized; however, seed protein levels are modest and resources for genetic improvement are scarce. The wild tetraploid, A. magna Murphy et Terrell (syn A. maroccana Gdgr., 2n = 4x = 28, CCDD), which contains approximately 31% seed protein, was hybridized with cultivated oat to produce a domesticated A. magna. Wild and cultivated accessions were crossed to generate a recombinant inbred line (RIL) population. Although these materials could be used to develop domesticated, high-protein oat, mapping and quantitative trait loci introgression is hindered by a near absence of genetic markers. Objectives of this study were to develop high-throughput, A. magna-specific markers; generate a genetic linkage map based on the A. magna RIL population; and map genes controlling oat domestication. A Diversity Arrays Technology (DArT) array derived from 10 A. magna genotypes was used to generate 2,688 genome-specific probes. These, with 12,672 additional oat clones, produced 2,349 polymorphic markers, including 498 (21.2%) from A. magna arrays and 1,851 (78.8%) from other Avena libraries. Linkage analysis included 974 DArT markers, 26 microsatellites, 13 SNPs, and 4 phenotypic markers, and resulted in a 14-linkage-group map. Marker-to-marker correlation coefficient analysis allowed classification of shared markers as unique or redundant, and putative linkage-group-to-genome anchoring. Results of this study provide for the first time a collection of high-throughput tetraploid oat markers and a comprehensive map of the genome, providing insights to the genome ancestry of oat and affording a resource for study of oat domestication, gene transfer, and comparative genomics.

Genetic analysis of wheat domestication and evolution under domestication

Wheat is undoubtedly one of the world's major food sources since the dawn of Near Eastern agriculture and up to the present day. Morphological, physiological, and genetic modifications involved in domestication and subsequent evolution under domestication were investigated in a tetraploid recombinant inbred line population, derived from a cross between durum wheat and its immediate progenitor wild emmer wheat. Experimental data were used to test previous assumptions regarding a protracted domestication process. The brittle rachis (Br) spike, thought to be a primary characteristic of domestication, was mapped to chromosome 2A as a single gene, suggesting, in light of previously reported Br loci (homoeologous group 3), a complex genetic model involved in spike brittleness. Twenty-seven quantitative trait loci (QTLs) conferring threshability and yield components (kernel size and number of kernels per spike) were mapped. The large number of QTLs detected in this and other studies suggests that following domestication, wheat evolutionary processes involved many genomic changes. The Br gene did not show either genetic (co-localization with QTLs) or phenotypic association with threshability or yield components, suggesting independence of the respective loci. It is argued here that changes in spike threshability and agronomic traits (e.g. yield and its components) are the outcome of plant evolution under domestication, rather than the result of a protracted domestication process. Revealing the genomic basis of wheat domestication and evolution under domestication, and clarifying their inter-relationships, will improve our understanding of wheat biology and contribute to further crop improvement.

Alteration in expression of hormone-related genes in wild emmer wheat roots associated with drought adaptation mechanisms

Transcriptomic and metabolomic profiles were used to unravel drought adaptation mechanisms in wild emmer wheat (Triticum turgidum ssp. dicoccoides), the progenitor of cultivated wheat, by comparing the response to drought stress in roots of genotypes contrasting in drought tolerance. The differences between the drought resistant (R) and drought susceptible (S) genotypes were characterized mainly by shifts in expression of hormone-related genes (e.g., gibberellins, abscisic acid (ABA) and auxin), including biosynthesis, signalling and response; RNA binding; calcium (calmodulin, caleosin and annexin) and phosphatidylinositol signalling, in the R genotype. ABA content in the roots of the R genotype was higher in the well-watered treatment and increased in response to drought, while in the S genotype ABA was invariant. The metabolomic profiling revealed in the R genotype a higher accumulation of tricarboxylic acid cycle intermediates and drought-related metabolites, including glucose, trehalose, proline and glycine. The integration of transcriptomics and metabolomics results indicated that adaptation to drought included efficient regulation and signalling pathways leading to effective bio-energetic processes, carbon metabolism and cell homeostasis. In conclusion, mechanisms of drought tolerance were identified in roots of wild emmer wheat, supporting our previous studies on the potential of this genepool as a valuable source for novel candidate genes to improve drought tolerance in cultivated wheat.

Multilevel regulation and signalling processes associated with adaptation to terminal drought in wild emmer wheat

Low water availability is the major environmental factor limiting crop productivity. Transcriptome analysis was used to study terminal drought response in wild emmer wheat, Triticum dicoccoides, genotypes contrasting in their productivity and yield stability under drought stress. A total of 5,892 differentially regulated transcripts were identified between drought and well-watered control and/or between drought resistant (R) and drought susceptible (S) genotypes. Functional enrichment analyses revealed that multilevel regulatory and signalling processes were significantly enriched among the drought-induced transcripts, in particular in the R genotype. Therefore, further analyses were focused on selected 221 uniquely expressed or highly abundant transcripts in the R genotype, as potential candidates for drought resistance genes. Annotation of the 221 genes revealed that 26% of them are involved in multilevel regulation, including: transcriptional regulation, RNA binding, kinase activity and calcium and abscisic acid signalling implicated in stomatal closure. Differential expression patterns were also identified in genes known to be involved in drought adaptation pathways, such as: cell wall adjustment, cuticular wax deposition, lignification, osmoregulation, redox homeostasis, dehydration protection and drought-induced senescence. These results demonstrate the potential of wild emmer wheat as a source for candidate genes for improving drought resistance.

Quantitative trait loci conferring grain mineral nutrient concentrations in durum wheat x wild emmer wheat RIL population

Mineral nutrient malnutrition, and particularly deficiency in zinc and iron, afflicts over 3 billion people worldwide. Wild emmer wheat, Triticum turgidum ssp. dicoccoides, genepool harbors a rich allelic repertoire for mineral nutrients in the grain. The genetic and physiological basis of grain protein, micronutrients (zinc, iron, copper and manganese) and macronutrients (calcium, magnesium, potassium, phosphorus and sulfur) concentration was studied in tetraploid wheat population of 152 recombinant inbred lines (RILs), derived from a cross between durum wheat (cv. Langdon) and wild emmer (accession G18-16). Wide genetic variation was found among the RILs for all grain minerals, with considerable transgressive effect. A total of 82 QTLs were mapped for 10 minerals with LOD score range of 3.2-16.7. Most QTLs were in favor of the wild allele (50 QTLs). Fourteen pairs of QTLs for the same trait were mapped to seemingly homoeologous positions, reflecting synteny between the A and B genomes. Significant positive correlation was found between grain protein concentration (GPC), Zn, Fe and Cu, which was supported by significant overlap between the respective QTLs, suggesting common physiological and/or genetic factors controlling the concentrations of these mineral nutrients. Few genomic regions (chromosomes 2A, 5A, 6B and 7A) were found to harbor clusters of QTLs for GPC and other nutrients. These identified QTLs may facilitate the use of wild alleles for improving grain nutritional quality of elite wheat cultivars, especially in terms of protein, Zn and Fe.

Genomic dissection of drought resistance in durum wheat x wild emmer wheat recombinant inbreed line population

Drought is the major factor limiting wheat productivity worldwide. The gene pool of wild emmer wheat, Triticum turgidum ssp. dicoccoides, harbours a rich allelic repertoire for morpho-physiological traits conferring drought resistance. The genetic and physiological bases of drought responses were studied here in a tetraploid wheat population of 152 recombinant inbreed lines (RILs), derived from a cross between durum wheat (cv. Langdon) and wild emmer (acc# G18-16), under contrasting water availabilities. Wide genetic variation was found among RILs for all studied traits. A total of 110 quantitative trait loci (QTLs) were mapped for 11 traits, with LOD score range of 3.0-35.4. Several QTLs showed environmental specificity, accounting for productivity and related traits under water-limited (20 QTLs) or well-watered conditions (15 QTLs), and in terms of drought susceptibility index (22 QTLs). Major genomic regions controlling productivity and related traits were identified on chromosomes 2B, 4A, 5A and 7B. QTLs for productivity were associated with QTLs for drought-adaptive traits, suggesting the involvement of several strategies in wheat adaptation to drought stress. Fifteen pairs of QTLs for the same trait were mapped to seemingly homoeologous positions, reflecting synteny between the A and B genomes. The identified QTLs may facilitate the use of wild alleles for improvement of drought resistance in elite wheat cultivars.

High-density genetic map of durum wheat x wild emmer wheat based on SSR and DArT markers

A genetic linkage map of tetraploid wheat was constructed based on a cross between durum wheat [Triticum turgidum ssp. durum (Desf.) MacKey] cultivar Langdon and wild emmer wheat [T. turgidum ssp. dicoccoides (Körn.) Thell.] accession G18-16. One hundred and fifty-two single-seed descent derived F(6) recombinant inbred lines (RILs) were analyzed with a total of 690 loci, including 197 microsatellite and 493 DArT markers. Linkage analysis defined 14 linkage groups. Most markers were mapped to the B-genome (60%), with an average of 57 markers per chromosome and the remaining 40% mapped to the A-genome, with an average of 39 markers per chromosome. To construct a stabilized (skeleton) map, markers interfering with map stability were removed. The skeleton map consisted of 307 markers with a total length of 2,317 cM and average distance of 7.5 cM between adjacent markers. The length of individual chromosomes ranged between 112 cM for chromosome 4B to 217 cM for chromosome 3B. A fraction (30.1%) of the markers deviated significantly from the expected Mendelian ratios; clusters of loci showing distorted segregation were found on chromosomes 1A, 1BL, 2BS, 3B, and 4B. DArT markers showed high proportion of clustering, which may be indicative of gene-rich regions. Three hundred and fifty-two new DArT markers were mapped for the first time on the current map. This map provides a useful groundwork for further genetic analyses of important quantitative traits, positional cloning, and marker-assisted selection, as well as for genome comparative genomics and genome organization studies in wheat and other cereals.

Genetic dissection of intermated recombinant inbred lines using a new genetic map of maize

A new genetic map of maize, ISU-IBM Map4, that integrates 2029 existing markers with 1329 new indel polymorphism (IDP) markers has been developed using intermated recombinant inbred lines (IRILs) from the intermated B73xMo17 (IBM) population. The website http://magi.plantgenomics.iastate.edu provides access to IDP primer sequences, sequences from which IDP primers were designed, optimized marker-specific PCR conditions, and polymorphism data for all IDP markers. This new gene-based genetic map will facilitate a wide variety of genetic and genomic research projects, including map-based genome sequencing and gene cloning. The mosaic structures of the genomes of 91 IRILs, an important resource for identifying and mapping QTL and eQTL, were defined. Analyses of segregation data associated with markers genotyped in three B73/Mo17-derived mapping populations (F2, Syn5, and IBM) demonstrate that allele frequencies were significantly altered during the development of the IBM IRILs. The observations that two segregation distortion regions overlap with maize flowering-time QTL suggest that the altered allele frequencies were a consequence of inadvertent selection. Detection of two-locus gamete disequilibrium provides another means to extract functional genomic data from well-characterized plant RILs.

The dynamics of the rec-system in variable environments: haploid selection in a cyclical two-state environment

The dynamics of a 3-locus infinite population with non-overlapping generations and panmixia was studied. Loci are di-allelic: two loci affect fitness under cyclical symmetric haploid selection while the third one is a modifier of recombination (rec-modifier). Selection favors alternatively haplotypes AB and ab or Ab and aB. It has been proven that under alternating selection (when period of selection consists of two generations) a dominant suppressor of recombination is displaced and the allele for non-zero recombination becomes fixed within the population. For populations with inversion heterozygosity within the selective system (i.e. with zero recombination in heterozygote for rec-modifier and non-zero for homozygotes) fixation of one of the alleles (depending on the initial point) at the rec-modifier locus is predicted. For other values of recombination parameters, the behavior of the system was studied numerically. A full bifurcation picture of parameters was obtained. Many of the results related to the case of a two-generation period hold also in the case of longer period lengths.

Multilocus consensus genetic maps (MCGM): formulation, algorithms, and results

In process of creating genetic maps different labs/research groups obtain overlapping parts of the map. Merging these parts into one integrative map is based on looking for maximum shared marker orders among the maps. Really, not all shared markers of such maps have consensus order that obstructs building of the integrative maps. In this paper we propose a new approach to build verified multilocus consensus genetic maps in which shared markers always are integrated in stable consensus order. The approach is based on combined analysis of initial mapping data rather than manipulating with previously constructed maps. We show that more effective and reliable solutions may be obtained based on "synchronized ordering" facilitated by cycles of "re-sampling-->ordering-->removing unstable markers". The proposed formulation of consensus genetic mapping can be considered as a version of traveling salesperson problem (TSP) that we refer to as synchronized-TSP. From the viewpoint of optimization, synchronized-TSP belongs to discrete constrained optimization problems. Earlier we developed new powerful and fast guided evolution strategy algorithms for some types of discrete constrained optimization. These algorithms were used here as a basis for solving more challenging problems of consensual marker ordering.

Evolutionary regulation of the blind subterranean mole rat, Spalax, revealed by genome-wide gene expression

We applied genome-wide gene expression analysis to the evolutionary processes of adaptive speciation of the Israeli blind subterranean mole rats of the Spalax ehrenbergi superspecies. The four Israeli allospecies climatically and adaptively radiated into the cooler, mesic northern domain (N) and warmer, xeric southern domain (S). The kidney and brain mRNAs of two N and two S animals were examined through cross-species hybridizations with two types of Affymetrix arrays (mouse and rat) and muscle mRNA of six N and six S animals with spotted cDNA mouse arrays. The initial microarray analysis was hypothesis-free, i.e., conducted without reference to the origin of animals. Principal component analysis revealed that 20-30% of the expression signal variability could be explained by the differentiation of N-S species. Similar N-S effects were obtained for all tissues and types of arrays: two Affymetrix microarrays using probe oligomer signals and the spotted array. Likewise, ANOVA and t test statistics demonstrated significant N-S ecogeographic divergence and region-tissue specificity in gene expression. Analysis of differential gene expression between species corroborates previous results deduced by allozymes and DNA molecular polymorphisms. Functional categories show significant N-S ecologic putative adaptive divergent up-regulation of genes highlighting a higher metabolism in N, and potential adaptive brain activity and kidney urine cycle pathways in S. The present results confirm ecologic-genomic separation of blind mole rats into N and S. Gene expression regulation appears to be central to the evolution of blind mole rats.

Fast and high precision algorithms for optimization in large-scale genomic problems

There are several very difficult problems related to genetic or genomic analysis that belong to the field of discrete optimization in a set of all possible orders. With n elements (points, markers, clones, sequences, etc.), the number of all possible orders is n!/2 and only one of these is considered to be the true order. A classical formulation of a similar mathematical problem is the well-known traveling salesperson problem model (TSP). Genetic analogues of this problem include: ordering in multilocus genetic mapping, evolutionary tree reconstruction, building physical maps (contig assembling for overlapping clones and radiation hybrid mapping), and others. A novel, fast and reliable hybrid algorithm based on evolution strategy and guided local search discrete optimization was developed for TSP formulation of the multilocus mapping problems. High performance and high precision of the employed algorithm named guided evolution strategy (GES) allows verification of the obtained multilocus orders based on different computing-intensive approaches (e.g., bootstrap or jackknife) for detection and removing unreliable marker loci, hence, stabilizing the resulting paths. The efficiency of the proposed algorithm is demonstrated on standard TSP problems and on simulated data of multilocus genetic maps up to 1000 points per linkage group.

Genetic basis of barley caryopsis dormancy and seedling desiccation tolerance at the germination stage

The genomic regions controlling caryopsis dormancy and seedling desiccation tolerance were identified using 152 F4 lines derived from a cross between Mona, a Swedish cultivar, and an Israeli xeric wild barley Hordeum spontaneum genotype collected at Wadi Qilt, Israel. Dormancy, the inability of a viable seed to germinate, and desiccation tolerance, the ability of the desiccated seedlings to revive after rehydration, were characterized by fitting the germination and revival data with growth curves, using three parameters: minimum, maximum, and slope of germination or revival rate derived by the least square method. The genetic map was constructed with 85 genetic markers (SSRs, AFLPs, STSs, and Dhn genes) using the MULTIPOINT: -mapping algorithm. Quantitative trait loci (QTLs) mapping was conducted with the MULTIQTL: package. Ten genomic regions were detected that affected the target traits, seven of which affected both dormancy and desiccation tolerance traits. Both the wild barley genotype and the Swedish cultivar contributed the favorite alleles for caryopsis dormancy, whereas seedling desiccation tolerance was attributed to alleles descending from the cultivar. The results indicate that some barley dormancy genes are lost during domestication and that dormancy QTLs are associated with abiotic stress tolerance.

Genetic basis of barley caryopsis dormancy and seedling desiccation tolerance at the germination stage

The genomic regions controlling caryopsis dormancy and seedling desiccation tolerance were identified using 152 F4 lines derived from a cross between Mona, a Swedish cultivar, and an Israeli xeric wild barley Hordeum spontaneum genotype collected at Wadi Qilt, Israel. Dormancy, the inability of a viable seed to germinate, and desiccation tolerance, the ability of the desiccated seedlings to revive after rehydration, were characterized by fitting the germination and revival data with growth curves, using three parameters: minimum, maximum, and slope of germination or revival rate derived by the least square method. The genetic map was constructed with 85 genetic markers (SSRs, AFLPs, STSs, and Dhn genes) using the MULTIPOINT: -mapping algorithm. Quantitative trait loci (QTLs) mapping was conducted with the MULTIQTL: package. Ten genomic regions were detected that affected the target traits, seven of which affected both dormancy and desiccation tolerance traits. Both the wild barley genotype and the Swedish cultivar contributed the favorite alleles for caryopsis dormancy, whereas seedling desiccation tolerance was attributed to alleles descending from the cultivar. The results indicate that some barley dormancy genes are lost during domestication and that dormancy QTLs are associated with abiotic stress tolerance.

Wild barley eibi1 mutation identifies a gene essential for leaf water conservation

Drought is a major abiotic stress that limits plant growth and crop productivity. A spontaneous wilty mutant (eibi1) hypersensitive to drought was identified from wild barley (Hordeum spontaneum Koch). eibi1 showed the highest relative water loss rate among the known wilty mutants, which indicates that eibi1 is one of the most drought-sensitive mutants. eibi1 had the same abscisic acid (ABA) level, the same ability to accumulate stress-induced ABA, and the same stomatal movement in response to light, dark, drought, and exogenous ABA as the wild type, revealing that eibi1 was neither an ABA-deficient nor an ABA-insensitive mutant. The eibi1 leaves had a larger chlorophyll efflux rate in 80% ethanol than the wild-type leaves; and the transpiration rate of eibi1 was more closely related to chlorophyll efflux rate than to stomatal density, demonstrating that the cuticle of eibi1 was defective. eibi1 will be a promising candidate to study the actual barrier layer in the cuticle that limits water loss of the plant. Exogenous ABA reduced leaf length growth in eibi1 more than in the wild type, implying an interaction on plant growth of ABA signal transduction and the eibi1 product. One may infer that the eibi1 product may reverse the growth inhibition induced by ABA.

Microsatellites within genes: structure, function, and evolution

Recently, increasingly more microsatellites, or simple sequence repeats (SSRs) have been found and characterized within protein-coding genes and their untranslated regions (UTRs). These data provide useful information to study possible SSR functions. Here, we review SSR distributions within expressed sequence tags (ESTs) and genes including protein-coding, 3'-UTRs and 5'-UTRs, and introns; and discuss the consequences of SSR repeat-number changes in those regions of both prokaryotes and eukaryotes. Strong evidence shows that SSRs are nonrandomly distributed across protein-coding regions, UTRs, and introns. Substantial data indicates that SSR expansions and/or contractions in protein-coding regions can lead to a gain or loss of gene function via frameshift mutation or expanded toxic mRNA. SSR variations in 5'-UTRs could regulate gene expression by affecting transcription and translation. The SSR expansions in the 3'-UTRs cause transcription slippage and produce expanded mRNA, which can be accumulated as nuclear foci, and which can disrupt splicing and, possibly, disrupt other cellular function. Intronic SSRs can affect gene transcription, mRNA splicing, or export to cytoplasm. Triplet SSRs located in the UTRs or intron can also induce heterochromatin-mediated-like gene silencing. All these effects caused by SSR expansions or contractions within genes can eventually lead to phenotypic changes. SSRs within genes evolve through mutational processes similar to those for SSRs located in other genomic regions including replication slippage, point mutation, and recombination. These mutational processes generate DNA changes that should be connected by DNA mismatch repair (MMR) system. Mutation that has escaped from the MMR system correction would become new alleles at the SSR loci, and then regulate and/or change gene products, and eventually lead to phenotype changes. Therefore, SSRs within genes should be subjected to stronger selective pressure than other genomic regions because of their functional importance. These SSRs may provide a molecular basis for fast adaptation to environmental changes in both prokaryotes and eukaryotes.

Efficient multipoint mapping: making use of dominant repulsion-phase markers

The paper is devoted to the problem of multipoint gene ordering with a particular focus on "dominance" complication that acts differently in conditions of coupling-phase and repulsion-phase markers. To solve the problem we split the dataset into two complementary subsets each containing shared codominant markers and dominant markers in the coupling-phase only. Multilocus ordering in the proposed algorithm is based on pairwise recombination frequencies and using the well-known travelling salesman problem (TSP) formalization. To obtain accurate results, we developed a multiphase algorithm that includes synchronized-marker ordering of two subsets assisted by re-sampling-based map verification, combining the resulting maps into an integrated map followed by verification of the integrated map. A new synchronized Evolution-Strategy discrete optimization algorithm was developed here for the proposed multilocus ordering approach in which common codominant markers facilitate stabilization of the marker order of the two complementary maps. High performance of the employed algorithm allows systematic treatment for the problem of verification of the obtained multilocus orders, based on computing-intensive bootstrap and jackknife technologies for detection and removing unreliable marker scores. The efficiency of the proposed algorithm was demonstrated on simulated and real data.

Biodiversity of 52 chicken populations assessed by microsatellite typing of DNA pools

In a project on the biodiversity of chickens funded by the European Commission (EC), eight laboratories collaborated to assess the genetic variation within and between 52 populations from a wide range of chicken types. Twenty-two di-nucleotide microsatellite markers were used to genotype DNA pools of 50 birds from each population. The polymorphism measures for the average, the least polymorphic population (inbred C line) and the most polymorphic population (Gallus gallus spadiceus) were, respectively, as follows: number of alleles per locus, per population: 3.5, 1.3 and 5.2; average gene diversity across markers: 0.47, 0.05 and 0.64; and proportion of polymorphic markers: 0.91, 0.25 and 1.0. These were in good agreement with the breeding history of the populations. For instance, unselected populations were found to be more polymorphic than selected breeds such as layers. Thus DNA pools are effective in the preliminary assessment of genetic variation of populations and markers. Mean genetic distance indicates the extent to which a given population shares its genetic diversity with that of the whole tested gene pool and is a useful criterion for conservation of diversity. The distribution of population-specific (private) alleles and the amount of genetic variation shared among populations supports the hypothesis that the red jungle fowl is the main progenitor of the domesticated chicken.

Genetic effects on microsatellite diversity in wild emmer wheat (Triticum dicoccoides) at the Yehudiyya microsite, Israel

This study investigated allele size constraints and clustering, and genetic effects on microsatellite (simple sequence repeat, SSR) diversity at 28 loci comprising seven types of tandem repeated dinucleotide motifs in a natural population of wild emmer wheat, Triticum dicoccoides, from a shade vs sun microsite in Yehudiyya, northeast of the Sea of Galilee, Israel. It was found that allele distribution at SSR loci is clustered and constrained with lower or higher boundary. This may imply that SSR have functional significance and natural constraints. Genetic factors, involving genome, chromosome, motif, and locus significantly affected SSR diversity. Genome B appeared to have a larger average repeat number (ARN), but lower variance in repeat number (sigma(ARN)(2)), and smaller number of alleles per locus than genome A. SSRs with compound motifs showed larger ARN than those with perfect motifs. The effects of replication slippage and recombinational effects (eg, unequal crossing over) on SSR diversity varied with SSR motifs. Ecological stresses (sun vs shade) may affect mutational mechanisms, influencing the level of SSR diversity by both processes.

Ecological stress and sex evolution in soil microfungi

The elucidation of the origin and maintenance of sex is a major unsolved problem in evolutionary biology. A number of hypotheses have been elaborated, but the scarcity of empirical data limits further progress. During recent years, the general inclination has changed towards pluralistic models of sex evolution, due partly to an increased diversity of studied organisms. Fungi are among the most promising organisms for testing sexual causation, as demonstrated in recent laboratory experiments. However, reconciling theory and evidence necessitates critical field observations. Here, we report new estimates of the distribution of morphologically sexual and asexual soil microfungi in nature, which indicate a remarkable trend towards increased sexuality with increasing climatic stress.

Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review

Microsatellites, or tandem simple sequence repeats (SSR), are abundant across genomes and show high levels of polymorphism. SSR genetic and evolutionary mechanisms remain controversial. Here we attempt to summarize the available data related to SSR distribution in coding and noncoding regions of genomes and SSR functional importance. Numerous lines of evidence demonstrate that SSR genomic distribution is nonrandom. Random expansions or contractions appear to be selected against for at least part of SSR loci, presumably because of their effect on chromatin organization, regulation of gene activity, recombination, DNA replication, cell cycle, mismatch repair system, etc. This review also discusses the role of two putative mutational mechanisms, replication slippage and recombination, and their interaction in SSR variation.

Climatic effects on microsatellite diversity in wild emmer wheat (Triticum dicoccoides) at the Yehudiyya microsite, Israel

Microsatellite (SSR) diversity at 28 loci comprising seven types of tandem dinucleotide repeated motifs was analyzed in 105 individual plants of wild emmer wheat, Triticum dicoccoides, from a microsite in Yehudiyya, northeast of the Sea of Galilee, Israel. The study area was less than 1000 m(2) and involved 12 paired plots distributed in a mosaic pattern. Each experiment involved very close (a few meters apart), but sharply divergent, microclimatic niches in the open park forest of Tabor oak: (1) sun, between trees, and (2) shade, under tree canopy. Significant microclimatic divergence characterized many loci displaying asymmetric and non-random distribution of repeat numbers. Niche-specific and niche-unique alleles and linkage disequilibria were found in the two sub-populations. Microsatellite diversity at both single- and two-locus levels is affected by microclimatic environment. The evidence reflects effects of ecological stresses and natural selection on SSR diversity, resulting presumably in adaptive structures.