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

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.

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.

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.

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.

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.

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.

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.

Whole-genome profiling and shotgun sequencing delivers an anchored, gene-decorated, physical map assembly of bread wheat chromosome 6A.

Bread wheat (Triticum aestivum L.) is the most important staple food crop for 35% of the world's population. International efforts are underway to facilitate an increase in wheat production, of which the International Wheat Genome Sequencing Consortium (IWGSC) plays an important role. As part of this effort, we have developed a sequence-based physical map of wheat chromosome 6A using whole-genome profiling (WGP™). The bacterial artificial chromosome (BAC) contig assembly tools fingerprinted contig (fpc) and linear topological contig (ltc) were used and their contig assemblies were compared. A detailed investigation of the contigs structure revealed that ltc created a highly robust assembly compared with those formed by fpc. The ltc assemblies contained 1217 contigs for the short arm and 1113 contigs for the long arm, with an L50 of 1 Mb. To facilitate in silico anchoring, WGP™ tags underlying BAC contigs were extended by wheat and wheat progenitor genome sequence information. Sequence data were used for in silico anchoring against genetic markers with known sequences, of which almost 79% of the physical map could be anchored. Moreover, the assigned sequence information led to the 'decoration' of the respective physical map with 3359 anchored genes. Thus, this robust and genetically anchored physical map will serve as a framework for the sequencing of wheat chromosome 6A, and is of immediate use for map-based isolation of agronomically important genes/quantitative trait loci located on this chromosome.

The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements.

BACKGROUND: The substantially large bread wheat genome, organized into highly similar three sub-genomes, renders genomic research challenging. The construction of BAC-based physical maps of individual chromosomes reduces the complexity of this allohexaploid genome, enables elucidation of gene space and evolutionary relationships, provides tools for map-based cloning, and serves as a framework for reference sequencing efforts. In this study, we constructed the first comprehensive physical map of wheat chromosome arm 5DS, thereby exploring its gene space organization and evolution. RESULTS: The physical map of 5DS was comprised of 164 contigs, of which 45 were organized into 21 supercontigs, covering 176 Mb with an N50 value of 2,173 kb. Fifty-eight of the contigs were larger than 1 Mb, with the largest contig spanning 6,649 kb. A total of 1,864 molecular markers were assigned to the map at a density of 10.5 markers/Mb, anchoring 100 of the 120 contigs (>5 clones) that constitute ~95 % of the cumulative length of the map. Ordering of 80 contigs along the deletion bins of chromosome arm 5DS revealed small-scale breaks in syntenic blocks. Analysis of the gene space of 5DS suggested an increasing gradient of genes organized in islands towards the telomere, with the highest gene density of 5.17 genes/Mb in the 0.67-0.78 deletion bin, 1.4 to 1.6 times that of all other bins. CONCLUSIONS: Here, we provide a chromosome-specific view into the organization and evolution of the D genome of bread wheat, in comparison to one of its ancestors, revealing recent genome rearrangements. The high-quality physical map constructed in this study paves the way for the assembly of a reference sequence, from which breeding efforts will greatly benefit.

A high-resolution physical map integrating an anchored chromosome with the BAC physical maps of wheat chromosome 6B.

BACKGROUND: A complete genome sequence is an essential tool for the genetic improvement of wheat. Because the wheat genome is large, highly repetitive and complex due to its allohexaploid nature, the International Wheat Genome Sequencing Consortium (IWGSC) chose a strategy that involves constructing bacterial artificial chromosome (BAC)-based physical maps of individual chromosomes and performing BAC-by-BAC sequencing. Here, we report the construction of a physical map of chromosome 6B with the goal of revealing the structural features of the third largest chromosome in wheat. RESULTS: We assembled 689 informative BAC contigs (hereafter reffered to as contigs) representing 91% of the entire physical length of wheat chromosome 6B. The contigs were integrated into a radiation hybrid (RH) map of chromosome 6B, with one linkage group consisting of 448 loci with 653 markers. The order and direction of 480 contigs, corresponding to 87% of the total length of 6B, were determined. We also characterized the contigs that contained a part of the nucleolus organizer region or centromere based on their positions on the RH map and the assembled BAC clone sequences. Analysis of the virtual gene order along 6B using the information collected for the integrated map revealed the presence of several chromosomal rearrangements, indicating evolutionary events that occurred on chromosome 6B. CONCLUSIONS: We constructed a reliable physical map of chromosome 6B, enabling us to analyze its genomic structure and evolutionary progression. More importantly, the physical map should provide a high-quality and map-based reference sequence that will serve as a resource for wheat chromosome 6B.

Flow sorting and sequencing meadow fescue chromosome 4F.

The analysis of large genomes is hampered by a high proportion of repetitive DNA, which makes the assembly of short sequence reads difficult. This is also the case in meadow fescue (Festuca pratensis), which is known for good abiotic stress resistance and has been used in intergeneric hybridization with ryegrasses (Lolium spp.) to produce Festulolium cultivars. In this work, we describe a new approach to analyze the large genome of meadow fescue, which involves the reduction of sample complexity without compromising information content. This is achieved by dissecting the genome to smaller parts: individual chromosomes and groups of chromosomes. As the first step, we flow sorted chromosome 4F and sequenced it by Illumina with approximately 50× coverage. This provided, to our knowledge, the first insight into the composition of the fescue genome, enabled the construction of the virtual gene order of the chromosome, and facilitated detailed comparative analysis with the sequenced genomes of rice (Oryza sativa), Brachypodium distachyon, sorghum (Sorghum bicolor), and barley (Hordeum vulgare). Using GenomeZipper, we were able to confirm the collinearity of chromosome 4F with barley chromosome 4H and the long arm of chromosome 5H. Several new tandem repeats were identified and physically mapped using fluorescence in situ hybridization. They were found as robust cytogenetic markers for karyotyping of meadow fescue and ryegrass species and their hybrids. The ability to purify chromosome 4F opens the way for more efficient analysis of genomic loci on this chromosome underlying important traits, including freezing tolerance. Our results confirm that next-generation sequencing of flow-sorted chromosomes enables an overview of chromosome structure and evolution at a resolution never achieved before.

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.

Chromosomes in the flow to simplify genome analysis.

Nuclear genomes of human, animals, and plants are organized into subunits called chromosomes. When isolated into aqueous suspension, mitotic chromosomes can be classified using flow cytometry according to light scatter and fluorescence parameters. Chromosomes of interest can be purified by flow sorting if they can be resolved from other chromosomes in a karyotype. The analysis and sorting are carried out at rates of 10(2)-10(4) chromosomes per second, and for complex genomes such as wheat the flow sorting technology has been ground-breaking in reducing genome complexity for genome sequencing. The high sample rate provides an attractive approach for karyotype analysis (flow karyotyping) and the purification of chromosomes in large numbers. In characterizing the chromosome complement of an organism, the high number that can be studied using flow cytometry allows for a statistically accurate analysis. Chromosome sorting plays a particularly important role in the analysis of nuclear genome structure and the analysis of particular and aberrant chromosomes. Other attractive but not well-explored features include the analysis of chromosomal proteins, chromosome ultrastructure, and high-resolution mapping using FISH. Recent results demonstrate that chromosome flow sorting can be coupled seamlessly with DNA array and next-generation sequencing technologies for high-throughput analyses. The main advantages are targeting the analysis to a genome region of interest and a significant reduction in sample complexity. As flow sorters can also sort single copies of chromosomes, shotgun sequencing DNA amplified from them enables the production of haplotype-resolved genome sequences. This review explains the principles of flow cytometric chromosome analysis and sorting (flow cytogenetics), discusses the major uses of this technology in genome analysis, and outlines future directions.

Functional features of a single chromosome arm in wheat (1AL) determined from its structure.

Bread wheat (Triticum aestivum L.) is one of the most important crops globally and a high priority for genetic improvement, but its large and complex genome has been seen as intractable to whole genome sequencing. Isolation of individual wheat chromosome arms has facilitated large-scale sequence analyses. However, so far there is no such survey of sequences from the A genome of wheat. Greater understanding of an A chromosome could facilitate wheat improvement and future sequencing of the entire genome. We have constructed BAC library from the long arm of T. aestivum chromosome 1A (1AL) and obtained BAC end sequences from 7,470 clones encompassing the arm. We obtained 13,445 (89.99%) useful sequences with a cumulative length of 7.57 Mb, representing 1.43% of 1AL and about 0.14% of the entire A genome. The GC content of the sequences was 44.7%, and 90% of the chromosome was estimated to comprise repeat sequences, while just over 1% encoded expressed genes. From the sequence data, we identified a large number of sites suitable for development of molecular markers (362 SSR and 6,948 ISBP) which will have utility for mapping this chromosome and for marker assisted breeding. From 44 putative ISBP markers tested 23 (52.3%) were found to be useful. The BAC end sequence data also enabled the identification of genes and syntenic blocks specific to chromosome 1AL, suggesting regions of particular functional interest and targets for future research.

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat.

BACKGROUND: Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat. RESULTS: We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus. CONCLUSION: Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.

A 3,000-loci transcription map of chromosome 3B unravels the structural and functional features of gene islands in hexaploid wheat.

To improve our understanding of the organization and regulation of the wheat (Triticum aestivum) gene space, we established a transcription map of a wheat chromosome (3B) by hybridizing a newly developed wheat expression microarray with bacterial artificial chromosome pools from a new version of the 3B physical map as well as with cDNA probes derived from 15 RNA samples. Mapping data for almost 3,000 genes showed that the gene space spans the whole chromosome 3B with a 2-fold increase of gene density toward the telomeres due to an increase in the number of genes in islands. Comparative analyses with rice (Oryza sativa) and Brachypodium distachyon revealed that these gene islands are composed mainly of genes likely originating from interchromosomal gene duplications. Gene Ontology and expression profile analyses for the 3,000 genes located along the chromosome revealed that the gene islands are enriched significantly in genes sharing the same function or expression profile, thereby suggesting that genes in islands acquired shared regulation during evolution. Only a small fraction of these clusters of cofunctional and coexpressed genes was conserved with rice and B. distachyon, indicating a recent origin. Finally, genes with the same expression profiles in remote islands (coregulation islands) were identified suggesting long-distance regulation of gene expression along the chromosomes in wheat.

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.