Manipulation and prediction of spike morphology traits for the improvement of grain yield in wheat

Discovering novel genomic regions explaining adaptation of bread wheat to conservation agriculture through GWAS

To sustainably increase wheat yield to meet the growing world population's food demand in the face of climate change, Conservation Agriculture (CA) is a promising approach. Still, there is a lack of genomic studies investigating the genetic basis of crop adaptation to CA. To dissect the genetic architecture of 19 morpho-physiological traits that could be involved in the enhanced adaptation and performance of genotypes under CA, we performed GWAS to identify MTAs under four contrasting production regimes viz., conventional tillage timely sown (CTTS), conservation agriculture timely sown (CATS), conventional tillage late sown (CTLS) and conservation agriculture late sown (CALS) using an association panel of 183 advanced wheat breeding lines along with 5 checks. Traits like Phi2 (Quantum yield of photosystem II; CATS:0.37, CALS: 0.31), RC (Relative chlorophyll content; CATS:55.51, CALS: 54.47) and PS1 (Active photosystem I centers; CATS:2.45, CALS: 2.23) have higher mean values in CA compared to CT under both sowing times. GWAS identified 80 MTAs for the studied traits across four production environments. The phenotypic variation explained (PVE) by these QTNs ranged from 2.15 to 40.22%. Gene annotation provided highly informative SNPs associated with Phi2, NPQ (Quantum yield of non-photochemical quenching), PS1, and RC which were linked with genes that play crucial roles in the physiological adaptation under both CA and CT. A highly significant SNP AX94651261 (9.43% PVE) was identified to be associated with Phi2, while two SNP markers AX94730536 (30.90% PVE) and AX94683305 (16.99% PVE) were associated with NPQ. Identified QTNs upon validation can be used in marker-assisted breeding programs to develop CA adaptive genotypes.

Grain yield trade-offs in spike-branching wheat can be mitigated by elite alleles affecting sink capacity and post-anthesis source activity

Introducing variations in inflorescence architecture, such as the 'Miracle-Wheat' (Triticum turgidum convar. compositum (L.f.) Filat.) with a branching spike, has relevance for enhancing wheat grain yield. However, in the spike-branching genotypes, the increase in spikelet number is generally not translated into grain yield advantage because of reduced grains per spikelet and grain weight. Here, we investigated if such trade-offs might be a function of source-sink strength by using 385 recombinant inbred lines developed by intercrossing the spike-branching landrace TRI 984 and CIRNO C2008, an elite durum (T. durum L.) cultivar; they were genotyped using the 25K array. Various plant and spike architectural traits, including flag leaf, peduncle, and spike senescence rate, were phenotyped under field conditions for 2 consecutive years. On chromosome 5AL, we found a new modifier QTL for spike branching, branched headt3 (bht-A3), which was epistatic to the previously known bht-A1 locus. Besides, bht-A3 was associated with more grains per spikelet and a delay in flag leaf senescence rate. Importantly, favourable alleles, viz. bht-A3 and grain protein content (gpc-B1) that delayed senescence, are required to improve grain number and grain weight in the spike-branching genotypes. In summary, achieving a balanced source-sink relationship might minimize grain yield trade-offs in Miracle-Wheat.

Identification of a novel SNP in the miR172 binding site of Q homoeolog AP2L-D5 is associated with spike compactness and agronomic traits in wheat (Triticum aestivum L.)

KEY MESSAGE

This study found that the compact spike locus of ANK-15 is on chromosome 5D instead of 2B. We have identified a new allele of AP2L-D5 as the candidate causal polymorphism. Spike architecture is a key determinant of wheat yield, a crop which supports much of the human diet but whose yield gains are stagnating. Spike architecture mutants offer opportunities to identify genetic factors contributing to inflorescence development. Here, we investigate the locus underlying the compact spike phenotype of mutant line ANK-15 by conducting mRNA-sequencing and genetic mapping using ANK-15 and its non-compact spike near-isogenic line Novosibirskaya 67 (N67). Previous literature has placed the compact spike locus of ANK-15 to chromosome 2B. However, based on the single nucleotide polymorphisms (SNPs) identified using mRNA-seq data, we were unable to detect polymorphisms between N67 and ANK-15 in the putative chromosome 2B region. We performed differential expression analysis of developing rachis and found that AP2L-D5, the D homoeolog of the domestication Q gene, is upregulated in ANK-15 in comparison to N67. ANK-15 carries a SNP in the microRNA172 binding site of AP2L-D5, which is predicted to lead to higher expression of AP2L-D5 due to decreased miRNA172-mediated degradation. Furthermore, we performed genetic mapping using an ANK-15 × N67 F2 population and found a single quantitative trait locus on chromosome 5D coinciding with the position of AP2L-D5. This result suggests that AP2L-D5 is likely the underlying causal gene for the compact spike phenotype in ANK-15. We performed a field trial to investigate the effect of the AP2L-D5 allele on agronomic traits and found that the AP2L-D5 allele from ANK-15 is associated with a significant reduction in height, increased thousand grain weight (TGW), and increased grain width.

Identification and validation of a major quantitative trait locus for spike length and compactness in the wheat (Triticum aestivum L.) line Chuanyu12D7

Spike length (SL) and spike compactness (SC) are crucial traits related to wheat (Triticum aestivum L.) yield potential. In this study, a backcrossed inbred lines (BILs) population segregating for SL/SC was developed by using a commercial variety chuanyu25 as recurrent parent and a backbone parent Chuanyu12D7. Bulked segregant analysis (BSA) combined with the Wheat 660K SNP array was performed to conduct quantitative trait locus (QTL) mapping. A major and stable SL/SC QTL (designated as QSl/Sc.cib-2D.1) was identified on chromosome 2DS, explaining 45.63-59.72% of the phenotypic variation. QSl/Sc.cib-2D.1 was mapped to a 102.29-Kb interval by flanking SNPs AX-110276364 and AX-111593853 using a BC₄F_2:3 population. Since QSl/Sc.cib-2D.1 is linked to the Rht8 gene, their additive effects on plant type and spike type were analysed. Remarkably, the superior allele of QSl/Sc.cib-2D.1 combined with Rht8 can increase SL and TGW, and decrese SC without any apparent trade-offs in other yield-related traits. In addition, the closely linked kompetitive allele-specific PCR (KASP) markers of this locus were developed for marker-assisted selection (MAS) breeding. Four genes within the physical interval were considered as potential candidates based on expression patterns as well as orthologous gene functions. These results laid the foundation for map-based cloning of the gene(s) underlying QSl/Sc.cib-2D.1 and its potential application in wheat ideotype breeding.

Multivariate analysis for yield and yield-related traits of amaranth genotypes from Ethiopia

The genus Amaranthus is one of the few dicotyledonous, non-grass mesophytes that use specialized C4 annuals or short-lived perennials to produce significant amounts of edible small-seeded pseudo cereals. In this study, we characterized the genetic diversity of 120 genotypes of amaranths collected from diverse amaranth-growing regions of Ethiopia using multivariate analysis of yield and yield-related traits. The experiments were carried out at Hawassa University, in the years 2020 and 2021. The experimental design was set up using an alpha lattice design and replicated two times. The collected data were examined for 24 descriptors. Principal component analysis showed that the first six principal components with eigenvalues greater than one contributed 80.41% of the variability. However, the first two principal components explained 52.42% of the total variation. The highest contributing traits in the first component were days to flowering, basal stem diameter, plant height at flowering, plant height at maturity, auxiliary inflorescence length, number of branches, terminal inflorescence lateral length, days to maturity, terminal inflorescence stalk length, leaf number, leaf length, top lateral branch length. The traits with the greatest weight on the second component were leaf area, basal lateral branch length, leaf length, and leaf width, grain filling period, grain sinking filling rate, and grain yield. Therefore, selection based on these traits would be effective for yield improvement in amaranth genotypes. Additionally, the hierarchical clustering grouped all the genotypes into five clusters. The pairwise generalized squared distance (D²) among the five clusters based on Mahalanobis's D² statistics revealed the maximum and highly significant genetic distance was observed between II and III (277.79), while the minimum inter-cluster distance observed between clusters I and II (39.50). The findings suggest that amaranth genotypes in Ethiopia have a lot of genetic variation, which might be used for future breeding and ought to be conserved.

Physical mapping of the wheat genes in low-recombination regions: radiation hybrid mapping of the C-locus

KEY MESSAGE

This work reports the physical mapping of an important gene affecting spike compactness located in a low-recombination region of hexaploid wheat. This work paves the way for the eventual isolation and characterization of the factor involved but also opens up possibilities to use this approach to precisely map other wheat genes located on proximal parts of wheat chromosomes that show highly reduced recombination. Mapping wheat genes, in the centromeric and pericentromeric regions (~ 2/3rd of a given chromosome), poses a formidable challenge due to highly suppressed recombination. Using an example of compact spike locus (C-locus), this study provides an approach to precisely map wheat genes in the pericentromeric and centromeric regions that house ~ 30% of wheat genes. In club-wheat, spike compactness is controlled by the dominant C-locus, but previous efforts have failed to localize it, on a particular arm of chromosome 2D. We integrated radiation hybrid (RH) and high-resolution genetic mapping to locate C-locus on the short arm of chromosome 2D. Flanking markers of the C-locus span a physical distance of 11.0 Mb (231.0-242 Mb interval) and contain only 11 high-confidence annotated genes. This work demonstrates the value of this integrated strategy in mapping dominant genes in the low-recombination regions of the wheat genome. A comparison of the mapping resolutions of the RH and genetic maps using common anchored markers indicated that the RH map provides ~ 9 times better resolution that the genetic map even with much smaller population size. This study provides a broadly applicable approach to fine map wheat genes in regions of suppressed recombination.

Identification of novel putative alleles related to important agronomic traits of wheat using robust strategies in GWAS

Principal component analysis (PCA) is widely used in various genetics studies. In this study, the role of classical PCA (cPCA) and robust PCA (rPCA) was evaluated explicitly in genome-wide association studies (GWAS). We evaluated 294 wheat genotypes under well-watered and rain-fed, focusing on spike traits. First, we showed that some phenotypic and genotypic observations could be outliers based on cPCA and different rPCA algorithms (Proj, Grid, Hubert, and Locantore). Hubert's method provided a better approach to identifying outliers, which helped to understand the nature of these samples. These outliers led to the deviation of the heritability of traits from the actual value. Then, we performed GWAS with 36,000 single nucleotide polymorphisms (SNPs) based on the traditional approach and two robust strategies. In the conventional approach and using the first three components of cPCA as population structure, 184 and 139 marker-trait associations (MTAs) were identified for five traits in well-watered and rain-fed environments, respectively. In the first robust strategy and when rPCA was used as population structure in GWAS, we observed that the Hubert and Grid methods identified new MTAs, especially for yield and spike weight on chromosomes 7A and 6B. In the second strategy, we followed the classical and robust principal component-based GWAS, where the first two PCs obtained from phenotypic variables were used instead of traits. In the recent strategy, despite the similarity between the methods, some new MTAs were identified that can be considered pleiotropic. Hubert's method provided a better linear combination of traits because it had the most MTAs in common with the traditional approach. Newly identified SNPs, including rs19833 (5B) and rs48316 (2B), were annotated with important genes with vital biological processes and molecular functions. The approaches presented in this study can reduce the misleading GWAS results caused by the adverse effect of outlier observations.

A framework for improving wheat spike development and yield based on the master regulatory TOR and SnRK gene systems

The low rates of yield gain in wheat breeding programs create an ominous situation for the world. Amongst the reasons for this low rate are issues manifested in spike development that result in too few spikelets, fertile florets, and therefore grains being produced. Phases in spike development are particularly sensitive to stresses of various kinds and origins, and these are partly responsible for the deficiencies in grain production and slow rates of gain in yield. The diversity of developmental processes, stresses, and the large numbers of genes involved make it particularly difficult to prioritize approaches in breeding programs without an overarching, mechanistic framework. Such a framework, introduced here, is provided around the master regulator target of rapamycin and sucrose non-fermenting-1 (SNF1)-related protein kinase complexes and their control by trehalose-6-phosphate and other molecules. Being master regulators of the balance between growth and growth inhibition under stress, these provide genetic targets for creating breakthroughs in yield enhancement. Examples of potential targets and experimental approaches are described.

Genome-wide association mapping of bread wheat genotypes using yield and grain morphology-related traits under different environments

Depleting water resources and increasing global temperature due to climate change are major challenges to agriculture and food security worldwide. Deciphering the underlying mechanisms of traits contributing to grain development and yield is essential for the development of climate-resilient cultivars. Therefore, this study assessed 105 bread wheat genotypes grown under control, drought, and heat-stress conditions for two crop seasons and performed a genome-wide association study (GWAS) using a 90k SNP array. The genotypes showed significant trait differences under all environmental conditions. Highly significant variation was observed, with moderate (50.09%) to high (76.19%) heritability in the studied germplasms. The studied traits were all also significantly positively correlated. A total of 541 significant associations (p ≤ 10^-3) between marker and trait (MTAs) were observed after crossing the FDR <0.05 threshold for all traits. Among these, 195, 179, and 167 significant MTAs were detected under control, drought, and heat-stress conditions, respectively. Under the control and drought conditions, pleiotropic loci BS00010616_51 and BS00010868_51 were observed on chromosomes 7B and 1B situated at 186.24 cM and 35.47 cM, respectively. Pleiotropic loci BS00010868_51, Kukri_c11154_1723, and Ex_c10068_1509 were identified on chromosomes 1B, 5B, and 2A, respectively, under control and heat stress conditions. A stable and consistent locus (Excalibur_c20796_395) on chromosome 7A, located at 372.34 cM, was also linked to grain morphology and yield-related attributes in control, drought, and heat-stress conditions. The results of the current study confirmed several previously reported MTAs for the traits under consideration and identified new MTAs under harsh climatic conditions. These SNPs will aid in the discovery of novel genes in wheat. SNPs showing significant associations may be used in marker-assisted selection and allow the development of drought- and heat-tolerant genotypes with high yields to address global food security concerns.

Genetic dissection of the roles of β-hydroxylases in carotenoid metabolism, photosynthesis, and plant growth in tetraploid wheat (Triticum turgidum L.)

Functional redundancy and subfunctionalization of β-hydroxylases in tetraploid wheat tissues open up opportunities for manipulation of carotenoid metabolism for trait improvement. The genetic diversity provided by subgenome homoeologs in allopolyploid wheat can be leveraged for developing improved wheat varieties with modified chemical traits, including profiles of carotenoids, which play critical roles in photosynthesis, photoprotection, and growth regulation. Carotenoid profiles are greatly influenced by hydroxylation catalyzed by β-hydroxylases (HYDs). To genetically dissect the contribution of HYDs to carotenoid metabolism and wheat growth and yield, we isolated loss-of-function mutants of the two homoeologs of HYD1 (HYD-A1 and HYD-B1) and HYD2 (HYD-A2 and HYD-B2) from the sequenced ethyl methanesulfonate mutant population of the tetraploid wheat cultivar Kronos, and generated various mutant combinations. Although functional redundancy between HYD1 and HYD2 paralogs was observed in leaves, HYD1 homoeologs contributed more than HYD2 homoeologs to carotenoid β-ring hydroxylation in this tissue. By contrast, the HYD2 homoeologs functioned toward production of lutein, the major carotenoid in mature grains, whereas HYD1 homoeologs had a limited role. These results collectively suggested subfunctionalization of HYD genes and homoeologs in different tissues of tetraploid wheat. Despite reduced photoprotective responses observed in the triple hyd-A1 hyd-B1 hyd-A2 and the quadruple hyd-A1 hyd-B1 hyd-A2 hyd-B2 combinatorial mutants, comprehensive plant phenotyping analysis revealed that all mutants analyzed were comparable to the control for growth, yield, and fertility, except for a slight delay in anthesis and senescence as well as accelerated germination in the quadruple mutant. Overall, this research takes steps toward untangling the functions of HYDs in wheat and has implications for improving performance and consumer traits of this economically important global crop.

Genome-Wide Association Analysis for Hybrid Breeding in Wheat

Disclosure of markers that are significantly associated with plant traits can help develop new varieties with desirable properties. This study determined the genome-wide associations based on DArTseq markers for six agronomic traits assessed in eight environments for wheat. Moreover, the association study for heterosis and analysis of the effects of markers grouped by linkage disequilibrium were performed based on mean values over all experiments. All results were validated using data from post-registration trials. GWAS revealed 1273 single nucleotide polymorphisms with biologically significant effects. Most polymorphisms were predicted to be modifiers of protein translation, with only two having a more pronounced effect. Markers significantly associated with the considered set of features were clustered within chromosomes based on linkage disequilibrium in 327 LD blocks. A GWAS for heterosis revealed 1261 markers with significant effects.

Automatic Detection and Counting of Wheat Spikelet Using Semi-Automatic Labeling and Deep Learning

The number of wheat spikelets is an important phenotypic trait and can be used to assess the grain yield of the wheat crop. However, manual counting of spikelets is time-consuming and labor-intensive. To develop a cost-effective and highly efficient phenotyping system for counting the number of spikelets under laboratory conditions, methods based on imaging processing techniques and deep learning were proposed to accurately detect and count spikelets from color images of wheat spikes captured at the grain filling stage. An unsupervised learning-based method was first developed to automatically detect and label spikelets from spike color images and build the datasets for the model training. Based on the constructed datasets, a deep convolutional neural network model was retrained using transfer learning to detect the spikelets. Testing results showed that the root mean squared errors, relative root mean squared errors, and the coefficients of determination between the automatic and manual counted spikelets for four wheat lines were 0.62, 0.58, 0.54, and 0.77; 3.96, 3.73, 3.34, and 4.94%; and 0.73, 0.78, 0.84, and 0.67, respectively. We demonstrated that the proposed methods can effectively estimate the number of wheat spikelets, which improves the counting efficiency of wheat spikelets and contributes to the analysis of the developmental characteristics of wheat spikes.

InDels Identification and Association Analysis with Spike and Awn Length in Chinese Wheat Mini-Core Collection

Diversity surveys of germplasm are important for gaining insight into the genomic basis for crop improvement; especially InDels, which are poorly understood in hexaploid common wheat. Here, we describe a map of 89,923 InDels from exome sequencing of 262 accessions of a Chinese wheat mini-core collection. Population structure analysis, principal component analysis and selective sweep analysis between landraces and cultivars were performed. Further genome-wide association study (GWAS) identified five QTL (Quantitative Trait Loci) that were associated with spike length, two of them, on chromosomes 2B and 6A, were detected in 10 phenotypic data sets. Assisted with RNA-seq data, we identified 14 and 21 genes, respectively that expressed in spike and rachis within the two QTL regions that can be further investigated for candidate genes discovery. Moreover, InDels were found to be associated with awn length on chromosomes 5A, 6B and 4A, which overlapped with previously reported genetic loci B1 (Tipped 1), B2 (Tipped 2) and Hd (Hooded). One of the genes TaAGL6 that was previously shown to affect floral organ development was found at the B2 locus to affect awn length development. Our study shows that trait-associated InDels may contribute to wheat improvement and may be valuable molecular markers for future wheat breeding.

A Multi-Omics Approach for Rapid Identification of Large Genomic Lesions at the Wheat Dense Spike (wds) Locus

Optimal spike architecture provides a favorable structure for grain development and yield improvement. However, the number of genes cloned to underlie wheat spike architecture is extremely limited. Here, we obtained a wheat dense spike mutant (wds) induced by 60Co treatment of a common wheat landrace Huangfangzhu that exhibited significantly reduced spike and grain lengths. The shortened spike length was caused by longitudinal reduction in number and length of rachis cells. We adopted a multi-omics approach to identify the genomic locus underlying the wds mutant. We performed Exome Capture Sequencing (ECS) and identified two large deletion segments, named 6BL.1 at 334.8∼424.3 Mb and 6BL.2, 579.4∼717.8 Mb in the wds mutant. RNA-seq analysis confirmed that genes located in these regions lost their RNA expression. We then found that the 6BL.2 locus was overlapping with a known spike length QTL, qSL6B.2. Totally, 499 genes were located within the deleted region and two of them were found to be positively correlated with long spike accessions but not the ones with short spike. One of them, TraesCS6B01G334600, a well-matched homolog of the rice OsBUL1 gene that works in the Brassinosteroids (BR) pathway, was identified to be involved in cell size and number regulation. Further transcriptome analysis of young spikes showed that hormone-related genes were enriched among differentially expressed genes, supporting TraesCS6B01G334600 as a candidate gene. Our work provides a strategy to rapid locate genetic loci with large genomic lesions in wheat and useful resources for future wheat study.

TaCol-B5 modifies spike architecture and enhances grain yield in wheat

Spike architecture influences grain yield in wheat. We report the map-based cloning of a gene determining the number of spikelet nodes per spike in common wheat. The cloned gene is named TaCOL-B5 and encodes a CONSTANS-like protein that is orthologous to COL5 in plant species. Constitutive overexpression of the dominant TaCol-B5 allele but without the region encoding B-boxes in a common wheat cultivar increases the number of spikelet nodes per spike and produces more tillers and spikes, thereby enhancing grain yield in transgenic plants under field conditions. Allelic variation in TaCOL-B5 results in amino acid substitutions leading to differential protein phosphorylation by the protein kinase TaK4. The TaCol-B5 allele is present in emmer wheat but is rare in a global collection of modern wheat cultivars.

Raise and characterization of a bread wheat hybrid line (Tulaykovskaya 10 × Saratovskaya 29) with chromosome 6Agi2 introgressed from Thinopyrum intermedium

Wheatgrass Thinopyrum intermedium is a source of agronomically valuable traits for common wheat. Partial wheat–wheatgrass amphidiploids and lines with wheatgrass chromosome substitutions are extensively used as intermediates in breeding programs. Line Agis 1 (6Agi2/6D) is present in the cultivar Tulaykovskaya 10 pedigree. Wheatgrass chromosome 6Agi2 carries multiple resistance to fungal diseases in various ecogeographical zones. In this work, we studied the transfer of chromosome 6Agi2 in hybrid populations Saratovskaya 29 × skaya 10 (S29 × T10) and Tulaykovskaya 10 × Saratovskaya 29 (T10 × S29). Chromosome 6Agi2 was identif ied by PCR with chromosome-specif ic primers and by genomic in situ hybridization (GISH). According to molecular data, 6Agi2 was transmitted to nearly half of the plants tested in the F2 and F3 generations. A new breeding line 49-14 (2n = 42) with chromosome pair 6Agi2 was isolated and characterized in T10 × S29 F5 by GISH. According to the results of our f ield experiment in 2020, the line had high productivity traits. The grain weights per plant (10.04 ± 0.93 g) and the number of grains per plant (259.36 ± 22.49) did not differ signif icantly from the parent varieties. The number of grains per spikelet in the main spike was signif icantly higher than in S29 ( p ≤ 0.001) or T10 ( p ≤ 0.05). Plants were characterized by the ability to set 3.77 ± 0.1 grains per spikelet, and this trait varied among individuals from 2.93 to 4.62. The grain protein content was 17.91 %, and the gluten content, 40.55 %. According to the screening for fungal disease resistance carried out in the f ield in 2018 and 2020, chromosome 6Agi2 makes plants retain immunity to the West Siberian population of brown rust and to dominant races of stem rust. It also provides medium resistant and medium susceptible types of response to yellow rust. The possibility of using lines/varieties of bread wheat with wheatgrass chromosomes 6Agi2 in breeding in order to increase protein content in the grain, to confer resistance to leaf diseases on plants and to create multif lowered forms is discussed.

Strategies of grain number determination differentiate barley row types

Gaining knowledge on fundamental interactions of various yield components is crucial to improve yield potential in small grain cereals. It is well known in barley that increasing grain number greatly improves yield potential; however, the yield components determining grain number and their association in barley row types are less explored. In this study, we assessed different yield components such as potential spikelet number (PSN), spikelet survival (SSL), spikelet number (SN), grain set (GS), and grain survival (GSL), as well as their interactions with grain number by using a selected panel of two- and six-rowed barley types. Also, to analyze the stability of these interactions, we performed the study in the greenhouse and the field. From this study, we found that in two-rowed barley, grain number determination is strongly influenced by PSN rather than SSL and/or GS in both growth conditions. Conversely, in six-rowed barley, grain number is associated with SSL instead of PSN and/or GS. Thus, our study showed that increasing grain number might be possible by augmenting PSN in two-rowed genotypes, while for six-rowed genotypes SSL needs to be improved. We speculate that this disparity of grain number determination in barley row types might be due to the fertility of lateral spikelets. Collectively, this study revealed that grain number in two-rowed barley largely depends on the developmental trait, PSN, while in six-rowed barley, it mainly follows the ability for SSL.

Towards Automated Analysis of Grain Spikes in Greenhouse Images Using Neural Network Approaches: A Comparative Investigation of Six Methods

Automated analysis of small and optically variable plant organs, such as grain spikes, is highly demanded in quantitative plant science and breeding. Previous works primarily focused on the detection of prominently visible spikes emerging on the top of the grain plants growing in field conditions. However, accurate and automated analysis of all fully and partially visible spikes in greenhouse images renders a more challenging task, which was rarely addressed in the past. A particular difficulty for image analysis is represented by leaf-covered, occluded but also matured spikes of bushy crop cultivars that can hardly be differentiated from the remaining plant biomass. To address the challenge of automated analysis of arbitrary spike phenotypes in different grain crops and optical setups, here, we performed a comparative investigation of six neural network methods for pattern detection and segmentation in RGB images, including five deep and one shallow neural network. Our experimental results demonstrate that advanced deep learning methods show superior performance, achieving over 90% accuracy by detection and segmentation of spikes in wheat, barley and rye images. However, spike detection in new crop phenotypes can be performed more accurately than segmentation. Furthermore, the detection and segmentation of matured, partially visible and occluded spikes, for which phenotypes substantially deviate from the training set of regular spikes, still represent a challenge to neural network models trained on a limited set of a few hundreds of manually labeled ground truth images. Limitations and further potential improvements of the presented algorithmic frameworks for spike image analysis are discussed. Besides theoretical and experimental investigations, we provide a GUI-based tool (SpikeApp), which shows the application of pre-trained neural networks to fully automate spike detection, segmentation and phenotyping in images of greenhouse-grown plants.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems

Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.

Genetic analysis of morphological traits of the spike and reproductivity elements of speltoid chemomutant Triticum aestivum

Long use of the methods of direct intraspecies hybridization in the selective breeding of varieties of Triticum aestivum L. has led to narrowing of their gene fund and close similarity of the genetic potential of their selective breeding traits. Using the method of experimental mutagenesis, one can induce systemic mutants with features of other hexaploid Triticum species which can freely cross-breed with other hybridization offspring, contributing to extension of genetic potential of cultivated wheat and creation of varieties with new levels of manifestation of economically beneficial traits. We studied the pattern of inheritance of morphological traits of the ear, length of the stem and the elements of productivity of speltoid chemomutant of T. aestivum. For the genetic analysis we used hybrids F1 and F2 of soft wheat obtained by cross-breeding speltoid macromutant (Smuhlianka speltoid), induced by the impact of aqueous solution of N-nitroso-N-methylurea (NMU) in the concentration of 0.025% on the seeds of Smuhlianka variety, with plants of Smuhlianka variety (Erythrospermum variety) and Podoloanka (Lutescens variety). To determine the pattern of inheritance the spike morphology, length of the stem and the elements of productivity in F1 hybrids, we calculated the extent of phenotype domination. In populations of F2, we examined plants with different combinations of phenotype manifestation of ear morphology. In F1 hybrids, the speltoid shape of the ear, absence of awns and red colour of the glumes indicated the dominant pattern of inheritance. The high level of phenotype domination of length of the stem and ear, number of spikelets in the main ear indicates the inheritance of the features according to intermediate, partly dominant and over-dominant types. Taking into account the segregation according to the features of spike morphology, awnedness and colour of glumes, the plants of F2 population were divided into phenotype classes, that is 12 and 6 in the combinations of respectively Smuhlianka speltoid × Smuhlianka and Smuhlianka speltoid × Podolianka. We determined that the obtained results are the consequences of dihybrid linkage which corresponds to the theoretical proportion of 12 : 3 : 1. Segregation into non-aristate and aristate plants corresponds to the proportion of monogene segregation of 3 : 1. Within separately distinguished phenotype classes, no independent inheritance of the shape of the ear and awnedness was observed. Dihybrid segregation of F2 plants into speltoid, squarehead and varieties Lutescence/Erythrospermum with quantitative superiority of speltoid plants suggests the control of the trait by two non-allele genes with epistatic interaction. The red colour of the glumes indicates the dominant monogenic pattern of inheritance. Absence of independent inheritance of the shape of the ear and awnedness indicates localization of genes which determine these features in one chromosome.

Detection of genomic regions associated with tiller number in Iranian bread wheat under different water regimes using genome-wide association study

Two of the important traits for wheat yield are tiller and fertile tiller number, both of which have been thought to increase cereal yield in favorable and unfavorable environments. A total of 6,349 single nucleotide polymorphism (SNP) markers from the 15 K wheat Infinium array were employed for genome-wide association study (GWAS) of tillering number traits, generating a physical distance of 14,041.6 Mb based on the IWGSC wheat genome sequence. GWAS analysis using Fixed and random model Circulating Probability Unification (FarmCPU) identified a total of 47 significant marker-trait associations (MTAs) for total tiller number (TTN) and fertile tiller number (FTN) in Iranian bread wheat under different water regimes. After applying a 5% false discovery rate (FDR) threshold, a total of 13 and 11 MTAs distributed on 10 chromosomes were found to be significantly associated with TTN and FTN, respectively. Linked single nucleotide polymorphisms for IWB39005 (2A) and IWB44377 (7A) were highly significantly associated (FDR < 0.01) with TTN and FTN traits. Moreover, to validate GWAS results, meta-analysis was performed and 30 meta-QTL regions were identified on 11 chromosomes. The integration of GWAS and meta-QTLs revealed that tillering trait in wheat is a complex trait which is conditioned by the combined effects of minor changes in multiple genes. The information provided by this study can enrich the currently available candidate genes and genetic resources pools, offering evidence for subsequent analysis of genetic adaptation of wheat to different climatic conditions of Iran and other countries.

Quantitative trait loci for agronomic traits in tetraploid wheat for enhancing grain yield in Kazakhstan environments

Durum wheat (Triticum turgidum L. ssp. durum) is one of the top crops in Kazakhstan, where it is cultivated in different ecological niches, mainly at higher latitudes in the steppe zone of the northern region. Therefore, local breeding programs for durum wheat are primarily focused on selection for high productivity in Northern Kazakhstan based on the introduction of promising foreign germplasm and the adoption of marker-assisted selection. In this study, a world tetraploid wheat collection consisted of 184 primitive and domesticated accessions, which were previously genotyped using 16,425 polymorphic SNP markers, was field-tested in Northern and South-eastern Kazakhstan. The field tests have allowed the identification of 80 durum wheat promising lines in Northern Kazakhstan in comparison with a local standard cultivar. Also, GGE (Genotype and Genotype by Environment) biplot analyses for yield performance revealed that accessions of T. dicoccum, T. carthlicum, and T. turanicum also have potential to improve durum wheat yield in the region. The genome-wide association study (GWAS) has allowed the identification of 83 MTAs (marker-trait associations) for heading date, seed maturation time, plant height, spike length, number of fertile spikes, number of kernels per spike, and thousand kernel weight. The comparison of the 83 identified MTAs with those previously reported in GWAS for durum wheat suggests that 38 MTAs are presumably novel, while the co-localization of a large number of MTAs with those previously published confirms the validity of the results of this study. The MTAs reported herewith will provide the opportunity to implement marker-assisted selection in ongoing durum wheat breeding projects targeting higher productivity in the region.

Molecular Markers Associated with Agro-Physiological Traits under Terminal Drought Conditions in Bread Wheat

Terminal drought stress poses a big challenge to sustain wheat grain production in rain-fed environments. This study aimed to utilize the genetically diverse pre-breeding lines for identification of genomic regions associated with agro-physiological traits at terminal stage drought stress in wheat. A total of 339 pre-breeding lines panel derived from three-way crosses of 'exotics × elite × elite' lines were evaluated in field conditions at Obregon, Mexico for two years under well irrigated as well as drought stress environments. Drought stress was imposed at flowering by skipping the irrigations at pre and post anthesis stage. Results revealed that drought significantly reduced grain yield (Y), spike length (SL), number of grains spikes-1 (NGS) and thousand kernel weight (TKW), while kernel abortion (KA) was increased. Population structure analysis in this panel uncovered three sub-populations. Genome wide linkage disequilibrium (LD) decay was observed at 2.5 centimorgan (cM). The haplotypes-based genome wide association study (GWAS) identified significant associations of Y, SL, and TKW on three chromosomes; 4A (HB10.7), 2D (HB6.10) and 3B (HB8.12), respectively. Likewise, associations on chromosomes 6B (HB17.1) and 3A (HB7.11) were found for NGS while on chromosome 3A (HB7.12) for KA. The genomic analysis information generated in the study can be efficiently utilized to improve Y and/or related parameters under terminal stage drought stress through marker-assisted breeding.

Plant Inflorescence Architecture: The Formation, Activity, and Fate of Axillary Meristems

The above-ground plant body in different plant species can have very distinct forms or architectures that arise by recurrent redeployment of a finite set of building blocks-leaves with axillary meristems, stems or branches, and flowers. The unique architectures of plant inflorescences in different plant families and species, on which this review focuses, determine the reproductive success and yield of wild and cultivated plants. Major contributors to the inflorescence architecture are the activity and developmental trajectories adopted by axillary meristems, which determine the degree of branching and the number of flowers formed. Recent advances in genetic and molecular analyses in diverse flowering plants have uncovered both common regulatory principles and unique players and/or regulatory interactions that underlie inflorescence architecture. Modulating activity of these regulators has already led to yield increases in the field. Additional insight into the underlying regulatory interactions and principles will not only uncover how their rewiring resulted in altered plant form, but will also enhance efforts at optimizing plant architecture in desirable ways in crop species.

TaAPO-A1, an ortholog of rice ABERRANT PANICLE ORGANIZATION 1, is associated with total spikelet number per spike in elite European hexaploid winter wheat (Triticum aestivum L.) varieties

We dissected the genetic basis of total spikelet number (TSN) along with other traits, viz. spike length (SL) and flowering time (FT) in a panel of 518 elite European winter wheat varieties. Genome-wide association studies (GWAS) based on 39,908 SNP markers revealed highly significant quantitative trait loci (QTL) for TSN on chromosomes 2D, 7A, and 7B, for SL on 5A, and FT on 2D, with 2D-QTL being the functional marker for the gene Ppd-D1. The physical region of the 7A-QTL for TSN revealed the presence of a wheat ortholog (TaAPO-A1) to APO1–a rice gene that positively controls the spikelet number on the panicles. Interspecific analyses of the TaAPO-A1 orthologs showed that it is a highly conserved gene important for floral development and present in a wide range of terrestrial plants. Intraspecific studies of the TaAPO-A1 across wheat genotypes revealed a polymorphism in the conserved F-box domain, defining two haplotypes. A KASP marker developed on the polymorphic site showed a highly significant association of TaAPO-A1 with TSN, explaining 23.2% of the total genotypic variance. Also, the TaAPO-A1 alleles showed weak but significant differences for SL and grain yield. Our results demonstrate the importance of wheat sequence resources to identify candidate genes for important traits based on genetic analyses.

Gene Expression Profiles and microRNA Regulation Networks in Tiller Primordia, Stem Tips, and Young Spikes of Wheat Guomai 301

Tillering and spike differentiation are two key events for wheat (Triticum aestivum L.). A study on the transcriptomes and microRNA group profiles of wheat at the two key developmental stages will bring insight into the molecular regulation mechanisms. Guomai 301 is a representative excellent new high yield wheat cultivar in the Henan province in China. The transcriptomes and microRNA (miRNA) groups of tiller primordia (TPs), stem tips (STs), and young spikes (YSs) in Guomai 301 were compared to each other. A total of 1741 tillering specifically expressed and 281 early spikes differentiating specifically expressed differentially expressed genes (DEGs) were identified. Six major expression profile clusters of tissue-specific DEGs for the three tissues were classified by gene co-expression analysis using K-means cluster. The ribosome (ko03010), photosynthesis-antenna proteins (ko00196), and plant hormone signal transduction (ko04075) were the main metabolic pathways in TPs, STs, and YSs, respectively. Similarly, 67 TP specifically expressed and 19 YS specifically expressed differentially expressed miRNAs were identified, 65 of them were novel. The roles of 3 well known miRNAs, tae-miR156, tae-miR164, and tae-miR167a, in post-transcriptional regulation were similar to that of other researches. There were 651 significant negative miRNA-mRNA interaction pairs in TPs and YSs, involving 63 differentially expressed miRNAs (fold change > 4) and 416 differentially expressed mRNAs. Among them 12 key known miRNAs and 16 novel miRNAs were further analyzed, and miRNA-mRNA regulatory networks during tillering and early spike differentiating were established.