MapChart: software for the graphical presentation of linkage maps and QTLs

Genome-wide identification and functional analysis of the cellulose synthase-like gene superfamily in common oat (Avena sativa L.)

Hemicelluloses constitute approximately one-third of the plant cell wall and can be used as a dietary fiber and food additive, and as raw materials for biofuels. Although genes involved in hemicelluloses synthesis have been investigated in some model plants, no comprehensive analysis has been conducted in common oat at present. In this study, we identified and systematically analyzed the cellulose synthase-like gene (Csl) family members in common oat and investigated them using various bioinformatics tools. The results showed that there are 76 members of the oat Csl gene family distributed on 17 chromosomes, and phylogenetic analysis indicated that the 76 Csl genes belong to the CslA, CslC, CslD, CslE, CslF, CslH, and CslJ subfamilies. A total of 14 classes of cis-acting elements were identified in the promoter regions, including hormone response, light response, cell development, and defense stress elements. The collinearity analysis identified 28 pairs of segmentally duplicated genes, most of which were found on chromosomes 2D and 6A. Expression pattern analysis showed that oat Csl genes display strong tissue-specific expression; of the 76 Csl genes, 33 were significantly up-regulated in stems and 30 were up-regulated in immature seeds. The expression of most members of the AsCsl gene family is repressed by abiotic stress, while the expression of some members is up-regulated by light. Immunoelectron microscopy shows that the product of AsCsl61, a member of CslF subfamily, mediates (1,3; 1,4)-β-D-glucan synthesis in transgenic Arabidopsis. These findings provide a fundamental understanding of the structural, functional, and evolutionary features of the oat Csl genes and may contribute to our general understanding of hemicellulose biosynthesis. Moreover, this information will be helpful in designing experiments for genetic manipulation of mixed-linkage glucan (MLG) synthesis with the goal of quality improvement in oat.

Genotyping by sequencing for the construction of oil palm (Elaeis guineensis Jacq.) genetic linkage map and mapping of yield related quantitative trait loci

Background

Oil palm (Elaeis guineensis Jacq.) is one of the major oil-producing crops. Improving the quality and increasing the production yield of oil palm have been the primary focuses of both conventional and modern breeding approaches. However, the conventional breeding approach for oil palm is very challenging due to its longevity, which results in a long breeding cycle. Thus, the establishment of marker assisted selection (MAS) for oil palm breeding programs would speed up the breeding pipeline by generating new oil palm varieties that possess high commercial traits. With the decreasing cost of sequencing, Genotyping-by-sequencing (GBS) is currently feasible to many researchers and it provides a platform to accelerate the discovery of single nucleotide polymorphism (SNP) as well as insertion and deletion (InDel) markers for the construction of a genetic linkage map. A genetic linkage map facilitates the identification of significant DNA regions associated with the trait of interest via quantitative trait loci (QTL) analysis.

Methods

A mapping population of 112 F1 individuals from a cross of Deli dura and Serdang pisifera was used in this study. GBS libraries were constructed using the double digestion method with HindIII and TaqI enzymes. Reduced representation libraries (RRL) of 112 F1 progeny and their parents were sequenced and the reads were mapped against the E. guineensis reference genome. To construct the oil palm genetic linkage map, informative SNP and InDel markers were used to discover significant DNA regions associated with the traits of interest. The nine traits of interest in this study were fresh fruit bunch (FFB) yield, oil yield (OY), oil to bunch ratio (O/B), oil to dry mesocarp ratio (O/DM) ratio, oil to wet mesocarp ratio (O/WM), mesocarp to fruit ratio (M/F), kernel to fruit ratio (K/F), shell to fruit ratio (S/F), and fruit to bunch ratio (F/B).

Results

A total of 2.5 million SNP and 153,547 InDel markers were identified. However, only a subset of 5,278 markers comprising of 4,838 SNPs and 440 InDels were informative for the construction of a genetic linkage map. Sixteen linkage groups were produced, spanning 2,737.6 cM for the maternal map and 4,571.6 cM for the paternal map, with average marker densities of one marker per 2.9 cM and one per 2.0 cM respectively, were produced. A QTL analysis was performed on nine traits; however, only QTL regions linked to M/F, K/F and S/F were declared to be significant. Of those QTLs were detected: two for M/F, four for K/F and one for S/F. These QTLs explained 18.1-25.6% of the phenotypic variance and were located near putative genes, such as casein kinase II and the zinc finger CCCH domain, which are involved in seed germination and growth. The identified QTL regions for M/F, K/F and S/F from this study could be applied in an oil palm breeding program and used to screen palms with desired traits via marker assisted selection (MAS).

Fine mapping of a new common bean anthracnose resistance gene (Co-18) to the proximal end of Pv10 in Indian landrace KRC-5

KEY MESSAGE

Mapping and fine mapping of bean anthracnose resistance genes is a continuous process. We report fine mapping of anthracnose resistance gene Co-18 which is the first anthracnose gene mapped to Pv10. The discovery of resistance gene is a major gain in the bean anthracnose pathosystem research. Among the Indian common bean landraces, KRC-5 exhibit high levels of resistance to the bean anthracnose pathogen Colletotrichum lindemuthianum. To precisely map the anthracnose resistance gene, we used a Recombinant Inbred Line (F2:9 RIL) population (KRC-5 × Jawala). The inheritance test revealed that KRC-5 carries a dominant resistance gene temporarily designated as Co-18. We discovered two RAPD markers linked to Co-18 among 287 RAPD markers. These RAPD markers were eventually developed into SCARs (Sc-OPR15 and Sc-OPF6) and flank Co-18 on chromosome Pv10 at a distance of 5.3 and 4.2 cM, respectively. At 4.0-4.1 Mb on Pv10, we detected a SNP (single-nucleotide polymorphism) signal. We synthesized 58 SSRs and 83 InDels from a pool of 135 SSRs and 1134 InDels, respectively. Five SSRs, four InDels, and two SCARs were used to generate the high-density linkage map, which led to the identification of two SSRs (SSR24 and SSR36) that are tightly linked to Co-18. These two SSRs flank the Co-18 to 178 kb genomic region with 13 candidate genes including five NLR (nucleotide-binding and leucine-rich repeat) genes. The closely linked markers SSR24 and SSR36 will be used in cloning and pyramiding of the Co-18 gene with other R genes to develop durable resistant bean varieties.

Genome-wide identification of the expansin gene family in netted melon and their transcriptional responses to fruit peel cracking

Introduction

Fruit cracking not only affects the appearance of netted melons (Cucumis melo L. var. reticulatus Naud.) but also decreases their marketability.

Methods

Herein, to comprehensively understand the role of expansin (EXP) proteins in netted melon, bioinformatics methods were employed to discover the EXP gene family in the melon genome and analyze its characteristic features. Furthermore, transcriptomics analysis was performed to determine the expression patterns of melon EXP (CmEXP) genes in crack-tolerant and crack-susceptible netted melon varieties.

Discussion

Thirty-three CmEXP genes were identified. Chromosomal location analysis revealed that CmEXP gene distribution was uneven on 12 chromosomes. In addition, phylogenetic tree analysis revealed that CmEXP genes could be categorized into four subgroups, among which the EXPA subgroup had the most members. The same subgroup members shared similar protein motifs and gene structures. Thirteen duplicate events were identified in the 33 CmEXP genes. Collinearity analysis revealed that the CmEXP genes had 50, 50, and 44 orthologous genes with EXP genes in cucumber, watermelon, and Arabidopsis, respectively. However, only nine orthologous EXP genes were observed in rice. Promoter cis-acting element analysis demonstrated that numerous cis-acting elements in the upstream promoter region of CmEXP genes participate in plant growth, development, and environmental stress responses. Transcriptomics analysis revealed 14 differentially expressed genes (DEGs) in the non-cracked fruit peels between the crack-tolerant variety 'Xizhoumi 17' (N17) and the crack-susceptible variety 'Xizhoumi 25' (N25). Among the 14 genes, 11 were upregulated, whereas the remaining three were downregulated in N17. In the non-cracked (N25) and cracked (C25) fruit peels of 'Xizhoumi 25', 24 DEGs were identified, and 4 of them were upregulated, whereas the remaining 20 were downregulated in N25. In the two datasets, only CmEXPB1 exhibited consistently upregulated expression, indicating its importance in the fruit peel crack resistance of netted melon. Transcription factor prediction revealed 56 potential transcription factors that regulate CmEXPB1 expression.

Results

Our study findings enrich the understanding of the CmEXP gene family and present candidate genes for the molecular breeding of fruit peel crack resistance of netted melon.

Characterization of the pyruvate kinase gene family in soybean and identification of a putative salt responsive gene GmPK21

BACKGROUND

As a key regulatory enzyme in the glycolysis pathway, pyruvate kinase (PK) plays crucial roles in multiple physiological processes during plant growth and is also involved in the abiotic stress response. However, little information is known about PKs in soybean.

RESULTS

In this study, we identified 27 PK family genes against the genome of soybean cultivar Zhonghuang13. They were classified into 2 subfamilies including PKc and PKp. 22 segmental duplicated gene pairs and 1 tandem duplicated gene pair were identified and all of them experienced a strong purifying selective pressure during evolution. Furthermore, the abiotic stresses (especially salt stress) and hormone responsive cis-elements were present in the promoters of GmPK genes, suggesting their potential roles in abiotic stress tolerance. By performing the qRT-PCR, 6 GmPK genes that continuously respond to both NaCl and ABA were identified. Subsequently, GmPK21, which represented the most significant change under NaCl treatment was chosen for further study. Its encoded protein GmPK21 was localized in the cytoplasm and plasma membrane. The transgenic Arabidopsis overexpressing GmPK21 exhibited weakened salinity tolerance.

CONCLUSIONS

This study provides genomic information of soybean PK genes and a molecular basis for mining salt tolerance function of PKs in the future.

Identification of two major QTLs for pod shell thickness in peanut (Arachis hypogaea L.) using BSA-seq analysis

BACKGROUND

Pod shell thickness (PST) is an important agronomic trait of peanut because it affects the ability of shells to resist pest infestations and pathogen attacks, while also influencing the peanut shelling process. However, very few studies have explored the genetic basis of PST.

RESULTS

An F2 segregating population derived from a cross between the thick-shelled cultivar Yueyou 18 (YY18) and the thin-shelled cultivar Weihua 8 (WH8) was used to identify the quantitative trait loci (QTLs) for PST. On the basis of a bulked segregant analysis sequencing (BSA-seq), four QTLs were preliminarily mapped to chromosomes 3, 8, 13, and 18. Using the genome resequencing data of YY18 and WH8, 22 kompetitive allele-specific PCR (KASP) markers were designed for the genotyping of the F2 population. Two major QTLs (qPSTA08 and qPSTA18) were identified and finely mapped, with qPSTA08 detected on chromosome 8 (0.69-Mb physical genomic region) and qPSTA18 detected on chromosome 18 (0.15-Mb physical genomic region). Moreover, qPSTA08 and qPSTA18 explained 31.1-32.3% and 16.7-16.8% of the phenotypic variation, respectively. Fifteen genes were detected in the two candidate regions, including three genes with nonsynonymous mutations in the exon region. Two molecular markers (Tif2_A08_31713024 and Tif2_A18_7198124) that were developed for the two major QTL regions effectively distinguished between thick-shelled and thin-shelled materials. Subsequently, the two markers were validated in four F2:3 lines selected.

CONCLUSIONS

The QTLs identified and molecular markers developed in this study may lay the foundation for breeding cultivars with a shell thickness suitable for mechanized peanut shelling.

LsMybW-encoding R2R3-MYB transcription factor is responsible for a shift from black to white in lettuce seed

KEY MESSAGE

We identified LsMybW as the allele responsible for the shift in color from black to white seeds in wild ancestors of lettuce to modern cultivars. Successfully selected white seeds are a key agronomic trait for lettuce cultivation and breeding; however, the mechanism underlying the shift from black-in its wild ancestor-to white seeds remains uncertain. We aimed to identify the gene/s responsible for white seed trait in lettuce. White seeds accumulated less proanthocyanidins than black seeds, similar to the phenotype observed in Arabidopsis TT2 mutants. Genetic mapping of a candidate gene was performed with double-digest RAD sequencing using an F2 population derived from a cross between "ShinanoPower" (white) and "Escort" (black). The white seed trait was controlled by a single recessive locus (48.055-50.197 Mbp) in linkage group 7. Using five PCR-based markers and numerous cultivars, eight candidate genes were mapped in the locus. Only the LG7_v8_49.251Mbp_HinfI marker, employing a single-nucleotide mutation in the stop codon of Lsat_1_v5_gn_7_35020.1, was completely linked to seed color phenotype. In addition, the coding region sequences for other candidate genes were identical in the resequence analysis of "ShinanoPower" and "Escort." Therefore, we proposed Lsat_1_v5_gn_7_35020.1 as the candidate gene and designated it as LsMybW (Lactuca sativa Myb White seeds), an ortholog encoding the R2R3-MYB transcription factor in Arabidopsis. When we validated the role of LsMybW through genome editing, LsMybW knockout mutants harboring an early termination codon showed a change in seed color from black to white. Therefore, LsMybW was the allele responsible for the shift in seed color. The development of a robust marker for marker-assisted selection and identification of the gene responsible for white seeds have implications for future breeding technology and physiological analysis.

Identification of Putative Quantitative Trait Loci for Improved Seed Oil Quality in Peanuts

Improving seed oil quality in peanut (Arachis hypogaea) has long been an aim of breeding programs worldwide. The genetic resources to achieve this goal are limited. We used an advanced recombinant inbred line (RIL) population derived from JH5 × KX01-6 to explore quantitative trait loci (QTL) affecting peanut oil quality and their additive effects, epistatic effects, and QTL × environment interactions. Gas chromatography (GC) analysis suggested seven fatty acids components were obviously detected in both parents and analyzed in a follow-up QTL analysis. The major components, palmitic acid (C16:0), oleic acid (C18:1), and linoleic acid (C18:2), exhibited considerable phenotypic variation and fit the two major gene and minor gene mixed-inheritance model. Seventeen QTL explained 2.57-38.72% of the phenotypic variation in these major components, with LOD values of 4.12-37.56 in six environments, and thirty-five QTL explained 0.94-32.21% of the phenotypic variation, with LOD values of 5.99-150.38 in multiple environments. Sixteen of these QTL were detected in both individual and multiple environments. Among these, qFA_08_1 was a novel QTL with stable, valuable and major effect. Two other major-effect QTL, qFA_09_2 and qFA_19_3, share the same physical position as FAD2A and FAD2B, respectively. Eleven stable epistatic QTL involving nine loci explained 1.30-34.97% of the phenotypic variation, with epistatic effects ranging from 0.09 to 6.13. These QTL could be valuable for breeding varieties with improved oil quality.

Mapping yield and yield-related traits using diverse common bean germplasm

Common bean (bean) is one of the most important legume crops, and mapping genes for yield and yield-related traits is essential for its improvement. However, yield is a complex trait that is typically controlled by many loci in crop genomes. The objective of this research was to identify regions in the bean genome associated with yield and a number of yield-related traits using a collection of 121 diverse bean genotypes with different yields. The beans were evaluated in replicated trials at two locations, over two years. Significant variation among genotypes was identified for all traits analyzed in the four environments. The collection was genotyped with the BARCBean6K_3 chip (5,398 SNPs), two yield/antiyield gene-based markers, and seven markers previously associated with resistance to common bacterial blight (CBB), including a Niemann-Pick polymorphism (NPP) gene-based marker. Over 90% of the single-nucleotide polymorphisms (SNPs) were polymorphic and separated the panel into two main groups of small-seeded and large-seeded beans, reflecting their Mesoamerican and Andean origins. Thirty-nine significant marker-trait associations (MTAs) were identified between 31 SNPs and 15 analyzed traits on all 11 bean chromosomes. Some of these MTAs confirmed genome regions previously associated with the yield and yield-related traits in bean, but a number of associations were not reported previously, especially those with derived traits. Over 600 candidate genes with different functional annotations were identified for the analyzed traits in the 200-Kb region centered on significant SNPs. Fourteen SNPs were identified within the gene model sequences, and five additional SNPs significantly associated with five different traits were located at less than 0.6 Kb from the candidate genes. The work confirmed associations between two yield/antiyield gene-based markers (AYD1m and AYD2m) on chromosome Pv09 with yield and identified their association with a number of yield-related traits, including seed weight. The results also confirmed the usefulness of the NPP marker in screening for CBB resistance. Since disease resistance and yield measurements are environmentally dependent and labor-intensive, the three gene-based markers (CBB- and two yield-related) and quantitative trait loci (QTL) that were validated in this work may be useful tools for simplifying and accelerating the selection of high-yielding and CBB-resistant bean cultivars.

Genome-wide analysis of the U-box E3 ligases gene family in potato (Solanum tuberosum L.) and overexpress StPUB25 enhance drought tolerance in transgenic Arabidopsis

BACKGROUND

Plant U-box (PUB) E3 ubiquitin ligases have vital effects on various biological processes. Therefore, a comprehensive and systematic identification of the members of the U-box gene family in potato will help to understand the evolution and function of U-box E3 ubiquitin ligases in plants.

RESULTS

This work identified altogether 74 PUBs in the potato (StPUBs) and examined their gene structures, chromosomal distributions, and conserved motifs. There were seventy-four StPUB genes on ten chromosomes with diverse densities. As revealed by phylogenetic analysis on PUBs within potato, Arabidopsis, tomato (Solanum lycopersicum), cabbage (Brassica oleracea), rice (Oryza sativa), and corn (Zea mays), were clustered into eight subclasses (C1-C8). According to synteny analysis, there were 40 orthologous StPUB genes to Arabidopsis, 58 to tomato, 28 to cabbage, 7 to rice, and 8 to corn. In addition, RNA-seq data downloaded from PGSC were utilized to reveal StPUBs' abiotic stress responses and tissue-specific expression in the doubled-monoploid potato (DM). Inaddition, we performed RNA-seq on the 'Atlantic' (drought-sensitive cultivar, DS) and the 'Qingshu NO.9' (drought-tolerant cultivar, DT) in early flowering, full-blooming, along with flower-falling stages to detect genes that might be involved in response to drought stress. Finally, quantitative real-time PCR (qPCR) was carried out to analyze three candidate genes for their expression levels within 100 mM NaCl- and 10% PEG 6000 (w/v)-treated potato plantlets for a 24-h period. Furthermore, we analyzed the drought tolerance of StPUB25 transgenic plants and found that overexpression of StPUB25 significantly increased peroxidase (POD) activity, reduced ROS (reactive oxygen species) and MDA (malondialdehyde) accumulation compared with wild-type (WT) plants, and enhancing drought tolerance of the transgenic plants.

CONCLUSION

In this study, three candidate genes related to drought tolerance in potato were excavated, and the function of StPUB25 under drought stress was verified. These results should provide valuable information to understand the potato StPUB gene family and investigate the molecular mechanisms of StPUBs regulating potato drought tolerance.

Functional allele of a MATE gene selected during domestication modulates seed color in chickpea

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

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

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

Genome-wide survey, molecular evolution and expression analysis of Auxin Response Factor (ARF) gene family indicating their key role in seed number per pod in pigeonpea (C. cajan L. Millsp.)

Auxin Response Factors (ARF) are a family of transcription factors that mediate auxin signalling and regulate multiple biological processes. Their crucial role in increasing plant biomass/yield influenced this study, where a systematic analysis of ARF gene family was carried out to identify the key proteins controlling embryo/seed developmental pathways in pigeonpea. A genome-wide scan revealed the presence of 12 ARF genes in pigeonpea, distributed across the chromosomes 1, 3, 4, 8 and 11. Domain analysis of ARF proteins showed the presence of B3 DNA binding, AUX response, and IAA domains. Majority of them are of nuclear origin, and do not exhibit the level of genomic expansion as observed in Glycine max (51 members). The duplication events seem to range from 31.6 to 42.3 million years ago (mya). Promoter analysis revealed the presence of multiple cis-acting elements related to stress responses, hormone signalling and other development processes. The expression atlas data highlighted the expression of CcARF8 in hypocotyl, bud and flower whereas, CcARF7 expression was significantly high in pod. The real-time expression of CcARF2, CcARF3 and CcARF18 was highest in genotypes with high seed number indicating their key role in regulating embryo development and determining seed set in pigeonpea.

Navigating rice seedling cold resilience: QTL mapping in two inbred line populations and the search for genes

Due to global climate change resulting in extreme temperature fluctuations, it becomes increasingly necessary to explore the natural genetic variation in model crops such as rice to facilitate the breeding of climate-resilient cultivars. To uncover genomic regions in rice involved in managing cold stress tolerance responses and to identify associated cold tolerance genes, two inbred line populations developed from crosses between cold-tolerant and cold-sensitive parents were used for quantitative trait locus (QTL) mapping of two traits: degree of membrane damage after 1 week of cold exposure quantified as percent electrolyte leakage (EL) and percent low-temperature seedling survivability (LTSS) after 1 week of recovery growth. This revealed four EL QTL and 12 LTSS QTL, all overlapping with larger QTL regions previously uncovered by genome-wide association study (GWAS) mapping approaches. Within the QTL regions, 25 cold-tolerant candidate genes were identified based on genomic differences between the cold-tolerant and cold-sensitive parents. Of those genes, 20% coded for receptor-like kinases potentially involved in signal transduction of cold tolerance responses; 16% coded for transcription factors or factors potentially involved in regulating cold tolerance response effector genes; and 64% coded for protein chaperons or enzymes potentially serving as cold tolerance effector proteins. Most of the 25 genes were cold temperature regulated and had deleterious nucleotide variants in the cold-sensitive parent, which might contribute to its cold-sensitive phenotype.

Genetic control of pod morphological traits and pod edibility in a common bean RIL population

KEY MESSAGE

QTL mapping, association analysis, and colocation study with previously reported QTL revealed three main regions controlling pod morphological traits and two loci for edible pod characteristics on the common bean chromosomes Pv01 and Pv06. Bean pod phenotype is a complex characteristic defined by the combination of different traits that determine the potential use of a genotype as a snap bean. In this study, the TUM RIL population derived from a cross between 'TU' (dry) and 'Musica' (snap) was used to investigate the genetic control of pod phenotype. The character was dissected into pod morphological traits (PMTs) and edible pod characteristics (EPC). The results revealed 35 QTL for PMTs located on seven chromosomes, suggesting a strong QTL colocation on chromosomes Pv01 and Pv06. Some QTL were colocated with previously reported QTL, leading to the mapping of 15 consensus regions associated with bean PMTs. Analysis of EPC of cooked beans revealed that two major loci with epistatic effect, located on chromosomes Pv01 and Pv06, are involved in the genetic control of this trait. An association study using a subset of the Spanish Diversity Panel (snap vs. non-snap) detected 23 genomic regions, with three regions being mapped at a position similar to those of two loci identified in the TUM population. The results demonstrated the relevant roles of Pv01 and Pv06 in the modulation of bean pod phenotype. Gene ontology enrichment analysis revealed a significant overrepresentation of genes regulating the phenylpropanoid metabolic process and auxin response in regions associated with PMTs and EPC, respectively. Both biological functions converged in the lignin biosynthetic pathway, suggesting the key role of the pathway in the genetic control of bean pod phenotype.

Genome-wide family prediction unveils molecular mechanisms underlying the regulation of agronomic traits in Urochloa ruziziensis

Tropical forage grasses, particularly those belonging to the Urochloa genus, play a crucial role in cattle production and serve as the main food source for animals in tropical and subtropical regions. The majority of these species are apomictic and tetraploid, highlighting the significance of U. ruziziensis, a sexual diploid species that can be tetraploidized for use in interspecific crosses with apomictic species. As a means to support breeding programs, our study investigates the feasibility of genome-wide family prediction in U. ruziziensis families to predict agronomic traits. Fifty half-sibling families were assessed for green matter yield, dry matter yield, regrowth capacity, leaf dry matter, and stem dry matter across different clippings established in contrasting seasons with varying available water capacity. Genotyping was performed using a genotyping-by-sequencing approach based on DNA samples from family pools. In addition to conventional genomic prediction methods, machine learning and feature selection algorithms were employed to reduce the necessary number of markers for prediction and enhance predictive accuracy across phenotypes. To explore the regulation of agronomic traits, our study evaluated the significance of selected markers for prediction using a tree-based approach, potentially linking these regions to quantitative trait loci (QTLs). In a multiomic approach, genes from the species transcriptome were mapped and correlated to those markers. A gene coexpression network was modeled with gene expression estimates from a diverse set of U. ruziziensis genotypes, enabling a comprehensive investigation of molecular mechanisms associated with these regions. The heritabilities of the evaluated traits ranged from 0.44 to 0.92. A total of 28,106 filtered SNPs were used to predict phenotypic measurements, achieving a mean predictive ability of 0.762. By employing feature selection techniques, we could reduce the dimensionality of SNP datasets, revealing potential genotype-phenotype associations. The functional annotation of genes near these markers revealed associations with auxin transport and biosynthesis of lignin, flavonol, and folic acid. Further exploration with the gene coexpression network uncovered associations with DNA metabolism, stress response, and circadian rhythm. These genes and regions represent important targets for expanding our understanding of the metabolic regulation of agronomic traits and offer valuable insights applicable to species breeding. Our work represents an innovative contribution to molecular breeding techniques for tropical forages, presenting a viable marker-assisted breeding approach and identifying target regions for future molecular studies on these agronomic traits.

Genome- and transcriptome-wide identification of trehalose-6-phosphate phosphatases (TPP) gene family and their expression patterns under abiotic stress and exogenous trehalose in soybean

BACKGROUND

Trehalose-6-phosphate phosphatase (TPP) is an essential enzyme catalyzing trehalose synthesis, an important regulatory factor for plant development and stress response in higher plants. However, the TPP gene family in soybean has not been reported.

RESULTS

A comprehensive analysis of the TPP gene family identified 18 GmTPPs classified into eight groups based on the phylogenetic relationships and the conservation of protein in six monocot and eudicot plants. The closely linked subfamilies had similar motifs and intron/exon numbers. Segmental duplication was the main driving force of soybean GmTPPs expansion. In addition, analysis of the cis-regulatory elements and promoter regions of GmTPPs revealed that GmTPPs regulated the response to several abiotic stresses. Moreover, RNA-seq and qRT-PCR analysis of the tissue-specific GmTPPs under different abiotic stresses revealed that most GmTPPs were associated with response to different stresses, including cold, drought, saline-alkali, and exogenous trehalose. Notably, exogenous trehalose treatment up-regulated the expression of most TPP genes under saline-alkali conditions while increasing the carbohydrate and trehalose levels and reducing reactive oxygen species (ROS) accumulation in soybean sprouts, especially in the saline-alkali tolerant genotype. Furthermore, the interaction network and miRNA target prediction revealed that GmTPPs interacted with abiotic stress response-related transcription factors.

CONCLUSIONS

The findings in this study lay a foundation for further functional studies on TPP-based breeding to improve soybean development and stress tolerance.

Sr65: a widely effective gene for stem rust resistance in wheat

KEY MESSAGE

Sr65 in chromosome 1A of Indian wheat landrace Hango-2 is a potentially useful all-stage resistance gene that currently protects wheat from stem rust in Australia, India, Africa and Europe. Stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), threatened global wheat production with the appearance of widely virulent races that included TTKSK and TTRTF. Indian landrace Hango-2 showed resistance to Pgt races in India and Australia. Screening of a Hango-2/Avocet 'S' (AvS) recombinant inbred line population identified two stem rust resistance genes, a novel gene (temporarily named as SrH2) from Hango-2 and Sr26 from AvS. A mapping population segregating for SrH2 alone was developed from two recombinant lines. SrH2 was mapped on the short arm of chromosome 1A, where it was flanked by KASP markers KASP_7944 (proximal) and KASP_12147 (distal). SrH2 was delimited to an interval of 1.8-2.3 Mb on chromosome arm 1AS. The failure to detect candidate genes through MutRenSeq and comparative genomic analysis with the pan-genome dataset indicated the necessity to generate a Hango-2 specific assembly for detecting the gene sequence linked with SrH2 resistance. MutRenSeq however enabled identification of SrH2-linked KASP marker sunCS_265. Markers KASP_12147 and sunCS_265 showed 92% and 85% polymorphism among an Australian cereal cultivar diversity panel and can be used for marker-assisted selection of SrH2 in breeding programs. The effectiveness of SrH2 against Pgt races from Europe, Africa, India, and Australia makes it a valuable resource for breeding stem rust-resistant wheat cultivars. Since no wheat-derived gene was previously located in chromosome arm 1AS, SrH2 represents a new locus and named as SR65.

Genome-wide identification and expression analysis of WRKY gene family members in red clover (Trifolium pratense L.)

Trifolium pratense is an important legume forage grass and a key component of sustainable livestock development. Serving as an essential component, the WRKY gene family, a crucial group of regulatory transcription factors in plants, holds significant importance in their response to abiotic stresses. However, there has been no systematic analysis conducted on the WRKY gene family in Trifolium pratense. This study conducted a comprehensive genomic characterization of the WRKY gene family in Trifolium pratense, utilizing the latest genomic data, resulting in the identification of 59 TpWRKY genes. Based on their structural features, phylogenetic characteristics, and conserved motif composition, the WRKY proteins were classified into three groups, with group II further subdivided into five subgroups (II-a, II-b, II-c, II-d, and II-e). The majority of the TpWRKYs in a group share a similar structure and motif composition. Intra-group syntenic analysis revealed eight pairs of duplicate segments. The expression patterns of 59 TpWRKY genes in roots, stems, leaves, and flowers were examined by analyzing RNA-seq data. The expression of 12 TpWRKY genes under drought, low-temperature (4°C), methyl jasmonate (MeJA) and abscisic acid (ABA) stresses was analyzed by RT-qPCR. The findings indicated that TpWRKY46 was highly induced by drought stress, and TpWRKY26 and TpWRKY41 were significantly induced by low temperature stress. In addition, TpWRKY29 and TpWRKY36 were greatly induced by MeJA stress treatment, and TpWRKY17 was significantly upregulated by ABA stress treatment. In this research, we identified and comprehensively analyzed the structural features of the WRKY gene family in T.pratense, along with determined the possible roles of WRKY candidate genes in abiotic stress. These discoveries deepen our understandings of how WRKY transcription factors contribute to species evolution and functional divergence, laying a solid molecular foundation for future exploration and study of stress resistance mechanisms in T.pratense.

Genome-wide analysis of UDP-glycosyltransferase gene family and identification of a flavonoid 7-O-UGT (AhUGT75A) enhancing abiotic stress in peanut (Arachis hypogaea L.)

BACKGROUND

Glycosylation, catalyzed by UDP-glycosyltransferase (UGT), was important for enhancing solubility, bioactivity, and diversity of flavonoids. Peanut (Arachis hypogaea L.) is an important oilseed and cash crop worldwide. In addition to provide high quality of edible oils and proteins, peanut seeds contain a rich source of flavonoid glycosides that benefit human health. However, information of UGT gene family was quite limited in peanut.

RESULTS

In present study, a total of 267 AhUGTs clustered into 15 phylogenetic groups were identified in peanut genome. Group I has greatly expanded to contain the largest number of AhUGT genes. Segmental duplication was the major driving force for AhUGT gene family expansion. Transcriptomic analysis of gene expression profiles in various tissues and under different abiotic stress treatments indicated AhUGTs were involved in peanut growth and abiotic stress response. AhUGT75A (UGT73CG33), located in mitochondria, was characterized as a flavonoid 7-O-UGT by in vitro enzyme assays. The transcript level of AhUGT75A was strongly induced by abiotic stress. Overexpression of AhUGT75A resulted in accumulating less amount of malondialdehyde (MDA) and superoxide, and enhancing tolerance against drought and/or salt stress in transgenic Arabidopsis. These results indicated AhUGT75A played important roles in conferring abiotic stress tolerance through reactive oxygen species scavenging.

CONCLUSIONS

Our research only not provides valuable information for functional characterization of UGTs in peanut, but also gives new insights into potential applications in breeding new cultivars with both desirable stress tolerance and health benefits.

Genome-wide identification of RNA recognition motif (RRM1) in Brassica rapa and functional analysis of RNA-binding protein (BrRBP) under low-temperature stress

BACKGROUND

The RNA recognition motif (RRM) is primarily engaged in the processing of mRNA and rRNA following gene transcription as well as the regulation of RNA transport; it is critical in preserving RNA stability.

RESULTS

In this study, we identified 102 members of the RRM1 gene family in Brassica rapa, which were dispersed across 10 chromosomes with the ninth chromosome being the most extensively distributed. The RRM1 gene family members of Brassica rapa and Arabidopsis thaliana were grouped into 14 subclades (I-XIV) using phylogenetic analysis. Moreover, the results of transcriptome analysis and RT-qPCR indicated that the expression of Brapa05T000840 was upregulated in the cultivars 'Longyou 7' and 'Longyou 99' following exposure to cold stress at a temperature of 4 °C for 24 h. The levels of expression in the leaves and growth cones of the 'Longyou 7' variety were found to be significantly higher than those observed in the 'Longyou 99' variety under conditions of low temperature and NaCl stress. It illustrates the involvement of the RRM1 gene in the physiological response to both low temperature and salt stress. In addition, it was observed that the survival rate of transgenic BrRBP (Brapa05T000840) Arabidopsis thaliana plants was notably higher compared to that of wild-type plants when subjected to varying durations of low temperature treatment. Furthermore, the expression of the BrRBP gene in transgenic plants exhibited an upward trend as the duration of low temperature treatment increased, reaching its peak at 24 h. The in-vivo enzymatic activity of reactive oxygen species-scavenging enzymes were found to be significantly elevated in comparison to wild-type plants, suggesting that the BrRBP gene may enhance the cold tolerance of Arabidopsis thaliana.

CONCLUSIONS

This study offers a significant foundation for comprehending the regulation mechanism of the RRM1 gene family in winter Brassica rapa subjected to cold stress, as well as for finding key genes associated with cold resistance.

Genome-wide identification and expression analysis of the NHX gene family under salt stress in wheat (Triticum aestivum L)

Salt stress affects plant growth and development, resulting in the loss of crop yield across the world, and sodium-proton antiporters (NHXs) are one of the genes known to promote salt tolerance in transgenic plants. In this study, we conducted a comprehensive genome-wide analysis and expression profile of NHX genes in wheat under salinity stress. We identified 30 TaNHX genes in wheat based on the Na+/H+ exchanger domain, with all genes containing an amiloride motif except one, a known for inhibiting Na+ ions in plants. Phylogenetic analysis classified these genes into three classes with subfamilies: 12 were localized in vacuoles, while 18 were in the endoplasmic reticulum and plasma membrane. Promoter analysis revealed stress-related cis-acting elements, indicating their potential role in abiotic stress tolerance. The non-synonymous (Ka)/synonymous (Ks) ratios highlighted that the majority of TaNHX genes experienced robust purifying selection throughout their evolutionary history. Transcriptomis data analysis and qRT-PCR demonstrated distinct expression patterns for TaNHX genes across various tissues when subjected to salt stress. Additionally, we predicted 20 different miRNA candidates targeting the identified TaNHX genes. Protein-protein interaction prediction revealed NHX6's involvement in the SOS1 pathway, while NHX1 gene exhibit proton antiporter activity. Molecular dynamics (MD) simulations were also conducted to examine the interactions of TaNHX1, TaNHX2, and TaNHX3. These results represent a significant advancement in our understanding of the molecular mechanisms governing Na+ transporters. This may also offer promising avenues for future studies aimed at unraveling the intricate details of their biological roles and applications.

Genome wide identification and characterization of the amino acid transporter (AAT) genes regulating seed protein content in chickpea (Cicer arietinum L.)

Amino acid transporters (AATs), besides, being a crucial component for nutrient partitioning system are also vital for growth and development of the plants and stress resilience. In order to understand the role of AAT genes in seed quality proteins, a comprehensive analysis of AAT gene family was carried out in chickpea leading to identification of 109 AAT genes, representing 10 subfamilies with random distribution across the chickpea genome. Several important stress responsive cis-regulatory elements like Myb, ABRE, ERE were detected in the promoter region of these CaAAT genes. Most of the genes belonging to the same sub-families shared the intron-exon distribution pattern owing to their conserved nature. Random distribution of these CaAAT genes was observed on plasma membrane, vacuolar membrane, Endoplasmic reticulum and Golgi membranes, which may be associated to distinct biochemical pathways. In total 92 out 109 CaAAT genes arise as result of duplication, among which segmental duplication was more prominent over tandem duplication. As expected, the phylogenetic tree was divided into 2 major clades, and further sub-divided into different sub-families. Among the 109 CaAAT genes, 25 were found to be interacting with 25 miRNAs, many miRNAs like miR156, miR159 and miR164 were interacting only with single AAT genes. Tissues specific expression pattern of many CaAAT genes was observed like CaAAP7 and CaAVT18 in nodules, CaAAP17, CaAVT5 and CaCAT9 in vegetative tissues while CaCAT10 and CaAAP23 in seed related tissues as per the expression analysis. Mature seed transcriptome data revealed that genotypes having high protein content (ICC 8397, ICC 13461) showed low CaAATs expression as compared to the genotypes having low protein content (FG 212, BG 3054). Amino acid profiling of these genotypes revealed a significant difference in amount of essential and non-essential amino acids, probably due to differential expression of CaAATs. Thus, the present study provides insights into the biological role of AAT genes in chickpea, which will facilitate their functional characterization and role in various developmental stages, stress responses and involvement in nutritional quality enhancement.

Mapping QTL for vernalization requirement identified adaptive divergence of the candidate gene Flowering Locus C in polyploid Camelina sativa

Vernalization requirement is an integral component of flowering in winter-type plants. The availability of winter ecotypes among Camelina species facilitated the mapping of quantitative trait loci (QTL) for vernalization requirement in Camelina sativa. An inter and intraspecific crossing scheme between related Camelina species, where one spring and two different sources of winter-type habit were used, resulted in the development of two segregating populations. Linkage maps generated with sequence-based markers identified three QTLs associated with vernalization requirement in C. sativa; two from the interspecific (chromosomes 13 and 20) and one from the intraspecific cross (chromosome 8). Notably, the three loci were mapped to different homologous regions of the hexaploid C. sativa genome. All three QTLs were found in proximity to Flowering Locus C (FLC), variants of which have been reported to affect the vernalization requirement in plants. Temporal transcriptome analysis for winter-type Camelina alyssum demonstrated reduction in expression of FLC on chromosomes 13 and 20 during cold treatment, which would trigger flowering, since FLC would be expected to suppress floral initiation. FLC on chromosome 8 also showed reduced expression in the C. sativa ssp. pilosa winter parent upon cold treatment, but was expressed at very high levels across all time points in the spring-type C. sativa. The chromosome 8 copy carried a deletion in the spring-type line, which could impact its functionality. Contrary to previous reports, all three FLC loci can contribute to controlling the vernalization response in C. sativa and provide opportunities for manipulating this requirement in the crop.

Quantitative Trait Locus Mapping for Fire Blight Resistance in an F2 Population of Malus fusca MAL0045 Uncovers Novel Resistance Loci

Several fire blight resistance loci in Malus genotypes map on different linkage groups (LGs) representing chromosomes of the domesticated apple. Prior genetics studies primarily focused on F1 populations. A strong resistance quantitative trait locus (QTL) explained up to 66% of phenotypic variance in an F1 progeny derived from crossing the highly resistant wild apple genotype Malus fusca MAL0045 and the highly susceptible apple cultivar 'Idared', which was previously mapped on LG10 (Mfu10) of MAL0045. Strains of the causative bacterial pathogen Erwinia amylovora, notably those that show a single nucleotide polymorphism in the avrRpt2EA effector protein sequence at position 156 (e.g., Ea3049), are more virulent and overcome some known fire blight resistance donors and their QTLs. However, MAL0045 is resistant to Ea3049 and Mfu10 is not overcome, but most of the F1 progeny were highly susceptible to this strain. This phenomenon led to the assumption that other putative resistance factors not segregating in the F1 progeny might be present in the genome of MAL0045. Here, we crossed F1 progeny together to obtain 135 F2 individuals. Facilitated by genotyping-by-sequencing and phenotypic assessments, we identified and mapped two novel resistance QTLs in these F2 individuals on LGs 4 and 15, which were not identified in the F1. To our knowledge, these are the first resistance QTLs mapped in F2 progeny in Malus. In addition, we report that neither MAL0045 nor Mfu10 is broken down by a highly aggressive U.S. strain, LA635, after analyses in the original F1 individuals. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Natural alleles of Mediator subunit genes modulate plant height in chickpea

Plant height (PH) is an important plant architectural trait targeted during Green Revolution to enhance crop yields. Identification of genes and natural alleles governing plant height without compromising agronomic performance can fill the lacuna of knowledge connecting ideal plant architecture with maximum achievable yield in chickpea. Through coherent strategy involving genome‐wide association study, QTL/fine mapping, map‐based cloning, molecular haplotyping, and downstream functional genomics, the current study identified two Mediator subunit genes namely, CaMED23 and CaMED5b and their derived natural alleles/haplotypes underlying the major QTLs and trans‐acting eQTLs regulating plant height in chickpea. Differential accumulation of haplotype‐specific transcripts of these two Mediator genes in corresponding haplotype‐introgressed near‐isogenic lines (NILs) correlates negatively with the plant height trait. Quantitative as well as qualitative estimation based on histology, scanning electron microscopy, and histochemical assay unraveled the reduced lengths and cell sizes of internodes along with compromised lignin levels in dwarf/semi‐dwarf chickpea NILs introgressed with superior CaMED23 and CaMED5b gene haplotypes. This observation, supported by global transcriptome profiling‐based diminished expression of various phenylpropanoid pathway genes upstream of lignin biosynthesis in dwarf/semi‐dwarf NILs, essentially links plant height with lignin accumulation. The identified molecular signatures in the Mediator subunit genes can be efficiently utilized to develop desirable dwarf/semi‐dwarf‐type chickpea cultivars without affecting their yield per plant via modulating lignin/phenylpropanoid biosynthesis.

Berry texture QTL and candidate gene analysis in grape (Vitis vinifera L.)

Berry texture is a noteworthy economic trait for grape; however, the genetic bases and the complex gene expression and regulatory mechanism for the diverse changes in berry texture are still poorly understood. In this study, the results suggest that it is difficult to obtain high-mesocarp firmness (MesF) and high-pericarp puncture hardness (PPH) grape cultivars with high pericarp brittleness (PerB). The high-density linkage map was constructed using whole-genome resequencing based on 151 F1 individuals originating from intraspecific hybridization between the firm-flesh cultivar 'Red Globe' and soft-flesh cultivar 'Muscat Hamburg'. The total length of the consensus map was 1613.17 cM, with a mean genetic distance between adjacent bin markers of 0.59 cM. Twenty-seven quantitative trait loci (QTLs) for berry MesF, PPH, and PerB were identified in linkage groups (LGs) 1, 3, 4, 6, 8, 9, 10, 11, 14, 16, and 17, including twelve QTLs that were firstly detected in LGs 6, 11, and 14. Fourteen promising candidate genes were identified from the stable QTL regions in LGs 10, 11, 14, and 17. In particular, VvWARK2 and VvWARK8 refer to chromosome 17 and are two promising candidate genes for MesF and PPH, as the VvWARK8 gene may increase pectin residue binding with WARK for high berry firmness maintenance and the allele for VvWARK2 carrying the 'CC' and 'GA' genotypes at Chr17:1836764 and Chr17:1836770 may be associated with non-hard texture grape cultivars. In addition, real-time quantitative polymerase chain reaction (RT-qPCR) verification revealed that the promising candidate transcription factor genes VvMYB4-like, VvERF113, VvWRKY31, VvWRKY1, and VvNAC83 may regulate cell wall metabolism candidate gene expression for grape berry texture changes.

In silico identification, characterization, and expression analysis of RNA recognition motif (RRM) containing RNA-binding proteins in Aedes aegypti

RNA-binding proteins (RBPs) are the proteins that bind RNAs and regulate their functioning. RBPs in mosquitoes are gaining attention due to their ability to bind flaviviruses and regulate their replication and transmission. Despite their relevance, RBPs in mosquitoes are not explored much. In this study, we screened the whole genome of Aedes aegypti, the primary vector of several pathogenic viruses, and identified the proteins containing RNA recognition motif (RRM), the most abundant protein domain in eukaryotes. Using several in silico strategies, a total of 135 RRM-containing RBPs were identified in Ae. aegypti. The proteins were characterized based on their available annotations and the sequence similarity with Drosophila melanogaster. Ae. aegypti RRM-containing RBPs included serine/arginine-rich (SR) proteins, polyadenylate-binding proteins (PABP), heteronuclear ribonucleoproteins (hnRNP), small nuclear ribonucleoproteins (snRNP), splicing factors, eukaryotic initiation factors, transformers, and nucleolysins. Phylogenetic analysis revealed that the proteins and the domain organization are conserved among Ae. aegypti, Bombyx mori, and Drosophila melanogaster. However, the gene length and the intron-exon organization varied across the insect species. Expression analysis of the genes encoding RBPs using publicly available RNA sequencing data for different developmental time points of the mosquito life cycle starting from the ovary and eggs up to the adults revealed stage-specific expression with several genes preferentially expressed in early embryonic stages and blood-fed female ovaries. This is the first database for the Ae. aegypti RBPs that can serve as the reference base for future investigations. Stage-specific genes can be further explored to determine their role in mosquito growth and development with a focus on developing novel mosquito control strategies.

Genome-wide identification of the C2H2-Zinc finger gene family and functional validation of CsZFP7 in citrus nucellar embryogenesis

KEY MESSAGE

Genome-wide identification of C2H2-ZF gene family in the poly- and mono-embryonic citrus species and validation of the positive role of CsZFP7 in sporophytic apomixis. The C2H2 zinc finger (C2H2-ZF) gene family is involved in plant vegetative and reproductive development. Although a large number of C2H2 zinc-finger proteins (C2H2-ZFPs) have been well characterized in some horticultural plants, little is known about the C2H2-ZFPs and their function in citrus. In this work, we performed a genome-wide sequence analysis and identified 97 and 101 putative C2H2-ZF gene family members in the genomes of sweet orange (C. sinensis, poly-embryonic) and pummelo (C. grandis, mono-embryonic), respectively. Phylogenetic analysis categorized citrus C2H2-ZF gene family into four clades, and their possible functions were inferred. According to the numerous regulatory elements on promoter, citrus C2H2-ZFPs can be divided into five different regulatory function types that indicate functional differentiation. RNA-seq data revealed 20 differentially expressed C2H2-ZF genes between poly-embryonic and mono-embryonic ovules at two stages of citrus nucellar embryogenesis, among them CsZFP52 specifically expressed in mono-embryonic pummelo ovules, while CsZFP7, 37, 44, 45, 67 and 68 specifically expressed in poly-embryonic sweet orange ovules. RT-qPCR further validated that CsZFP7 specifically expressed at higher levels in poly-embryonic ovules, and down-regulation of CsZFP7 in the poly-embryonic mini citrus (Fortunella hindsii) increased rate of mono-embryonic seeds compared with the wild type, indicating the regulatory potential of CsZFP7 in nucellar embryogenesis of citrus. This work provided a comprehensive analysis of C2H2-ZF gene family in citrus, including genome organization and gene structure, phylogenetic relationships, gene duplications, possible cis-elements on promoter regions and expression profiles, especially in the poly- and mono-embryogenic ovules, and suggested that CsZFP7 is involved in nucellar embryogenesis.

Two haplotype-resolved genomes reveal important flower traits in bigleaf hydrangea (Hydrangea macrophylla) and insights into Asterid evolution

The Hydrangea genus belongs to the Hydrangeaceae family, in the Cornales order of flowering plants, which early diverged among the Asterids, and includes several species that are commonly used ornamental plants. Of them, Hydrangea macrophylla is one of the most valuable species in the nursery trade, yet few genomic resources are available for this crop or closely related Asterid species. Two high-quality haplotype-resolved reference genomes of hydrangea cultivars 'Veitchii' and 'Endless Summer' [highest quality at 2.22 gigabase pairs (Gb), 396 contigs, N50 22.8 megabase pairs (Mb)] were assembled and scaffolded into the expected 18 pseudochromosomes. Utilizing the newly developed high-quality reference genomes along with high-quality genomes of other related flowering plants, nuclear data were found to support a single divergence point in the Asterids clade where both the Cornales and Ericales diverged from the euasterids. Genetic mapping with an F1 hybrid population demonstrated the power of linkage mapping combined with the new genomic resources to identify the gene for inflorescence shape, CYP78A5 located on chromosome 4, and a novel gene, BAM3 located on chromosome 17, for causing double flower. Resources developed in this study will not only help to accelerate hydrangea genetic improvement but also contribute to understanding the largest group of flowering plants, the Asterids.

Wheat stripe rust resistance locus YR63 is a hot spot for evolution of defence genes - a pangenome discovery

BACKGROUND

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), poses a threat to global wheat production. Deployment of widely effective resistance genes underpins management of this ongoing threat. This study focused on the mapping of stripe rust resistance gene YR63 from a Portuguese hexaploid wheat landrace AUS27955 of the Watkins Collection.

RESULTS

YR63 exhibits resistance to a broad spectrum of Pst races from Australia, Africa, Asia, Europe, Middle East and South America. It was mapped to the short arm of chromosome 7B, between two single nucleotide polymorphic (SNP) markers sunCS_YR63 and sunCS_67, positioned at 0.8 and 3.7 Mb, respectively, in the Chinese Spring genome assembly v2.1. We characterised YR63 locus using an integrated approach engaging targeted genotyping-by-sequencing (tGBS), mutagenesis, resistance gene enrichment and sequencing (MutRenSeq), RNA sequencing (RNASeq) and comparative genomic analysis with tetraploid (Zavitan and Svevo) and hexaploid (Chinese Spring) wheat genome references and 10+ hexaploid wheat genomes. YR63 is positioned at a hot spot enriched with multiple nucleotide-binding and leucine rich repeat (NLR) and kinase domain encoding genes, known widely for defence against pests and diseases in plants and animals. Detection of YR63 within these gene clusters is not possible through short-read sequencing due to high homology between members. However, using the sequence of a NLR member we were successful in detecting a closely linked SNP marker for YR63 and validated on a panel of Australian bread wheat, durum and triticale cultivars.

CONCLUSIONS

This study highlights YR63 as a valuable source for resistance against Pst in Australia and elsewhere. The closely linked SNP marker will facilitate rapid introgression of YR63 into elite cultivars through marker-assisted selection. The bottleneck of this study reinforces the necessity for a long-read sequencing such as PacBio or Oxford Nanopore based techniques for accurate detection of the underlying resistance gene when it is part of a large gene cluster.

Genome-Wide Identification and Expression Pattern of MYB Family Transcription Factors in Erianthus fulvus

MYB family genes have many functions and are widely involved in plant abiotic-stress responses. Erianthus fulvus is an important donor material for stress-resistance genes in sugarcane breeding. However, the MYB family genes in E. fulvus have not been systematically investigated. In this study, 133 EfMYB genes, including 48 Ef1R-MYB, 84 EfR2R3-MYB and 1 Ef3R-MYB genes, were identified in the E. fulvus genome. Among them, the EfR2R3-MYB genes were classified into 20 subgroups. In addition, these EfMYB genes were unevenly distributed across 10 chromosomes. A total of 4 pairs of tandemly duplicated EfMYB genes and 21 pairs of segmentally duplicated EfMYB genes were identified in the E. fulvus genome. Protein-interaction analysis predicted that 24 EfMYB proteins had potential interactions with 14 other family proteins. The EfMYB promoter mainly contains cis-acting elements related to the hormone response, stress response, and light response. Expression analysis showed that EfMYB39, EfMYB84, and EfMYB124 could be significantly induced using low-temperature stress. EfMYB30, EfMYB70, EfMYB81, and EfMYB101 responded positively to drought stress. ABA treatment significantly induced EfMYB1, EfMYB30, EfMYB39, EfMYB84, and EfMYB130. All nine genes were induced using MeJA treatment. These results provide comprehensive information on EfMYB genes and can serve as a reference for further studies of gene function.

Identification of Genomic Regions Associated with Seedling Frost Tolerance in Sorghum

Sorghum bicolor (L.) Moench is the fifth most valuable cereal crop globally. Although sorghum is tolerant to drought and elevated temperatures, it is susceptible to chilling, frost, and freezing stresses. Sorghum seeds planted in April may encounter frequent frost during late April and early May. Early spring freezing temperatures adversely affect crop development and yield. This study aims to identify genomic regions associated with frost tolerance at the seedlings stage. Breeding freeze-tolerant cultivars require selection for freeze tolerance in nurseries. However, the unpredictability of environmental conditions complicates the identification of freeze-tolerant genotypes. An indoor selection protocol has been developed to investigate the genetic determinism of freeze tolerance at the seedling stages and its correlation with several developmental traits. To accomplish this, we used two populations of recombinant inbred lines (RIL) developed from crosses between cold-tolerant (CT19, ICSV700) and cold-sensitive (TX430, M81E) parents. The derived RIL populations were evaluated for single nucleotide polymorphism (SNP) using genotype-by-sequencing (GBS) under controlled environments for their response to freezing stress. Linkage maps were constructed with 464 and 875 SNPs for the CT19 X TX430 (C1) and ICSV700 X M81E(C2) populations. Using quantitative trait loci (QTL) mapping, we identified six QTLs conferring tolerance to freezing temperatures. One QTL in the C1 population and four QTLs in the C2 population, explain 17.75-98% of the phenotypic variance of traits measured. Proline leaf content was increased in response to exposing the seedlings to low temperatures. Candidate QTLs identified in this study could be further exploited to develop frost-tolerant cultivars as proxies in marker-assisted breeding, genomic selection, and genetic engineering.

Glyceraldehyde-3-phosphate dehydrogenase Gh_GAPDH9 is associated with drought resistance in Gossypium hirsutum

Background

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is the central enzyme of glycolysis and plays important regulatory roles in plant growth and development and responses to adverse stress conditions. However, studies on the characteristics and functions of cotton GAPDH family genes are still lacking.

Methods

In this study, genome-wide identification of the cotton GAPDH gene family was performed, and the phylogeny, gene structures, promoter progenitors and expression profiles of upland cotton GAPDH gene family members were explored by bioinformatics analysis to highlight potential functions. The functions of GhGAPDH9 in response to drought stress were initially validated based on RNA-seq, qRT‒PCR, VIGS techniques and overexpression laying a foundation for further studies on the functions of GAPDH genes.

Results

This study is the first systematic analysis of the cotton GAPDH gene family, which contains a total of 84 GAPDH genes, among which upland cotton contains 27 members. Quantitative, phylogenetic and covariance analyses of the genes revealed that the GAPDH gene family has been conserved during the evolution of cotton. Promoter analysis revealed that most cis-acting elements were related to MeJA and ABA. Based on the identified promoter cis-acting elements and RNA-seq data, it was hypothesized that Gh_GAPDH9, Gh_GAPDH11, Gh_GAPDH19 and Gh_GAPDH21 are involved in the response of cotton to abiotic stress. The expression levels of the Gh_GAPDH9 gene in two drought-resistant and two drought-sensitive materials were analyzed by qRT‒PCR and found to be high early in the treatment period in the drought-resistant material. The silencing of Gh_GAPDH9 based on virus-induced gene silencing (VIGS) technology resulted in significant leaf wilting or whole-plant dieback in silenced plants after drought stress compared to the control. The content of-malondialdehyde (MDA) in cotton leaves was significantly increased, and the content of proline (Pro) and chlorophyll (Chl) was reduced. In addition, the leaf wilting and dryness of transgenic lines under drought stress were lower than those of wild-type Arabidopsis, indicating that Gh_GAPDH9 is a positive regulator of drought resistance. In conclusion, our results demonstrate that GAPDH genes play an important role in the response of cotton to abiotic stresses and provide preliminary validation of the function of the Gh_GAPDH9 gene under drought stress. These findings provide an important theoretical basis for further studies on the function of the Gh_GAPDH9 gene and the molecular mechanism of the drought response in cotton.

Loss-of-function of an α-SNAP gene confers resistance to soybean cyst nematode

Plant-parasitic nematodes are one of the most economically impactful pests in agriculture resulting in billions of dollars in realized annual losses worldwide. Soybean cyst nematode (SCN) is the number one biotic constraint on soybean production making it a priority for the discovery, validation and functional characterization of native plant resistance genes and genetic modes of action that can be deployed to improve soybean yield across the globe. Here, we present the discovery and functional characterization of a soybean resistance gene, GmSNAP02. We use unique bi-parental populations to fine-map the precise genomic location, and a combination of whole genome resequencing and gene fragment PCR amplifications to identify and confirm causal haplotypes. Lastly, we validate our candidate gene using CRISPR-Cas9 genome editing and observe a gain of resistance in edited plants. This demonstrates that the GmSNAP02 gene confers a unique mode of resistance to SCN through loss-of-function mutations that implicate GmSNAP02 as a nematode virulence target. We highlight the immediate impact of utilizing GmSNAP02 as a genome-editing-amenable target to diversify nematode resistance in commercially available cultivars.

Genome-wide identification of peanut IGT family genes and their potential roles in the development of plant architecture

IGT family genes play essential roles in shaping plant architecture. However, limited amount of information is available about IGT family genes in peanuts (Arachis hypogaea). In the current study, 13 AhIGT genes were identified and classified into three groups based on their phylogenetic relationship. Gene structure, conserved domain analyses indicated all AhIGTs were observed to share a similar exon-intron distribution pattern. AhIGTs within the same subfamily maintained a consistent motif composition. Chromosomal localization and synteny analyses showed that AhIGTs were unevenly localized on 9 chromosomes and that segmental duplication and purifying selection may have played important roles in the evolution of AhIGT genes. The analysis of conserved motifs, GO annotation, and transcript profile suggested that AhLAZY1-3 may play roles in gravity sensing and shaping peanut plant architecture. Transcript profile analysis suggested that AhTAC1 could potentially be involved gynophore ('peg') penetration into the soil. The cis-element analysis revealed that the light-responsive elements accounted for most of all cis-acting elements. Furthermore, qRT-PCR analysis showed that the expression of several AhIGT genes, like AhTAC1-2/4, was light-dependent, indicating that these genes may regulate plant architecture in response to light signals. This study may facilitate functional studies of the IGT genes in peanut.

Rht12b, a widely used ancient allele of TaGA2oxA13, reduces plant height and enhances yield potential in wheat

KEY MESSAGE

We identified a new wheat dwarfing allele Rht12b conferring reduced height and higher grain yield, pinpointed its causal variations, developed a breeding-applicable marker, and traced its origin and worldwide distribution. Plant height control is essential to optimize lodging resistance and yield gain in crops. RHT12 is a reduced height (Rht) locus that is identified in a mutationally induced dwarfing mutant and encodes a gibberellin 2-oxidase TaGA2oxA13. However, the artificial dwarfing allele is not used in wheat breeding due to excessive height reduction. Here, we confirmed a stable Rht locus, overlapping with RHT12, in a panel of wheat cultivars and its dwarfing allele reduced plant height by 5.4-8.2 cm, equivalent to Rht12b, a new allele of RHT12. We validated the effect of Rht12b on plant height in a bi-parent mapping population. Importantly, wheat cultivars carrying Rht12b had higher grain yield than those with the contrasting Rht12a allele. Rht12b conferred higher expression level of TaGA2oxA13. Transient activation assays defined SNP-390(C/A) in the promoter of TaGA2oxA13 as the causal variation. An efficient kompetitive allele-specific PCR marker was developed to diagnose Rht12b. Conjoint analysis showed that Rht12b plus the widely used Rht-D1b, Rht8 and Rht24b was the predominant Rht combination and conferred a moderate plant height in tested wheat cultivars. Evolutionary tracking uncovered that RHT12 locus arose from a tandem duplication event with Rht12b firstly appearing in wild emmer. The frequency of Rht12b was approximately 70% (700/1005) in a worldwide wheat panel and comparable to or higher than those of other widely used Rht genes, suggesting it had been subjected to positive selection. These findings not only identify a valuable Rht gene for wheat improvement but also develop a functionally diagnostic tool for marker-assisted breeding.

Identification of QTLs for resistance to 10 pathotypes of Plasmodiophora brassicae in Brassica oleracea cultivar ECD11 through genotyping-by-sequencing

KEY MESSAGE

Two major quantitative trait loci (QTLs) and five minor QTLs for 10 pathotypes were identified on chromosomes C01, C03, C04 and C08 through genotyping-by-sequencing from Brassica oleracea. Clubroot caused by Plasmodiophora brassicae is an important disease in brassica crops. Managing clubroot disease of canola on the Canadian prairie is challenging due to the continuous emergence of new pathotypes. Brassica oleracea is considered a major source of quantitative resistance to clubroot. Genotyping-by-sequencing (GBS) was performed in the parental lines; T010000DH3 (susceptible), ECD11 (resistant) and 124 BC1 plants. A total of 4769 high-quality polymorphic SNP loci were obtained and distributed on 9 chromosomes of B. oleracea. Evaluation of 124 BC1S1 lines for resistance to 10 pathotypes: 3A, 2B, 5C, 3D, 5G, 3H, 8J, 5K, 5L and 3O of P. brassicae, was carried out. Seven QTLs, 5 originating from ECD11 and 2 from T010000DH3, were detected. One major QTL designated as Rcr_C03-1 on C03 contributed 16.0-65.6% of phenotypic variation explained (PVE) for 8 pathotypes: 2B, 5C, 5G, 3H, 8J, 5K, 5L and 3O. Another major QTL designated as Rcr_C08-1 on C08 contributed 8.3 and 23.5% PVE for resistance to 8J and 5K, respectively. Five minor QTLs designated as Rcr_C01-1, Rcr_C03-2, Rcr_C03-3, Rcr_C04-1 and Rcr_C08-2 were detected on chromosomes C01, C03, C04 and C08 that contributed 8.3-23.5% PVE for 5 pathotypes each of 3A, 2B, 3D, 8J and 5K. There were 1, 10 and 4 genes encoding TIR-NBS-LRR/CC-NBS-LRR class disease resistance proteins in the Rcr_C01-1, Rcr_C03-1 and Rcr_C08-1 flanking regions. The syntenic regions of the two major QTLs Rcr_C03-1 and Rcr_C08-1 in the B. rapa genome 'Chiifu' were searched.

Unraveling the genetic basis of grain number-related traits in a wheat-Agropyron cristatum introgressed line through high-resolution linkage mapping

BACKGROUND

Grain number per spike (GNS) is a pivotal determinant of grain yield in wheat. Pubing 3228 (PB3228), a wheat-Agropyron cristatum germplasm, exhibits a notably higher GNS.

RESULTS

In this study, we developed a recombinant inbred line (RIL) population derived from PB3228/Gao8901 (PG-RIL) and constructed a high-density genetic map comprising 101,136 loci, spanning 4357.3 cM using the Wheat 660 K SNP array. The genetic map demonstrated high collinearity with the wheat assembly IWGSC RefSeq v1.0. Traits related to grain number and spikelet number per spike were evaluated in seven environments for quantitative trait locus (QTL) analysis. Five environmentally stable QTLs were detected in at least three environments. Among these, two major QTLs, QGns-4A.2 and QGns-1A.1, associated with GNS, exhibited positive alleles contributed by PB3228. Further, the conditional QTL analysis revealed a predominant contribution of PB3228 to the GNS QTLs, with both grain number per spikelet (GNSL) and spikelet number per spike (SNS) contributing to the overall GNS trait. Four kompetitive allele-specific PCR (KASP) markers that linked to QGns-4A.2 and QGns-1A.1 were developed and found to be effective in verifying the QTL effect within a diversity panel. Compared to previous studies, QGns-4A.2 exhibited stability across different trials, while QGns-1A.1 represents a novel QTL. The results from unconditional and conditional QTL analyses are valuable for dissecting the genetic contribution of the component traits to GNS at the individual QTL level and for understanding the genetic basis of the superior grain number character in PB3228. The KASP markers can be utilized in marker-assisted selection for enhancing GNS.

CONCLUSIONS

Five environmentally stable QTLs related to grain number and spikelet number per spike were identified. PB3228 contributed to the majority of the QTLs associated with GNS.

The Basic/Helix-Loop-Helix Transcription Factor Family Gene RcbHLH112 Is a Susceptibility Gene in Gray Mould Resistance of Rose (Rosa Chinensis)

The basic/helix-loop-helix (bHLH) family is a major family of transcription factors in plants. Although it has been reported that bHLH plays a defensive role against pathogen infection in plants, there is no comprehensive study on the bHLH-related defence response in rose (Rosa sp.). In this study, a genome-wide analysis of bHLH family genes (RcbHLHs) in rose was carried out, including their phylogenetic relationships, gene structure, chromosome localization and collinearity analysis. Via phylogenetic analysis, a total of 121 RcbHLH genes in the rose genome were divided into 21 sub-groups. These RcbHLHs are unevenly distributed in all 7 chromosomes of rose. The occurrence of gene duplication events indicates that whole-genome duplication and segmental duplication may play a key role in gene duplication. Ratios of non-synonymous to synonymous mutation frequency (Ka/Ks) analysis showed that the replicated RcbHLH genes mainly underwent purification selection, and their functional differentiation was limited. Gene expression analysis showed that 46 RcbHLHs were differentially expressed in rose petals upon B. cinerea infection. It is speculated that these RcbHLHs are candidate genes that regulate the response of rose plants to B. cinerea infection. Virus-induced gene silencing (VIGS) confirmed that RcbHLH112 in rose is a susceptibility factor for infection with B. cinerea. This study provides useful information for further study of the functions of the rose bHLH gene family.

Heat-responsive microRNAs participate in regulating the pollen fertility stability of CMS-D2 restorer line under high-temperature stress

Anther development and pollen fertility of cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) restorer lines are susceptible to continuous high-temperature (HT) stress in summer, which seriously hinders the large-scale application of "three-line" hybrids in production. Here, integrated small RNA, transcriptome, degradome, and hormone profiling was performed to explore the roles of microRNAs (miRNAs) in regulating fertility stability in mature pollens of isonuclear alloplasmic near-isogenic restorer lines NH and SH under HT stress at two environments. A total of 211 known and 248 novel miRNAs were identified, of which 159 were differentially expressed miRNAs (DEMs). Additionally, 45 DEMs in 39 miRNA clusters (PmCs) were also identified, and most highly expressed miRNAs were significantly induced in SH under extreme HT, especially four MIR482 and six MIR6300 family miRNAs. PmC28 was located in the fine-mapped interval of the Rf1 gene and contained two DEMs, gra-miR482_L-2R + 2 and gma-miR2118a-3p_R + 1_1ss18TG. Transcriptome sequencing identified 6281 differentially expressed genes, of which heat shock protein (HSP)-related genes, such as HSP70, HSP22, HSP18.5-C, HSP18.2 and HSP17.3-B, presented significantly reduced expression levels in SH under HT stress. Through integrating multi-omics data, we constructed a comprehensive molecular network of miRNA-mRNA-gene-KEGG containing 35 pairs of miRNA/target genes involved in regulating the pollen development in response to HT, among which the mtr-miR167a_R + 1, tcc-miR167c and ghr-miR390a, tcc-miR396c_L-1 and ghr-MIR169b-p3_1ss6AG regulated the pollen fertility by influencing ARF8 responsible for the auxin signal transduction, ascorbate and aldarate metabolism, and the sugar and lipid metabolism and transport pathways, respectively. Further combination with hormone analysis revealed that HT-induced jasmonic acid signaling could activate the expression of downstream auxin synthesis-related genes and cause excessive auxin accumulation, followed by a cascade of auxin signal transduction, ultimately resulting in pollen abortion. The results provide a new understanding of how heat-responsive miRNAs regulate the stability of fertility restoration for CMS-D2 cotton under heat stress.

Identification of the High-Affinity Potassium Transporter Gene Family (HKT) in Brassica U-Triangle Species and Its Potential Roles in Abiotic Stress in Brassica napus L

Members of the high-affinity potassium transporter (HKT) protein family regulate the uptake and homeostasis of sodium and potassium ions, but little research describes their roles in response to abiotic stresses in rapeseed (Brassica napus L.). In this study, we identified and characterized a total of 36 HKT genes from the species comprising the triangle of U model (U-triangle species): B. rapa, B. nigra, B. oleracea, B. juncea, B. napus, and B. carinata. We analyzed the phylogenetic relationships, gene structures, motif compositions, and chromosomal distributions of the HKT family members of rapeseed. Based on their phylogenetic relationships and assemblage of functional domains, we classified the HKT members into four subgroups, HKT1;1 to HKT1;4. Analysis of the nonsynonymous substitutions (Ka), synonymous substitutions (Ks), and the Ka/Ks ratios of HKT gene pairs suggested that these genes have experienced strong purifying selective pressure after duplication, with their evolutionary relationships supporting the U-triangle theory. Furthermore, the expression profiles of BnaHKT genes varies among potassium, phytohormone and heavy-metal treatment. Their repression provides resistance to heavy-metal stress, possibly by limiting uptake. Our results systematically reveal the characteristics of HKT family proteins and their encoding genes in six Brassica species and lay a foundation for further exploration of the role of HKT family genes in heavy-metal tolerance.

QTL analysis of traits related to seed size and shape in sesame (Sesamum indicum L.)

Seed size and shape are important traits that determine seed yield in sesame. Understanding the genetic basis of seed size and shape is essential for improving the yield of sesame. In this study, F2 and BC1 populations were developed by crossing the Yuzhi 4 and Bengal small-seed (BS) lines for detecting the quantitative trait loci (QTLs) of traits related to seed size and shape. A total of 52 QTLs, including 13 in F2 and 39 in BC1 populations, for seed length (SL), seed width (SW), and length to width ratio (L/W) were identified, explaining phenotypic variations from 3.68 to 21.64%. Of these QTLs, nine stable major QTLs were identified in the two populations. Notably, three major QTLs qSL-LG3-2, qSW-LG3-2, and qSW-LG3-F2 that accounted for 4.94-16.34% of the phenotypic variations were co-localized in a 2.08 Mb interval on chromosome 1 (chr1) with 279 candidate genes. Three stable major QTLs qSL-LG6-2, qLW-LG6, and qLW-LG6-F2 that explained 8.14-33.74% of the phenotypic variations were co-localized in a 3.27 Mb region on chr9 with 398 candidate genes. In addition, the stable major QTL qSL-LG5 was co-localized with minor QTLs qLW-LG5-3 and qSW-LG5 to a 1.82 Mb region on chr3 with 195 candidate genes. Gene annotation, orthologous gene analysis, and sequence analysis indicated that three genes are likely involved in sesame seed development. These results obtained herein provide valuable in-formation for functional gene cloning and improving the seed yield of sesame.

Transcriptome analysis reveals ZmERF055 contributes to waterlogging tolerance in sweetcorn

Waterlogging is a major disaster damaging crop production. However, most sweetcorn cultivars are not tolerant to waterlogging, which severely threatens their production. In order to understand the genetic mechanisms underlying waterlogging tolerance in sweetcorn, this study conducted a comprehensive investigation of sweetcorn waterlogging tolerance at the levels of physiology, biochemistry, and transcriptome in two sweetcorn CSSLs (chromosome segment substitution lines), D120 and D81. We found that D120 showed increased plant height, root length, root area, adventitious root numbers, antioxidant enzyme activities, and aerenchyma area ratio compared to D81. The transcriptome results showed that 2492 and 2351 differentially expressed genes (DEGs) were obtained at 4 h and 8 h of waterlogging treatment, respectively. Genes involved in reactive oxygen species (ROS) homeostasis, photosynthesis, and alcohol fermentation are sensitive in the waterlogging tolerant genotype D120, resulting in enhanced ROS scavenging ability, adventitious roots, and aerenchyma formation. Additionally, ethylene-, auxin-, and ABA-related genes exhibited different responses to waterlogging stress in sweetcorn. We integrated transcriptome and differential chromosomal fragments data and identified that ZmERF055 on chromosome 9 was directly involved in waterlogging stress. ZmERF055-overexpressing plants consistently exhibited significantly increased waterlogging tolerance and ROS homeostasis in Arabidopsis. These results offer a network of plant hormone signaling, ROS homeostasis, and energy metabolism co-modulating waterlogging tolerance in sweetcorn. Additionally, the findings support ZmERF055 as a potential ideal target gene in crop breeding to improve plant waterlogging tolerance.

Key processes required for the different stages of fungal carnivory by a nematode-trapping fungus

Nutritional deprivation triggers a switch from a saprotrophic to predatory lifestyle in soil-dwelling nematode-trapping fungi (NTF). In particular, the NTF Arthrobotrys oligospora secretes food and sex cues to lure nematodes to its mycelium and is triggered to develop specialized trapping devices. Captured nematodes are then invaded and digested by the fungus, thus serving as a food source. In this study, we examined the transcriptomic response of A. oligospora across the stages of sensing, trap development, and digestion upon exposure to the model nematode Caenorhabditis elegans. A. oligospora enacts a dynamic transcriptomic response, especially of protein secretion-related genes, in the presence of prey. Two-thirds of the predicted secretome of A. oligospora was up-regulated in the presence of C. elegans at all time points examined, and among these secreted proteins, 38.5% are predicted to be effector proteins. Furthermore, functional studies disrupting the t-SNARE protein Sso2 resulted in impaired ability to capture nematodes. Additionally, genes of the DUF3129 family, which are expanded in the genomes of several NTF, were highly up-regulated upon nematode exposure. We observed the accumulation of highly expressed DUF3129 proteins in trap cells, leading us to name members of this gene family as Trap Enriched Proteins (TEPs). Gene deletion of the most highly expressed TEP gene, TEP1, impairs the function of traps and prevents the fungus from capturing prey efficiently. In late stages of predation, we observed up-regulation of a variety of proteases, including metalloproteases. Following penetration of nematodes, these metalloproteases facilitate hyphal growth required for colonization of prey. These findings provide insights into the biology of the predatory lifestyle switch in a carnivorous fungus and provide frameworks for other fungal-nematode predator-prey systems.

HSF and Hsp Gene Families in sunflower: a comprehensive genome-wide determination survey and expression patterns under abiotic stress conditions

Sunflowers belong to the Asteraceae family, which comprises nutrimental and economic oilseed plants. Heat shock proteins (Hsps) are protein families vital for all organisms' growth and survival. Besides the ordinary conditions, the expression of these proteins ascends during abiotic stress factors such as high temperature, salinity, and drought. Using bioinformatics approaches, the current study identified and analyzed HSF and Hsp gene family members in the sunflower (Helianthus annuus L.) plant. HSF, sHsp, Hsp40, Hsp60, Hsp70, Hsp90, and Hsp100 domains were analyzed in the sunflower genome, and 88, 72, 192, 52, 85, 49, and 148 genes were identified, respectively. The motif structures of the proteins in the same phylogenetic tree were similar, and the α-helical form was dominant in all the protein families except for sHsp. The estimated three-dimensional structure of 28 sHsp proteins was determined as β-sheets. Considering protein-protein interactions, the Hsp60-09 protein (38 interactions) was found to be the most interacting protein. The most orthologous gene pairs (58 genes) were identified between Hsp70 genes and Arabidopsis genes. The expression analysis of selected genes was performed under high temperature, drought, and high temperature-drought combined stress conditions in two sunflower cultivars. In stress conditions, gene expressions were upregulated for almost all genes in the first half and first hours at large. The expressions of HanHSF-45 and HanHsp70-29 genes were raised in two cultivars under high temperature and high temperature-drought combined stress conditions. This study presents a blueprint for subsequent research and delivers comprehensive knowledge of this vital protein domain.

Development of an SNP marker set for marker-assisted backcrossing using genotyping-by-sequencing in tetraploid perilla

High-quality molecular markers are essential for marker-assisted selection to accelerate breeding progress. Compared with diploid species, recently diverged polyploid crop species tend to have highly similar homeologous subgenomes, which is expected to limit the development of broadly applicable locus-specific single-nucleotide polymorphism (SNP) assays. Furthermore, it is particularly challenging to make genome-wide marker sets for species that lack a reference genome. Here, we report the development of a genome-wide set of kompetitive allele specific PCR (KASP) markers for marker-assisted recurrent selection (MARS) in the tetraploid minor crop perilla. To find locus-specific SNP markers across the perilla genome, we used genotyping-by-sequencing (GBS) to construct linkage maps of two F2 populations. The two resulting high-resolution linkage maps comprised 2326 and 2454 SNP markers that spanned a total genetic distance of 2133 cM across 16 linkage groups and 2169 cM across 21 linkage groups, respectively. We then obtained a final genetic map consisting of 22 linkage groups with 1123 common markers from the two genetic maps. We selected 96 genome-wide markers for MARS and confirmed the accuracy of markers in the two F2 populations using a high-throughput Fluidigm system. We confirmed that 91.8% of the SNP genotyping results from the Fluidigm assay were the same as the results obtained through GBS. These results provide a foundation for marker-assisted backcrossing and the development of new varieties of perilla.

Two QTL regions for spike length showing pleiotropic effects on Fusarium head blight resistance and thousand-grain weight in bread wheat (Triticum aestivum L.)

Spike length (SL) plays an important role in the yield improvement of wheat and is significantly associated with other traits. Here, we used a recombinant inbred line (RIL) population derived from a cross between Yangmai 12 (YM12) and Yanzhan 1 (YZ1) to construct a genetic linkage map and identify quantitative trait loci (QTL) for SL. A total of 5 QTL were identified for SL, among which QSl.yaas-3A and QSl.yaas-5B are two novel QTL for SL. The YZ1 alleles at QSl.yaas-2D and QSl.yaas-5A, and the YM12 alleles at QSl.yaas-2A, QSl.yaas-3A, and QSl.yaas-5B conferred increasing SL effects. Two major QTL QSl.yaas-5A and QSl.yaas-5B explained 9.11-15.85% and 9.01-12.85% of the phenotypic variations, respectively. Moreover, the positive alleles of QSl.yaas-5A and QSl.yaas-5B could significantly increase Fusarium head blight (FHB) resistance (soil surface inoculation and spray inoculation were used) and thousand-grain weight (TGW) in the RIL population. Kompetitive allele-specific PCR (KASP) markers for QSl.yaas-5A and QSl.yaas-5B were developed and validated in an additional panel of 180 wheat cultivars/lines. The cultivars/lines harboring both the positive alleles of QSl.yaas-5A and QSl.yaas-5B accounted for only 28.33% of the validation populations and had the longest SL, best FHB resistance (using spray inoculation), and highest TGW. A total of 358 and 200 high-confidence annotated genes in QSl.yaas-5A and QSl.yaas-5B were identified, respectively. Some of the genes in these two regions were involved in cell development, disease resistance, and so on. The results of this study will provide a basis for directional breeding of longer SL, higher TGW, and better FHB resistance varieties and a solid foundation for fine-mapping QSl.yaas-5A and QSl.yaas-5B in future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01427-8.

Genome-wide analysis of the LAZ1 gene family in Gossypium hirsutum

BACKGROUND

As the world's leading fiber crop and a major oil-producing crop, cotton fiber yield and fiber quality are affected by environmental stresses, especially heat, drought and salinity. The LAZ1 (Lazarus 1) family genes are responsive to abscisic acid, drought, and salt treatments. Currently, mining and functional analyses of LAZ1 family genes in cotton have not been reported.

METHODS AND RESULTS

In this study, 20 GhLAZ1 genes, designated GhLAZ1-1 - GhLAZ1-20, were identified in the genome of Gossypium hirsutum through the construction of an HMM model, and their molecular properties, chromosomal localization, phylogeny, gene structure, evolutionary selection pressure, promoter cis elements and gene expression under salt stress were analyzed. With the exception of GhLAZ1-17 and GhLAZ1-20, the remaining 18 GhLAZ1 genes were unevenly localized on 13 chromosomes in G. hirsutum; evolutionary analysis showed that these genes could be divided into three subfamilies; and evolutionary selection pressure analysis demonstrated that the GhLAZ1 genes were all under purifying selection. Many elements related to light responses, hormone responses, and abiotic stresses were predicted on the GhLAZ1 family gene promoters, and real-time quantitative PCR results showed that GhLAZ1-2, GhLAZ1-8, and GhLAZ1-18 were upregulated significantly in salt-treated cotton leaves.

CONCLUSIONS

Our results suggested that GhLAZ1 genes were involved in the salt tolerance mechanism in G. hirsutum and provided a reference for further exploring the function and molecular mechanism of LAZ1 genes.

Molecular Mapping of Yr85 and Comparison with Other Genes for Resistance to Stripe Rust on Wheat Chromosome 1B

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici, is one of the most serious plant diseases worldwide. Resistant cultivars are the most effective way to control the disease. YrTr1 is an important stripe rust resistance gene that has been used in wheat breeding programs and is represented in the host differential set to identify P. striiformis f. sp. tritici races in the United States. To map YrTr1, AvSYrTr1NIL was backcrossed to its recurrent parent Avocet S (AvS). Seedlings of BC7F2, BC7F3, and BC8F1 populations were tested with YrTr1-avirulent races under controlled conditions, and BC7F2 plants were genotyped using simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers. YrTr1 was mapped to the short arm of chromosome 1B using four SSR and seven SNP markers. The genetic distances of YrTr1 from the nearest flanking markers IWA2583 and IWA7480 were 1.8 and 1.3 centimorgans (cM), respectively. DNA amplification of a set of 21 Chinese Spring (CS) nulli-tetrasomic lines and seven CS 1B deletion lines with three SSR markers confirmed the chromosome arm location and further placed the gene in chromosomal bin region 1BS18 (0.5). The gene was determined to be about 7.4 cM proximal to Yr10. Based on multirace response array and chromosomal location, YrTr1 was determined to be different from other permanently named stripe rust resistance genes in chromosome arm 1BS and was named Yr85.