Unveiling kiwifruit TCP genes: evolution, functions, and expression insights

The TEOSINTE-BRANCHED1/CYCLOIDEA/PROLEFERATING-CELL-FACTORS (TCP) gene family is a plant-specific transcriptional factor family involved in leaf morphogenesis and senescence, lateral branching, hormone crosstalk, and stress responses. To date, a systematic study on the identification and characterization of the TCP gene family in kiwifruit has not been reported. Additionally, the function of kiwifruit TCPs in regulating kiwifruit responses to the ethylene treatment and bacterial canker disease pathogen (Pseudomonas syringae pv. actinidiae, Psa) has not been investigated. Here, we identified 40 and 26 TCP genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes, respectively. The synteny analysis of AcTCPs illustrated that whole-genome duplication accounted for the expansion of the TCP family in Ac. Phylogenetic, conserved domain, and selection pressure analysis indicated that TCP family genes in Ac and Ae had undergone different evolutionary patterns after whole-genome duplication (WGD) events, causing differences in TCP gene number and distribution. Our results also suggested that protein structure and cis-element architecture in promoter regions of TCP genes have driven the function divergence of duplicated gene pairs. Three and four AcTCP genes significantly affected kiwifruit responses to the ethylene treatment and Psa invasion, respectively. Our results provided insight into general characters, evolutionary patterns, and functional diversity of kiwifruit TCPs.

Comprehensive analysis of circRNA expression profiles in postmenopausal women differing in bone mineral density

Postmenopausal osteoporosis (PMOP) seriously endangers the bone health of older women. Although there are currently indicators to diagnose PMOP, early diagnostic biomarkers are lacking. Circular ribonucleic acid (circRNA) has a stable structure, regulates gene expression, participates in the pathological process of disease, and has the potential to become a biomarker. The purpose of this study was to investigate circRNAs that could be used to predict patients with early PMOP. Ribonucleic acid (RNA) sequencing was performed on peripheral blood leukocytes from 15 female patients to identify differential circRNAs between different groups. Using bioinformatics analysis, enrichment analysis was performed to discover relevant functions and pathways. CircRNA-micro ribonucleic acid (miRNA) interaction analysis and messenger ribonucleic acid (mRNA) prediction and network construction help us to understand the relationship between circRNA, miRNA, and mRNA. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the gene expression of candidate circRNAs. We screened out 2 co-expressed differential circRNAs, namely hsa_circ_0060849 and hsa_circ_0001394. By analyzing the regulatory network, a total of 54 miRNAs and 57 osteoporosis-related mRNAs were identified, which, as potential downstream target genes of hsa_circ_0060849 and hsa_circ_0001394, may play a key role in the occurrence and development of PMOP. The occurrence and development of PMOP is regulated by circRNAs, and hsa_circ_0060849 and hsa_circ_0001394 can be used as new diagnostic markers and therapeutic targets for early PMOP.

Bioinformatics study and experimental evaluation of miR-182, and miR-34 expression profiles in Tuberculosis and lung cancer

Background

Lung cancer is one of the most dangerous diseases among cancers and tuberculosis is one of the deadliest infectious diseases in the world. Many studies have mentioned the connection between lung cancer and tuberculosis, and also the microRNAs that play a significant role in the development of these two diseases. This study aims to use different databases to find effective miRNAs and their role on different genes on lung and tuberculosis diseases. Also determining the role of miR-34a and miR-182 in lung cancer and tuberculosis.

Methods

Using the GEO database, the influential microRNA databases were studied in two diseases. Finally, regarding bioinformatics results and literature studies, two miR-34a and miR-182 were selected. The role of these microRNAs and their target genes was carefully evaluated using bioinformatics. The expression of microRNAs in the blood plasma of patients with lung cancer and tuberculosis and healthy people were investigated.

Results

According to the GEO database, miR-34a and miR-182 are microRNAs that affect tuberculosis and lung cancer. By checking the miRBase, miRcode, Diana, miRDB, galaxy, KEGG databases, the role of these microRNAs on genes and different molecular pathways and their effect on these microRNAs were mentioned. The results of the present study showed that the expression of miR-34a and miR-182 was lower than that of healthy people. The P value amount for miR-182 was <0.0001 and for miR-34a was 0.3380.

Conclusion

Reducing the expression pattern of these microRNAs indicates their role in lung cancer and tuberculosis occurrence. Therefore, these microRNAs can be used as a biomarker for prognosis, diagnosis, and treatment methods.

The Evolution, Expression Patterns, and Domestication Selection Analysis of the Annexin Gene Family in the Barley Pan-Genome

Plant annexins constitute a conserved protein family that plays crucial roles in regulating plant growth and development, as well as in responses to both biotic and abiotic stresses. In this study, a total of 144 annexin genes were identified in the barley pan-genome, comprising 12 reference genomes, including cultivated barley, landraces, and wild barley. Their chromosomal locations, physical-chemical characteristics, gene structures, conserved domains, and subcellular localizations were systematically analyzed to reveal the certain differences between wild and cultivated populations. Through a cis-acting element analysis, co-expression network, and large-scale transcriptome analysis, their involvement in growth, development, and responses to various stressors was highlighted. It is worth noting that HvMOREXann5 is only expressed in pistils and anthers, indicating its crucial role in reproductive development. Based on the resequencing data from 282 barley accessions worldwide, genetic variations in thefamily were investigated, and the results showed that 5 out of the 12 identified HvMOREXanns were affected by selection pressure. Genetic diversity and haplotype frequency showed notable reductions between wild and domesticated barley, suggesting that a genetic bottleneck occurred on the annexin family during the barley domestication process. Finally, qRT-PCR analysis confirmed the up-regulation of HvMOREXann7 under drought stress, along with significant differences between wild accessions and varieties. This study provides some insights into the genome organization and genetic characteristics of the annexin gene family in barley at the pan-genome level, which will contribute to better understanding its evolution and function in barley and other crops.

Genome-Wide Identification, Expression Analysis under Abiotic Stress and Co-Expression Analysis of MATE Gene Family in Torreya grandis

The multidrug and toxin efflux (MATE) family participates in numerous biological processes and plays important roles in abiotic stress responses. However, information about the MATE family genes in Torreya grandis remains unclear. In this study, our genome-wide investigation identified ninety MATE genes in Torreya grandis, which were divided into five evolutionary clades. TgMATE family members are located on eleven chromosomes, and a total of thirty TgMATEs exist in tandem duplication. The promoter analysis showed that most TgMATEs contain the cis-regulatory elements associated with stress and hormonal responses. In addition, we discovered that most TgMATE genes responded to abiotic stresses (aluminum, drought, high temperatures, and low temperatures). Weighted correlation network analysis showed that 147 candidate transcription factor genes regulated the expression of 14 TgMATE genes, and it was verified through a double-luciferase assay. Overall, our findings offer valuable information for the characterization of the TgMATE gene mechanism in responding to abiotic stress and exhibit promising prospects for the stress tolerance breeding of Torreya grandis.

Genome-Wide Characterization and Expression Analyses of Major Latex Protein Gene Family in Populus simonii × P. nigra

Major latex proteins, or MLPs, are crucial to plants' capacity to grow, develop, and endure biotic and abiotic stresses. The MLP gene family has been found in numerous plants, but little is known about its role in Populus simonii × P. nigra. This study discovered and assessed 43 PtMLP genes that were unevenly dispersed throughout 12 chromosomes in terms of their physicochemical characteristics, gene structure, conserved motifs, and protein localization. Based on their phylogeny and protein structural characteristics, three separate subclasses of PtMLP family were identified. Segmental and tandem duplication were found to be essential variables in the expansion of the PtMLP genes. The involvement of the PtMLP genes in growth and development, as well as in the responses to different hormones and stresses, was demonstrated by cis-regulatory element prediction. The PtMLP genes showed varying expression patterns in various tissues and under different conditions (cold, salt, and drought stress), as demonstrated in RNA-Seq databases, suggesting that PsnMLP may have different functions. Following the further investigation of the genes demonstrating notable variations in expression before and after the application of three stresses, PsnMLP5 was identified as a candidate gene. Subsequent studies revealed that PsnMLP5 could be induced by ABA treatment. This study paves the way for further investigations into the MLP genes' functional mechanisms in response to abiotic stressors, as well as the ways in which they can be utilized in poplar breeding for improved stress tolerance.

Genome-Wide Identification and Characterization of the Salvia miltiorrhiza Histone Deacetylase (HDAC) Family in Response to Multiple Abiotic Stresses

Salvia miltiorrhiza is a plant commonly used in traditional Chinese medicine. Its material bases for treating diseases are tanshinones and phenolic acids, including salvianolic acids. Histone deacetylase proteins (HDACs) are a class of specific functional enzymes that interact with acetylation groups on the N-terminal lysine of histone proteins further regulate gene transcription through structural changes at the chromatin level. HDACs involved in the growth and development of various plants, and induced by plant hormones to regulate the internal environment of plants to resist stress, at the same time affect the accumulation of some secondary metabolites. However, the role of SmHDACs on the accumulation of salvianolic acid in S. miltiorrhiza remains unclear. In this study, 16 SmHDACs genes were identified from the high-quality S. miltiorrhiza genome, their physicochemical properties were predicted. In phylogenetic trees co-constructed with HDACs proteins from other plants, SmHDACs was divided into three subfamilies, each with similar motif and conserved domain composition. The distribution of the three subfamilies is similar to that of dicotyledonous plants. Chromosome localization analysis showed that SmHDACs genes were randomly located. Cis-acting element analysis predicted that SmHDACs gene expression may be related to and induced by various phytohormones, such as MeJA and ABA. By combining the expression pattern and co-expression network induced by phytohormones, we speculate that SmHDACs may further influence the synthesis of salvianolic acid, and identified SmHDA5, a potential functional gene, then speculate its downstream target based on the co-expression network. In summary, we analyzed the SmHDACs gene family of S. miltiorrhiza and screened out the potential functional gene SmHDA5. From the perspective of epigenetics, we proposed the molecular mechanism of plant hormone promoting salvianolic acid synthesis, which filled the gap in the subdivision of histone deacetylase in S. miltiorrhiza research, provided a theoretical basis for the culture and transformation of S. miltiorrhiza germplasm resources.

Comprehensive analysis of sesame LRR-RLKs: structure, evolution and dynamic expression profiles under Macrophomina phaseolina stress

Leucine-rich repeat receptor-like kinases (LRR-RLKs) can participate in the regulation of plant growth and development, immunity and signal transduction. Sesamum indicum, one of the most important oil crops, has a significant role in promoting human health. In this study, 175 SiLRR-RLK genes were identified in S. indicum, and they were subdivided into 12 subfamilies by phylogenetic analysis. Gene duplication analysis showed that the expansion of the SiLRR-RLK family members in the sesame was mainly due to segmental duplication. Moreover, the gene expansion of subfamilies IV and III contributed to the perception of stimuli under M. phaseolina stress in the sesame. The collinearity analysis with other plant species revealed that the duplication of SiLRR-RLK genes occurred after the differentiation of dicotyledons and monocotyledons. The expression profile analysis and functional annotation of SiLRR-RLK genes indicated that they play a vital role in biotic stress. Furthermore, the protein-protein interaction and coexpression networks suggested that SiLRR-RLKs contributed to sesame resistance to Macrophomina phaseolina by acting alone or as a polymer with other SiLRR-RLKs. In conclusion, the comprehensive analysis of the SiLRR-RLK gene family provided a framework for further functional studies on SiLRR-RLK genes.

Comprehensive Characterization of the C3HC4 RING Finger Gene Family in Potato (Solanum tuberosum L.): Insights into Their Involvement in Anthocyanin Biosynthesis

The C3HC4 RING finger gene (RING-HC) family is a zinc finger protein crucial to plant growth. However, there have been no studies on the RING-HC gene family in potato. In this study, 77 putative StRING-HCs were identified in the potato genome and grouped into three clusters based on phylogenetic relationships, the chromosome distribution, gene structure, conserved motif, gene duplication events, and synteny relationships, and cis-acting elements were systematically analyzed. By analyzing RNA-seq data of potato cultivars, the candidate StRING-HC genes that might participate in tissue development, abiotic stress, especially drought stress, and anthocyanin biosynthesis were further determined. Finally, a StRING-HC gene (Soltu.DM.09G017280 annotated as StRNF4-like), which was highly expressed in pigmented potato tubers was focused on. StRNF4-like localized in the nucleus, and Y2H assays showed that it could interact with the anthocyanin-regulating transcription factors (TFs) StbHLH1 of potato tubers, which is localized in the nucleus and membrane. Transient assays showed that StRNF4-like repressed anthocyanin accumulation in the leaves of Nicotiana tabacum and Nicotiana benthamiana by directly suppressing the activity of the dihydroflavonol reductase (DFR) promoter activated by StAN1 and StbHLH1. The results suggest that StRNF4-like might repress anthocyanin accumulation in potato tubers by interacting with StbHLH1. Our comprehensive analysis of the potato StRING-HCs family contributes valuable knowledge to the understanding of their functions in potato development, abiotic stress, hormone signaling, and anthocyanin biosynthesis.

Genome-wide characterization and expression analysis of MADS-box transcription factor gene family in Perilla frutescens

MADS-box transcription factors are widely involved in the regulation of plant growth, developmental processes, and response to abiotic stresses. Perilla frutescens, a versatile plant, is not only used for food and medicine but also serves as an economical oil crop. However, the MADS-box transcription factor family in P. frutescens is still largely unexplored. In this study, a total of 93 PfMADS genes were identified in P. frutescens genome. These genes, including 37 Type I and 56 Type II members, were randomly distributed across 20 chromosomes and 2 scaffold regions. Type II PfMADS proteins were found to contain a greater number of motifs, indicating more complex structures and diverse functions. Expression analysis revealed that most PfMADS genes (more than 76 members) exhibited widely expression model in almost all tissues. The further analysis indicated that there was strong correlation between some MIKCC-type PfMADS genes and key genes involved in lipid synthesis and flavonoid metabolism, which implied that these PfMADS genes might play important regulatory role in the above two pathways. It was further verified that PfMADS47 can effectively mediate the regulation of lipid synthesis in Chlamydomonas reinhardtii transformants. Using cis-acting element analysis and qRT-PCR technology, the potential functions of six MIKCC-type PfMADS genes in response to abiotic stresses, especially cold and drought, were studied. Altogether, this study is the first genome-wide analysis of PfMADS. This result further supports functional and evolutionary studies of PfMADS gene family and serves as a benchmark for related P. frutescens breeding studies.

The Identification and Expression Analysis of the Liriodendron chinense F-Box Gene Family

The F-box gene family is one of the largest gene families in plants, and it plays a crucial role in regulating plant development, reproduction, cellular protein degradation, and response to biotic and abiotic stresses. Despite their significance, a comprehensive analysis of the F-box gene family in Liriodendron chinense and other magnoliaceae species has not been reported. In this study, we report for the first time the identification of 144 full-length F-box genes in L. chinense. Based on specific domains and phylogenetic analyses, these genes were divided into 10 distinct subfamilies. We further analyzed their gene structure, conserved domain and chromosome distribution, genome-wide replication events, and collinearity. Additionally, based on GO analysis, we found that F-box genes exhibit functional specificity, with a significant proportion of them being involved in protein binding (GO:0005515), suggesting that F-box genes may play an important role in gene regulation in L. chinense. Transcriptome data and q-PCR results also showed that F-box genes are involved in the development of multiple tissues in L. chinense, regulate the somatic embryogenesis of Liriodendron hybrids, and play a pivotal role in abiotic stress. Altogether, these findings provide a foundation for understanding the biological function of F-box genes in L. chinense and other plant species.

Identification of Differential Circular RNA Expression Profiles and Functional Networks in Human Macrophages Induced by Virulent and Avirulent Mycobacterium tuberculosis Strains

Circular RNAs (circRNAs) are noncoding RNAs with diverse functions. However, most Mycobacterium tuberculosis (M.tb)-related circRNAs remain undiscovered. In this study, we infected THP-1 cells with virulent and avirulent M.tb strains and then sequenced the cellular circRNAs. Bioinformatic analysis predicted 58,009 circRNAs in all the cells. In total, 2035 differentially expressed circRNAs were identified between the M.tb-infected and uninfected THP-1 cells and 1258 circRNAs were identified in the virulent and avirulent M.tb strains. Further, the top 10 circRNAs were confirmed by Sanger sequencing, among which four circRNAs, namely circSOD2, circCHSY1, circTNFRSF21, and circDHTKD1, which were highly differentially expressed in infected cells compared with those in uninfected cells, were further confirmed by ring formation, specific primers, and RNase R digestion. Next, circRNA-miRNA-mRNA subnetworks were constructed, such as circDHTKD1/miR-660-3p/IL-12B axis. Some of the individual downstream genes, such as miR-660-3p and IL-12B, were previously reported to be associated with cellular defense against pathological processes induced by M.tb infection. Because macrophages are important immune cells and the major host cells of M.tb, these findings provide novel ideas for exploring the M.tb pathogenesis and host defense by focusing on the regulation of circRNAs during M.tb infection.

Genome-wide identification of the C2H2 zinc finger gene family and expression analysis under salt stress in sweetpotato

Introduction

The higher plant transcription factor C2H2 zinc finger protein (C2H2-ZFP) is essential for plant growth, development, and stress response. There are limited studies on C2H2-ZFP genes in sweetpotato, despite a substantial number of C2H2-ZFP genes having been systematically found in plants.

Methods

In this work, 178 C2H2-ZFP genes were found in sweetpotato, distributed randomly on 15 chromosomes, and given new names according to where they were located. These members of the zinc finger gene family are separated into six branches, as shown by the phylogenetic tree. 24 tandem repeats of IbZFP genes and 46 fragment repeats were identified, and a homology study revealed that IbZFP genes linked more regions with wild relative species of sweetpotato as well as rhizome plants like potato and cassava. And we analyzed the expression patterns of IbZFP genes during the early development of sweetpotato storage roots (SRs) and salt stress using transcriptome data, and identified 44 IbZFP genes that exhibited differences in expression levels during the early expansion of sweetpotato SRs in different varieties, and 92 IbZFP genes that exhibited differences in expression levels under salt stress in salt tolerant and salt sensitive sweetpotato varieties. Additionally, we cloned six IbZFP genes in sweetpotato and analyzed their expression patterns in different tissues, their expression patterns under abiotic stress and hormone treatment, and subcellular localization.

Results and discussion

The results showed that the IbZFP genes had tissue specificity in sweetpotato and were induced to varying degrees by drought and salt stress. ABA and GA3 treatments also affected the expression of the IbZFP genes. We selected IbZFP105, which showed significant differences in expression levels under salt stress and ABA treatment, to be heterologously expressed in Arabidopsis thaliana. We found that IbZFP105 OE lines exhibited higher tolerance to salt stress and ABA stress. This indicates that IbZFP105 can enhance the salt tolerance of plants. These results systematically identified the evolution and expression patterns of members of the C2H2-ZFP gene family in sweetpotato, providing a theoretical basis for studying the role of IbZFP genes in the development of sweetpotato SRs and in resistance to stress.

A novel splice variant of goat CPT1a gene and their diverse mRNA expression profiles

The Alternative splicing (AS) of Carnitine palmitoyltransferase 1a (CPT1a) and their expression profiles had never been illuminated in goats until now. Herein, a novel splice transcript in the CPT1a gene that is predicted to result in the skipping of exons 6-19 (CPT1a-sv1) has been isolated in addition to the full-length transcript in goats. The result of RT-PCR showed that CPT1a-sv1 is 606 bp in length and consists of 6 exons. A novel exon 6 was consisted of partial exon 5 and partial exon 19, compared to that in CPT1a. RT-qPCR analysis showed that the expression patterns of CPT1a and CPT1a-sv1 are spatially different. In both kid and adult goats, the CPT1a transcript is strongly expressed in the liver, spleen, lung, kidney, and brain tissues. However, CPT1a-sv1 has a strong tissue-specific expression pattern, with moderate RNA levels in the liver and brain of kids, while highly expressed in the liver and minimally expressed in the brain of adults. We observed two transcripts to be involved in brain development. These findings improve our understanding of the function of the CPT1a gene in goats and provide information on the molecular mechanism of AS events.

Characterization of CrufCSP1 and Its Potential Involvement in Host Location by Cotesia ruficrus (Hymenoptera: Braconidae), an Indigenous Parasitoid of Spodoptera frugiperda (Lepidoptera: Noctuidae) in China

Chemosensory proteins (CSPs) are a class of soluble proteins that facilitate the recognition of chemical signals in insects. While CSP genes have been identified in many insect species, studies investigating their function remain limited. Cotesia ruficrus (Hymenoptera: Braconidae) holds promise as an indigenous biological control agent for managing the invasive pest Spodoptera frugiperda (Lepidoptera: Noctuidae) in China. This study aimed to shed light on the gene expression, ligand binding, and molecular docking of CrufCSP1 in C. ruficrus. A RT-qPCR analysis revealed that the expression of CrufCSP1 was higher in the wings, with male adults exhibiting significantly higher relative expression levels than other developmental stages. A fluorescence competitive binding analysis further demonstrated that CrufCSP1 has a high binding ability with several host-related volatiles, with trans-2-hexenal, octanal, and benzaldehyde showing the strongest affinity to CrufCSP1. A molecular docking analysis indicated that specific amino acid residues (Phe24, Asp25, Thr53, and Lys81) of CrufCSP1 can bind to these specific ligands. Together, these findings suggest that CrufCSP1 may play a crucial role in the process of C. ruficrus locating hosts. This knowledge can contribute to the development of more efficient and eco-friendly strategies for protecting crops and managing pests.

Differential Expression of miRNAs, lncRNAs, and circRNAs between Ovaries and Testes in Common Carp (Cyprinus carpio)

Female common carp grow faster than male individuals, implying that rearing females could be more profitable in aquaculture. Non-coding RNAs (ncRNAs) serve as versatile regulators with multiple functions in diverse biological processes. However, the roles of ncRNAs in the sex differentiation of common carp are less studied. In this study, we investigated the expression profiles of ncRNAs, including miRNAs, lncRNAs, and circRNAs, in the gonads to comprehend the roles of ncRNAs in sex differentiation in common carp. A substantial number of differentially expressed (DE) ncRNAs in ovaries and testes were identified. Some miRNAs, notably miR-205, miR-214, and miR-460-5p, might modulate hormone synthesis and thus maintain sex. A novel miRNA, novel_158, was predicted to suppress the expression of foxl3. DE lncRNAs were associated with oocyte meiosis, GnRH signaling pathways, and steroid biosynthesis, while DE circRNA target genes were enriched in the WNT signaling pathway and MAPK signaling pathway. We also analyzed ncRNA-mRNA interactions to shed light on the crosstalk between competing endogenous RNAs (ceRNAs), which is the critical mechanism by which lncRNAs and circRNAs function. Some lncRNAs and circRNAs may be able to competitively bind novel_313, a new miRNA, and thus regulate hsd17β3. Our research will provide a valuable resource for understanding the genetic basis of gonadal differentiation and development in common carp.

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

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

Comparative transcriptomics analysis of testicular miRNA from indicine and taurine cattle

Numerous studies have shown that several microRNAs (miRNAs) are specifically expressed in testis, play an essential role in regulating testicular spermatogenesis. Hainan and Mongolian cattle are two representative Chinese native cattle breeds representing Bos indicus (indicine cattle) and Bos taurus (taurine cattle), respectively, which are distributed in hot Hainan and cold Inner Mongolia province. To study the functional differences of miRNA in spermatogenesis between indicine and taurine cattle, six mature testes samples from indicine cattle (n = 3) and taurine cattle (n = 3) were collected, respectively. We detected miRNA expression using small RNA sequencing technology following bioinformatic analysis. A total of 578 known miRNAs and 132 novel miRNAs were detected in the six libraries. Among the 710 miRNAs, 564 miRNAs were expressed in both indicine and taurine cattle, 73 miRNAs were found solely in indicine cattle and 73 miRNAs were found solely in taurine cattle. After further analysis, among the miRNAs were identified in both indicine and taurine cattle, 184 miRNAs were differentially expressed (|log2 fold change| ≥ 1 and corrected p-value <0.05). Among the miRNAs that were only expressed in indicine cattle, 10 miRNAs were differentially expressed, whereas, among the miRNAs that were only expressed in taurine cattle, six miRNAs were differentially expressed. The enrichment analysis result showed that predicted target genes of a total of 200 differentially expressed miRNAs were enriched on some testicular spermatogenesis-related Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, especially mitogen-activated protein kinase (MAPK) signaling pathway. These findings identify miRNAs as key factors to regulate spermatogenesis in both indicine and taurine cattle, which may also be helpful for improving cattle reproductive performance in future studies.

Genome-wide identification of PEBP gene family in pineapple reveal its potential functions in flowering

Phosphatidylethanolamine binding protein (PEBP) plays an important role in regulating flowering time and morphogenesis of plants. However, the identification and functional analysis of PEBP gene in pineapple (AcPEBP) have not been systematically studied. The pineapple genome contained 11 PEBP family members, which were subsequently classified into three subfamilies (FT-like, TFL-like and MFT-like) based on phylogenetic relationships. The arrangement of these 11 shows an unequal pattern across the six chromosomes of pineapple the pineapple genome. The anticipated outcomes of the promoter cis-acting elements indicate that the PEBP gene is subject to regulation by diverse light signals and endogenous hormones such as ethylene. The findings from transcriptome examination and quantitative real-time polymerase chain reaction (qRT-PCR) indicate that FT-like members AcFT3 and AcFT4 display a heightened expression level, specifically within the floral structures. The expression of AcFT3 and AcFT4 increases sharply and remains at a high level after 4 days of ethylene induction, while the expression of AcFT7 and AcMFT1 decreases gradually during the flowering process. Additionally, AcFT3, AcFT4 and AcFT7 show specific expression in different floral organs of pineapple. These outcomes imply that members belonging to the FT-like subfamily may have a significant impact on the process of bud differentiation and flower development. Through transcriptional activation analysis, it was determined that AcFT4 possesses transcriptional activation capability and is situated in the nucleus and peripheral cytoplasm. Overexpression of AcFT4 in Arabidopsis resulted in the promotion of early flowering by 6-7 days. The protein interaction prediction network identified potential flower regulators, including CO, AP1, LFY and SOC1, that may interact with PEBP proteins. This study explores flower development in pineapple, thereby serving as a valuable reference for future research endeavors in this domain.

Genome-Wide Identification and Characterization of the NAC Gene Family and Its Involvement in Cold Response in Dendrobium officinale

The NAC (NAM, ATAF1/2 and CUC2) gene family is one of the largest plant-specific transcription factor families, functioning as crucial regulators in diverse biological processes such as plant growth and development as well as biotic and abiotic stress responses. Although it has been widely characterized in many plants, the significance of the NAC family in Dendrobium officinale remained elusive up to now. In this study, a genome-wide search method was conducted to identify NAC genes in Dendrobium officinale (DoNACs) and a total of 110 putative DoNACs were obtained. Phylogenetic analysis classified them into 15 subfamilies according to the nomenclature in Arabidopsis and rice. The members in the subfamilies shared more similar gene structures and conversed protein domain compositions. Furthermore, the expression profiles of these DoNACs were investigated in diverse tissues and under cold stress by RNA-seq data. Then, a total of five up-regulated and five down-regulated, cold-responsive DoNACs were validated through QRT-PCR analysis, demonstrating they were involved in regulating cold stress response. Additionally, the subcellular localization of two down-regulated candidates (DoNAC39 and DoNAC58) was demonstrated to be localized in the nuclei. This study reported the genomic organization, protein domain compositions and expression patterns of the NAC family in Dendrobium officinale, which provided targets for further functional studies of DoNACs and also contributed to the dissection of the role of NAC in regulating cold tolerance in Dendrobium officinale.

Spatiotemporal Patterns of Five Small Heat Shock Protein Genes in Hyphantria cunea in Response to Thermal Stress

Hyphantria cunea (Drury), a destructive polyphagous pest, has been spreading southward after invading northern China, which indicates that this insect species is facing a huge thermal challenge. Small heat shock proteins (sHSPs) function as ATP-independent molecular chaperones that protect insects from heat stress damage. In order to explore the role of sHSPs in the thermotolerance of H. cunea, five novel sHSP genes of H. cunea were cloned, including an orthologous gene (HcHSP21.4) and four species-specific sHSP genes (HcHSP18.9, HcHSP20.1, HcHSP21.5, and HcHSP29.8). Bioinformatics analysis showed that the proteins encoded by these five HcHSPs contained typical α-crystallin domains. Quantitative real-time PCR analysis revealed the ubiquitous expression of all HcHSPs across all developmental stages of H. cunea, with the highest expression levels in pupae and adults. Four species-specific HcHSPs were sensitive to high temperatures. The expression levels of HcHSPs were significantly up-regulated under heat stress and increased with increasing temperature. The expression levels of HcHSPs in eggs exhibited an initial up-regulation in response to a temperature of 40 °C. In other developmental stages, the transcription of HcHSPs was immediately up-regulated at 30 °C or 35 °C. HcHSPs transcripts were abundant in the cuticle before and after heat shock. The expression of HcHSP21.4 showed weak responses to heat stress and constitutive expression in the tissues tested. These results suggest that most of the HcHSPs are involved in high-temperature response and may also have functions in the normal development and reproduction of H. cunea.

Transcription factors NF-YB involved in embryogenesis and hormones responses in Dimocarpus Longan Lour

Introduction

NF-YB transcription factor is an important regulatory factor in plant embryonic development.

Results

In this study, 15 longan NF-YB (DlNF-YB) family genes were systematically identified in the whole genome of longan, and a comprehensive bioinformatics analysis of DlNF-YB family was performed. Comparative transcriptome analysis of DlNF-YBs expression in different tissues, early somatic embryogenesis (SE), and under different light and temperature treatments revealed its specific expression profiles and potential biological functions in longan SE. The qRT-PCR results implied that the expression patterns of DlNF-YBs were different during SE and the zygotic embryo development of longan. Supplementary 2,4-D, NPA, and PP333 in longan EC notably inhibited the expression of DlNF-YBs; ABA, IAA, and GA3 suppressed the expressions of DlNF-YB6 and DlNF-YB9, but IAA and GA3 induced the other DlNF-YBs. Subcellular localization indicated that DlNF-YB6 and DlNF-YB9 were located in the nucleus. Furthermore, verification by the modified 5'RNA Ligase Mediated Rapid Amplification of cDNA Ends (5' RLM-RACE) method demonstrated that DlNF-YB6 was targeted by dlo-miR2118e, and dlo-miR2118e regulated longan somatic embryogenesis (SE) by targeting DlNF-YB6. Compared with CaMV35S- actuated GUS expression, DlNF-YB6 and DlNF-YB9 promoters significantly drove GUS expression. Meanwhile, promoter activities were induced to the highest by GA3 but suppressed by IAA. ABA induced the activities of the promoter of DlNF-YB9, whereas it inhibited the promoter of DlNF-YB6.

Discussion

Hence, DlNF-YB might play a prominent role in longan somatic and zygotic embryo development, and it is involved in complex plant hormones signaling pathways.

Genome-Wide Identification and Analysis of TCP Gene Family among Three Dendrobium Species

Dendrobium orchids, which are among the most well-known species of orchids, are appreciated for their aesthetic appeal across the globe. Furthermore, due to their strict living conditions, they have accumulated high levels of active ingredients, resulting not only in their medicinal value but also in their strong ability to respond to harsh environments. The TCP gene family plays an important role in plant growth and development, and signal transduction. However, these genes have not been systematically investigated in Dendrobium species. In this study, we detected a total of 24, 23, and 14 candidate TCP members in the genome sequences of D. officinale, D. nobile, and D. chrysotoxum, respectively. These genes were classified into three clades on the basis of a phylogenetic analysis. The TCP gene numbers among Dendrobium species were still highly variable due to the independent loss of genes in the CIN clade. However, only three gene duplication events were detected, with only one tandem duplication event (DcTCP9/DcTCP10) in D. chrysotoxum and two pairs of paralogous DoTCP gene duplication events (DoTCP1/DoTCP23 and DoTCP16/DoTCP24) in D. officinale. A total of 25 cis-acting elements of TCPs related to hormone/stress and light responses were detected. Among them, the proportions of hormone response, light response, and stress response elements in D. officinale (100/421, 127/421, and 171/421) were similar to those in D. nobile (83/352, 87/352, and 161/352). Using qRT-PCR to determine their expression patterns under MeJA treatment, four DoTCPs (DoTCP2, DoTCP4, DoTCP6, and DoTCP14) were significantly upregulated under MeJA treatment, which indicates that TCP genes may play important roles in responding to stress. Under ABA treatment, seven DoTCPs (DoTCP3, DoTCP7, DoTCP9, DoTCP11, DoTCP14, DoTCP15, and DoTCP21) were significantly upregulated, indicating that TCP genes may also play an important role in hormone response. Therefore, these results can provide useful information for studying the evolution and function of TCP genes in Dendrobium species.

Phylogenetic Analysis of the PR-4 Gene Family in Euphorbiaceae and Its Expression Profiles in Tung Tree (Vernicia fordii)

Pathogenesis-related protein-4 (PR-4) is generally believed to be involved in physiological processes. However, a comprehensive investigation of this protein in tung tree (Vernicia fordii) has yet to be conducted. In this study, we identified 30 PR-4 genes in the genomes of Euphorbiaceae species and investigated their domain organization, evolution, promoter cis-elements, expression profiles, and expression profiles in the tung tree. Sequence and structural analyses indicated that VF16136 and VF16135 in the tung tree could be classified as belonging to Class II and I, respectively. Phylogenetic and Ka/Ks analyses revealed that Hevea brasiliensis exhibited a significantly expanded number of PR-4 genes. Additionally, the analysis of promoter cis-elements suggested that two VfPR-4 genes may play a role in the response to hormones and biotic and abiotic stress of tung trees. Furthermore, the expression patterns of VfPR-4 genes and their responses to 6-BA, salicylic acid, and silver nitrate in inflorescence buds of tung trees were evaluated using qRT-PCR. Notably, the expression of two VfPR-4 genes was found to be particularly high in leaves and early stages of tung seeds. These results suggest that VF16136 and VF16135 may have significant roles in the development of leaves and seeds in tung trees. Furthermore, these genes were found to be responsive to 6-BA, salicylic acid, and silver nitrate in the development of inflorescence buds. This research provides valuable insights for future investigation into the functions of PR-4 genes in tung trees.

Identification of GA2ox Family Genes and Expression Analysis under Gibberellin Treatment in Pineapple (Ananas comosus (L.) Merr.)

Gibberellin (GAs) plays an important regulatory role in the development and growth of pineapple (Ananas comosus (L.) Merr.). Bioinformatics was used to confirm the differential expression of GA2 gibberellin oxidase gene AcGA2oxs in the pineapple genome, which laid the foundation for exploring its role in pineapple. In this study, 42 GA2ox genes (AcGA2oxs) were identified in the pineapple genome, named from AcGA2ox1 to AcGA2ox42, and divided into four groups according to phylogenetic analysis. We also analyzed the gene structure, conserved motifs and chromosome localization of AcGA2oxs. AcGA2oxs within the same group had similar gene structure and motifs composition. Collinear analysis and cis-element analysis provided the basis for understanding the evolution and function of GA2ox genes in pineapple. In addition, we selected different tissue parts to analyze the expression profile of AcGA2oxs, and the results show that 41 genes were expressed, except for AcGA2ox18. AcGA2ox18 may not be expressed in these sites or may be pseudogenes. qRT-PCR (real-time fluorescence quantitative PCR) was used to detect the relative expression levels of the GA2ox gene family under different concentrations of GA3 treatment, and it was found that AcGA2ox gene expression was upregulated in different degrees under GA3 treatment. These results provide useful information for further study on the evolution and function of the GA2ox family in pineapple.

Genome-Wide Detection of SPX Family and Profiling of CoSPX-MFS3 in Regulating Low-Phosphate Stress in Tea-Oil Camellia

Camellia oleifera a member of the family Theaceae, is a phosphorus (P) tolerator native to southern China. The SPX gene family critically regulates plant growth and development and maintains phosphate (Pi) homeostasis. However, the involvement of SPX genes in Pi signaling in Tea-Oil Camellia remains unknown. In this work, 20 SPX genes were identified and categorized into four subgroups. Conserved domains, motifs, gene structure, chromosomal location and gene duplication events were also investigated in the SPX gene family. Defense and stress responsiveness cis-elements were identified in the SPX gene promoters, which participated in low-Pi stress responses. Based on transcriptome data and qRT-PCR results, nine CoSPX genes had similar expression patterns and eight genes (except CoPHO1H3) were up-regulated at 30 days after exposure to low-Pi stress. CoSPX-MFS3 was selected as a key candidate gene by WGCNA analysis. CoSPX-MFS3 was a tonoplast protein. Overexpression of CoSPX-MFS3 in Arabidopsis promoted the accumulation of total P content and decreased the anthocyanin content. Overexpression of CoSPX-MFS3 could enhance low-Pi tolerance by increased biomass and organic acid contents in transgenic Arabidopsis lines. Furthermore, the expression patterns of seven phosphate starvation genes were higher in transgenic Arabidopsis than those in the wild type. These results highlight novel physiological roles of the SPX family genes in C. oleifera under low-Pi stress, and lays the foundation for a deeper knowledge of the response mechanism of C. oleifera to low-Pi stress.

Expression profiles and potential roles of serum tRNA‑derived fragments in diabetic nephropathy

Diabetic nephropathy (DN) is one of the most important causes of end-stage renal disease and current treatments are ineffective in preventing its progression. Transfer RNA (tRNA)-derived fragments (tRFs), which are small non-coding fragments derived from tRNA precursors or mature tRNAs, have a critical role in various human diseases. The present study aimed to investigate the expression profile and potential functions of tRFs in DN. High-throughput sequencing technology was employed to detect the differential serum levels of tRFs between DN and diabetes mellitus and to validate the reliability of the sequencing results using reverse transcription-quantitative PCR. Ultimately, six differentially expressed (DE) tRFs were identified (P<0.05; |log2fold change| ≥1), including three upregulated (tRF5-GluCTC, tRF5-AlaCGC and tRF5-ValCAC) and three downregulated tRFs (tRF5-GlyCCC, tRF3-GlyGCC and tRF3-IleAAT). Potential functions and regulatory mechanisms of these DE tRFs were further evaluated using an applied bioinformatics-based analysis. Gene ontology analysis revealed that the DE tRFs are mainly enriched in biological processes, including axon guidance, Rad51 paralog (Rad51)B-Rad51C-Rad51D-X-Ray repair cross-complementing 2 complex, nuclear factor of activated T-cells protein binding and fibroblast growth factor-activated receptor activity. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that they are associated with axon guidance, neurotrophin signaling, mTOR signaling, AMPK signaling and epidermal growth factor receptor family signaling pathways. In conclusion, the present findings indicated that tRFs were DE in DN and may be involved in the regulation of DN pathology through multiple pathways, thereby providing a new perspective for the study of DN therapeutic targets.

Genome-wide characterization of SOS1 gene family in potato (Solanum tuberosum) and expression analyses under salt and hormone stress

Salt Overly Sensitive 1 (SOS1) is one of the members of the Salt Overly Sensitive (SOS) signaling pathway and plays critical salt tolerance determinant in plants, while the characterization of the SOS1 family in potato (Solanum tuberosum) is lacking. In this study, 37 StSOS1s were identified and found to be unevenly distributed across 10 chromosomes, with most of them located on the plasma membrane. Promoter analysis revealed that the majority of these StSOS1 genes contain abundant cis-elements involved in various abiotic stress responses. Tissue specific expression showed that 21 of the 37 StSOS1s were widely expressed in various tissues or organs of the potato. Molecular interaction network analysis suggests that 25 StSOS1s may interact with other proteins involved in potassium ion transmembrane transport, response to salt stress, and cellular processes. In addition, collinearity analysis showed that 17, 8, 1 and 5 of orthologous StSOS1 genes were paired with those in tomato, pepper, tobacco, and Arabidopsis, respectively. Furthermore, RT-qPCR results revealed that the expression of StSOS1s were significant modulated by various abiotic stresses, in particular salt and abscisic acid stress. Furthermore, subcellular localization in Nicotiana benthamiana suggested that StSOS1-13 was located on the plasma membrane. These results extend the comprehensive overview of the StSOS1 gene family and set the stage for further analysis of the function of genes in SOS and hormone signaling pathways.

Molecular Evolution and Genetic Variation of G2-Like Transcription Factor Genes in Wheat (Triticum aestivum L.)

The GOLDEN2-LIKE (G2-like) gene family members provide significant contributions to the growth and development of plants. In this study, a total of 76 wheat G2-like gene family members (TaG1-TaG76) were detected in the wheat genome and were categorized into three groups (including six subgroups) based on the gene structure and protein motif analyses. These genes were unevenly distributed in 19 of 21 wheat chromosomes. A total of 63 segmental duplication pairs of TaG2-like genes were identified in the wheat genome. The expression levels of all the TaG2-like genes indicated that TaG2-like genes showed different expression patterns in various organs and tissues. Moreover, the transcriptions of TaG2-like genes were significantly affected under abiotic stress (cold, ABA, NaCl, and PEG). This study offered valuable insights into the functional characterization of TaG2-like genes in wheat.

Genome-wide investigation and expression profiling of LOR gene family in rapeseed under salinity and ABA stress

The Brassica napus (B. napus) LOR (Lurp-One-Related) gene family is a little-known gene family characterized by a conserved LOR domain in the proteins. Limited research in Arabidopsis showed that LOR family members played important roles in Hyaloperonospora parasitica (Hpa) defense. Nevertheless, there is a paucity of research investigating the role of the LOR gene family towards their responses to abiotic stresses and hormone treatments. This study encompassed a comprehensive survey of 56 LOR genes in B. napus, which is a prominent oilseed crop that holds substantial economic significance in China, Europe, and North America. Additionally, the study evaluated the expression profiles of these genes in response to salinity and ABA stress. Phylogenetic analysis showed that 56 BnLORs could be divided into 3 subgroups (8 clades) with uneven distribution on 19 chromosomes. 37 out of 56 BnLOR members have experienced segmental duplication and 5 of them have undergone tandem repeats events with strong evidence of purifying selection. Cis-regulatory elements (CREs) analysis indicated that BnLORs involved in process such as light response, hormone response, low temperature response, heat stress response, and dehydration response. The expression pattern of BnLOR family members revealed tissue specificity. RNA-Seq and qRT-PCR were used to validate BnLOR gene expression under temperature, salinity and ABA stress, revealing that most BnLORs showed inducibility. This study enhanced our comprehension of the B. napus LOR gene family and could provide valuable information for identifying and selecting genes for stress resistant breeding.

Comparative expression profiles of carboxylesterase orthologous CXE14 in two closely related tea geometrid species, Ectropis obliqua Prout and Ectropis grisescens Warren

Insect carboxylesterases (CXEs) can be expressed in multiple tissues and play crucial roles in detoxifying xenobiotic insecticides and degrading olfactory cues. Therefore, they have been considered as an important target for development of eco-friendly insect pest management strategies. Despite extensive investigation in most insect species, limited information on CXEs in sibling moth species is currently available. The Ectropis obliqua Prout and Ectropis grisescens Warren are two closely related tea geometrid species, which share the same host of tea plant but differ in geographical distribution, sex pheromone composition, and symbiotic bacteria abundance, providing an excellent mode species for studies of functional diversity of orthologous CXEs. In this study, we focused on EoblCXE14 due to its previously reported non-chemosensory organs-biased expression. First, the EoblCXE14 orthologous gene EgriCXE14 was cloned and sequence characteristics analysis showed that they share a conserved motif and phylogenetic relationship. Quantitative real-time polymerase chain reaction (qRT-PCR) was then used to compare the expression profiles between two Ectropis spp. The results showed that EoblCXE14 was predominately expressed in E. obliqua larvae, whereas EgriCXE14 was abundant in E. grisescens at multiple developmental stages. Interestingly, both orthologous CXEs were highly expressed in larval midgut, but the expression level of EoblCXE14 in E. obliqua midgut was significantly higher than that of EgriCXE14 in E. grisescens midgut. In addition, the potential effect of symbiotic bacteria Wolbachia on the CXE14 was examined. This study is the first to provide comparative expression profiles of orthologous CXE genes in two sibling geometrid moth species and the results will help further elucidate CXEs functions and identify a potential target for tea geometrid pest control.

[Identification of heat shock protein hsp70 family genes from Rana amurensis and its expression profiles upon infection]

Heat shock proteins (HSPs) widely exist in all organisms, the structures of which are usually extraordinarily conservative. They are also well-known stress proteins that are involved in response to physical, chemical and biological stresses. HSP70 is an important member of the HSPs family. In order to study the roles of amphibians HSP70 during infection, the cDNA sequence of Rana amurensis hsp70 family genes were cloned by homologous cloning method. The sequence characteristics, three-dimensional structure and genetic relationship of Ra-hsp70s were analyzed by bioinformatics methods. The expression profiles under bacterial infection were also analyzed by real-time quantitative PCR (qRT-PCR). Expression and localization of HSP70 protein were tested by immunohistochemical techniques. The results showed that three conservative tag sequences of HSP70 family, HSPA5, HSPA8 and HSPA13, were found in HSP70. Phylogenetic tree analysis indicated four members are distributed in four different branches, and members with the same subcellular localization motif are distributed in the same branch. The relative expression levels of the mRNA of four members were all significantly upregulated (P < 0.01) upon infection, but the time for up-regulating the expression levels were diverse in different tissues. The immunohistochemical analysis showed that HSP70 was expressed to different degrees in the cytoplasm of liver, kidney, skin and stomach tissue. The four members of Ra-hsp70 family have ability to respond bacterial infection to varying degrees. Therefore, it was proposed that they are involved in biological processes against pathogen and play different biological functions. The study provides a theoretical basis for functional studies of HSP70 gene in amphibians.

Genome-Wide Identification and Expression Pattern Analysis of Dirigent Members in the Genus Oryza

Dirigent (DIR) members have been shown to play essential roles in plant growth, development and adaptation to environmental changes. However, to date, there has been no systematic analysis of the DIR members in the genus Oryza. Here, 420 genes were identified from nine rice species to have the conserved DIR domain. Importantly, the cultivated rice species Oryza sativa has more DIR family members than the wild rice species. DIR proteins in rice could be classified into six subfamilies based on phylogeny analysis. Gene duplication event analysis suggests that whole genome/segmental duplication and tandem duplication are the primary drivers for DIR genes' evolution in Oryza, while tandem duplication is the main mechanism of gene family expansion in the DIR-b/d and DIR-c subfamilies. Analysis of the RNA sequencing data indicates that OsjDIR genes respond to a wide range of environmental factors, and most OsjDIR genes have a high expression level in roots. Qualitative reverse transcription PCR assays confirmed the responsiveness of OsjDIR genes to the undersupply of mineral elements, the excess of heavy metals and the infection of Rhizoctonia solani. Furthermore, there exist extensive interactions between DIR family members. Taken together, our results shed light on and provide a research foundation for the further exploration of DIR genes in rice.

Genome-Wide Identification and Chilling Stress Analysis of the NF-Y Gene Family in Melon

The nuclear factor Y (NF-Y) transcription factor contains three subfamilies: NF-YA, NF-YB, and NF-YC. The NF-Y family have been reported to be key regulators in plant growth and stress responses. However, little attention has been given to these genes in melon (Cucumis melo L.). In this study, twenty-five NF-Ys were identified in the melon genome, including six CmNF-YAs, eleven CmNF-YBs, and eight CmNF-YCs. Their basic information (gene location, protein characteristics, and subcellular localization), conserved domains and motifs, and phylogeny and gene structure were subsequently analyzed. Results showed highly conserved motifs exist in each subfamily, which are distinct between subfamilies. Most CmNF-Ys were expressed in five tissues and exhibited distinct expression patterns. However, CmNF-YA6, CmNF-YB1/B2/B3/B8, and CmNF-YC6 were not expressed and might be pseudogenes. Twelve CmNF-Ys were induced by cold stress, indicating the NF-Y family plays a key role in melon cold tolerance. Taken together, our findings provide a comprehensive understanding of CmNF-Y genes in the development and stress response of melon and provide genetic resources for solving the practical problems of melon production.

Genome-Wide Analysis of the LATERAL ORGAN BOUNDARIES Domain (LBD) Members in Alfalfa and the Involvement of MsLBD48 in Nitrogen Assimilation

The LATERAL ORGAN BOUNDARIES DOMAIN (LBD) proteins, a transcription factor family specific to the land plants, have been implicated in multiple biological processes including organ development, pathogen response and the uptake of inorganic nitrogen. The study focused on LBDs in legume forage Alfalfa. The genome-wide analysis revealed that in Alfalfa 178 loci across 31 allelic chromosomes encoded 48 unique LBDs (MsLBDs), and the genome of its diploid progenitor M. sativa spp. Caerulea encoded 46 LBDs. Synteny analysis indicated that the expansion of AlfalfaLBDs was attributed to the whole genome duplication event. The MsLBDs were divided into two major phylogenetic classes, and the LOB domain of the Class I members was highly conserved relative to that of the Class II. The transcriptomic data demonstrated that 87.5% of MsLBDs were expressed in at least one of the six test tissues, and Class II members were preferentially expressed in nodules. Moreover, the expression of Class II LBDs in roots was upregulated by the treatment of inorganic nitrogen such as KNO₃ and NH₄Cl (0.3 mM). The overexpression of MsLBD48, a Class II member, in Arabidopsis resulted in growth retardance with significantly declined biomass compared with the non-transgenic plants, and the transcription level of the genes involved in nitrogen uptake or assimilation, including NRT1.1, NRT2.1, NIA1 and NIA2 was repressed. Therefore, the LBDs in Alfalfa are highly conserved with their orthologs in embryophytes. Our observations that ectopic expression of MsLBD48 inhibited Arabidopsis growth by repressing nitrogen adaption suggest the negative role of the transcription factor in plant uptake of inorganic nitrogen. The findings imply the potential application of MsLBD48 in Alfalfa yield improvement via gene editing.

New Low Morphine Opium Poppy Genotype Obtained by TILLING Approach

The opium poppy's ability to produce various alkaloids is both useful and problematic. Breeding of new varieties with varying alkaloid content is therefore an important task. In this paper, the breeding technology of new low morphine poppy genotypes, based on a combination of a TILLING approach and single-molecule real-time NGS sequencing, is presented. Verification of the mutants in the TILLING population was obtained using RT-PCR and HPLC methods. Only three of the single-copy genes of the morphine pathway among the eleven genes were used for the identification of mutant genotypes. Point mutations were obtained only in one gene (CNMT) while an insertion was obtained in the other (SalAT). Only a few expected transition SNPs from G:C to A:T were obtained. In the low morphine mutant genotype, the production of morphine was decreased to 0.1% from 1.4% in the original variety. A comprehensive description of the breeding process, a basic characterization of the main alkaloid content, and a gene expression profile for the main alkaloid-producing genes is provided. Difficulties with the TILLING approach are also described and discussed.

Evolutionary Analysis of Respiratory Burst Oxidase Homolog (RBOH) Genes in Plants and Characterization of ZmRBOHs

The respiratory burst oxidase homolog (RBOH), as the key producer of reactive oxygen species (ROS), plays an essential role in plant development. In this study, a bioinformatic analysis was performed on 22 plant species, and 181 RBOH homologues were identified. A typical RBOH family was identified only in terrestrial plants, and the number of RBOHs increased from non-angiosperms to angiosperms. Whole genome duplication (WGD)/segmental duplication played a key role in RBOH gene family expansion. Amino acid numbers of 181 RBOHs ranged from 98 to 1461, and the encoded proteins had molecular weights from 11.1 to 163.6 kDa, respectively. All plant RBOHs contained a conserved NADPH_Ox domain, while some of them lacked the FAD_binding_8 domain. Plant RBOHs were classified into five main subgroups by phylogenetic analysis. Most RBOH members in the same subgroup showed conservation in both motif distribution and gene structure composition. Fifteen ZmRBOHs were identified in maize genome and were positioned in eight maize chromosomes. A total of three pairs of orthologous genes were found in maize, including ZmRBOH6/ZmRBOH8, ZmRBOH4/ZmRBOH10 and ZmRBOH15/ZmRBOH2. A Ka/Ks calculation confirmed that purifying selection was the main driving force in their evolution. ZmRBOHs had typical conserved domains and similar protein structures. cis-element analyses together with the expression profiles of the ZmRBOH genes in various tissues and stages of development suggested that ZmRBOH was involved in distinct biological processes and stress responses. Based on the RNA-Seq data and qRT-PCR analysis, the transcriptional response of ZmRBOH genes was examined under various abiotic stresses, and most of ZmRBOH genes were up-regulated by cold stress. These findings provide valuable information for further revealing the biological roles of ZmRBOH genes in plant development and abiotic stress responses.

Identification of QTLs and critical genes related to sugarcane mosaic disease resistance

Mosaic viral diseases affect sugarcane productivity worldwide. Mining disease resistance-associated molecular markers or genes is a key component of disease resistance breeding programs. In the present study, 285 F₁ progeny were produced from a cross between Yuetang 93-159, a moderately resistant variety, and ROC22, a highly susceptible variety. The mosaic disease symptoms of these progenies, with ROC22 as the control, were surveyed by natural infection under 11 different environmental conditions in the field and by artificial infections with a mixed sugarcane mosaic virus (SCMV) and sorghum mosaic virus (SrMV) inoculum. Analysis of consolidated survey data enabled the identification of 29 immune, 55 highly resistant, 70 moderately resistant, 62 susceptible, and 40 highly susceptible progenies. The disease response data and a high-quality SNP genetic map were used in quantitative trait locus (QTL) mapping. The results showed that the correlation coefficients (0.26~0.91) between mosaic disease resistance and test environments were significant (p< 0.001), and that mosaic disease resistance was a highly heritable quantitative trait (H² = 0.85). Seven mosaic resistance QTLs were located to the SNP genetic map, each QTL accounted for 3.57% ~ 17.10% of the phenotypic variation explained (PVE). Furthermore, 110 pathogen response genes and 69 transcription factors were identified in the QTLs interval. The expression levels of nine genes (Soffic.07G0015370-1P, Soffic.09G0015410-2T, Soffic.09G0016460-1T, Soffic.09G0016460-1P, Soffic.09G0017080-3C, Soffic.09G0018730-3P, Soffic.09G0018730-3C, Soffic.09G0019920-3C and Soffic.03G0019710-2C) were significantly different between resistant and susceptible progenies, indicating their key roles in sugarcane resistance to SCMV and SrMV infection. The seven QTLs and nine genes can provide a certain scientific reference to help sugarcane breeders develop varieties resistant to mosaic diseases.

Identification and characterization of UDP-glycosyltransferase genes and the potential role in response to insecticides exposure in Bactrocera dorsalis

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis (Hendel) is a worldwide pest damaging a wide range of hosts. Due to the long-term indiscriminate use of insecticides, B. dorsalis has developed serious resistance to several insecticides. UDP-glycosyltransferases (UGTs) are secondary metabolic enzymes involved in biotransformation and play an important role in the metabolism of plant secondary metabolites and synthetic insecticides in insects. Thus, we suspect that UGTs in B. dorsalis play an important role in insecticide tolerance.

RESULTS

In this study, 31 UGT genes were identified in the genome of B. dorsalis, belonging to 13 subfamilies. Real-time quantitative polymerase chain reaction (RT-qPCR) results revealed that 12 UGT genes were highly expressed in the antennae, midgut, Malpighian tubule and fat body. The mRNA expressions of 17 UGT genes were up-regulated upon exposure to λ-cyhalothrin, imidacloprid, abamectin and chlorpyrifos. Knockdown of the selected five UGT genes (BdUGT301D2, BdUGT35F2, BdUGT36K2, BdUGT49D2, BdUGT50B5) by RNA interference increased the mortality of B. dorsalis from 9.29% to 27.22% upon exposure to four insecticides.

CONCLUSION

The abundance of UGTs in B. dorsalis is similar to other insect species, and 12 out of 31 UGTs were specifically expressed in metabolic tissues, suggesting a key role in detoxification. Down-regulation of five selected UGT genes increased the susceptibility of B. dorsalis to various insecticides, indicating that UGTs may play an important role in tolerance of B. dorsalis to multiple insecticides. © 2022 Society of Chemical Industry.

Genome-Wide Identification, Expression Analysis, and Potential Roles under Abiotic Stress of the YUCCA Gene Family in Mungbean (Vigna radiata L.)

YUCCA, belonging to the class B flavin-dependent monooxygenases, catalyzes the rate-limiting step for endogenous auxin synthesis and is implicated in plant-growth regulation and stress response. Systematic analysis of the YUCCA gene family and its stress response benefits the dissection of regulation mechanisms and breeding applications. In this study, 12 YUCCA genes were identified from the mungbean (Vigna radiata L.) genome and were named based on their similarity to AtYUCCAs. Phylogenetic analysis revealed that the 12 VrYUCCAs could be divided into 4 subfamilies. The evidence from enzymatic assays in vitro and transgenetic Arabidopsis in vivo indicated that all the isolated VrYUCCAs had biological activity in response to IAA synthesis. Expression pattern analysis showed that functional redundancy and divergence existed in the VrYUCCA gene family. Four VrYUCCAs were expressed in most tissues, and five VrYUCCAs were specifically highly expressed in the floral organs. The response toward five stresses, namely, auxin (indole-3-acetic acid, IAA), salinity, drought, high temperatures, and cold, was also investigated here. Five VrYUCCAs responded to IAA in the root, while only VrYUCCA8a was induced in the leaf. VrYUCCA2a, VrYUCCA6a, VrYUCCA8a, VrYUCCA8b, and VrYUCCA10 seemed to dominate under abiotic stresses, due to their sensitivity to the other four treatments. However, the response modes of the VrYUCCAs varied, indicating that they may regulate different stresses in distinct ways to finely adjust IAA content. The comprehensive analysis of the VrYUCCAs in this study lays a solid foundation for further investigation of VrYUCCA genes' mechanisms and applications in breeding.

Genome-wide identification and expression analysis of two-component system genes in sweet potato (Ipomoea batatas L.)

Two-component system (TCS), which comprises histidine kinases (HKs), histidine phosphotransfer proteins (HPs), and response regulators (RRs), plays essential roles in regulating plant growth, development, and response to various environmental stimuli. TCS genes have been comprehensively identified in various plants, while studies on the genome-wide identification and analysis of TCS in sweet potato were still not reported. Therefore, in this study, a total of 90 TCS members consisting of 20 HK(L)s, 11 HPs, and 59 RRs were identified in the genome of Ipomoea batatas. Furthermore, their gene structures, conserved domains, and phylogenetic relationships were analyzed in detail. Additionally, the gene expression profiles in various organs were analyzed, and response patterns to adverse environmental stresses were investigated. The results showed that these 90 TCS genes were mapped on 15 chromosomes with a notably uneven distribution, and the expansion of TCS genes in sweet potato was attributed to both segmental and tandem duplications. The majority of the TCS genes showed distinct organ-specific expression profiles, especially in three types of roots (stem roots, fibrous roots, tuberous roots). Moreover, most of the TCS genes were either induced or suppressed upon treatment with abiotic stresses (drought, salinity, cold, heat) and exogenous phytohormone abscisic acid (ABA). In addition, the yeast-two hybrid system was used to reveal the HK-HP-RR protein-protein interactions. IbHP1, IbHP2, IbHP4, and IbHP5 could interact with three HKs (IbHK1a, IbHK1b, and IbHK5), and also interact with majority of the type-B RRs (IbRR20-IbRR28), while no interaction affinity was detected for IbHP3. Our systematic analyses could provide insights into the characterization of the TCS genes, and further the development of functional studies in sweet potato.

Comprehensive analysis of aberrantly expressed circRNAs, mRNAs and lncRNAs in patients with nasopharyngeal carcinoma

BACKGROUND

The location of nasopharyngeal cancer is hidden, so it is difficult to diagnose at an early stage. In this study, we aimed to investigate the expression profiles of circRNAs, mRNAs and IncRNAs and to provide some basis for further studies.

METHODS

Expression profiles of circRNAs, mRNAs, and lncRNAs were analyzed using microarray techniques. The differentially expressed ncRNA was calculated by bioinformatics.

RESULTS

A total of 3048 circRNAs, 2179 lncRNAs, and 2015 mRNAs were detected to be significantly differentially expressed in NPC. The most upregulated circRNAs, lncRNAs, and mRNAs were hsa-circ-0067562, NONHSAT232922.1, and HOXB13, respectively. And, the most downregulated circRNAs, lncRNAs, and mRNAs were hsa_circ_0078837, lnc-TTC8-4:3, and LTF, respectively. The number of upregulated DE lncRNAs was more than twice than those downregulated. Our data showed that 80.44% of pairs of lncRNAs and cis-mRNAs demonstrated positive correlations. For lncRNAs and trans-mRNAs pairs, 53.7% of pairs showed positive correlation. LncRNA-mediated cis regulation is a prevalent regulatory mode in the development of nasopharyngeal carcinoma. CR1, LRMP and SORBS2 are predicted to be mediated not only by cis-acting lncRNA modes of action, but also by trans-acting lncRNA mechanisms. Additionally, we constructed a diagnostic prediction model with a high sensitivity and specificity.

CONCLUSION

Our study characterized the landscape of circRNAs, mRNAs and lncRNAs in NPC tissue and provided novel insights into the molecular mechanisms of NPC.

Identification and expression profile analysis of chemosensory genes in pine needle gall midge, Thecodiplosis japonensis (Diptera: Cecidomyiidae)

Insects have highly specialized and sensitive olfactory systems involving several chemosensory genes to locate their mates and hosts or escape from predators. Pine needle gall midge, Thecodiplosis japonensis (Diptera: Cecidomyiidae), has invaded China since 2016 and caused serious damage. Till now, there is no environmentally friendly measure to control this gall midge. Screening molecules with high affinity to target odorant-binding protein to develop highly efficient attractants is a potential pest management method. However, the chemosensory genes in T. japonensis are still unclear. We identified 67 chemosensory-related genes in the transcriptomes of antennae, including 26 OBPs, 2 CSPs, 17 ORs, 3 SNMPs, 6 GRs, and 13 IRs, using high throughput sequencing. Phylogenetic analysis of these six chemosensory gene families among Dipteran was performed to classify and predict the functions. The expression profiles of OBPs, CSPs and ORs were validated by quantitative real-time PCR. 16 of the 26 OBPs were biased expressed in antennae. TjapORco and TjapOR5 were highly expressed in the antenna of unmated male and female adults. The functions of related OBPs and ORs genes were also discussed. These results provide a basis for the functional research on chemosensory genes at the molecular level.

Genome-wide analysis of R2R3-MYB transcription factors in Boehmeria nivea (L.) gaudich revealed potential cadmium tolerance and anthocyanin biosynthesis genes

Gene family, especially MYB as one of the largest transcription factor family in plants, the study of its subfunctional characteristics is a key step in the study of plant gene function. The sequencing of ramie genome provides a good opportunity to study the organization and evolutionary characters of the ramie MYB gene at the whole genome level. In this study, a total of 105 BnGR2R3-MYB genes were identified from ramie genome and subsequently grouped into 35 subfamilies according to phylogeny divergence and sequences similarity. Chromosomal localization, gene structure, synteny analysis, gene duplication, promoter analysis, molecular characteristics and subcellular localization were accomplished using several bioinformatics tools. Collinearity analysis showed that the segmental and tandem duplication events is the dominant form of the gene family expansion, and duplications prominent in distal telomeric regions. Highest syntenic relationship was obtained between BnGR2R3-MYB genes and that of Apocynum venetum (88). Furthermore, transcriptomic data and phylogenetic analysis revealed that BnGMYB60, BnGMYB79/80 and BnGMYB70 might inhibit the biosynthesis of anthocyanins, and UPLC-QTOF-MS data further supported the results. qPCR and phylogenetic analysis revealed that the six genes (BnGMYB9, BnGMYB10, BnGMYB12, BnGMYB28, BnGMYB41, and BnGMYB78) were cadmium stress responsive genes. Especially, the expression of BnGMYB10/12/41 in roots, stems and leaves all increased more than 10-fold after cadmium stress, and in addition they may interact with key genes regulating flavonoid biosynthesis. Thus, a potential link between cadmium stress response and flavonoid synthesis was identified through protein interaction network analysis. The study thus provided significant information into MYB regulatory genes in ramie and may serve as a foundation for genetic enhancement and increased productivity.

Identification and Pharmacological Characterization of Two Serotonin Type 7 Receptor Isoforms from Mythimna separata

Serotonin (5-hydroxytryptamine, 5-HT) is an important neuroactive molecule, as neurotransmitters regulate various biological functions in vertebrates and invertebrates by binding and activating specific 5-HT receptors. The pharmacology and tissue distribution of 5-HT receptors have been investigated in several model insects, and these receptors are recognized as potential insecticide targets. However, little is known about the pharmacological characterization of the 5-HT receptors in important agricultural pests. In this study, we investigated the sequence, pharmacology, and tissue distribution of 5-HT7 receptors from oriental armyworm Mythimna separata (Walker) (Lepidoptera: Noctuidae), an important migratory and polyphagous pest species. We found that the 5-HT7 receptor gene encodes two molecularly distinct transcripts, Msep5-HT7L and Msep5-HT7S, by the mechanism of alternative splicing in M. separata. Msep5-HT7S differs from Msep5-HT7L based on the deletion of 95 amino acids within the third intracellular loop. Two Msep5-HT7 receptor isoforms were activated by 5-HT and synthetic agonists α-methylserotonin, 8-hydroxy-DPAT, and 5-methoxytryptamine, resulting in increased intracellular cAMP levels in a dose-dependent manner, although these agonists showed much poorer potency and efficacy than 5-HT. The maximum efficacy of 5-HT compared to the two 5-HT isoforms was equivalent, but 5-HT exhibited 2.63-fold higher potency against the Msep5-HT7S than the Msep5-HT7L receptor. These two isoforms were also blocked by the non-selective antagonist methiothepin and the selective antagonists WAY-100635, ketanserin, SB-258719, and SB-269970. Moreover, two distinct mRNA transcripts were expressed preferentially in the brain and chemosensory organs of M. separata adults, as determined by qPCR assay. This study is the first comprehensive characterization of two splicing isoforms of 5-HT7 receptors in M. separata, and the first to demonstrate that alternative splicing is also the mechanism for producing multiple 5-HT7 isoforms in insects. Pharmacological and gene expression profiles offer important information that could facilitate further exploration of their function in the central nervous system and peripheral chemosensory organs, and may even contribute to the development of new selective pesticides.

Characteristics and Expression Analysis of Invertase Gene Family in Common Wheat (Triticum aestivum L.)

Invertase (INV) irreversibly catalyzes the conversion of sucrose into glucose and fructose, playing important role in plant development and stress tolerance. However, the functions of INV genes in wheat have been less studied. In this study, a total of 126 TaINV genes were identified using a genome-wide search method, which could be classified into five classes (TaCWI-α, TaCWI-β, TaCI-α, TaCI-β, and TaVI) based on phylogenetic relationship. A total of 101 TaINVs were collinear with their ancestors in the synteny analysis, and we speculated that polyploidy events were the main force in the expansion of the TaINV gene family. Compared with TaCI, TaCWI and TaVI are more similar in gene structure and protein properties. Transcriptome sequencing analysis showed that TaINVs expressed in multiple tissues with different expression levels. Among 19 tissue-specific expressed TaINVs, 12 TaINVs showed grain-specific expression pattern and might play an important role in wheat grain development. In addition, qRT-PCR results further confirmed that TaCWI50 and TaVI27 show different expression in grain weight NILs. Our results demonstrated that the high expression of TaCWI50 and TaVI27 may be associated with a larger TGW phenotype. This work provides the foundations for understanding the grain development mechanism.

Comprehensive Analysis of Potential Prognostic Values of ANGPTLs in Colorectal Cancer

Colorectal cancer (CRC) is one of the most common malignant tumors in the world. CRC recurrence and metastasis cause poor prognosis. ANGPTLs (angiopoietin-like proteins) are a family of proteins that are widely involved in metabolic disease and tumorigenesis. The roles of ANGPTLs in CRC are still controversial and deserve further research. In this study, several databases were employed to explore the expression profiles, prognostic values, genetic alterations, potential biological function, and immune infiltration correlation of ANGPTLs in CRC. The expression of ANGPTL4 was significantly positively correlated with the stage of CRC. Therefore, cell and molecular experiments were further performed to explore the roles of ANGPTL4. Our results showed that the transcriptions of ANGPTLs in colon cancer and rectal cancer tissues were lower than those in normal tissues, but the protein expression varied among different ANGPTLs. In addition, the high expression of ANGPTLs led to a relatively poor oncological outcome. Specifically, the expression of ANGPTL4 is significantly positively correlated with the stage of CRC. Further investigation revealed that ANGPTLs are mainly involved in signal transduction and the regulation of transcription, while KEGG pathway analyses demonstrated pathways in cancer. Additionally, we also observed that ANGPTL4 could promote the proliferation and migration of CRC cells, and four specific small molecule compounds had potential ANGPTL4-binding capabilities, suggesting the clinical application of these small molecule compounds on CRC treatment. Our findings imply the prognostic values and potential therapeutic targets of ANGPTLs in CRC.

Effects of abiotic stresses on the expression of chitinase-like genes in Acyrthosiphon pisum

Insect chitinases play a crucial part to digest chitin in the exoskeleton during the molting process. However, research on insect chitinase related to the environmental stress response is very limited. This study was the first conducted to expression analysis of chitinase- related genes in A. pisum under abiotic stresses. Here, we identified five chitinase-like proteins (ApIDGF, ApCht3, ApCht7, ApCht10 and ApENGase), and clustered them into five groups (group II, III, V, Ⅹ, and ENGase). Developmental expression analysis revealed that the five A. pisum chitinase-related genes were expressed at whole developmental stages with different relative expression patterns. When aphids were exposed to various abiotic stresses including temperature, insecticide and the stress 20-hydroxyecdysone (20E), all five chitinase genes were differentially expressed in A. pisum. The results showed that insecticide such as imidacloprid down-regulated the expression of these five Cht-related genes. Analysis of temperature stress of A. pisum chitinase suggested that ApCht7 expression was high at 10°C, which demonstrates its important role in pea aphids under low temperature. Conversely, ApCht10 was more active under high temperature stress, as it was significantly up-regulated at 30°C. Besides, 20E enhanced ApCht3 and ApCht10 expression in A. pisum, but reduced ApCht7 expression. These findings provide basic information and insights for the study of the role of these genes under abiotic stress, which advances our knowledge in the management of pea aphids under multiple stresses.

Molecular and Pharmacological Characterization of β-Adrenergic-like Octopamine Receptors in the Endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae)

Octopamine (OA) is structurally and functionally similar to adrenaline/noradrenaline in vertebrates, and OA modulates diverse physiological and behavioral processes in invertebrates. OA exerts its actions by binding to specific octopamine receptors (OARs). Functional and pharmacological characterization of OARs have been investigated in several insects. However, the literature on OARs is scarce for parasitoids. Here we cloned three β-adrenergic-like OARs (CcOctβRs) from Cotesia chilonis. CcOctβRs share high similarity with their own orthologous receptors. The transcript levels of CcOctβRs were varied in different tissues. When heterologously expressed in CHO-K1 cells, CcOctβRs induced cAMP production, and were dose-dependently activated by OA, TA and putative octopaminergic agonists. Their activities were inhibited by potential antagonists and were most efficiently blocked by epinastine. Our study offers important information about the molecular and pharmacological properties of β-adrenergic-like OARs from C. chilonis that will provide the basis to reveal the contribution of individual receptors to the physiological processes and behaviors in parasitoids.

Comprehensive analysis of the NAC transcription factor gene family in Kandelia obovata reveals potential members related to chilling tolerance

BACKGROUND

Kandelia obovata is an important mangrove species extensively distributed in Eastern Asia that is susceptible to low-temperature stress. NAC (NAM, ATAF1/2 and CUC2) domain proteins are transcription factors (TFs) that play various roles in plant growth and development and in the plant response to environmental stresses. Nevertheless, genome-wide analyses of K. obovata NAC genes (KoNACs) and their responses to chilling stress have rarely been studied.

METHODS

The KoNAC gene family was identified and characterized using bioinformatic analysis, the subcellular location of some NAC proteins was confirmed using confocal microscopy analysis, and the KoNACs that responded to chilling stress were screened using RNA-seq and qRT-PCR analysis.

RESULTS

A total of 79 KoNACs were identified, and they were unequally distributed across all 18 chromosomes of K. obovata. The KoNAC proteins could be divided into 16 subgroups according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana. The KoNACs exhibited greater synteny with A. thaliana sequences than with Oryza sativa sequences, indicating that KoNACs underwent extensive evolution after the divergence of dicotyledons and monocotyledons. Segmental duplication was the main driving force of the expansions of KoNAC genes. Confocal microscopy analysis verified that the four randomly selected KoNACs localized to the nucleus, indicating the accuracy of the bioinformatic predictions. Tissue expression pattern analysis demonstrated that some KoNAC genes showed tissue-specific expression, suggesting that these KoNACs might be important for plant development and growth. Additionally, the expression levels of 19 KoNACs were significantly (15 positively and 4 negatively) induced by cold treatment, demonstrating that these KoNACs might play important roles during cold stress responses and might be candidate genes for the genetic engineering of K. obovata with enhanced chilling stress tolerance. Coexpression network analysis revealed that 381 coexpressed pairs (between 13 KoNACs and 284 other genes) were significantly correlated.

CONCLUSIONS

Seventy-nine KoNACs were identified in K. obovata, nineteen of which displayed chilling-induced expression patterns. These genes may serve as candidates for functional analyses of KoNACs engaged in chilling stress. Our results lay the foundation for evolutionary analyses of KoNACs and their molecular mechanisms in response to environmental stress.