Identification and characterization of the CRK gene family in the wheat genome and analysis of their expression profile in response to high temperature-induced male sterility.

Cysteine-rich receptor-like kinases (CRKs) play many important roles during plant development, including defense responses under both biotic and abiotic stress, reactive oxygen species (ROS) homeostasis, callose deposition and programmed cell death (PCD). However, there are few studies on the involvement of the CRK family in male sterility due to heat stress in wheat (Triticum aestivum L.). In this study, a genome-wide characterization of the CRK family was performed to investigate the structural and functional attributes of the wheat CRKs in anther sterility caused by heat stress. A total of 95 CRK genes were unevenly distributed on 18 chromosomes, with the most genes distributed on chromosome 2B. Paralogous homologous genes with Ka/Ks ratios less than 1 may have undergone strong purifying selection during evolution and are more functionally conserved. The collinearity analysis results of CRK genes showed that wheat and Arabidopsis (A. thaliana), foxtail millet, Brachypodium distachyon (B. distachyon), and rice have three, 12, 15, and 11 pairs of orthologous genes, respectively. In addition, the results of the network interactions of genes and miRNAs showed that five miRNAs were in the hub of the interactions map, namely tae-miR9657b-5p, tae-miR9780, tae-miR9676-5p, tae-miR164, and tae-miR531. Furthermore, qRT-PCR validation of the six TaCRK genes showed that they play key roles in the development of the mononuclear stage anthers, as all six genes were expressed at highly significant levels in heat-stressed male sterile mononuclear stage anthers compared to normal anthers. We hypothesized that the TaCRK gene is significant in the process of high-temperature-induced sterility in wheat based on the combination of anther phenotypes, paraffin sections, and qRT-PCR data. These results improve our understanding of their relationship.

Genetic characterization and linkage analysis of spotted leaf 6,liguleless and lax panicle traits in mutant rice.

Phenotypic mutants are valuable resources for elucidating the function of genes responsible for their expression. This study examined mutant rice strains expressing three traits: spotted leaf 6 (spl6), lax panicle (lax), and liguleless (lg). In the mutant, the spl6 phenotype was a genetically programmed lesion-mimicking mutation (LMM) that displayed spontaneously scattered spots across the leaf surface. In the lg trait, the plant lacked a collar region, and there were no auricles and ligules at the junction of the leaf blade and leaf sheath. The lax panicle trait manifested as sparely arranged spikelets resulting from the terminal spikelet with no lateral spikelets, which caused a drastic reduction of the total seed number in the mutant. All three mutant genes were genetically recessive and had nuclear gene regulation. The dihybrid segregation of the lg gene was classified independently according to the Mendelian 9:3:3:1 dihybrid segregation ratio in the F2 generation, suggesting that the lg gene is not linked to the same chromosome as the lax and spl6 genes. On the other hand, spl6 and lax were not assorted independently, indicating that they are closely linked on chromosome 1 in rice. Additional linkage analysis from the recombination of spl6 and lax genes reconfirmed that the two genes were ~9.4 cM away from each other. The individual single-gene mutant plant from one plant with a three-gene mutation (spl6, lax, and lg) was isolated and characterized, which will be a crucial resource for the gene cloning and molecular characterization of these genes.

RNA-Seq Bulked Segregant Analysis of an Exotic B. napus ssp. napobrassica (Rutabaga) F2 Population Reveals Novel QTLs for Breeding Clubroot-Resistant Canola.

In this study, a rutabaga (Brassica napus ssp. napobrassica) donor parent FGRA106, which exhibited broad-spectrum resistance to 17 isolates representing 16 pathotypes of Plasmodiophora brassicae, was used in genetic crosses with the susceptible spring-type canola (B. napus ssp. napus) accession FG769. The F2 plants derived from a clubroot-resistant F1 plant were screened against three P. brassicae isolates representing pathotypes 3A, 3D, and 3H. Chi-square (χ2) goodness-of-fit tests indicated that the F2 plants inherited two major clubroot resistance genes from the CR donor FGRA106. The total RNA from plants resistant (R) and susceptible (S) to each pathotype were pooled and subjected to bulked segregant RNA-sequencing (BSR-Seq). The analysis of gene expression profiles identified 431, 67, and 98 differentially expressed genes (DEGs) between the R and S bulks. The variant calling method indicated a total of 12 (7 major + 5 minor) QTLs across seven chromosomes. The seven major QTLs included: BnaA5P3A.CRX1.1, BnaC1P3H.CRX1.2, and BnaC7P3A.CRX1.1 on chromosomes A05, C01, and C07, respectively; and BnaA8P3D.CRX1.1, BnaA8P3D.RCr91.2/BnaA8P3H.RCr91.2, BnaA8P3H.Crr11.3/BnaA8P3D.Crr11.3, and BnaA8P3D.qBrCR381.4 on chromosome A08. A total of 16 of the DEGs were located in the major QTL regions, 13 of which were on chromosome C07. The molecular data suggested that clubroot resistance in FGRA106 may be controlled by major and minor genes on both the A and C genomes, which are deployed in different combinations to confer resistance to the different isolates. This study provides valuable germplasm for the breeding of clubroot-resistant B. napus cultivars in Western Canada.

Identification of Quantitative Trait Loci Controlling Root Morphological Traits in an Interspecific Soybean Population Using 2D Imagery Data.

Roots are the hidden and most important part of plants. They serve as stabilizers and channels for uptaking water and nutrients and play a crucial role in the growth and development of plants. Here, two-dimensional image data were used to identify quantitative trait loci (QTL) controlling root traits in an interspecific mapping population derived from a cross between wild soybean 'PI366121' and cultivar 'Williams 82'. A total of 2830 single-nucleotide polymorphisms were used for genotyping, constructing genetic linkage maps, and analyzing QTLs. Forty-two QTLs were identified on twelve chromosomes, twelve of which were identified as major QTLs, with a phenotypic variation range of 36.12% to 39.11% and a logarithm of odds value range of 12.01 to 17.35. Two significant QTL regions for the average diameter, root volume, and link average diameter root traits were detected on chromosomes 3 and 13, and both wild and cultivated soybeans contributed positive alleles. Six candidate genes, Glyma.03G027500 (transketolase/glycoaldehyde transferase), Glyma.03G014500 (dehydrogenases), Glyma.13G341500 (leucine-rich repeat receptor-like protein kinase), Glyma.13G341400 (AGC kinase family protein), Glyma.13G331900 (60S ribosomal protein), and Glyma.13G333100 (aquaporin transporter) showed higher expression in root tissues based on publicly available transcriptome data. These results will help breeders improve soybean genetic components and enhance soybean root morphological traits using desirable alleles from wild soybeans.

Chromosome-level genome assembly of the threatened resource plant Cinnamomum chago.

Cinnamomum chago is a tree species endemic to Yunnan province, China, with potential economic value, phylogenetic importance, and conservation priority. We assembled the genome of C. chago using multiple sequencing technologies, resulting in a high-quality, chromosomal-level genome with annotation information. The assembled genome size is approximately 1.06 Gb, with a contig N50 length of 92.10 Mb. About 99.92% of the assembled sequences could be anchored to 12 pseudo-chromosomes, with only one gap, and 63.73% of the assembled genome consists of repeat sequences. In total, 30,497 genes were recognized according to annotation, including 28,681 protein-coding genes. This high-quality chromosome-level assembly and annotation of C. chago will assist us in the conservation and utilization of this valuable resource, while also providing crucial data for studying the evolutionary relationships within the Cinnamomum genus, offering opportunities for further research and exploration of its diverse applications.

Genome-wide identification and expression analysis of the SET domain-containing gene family in potato (Solanum tuberosum L.).

Genes containing the SET domain can catalyse histone lysine methylation, which in turn has the potential to cause changes to chromatin structure and regulation of the transcription of genes involved in diverse physiological and developmental processes. However, the functions of SET domain-containing (StSET) genes in potato still need to be studied. The objectives of our study can be summarized as in silico analysis to (i) identify StSET genes in the potato genome, (ii) systematically analyse gene structure, chromosomal distribution, gene duplication events, promoter sequences, and protein domains, (iii) perform phylogenetic analyses, (iv) compare the SET domain-containing genes of potato with other plant species with respect to protein domains and orthologous relationships, (v) analyse tissue-specific expression, and (vi) study the expression of StSET genes in response to drought and heat stresses. In this study, we identified 57 StSET genes in the potato genome, and the genes were physically mapped onto eleven chromosomes. The phylogenetic analysis grouped these StSET genes into six clades. We found that tandem duplication through sub-functionalisation has contributed only marginally to the expansion of the StSET gene family. The protein domain TDBD (PFAM ID: PF16135) was detected in StSET genes of potato while it was absent in all other previously studied species. This study described three pollen-specific StSET genes in the potato genome. Expression analysis of four StSET genes under heat and drought in three potato clones revealed that these genes might have non-overlapping roles under different abiotic stress conditions and durations. The present study provides a comprehensive analysis of StSET genes in potatoes, and it serves as a basis for further functional characterisation of StSET genes towards understanding their underpinning biological mechanisms in conferring stress tolerance.

Genome-wide characterization of cytochrome P450 genes reveals the potential roles in fruit ripening and response to cold stress in tomato.

Plant cytochrome P450 (CYP) superfamily, the largest enzyme metabolism family, has been identified in many species and plays a vital role in plant development and stress response via secondary metabolite biosynthesis. A comprehensive identification and functional investigation of CYPs in tomato plants would contribute to deeper understanding of their biological significance. In this study, 268 tomato CYP genes were identified and found to be unevenly located on 12 chromosomes. Based on the phylogenetic analysis, these 268 SlCYPs were classed into two distinct clades (A-type and non-A-type) and nine clans, including 48 families. Moreover, 67 tandem and 22 WGD (whole genome duplication)/segmental duplication events were detected, of which 12 SlCYP genes experienced both WGD/segmental and tandem duplication events, indicating that tandem duplication plays a major role in the expansion of the SlCYP family. Besides, 48 pairs containing 41 SlCYP and 44 AtCYP genes were orthologous, while 216 orthologous pairs were obtained between tomato and potato. The expression level of all SlCYP genes in tomato tissues at different development stages was analyzed, and most expressed SlCYPs showed a tissue-specific pattern. Meanwhile, 143 differentially expressed SlCYPs were identified under cold stress. Furthermore, the RT-qPCR results indicated that SlCYPs may be involved in fruit ripening and cold tolerance in tomato seedlings. These findings provide valuable insights into the evolutionary relationships and functional characteristics of SlCYPs, which can be utilized for further investigation of fruit metabolic pathways and cold tolerance in tomato.

Precision mapping and expression analysis of recessive bacterial blight resistance gene xa-45(t) from Oryza glaberrima.

BACKGROUND: Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is one of the most devastating diseases of rice leading to huge yield losses in Southeast Asia. The recessive resistance gene xa-45(t) from Oryza glaberrima IRGC102600B, mapped on rice chromosome 8, spans 80 Kb with 9 candidate genes on Nipponbare reference genome IRGSP-1.0. The xa-45(t) gene provides durable resistance against all the ten Xanthomonas pathotypes of Northern India, thus aiding in the expansion of recessive bacterial blight resistance gene pool. Punjab Rice PR127, carrying xa-45(t), was released for wider use in breeding programs. This study aims to precisely locate the target gene among the 9 candidates conferring resistance to bacterial blight disease. METHODS AND RESULTS: Sanger sequencing of all nine candidate genes revealed seven SNPs and an Indel between the susceptible parent Pusa 44 and the resistant introgression line IL274. The genotyping with polymorphic markers identified three recombinant breakpoints for LOC_Os08g42370, and LOC_Os08g42400, 15 recombinants for LOC_Os08g423420 and 26 for LOC_Os08g42440 out of 190 individuals. Relative expression analysis across six time intervals (0, 8, 24, 48, 72, and 96 h) after bacterial blight infection showed over expression of LOC_Os08g42410-specific transcripts in IL274 compared to Pusa 44, with a significant 4.46-fold increase observed at 72 h post-inoculation. CONCLUSIONS: The Indel marker at the locus LOC_Os08g42410 was found co-segregating with the phenotype, suggesting its candidacy towards xa-45(t). The transcript abundance assay provides strong evidence for the involvement of LOC_Os08g42410 in the resistance conferred by the bacterial blight gene xa-45(t).

Genome-wide prediction and functional analysis of WOX genes in blueberry.

BACKGROUND: WOX genes are a class of plant-specific transcription factors. The WUSCHEL-related homeobox (WOX) family is a member of the homeobox transcription factor superfamily. Previous studies have shown that WOX members play important roles in plant growth and development. However, studies of the WOX gene family in blueberry plants have not been reported. RESULTS: In order to understand the biological function of the WOX gene family in blueberries, bioinformatics were used methods to identify WOX gene family members in the blueberry genome, and analyzed the basic physical and chemical properties, gene structure, gene motifs, promoter cis-acting elements, chromosome location, evolutionary relationships, expression pattern of these family members and predicted their functions. Finally, 12 genes containing the WOX domain were identified and found to be distributed on eight chromosomes. Phylogenetic tree analysis showed that the blueberry WOX gene family had three major branches: ancient branch, middle branch, and WUS branch. Blueberry WOX gene family protein sequences differ in amino acid number, molecular weight, isoelectric point and hydrophobicity. Predictive analysis of promoter cis-acting elements showed that the promoters of the VdWOX genes contained abundant light response, hormone, and stress response elements. The VdWOX genes were induced to express in both stems and leaves in response to salt and drought stress. CONCLUSIONS: Our results provided comprehensive characteristics of the WOX gene family and important clues for further exploration of its role in the growth, development and resistance to various stress in blueberry plants.

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture.

KEY MESSAGE: A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Fine-mapping and validation of the major quantitative trait locus QFlANG-4B for flag leaf angle in wheat.

KEY MESSAGE: This study precisely mapped and validated a quantitative trait locus (QTL) located on chromosome 4B for flag leaf angle in wheat. Flag leaf angle (FLANG) is closely related to crop architecture and yield. We previously identified the quantitative trait locus (QTL) QFLANG-4B for FLANG on chromosome 4B, located within a 14-cM interval flanked by the markers Xbarc20 and Xzyh357, using a mapping population of recombinant inbred lines (RILs) derived from a cross between Nongda3331 (ND3331) and Zang1817. In this study, we fine-mapped QFLANG-4B and validated its associated genetic effect. We developed a BC3F3 population using ND3331 as the recurrent parent through marker-assisted selection, as well as near-isogenic lines (NILs) by selfing BC3F3 plants carrying different heterozygous segments for the QFLANG-4B region. We obtained eight recombinant types for QFLANG-4B, narrowing its location down to a 5.3-Mb region. This region contained 76 predicted genes, 7 of which we considered to be likely candidate genes for QFLANG-4B. Marker and phenotypic analyses of individual plants from the secondary mapping populations and their progeny revealed that the FLANG of the ND3331 allele is significantly higher than that of the Zang1817 allele in multiple environments. These results not only provide a basis for the map-based cloning of QFLANG-4B, but also indicate that QFLANG-4B has great potential for marker-assisted selection in wheat breeding programs designed to improve plant architecture and yield.

Chromosome stability of synthetic Triticum turgidum-Aegilops umbellulata hybrids.

BACKGROUND: Unreduced gamete formation during meiosis plays a critical role in natural polyploidization. However, the unreduced gamete formation mechanisms in Triticum turgidum-Aegilops umbellulata triploid F1 hybrid crosses and the chromsome numbers and compostions in T. turgidum-Ae. umbellulata F2 still not known. RESULTS: In this study, 11 T.turgidum-Ae. umbellulata triploid F1 hybrid crosses were produced by distant hybridization. All of the triploid F1 hybrids had 21 chromosomes and two basic pathways of meiotic restitution, namely first-division restitution (FDR) and single-division meiosis (SDM). Only FDR was found in six of the 11 crosses, while both FDR and SDM occurred in the remaining five crosses. The chromosome numbers in the 127 selfed F2 seeds from the triploid F1 hybrid plants of 10 crosses (no F2 seeds for STU 16) varied from 35 to 43, and the proportions of euploid and aneuploid F2 plants were 49.61% and 50.39%, respectively. In the aneuploid F2 plants, the frequency of chromosome loss/gain varied among genomes. The chromosome loss of the U genome was the highest (26.77%) among the three genomes, followed by that of the B (22.83%) and A (11.81%) genomes, and the chromosome gain for the A, B, and U genomes was 3.94%, 3.94%, and 1.57%, respectively. Of the 21 chromosomes, 7U (16.54%), 5 A (3.94%), and 1B (9.45%) had the highest loss frequency among the U, A, and B genomes. In addition to chromosome loss, seven chromosomes, namely 1 A, 3 A, 5 A, 6 A, 1B, 1U, and 6U, were gained in the aneuploids. CONCLUSION: In the aneuploid F2 plants, the frequency of chromosome loss/gain varied among genomes, chromsomes, and crosses. In addition to variations in chromosome numbers, three types of chromosome translocations including 3UL·2AS, 6UL·1AL, and 4US·6AL were identified in the F2 plants. Furthermore, polymorphic fluorescence in situ hybridization karyotypes for all the U chromosomes were also identified in the F2 plants when compared with the Ae. umbellulata parents. These results provide useful information for our understanding the naturally occurred T. turgidum-Ae. umbellulata amphidiploids.

The holocentricity in the dioecious nutmeg (Myristica fragrans) is not based on major satellite repeats.

Holocentric species are characterized by the presence of centromeres throughout the length of the chromosomes. We confirmed the holocentricity of the dioecious, small chromosome-size species Myristica fragrans based on the chromosome-wide distribution of the centromere-specific protein KNL1, α-tubulin fibers, and the cell cycle-dependent histone H3 serine 28 phosphorylation (H3S28ph) mark. Each holocentromere is likely composed of, on average, ten centromere units, but none of the identified and in situ hybridized high-copy satellite repeats is centromere-specific. No sex-specific major repeats are present in the high-copy repeat composition of male or female plants, or a significant difference in genome size was detected. Therefore, it is unlikely that M. fragrans possesses heteromorphic sex chromosomes.

Identification of candidate genes for Fusarium head blight resistance from QTLs using RIL population in wheat.

Fusarium head blight (FHB) stands out as one of the most devastating wheat diseases and leads to significantly grain yield losses and quality reductions in epidemic years. Exploring quantitative trait loci (QTL) for FHB resistance is a critical step for developing new FHB-resistant varieties. We previously constructed a genetic map of unigenes (UG-Map) according to the physical positions using a set of recombinant-inbred lines (RILs) derived from the cross of 'TN18 × LM6' (TL-RILs). Here, the number of diseased spikelets (NDS) and relative disease index (RDI) for FHB resistance were investigated under four environments using TL-RILs, which were distributed across 13 chromosomes. A number of 36 candidate genes for NDS and RDI from of 19 stable QTLs were identified. The average number of candidate genes per QTL was 1.89, with 14 (73.7%), two (10.5%), and three (15.8%) QTLs including one, two, and 3-10 candidate genes, respectively. Among the 24 candidate genes annotated in the reference genome RefSeq v1.1, the homologous genes of seven candidate genes, including TraesCS4B02G227300 for QNds/Rdi-4BL-4553, TraesCS5B02G303200, TraesCS5B02G303300, TraesCS5B02G303700, TraesCS5B02G303800 and TraesCS5B02G304000 for QNds/Rdi-5BL-9509, and TraesCS7A02G568400 for QNds/Rdi-7AL-14499, were previously reported to be related to FHB resistance in wheat, barely or Brachypodium distachyon. These genes should be closely associated with FHB resistance in wheat. In addition, the homologous genes of five genes, including TraesCS1A02G037600LC for QNds-1AS-2225, TraesCS1D02G017800 and TraesCS1D02G017900 for QNds-1DS-527, TraesCS1D02G018000 for QRdi-1DS-575, and TraesCS4B02G227400 for QNds/Rdi-4BL-4553, were involved in plant defense responses against pathogens. These genes should be likely associated with FHB resistance in wheat.

Identification of soybean mosaic virus strain SC7 resistance loci and candidate genes in soybean [Glycine max (L.) Merr.].

Soybean [Glycine max (L.) Merr.] is an important legume crop worldwide, which provides abundant plant protein and oil for human beings. Soybean mosaic virus (SMV) can cause serious damage to the yield and quality of soybean, but it is difficult to control SMV with chemicals, breeding SMV-resistant varieties has become the most effective way to control the disease. Therefore, it is important to identify SMV resistance genes from soybean resources and apply them to soybean breeding. In this study, the disease rates (DRs) of 219 soybean accessions to SMV strain SC7 in two environments were investigated. A high-density NJAU 355 K SoySNP array was used for genome-wide association study (GWAS) of DR. A 274 kb region on chromosome 15 (1,110,567 bp to 1,384,173 bp) was repeatedly detected in two environments. Six new significant single nucleotide polymorphisms (SNPs) on chromosome 15 were identified. Four of these six SNPs were located within two candidate genes, Glyma.15G015700 and Glyma.15G015800. The elite haplotype Glyma.15G015700Hap I with low DR exhibited strong resistance to SC7. The expression of Glyma.15G015700 in the SMV-resistant accession increased significantly after inoculation with SC7. Furthermore, most of the proteins predicted to interact with Glyma.15G015700 are heat shock proteins, which have been shown to be related to disease resistance. In summary, new SMV resistance loci and a new candidate gene, Glyma.15G015700, were identified and might be utilized in further soybean disease resistance breeding.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis.

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Genome-wide identification and expression analysis of the Eriobotrya japonica TIFY gene family reveals its functional diversity under abiotic stress conditions.

BACKGROUND: Plant-specific TIFY proteins are widely found in terrestrial plants and play important roles in plant adversity responses. Although the genome of loquat at the chromosome level has been published, studies on the TIFY family in loquat are lacking. Therefore, the EjTIFY gene family was bioinformatically analyzed by constructing a phylogenetic tree, chromosomal localization, gene structure, and adversity expression profiling in this study. RESULTS: Twenty-six EjTIFY genes were identified and categorized into four subfamilies (ZML, JAZ, PPD, and TIFY) based on their structural domains. Twenty-four EjTIFY genes were irregularly distributed on 11 of the 17 chromosomes, and the remaining two genes were distributed in fragments. We identified 15 covariate TIFY gene pairs in the loquat genome, 13 of which were involved in large-scale interchromosomal segmental duplication events, and two of which were involved in tandem duplication events. Many abiotic stress cis-elements were widely present in the promoter region. Analysis of the Ka/Ks ratio showed that the paralogous homologs of the EjTIFY family were mainly subjected to purifying selection. Analysis of the RNA-seq data revealed that a total of five differentially expressed genes (DEGs) were expressed in the shoots under gibberellin treatment, whereas only one gene was significantly differentially expressed in the leaves; under both low-temperature and high-temperature stresses, there were significantly differentially expressed genes, and the EjJAZ15 gene was significantly upregulated under both low- and high-temperature stress. RNA-seq and qRT-PCR expression analysis under salt stress conditions revealed that EjJAZ2, EjJAZ4, and EjJAZ9 responded to salt stress in loquat plants, which promoted resistance to salt stress through the JA pathway. The response model of the TIFY genes in the jasmonic acid pathway under salt stress in loquat was systematically summarized. CONCLUSIONS: These results provide a theoretical basis for exploring the characteristics and functions of additional EjTIFY genes in the future. This study also provides a theoretical basis for further research on breeding for salt stress resistance in loquat. RT-qPCR analysis revealed that the expression of one of the three EjTIFY genes increased and the expression of two decreased under salt stress conditions, suggesting that EjTIFY exhibited different expression patterns under salt stress conditions.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii.

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

Clpf encodes pentatricopeptide repeat protein (PPR5) and regulates pink flesh color in watermelon (Citrullus lanatus L.).

KEY MESSAGE: The gene controlling pink flesh in watermelon was finely mapped to a 55.26-kb region on chromosome 6. The prime candidate gene, Cla97C06G122120 (ClPPR5), was identified through forward genetics. Carotenoids offer numerous health benefits; while, they cannot be synthesized by the human body. Watermelon stands out as one of the richest sources of carotenoids. In this study, genetic generations derived from parental lines W15-059 (red flesh) and JQ13-3 (pink flesh) revealed the presence of the recessive gene Clpf responsible for the pink flesh (pf) trait in watermelon. Comparative analysis of pigment components and microstructure indicated that the disparity in flesh color between the parental lines primarily stemmed from variations in lycopene content, as well as differences in chromoplast number and size. Subsequent bulk segregant analysis (BSA-seq) and genetic mapping successfully narrowed down the Clpf locus to a 55.26-kb region on chromosome 6, harboring two candidate genes. Through sequence comparison and gene expression analysis, Cla97C06G122120 (annotated as a pentatricopeptide repeat, PPR) was predicted as the prime candidate gene related to pink flesh trait. To further investigate the role of the PPR gene, its homologous gene in tomato was silenced using a virus-induced system. The resulting silenced fruit lines displayed diminished carotenoid accumulation compared with the wild-type, indicating the potential regulatory function of the PPR gene in pigment accumulation. This study significantly contributes to our understanding of the forward genetics underlying watermelon flesh traits, particularly in relation to carotenoid accumulation. The findings lay essential groundwork for elucidating mechanisms governing pigment synthesis and deposition in watermelon flesh, thereby providing valuable insights for future breeding strategies aimed at enhancing fruit quality and nutritional value.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety.

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.