Genome-wide investigation of the TIFY transcription factors in alfalfa (Medicago sativa L.): identification, analysis, and expression.

BACKGROUND: Alfalfa (Medicago sativa L.) is an essential leguminous forage with high nutrition and strong adaptability. The TIFY family is a plant-specific transcription factor identified in many plants. However, few reports have been reported on the phylogenetic analysis and gene expression profiling of TIFY family genes in alfalfa. RESULT: A total of 84 TIFY genes belonging to 4 categories were identified in alfalfa, including 58 MsJAZs, 18 MsZMLs, 4 MsTIFYs and 4 MsPPDs, respectively. qRT-PCR data from 8 genes in different tissues revealed that most MsTIFY genes were highly expressed in roots. The expression of MsTIFY14 was up-regulated after different times in both thrips-resistant and susceptible alfalfa after thrips feeding, and the expression of the remaining MsTIFYs had a strong correlation with the time of thrips feeding. Different abiotic stresses, including drought, salt, and cold, could induce or inhibit the expression of MsTIFY genes to varying degrees. In addition, the eight genes were all significantly up-regulated by JA and/or SA. Interestingly, MsTIFY77 was induced considerably by all the biotic, abiotic, or plant hormones (JA or SA) except ABA. CONCLUSION: Our study identified members of the TIFY gene family in alfalfa and analyzed their structures and possible functions. It laid the foundation for further research on the molecular functions of TIFYs in alfalfa.

Genome-wide analysis of the passion fruit invertase gene family reveals involvement of PeCWINV5 in hexose accumulation.

BACKGROUND: Invertases (INVs) are key enzymes in sugar metabolism, cleaving sucrose into glucose and fructose and playing an important role in plant development and the stress response, however, the INV gene family in passion fruit has not been systematically reported. RESULTS: In this study, a total of 16 PeINV genes were identified from the passion fruit genome and named according to their subcellular location and chromosome position. These include six cell wall invertase (CWINV) genes, two vacuolar invertase (VINV) genes, and eight neutral/alkaline invertase (N/AINV) genes. The gene structures, phylogenetic tree, and cis-acting elements of PeINV gene family were predicted using bioinformatics methods. Results showed that the upstream promoter region of the PeINV genes contained various response elements; particularly, PeVINV2, PeN/AINV3, PeN/AINV5, PeN/AINV6, PeN/AINV7, and PeN/AINV8 had more response elements. Additionally, the expression profiles of PeINV genes under different abiotic stresses (drought, salt, cold temperature, and high temperature) indicated that PeCWINV5, PeCWINV6, PeVINV1, PeVINV2, PeN/AINV2, PeN/AINV3, PeN/AINV6, and PeN/AINV7 responded significantly to these abiotic stresses, which was consistent with cis-acting element prediction results. Sucrose, glucose, and fructose are main soluble components in passion fruit pulp. The contents of total soluble sugar, hexoses, and sweetness index increased significantly at early stages during fruit ripening. Transcriptome data showed that with an increase in fruit development and maturity, the expression levels of PeCWINV2, PeCWINV5, and PeN/AINV3 exhibited an up-regulated trend, especially for PeCWINV5 which showed highest abundance, this correlated with the accumulation of soluble sugar and sweetness index. Transient overexpression results demonstrated that the contents of fructose, glucose and sucrose increased in the pulp of PeCWINV5 overexpressing fruit. It is speculated that this cell wall invertase gene, PeCWINV5, may play an important role in sucrose unloading and hexose accumulation. CONCLUSION: In this study, we systematically identified INV genes in passion fruit for the first time and further investigated their physicochemical properties, evolution, and expression patterns. Furthermore, we screened out a key candidate gene involved in hexose accumulation. This study lays a foundation for further study on INV genes and will be beneficial on the genetic improvement of passion fruit breeding.

Genome-wide identification of kiwifruit K+ channel Shaker family members and their response to low-K+ stress.

BACKGROUND: 'Hongyang' kiwifruit (Actinidia chinensis cv 'Hongyang') is a high-quality variety of A. chinensis with the advantages of high yield, early ripening, and high stress tolerance. Studies have confirmed that the Shaker K+ genes family is involved in plant uptake and distribution of potassium (K+). RESULTS: Twenty-eight Shaker genes were identified and analyzed from the 'Hongyang' kiwifruit (A. chinensis cv 'Hongyang') genome. Subcellular localization results showed that more than one-third of the AcShaker genes were on the cell membrane. Phylogenetic analysis indicated that the AcShaker genes were divided into six subfamilies (I-VI). Conservative model, gene structure, and structural domain analyses showed that AcShaker genes of the same subfamily have similar sequence features and structure. The promoter cis-elements of the AcShaker genes were classified into hormone-associated cis-elements and environmentally stress-associated cis-elements. The results of chromosomal localization and duplicated gene analysis demonstrated that AcShaker genes were distributed on 18 chromosomes, and segmental duplication was the prime mode of gene duplication in the AcShaker family. GO enrichment analysis manifested that the ion-conducting pathway of the AcShaker family plays a crucial role in regulating plant growth and development and adversity stress. Compared with the transcriptome data of the control group, all AcShaker genes were expressed under low-K+stress, and the expression differences of three genes (AcShaker15, AcShaker17, and AcShaker22) were highly significant. The qRT-PCR results showed a high correlation with the transcriptome data, which indicated that these three differentially expressed genes could regulate low-K+ stress and reduce K+ damage in kiwifruit plants, thus improving the resistance to low-K+ stress. Comparison between the salt stress and control transcriptomic data revealed that the expression of AcShaker11 and AcShaker18 genes was significantly different and lower under salt stress, indicating that both genes could be involved in salt stress resistance in kiwifruit. CONCLUSION: The results showed that 28 AcShaker genes were identified and all expressed under low K+ stress, among which AcShaker22 was differentially and significantly upregulated. The AcShaker22 gene can be used as a candidate gene to cultivate new varieties of kiwifruit resistant to low K+ and provide a reference for exploring more properties and functions of the AcShaker genes.

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion.

BACKGROUND: Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions. RESULTS: Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP". CONCLUSION: The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

Haplotype-based pangenomes reveal genetic variations and climate adaptations in moso bamboo populations.

Moso bamboo (Phyllostachys edulis), an ecologically and economically important forest species in East Asia, plays vital roles in carbon sequestration and climate change mitigation. However, intensifying climate change threatens moso bamboo survival. Here we generate high-quality haplotype-based pangenome assemblies for 16 representative moso bamboo accessions and integrated these assemblies with 427 previously resequenced accessions. Characterization of the haplotype-based pangenome reveals extensive genetic variation, predominantly between haplotypes rather than within accessions. Many genes with allele-specific expression patterns are implicated in climate responses. Integrating spatiotemporal climate data reveals more than 1050 variations associated with pivotal climate factors, including temperature and precipitation. Climate-associated variations enable the prediction of increased genetic risk across the northern and western regions of China under future emissions scenarios, underscoring the threats posed by rising temperatures. Our integrated haplotype-based pangenome elucidates moso bamboo's local climate adaptation mechanisms and provides critical genomic resources for addressing intensifying climate pressures on this essential bamboo. More broadly, this study demonstrates the power of long-read sequencing in dissecting adaptive traits in climate-sensitive species, advancing evolutionary knowledge to support conservation.

An Improved Chromosome-scale Genome Assembly and Population Genetics resource for Populus tremula.

Aspen (Populus tremula L.) is a keystone species and a model system for forest tree genomics. We present an updated resource comprising a chromosome-scale assembly, population genetics and genomics data. Using the resource, we explore the genetic basis of natural variation in leaf size and shape, traits with complex genetic architecture. We generated the genome assembly using long-read sequencing, optical and high-density genetic maps. We conducted whole-genome resequencing of the Umeå Aspen (UmAsp) collection. Using the assembly and re-sequencing data from the UmAsp, Swedish Aspen (SwAsp) and Scottish Aspen (ScotAsp) collections we performed genome-wide association analyses (GWAS) using Single Nucleotide Polymorphisms (SNPs) for 26 leaf physiognomy phenotypes. We conducted Assay of Transposase Accessible Chromatin sequencing (ATAC-Seq), identified genomic regions of accessible chromatin, and subset SNPs to these regions, improving the GWAS detection rate. We identified candidate long non-coding RNAs in leaf samples, quantified their expression in an updated co-expression network, and used this to explore the functions of candidate genes identified from the GWAS. A GWAS found SNP associations for seven traits. The associated SNPs were in or near genes annotated with developmental functions, which represent candidates for further study. Of particular interest was a ~177-kbp region harbouring associations with several leaf phenotypes in ScotAsp. We have incorporated the assembly, population genetics, genomics, and GWAS data into the PlantGenIE.org web resource, including updating existing genomics data to the new genome version, to enable easy exploration and visualisation. We provide all raw and processed data to facilitate reuse in future studies.

Genome-Wide Analysis and Expression Profiling of Soybean RbcS Family in Response to Plant Hormones and Functional Identification of GmRbcS8 in Soybean Mosaic Virus.

Rubisco small subunit (RbcS), a core component with crucial effects on the structure and kinetic properties of the Rubisco enzyme, plays an important role in response to plant growth, development, and various stresses. Although Rbcs genes have been characterized in many plants, their muti-functions in soybeans remain elusive. In this study, a total of 11 GmRbcS genes were identified and subsequently divided into three subgroups based on a phylogenetic relationship. The evolutionary analysis revealed that whole-genome duplication has a profound effect on GmRbcSs. The cis-acting elements responsive to plant hormones, development, and stress-related were widely found in the promoter region. Expression patterns based on the RT-qPCR assay exhibited that GmRbcS genes are expressed in multiple tissues, and notably Glyma.19G046600 (GmRbcS8) exhibited the highest expression level compared to other members, especially in leaves. Moreover, differential expressions of GmRbcS genes were found to be significantly regulated by exogenous plant hormones, demonstrating their potential functions in diverse biology processes. Finally, the function of GmRbcS8 in enhancing soybean resistance to soybean mosaic virus (SMV) was further determined through the virus-induced gene silencing (VIGS) assay. All these findings establish a strong basis for further elucidating the biological functions of RbcS genes in soybeans.

Genome-Wide Identification of microRNAs Associated with Starch Biosynthesis and Endosperm Development in Foxtail Millet.

Foxtail millet is one of the oldest crops, and its endosperm contains up to 70% of starch. Grain filling is an important starch accumulation process associated with foxtail millet yield and quality. However, the molecular mechanisms of grain filling in foxtail millet are relatively unclear. Here, we investigate the genes and regulated miRNAs associated with starch synthesis and metabolism in foxtail millet using high-throughput small RNA, mRNA and degradome sequencing. The regulation of starch synthesis and quality is carried out mainly at the 15 DAA to 35 DAA stage during grain filling. The DEGs between waxy and non-waxy foxtail millet were significant, especially for GBSS. Additionally, ptc-miR169i_R+2_1ss21GA, fve-miR396e_L-1R+1, mtr-miR162 and PC-5p-221_23413 regulate the expression of genes associated with the starch synthesis pathway in foxtail millet. This study provides new insights into the molecular mechanisms of starch synthesis and quality formation in foxtail millet.

Genome-Wide Identification of NAC Gene Family Members of Tree Peony (Paeonia suffruticosa Andrews) and Their Expression under Heat and Waterlogging Stress.

An important family of transcription factors (TFs) in plants known as NAC (NAM, ATAF1/2, and CUC2) is crucial for the responses of plants to environmental stressors. In this study, we mined the NAC TF family members of tree peony (Paeonia suffruticosa Andrews) from genome-wide data and analyzed their response to heat and waterlogging stresses in conjunction with transcriptome data. Based on tree peony's genomic information, a total of 48 PsNAC genes were discovered. Based on how similar their protein sequences were, these PsNAC genes were divided into 14 branches. While the gene structures and conserved protein motifs of the PsNAC genes within each branch were largely the same, the cis-acting elements in the promoter region varied significantly. Transcriptome data revealed the presence of five PsNAC genes (PsNAC06, PsNAC23, PsNAC38, PsNAC41, PsNAC47) and one PsNAC gene (PsNAC37) in response to heat and waterlogging stresses, respectively. qRT-PCR analysis reconfirmed the response of these five PsNAC genes to heat stress and one PsNAC gene to waterlogging stress. This study lays a foundation for the study of the functions and regulatory mechanisms of NAC TFs in tree peony. Meanwhile, the NAC TFs of tree peony in response to heat and waterlogging stress were excavated, which is of great significance for the selection and breeding of new tree peony varieties with strong heat and waterlogging tolerance.

Exploring Evolutionary Pathways and Abiotic Stress Responses through Genome-Wide Identification and Analysis of the Alternative Oxidase (AOX) Gene Family in Common Oat (Avena sativa).

The alternative oxidase (AOX), a common terminal oxidase in the electron transfer chain (ETC) of plants, plays a crucial role in stress resilience and plant growth and development. Oat (Avena sativa), an important crop with high nutritional value, has not been comprehensively studied regarding the AsAOX gene family. Therefore, this study explored the responses and potential functions of the AsAOX gene family to various abiotic stresses and their potential evolutionary pathways. Additionally, we conducted a genome-wide analysis to explore the evolutionary conservation and divergence of AOX gene families among three Avena species (Avena sativa, Avena insularis, Avena longiglumis) and four Poaceae species (Avena sativa, Oryza sativa, Triticum aestivum, and Brachypodium distachyon). We identified 12 AsAOX, 9 AiAOX, and 4 AlAOX gene family members. Phylogenetic, motif, domain, gene structure, and selective pressure analyses revealed that most AsAOXs, AiAOXs, and AlAOXs are evolutionarily conserved. We also identified 16 AsAOX segmental duplication pairs, suggesting that segmental duplication may have contributed to the expansion of the AsAOX gene family, potentially preserving these genes through subfunctionalization. Chromosome polyploidization, gene structural variations, and gene fragment recombination likely contributed to the evolution and expansion of the AsAOX gene family as well. Additionally, we hypothesize that AsAOX2 may have potential function in resisting wounding and heat stresses, while AsAOX4 could be specifically involved in mitigating wounding stress. AsAOX11 might contribute to resistance against chromium and waterlogging stresses. AsAOX8 may have potential fuction in mitigating ABA-mediated stress. AsAOX12 and AsAOX5 are most likely to have potential function in mitigating salt and drought stresses, respectively. This study elucidates the potential evolutionary pathways of the AsAOXs gene family, explores their responses and potential functions to various abiotic stresses, identifies potential candidate genes for future functional studies, and facilitates molecular breeding applications in A. sativa.

Genome-Wide Identification of NAC Family Genes and Their Expression Analyses in Response to Osmotic Stress in Cannabis sativa L.

NAC (NAM, ATAF1/2, and CUC2) transcription factors are unique and essential for plant growth and development. Although the NAC gene family has been identified in a wide variety of plants, its chromosomal location and function in Cannabis sativa are still unknown. In this study, a total of 69 putative CsNACs were obtained, and chromosomal location analysis indicated that the CsNAC genes mapped unevenly to 10 chromosomes. Phylogenetic analyses showed that the 69 CsNACs could be divided into six subfamilies. Additionally, the CsNAC genes in group IV-a are specific to Cannabis sativa and contain a relatively large number of exons. Promoter analysis revealed that most CsNAC promoters contained cis-elements related to plant hormones, the light response, and abiotic stress. Furthermore, transcriptome expression profiling revealed that 24 CsNAC genes in two Cannabis sativa cultivars (YM1 and YM7) were significantly differentially expressed under osmotic stress, and these 12 genes presented differential expression patterns across different cultivars according to quantitative real-time PCR (RT-qPCR) analysis. Among these, the genes homologous to the CsNAC18, CsNAC24, and CsNAC61 genes have been proven to be involved in the response to abiotic stress and might be candidate genes for further exploration to determine their functions. The present study provides a comprehensive insight into the sequence characteristics, structural properties, evolutionary relationships, and expression patterns of NAC family genes under osmotic stress in Cannabis sativa and provides a basis for further functional characterization of CsNAC genes under osmotic stress to improve agricultural traits in Cannabis sativa.

Pangenome Identification and Analysis of Terpene Synthase Gene Family Members in Gossypium.

Terpene synthases (TPSs), key gatekeepers in the biosynthesis of herbivore-induced terpenes, are pivotal in the diversity of terpene chemotypes across and within plant species. Here, we constructed a gene-based pangenome of the Gossypium genus by integrating the genomes of 17 diploid and 10 tetraploid species. Within this pangenome, 208 TPS syntelog groups (SGs) were identified, comprising 2 core SGs (TPS5 and TPS42) present in all 27 analyzed genomes, 6 softcore SGs (TPS11, TPS12, TPS13, TPS35, TPS37, and TPS47) found in 25 to 26 genomes, 131 dispensable SGs identified in 2 to 24 genomes, and 69 private SGs exclusive to a single genome. The mutational load analysis of these identified TPS genes across 216 cotton accessions revealed a great number of splicing variants and complex splicing patterns. The nonsynonymous/synonymous Ka/Ks value for all 52 analyzed TPS SGs was less than one, indicating that these genes were subject to purifying selection. Of 208 TPS SGs encompassing 1795 genes, 362 genes derived from 102 SGs were identified as atypical and truncated. The structural analysis of TPS genes revealed that gene truncation is a major mechanism contributing to the formation of atypical genes. An integrated analysis of three RNA-seq datasets from cotton plants subjected to herbivore infestation highlighted nine upregulated TPSs, which included six previously characterized TPSs in G. hirsutum (AD1_TPS10, AD1_TPS12, AD1_TPS40, AD1_TPS42, AD1_TPS89, and AD1_TPS104), two private TPSs (AD1_TPS100 and AD2_TPS125), and one atypical TPS (AD2_TPS41). Also, a TPS-associated coexpression module of eight genes involved in the terpenoid biosynthesis pathway was identified in the transcriptomic data of herbivore-infested G. hirsutum. These findings will help us understand the contributions of TPS family members to interspecific terpene chemotypes within Gossypium and offer valuable resources for breeding insect-resistant cotton cultivars.

Genome-wide identification and expression analysis of the SPL gene family and its response to abiotic stress in barley (Hordeum vulgare L.).

BACKGROUND: Squamosa promoter-binding protein-like (SPL) is a plant-specific transcription factor that is widely involved in the regulation of plant growth and development, including flower and grain development, stress responses, and secondary metabolite synthesis. However, this gene family has not been comprehensively evaluated in barley, the most adaptable cereal crop with a high nutritional value. RESULTS: In this study, a total of 15 HvSPL genes were identified based on the Hordeum vulgare genome. These genes were named HvSPL1 to HvSPL15 based on the chromosomal distribution of the HvSPL genes and were divided into seven groups (I, II, III, V, VI, VII, and VIII) based on the phylogenetic tree analysis. Chromosomal localization revealed one pair of tandem duplicated genes and one pair of segmental duplicated genes. The HvSPL genes exhibited the highest collinearity with the monocotyledonous plant, Zea mays (27 pairs), followed by Oryza sativa (18 pairs), Sorghum bicolor (16 pairs), and Arabidopsis thaliana (3 pairs), and the fewest homologous genes with Solanum lycopersicum (1 pair). The distribution of the HvSPL genes in the evolutionary tree was relatively scattered, and HvSPL proteins tended to cluster with SPL proteins from Z. mays and O. sativa, indicating a close relationship between HvSPL and SPL proteins from monocotyledonous plants. Finally, the spatial and temporal expression patterns of the 14 HvSPL genes from different subfamilies were determined using quantitative real-time polymerase chain reaction (qRT-PCR). Based on the results, the HvSPL gene family exhibited tissue-specific expression and played a regulatory role in grain development and abiotic stress. HvSPL genes are highly expressed in various tissues during seed development. The expression levels of HvSPL genes under the six abiotic stress conditions indicated that many genes responded to stress, especially HvSPL8, which exhibited high expression under multiple stress conditions, thereby warranting further attention. CONCLUSION: In this study, 15 SPL gene family members were identified in the genome of Hordeum vulgare, and the phylogenetic relationships, gene structure, replication events, gene expression, and potential roles of these genes in millet development were studied. Our findings lay the foundation for exploring the HvSPL genes and performing molecular breeding of barley.

Genome-wide association study of fiber quality traits in US upland cotton (Gossypium hirsutum L.).

KEY MESSAGE: A GWAS in an elite diversity panel, evaluated across 10 environments, identified genomic regions regulating six fiber quality traits, facilitating genomics-assisted breeding and gene discovery in upland cotton. In this study, an elite diversity panel of 348 upland cotton accessions was evaluated in 10 environments across the US Cotton Belt and genotyped with the cottonSNP63K array, for a genome-wide association study of six fiber quality traits. All fiber quality traits, upper half mean length (UHML: mm), fiber strength (FS: g tex-1), fiber uniformity (FU: %), fiber elongation (FE: %), micronaire (MIC) and short fiber content (SFC: %), showed high broad-sense heritability (> 60%). All traits except FE showed high genomic heritability. UHML, FS and FU were all positively correlated with each other and negatively correlated with FE, MIC and SFC. GWAS of these six traits identified 380 significant marker-trait associations (MTAs) including 143 MTAs on 30 genomic regions. These 30 genomic regions included MTAs identified in at least three environments, and 23 of them were novel associations. Phenotypic variation explained for the MTAs in these 30 genomic regions ranged from 6.68 to 11.42%. Most of the fiber quality-associated genomic regions were mapped in the D-subgenome. Further, this study confirmed the pleiotropic region on chromosome D11 (UHML, FS and FU) and identified novel co-localized regions on D04 (FU, SFC), D05 (UHML, FU, and D06 UHML, FU). Marker haplotype analysis identified superior combinations of fiber quality-associated genomic regions with high trait values (UHML = 32.34 mm; FS = 32.73 g tex-1; FE = 6.75%). Genomic analyses of traits, haplotype combinations and candidate gene information described in the current study could help leverage genetic diversity for targeted genetic improvement and gene discovery for fiber quality traits in cotton.

Orthogonal genome editing and transcriptional activation in tomato using CRISPR-Combo systems.

KEY MESSAGE: The CRISPR-Combo systems (Cas9-Combo and CBE-Combo) are designed for comprehensive genetic manipulation, enabling Cas9-based targeted mutagenesis or cytosine base editing with simultaneous gene activation in tomato stable lines.

Genomes of diverse Actinidia species provide insights into cis-regulatory motifs and genes associated with critical traits.

BACKGROUND: Kiwifruit, belonging to the genus Actinidia, represents a unique fruit crop characterized by its modern cultivars being genetically diverse and exhibiting remarkable variations in morphological traits and adaptability to harsh environments. However, the genetic mechanisms underlying such morphological diversity remain largely elusive. RESULTS: We report the high-quality genomes of five Actinidia species, including Actinidia longicarpa, A. macrosperma, A. polygama, A. reticulata, and A. rufa. Through comparative genomics analyses, we identified three whole genome duplication events shared by the Actinidia genus and uncovered rapidly evolving gene families implicated in the development of characteristic kiwifruit traits, including vitamin C (VC) content and fruit hairiness. A range of structural variations were identified, potentially contributing to the phenotypic diversity in kiwifruit. Notably, phylogenomic analyses revealed 76 cis-regulatory elements within the Actinidia genus, predominantly associated with stress responses, metabolic processes, and development. Among these, five motifs did not exhibit similarity to known plant motifs, suggesting the presence of possible novel cis-regulatory elements in kiwifruit. Construction of a pan-genome encompassing the nine Actinidia species facilitated the identification of gene DTZ79_23g14810 specific to species exhibiting extraordinarily high VC content. Expression of DTZ79_23g14810 is significantly correlated with the dynamics of VC concentration, and its overexpression in the transgenic roots of kiwifruit plants resulted in increased VC content. CONCLUSIONS: Collectively, the genomes and pan-genome of diverse Actinidia species not only enhance our understanding of fruit development but also provide a valuable genomic resource for facilitating the genome-based breeding of kiwifruit.

Genome-wide identification of CaWD40 proteins reveals the involvement of a novel complex (CaAN1-CaDYT1-CaWD40-91) in anthocyanin biosynthesis and genic male sterility in Capsicum annuum.

BACKGROUND: The WD40 domain, one of the most abundant in eukaryotic genomes, is widely involved in plant growth and development, secondary metabolic biosynthesis, and mediating responses to biotic and abiotic stresses. WD40 repeat (WD40) protein has been systematically studied in several model plants but has not been reported in the Capsicum annuum (pepper) genome. RESULTS: Herein, 269, 237, and 257 CaWD40 genes were identified in the Zunla, CM334, and Zhangshugang genomes, respectively. CaWD40 sequences from the Zunla genome were selected for subsequent analysis, including chromosomal localization, phylogenetic relationships, sequence characteristics, motif compositions, and expression profiling. CaWD40 proteins were unevenly distributed on 12 chromosomes, encompassing 19 tandem duplicate gene pairs. The 269 CaWD40s were divided into six main branches (A to F) with 17 different types of domain distribution. The CaWD40 gene family exhibited diverse expression patterns, and several genes were specifically expressed in flowers and seeds. Yeast two-hybrid (Y2H) and dual-luciferase assay indicated that CaWD40-91 could interact with CaAN1 and CaDYT1, suggesting its involvement in anthocyanin biosynthesis and male sterility in pepper. CONCLUSIONS: In summary, we systematically characterized the phylogeny, classification, structure, and expression of the CaWD40 gene family in pepper. Our findings provide a valuable foundation for further functional investigations on WD40 genes in pepper.

First insight into the genomes of the Pulmonaria officinalis group (Boraginaceae) provided by repeatome analysis and comparative karyotyping.

BACKGROUND: The genus Pulmonaria (Boraginaceae) represents a taxonomically complex group of species in which morphological similarity contrasts with striking karyological variation. The presence of different numbers of chromosomes in the diploid state suggests multiple hybridization/polyploidization events followed by chromosome rearrangements (dysploidy). Unfortunately, the phylogenetic relationships and evolution of the genome, have not yet been elucidated. Our study focused on the P. officinalis group, the most widespread species complex, which includes two morphologically similar species that differ in chromosome number, i.e. P. obscura (2n = 14) and P. officinalis (2n = 16). Ornamental cultivars, morphologically similar to P. officinalis (garden escapes), whose origin is unclear, were also studied. Here, we present a pilot study on genome size and repeatome dynamics of these closely related species in order to gain new information on their genome and chromosome structure. RESULTS: Flow cytometry confirmed a significant difference in genome size between P. obscura and P. officinalis, corresponding to the number of chromosomes. Genome-wide repeatome analysis performed on genome skimming data showed that retrotransposons were the most abundant repeat type, with a higher proportion of Ty3/Gypsy elements, mainly represented by the Tekay lineage. Comparative analysis revealed no species-specific retrotransposons or striking differences in their copy number between the species. A new set of chromosome-specific cytogenetic markers, represented by satellite DNAs, showed that the chromosome structure in P. officinalis was more variable compared to that of P. obscura. Comparative karyotyping supported the hybrid origin of putative hybrids with 2n = 15 collected from a mixed population of both species and outlined the origin of ornamental garden escapes, presumably derived from the P. officinalis complex. CONCLUSIONS: Large-scale genome size analysis and repeatome characterization of the two morphologically similar species of the P. officinalis group improved our knowledge of the genome dynamics and differences in the karyotype structure. A new set of chromosome-specific cytogenetic landmarks was identified and used to reveal the origin of putative hybrids and ornamental cultivars morphologically similar to P. officinalis.

Genetic resources of African mahogany in Brazil: genomic diversity and structure of forest plantations.

BACKGROUND: African mahogany species (Khaya sp.) have been introduced to Brazil gaining increasing economic interest over the last years, as they produce high quality wood for industrial applications. To this date, however, the knowledge available on the genetic basis of African mahogany plantations in Brazil is limited, which has driven this study to examine the extent of genetic diversity and structure of three cultivated species (Khaya grandifoliola, Khaya senegalensis and Khaya ivorensis) and their prospects for forest breeding. RESULTS: In total, 115 individuals were genotyped (48 of K. grandifoliola, 34 of K. senegalensis and 33 of K. ivorensis) for 3,330 filtered neutral loci obtained from genotyping-by-sequencing for the three species. The number of SNPs varied from 2,951 in K. ivorensis to 4,754 in K. senegalensis. Multiloci clustering, principal component analysis, Bayesian structure and network analyses showed a clear genetic separation among the three species. Structure analysis also showed internal structure within each species, highlighting genetic subgroups that could be sampled for selecting distinct genotypes for further breeding, although the genetic distances are moderate to low. CONCLUSION: In our study, SNP markers efficiently assessed the genomic diversity of African mahogany forest plantations in Brazil. Our genetic data clearly separated the three Khaya species. Moreover, pairwise estimates of genetic distances among individuals within each species showed considerable genetic divergence among individuals. By genotyping 115 pre-selected individuals with desirable growth traits, allowed us not only to recommend superior genotypes but also to identify genetically distinct individuals for use in breeding crosses.

Genome-wide association analysis identifies candidates of three bulb traits in garlic.

Garlic bulbs generally possess several swelling cloves, and the swelling degree of the bulbs determines its yield and appearance quality. However, the genetic basis underlying bulb traits remains poorly known. To address this issue, we performed a genome-wide association analysis for three bulb traits: bulb weight, diameter, and height. It resulted in the identification of 51 significant associated signals from 38 genomic regions. Twelve genes from the associated regions, whose transcript abundances in the developmental bulb showed significant correlations with the investigated traits in 81 garlic accessions, were considered the candidates of the corresponding locus. We focused on five of these candidates and their variations and revealed that the promoter variations of fructose-bisphosphate aldolase-encoding Asa8G05696.1 and beta-fructofuranosidase-encoding Asa6G01167.1 are responsible for the functional diversity of these two genes in garlic population. Interestingly, our results revealed that all candidates we focused on experienced a degree of selection during garlic evolutionary history, and different genotypes of them were retained in two China-cultivated garlic groups. Taken together, these results suggest a potential involvement of those candidates in the parallel evolution of garlic bulb organs in two China-cultivated garlic groups. This study provides important insights into the genetic basis of garlic bulb traits and their evolution.