Rapid identification of a stripe rust resistant gene in a space-induced wheat mutant using specific locus amplified fragment (SLAF) sequencing

TavWA1 is critical for wheat growth by modulating cell morphology and arrangement

Plant growth is determined by the production of cells and initiation of new organs. Exploring genes that control cell number and cell size is of great significance for understanding plant growth regulation. In this study, we characterized two wheat mutants, ah and dl, with abnormal growth. The ah mutant is a naturally occurring variant characterized by severe dwarfism, increased tiller number, and reduced grain length, while the dl mutant is derived from an ethyl methane sulfonate (EMS)-mutagenized population and exhibits smaller grain size and slightly reduced plant height. Cytological analyses revealed abnormal cell number, cell morphology and arrangement in the stems and leaves of the ah mutant, along with reduced cell length in the grains of the dl mutant. Map-based cloning identified that both mutants carry mutations in the same gene TavWA1-7D, which encodes a protein with a von Willebrand factor A (vWA) domain. The ah mutant harbors a 174-bp insertion in the 1,402-bp coding sequence (CDS) of TavWA1-7D, causing premature termination of protein translation, while the dl mutant contains a Glu420Lys substitution. Mimicking the TavWA1-7Dah through clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease 9-mediated genome editing leads to a severe dwarfism phenotype. The C-terminus of the protein is crucial for its correct subcellular localization and interaction, supporting its critical role for TavWA1-7D function. Proteomic analysis showed that the dwarf phenotype of the ah mutant is associated with impaired photosynthesis, ribosome function, and nucleosome formation. Additionally, TavWA1-7D interacts with an E3 ligase, TaVIP1-3B, the expression levels of which are elevated in both mutants. Overexpression and knockout studies of TaVIP1-3B demonstrated its negative regulatory role in cell length and grain size. Together, our findings suggest that TavWA1-7D plays a vital role in regulating wheat growth and yield-related traits, with the dl mutant's short grain phenotype being associated with TaVIP1-3B expression levels.

Integrated genome-wide association and transcriptomic analysis to identify receptor kinase genes to stripe rust resistance in wheat germplasm from southwestern China

Stripe rust of wheat, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat worldwide. Identification of new and elite Pst-resistance loci or genes has the potential to enhance overall resistance to this pathogen. Here, we conducted an integrated genome-wide association study (GWAS) and transcriptomic analysis to screen for loci associated with resistance to stripe rust in 335 accessions from Yunnan, including 311 landraces and 24 cultivars. Based on the environmental phenotype, we identified 113 protein kinases significantly associated with Pst resistance using mixed linear model (MLM) and generalized linear model (GLM) models. Transcriptomic analysis revealed that 52 of 113 protein kinases identified by GWAS were up and down regulated in response to Pst infection. Among these genes, a total of 15 receptor kinase genes were identified associated with Pst resistance. 11 candidate genes were newly discovered in Yunnan wheat germplasm. Our results revealed that resistance alleles to stripe rust were accumulated in Yunnan wheat germplasm, implying direct or indirect selection for improving stripe rust resistance in elite wheat breeding programs.

A diverse panel of 755 bread wheat accessions harbors untapped genetic diversity in landraces and reveals novel genetic regions conferring powdery mildew resistance

KEY MESSAGE

A bread wheat panel reveals rich genetic diversity in Turkish, Pakistani and Iranian landraces and novel resistance loci to diverse powdery mildew isolates via subsetting approaches in association studies. Wheat breeding for disease resistance relies on the availability and use of diverse genetic resources. More than 800,000 wheat accessions are globally conserved in gene banks, but they are mostly uncharacterized for the presence of resistance genes and their potential for agriculture. Based on the selective reduction of previously assembled collections for allele mining for disease resistance, we assembled a trait-customized panel of 755 geographically diverse bread wheat accessions with a focus on landraces, called the LandracePLUS panel. Population structure analysis of this panel based on the TaBW35K SNP array revealed an increased genetic diversity compared to 632 landraces genotyped in an earlier study and 17 high-quality sequenced wheat accessions. The additional genetic diversity found here mostly originated from Turkish, Iranian and Pakistani landraces. We characterized the LandracePLUS panel for resistance to ten diverse isolates of the fungal pathogen powdery mildew. Performing genome-wide association studies and dividing the panel further by a targeted subsetting approach for accessions of distinct geographical origin, we detected several known and already cloned genes, including the Pm2a gene. In addition, we identified 22 putatively novel powdery mildew resistance loci that represent useful sources for resistance breeding and for research on the mildew-wheat pathosystem. Our study shows the value of assembling trait-customized collections and utilizing a diverse range of pathogen races to detect novel loci. It further highlights the importance of integrating landraces of different geographical origins into future diversity studies.

Mapping and identification of a reverse mutation of Rht2 that enhances plant height and thousand grain weight in an elite wheat mutant induced by spaceflight

As climate change continues to negatively impact our farmlands, abiotic factors like salinity and drought stress increasingly threaten global food security. The development of elite germplasms with resistance to multiple abiotic stresses is essential for breeding climate-resilient wheat cultivars. In this study, we determined that the previously reported salt-tolerant st1 mutant, obtained via spaceflight mutagenesis, may also resist to drought stress at the seedling stage. Moreover, our field trial revealed that yield-related traits including plant height, 1000-grain weight, and spike number per plant were significantly increased in st1 compared to the wild type. An F2 population of 334 individuals derived from a cross between the wild type and st1 displayed a bimodal distribution indicating that st1 plant height is controlled by a single major gene. Our Bulked Segregant Analysis and exome capture sequencing indicate that this gene is located on chromosome 4D. Further genetic linkage and gene sequence analysis suggests that a reverse mutation of Rht2 is putatively responsible for plant height variation in st1. Our genotypic and phenotypic analysis of the F2 population and F3 lines indicate that this reverse mutation significantly increases plant height and thousand grain weight but slightly decreases spike number per plant. Together, these results supply helpful information for the utilization of Rht2 in wheat breeding and provide an important material for breeding environmentally resilient, high-yield wheat varieties.

Phenotyping and Exploitation of Kompetitive Allele-Specific PCR Assays for Genes Underpinning Leaf Rust Resistance in New Spring Wheat Mutant Lines

Leaf rust (Puccinia triticina Eriks) is a wheat disease causing substantial yield losses in wheat production globally. The identification of genetic resources with permanently effective resistance genes and the generation of mutant lines showing increased levels of resistance allow the efficient incorporation of these target genes into germplasm pools by marker-assisted breeding. In this study, new mutant (M3 generation) lines generated from the rust-resistant variety Kazakhstanskaya-19 were developed using gamma-induced mutagenesis through 300-, 350-, and 400-Gy doses. In field trials after leaf rust inoculation, 75 mutant lines showed adult plant resistance. These lines were evaluated for resistance at the seedling stage via microscopy in greenhouse experiments. Most of these lines (89.33%) were characterized as resistant at both developmental stages. Hyperspectral imaging analysis indicated that infected leaves of wheat genotypes showed increased relative reflectance in visible and near-infrared light compared to the non-infected genotypes, with peak means at 462 and 644 nm, and 1936 and 2392 nm, respectively. Five spectral indexes, including red edge normalized difference vegetation index (RNDVI), structure-insensitive pigment index (SIPI), ratio vegetation index (RVSI), water index (WI), and normalized difference water index (NDWI), demonstrated significant potential for determining disease severity at the seedling stage. The most significant differences in reflectance between susceptible and resistant mutant lines appeared at 694.57 and 987.51 nm. The mutant lines developed were also used for the development and validation of KASP markers for leaf rust resistance genes Lr1, Lr2a, Lr3, Lr9, Lr10, and Lr17. The mutant lines had high frequencies of "a" resistance alleles (0.88) in all six Lr genes, which were significantly associated with seedling resistance and suggest the potential of favorable haplotype introgression through functional markers. Nine mutant lines characterized by the presence of "b" alleles in Lr9 and Lr10-except for one line with allele "a" in Lr9 and three mutant lines with allele "a" in Lr10-showed the progressive development of fungal haustorial mother cells 72 h after inoculation. One line from 300-Gy-dosed mutant germplasm with "b" alleles in Lr1, Lr2a, Lr10, and Lr17 and "a" alleles in Lr3 and Lr9 was characterized as resistant based on the low number of haustorial mother cells, suggesting the contribution of the "a" alleles of Lr3 and Lr9.

Genome-Wide Identification, Characterization and Expression Analysis of the TaDUF724 Gene Family in Wheat (Triticum aestivum)

Unknown functional domain (DUF) proteins constitute a large number of functionally uncharacterized protein families in eukaryotes. DUF724s play crucial roles in plants. However, the insight understanding of wheat TaDUF724s is currently lacking. To explore the possible function of TaDUF724s in wheat growth and development and stress response, the family members were systematically identified and characterized. In total, 14 TaDUF724s were detected from a wheat reference genome; they are unevenly distributed across the 11 chromosomes, and, according to chromosome location, they were named TaDUF724-1 to TaDUF724-14. Evolution analysis revealed that TaDUF724s were under negative selection, and fragment replication was the main reason for family expansion. All TaDUF724s are unstable proteins; most TaDUF724s are acidic and hydrophilic. They were predicted to be located in the nucleus and chloroplast. The promoter regions of TaDUF724s were enriched with the cis-elements functionally associated with growth and development, as well as being hormone-responsive. Expression profiling showed that TaDUF724-9 was highly expressed in seedings, roots, leaves, stems, spikes and grains, and strongly expressed throughout the whole growth period. The 12 TaDUF724 were post-transcription regulated by 12 wheat MicroRNA (miRNA) through cleavage and translation. RT-qPCR showed that six TaDUF724s were regulated by biological and abiotic stresses. Conclusively, TaDUF724s were systematically analyzed using bioinformatics methods, which laid a theoretical foundation for clarifying the function of TaDUF724s in wheat.

Single-step genome-wide association study for susceptibility to Teratosphaeria nubilosa and precocity of vegetative phase change in Eucalyptus globulus

Introduction

Mycosphaerella leaf disease (MLD) is one of the most prevalent foliar diseases of Eucalyptus globulus plantations around the world. Since resistance management strategies have not been effective in commercial plantations, breeding to develop more resistant genotypes is the most promising strategy. Available genomic information can be used to detect genomic regions associated with resistance to MLD, which could significantly speed up the process of genetic improvement.

Methods

We investigated the genetic basis of MLD resistance in a breeding population of E. globulus which was genotyped with the EUChip60K SNP array. Resistance to MLD was evaluated through resistance of the juvenile foliage, as defoliation and leaf spot severity, and through precocity of change to resistant adult foliage. Genome-wide association studies (GWAS) were carried out applying four Single-SNP models, a Genomic Best Linear Unbiased Prediction (GBLUP-GWAS) approach, and a Single-step genome-wide association study (ssGWAS).

Results

The Single-SNP (model K) and GBLUP-GWAS models detected 13 and 16 SNP-trait associations in chromosomes 2, 3 y 11; whereas the ssGWAS detected 66 SNP-trait associations in the same chromosomes, and additional significant SNP-trait associations in chromosomes 5 to 9 for the precocity of phase change (proportion of adult foliage). For this trait, the two main regions in chromosomes 3 and 11 were identified for the three approaches. The SNPs identified in these regions were positioned near the key miRNA genes, miR156.5 and miR157.4, which have a main role in the regulation of the timing of vegetative change, and also in the response to environmental stresses in plants.

Discussion

Our results demonstrated that ssGWAS was more powerful in detecting regions that affect resistance than conventional GWAS approaches. Additionally, the results suggest a polygenic genetic architecture for the heteroblastic transition in E. globulus and identified useful SNP markers for the development of marker-assisted selection strategies for resistance to MLD.

Genome-wide identification and functional characterization of wheat Brassinazole-resistant transcription factors in response to abiotic stresses and stripe rust infection

Brassinazole-resistant (BZR) transcription factors (TFs) are key players in brassinolides (BRs) signaling pathway, which is widely involved in regulating plant growth and development, as well as in plant responding to a variety stresses. Despite their critical roles, little is known about BZR TFs in wheat. In this study, we performed genome-wide analysis of BZR gene family from wheat genome, and 20 TaBZRs were identified. Based on the phylogenetic relationships of TaBZR and BZRs from rice and Arabidopsis, all BZR genes were clustered into four groups. The intron-exon structural patterns and conserved protein motifs of TaBZRs showed high group specificity. TaBZR5, 7, and 9 were significantly induced after salt, drought treatment, and stripe rust infection. However, TaBZR16, which was significantly upregulated under NaCl application, was not expressed during wheat-stripe rust fungus interaction. These results indicated that BZR genes in wheat play different roles in response to various stresses. The results of this study will lay a foundation for further in-depth functional studies of TaBZRs and will provide information for the breeding and genetic improvement of wheat against drought and salt stresses.

Phenotypic characterization of the wheat temperature-sensitive leaf color mutant and physical mapping of mutant gene by reduced-representation sequencing

Few available leaf color mutants in crops have greatly limited the understanding of photosynthesis mechanisms, leading to few accomplishments in crop yield improvement via enhanced photosynthetic efficiency. Here, a noticeable albino mutant, CN19M06, was identified. A comparison between CN19M06 and the wild type CN19 at different temperatures showed that the albino mutant was temperature-sensitive and produced leaves with a decreased chlorophyll content at temperatures below 10 °C. Genetic analysis suggested that the albinism was controlled by one recessive nuclear gene named TSCA1, which was putatively assigned to the region of 718.1-729.8 Mb on chromosome 2AL using bulked-segregant analysis and double-digest restriction site-associated DNA. Finally, molecular linkage analysis physically anchored TSCA1 to a narrowed region of 718.8-725.3 Mb with a 6.5 Mb length on 2AL flanked by InDel 18 and InDel 25 with 0.7 cM genetic interval. Among the 111 annotated functional genes in the corresponding chromosomal region, only TraesCS2A01G487900 of the PAP fibrillin family was both related to chlorophyll metabolism and temperature sensitivity; therefore, it was considered the putative candidate gene of TSCA1. Overall, CN19M06 has great potential for exploring the molecular mechanism of photosynthesis and monitoring temperature changes in wheat production.

Quantitative trait locus mapping and identification of candidate genes for resistance to Verticillium wilt in four recombinant inbred line populations of Gossypium hirsutum

Improving resistance to Verticillium wilt is of great significance for achieving high and stable yields of Upland cotton (Gossypium hirsutum). To deeply understand the genetic basis of cotton resistance to Verticillium wilt, Verticillium wilt-resistant Upland Lumianyan 28 and four Verticillium wilt-susceptible Acala cotton cultivars were used to create four recombinant inbred line (RIL) populations of 469 families through nested hybridization. Phenotypic data collected in five stressful environments were used to select resistant and sensitive lines and create a mixed pool of extreme phenotypes for BSA-seq. A total of 8 QTLs associated with Verticillium wilt resistance were identified on 4 chromosomes, of which qVW-A12-5 was detected simultaneously in the RIL populations and in one of the RIL populations and was identified for the first time. According to the sequence comparison and transcriptome analysis of candidate genes in the QTL interval between parents and pools, 4 genes were identified in the qVW-A12-5 interval. qRT-PCR of parental and phenotypically extreme lines revealed that Gh_CPR30 was induced by and may be a candidate gene for resistance to Verticillium wilt in G. hirsutum. Furthermore, VIGS technology revealed that the disease severity index (DSI) of the Gh_CPR30-silenced plants was significantly higher than that of the control. These results indicate that the Gh_CPR30 gene plays an important role in the resistance of G. hirsutum to Verticillium wilt, and the study provides a molecular basis for analyzing the molecular mechanism underlying G. hirsutum resistance to Verticillium wilt.

A pleiotropic QTL increased economic water use efficiency in bread wheat (Triticum aestivum L.)

Wheat is one of the most important food crops in the world and drought can severely impact on wheat productivity. The identification and deployment of genes for improved water use efficiency (WUE) can help alleviate yield loss under water limitation. In this study, a high-density genetic linkage map of wheat recombinant inbred lines (Ningchun 4 x Ningchun 27) containing 8751 specific locus amplified fragment (SLAF) tags (including 14757 SNPs), with a total map distance of 1685 cM and an average inter-marker map distance of 0.19 cM was constructed by SLAF-seq technology. The economic yield WUE and nine related traits under three water treatments was monitored over four years. The results showed that loci conditioning WUE were also associated with grain carbon isotope discrimination (CID), flag leaf chlorophyll content, plant height, 1000-grain weight, grain weight per spike and grain number per spike. One locus on chromosome 2B explained 26.3% WUE variation in multiple environments. Under good soil moisture conditions before flowering, the high CID genotype QWue.acn-2B Ningchun 27, was associated with WUE, high grain weight per spike, and kilo-grain weight. Under rain-fed conditions, the low CID genotype QWue.acn-2B Ningchun 4 tended to maintain more spike number and was associated with improved WUE and yield. The introduction of good chromosome fragments of QWue.acn-2B into elite lines by molecular marker assisted selection will boost up the cultivation of high-yield and water-saving wheat varieties.

Fine mapping of the reduced height gene Rht22 in tetraploid wheat landrace Jianyangailanmai (Triticum turgidum L.)

Rht22 was fine mapped in the interval of 0.53-1.48 Mb on 7AS, which reduces cell number of internode to cause semi-dwarfism in Jianyangailanmai. As a valuable germplasm resource for wheat genetic improvement, tetraploid wheat has several reduced height (Rht) and enhanced harvest index genes. Rht22, discovered in Jianyangailanmai (JAM, Triticum turgidum L., 2n = 4x = 28, AABB), significantly increases the spikelet number per spike, but its accurate chromosomal position is still unknown. In this study, a high-density genetic map was constructed using specific-length amplified fragment sequencing in an F₇ RIL_DJ population, which was derived from a cross between dwarf Polish wheat (T. polonicum L., 2n = 4x = 28, AABB) and JAM. Two plant height loci, Qph.sicau-4B and Qph.sicau-7A, were mapped on chromosomes 4BS and 7AS, respectively. Qph.sicau-7A was mapped to the 0.33-4.46 Mb interval on 7AS and likely represents the candidate region of Rht22. Fine mapping confirmed and narrowed Rht22 on chromosome arm 7AS between Xbag295.s53 and Xb295.191 in three different populations. The physical region ranged from 0.53 to 1.48 Mb and included 18 candidate genes. Transcriptome analysis of two pairs of near-isogenic lines revealed that 135 differentially expressed genes (DEGs) were associated with semi-dwarfism. Of these, the expression of 83 annotated DEGs involved in hormones synthesis and signal transduction, cell wall composition, DNA replication, microtubule and phragmoplast arrays was significantly down-regulated in the semi-dwarf line. Therefore, Rht22 causes semi-dwarfism in JAM by disrupting these cellular processes, which impairs cell proliferation and reduces internode cell number.

Meta-QTLs for multiple disease resistance involving three rusts in common wheat (Triticum aestivum L.)

In wheat, multiple disease resistance meta-QTLs (MDR-MQTLs) and underlying candidate genes for the three rusts were identified which may prove useful for development of resistant cultivars. Rust diseases in wheat are a major threat to global food security. Therefore, development of multiple disease-resistant cultivars (resistant to all three rusts) is a major goal in all wheat breeding programs worldwide. In the present study, meta-QTLs and candidate genes for multiple disease resistance (MDR) involving all three rusts were identified using 152 individual QTL mapping studies for resistance to leaf rust (LR), stem rust (SR), and yellow rust (YR). From these 152 studies, a total of 1,146 QTLs for resistance to three rusts were retrieved, which included 368 QTLs for LR, 291 QTLs for SR, and 487 QTLs for YR. Of these 1,146 QTLs, only 718 QTLs could be projected onto the consensus map saturated with 2, 34,619 markers. Meta-analysis of the projected QTLs resulted in the identification of 86 MQTLs, which included 71 MDR-MQTLs. Ten of these MDR-MQTLs were referred to as the 'Breeders' MQTLs'. Seventy-eight of the 86 MQTLs could also be anchored to the physical map of the wheat genome, and 54 MQTLs were validated by marker-trait associations identified during earlier genome-wide association studies. Twenty MQTLs (including 17 MDR-MQTLs) identified in the present study were co-localized with 44 known R genes. In silico expression analysis allowed identification of several differentially expressed candidate genes (DECGs) encoding proteins carrying different domains including the following: NBS-LRR, WRKY domains, F-box domains, sugar transporters, transferases, etc. The introgression of these MDR loci into high-yielding cultivars should prove useful for developing high yielding cultivars with resistance to all the three rusts.

Expression Profiles and Characteristics of Apple lncRNAs in Roots, Phloem, Leaves, Flowers, and Fruit

LncRNAs impart crucial effects on various biological processes, including biotic stress responses, abiotic stress responses, fertility and development. The apple tree is one of the four major fruit trees in the world. However, lncRNAs's roles in different tissues of apple are unknown. We identified the lncRNAs in five tissues of apples including the roots, phloem, leaves, flowers, and fruit, and predicted the intricate regulatory networks. A total of 9440 lncRNAs were obtained. LncRNA target prediction revealed 10,628 potential lncRNA-messenger RNA (mRNA) pairs, 9410 pairs functioning in a cis-acting fashion, and 1218 acting in a trans-acting fashion. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to photosynthesis-antenna proteins, single-organism metabolic process and glutathione metabolism. Additionally, a total of 88 lncRNAs have various functions related to microRNAs (miRNAs) as miRNA precursors. Interactions between lncRNAs and miRNAs were predicted, 1341 possible interrelations between 187 mdm-miRNAs and 174 lncRNAs (1.84%) were identified. MSTRG.121644.5, MSTRG.121644.8, MSTRG.2929.2, MSTRG.3953.2, MSTRG.63448.2, MSTRG.9870.2, and MSTRG.9870.3 could participate in the functions in roots as competing endogenous RNAs (ceRNAs). MSTRG.11457.2, MSTRG.138614.2, and MSTRG.60895.2 could adopt special functions in the fruit by working with miRNAs. A further analysis showed that different tissues formed special lncRNA-miRNA-mRNA networks. MSTRG.60895.2-mdm-miR393-MD17G1009000 may participate in the anthocyanin metabolism in the fruit. These findings provide a comprehensive view of potential functions for lncRNAs, corresponding target genes, and related lncRNA-miRNA-mRNA networks, which will increase our knowledge of the underlying development mechanism in apple.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Candidate powdery mildew resistance gene in wheat landrace cultivar Hongyoumai discovered using SLAF and BSR-seq

BACKGROUND

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is an important disease affecting wheat production. Planting resistant cultivars is an effective, safe, and economical method to control the disease. Map construction using next-generation sequencing facilitates gene cloning based on genetic maps and high-throughput gene expression studies. In this study, specific-locus amplified fragment sequencing (SLAF) was used to analyze Huixianhong (female parent), Hongyoumai (male parent) and two bulks (50 homozygous resistant and 50 susceptible F_2:3 segregating population derived from Huixianhong × Hongyoumai to determine a candidate gene region for resistance to powdery mildew on the long arm of chromosome 7B in wheat landrace Hongyoumai. Gene expressions of candidate regions were obtained using bulked segregant RNA-seq in 10 homozygous resistant and 10 susceptible progeny inoculated by Bgt.. Candidate genes were obtained using homology-based cloning in two parents.

RESULTS

A 12.95 Mb long candidate region in chromosome 7BL was identified, and five blocks in SLAF matched the scaffold of the existing co-segregation marker Xmp1207. In the candidate region, 39 differentially expressed genes were identified using RNA-seq, including RGA4 (Wheat_Chr_Trans_newGene_16173)-a disease resistance protein whose expression was upregulated in the resistant pool at 16 h post inoculation with Bgt. Quantitative reverse transcription (qRT)-PCR was used to further verify the expression patterns in Wheat_Chr_Trans_newGene_16173 that were significantly different in the two parents Hongyoumai and Huixianhong. Two RGA4 genes were cloned based on the sequence of Wheat_Chr_Trans_newGene_16173, respectively from two parent and there was one amino acid mutation: S to G in Huixianhong on 510 loci.

CONCLUSION

The combination of SLAF and BSR-seq methods identified a candidate region of pmHYM in the chromosome 7BL of wheat landrace cultivar Hongyoumai. Comparative analysis between the scaffold of co-segregating marker Xmp1207 and SLAF-seq showed five matching blocks. qRT-PCR showed that only the resistant gene Wheat_Chr_Trans_newGene_16173 was significantly upregulated in the resistant parent Hongyoumai after inoculation with Bgt, and gene cloning revealed a difference in one amino acid between the two parent genes, indicating it was involved in the resistance response and may be the candidate resistance gene pmHYM.

VaAPRT3 Gene is Associated With Sex Determination in Vitis amurensis

In the past decade, progress has been made in sex determination mechanism in Vitis. However, genes responsible for sexual differentiation and its mechanism in V. amurensis remain unknown. Here, we identify a sex determination candidate gene coding adenine phosphoribosyl transferase 3 (VaAPRT3) in V. amurensis. Cloning and sequencing of the VaAPRT3 gene allowed us to develop a molecular marker able to discriminate female individuals from males or hermaphrodites based on a 22-bp InDel. Gene expression and endogenous cytokinin content analysis revealed that the VaAPRT3 gene is involved in sex determination or, to be precise, in female organ differentiation, through regulating cytokinin metabolism in V. amurensis. This study enlarged the understanding of sex determination mechanism in the genus Vitis, and the sex marker could be used as a helpful tool for sexual identification in breeding programs as well as in investigation and collection of V. amurensis germplasms.

Rapid Identification of a Stripe Rust Resistance Gene YrXK in Chinese Wheat Line Xike01015 Using Specific Locus Amplified Fragment (SLAF) Sequencing

Wheat stripe rust, an airborne fungal disease caused by Puccinia striiformis Westend. f. sp. tritici, is one of the most devastating diseases of wheat. Chinese wheat cultivar Xike01015 displays high levels of all-stage resistance (ASR) to the current predominant P. striiformis f. sp. tritici race CYR33. In this study, a single dominant gene, designated YrXk, was identified in Xike01015 conferring resistance to CYR33 with genetic analysis of F₂ and BC₁ populations from a cross of Mingxian169 (susceptible) and Xike01015. The specific length amplified fragment sequencing (SLAF-seq) strategy was used to construct a linkage map in the F₂ population. Quantitative trait loci (QTL) analysis mapped YrXk to a 12.4-Mb segment on chromosome1 BS, explaining >86.96% of the phenotypic variance. Gene annotation in the QTL region identified three differential expressed candidate genes, TraesCS1B02G168600.1, TraesCS1B02G170200.1, and TraesCS1B02G172400.1. The qRT-PCR results showed that TraesCS1B02G172400.1 and TraesCS1B02G168600.1 are upregulated and that TraesCS1B02G170200.1 is slightly downregulated after inoculation with CYR33 in the seedling stage, which indicates that these genes may function in wheat resistance to stripe rust. The results of this study can be used in wheat breeding for improving resistance to stripe rust.

Mining for New Sources of Resistance to Powdery Mildew in Genetic Resources of Winter Wheat

Genetic pathogen control is an economical and sustainable alternative to the use of chemicals. In order to breed resistant varieties, information about potentially unused genetic resistance mechanisms is of high value. We phenotyped 8,316 genotypes of the winter wheat collection of the German Federal ex situ gene bank for Agricultural and Horticultural Crops, Germany, for resistance to powdery mildew (PM), Blumeria graminis f. sp. tritici, one of the most important biotrophic pathogens in wheat. To achieve this, we used a semi-automatic phenotyping facility to perform high-throughput detached leaf assays. This data set, combined with genotyping-by-sequencing (GBS) marker data, was used to perform a genome-wide association study (GWAS). Alleles of significantly associated markers were compared with SNP profiles of 171 widely grown wheat varieties in Germany to identify currently unexploited resistance conferring genes. We also used the Chinese Spring reference genome annotation and various domain prediction algorithms to perform a domain enrichment analysis and produced a list of candidate genes for further investigation. We identified 51 significantly associated regions. In most of these, the susceptible allele was fixed in the tested commonly grown wheat varieties. Eleven of these were located on chromosomes for which no resistance conferring genes have been previously reported. In addition to enrichment of leucine-rich repeats (LRR), we saw enrichment of several domain types so far not reported as relevant to PM resistance, thus, indicating potentially novel candidate genes for the disease resistance research and prebreeding in wheat.

Elucidating the Molecular Mechanisms by which Seed-Borne Endophytic Fungi, Epichloë gansuensis, Increases the Tolerance of Achnatherum inebrians to NaCl Stress

Seed-borne endophyte Epichloë gansuensis enhance NaCl tolerance in Achnatherum inebrians and increase its biomass. However, the molecular mechanism by which E. gansuensis increases the tolerance of host grasses to NaCl stress is unclear. Hence, we firstly explored the full-length transcriptome information of A. inebrians by PacBio RS II. In this work, we obtained 738,588 full-length non-chimeric reads, 36,105 transcript sequences and 27,202 complete CDSs from A. inebrians. We identified 3558 transcription factors (TFs), 15,945 simple sequence repeats and 963 long non-coding RNAs of A. inebrians. The present results show that 2464 and 1817 genes were differentially expressed by E. gansuensis in the leaves of E+ and E− plants at 0 mM and 200 mM NaCl concentrations, respectively. In addition, NaCl stress significantly regulated 4919 DEGs and 502 DEGs in the leaves of E+ and E− plants, respectively. Transcripts associated with photosynthesis, plant hormone signal transduction, amino acids metabolism, flavonoid biosynthetic process and WRKY TFs were differentially expressed by E. gansuensis; importantly, E. gansuensis up-regulated biology processes (brassinosteroid biosynthesis, oxidation–reduction, cellular calcium ion homeostasis, carotene biosynthesis, positive regulation of proteasomal ubiquitin-dependent protein catabolism and proanthocyanidin biosynthesis) of host grass under NaCl stress, which indicated an increase in the ability of host grasses’ adaptation to NaCl stress. In conclusion, our study demonstrates the molecular mechanism for E. gansuensis to increase the tolerance to salt stress in the host, which provides a theoretical basis for the molecular breed to create salt-tolerant forage with endophytes.

Meta-QTLs and candidate genes for stripe rust resistance in wheat

In bread wheat, meta-QTL analysis was conducted using 353 QTLs that were available from earlier studies. When projected onto a dense consensus map comprising 76,753 markers, only 184 QTLs with the required information, could be utilized leading to identification of 61 MQTLs spread over 18 of the 21 chromosomes (barring 5D, 6D and 7D). The range for mean R2 (PVE %) was 1.9% to 48.1%, and that of CI was 0.02 to 11.47 cM; these CIs also carried 37 Yr genes. Using these MQTLs, 385 candidate genes (CGs) were also identified. Out of these CGs, 241 encoded known R proteins and 120 showed differential expression due to stripe rust infection at the seedling stage; the remaining 24 CGs were common in the sense that they encoded R proteins as well as showed differential expression. The proteins encoded by CGs carried the following widely known domains: NBS-LRR domain, WRKY domains, ankyrin repeat domains, sugar transport domains, etc. Thirteen breeders' MQTLs (PVE > 20%) including four pairs of closely linked MQTLs are recommended for use in wheat molecular breeding, for future studies to understand the molecular mechanism of stripe rust resistance and for gene cloning.

In Silico and Transcription Analysis of Trehalose-6-phosphate Phosphatase Gene Family of Wheat: Trehalose Synthesis Genes Contribute to Salinity, Drought Stress and Leaf Senescence

Trehalose-6-phosphate phosphatase (TPP) genes take part in trehalose metabolism and also in stress tolerance, which has been well documented in many species but poorly understood in wheat. The present research has identified a family of 31 TPP genes in Triticum aestivum L. through homology searches and classified them into five clades by phylogenetic tree analysis, providing evidence of an evolutionary status with Hordeum vulgare, Brachypodium distachyon and Oryza sativa. The exon-intron distribution revealed a discrete evolutionary history and projected possible gene duplication occurrences. Furthermore, different computational approaches were used to analyze the physical and chemical properties, conserved domains and motifs, subcellular and chromosomal localization, and three-dimensional (3-D) protein structures. Cis-regulatory elements (CREs) analysis predicted that TaTPP promoters consist of CREs related to plant growth and development, hormones, and stress. Transcriptional analysis revealed that the transcription levels of TaTPPs were variable in different developmental stages and organs. In addition, qRT-PCR analysis showed that different TaTPPs were induced under salt and drought stresses and during leaf senescence. Therefore, the findings of the present study give fundamental genomic information and possible biological functions of the TaTPP gene family in wheat and will provide the path for a better understanding of TaTPPs involvement in wheat developmental processes, stress tolerance, and leaf senescence.

Revealing the coexistence of differentiation and communication in an endemic hare, Lepus yarkandensis (Mammalia, Leporidae) using specific-length amplified fragment sequencing

BACKGROUND

The Yarkand hare (Lepus yarkandensis Günther, 1875) is endemic to oasis and desert areas around the Tarim Basin in the Xinjiang Uyghur Autonomous Region of northwest China; however, genome-wide information for this species remains limited. Moreover, the genetic variation, genetic structure, and phylogenetic relationships of Yarkand hare from the plateau mountain regions have not been reported. Thus, we used specific-length amplified fragment sequencing (SLAF-seq) technology to evaluate the genetic diversity of 76 Yarkand hares from seven geographic populations in the northern and southwestern parts of the Tarim Basin to investigate single-nucleotide polymorphism (SNP) marker-based population differentiation and evolutionary processes. Selective sweep analysis was conducted to identify genetic differences between populations.

RESULTS

Using SLAF-seq, a total of 1,835,504 SNPs were initially obtained, of which 308,942 high-confidence SNPs were selected for further analysis. Yarkand hares exhibited a relatively high degree of genetic diversity at the SNP level. Based on pairwise FST estimates, the north and southwest groups showed a moderate level of genetic differentiation. Phylogenetic tree and population structure analyses demonstrated evident systematic phylogeographical structure patterns consistent with the geographical distribution of the hares. Hierarchical analysis of molecular variation further indicated that genetic variation was mainly observed within populations. Low to moderate genetic differentiation also occurred among populations despite a common genomic background, likely due to geographical barriers, genetic drift, and differential selection pressure of distinct environments. Nevertheless, the observed lineage-mixing pattern, as indicated by the evolutionary tree, principal component analysis, population structure, and TreeMix analyses, suggests a certain degree of gene flow between the north and southwest groups. This may be related to the migration of hares to high-altitude water sources southwest of the basin during glacial climatic oscillations, as well as river re-diffusion and oasis restoration in the basin following the glacial period. We also identified candidate genes, and their associated gene ontology terms and pathways, related to the adaptation of Yarkand hares to different environmental habitats.

CONCLUSIONS

The identified genome-wide SNPs, genetic diversity, and population structure of Yarkand hares expand our understanding of the genetic background of this endemic species and provide valuable insights into its environmental adaptation, allowing for further exploration of the underlying mechanisms.

Genome-wide identification of ZF-HD gene family in Triticum aestivum: Molecular evolution mechanism and function analysis

ZF-HD family genes play important roles in plant growth and development. Studies about the whole genome analysis of ZF-HD gene family have been reported in some plant species. In this study, the whole genome identification and expression profile of the ZF-HD gene family were analyzed for the first time in wheat. A total of 37 TaZF-HD genes were identified and divided into TaMIF and TaZHD subfamilies according to the conserved domain. The phylogeny tree of the TaZF-HD proteins was further divided into six groups based on the phylogenetic relationship. The 37 TaZF-HDs were distributed on 18 of 21 chromosomes, and almost all the genes had no introns. Gene duplication and Ka/Ks analysis showed that the gene family may have experienced powerful purification selection pressure during wheat evolution. The qRT-PCR analysis showed that TaZF-HD genes had significant expression patterns in different biotic stress and abiotic stress. Through subcellular localization experiments, we found that TaZHD6-3B was located in the nucleus, while TaMIF4-5D was located in the cell membrane and nucleus. Our research contributes to a comprehensive understanding of the TaZF-HD family, provides a new perspective for further research on the biological functions of TaZF-HD genes in wheat.

Genome-wide identification and expression analysis of ADP-ribosylation factors associated with biotic and abiotic stress in wheat (Triticum aestivum L.)

The ARF gene family plays important roles in intracellular transport in eukaryotes and is involved in conferring tolerance to biotic and abiotic stresses in plants. To explore the role of these genes in the development of wheat (Triticum aestivum L.), 74 wheat ARF genes (TaARFs; including 18 alternate transcripts) were identified and clustered into seven sub-groups. Phylogenetic analysis revealed that TaARFA1 sub-group genes were strongly conserved. Numerous cis-elements functionally associated with the stress response and hormones were identified in the TaARFA1 sub-group, implying that these TaARFs are induced in response to abiotic and biotic stresses in wheat. According to available transcriptome data and qRT-PCR analysis, the TaARFA1 genes displayed tissue-specific expression patterns and were regulated by biotic stress (powdery mildew and stripe rust) and abiotic stress (cold, heat, ABA, drought and NaCl). Protein interaction network analysis further indicated that TaARFA1 proteins may interact with protein phosphatase 2C (PP2C), which is a key protein in the ABA signaling pathway. This comprehensive analysis will be useful for further functional characterization of TaARF genes and the development of high-quality wheat varieties.

Genome-Wide Identification and Characterization of Wheat 14-3-3 Genes Unravels the Role of TaGRF6-A in Salt Stress Tolerance by Binding MYB Transcription Factor

14-3-3 proteins are a large multigenic family of general regulatory factors (GRF) ubiquitously found in eukaryotes and play vital roles in the regulation of plant growth, development, and response to stress stimuli. However, so far, no comprehensive investigation has been performed in the hexaploid wheat. In the present study, A total of 17 potential 14-3-3 gene family members were identified from the Chinese Spring whole-genome sequencing database. The phylogenetic comparison with six 14-3-3 families revealed that the majority of wheat 14-3-3 genes might have evolved as an independent branch and grouped into ε and non-ε group using the phylogenetic comparison. Analysis of gene structure and motif indicated that 14-3-3 protein family members have relatively conserved exon/intron arrangement and motif composition. Physical mapping showed that wheat 14-3-3 genes are mainly distributed on chromosomes 2, 3, 4, and 7. Moreover, most 14-3-3 members in wheat exhibited significantly down-regulated expression in response to alkaline stress. VIGS assay and protein-protein interaction analysis further confirmed that TaGRF6-A positively regulated slat stress tolerance by interacting with a MYB transcription factor, TaMYB64. Taken together, our findings provide fundamental information on the involvement of the wheat 14-3-3 family in salt stress and further investigating their molecular mechanism.

Genome-wide analysis of growth-regulating factors (GRFs) in Triticum aestivum

The Growth-Regulating Factor (GRF) family encodes a type of plant-specific transcription factor (TF). GRF members play vital roles in plant development and stress response. Although GRF family genes have been investigated in a variety of plants, they remain largely unstudied in bread wheat (Triticum aestivum L.). The present study was conducted to comprehensively identify and characterize the T. aestivum GRF (TaGRF) gene family members. We identified 30 TaGRF genes, which were divided into four groups based on phylogenetic relationship. TaGRF members within the same subgroup shared similar motif composition and gene structure. Synteny analysis suggested that duplication was the dominant reason for family member expansion. Expression pattern profiling showed that most TaGRF genes were highly expressed in growing tissues, including shoot tip meristems, stigmas and ovaries, suggesting their key roles in wheat growth and development. Further qRT-PCR analysis revealed that all 14 tested TaGRFs were significantly differentially expressed in responding to drought or salt stresses, implying their additional involvement in stress tolerance of wheat. Our research lays a foundation for functional determination of TaGRFs, and will help to promote further scrutiny of their regulatory network in wheat development and stress response.

Rapid Mining of Candidate Genes for Verticillium Wilt Resistance in Cotton Based on BSA-Seq Analysis

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as "cancer" of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F_2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F₂ segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

Role of Genetics, Genomics, and Breeding Approaches to Combat Stripe Rust of Wheat

Puccinia striiformis (Pst) is a devastating biotrophic fungal pathogen that causes wheat stripe rust. It usually loves cool and moist places and can cause 100% crop yield losses in a single field when ideal conditions for disease incidence prevails. Billions of dollars are lost due to fungicide application to reduce stripe rust damage worldwide. Pst is a macrocyclic, heteroecious fungus that requires primary (wheat or grasses) as well as secondary host (Berberis or Mahonia spp.) for completion of life cycle. In this review, we have summarized the knowledge about pathogen life cycle, genes responsible for stripe rust resistance, and susceptibility in wheat. In the end, we discussed the importance of conventional and modern breeding tools for the development of Pst-resistant wheat varieties. According to our findings, genetic engineering and genome editing are less explored tools for the development of Pst-resistant wheat varieties; hence, we highlighted the putative use of advanced genome-modifying tools, i.e., base editing and prime editing, for the development of Pst-resistant wheat.

Application of Genomics Tools in Wheat Breeding to Attain Durable Rust Resistance

Wheat is an important source of dietary protein and calories for the majority of the world's population. It is one of the largest grown cereal in the world occupying over 215 M ha. Wheat production globally is challenged by biotic stresses such as pests and diseases. Of the 50 diseases of wheat that are of economic importance, the three rust diseases are the most ubiquitous causing significant yield losses in the majority of wheat production environments. Under severe epidemics they can lead to food insecurity threats amid the continuous evolution of new races of the pathogens, shifts in population dynamics and their virulence patterns, thereby rendering several effective resistance genes deployed in wheat breeding programs vulnerable. This emphasizes the need to identify, characterize, and deploy effective rust-resistant genes from diverse sources into pre-breeding lines and future wheat varieties. The use of genetic resistance has been marked as eco-friendly and to curb the further evolution of rust pathogens. Deployment of multiple rust resistance genes including major and minor genes in wheat lines could enhance the durability of resistance thereby reducing pathogen evolution. Advances in next-generation sequencing (NGS) platforms and associated bioinformatics tools have revolutionized wheat genomics. The sequence alignment of the wheat genome is the most important landmark which will enable genomics to identify marker-trait associations, candidate genes and enhanced breeding values in genomic selection (GS) studies. High throughput genotyping platforms have demonstrated their role in the estimation of genetic diversity, construction of the high-density genetic maps, dissecting polygenic traits, and better understanding their interactions through GWAS (genome-wide association studies) and QTL mapping, and isolation of R genes. Application of breeder's friendly KASP assays in the wheat breeding program has expedited the identification and pyramiding of rust resistance alleles/genes in elite lines. The present review covers the evolutionary trends of the rust pathogen and contemporary wheat varieties, and how these research strategies galvanized to control the wheat killer genus Puccinia. It will also highlight the outcome and research impact of cost-effective NGS technologies and cloning of rust resistance genes amid the public availability of common and tetraploid wheat reference genomes.

Genome-Wide Mining of Wheat DUF966 Gene Family Provides New Insights Into Salt Stress Responses

Domain of unknown function (DUF) proteins constitute a great deal of families of functionally uncharacterized proteins in eukaryotes. The DUF966 gene family is found in monocotyledons, dicotyledons, mosses, and other species. However, little is known about the functions of DUF966 genes in wheat (Triticum aestivum L.). In this study, we identified and characterized the TaDUF966 gene family members in wheat by in silico analysis. A total of 28 TaDUF966 proteins were identified in wheat. Phylogenetic analysis divided these proteins into two groups (Groups I and II). Proteins in each group showed a highly conserved DUF966 domain and conserved motif distribution, implying their functional conservation. Analysis of gene expression profiling data showed that some TaDUF966 genes were induced by salt stress. We further confirmed the role of TaDUF966-9B in salt stress using virus induced gene silencing (VIGS) assay. Compared with the empty vector control, the TaDUF966-9B knockdown plants exhibited severe leaf curling at 10 days post-inoculation with BSMV under salt stress, suggesting that TaDUF966 genes play a vital role in salt stress tolerance in wheat. Taken together, these results expand our knowledge of the evolution of the DUF966 gene family in wheat and promote the potential application of these genes in wheat genetic improvement.

Fine-Mapping of the Powdery Mildew Resistance Gene mlxbd in the Common Wheat Landrace Xiaobaidong

Powdery mildew, which is caused by Blumeria graminis f. sp. tritici (Bgt), is a disease of wheat worldwide. Xiaobaidong is a Chinese wheat landrace, which still maintains good resistance against powdery mildew. To obtain more genetic markers closely linked to the powdery mildew resistance gene mlxbd and narrow the candidate region for its isolation, new simple sequence repeats and cross intron-spanning markers were designed based on the genome sequence of Triticum aestivum cultivar Chinese Spring chromosome 7BL. The flanking markers 7BLSSR49 and WGGC5746 were found to be tightly linked to mlxbd at genetic distances of 0.4 cM and 0.3 cM, respectively. The resistance locus was mapped to a 63.40 kb and 0.29 Mb region of the Chinese Spring genome and Zavitan genome, respectively. The linked markers of mlxbd could be used as diagnostic markers for mlxbd. The linked molecular markers and delineated genomic region in the sequenced Chinese Spring genome will assist the future map-based cloning of mlxbd.

Genome-Wide Identification, Structure Characterization, and Expression Profiling of Dof Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

DNA binding with one finger (Dof) proteins are plant-specific transcription factors with crucial roles in plant growth and stress response. Even so, little is known about them in wheat. In this study, 108 wheat Dof (TaDof) genes across 21 chromosomes were detected. Although variable in sequence length, molecular weight, and isoelectric point, all TaDof proteins contained conserved zinc-finger structures and were phylogenetically divided into 7 sub-groups. Exon/intron and motif analyses suggested that TaDof structures and conserved motifs were similar within sub-groups but diverse among sub-groups. Many segmental duplications were identified and Ka/Ks and inter-species synthetic analyses indicated that polyploidization was main reason for increased number of TaDofs. Prediction and experimental confirmation revealed that TaDofs functioned as transcription factors in the nucleus. Expression pattern profiling showed that TaDofs specifically affected growth and development, and biotic and abiotic stress responses. Wheat miRNAs and cis-regulator were predicted as essential players in molding TaDofs expression patterns. qRT-PCR analysis revealed that TaDofs were induced by salt and drought stresses. Customized annotation revealed that TaDofs were widely involved in phytohormone response, defense, growth and development, and metabolism. Our study provided a comprehensive understanding to wheat TaDofs.

Quantitative Trait Locus (QTLs) Mapping for Quality Traits of Wheat Based on High Density Genetic Map Combined With Bulked Segregant Analysis RNA-seq (BSR-Seq) Indicates That the Basic 7S Globulin Gene Is Related to Falling Number

Numerous quantitative trait loci (QTLs) have been identified for wheat quality; however, most are confined to low-density genetic maps. In this study, based on specific-locus amplified fragment sequencing (SLAF-seq), a high-density genetic map was constructed with 193 recombinant inbred lines derived from Chuanmai 42 and Chuanmai 39. In total, 30 QTLs with phenotypic variance explained (PVE) up to 47.99% were identified for falling number (FN), grain protein content (GPC), grain hardness (GH), and starch pasting properties across three environments. Five NAM genes closely adjacent to QGPC.cib-4A probably have effects on GPC. QGH.cib-5D was the only one detected for GH with high PVE of 33.31-47.99% across the three environments and was assumed to be related to the nearest pina-D1 and pinb-D1genes. Three QTLs were identified for FN in at least two environments, of which QFN.cib-3D had relatively higher PVE of 16.58-25.74%. The positive effect of QFN.cib-3D for high FN was verified in a double-haploid population derived from Chuanmai 42 × Kechengmai 4. The combination of these QTLs has a considerable effect on increasing FN. The transcript levels of Basic 7S globulin and Basic 7S globulin 2 in QFN.cib-3D were significantly different between low FN and high FN bulks, as observed through bulk segregant RNA-seq (BSR). These QTLs and candidate genes based on the high-density genetic map would be beneficial for further understanding of the genetic mechanism of quality traits and molecular breeding of wheat.

Construction of a High-Density Genetic Map Based on SLAF Markers and QTL Analysis of Leaf Size in Rice

Leaf shape is an important agronomic trait for constructing an ideal plant type in rice, and high-density genetic map is facilitative in improving accuracy and efficiency for quantitative trait loci (QTL) analysis of leaf trait. In this study, a high-density genetic map contained 10,760 specific length amplified fragment sequencing (SLAF) markers was established based on 149 recombinant inbred lines (RILs) derived from the cross between Rekuangeng (RKG) and Taizhong1 (TN1), which exhibited 1,613.59 cM map distance with an average interval of 0.17 cM. A total of 24 QTLs were detected and explained the phenotypic variance ranged from 9% to 33.8% related to the leaf morphology across two areas. Among them, one uncloned major QTL qTLLW1 (qTLL1 and qTLLW1) involved in regulating leaf length and leaf width with max 33.8% and 22.5% phenotypic variance respectively was located on chromosome 1, and another major locus qTLW4 affecting leaf width accounted for max 25.3% phenotypic variance was mapped on chromosome 4. Fine mapping and qRT-PCR expression analysis indicated that qTLW4 may be allelic to NAL1 (Narrow leaf 1) gene.

Genome-wide identification and transcriptional expression analysis of superoxide dismutase (SOD) family in wheat (Triticum aestivum)

Superoxide dismutases (SODs) are a family of key antioxidant enzymes that play a crucial role in plant growth and development. Previously, this gene family has been investigated in Arabidopsis and rice. In the present study, a genome-wide analysis of the SOD gene family in wheat were performed. Twenty-six SOD genes were identified from the whole genome of wheat, including 17 Cu/Zn-SODs, six Fe-SODs, and three Mn-SODs. The chromosomal location mapping analysis indicated that these three types of SOD genes were only distributed on 2, 4, and 7 chromosomes, respectively. Phylogenetic analyses of wheat SODs and several other species revealed that these SOD proteins can be assigned to two major categories. SOD1 mainly comprises of Cu/Zn-SODs, and SOD2 mainly comprises of Fe-SODs and Mn-SODs. Gene structure and motif analyses indicated that most of the SOD genes showed a relatively conserved exon/intron arrangement and motif composition. Analyses of transcriptional data indicated that most of the wheat SOD genes were expressed in almost all of the examined tissues and had important functions in abiotic stress resistance. Finally, quantitative real-time polymerase chain reaction (qRT-PCR) analysis was used to reveal the regulating roles of wheat SOD gene family in response to NaCl, mannitol, and polyethylene glycol stresses. qRT-PCR showed that eight randomly selected genes with relatively high expression levels responded to all three stresses based on released transcriptome data. However, their degree of response and response patterns were different. Interestingly, among these genes, TaSOD1.7, TaSOD1.9, TaSOD2.1, and TaSOD2.3 feature research value owing to their remarkable expression-fold change in leaves or roots under different stresses. Overall, our results provide a basis of further functional research on the SOD gene family in wheat and facilitate their potential use for applications in the genetic improvement on wheat in drought and salt stress environments.

Genome-Wide Characterization and Expression Profiling of Squamosa Promoter Binding Protein-like (SBP) Transcription Factors in Wheat (Triticum aestivum L.)

Transcription factors (TFs) play fundamental roles in the developmental processes of all living organisms. Squamosa Promoter Binding Protein-like (SBP/SBP-Box) is a major family of plant-specific TFs, which plays important roles in multiple processes involving plant growth and development. While some work has been done, there is a lot more that is yet to be discovered in the hexaploid wheat SBP (TaSBP) family. With the completion of whole genome sequencing, genome-wide analysis of SBPs in common hexaploid wheat is now possible. In this study, we used protein–protein Basic Local Alignment Search Tool (BLASTp) to hunt the newly released reference genome sequence of hexaploid wheat (Chinese spring). Seventy-four TaSBP proteins (belonging to 56 genes) were identified and clustered into five groups. Gene structure and motif analysis indicated that most TaSBPs have relatively conserved exon–intron arrangements and motif composition. Analysis of transcriptional data showed that many TaSBP genes responded to some biological and abiotic stresses with different expression patterns. Moreover, three TaSBP genes were generally expressed in the majority of tissues throughout the wheat growth and also responded to many environmental biotic and abiotic stresses. Collectively, the detailed analyses presented here will help in understanding the roles of the TaSBP and also provide a reference for the further study of its biological function in wheat.

Identification of QTLs for Stripe Rust Resistance in a Recombinant Inbred Line Population

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating fungal diseases of wheat worldwide. It is essential to discover more sources of stripe rust resistance genes for wheat breeding programs. Specific locus amplified fragment sequencing (SLAF-seq) is a powerful tool for the construction of high-density genetic maps. In this study, a set of 200 recombinant inbred lines (RILs) derived from a cross between wheat cultivars Chuanmai 42 (CH42) and Chuanmai 55 (CH55) was used to construct a high-density genetic map and to identify quantitative trait loci (QTLs) for stripe rust resistance using SLAF-seq technology. A genetic map of 2828.51 cM, including 21 linkage groups, contained 6732 single nucleotide polymorphism markers (SNP). Resistance QTLs were identified on chromosomes 1B, 2A, and 7B; Qyr.saas-7B was derived from CH42, whereas Qyr.saas-1B and Qyr.saas-2A were from CH55. The physical location of Qyr.saas-1B, which explained 6.24–34.22% of the phenotypic variation, overlapped with the resistance gene Yr29. Qyr.saas-7B accounted for up to 20.64% of the phenotypic variation. Qyr.saas-2A, a minor QTL, was found to be a likely new stripe rust resistance locus. A significant additive effect was observed when all three QTLs were combined. The combined resistance genes could be of value in breeding wheat for stripe rust resistance.

QTL analysis of durable stripe rust resistance in the North American winter wheat cultivar Skiles

This study determined the effects of growth stage and temperature on expression of high-temperature adult-plant resistance to stripe rust, mapped six QTL for durable resistance in winter wheat Skiles using a doubled haploid population, and selected breeding lines with different combinations of the QTL using marker-assisted selection. The winter wheat cultivar Skiles has a high level of high-temperature adult-plant (HTAP) resistance to stripe rust caused by Puccinia striiformis f. sp. tritici (Pst). The Skiles HTAP resistance was highly effective at the adult-plant stage even under low temperatures, but high temperatures induced earlier expression and increased levels of resistance. To map resistance genes, Skiles was crossed with the susceptible cultivar Avocet S and a doubled haploid (DH) population was developed. The DH population was tested in fields at Pullman, WA, in 2016, 2017 and 2018, Mount Vernon, WA, in 2017 and 2018 under natural infection, and an environmentally controlled greenhouse at the adult-plant stage with the currently predominant race PSTv-37. The population was genotyped using the 90 K Illumina iSelect wheat SNP chip and selected SSR markers on specific chromosomes. In total, 2526 polymorphic markers were used for QTL mapping and six QTL were detected. Two of the six QTL had major effects across all environments, with one mapped on chromosome 3BS, explaining up to 28.2% of the phenotypic variation and the other on chromosome 4BL, explaining up to 41.8%. Minor QTL were mapped on chromosomes 1BL, 5AL, 6B and 7DL. Genotyping 140 wheat cultivars from the US Pacific Northwest revealed high polymorphism of markers for five of the QTL, and five highly resistant lines with the five QTL were selected from Skiles-derived breeding lines using the markers. This study demonstrated that multiple QTL with mostly additive effects contributed to the high-level HTAP resistance in Skiles.

Molecular Mapping of a Recessive Powdery Mildew Resistance Gene in Wheat Cultivar Tian Xuan 45 Using Bulked Segregant Analysis with Polymorphic Single Nucleotide Polymorphism Relative Ratio Distribution

Powdery mildew is a destructive foliar disease of wheat worldwide. Wheat cultivar Tian Xuan 45 exhibits resistance to the highly virulent isolate HY5. Genetic analysis of the F₂ and F_2:3 populations of a cultivar Ming Xian 169/Tian Xuan 45 cross revealed that the resistance to HY5 was controlled by a single recessive gene, temporarily designated as PmTx45. A Manhattan plot with the relative frequency distribution of single nucleotide polymorphisms (SNPs) was used to rapidly narrow down the possible chromosomal regions of the associated genes. This microarray-based bulked segregant analysis (BSA) largely improved traditional analytical methods. PmTx45 was located in chromosomal bin 4BL5-0.86-1.00 and was flanked by SNP marker AX-110673642 and intron length polymorphism (ILP) marker ILP-4B01G269900 with genetic distances of 3.0 and 2.6 cM, respectively. Molecular detection in a panel of wheat cultivars using the markers linked to PmTx45 showed that the presence of PmTx45 in commercial wheat cultivars was rare. Resistance spectrum and chromosomal position analyses indicated that PmTx45 may be a novel recessive gene with moderate powdery mildew resistance. This new microarray-based BSA method is feasible and effective and has the potential application for mapping genes in wheat in marker-assisted breeding.

Identification of cucumber circular RNAs responsive to salt stress

BACKGROUND

Circular RNAs (circRNAs) are 3'-5' head-to-tail covalently closed non-coding RNA that have been proved to play essential roles in many cellular and developmental processes. However, no information relate to cucumber circRNAs is available currently, especially under salt stress condition.

RESULTS

In this study, we sequenced circRNAs in cucumber and a total of 2787 were identified, with 1934 in root and 44 in leaf being differentially regulated under salt stress. Characteristics analysis of these circRNAs revealed following features: most of them are exon circRNAs (79.51%) and they prefer to arise from middle exon(s) of parent genes (2035/2516); moreover, most of circularization events (88.3%) use non-canonical-GT/AG splicing signals; last but not least, pairing-driven circularization is not the major way to generate cucumber circRNAs since very few circRNAs (18) contain sufficient flanking complementary sequences. Annotation and enrichment analysis of both parental genes and target mRNAs were launched to uncover the functions of differentially expressed circRNAs induced by salt stress. The results showed that circRNAs may be paly roles in salt stress response by mediating transcription, signal transcription, cell cycle, metabolism adaptation, and ion homeostasis related pathways. Moreover, circRNAs may function to regulate proline metabolisms through regulating associated biosynthesis and degradation genes.

CONCLUSIONS

The present study identified large number of cucumber circRNAs and function annotation revealed their possible biological roles in response to salt stress. Our findings will lay a solid foundation for further structure and function studies of cucumber circRNAs.

Transcriptome Analysis Reveals the Accumulation Mechanism of Anthocyanins in Buckwheat (Fagopyrum esculentum Moench) Cotyledons and Flowers

Buckwheat (Fagopyrum esculentum) is a valuable crop which can produce multiple human beneficial secondary metabolites, for example, the anthocyanins in sprouts and flowers. However, as the predominant group of visible polyphenols in pigmentation, little is known about the molecular mechanisms underlying the anthocyanin biosynthesis within buckwheat. In this study, a comparative transcriptome analysis of green and red common buckwheat cultivars was carried out through RNA sequencing. Overall, 3727 and 5323 differently expressed genes (DEGs) were identified in flowers and cotyledons, respectively. Through GO and KEGG analysis, we revealed that DEGs in flowers and cotyledons are predominately involved in biosynthesis of anthocyanin. A total of 42 unigenes encoding 11 structural enzymes of the anthocyanin biosynthesis were identified as DEGs. We also identified some transcription factor families involved in the regulation of anthocyanin biosynthesis. Real-time qPCR validation of candidate genes was performed in flowers and cotyledons, and the results suggested that the high expression level of structural genes involved in anthocyanin biosynthetic pathway promotes anthocyanin accumulation. Our results provide the insight understanding for coloration of red common buckwheat.

Genome-Wide Assessment of Avocado Germplasm Determined from Specific Length Amplified Fragment Sequencing and Transcriptomes: Population Structure, Genetic Diversity, Identification, and Application of Race-Specific Markers

Genomic data is a powerful tool. However, the phylogenetic relationships among different ecological races of avocado remain unclear. Here, we used the results from specific length amplified fragment sequencing (SLAF-seq) and transcriptome data to infer the population structure and genetic diversity of 21 avocado cultivars and reconstructed the phylogeny of three ecological races and two interracial hybrids. The results of the three analyses performed (unweighted pair-group methods with arithmetic means (UPGMA) cluster, Principal coordinate analysis (PCoA), and STRUCTURE) based on single nucleotide polymorphisms (SNPs) from SLAF-seq all indicated the existence of two populations based on botanical race: Mexican–Guatemalan and West Indian genotype populations. Our results based on SNPs from SLAF-seq indicated that the Mexican and Guatemalan races were more closely related to each other than either was to the West Indian race, which also was confirmed in the UPGMA cluster results based on SNPs from transcriptomic data. SNPs from SLAF-seq provided strong evidence that the Guatemalan, Mexican, and Guatemalan × Mexican hybrid accession possessed higher genetic diversity than the West Indian races and Guatemalan × West Indian hybrid accessions. Six race-specific Kompetitive allele specific PCR (KASP) markers based on SNPs from SLAF-seq were then developed and validated.