Jan Šafář

The genus Hordeum, which includes cultivated barley and its wild relatives, represents a diverse group of grasses with agricultural, ecological, and evolutionary importance. Here, we present a genus-wide pangenome based on 25 chromosome-scale genome sequences encompassing all known diploid Hordeum species and selected tetraploids. Comparative analysis reveals extensive structural variations, including large-scale inversions and chromosomal rearrangements, which disrupt proximal chromosomal collinearity and reflect lineage-specific genome remodeling. Using single-copy gene phylogenies, we identify multiple deep interspecific hybridization events, segmental exchanges, and introgressions, particularly in regions enriched for conserved inversions. We also uncover widespread horizontal gene transfer from Panicoideae grasses into Hordeum, with functional integration of transferred genes. Transposable element dynamics, especially retrotransposon proliferation, account for much of the observed genome size variation and offer insights into polyploid genome evolution. Whole-genome alignments reveal contrasting evolutionary patterns: while coding genes remain highly conserved, regulatory elements show rapid turnover, providing a substrate for phenotypic diversification. Case studies of vernalization and glume development demonstrate how natural variation in regulatory and coding sequences contributes to flowering time and floral morphology. This comprehensive pangenomic framework establishes Hordeum as a model for studying genome evolution, structural variation, and the genetic basis of adaptive traits, offering new avenues for barley improvement.

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.

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.

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

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.

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

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 based on a BAC-based physical map, BAC pool paired end sequencing, chromosome arm specific mate-pair sequencing and Bionano optical mapping with the IWGSC 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 be now 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 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.

Switches in heterogamety occasionally occur both in animals and plants, although plant sex determination systems are mostly more recently evolved than those of animals, and have had less time for switches to occur. However, our previous research revealed a switch in heterogamety in section Otites of the plant genus Silene . Here we analyse in detail the evolution of genetic sex determination in section Otites , which is estimated to have evolved 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 , in a different section of the genus. Our phylogenetic analysis suggests a possible change from female to male heterogamety within Silene section Otites, making these species suitable for detailed studies of the events involved.

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.

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