Genome-wide impacts of alien chromatin introgression on wheat gene transcriptions

Disomic chromosome 3R(3B) substitution causes a complex of meiotic abnormalities in bread wheat Triticum aestivum L

Triticum aestivum L. lines introgressed with alien chromosomes create a new genetic background that changes the gene expression of both wheat and donor chromosomes. The genes involved in meiosis regulation are localized on wheat chromosome 3B. The purpose of the present study was to investigate the effect of wheat chromosome 3B substituted with homoeologous rye chromosome 3R on meiosis regulation in disomically substituted wheat line 3R(3B). Employing immunostaining with antibodies against microtubule protein, α-tubulin, and the centromere-specific histone (CENH3), as well as FISH, we analyzed microtubule cytoskeleton dynamics and wheat and rye 3R chromosomes behavior in 3R(3B) (Triticum aestivum L. variety Saratovskaya 29 × Secale cereale L. variety Onokhoiskaya) meiosis. The results revealed a set of abnormalities in the microtubule dynamics and chromosome behavior in both first and second divisions. A feature of metaphase I in 3R(3B) was a decrease in the chiasmata number compared with variety Saratovskaya 29, 34.9 ± 0.62 and 41.92 ± 0.38, respectively. Rye homologs 3R in 13.18 % of meiocytes did not form bivalents. Chromosomes were characterized by varying degrees of compaction; 53.33 ± 14.62 cells lacked a metaphase plate. Disturbances were found in microtubule nucleation at the bivalent kinetochores and in their convergence at the spindle division poles. An important feature of meiosis was the asynchronous chromosome behavior in the second division and dyads at the telophase II in 8-13 % of meiocytes, depending on the anther studied. Considering the 3R(3B) meiotic phenotype, chromosome 3B contains the genes involved in the regulation of meiotic division, and substituting 3B3B chromosomes with rye 3R3R does not compensate for their absence.

Wheat powdery mildew resistance gene Pm13 encodes a mixed lineage kinase domain-like protein

Wheat powdery mildew is one of the most destructive diseases threatening global wheat production. The wild relatives of wheat constitute rich sources of diversity for powdery mildew resistance. Here, we report the map-based cloning of the powdery mildew resistance gene Pm13 from the wild wheat species Aegilops longissima. Pm13 encodes a mixed lineage kinase domain-like (MLKL) protein that contains an N-terminal-domain of MLKL (MLKL_NTD) domain in its N-terminus and a C-terminal serine/threonine kinase (STK) domain. The resistance function of Pm13 is validated by mutagenesis, gene silencing, transgenic assay, and allelic association analyses. The development of introgression lines with significantly reduced chromosome segments of Ae. longissima encompassing Pm13 enables widespread deployment of this gene into wheat cultivars. The cloning of Pm13 may provide valuable insights into the molecular mechanisms underlying Pm13-mediated powdery mildew resistance and highlight the important roles of kinase fusion proteins (KFPs) in wheat immunity.

Introgression of the bread wheat D genome encoded Lr34/Yr18/Sr57/Pm38/Ltn1 adult plant resistance gene into Triticum turgidum (durum wheat)

KEY MESSAGE

Lack of function of a D-genome adult plant resistance gene upon introgression into durum wheat. The wheat Lr34/Yr18/Sr57/Pm38/Ltn1 adult plant resistance gene (Lr34), located on chromosome arm 7DS, provides broad spectrum, partial, adult plant resistance to leaf rust, stripe rust, stem rust and powdery mildew. It has been used extensively in hexaploid bread wheat (AABBDD) and conferred durable resistance for many decades. These same diseases also occur on cultivated tetraploid durum wheat and emmer wheat but transfer of D genome sequences to those subspecies is restricted due to very limited intergenomic recombination. Herein we have introgressed the Lr34 gene into chromosome 7A of durum wheat. Durum chromosome substitution line Langdon 7D(7A) was crossed to Cappelli ph1c, a mutant derivative of durum cultivar Cappelli homozygous for a deletion of the chromosome pairing locus Ph1. Screening of BC1F2 plants and their progeny by KASP and PCR markers, 90 K SNP genotyping and cytology identified 7A chromosomes containing small chromosome 7D fragments encoding Lr34. However, in contrast to previous transgenesis experiments in durum wheat, resistance to wheat stripe rust was not observed in either Cappelli/Langdon 7D(7A) or Bansi durum plants carrying this Lr34 encoding segment due to low levels of Lr34 gene expression.

KEY MESSAGE

Global transcriptome analysis of allopolyploidization reveals large-scale repression of the D-subgenome in synthetic hexaploid wheat

Synthetic hexaploid wheat (SHW) lines are created as pre-breeding germplasm to diversify the D subgenome of hexaploid wheat and capitalize upon the untapped genetic diversity of the Aegilops tauschii gene pool. However, the phenotypes observed in the Ae. tauschii parents are not always recovered in the SHW lines, possibly due to inter-subgenome interactions. To elucidate this post-polyploidization genome reprogramming phenomenon, we performed RNA-seq of four SHW lines and their corresponding tetraploid and diploid parents, across ten tissues and three biological replicates. Homoeologue expression bias (HEB) analysis using more than 18,000 triads suggests massive suppression of homoeoalleles of the D subgenome in SHWs. Comparative transcriptome analysis of the whole-genome gene set further corroborated this finding. Alternative splicing analysis of the high-confidence genes indicates an additional layer of complexity where all five splice events are identified, and retained intron is predominant. Homoeologue expression upon resynthesis of hexaploid wheat has implications to the usage and handling of this germplasm in breeding as it relates to capturing the effects of epistatic interaction across subgenomes upon polyploidization. Special considerations must be given to this germplasm in pre-breeding activities to consider the extent of the inter-subgenome interactions on gene expression and their impact on traits for crop improvement.

Whole-genome sequencing uncovers the structural and transcriptomic landscape of hexaploid wheat/Ambylopyrum muticum introgression lines

Wheat is a globally vital crop, but its limited genetic variation creates a challenge for breeders aiming to maintain or accelerate agricultural improvements over time. Introducing novel genes and alleles from wheat's wild relatives into the wheat breeding pool via introgression lines is an important component of overcoming this low variation but is constrained by poor genomic resolution and limited understanding of the genomic impact of introgression breeding programmes. By sequencing 17 hexaploid wheat/Ambylopyrum muticum introgression lines and the parent lines, we have precisely pinpointed the borders of introgressed segments, most of which occur within genes. We report a genome assembly and annotation of Am. muticum that has facilitated the identification of Am. muticum resistance genes commonly introgressed in lines resistant to stripe rust. Our analysis has identified an abundance of structural disruption and homoeologous pairing across the introgression lines, likely caused by the suppressed Ph1 locus. mRNAseq analysis of six of these introgression lines revealed that novel introgressed genes are rarely expressed and those that directly replace a wheat orthologue have a tendency towards downregulation, with no discernible compensation in the expression of homoeologous copies. This study explores the genomic impact of introgression breeding and provides a schematic that can be followed to characterize introgression lines and identify segments and candidate genes underlying the phenotype. This will facilitate more effective utilization of introgression pre-breeding material in wheat breeding programmes.

Genomic sequencing of Thinopyrum elongatum chromosome arm 7EL, carrying fusarium head blight resistance, and characterization of its impact on the transcriptome of the introgressed line CS-7EL

Background

The tall wheatgrass species Thinopyrum elongatum carries a strong fusarium head blight (FHB) resistance locus located on the long arm of chromosome 7 (7EL) as well as resistance to leaf and stem rusts, all diseases with a significant impact on wheat production. Towards understanding the contribution of Th. elongatum 7EL to improvement of disease resistance in wheat, the genomic sequence of the 7EL fragment present in the wheat Chinese Spring (CS) telosomic addition line CS-7EL was determined and the contribution and impact of 7EL on the rachis transcriptome during FHB infection was compared between CS and CS-7EL.

Results

We assembled the Th. elongatum 7EL chromosome arm using a reference-guided approach. Combining this assembly with the available reference sequence for CS hexaploid wheat provided a reliable reference for interrogating the transcriptomic differences in response to infection conferred by the 7EL fragment. Comparison of the transcriptomes of rachis tissues from CS and CS-7EL showed expression of Th. elongatum transcripts as well as modulation of wheat transcript expression profiles in the CS-7EL line. Expression profiles at 4 days after infection with Fusarium graminearum, the causal agent of FHB, showed an increased in expression of genes associated with an effective defense response, in particular glucan endo-1,3-beta-glucosidases and chitinases, in the FHB-resistant line CS-7EL while there was a larger increase in differential expression for genes associated with the level of fungal infection in the FHB-susceptible line CS. One hundred and seven 7EL transcripts were expressed in the smallest 7EL region defined to carry FHB resistance.

Conclusion

7EL contributed to CS-7EL transcriptome by direct expression and through alteration of wheat transcript profiles. FHB resistance in CS-7EL was associated with transcriptome changes suggesting a more effective defense response. A list of candidate genes for the FHB resistance locus on 7EL has been established.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08433-8.

Rye B chromosomes differently influence the expression of A chromosome-encoded genes depending on the host species

The B chromosome (B) is a dispensable component of the genome in many species. To evaluate the impact of Bs on the transcriptome of the standard A chromosomes (A), comparative RNA-seq analyses of rye and wheat anthers with and without additional rye Bs were conducted. In both species, 5-6% of the A-derived transcripts across the entire genomes were differentially expressed in the presence of  2Bs. The GO term enrichment analysis revealed that Bs influence A chromosome encoded processes like "gene silencing"; "DNA methylation or demethylation"; "chromatin silencing"; "negative regulation of gene expression, epigenetic"; "post-embryonic development"; and "chromosome organization." 244 B chromosome responsive A-located genes in + 2B rye and + B wheat shared the same biological function. Positively correlated with the number of Bs, 939 and 1391 B-specific transcripts were identified in + 2B and + 4B wheat samples, respectively. 85% of B-transcripts in + 2B were also found in + 4B transcriptomes. 297 B-specific transcripts were identified in + 2B rye, and 27% were common to the B-derived transcripts identified in + B wheat. Bs encode mobile elements and housekeeping genes, but most B-transcripts were without detectable similarity to known genes. Some of these genes are involved in cell division-related functions like Nuf2 and might indicate their importance in maintaining Bs. The transcriptome analysis provides new insights into the complex interrelationship between standard A chromosomes and supernumerary B chromosomes.

The decreased expression of GW2 homologous genes contributed to the increased grain width and thousand‑grain weight in wheat-Dasypyrum villosum 6VS·6DL translocation lines

This study demonstrated that the aberrant transcription of DvGW2 contributed to the increased grain width and thousand-grain weight in wheat-Dasypyrum villosum T6VS·6DL translocation lines. Due to the high immunity to powdery mildew, Dasypyrum villosum 6VS has been one of the most successful applications of the wild relatives in modern wheat breeding. Along with the desired traits, side-effects could be brought when large alien chromosome fragments are introduced into wheat, but little is known about effects of 6VS on agronomic traits. Here, we found that T6VS·6DL translocation had significantly positive effects on grain weight, plant heightand spike length, and small negative effects on total spikelet number and spikelet compactness using recipient and wheat-D. villosum T6VS·6DL allohexaploid wheats, Wan7107 and Pm97033. Further analysis showed that the 6VS segment might exert direct genetic effect on grain width, then driving the increase of thousand-grain weight. Furthermore, comparative transcriptome analysis identified 2549 and 1282 differentially expressed genes (DEGs) and 2220 and 1496 specifically expressed genes (SEGs) at 6 days after pollination (DAP) grains and 15 DAP endosperms, respectively. Enrichment analysis indicated that the process of cell proliferation category was over-represented in the DEGs. Notably, two homologous genes, TaGW2-D1 and DvGW2, were identified as putative candidate genes associated with grain weight and yield. The expression analysis showed that DvGW2 had an aberrant expression in Pm97033, resulting in significantly lower total expression level of GW2 than Wan7107, which drives the increase of grain weight and width in Pm97033. Collectively, our data indicated that the compromised expression of DvGW2 is critical for increased grain width and weight in T6VS·6DL translocation lines.

Development of Novel Wheat-Aegilops longissima 3S1 Translocations Conferring Powdery Mildew Resistance and Specific Molecular Markers for Chromosome 3S1

Powdery mildew of wheat, caused by Blumeria graminis f. sp. tritici, is a destructive disease of common wheat. Cultivation of resistant varieties is the most cost-effective disease management strategy. Previous studies reported that chromosome 3S^l#2 present in Chinese Spring (CS)-Aegilops longissima 3S^l#2(3B) disomic substitution line TA3575 conferred resistance to powdery mildew. In this study, we further located the powdery mildew resistance gene(s) to the short arm of chromosome 3S^l#2 (3S^l#2S) by evaluating for B. graminis f. sp. tritici resistance of newly developed CS-Ae. longissima 3S^l#2 translocation lines. Meanwhile, TA7545, a previously designated CS-Ae. longissima 3S^l#3 disomic addition line, was reidentified as an isochromosome 3S^l#3S addition line and evaluated to confer resistance to powdery mildew, thus locating the resistance gene(s) to the short arm of chromosome 3S^l#3 (3S^l#3S). Based on transcriptome sequences of TA3575, 10 novel chromosome 3S^lS-specific markers were developed, of which 5 could be used to distinguish between 3S^l#2S and 3S^l#3S derived from Ae. longissima accessions TL20 and TA1910 (TAM4) and the remaining 5 could identify both 3S^l#2S and 3S^l#3S. Also, CL897, one of five markers specific to both 3S^l#2S and 3S^l#3S, could be used to detect Pm13 located at chromosome 3S^l#1S from Ae. longissima accession TL01 in diverse wheat genetic backgrounds. The powdery mildew resistance genes on chromosomes 3S^l#2S and 3S^l#3S, the CS-Ae. longissima 3S^l#2 translocation lines, and the 3S^lS-specific markers developed in this study will facilitate the transfer of B. graminis f. sp. tritici resistance genes into common wheat and provide new germplasm resources for powdery mildew resistance breeding.

Quality evaluation of near-isogenic line of the wheat variety HD2733 carrying the Lr24/Sr24 genomic region

A near-isogenic line (NIL) of the Indian wheat variety HD2733, carrying an introgressed Lr24/Sr24 genomic region was used for studying the effect of this introgression on quality traits. Data on the grain yield and 21 quality traits were recorded in this NIL and its recurrent parent (RP), both of which were grown in a randomized block design for two consecutive years. The statistical analysis revealed that grain yield was on par between the NIL and the RP. The NIL and its RP were both hard grained but the NIL showed a grain hardness index reduced by 9.7%. However, quality traits such as grain weight, protein content, sedimentation value, gluten traits, and solvent retention capacity were significantly higher in the NIL. The NIL also showed an increase in dough stability, a lower degree of softening and a higher farinograph quality number. These results indicated that the NIL could be utilized for hard grain, high protein and strong gluten-based products. An overall improvement in the quality of the NIL over its recurrent parent and without any yield penalty suggests that the Lr24/Sr24 genomic region could be gainfully utilized in wheat breeding for improving the industrial quality of wheat without jeopardising grain yield. The authors suggest that the improved quality of the NIL may be due to the genomic segment carried along with the Lr24/Sr24 genes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02679-x.

3-D Nucleus Architecture in Oat × Maize Addition Lines

The nucleus architecture of hybrid crop plants is not a well-researched topic, yet it can have important implications for their genetic stability and usefulness in the successful expression of agronomically desired traits. In this work we studied the spatial distribution of introgressed maize chromatin in oat × maize addition lines with the number of added maize chromosomes varying from one to four. The number of chromosome additions was confirmed by genomic in situ hybridization (GISH). Maize chromosome-specific simple sequence repeat (SSR) markers were used to identify the added chromosomes. GISH on 3-D root and leaf nuclei was performed to assess the number, volume, and position of the maize-chromatin occupied regions. We revealed that the maize chromosome territory (CT) associations of varying degree prevailed in the double disomic lines, while CT separation was the most common distribution pattern in the double monosomic line. In all analyzed lines, the regions occupied by maize CTs were located preferentially at the nuclear periphery. A comparison between the tissues showed that the maize CTs in the leaf nuclei are positioned closer to the center of the nucleus than in the root nuclei. These findings shed more light on the processes that shape the nucleus architecture in hybrids.