Population structure and association analyses of the core collection of hexaploid accessions conserved ex situ in the Japanese gene bank NBRP-Wheat

Subgenome-informed statistical modeling of transcriptomes in 25 common wheat accessions reveals cis- and trans-regulation architectures

Common wheat is an allohexaploid, where it is difficult to obtain homoeolog-distinguished transcriptome data. Lasy-Seq, a type of 3' RNA-seq, is technology efficient at obtaining homoeolog-distinguished transcriptomes. Here, we applied Lasy-Seq to obtain transcriptome data from the seedlings, second leaves, and root tips of 25 common wheat lines mainly from East Asia. Roots and seedlings were similar to each other in transcriptome profiles, but they were different from the leaves. We then asked how three homoeologous genes from different subgenomes (i.e. triads) show different levels of expression. Specifically, we examined the effects of subgenomes, lines, and their interaction on the expression levels of each homoeolog triad, separately in each tissue. Of the 19 805 homoeolog triads, 51-55% showed significant effect of subgenome, suggesting cis-regulation, whereas 24-30% showed significant effect line, suggesting trans-regulation. We also found that 7.7-9.0% triads showed significant effects of the interaction. Hierarchical clustering and co-trans regulation network analysis of homoeolog triads revealed that the patterns of expression polymorphisms among the lines were shared in different genes. Our results also implied that expression variation between lines is caused by changes in a smaller number of common trans-factors. We performed gene ontology (GO)-term enrichment analysis using newly annotated and substantially improved GO annotations, which revealed that GO terms related to each tissue-type function were enriched in genes expressed in the leaves and roots. Our information provides fundamental knowledge for the future breeding of plants possessing complex gene regulatory networks such as common wheat.

A Catalog of GNI-A1 Genes That Regulate Floret Fertility in a Diverse Bread Wheat Collection

Modifying inflorescence architecture improves grain number and grain weight in bread wheat (Triticum aestivum). Allelic variation in Grain Number Increase 1 (GNI-A1) genes, encoding a homeodomain leucine zipper class I transcription factor, influences grain number and yield. However, allelic information about GNI-A1 in diverse germplasms remains limited. Here, we investigated GNI-A1 alleles in a panel of 252 diverse bread wheat accessions (NBRP core collection and HRO breeder's panel) by target resequencing. Cultivars carrying the reduced-function allele (105Y) were predominant in the NBRP panel, whereas the 105N functional allele was the major type in the HRO panel. Cultivars with the 105Y allele were distributed in Asian landraces but not in European genotypes. Association analysis demonstrated that floret fertility, together with grain size, were improved in cultivars in the NBRP core collection carrying the 105Y allele. These results imply that different alleles of GNI-A1 have been locally selected, with the 105Y allele selected in East Asia and the 105N allele selected in Europe.

Evaluation of the Allelic Variations in Vernalisation (VRN1) and Photoperiod (PPD1) Genes and Genetic Diversity in a Spanish Spelt Wheat Collection

Allelic variation within genes controlling the vernalisation requirement (VRN1) and photoperiod response (PPD1) determines the adaptation of wheat to different environmental growing conditions as well as influences other traits related to grain yield. This study aimed to screen a Spanish spelt wheat collection using gene-specific molecular markers for VRN-A1, VRN-B1, VRN-D1, and PPD-D1 loci and to phenotype for heading date (HD) in both field and greenhouse experiments under a long photoperiod and without vernalisation. Fifty-five spelt genotypes (91.7%) exhibited a spring growth habit, and all of them carried at least one dominant VRN1 allele, whereas five (8.3%) genotypes had a winter growth habit, and they carried the triple recessive allele combination. The Vrn-D1s was the most frequent allele in the studied set of spelt accessions, and it was found in combination with both the dominant Vrn-A1b and/or Vrn-B1a alleles in 88.3% of the spelt accessions tested. All spelt accessions carried the photoperiod-sensitive Ppd-D1b allele, which may explain the late heading of spelt germplasm compared to the commercial spring bread wheat Setenil used as a control. The least significant difference test showed significant differences between allelic combinations, the earliest accessions being those carrying two or three dominant alleles, followed by the one-gene combinations. In addition, the genetic diversity was evaluated through capillary electrophoresis using 15 wheat simple sequence repeat (SSR) markers. Most markers had high levels of polymorphism, producing 95 different alleles which ranged between 53 and 279 bp in size. Based on the polymorphic information content values obtained (from 0.51 to 0.97), 12 out of the 15 SSRs were catalogued as informative markers (values > 0.5). According to the dendrogram generated, the spelt accessions clustered as a separate group from the commercial bread wheat Setenil. Knowledge of VRN1 and PPD1 alleles, heading time, and genetic variability using SSR markers is valuable for spelt wheat breeding programs.

Natural variations of wheat EARLY FLOWERING 3 highlight their contributions to local adaptation through fine-tuning of heading time

KEY MESSAGE

We identified a large chromosomal deletion containing TaELF-B3 that confers early flowering in wheat. This allele has been preferred in recent wheat breeding in Japan to adapt to the environment. Heading at the appropriate time in each cultivation region can greatly contribute to stabilizing and maximizing yield. Vrn-1 and Ppd-1 are known as the major genes for vernalization requirement and photoperiod sensitivity in wheat. Genotype combinations of Vrn-1 and Ppd-1 can explain the variation in heading time. However, the genes that can explain the remaining variations in heading time are largely unknown. In this study, we aimed to identify the genes conferring early heading using doubled haploid lines derived from Japanese wheat varieties. Quantitative trait locus (QTL) analysis revealed a significant QTL on the long arm of chromosome 1B in multiple growing seasons. Genome sequencing using Illumina short reads and Pacbio HiFi reads revealed a large deletion of a ~ 500 kb region containing TaELF-B3, an orthologue of Arabidopsis clock gene EARLY FLOWERING 3 (ELF3). Plants with the deleted allele of TaELF-B3 (ΔTaELF-B3 allele) headed earlier only under short-day vernalization conditions. Higher expression levels of clock- and clock-output genes, such as Ppd-1 and TaGI, were observed in plants with the ΔTaELF-B3 allele. These results suggest that the deletion of TaELF-B3 causes early heading. Of the TaELF-3 homoeoalleles conferring early heading, the ΔTaELF-B3 allele showed the greatest effect on the early heading phenotype in Japan. The higher allele frequency of the ΔTaELF-B3 allele in western Japan suggests that the ΔTaELF-B3 allele was preferred during recent breeding to adapt to the environment. TaELF-3 homoeologs will help to expand the cultivated area by fine-tuning the optimal timing of heading in each environment.

Global genomic analyses of wheat powdery mildew reveal association of pathogen spread with historical human migration and trade

The fungus Blumeria graminis f. sp. tritici causes wheat powdery mildew disease. Here, we study its spread and evolution by analyzing a global sample of 172 mildew genomes. Our analyses show that B.g. tritici emerged in the Fertile Crescent during wheat domestication. After it spread throughout Eurasia, colonization brought it to America, where it hybridized with unknown grass mildew species. Recent trade brought USA strains to Japan, and European strains to China. In both places, they hybridized with local ancestral strains. Thus, although mildew spreads by wind regionally, our results indicate that humans drove its global spread throughout history and that mildew rapidly evolved through hybridization.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Diversity analysis of 80,000 wheat accessions reveals consequences and opportunities of selection footprints

Undomesticated wild species, crop wild relatives, and landraces represent sources of variation for wheat improvement to address challenges from climate change and the growing human population. Here, we study 56,342 domesticated hexaploid, 18,946 domesticated tetraploid and 3,903 crop wild relatives in a massive-scale genotyping and diversity analysis. Using DArTseq^TM technology, we identify more than 300,000 high-quality SNPs and SilicoDArT markers and align them to three reference maps: the IWGSC RefSeq v1.0 genome assembly, the durum wheat genome assembly (cv. Svevo), and the DArT genetic map. On average, 72% of the markers are uniquely placed on these maps and 50% are linked to genes. The analysis reveals landraces with unexplored diversity and genetic footprints defined by regions under selection. This provides fertile ground to develop wheat varieties of the future by exploring specific gene or chromosome regions and identifying germplasm conserving allelic diversity missing in current breeding programs.

Genebanks hold comprehensive collections of wild species, wild relatives, and landraces that are useful for genetic improvement. Here, the authors report the genotype of nearly 80,000 wheat accessions using DArTseq technology to show the less explored genetic diversity.

Variations in radioactive cesium accumulation in wheat germplasm from fields affected by the 2011 Fukushima nuclear power plant accident

Decreasing the transfer of radioactive cesium (RCs) from soil to crops has been important since the deposition of RCs in agricultural soil owing to the Fukushima nuclear power plant accident of 2011. We investigated the genotypic variation in RCs accumulation in 234 and 198 hexaploid wheat (Triticum spp.) varieties in an affected field in 2012 and 2013, respectively. The effects of soil exchangeable potassium (ExK) content to RCs accumulation in wheat varieties were also evaluated. A test field showed fourfold differences in soil ExK contents based on location, and the wheat varieties grown in areas with lower soil ExK contents tended to have higher grain RCs concentrations. RCs concentrations of shoots, when corrected by the soil ExK content, were positively significantly correlated between years, and RCs concentrations of shoots were significantly correlated with the grain RCs concentration corrected by the soil ExK content. These results indicated that there were genotypic variations in RCs accumulation. The grain to shoot ratio of RCs also showed significant genotypic variation. Wheat varieties with low RCs accumulations were identified. They could contribute to the research and breeding of low RCs accumulating wheat and to agricultural production in the area affected by RCs deposition.