Mapping a leaf senescence gene els1 by BSR-Seq in common wheat

Identification of the early leaf senescence gene ELS3 in bread wheat (Triticum aestivum L.)

MAIN CONCLUSION

Characterization of the early leaf senescence mutant els3 and identification of its causal gene ELS3, which encodes an LRR-RLK protein in wheat. Leaf senescence is an important agronomic trait that affects both crop yield and quality. However, few senescence-related genes in wheat have been cloned and functionally analyzed. Here, we report the characterization of the early leaf senescence mutant els3 and fine mapping of its causal gene ELS3 in wheat. Compared with wild-type Yanzhan4110 (YZ4110), the els3 mutant had a decreased chlorophyll content and a degraded chloroplast structure after the flowering stage. Further biochemical assays in flag leaves showed that the superoxide anion and hydrogen peroxide contents increased, while the activities of antioxidant enzymes, including catalase, superoxide dismutase and glutathione reductase, decreased gradually after the flowering stage in the els3 mutant. To clone the causal gene underlying the phenotype of leaf senescence, a genetic map was constructed using 10,133 individuals of F2:3 populations, and ELS3 was located in a 2.52 Mb region on chromosome 2DL containing 16 putative genes. Subsequent sequence analysis and gene annotation identified only one SNP (C to T) in the first exon of TraesCS2D02G332700, resulting in an amino acid substitution (Pro329Ser), and TraesCS2D02G332700 was preliminarily considered as the candidate gene of ELS3. ELS3 encodes a leucine-rich repeat receptor-like kinase (LRR-RLK) protein that is localized on the cell membrane. We also found that the transient expression of mutant TraesCS2D02G332700 can induce leaf senescence in N. benthamiana. Taken together, TraesCS2D02G332700 is likely to be the candidate gene of ELS3 and may have a function in regulating leaf senescence.

TaACTIN7-D regulates plant height and grain shape in bread wheat

Exploitation of new gene resources and genetic networks contributing to the control of crop yield-related traits, such as plant height, grain size, and shape, may enable us to breed modern high-yielding wheat varieties through molecular methods. In this study, via ethylmethanesulfonate mutagenesis, we identify a wheat mutant plant, mu-597, that shows semi-dwarf plant architecture and round grain shape. Through bulked segregant RNA-seq and map-based cloning, the causal gene for the semi-dwarf phenotype of mu-597 is located. We find that a single-base mutation in the coding region of TaACTIN7-D (TaACT7-D), leading to a Gly-to-Ser (G65S) amino acid mutation at the 65th residue of the deduced TaACT7-D protein, can explain the semi-dwarfism and round grain shape of mu-597. Further evidence shows that the G65S mutation in TaACT7-D hinders the polymerization of actin from monomeric (G-actin) to filamentous (F-actin) status while attenuates wheat responses to multiple phytohormones, including brassinosteroids, auxin, and gibberellin. Together, these findings not only define a new semi-dwarfing gene resource that can be potentially used to design plant height and grain shape of bread wheat but also establish a direct link between actin structure modulation and phytohormone signal transduction.

Identification and validation of quantitative trait loci for chlorophyll content of flag leaf in wheat under different phosphorus treatments

In wheat, the leaf chlorophyll content in flag leaves is closely related to the degree of phosphorus stress. Identifying major genes/loci associated with chlorophyll content in flag leaves under different phosphorus conditions is critical for breeding wheat varieties resistant to low phosphorus (P). Under normal, medium, and low phosphorus conditions, the chlorophyll content of flag leaves was investigated by a double haploid (DH) population derived from a cross between two popular wheat varieties Jinmai 47 and Jinmai 84, at different grain filling stages. Chlorophyll content of the DH population and parents decreased gradually during the S1 to the S3 stages and rapidly at the S4 stage. At the S4 stage, the chlorophyll content of the DH population under low phosphorus conditions was significantly lower than under normal phosphate conditions. Using a wheat 15K single-nucleotide polymorphism (SNP) panel, a total of 157 QTLs were found to be associated with chlorophyll content in flag leaf and were identified under three phosphorus conditions. The phenotypic variation explained (PVE) ranged from 3.07 to 31.66%. Under three different phosphorus conditions, 36, 30, and 48 QTLs for chlorophyll content were identified, respectively. Six major QTLs Qchl.saw-2B.1, Qchl.saw-3B.1, Qchl.saw-4D.1, Qchl.saw-4D.2, Qchl.saw-5A.9 and Qchl.saw-6A.4 could be detected under multiple phosphorus conditions in which Qchl.saw-4D.1, Qchl.saw-4D.2, and Qchl.saw-6A.4 were revealed to be novel major QTLs. Moreover, the closely linked SNP markers of Qchl.saw-4D.1 and Qchl.saw-4D.2 were validated as KASP markers in a DH population sharing the common parent Jinmai 84, showed extreme significance (P <0.01) in more than three environments under different phosphorus conditions, which has the potential to be utilized in molecular marker-assisted breeding for low phosphorus tolerance in wheat.

Capturing Wheat Phenotypes at the Genome Level

Recent technological advances in next-generation sequencing (NGS) technologies have dramatically reduced the cost of DNA sequencing, allowing species with large and complex genomes to be sequenced. Although bread wheat (Triticum aestivum L.) is one of the world's most important food crops, efficient exploitation of molecular marker-assisted breeding approaches has lagged behind that achieved in other crop species, due to its large polyploid genome. However, an international public-private effort spanning 9 years reported over 65% draft genome of bread wheat in 2014, and finally, after more than a decade culminated in the release of a gold-standard, fully annotated reference wheat-genome assembly in 2018. Shortly thereafter, in 2020, the genome of assemblies of additional 15 global wheat accessions was released. As a result, wheat has now entered into the pan-genomic era, where basic resources can be efficiently exploited. Wheat genotyping with a few hundred markers has been replaced by genotyping arrays, capable of characterizing hundreds of wheat lines, using thousands of markers, providing fast, relatively inexpensive, and reliable data for exploitation in wheat breeding. These advances have opened up new opportunities for marker-assisted selection (MAS) and genomic selection (GS) in wheat. Herein, we review the advances and perspectives in wheat genetics and genomics, with a focus on key traits, including grain yield, yield-related traits, end-use quality, and resistance to biotic and abiotic stresses. We also focus on reported candidate genes cloned and linked to traits of interest. Furthermore, we report on the improvement in the aforementioned quantitative traits, through the use of (i) clustered regularly interspaced short-palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated gene-editing and (ii) positional cloning methods, and of genomic selection. Finally, we examine the utilization of genomics for the next-generation wheat breeding, providing a practical example of using in silico bioinformatics tools that are based on the wheat reference-genome sequence.

Identification of Genetic Loci Affecting Flag Leaf Chlorophyll in Wheat Grown under Different Water Regimes

Chlorophyll content of the flag leaf is an important trait for drought resistance in wheat under drought stress. Understanding the regulatory mechanism of flag leaf chlorophyll content could accelerate breeding for drought resistance. In this study, we constructed a recombinant inbred line (RIL) population from a cross of drought-sensitive variety DH118 and drought-resistant variety Jinmai 919, and analyzed the chlorophyll contents of flag leaves in six experimental locations/years using the Wheat90K single-nucleotide polymorphism array. A total of 29 quantitative trait loci (QTLs) controlling flag leaf chlorophyll were detected with contributions to phenotypic variation ranging from 4.67 to 23.25%. Twelve QTLs were detected under irrigated conditions and 18 were detected under dryland (drought) conditions. Most of the QTLs detected under the different water regimes were different. Four major QTLs (Qchl.saw-3B.2, Qchl.saw-5A.2, Qchl.saw-5A.3, and Qchl.saw-5B.2) were detected in the RIL population. Qchl.saw-3B.2, possibly more suitable for marker-assisted selection of genotypes adapted to irrigated conditions, was validated by a tightly linked kompetitive allele specific PCR (KASP) marker in a doubled haploid population derived from a different cross. Qchl.saw-5A.3, a novel stably expressed QTL, was detected in the dryland environments and explained up to 23.25% of the phenotypic variation, and has potential for marker-assisted breeding of genotypes adapted to dryland conditions. The stable and major QTLs identified here add valuable information for understanding the genetic mechanism underlying chlorophyll content and provide a basis for molecular marker–assisted breeding.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

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.

Ne2, a typical CC-NBS-LRR-type gene, is responsible for hybrid necrosis in wheat

Hybrid necrosis, caused by complementary genes Ne1 and Ne2, is a serious barrier for combining desirable traits from different genotypes of wheat, affecting the full utilisation of heterosis. To date, both Ne1 and Ne2 are still not isolated although they were documented decades ago. We report here the map-based cloning and functional characterisation of Ne2, encoding a coiled coil-nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR) protein. Homozygous frameshift mutations generated using the CRISPR/Cas9 approach confirmed the Ne2-inducing hybrid necrosis in wheat. Upregulated expression of Ne2 induced by Ne1 and excess hydrogen peroxide accumulation are associated with the necrosis formation. Genetic analyses of a Ne2 allele (Ne2^m ) and leaf rust resistance gene LrLC10/Lr13 revealed that they might be the same gene. Furthermore, we demonstrated that the frequency of the Ne2 allele was much lower in landraces (2.00%) compared with that in modern cultivars (13.62%), suggesting that Ne2 allele has been partially applied in wheat genetic improvement. Our findings open opportunities of thoroughly investigating the molecular mechanism of hybrid necrosis, selecting Lr13 and simultaneously avoiding hybrid necrosis in wheat breeding through marker-assisted selection.

Fine mapping of Ne1, the hybrid necrosis gene complementary to Ne2 in common wheat (Triticum aestivum L.)

Apart from confinement of Ne1 to a 4.45 Mb genomic segment, markers closely linked to Ne2 were identified and incomplete dominance of both genes in conditioning necrosis severity was shown. Hybrid necrosis in plants is characterized by premature death of leaves or plants in F₁ hybrids. Interaction of two complementary dominant genes Ne1 and Ne2 in wheat (Triticum aestivum L.) is known to cause hybrid necrosis. However, the mechanism underlying this necrosis is still elusive. To obtain markers closely-linked to these two genes, Ne1-carrying cultivar Zheng891 was crossed with Ne2-carrying cultivar Pan555. Using BC₁F₁ plants derived from crosses of the F₁ plants with the two parental lines, Ne1 and Ne2 were mapped to a 2.2 cM interval and a 2.3 cM interval with newly developed markers, respectively. Ne1 was further delimited to a 0.19 cM interval using 2015 Ne2-carrying F₂ plants. Xwgrc3146, Xwgrc3147 and Xwgrc3150, three of the four markers co-segregating with Ne1, were all Zheng891-dominant, suggesting that, compared with Pan555, Ne1 is located in a region with substantial sequence diversity. The Ne1 interval is syntenic to chromosomes 5H, 4, 9 and 2 of barley, Brachypodium distachyon, rice and sorghum, respectively, and corresponds to a 4.45 Mb Chinese Spring sequence. Variations in necrosis severity of the F₂ plants differing in Ne1 and Ne2 genotypes implied that these two genes are incompletely dominant in determining the timing and severity of necrosis.

A rare single nucleotide variant in Pm5e confers powdery mildew resistance in common wheat

Powdery mildew poses severe threats to wheat production. The most sustainable way to control this disease is through planting resistant cultivars. We report the map-based cloning of the powdery mildew resistance allele Pm5e from a Chinese wheat landrace. We applied a two-step bulked segregant RNA sequencing (BSR-Seq) approach in developing tightly linked or co-segregating markers to Pm5e. The first BSR-Seq used phenotypically contrasting bulks of recombinant inbred lines (RILs) to identify Pm5e-linked markers. The second BSR-Seq utilized bulks of genetic recombinants screened from a fine-mapping population to precisely quantify the associated genomic variation in the mapping interval, and identified the Pm5e candidate genes. The function of Pm5e was validated by transgenic assay, loss-of-function mutants and haplotype association analysis. Pm5e encodes a nucleotide-binding domain leucine-rich-repeat-containing (NLR) protein. A rare nonsynonymous single nucleotide variant (SNV) within the C-terminal leucine rich repeat (LRR) domain is responsible for the gain of powdery mildew resistance function of Pm5e, an allele endemic to wheat landraces of Shaanxi province of China. Results from this study demonstrate the value of landraces in discovering useful genes for modern wheat breeding. The key SNV associated with powdery mildew resistance will be useful for marker-assisted selection of Pm5e in wheat breeding programs.

Identification and genetic analysis of EMS-mutagenized wheat mutants conferring lesion-mimic premature aging

BACKGROUND

Lesion-mimic and premature aging (lmpa) mutant lmpa1 was identified from the ethyl methane sulfonate (EMS) mutant library in the bread wheat variety Keda 527 (KD527) background. To reveal the genetic basis of lmpa1 mutant, phenotypic observations and analyses of chlorophyll content and photosynthesis were carried out in lmpa1, KD527 and their F1 and F2 derivatives. Further, bulked segregation analysis (BSA) in combination with a 660 K SNP array were conducted on the F2 segregation population of lmpa1/Chinese spring (CS) to locate the lmpa1 gene.

RESULTS

Most agronomic traits of lmpa1 were similar to those of KD527 before lesion-like spots appeared. Genetic analysis indicated that the F1 plants from the crossing of lmpa1 and KD527 exhibited the lmpa phenotype and the F2 progenies showed a segregation of normal (wild type, WT) and lmpa, with the ratios of lmpa: WT = 124:36(χ2 = 1.008 < =3.841), indicating that lmpa is a dominant mutation. The combination of BSA and the SNP array analysis of CS, lmpa1 and lmpa1/CS F2 WT pool (50 plants) and lmpa pool (50 plants) showed that polymorphic SNPs were enriched on chromosome 5A, within a region of 30-40 Mb, indicating that the wheat premature aging gene Lmpa1 was probably located on the short arm of chromosome 5A.

CONCLUSIONS

EMS-mutagenized mutant lmpa1 deriving from elite wheat line KD527 conferred lmpa. Lmpa phenotype of lmpa1 mutant is controlled by a single dominant allele designated as Lmpa1, which affected wheat growth and development and reduced the thousand grain weight (tgw) of single plant in wheat. The gene Lmpa1 was tentatively located within the region of 30-40 Mb near to the short arm of chromosome 5A.

The Battle to Sequence the Bread Wheat Genome: A Tale of the Three Kingdoms

In the year 2018, the world witnessed the finale of the race to sequence the genome of the world's most widely grown crop, the common wheat. Wheat has been known to bear a notoriously large and complicated genome of a polyploidy nature. A decade competition to sequence the wheat genome initiated with a single consortium of multiple countries, taking a conventional strategy similar to that for sequencing Arabidopsis and rice, became ferocious over time as both sequencing technologies and genome assembling methodologies advanced. At different stages, multiple versions of genome sequences of the same variety (e.g., Chinese Spring) were produced by several groups with their special strategies. Finally, 16 years after the rice genome was finished and 9 years after that of maize, the wheat research community now possesses its own reference genome. Armed with these genomics tools, wheat will reestablish itself as a model for polyploid plants in studying the mechanisms of polyploidy evolution, domestication, genetic and epigenetic regulation of homoeolog expression, as well as defining its genetic diversity and breeding on the genome level. The enhanced resolution of the wheat genome should also help accelerate development of wheat cultivars that are more tolerant to biotic and/or abiotic stresses with better quality and higher yield.

High-Resolution Mapping of the Novel Early Leaf Senescence Gene Els2 in Common Wheat

Early leaf senescence negatively impacts the grain yield in wheat (Triticum aestivum L.). Induced mutants provide an important resource for mapping and cloning of genes for early leaf senescence. In our previous study, Els2, a single incomplete dominance gene, that caused early leaf senescence phenotype in the wheat mutant LF2099, had been mapped on the long arm of chromosome 2B. The objective of this study was to develop molecular markers tightly linked to the Els2 gene and construct a high-resolution map surrounding the Els2 gene. Three tightly linked single-nucleotide polymorphism (SNP) markers were obtained from the Illumina Wheat 90K iSelect SNP genotyping array and converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers. To saturate the Els2 region, the Axiom® Wheat 660K SNP array was used to screen bulked extreme phenotype DNA pools, and 9 KASP markers were developed. For fine mapping of the Els2 gene, these KASP markers and previously identified polymorphic markers were analyzed in a large F₂ population of the LF2099 × Chinese Spring cross. The Els2 gene was located in a 0.24-cM genetic region flanked by the KASP markers AX-111643885 and AX-111128667, which corresponded to a physical interval of 1.61 Mb in the Chinese Spring chromosome 2BL containing 27 predicted genes with high confidence. The study laid a foundation for a map-based clone of the Els2 gene controlling the mutation phenotype and revealing the molecular regulatory mechanism of wheat leaf senescence.

Appraisal of wheat genomics for gene discovery and breeding applications: a special emphasis on advances in Asia

We discussed the most recent efforts in wheat functional genomics to discover new genes and their deployment in breeding with special emphasis on advances in Asian countries. Wheat research community is making significant progress to bridge genotype-to-phenotype gap and then applying this knowledge in genetic improvement. The advances in genomics and phenomics have intrigued wheat researchers in Asia to make best use of this knowledge in gene and trait discovery. These advancements include, but not limited to, map-based gene cloning, translational genomics, gene mapping, association genetics, gene editing and genomic selection. We reviewed more than 57 homeologous genes discovered underpinning important traits and multiple strategies used for their discovery. Further, the complementary advancements in wheat phenomics and analytical approaches to understand the genetics of wheat adaptability, resilience to climate extremes and resistance to pest and diseases were discussed. The challenge to build a gold standard reference genome sequence of bread wheat is now achieved and several de novo reference sequences from the cultivars representing different gene pools will be available soon. New pan-genome sequencing resources of wheat will strengthen the foundation required for accelerated gene discovery and provide more opportunities to practice the knowledge-based breeding.

Alternative splicing of a barley gene results in an excess-tillering and semi-dwarf mutant

An excess-tillering semi-dwarf gene Hvhtd was identified from an EMS-induced mutant in barley and alternative splicing results in excess-tillering semi-dwarf traits. Tillering and plant height are important traits determining plant architecture and grain production in cereal crops. This study identified an excess-tillering semi-dwarf mutant (htd) from an EMS-treated barley population. Genetic analysis of the F₁, F₂, and F_2:3 populations showed that a single recessive gene controlled the excess-tillering semi-dwarf in htd. Using BSR-Seq and gene mapping, the Hvhtd gene was delimited within a 1.8 Mb interval on chromosome 2HL. Alignment of the RNA-Seq data with the functional genes in the interval identified a gene HORVU2Hr1G098820 with alternative splicing between exon2 and exon3 in the mutant, due to a G to A single-nucleotide substitution at the exon and intron junction. An independent mutant with a similar phenotype confirmed the result, with alternative splicing between exon3 and exon4. In both cases, the alternative splicing resulted in a non-functional protein. And the gene HORVU2Hr1G098820 encodes a trypsin family protein and may be involved in the IAA signaling pathway and differs from the mechanism of Green Revolution genes in the gibberellic acid metabolic pathway.

Current achievements and future prospects in the genetic breeding of chrysanthemum: a review

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a leading flower with applied value worldwide. Developing new chrysanthemum cultivars with novel characteristics such as new flower colors and shapes, plant architectures, flowering times, postharvest quality, and biotic and abiotic stress tolerance in a time- and cost-efficient manner is the ultimate goal for breeders. Various breeding strategies have been employed to improve the aforementioned traits, ranging from conventional techniques, including crossbreeding and mutation breeding, to a series of molecular breeding methods, including transgenic technology, genome editing, and marker-assisted selection (MAS). In addition, the recent extensive advances in high-throughput technologies, especially genomics, transcriptomics, proteomics, metabolomics, and microbiomics, which are collectively referred to as omics platforms, have led to the collection of substantial amounts of data. Integration of these omics data with phenotypic information will enable the identification of genes/pathways responsible for important traits. Several attempts have been made to use emerging molecular and omics methods with the aim of accelerating the breeding of chrysanthemum. However, applying the findings of such studies to practical chrysanthemum breeding remains a considerable challenge, primarily due to the high heterozygosity and polyploidy of the species. This review summarizes the recent achievements in conventional and modern molecular breeding methods and emerging omics technologies and discusses their future applications for improving the agronomic and horticultural characteristics of chrysanthemum.

Rapid identification of a candidate nicosulfuron sensitivity gene (Nss) in maize (Zea mays L.) via combining bulked segregant analysis and RNA-seq

A candidate nicosulfuron sensitivity gene Nss was identified by combining bulked segregant analysis and RNA-seq. Multiple mutations of this gene were discovered in nicosulfuron-sensitive maize compared with the tolerant. It has been demonstrated that variabilities exist in maize response to nicosulfuron. Two nicosulfuron-sensitive inbred lines (HB39, HB41) and two tolerant inbred lines (HB05, HB09) were identified via greenhouse and field trials. Genetic analysis indicated that the sensitivity to nicosulfuron in maize was controlled by a single, recessive gene. To precisely and rapidly map the nicosulfuron sensitivity gene (Nss), two independent F₂ segregating populations, Population A (HB41 × HB09) and Population B (HB39 × HB05), were constructed. By applying bulked segregant RNA-Seq (BSR-Seq), the Nss gene was, respectively, mapped on the short arm of chromosome 5 (chr5: 1.1-15.3 Mb) and (chr5: 0.5-18.2 Mb) using two populations, with 14.2 Mb region in common. Further analysis revealed that there were 43 and 119 differentially expressed genes in the mapping intervals, with 18 genes in common. Gene annotation results showed that a cytochrome P450 gene (CYP81A9) appeared to be the candidate gene of Nss associated with nicosulfuron sensitivity in maize. Sequence analysis demonstrated that two common deletion mutations existed in the sensitive maize, which might lead to the nicosulfuron sensitivity in maize. The results will make valuable contributions to the understanding of molecular mechanism of herbicide sensitivity in maize.