The transcriptome landscape of developing barley seeds.

Cereal grains are an important source of food and feed. To provide comprehensive spatiotemporal information about biological processes in developing seeds of cultivated barley (Hordeum vulgare L. subsp. vulgare), we performed a transcriptomic study of the embryo, endosperm, and seed maternal tissues collected from grains 4-32 days after pollination. Weighted gene co-expression network and motif enrichment analyses identified specific groups of genes and transcription factors (TFs) potentially regulating barley seed tissue development. We defined a set of tissue-specific marker genes and families of TFs for functional studies of the pathways controlling barley grain development. Assessing selected groups of chromatin regulators revealed that epigenetic processes are highly dynamic and likely play a major role during barley endosperm development. The repressive H3K27me3 modification is globally reduced in endosperm tissues and at specific genes related to development and storage compounds. Altogether, this atlas uncovers the complexity of developmentally regulated gene expression in developing barley grains.

A new wild emmer wheat panel allows to map new loci associated with resistance to stem rust at seedling stage.

Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is a major wheat disease worldwide. A collection of 283 wild emmer wheat [Triticum turgidum L. subsp. dicoccoides (Körn. ex Asch. & Graebn.) Thell] accessions, representative of the entire Fertile Crescent region where wild emmer naturally occurs, was assembled, genotyped, and characterized for population structure, genetic diversity, and rate of linkage disequilibrium (LD) decay. Then, the collection was employed for mapping Pgt resistance genes, as a proof of concept of the effectiveness of genome-wide association studies in wild emmer. The collection was evaluated in controlled conditions for reaction to six common Pgt pathotypes (TPMKC, TTTTF, JRCQC, TRTTF, TTKSK/Ug99, and TKTTF). Most resistant accessions originated from the Southern Levant wild emmer lineage, with some showing a resistance reaction toward three to six tested races. Association analysis was conducted considering a 12K polymorphic single-nucleotide polymorphisms dataset, kinship relatedness between accessions, and population structure. Eleven significant marker-trait associations (MTA) were identified across the genome, which explained from 17% to up to 49% of phenotypic variance with an average 1.5 additive effect (based on the 1-9 scoring scale). The identified loci were either effective against single or multiple races. Some MTAs colocalized with known Pgt resistance genes, while others represent novel resistance loci useful for durum and bread wheat prebreeding. Candidate genes with an annotated function related to plant response to pathogens were identified at the regions linked to the resistance and defined according to the estimated small LD (about 126 kb), as typical of wild species.

Intact ribosomal DNA arrays of Potentilla origin detected in Erythronium nucleus suggest recent eudicot-to-monocot horizontal transfer.

During our initial phylogenetic study of the monocot genus Erythronium (Liliaceae), we observed peculiar eudicot-type internal transcribed spacer (ITS) sequences in a dataset derived from genomic DNA of Erythronium dens-canis. This raised the possibility of horizontal transfer of a eudicot alien ribosomal DNA (rDNA) into the Erythronium genome. In this work we aimed to support this hypothesis by carrying out genomic, molecular, and cytogenetic analyses. Genome skimming coupled by PacBio HiFi sequencing of a bacterial artificial chromosome clone derived from flow-sorted nuclei was used to characterise the alien 45S rDNA. Integration of alien rDNA in the recipient genome was further proved by Southern blotting and fluorescence in situ hybridization using specific probes. Alien rDNA, nested among Potentilla species in phylogenetic analysis, likely entered the Erythronium lineage in the common ancestor of E. dens-canis and E. caucasicum. Transferred eudicot-type rDNA preserved its tandemly arrayed feature on a single chromosome and was found to be transcribed in the monocot host, albeit much less efficiently than the native counterpart. This study adds a new example to the rarely documented nuclear-to-nuclear jumps of DNA between eudicots and monocots while holding the scientific community continually in suspense about the mode of DNA transfer.

Wild emmer wheat, the progenitor of modern bread wheat, exhibits great diversity in the VERNALIZATION1 gene.

Wild emmer wheat is an excellent reservoir of genetic variability that can be utilized to improve cultivated wheat to address the challenges of the expanding world population and climate change. Bearing this in mind, we have collected a panel of 263 wild emmer wheat (WEW) genotypes across the Fertile Crescent. The genotypes were grown in different locations and phenotyped for heading date. Genome-wide association mapping (GWAS) was carried out, and 16 SNPs were associated with the heading date. As the flowering time is controlled by photoperiod and vernalization, we sequenced the VRN1 gene, the most important of the vernalization response genes, to discover new alleles. Unlike most earlier attempts, which characterized known VRN1 alleles according to a partial promoter or intron sequences, we obtained full-length sequences of VRN-A1 and VRN-B1 genes in a panel of 95 wild emmer wheat from the Fertile Crescent and uncovered a significant sequence variation. Phylogenetic analysis of VRN-A1 and VRN-B1 haplotypes revealed their evolutionary relationships and geographic distribution in the Fertile Crescent region. The newly described alleles represent an attractive resource for durum and bread wheat improvement programs.

Contemplation on wheat vernalization.

Vernalization is a period of low non-freezing temperatures, which provides the competence to flower. This mechanism ensures that plants sown before winter develop reproductive organs in more favourable conditions during spring. Such an evolutionary mechanism has evolved in both monocot and eudicot plants. Studies in monocots, represented by temperate cereals like wheat and barley, have identified and proposed the VERNALIZATION1 (VRN1) gene as a key player in the vernalization response. VRN1 belongs to MADS-box transcription factors and is expressed in the leaves and the apical meristem, where it subsequently promotes flowering. Despite substantial research advancement in the last two decades, there are still gaps in our understanding of the vernalization mechanism. Here we summarise the present knowledge of wheat vernalization. We discuss VRN1 allelic variation, review vernalization models, talk VRN1 copy number variation and devernalization phenomenon. Finally, we suggest possible future directions of the vernalization research in wheat.

In-Depth Sequence Analysis of Bread Wheat VRN1 Genes.

The VERNALIZATION1 (VRN1) gene encodes a MADS-box transcription factor and plays an important role in the cold-induced transition from the vegetative to reproductive stage. Allelic variability of VRN1 homoeologs has been associated with large differences in flowering time. The aim of this study was to investigate the genetic variability of VRN1 homoeologs (VRN-A1, VRN-B1 and VRN-D1). We performed an in-depth sequence analysis of VRN1 homoeologs in a panel of 105 winter and spring varieties of hexaploid wheat. We describe the novel allele Vrn-B1f with an 836 bp insertion within intron 1 and show its specific expression pattern associated with reduced heading time. We further provide the complete sequence of the Vrn-A1b allele, revealing a 177 bp insertion in intron 1, which is transcribed into an alternative splice variant. Copy number variation (CNV) analysis of VRN1 homoeologs showed that VRN-B1 and VRN-D1 are present in only one copy. The copy number of recessive vrn-A1 ranged from one to four, while that of dominant Vrn-A1 was one or two. Different numbers of Vrn-A1a copies in the spring cultivars Branisovicka IX/49 and Bastion did not significantly affect heading time. We also report on the deletion of secondary structures (G-quadruplex) in promoter sequences of cultivars with more vrn-A1 copies.

A Panicum-derived chromosomal segment captured by Hordeum a few million years ago preserves a set of stress-related genes.

Intra-specific variability is a cornerstone of evolutionary success of species. Acquiring genetic material from distant sources is an important adaptive mechanism in bacteria, but it can also play a role in eukaryotes. In this paper, we investigate the nature and evolution of a chromosomal segment of panicoid (Poaceae, Panicoideae) origin occurring in the nuclear genomes of species of the barley genus Hordeum (Pooideae). The segment, spanning over 440 kb in the Asian Hordeum bogdanii and 219 kb in the South American Hordeum pubiflorum, resides on a pair of nucleolar organizer region (NOR)-bearing chromosomes. Conserved synteny and micro-collinearity of the segment in both species indicate a common origin of the segment, which was acquired before the split of the respective barley lineages 5-1.7 million years ago. A major part of the foreign DNA consists of several approximately 68 kb long repeated blocks containing five stress-related protein-coding genes and transposable elements (TEs). Whereas outside these repeats, the locus was invaded by multiple TEs from the host genome, the repeated blocks are rather intact and appear to be preserved. The protein-coding genes remained partly functional, as indicated by conserved reading frames, a low amount of non-synonymous mutations, and expression of mRNA. A screen across Hordeum species targeting the panicoid protein-coding genes revealed the presence of the genes in all species of the section Stenostachys. In summary, our study shows that grass genomes can contain large genomic segments obtained from distantly related species. These segments usually remain undetected, but they may play an important role in the evolution and adaptation of species.

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa.

BACKGROUND AND AIMS: Dioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY1Y2 sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes. METHODS: We flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons. KEY RESULTS: We identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes. CONCLUSIONS: The X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.

Identification of polycomb repressive complex 1 and 2 core components in hexaploid bread wheat.

BACKGROUND: Polycomb repressive complexes 1 and 2 play important roles in epigenetic gene regulation by posttranslationally modifying specific histone residues. Polycomb repressive complex 2 is responsible for the trimethylation of lysine 27 on histone H3; Polycomb repressive complex 1 catalyzes the monoubiquitination of histone H2A at lysine 119. Both complexes have been thoroughly studied in Arabidopsis, but the evolution of polycomb group gene families in monocots, particularly those with complex allopolyploid origins, is unknown. RESULTS: Here, we present the in silico identification of the Polycomb repressive complex 1 and 2 (PRC2, PRC1) subunits in allohexaploid bread wheat, the reconstruction of their evolutionary history and a transcriptional analysis over a series of 33 developmental stages. We identified four main subunits of PRC2 [E(z), Su(z), FIE and MSI] and three main subunits of PRC1 (Pc, Psc and Sce) and determined their chromosomal locations. We found that most of the genes coding for subunit proteins are present as paralogs in bread wheat. Using bread wheat RNA-seq data from different tissues and developmental stages throughout plant ontogenesis revealed variable transcriptional activity for individual paralogs. Phylogenetic analysis showed a high level of protein conservation among temperate cereals. CONCLUSIONS: The identification and chromosomal location of the Polycomb repressive complex 1 and 2 core components in bread wheat may enable a deeper understanding of developmental processes, including vernalization, in commonly grown winter wheat.

Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies.

Sex-chromosome systems tend to be highly conserved and knowledge about their evolution typically comes from macroevolutionary inference. Rapidly evolving complex sex-chromosome systems represent a rare opportunity to study the mechanisms of sex-chromosome evolution at unprecedented resolution. Three cryptic species of wood-white butterflies-Leptidea juvernica, L. sinapis and L. reali-have each a unique set of multiple sex-chromosomes with 3-4 W and 3-4 Z chromosomes. Using a transcriptome-based microarray for comparative genomic hybridisation (CGH) and a library of bacterial artificial chromosome (BAC) clones, both developed in L. juvernica, we identified Z-linked Leptidea orthologs of Bombyx mori genes and mapped them by fluorescence in situ hybridisation (FISH) with BAC probes on multiple Z chromosomes. In all three species, we determined synteny blocks of autosomal origin and reconstructed the evolution of multiple sex-chromosomes. In addition, we identified W homologues of Z-linked orthologs and characterised their molecular differentiation. Our results suggest that the multiple sex-chromosome system evolved in a common ancestor as a result of dynamic genome reshuffling through repeated rearrangements between the sex chromosomes and autosomes, including translocations, fusions and fissions. Thus, the initial formation of neo-sex chromosomes could not have played a role in reproductive isolation between these Leptidea species. However, the subsequent species-specific fissions of several neo-sex chromosomes could have contributed to their reproductive isolation. Then, significantly increased numbers of Z-linked genes and independent neo-W chromosome degeneration could accelerate the accumulation of genetic incompatibilities between populations and promote their divergence resulting in speciation.

Evolution of sex determination and heterogamety changes in section Otites of the genus Silene.

Switches in heterogamety are known to occur in both animals and plants. Although plant sex determination systems probably often evolved more recently than those in several well-studied animals, including mammals, and have had less time for switches to occur, we previously detected a switch in heterogamety in the plant genus Silene: section Otites has both female and male heterogamety, whereas S. latifolia and its close relatives, in a different section of the genus, Melandrium (subgenus Behenantha), all have male heterogamety. Here we analyse the evolution of sex chromosomes in section Otites, which is estimated to have evolved only about 0.55 MYA. Our study confirms female heterogamety in S. otites and newly reveals female heterogamety in S. borysthenica. Sequence analyses and genetic mapping show that the sex-linked regions of these two species are the same, but the region in S. colpophylla, a close relative with male heterogamety, is different. The sex chromosome pairs of S. colpophylla and S. otites each correspond to an autosome of the other species, and both differ from the XY pair in S. latifolia. Silene section Otites species are suitable for detailed studies of the events involved in such changes, and our phylogenetic analysis suggests a possible change from female to male heterogamety within this section. Our analyses suggest a possibility that has so far not been considered, change in heterogamety through hybridization, in which a male-determining chromosome from one species is introgressed into another one, and over-rides its previous sex-determining system.

Divergence between bread wheat and Triticum militinae in the powdery mildew resistance QPm.tut-4A locus and its implications for cloning of the resistance gene.

A segment of Triticum militinae chromosome 7G harbors a gene(s) conferring powdery mildew resistance which is effective at both the seedling and the adult plant stages when transferred into bread wheat (T. aestivum). The introgressed segment replaces a piece of wheat chromosome arm 4AL. An analysis of segregating materials generated to positionally clone the gene highlighted that in a plant heterozygous for the introgression segment, only limited recombination occurs between the introgressed region and bread wheat 4A. Nevertheless, 75 genetic markers were successfully placed within the region, thereby confining the gene to a 0.012 cM window along the 4AL arm. In a background lacking the Ph1 locus, the localized rate of recombination was raised 33-fold, enabling the reduction in the length of the region containing the resistance gene to a 480 kbp stretch harboring 12 predicted genes. The substituted segment in the reference sequence of bread wheat cv. Chinese Spring is longer (640 kbp) and harbors 16 genes. A comparison of the segments' sequences revealed a high degree of divergence with respect to both their gene content and nucleotide sequence. Of the 12 T. militinae genes, only four have a homolog in cv. Chinese Spring. Possible candidate genes for the resistance have been identified based on function predicted from their sequence.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome.

BACKGROUND: Numerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome. RESULTS: Using chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region. CONCLUSIONS: Sufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

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

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

Correction: Heritable heading time variation in wheat lines with the same number of Ppd-B1 gene copies.

[This corrects the article DOI: 10.1371/journal.pone.0183745.].

Heritable heading time variation in wheat lines with the same number of Ppd-B1 gene copies.

The ability of plants to identify an optimal flowering time is critical for ensuring the production of viable seeds. The main environmental factors that influence the flowering time include the ambient temperature and day length. In wheat, the ability to assess the day length is controlled by photoperiod (Ppd) genes. Due to its allohexaploid nature, bread wheat carries the following three Ppd-1 genes: Ppd-A1, Ppd-B1 and Ppd-D1. While photoperiod (in)sensitivity controlled by Ppd-A1 and Ppd-D1 is mainly determined by sequence changes in the promoter region, the impact of the Ppd-B1 alleles on the heading time has been linked to changes in the copy numbers (and possibly their methylation status) and sequence changes in the promoter region. Here, we report that plants with the same number of Ppd-B1 copies may have different heading times. Differences were observed among F7 lines derived from crossing two spring hexaploid wheat varieties. Several lines carrying three copies of Ppd-B1 headed 16 days later than other plants in the population with the same number of gene copies. This effect was associated with changes in the gene expression level and methylation of the Ppd-B1 gene.

Physical Map of the Short Arm of Bread Wheat Chromosome 3D.

Bread wheat ( L.) is one of the most important crops worldwide. Although a reference genome sequence would represent a valuable resource for wheat improvement through genomics-assisted breeding and gene cloning, its generation has long been hampered by its allohexaploidy, high repeat content, and large size. As a part of a project coordinated by the International Wheat Genome Sequencing Consortium (IWGSC), a physical map of the short arm of wheat chromosome 3D (3DS) was prepared to facilitate reference genome assembly and positional gene cloning. It comprises 869 contigs with a cumulative length of 274.5 Mbp and represents 85.5% of the estimated chromosome arm size. Eighty-six Mbp of survey sequences from chromosome arm 3DS were assigned in silico to physical map contigs via next-generation sequencing of bacterial artificial chromosome pools, thus providing a high-density framework for physical map ordering along the chromosome arm. About 60% of the physical map was anchored in this single experiment. Finally, 1393 high-confidence genes were anchored to the physical map. Comparisons of gene space of the chromosome arm 3DS with genomes of closely related species [ (L.) P.Beauv., rice ( L.), and sorghum [ (L.) Moench] and homeologous wheat chromosomes provided information about gene movement on the chromosome arm.

Multiple horizontal transfers of nuclear ribosomal genes between phylogenetically distinct grass lineages.

The movement of nuclear DNA from one vascular plant species to another in the absence of fertilization is thought to be rare. Here, nonnative rRNA gene [ribosomal DNA (rDNA)] copies were identified in a set of 16 diploid barley (Hordeum) species; their origin was traceable via their internal transcribed spacer (ITS) sequence to five distinct Panicoideae genera, a lineage that split from the Pooideae about 60 Mya. Phylogenetic, cytogenetic, and genomic analyses implied that the nonnative sequences were acquired between 1 and 5 Mya after a series of multiple events, with the result that some current Hordeum sp. individuals harbor up to five different panicoid rDNA units in addition to the native Hordeum rDNA copies. There was no evidence that any of the nonnative rDNA units were transcribed; some showed indications of having been silenced via pseudogenization. A single copy of a Panicum sp. rDNA unit present in H. bogdanii had been interrupted by a native transposable element and was surrounded by about 70 kbp of mostly noncoding sequence of panicoid origin. The data suggest that horizontal gene transfer between vascular plants is not a rare event, that it is not necessarily restricted to one or a few genes only, and that it can be selectively neutral.

Construction of BAC Libraries from Flow-Sorted Chromosomes.

Cloned DNA libraries in bacterial artificial chromosome (BAC) are the most widely used form of large-insert DNA libraries. BAC libraries are typically represented by ordered clones derived from genomic DNA of a particular organism. In the case of large eukaryotic genomes, whole-genome libraries consist of a hundred thousand to a million clones, which make their handling and screening a daunting task. The labor and cost of working with whole-genome libraries can be greatly reduced by constructing a library derived from a smaller part of the genome. Here we describe construction of BAC libraries from mitotic chromosomes purified by flow cytometric sorting. Chromosome-specific BAC libraries facilitate positional gene cloning, physical mapping, and sequencing in complex plant genomes.

Characterization of new allele influencing flowering time in bread wheat introgressed from Triticum militinae.

Flowering time variation was identified within a mapping population of doubled haploid lines developed from a cross between the introgressive line 8.1 and spring bread wheat cv. Tähti. The line 8.1 carried introgressions from tetraploid Triticum militinae in the cv. Tähti genetic background on chromosomes 1A, 2A, 4A, 5A, 7A, 1B and 5B. The most significant QTL for the flowering time variation was identified within the introgressed region on chromosome 5A and its largest effect was associated with the VRN-A1 locus, accounting for up to 70% of phenotypic variance. The allele of T. militinae origin was designated as VRN-A1f-like. The effect of the VRN-A1f-like allele was verified in two other mapping populations. QTL analysis identified that in cv. Tähti and cv. Mooni genetic background, VRN-A1f-like allele incurred a delay of 1.9-18.6 days in flowering time, depending on growing conditions. Sequence comparison of the VRN-A1f-like and VRN-A1a alleles from the parental lines of the mapping populations revealed major mutations in the promoter region as well as in the first intron, including insertion of a MITE element and a large deletion. The sequence variation allowed construction of specific diagnostic PCR markers for VRN-A1f-like allele determination. Identification and quantification of the effect of the VRN-A1f-like allele offers a useful tool for wheat breeding and for studying fine-scale regulation of flowering pathways in wheat.