Genome Duplication Events Research Articles

Dipterocarpoidae genomics reveal their demography and adaptations to Asian rainforests

Dipterocarpoideae species form the emergent layer of Asian rainforests. They are the indicator species for Asian rainforest distribution, but they are severely threatened. Here, to understand their adaptation and population decline, we assemble high-quality genomes of seven Dipterocarpoideae species including two autotetraploid species. We estimate the divergence time between Dipterocarpoideae and Malvaceae and within Dipterocarpoideae to be 108.2 (97.8‒118.2) and 88.4 (77.7‒102.9) million years ago, and we identify a whole genome duplication event preceding dipterocarp lineage diversification. We find several genes that showed a signature of selection, likely associated with the adaptation to Asian rainforests. By resequencing of two endangered species, we detect an expansion of effective population size after the last glacial period and a recent sharp decline coinciding with the history of local human activities. Our findings contribute to understanding the diversification and adaptation of dipterocarps and highlight anthropogenic disturbances as a major factor in their endangered status.

Unifying duplication episode clustering and gene-species mapping inference

We present a novel problem, called MetaEC, which aims to infer gene-species assignments in a collection of partially leaf-labeled gene trees labels by minimizing the size of duplication episode clustering (EC). This problem is particularly relevant in metagenomics, where incomplete data often poses a challenge in the accurate reconstruction of gene histories. To solve MetaEC, we propose a polynomial time dynamic programming (DP) formulation that verifies the existence of a set of duplication episodes from a predefined set of episode candidates. In addition, we design a method to infer distributions of gene-species mappings. We then demonstrate how to use DP to design an algorithm that solves MetaEC. Although the algorithm is exponential in the worst case, we introduce a heuristic modification of the algorithm that provides a solution with the knowledge that it is exact. To evaluate our method, we perform two computational experiments on simulated and empirical data containing whole genome duplication events, showing that our algorithm is able to accurately infer the corresponding events.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae.

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

The functional divergence of homologous GPAT9 genes contributes to the erucic acid content of Brassica napus seeds

BackgroundThe early allopolyploid Brassica napus was a hybrid of two Brassica species, that had undergone a whole genome duplication event followed by genome restructuring, including deletions and small scale duplications. A large number of homologous genes appeared functional divergence during species domestication. Due to the high conservation of de novo glycerolipid biosynthesis, multiple homologues of glycerol-3-phosphate acyltransferases (GPATs) have been found in B. napus. Moreover, the functional variances among these homologous GPAT-encoding genes are unclear.ResultsIn this study, four B. napus homologous genes encoding glycerol-3-phosphate acyltransferase 9 (BnaGPAT9) were characterized. Although a bioinformatics analysis indicated high protein sequence similarity, the homologues demonstrated tissue-specific expression patterns and functional divergence. Yeast genetic complementation assays revealed that BnaGPAT9-A1/C1 homologues but not BnaGPAT9-A10/C9 homologues encoded functional GPAT enzymes. Furthermore, a single nucleotide polymorphism of BnaGPAT9-C1 that occurred during the domestication process was associated with enzyme activity and contributed to the fatty acid composition. The seed-specific expression of BnGPAT9-C11124A increased the erucic acid content in the transformant seeds.ConclusionsThis study revealed that BnaGPAT9 gene homologues evolved into functionally divergent forms with important roles in erucic acid biosynthesis.

Molecular cloning and characterization of BNIP3 and NIX1/2 and their role in DHA-induced mitophagy and apoptosis in grass carp (Ctenopharyngodon idellus) adipocytes

B-cell lymphoma-2 and adenovirus E1B 19-kDa-interacting protein 3 (BNIP3) and BNIP3 like (BNIP3L or NIX) play a vital role in regulating mitophagy and the intrinsic apoptosis in mammals, but their gene characterizations remain unclear in fish. Herein, bnip3, nix1 and nix2 were isolated and characterized from grass carp (Ctenopharyngodon idellus), which encode peptides of 194, 233 and 222 amino acids, respectively. As typical BH3-only proteins, grass carp BNIP3, NIX1 and NIX2 proteins contain BH3 and C-terminal transmembrane domains for inducing apoptosis. Moreover, the LC3-interacting region motif of BNIP3, NIX1 and NIX2 is also conserved in grass carp. Phylogenetic analyses also demonstrated that nix1 and nix2 may have originated from the genome duplication event. Expression pattern analysis indicated that bnip3, nix1 and nix2 were highest expressed in brain, followed by eye (bnip3) and liver (nix1 and nix2). BNIP3, NIX1 and NIX2 localized to the nucleus and the cytoplasm, with a predominant localization to mitochondria within the cytoplasm. In the present study, we found that 200 μM DHA impaired the mitochondrial function, manifested as the decreased antioxidant ability, cellular ATP content and mitochondrial membrane potential in grass carp adipocytes. In addition, the gene expression and enzyme activities of caspase family were significantly increased in 200 μM DHA group, indicating that adipocyte apoptosis was induced. Meanwhile, DHA increased the gene expression of bnip3, nix1 and nix2 in a dose-dependent manner in grass carp adipocytes. The colocalization of mitochondria and lysosomes was promoted by 200 μM DHA treatment, implying that BNIP3/NIX-related mitophagy was activated in adipocytes. Based on these findings, it can be inferred that BNIP3/NIX-related mitophagy may be involved in the adipocyte apoptosis induced by DHA in grass carp.

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

Impact of Whole Genome Duplication Events on the Diversification of Legumes

Impact of Whole Genome Duplication Events on the Diversification of Legumes

Structure-function analyses reveal Arabidopsis thaliana HDA7 to be an inactive histone deacetylase

Histone deacetylases (HDACs), responsible for the removal of acetyl groups from histone tails, are important epigenetic factors. They play a critical role in the regulation of gene expression and are significant in the context of plant growth and development. The Rpd3/Hda1 family of HDACs is reported to regulate key biological processes in plants, such as stress response, seed, embryonic, and floral development. Here, we characterized Arabidopsis thaliana HDA7, a Class I, Rpd3/Hda1 family HDAC. SAXS and AUC results show that the recombinantly expressed and purified histone deacetylase domain of AtHDA7 exists as a monomer in solution. Further, the crystal structure showed AtHDA7 to fold into the typical α/β arginase fold, characteristic of Rpd3/Hda1 family HDACs. Sequence analysis revealed that the Asp and His residues of the catalytic ‘XDXH’ motif present in functional Rpd3/Hda1 family HDACs are mutated to Gly and Pro, respectively, in AtHDA7, suggesting that it might be catalytically inactive. The Asp and His residues are important for Zn2+-binding. Not surprisingly, the crystal structure did not have Zn2+ bound in the catalytic pocket, which is essential for the HDAC activity. Further, our in vitro activity assay revealed AtHDA7 to be inactive as an HDAC. A search in the sequence databases suggested that homologs of AtHDA7 are found exclusively in the Brassicaceae family to which Arabidopsis belongs. It is possible that HDA7 descended from HDA6 through whole genome duplication and triplication events during evolution, as suggested in a previous phylogenetic study.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D.

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B.racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

Comparative genomic analysis of chemosensory-related gene families in gastropods

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann–Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

Abstract The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Genome-wide analysis of SET domain genes and the function of GhSDG51 during salt stress in upland cotton (Gossypium hirsutum L.)

BackgroundCotton, being extensively cultivated, holds immense economic significance as one of the most prominent crops globally. The SET (Su(var), E, and Trithorax) domain-containing protein is of significant importance in plant development, growth, and response to abiotic stress by modifying the lysine methylation status of histone. However, the comprehensive identification of SET domain genes (SDG) have not been conducted in upland cotton (Gossypium hirsutum L.).ResultsA total of 229 SDGs were identified in four Gossypium species, including G. arboretum, G. raimondii, G. hirsutum, and G. barbadense. These genes could distinctly be divided into eight groups. The analysis of gene structure and protein motif revealed a high degree of conservation among the SDGs within the same group. Collinearity analysis suggested that the SDGs of Gossypium species and most of the other selected plants were mainly expanded by dispersed duplication events and whole genome duplication (WGD) events. The allopolyploidization event also has a significant impact on the expansion of SDGs in tetraploid Gossypium species. Furthermore, the characteristics of these genes have been relatively conserved during the evolution. Cis-element analysis revealed that GhSDGs play a role in resistance to abiotic stresses and growth development. Furthermore, the qRT-PCR results have indicated the ability of GhSDGs to respond to salt stress. Co-expression analysis revealed that GhSDG51 might co-express with genes associated with salt stress. In addition, the silencing of GhSDG51 in cotton by the virus-induced gene silencing (VIGS) method suggested a potential positive regulatory role of GhSDG51 in salt stress.ConclusionsThe results of this study comprehensively analyze the SDGs in cotton and provide a basis for understanding the biological role of SDGs in the stress resistance in upland cotton.

Expectations of duplicate gene retention under the gene duplicability hypothesis

BackgroundGene duplication is an important process in evolution. What causes some genes to be retained after duplication and others to be lost is a process not well understood. The most prevalent theory is the gene duplicability hypothesis, that something about the function and number of interacting partners (number of subunits of protein complex, etc.), determines whether copies have more opportunity to be retained for long evolutionary periods. Some genes are also more susceptible to dosage balance effects following WGD events, making them more likely to be retained for longer periods of time. One would expect these processes that affect the retention of duplicate copies to affect the conditional probability ratio after consecutive whole genome duplication events. The probability that a gene will be retained after a second whole genome duplication event (WGD2), given that it was retained after the first whole genome duplication event (WGD1) versus the probability a gene will be retained after WGD2, given it was lost after WGD1 defines the probability ratio that is calculated.ResultsSince duplicate gene retention is a time heterogeneous process, the time between the events (t1) and the time since the most recent event (t2) are relevant factors in calculating the expectation for observation in any genome. Here, we use a survival analysis framework to predict the probability ratio for genomes with different values of t1 and t2 under the gene duplicability hypothesis, that some genes are more susceptible to selectable functional shifts, some more susceptible to dosage compensation, and others only drifting. We also predict the probability ratio with different values of t1 and t2 under the mutational opportunity hypothesis, that probability of retention for certain genes changes in subsequent events depending upon how they were previously retained. These models are nested such that the mutational opportunity model encompasses the gene duplicability model with shifting duplicability over time. Here we present a formalization of the gene duplicability and mutational opportunity hypotheses to characterize evolutionary dynamics and explanatory power in a recently developed statistical framework.ConclusionsThis work presents expectations of the gene duplicability and mutational opportunity hypotheses over time under different sets of assumptions. This expectation will enable formal testing of processes leading to duplicate gene retention.

Comparative genomic analysis of the RabGAP gene family in seven Rosaceae species, and functional identification of PbrRabGAP10 in controlling pollen tube growth by mediating cellulose deposition in pear

Rab GTPase-activating proteins (RabGAPs), serving as crucial signaling switches, play essential roles in several physiological processes related to plant growth and development. However, despite their importance, information regarding the RabGAP gene family and their biological functions remains unknown in the Rosaceae. In this study, we identified a total of 127 RabGAP genes in seven Rosaceae species, which were divided into five subfamilies. Our findings indicate that whole genome duplication (WGD) events or dispersed duplication events largely contributed to the expansion of RabGAP family members within Rosaceae species. Through tissue-specific expression analyses, we revealed that the PbrRabGAP genes exhibited distinct expression patterns in different pear tissues. Furthermore, by examining the expression pattern during pollen development and employing an antisense oligonucleotide approach, we demonstrated that PbrRabGAP10, located in the cytoplasm, mediates the imbalance of cellulose distribution, thus regulating pollen tube elongation. In conclusion, the present study offers an overview of the RabGAP family in Rosaceae genomes and serves as the basis for further functional studies.

Origin and early divergence of tandem duplicated sorbitol transporter genes in Rosaceae: insights from evolutionary analysis of the SOT gene family in angiosperms.

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol transporters (SOTs) play a vital role in facilitating sorbitol allocation from source to sink organs and sugar accumulation in sink organs. While prior research has addressed gene duplications within the SOT gene family in Rosaceae, the precise origin and evolutionary dynamics of these duplications remain unclear, largely due to the complicated interplay of whole genome duplications and tandem duplications. Here, we investigated the synteny relationships among all identified Polyol/Monosaccharide Transporter (PLT) genes in 61 angiosperm genomes and SOT genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of PLTs and SOTs across diverse plant lineages. We found that Rosaceae SOTs, as PLT family members, originated from a pair of tandemly duplicated PLT genes within Class III-A. Furthermore, our investigation highlights the role of lineage-specific and synergistic duplications in Amygdaloideae in contributing to the expansion of SOTs in Rosaceae plants. Collectively, our findings provide insights into the genomic origins, duplication events, and subsequent divergence of SOT gene family members. Such insights lay a crucial foundation for comprehensive functional characterizations in future studies.

The evolution of NLRC3 subfamily genes in Sebastidae teleost fishes

BackgroundWith more than 36,000 valid fish species, teleost fishes constitute the most species-rich vertebrate clade and exhibit extensive genetic and phenotypic variation, including diverse immune defense strategies. NLRC3 subfamily genes, which are specific to fishes, play vital roles in the immune system of teleosts. The evolution of teleosts has been impacted by several whole-genome duplication (WGD) events, which might be a key reason for the expansions of the NLRC3 subfamily, but detailed knowledge of NLRC3 subfamily evolution in the family Sebastidae is still limited.ResultsPhylogenetic inference of NLRC3 subfamily protein sequences were conducted to evaluate the orthology of NLRC3 subfamily genes in black rockfish (Sebastes schlegilii), 13 other fish species from the families Sebastidae, Serranidae, Gasterosteidae and Cyclopteridae, and three species of high vertebrates (bird, reptile and amphibian). WGD analyses were used to estimate expansions and contractions of the NLRC3 subfamily, and patterns of expression of NLRC3 subfamily genes in black rockfish following bacterial infections were used to investigate the functional roles of these genes in the traditional and mucosal immune system of the Sebastidae. Different patterns of gene expansions and contractions were observed in 17 fish and other species examined, and one and two whole-genome duplication events were observed in two members of family Sebastidae (black rockfish and honeycomb rockfish, Sebastes umbrosus), respectively. Subsequently, 179 copy numbers of NLRC3 genes were found in black rockfish and 166 in honeycomb rockfish. Phylogenetic analyses corroborated the conservation and evolution of NLRC3 orthologues between Sebastidae and other fish species. Finally, differential expression analyses provided evidence of the immune roles of NLRC3 genes in black rockfish during bacterial infections and gene ontology analysis also indicated other functional roles.ConclusionsWe hypothesize that NLRC3 genes have evolved a variety of different functions, in addition to their role in the immune response, as a result of whole genome duplication events during teleost diversification. Importantly, this study had underscored the importance of sampling across taxonomic groups, to better understand the evolutionary patterns of the innate immunity system on which complex immunological novelties arose. Moreover, the results in this study could extend current knowledge of the plasticity of the immune system.

Characterization and transcript expression analyses of four Atlantic salmon (Salmo salar) serpinh1 paralogues provide evidence of evolutionary divergence

Atlantic salmon (Salmo salar) are not only the world’s most economically important farmed fish in terms of total value, but also a salmonid, which means that they are invaluable for studies of the evolutionary fate of genes following multiple whole-genome duplication (WGD) events. In this study, four paralogues of the molecular chaperone serpinh1 were characterized in Atlantic salmon, as while this gene is considered to be a sensitive biomarker of heat stress in salmonids, mammalian studies have also identified it as being essential for collagen structural assembly and integrity. The four salmon paralogues were cloned and sequenced so that in silico analyses at the nucleotide and deduced amino acid levels could be performed. In addition, qPCR was used to measure: paralogue- and sex-specific constitutive serpinh1 expression across 17 adult tissues; and their expression in the liver and head kidney of male Atlantic salmon as affected by stress phenotype (high vs. low responder), increased temperature, and injection with a multi-valent vaccine. Compared to the other three paralogues, serpinh1a-2 had a unique constitutive expression profile across the 17 tissues. Although stress phenotype had minimal impact on the transcript expression of the four paralogues, injection with a commercial vaccine containing several formalin inactivated bacterins increased the expression of most paralogues (by 1.1 to 4.5-fold) across both tissues. At 20 °C, the expression levels of serpinh1a-1 and serpinh1a-2 were generally lower (by −1.1- to −1.6-fold), and serpinh1b-1 and serpinh1b-2 were 10.2- to 19.0-fold greater, in comparison to salmon held at 12 °C. With recent studies suggesting a putative link between serpinh1 and upper thermal tolerance in salmonids, the current research is a valuable first step in elucidating the potential mechanisms involved. This research: supports the use of serpinh1b-1 and serpinh1b-2 as a biomarkers of heat stress in salmon; and provides evidence of neo- and/or subfunctionalization between the paralogues, and important insights into how multiple genome duplication events can potentially lead to evolutionary divergence.

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

BackgroundEpigenetic 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.ResultsIn 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.ConclusionsHerein, 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.

Characterization of caspase gene family in Sebastes schlegelii and their expression profiles under Aeromonas salmonicida and Vibrio anguillarum infection

The caspase, functioning as a proteinase, plays a crucial role in eukaryotic cell apoptosis, regulation of apoptosis, cellular growth, differentiation, and immunity. The identification of caspase gene family in Sebastes schlegelii is of great help to understand its antimicrobial research. In S. schlegelii, we totally identified nine caspase genes, including four apoptosis initiator caspases (caspase 2, caspase 8, caspase 9 and caspase 10), four apoptosis executioners (caspase 3a, caspase 3b, caspase 6, and caspase 7) and one inflammatory executioner (caspase 1). The duplication of caspase 3 genes on chr3 and chr8 may have been facilitated by whole genome duplication (WGD) events or other complex evolutionary processes. In general, the number of caspase genes relatively conserved in high vertebrates, while exhibiting variation in teleosts. Furthermore, syntenic analysis and phylogenetic relationships analysis supported the correct classification of these caspase gene family in S. schlegelii, especially for genes with duplicated copies. Additionally, the expression patterns of these caspase genes in different tissues of S. schlegelii under healthy conditions were assessed. The results revealed that the expression levels of most caspase genes were significantly elevated in the intestine, spleen, and liver. To further investigate the potential immune functions of these caspase genes in S. schlegelii, we challenged individuals with A. salmonicida and V. anguillarum, respectively. After infection with A. salmonicida, the expression levels of caspase 1 in the liver and spleen of S. schlegelii remained consistently elevated throughout the infection time points. The expression levels of most caspase family members in the intestine exhibited significant divergence following V. anguillarum infection. This study provides a comprehensive understanding of the caspase gene families in S. schlegelii, thereby establishing a solid foundation for further investigations into the functional roles of these caspase genes.

Systematic analysis of the NDR1/HIN1-like (NHL) family in Gossypium hirsutum reveals a role of GhNHL69 in responding to cold stress

NDR1/HIN1-like (NHL) proteins widely exist in the plant kingdom, and play an essential role in resistance to both biotic and abiotic stress. However, the study of NHL family in cotton is still lacking. In the present study, a total of 40, 46, 79, and 78 NHL genes were detected by genome-wide identification in G.arboreum, G.raimondii, G.barbadense, and G.hirsutum, respectively. By analyzing the duplication pattern, the whole genome duplication (WGD) and allotetraploidization events play a vital role in the expansion of GhNHLs. The analysis of gene structure, conserved protein motif, and Ka/Ks values suggest that the structure and function might remain conserved during the evolution. The cis-acting elements analysis revealed that the GhNHL genes could respond to the abiotic stress, and the further transcriptome analysis and qRT-PCR assay indicated that GhNHLs might play a role in the cold stress tolerance. The WGCNA (Weighted correlation network analysis) analysis showed that GhNHL69 co-express with many transcription factors that are related to cold stress. By performing virus-induced gene silencing (VIGS) assays, the GhNHL69 silencing plant exhibited more severe dehydration and injury compared to the control groups. Additionally, we further found the expression of many abiotic genes was inhibited in the GhNHL69 silenced plants. This study gained insight into NHL members in cotton and mined a gene related to cold stress and lays the important basis for further study on NHL genes related to cold stress tolerance.