Structurally divergent and recurrently mutated regions of primate genomes.

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

The Penium margaritaceum Genome: Hallmarks of the Origins of Land Plants.

The evolutionary features and molecular innovations that enabled plants to first colonize land are not well understood. Here, insights are provided through our report of the genome sequence of the unicellular alga Penium margaritaceum, a member of the Zygnematophyceae, the sister lineage to land plants. The genome has a high proportion of repeat sequences that are associated with massive segmental gene duplications, likely facilitating neofunctionalization. Compared with representatives of earlier diverging algal lineages, P. margaritaceum has expanded repertoires of gene families, signaling networks, and adaptive responses that highlight the evolutionary trajectory toward terrestrialization. These encompass a broad range of physiological processes and protective cellular features, such as flavonoid compounds and large families of modifying enzymes involved in cell wall biosynthesis, assembly, and remodeling. Transcriptome profiling further elucidated adaptations, responses, and selective pressures associated with the semi-terrestrial ecosystems of P. margaritaceum, where a simple body plan would be an advantage.

Transcriptional Architecture of Synaptic Communication Delineates GABAergic Neuron Identity.

Understanding the organizational logic of neural circuits requires deciphering the biological basis of neuronal diversity and identity, but there is no consensus on how neuron types should be defined. We analyzed single-cell transcriptomes of a set of anatomically and physiologically characterized cortical GABAergic neurons and conducted a computational genomic screen for transcriptional profiles that distinguish them from one another. We discovered that cardinal GABAergic neuron types are delineated by a transcriptional architecture that encodes their synaptic communication patterns. This architecture comprises 6 categories of ∼40 gene families, including cell-adhesion molecules, transmitter-modulator receptors, ion channels, signaling proteins, neuropeptides and vesicular release components, and transcription factors. Combinatorial expression of select members across families shapes a multi-layered molecular scaffold along the cell membrane that may customize synaptic connectivity patterns and input-output signaling properties. This molecular genetic framework of neuronal identity integrates cell phenotypes along multiple axes and provides a foundation for discovering and classifying neuron types.

The immunoglobulin J chain is an evolutionarily co-opted chemokine.

The joining (J) chain regulates polymerization of multimeric Immunoglobulin(Ig)M and IgA, forming a disulfide bond to the C termini of their Ig heavy chains, and it controls IgM/IgA transport across mucosal epithelia. Like Ig itself and human-like adaptive immunity, J chain emerged in jawed vertebrates (gnathostomes), but its origin has remained mysterious since its discovery over 50 y ago. Here, we show unexpectedly that J chain is a member of the CXCL chemokine family. The J chain gene (JCHAIN) is linked to clustered CXCL chemokine loci in all gnathostomes except actinopterygians that lost JCHAIN. JCHAIN and most CXCL genes have four exons with the same intron phases, including the same cleavage site for the signal peptide/mature protein. The second exon of both genes encodes a CXC motif at the same position, and the lengths of exons 1 to 3 are similar. No other gene in the human secretome shares all of these characteristics. In contrast, intrachain disulfide bonds of the two proteins are completely different, likely due to modifications in J chain to direct Ig polymerization and mucosal transport. Crystal structures of CXCL8 and J chain share a conserved beta-strand core but diverge otherwise due to different intrachain disulfide bonds and extension of the J chain C terminus. Identification of this ancestral affiliation between J chain and CXCL chemokines addresses an age-old problem in immunology.

Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae.

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has since followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appears to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, Sex Peptide Receptor, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Insight into the mechanism of H+-coupled nucleobase transport.

Members of the nucleobase/ascorbic acid transporter (NAT) gene family are found in all kingdoms of life. In mammals, the concentrative uptake of ascorbic acid (vitamin C) by members of the NAT family is driven by the Na+ gradient, while the uptake of nucleobases in bacteria is powered by the H+ gradient. Here, we report the structure and function of PurTCp, a NAT family member from Colwellia psychrerythraea. The structure of PurTCp was determined to 2.80 Å resolution by X-ray crystallography. PurTCp forms a homodimer, and each protomer has 14 transmembrane segments folded into a transport domain (core domain) and a scaffold domain (gate domain). A purine base is present in the structure and defines the location of the substrate binding site. Functional studies reveal that PurTCp transports purines but not pyrimidines and that purine binding and transport is dependent on the pH. Mutation of a conserved aspartate residue close to the substrate binding site reveals the critical role of this residue in H+-dependent transport of purines. Comparison of the PurTCp structure with transporters of the same structural fold suggests that rigid-body motions of the substrate-binding domain are central for substrate translocation across the membrane.

AT-hook motif nuclear localized transcription factors function redundantly in promoting root growth through modulation of redox homeostasis.

Maintaining an optimal redox status is essential for plant growth and development, particularly when the plants are under stress. AT-hook motif nuclear localized (AHL) proteins are evolutionarily conserved transcription factors in plants. Much of our understanding about this gene family has been derived from studies on clade A members. To elucidate the functions of clade B genes, we first analyzed their spatial expression patterns using transgenic plants expressing a nuclear localized GFP under the control of their promoter sequences. AHL1, 2, 6, 7, and 10 were further functionally characterized owing to their high expression in the root apical meristem. Through mutant analyses and transgenic studies, we showed that these genes have the ability to promote root growth. Using yeast one-hybrid and dual luciferase assays, we demonstrated that AHL1, 2, 6, 7, and 10 are transcription regulators and this activity is required for their roles in root growth. Although mutants for these genes did not showed obvious defects in root growth, transgenic plants expressing their fusion proteins with the SRDX repressor motif exhibited a short-root phenotype. Through transcriptome analysis, histochemical staining and molecular genetics experiments, we found that AHL10 maintains redox homeostasis via direct regulation of glutathione transferase (GST) genes. When the transcript level of GSTF2, a top-ranked target of AHL10, was reduced by RNAi, the short-root phenotype in the AHL10-SRDX expressing plant was largely rescued. These results together suggest that AHL genes function redundantly in promoting root growth through direct regulation of redox homeostasis.

The Evolution of Transglutaminases Underlies the Origin and Loss of Cornified Skin Appendages in Vertebrates.

Transglutaminases (TGMs) cross-link proteins by introducing covalent bonds between glutamine and lysine residues. These cross-links are essential for epithelial cornification which enables tetrapods to live on land. Here, we investigated which evolutionary adaptations of vertebrates were associated with specific changes in the family of TGM genes. We determined the catalog of TGMs in the main clades of vertebrates, performed a comprehensive phylogenetic analysis of TGMs, and localized the distribution of selected TGMs in tissues. Our data suggest that TGM1 is the phylogenetically oldest epithelial TGM, with orthologs being expressed in the cornified teeth of the lamprey, a basal vertebrate. Gene duplications led to the origin of TGM10 in stem vertebrates, the origin of TGM2 in jawed vertebrates, and an increasing number of epithelium-associated TGM genes in the lineage leading to terrestrial vertebrates. TGM9 is expressed in the epithelial egg tooth, and its evolutionary origin in stem amniotes coincided with the evolution of embryonic development in eggs that are surrounded by a protective shell. Conversely, viviparous mammals have lost both the epithelial egg tooth and TGM9. TGM3 and TGM6 evolved as regulators of cornification in hair follicles and underwent pseudogenization upon the evolutionary loss of hair in cetaceans. Taken together, this study reveals the gain and loss of vertebrate TGM genes in association with the evolution of cornified skin appendages and suggests an important role of TGM9 in the evolution of amniotes.

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles.

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

A novel broad spectrum venom metalloproteinase autoinhibitor in the rattlesnake Crotalus atrox evolved via a shift in paralog function.

The complexity of snake venom composition reflects adaptation to the diversity of prey and may be driven at times by a coevolutionary arms race between snakes and venom-resistant prey. However, many snakes are also resistant to their own venom due to serum-borne inhibitors of venom toxins, which raises the question of how snake autoinhibitors maintain their efficacy as venom proteins evolve. To investigate this potential three-way arms race among venom, prey, and autoinhibitors, we have identified and traced the evolutionary origin of serum inhibitors of snake venom metalloproteinases (SVMPs) in the Western Diamondback rattlesnake Crotalus atrox which possesses the largest known battery of SVMP genes among crotalids examined. We found that C. atrox expresses five members of a Fetuin A-related metalloproteinase inhibitor family but that one family member, FETUA-3, is the major SVMP inhibitor that binds to approximately 20 different C. atrox SVMPs and inhibits activities of all three SVMP classes. We show that the fetua-3 gene arose deep within crotalid evolution before the origin of New World species but, surprisingly, fetua-3 belongs to a different paralog group than previously identified SVMP inhibitors in Asian and South American crotalids. Conversely, the C. atrox FETUA-2 ortholog of previously characterized crotalid SVMP inhibitors shows limited activity against C. atrox SVMPs. These results reveal that there has been a functional evolutionary shift in the major SVMP inhibitor in the C. atrox lineage as the SVMP family expanded and diversified in the Crotalus lineage. This broad-spectrum inhibitor may be of potential therapeutic interest.

Complex evolutionary patterns within the tubulin gene family of ciliates, unicellular eukaryotes with diverse microtubular structures.

BACKGROUND: Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear. RESULTS: We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, ß_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, ß_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and ß_Tub varied among different classes. (3) Ciliate α- and ß-tubulin isoforms could be classified into an "ancestral group" present in LECA and a "divergent group" containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the "ancestral group" of α_Tub and ß_Tub. The "divergent group" might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events. CONCLUSIONS: Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution.

Whole genome sequencing of Crassostrea ariakensis (Mollusca: Ostreidae) and C. hongkongensis expands understandings of stress resistance in sessile oysters.

To understand the environmental adaptations among sessile bivalves lacking adaptive immunity, a series of analyses were conducted, with special emphasis on the widely distributed C. ariakensis. Employing Pacbio sequencing and Hi-C technologies, whole genome for each of a C. ariakensis (southern China) and C. hongkongensis individual was generated, with the contig N50 reaching 6.2 and 13.0 Mb, respectively. Each genome harbored over 30,000 protein-coding genes, with approximately half of each genome consisting of repeats. Genome alignment suggested possible introgression between C. gigas and C. ariakensis (northern China), and re-sequencing data corroborated this result and indicated significant gene flow between C. gigas and C. ariakensis. These introgressed candidates, well-represented by genes related to immunity and osmotic pressure, may be associated with environmental stresses. Gene family dynamics modeling suggested immune-related genes were well represented among the expanded genes in C. ariakensis. These outcomes could be attributed to the spread of C. ariakensis.

Genome-wide identification of the TIFY gene family in tobacco and expression analysis in response to Ralstonia solanacearum infection.

The TIFY gene family plays an essential role in plant development and abiotic and biotic stress responses. In this study, genome-wide identification of TIFY members in tobacco and their expression pattern analysis in response to Ralstonia solanacearum infection were performed. A total of 33 TIFY genes were identified, including the TIFY, PPD, ZIM&ZML and JAZ subfamilies. Promoter analysis results indicated that a quantity of light-response, drought-response, SA-response and JA-response cis-elements exist in promoter regions. The TIFY gene family exhibited expansion and possessed gene redundancy resulting from tobacco ploidy change. In addition, most NtTIFYs equivalently expressed in roots, stems and leaves, while NtTIFY1, NtTIFY4, NtTIFY18 and NtTIFY30 preferentially expressed in roots. The JAZ III clade showed significant expression changes after inoculation with R. solanacearum, and the expression of NtTIFY7 in resistant varieties, compared with susceptible varieties, was more stably induced. Furthermore, NtTIFY7-silenced plants, compared with the control plants, were more susceptible to bacterial wilt. These results lay a foundation for exploring the evolutionary history of TIFY gene family and revealing gene function of NtTIFYs in tobacco bacterial wilt resistance.

RhoA rescues cardiac senescence by regulating Parkin-mediated mitophagy.

Heart failure is one of the leading causes of death worldwide. RhoA, a small GTPase, governs actin dynamics in various tissue and cell types, including cardiomyocytes; however, its involvement in cardiac function has not been fully elucidated. Here, we generated cardiomyocyte-specific RhoA conditional knockout (cKO) mice, which demonstrated a significantly shorter lifespan with left ventricular dilation and severely impaired ejection fraction. We found that the cardiac tissues of the cKO mice exhibited structural disorganization with fibrosis and also exhibited enhanced senescence compared with control mice. In addition, we show that cardiomyocyte mitochondria were structurally abnormal in the aged cKO hearts. Clearance of damaged mitochondria by mitophagy was remarkably inhibited in both cKO cardiomyocytes and RhoA-knockdown HL-1 cultured cardiomyocytes. In RhoA-depleted cardiomyocytes, we reveal that the expression of Parkin, an E3 ubiquitin ligase that plays a crucial role in mitophagy, was reduced, and expression of N-Myc, a negative regulator of Parkin, was increased. We further reveal that the RhoA-Rho kinase axis induced N-Myc phosphorylation, which led to N-Myc degradation and Parkin upregulation. Re-expression of Parkin in RhoA-depleted cardiomyocytes restored mitophagy, reduced mitochondrial damage, attenuated cardiomyocyte senescence, and rescued cardiac function both in vitro and in vivo. Finally, we found that patients with idiopathic dilated cardiomyopathy without causal mutations for dilated cardiomyopathy showed reduced cardiac expression of RhoA and Parkin. These results suggest that RhoA promotes Parkin-mediated mitophagy as an indispensable mechanism contributing to cardioprotection in the aging heart.

Palmitoylation-dependent regulation of cardiomyocyte Rac1 signaling activity and minor effects on cardiac hypertrophy.

S-palmitoylation is a reversible lipid modification catalyzed by 23 S-acyltransferases with a conserved zinc finger aspartate-histidine-histidine-cysteine (zDHHC) domain that facilitates targeting of proteins to specific intracellular membranes. Here we performed a gain-of-function screen in the mouse and identified the Golgi-localized enzymes zDHHC3 and zDHHC7 as regulators of cardiac hypertrophy. Cardiomyocyte-specific transgenic mice overexpressing zDHHC3 show cardiac disease, and S-acyl proteomics identified the small GTPase Rac1 as a novel substrate of zDHHC3. Notably, cardiomyopathy and congestive heart failure in zDHHC3 transgenic mice is preceded by enhanced Rac1 S-palmitoylation, membrane localization, activity, downstream hypertrophic signaling, and concomitant induction of all Rho family small GTPases whereas mice overexpressing an enzymatically dead zDHHC3 mutant show no discernible effect. However, loss of Rac1 or other identified zDHHC3 targets Gαq/11 or galectin-1 does not diminish zDHHC3-induced cardiomyopathy, suggesting multiple effectors and pathways promoting decompensation with sustained zDHHC3 activity. Genetic deletion of Zdhhc3 in combination with Zdhhc7 reduces cardiac hypertrophy during the early response to pressure overload stimulation but not over longer time periods. Indeed, cardiac hypertrophy in response to 2 weeks of angiotensin-II infusion is not diminished by Zdhhc3/7 deletion, again suggesting other S-acyltransferases or signaling mechanisms compensate to promote hypertrophic signaling. Taken together, these data indicate that the activity of zDHHC3 and zDHHC7 at the cardiomyocyte Golgi promote Rac1 signaling and maladaptive cardiac remodeling, but redundant signaling effectors compensate to maintain cardiac hypertrophy with sustained pathological stimulation in the absence of zDHHC3/7.

Comprehensive genome-wide analysis of the DREB gene family in Moso bamboo (Phyllostachys edulis): evidence for the role of PeDREB28 in plant abiotic stress response.

Dehydration response element binding (DREB) proteins are vital for plant abiotic stress responses, but the understanding of DREBs in bamboo, an important sustainable non-timber forest product, is limited. Here we conducted a comprehensive genome-wide analysis of the DREB gene family in Moso bamboo, representing the most important running bamboo species in Asia. In total, 44 PeDREBs were identified, and information on their gene structures, protein motifs, phylogenetic relationships, and stress-related cis-regulatory elements (CREs) was provided. Based on the bioinformatical analysis, we further analyzed PeDREBs from the A5 group and found that four of five PeDREB transcripts were induced by salt, drought, and cold stresses, and their proteins could bind to stress-related CREs. Among these, PeDREB28 was selected as a promising candidate for further functional characterization. PeDREB28 is localized in nucleus, has transcriptional activation activity, and could bind to the DRE- and coupling element 1- (CE1) CREs. Overexpression of PeDREB28 in Arabidopsis and bamboo improved plant abiotic stress tolerance. Transcriptomic analysis showed that broad changes due to the overexpression of PeDREB28. Furthermore, 628 genes that may act as the direct PeDREB28 downstream genes were identified by combining DAP-seq and RNA-seq analysis. Moreover, we confirmed that PeDREB28 could bind to the promoter of pyrabactin-resistance-like gene (DlaPYL3), which is a homolog of abscisic acid receptor in Arabidopsis, and activates its expression. In summary, our study provides important insights into the DREB gene family in Moso bamboo, and contributes to their functional verification and genetic engineering applications in the future.

Genome-wide molecular evolution analysis of the GRF and GIF gene families in Plantae (Archaeplastida).

BACKGROUND: Plant growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) interact with each other and collectively have important regulatory roles in plant growth, development, and stress responses. Therefore, it is of great significance to explore the systematic evolution of GRF and GIF gene families. However, our knowledge and understanding of the role of GRF and GIF genes during plant evolution has been fragmentary. RESULTS: In this study, a large number of genomic and transcriptomic datasets of algae, mosses, ferns, gymnosperms and angiosperms were used to systematically analyze the evolution of GRF and GIF genes during the evolution of plants. The results showed that GRF gene first appeared in the charophyte Klebsormidium nitens, whereas the GIF genes originated relatively early, and these two gene families were mainly expanded by segmental duplication events after plant terrestrialization. During the process of evolution, the protein sequences and functions of GRF and GIF family genes are relatively conservative. As cooperative partner, GRF and GIF genes contain the similar types of cis-acting elements in their promoter regions, which enables them to have similar transcriptional response patterns, and both show higher levels of expression in reproductive organs and tissues and organs with strong capacity for cell division. Based on protein-protein interaction analysis and verification, we found that the GRF-GIF protein partnership began to be established in pteridophytes and is highly conserved across different terrestrial plants. CONCLUSIONS: These results provide a foundation for further exploration of the molecular evolution and biological functions of GRF and GIF genes.

Genome-wide identification, phylogeny and expression analysis of Hsf gene family in Verbena bonariensis under low-temperature stress.

BACKGROUND: The heat shock transcription factor (Hsf) is a crucial regulator of plant stress resistance, playing a key role in plant stress response, growth, and development regulation. RESULTS: In this study, we utilized bioinformatics tools to screen 25 VbHsf members, which were named VbHsf1-VbHsf25. We used bioinformatics methods to analyze the sequence structure, physicochemical properties, conserved motifs, phylogenetic evolution, chromosome localization, promoter cis-acting elements, collinearity, and gene expression of Hsf heat shock transcription factor family members under low-temperature stress. The results revealed that the majority of the Hsf genes contained motif1, motif2, and motif3, signifying that these three motifs were highly conserved in the Hsf protein sequence of Verbena bonariensis. Although there were some variations in motif deletion among the members, the domain remained highly conserved. The theoretical isoelectric point ranged from 4.17 to 9.71, with 21 members being unstable proteins and the remainder being stable proteins. Subcellular localization predictions indicated that all members were located in the nucleus. Phylogenetic analysis of the Hsf gene family in V. bonariensis and Arabidopsis thaliana revealed that the Hsf gene family of V. bonariensis could be categorized into three groups, with group A comprising 17 members and group C having at least two members. Among the 25 Hsf members, there were 1-3 exons located on seven chromosome fragments, which were unevenly distributed. Collinearity analysis demonstrated the presence of seven pairs of homologous genes in the VbHsf gene family. The Ka/Ks ratios were less than one, indicating that the VbHsf gene underwent purification selection pressure. Additionally, nine genes in V. bonariensis were found to have collinearity with A. thaliana. Promoter analysis revealed that the promoters of all VbHsf genes contained various types of cis-acting elements related to hormones and stress. Based on RNA-seq data, qRT-PCR analysis of six highly expressed genes was performed, and it was found that VbHsf5, VbHsf14, VbHsf17, VbHsf18, VbHsf20 and VbHsf21 genes were highly expressed at 12 h of low-temperature treatment, and the expression decreased after 24 h, among which VbHsf14 was up-regulated at 12 h of low-temperature by 70-fold. CONCLUSIONS: Our study may help reveal the important roles of Hsf in plant development and show insight for the further molecular breeding of V. bonariensis.

Comprehensive analysis of pepper (Capsicum annuum) RAV genes family and functional identification of CaRAV1 under chilling stress.

BACKGROUND: Despite its known significance in plant abiotic stress responses, the role of the RAV gene family in the response of Capsicum annuum to chilling stress remains largely unexplored. RESULTS: In this study, we identified and characterized six members of the CaRAV gene subfamily in pepper plants through genome-wide analysis. Subsequently, the CaRAV subfamily was classified into four branches based on homology with Arabidopsis thaliana, each exhibiting relatively conserved domains within the branch. We discovered that light response elements accounted for the majority of CaRAVs, whereas low-temperature response elements were specific to the NGA gene subfamily. After pepper plants were subjected to chilling stress, qRT‒PCR analysis revealed that CaRAV1, CaRAV2 and CaNGA1 were significantly induced in response to chilling stress, indicating that CaRAVs play a role in the response to chilling stress. Using virus-induced gene silencing (VIGS) vectors, we targeted key members of the CaRAV gene family. Under normal growth conditions, the MDA content and SOD enzyme activity of the silenced plants were slightly greater than those of the control plants, and the REC activity was significantly greater than that of the control plants. The levels of MDA and electrolyte leakage were greater in the silenced plants after they were exposed to chilling stress, and the POD and CAT enzyme activities were significantly lower than those in the control, which was particularly evident under repeated chilling stress. In addition, the relative expression of CaPOD and CaCAT was greater in V2 plants upon repeated chilling stress, especially CaCAT was significantly greater in V2 plants than in the other two silenced plants, with 3.29 and 1.10 increases within 12 and 24 h. These findings suggest that CaRAV1 and CaNGA1 positively regulate the response to chilling stress. CONCLUSIONS: Silencing of key members of the CaRAV gene family results in increased susceptibility to chilling damage and reduced antioxidant enzyme activity in plants, particularly under repeated chilling stress. This study provides valuable information for understanding the classification and putative functions of RAV transcription factors in pepper plants.

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

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