Identification of CEP peptides encoded by the tobacco (Nicotiana tabacum) genome and characterization of their roles in osmotic and salt stress responses

Members of the CEP (C-terminally Encoded Peptide) gene family have been shown to be involved in various developmental processes and stress responses in plants. In order to understand the roles of CEP peptides in stress response, a comprehensive bioinformatics approach was employed to identify NtCEP genes in tobacco (Nicotiana tabacum L.) and to analyze their potential roles in stress responses. Totally 21 NtCEP proteins were identified and categorized into two subgroups based on their CEP domains. Expression changes of the NtCEP genes in response to various abiotic stresses were analyzed via qRT-PCR and the results showed that a number of NtCEPs were significant up-regulated under drought, salinity, or temperature stress conditions. Furthermore, application of synthesized peptides derived from NtCEP5, NtCEP13, NtCEP14, and NtCEP17 enhanced plant tolerance to different salt stress treatments. NtCEP5, NtCEP9 and NtCEP14, and NtCEP17 peptides were able to promote osmotic tolerance of tobacco plants. The results from this study suggest that NtCEP peptides may serve as important signaling molecules in tobacco's response to abiotic stresses.

Clearing of Vascular Tissue in Arabidopsis thaliana for Reporter Analysis of Gene Expression

To study the gene regulatory mechanisms modulating development is essential to visualize gene expression patterns at cellular resolution. However, this kind of analysis has been limited as a consequence of the plant tissues' opacity. In the last years, ClearSee has been increasingly used to obtain high-quality imaging of plant tissue anatomy combined with the visualization of gene expression patterns. ClearSee is established as a major tissue clearing technique due to its simplicity and versatility.In this chapter, we outline an easy-to-follow ClearSee protocol to analyze gene expression of reporters using either β-glucuronidase (GUS) or fluorescent protein (FP) tags, compatible with different dyes to stain cell walls. We detail materials, equipment, solutions, and procedures to easily implement ClearSee for the study of vascular development in Arabidopsis thaliana, but the protocol can be easily adapted to a variety of plant tissues in a wide range of plant species.

Genome-wide identification and analysis of MIKC-type MADS-box genes expression in Chimonanthus salicifolius

BACKGROUND

MIKC type MADS-box transcription factors are one of the largest gene families and play a pivotal role in flowering time and flower development. Chimonanthus salicifolius belongs to the family Calycanthaceae and has a unique flowering time and flowering morphology compared to other Chimonanthus species, but the research on MIKC type MADS-box gene family of C. salicifolius has not been reported.

OBJECTIVE

Identification, comprehensive bioinformatic analysis, the expression pattern of MIKC-type MADS-box gene family from different tissues of C. salicifolius.

METHODS

Genome-wide investigation and expression pattern under different tissues of the MIKC-type MADS-box gene family in C. salicifolius, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, promoter cis-acting element were performed.

RESULTS

A total of 29 MIKC-type MADS-box genes were identified from the whole genome sequencing. Interspecies synteny analysis revealed more significant collinearity between C. salicifolius and the magnoliids species compared to eudicots and monocots. MIKC-type MADS-box genes from the same subfamily share similar distribution patterns, gene structure, and expression patterns. Compared with Arabidopsis thaliana, Nymphaea colorata, and Chimonanthus praecox, the FLC genes were absent in C. salicifolius, while the AGL6 subfamily was expanded in C. salicifolius. The selectively expanded promoter (AGL6) and lack of repressor (FLC) genes may explain the earlier flowering in C. salicifolius. The loss of the AP3 homologous gene in C. salicifolius is probably the primary cause of the morphological distinction between C. salicifolius and C. praecox. The csAGL6a gene is specifically expressed in the flowering process and indicates the potential function of promoting flowering.

CONCLUSION

This study offers a genome-wide identification and expression profiling of the MIKC-types MADS-box genes in the C. salicifolius, and establishes the foundation for screening flowering development genes and understanding the potential function of the MIKC-types MADS-box genes in the C. salicifolius.

Molecular characterization and expression analysis of immunoglobulins (IgM and IgT) heavy chains in black rockfish (Sebastes schlegelii) that response to bacterial challenge

Sebastes schlegelii is a kind of fish with great economic values. Recently, with the continuous expansion of aquaculture scale and the continuous improvement of aquaculture density, outbreak of various diseases has caused huge economic losses to its aquaculture industry. Study of fish immune system can help to understand the mechanism of immune response to external pathogens and can promote the development of immune prevention and control methods. Immunoglobulins (Igs) are complex glycoproteins that appear to be unique to the vertebrates that can recognize a wide variety of pathogens and recruit immune cells and molecules to destroy pathogens, which are generated by a series of rearrangement and somatic mutations. We therefore studied the immunoglobulin genes of S. schlegelii in view of their important roles in resisting to external pathogen infections. In this study, the immunoglobulin heavy chain genes (sIgM, mIgM, sIgT, and mIgT) of S. schlegelii were successfully identified and cloned. Phylogenetic analysis showed that the IgM and IgT genes of S. schlegelii were clustered together with homologous genes of other species, indicating that they were highly conserved during the evolutionary process. Collinearity analysis showed that the immunoglobulin genes and their adjacent genes were aligned with zebrafish, Atlantic salmon and tilapia, which further confirmed the conserved immunoglobulin gene of teleost. Expression analysis of healthy tissues showed that the expression levels of sIgM, sIgT and mIgT were the highest in the skin, while mIgM was the highest in spleen. After different bacterial infection, IgM and IgT were significantly expressed in skin and gill, which may be because skin and gill are the first line of defense against the infection pathogens. Subcellular localization showed that the mIgT protein was expressed in both the cell membrane and cytoplasm. Meanwhile, recombinant protein of mIgT was obtained in vitro, which laid a foundation for subsequent protein function studies. These results provide a theoretical basis for understanding the immunity role of immunoglobulin in S. schlegelii.

"In the light of evolution:" keratins as exceptional tumor biomarkers

Keratins (KRTs) are the intermediate filament-forming proteins of epithelial cells, classified, according to their physicochemical properties, into "soft" and "hard" keratins. They have a key role in several aspects of cancer pathophysiology, including cancer cell invasion and metastasis, and several members of the KRT family serve as diagnostic or prognostic markers. The human genome contains both, functional KRT genes and non-functional KRT pseudogenes, arranged in two uninterrupted clusters on chromosomes 12 and 17. This characteristic renders KRTs ideal for evolutionary studies. Herein, comprehensive phylogenetic analyses of KRT homologous proteins in the genomes of major taxonomic divisions were performed, so as to fill a gap in knowledge regarding the functional implications of keratins in cancer biology among tumor-bearing species. The differential expression profiles of KRTs in diverse types of cancers were investigated by analyzing high-throughput data, as well. Several KRT genes, including the phylogenetically conserved ones, were found to be deregulated across several cancer types and to participate in a common protein-protein interaction network. This indicates that, at least in cancer-bearing species, these genes might have been under similar evolutionary pressure, perhaps to support the same important function(s). In addition, semantic relations between KRTs and cancer were detected through extensive text mining. Therefore, by applying an integrative in silico pipeline, the evolutionary history of KRTs was reconstructed in the context of cancer, and the potential of using non-mammalian species as model organisms in functional studies on human cancer-associated KRT genes was uncovered.

Plant growth promoting rhizobacteria modulates the antioxidant defense and the expression of stress-responsive genes providing Pb accumulation and tolerance of grass pea

To ensure the success of phytoremediation, it is important to consider the appropriate combination of plants and microorganisms. This study was conducted to get a better insight into the underlying molecular and biochemical mechanism of grass pea (Lathyrus sativus L.) induced by plant growth promoting rhizobacteria (PGPR), when exposed for 3, 6, 9, and 14 days to 1 mM Pb in a hydroponic system. The significant positive effect of bacterial inoculation was reproduced in various parameters. Results indicated that inoculation of PGPR significantly increased the accumulation of Pb by 20%, 66%, 43%, and 36% in roots and by 46%, 55%, 37%, and 46% in shoots, respectively after 3, 6, 9, and 14 days of metal exposure compared to the uninoculated plants. The metal accumulation in grass pea plants triggered a significant elevation in the synthesis of non-protein thiols (NPT), particularly in inoculated plant leaves where it was about 3 and 2-fold higher than the uninoculated set on the 6th and the 9th day. Nevertheless, Pb treatment significantly increased oxidative stress and membrane damage in leaves with the highest hydrogen peroxide (H2O2) production and tissue malondialdehyde (MDA) concentration recorded in uninoculated plants. Furthermore, the PGPR inoculation alleviated the oxidative stress, improved significantly plant tolerance, and modulated the activities of antioxidant enzymes (SOD, CAT, APX, GR, DHAR, and MDHAR). Similarly, the expression patterns of LsPCS, LsGCN, LsCNGC, LsGR, and LsGST through qRT-PCR demonstrated that bacterial inoculation significantly induced gene expression levels in leaves 6 days after Pb treatment, indicating that PGPR act as regulators of stress-responsive genes. The findings suggest the key role of PGPR (R. leguminosarum (M5) + Pseudomonas fluorescens (K23) + Luteibacter sp. + Variovorax sp.) in enhancing Pb accumulation, reducing metal toxicity, strengthening of the antioxidant system, and conferring higher Pb tolerance to grass pea plants. Hence, the association Lathyrus sativus-PGPR is an effective tool to achieve the goal of remediation of Pb contaminated sites.

Comparison of meiotic transcriptomes of three maize inbreds with different origins reveals differences in cell cycle and recombination

Background

Cellular events during meiosis can differ between inbred lines in maize. Substantial differences in the average numbers of chiasmata and double-strand breaks (DSBs) per meiotic cell have been documented among diverse inbred lines of maize: CML228, a tropical maize inbred line, B73 and Mo17, temperate maize lines. To determine if gene expression might explain these observed differences, an RNA-Seq experiment was performed on CML228 male meiocytes which was compared to B73 and Mo17 male meiocytes, where plants were grown in the same controlled environment.

Results

We found that a few DSB-repair/meiotic genes which promote class I crossovers (COs) and the Zyp1 gene which limits newly formed class I COs were up-regulated, whereas Mus81 homolog 2 which promotes class II COs was down-regulated in CML228. Although we did not find enriched gene ontology (GO) categories directly related to meiosis, we found that GO categories in membrane, localization, proteolysis, energy processes were up-regulated in CML228, while chromatin remodeling, epigenetic regulation, and cell cycle related processes including meiosis related cell cycle processes were down-regulated in CML228. The degree of similarity in expression patterns between the three maize lines reflect their genetic relatedness: B73 and Mo17 had similar meiotic expressions and CML228 had a more distinct expression profile.

Conclusions

We found that meiotic related genes were mostly conserved among the three maize inbreds except for a few DSB-repair/meiotic genes. The findings that the molecular players in limiting class I CO formation (once CO assurance is achieved) were up-regulated and those involved in promoting class II CO formation were down-regulated in CML228 agree with the lower chiasmata number observed in CML228 previously. In addition, epigenetics such as chromatin remodeling and histone modification might play a role. Transport and energy-related processes was up-regulated and Cyclin13 was down-regulated in CML228. The direction of gene expression of these processes agree with that previously found in meiotic tissues compared with vegetative tissues. In summary, we used different natural maize inbred lines from different climatic conditions and have shown their differences in expression landscape in male meiocytes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08922-w.

Genome-wide identification and expression analysis of the ftsH protein family and its response to abiotic stress in Nicotiana tabacum L

BACKGROUND

The filamentous temperature-sensitive H protease (ftsH) gene family plays an important role in plant growth and development. FtsH proteins belong to the AAA protease family. Studies have shown that it is a key gene for plant chloroplast development and photosynthesis regulation. In addition, the ftsH gene is also involved in plant response to stress. At present, the research and analysis of the ftsH gene family are conducted in microorganisms such as Escherichia coli and Oenococcus and various plants such as Arabidopsis, pear, rice, and corn. However, analysis reports on ftsH genes from tobacco (Nicotiana tabacum L.), an important model plant, are still lacking. Since ftsH genes regulate plant growth and development, it has become necessary to systematically study this gene in an economically important plant like tobacco.

RESULTS

This is the first study to analyze the ftsH gene from Nicotiana tabacum L. K326 (NtftsH). We identified 20 ftsH genes from the whole genome sequence, renamed them according to their chromosomal locations, and divided them into eight subfamilies. These 20 NtftsH genes were unevenly distributed across the 24 chromosomes. We found four pairs of fragment duplications. We further investigated the collinearity between these genes and related genes in five other species. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis identified differential expression patterns of NtftsH in different tissues and under various abiotic stress conditions.

CONCLUSIONS

This study provides a comprehensive analysis of the NtftsH gene family. The exon-intron structure and motif composition are highly similar in NtftsH genes that belong to the same evolutionary tree branch. Homology analysis and phylogenetic comparison of ftsH genes from several different plants provide valuable clues for studying the evolutionary characteristics of NtftsH genes. The NtftsH genes play important roles in plant growth and development, revealed by their expression levels in different tissues as well as under different stress conditions. Gene expression and phylogenetic analyses will provide the basis for the functional analysis of NtftsH genes. These results provide a valuable resource for a better understanding of the biological role of the ftsH genes in the tobacco plant.

Genome-wide identification and expression analysis of the bZIP transcription factor family genes in response to abiotic stress in Nicotiana tabacum L

BACKGROUND

The basic leucine zipper (bZIP) transcription factor (TF) is one of the largest families of transcription factors (TFs). It is widely distributed and highly conserved in animals, plants, and microorganisms. Previous studies have shown that the bZIP TF family is involved in plant growth, development, and stress responses. The bZIP family has been studied in many plants; however, there is little research on the bZIP gene family in tobacco.

RESULTS

In this study, 77 bZIPs were identified in tobacco and named NtbZIP01 through to NtbZIP77. These 77 genes were then divided into eleven subfamilies according to their homology with Arabidopsis thaliana. NtbZIPs were unevenly distributed across twenty-two tobacco chromosomes, and we found sixteen pairs of segmental duplication. We further studied the collinearity between these genes and related genes of six other species. Quantitative real-time polymerase chain reaction analysis identified that expression patterns of bZIPs differed, including in different organs and under various abiotic stresses. NtbZIP49 might be important in the development of flowers and fruits; NtbZIP18 might be an important regulator in abiotic stress.

CONCLUSIONS

In this study, the structures and functions of the bZIP family in tobacco were systematically explored. Many bZIPs may play vital roles in the regulation of organ development, growth, and responses to abiotic stresses. This research has great significance for the functional characterisation of the tobacco bZIP family and our understanding of the bZIP family in higher plants.

General gene expression patterns and stemness of the gingiva and dental pulp

Background/purpose

Due to the unique properties of healing processes and cellular differentiation, the gingiva and dental pulp have attracted attention as a potential source of mesenchymal stem cells (MSCs). The purpose of this study was to obtain molecular-level information on these tissues in terms of their function and differentiation processes and investigate stemness.

Materials and methods

Healthy gingival tissues were collected from patients (n = 9; aged 7–12 years) who underwent simple surgical procedures, and normal dental pulp tissues were obtained from patients (n = 25; aged 11–25 years) undergoing tooth extraction for orthodontic reasons. Complementary DNA microarray, qRT-qPCR, and immunohistochemical staining were performed to assess general and MSC gene expression patterns.

Results

In the gingival tissue, genes related to keratinization, the formation of epithelial cells and ectoderm, and immune and/or inflammatory responses were highly expressed. Meanwhile, in the dental pulp tissue, genes related to ion transport, neuronal development and axon guidance, bone and enamel mineralization, extracellular matrix organization, and angiogenesis were highly expressed. When focusing on the expression of MSC genes, induced pluripotent stem (iPS) cell genes, such as Sox2, c-Myc, and KLF4, were expressed at higher levels in the gingival tissue, whereas dental stem cell genes, such as NT5E and VCAM1, were expressed in dental pulp tissue.

Conclusion

We found different general and MSC gene expression patterns between the gingival and dental pulp tissue. These results have implications for future regenerative medicine, considering the application of gingival tissue as a potential source of iPS cells.

Genome-wide and expression analysis of B-box gene family in pepper

BACKGROUND

BBX transcription factors are a kind of zinc finger transcription factors with one or two B-box domains, which partilant in plant growth, development and response to abiotic or biotic stress. The BBX family has been identified in Arabidopsis, rice, tomato and some other model plant genomes.

RESULTS

Here, 24 CaBBX genes were identified in pepper (Capsicum annuum L.), and the phylogenic analysis, structures, chromosomal location, gene expression patterns and subcellular localizations were also carried out to understand the evolution and function of CaBBX genes. All these CaBBXs were divided into five classes, and 20 of them distributed in 11 of 12 pepper chromosomes unevenly. Most duplication events occurred in subgroup I. Quantitative RT-PCR indicated that several CaBBX genes were induced by abiotic stress and hormones, some had tissue-specific expression profiles or differentially expressed at developmental stages. Most of CaBBX members were predicated to be nucleus-localized in consistent with the transient expression assay by onion inner epidermis of the three tested CaBBX members (CaBBX5, 6 and 20).

CONCLUSION

Several CaBBX genes were induced by abiotic stress and exogenous phytohormones, some expressed tissue-specific and variously at different developmental stage. The detected CaBBXs act as nucleus-localized transcription factors. Our data might be a foundation in the identification of CaBBX genes, and a further understanding of their biological function in future studies.

Complement genes in black rockfish (Sebastods schlegelii): genome-wide identification, evolution and their potential functions in response to Vibrio anguillarum infection

As complex components of innate immune system, members of complement system play crucial roles during the process of defensing against pathogens. Black rockfish (Sebastes schlegelii) is one of the important aquaculture species in East Asian. However, studies of complement genes in black rockfish and its related immune activities are still lacking. Therefore, a total of 112 members of the complement genes were identified from the genome of black rockfish and were classified into five subgroups. According to their functional annotations, 30 genes belonged to pattern recognition, 6 genes belonged to proteases, 14 genes belonged to complement components, 36 genes belonged to receptors, and 26 genes belonged to regulators. It can be found that many complement genes evolved into multi-copies, especially in teleost, which may be influenced by whole-genome duplication or tandem duplication events. Complement genes were randomly distributed on 22 chromosomes. The number of introns of complement genes varied from 1 to 70. Results of the expression patterns of 10 randomly selected genes from 5 subtypes response to Vibrio anguillarum infection revealed that most of the members of the complement genes were induced in gill and skin. In contrast, most genes in intestine showed downregulation. This study systematically characterized and analyzed the complement genes in black rockfish and provided new insights into their functions responding to bacterial infection.

Genetic basis of negative heterosis for growth traits in chickens revealed by genome-wide gene expression pattern analysis

Background

Heterosis is an important biological phenomenon that has been extensively utilized in agricultural breeding. However, negative heterosis is also pervasively observed in nature, which can cause unfavorable impacts on production performance. Compared with systematic studies of positive heterosis, the phenomenon of negative heterosis has been largely ignored in genetic studies and breeding programs, and the genetic mechanism of this phenomenon has not been thoroughly elucidated to date. Here, we used chickens, the most common agricultural animals worldwide, to determine the genetic and molecular mechanisms of negative heterosis.

Results

We performed reciprocal crossing experiments with two distinct chicken lines and found that the body weight presented widely negative heterosis in the early growth of chickens. Negative heterosis of carcass traits was more common than positive heterosis, especially breast muscle mass, which was over − 40% in reciprocal progenies. Genome-wide gene expression pattern analyses of breast muscle tissues revealed that nonadditivity, including dominance and overdominace, was the major gene inheritance pattern. Nonadditive genes, including a substantial number of genes encoding ATPase and NADH dehydrogenase, accounted for more than 68% of differentially expressed genes in reciprocal crosses (4257 of 5587 and 3617 of 5243, respectively). Moreover, nonadditive genes were significantly associated with the biological process of oxidative phosphorylation, which is the major metabolic pathway for energy release and animal growth and development. The detection of ATP content and ATPase activity for purebred and crossbred progenies further confirmed that chickens with lower muscle yield had lower ATP concentrations but higher hydrolysis activity, which supported the important role of oxidative phosphorylation in negative heterosis for growth traits in chickens.

Conclusions

These findings revealed that nonadditive genes and their related oxidative phosphorylation were the major genetic and molecular factors in the negative heterosis of growth in chickens, which would be beneficial to future breeding strategies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00574-2.

Genome-wide survey, characterization, and expression analysis of bZIP transcription factors in Chenopodium quinoa

BACKGROUND

Chenopodium quinoa Willd. (quinoa) is a pseudocereal crop of the Amaranthaceae family and represents a promising species with the nutritional content and high tolerance to stressful environments, such as soils affected by high salinity. The basic leucine zipper (bZIP) transcription factor represents exclusively in eukaryotes and can be related to many biological processes. So far, the genomes of quinoa and 3 other Amaranthaceae crops (Spinacia oleracea, Beta vulgaris, and Amaranthus hypochondriacus) have been fully sequenced. However, information about the bZIPs in these Amaranthaceae species is limited, and genome-wide analysis of the bZIP family is lacking in quinoa.

RESULTS

We identified 94 bZIPs in quinoa (named as CqbZIP1-CqbZIP94). All the CqbZIPs were phylogenetically splitted into 12 distinct subfamilies. The proportion of CqbZIPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures and protein motifs. Besides, 32 duplicated CqbZIP gene pairs were investigated, and the duplicated CqbZIPs had mainly undergone purifying selection pressure, which suggested that the functions of the duplicated CqbZIPs might not diverge much. Moreover, we identified the bZIP members in 3 other Amaranthaceae species, and 41, 32, and 16 orthologous gene pairs were identified between quinoa and S. oleracea, B. vulgaris, and A. hypochondriacus, respectively. Among them, most were a single copy being present in S. oleracea, B. vulgaris, and A. hypochondriacus, and two copies being present in allotetraploid quinoa. The function divergence within the bZIP orthologous genes might be limited. Additionally, 11 selected CqbZIPs had specific spatial expression patterns, and 6 of 11 CqbZIPs were up-regulated in response to salt stress. Among the selected CqbZIPs, 3 of 4 duplicated gene pairs shared similar expression patterns, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, motif and gene structure, expansion pattern, and expression profile of the bZIP family in quinoa. Our results would lay an important foundation for functional and evolutionary analysis of CqbZIPs, and provide promising candidate genes for further investigation in tissue specificity and their functional involvement in quinoa's resistance to salt stress.

SQUAMOSA Promoter Binding Protein-Like (SPL) Gene Family: TRANSCRIPTOME-Wide Identification, Phylogenetic Relationship, Expression Patterns and Network Interaction Analysis in Panax ginseng C. A. Meyer

: SPL (SQUAMOSA promoter binding protein-like) gene family is specific transcription factor in the plant that have an important function for plant growth and development. Although the SPL gene family has been widely studied and reported in many various plant species from gymnosperm to angiosperm, there are no systematic studies and reports about the SPL gene family in Panax ginseng C. A. Meyer. In this study, we conducted transcriptome-wide identification, evolutionary analysis, structure analysis, and expression characteristics analysis of SPL gene family in Panax ginseng by bioinformatics. We annotated the PgSPL gene family and found that they might involve in multiple functions including encoding structural proteins, but the main function were still focused on the binding function. The result showed that 106 PgSPL transcripts were classified into two clades - A and B, both of which respectively consisted of three groups. Besides, we profiled PgSPL transcripts' genotypic, temporal, and spatial expression characteristics. Furthermore, we calculated the correlation of PgSPL transcripts in the 14 tissues of a 4 years old ginseng and 42 farmers' cultivars farmers' cultivars of 4 years old ginsengs' roots with both results showing that SPL transcripts formed a single network, which indicated that PgSPLs inter-coordinated when performing their functions. What's more, we found that most PgSPL transcripts tended to express in older ginseng instead of younger ginseng, which was not only reflected in the expression of more types of SPL transcripts in older ginseng, but also in the higher expression of SPL transcripts in older ginseng. Additionally, we found that four PgSPL transcripts were only massively expressed in roots. According to PgSPL transcripts' expression characteristics, we found that PgSPL23-35 and PgSPL24-09 were most proper two transcripts to further study as ginseng age's molecular marker. These results provide the basis for further elucidation of the PgSPL transcripts' biological function in ginseng and ginseng genetics improvement and gene breeding in the future.

Evolution and expression of genes encoding TCP transcription factors in Solanum tuberosum reveal the involvement of StTCP23 in plant defence

Background

The plant-specific Teosinte branched1/Cycloidea/Proliferating cell factor (TCP) family of transcription factors is involved in the regulation of cell growth and proliferation, performing diverse functions in plant growth and development. In addition, TCP transcription factors have recently been shown to be targets of pathogenic effectors and are likely to play a vital role in plant immunity. No comprehensive analysis of the TCP family members in potato (Solanum tuberosum L.) has been undertaken, however, and whether their functions are conserved in potato remains unknown.

Results

To assess TCP gene evolution in potato, we identified TCP-like genes in several publicly available databases. A total of 23 non-redundant TCP transcription factor-encoding genes were identified in the potato genome and subsequently subjected to a systematic analysis that included determination of their phylogenetic relationships, gene structures and expression profiles in different potato tissues under basal conditions and after hormone treatments. These assays also confirmed the function of the class I TCP StTCP23 in the regulation of plant growth and defence.

Conclusions

This is the first genome-wide study including a systematic analysis of the StTCP gene family in potato. Identification of the possible functions of StTCPs in potato growth and defence provides valuable information for our understanding of the classification and functions of the TCP genes in potato.

larvalign: Aligning Gene Expression Patterns from the Larval Brain of Drosophila melanogaster

The larval brain of the fruit fly Drosophila melanogaster is a small, tractable model system for neuroscience. Genes for fluorescent marker proteins can be expressed in defined, spatially restricted neuron populations. Here, we introduce the methods for 1) generating a standard template of the larval central nervous system (CNS), 2) spatial mapping of expression patterns from different larvae into a reference space defined by the standard template. We provide a manually annotated gold standard that serves for evaluation of the registration framework involved in template generation and mapping. A method for registration quality assessment enables the automatic detection of registration errors, and a semi-automatic registration method allows one to correct registrations, which is a prerequisite for a high-quality, curated database of expression patterns. All computational methods are available within the larvalign software package: https://github.com/larvalign/larvalign/releases/tag/v1.0

Quality Control of Single-Cell RNA-seq

Single-cell RNA-seq (scRNA-seq) is emerging as a promising technology to characterize and dissect the cell-to-cell variability. However, the mixture of technical noise and intrinsic biological variability makes separating technical artifacts from real biological variation cells particularly challenging. Proper detection and filtering out technical artifacts before downstream analysis are critical. Here, we present a protocol that integrates both gene expression patterns and data quality to detect technical artifacts in scRNA-seq samples.

Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis)

Circular RNA (circRNA) discovery, expression patterns and experimental validation in developing tea leaves indicates its correlation with circRNA-parental genes and potential roles in ceRNA interaction network. Circular RNAs (circRNAs) have recently emerged as a novel class of abundant endogenous stable RNAs produced by circularization with regulatory potential. However, identification of circRNAs in plants, especially in non-model plants with large genomes, is challenging. In this study, we undertook a systematic identification of circRNAs from different stage tissues of tea plant (Camellia sinensis) leaf development using rRNA-depleted circular RNA-seq. By combining two state-of-the-art detecting tools, we characterized 3174 circRNAs, of which 342 were shared by each approach, and thus considered high-confidence circRNAs. A few predicted circRNAs were randomly chosen, and 20 out of 24 were experimental confirmed by PCR and Sanger sequencing. Similar in other plants, tissue-specific expression was also observed for many C. sinensis circRNAs. In addition, we found that circRNA abundances were positively correlated with the mRNA transcript abundances of their parental genes. qRT-PCR validated the differential expression patterns of circRNAs between leaf bud and young leaf, which also indicated the low expression abundance of circRNAs compared to the standard mRNAs from the parental genes. We predicted the circRNA-microRNA interaction networks, and 54 of the differentially expressed circRNAs were found to have potential tea plant miRNA binding sites. The gene sets encoding circRNAs were significantly enriched in chloroplasts related GO terms and photosynthesis/metabolites biosynthesis related KEGG pathways, suggesting the candidate roles of circRNAs in photosynthetic machinery and metabolites biosynthesis during leaf development.

Expression patterns of heat shock protein genes in Rita rita from natural riverine habitat as biomarker response against environmental pollution

River pollution is one of the principal environmental concerns and biomonitoring tools can play an important role in pollution assessment in the riverine environment. Heat shock proteins (Hsps) have been found to be suitable tools for monitoring stress response. In the present study, expression analyses of hsp genes (hsp27, hsp47, hsp60, hsp70, hsc70, and hsp90) and selected hsp-regulatory genes (hsf1, hyou1, ask1, jnk) were carried out by RT-qPCR in catfish Rita rita collected from selected stretches of river Ganga to investigate changes in their expression patterns as biomarker response. Water quality characteristics were measured in terms of physico-chemical characteristics (DO, BOD, COD, pH, conductivity), element profile (arsenic, mercury, cadmium, lead, chromium, zinc, copper) and persistent organic pollutants (POPs; HCH, DDT, aldrin, endosulphan, heptachlor). Water quality index was calculated and sampling sites were categorized as good/medium/bad. Multivariate analysis was carried out taking the water quality parameters and the fold changes in hsp gene expression as variables, which showed that hsp47 and hsp70b correlated well with BOD, an indicator of organic pollution. To identify the organic pollutant(s) which could be influencing the expression of hsps, again multivariate analysis was employed taking concentration of POPs and fold changes of hsps, which showed up-regulation of hsp47 and hsp70b (HSP72i) correlated well with concentrations of aldrin and HCH. Synergistic effects of these POPs could be responsible for the up-regulation of said hsps, although individually present in low concentration; thus, indicating synergistic effect of the POPs on hsp47 and hsp70b up-regulation as biomarker response.

Genome-wide identification and expression analyses of the pectate lyase (PEL) gene family in cotton (Gossypium hirsutum L.)

Background

Pectin is a major component and structural polysaccharide of the primary cell walls and middle lamella of higher plants. Pectate lyase (PEL, EC 4.2.2.2), a cell wall modification enzyme, degrades de-esterified pectin for cell wall loosening, remodeling and rearrangement. Nevertheless, there have been few studies on PEL genes and no comprehensive analysis of the PEL gene family in cotton.

Results

We identified 53, 42 and 83 putative PEL genes in Gossypium raimondii (D5), Gossypium arboreum (A2), and Gossypium hirsutum (AD1), respectively. These PEL genes were classified into five subfamilies (I-V). Members from the same subfamilies showed relatively conserved gene structures, motifs and protein domains. An analysis of gene chromosomal locations and gene duplication revealed that segmental duplication likely contributed to the expansion of the GhPELs. The 2000 bp upstream sequences of all the GhPELs contained auxin response elements. A transcriptomic data analysis showed that 62 GhPELs were expressed in various tissues. Notably, most (29/32) GhPELs of subfamily IV were preferentially expressed in the stamen, and five GhPELs of subfamily V were prominently expressed at the fiber elongation stage. In addition, qRT-PCR analysis revealed the expression characteristics of 24 GhPELs in four pollen developmental stages and significantly different expression of some GhPELs between long- and short-fiber cultivars. Moreover, some members were responsive to IAA treatment. The results indicate that GhPELs play significant and functionally diverse roles in the development of different tissues.

Conclusions

In this study, we comprehensively analyzed PELs in G. hirsutum, providing a foundation to better understand the functions of GhPELs in different tissues and pathways, especially in pollen, fiber and the auxin signaling pathway.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5047-5) contains supplementary material, which is available to authorized users.

Homoeolog expression bias and expression level dominance in resynthesized allopolyploid Brassica napus

Background

Allopolyploids require rapid genetic and epigenetic modifications to reconcile two or more sets of divergent genomes. To better understand the fate of duplicate genes following genomic mergers and doubling during allopolyploid formation, in this study, we explored the global gene expression patterns in resynthesized allotetraploid Brassica napus (AACC) and its diploid parents B. rapa (AA) and B. oleracea (CC) using RNA sequencing of leaf transcriptomes.

Results

We found that allopolyploid B. napus formation was accompanied by extensive changes (approximately one-third of the expressed genes) in the parental gene expression patterns (‘transcriptome shock’). Interestingly, the majority (85%) of differentially expressed genes (DEGs) were downregulated in the allotetraploid. Moreover, the homoeolog expression bias (relative contribution of homoeologs to the transcriptome) and expression level dominance (total expression level of both homoeologs) were thoroughly investigated by monitoring the expression of 23,766 B. oleracea-B. rapa orthologous gene pairs. Approximately 36.5% of the expressed gene pairs displayed expression bias with a slight preference toward the A-genome. In addition, 39.6, 4.9 and 9.0% of the expressed gene pairs exhibited expression level dominance (ELD), additivity expression and transgressive expression, respectively. The genome-wide ELD was also biased toward the A-genome in the resynthesized B. napus. To explain the ELD phenomenon, we compared the individual homoeolog expression levels relative to those of the diploid parents and found that ELD in the direction of the higher-expression parent can be explained by the downregulation of homoeologs from the dominant parent or upregulation of homoeologs from the nondominant parent; however, ELD in the direction of the lower-expression parent can be explained only by the downregulation of the nondominant parent or both homoeologs. Furthermore, Gene Ontology (GO) enrichment analysis suggested that the alteration in the gene expression patterns could be a prominent cause of the phenotypic variation between the newly formed B. napus and its parental species.

Conclusions

Collectively, our data provide insight into the rapid repatterning of gene expression at the beginning of Brassica allopolyploidization and enhance our knowledge of allopolyploidization processes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4966-5) contains supplementary material, which is available to authorized users.

Expression patterns and mutation analysis of p53 in fish Rita rita from polluted riverine environment

The present study was undertaken to investigate the alterations in gene expression patterns and for mutation analysis of p53 in the riverine catfish Rita rita collected from polluted riverine habitat. The partial p53 gene sequence of Rita rita generated showed a high degree of similarities with the DNA binding domains of fishes, mice and human. Transcriptomic analysis, carried out by quantitative real-time Polymerase Chain Reaction (RT-qPCR), showed significant down-regulation of p53 in fishes collected from most of the polluted stretches. Similar trend in protein abundance was observed by western blot analysis. Down-regulation of p53 was more pronounced in gill than liver. Expression patterns of p53 suggest that exposure to a multitude of contaminants in the natural riverine ecosystem could suppress the expression of p53. Genomic DNA showed a low stained smear pattern upon electrophoresis, with no evidence of DNA fragmentation. For mutation analysis PCR-SSCP followed by sequence analysis was carried out, which identified eight mutations; two at codon level and six missense mutations in the DNA binding domain IV and V. Secondary structure prediction showed that these mutations could lead to impairment of protein structure. Thus, the present study indicated that aquatic pollution has impacted these lower vertebrates which are reflected by the down-regulation of tumor suppressor protein (p53) in majority of the stretches studied.

Visualization of Gene Expression Patterns by In Situ Hybridization on Early Stages of Development of Xenopus laevis

In situ hybridization performed using whole fixed embryos provides accurate and detailed visualization of gene expression patterns. These patterns are useful for investigating spatial patterns of gene expression in normally developing embryos but can also be useful in investigating the effects of genetic or environmental changes on expression of genetic markers characteristic of particular tissues, organs, or genetic pathways. Our lab's protocol for whole-mount in situ hybridization is presented.

RNA In Situ Hybridization on Planarian Paraffin Sections

RNA in situ hybridization techniques are an important tool for the study of gene expression patterns in freshwater planarians. Here I describe a RNA in situ hybridization method on histological paraffin sections of planarian tissue. This protocol allows the visualization of gene expression at cellular or subcellular resolution. Following paraffin-embedding and sectioning of planarians, the resulting sections are hybridized with hapten-labeled RNA probes. Subsequent immunological probe detection is carried out with either chromogenic or fluorescent development. This protocol can be performed alone, or in combination with other immunodetection techniques, and represents a useful alternative to whole-mount protocols more commonly used in the community.

A genome-wide analysis of the small auxin-up RNA (SAUR) gene family in cotton

Background

Small auxin-up RNA (SAUR) gene family is the largest family of early auxin response genes in higher plants, which have been implicated in the regulation of multiple biological processes. However, no comprehensive analysis of SAUR genes has been reported in cotton (Gossypium spp.).

Results

In the study, we identified 145, 97, 214, and 176 SAUR homologous genes in the sequenced genomes of G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. A phylogenetic analysis revealed that the SAUR genes can be classified into 10 groups. A further analysis of chromosomal locations and gene duplications showed that tandem duplication and segmental duplication events contributed to the expansion of the SAUR gene family in cotton. An exon-intron organization and motif analysis revealed the conservation of SAUR-specific domains, and the auxin responsive elements existed in most of the upstream sequences. The expression levels of 16 GhSAUR genes in response to an exogenous application of IAA were determined by a quantitative RT-PCR analysis. The genome-wide RNA-seq data and qRT-PCR analysis of selected SAUR genes in developing fibers revealed their differential expressions. The physical mapping showed that 20 SAUR genes were co-localized with fiber length quantitative trait locus (QTL) hotspots. Single nucleotide polymorphisms (SNPs) were detected for 12 of these 20 genes between G. hirsutum and G. barbadense, but no SNPs were identified between two backcross inbred lines with differing fiber lengths derived from a cross between the two cultivated tetraploids.

Conclusions

This study provides an important piece of genomic information for the SAUR genes in cotton and lays a solid foundation for elucidating the functions of SAUR genes in auxin signaling pathways to regulate cotton growth.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4224-2) contains supplementary material, which is available to authorized users.

Novel messenger RNAs for body fluid identification

In forensic investigations, the identification of the cellular or body fluid source of biological evidence can provide crucial probative information for the court. Messenger RNA (mRNA) profiling has become a valuable tool for body fluid and cell type identification due to its high sensitivity and compatibility with DNA analysis. However, using a single marker to determine the somatic origin of a sample can lead to misinterpretation as a result of cross-reactions. While false positives may be avoided through the simultaneous detection of multiple markers per body fluid, this approach is currently limited by the small number of known differentially expressed mRNAs. Here we characterise six novel mRNAs, partly identified from RNA-Seq, which can supplement existing markers for the detection of circulatory blood, semen (with and without spermatozoa), and menstrual fluid: HBD and SLC4A1 for blood, TNP1 for spermatozoa, KLK2 for seminal fluid, and MMP3 and STC1 for menstrual fluid. Respective expression profiles were evaluated by singleplex endpoint reverse transcription polymerase chain reaction (RT-PCR). HBD, SLC4A1, and KLK2 were specific to their target body fluids. TNP1, MMP3, and STC1 each cross-reacted with two non-target samples; however, these signals were below 350RFU, not reproducible, and likely resulted from large body fluid inputs. All candidates were more sensitive for the detection of their target body fluids than corresponding well-known mRNAs, in particular those for menstrual fluid. The increased sensitivities were statistically significant, except for KLK2. Thus, the new mRNAs introduced here are promising new targets for improved body fluid profiling.

Genome-wide identification, phylogeny and expression analyses of SCARECROW-LIKE(SCL) genes in millet (Setaria italica)

As a member of the GRAS gene family, SCARECROW-LIKE (SCL) genes encode transcriptional regulators that are involved in plant information transmission and signal transduction. In this study, 44 SCL genes including two SCARECROW genes in millet were identified to be distributed on eight chromosomes, except chromosome 6. All the millet genes contain motifs 6-8, indicating that these motifs are conserved during the evolution. SCL genes of millet were divided into eight groups based on the phylogenetic relationship and classification of Arabidopsis SCL genes. Several putative millet orthologous genes in Arabidopsis, maize and rice were identified. High throughput RNA sequencing revealed that the expressions of millet SCL genes in root, stem, leaf, spica, and along leaf gradient varied greatly. Analyses combining the gene expression patterns, gene structures, motif compositions, promoter cis-elements identification, alternative splicing of transcripts and phylogenetic relationship of SCL genes indicate that the these genes may play diverse functions. Functionally characterized SCL genes in maize, rice and Arabidopsis would provide us some clues for future characterization of their homologues in millet. To the best of our knowledge, this is the first study of millet SCL genes at the genome wide level. Our work provides a useful platform for functional analysis of SCL genes in millet, a model crop for C₄ photosynthesis and bioenergy studies.

A genome-wide analysis of the lysophosphatidate acyltransferase (LPAAT) gene family in cotton: organization, expression, sequence variation, and association with seed oil content and fiber quality

BACKGROUND

Lysophosphatidic acid acyltransferase (LPAAT) encoded by a multigene family is a rate-limiting enzyme in the Kennedy pathway in higher plants. Cotton is the most important natural fiber crop and one of the most important oilseed crops. However, little is known on genes coding for LPAATs involved in oil biosynthesis with regard to its genome organization, diversity, expression, natural genetic variation, and association with fiber development and oil content in cotton.

RESULTS

In this study, a comprehensive genome-wide analysis in four Gossypium species with genome sequences, i.e., tetraploid G. hirsutum- AD₁ and G. barbadense- AD₂ and its possible ancestral diploids G. raimondii- D₅ and G. arboreum- A₂, identified 13, 10, 8, and 9 LPAAT genes, respectively, that were divided into four subfamilies. RNA-seq analyses of the LPAAT genes in the widely grown G. hirsutum suggest their differential expression at the transcriptional level in developing cottonseeds and fibers. Although 10 LPAAT genes were co-localised with quantitative trait loci (QTL) for cottonseed oil or protein content within a 25-cM region, only one single strand conformation polymorphic (SSCP) marker developed from a synonymous single nucleotide polymorphism (SNP) of the At-Gh13LPAAT5 gene was significantly correlated with cottonseed oil and protein contents in one of the three field tests. Moreover, transformed yeasts using the At-Gh13LPAAT5 gene with the two sequences for the SNP led to similar results, i.e., a 25-31% increase in palmitic acid and oleic acid, and a 16-29% increase in total triacylglycerol (TAG).

CONCLUSIONS

The results in this study demonstrated that the natural variation in the LPAAT genes to improving cottonseed oil content and fiber quality is limited; therefore, traditional cross breeding should not expect much progress in improving cottonseed oil content or fiber quality through a marker-assisted selection for the LPAAT genes. However, enhancing the expression of one of the LPAAT genes such as At-Gh13LPAAT5 can significantly increase the production of total TAG and other fatty acids, providing an incentive for further studies into the use of LPAAT genes to increase cottonseed oil content through biotechnology.

Identification and characterization of a novel Toll-like receptor 2 homologue in the large yellow croaker Larimichthys crocea

Toll-like receptors (TLRs) are key components of innate immunity that play significant roles in immune defence against pathogen invasion. In the present study, we identified a novel TLR2 homologue (LycTLR2b) in large yellow croaker (Larimichthys crocea) that shared low sequence identity with the previously reported large yellow croaker TLR2 (tentatively named LycTLR2a). The full-length cDNA of LycTLR2b was 2926 nucleotides (nt) long and encoded a protein consisting of 797 amino acids (aa). The deduced LycTLR2b protein exhibited a typical TLR domain architecture including a signal peptide, seven leucine-rich repeats (LRRs) in the extracellular region, a transmembrane domain, and a Toll-Interleukin 1 receptor (TIR) domain in the cytoplasmic region. Phylogenetic analysis showed that both LycTLR2a and LycTLR2b fall into a major clade formed by all TLR2 sequences, and are divided into two distinct branches. Genomic organization revealed that the LycTLR2b gene lacks intron, which is similar to zebrafish and human TLR2 genes, whereas the LycTLR2a gene contains multiple introns, as found in damselfish TLR2a and Fugu TLR2 genes. Syntenic analysis suggested that the occurrence of LycTLR2a and LycTLR2b may result from a relatively recent genome duplication event. LycTLR2b mRNA was constitutively expressed in all tissues examined although at different levels. Following bacterial vaccine challenge, LycTLR2b expression levels were significantly up-regulated in both spleen and head kidney tissues. Taken together, these results indicated that two different TLR2 homologues, which may play roles in antibacterial immunity, exist in large yellow croaker.

Global analysis of gene expression and projection target correlations in the mouse brain

Recent studies have shown that projection targets in the mouse neocortex are correlated with their gene expression patterns. However, a brain-wide quantitative analysis of the relationship between voxel genetic composition and their projection targets is lacking to date. Here we extended those studies to perform a global, integrative analysis of gene expression and projection target correlations in the mouse brain. By using the Allen Brain Atlas data, we analyzed the relationship between gene expression and projection targets. We first visualized and clustered the two data sets separately and showed that they both exhibit strong spatial autocorrelation. Building upon this initial analysis, we conducted an integrative correlation analysis of the two data sets while correcting for their spatial autocorrelation. This resulted in a correlation of 0.19 with significant p value. We further identified the top genes responsible for this correlation using two greedy gene ranking techniques. Using only the top genes identified by those techniques, we recomputed the correlation between these two data sets. This led to correlation values up to 0.49 with significant p values. Our results illustrated that although the target specificity of neurons is in fact complex and diverse, yet they are strongly affected by their genetic and molecular compositions.

High-resolution prediction of mouse brain connectivity using gene expression patterns

The brain is a multi-level system in which the high-level functions are generated by low-level genetic mechanisms. Thus, elucidating the relationship among multiple brain levels via correlative and predictive analytics is an important area in brain research. Currently, studies in multiple species have indicated that the spatiotemporal gene expression patterns are predictive of brain wiring. Specifically, results on the worm Caenorhabditis elegans have shown that the prediction of neuronal connectivity using gene expression signatures yielded statistically significant results. Recent studies on the mammalian brain produced similar results at the coarse regional level. In this study, we provide the first high-resolution, large-scale integrative analysis of the transcriptome and connectome in a single mammalian brain at a fine voxel level. By using the Allen Brain Atlas data, we predict voxel-level brain connectivity based on the gene expressions in the adult mouse brain. We employ regularized models to show that gene expression is predictive of connectivity at the voxel-level with an accuracy of 93%. We also identify a set of genes playing the most important role in connectivity prediction. We use only this small number of genes to predict the brain wiring with an accuracy over 80%. We discover that these important genes are enriched in neurons as compared to glia, and they perform connectivity-related functions. We perform several interesting correlative studies to further elucidate the transcriptome-connectome relationship.