Tissue-specific Expression Patterns Research Articles

Abstract 1388: Spatial transcriptome analysis of EMT trait and corresponding immune composition study in HNSCC

Abstract Head and neck squamous cell carcinoma (HNSCC) is a prevalent malignant cancer worldwide. Conventional treatments for this type of cancer involve surgical resection, radiation therapy, and chemotherapy. However, in cases where the tumor recurs or metastasizes distantly, rendering surgical resection or radiation therapy impractical, the only viable option is the utilization of targeted drugs in conjunction with chemotherapy. Despite such interventions, the prognosis remains grim. Hence, the identification of genetic characteristics or biomarkers associated with easy metastasis or recurrence in the early stages of the tumor has become imperative for effective head and neck cancer treatment. The objective of this study is to analyze tumor specimens from the primary site and lymph node metastasis post-surgical resection in HNSCC patients. We employed the GeoMx Digital Spatial Profiler (DSP) for imaging-based digitized gene expression analysis. The DSP provides gene expression information at specific locations we selected, such as the primary tumor center core, primary tumor invasion borders, and lymph node metastasis specimens. The GeoMx DSP spatial analysis tool illustrates the cancer transcriptome atlas (CTA), encompassing the expression of 1825 genes related to cancer and immunity from the primary site to the metastatic site. We identified tissue-specific EMT gene expression patterns in paired tumor regions of interest (ROIs). ITGB1 emerged as the most significant co-expression marker in the primary outer tumor and metastatic tumor. Distinct EMT signature genes expressed in the HNSCC primary site and lymph node metastatic site were also discovered. Additionally, we conducted CTA profiling in correlated stroma ROIs. Through myeloid activity analysis in HNSCC paired samples, we identified distinct myeloid activation signatures in primary and metastatic site stroma cells. These results will provide crucial insights into head and neck cancer treatment strategies based on different tumor site characteristics. Citation Format: Chih-Hung Chung, Muh-Hwa Yang. Spatial transcriptome analysis of EMT trait and corresponding immune composition study in HNSCC [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1388.

Genome-Wide Identification of PYL/RCAR ABA Receptors and Functional Analysis of LbPYL10 in Heat Tolerance in Goji (Lycium barbarum).

The characterization of the PYL/RCAR ABA receptors in a great deal of plant species has dramatically advanced the study of ABA functions involved in key physiological processes. However, the genes in this family are still unclear in Lycium (Goji) plants, one of the well-known economically, medicinally, and ecologically valuable fruit crops. In the present work, 12 homologs of Arabidopsis PYL/RCAR ABA receptors were first identified and characterized from Lycium (L.) barbarum (LbPYLs). The quantitative real-time PCR (qRT-PCR) analysis showed that these genes had clear tissue-specific expression patterns, and most of them were transcribed in the root with the largest amount. Among the three subfamilies, while the Group I and Group III members were down-regulated by extraneous ABA, the Group II members were up-regulated. At 42 °C, most transcripts showed a rapid and violent up-regulation response to higher temperature, especially members of Group II. One of the genes in the Group II members, LbPYL10, was further functionally validated by virus-induced gene silencing (VIGS) technology. LbPYL10 positively regulates heat stress tolerance in L. barbarum by alleviating chlorophyll degradation, thus maintaining chlorophyll stability. Integrating the endogenous ABA level increase following heat stress, it may be concluded that LbPYL-mediated ABA signaling plays a vital role in the thermotolerance of L. barbarum plants. Our results highlight the strong potential of LbPYL genes in breeding genetically modified L. barbarum crops that acclimate to climate change.

The comprehensive landscape of prognosis, immunity, and function of the GLI family by pan-cancer and single-cell analysis.

The Hedgehog (Hh) signaling pathway has been implicated in the pathogenesis of various cancers. However, the roles of the downstream GLI family (GLI1, GLI2, and GLI3) in tumorigenesis remain elusive. This study aimed to unravel the genetic alterations of GLI1/2/3 in cancer and their association with the immune microenvironment and related signaling pathways. Firstly, we evaluated the expression profiles of GLI1/2/3 in different cancer types, analyzed their prognostic and predictive values, and assessed their correlation with tumor-infiltrating immune cells. Secondly, we explored the relationships between GLI1/2/3 and genetic mutations, epigenetic modifications, and clinically relevant drugs. Finally, we performed enrichment analysis to decipher the underlying mechanisms of GLI1/2/3 in cancer initiation and progression. Our results revealed that the expression levels of GLI1/2/3 were positively correlated in most cancer tissues, suggesting a cooperative role of these factors in tumorigenesis. We also identified tissue-specific expression patterns of GLI1/2/3, which may reflect the distinct functions of these factors in different cell types. Furthermore, GLI1/2/3 expression displayed significant associations with poor prognosis in several cancers, indicating their potential as prognostic biomarkers and therapeutic targets. Importantly, we found that GLI1/2/3 modulated the immune microenvironment by regulating the recruitment, activation, and polarization of cancer-associated fibroblasts, endothelial cells, and macrophages. Additionally, functional enrichment analyses indicated that GLI1/2/3 are involved in the regulation of epithelial-mesenchymal transition (EMT). Together, our findings shed new light on the roles of GLI1/2/3 in tumorigenesis and provide a potential basis for the development of novel therapeutic strategies targeting GLI-mediated signaling pathways in cancer.

Cloning and functional analysis of Gb4CL1 and Gb4CL2 from Ginkgo biloba.

4-Coumarate-CoA ligase (4CL) gene plays vital roles in plant growth and development, especially the regulation of lignin metabolism and flavonoid synthesis. To investigate the potential function of 4CL in the lignin biosynthesis of Ginkgo biloba, this study identified two 4CL genes, Gb4CL1 and Gb4CL2, from G. biloba genome. Based on the phylogenetic tree analysis, Gb4CL1 and Gb4CL2 protein were classified into Class I, which has been confirmed to be involved in lignin biosynthesis. Therefore, it can be inferred that these two genes may also participate in lignin metabolism. The tissue-specific expression patterns of these two genes revealed that Gb4CL1 was highly expressed in microstrobilus, whereas Gb4CL2 was abundant in immature leaves. The onion transient expression assay indicated that Gb4CL1 was predominantly localized in the nucleus, indicating its potential involvement in nuclear functions, while Gb4CL2 was observed in the cell wall, suggesting its role in cell wall-related processes. Phytohormone response analysis revealed that the expression of both genes was upregulated in response to indole acetic acid, while methyl jasmonate suppressed it, gibberellin exhibited opposite effects on these genes. Furthermore, Gb4CL1 and Gb4CL2 expressed in all tissues containing lignin that showed a positive correlation with lignin content. Thus, these findings suggest that Gb4CL1 and Gb4CL2 are likely involved in lignin biosynthesis. Gb4CL1 and Gb4CL2 target proteins were successfully induced in Escherichia coli BL21 with molecular weights of 85.5 and 89.2kDa, proving the integrity of target proteins. Our findings provided a basis for revealing that Gb4CL participated in lignin synthesis in G. biloba.

Genome-wide identification and expression analysis of Autophagy-Related Genes in eggplant (Solanum melongena L.)

Autophagy, a highly conserved cellular process in plants, plays a vital role in regulating abiotic stress responses and senescence. Eggplant is a widely cultivated vegetable crop, its fruit is susceptible to cold damage, and its molecular mechanisms underlying stress tolerance, fruit ripening, and senescence remain largely uncharacterized. In this study, we systematically identified 41 Autophagy-Related Genes (ATGs) in eggplant and through bioinformatics methods, we analyzed their gene structures, evolutionary features, chromosomal locations, and promoter elements. Additionally, we conducted RNA-seq transcriptome analysis of various eggplant organs as well as under different stress conditions, and assessed the expression patterns of the SmATG8 under cold stress in eggplant fruit. Our analysis results show that the SmATGs promoters have diverse cis-regulatory elements, with most SmATGs exhibiting low and stable expression levels in different tissues. Some genes displayed tissue-specific expression patterns. Under 4 °C low-temperature stress, the expression of the SmATG8 family was continuously upregulated, indicating that autophagy plays a crucial role in defending against cold damage in eggplant. Furthermore, we found that ATGs in dicots and monocots exhibit distinct homology patterns. Within the dicots, eggplant and pepper, both being Solanaceae plants, their ATGs exhibit high homology, suggesting that SmATGs in Solanaceae plants have unique homologous relationships. The analysis of SmATGs unveiled their potential contributions to plant stress responses, fruit ripening, and senescence, thus offering a theoretical foundation for further exploration of stress resistance and senescence mechanisms in eggplant.

Identification and Characterization of the AREB/ABF Gene Family in Three Orchid Species and Functional Analysis of DcaABI5 in Arabidopsis

AREB/ABF (ABA response element binding) proteins in plants are essential for stress responses, while our understanding of AREB/ABFs from orchid species, important traditional medicinal and ornamental plants, is limited. Here, twelve AREB/ABF genes were identified within three orchids’ complete genomes and classified into three groups through phylogenetic analysis, which was further supported with a combined analysis of their conserved motifs and gene structures. The cis-element analysis revealed that hormone response elements as well as light and stress response elements were widely rich in the AREB/ABFs. A prediction analysis of the orchid ABRE/ABF-mediated regulatory network was further constructed through cis-regulatory element (CRE) analysis of their promoter regions. And it revealed that several dominant transcriptional factor (TF) gene families were abundant as potential regulators of these orchid AREB/ABFs. Expression profile analysis using public transcriptomic data suggested that most AREB/ABF genes have distinct tissue-specific expression patterns in orchid plants. Additionally, DcaABI5 as a homolog of ABA INSENSITIVE 5 (ABI5) from Arabidopsis was selected for further analysis. The results showed that transgenic Arabidopsis overexpressing DcaABI5 could rescue the ABA-insensitive phenotype in the mutant abi5. Collectively, these findings will provide valuable information on AREB/ABF genes in orchids.

CHARACTERIZATION OF GALACTINOL SYNTHASE (GOLS) GENES IN IPOMOEA TRILOBA AND I. TRIFIDA: AN APPROACH USING BIOINFORMATICS

Galactinol synthase (GolS - EC 2.4.1.123) is classified as a key enzyme that catalyzes the first step in the synthesis pathway of the raffinose family (RFOs). Although GolS genes have been characterized in several important species, their characterization in sweet potatoes has yet to be explored. Ipomoea trifida (Kunth) G. Don (2n = 2x = 30) is currently described as one of the closest ancestors of the sweet potato and considered an excellent crossbreeding species, allowing the introgression of important genes such as GolS. This study aimed to identify and characterize in silico the GolS genes in I. trifida and compare with I. triloba. We identified nine GolS genes, five in I. triloba and four in I. trifida. Our study encompassed various aspects, including gene structure analysis, motif identification, chromosomal distribution, synteny analysis, and gene expression. The presence of gene duplications and purifying selection were highlighted, suggesting the evolutionary significance of GolS genes in these species. Phylogenetic analysis categorized GolS proteins into three groups, potentially reflecting distinct functional roles. Furthermore, synteny analysis revealed orthologous relationships between GolS genes in the studied species and related plants, contributing to our understanding of their evolutionary history. In silico expression analysis across diverse tissues unveiled tissue-specific expression patterns, hinting at specialized roles for GolS genes in different plant organs. These findings contribute to the broader field of plant genetics, carbohydrate metabolism, and agriculture, offering opportunities for crop improvement and sustainable food production.

Genome-wide identification, expression profiling, and protein interaction analysis of the CCoAOMT gene family in the tea plant (Camellia sinensis)

BackgroundThe caffeoyl-CoA-O methyltransferase (CCoAOMT) family plays a crucial role in the oxidative methylation of phenolic substances and is involved in various plant processes, including growth, development, and stress response. However, there is a limited understanding of the interactions among CCoAOMT protein members in tea plants.ResultsIn this study, we identified 10 members of the CsCCoAOMT family in the genome of Camellia sinensis (cultivar ‘HuangDan’), characterized by conserved gene structures and motifs. These CsCCoAOMT members were located on six different chromosomes (1, 2, 3, 4, 6, and 14). Based on phylogenetic analysis, CsCCoAOMT can be divided into two groups: I and II. Notably, the CsCCoAOMT members of group Ia are likely to be candidate genes involved in lignin biosynthesis. Moreover, through the yeast two-hybrid (Y2H) assay, we established protein interaction networks for the CsCCoAOMT family, revealing 9 pairs of members with interaction relationships.ConclusionsWe identified the CCoAOMT gene family in Camellia sinensis and conducted a comprehensive analysis of their classifications, phylogenetic and synteny relationships, gene structures, protein interactions, tissue-specific expression patterns, and responses to various stresses. Our findings shed light on the evolution and composition of CsCCoAOMT. Notably, the observed interaction among CCoAOMT proteins suggests the potential formation of the O-methyltransferase (OMT) complex during the methylation modification process, expanding our understanding of the functional roles of this gene family in diverse biological processes.

Genome-wide identification and comparative analysis of the Amino Acid Transporter (AAT) gene family and their roles during Phaseolus vulgaris symbioses

Amino acid transporters (AATs) are essential integral membrane proteins that serve multiple roles, such as facilitating the transport of amino acids across cell membranes. They play a crucial role in the growth and development of plants. Phaseolus vulgaris, a significant legume crop, serves as a valuable model for studying root symbiosis. In this study, we have conducted an exploration of the AAT gene family in P. vulgaris. In this research, we identified 84 AAT genes within the P. vulgaris genome sequence and categorized them into 12 subfamilies based on their similarity and phylogenetic relationships with AATs found in Arabidopsis and rice. Interestingly, these AAT genes were not evenly distributed across the chromosomes of P. vulgaris . Instead, there was an unusual concentration of these genes located toward the outer edges of chromosomal arms. Upon conducting motif analysis and gene structural analysis, we observed a consistent presence of similar motifs and an intron-exon distribution pattern among the subfamilies. When we analyzed the expression profiles of PvAAT genes, we noted tissue-specific expression patterns. Furthermore, our investigation into AAT gene expression under rhizobial and mycorrhizal symbiotic conditions revealed that certain genes exhibited high levels of expression. Specifically, ATLa5 and LHT2 was notably upregulated under both symbiotic conditions. These findings point towards a potential role of AATs in the context of rhizobial and mycorrhizal symbiosis in P. vulgaris, in addition to their well-established regulatory functions.

Molecular cloning and characterization of farnesyl diphosphate synthase from Rosa rugosa Thunb associated with salinity stress.

Rosa rugosa, a renowned ornamental plant, is cultivated for its essential oil containing valuable monoterpenes, sesquiterpenes, and other compounds widely used in the floriculture industry. Farnesyl diphosphate synthase (FPPS) is a key enzyme involved in the biosynthesis of sesquiterpenes and triterpenes for abiotic or biotic stress. In this study, we successfully cloned and characterized a full-length FPPS- encoding cDNA identified as RrFPPS1 using RT-PCR from R. rugosa. Phylogenetic analysis showed that RrFPPS1 belonged to the angiosperm-FPPS clade. Transcriptomic and RT-qPCR analyses revealed that the RrFPPS1 gene had tissue-specific expression patterns. Subcellular localization analysis using Nicotiana benthamiana leaves showed that RrFPPS1 was a cytoplasmic protein. In vitro enzymatic assays combined with GC-MS analysis showed that RrFPPS1 produced farnesyl diphosphate (FPP) using isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) as substrates to provide a precursor for sesquiterpene and triterpene biosynthesis in the plant. Additionally, our research found that RrFPPS1 was upregulated under salt treatment. These substantial findings contribute to an improved understanding of terpene biosynthesis in R. rugosa and open new opportunities for advancements in horticultural practices and fragrance industries by overexpression of the RrFPPS1 gene in vivo increased FPP production and subsequently led to elevated sesquiterpene yields in the future. The knowledge gained from this study can potentially lead to the development of enhanced varieties of R. rugosa with improved aroma, medicinal properties, and resilience to environmental stressors.

Long Intergenic Non-Coding RNAs of Human Chromosome 18: Focus on Cancers.

Malignant neoplasms are characterized by high molecular heterogeneity due to multilevel deregulation of gene expression and cellular functions. It is known that non-coding RNAs, including long intergenic non-coding RNAs (lincRNAs), can play significant roles in cancer biology. The current review focuses on a systematical analysis of genomic, transcriptomic, epigenomic, interactomic, and literature data on 65 lincRNAs of human chromosome 18 in the context of pan-cancer studies. The entire group of lincRNAs can be conditionally divided into 4 subgroups depending on experimental evidence on direct or indirect involvement in cancers and the biological associations with cancers, which we found during the data-mining process: the most studied (5 lincRNAs), moderately or poorly studied (11 lincRNAs), and understudied (31 lincRNAs). For the remaining 18 lincRNAs, data for analysis were fragmentary or missing. Among the key findings were the following: Of the lincRNAs of human chromosome 18, 40% have tissue-specific expression patterns, 22% of lincRNAs are known to have gene fusions, 40% of lincRNAs are prone to gene amplifications and/or deletions in cancers at a frequency greater than 3%, and 23% of lincRNAs are differentially expressed across cancer types, whereas 7% have subtype-specific expression patterns. LincRNAs' interactomes consist of 'master' microRNAs and 47 proteins (including cancer-associated proteins and microRNAs) that can interact with 3 or more lincRNAs. Functional enrichment analysis of a set of highly co-expressed genes retrieved for 17 lincRNAs in different cancer types indicated the potential associations of these lincRNAs with cellular signaling pathways. Six lincRNAs encoded small open-reading frame (smORF) proteins with emerging roles in cancers, and microRNAs as well as proteins with known functions in molecular carcinogenesis can bind to coding regions of smORFs. We identified seven transcriptomic signatures with potential prognostic value, consisting of two to seven different lincRNAs only. Taken together, the literature, biomedical, and molecular biology data analyzed indicated that only five of all lincRNAs of human chromosome 18 are cancer-associated, while eleven other lincRNAs have the tendency to be associated with cancers.

Genome-Wide Analysis of Aux/IAA Gene Family in Artemisia argyi: Identification, Phylogenetic Analysis, and Determination of Response to Various Phytohormones.

Artemisia argyi is a traditional herbal medicine plant, and its folium artemisia argyi is widely in demand due to moxibustion applications globally. The Auxin/indole-3-acetic acid (Aux/IAA, or IAA) gene family has critical roles in the primary auxin-response process, with extensive involvement in plant development and stresses, controlling various essential traits of plants. However, the systematic investigation of the Aux/IAA gene family in A. argyi remains limited. In this study, a total of 61 Aux/IAA genes were comprehensively identified and characterized. Gene structural analysis indicated that 46 Aux/IAA proteins contain the four typical domains, and 15 Aux/IAA proteins belong to non-canonical IAA proteins. Collinear prediction and phylogenetic relationship analyses suggested that Aux/IAA proteins were grouped into 13 distinct categories, and most Aux/IAA genes might experience gene loss during the tandem duplication process. Promoter cis-element investigation indicated that Aux/IAA promoters contain a variety of plant hormone response and stress response cis-elements. Protein interaction prediction analysis demonstrated that AaIAA26/29/7/34 proteins are possibly core members of the Aux/IAA family interaction. Expression analysis in roots and leaves via RNA-seq data indicated that the expression of some AaIAAs exhibited tissue-specific expression patterns, and some AaIAAs were involved in the regulation of salt and saline-alkali stresses. In addition, RT-qPCR results indicated that AaIAA genes have differential responses to auxin, with complex response patterns in response to other hormones, indicating that Aux/IAA may play a role in connecting auxin and other hormone signaling pathways. Overall, these findings shed more light on AaIAA genes and offer critical foundational knowledge toward the elucidation of their function during plant growth, stress response, and hormone networking of Aux/IAA family genes in A. argyi.

Gene co-expression network analysis reveal core responsive genes in Parascaris univalens tissues following ivermectin exposure.

Anthelmintic resistance in equine parasite Parascaris univalens, compromises ivermectin (IVM) effectiveness and necessitates an in-depth understanding of its resistance mechanisms. Most research, primarily focused on holistic gene expression analyses, may overlook vital tissue-specific responses and often limit the scope of novel genes. This study leveraged gene co-expression network analysis to elucidate tissue-specific transcriptional responses and to identify core genes implicated in the IVM response in P. univalens. Adult worms (n = 28) were exposed to 10-11 M and 10-9 M IVM in vitro for 24 hours. RNA-sequencing examined transcriptional changes in the anterior end and intestine. Differential expression analysis revealed pronounced tissue differences, with the intestine exhibiting substantially more IVM-induced transcriptional activity. Gene co-expression network analysis identified seven modules significantly associated with the response to IVM. Within these, 219 core genes were detected, largely expressed in the intestinal tissue and spanning diverse biological processes with unspecific patterns. After 10-11 M IVM, intestinal tissue core genes showed transcriptional suppression, cell cycle inhibition, and ribosomal alterations. Interestingly, genes PgR028_g047 (sorb-1), PgB01_g200 (gmap-1) and PgR046_g017 (col-37 & col-102) switched from downregulation at 10-11 M to upregulation at 10-9 M IVM. The 10-9 M concentration induced expression of cuticle and membrane integrity core genes in the intestinal tissue. No clear core gene patterns were visible in the anterior end after 10-11 M IVM. However, after 10-9 M IVM, the anterior end mostly displayed downregulation, indicating disrupted transcriptional regulation. One interesting finding was the non-modular calcium-signaling gene, PgR047_g066 (gegf-1), which uniquely connected 71 genes across four modules. These genes were enriched for transmembrane signaling activity, suggesting that PgR047_g066 (gegf-1) could have a key signaling role. By unveiling tissue-specific expression patterns and highlighting biological processes through unbiased core gene detection, this study reveals intricate IVM responses in P. univalens. These findings suggest alternative drug uptake of IVM and can guide functional validations to further IVM resistance mechanism understanding.

In Silico Identification and Characterization of Rare Cold Inducible 2 (RCI2) Gene Family in Cotton.

RCI2/PMP3s are involved in biotic and abiotic stresses and have an influence on the regulation of many genes. RCI2/PMP3 genes, which particularly encode small membrane proteins of the PMP3 family, are involved in abiotic stress responses in plants. In this work, in silico studies were used to investigate RCI2's potential function in stress tolerance and organogenesis. We conducted an extensive study of the RCI2 gene family and revealed 36 RCI2 genes from cotton species that were distributed over 36 chromosomes of the cotton genome. Functional and phylogenetic examination of the RCI2/PMP3 gene family has been studied in Arabidopsis, but in cotton, the RCI2/PMP3 genes have not yet been studied. Phylogenetic and sequencing studies revealed that cotton RCI2s are conserved, with most of them categorized into six distinct clades. A chromosome distribution and localization study indicated that cotton RCI2 genes were distributed unevenly on 36 chromosomes with segmental duplications, suggesting that the cotton RCI2 family is evolutionarily conserved. Many cis-elements related to stress responsiveness, development, and hormone responsiveness were detected in the promoter regions of the cotton RCI2. Moreover, the 36 cotton RCI2s revealed tissue-specific expression patterns in the development of cotton performed by transcriptome analysis. Gene structure analysis indicated that nearly all RCI2 genes have two exons and one intron. All of the cotton RCI2 genes were highly sensitive to drought, abscisic acid, salt, and cold treatments, demonstrating that they may be employed as genetic objects to produce stress-resistant plants.

Identification of RSK substrates using an analog-sensitive kinase approach

The p90 ribosomal S6 kinases (RSK) family of serine/threonine kinases comprises four isoforms (RSK1-4) that lie downstream of the ERK1/2 mitogen-activated protein kinase (MAPK) pathway. RSKs are implicated in fine tuning of cellular processes such as translation, transcription, proliferation, and motility. Previous work showed that pathogens such as Cardioviruses could hijack any of the four RSK isoforms to inhibit PKR activation or to disrupt cellular nucleocytoplasmic trafficking. In contrast, some reports suggest non-redundant functions for distinct RSK isoforms, whereas Coffin-Lowry syndrome has only been associated with mutations in the gene encoding RSK2. In this work, we used the analog-sensitive kinase strategy to ask whether the cellular substrates of distinct RSK isoforms differ. We compared the substrates of two of the most distant RSK isoforms: RSK1 and RSK4. We identified a series of potential substrates for both RSKs in cells, and validated RanBP3, PDCD4, IRS2 and ZC3H11A as substrates of both RSK1 and RSK4, and SORBS2 as an RSK1 substrate. In addition, using mutagenesis and inhibitors, we confirmed analog-sensitive kinase data showing that endogenous RSKs phosphorylate TRIM33 at S1119. Our data thus identify a series of potential RSK substrates and suggest that the substrates of RSK1 and RSK4 largely overlap and that the specificity of the various RSK isoforms likely depends on their cell- or tissue-specific expression pattern.

Molecular characterization and expression pattern of Rubisco activase gene GhRCAβ2 in upland cotton (Gossypium hirsutum L.).

Rubisco activase (RCA) is a pivotal enzyme that can catalyse the activation of Rubisco in carbon assimilation pathway. Many studies have shown that RCA may be a potential target for genetic manipulation aimed at enhancing photosynthetic efficiency and crop yield. To understand the biological function of the GhRCAβ2 gene in upland cotton, we cloned the coding sequence (CDS) of the GhRCAβ2 gene and investigated its sequence features, evolutionary relationship, subcellular localization, promoter sequence and expression pattern. The bioinformatics tools were used to analyze the sequence features of GhRCAβ2 protein. Transient transformation of Arabidopsis mesophyll protoplasts was performed to determine the subcellular localization of the GhRCAβ2 protein. The expression pattern of the GhRCAβ2 gene was examined by analyzing transcriptome data and using the quantitative real-time PCR (qRT-PCR). The full-length CDS of GhRCAβ2 was 1317bp, and it encoded a protein with a chloroplast transit peptide. The GhRCAβ2 had two conserved ATP-binding domains, and did not have the C-terminal extension (CTE) domain that was unique to the RCA α-isoform in plants. Evolutionarily, GhRCAβ2 was clustered in Group A, and had a close evolutionary relationship with the soybean RCA. Western blot analysis demonstrated that GhRCAβ2 was immunoreactive to the RCA antibody displaying a molecular weight similar to that of the RCA β-isoform. The GhRCAβ2 protein was found in chloroplast, aligning with its role as a vital enzyme in the process of photosynthesis. The GhRCAβ2 gene had a leaf tissue-specific expression pattern, and the yellow-green leaf mutant exhibited a decreased expression of GhRCAβ2 in comparison to the wild-type cotton plants. The GhRCAβ2 promoter contained several cis-acting elements that respond to light, phytohormones and stress, suggesting that the expression of GhRCAβ2 may be regulated by these factors. An additional examination of stress response indicated that GhRCAβ2 expression was influenced by cold, heat, salt, and drought stress. Notably, diverse expression pattern was observed across different stress conditions. Additionally, low phosphorus and low potassium stress may result in a notable reduction in the expression of GhRCAβ2 gene. Our findings will establish a basis for further understanding the function of the GhRCAβ2 gene, as well as providing valuable genetic knowledge to improve cotton photosynthetic efficiency and yield under challenging environmental circumstances.

Genome-wide identification and analysis of WD40 proteins reveal that NtTTG1 enhances drought tolerance in tobacco (Nicotiana tabacum)

BackgroundWD40 proteins, which are highly prevalent in eukaryotes, play important roles in plant development and stress responses. However, systematic identification and exploration of WD40 proteins in tobacco have not yet been conducted.ResultsIn this study, a total of 399 WD40 regulatory genes were identified in common tobacco (Nicotiana tabacum). Gene structure and motif analysis revealed structural and functional diversity among different clades of tobacco WD40 regulatory genes. The expansion of tobacco WD40 regulatory genes was mainly driven by segmental duplication and purifying selection. A potential regulatory network of NtWD40s suggested that NtWD40s might be regulated by miRNAs and transcription factors in various biological processes. Expression pattern analysis via transcriptome analysis and qRT-PCR revealed that many NtWD40s exhibited tissue-specific expression patterns and might be involved in various biotic and abiotic stresses. Furthermore, we have validated the critical role of NtTTG1, which was located in the nuclei of trichome cells, in enhancing the drought tolerance of tobacco plants.ConclusionsOur study provides comprehensive information to better understand the evolution of WD40 regulatory genes and their roles in different stress responses in tobacco.

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

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

Transcriptome analysis reveals the molecular mechanisms of rubber biosynthesis and laticifer differentiation during rubber seed germination.

The molecular mechanisms underlying the initiation of natural rubber synthesis and laticifer differentiation have not been fully elucidated. In this study, we conducted a time-series transcriptome analysis of five rubber tree tissues at four stages of seed germination. A total of 161,199 DEGs were identified between the two groups, including most 16,673 DEGs (A3 vs B3 and A3 vs C3) and lest 1,210 DEGs (C2 vs D2). We found that the maturation of the seed is accompanied by the formation of laticifer cells in cotyledon. Meanwhile, the analysis of hormones related genes expression may provide effective clues for us to promote the differentiation of laticifer cells in seeds by hormones in the future. In this study, hormone-related gene enrichment analyses revealed that IAA, GA, and CTK were activated in laticifer containing tissues. Similarly, GO and GEGG analysis showed that hormone pathways, especially the auxin pathway, are enriched. Gene expression clustering was analyzed using the short time-series expression miner (STEM), and the analysis revealed four distinct trends in the gene expression profiles. Moreover, we enriched transcription factor (TF) enrichment in cotyledon and embryonic axis tissues, and the MYB type exhibited the most significant difference. Furthermore, our findings revealed that genes related to rubber synthesis exhibited tissue-specific expression patterns during seed germination. Notably, key genes associated with rubber biosynthesis, specifically small rubber particle protein (SRPP) and cis-prenyltransferase (CPT), exhibited significant changes in expression in cotyledon and embryonic axis tissues, suggesting synchronous rubber synthesis with seed germination. Our staining results reveled that laticifer cells were exits in the cotyledon before seed imbibition stage. In conclusion, these results lay the foundation for exploring the molecular mechanisms underlying laticifer differentiation and rubber synthesis during seed germination, deepening our understanding of the initiation stages of rubber biosynthesis and laticifer differentiation.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster.

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.