Vitellogenin and vitellogenin receptor gene expression is associated with male and female parenting in a subsocial insect

Complex social behaviour in Hymenoptera has been hypothesized to evolve by co-opting reproductive pathways (the ovarian ground plan hypothesis, OGPH) and gene networks (the reproductive ground plan hypothesis, RGPH). In support of these hypotheses, in eusocial Hymenoptera where there is reproductive division of labour, the yolk precursor protein vitellogenin (Vg) influences the expression of worker social behaviour. We suggest that co-opting genes involved in reproduction may occur more generally than just in the evolution of eusociality; i.e. underlie earlier stages of social evolution such as the evolution of parental care, given that reproduction and parental care rarely overlap. We therefore examined vitellogenin (vg) gene expression associated with parental care in the subsocial beetle Nicrophorus vespilloides. We found a significant reduction in the expression of vg and its receptor, vgr, in head tissue during active parental care, and confirmed that the receptor is expressed in the brains of both sexes. Ours is the first study to show that vgr is expressed in the brain of a non-eusocial insect. Given the association between behaviour and gene expression in both sexes, and the presence of vitellogenin receptors in the brain, we suggest that Vg was co-opted early in the evolution of sociality to have a regulatory function. This extends the association of Vg in parenting to subsocial species and outside of the Hymenoptera, and supports the hypothesis that the OGPH is general and that heterochrony in gene expression is important in the evolution of social behaviour and precedes subsequent evolutionary specialization of social roles.     

Hypervariable genes—experimental error or hidden dynamics

In a homogeneous group of samples, not all genes of high variability stem from experimental errors in microarray experiments. These expression variations can be attributed to many factors including natural biological oscillations or metabolic processes. The behavior of these genes can tease out important clues about naturally occurring dynamic processes in the organism or experimental system under study. We developed a statistical procedure for the selection of genes with high variability denoted hypervariable (HV) genes. After the exclusion of low expressed genes and a stabilizing log-transformation, the majority of genes have comparable residual variability. Based on an F-test, HV genes are selected as having a statistically significant difference from the majority of variability stabilized genes measured by the ‘reference group’. A novel F-test clustering technique, further noted as ‘F-means clustering’, groups HV genes with similar variability patterns, presumably from their participation in a common dynamic biological process. F-means clustering establishes, for the first time, groups of co-expressed HV genes and is illustrated with microarray data from patients with juvenile rheumatoid arthritis and healthy controls.     

Adaptation of Musca domestica L. Field Population to Laboratory Breeding Causes Transcriptional Alterations

Background

The housefly, Musca domestica, has developed resistance to most insecticides applied for its control. Expression of genes coding for detoxification enzymes play a role in the response of the housefly when encountered by a xenobiotic. The highest level of constitutive gene expression of nine P450 genes was previously found in a newly-collected susceptible field population in comparison to three insecticide-resistant laboratory strains and a laboratory reference strain.

Results

We compared gene expression of five P450s by qPCR as well as global gene expression by RNAseq in the newly-acquired field population (845b) in generation F₁, F₁₃ and F₂₉ to test how gene expression changes following laboratory adaption. Four (CYP6A1, CYP6A36, CYP6D3, CYP6G4) of five investigated P450 genes adapted to breeding by decreasing expression. CYP6D1 showed higher female expression in F₂₉ than in F₁. For males, about half of the genes accessed in the global gene expression were up-regulated in F₁₃ and F₂₉ in comparison with the F₁ population. In females, 60% of the genes were up-regulated in F₁₃ in comparison with F₁, while 33% were up-regulated in F₂₉. Forty potential P450 genes were identified. In most cases, P450 gene expression was decreased in F₁₃ flies in comparison with F₁. Gene expression then increased from F₁₃ to F₂₉ in males and decreased further in females.

Conclusion

The global gene expression changes massively during adaptation to laboratory breeding. In general, global expression decreased as a result of laboratory adaption in males, while female expression was not unidirectional. Expression of P450 genes was in general down-regulated as a result of laboratory adaption. Expression of hexamerin, coding for a storage protein was increased, while gene expression of genes coding for amylases decreased. This suggests a major impact of the surrounding environment on gene response to xenobiotics and genetic composition of housefly strains.     

Multi-Tissue Microarray Analysis Identifies a Molecular Signature of Regeneration

The inability to functionally repair tissues that are lost as a consequence of disease or injury remains a significant challenge for regenerative medicine. The molecular and cellular processes involved in complete restoration of tissue architecture and function are expected to be complex and remain largely unknown. Unlike humans, certain salamanders can completely regenerate injured tissues and lost appendages without scar formation. A parsimonious hypothesis would predict that all of these regenerative activities are regulated, at least in part, by a common set of genes. To test this hypothesis and identify genes that might control conserved regenerative processes, we performed a comprehensive microarray analysis of the early regenerative response in five regeneration-competent tissues from the newt Notophthalmus viridescens. Consistent with this hypothesis, we established a molecular signature for regeneration that consists of common genes or gene family members that exhibit dynamic differential regulation during regeneration in multiple tissue types. These genes include members of the matrix metalloproteinase family and its regulators, extracellular matrix components, genes involved in controlling cytoskeleton dynamics, and a variety of immune response factors. Gene Ontology term enrichment analysis validated and supported their functional activities in conserved regenerative processes. Surprisingly, dendrogram clustering and RadViz classification also revealed that each regenerative tissue had its own unique temporal expression profile, pointing to an inherent tissue-specific regenerative gene program. These new findings demand a reconsideration of how we conceptualize regenerative processes and how we devise new strategies for regenerative medicine.     

Expression patterns of WRKY genes in di-haploid Populus simonii × P. nigra in response to salinity stress revealed by quantitative real-time PCR and RNA sequencing

Key message

Spatio-temporal expression patterns of 13 out of 119 poplar WRKY genes indicated dynamic and tissue-specific roles of WRKY family proteins in salinity stress tolerance.

Abstract

To understand the expression patterns of poplar WRKY genes under salinity stress, 51 of the 119 WRKY genes were selected from di-haploid Populus simonii × P. nigra by quantitative real-time PCR (qRT-PCR). We used qRT-PCR to profile the expression of the top 13 genes under salinity stress across seven time points, and employed RNA-Seq platforms to cross-validate it. Results demonstrated that all the 13 WRKY genes were expressed in root, stem, and leaf tissues, but their expression levels and overall patterns varied notably in these tissues. Regarding overall gene expression in roots, the 13 genes were significantly highly expressed at all six time points after the treatment, reaching the plateau of expression at hour 9. In leaves, the 13 genes were similarly up-regulated from 3 to 12 h in response to NaCl treatment. In stems, however, expression levels of the 13 genes did not show significant changes after the NaCl treatment. Regarding individual gene expression across the time points and the three tissues, the 13 genes can be classified into three clusters: the lowly expressed Cluster 1 containing PthWRKY28, 45 and 105; intermediately expressed Clusters 2 including PthWRKY56, 88 and 116; and highly expressed Cluster 3 consisting of PthWRKY41, 44, 51, 61, 62, 75 and 106. In general, genes in Cluster 2 and 3 displayed a dynamic pattern of “induced amplification—recovering”, suggesting that these WRKY genes and corresponding pathways may play a critical role in mediating salt response and tolerance in a dynamic and tissue-specific manner.     

Expression of cell cycle regulatory factors hus1, gadd45a, rb1, cdkn2a and mre11a correlates with expression of clock gene per2 in human colorectal carcinoma tissue

Deregulated expression of clock gene per2 has previously been associated with progression of cancer. The aim of the present study was to identify genes related to per2 expression and involved in cell cycle control. Patients surgically treated for colorectal carcinoma with up-regulated and down-regulated per2 expression in cancer versus adjacent tissue were studied. Total RNA from cancer tissue of these patients was used to specify genes associated with altered per2 expression using the Human Cell Cycle RT² profiler PCR array system. We identified seven genes positively correlated (hus1, gadd45α, rb1, cdkn2a, cdk5rp1, mre11a, sumo1) and two genes negatively correlated (cdc20, birc5) with per2 expression. Expression of these seven genes was subsequently measured by real time PCR in all patients of the cohort. Patients were divided into three groups according to TNM classification. We observed an increase in gene expression in cancer tissue compared to adjacent tissue in the first group of patients in all genes measured. Expression of genes positively associated with per2 gene expression was dependent on tumor staging and changes were observed preferentially in cancer tissue. For genes negatively associated with per2 expression we also detected changes in expression dependent on tumor staging. Expression of cdc20 and birc5 was increasing in the proximal tissue and decreasing in the cancer tissue. These results implicate functional involvement of per2 in the process of carcinogenesis via newly uncovered genes. The relevancy of gene expression for determination of diagnosis and prognosis should be considered in relation to tumor staging.     

Independence of genetic variation between circadian rhythm and development time in the seed beetle, Callosobruchus chinensis

A positive genetic correlation between periods of circadian rhythm and developmental time supports the hypothesis that circadian clocks are implicated in the timing of development. Empirical evidence for this genetic correlation in insects has been documented in two fly species. In contrast, here we show that there is no evidence of genetic correlation between circadian rhythm and development time in the adzuki bean beetle, Callosobruchus chinensis. This species has variation that is explained by a major gene in the expression and period length of circadian rhythm between strains. In this study, we found genetic variation in development time between the strains. The development time was not covaried with either the incidence or the period length of circadian rhythm among the strains. Crosses between strains suggest that development time is controlled by a polygene. In the F₂ individuals from the crosses, the circadian rhythm is attributable to allelic variation in the major gene. Across the F₂ individuals, development time was not correlated with either the expression or the period length of circadian rhythm. Thus, we found no effects of major genes responsible for variation in the circadian rhythm on development time in C. chinensis. Our findings collectively give no support to the hypothesis that the circadian clock is involved in the regulation of development time in this species.     

Effect of dietary fish oil on selected inflammatory markers in pigs

The present study tested a hypothesis that dietary fish oil (eicosapentaenoic acid+docosahexaenoic acid) in a commonly achievable dose ameliorates a systemic inflammation in pigs. Two groups of pigs of 16 animals each were fed a diet with either 2.5% of fish oil (F) or a control diet with 2.5% of palm oil (P). After 70 days of fattening, eight F and eight P pigs were challenged (F⁺; P⁺) i.v. by lipopolysaccharide. After 3 h, all pigs were sacrificed and blood, liver and visceral adipose tissue (VAT) samples were taken. No significant effect (P>0.05) of dietary oil on the feed intake and daily weight gain was found out. Less neutrophils (16.8% v. 28.8%; P<0.05) were found in the F⁺-leukocytes of the peripheral blood; F⁺ pigs had lower (P<0.05) percentage of the swine leukocyte antigen-D-related CD163⁺ (SLA-DR⁺ CD163⁺) macrophages in the VAT (15.4% v. 21.8%) and lower expression of the SLA-DR-CD163⁺ surface molecules of the VAT macrophages. No difference (P>0.05) between F⁺ and P⁺ pigs in the peroxisome proliferator-activated receptor γ, GPR120, Adipor1 and Adipor2 (adiponectin receptor) gene expression, respectively, was established; plasma adiponectin was the same (21.1 ng/ml) in F⁺ and P⁺ pigs. In comparison with the P⁺ pigs, increased expression of the lipopolysaccharide-binding protein (LBP) gene and intercellular adhesion molecule 1 (ICAM1) gene was found out in the liver of the F⁺ pigs; expression of the tumor necrosis factor α (TNFα) gene was higher in the liver but lower in the VAT of the F⁺ pigs (P<0.05). The F⁺ pigs had higher (P<0.05) plasma concentration of both anti-inflammatory cytokine interleukin-4 (0.46 v. 0.04 ng/ml) and pro-inflammatory TNF-α (13.41 v. 7.72 ng/ml). It was concluded that dietary fish oil at the tested amount had a negligible effect on expression of the evaluated receptor genes and plasma adiponectin, and had an ambiguous effect on expression of cytokine genes and plasma cytokine levels.     

Correlations of expression of cell wall biosynthesis genes with variation in biomass composition in shrub willow (Salix spp.) biomass crops

We have measured significant genetically determined variation in biomass composition among breeding populations of shrub willow, a biomass feedstock crop. This project was aimed to ask whether patterns of cell wall gene expression can be correlated with genetic variation in biomass composition at harvest, in order to develop assays of early differences in gene expression as indicators of harvestable biomass chemical composition and potentially reduce the time of selection for new willow genotypes. Previous studies have demonstrated that manipulation of expression of cell wall biosynthetic genes results in altered biomass chemical composition. We analyzed genes encoding enzymes involved in lignin biosynthesis and carbohydrate active enzymes selected based on their functional characterization and conservation in Populus trichocarpa and Arabidopsis thaliana. Fragments of 20 genes were cloned from young stem cDNA of Salix sachalinensis and Salix miyabeana. Expression profiling in willow stem apical tissue and developing stem tissue was performed for each isolated gene using probe-based quantitative real-time PCR. Two willow parental genotypes and six progeny within a hybrid family were selected for analysis, and significant differences in expression among the individuals and between tissue types were observed for most of the genes. Significant correlations between patterns of gene expression and variation in the biomass chemical composition of those genotypes provide insight into the genetic regulation of lignocellulosic deposition in this important bioenergy crop and could be utilized as a tool for early selection of new genotypes.     

Zonation of Nitrogen and Glucose Metabolism Gene Expression upon Acute Liver Damage in Mouse

Zonation of metabolic activities within specific structures and cell types is a phenomenon of liver organization and ensures complementarity of variant liver functions like protein production, glucose homeostasis and detoxification. To analyze damage and regeneration of liver tissue in response to a toxic agent, expression of liver specific enzymes was analyzed by in situ hybridization in mouse over a 6 days time course following carbon tetrachloride (CCl₄) injection. CCl₄ mixed with mineral oil was administered to BALB/c mice by intraperitoneal injection, and mice were sacrificed at different time points post injection. Changes in the expression of albumin (Alb), arginase (Arg1), glutaminase 2 (Gls2), Glutamine synthetase (Gs), glucose-6-phosphatase (G6pc), glycogen synthase 2 (Gys2), Glycerinaldehyd-3-phosphat-Dehydrogenase (Gapdh), Cytochrom p450 2E1 (Cyp2e1) and glucagon receptor (Gcgr) genes in the liver were studied by in situ hybridization and qPCR. We observed significant changes in gene expression of enzymes involved in nitrogen and glucose metabolism and their local distribution following CCl₄ injury. We also found that Cyp2e1, the primary metabolizing enzyme for CCl₄, was strongly expressed in the pericentral zone during recovery. Furthermore, cells in the damaged area displayed distinct gene expression profiles during the analyzed time course and showed complete recovery with strong albumin production 6 days after CCl₄ injection. Our results indicate that despite severe damage, liver cells in the damaged area do not simply die but instead display locally adjusted gene expression supporting damage response and recovery.     

Analysis of the real EADGENE data set: Multivariate approaches and post analysis (Open Access publication)

The aim of this paper was to describe, and when possible compare, the multivariate methods used by the participants in the EADGENE WP1.4 workshop. The first approach was for class discovery and class prediction using evidence from the data at hand. Several teams used hierarchical clustering (HC) or principal component analysis (PCA) to identify groups of differentially expressed genes with a similar expression pattern over time points and infective agent (E. coli or S. aureus). The main result from these analyses was that HC and PCA were able to separate tissue samples taken at 24 h following E. coli infection from the other samples. The second approach identified groups of differentially co-expressed genes, by identifying clusters of genes highly correlated when animals were infected with E. coli but not correlated more than expected by chance when the infective pathogen was S. aureus. The third approach looked at differential expression of predefined gene sets. Gene sets were defined based on information retrieved from biological databases such as Gene Ontology. Based on these annotation sources the teams used either the GlobalTest or the Fisher exact test to identify differentially expressed gene sets. The main result from these analyses was that gene sets involved in immune defence responses were differentially expressed.     

Temporal gene expression profile after acute electroconvulsive stimulation in the rat

Electroconvulsive therapy (ECT) remains one of the most effective treatments of major depression. It has been suggested that the mechanisms of action involve gene expression. In recent decades there have been several investigations of gene expression following both acute and chronic electroconvulsive stimulation (ECS). These studies have focused on several distinct gene targets but have generally included only few time points after ECS for measuring gene expression. Here we measured gene expression of three types of genes: Immediate early genes, synaptic proteins, and neuropeptides at six time points following an acute ECS. We find significant increases for c-Fos, Egr1, Neuritin 1 (Nrn 1), Bdnf, Snap29, Synaptotagmin III (Syt 3), Synapsin I (Syn 1), and Psd95 at differing time points after ECS. For some genes these changes are prolonged whereas for others they are transient. Npy expression significantly increases whereas the gene expression of its receptors Npy1r, Npy2r, and Npy5r initially decreases. These decreases are followed by a significant increase for Npy2r, suggesting anticonvulsive adaptations following seizures. In summary, we find distinct changes in mRNA quantities that are characteristic for each gene. Considering the observed transitory and inverse changes in expression patterns, these data underline the importance of conducting measurements at several time points post-ECS.     

Effects of Warm Ischemic Time on Gene Expression Profiling in Colorectal Cancer Tissues and Normal Mucosa

Background

Genome-wide gene expression analyses of tumors are a powerful tool to identify gene signatures associated with biologically and clinically relevant characteristics and for several tumor types are under clinical validation by prospective trials. However, handling and processing of clinical specimens may significantly affect the molecular data obtained from their analysis. We studied the effects of tissue handling time on gene expression in human normal and tumor colon tissues undergoing routine surgical procedures.

Methods

RNA extracted from specimens of 15 patients at four time points (for a total of 180 samples) after surgery was analyzed for gene expression on high-density oligonucleotide microarrays. A mixed-effects model was used to identify probes with different expression means across the four different time points. The p-values of the model were adjusted with the Bonferroni method.

Results

Thirty-two probe sets associated with tissue handling time in the tumor specimens, and thirty-one in the normal tissues, were identified. Most genes exhibited moderate changes in expression over the time points analyzed; however four of them were oncogenes, and two confirmed the effect of tissue handling by independent validation.

Conclusions

Our results suggest that a critical time point for tissue handling in colon seems to be 60 minutes at room temperature. Although the number of time-dependent genes we identified was low, the three genes that already showed changes at this time point in tumor samples were all oncogenes, hence recommending standardization of tissue-handling protocols and effort to reduce the time from specimen removal to snap freezing accounting for warm ischemia in this tumor type.     

Isolation of mRNA from specific tissues of Drosophila by mRNA tagging

To study the function of specific cells or tissues using genomic tools like microarray analyses, it is highly desirable to obtain mRNA from a homogeneous source. However, this is particularly challenging for small organisms, like Caenorhabditis elegans and Drosophila melanogaster. We have optimized and applied a new technique, mRNA tagging, to isolate mRNA from specific tissues of D.melanogaster. A FLAG-tagged poly(A)-binding protein (PABP) is expressed in a specific tissue and mRNA from that tissue is thus tagged by the recombinant PABP and separated from mRNA in other tissues by co-immunoprecipitation with a FLAG-tag specific antibody. The fractionated mRNA is then amplified and used as probe in microarray experiments. As a test system, we employed the procedures to identify genes expressed in Drosophila photoreceptor cells. We found that most known photoreceptor cell-specific mRNAs were identified by mRNA tagging. Furthermore, at least 11 novel genes have been identified as enriched in photoreceptor cells. mRNA tagging is a powerful general method for profiling gene expression in specific tissues and for identifying tissue-specific genes.