Gene duplication as a major force in evolution

Genome-wide identification of the Hsp70 gene family in Penaeus chinensis and their response to environmental stress

The heat shock protein 70 (HSP70) gene family plays a crucial role in the response of organisms to environmental stress. However, it has not been systematically characterized in shrimp. In this study, we identified 25 PcHsp70 genes in the Penaeus chinensis genome. The encoded proteins were categorized into six subgroups based on phylogenetic relationships. Tandem duplication was the main driver of amplification in the PcHsp70 family, and the genes have experienced strong purifying selection during evolution. Transcriptome data analysis revealed that the 25 PcHsp70 members have different expression patterns in shrimp under conditions of low temperature, low salinity, and white spot syndrome virus infection. Among them, PcHsp70.11 was significantly induced under all three stress conditions, suggesting that this gene plays an important role in response to environmental stress in P. chinensis. To the best of our knowledge, this is the first study to systematically analyze the Hsp70 gene family in shrimp. The results provide important information on shrimp Hsp70s, contributing to a better understanding of the role of these genes in environmental stress and providing a basis for further functional studies.

The Molecular Basis of Differentiation Wave Activity in Embryogenesis

As development varies greatly across the tree of life, it may seem difficult to suggest a model that proposes a single mechanism for understanding collective cell behaviors and the coordination of tissue formation. Here we propose a mechanism called differentiation waves, which unify many disparate results involving developmental systems from across the tree of life. We demonstrate how a relatively simple model of differentiation proceeds not from function-related molecular mechanisms, but from so-called differentiation waves. A phenotypic model of differentiation waves is introduced, and its relation to molecular mechanisms is proposed. These waves contribute to a differentiation tree, which is an alternate way of viewing cell lineage and local action of the molecular factors. We construct a model of differentiation wave-related molecular mechanisms (genome, epigenome, and proteome) based on bioinformatic data from the nematode Caenorhabditis elegans. To validate this approach across different modes of development, we evaluate protein expression across different types of development by comparing Caenorhabditis elegans with several model organisms: fruit flies (Drosophila melanogaster), yeast (Saccharomyces cerevisiae), and mouse (Mus musculus). Inspired by gene regulatory networks, two Models of Interactive Contributions (fully-connected MICs and ordered MICs) are used to suggest potential genomic contributions to differentiation wave-related proteins. This, in turn, provides a framework for understanding differentiation and development.

New genes driven by segmental duplications share a testis-specific expression pattern in the chromosome-level genome assembly of tree sparrow

Based on a chromosome-level genome assembly, a burst of new genes with different structures but a similar testis-specific expression pattern was detected in tree sparrow.

Genome-wide characterization and expression profiling of FARL (FHY3/FAR1) family genes in Zea mays

A significant role of the plant is played by the transcription factor FARL, which is light signal transduction as well as plant growth and development. Despite being transposases, FARL has developed a variety of dominant biological actions in evolution and speciation. On the other hand, little is known about the Zea mays FARL protein family. This study identifies and characterizes fifteen ZmFARL genes genome-wide, and RNA sequencing data was used to profile their expression. 105 FARL proteins from five plant species were classified into five groups based on sequence alignment and phylogeny. The ZmFARL genes' exon-intron and motif distribution were conserved based on their evolutionary group. The fifteen ZmFARL genes were distributed over seven of the ten Z. mays chromosomes, although no duplication was discovered. Cis-element analysis reveals that ZmFARL genes play a variety of activities, including tissue-specific, stress- and hormone-responsive expressions. Furthermore, the results of the RNA sequencing used to profile expression showed that the genes ZmFARL2 and ZmFARL5 were much more expressed than other genes in various tissues, particularly in leaf characteristics. The identification of likely genes involved in cellular activity in Z. mays and related species will be aided by the characterization of the FARL genes.

Functional investigation of the two ClpPs and three ClpXs in Myxococcus xanthus DK1622

ClpXP is a protease complex that plays important roles in protein quality control and cell cycle regulation, but the functions of multiple ClpXs and multiple ClpPs in M. xanthus remain unknown. The genome of Myxococcus xanthus DK1622 contains two clpPs and three clpXs. The clpP1 and clpX1 genes are cotranscribed and are both essential, while the other copies are isolated in the genome and are deletable. The deletion of clpX2 caused the mutant to be deficient in fruiting body development, while the clpX3 gene is involved in resistance to thermal stress. Both ClpPs possess catalytic active sites, but only ClpP1 shows in vitro peptidase activity on the typical substrate Suc-LY-AMC. All of these clpP and clpX genes exhibit strong transcriptional upregulation in the stationary phase, and the transcription of the three clpX genes appears to be coordinated. Our results demonstrated that multiple ClpPs and multiple ClpXs are functionally divergent and may assist in the environmental adaptation and functional diversification of M. xanthus.IMPORTANCEClpXP is an important protease complex of bacteria and is involved in various physiological processes. Myxococcus xanthus DK1622 possesses two ClpPs and three ClpXs with unclear functions. We investigated the functions of these genes and demonstrated the essential roles of clpP1 and clpX1. Only ClpP1 has in vitro peptidase activity on Suc-LY-AMC, and the isolated clpX copies participate in distinct cellular processes. All of these genes exhibited significant transcriptional upregulation in the stationary phase. Divergent functions appear in multiple ClpPs and multiple ClpXs in M. xanthus DK1622.

Genome-wide identification of Gα family in grass carp (Ctenopharyngodon idella) and reproductive regulation functional characteristics of Cignaq

BACKGROUND

The Gα family plays a crucial role in the complex reproductive regulatory network of teleosts. However, the characterization and function of Gα family members, especially Gαq, remain poorly understood in teleosts. To analyze the characterization, expression, and function of grass carp (Ctenopharyngodon idella) Gαq, we identified the Gα family members in grass carp genome, and analyzed the expression, distribution, and signal transduction of Gαq/gnaq. We also explored the role of Gαq in the reproductive regulation of grass carp.

RESULTS

Our results showed that the grass carp genome contains 27 Gα genes with 46 isoforms, which are divided into four subfamilies: Gαs, Gαi/o, Gαq/11, and Gα12/13. The expression level of Cignaq in the testis was the highest and significantly higher than in other tissues, followed by the hypothalamus and brain. The luteinizing hormone receptor (LHR) was mainly localized to the nucleus in grass carp oocytes, with signals also present in follicular cells. In contrast, Gαq signal was mainly found in the cytoplasm of oocytes, with no signal in follicular cells. In the testis, Gαq and LHR were co-localized in the cytoplasm. Furthermore, the grass carp Gαq recombinant protein significantly promoted Cipgr expression.

CONCLUSIONS

These results provided preliminary evidence for understanding the role of Gαq in the reproductive regulation of teleosts.

The GC% landscape of the Nucleocytoviricota

Genomic studies on sequence composition employ various approaches, such as calculating the proportion of guanine and cytosine within a given sequence (GC% content), which can shed light on various aspects of the organism's biology. In this context, GC% can provide insights into virus-host relationships and evolution. Here, we present a comprehensive gene-by-gene analysis of 61 representatives belonging to the phylum Nucleocytoviricota, which comprises viruses with the largest genomes known in the virosphere. Parameters were evaluated not only based on the average GC% of a given viral species compared to the entire phylum but also considering gene position and phylogenetic history. Our results reveal that while some families exhibit similar GC% among their representatives (e.g., Marseilleviridae), others such as Poxviridae, Phycodnaviridae, and Mimiviridae have members with discrepant GC% values, likely reflecting adaptation to specific biological cycles and hosts. Interestingly, certain genes located at terminal regions or within specific genomic clusters show GC% values distinct from the average, suggesting recent acquisition or unique evolutionary pressures. Horizontal gene transfer and the presence of potential paralogs were also assessed in genes with the most discrepant GC% values, indicating multiple evolutionary histories. Taken together, to the best of our knowledge, this study represents the first global and gene-by-gene analysis of GC% distribution and profiles within genomes of Nucleocytoviricota members, highlighting their diversity and identifying potential new targets for future studies.

Rp3: Ribosome profiling-assisted proteogenomics improves coverage and confidence during microprotein discovery

There has been a dramatic increase in the identification of non-canonical translation and a significant expansion of the protein-coding genome. Among the strategies used to identify unannotated small Open Reading Frames (smORFs) that encode microproteins, Ribosome profiling (Ribo-Seq) is the gold standard for the annotation of novel coding sequences by reporting on smORF translation. In Ribo-Seq, ribosome-protected footprints (RPFs) that map to multiple genomic sites are removed since they cannot be unambiguously assigned to a specific genomic location. Furthermore, RPFs necessarily result in short (25-34 nucleotides) reads, increasing the chance of multi-mapping alignments, such that smORFs residing in these regions cannot be identified by Ribo-Seq. Moreover, it has been challenging to identify protein evidence for Ribo-Seq. To solve this, we developed Rp3, a pipeline that integrates proteogenomics and Ribosome profiling to provide unambiguous evidence for a subset of microproteins missed by current Ribo-Seq pipelines. Here, we show that Rp3 maximizes proteomics detection and confidence of microprotein-encoding smORFs.

A Rigorous Framework to Classify the Postduplication Fate of Paralogous Genes

Gene duplication has a central role in evolution; still, little is known on the fates of the duplicated copies, their relative frequency, and on how environmental conditions affect them. Moreover, the lack of rigorous definitions concerning the fate of duplicated genes hinders the development of a global vision of this process. In this paper, we present a new framework aiming at characterizing and formally differentiating the fate of duplicated genes. Our framework has been tested via simulations, where the evolution of populations has been simulated using aevol, an in silico experimental evolution platform. Our results show several patterns that confirm some of the conclusions from previous studies, while also exhibiting new tendencies; this may open up new avenues to better understand the role of duplications as a driver of evolution.

Terpene synthases GhTPS6 and GhTPS47 participate in resistance to Verticillium dahliae in upland cotton

Terpene synthases (TPSs) are enzymes responsible for catalyzing the production of diverse terpenes, the largest class of secondary metabolites in plants. Here, we identified 107 TPS gene loci encompassing 92 full-length TPS genes in upland cotton (Gossypium hirsutum L.). Phylogenetic analysis showed they were divided into six subfamilies. Segmental duplication and tandem duplication events contributed greatly to the expansion of TPS gene family, particularly the TPS-a and TPS-b subfamilies. Expression profile analysis screened out that GhTPSs may mediate the interaction between cotton and Verticillium dahliae. Three-dimensional structures and subcellular localizations of the two selected GhTPSs, GhTPS6 and GhTPS47, which belong to the TPS-a subfamily, demonstrated similarity in protein structures and nucleus and cytoplasm localization. Virus-induced gene silencing (VIGS) of the two GhTPSs yielded plants characterized by increased wilting and chlorosis, more severe vascular browning, and higher disease index than control plants. Additionally, knockdown of GhTPS6 and GhTPS47 led to the down-regulation of cotton terpene synthesis following V. dahliae infection, indicating that these two genes may positively regulate resistance to V. dahliae through the modulation of disease-resistant terpene biosynthesis. Overall, our study represents a comprehensive analysis of the G. hirsutum TPS gene family, revealing their potential roles in defense responses against Verticillium wilt.

Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

Link prediction in protein-protein interaction network: A similarity multiplied similarity algorithm with paths of length three

Protein-protein interactions (PPIs) are crucial for various biological processes, and predicting PPIs is a major challenge. To solve this issue, the most common method is link prediction. Currently, the link prediction methods based on network Paths of Length Three (L3) have been proven to be highly effective. In this paper, we propose a novel link prediction algorithm, named SMS, which is based on L3 and protein similarities. We first design a mixed similarity that combines the topological structure and attribute features of nodes. Then, we compute the predicted value by summing the product of all similarities along the L3. Furthermore, we propose the Max Similarity Multiplied Similarity (maxSMS) algorithm from the perspective of maximum impact. Our computational prediction results show that on six datasets, including S. cerevisiae, H. sapiens, and others, the maxSMS algorithm improves the precision of the top 500, area under the precision-recall curve, and normalized discounted cumulative gain by an average of 26.99%, 53.67%, and 6.7%, respectively, compared to other optimal methods.

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

Evolutionary and Expression Analysis of the Pig MAGE Gene Family

The melanoma-associated antigen (MAGE) family found in eukaryotes plays a crucial role in cell proliferation and differentiation, spermatogenesis, neural development, etc. This study explored the validation and evolution of MAGE genes in eukaryotic genomes and their distribution and expression patterns in pigs. In total, 249 MAGE genes were found on 13 eukaryotic species. In total, 33, 25, and 18 genes were located on human, mouse, and pig genomes, respectively. We found eight, four, and three tandemly duplicated gene clusters on the human, mouse, and pig genomes, respectively. The majority of MAGE genes in mammals are located on the X chromosome. According to the phylogenetic analysis, the MAGE family genes were classified into 11 subfamilies. The NDN gene in zebrafish (DreNDN) was the root of this evolutionary tree. In total, 10 and 11 MAGE genes on human and mouse genomes, respectively, exhibited a collinearity relationship with the MAGE genes on pig genomes. Taking the MAGE family genes in pigs, the MAGE subfamilies had similar gene structures, protein motifs, and biochemical attributes. Using the RNA-seq data of Duroc pigs and Rongchang pigs, we detected that the expression of type I MAGE genes was higher in reproductive tissues, but type II MAGE genes were predominantly expressed in the brain tissue. These findings are a valuable resource for gaining insight into the evolution and expression of the MAGE family genes.

Evolutionary Patterns of Maternal Recognition of Pregnancy and Implantation in Eutherian Mammals

The implantation of the embryo into the maternal endometrium is a complex process associated with the evolution of viviparity and placentation in mammals. In this review, we provide an overview of maternal recognition of pregnancy signals and implantation modes in eutherians, focusing on their diverse mechanisms and evolutionary patterns. Different pregnancy recognition signals and implantation modes have evolved in eutherian mammals, reflecting the remarkable diversity of specializations in mammals following the evolution of viviparity. Superficial implantation is the ancestral implantation mode in Eutheria and its major clades. The other modes, secondary, partially, and primary interstitial implantation have each independently evolved multiple times in the evolutionary history of eutherians. Although significant progress has been made in understanding pregnancy recognition signals and implantation modes, there is still much to uncover. Rodents and chiropterans (especially Phyllostomidae) offer valuable opportunities for studying the transitions among implantation modes, but data is still scarce for these diverse orders. Further research should focus on unstudied taxa so we can establish robust patterns of evolutionary changes in pregnancy recognition signaling and implantation modes.

Genome-wide identification of MATE and ALMT genes and their expression profiling in mungbean (Vigna radiata L.) under aluminium stress

The Multidrug and toxic compound extrusion (MATE) and aluminium activated malate transporter (ALMT) gene families are involved in response to aluminium (Al) stress. In this study, we identified 48 MATE and 14 ALMT gene families in Vigna radiata genome and classified into 5 (MATE) and 3 (ALMT) clades by phylogenetic analysis. All the VrMATE and VrALMT genes were distributed across mungbean chromosomes. Tandem duplication was the main driving force for evolution and expansion of MATE gene family. Collinearity of mungbean with soybean indicated that MATE gene family is closely linked to Glycine max. Eight MATE transporters in clade 2 were found to be associated with previously characterized Al tolerance related MATEs in various plant species. Citrate exuding motif (CEM) was present in seven VrMATEs of clade 2. Promoter analysis revealed abundant plant hormone and stress responsive cis-elements. Results from quantitative real time-polymerase chain reaction (qRT-PCR) revealed that VrMATE19, VrMATE30 and VrALMT13 genes were markedly up-regulated at different time points under Al stress. Overall, this study offers a new direction for further molecular characterization of the MATE and ALMT genes in mungbean for Al tolerance.

Genome-Wide Identification and Expression Pattern Analysis of GATA Gene Family in Orchidaceae

The GATA transcription factors play crucial roles in plant growth, development, and responses to environmental stress. Despite extensive studies of GATA genes in many plants, their specific functions and mechanisms in orchids remain unexplored. In our study, a total of 149 GATA genes were identified in the genomes of seven sequenced orchid species (20 PeqGATAs, 23 CgGATAs, 24 CeGATAs, 23 DcaGATAs, 20 DchGATAs, 27 DnoGATAs, and 12 GelGATAs), classified into four subfamilies. Subfamily I typically contains genes with two exons, while subfamily II contains genes with two or three exons. Most members of subfamilies III and IV have seven or eight exons, with longer introns compared to subfamilies I and II. In total, 24 pairs (CgGATAs-DchGATAs), 27 pairs (DchGATAs-DnoGATAs), and 14 pairs (DnoGATAs-GelGATAs) of collinear relationships were identified. Cis-acting elements in GATA promoters were mainly enriched in abscisic acid (ABA) response elements and methyl jasmonate (MeJA) elements. Expression patterns and RT-qPCR analysis revealed that GATAs are involved in the regulation of floral development in orchids. Furthermore, under high-temperature treatment, GL17420 showed an initial increase followed by a decrease, GL18180 and GL17341 exhibited a downregulation followed by upregulation and then a decrease, while GL30286 and GL20810 displayed an initial increase followed by slight inhibition and then another increase, indicating diverse regulatory mechanisms of different GATA genes under heat stress. This study explores the function of GATA genes in orchids, providing a theoretical basis and potential genetic resources for orchid breeding and stress resistance improvement.

Genome-Based Identification of the Dof Gene Family in Three Cymbidium Species and Their Responses to Heat Stress in Cymbidium goeringii

As an important genus in Orchidaceae, Cymbidium has rich ecological diversity and significant economic value. DNA binding with one zinc finger (Dof) proteins are pivotal plant-specific transcription factors that play crucial roles in the growth, development, and stress response of plants. Although the Dof genes have been identified and functionally analyzed in numerous plants, exploration in Orchidaceae remains limited. We conducted a thorough analysis of the Dof gene family in Cymbidium goeringii, C. ensifolium, and C. sinensis. In total, 91 Dof genes (27 CgDofs, 34 CeDofs, 30 CsDofs) were identified, and Dof genes were divided into five groups (I-V) based on phylogenetic analysis. All Dof proteins have motif 1 and motif 2 conserved domains and over half of the genes contained introns. Chromosomal localization and collinearity analysis of Dof genes revealed their evolutionary relationships and potential gene duplication events. Analysis of cis-elements in CgDofs, CeDofs, and CsDofs promoters showed that light-responsive cis-elements were the most common, followed by hormone-responsive elements, plant growth-related elements, and abiotic stress response elements. Dof proteins in three Cymbidium species primarily exhibit a random coil structure, while homology modeling exhibited significant similarity. In addition, RT-qPCR analysis showed that the expression levels of nine CgDofs changed greatly under heat stress. CgDof03, CgDof22, CgDof27, CgDof08, and CgDof23 showed varying degrees of upregulation. Most upregulated genes under heat stress belong to group I, indicating that the Dof genes in group I have great potential for high-temperature resistance. In conclusion, our study systematically demonstrated the molecular characteristics of Dof genes in different Cymbidium species, preliminarily revealed the patterns of heat stress, and provided a reference for further exploration of stress breeding in orchids.

Comprehensive genome-wide identification of Snf2 gene family and their expression profile under salt stress in six Brassica species of U's triangle model

MAIN CONCLUSION

We comprehensively identified and analyzed the Snf2 gene family. Some Snf2 genes were involved in responding to salt stress based on the RNA-seq and qRT-PCR analysis. Sucrose nonfermenting 2 (Snf2) proteins are core components of chromatin remodeling complexes that not only alter DNA accessibility using the energy of ATP hydrolysis, but also play a critical regulatory role in growth, development, and stress response in eukaryotes. However, the comparative study of Snf2 gene family in the six Brassica species in U's triangle model remains unclear. Here, a total of 405 Snf2 genes were identified, comprising 53, 50, and 46 in the diploid progenitors: Brassica rapa (AA, 2n = 20), Brassica nigra (BB, 2n = 16), and Brassica oleracea (CC, 2n = 18), and 93, 91, and 72 in the allotetraploid: Brassica juncea (AABB, 2n = 36), Brassica napus (AACC, 2n = 38), and Brassica carinata (BBCC, 2n = 34), respectively. These genes were classified into six clades and further divided into 18 subfamilies based on their conserved motifs and domains. Intriguingly, these genes showed highly conserved chromosomal distributions and gene structures, indicating that few dynamic changes occurred during the polyploidization. The duplication modes of the six Brassica species were diverse, and the expansion of most Snf2 in Brassica occurred primarily through dispersed duplication (DSD) events. Additionally, the majority of Snf2 genes were under purifying selection during polyploidization, and some Snf2 genes were associated with various abiotic stresses. Both RNA-seq and qRT-PCR analysis showed that the expression of BnaSnf2 genes was significantly induced under salt stress, implying their involvement in salt tolerance response in Brassica species. The results provide a comprehensive understanding of the Snf2 genes in U's triangle model species, which will facilitate further functional analysis of the Snf2 genes in Brassica plants.

Omics approaches to understand the MADS-box gene family in common bean (Phaseolus vulgaris L.) against drought stress

MADS-box genes are known to play important roles in diverse aspects of growth/devolopment and stress response in several plant species. However, no study has yet examined about MADS-box genes in P. vulgaris. In this study, a total of 79 PvMADS genes were identified and classified as type I and type II according to the phylogenetic analysis. While both type I and type II PvMADS classes were found to contain the MADS domain, the K domain was found to be present only in type II PvMADS proteins, in agreement with the literature. All chromosomes of the common bean were discovered to contain PvMADS genes and 17 paralogous gene pairs were identified. Only two of them were tandemly duplicated gene pairs (PvMADS-19/PvMADS-23 and PvMADS-20/PvMADS-24), and the remaining 15 paralogous gene pairs were segmentally duplicated genes. These duplications were found to play an important role in the expansion of type II PvMADS genes. Moreover, the RNAseq and RT-qPCR analyses showed the importance of PvMADS genes in response to drought stress in P. vulgaris.

A novel C-type lectin protein (BjCTL5) interacts with apoptosis stimulating proteins of p53 (ASPP) to activate NF-κB signaling pathway in primitive chordate

C-type lectin proteins (CTLs), a group of pattern recognition receptors (PRRs), play pivotal roles in immune responses. However, the signal transduction and regulation of CTLs in cephalochordates have yet to be explored. In this study, we examined the composition of CTLs in Branchiostoma japonicum, identifying a total of 272 CTLs. These CTLs underwent further analysis concerning domain arrangement, tandem and segmental duplication events. A multidomain C-type lectin gene, designated as BjCTL5, encompassing CLECT, KR, CUB, MAM, and SR domains, was the focal point of our investigation. BjCTL5 exhibits ubiquitous expression across all detected tissues and is responsive to stimulation by LPS, mannose, and poly (I:C). The recombinant protein of BjCTL5 can bind to Escherichia coli and Staphylococcus aureus, inducing their agglutination and inhibiting the proliferation of S. aureus. Yeast two-hybrid, CoIP, and confocal immunofluorescence experiments revealed the interaction between BjCTL5 and apoptosis-stimulating proteins of p53, BjASPP. Intriguingly, BjCTL5 was observed to induce the luciferase activity of the NF-κB promoter in HEK293T cells. These results suggested a potential interaction between BjCTL5 and BjASPP, implicating that they involve in the activation of the NF-κB signaling pathway, which provides an evolutionary viewpoint on NF-κB signaling pathway in primitive chordate.

Genome-scale identification, expression and evolution analysis of B-box members in Dendrobium huoshanense

B-box (BBX) proteins have been recognized as vital determinants in plant development, morphogenesis, and adaptive responses to a myriad of environmental stresses. These zinc-finger proteins play a pivotal role in various biological processes. Their influence spans photomorphogenesis, the regulation of flowering, and imparting resilience to a wide array of challenges, encompassing both biotic and abiotic factors. Chromosome localization, gene structure and conserved motifs, phylogenetic analysis, collinearity analysis, expression profiling, fluorescence quantitative analysis, and tobacco transient transformation methods were used for functional localization and expression pattern analysis of the DhBBX gene. A total of 23 DhBBX members were identified from Dendrobium huoshanense. Subsequent phylogenetic evaluations effectively segregated these genes into five discrete evolutionary subsets. The predictions of subcellular localizations revealed that all these proteins were localized in the nucleus. The genetic composition and patterns showed that the majority of these genes consisted of several exons, with a few variations that could be attributed to transposon insertion. A comprehensive analysis using qRT-PCR was conducted to unravel the expression patterns of these genes in D. huoshanense, with a specific concentration on their responses to various hormone treatments and cold stress. Subcellular localization reveals that DhBBX21 and DhBBX9 are located in the nucleus. Our results provide a deep comprehension of the complex regulatory mechanisms of BBXs in response to various environmental and hormonal stimuli. These discoveries encourage further detailed and focused investigations into the operational dynamics of the BBX gene family in a wider range of plant species.

Copy number variation and population-specific immune genes in the model vertebrate zebrafish

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Genetic characterization, structural analysis, and detection of positive selection in small heat shock proteins of Cypriniformes and Clupeiformes

In the current investigation, a total of 42 full-length, non-redundant small heat shock proteins (sHsp) were detected in Cyprinus carpio, Labeo rohita, Danio rerio, Salmo salar, Oncorhynchus mykiss, and Clupea harengus. The sHsp genes were classified into three groups based on phylogenetic analysis. All the sHsps were shown to have higher aliphatic index values, which is an indication that these proteins are more thermally stable. The hydrophilic nature of sHsps was deduced from the fact that all fish species had negative GRAVY scores. In all of the representative fish species, sHsp genes were assigned to distinct chromosomes in an inconsistent and unequal manner. Segmental duplications are the main events that have contributed to the expansion of the sHsp genes in all species. We were also able to determine the selective pressure that was placed on particular codons and discovered several significant coding sites within the coding region of sHsps. Eventually, diversifying positive selection was found to be connected with evolutionary changes in sHsp proteins, which showed that gene evolution controlled the fish adaption event in response to environmental conditions. Clarification of the links between sHsps and environmental stress in fish will be achieved through rigorous genomic comparison, which will also yield substantial new insights.

Pleiotropy alleviates the fitness costs associated with resource allocation trade-offs in immune signalling networks

Many genes and signalling pathways within plant and animal taxa drive the expression of multiple organismal traits. This form of genetic pleiotropy instigates trade-offs among life-history traits if a mutation in the pleiotropic gene improves the fitness contribution of one trait at the expense of another. Whether or not pleiotropy gives rise to conflict among traits, however, likely depends on the resource costs and timing of trait deployment during organismal development. To investigate factors that could influence the evolutionary maintenance of pleiotropy in gene networks, we developed an agent-based model of co-evolution between parasites and hosts. Hosts comprise signalling networks that must faithfully complete a developmental programme while also defending against parasites, and trait signalling networks could be independent or share a pleiotropic component as they evolved to improve host fitness. We found that hosts with independent developmental and immune networks were significantly more fit than hosts with pleiotropic networks when traits were deployed asynchronously during development. When host genotypes directly competed against each other, however, pleiotropic hosts were victorious regardless of trait synchrony because the pleiotropic networks were more robust to parasite manipulation, potentially explaining the abundance of pleiotropy in immune systems despite its contribution to life history trade-offs.

Genome-wide identification and expression analysis of the SPL transcription factor family and its response to abiotic stress in Pisum sativum L

Squamous promoter binding protein-like (SPL) genes encode plant-specific transcription factors (TFs) that play essential roles in modulating plant growth, development, and stress response. Pea (Pisum sativum L.) is a coarse grain crop of great importance in food production, biodiversity conservation and molecular genetic research, providing genetic information and nutritional resources for improving agricultural production and promoting human health. However, only limited researches on the structure and functions of SPL genes exist in pea (PsSPLs). In this study, we identified 22 PsSPLs and conducted a genome-wide analysis of their physical characteristics, chromosome distribution, gene structure, phylogenetic evolution and gene expression patterns. As a result, the PsSPLs were unevenly distributed on the seven chromosomes of pea and harbored the SBP domain, which is composed of approximately 76 amino acid residues. The phylogenetic analysis revealed that the PsSPLs clustered into eight subfamilies and showed high homology with SPL genes in soybean. Further analysis showed the presence of segmental duplications in the PsSPLs. The expression patterns of 22 PsSPLs at different tissues, developmental stages and under various stimulus conditions were evaluated by qRT-PCR method. It was found that the expression patterns of PsSPLs from the same subfamily were similar in different tissues, the transcripts of most PsSPLs reached the maximum peak value at 14 days after anthesis in the pod. Abiotic stresses can cause significantly up-regulated PsSPL19 expression with spatiotemporal specificity, in addition, four plant hormones can cause the up-regulated expression of most PsSPLs including PsSPL19 in a time-dependent manner. Therefore, PsSPL19 could be a key candidate gene for signal transduction during pea growth and development, pod formation, abiotic stress and plant hormone response. Our findings should provide insights for the elucidating of development regulation mechanism and breeding for resistance to abiotic stress pea.

Jasmonate ZIM Domain Protein (JAZ) Gene SLJAZ15 Increases Resistance to Orobanche aegyptiaca in Tomato

Orobanche aegyptiaca Pers. is a holoparasitic plant that severely reduces tomato (Solanum lycopersicum L.) production in China. However, there is a lack of effective control methods and few known sources of genetic resistance. In this study, we focused on key genes in the JAZ family, comparing the JAZ family in Arabidopsis thaliana (L. Heynh.) to the tomato genome. After identifying the JAZ family members in S. lycopersicum, we performed chromosomal localization and linear analysis with phylogenetic relationship analysis of the JAZ family. We also analyzed the gene structure of the JAZ gene family members in tomato and the homology of the JAZ genes among the different species to study their relatedness. The key genes for O. aegyptiaca resistance were identified using VIGS (virus-induced gene silencing), and the parasitization rate of silenced tomato plants against O. aegyptiaca increased by 47.23-91.13%. The genes were localized in the nucleus by subcellular localization. Heterologous overexpression in A. thaliana showed that the key gene had a strong effect on the parasitization process of O. aegyptiaca, and the overexpression of the key gene reduced the parasitization rate of O. aegyptiaca 1.69-fold. Finally, it was found that the SLJAZ15 gene can positively regulate the hormone content in tomato plants and affect plant growth and development, further elucidating the function of this gene.

Genome-Wide Identification Analysis of GST Gene Family in Wild Blueberry Vaccinium duclouxii and Their Impact on Anthocyanin Accumulation

Vaccinium duclouxii, a wild blueberry species native to the mountainous regions of southwestern China, is notable for its exceptionally high anthocyanin content, surpassing that of many cultivated varieties and offering significant research potential. Glutathione S-transferases (GSTs) are versatile enzymes crucial for anthocyanin transport in plants. Yet, the GST gene family had not been previously identified in V. duclouxii. This study utilized a genome-wide approach to identify and characterize the GST gene family in V. duclouxii, revealing 88 GST genes grouped into seven distinct subfamilies. This number is significantly higher than that found in closely related species, with these genes distributed across 12 chromosomes and exhibiting gene clustering. A total of 46 members are classified as tandem duplicates. The gene structure of VdGST is relatively conserved among related species, showing closer phylogenetic relations to V. bracteatum and evidence of purifying selection. Transcriptomic analysis and qRT-PCR indicated that VdGSTU22 and VdGSTU38 were highly expressed in flowers, VdGSTU29 in leaves, and VdGSTF11 showed significant expression in ripe and fully mature fruits, paralleling trends seen with anthocyanin accumulation. Subcellular localization identified VdGSTF11 primarily in the plasma membrane, suggesting a potential role in anthocyanin accumulation in V. duclouxii fruits. This study provides a foundational basis for further molecular-level functional analysis of the transport and accumulation of anthocyanins in V. duclouxii, enhancing our understanding of the molecular mechanisms underlying anthocyanin metabolism in this valuable species.

Identification and characterization of the plasma membrane H+-ATPase genes in Brassica napus and functional analysis of BnHA9 in salt tolerance

As a primary proton pump, plasma membrane (PM) H+-ATPase plays critical roles in regulating plant growth, development, and stress responses. PM H+-ATPases have been well characterized in many plant species. However, no comprehensive study of PM H+-ATPase genes has been performed in Brassica napus (rapeseed). In this study, we identified 32 PM H+-ATPase genes (BnHAs) in the rapeseed genome, and they were distributed on 16 chromosomes. Phylogenetical and gene duplication analyses showed that the BnHA genes were classified into five subfamilies, and the segmental duplication mainly contributed to the expansion of the rapeseed PM H+-ATPase gene family. The conserved domain and subcellular analyses indicated that BnHAs encoded canonical PM H+-ATPase proteins with 14 highly conserved domains and localized on PM. Cis-acting regulatory element and expression pattern analyses indicated that the expression of BnHAs possessed tissue developmental stage specificity. The 25 upstream open reading frames with the canonical initiation codon ATG were predicted in the 5' untranslated regions of 11 BnHA genes and could be used as potential target sites for improving rapeseed traits. Protein interaction analysis showed that BnBRI1.c associated with BnHA2 and BnHA17, indicating that the conserved activity regulation mechanism of BnHAs may be present in rapeseed. BnHA9 overexpression in Arabidopsis enhanced the salt tolerance of the transgenic plants. Thus, our results lay a foundation for further research exploring the biological functions of PM H+-ATPases in rapeseed.

Genome-Wide Identification, Phylogenetic and Expression Analysis of Expansin Gene Family in Medicago sativa L

Expansins, a class of cell-wall-loosening proteins that regulate plant growth and stress resistance, have been studied in a variety of plant species. However, little is known about the Expansins present in alfalfa (Medicago sativa L.) due to the complexity of its tetraploidy. Based on the alfalfa (cultivar "XinjiangDaye") reference genome, we identified 168 Expansin members (MsEXPs). Phylogenetic analysis showed that MsEXPs consist of four subfamilies: MsEXPAs (123), MsEXPBs (25), MsEXLAs (2), and MsEXLBs (18). MsEXPAs, which account for 73.2% of MsEXPs, and are divided into twelve groups (EXPA-I-EXPA-XII). Of these, EXPA-XI members are specific to Medicago trunctula and alfalfa. Gene composition analysis revealed that the members of each individual subfamily shared a similar structure. Interestingly, about 56.3% of the cis-acting elements were predicted to be associated with abiotic stress, and the majority were MYB- and MYC-binding motifs, accounting for 33.9% and 36.0%, respectively. Our short-term treatment (≤24 h) with NaCl (200 mM) or PEG (polyethylene glycol, 15%) showed that the transcriptional levels of 12 MsEXPs in seedlings were significantly altered at the tested time point(s), indicating that MsEXPs are osmotic-responsive. These findings imply the potential functions of MsEXPs in alfalfa adaptation to high salinity and/or drought. Future studies on MsEXP expression profiles under long-term (>24 h) stress treatment would provide valuable information on their involvement in the response of alfalfa to abiotic stress.

Characterization and expression profiles of the ZmLBD gene family in Zea mays

BACKGROUND

The Lateral Organ Boundaries Domain (LBD) gene family is a family of plant-specific transcription factors (TFs) that are widely involved in processes such as lateral organ formation, stress response, and nutrient metabolism. However, the function of LBD genes in maize remains poorly understood.

METHODS AND RESULTS

In this study, a total of 49 ZmLBD genes were identified at the genome-wide level of maize, they were classified into nine branches based on phylogenetic relationships, and all of them were predicted to be nuclear localized. The 49 ZmLBD genes formed eight pairs of segmental duplicates, and members of the same branches' members had similar gene structure and conserved motif composition. The promoters of ZmLBD genes contain multiple types of cis-acting elements. In addition, by constructing the regulatory network of ZmLBD and other genes and miRNAs, 12 and 22 ZmLBDs were found to be involved in the gene regulatory network and miRNA regulatory network, respectively. The expression pattern analysis suggests that ZmLBD genes may be involved in different biological pathways, and drought stress induced the expressions of two inbred lines.

CONCLUSIONS

The findings enhance our comprehension of the potential roles of the ZmLBD gene family in maize growth and development, which is pivotal for genetic enhancement and breeding efforts pertaining to this significant crop.

Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations

Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43 years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.

Conserved switch genes that arose via whole-genome duplication regulate a cannibalistic nematode morph

Experimental genetics in a nematode reveals a key role for developmental plasticity in the evolution of nutritional diversity.

Genome-Wide Identification of NAC Family Genes in Oat and Functional Characterization of AsNAC109 in Abiotic Stress Tolerance

The plant-specific NAC gene family is one of the largest transcription factor families, participating in plant growth regulation and stress response. Despite extensive characterization in various plants, our knowledge of the NAC family in oat is lacking. Herein, we identified 333 NAC genes from the latest release of the common oat genome. We provide a comprehensive overview of the oat NAC gene family, covering gene structure, chromosomal localization, phylogenetic characteristics, conserved motif compositions, and gene duplications. AsNAC gene expression in different tissues and the response to various abiotic stresses were characterized using RT-qPCR. The main driver of oat NAC gene family expansion was identified as segmental duplication using collinearity analysis. In addition, the functions of AsNAC109 in regulating abiotic stress tolerance in Arabidopsis were clarified. This is the first genome-wide investigation of the NAC gene family in cultivated oat, which provided a unique resource for subsequent research to elucidate the mechanisms responsible for oat stress tolerance and provides valuable clues for the improvement of stress resistance in cultivated oat.

Insights into cucumber (Cucumis sativus) genetics: Genome-wide discovery and computational analysis of the Calreticulin Domain-Encoding gene (CDEG) family

Cucumber is an essential vegetable crop throughout the world. Cucumber development is vital for accomplishing both quality and productivity requirements. Meanwhile, numerous factors have resulted in substantial cucumber losses. However, the calreticulin domain-encoding genes (CDEGs) in cucumber were not well-characterized and had little function. In the genome-wide association study (GWAS), we recognized and characterized the CDEGs in Cucumis sativus (cucumber). Through a comprehensive study of C. sativus, our research has unveiled the presence of three unique genes, denoted as CsCRTb, CsCRT3, and CsCNX1, unevenly distributed on three chromosomes in the genome of C. sativus. In accordance to the phylogenetic investigation, these genes may be categorized into three subfamilies. Based on the resemblance with AtCDE genes, we reorganized the all CsCDE genes in accordance with international nomenclature. The expression analysis and cis-acting components revealed that each of CsCDE gene promoter region enclosed number of cis-elements connected with hormone and stress response. According to subcellular localization studies demonstrated that, they were found in deferent locations of the cell such as endoplasmic reticulum, plasma membrane, golgi apparatus, and vacuole, according to subcellular localization studies. Chromosomal distribution analysis and synteny analysis demonstrated the probability of segmental or tandem duplications within the cucumber CDEG gene family. Additionally, miRNAs displayed diverse modes of action, including mRNA cleavage and translational inhibition. We used the RNA seq data to analyze the expression of CDEG genes in response to cold stress and also improved cold tolerance, which was brought on by treating cucumber plants to an exogenous chitosan oligosaccharide spray. Our investigation revealed that these genes responded to this stress in a variety of ways, demonstrating that they may adapt quickly to environmental changes in cucumber plants. This study provides a base for further understanding in reference to CDE gene family and reveals that genes play significant functions in cucumber stress responses.

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

Genome-wide analysis of the KNOX gene family in Moso bamboo: insights into their role in promoting the rapid shoot growth

BACKGROUND

KNOTTED1-like homeobox (KNOX) genes, plant-specific homologous box transcription factors (TFs), play a central role in regulating plant growth, development, organ formation, and response to biotic and abiotic stresses. However, a comprehensive genome-wide identification of the KNOX genes in Moso bamboo (Phyllostachys edulis), the fastest growing plant, has not yet been conducted, and the specific biological functions of this family remain unknown.

RESULTS

The expression profiles of 24 KNOX genes, divided into two subfamilies, were determined by integrating Moso bamboo genome and its transcriptional data. The KNOX gene promoters were found to contain several light and stress-related cis-acting elements. Synteny analysis revealed stronger similarity with rice KNOX genes than with Arabidopsis KNOX genes. Additionally, several conserved structural domains and motifs were identified in the KNOX proteins. The expansion of the KNOX gene family was primarily regulated by tandem duplications. Furthermore, the KNOX genes were responsive to naphthaleneacetic acid (NAA) and gibberellin (GA) hormones, exhibiting distinct temporal expression patterns in four different organs of Moso bamboo. Short Time-series Expression Miner (STEM) analysis and quantitative real-time PCR (qRT-PCR) assays demonstrated that PeKNOX genes may play a role in promoting rapid shoot growth. Additionally, Gene Ontology (GO) and Protein-Protein Interaction (PPI) network enrichment analyses revealed several functional annotations for PeKNOXs. By regulating downstream target genes, PeKNOXs are involved in the synthesis of AUX /IAA, ultimately affecting cell division and elongation.

CONCLUSIONS

In the present study, we identified and characterized a total of 24 KNOX genes in Moso bamboo and investigated their physiological properties and conserved structural domains. To understand their functional roles, we conducted an analysis of gene expression profiles using STEM and RNA-seq data. This analysis successfully revealed regulatory networks of the KNOX genes, involving both upstream and downstream genes. Furthermore, the KNOX genes are involved in the AUX/IAA metabolic pathway, which accelerates shoot growth by influencing downstream target genes. These results provide a theoretical foundation for studying the molecular mechanisms underlying the rapid growth and establish the groundwork for future research into the functions and transcriptional regulatory networks of the KNOX gene family.

Genome-wide identification and characterization of protein phosphatase 2C (PP2C) gene family in sunflower (Helianthus annuus L.) and their expression profiles in response to multiple abiotic stresses

Plant protein phosphatase 2C (PP2C) plays vital roles in responding to various stresses, stimulating growth factors, phytohormones, and metabolic activities in many important plant species. However, the PP2C gene family has not been investigated in the economically valuable plant species sunflower (Helianthus annuus L.). This study used comprehensive bioinformatics tools to identify and characterize the PP2C gene family members in the sunflower genome (H. annuus r1.2). Additionally, we analyzed the expression profiles of these genes using RNA-seq data under four different stress conditions in both leaf and root tissues. A total of 121 PP2C genes were identified in the sunflower genome distributed unevenly across the 17 chromosomes, all containing the Type-2C phosphatase domain. HanPP2C genes are divided into 15 subgroups (A-L) based on phylogenetic tree analysis. Analyses of conserved domains, gene structures, and motifs revealed higher structural and functional similarities within various subgroups. Gene duplication and collinearity analysis showed that among the 53 HanPP2C gene pairs, 48 demonstrated segmental duplications under strong purifying selection pressure, with only five gene pairs showing tandem duplications. The abundant segmental duplication was observed compared to tandem duplication, which was the major factor underlying the dispersion of the PP2C gene family in sunflowers. Most HanPP2C proteins were localized in the nucleus, cytoplasm, and chloroplast. Among the 121 HanPP2C genes, we identified 71 miRNAs targeting 86 HanPP2C genes involved in plant developmental processes and response to abiotic stresses. By analyzing cis-elements, we identified 63 cis-regulatory elements in the promoter regions of HanPP2C genes associated with light responsiveness, tissue-specificity, phytohormone, and stress responses. Based on RNA-seq data from two sunflower tissues (leaf and root), 47 HanPP2C genes exhibited varying expression levels in leaf tissue, while 49 HanPP2C genes showed differential expression patterns in root tissue across all stress conditions. Transcriptome profiling revealed that nine HanPP2C genes (HanPP2C12, HanPP2C36, HanPP2C38, HanPP2C47, HanPP2C48, HanPP2C53, HanPP2C54, HanPP2C59, and HanPP2C73) exhibited higher expression in leaf tissue, and five HanPP2C genes (HanPP2C13, HanPP2C47, HanPP2C48, HanPP2C54, and HanPP2C95) showed enhanced expression in root tissue in response to the four stress treatments, compared to the control conditions. These results suggest that these HanPP2C genes may be potential candidates for conferring tolerance to multiple stresses and further detailed characterization to elucidate their functions. From these candidates, 3D structures were predicted for six HanPP2C proteins (HanPP2C47, HanPP2C48, HanPP2C53, HanPP2C54, HanPP2C59, and HanPP2C73), which provided satisfactory models. Our findings provide valuable insights into the PP2C gene family in the sunflower genome, which could play a crucial role in responding to various stresses. This information can be exploited in sunflower breeding programs to develop improved cultivars with increased abiotic stress tolerance.

Enzymes Involved in Antioxidant and Detoxification Processes Present Changes in the Expression Levels of Their Coding Genes under the Stress Caused by the Presence of Antimony in Tomato

Currently, there is an increasing presence of heavy metals and metalloids in soils and water due to anthropogenic activities. However, the biggest problem caused by this increase is the difficulty in recycling these elements and their high permanence in soils. There are plants with great capacity to assimilate these elements or make them less accessible to other organisms. We analyzed the behavior of Solanum lycopersicum L., a crop with great agronomic interest, under the stress caused by antimony (Sb). We evaluated the antioxidant response throughout different exposure times to the metalloid. Our results showed that the enzymes involved in the AsA-GSH cycle show changes in their expression level under the stress caused by Sb but could not find a relationship between the NITROSOGLUTATHIONE REDUCTASE (GSNOR) expression data and nitric oxide (NO) content in tomato roots exposed to Sb. We hypothesize that a better understanding of how these enzymes work could be key to develop more tolerant varieties to this kind of abiotic stress and could explain a greater or lesser phytoremediation capacity. Moreover, we deepened our knowledge about Glutathione S-transferase (GST) and Glutathione Reductase (GR) due to their involvement in the elimination of the xenobiotic component.

Mitogenomic architecture and evolution of the soil ciliates Colpoda

Colpoda are cosmopolitan unicellular eukaryotes primarily inhabiting soil and benefiting plant growth, but they remain one of the least understood taxa in genetics and genomics within the realm of ciliated protozoa. Here, we investigate the architecture of de novo assembled mitogenomes of six Colpoda species, using long-read sequencing and involving 36 newly isolated natural strains in total. The mitogenome sizes span from 43 to 63 kbp and typically contain 28-33 protein-coding genes. They possess a linear structure with variable telomeres and central repeats, with one Colpoda elliotti strain isolated from Tibet harboring the longest telomeres among all studied ciliates. Phylogenomic analyses reveal that Colpoda species started to diverge more than 326 million years ago, eventually evolving into two distinct groups. Collinearity analyses also reveal significant genomic divergences and a lack of long collinear blocks. One of the most notable features is the exceptionally high level of gene rearrangements between mitochondrial genomes of different Colpoda species, dominated by gene loss events. Population-level mitogenomic analysis on natural strains also demonstrates high sequence divergence, regardless of geographic distance, but the gene order remains highly conserved within species, offering a new species identification criterion for Colpoda species. Furthermore, we identified underlying heteroplasmic sites in the majority of strains of three Colpoda species, albeit without a discernible recombination signal to account for this heteroplasmy. This comprehensive study systematically unveils the mitogenomic structure and evolution of these ancient and ecologically significant Colpoda ciliates, thus laying the groundwork for a deeper understanding of the evolution of unicellular eukaryotes.IMPORTANCEColpoda, one of the most widespread ciliated protozoa in soil, are poorly understood in regard to their genetics and evolution. Our research revealed extreme mitochondrial gene rearrangements dominated by gene loss events, potentially leading to the streamlining of Colpoda mitogenomes. Surprisingly, while interspecific rearrangements abound, our population-level mitogenomic study revealed a conserved gene order within species, offering a potential new identification criterion. Phylogenomic analysis traced their lineage over 326 million years, revealing two distinct groups. Substantial genomic divergence might be associated with the lack of extended collinear blocks and relaxed purifying selection. This study systematically reveals Colpoda ciliate mitogenome structures and evolution, providing insights into the survival and evolution of these vital soil microorganisms.

ZjHXK5 and ZjHXK6 negatively regulate the sugar metabolism of Ziziphus jujuba Mill

Hexokinase (HXK) plays a crucial role in plants, catalyzing the phosphorylation of hexose substances, which is one of the key steps in sugar metabolism and energy production. While HXK genes have been well-studied in model plants, the evolutionary and functional characteristics of HXK gene family in jujube is unknow. In this study, the HXK gene family members were identified by bioinformatics methods, the key members regulating glucose metabolism were identified by transcriptome data, and finally the function of the key genes was verified by instantaneous and stable genetic transformation. Our results showed that seven HXK genes were identified in the jujube genome, all of which were predict located in the chloroplast and contain Hexokinase-1 (PF00349) and Hexokinase-2 (PF03727) conserved domains. Most of HXK proteins were transmembrane protein with stable, lipid-soluble, hydrophilic. The secondary structure of ZjHXK proteins main α-helix, and contains two distinct tertiary structure. All ZjHXK genes contain nine exons and eight introns. Predictions of cis-regulatory elements indicate that the promoter region of ZjHXK contains a large number of MeJA responsive elements. Finally, combined with the analysis of the relationship between the expression and glucose metabolism, found that ZjHXK5 and ZjHXK6 may the key genes regulating sugar metabolism. Transient overexpression of ZjHXK5 and ZjHXK6 on jujube, or allogeneic overexpression of ZjHXK5 and ZjHXK6 on tomato would significantly reduce the content of total sugar and various sugar components. Transient silencing of ZjHXK5 and ZjHXK6 genes results in a significant increase in sucrose and total sugar content. Interestingly, the expression of ZjHXK5 and ZjHXK6 were also affected by methyl jasmonate.

Genomic investigation of duplication, functional conservation, and divergence in the LRR-RLK Family of Saccharum

BACKGROUND

Sugarcane (Saccharum spp.) holds exceptional global significance as a vital crop, serving as a primary source of sucrose, bioenergy, and various by-products. The optimization of sugarcane breeding by fine-tuning essential traits has become crucial for enhancing crop productivity and stress resilience. Leucine-rich repeat receptor-like kinases (LRR-RLK) genes present promising targets for this purpose, as they are involved in various aspects of plant development and defense processes.

RESULTS

Here, we present a detailed overview of phylogeny and expression of 288 (495 alleles) and 312 (1365 alleles) LRR-RLK genes from two founding Saccharum species, respectively. Phylogenetic analysis categorized these genes into 15 subfamilies, revealing considerable expansion or reduction in certain LRR-type subfamilies. Compared to other plant species, both Saccharum species had more significant LRR-RLK genes. Examination of cis-acting elements demonstrated that SsLRR-RLK and SoLRR-RLK genes exhibited no significant difference in the types of elements included, primarily involved in four physiological processes. This suggests a broad conservation of LRR-RLK gene function during Saccharum evolution. Synteny analysis indicated that all LRR-RLK genes in both Saccharum species underwent gene duplication, primarily through whole-genome duplication (WGD) or segmental duplication. We identified 28 LRR-RLK genes exhibiting novel expression patterns in response to different tissues, gradient development leaves, and circadian rhythm in the two Saccharum species. Additionally, SoLRR-RLK104, SoLRR-RLK7, SoLRR-RLK113, and SsLRR-RLK134 were identified as candidate genes for sugarcane disease defense response regulators through transcriptome data analysis of two disease stresses. This suggests LRR-RLK genes of sugarcane involvement in regulating various biological processes, including leaf development, plant morphology, photosynthesis, maintenance of circadian rhythm stability, and defense against sugarcane diseases.

CONCLUSIONS

This investigation into gene duplication, functional conservation, and divergence of LRR-RLK genes in two founding Saccharum species lays the groundwork for a comprehensive genomic analysis of the entire LRR-RLK gene family in Saccharum. The results reveal LRR-RLK gene played a critical role in Saccharum adaptation to diverse conditions, offering valuable insights for targeted breeding and precise phenotypic adjustments.

Genome-wide identification of GATA transcription factor family and the effect of different light quality on the accumulation of terpenoid indole alkaloids in Uncaria rhynchophylla

Uncaria rhynchophylla is an evergreen vine plant, belonging to the Rubiaceae family, that is rich in terpenoid indole alkaloids (TIAs) that have therapeutic effects on hypertension and Alzheimer's disease. GATA transcription factors (TF) are a class of transcription regulators that participate in the light response regulation, chlorophyll synthesis, and metabolism, with the capability to bind to GATA cis-acting elements in the promoter region of target genes. Currently the charactertics of GATA TFs in U. rhynchophylla and how different light qualities affect the expression of GATA and key enzyme genes, thereby affecting the changes in U. rhynchophylla alkaloids have not been investigated. In this study, 25 UrGATA genes belonging to four subgroups were identified based on genome-wide analysis. Intraspecific collinearity analysis revealed that only segmental duplications were identified among the UrGATA gene family. Collinearity analysis of GATA genes between U. rhynchophylla and four representative plant species, Arabidopsis thaliana, Oryza sativa, Coffea Canephora, and Catharanthus roseus was also performed. U. rhynchophylla seedlings grown in either red lights or under reduced light intensity had altered TIAs content after 21 days. Gene expression analysis reveal a complex pattern of expression from the 25 UrGATA genes as well as a number of key TIA enzyme genes. UrGATA7 and UrGATA8 were found to have similar expression profiles to key enzyme TIA genes in response to altered light treatments, implying that they may be involved in the regulation TIA content. In this research, we comprehensively analyzed the UrGATA TFs, and offered insight into the involvement of UrGATA TFs from U. rhynchophylla in TIAs biosynthesis.

Genome-Wide Identification and Expression Profile Analysis of Sugars Will Eventually Be Exported Transporter (SWEET) Genes in Zantedeschia elliottiana and Their Responsiveness to Pectobacterium carotovora subspecies Carotovora (Pcc) Infection

SWEET, sugars will eventually be exported transporter, is a novel class of sugar transporter proteins that can transport sugars across membranes down a concentration gradient. It plays a key role in plant photosynthetic assimilates, phloem loading, nectar secretion from nectar glands, seed grouting, pollen development, pathogen interactions, and adversity regulation, and has received widespread attention in recent years. To date, systematic analysis of the SWEET family in Zantedeschia has not been documented, although the genome has been reported in Zantedeschia elliottiana. In this study, 19 ZeSWEET genes were genome-wide identified in Z. elliottiana, and unevenly located in 10 chromosomes. They were further clustered into four clades by a phylogenetic tree, and almost every clade has its own unique motifs. Synthetic analysis confirmed two pairs of segmental duplication events of ZeSWEET genes. Heatmaps of tissue-specific and Pectobacterium carotovora subsp. Carotovora (Pcc) infection showed that ZeSWEET genes had different expression patterns, so SWEETs may play widely varying roles in development and stress tolerance in Zantedeschia. Moreover, quantitative reverse transcription-PCR (qRT-PCR) analysis revealed that some of the ZeSWEETs responded to Pcc infection, among which eight genes were significantly upregulated and six genes were significantly downregulated, revealing their potential functions in response to Pcc infection. The promoter sequences of ZeSWEETs contained 51 different types of the 1380 cis-regulatory elements, and each ZeSWEET gene contained at least two phytohormone responsive elements and one stress response element. In addition, a subcellular localization study indicated that ZeSWEET07 and ZeSWEET18 were found to be localized to the plasma membrane. These findings provide insights into the characteristics of SWEET genes and contribute to future studies on the functional characteristics of ZeSWEET genes, and then improve Pcc infection tolerance in Zantedeschia through molecular breeding.

DOF gene family expansion and diversification

DOF (DNA binding with one finger) proteins are part of a plant-specific transcription factor (TF) gene family widely involved in plant development and stress responses. Many studies have uncovered their structural and functional characteristics in recent years, leading to a rising number of genome-wide identification study approaches, unveiling the DOF family expansion in angiosperm species. Nonetheless, these studies primarily concentrate on particular taxonomic groups. Identifying DOF TFs within less-represented groups is equally crucial, as it enhances our comprehension of their evolutionary history, contributions to plant phenotypic diversity, and role in adaptation. This review summarizes the main findings and progress of genome-wide identification and characterization studies of DOF TFs in Viridiplantae, exposing their roles as players in plant adaptation and a glimpse of their evolutionary history. We also present updated data on the identification and number of DOF genes in native and wild species. Altogether, these data, comprising a phylogenetic analysis of 2124 DOF homologs spanning 83 different species, will contribute to identifying new functional DOF groups, adding to our understanding of the mechanisms driving plant evolution and offering valuable insights into their potential applications.

The evolution and functional divergence of 10 Apolipoprotein D-like genes in Nilaparvata lugens

Apolipoprotein D (ApoD), a member of the lipocalin superfamily of proteins, is involved in lipid transport and stress resistance. Whereas only a single copy of the ApoD gene is found in humans and some other vertebrates, there are typically several ApoD-like genes in insects. To date, there have been relatively few studies that have examined the evolution and functional differentiation of ApoD-like genes in insects, particularly hemi-metabolous insects. In this study, we identified 10 ApoD-like genes (NlApoD1-10) with distinct spatiotemporal expression patterns in Nilaparvata lugens (BPH), which is an important pest of rice. NlApoD1-10 were found to be distributed on 3 chromosomes in a tandem array of NlApoD1/2, NlApoD3-5, and NlApoD7/8, and show sequence and gene structural divergence in the coding regions, indicating that multiple gene duplication events occurred during evolution. Phylogenetic analysis revealed that NlApoD1-10 can be clustered into 5 clades, with NlApoD3-5 and NlApoD7/8 potentially evolving exclusively in the Delphacidae family. Functional screening using an RNA interference approach revealed that only NlApoD2 was essential for BPH development and survival, whereas NlApoD4/5 are highly expressed in testes, and might play roles in reproduction. Moreover, stress response analysis revealed that NlApoD3-5/9, NlApoD3-5, and NlApoD9 were up-regulated after treatment with lipopolysaccharide, H2 O2 , and ultraviolet-C, respectively, indicating their potential roles in stress resistance.

Comparative genomics reveals the dynamic evolutionary history of cement protein genes of barnacles from intertidal to deep-sea hydrothermal vents

Thoracican barnacles are a diverse group of marine organisms for which the availability of genome assemblies is currently limited. In this study, we sequenced the genomes of two neolepadoid species (Ashinkailepas kermadecensis, Imbricaverruca yamaguchii) from hydrothermal vents, in addition to two intertidal species. Genome sizes ranged from 481 to 1054 Mb, with repetitive sequence contents of 21.2% to 50.7%. Concordance rates of orthologs and heterozygosity rates were between 82.4% and 91.7% and between 1.0% and 2.1%, respectively, indicating high genetic diversity and heterozygosity. Based on phylogenomic analyses, we revised the nomenclature of cement genes encoding cement proteins that are not homologous to any known proteins. The major cement gene, CP100A, was found in all thoracican species, including vent-associated neolepadoids, and was hypothesised to be essential for thoracican settlement. Duplicated genes, CP100B and CP100C, were found only in balanids, suggesting potential functional redundancy or acquisition of new functions associated with the calcareous base. An ancestor of CP52 genes was duplicated dynamically among lepadids, pollicipedids with multiple copies on a single scaffold, and balanids with multiple sequential repeats of the conserved regions, but no CP52 genes were found in neolepadoids, providing insights into cement gene evolution among thoracican lineages. This study enhances our understanding of the adhesion mechanisms of thoracicans in underwater environments. The newly sequenced genomes provide opportunities for studying their evolution and ecology, shedding light on their adaptation to diverse marine environments, and contributing to our knowledge of barnacle biology with valuable genomic resources for further studies in this field.

ZmNRT1.1B (ZmNPF6.6) determines nitrogen use efficiency via regulation of nitrate transport and signalling in maize

Nitrate (NO3 - ) is crucial for optimal plant growth and development and often limits crop productivity under low availability. In comparison with model plant Arabidopsis, the molecular mechanisms underlying NO3 - acquisition and utilization remain largely unclear in maize. In particular, only a few genes have been exploited to improve nitrogen use efficiency (NUE). Here, we demonstrated that NO3 - -inducible ZmNRT1.1B (ZmNPF6.6) positively regulated NO3 - -dependent growth and NUE in maize. We showed that the tandem duplicated proteoform ZmNRT1.1C is irrelevant to maize seedling growth under NO3 - supply; however, the loss of function of ZmNRT1.1B significantly weakened plant growth under adequate NO3 - supply under both hydroponic and field conditions. The 15 N-labelled NO3 - absorption assay indicated that ZmNRT1.1B mediated the high-affinity NO3 - -transport and root-to-shoot NO3 - translocation. Transcriptome analysis further showed, upon NO3 - supply, ZmNRT1.1B promotes cytoplasmic-to-nuclear shuttling of ZmNLP3.1 (ZmNLP8), which co-regulates the expression of genes involved in NO3 - response, cytokinin biosynthesis and carbon metabolism. Remarkably, overexpression of ZmNRT1.1B in modern maize hybrids improved grain yield under N-limiting fields. Taken together, our study revealed a crucial role of ZmNRT1.1B in high-affinity NO3 - transport and signalling and offers valuable genetic resource for breeding N use efficient high-yield cultivars.

Uncovering Evolutionary Adaptations in Common Warthogs through Genomic Analyses

In the Suidae family, warthogs show significant survival adaptability and trait specificity. This study offers a comparative genomic analysis between the warthog and other Suidae species, including the Luchuan pig, Duroc pig, and Red River hog. By integrating the four genomes with sequences from the other four species, we identified 8868 single-copy orthologous genes. Based on 8868 orthologous protein sequences, phylogenetic assessments highlighted divergence timelines and unique evolutionary branches within suid species. Warthogs exist on different evolutionary branches compared to DRCs and LCs, with a divergence time preceding that of DRC and LC. Contraction and expansion analyses of warthog gene families have been conducted to elucidate the mechanisms of their evolutionary adaptations. Using GO, KEGG, and MGI databases, warthogs showed a preference for expansion in sensory genes and contraction in metabolic genes, underscoring phenotypic diversity and adaptive evolution direction. Associating genes with the QTLdb-pigSS11 database revealed links between gene families and immunity traits. The overlap of olfactory genes in immune-related QTL regions highlighted their importance in evolutionary adaptations. This work highlights the unique evolutionary strategies and adaptive mechanisms of warthogs, guiding future research into the distinct adaptability and disease resistance in pigs, particularly focusing on traits such as resistance to African Swine Fever Virus.

Genome-Wide Comparative Analysis of SRCR Gene Superfamily in Invertebrates Reveals Massive and Independent Gene Expansions in the Sponge and Sea Urchin

Without general adaptative immunity, invertebrates evolved a vast number of heterogeneous non-self recognition strategies. One of those well-known adaptations is the expansion of the immune receptor gene superfamily coding for scavenger receptor cysteine-rich domain containing proteins (SRCR) in a few invertebrates. Here, we investigated the evolutionary history of the SRCR gene superfamily (SRCR-SF) across 29 metazoan species with an emphasis on invertebrates. We analyzed their domain architectures, genome locations and phylogenetic distribution. Our analysis shows extensive genome-wide duplications of the SRCR-SFs in Amphimedon queenslandica and Strongylocentrotus purpuratus. Further molecular evolution study reveals various patterns of conserved cysteines in the sponge and sea urchin SRCR-SFs, indicating independent and convergent evolution of SRCR-SF expansion during invertebrate evolution. In the case of the sponge SRCR-SFs, a novel motif with seven conserved cysteines was identified. Exon-intron structure analysis suggests the rapid evolution of SRCR-SFs during gene duplications in both the sponge and the sea urchin. Our findings across nine representative metazoans also underscore a heightened expression of SRCR-SFs in immune-related tissues, notably the digestive glands. This observation indicates the potential role of SRCR-SFs in reinforcing distinct immune functions in these invertebrates. Collectively, our results reveal that gene duplication, motif structure variation, and exon-intron divergence might lead to the convergent evolution of SRCR-SF expansions in the genomes of the sponge and sea urchin. Our study also suggests that the utilization of SRCR-SF receptor duplication may be a general and basal strategy to increase immune diversity and tissue specificity for the invertebrates.