Bulked Segregant RNA-seq Reveals Differential Expression and SNPs of Candidate Genes Associated with Waterlogging Tolerance in Maize

The putative vacuolar processing enzyme gene TaVPE3cB is a candidate gene for wheat stem pith-thickness

KEY MESSAGE

The vacuolar processing enzyme gene TaVPE3cB is identified as a candidate gene for a QTL of wheat pith-thickness on chromosome 3B by BSR-seq and differential expression analyses. The high pith-thickness (PT) of the wheat stem could greatly enhance stem mechanical strength, especially the basal internodes which support the heavier upper part, such as upper stems, leaves and spikes. A QTL for PT in wheat was previously discovered on 3BL in a double haploid population of 'Westonia' × 'Kauz'. Here, a bulked segregant RNA-seq analysis was applied to identify candidate genes and develop associated SNP markers for PT. In this study, we aimed at screening differentially expressed genes (DEGs) and SNPs in the 3BL QTL interval. Sixteen DEGs were obtained based on BSR-seq and differential expression analyses. Twenty-four high-probability SNPs in eight genes were identified by comparing the allelic polymorphism in mRNA sequences between the high PT and low PT samples. Among them, six genes were confirmed to be associated with PT by qRT-PCR and sequencing. A putative vacuolar processing enzyme gene TaVPE3cB was screened out as a potential PT candidate gene in Australian wheat 'Westonia'. A robust SNP marker associated with TaVPE3cB was developed, which can assist in the introgression of TaVPE3cB.b in wheat breeding programs. In addition, we also discussed the function of other DEGs which may be related to pith development and programmed cell death (PCD). A five-level hierarchical regulation mechanism of stem pith PCD in wheat was proposed.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18 Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

QTL and Candidate Genes: Techniques and Advancement in Abiotic Stress Resistance Breeding of Major Cereals

At least 75% of the world's grain production comes from the three most important cereal crops: rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays). However, abiotic stressors such as heavy metal toxicity, salinity, low temperatures, and drought are all significant hazards to the growth and development of these grains. Quantitative trait locus (QTL) discovery and mapping have enhanced agricultural production and output by enabling plant breeders to better comprehend abiotic stress tolerance processes in cereals. Molecular markers and stable QTL are important for molecular breeding and candidate gene discovery, which may be utilized in transgenic or molecular introgression. Researchers can now study synteny between rice, maize, and wheat to gain a better understanding of the relationships between the QTL or genes that are important for a particular stress adaptation and phenotypic improvement in these cereals from analyzing reports on QTL and candidate genes. An overview of constitutive QTL, adaptive QTL, and significant stable multi-environment and multi-trait QTL is provided in this article as a solid framework for use and knowledge in genetic enhancement. Several QTL, such as DRO1 and Saltol, and other significant success cases are discussed in this review. We have highlighted techniques and advancements for abiotic stress tolerance breeding programs in cereals, the challenges encountered in introgressing beneficial QTL using traditional breeding techniques such as mutation breeding and marker-assisted selection (MAS), and the in roads made by new breeding methods such as genome-wide association studies (GWASs), the clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system, and meta-QTL (MQTL) analysis. A combination of these conventional and modern breeding approaches can be used to apply the QTL and candidate gene information in genetic improvement of cereals against abiotic stresses.

Screening of Candidate Genes Associated with Brown Stripe Resistance in Sugarcane via BSR-seq Analysis

Sugarcane brown stripe (SBS), caused by the fungal pathogen Helminthosporium stenospilum, is one of the most serious threats to sugarcane production. However, its outbreaks and epidemics require suitable climatic conditions, resulting in the inefficient improvement of the SBS resistance by phenotype selection. The sugarcane F₁ population of SBS-resistant YT93-159 × SBS-susceptible ROC22 was used for constructing the bulks. Bulked segregant RNA-seq (BSR-seq) was then performed on the parents YT93-159 (T01) and ROC22 (T02), and the opposite bulks of 30 SBS-susceptible individuals mixed bulk (T03) and 30 SBS-resistant individuals mixed bulk (T04) collected from 287 F₁ individuals. A total of 170.00 Gb of clean data containing 297,921 SNPs and 70,426 genes were obtained. Differentially expressed genes (DEGs) analysis suggested that 7787 and 5911 DEGs were identified in the parents (T01 vs. T02) and two mixed bulks (T03 vs. T04), respectively. In addition, 25,363 high-quality and credible SNPs were obtained using the genome analysis toolkit GATK for SNP calling. Subsequently, six candidate regions with a total length of 8.72 Mb, which were located in the chromosomes 4B and 7C of sugarcane wild species Saccharum spontaneum, were identified, and 279 genes associated with SBS-resistance were annotated by ED algorithm and ΔSNP-index. Furthermore, the expression profiles of candidate genes were verified by quantitative real-time PCR (qRT-PCR) analysis, and the results showed that eight genes (LRR-RLK, DHAR1, WRKY7, RLK1, BLH4, AK3, CRK34, and NDA2) and seven genes (WRKY31, CIPK2, CKA1, CDPK6, PFK4, CBL2, and PR2) of the 20 tested genes were significantly up-regulated in YT93-159 and ROC22, respectively. Finally, a potential molecular mechanism of sugarcane response to H. stenospilum infection is illustrate that the activations of ROS signaling, MAPK cascade signaling, Ca²⁺ signaling, ABA signaling, and the ASA-GSH cycle jointly promote the SBS resistance in sugarcane. This study provides abundant gene resources for the SBS resistance breeding in sugarcane.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Expressing the sunflower transcription factor HaHB11 in maize improves waterlogging and defoliation tolerance

The sunflower (Helianthus annuus) transcription factor HaHB11 (H. annuus  Homeobox 11) belongs to the homeodomain-leucine zipper family and confers improved yield to maize (Zea mays) hybrids (HiII × B73) and lines. Here we report that transgenic maize lines expressing HaHB11 exhibited better performance under waterlogging, both in greenhouse and field trials carried out during three growth cycles. Transgenic plants had increased chlorophyll content, wider stems, more nodal roots, greater total aerial biomass, a higher harvest index, and increased plant grain yield. Under severe defoliation caused by a windstorm during flowering, transgenic genotypes were able to set more grains than controls. This response was confirmed in controlled defoliation assays. Hybrids generated by crossing B73 HaHB11 lines with the contrasting Mo17 lines were also tested in the field and exhibited the same beneficial traits as the parental lines, compared with their respective controls. Moreover, they were less penalized by stress than commercial hybrids. Waterlogging tolerance increased via improvement of the root system, including more xylem vessels, reduced tissue damage, less superoxide accumulation, and altered carbohydrate metabolism. Multivariate analyses corroborated the robustness of the differential traits observed. Furthermore, canopy spectral reflectance data, computing 29 vegetation indices associated with biomass, chlorophyll, and abiotic stress, helped to distinguish genotypes as well as their growing conditions. Altogether the results reported here indicate that this sunflower gene constitutes a suitable tool to improve maize plants for environments prone to waterlogging and/or wind defoliation.

Genomics and Pathways Involved in Maize Resistance to Fusarium Ear Rot and Kernel Contamination With Fumonisins

Fusarium verticillioides is a causal agent of maize ear rot and produces fumonisins, which are mycotoxins that are toxic to animals and humans. In this study, quantitative trait loci (QTLs) and bulk-segregant RNA-seq approaches were used to uncover genomic regions and pathways involved in resistance to Fusarium ear rot (FER) and to fumonisin accumulation in maize kernels. Genomic regions at bins 4.07-4.1, 6-6.01, 6.04-6.05, and 8.05-8.08 were related to FER resistance and/or reduced fumonisin levels in kernels. A comparison of transcriptomes between resistant and susceptible inbred bulks 10 days after inoculation with F. verticillioides revealed 364 differentially expressed genes (DEGs). In the resistant inbred bulks, genes involved in sink metabolic processes such as fatty acid and starch biosynthesis were downregulated, as well as those involved in phytosulfokine signaling and many other genes involved in cell division; while genes involved in secondary metabolism and compounds/processes related to resistance were upregulated, especially those related to cell wall biosynthesis/rearrangement and flavonoid biosynthesis. These trends are indicative of a growth-defense trade-off. Among the DEGs, Zm00001d053603, Zm00001d035562, Zm00001d037810, Zm00001d037921, and Zm00001d010840 were polymorphic between resistant and susceptible bulks, were located in the confidence intervals of detected QTLs, and showed large differences in transcript levels between the resistant and susceptible bulks. Thus, they were identified as candidate genes involved in resistance to FER and/or reduced fumonisin accumulation.

Bulked Segregant RNA Sequencing Revealed Difference Between Virulent and Avirulent Brown Planthoppers

The brown planthopper (Nilaparvata lugens Stål, BPH) is one of the most devastating insect pests of rice (Oryza sativa L.), but BPH populations have varying degrees of virulence to rice varieties carrying different resistance genes. To help efforts to characterize these variations we applied bulked segregant RNA sequencing (BSR-seq) to identify differentially expressed genes (DEGs) and genetic loci associated with BPH virulence to YHY15 rice plants carrying the resistance gene Bph15. BPHs that are highly virulent or avirulent to these plants were selected from an F2 population to form two contrasting bulks, and BSR-seq identified 751 DEGs between the bulks. Genes associated with carbohydrate, amino acid and nucleotide metabolism, the endocrine system, and signal transduction were upregulated in the avirulent insects when they fed on these plants. The results also indicated that shifts in lipid metabolism and digestive system pathways were crucial for the virulent BPHs' adaptation to the resistant rice. We identified 24 single-nucleotide polymorphisms (SNPs) in 21 genes linked with BPH virulence. Possible roles of genes apparently linked to BPH virulence are discussed. Our results provide potentially valuable information for further studies of BPH virulence mechanisms and development of robust control strategies.

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus. The yellow-seed trait (YST) has always been a main breeding objective for rapeseed because yellow-seeded B. napus generally contains higher oil contents, fewer pigments and polyphenols and lower fiber content than black-seeded B. napus, although the mechanism controlling this correlation remains unclear. In this study, QTL mapping was implemented for YST based on a KN double haploid population derived from the hybridization of yellow-seeded B. napus N53-2 with a high oil content and black-seeded Ken-C8 with a relatively low oil content. Ten QTLs were identified, including four stable QTLs that could be detected in multiple environments. A major QTL, cqSC-A09, on chromosome A09 was identified by both QTL mapping and BSR-Seq technology, and explained more than 41% of the phenotypic variance. The major QTL cqSC-A09 for YST not only controls the seed color but also affects the oil and fiber contents in seeds. More importantly, the advantageous allele could increase the oil content and reduce the pigment and fiber content at the same time. This is the first QTL reported to control seed color, oil content and fiber content simultaneously with a large effect and has great application value for breeding high oil varieties with high seed quality. Important candidate genes, including BnaA09. JAZ1, BnaA09. GH3.3 and BnaA09. LOX3, were identified for cqSC-A09 by combining sequence variation annotation, expression differences and an interaction network, which lays a foundation for further cloning and breeding applications in the future.

Poplar's Waterlogging Resistance Modeling and Evaluating: Exploring and Perfecting the Feasibility of Machine Learning Methods in Plant Science

Floods, as one of the most common disasters in the natural environment, have caused huge losses to human life and property. Predicting the flood resistance of poplar can effectively help researchers select seedlings scientifically and resist floods precisely. Using machine learning algorithms, models of poplar's waterlogging tolerance were established and evaluated. First of all, the evaluation indexes of poplar's waterlogging tolerance were analyzed and determined. Then, significance testing, correlation analysis, and three feature selection algorithms (Hierarchical clustering, Lasso, and Stepwise regression) were used to screen photosynthesis, chlorophyll fluorescence, and environmental parameters. Based on this, four machine learning methods, BP neural network regression (BPR), extreme learning machine regression (ELMR), support vector regression (SVR), and random forest regression (RFR) were used to predict the flood resistance of poplar. The results show that random forest regression (RFR) and support vector regression (SVR) have high precision. On the test set, the coefficient of determination (R²) is 0.8351 and 0.6864, the root mean square error (RMSE) is 0.2016 and 0.2780, and the mean absolute error (MAE) is 0.1782 and 0.2031, respectively. Therefore, random forest regression (RFR) and support vector regression (SVR) can be given priority to predict poplar flood resistance.

iTRAQ Proteomic Analysis of Wheat (Triticum aestivum L.) Genotypes Differing in Waterlogging Tolerance

Transient and chronic waterlogging constrains crop production in many regions of the world. Here, we invoke a novel iTRAQ-based proteomic strategy to elicit protein synthesis and regulation responses to waterlogging in tolerant (XM 55) and sensitive genotypes (YM 158). Of the 7,710 proteins identified, 16 were distinct between the two genotypes under waterlogging, partially defining a proteomic basis for waterlogging tolerance (and sensitivity). We found that 11 proteins were up-regulated and 5 proteins were down-regulated; the former included an Fe-S cluster assembly factor, heat shock cognate 70, GTP-binding protein SAR1A-like and CBS domain-containing protein. Down-regulated proteins contained photosystem II reaction center protein H, carotenoid 9, 10 (9', 10')-cleavage dioxygenase-like, psbP-like protein 1 and mitochondrial ATPase inhibitor. We showed that nine proteins responded to waterlogging with non-cultivar specificity: these included 3-isopropylmalate dehydratase large subunit, solanesyl-diphosphate synthase 2, DEAD-box ATP-dependent RNA helicase 3, and 3 predicted or uncharacterized proteins. Sixteen of the 28 selected proteins showed consistent expression patterns between mRNA and protein levels. We conclude that waterlogging stress may redirect protein synthesis, reduce chlorophyll synthesis and enzyme abundance involved in photorespiration, thus influencing synthesis of other metabolic enzymes. Collectively, these factors accelerate the accumulation of harmful metabolites in leaves in waterlogging-susceptible genotypes. The differentially expressed proteins enumerated here could be used as biological markers for enhancing waterlogging tolerance as part of future crop breeding programs.

Identification of Candidate Genes Associated With Hypoxia Tolerance in Trachinotus blochii Using Bulked Segregant Analysis and RNA-Seq

Golden Pompano (Trachinotus blochii) has rapidly developed into the one of the main valuable fish species in Chinese marine aquaculture. Due to its rapid growth, active metabolism, and high oxygen consumption, hypoxia will increase its mortality and cause serious economic losses. We constructed two experimental groups of fish with different degrees of tolerance to hypoxia, used BSR-Seq analysis based on genome and genetic linkage groups to locate SNPs and genes that were related to the differences in hypoxia tolerance. The results showed that hypoxia tolerance SNPs of golden pompano may be jointly determined by multiple linkage groups, especially linkage groups 18 and 22. There were 768 and 348 candidate genes located in the candidate regions of the brain and liver, respectively. These genes were mainly involved in anaerobic energy metabolism, stress response, immune response, waste discharge, and cell death. The prostaglandin-endoperoxide synthase 2 (PTGS2) on LG8, which is involved in the metabolism of arachidonic acid, has a G/A nonsynonymous mutation at position 20641628, and the encoded amino acid was changed from hydrophobic aspartic acid to asparaginate. The specific pathway of the RIG-I-like receptor signaling pathway in the liver may mediate the metabolic system and the immune system, linking glucose metabolism with immune regulation. The death of the hypoxia-intolerant group may be due to the accumulation of lactic acid caused by the activation of anaerobic glycolysis during the early stage of hypoxia stress, and the activation of type I interferon was inhibited, which resulted in decreased immunity. Among the genes involved in the RIG-I-like receptor signaling pathway, the CYLD Lysine 63 Deubiquitinase (CYLD) located on LG16 had a G/T nonsynonymous mutation at position 13629651, and the encoded amino acid was changed from alanine acid to valine. The interferon induced with helicase C domain 1 (Ifih1) located on LG18 has a G/C nonsynonymous mutation at position 16153700, and the encoded hydrophilic glycine was changed to hydrophobic alanine. Our findings suggest these SNPs may assist in the molecular breeding of hypoxia-tolerant golden pompano, and speculate that the balance of glucose and lipid metabolism plays a key role in Trachinotus blochii under acute hypoxia.

Opportunities for Improving Waterlogging Tolerance in Cereal Crops-Physiological Traits and Genetic Mechanisms

Waterlogging occurs when soil is saturated with water, leading to anaerobic conditions in the root zone of plants. Climate change is increasing the frequency of waterlogging events, resulting in considerable crop losses. Plants respond to waterlogging stress by adventitious root growth, aerenchyma formation, energy metabolism, and phytohormone signalling. Genotypes differ in biomass reduction, photosynthesis rate, adventitious roots development, and aerenchyma formation in response to waterlogging. We reviewed the detrimental effects of waterlogging on physiological and genetic mechanisms in four major cereal crops (rice, maize, wheat, and barley). The review covers current knowledge on waterlogging tolerance mechanism, genes, and quantitative trait loci (QTL) associated with waterlogging tolerance-related traits, the conventional and modern breeding methods used in developing waterlogging tolerant germplasm. Lastly, we describe candidate genes controlling waterlogging tolerance identified in model plants Arabidopsis and rice to identify homologous genes in the less waterlogging-tolerant maize, wheat, and barley.

Identification of major QTL for waterlogging tolerance in maize using genome-wide association study and bulked sample analysis

Waterlogging has increasingly become one of the major constraints to maize (Zea mays L.) production in some maize growing areas as it seriously decreases the yield. Waterlogging tolerance in maize germplasm provides a basis for maize waterlogging improvement. In this study, nine seedling traits, plant height (PH), root length (RL), shoot dry weight (SDW), root dry weight (RDW), adventitious root number (ARN), node number of brace root (BRNN), brace root number (BRN), brace root dry weigh (BRDW), survival rate (SR), and the secondary traits that were defined as relative phenotypic value of seedling traits under waterlogging and control treatments were used in a natural population that contain 365 inbred lines to evaluate the waterlogging tolerance of tropical maize. The result showed that maize waterlogging tolerance was genetically controlled and seedling traits were significantly different between the control and waterlogging treatments. PH, RL, SDW, and RDW are important seedling traits for waterlogging tolerance identification. Some tropical maize inbred lines were identified with extreme waterlogging tolerance that can provide an important germplasm resource for breeding. Population structure analysis showed that two major phylogenetic subgroups in tropical maize could be identified. Genome-wide association study (GWAS) using 39,266 single nucleotide polymorphisms (SNPs) across the whole genome identified 49 trait-SNPs distributed on over all 10 chromosomes excluding chromosome 10. Seventy-one significant SNPs, distributed on all 10 chromosomes excluding chromosome 5, were identified by extend bulked sample analysis (Ext-BSA) based on the inbred lines with extreme phenotypes. GWAS and Ext-BSA identified the same loci on bin1.07, bin6.01, bin2.09, bin6.04, bin7.02, and bin7.03. Nine genes were proposed as potential candidate genes. Cloning and functional validation of these genes would be helpful for understanding the molecular mechanism of waterlogging tolerance in maize.

Identification and validation of genetic locus Rfk1 for wheat fertility restoration in the presence of Aegilops kotschyi cytoplasm

Major fertility restorer locus for Aegilops kotschyi cytoplasm in wheat, Rfk1, was mapped to chromosome arm 1BS. Most likely candidate gene is TraesCS1B02G197400LC, which is predicted to encode a pectinesterase/pectinesterase inhibitor. Cytoplasmic male sterility (CMS) is widely used for heterosis and hybrid seed production in wheat. Genes related to male fertility restoration in the presence of Aegilops kotschyi cytoplasm have been reported, but the fertility restoration-associated gene loci have not been investigated systematically. In this study, a BC₁F₁ population derived from a backcross between KTP116A, its maintainer line TP116B, and its restorer line LK783 was employed to map fertility restoration by bulked segregant RNA-Seq (BSR-Seq). A major fertility allele restorer locus for Ae. kotschyi cytoplasm in wheat, Rfk1, was mapped to chromosome arm 1BS, and it was contributed by LK783. Morphological and cytological studies showed that male fertility restoration occurred mainly after the late uninucleate stage. Based on simple sequence repeat and single-nucleotide polymorphism genotyping, the gene locus was located between Xnwafu_6 and Xbarc137 on chromosome arm 1BS. To further isolate the specific region, six Kompetitive allele-specific polymerase chain reaction markers derived from BSR-Seq were developed to delimit Rfk1 within physical intervals of 26.0 Mb. After searching for differentially expressed genes within the candidate interval in the anthers and sequencing analysis, TraesCS1B02G197400LC was identified as a candidate gene for Rfk1 and it was predicted to encode a pectinesterase/pectinesterase inhibitor. Expression analysis also confirmed that it was specifically expressed in the anthers, and its expression level was higher in fertile lines compared with sterile lines. Thus, TraesCS1B02G197400LC was identified as the most likely candidate gene for Rfk1, thereby providing insights into the fertility restoration mechanism for K-type CMS in wheat.

BSR-Seq analysis provides insights into the cold stress response of Actinidia arguta F1 populations

BACKGROUND

Freezing injury, which is an important abiotic stress in horticultural crops, influences the growth and development and the production area of kiwifruit (Actinidia Lind1). Among Actinidia species, Actinidia arguta has excellent cold resistance, but knowledge relevant to molecular mechanisms is still limited. Understanding the mechanism underlying cold resistance in kiwifruit is important for breeding cold resistance.

RESULTS

In our study, a population resulting from the cross of A. arguta 'Ruby-3' × 'Kuilv' male was generated for kiwifruit hardiness study, and 20 cold-tolerant and 20 cold-sensitive populations were selected from 492 populations according to their LT50. Then, we performed bulked segregant RNA-seq combined with single-molecule real-time sequencing to identify differentially expressed genes that provide cold hardiness. We found that the content of soluble sucrose and the activity of β-amylase were higher in the cold-tolerant population than in the cold-sensitive population. Upon - 30 °C low-temperature treatment, 126 differentially expressed genes were identify; the expression of 59 genes was up-regulated and that of 67 genes was down-regulated between the tolerant and sensitive pools, respectively. KEGG pathway analysis showed that the DEGs were primarily related to starch and sucrose metabolism, amino sugar and nucleotide sugar metabolism. Ten major key enzyme-encoding genes and two regulatory genes were up-regulated in the tolerant pool, and regulatory genes of the CBF pathway were found to be differentially expressed. In particular, a 14-3-3 gene was down-regulated and an EBF gene was up-regulated. To validate the BSR-Seq results, 24 DEGs were assessed via qRT-PCR, and the results were consistent with those obtained by BSR-Seq.

CONCLUSION

Our research provides valuable insights into the mechanism related to cold resistance in Actinidia and identified potential genes that are important for cold resistance in kiwifruit.

A temporal hierarchy underpins the transcription factor-DNA interactome of the maize UPR

Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large-scale systems-level analysis of the protein-DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high-throughput enhanced yeast one-hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis-regulatory elements (CREs) in ER stress-responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Finally, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription-regulatory mechanisms in the stress responses and crop improvement.

Repressed Gene Expression of Photosynthetic Antenna Proteins Associated with Yellow Leaf Variation as Revealed by Bulked Segregant RNA-seq in Tea Plant Camellia sinensis

The young leaves and shoots of albino tea cultivars are usually characterized as having a yellow or pale color, high amino acid, and low catechin. Increasing attention has been paid to albino tea cultivars in recent years because their tea generally shows high umami and reduced astringency. However, the genetic mechanism of yellow-leaf variation in albino tea cultivar has not been elucidated clearly. In this study, bulked segregant RNA-seq (BSR-seq) was performed on bulked yellow- and green-leaf hybrid progenies from a leaf color variation population. A total of 359 and 1134 differentially expressed genes (DEGs) were identified in the yellow and green hybrid bulked groups (Y_f vs G_f) and parent plants (Y_p vs G_p), respectively. The significantly smaller number of DEGs in Y_f versus G_f than in Y_p versus G_p indicated that individual differences could be reduced within the same hybrid progeny. Analysis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes revealed that the photosynthetic antenna protein was most significantly enriched in either the bulked groups or their parents. Interaction was found among light-harvesting chlorophyll a/b -binding proteins (LHC), heat shock proteins (HSPs), and enzymes involved in cuticle formation. Combined with the transcriptomic expression profile, results showed that the repressed genes encoding LHC were closely linked to aberrant chloroplast development in yellow-leaf tea plants. Furthermore, the photoprotection and light stress response possessed by genes involved in HSP protein interaction and cuticle formation were discussed. The expression profile of DEGs was verified via quantitative real-time PCR analysis of the bulked samples and other F₁ individuals. In summary, using BSR-seq on a hybrid population eliminated certain disturbing effects of genetic background and individual discrepancy, thereby helping this study to intensively focus on the key genes controlling leaf color variation in yellow-leaf tea plants.

Dissecting the genetics of cold tolerance in a multiparental maize population

We identify the largest amount of QTLs for cold tolerance in maize; mainly associated with photosynthetic efficiency, which opens new possibilities for genomic selection for cold tolerance in maize. Breeding for cold tolerance in maize is an important objective in temperate areas. The objective was to carry out a highly efficient study of quantitative trait loci (QTLs) for cold tolerance in maize. We evaluated 406 recombinant inbred lines from a multi-parent advanced generation intercross (MAGIC) population in a growth chamber under cold and control conditions, and in the field at early and normal sowing. We recorded cold tolerance-related traits, including the number of days from sowing to emergence, chlorophyll content and maximum quantum efficiency of photosystem II (F_v/F_m). Association mapping was based on genotyping with near one million single nucleotide polymorphism (SNP) markers. We found 858 SNPs significantly associated with all traits, most of them under cold conditions and early sowing. Most QTLs were associated with chlorophyll and F_v/F_m. Many candidate genes coincided between the current research and previous reports. These results suggest that (1) the MAGIC population is an efficient tool for identifying QTLs for cold tolerance; (2) most QTLs for cold tolerance were associated with F_v/F_m; (3) most of these QTLs were located in specific genomic regions, particularly bin 10.04; (4) the current study allows genetically improving cold tolerance with genome-wide selection.

A predicted NEDD8 conjugating enzyme gene identified as a Capsicum candidate Rf gene using bulk segregant RNA sequencing

Cytoplasmic male sterility (CMS) is an important tool for producing F₁ hybrids, which can exhibit heterosis. The companion system, restorer-of-fertility (Rf), is poorly understood at the molecular level and would be valuable in producing restorer lines for hybrid seed production. The identity of the Rf gene in Capsicum (pepper) is currently unclear. In this study, using bulked segregant RNA sequencing (BSR-seq), a strong candidate Rf gene, Capana06g002866, which is annotated as a NEDD8 conjugating enzyme E2, was identified. Capana06g002866 has an ORF of 555 bp in length encoding 184 amino acids; it can be cloned from F₁ plants from the hybridization of the CMS line 8A and restorer line R1 but is not found in CMS line 8A. With qRT-PCR validation, Capana06g002866 was found to be upregulated in restorer accessions compared to sterile accessions. The relative expression in flower buds increased with the developmental stage in F₁ plants, while the expression was very low in all flower bud stages of the CMS lines. These results provide new insights into the Rf gene in pepper and will be useful for other crops utilizing the CMS system.

Comparative transcriptional analysis of Capsicum flower buds between a sterile flower pool and a restorer flower pool provides insight into the regulation of fertility restoration

BACKGROUND

Cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system is an important mechanism to produce F₁ hybrid seeds. Understanding the interaction that controls restoration at a molecular level will benefit plant breeders. The CMS is caused by the interaction between mitochondrial and nuclear genes, with the CMS phenotype failing to produce functional anthers, pollen, or male gametes. Thus, understanding the complex processes of anther and pollen development is a prerequisite for understanding the CMS system. Currently it is accepted that the Rf gene in the nucleus restores the fertility of CMS, however the Rf gene has not been cloned. In this study, CMS line 8A and the Rf line R1, as well as a sterile pool (SP) of accessions and a restorer pool (RP) of accessions analyzed the differentially expressed genes (DEGs) between CMS and its fertility restorer using the conjunction of RNA sequencing and bulk segregation analysis.

RESULTS

A total of 2274 genes were up-regulated in R1 as compared to 8A, and 1490 genes were up-regulated in RP as compared to SP. There were 891 genes up-regulated in both restorer accessions, R1 and RP, as compared to both sterile accessions, 8A and SP. Through annotation and expression analysis of co-up-regulated expressed genes, eight genes related to fertility restoration were selected. These genes encode putative fructokinase, phosphatidylinositol 4-phosphate 5-kinase, pectate lyase, exopolygalacturonase, pectinesterase, cellulose synthase, fasciclin-like arabinogalactan protein and phosphoinositide phospholipase C. In addition, a phosphatidylinositol signaling system and an inositol phosphate metabolism related to the fertility restorer of CMS were ranked as the most likely pathway for affecting the restoration of fertility in pepper.

CONCLUSIONS

Our study revealed that eight genes were related to the restoration of fertility, which provides new insight into understanding the molecular mechanism of fertility restoration of CMS in Capsicum.

Single-plant GWAS coupled with bulk segregant analysis allows rapid identification and corroboration of plant-height candidate SNPs

BACKGROUND

Genome wide association studies (GWAS) are a powerful tool for identifying quantitative trait loci (QTL) and causal single nucleotide polymorphisms (SNPs)/genes associated with various important traits in crop species. Typically, GWAS in crops are performed using a panel of inbred lines, where multiple replicates of the same inbred are measured and the average phenotype is taken as the response variable. Here we describe and evaluate single plant GWAS (sp-GWAS) for performing a GWAS on individual plants, which does not require an association panel of inbreds. Instead sp-GWAS relies on the phenotypes and genotypes from individual plants sampled from a randomly mating population. Importantly, we demonstrate how sp-GWAS can be efficiently combined with a bulk segregant analysis (BSA) experiment to rapidly corroborate evidence for significant SNPs.

RESULTS

In this study we used the Shoepeg maize landrace, collected as an open pollinating variety from a farm in Southern Missouri in the 1960's, to evaluate whether sp-GWAS coupled with BSA can efficiently and powerfully used to detect significant association of SNPs for plant height (PH). Plant were grown in 8 locations across two years and in total 768 individuals were genotyped and phenotyped for sp-GWAS. A total of 306 k polymorphic markers in 768 individuals evaluated via association analysis detected 25 significant SNPs (P ≤ 0.00001) for PH. The results from our single-plant GWAS were further validated by bulk segregant analysis (BSA) for PH. BSA sequencing was performed on the same population by selecting tall and short plants as separate bulks. This approach identified 37 genomic regions for plant height. Of the 25 significant SNPs from GWAS, the three most significant SNPs co-localize with regions identified by BSA.

CONCLUSION

Overall, this study demonstrates that sp-GWAS coupled with BSA can be a useful tool for detecting significant SNPs and identifying candidate genes. This result is particularly useful for species/populations where association panels are not readily available.

A SNP Mutation of SiCRC Regulates Seed Number Per Capsule and Capsule Length of cs1 Mutant in Sesame

Seed number per capsule (SNC) is a major factor influencing seed yield and is an important trait with complex gene interaction effects. We first performed genetic analysis, gene cloning, and molecular mechanism study for an EMS-induced sesame mutant cs1 with fewer SNC and shorter capsule length (CL). The mutant traits were due to the pleiotropism of a regressive gene (Sics1). Capsule hormone determination showed that five out of 12 hormones, including auxin indole-3-acetic acid (IAA), had significantly different levels between wild type (WT) and mutant type (MT). KEGG pathway analysis showed that plant hormone signal transduction, especially the auxin signal transduction pathway, was the most abundant differentially expressed signaling pathway. After the cross-population association and regional genome screening, we found that three homozygous loci were retained in cs1. Further analysis of these three loci resulted in the identification of SiCRC as the candidate gene for cs1. SiCRC consists of seven exons and six introns encoding 163 amino acids. The SiCRC in cs1 showed a point mutation at intron 5 and exon 6 junction, resulting in the splice site being frame-shifted eight nucleotides further downstream, causing incorrect splicing. Taken together, we assumed the SNP mutation in SiCRC disrupted the function of the transcription factor, which might act downstream of the CRC-auxin signal transduction pathway, resulting in a shorter CL and less SNC mutation of cs1 in sesame. Our results highlight the molecular framework underlying the transcription factor CRC-mediated role of auxin transduction in SNC and CL development.

Submergence and Waterlogging Stress in Plants: A Review Highlighting Research Opportunities and Understudied Aspects

Soil flooding creates composite and complex stress in plants known as either submergence or waterlogging stress depending on the depth of the water table. In nature, these stresses are important factors dictating the species composition of the ecosystem. On agricultural land, they cause economic damage associated with long-term social consequences. The understanding of the plant molecular responses to these two stresses has benefited from research studying individual components of the stress, in particular low-oxygen stress. To a lesser extent, other associated stresses and plant responses have been incorporated into the molecular framework, such as ion and ROS signaling, pathogen susceptibility, and organ-specific expression and development. In this review, we aim to highlight known or suspected components of submergence/waterlogging stress that have not yet been thoroughly studied at the molecular level in this context, such as miRNA and retrotransposon expression, the influence of light/dark cycles, protein isoforms, root architecture, sugar sensing and signaling, post-stress molecular events, heavy-metal and salinity stress, and mRNA dynamics (splicing, sequestering, and ribosome loading). Finally, we explore biotechnological strategies that have applied this molecular knowledge to develop cultivars resistant to flooding or to offer alternative uses of flooding-prone soils, like bioethanol and biomass production.

Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing

A male sterile mutant was created by ⁶⁰Co γ-rays of microspores isolated from Chinese cabbage DH line 'FT'. A candidate gene for the male sterile trait was identified as Bra010198. Male sterility is used for hybrid seed production in Chinese cabbage. In this study, we derived a male sterile mutant (ftms) from Chinese cabbage DH line 'FT' by irradiating microspores with ⁶⁰Co γ-rays and realized the rapid trait transformation from male fertility to sterility for creating valuable breeding materials. Genetic analysis indicated that the male sterile trait is controlled by a single recessive nuclear gene, ftms. Microspore development in mutant ftms was aborted at the tetrad stage and associated with severely retarded degeneration and vacuolation of tapetum. Using BSR-seq analysis, the candidate region for ftms was mapped on chromosome A05. A large F₂ population was created, and the region was narrowed to approximately 1.7-Mb between markers Indel20 and Indel14 via linkage analysis. The recombination frequency was extremely suppressed because the region was located on the chromosome A05 centromere. Whole-genome resequencing of mutant ftms and wild-type 'FT' aligned only one nonsynonymous SNP to Bra010198; this gene is a homolog of Arabidopsis KNS4/UPEX1, which encodes a putative β-(1,3)-galactosyltransferase that controls pollen exine development. Comparative sequencing verified the SNP position on the fifth exon of Bra010198 in mutant ftms. Further genotyping revealed that the male sterile phenotype was fully co-segregated with this SNP. Quantitative real-time PCR indicated that Bra0101918 specifically expressed in stamen. The data presented herein suggested that Bra010198 is a strong candidate gene for ftms. Hence, we developed a male sterile line for potential application in breeding and expanded the knowledge about the molecular mechanism underlying male sterility in Chinese cabbage.

Analysis of genome-wide variants through bulked segregant RNA sequencing reveals a major gene for resistance to Plasmodiophora brassicae in Brassica oleracea

Two cabbage (Brassica oleracea) cultivars 'Tekila' and 'Kilaherb' were identified as resistant to several pathotypes of Plasmodiophora brassicae. In this study, we identified a clubroot resistance gene (Rcr7) in 'Tekila' for resistance to pathotype 3 of P. brassicae from a segregating population derived from 'Tekila' crossed with the susceptible line T010000DH3. Genetic mapping was performed by identifying the percentage of polymorphic variants (PPV), a new method proposed in this study, through bulked segregant RNA sequencing. Chromosome C7 carried the highest PPV (42%) compared to the 30-34% in the remaining chromosomes. A peak with PPV (56-73%) was found within the physical interval 41-44 Mb, which indicated that Rcr7 might be located in this region. Kompetitive Allele-Specific PCR was used to confirm the association of Rcr7 with SNPs in the region. Rcr7 was flanked by two SNP markers and co-segregated with three SNP markers in the segregating population of 465 plants. Seven genes encoding TIR-NBS-LRR disease resistance proteins were identified in the target region, but only two genes, Bo7g108760 and Bo7g109000, were expressed. Resistance to pathotype 5X was also mapped to the same region as Rcr7. B. oleracea lines including 'Kilaherb' were tested with five SNP markers for Rcr7 and for resistance to pathotype 3; 11 of 25 lines were resistant, but 'Kilaherb' was the only line that carried the SNP alleles associated with Rcr7. The presence of Rcr7 in 'Kilaherb' for resistance to both pathotypes 3 and 5X was confirmed through linkage analysis.

RNA Sequencing-Based Bulked Segregant Analysis Facilitates Efficient D-genome Marker Development for a Specific Chromosomal Region of Synthetic Hexaploid Wheat

Common wheat originated from interspecific hybridization between cultivated tetraploid wheat and its wild diploid relative Aegilops tauschii followed by amphidiploidization. This evolutionary process can be reproduced artificially, resulting in synthetic hexaploid wheat lines. Here we performed RNA sequencing (RNA-seq)-based bulked segregant analysis (BSA) using a bi-parental mapping population of two synthetic hexaploid wheat lines that shared identical A and B genomes but included with D-genomes of distinct origins. This analysis permitted identification of D-genome-specific polymorphisms around the Net2 gene, a causative locus to hybrid necrosis. The resulting single nucleotide polymorphisms (SNPs) were classified into homoeologous polymorphisms and D-genome allelic variations, based on the RNA-seq results of a parental tetraploid and two Ae. tauschii accessions. The difference in allele frequency at the D-genome-specific SNP sites between the contrasting bulks (ΔSNP-index) was higher on the target chromosome than on the other chromosomes. Several SNPs with the highest ΔSNP-indices were converted into molecular markers and assigned to the Net2 chromosomal region. These results indicated that RNA-seq-based BSA can be applied efficiently to a synthetic hexaploid wheat population to permit molecular marker development in a specific chromosomal region of the D genome.

Improving Flooding Tolerance of Crop Plants

A major problem of climate change is the increasing duration and frequency of heavy rainfall events. This leads to soil flooding that negatively affects plant growth, eventually leading to death of plants if the flooding persists for several days. Most crop plants are very sensitive to flooding, and dramatic yield losses occur due to flooding each year. This review summarizes recent progress and approaches to enhance crop resistance to flooding. Most experiments have been done on maize, barley, and soybean. Work on other crops such as wheat and rape has only started. The most promising traits that might enhance crop flooding tolerance are anatomical adaptations such as aerenchyma formation, the formation of a barrier against radial oxygen loss, and the growth of adventitious roots. Metabolic adaptations might be able to improve waterlogging tolerance as well, but more studies are needed in this direction. Reasonable approaches for future studies are quantitative trait locus (QTL) analyses or genome-wide association (GWA) studies in combination with specific tolerance traits that can be easily assessed. The usage of flooding-tolerant relatives or ancestral cultivars of the crop of interest in these experiments might enhance the chances of finding useful tolerance traits to be used in breeding.