Grain dispersal mechanism in cereals arose from a genome duplication followed by changes in spatial expression of genes involved in pollen development.

KEY MESSAGE: Grain disarticulation in wild progenitor of wheat and barley evolved through a local duplication event followed by neo-functionalization resulting from changes in location of gene expression. One of the most critical events in the process of cereal domestication was the loss of the natural mode of grain dispersal. Grain dispersal in barley is controlled by two major genes, Btr1 and Btr2, which affect the thickness of cell walls around the disarticulation zone. The barley genome also encodes Btr1-like and Btr2-like genes, which have been shown to be the ancestral copies. While Btr and Btr-like genes are non-redundant, the biological function of Btr-like genes is unknown. We explored the potential biological role of the Btr-like genes by surveying their expression profile across 212 publicly available transcriptome datasets representing diverse organs, developmental stages and stress conditions. We found that Btr1-like and Btr2-like are expressed exclusively in immature anther samples throughout Prophase I of meiosis within the meiocyte. The similar and restricted expression profile of these two genes suggests they are involved in a common biological function. Further analysis revealed 141 genes co-expressed with Btr1-like and 122 genes co-expressed with Btr2-like, with 105 genes in common, supporting Btr-like genes involvement in a shared molecular pathway. We hypothesize that the Btr-like genes play a crucial role in pollen development by facilitating the formation of the callose wall around the meiocyte or in the secretion of callase by the tapetum. Our data suggest that Btr genes retained an ancestral function in cell wall modification and gained a new role in grain dispersal due to changes in their spatial expression becoming spike specific after gene duplication.

The effect of heat stress on sugar beet recombination.

Meiotic recombination plays a crucial role in plant breeding through the creation of new allelic combinations. Therefore, lack of recombination in some genomic regions constitutes a constraint for breeding programmes. In sugar beet, one of the major crops in Europe, recombination occurs mainly in the distal portions of the chromosomes, and so the development of simple approaches to change this pattern is of considerable interest for future breeding and genetics. In the present study, the effect of heat stress on recombination in sugar beet was studied by treating F1 plants at 28 °C/25 °C (day/night) and genotyping the progeny. F1 plants were reciprocally backcrossed allowing the study of male and female meiosis separately. Genotypic data indicated an overall increase in crossover frequency of approximately one extra crossover per meiosis, with an associated increase in pericentromeric recombination under heat treatment. Our data indicate that the changes were mainly induced by alterations in female meiosis only, showing that heterochiasmy in sugar beet is reduced under heat stress. Overall, despite the associated decrease in fertility, these data support the potential use of heat stress to foster recombination in sugar beet breeding programmes.

3D RNA-seq: a powerful and flexible tool for rapid and accurate differential expression and alternative splicing analysis of RNA-seq data for biologists.

RNA-sequencing (RNA-seq) analysis of gene expression and alternative splicing should be routine and robust but is often a bottleneck for biologists because of different and complex analysis programs and reliance on specialized bioinformatics skills. We have developed the '3D RNA-seq' App, an R shiny App and web-based pipeline for the comprehensive analysis of RNA-seq data from any organism. It represents an easy-to-use, flexible and powerful tool for analysis of both gene and transcript-level gene expression to identify differential gene/transcript expression, differential alternative splicing and differential transcript usage (3D) as well as isoform switching from RNA-seq data. 3D RNA-seq integrates state-of-the-art differential expression analysis tools and adopts best practice for RNA-seq analysis. The program is designed to be run by biologists with minimal bioinformatics experience (or by bioinformaticians) allowing lab scientists to analyse their RNA-seq data. It achieves this by operating through a user-friendly graphical interface which automates the data flow through the programs in the pipeline. The comprehensive analysis performed by 3D RNA-seq is extremely rapid and accurate, can handle complex experimental designs, allows user setting of statistical parameters, visualizes the results through graphics and tables, and generates publication quality figures such as heat-maps, expression profiles and GO enrichment plots. The utility of 3D RNA-seq is illustrated by analysis of data from a time-series of cold-treated Arabidopsis plants and from dexamethasone-treated male and female mouse cortex and hypothalamus data identifying dexamethasone-induced sex- and brain region-specific differential gene expression and alternative splicing.

Preparation of Barley Pollen Mother Cells for Confocal and Super Resolution Microscopy.

Recombination (crossover) drives the release of genetic diversity in plant breeding programs. However, in barley, recombination is skewed toward the telomeric ends of its seven chromosomes, restricting the re-assortment of about 30% of the genes located in the centromeric regions of its large 5.1 Gb genome. A better understanding of meiosis and recombination could provide ways of modulating crossover distribution and frequency in barley as well as in other grasses, including wheat. While most research on recombination has been carried out in the model plant Arabidopsis thaliana, recent studies in barley (Hordeum Vulgare) have provided new insights into the control of crossing over in large genome species. A major achievement in these studies has been the use of cytological procedures to follow meiotic events. This protocol provides detailed practical steps required to perform immunostaining of barley meiocytes (pollen mother cells) for confocal or structured illumination microscopy.

Evolutionary relationships among barley and Arabidopsis core circadian clock and clock-associated genes.

The circadian clock regulates a multitude of plant developmental and metabolic processes. In crop species, it contributes significantly to plant performance and productivity and to the adaptation and geographical range over which crops can be grown. To understand the clock in barley and how it relates to the components in the Arabidopsis thaliana clock, we have performed a systematic analysis of core circadian clock and clock-associated genes in barley, Arabidopsis and another eight species including tomato, potato, a range of monocotyledonous species and the moss, Physcomitrella patens. We have identified orthologues and paralogues of Arabidopsis genes which are conserved in all species, monocot/dicot differences, species-specific differences and variation in gene copy number (e.g. gene duplications among the various species). We propose that the common ancestor of barley and Arabidopsis had two-thirds of the key clock components identified in Arabidopsis prior to the separation of the monocot/dicot groups. After this separation, multiple independent gene duplication events took place in both monocot and dicot ancestors.

Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps.

The genome of potato, a major global food crop, was recently sequenced. The work presented here details the integration of the potato reference genome (DM) with a new sequence-tagged site marker-based linkage map and other physical and genetic maps of potato and the closely related species tomato. Primary anchoring of the DM genome assembly was accomplished by the use of a diploid segregating population, which was genotyped with several types of molecular genetic markers to construct a new ~936 cM linkage map comprising 2469 marker loci. In silico anchoring approaches used genetic and physical maps from the diploid potato genotype RH89-039-16 (RH) and tomato. This combined approach has allowed 951 superscaffolds to be ordered into pseudomolecules corresponding to the 12 potato chromosomes. These pseudomolecules represent 674 Mb (~93%) of the 723 Mb genome assembly and 37,482 (~96%) of the 39,031 predicted genes. The superscaffold order and orientation within the pseudomolecules are closely collinear with independently constructed high density linkage maps. Comparisons between marker distribution and physical location reveal regions of greater and lesser recombination, as well as regions exhibiting significant segregation distortion. The work presented here has led to a greatly improved ordering of the potato reference genome superscaffolds into chromosomal "pseudomolecules".

QTL mapping of yield, agronomic and quality traits in tetraploid potato (Solanum tuberosum subsp. tuberosum).

Interval mapping of quantitative trait loci (QTL) for 16 yield, agronomic and quality traits in potato was performed on a tetraploid full-sib family comprising 227 clones from a cross between processing clone 12601ab1 and table cultivar Stirling. Thirty-eight AFLP primer combinations and six SSRs provided 514 informative markers which formed a molecular marker map comprising 12 linkage groups (LGs) in 12601ab1 (nine with four homologous chromosomes) which were aligned with 12 in Stirling (11 with four homologous chromosomes), with four partial groups remaining in 12601ab1. Two LGs were identified unequivocally as chromosomes IV and V and eight others were tentatively assigned with chromosomes VII and X unidentified. All of the traits scored had moderately high heritabilities with 54-92% of the variation in clone means over 3 years and two replicates being due to genetic differences. A total of 39 QTLs were identified. A QTL for maturity was identified on chromosome V which explained 56% of the phenotypic variance, whereas the other QTLs individually explained between 5.4 and 16.5%. However, six QTLs were detected for after-cooking blackening and four for each of regularity of tuber shape, fry colour both after storage at 4 and 10 degrees C and sprouting. Just two QTLs were found for each of yield, the two 'overall' scores, crop emergence, tuber size and common scab and just one QTL was detected for each of dry matter content, keeping quality, growth cracks and internal condition. The implications for practical potato breeding and for practical linkage and QTL analysis in autotetraploids are discussed.

Construction of a 10,000-marker ultradense genetic recombination map of potato: providing a framework for accelerated gene isolation and a genomewide physical map.

An ultradense genetic linkage map with >10,000 AFLP loci was constructed from a heterozygous diploid potato population. To our knowledge, this is the densest meiotic recombination map ever constructed. A fast marker-ordering algorithm was used, based on the minimization of the total number of recombination events within a given marker order in combination with genotyping error-detection software. This resulted in "skeleton bin maps," which can be viewed as the most parsimonious marker order. The unit of distance is not expressed in centimorgans but in "bins." A bin is a position on the genetic map with a unique segregation pattern that is separated from adjacent bins by a single recombination event. Putative centromeres were identified by a strong clustering of markers, probably due to cold spots for recombination. Conversely, recombination hot spots resulted in large intervals of up to 15 cM without markers. The current level of marker saturation suggests that marker density is proportional to physical distance and independent of recombination frequency. Most chromatids (92%) recombined once or never, suggesting strong chiasma interference. Absolute chiasma interference within a chromosome arm could not be demonstrated. Two examples of contig construction and map-based cloning have demonstrated that the marker spacing was in accordance with the expected physical distance: approximately one marker per BAC length. Currently, the markers are used for genetic anchoring of a physical map of potato to deliver a sequence-ready minimal tiling path of BAC contigs of specific chromosomal regions for the potato genome sequencing consortium (http://www.potatogenome.net).

A single domestication for potato based on multilocus amplified fragment length polymorphism genotyping.

The cultivated potato, Solanum tuberosum, ultimately traces its origin to Andean and Chilean landraces developed by pre-Colombian cultivators. These Andean landraces exhibit tremendous morphological and genetic diversity, and are distributed throughout the Andes, from western Venezuela to northern Argentina, and in southern Chile. The wild species progenitors of these landraces have long been in dispute, but all hypotheses center on a group of approximately 20 morphologically very similar tuber-bearing (Solanum section Petota) wild taxa referred to as the S. brevicaule complex, distributed from central Peru to northern Argentina. We present phylogenetic analyses based on the representative cladistic diversity of 362 individual wild (261) and landrace (98) members of potato (all tuber-bearing) and three outgroup non-tuber-bearing members of Solanum section Etuberosum, genotyped with 438 robust amplified fragment length polymorphisms. Our analyses are consistent with a hypothesis of a "northern" (Peru) and "southern" (Bolivia and Argentina) cladistic split for members of the S. brevicaule complex, and with the need for considerable reduction of species in the complex. In contrast to all prior hypotheses, our data support a monophyletic origin of the landrace cultivars from the northern component of this complex in Peru, rather than from multiple independent origins from various northern and southern members.

Interval mapping of quantitative trait loci for resistance to late blight [Phytophthora infestans (Mont.) de Bary], height and maturity in a tetraploid population of potato (Solanum tuberosum subsp. tuberosum).

Interval mapping of quantitative trait loci (QTL) for resistance to late blight, height, and maturity was performed on a tetraploid full-sib family of potato comprising 227 clones from a cross between a susceptible parent, 12601ab1, and a resistant cultivar, Stirling, which were of similar height and main crop maturity. Thirty-eight AFLP primer combinations provided 585 informative markers, and 23 SSRs proved useful for identifying linkage groups (LGs). A simplex QTL allele was found on LGV of Stirling close to marker STM3179, which was associated with early maturity, short plants, and susceptibility to blight and explained 54.7, 26.5, 26.3, and 17.5% of the variation for maturity, height, tuber blight, and foliage blight. When the residuals from the regressions of foliage and tuber blight on maturity were analyzed, there was no significant effect of a QTL on LGV, but a duplex QTL allele for resistance was found on LGIV of Stirling, which explained 30.7 and 13.6% of the variation for foliage and tuber blight on an additive model. Partial dominance for resistance explained even more of the variation, up to 37.2% for foliage blight. A major gene for blight resistance in Stirling was also mapped to LGXI.

Toward a marker-dense meiotic map of the potato genome: lessons from linkage group I.

Segregation data were obtained for 1260 potato linkage group I-specific AFLP loci from a heterozygous diploid potato population. Analytical tools that identified potential typing errors and/or inconsistencies in the data and that assembled cosegregating markers into bins were applied. Bins contain multiple-marker data sets with an identical segregation pattern, which is defined as the bin signature. The bin signatures were used to construct a skeleton bin map that was based solely on observed recombination events. Markers that did not match any of the bin signatures exactly (and that were excluded from the calculation of the skeleton bin map) were placed on the map by maximum likelihood. The resulting maternal and paternal maps consisted of 95 and 101 bins, respectively. Markers derived from EcoRI/MseI, PstI/MseI, and SacI/MseI primer combinations showed different genetic distributions. Approximately three-fourths of the markers placed into a bin were considered to fit well on the basis of an estimated residual "error rate" of 0-3%. However, twice as many PstI-based markers fit badly, suggesting that parental PstI-site methylation patterns had changed in the population. Recombination frequencies were highly variable across the map. Inert, presumably centromeric, regions caused extensive marker clustering while recombination hotspots (or regions identical by descent) resulted in empty bins, despite the level of marker saturation.