Genome-Wide Identification and Expression Pattern Profiling of the Aquaporin Gene Family in Papaya (Carica papaya L.)

Aquaporins (AQPs) are mainly responsible for the transportation of water and other small molecules such as CO2 and H2O2, and they perform diverse functions in plant growth, in development, and under stress conditions. They are also active participants in cell signal transduction in plants. However, little is known about AQP diversity, biological functions, and protein characteristics in papaya. To better understand the structure and function of CpAQPs in papaya, a total of 29 CpAQPs were identified and classified into five subfamilies. Analysis of gene structure and conserved motifs revealed that CpAQPs exhibited a degree of conservation, with some differentiation among subfamilies. The predicted interaction network showed that the PIP subfamily had the strongest protein interactions within the subfamily, while the SIP subfamily showed extensive interaction with members of the PIP, TIP, NIP, and XIP subfamilies. Furthermore, the analysis of CpAQPs' promoters revealed a large number of cis-elements participating in light, hormone, and stress responses. CpAQPs exhibited different expression patterns in various tissues and under different stress conditions. Collectively, these results provided a foundation for further functional investigations of CpAQPs in ripening, as well as leaf, flower, fruit, and seed development. They also shed light on the potential roles of CpAQP genes in response to environmental factors, offering valuable insights into their biological functions in papaya.

Chromosome-scale genome assembly of sweet tea (Lithocarpus polystachyus Rehder)

Lithocarpus, with >320 species, is the second largest genus of Fagaceae. However, the lack of a reference genome limits the molecular biology and functional study of Lithocarpus species. Here, we report the chromosome-scale genome assembly of sweet tea (Lithocarpus polystachyus Rehder), the first Lithocarpus species to be sequenced to date. Sweet tea has a 952-Mb genome, with a 21.4-Mb contig N50 value and 98.6% complete BUSCO score. In addition, the per-base consensus accuracy and completeness of the genome were estimated at 60.6 and 81.4, respectively. Genome annotation predicted 37,396 protein-coding genes, with repetitive sequences accounting for 64.2% of the genome. The genome did not undergo whole-genome duplication after the gamma (γ) hexaploidy event. Phylogenetic analysis showed that sweet tea diverged from the genus Quercus approximately at 59 million years ago. The high-quality genome assembly and gene annotation resources enrich the genomics of sweet tea, and will facilitate functional genomic studies in sweet tea and other Fagaceae species.

Frequent allopolyploidy with distant progenitors in the moss genera Physcomitrium and Entosthodon (Funariaceae) identified via subgenome phasing of targeted nuclear genes

Allopolyploids represent a new frontier in species discovery among embryophytes. Within mosses, allopolyploid discovery is challenged by low morphological complexity. The rapid expansion of sequencing approaches in addition to computational developments to identifying genome merger and whole-genome duplication using variation among nuclear loci representing homeologs has allowed for increased allopolyploid discovery among mosses. Here, we test a novel approach to phasing homeologs within loci and phasing loci across subgenomes, or subgenome assignment, called Homologizer, in the family Funariaceae. We confirm the intergeneric hybrid nature of Entosthodon hungaricus, and the allopolyploid origin of Physcomitrium eurystomum and one population of Physcomitrium collenchymatum. We also reveal that hybridization gave rise to Physcomitrium immersum, as well as to yet unrecognized lineages sharing the phenotype of Physcomitrium pyriforme and Physcomitrium sphaericum. Our findings demonstrate the utility of our approach when working with polyploid genomes, and its value in identifying progenitor species using target capture data.

A genome assembly for Orinus kokonorica provides insights into the origin, adaptive evolution and further diversification of two closely related grass genera

Ancient whole-genome duplication (WGD) or polyploidization is prevalent in plants and has played a crucial role in plant adaptation. However, the underlying genomic basis of ecological adaptation and subsequent diversification after WGD are still poorly understood in most plants. Here, we report a chromosome-scale genome assembly for the genus Orinus (Orinus kokonorica as representative) and preform comparative genomics with its closely related genus Cleistogenes (Cleistogenes songorica as representative), both belonging to a newly named subtribe Orininae of the grass subfamily Chloridoideae. The two genera may share one paleo-allotetraploidy event before 10 million years ago, and the two subgenomes of O. kokonorica display neither fractionation bias nor global homoeolog expression dominance. We find substantial genome rearrangements and extensive structural variations (SVs) between the two species. With comparative transcriptomics, we demonstrate that functional innovations of orthologous genes may have played an important role in promoting adaptive evolution and diversification of the two genera after polyploidization. In addition, copy number variations and extensive SVs between orthologs of flower and rhizome related genes may contribute to the morphological differences between the two genera. Our results provide new insights into the adaptive evolution and subsequent diversification of the two genera after polyploidization.

Distribution and Location of BEVs in Different Genotypes of Bananas Reveal the Coevolution of BSVs and Bananas

Members of the family Caulimoviridae contain abundant endogenous pararetroviral sequences (EPRVs) integrated into the host genome. Banana streak virus (BSV), a member of the genus Badnavirus in this family, has two distinct badnaviral integrated sequences, endogenous BSV (eBSV) and banana endogenous badnavirus sequences (BEVs). BEVs are distributed widely across the genomes of different genotypes of bananas. To clarify the distribution and location of BEVs in different genotypes of bananas and their coevolutionary relationship with bananas and BSVs, BEVs and BSVs were identified in 102 collected banana samples, and a total of 327 BEVs were obtained and categorized into 26 BEVs species with different detection rates. However, the majority of BEVs were found in Clade II, and a few were clustered in Clade I. Additionally, BEVs and BSVs shared five common conserved motifs. However, BEVs had two unique amino acids, methionine and lysine, which differed from BSVs. BEVs were distributed unequally on most of chromosomes and formed hotspots. Interestingly, a colinear relationship of BEVs was found between AA and BB, as well as AA and SS genotypes of bananas. Notably, the chromosome integration time of different BEVs varied. Based on our findings, we propose that the coevolution of bananas and BSVs is driven by BSV Driving Force (BDF), a complex interaction between BSVs, eBSVs, and BEVs. This study provides the first clarification of the relationship between BEVs and the coevolution of BSVs and bananas in China.

Toward Transgene-Free Transposon-Mediated Biological Mutagenesis for Plant Breeding

Genetic diversity is a key factor for plant breeding. The birth of novel genic and genomic variants is also crucial for plant adaptation in nature. Therefore, the genomes of almost all living organisms possess natural mutagenic mechanisms. Transposable elements (TEs) are a major mutagenic force driving genetic diversity in wild plants and modern crops. The relatively rare TE transposition activity during the thousand-year crop domestication process has led to the phenotypic diversity of many cultivated species. The utilization of TE mutagenesis by artificial and transient acceleration of their activity in a controlled mode is an attractive foundation for a novel type of mutagenesis called TE-mediated biological mutagenesis. Here, I focus on TEs as mutagenic sources for plant breeding and discuss existing and emerging transgene-free approaches for TE activation in plants. Furthermore, I also review the non-randomness of TE insertions in a plant genome and the molecular and epigenetic factors involved in shaping TE insertion preferences. Additionally, I discuss the molecular mechanisms that prevent TE transpositions in germline plant cells (e.g., meiocytes, pollen, egg and embryo cells, and shoot apical meristem), thereby reducing the chances of TE insertion inheritance. Knowledge of these mechanisms can expand the TE activation toolbox using novel gene targeting approaches. Finally, the challenges and future perspectives of plant populations with induced novel TE insertions (iTE plant collections) are discussed.

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.

CRISPR/Cas9-mediated seamless gene replacement in protoplasts expands the resistance spectrum to TMV-U1 strain in regenerated Nicotiana tabacum

CRISPR/Cas-based genome editing is now extensively used in plant breeding and continues to evolve. Most CRISPR/Cas current applications in plants focus on gene knock-outs; however, there is a pressing need for new methods to achieve more efficient delivery of CRISPR components and gene knock-ins to improve agronomic traits of crop cultivars. We report here a genome editing system that combines the advantages of protoplast technologies with recent CRISPR/Cas advances to achieve seamless large fragment insertions in the model Solanaceae plant Nicotiana tabacum. With this system, two resistance-related regions of the N' gene were replaced with homologous fragments from the N'alata gene to confer TMV-U1 resistance in the T0 generation of GMO-free plants. Our study establishes a reliable genome-editing tool for efficient gene modifications and provides a detailed description of the optimization process to assist other researchers adapt this system for their needs.

Genome evolution in plants and the origins of innovation

Plant evolution has been characterised by a series of major novelties in their vegetative and reproductive traits that have led to greater complexity. Underpinning this diversification has been the evolution of the genome. When viewed at the scale of the plant kingdom, plant genome evolution has been punctuated by conspicuous instances of gene and whole-genome duplication, horizontal gene transfer and extensive gene loss. The periods of dynamic genome evolution often coincide with the evolution of key traits, demonstrating the coevolution of plant genomes and phenotypes at a macroevolutionary scale. Conventionally, plant complexity and diversity have been considered through the lens of gene duplication and the role of gene loss in plant evolution remains comparatively unexplored. However, in light of reductive evolution across multiple plant lineages, the association between gene loss and plant phenotypic diversity warrants greater attention.

Genome-wide association study of soluble solids content, flesh color, and fruit shape in citron watermelon

Fruit quality traits are crucial determinants of consumers' willingness to purchase watermelon produce, making them major goals for breeding programs. There is limited information on the genetic underpinnings of fruit quality traits in watermelon. A total of 125 citron watermelon (Citrullus amarus) accessions were genotyped using single nucleotide polymorphisms (SNPs) molecular markers generated via whole-genome resequencing. A total of 2,126,759 genome-wide SNP markers were used to uncover marker-trait associations using single and multi-locus GWAS models. High broad-sense heritability for fruit quality traits was detected. Correlation analysis among traits revealed positive relationships, with the exception of fruit diameter and fruit shape index (ratio of fruit length to fruit diameter), which was negative. A total of 37 significant SNP markers associated with soluble solids content, flesh color, fruit length, fruit diameter, and fruit shape index traits were uncovered. These peak SNPs accounted for 2.1%-23.4% of the phenotypic variation explained showing the quantitative inheritance nature of the evaluated traits. Candidate genes relevant to fruit quality traits were uncovered on chromosomes Ca01, Ca03, Ca06, and Ca07. These significant molecular markers and candidate genes will be useful in marker-assisted breeding of fruit quality traits in watermelon.

Comprehensive Evolution and Expression anaLysis of PHOSPHATE 1 Gene Family in Allotetraploid Brassica napus and Its Diploid Ancestors

The members of PHOSPHATE 1 (PHO1) family play important roles in plant phosphate (Pi) transport and adaptation to Pi deficiency. The functions of PHO1 family proteins have been reported in several plant species, with the exception of Brassica species. Here, we identified 23, 23, and 44 putative PHO1 family genes in Brassica rapa, Brassica oleracea, and Brassica napus by whole genome analysis, respectively. The phylogenetic analysis divided PHO1 family proteins into eight groups, which represented the orthologous relationships among PHO1 members. The gene structure and the conserved motif analysis indicated that the most PHO1 family genes had similar gene structures and the PHO1 proteins shared mutual conserved motifs. The chromosome distribution analysis showed that the majority of BnPHO1 family genes distributed analogously at chromosomes with BrPHO1 and BoPHO1 family genes. The data showed that PHO1 family genes were highly conserved during evolution from diploid to tetraploid. Furthermore, the expression analysis showed that PHO1 family genes had different expression patterns in plant tissues, suggesting the diversity of gene functions in Brassica species. Meanwhile, the expression analysis also revealed that some PHO1 family genes were significantly responsive to Pi deficiency, suggesting that PHO1 family genes play critical roles in Pi uptake and homeostasis under low Pi stress. Altogether, the characteristics of PHO1 family genes provide a reliable groundwork for further dissecting their functions in Brassica species.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Plant trait relationships are maintained within a major crop species: lack of artificial selection signal and potential for improved agronomic performance

The exploration of phenotypic spaces of large sets of plant species has considerably increased our understanding of diversification processes in the plant kingdom. Nevertheless, such advances have predominantly relied on interspecific comparisons that hold several limitations. Here, we grew in the field a unique set of 179 inbred lines of durum wheat, Triticum turgidum spp. durum, characterized by variable degrees of artificial selection. We measured aboveground and belowground traits as well as agronomic traits to explore the functional and agronomic trait spaces and to investigate trait-to-agronomic performance relationships. We showed that the wheat functional trait space shared commonalities with global cross-species spaces previously described, with two main axes of variation: a root foraging axis and a slow-fast trade-off axis. Moreover, we detected a clear signature of artificial selection on the variation of agronomic traits, unlike functional traits. Interestingly, we identified alternative phenotypic combinations that can optimize crop performance. Our work brings insightful knowledge about the structure of phenotypic spaces of domesticated plants and the maintenance of phenotypic trade-offs in response to artificial selection, with implications for trade-off-free and multi-criteria selection in plant breeding.

Developing an efficient and visible prime editing system to restore tobacco 8-hydroxy-copalyl diphosphate gene for labdane diterpene Z-abienol biosynthesis

Prime editing (PE) is a versatile CRISPR-Cas based precise genome-editing platform widely used to introduce a range of possible base conversions in various organisms. However, no PE systems have been shown to induce heritable mutations in tobacco, nor in any other dicot. In this study, we generated an efficient PE system in tobacco that not only introduced heritable mutations, but also enabled anthocyanin-based reporter selection of transgene-free T1 plants. This system was used to confer Z-abienol biosynthesis in the allotetraploid tobacco cultivar HHDJY by restoring a G>T conversion in the NtCPS2 gene. High levels of Z-abienol were detected in the leaves of homozygous T1 plants at two weeks after topping. This study describes an advance in PE systems and expands genome-editing toolbox in tobacco, even in dicots, for use in basic research and molecular breeding. And restoring biosynthesis of Z-abienol in tobacco might provide an efficient way to obtain Z-abienol in plants.

Genome-wide analysis of ALDH gene family in jujube and identification of ZjALDH3F3 for its important role in high-temperature tolerance

Aldehyde dehydrogenases (ALDHs) are NAD(P)-dependent enzymes that oxidize aliphatic and aromatic aldehydes. They play crucial roles in various biological processes and plant responses to stress. The impact of high temperatures on jujube quality and yield has been well documented. Nevertheless, the involvement of ALDHs in the response to heat stress remains poorly understood. This study aimed to identify ZjALDHs in the jujube genome (Ziziphus jujuba var. spinosa) and conducted in silico analyses. Phylogenetic analyses indicated that ALDHs in plants, including jujube, can be divided into ten families, and members from the same family share conserved gene and protein structures. Quantitative real-time PCR (qRT-PCR) and β-glucuronidase (GUS) histochemical staining were used to analyze the expression patterns of ZjALDHs in response to elevated temperatures. We identified a ZjALDH (ZjALDH3F3) gene displaying a significant upregulation and down-regulation, respectively in heat-resistant (HR) and heat-sensitive (HS) jujube in response to heat treatments. Such specific responses are probably attributed to the different heat-responsive cis-elements of ZjALDH3F3 in HR and HS jujubes. ZjALDH3F3 over-expressed in tobacco increased heat tolerance, as evidenced by the reduced accumulation of reactive oxygen species (ROS) and elevated activity of antioxidant enzymes. The qRT-PCR results indicated that the expression of antioxidant enzymes, abscisic acid (ABA), and stress-responsive genes was enhanced in transgenic tobacco. This study sheds novel light on the function of ZjALDHs in heat resistance of jujube.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Genome-Wide Identification and Comprehensive Analysis of the FtsH Gene Family in Soybean (Glycine max)

The filamentation temperature-sensitive H (FtsH) gene family is critical in regulating plant chloroplast development and photosynthesis. It plays a vital role in plant growth, development, and stress response. Although FtsH genes have been identified in a wide range of plants, there is no detailed study of the FtsH gene family in soybean (Glycine max). Here, we identified 34 GmFtsH genes, which could be categorized into eight groups, and GmFtsH genes in the same group had similar structures and conserved protein motifs. We also performed intraspecific and interspecific collinearity analysis and found that the GmFtsH family has large-scale gene duplication and is more closely related to Arabidopsis thaliana. Cis-acting elements analysis in the promoter region of the GmFtsH genes revealed that most genes contain developmental and stress response elements. Expression patterns based on transcriptome data and real-time reverse transcription quantitative PCR (qRT-PCR) showed that most of the GmFtsH genes were expressed at the highest levels in leaves. Then, GO enrichment analysis indicated that GmFtsH genes might function as a protein hydrolase. In addition, the GmFtsH13 protein was confirmed to be localized in chloroplasts by a transient expression experiment in tobacco. Taken together, the results of this study lay the foundation for the functional determination of GmFtsH genes and help researchers further understand the regulatory network in soybean leaf development.

Genome-wide identification of bZIP gene family and expression analysis of BhbZIP58 under heat stress in wax gourd

BACKGROUND

The basic leucine zipper (bZIP) transcription factor family is one of the most abundant and evolutionarily conserved gene families in plants. It assumes crucial functions in the life cycle of plants, including pathogen defense, secondary metabolism, stress response, seed maturation, and flower development. Although the genome of wax gourd has been published, little is known about the functions, evolutionary background, and gene expression patterns of the bZIP gene family, which limits its utilization.

RESULTS

A total of 61 bZIP genes (BhbZIPs) were identified from wax gourd (Benincasa hispida) genome and divided into 12 subgroups. Whole-genome duplication (WGD) and dispersed duplication (DSD) were the main driving forces of bZIP gene family expansion in wax gourd, and this family may have undergone intense purifying selection pressure during the evolutionary process. We selected BhbZIP58, only one in the member of subgroup B, to study its expression patterns under different stresses, including heat, salt, drought, cold stress, and ABA treatment. Surprisingly, BhbZIP58 had a dramatic response under heat stress. BhbZIP58 showed the highest expression level in the root compared with leaves, stem, stamen, pistil, and ovary. In addition, BhbZIP58 protein was located in the nucleus and had transcriptional activation activity. Overexpression of BhbZIP58 in Arabidopsis enhanced their heat tolerance.

CONCLUSIONS

In this study, bZIP gene family is systematically bioinformatically in wax gourd for the first time. Particularly, BhbZIP58 may have an important role in heat stress. It will facilitate further research on the bZIP gene family regarding their evolutionary history and biological functions.

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

CRISPR/Cas as a Genome-Editing Technique in Fruit Tree Breeding

CRISPR (short for "Clustered Regularly Interspaced Short Palindromic Repeats") is a technology that research scientists use to selectively modify the DNA of living organisms. CRISPR was adapted for use in the laboratory from the naturally occurring genome-editing systems found in bacteria. In this work, we reviewed the methods used to introduce CRISPR/Cas-mediated genome editing into fruit species, as well as the impacts of the application of this technology to activate and knock out target genes in different fruit tree species, including on tree development, yield, fruit quality, and tolerance to biotic and abiotic stresses. The application of this gene-editing technology could allow the development of new generations of fruit crops with improved traits by targeting different genetic segments or even could facilitate the introduction of traits into elite cultivars without changing other traits. However, currently, the scarcity of efficient regeneration and transformation protocols in some species, the fact that many of those procedures are genotype-dependent, and the convenience of segregating the transgenic parts of the CRISPR system represent the main handicaps limiting the potential of genetic editing techniques for fruit trees. Finally, the latest news on the legislation and regulations about the use of plants modified using CRISPR/Cas systems has been also discussed.

Genome-Wide Identification and Characterization of CDPK Gene Family in Cultivated Peanut (Arachis hypogaea L.) Reveal Their Potential Roles in Response to Ca Deficiency

This study identified 45 calcium-dependent protein kinase (CDPK) genes in cultivated peanut (Arachis hypogaea L.), which are integral in plant growth, development, and stress responses. These genes, classified into four subgroups based on phylogenetic relationships, are unevenly distributed across all twenty peanut chromosomes. The analysis of the genetic structure of AhCDPKs revealed significant similarity within subgroups, with their expansion primarily driven by whole-genome duplications. The upstream promoter sequences of AhCDPK genes contained 46 cis-acting regulatory elements, associated with various plant responses. Additionally, 13 microRNAs were identified that target 21 AhCDPK genes, suggesting potential post-transcriptional regulation. AhCDPK proteins interacted with respiratory burst oxidase homologs, suggesting their involvement in redox signaling. Gene ontology and KEGG enrichment analyses affirmed AhCDPK genes' roles in calcium ion binding, protein kinase activity, and environmental adaptation. RNA-seq data revealed diverse expression patterns under different stress conditions. Importantly, 26 AhCDPK genes were significantly induced when exposed to Ca deficiency during the pod stage. During the seedling stage, four AhCDPKs (AhCDPK2/-25/-28/-45) in roots peaked after three hours, suggesting early signaling roles in pod Ca nutrition. These findings provide insights into the roles of CDPK genes in plant development and stress responses, offering potential candidates for predicting calcium levels in peanut seeds.

CoreSNP: an efficient pipeline for core marker profile selection from genome-wide SNP datasets in crops

BACKGROUND

DNA marker profiles play a crucial role in the identification and registration of germplasm, as well as in the distinctness, uniformity, and stability (DUS) testing of new plant variety protection. However, selecting minimal marker sets from large-scale SNP dataset can be challenging to distinguish a maximum number of samples.

RESULTS

Here, we developed the CoreSNP pipeline using a "divide and conquer" strategy and a "greedy" algorithm. The pipeline offers adjustable parameters to guarantee the distinction of each sample pair with at least two markers. Additionally, it allows datasets with missing loci as input. The pipeline was tested in barley, soybean, wheat, rice and maize. A few dozen of core SNPs were efficiently selected in different crops with SNP array, GBS, and WGS dataset, which can differentiate thousands of individual samples. The core SNPs were distributed across all chromosomes, exhibiting lower pairwise linkage disequilibrium (LD) and higher polymorphism information content (PIC) and minor allele frequencies (MAF). It was shown that both the genetic diversity of the population and the characteristics of the original dataset can significantly influence the number of core markers. In addition, the core SNPs capture a certain level of the original population structure.

CONCLUSIONS

CoreSNP is an efficiency way of core marker sets selection based on Genome-wide SNP datasets of crops. Combined with low-density SNP chip or genotyping technologies, it can be a cost-effective way to simplify and expedite the evaluation of genetic resources and differentiate different crop varieties. This tool is expected to have great application prospects in the rapid comparison of germplasm and intellectual property protection of new varieties.

Genome-Wide Identification and Analysis of the ABCF Gene Family in Triticum aestivum

The ATP-binding cassette (ABC) superfamily of proteins is a group of evolutionarily conserved proteins. The ABCF subfamily is involved in ribosomal synthesis, antibiotic resistance, and transcriptional regulation. However, few studies have investigated the role of ABCF in wheat (Triticum aestivum) immunity. Here, we identified 18 TaABCFs and classified them into four categories based on their domain characteristics. Functional similarity between Arabidopsis and wheat ABCF genes was predicted using phylogenetic analysis. A comprehensive genome-wide analysis of gene structure, protein motifs, chromosomal location, and cis-acting elements was also performed. Tissue-specific analysis and expression profiling under temperature, hormonal, and viral stresses were performed using real-time quantitative reverse transcription polymerase chain reaction after randomly selecting one gene from each group. The results revealed that all TaABCF genes had the highest expression at 25 °C and responded to methyl jasmonate induction. Notably, TaABCF2 was highly expressed in all tissues except the roots, and silencing it significantly increased the accumulation of Chinese wheat mosaic virus or wheat yellow mosaic virus in wheat leaves. These results indicated that TaABCF may function in response to viral infection, laying the foundation for further studies on the mechanisms of this protein family in plant defence.

Genetic Analysis Based on CRISPR/Cas9 Technology in Plants

Genome-editing technology is a type of genetic engineering in which DNA is inserted into, replaced in, or deleted from the genome using artificially engineered nucleases or genetic scissors [...].

A chromosome-level genome assembly of the Knoxia roxburghii (Rubiaceae)

Knoxia roxburghii is a well-known medicinal plant that is widely distributed in southern China and Southeast Asia. Its dried roots, known as hongdaji in traditional Chinese medicine, are used to treat a range of diseases, including cancers, carbuncles, and ascites. In this study, we report a de novo chromosome-level genome sequence for this diploid plant, which has a length of approximately 446.30 Mb with a contig N50 size of 42.26 Mb and scaffold N50 size of 44.38 Mb. Approximately 99.78% of the assembled sequences were anchored to 10 pseudochromosomes and 3 gapless assembled chromosomes were included in this assembly. A total of 24,507 genes were annotated, along with 68.92% of repetitive elements. Overall, our results will facilitate further active component biosynthesis for K. roxburghii and provide insights for future functional genomic studies and DNA-informed breeding.

Bioinformatic Assessment and Expression Profiles of the AP2/ERF Superfamily in the Melastoma dodecandrum Genome

AP2/ERF transcription factors play crucial roles in various biological activities, including plant growth, development, and responses to biotic and abiotic stressors. However, limited research has been conducted on the AP2/ERF genes of Melastoma dodecandrum for breeding of this potential fruit crop. Leveraging the recently published whole genome sequence, we conducted a comprehensive assessment of this superfamily and explored the expression patterns of AP2/ERF genes at a genome-wide level. A significant number of genes, totaling 218, were discovered to possess the AP2 domain sequence and displayed notable structural variations among five subfamilies. An uneven distribution of these genes was observed on 12 pseudochromosomes as the result of gene expansion facilitated by segmental duplications. Analysis of cis-acting elements within promoter sites and 87.6% miRNA splicing genes predicted their involvement in multiple hormone responses and abiotic stresses through transcriptional and post-transcriptional regulations. Transcriptome analysis combined with qRT-PCR results indicated that certain candidate genes are involved in tissue formation and the response to developmental changes induced by IAA hormones. Overall, our study provides valuable insights into the evolution of ERF genes in angiosperms and lays a solid foundation for future breeding investigations aimed at improving fruit quality and enhancing adaptation to barren land environments.

CRISPR Variants for Gene Editing in Plants: Biosafety Risks and Future Directions

The CRISPR genome editing technology is a crucial tool for enabling revolutionary advancements in plant genetic improvement. This review shows the latest developments in CRISPR/Cas9 genome editing system variants, discussing their benefits and limitations for plant improvement. While this technology presents immense opportunities for plant breeding, it also raises serious biosafety concerns that require careful consideration, including potential off-target effects and the unintended transfer of modified genes to other organisms. This paper highlights strategies to mitigate biosafety risks and explores innovative plant gene editing detection methods. Our review investigates the international biosafety guidelines for gene-edited crops, analyzing their broad implications for agricultural and biotechnology research and advancement. We hope to provide illuminating and refined perspectives for industry practitioners and policymakers by evaluating CRISPR genome enhancement in plants.

Genome-wide identification, classification, and characterization of lectin gene superfamily in sweet orange (Citrus sinensis L.)

Lectins are sugar-binding proteins found abundantly in plants. Lectin superfamily members have diverse roles, including plant growth, development, cellular processes, stress responses, and defense against microbes. However, the genome-wide identification and functional analysis of lectin genes in sweet orange (Citrus sinensis L.) remain unexplored. Therefore, we used integrated bioinformatics approaches (IBA) for in-depth genome-wide identification, characterization, and regulatory factor analysis of sweet orange lectin genes. Through genome-wide comparative analysis, we identified a total of 141 lectin genes distributed across 10 distinct gene families such as 68 CsB-Lectin, 13 CsLysin Motif (LysM), 4 CsChitin-Bind1, 1 CsLec-C, 3 CsGal-B, 1 CsCalreticulin, 3 CsJacalin, 13 CsPhloem, 11 CsGal-Lec, and 24 CsLectinlegB.This classification relied on characteristic domain and phylogenetic analysis, showing significant homology with Arabidopsis thaliana's lectin gene families. A thorough analysis unveiled common similarities within specific groups and notable variations across different protein groups. Gene Ontology (GO) enrichment analysis highlighted the predicted genes' roles in diverse cellular components, metabolic processes, and stress-related regulation. Additionally, network analysis of lectin genes with transcription factors (TFs) identified pivotal regulators like ERF, MYB, NAC, WRKY, bHLH, bZIP, and TCP. The cis-acting regulatory elements (CAREs) found in sweet orange lectin genes showed their roles in crucial pathways, including light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and more. These findings will aid in the in-depth molecular examination of these potential genes and their regulatory elements, contributing to targeted enhancements of sweet orange species in breeding programs.

CRISPR-broad: combined design of multi-targeting gRNAs and broad, multiplex target finding

In CRISPR-Cas and related nuclease-mediated genome editing, target recognition is based on guide RNAs (gRNAs) that are complementary to selected DNA regions. While single site targeting is fundamental for localized genome editing, targeting to expanded and multiple chromosome elements is desirable for various biological applications such as genome mapping and epigenome editing that make use of different fusion proteins with enzymatically dead Cas9. The current gRNA design tools are not suitable for this task, as these are optimized for defining single gRNAs for unique loci. Here, we introduce CRISPR-broad, a standalone, open-source application that defines gRNAs with multiple but specific targets in large continuous or spread regions of the genome, as defined by the user. This ability to identify multi-targeting gRNAs and corresponding multiple targetable regions in genomes is based on a novel aggregate gRNA scoring derived from on-target windows and off-target sites. Applying the new tool to the genomes of two model species, C. elegans and H. sapiens, we verified its efficiency in determining multi-targeting gRNAs and ranking potential target regions optimized for broad targeting. Further, we demonstrated the general usability of CRISPR-broad by cellular mapping of a large human genome element using dCas9 fused to green fluorescent protein.

A chromosome-level genome assembly of Korean mint (Agastache rugosa)

Agastache rugosa, also known as Korean mint, is a perennial plant from the Lamiaceae family that is traditionally used for various ailments and contains antioxidant and antibacterial phenolic compounds. Molecular breeding of A. rugosa can enhance secondary metabolite production and improve agricultural traits, but progress in this field has been delayed due to the lack of chromosome-scale genome information. Herein, we constructed a chromosome-level reference genome using Nanopore sequencing and Hi-C technology, resulting in a final genome assembly with a scaffold N50 of 52.15 Mbp and a total size of 410.67 Mbp. Nine pseudochromosomes accounted for 89.1% of the predicted genome. The BUSCO analysis indicated a high level of completeness in the assembly. Repeat annotation revealed 561,061 repeat elements, accounting for 61.65% of the genome, with Copia and Gypsy long terminal repeats being the most abundant. A total of 26,430 protein-coding genes were predicted, with an average length of 1,184 bp. The availability of this chromosome-scale genome will advance our understanding of A. rugosa's genetic makeup and its potential applications in various industries.

PlantIntronDB: a database for plant introns that host functional elements

Although more and more attention has been focused on introns and the important role of plant introns in plant growth and development has been discovered, there is still a lack of an open and comprehensive database on plant introns with functional elements in current research. In order to make full use of large-scale sequencing data and help researchers in related fields to achieve high-throughput functional verification of identified plant introns with functional elements, we designed a database containing five plant species, PlantIntronDB and systematically analyzed 358, 59, 185, 210 and 141 RNA-seq samples from Arabidopsis thaliana (Arabidopsis), Gossypium raimondii (cotton), Zea mays (maize), Brassica napus (oilseed rape) and Oryza sativa Japonica Group (rice). In total, we found 100 126 introns that host functional elements in these five species. Specifically, we found that among all species, the number of introns with functional elements on the positive and negative strands is similar, with a length mostly smaller than 1500 bp, and the Adenine/Thymine (A/T) content is much higher than that of Guanine/Cytosine (G/C). In addition, the distribution of functional elements in introns varies among different species. All the above data can be downloaded for free in this database. This database provides a concise, comprehensive and user-friendly web interface, allowing users to easily retrieve target data based on their needs, using relevant organizational options. The database operation is simple and convenient, aiming to provide strong data support for researchers in related fields to study plant introns that host functional elements, including circular RNAs, lncRNAs, etc. Database URL:  http://deepbiology.cn/PlantIntronDB/.

Genome-wide identification and expression analysis of the GASA gene family in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

The Gibberellic Acid-Stimulated Arabidopsis (GASA) gene family is widely involved in the regulation of plant growth, development, and stress response. However, information on the GASA gene family has not been reported in Chinese cabbage (Brassica rapa L. ssp. pekinensis).

RESULTS

Here, we conducted genome-wide identification and analysis of the GASA genes in Chinese cabbage. In total, 15 GASA genes were identified in the Chinese cabbage genome, and the physicochemical property, subcellular location, and tertiary structure of the corresponding GASA proteins were elucidated. Phylogenetic analysis, conserved motif, and gene structure showed that the GASA proteins were divided into three well-conserved subfamilies. Synteny analysis proposed that the expansion of the GASA genes was influenced mainly by whole-genome duplication (WGD) and transposed duplication (TRD) and that duplication gene pairs were under negative selection. Cis-acting elements of the GASA promoters were involved in plant development, hormonal and stress responses. Expression profile analysis showed that the GASA genes were widely expressed in different tissues of Chinese cabbage, but their expression patterns appeared to diverse. The qRT-PCR analysis of nine GASA genes confirmed that they responded to salt stress, heat stress, and hormonal triggers.

CONCLUSIONS

Overall, this study provides a theoretical basis for further exploring the important role of the GASA gene family in the functional genome of Chinese cabbage.

Chromosome-level genome assembly and annotation of Zicaitai (Brassica rapa var. purpuraria)

Zicaitai is a seasonal vegetable known for its high anthocyanin content in both stalks and leaves, yet its reference genome has not been published to date. Here, we generated the first chromosome-level genome assembly of Zicaitai using a combination of PacBio long-reads, Illumina short-reads, and Hi-C sequencing techniques. The final genome length is 474.12 Mb with a scaffold N50 length of 43.82 Mb, a BUSCO score of 99.30% and the LAI score of 10.14. Repetitive elements accounted for 60.89% (288.72 Mb) of the genome, and Hi-C data enabled the allocation of 430.87 Mb of genome sequences to ten pseudochromosomes. A total of 42,051 protein-coding genes were successfully predicted using multiple methods, of which 99.74% were functionally annotated. Notably, comparing the genome of Zicaitai with seven other species in the Cruciferae family revealed strong conservation in terms of gene numbers and structures. Overall, the high-quality genome assembly provides a critical resource for studying the genetic basis of important agronomic traits in Zicaitai.

Genome-Wide Analysis of bHLH Family Genes and Identification of Members Associated with Cold/Drought-Induced Photoinhibition in Kandelia obovata

Plant basic helix-loop-helix (bHLH) transcription factors play pivotal roles in responding to stress, including cold and drought. However, it remains unclear how bHLH family genes respond to these stresses in Kandelia obovata. In this study, we identified 75 bHLH members in K. obovata, classified into 11 subfamilies and unevenly distributed across its 18 chromosomes. Collineation analysis revealed that segmental duplication primarily drove the expansion of KobHLH genes. The KobHLH promoters were enriched with elements associated with light response. Through RNA-seq, we identified several cold/drought-associated KobHLH genes. This correlated with decreased net photosynthetic rates (Pn) in the leaves of cold/drought-treated plants. Weighted gene co-expression network analysis (WGCNA) confirmed that 11 KobHLH genes were closely linked to photoinhibition in photosystem II (PS II). Among them, four Phytochrome Interacting Factors (PIFs) involved in chlorophyll metabolism were significantly down-regulated. Subcellular localization showed that KobHLH52 and KobHLH30 were located in the nucleus. Overall, we have comprehensively analyzed the KobHLH family and identified several members associated with photoinhibition under cold or drought stress, which may be helpfulfor further cold/drought-tolerance enhancement and molecular breeding through genetic engineering in K. obovata.

A comprehensive lettuce variation map reveals the impact of structural variations in agronomic traits

BACKGROUND

As an important vegetable crop, cultivated lettuce is grown worldwide and a great variety of agronomic traits have been preserved within germplasm collections. The mechanisms underlying these phenotypic variations remain to be elucidated in association with sequence variations. Compared with single nucleotide polymorphisms, structural variations (SVs) that have more impacts on gene functions remain largely uncharacterized in the lettuce genome.

RESULTS

Here, we produced a comprehensive SV set for 333 wild and cultivated lettuce accessions. Comparison of SV frequencies showed that the SVs prevalent in L. sativa affected the genes enriched in carbohydrate derivative catabolic and secondary metabolic processes. Genome-wide association analysis of seven agronomic traits uncovered potentially causal SVs associated with seed coat color and leaf anthocyanin content.

CONCLUSION

Our work characterized a great abundance of SVs in the lettuce genome, and provides a valuable genomic resource for future lettuce breeding.

Targeted large fragment deletion in plants using paired crRNAs with type I CRISPR system

The CRISPR-Cas systems have been widely used as genome editing tools, with type II and V systems typically introducing small indels, and type I system mediating long-range deletions. However, the precision of type I systems for large fragment deletion is still remained to be optimized. Here, we developed a compact Cascade-Cas3 Dvu I-C system with Cas11c for plant genome editing. The Dvu I-C system was efficient to introduce controllable large fragment deletion up to at least 20 kb using paired crRNAs. The paired-crRNAs design also improved the controllability of deletions for the type I-E system. Dvu I-C system was sensitive to spacer length and mismatch, which was benefit for target specificity. In addition, we showed that the Dvu I-C system was efficient for generating stable transgenic lines in maize and rice with the editing efficiency up to 86.67%. Overall, Dvu I-C system we developed here is powerful for achieving controllable large fragment deletions.

Genome-wide identification of Glutathione peroxidase (GPX) family genes and silencing TaGPX3.2A reduced disease resistance in wheat

Glutathione peroxidase (GPX) is a crucial enzyme that scavenges reactive oxygen species in plants, playing a vital role in enhancing plant stress resistance. In this study, we identified 14 glutathione peroxidase genes (TaGPXs) from common hexaploid wheat (Triticum aestivum L.). These genes were subsequently categorized into three distinct groups based on their phylogenetic relationships. Simultaneously, a preliminarily analysis was conducted on the protein characteristics, chromosome localization, gene structure, cis-regulatory elements and transcriptome. Using reverse transcription quantitative PCR to analyze the expression patterns of five GPX genes that were investigated under various exogenous hormone treatments. According to the qRT-PCR analysis, it indicated that TaGPX genes have the distinct expression patterns. The enzyme activities in transiently overexpressed Nicotiana benthamiana (TaGPX3.2A and TaGPX3.4A) leaves were measured under salt and drought stresses, showed that peroxidase (POD) exhibited higher enzyme activity under stresses. Silencing TaGPX3.2A by virus-induced gene silencing (VIGS) led to reduced resistance of wheat to Fusarium graminearum, indicating that TaGPX3.2A plays a crucial role in enhancing wheat resistance against F. graminearum. This research provides a foundational basis for further investigations on the functional characterization of TaGPXs family members. And in the future it is provides valuable resources for genetic improvement of wheat resistance.

Genetic analysis of iron, zinc and grain yield in wheat-Aegilops derivatives using multi-locus GWAS

BACKGROUND

Wheat is a major staple crop and helps to reduce worldwide micronutrient deficiency. Investigating the genetics that control the concentrations of iron (Fe) and zinc (Zn) in wheat is crucial. Hence, we undertook a comprehensive study aimed at elucidating the genomic regions linked to the contents of Fe and Zn in the grain.

METHODS AND RESULTS

We performed the multi-locus genome-wide association (ML-GWAS) using a panel of 161 wheat-Aegilops substitution and addition lines to dissect the genomic regions controlling grain iron (GFeC), and grain zinc (GZnC) contents. The wheat panel was genotyped using 10,825 high-quality SNPs and phenotyped in three different environments (E1-E3) during 2017-2019. A total of 111 marker-trait associations (MTAs) (at p-value < 0.001) were detected that belong to all three sub-genomes of wheat. The highest number of MTAs were identified for GFeC (58), followed by GZnC (44) and yield (9). Further, six stable MTAs were identified for these three traits and also two pleiotropic MTAs were identified for GFeC and GZnC. A total of 1291 putative candidate genes (CGs) were also identified for all three traits. These CGs encode a diverse set of proteins, including heavy metal-associated (HMA), bZIP family protein, AP2/ERF, and protein previously associated with GFeC, GZnC, and grain yield.

CONCLUSIONS

The significant MTAs and CGs pinpointed in this current study are poised to play a pivotal role in enhancing both the nutritional quality and yield of wheat, utilizing marker-assisted selection (MAS) techniques.

Genome-wide identification and characterization of microsatellite markers in Bactrian Camel

Simple sequence repeats (SSRs) have been widely used for parentage testing, marker-assisted selection, and evolutionary studies. The insufficient availability of SSR markers in Bactrian camels partially accounts for the lack of systematic breeding. Therefore, we aimed to establish a comprehensive SSR dataset for the Bactrian camel. Our approach involved genome searching to locate every SSR in the genome, SSR-enriched sequencing to acquire polymorphism information, and literature research to collect published data. The resulting dataset contains 213,711 SSRs and details their characteristics, including genome coordinates, motifs, lengths, annotations, PCR primers, and polymorphism information. The dataset reveals a biased distribution of SSRs in the Bactrian camel genome, reflecting the mutation mechanism and complex evolution of SSRs. In practice, we successfully demonstrated the utility of the dataset through parentage testing using 15 randomly selected SSRs. This comprehensive dataset can facilitate systematic breeding and enable QTL mapping and GWAS of the Bactrian camel.

Targeted knockout of a conserved plant mitochondrial gene by genome editing

Fusion proteins derived from transcription activator-like effectors (TALEs) have emerged as genome editing tools for mitochondria. TALE nucleases (TALENs) have been applied to delete chimaeric reading frames and duplicated (redundant) genes but produced complex genomic rearrangements due to the absence of non-homologous end-joining. Here we report the targeted deletion of a conserved mitochondrial gene, nad9, encoding a subunit of respiratory complex I. By generating a large number of TALEN-mediated mitochondrial deletion lines, we isolated, in addition to mutants with rearranged genomes, homochondriomic mutants harbouring clean nad9 deletions. Characterization of the knockout plants revealed impaired complex I biogenesis, male sterility and defects in leaf and flower development. We show that these defects can be restored by expressing a functional Nad9 protein from the nuclear genome, thus creating a synthetic cytoplasmic male sterility system. Our data (1) demonstrate the feasibility of using genome editing to study mitochondrial gene functions by reverse genetics, (2) highlight the role of complex I in plant development and (3) provide proof-of-concept for the construction of synthetic cytoplasmic male sterility systems for hybrid breeding by genome editing.

PloiDB: the plant ploidy database

See also the Commentary on this article by Spoelhof et al., 240: 909–911.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus

BACKGROUND

Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks.

RESULTS

In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications.

CONCLUSIONS

Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes' rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.