Genome-wide analysis of the chromatin sites targeted by HEX 70a storage protein in the honeybee brain and fat body

Hexamerins, the proteins massively stored in the larval haemolymph of insects, are gradually used throughout metamorphosis as a source of raw material and energy for the development of adult tissues. Such behaviour defined hexamerins as storage proteins. Immunofluorescence experiments coupled with confocal microscopy show a hexamerin, HEX 70a, in the nucleus of the brain and fat body cells from honeybee workers, an unexpected localization for a storage protein. HEX 70a colocalizes with fibrillarin, a nucleolar-specific protein and H3 histone, thus suggesting a potential role as a chromatin-binding protein. This was investigated through chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq). The significant HEX 70a-DNA binding sites were mainly localized at the intergenic, promoter and intronic regions. HEX 70a targeted DNA stretches mapped to the genomic regions encompassing genes with relevant functional attributes. Several HEX 70a targeted genes were associated with H3K27ac or/and H3K27me3, known as active and repressive histone marks. Brain and fat body tissues shared a fraction of the HEX 70 targeted genes, and tissue-specific targets were also detected. The presence of overrepresented DNA motifs in the binding sites is consistent with specific HEX 70a-chromatin association. In addition, a search for HEX 70a targets in RNA-seq public libraries of fat bodies from nurses and foragers revealed differentially expressed targets displaying hex 70a-correlated developmental expression, thus supporting a regulatory activity for HEX 70a. Our results support the premise that HEX 70a is a moonlighting protein that binds chromatin and has roles in the brain and fat body cell nuclei, apart from its canonical role as a storage protein.

The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae)

Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts' nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.

Worker bees (Apis mellifera) deprived of pollen in the first week of adulthood exhibit signs of premature aging

Pollinator populations, including bees, are in rapid decline in many parts of the world, raising concerns over the future of ecosystems and food production. Among the factors involved in these declines, poor nutrition deserves attention. The diet consumed by adult worker honeybees (Apis mellifera) is crucial for their behavioral maturation, i.e., the progressive division of labor they perform, such as nurse bees initially and later in life as foragers. Poor pollen nutrition is known to reduce the workers' lifespan, but the underlying physiological and genetic mechanisms are not fully understood. Here we investigate how the lack of pollen in the diet of workers during their first week of adult life can affect age-related phenotypes. During the first seven days of adult life, newly emerged workers were fed either a pollen-deprived (PD) diet mimicking that of an older bee, or a control pollen-rich (PR) diet, as typically consumed by young bees. The PD-fed bees showed alterations in their fat body transcriptome, such as a switch from a protein-lipid based metabolism to a carbohydrate-based metabolism, and a reduced expression of genes involved with immune response. The absence of pollen in the diet also led to an accumulation of oxidative stress markers in fat body tissue and alterations in the cuticular hydrocarbon profiles, which became similar to those of chronologically older bees. Together, our data indicate that the absence of pollen during first week of adulthood triggers the premature onset of an aging-related worker phenotype.

Identification of trailer heat zones and associated heat stress in weaner pigs transported by road in tropical climates

Heat stress is major welfare concern during transport of pigs in tropical climates, which can also lead to direct production costs. This study evaluated the dynamics of heat zones through the load and their relationship with heat stress of weaner pigs during road transport in a tropical climate. Both environmental (e.g. temperature and relative humidity) and physiological (e.g. respiratory frequency and lactate) measures were recorded from four vehicle journeys (70 km distance, 216 weaner pigs within each trailer load) within Ceará, northeastern Brazil. Geostatistics and fluid dynamics simulation techniques were applied to understand the dynamics of heat zones and ventilation patterns the truckload. Statistics based on canonical discriminant analysis and ANOVA were performed to verify the relationship between heat zones and heat stress in pigs. The results showed that, during transport, the generation of heat zones occurred with different magnitudes along the load (P < 0.05), which was harmonized by the ventilation dynamics. There was a heat core with high energy content, in the front region of the lower deck (LD) of the trailer. In this zone, weaners pigs had higher rectal temperature (+1.8 °C temperature difference), respiratory frequency (LD = 94 ± 1.3 breaths/min; UD = 86 ± 1.3 breaths/min), and blood cortisol concentration (LD = 32.9 ± 0.8 ng/mL; UD = 30.18 ± 0.6 ng/mL) (all P < 0.05). Weaners pigs transported in the upper deck (UD) compartments had the highest skin temperature (LD = 38.13 ± 0.3 °C; UD = 38.9 ± 0.22 °C) and the highest mean values of blood lactate (LD = 65.5 ± 1.11 m/M; UD = 71.60 ± 1.19 m/M) and Creatine kinase (LD = 3891.23 ± 69U/L; UD = 4107.43 ± 62U/L) (P < 0.05). Weaners transported in compartments of the LD of trailer were more susceptible to heat stress, while weaners in the UD compartments were more susceptible to physical stress and muscle exhaustion. These results provide additional evidence of heat zones within trailer compartments and highlight the requirement for the planning of pig transport operations in tropical climates to mitigate risks of heat stress.

Caste-specific gene expression underlying the differential adult brain development in the honeybee Apis mellifera

Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynβ). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories.

Spliced genes in muscle from Nelore Cattle and their association with carcass and meat quality

Transcript data obtained by RNA-Seq were used to identify differentially expressed alternatively spliced genes in ribeye muscle tissue between Nelore cattle that differed in their ribeye area (REA) or intramuscular fat content (IF). A total of 166 alternatively spliced transcripts from 125 genes were significantly differentially expressed in ribeye muscle between the highest and lowest REA groups (p ≤ 0.05). For animals selected on their IF content, 269 alternatively spliced transcripts from 219 genes were differentially expressed in ribeye muscle between the highest and lowest IF animals. Cassette exons and alternative 3′ splice sites were the most frequently found alternatively spliced transcripts for REA and IF content. For both traits, some differentially expressed alternatively spliced transcripts belonged to myosin and myotilin gene families. The hub transcripts were identified for REA (LRRFIP1, RCAN1 and RHOBTB1) and IF (TRIP12, HSPE1 and MAP2K6) have an important role to play in muscle cell degradation, development and motility. In general, transcripts were found for both traits with biological process GO terms that were involved in pathways related to protein ubiquitination, muscle differentiation, lipids and hormonal systems. Our results reinforce the biological importance of these known processes but also reveal new insights into the complexity of the whole cell muscle mRNA of Nelore cattle.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

Intraspecific Variation within the Utricularia amethystina Species Morphotypes Based on Chloroplast Genomes

Utricularia amethystina Salzm. ex A.St.-Hil. & Girard (Lentibulariaceae) is a highly polymorphic carnivorous plant taxonomically rearranged many times throughout history. Herein, the complete chloroplast genomes (cpDNA) of three U. amethystina morphotypes: purple-, white-, and yellow-flowered, were sequenced, compared, and putative markers for systematic, populations, and evolutionary studies were uncovered. In addition, RNA-Seq and RNA-editing analysis were employed for functional cpDNA evaluation. The cpDNA of three U. amethystina morphotypes exhibits typical quadripartite structure. Fine-grained sequence comparison revealed a high degree of intraspecific genetic variability in all morphotypes, including an exclusive inversion in the psbM and petN genes in U. amethystina yellow. Phylogenetic analyses indicate that U. amethystina morphotypes are monophyletic. Furthermore, in contrast to the terrestrial Utricularia reniformis cpDNA, the U. amethystina morphotypes retain all the plastid NAD(P)H-dehydrogenase (ndh) complex genes. This observation supports the hypothesis that the ndhs in terrestrial Utricularia were independently lost and regained, also suggesting that different habitats (aquatic and terrestrial) are not related to the absence of Utricularia ndhs gene repertoire as previously assumed. Moreover, RNA-Seq analyses recovered similar patterns, including nonsynonymous RNA-editing sites (e.g., rps14 and petB). Collectively, our results bring new insights into the chloroplast genome architecture and evolution of the photosynthesis machinery in the Lentibulariaceae.

Exploring integument transcriptomes, cuticle ultrastructure, and cuticular hydrocarbons profiles in eusocial and solitary bee species displaying heterochronic adult cuticle maturation

Differences in the timing of exoskeleton melanization and sclerotization are evident when comparing eusocial and solitary bees. This cuticular maturation heterochrony may be associated with life style, considering that eusocial bees remain protected inside the nest for many days after emergence, while the solitary bees immediately start outside activities. To address this issue, we characterized gene expression using large-scale RNA sequencing (RNA-seq), and quantified cuticular hydrocarbon (CHC) through gas chromatography-mass spectrometry in comparative studies of the integument (cuticle plus its underlying epidermis) of two eusocial and a solitary bee species. In addition, we used transmission electron microscopy (TEM) for studying the developing cuticle of these and other three bee species also differing in life style. We found 13,200, 55,209 and 30,161 transcript types in the integument of the eusocial Apis mellifera and Frieseomelitta varia, and the solitary Centris analis, respectively. In general, structural cuticle proteins and chitin-related genes were upregulated in pharate-adults and newly-emerged bees whereas transcripts for odorant binding proteins, cytochrome P450 and antioxidant proteins were overrepresented in foragers. Consistent with our hypothesis, a distance correlation analysis based on the differentially expressed genes suggested delayed cuticle maturation in A. mellifera in comparison to the solitary bee. However, this was not confirmed in the comparison with F. varia. The expression profiles of 27 of 119 genes displaying functional attributes related to cuticle formation/differentiation were positively correlated between A. mellifera and F. varia, and negatively or non-correlated with C. analis, suggesting roles in cuticular maturation heterochrony. However, we also found transcript profiles positively correlated between each one of the eusocial species and C. analis. Gene co-expression networks greatly differed between the bee species, but we identified common gene interactions exclusively between the eusocial species. Except for F. varia, the TEM analysis is consistent with cuticle development timing adapted to the social or solitary life style. In support to our hypothesis, the absolute quantities of n-alkanes and unsaturated CHCs were significantly higher in foragers than in the earlier developmental phases of the eusocial bees, but did not discriminate newly-emerged from foragers in C. analis. By highlighting differences in integument gene expression, cuticle ultrastructure, and CHC profiles between eusocial and solitary bees, our data provided insights into the process of heterochronic cuticle maturation associated to the way of life.

Draft genome of Thermomonospora sp. CIT 1 (Thermomonosporaceae) and in silico evidence of its functional role in filter cake biomass deconstruction

The filter cake from sugar cane processing is rich in organic matter and nutrients, which favors the proliferation of microorganisms with potential to deconstruct plant biomass. From the metagenomic data of this material, we assembled a draft genome that was phylogenetically related to Thermomonospora curvata DSM 43183, which shows the functional and ecological importance of this bacterium in the filter cake. Thermomonospora is a gram-positive bacterium that produces cellulases in compost, and it can survive temperatures of 60 ºC. We identified a complete set of biomass depolymerizing enzymes in the draft genome of Thermomonospora sp. CIT 1, such as α-amylase, catalase-peroxidases, β-mannanase, and arabinanase, demonstrating the potential of this bacterium to deconstruct the components of starch, lignin, and hemicellulose. In addition, the draft genome of Thermomonospora sp. CIT 1 contains 18 genes that do not share identity with five other species of Thermomonospora, suggesting that this bacterium has different genetic characteristics than those present in genomes reported so far for this genus. These findings add a new dimension to the current understanding of the functional profile of this microorganism that inhabits agro-industrial waste, which may boost new gene discoveries and be of importance for application in the production of bioethanol.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

Comparative genomic analysis of Genlisea (corkscrew plants-Lentibulariaceae) chloroplast genomes reveals an increasing loss of the ndh genes

In the carnivorous plant family Lentibulariaceae, all three genome compartments (nuclear, chloroplast, and mitochondria) have some of the highest rates of nucleotide substitutions across angiosperms. While the genera Genlisea and Utricularia have the smallest known flowering plant nuclear genomes, the chloroplast genomes (cpDNA) are mostly structurally conserved except for deletion and/or pseudogenization of the NAD(P)H-dehydrogenase complex (ndh) genes known to be involved in stress conditions of low light or CO2 concentrations. In order to determine how the cpDNA are changing, and to better understand the evolutionary history within the Genlisea genus, we sequenced, assembled and analyzed complete cpDNA from six species (G. aurea, G. filiformis, G. pygmaea, G. repens, G. tuberosa and G. violacea) together with the publicly available G. margaretae cpDNA. In general, the cpDNA structure among the analyzed Genlisea species is highly similar. However, we found that the plastidial ndh genes underwent a progressive process of degradation similar to the other terrestrial Lentibulariaceae cpDNA analyzed to date, but in contrast to the aquatic species. Contrary to current thinking that the terrestrial environment is a more stressful environment and thus requiring the ndh genes, we provide evidence that in the Lentibulariaceae the terrestrial forms have progressive loss while the aquatic forms have the eleven plastidial ndh genes intact. Therefore, the Lentibulariaceae system provides an important opportunity to understand the evolutionary forces that govern the transition to an aquatic environment and may provide insight into how plants manage water stress at a genome scale.

The mitochondrial genome of the terrestrial carnivorous plant Utricularia reniformis (Lentibulariaceae): Structure, comparative analysis and evolutionary landmarks

The carnivorous plants of the family Lentibulariaceae have attained recent attention not only because of their interesting lifestyle, but also because of their dynamic nuclear genome size. Lentibulariaceae genomes span an order of magnitude and include species with the smallest genomes in angiosperms, making them a powerful system to study the mechanisms of genome expansion and contraction. However, little is known about mitochondrial DNA (mtDNA) sequences of this family, and the evolutionary forces that shape this organellar genome. Here we report the sequencing and assembly of the complete mtDNA from the endemic terrestrial Brazilian species Utricularia reniformis. The 857,234bp master circle mitochondrial genome encodes 70 transcriptionaly active genes (42 protein-coding, 25 tRNAs and 3 rRNAs), covering up to 7% of the mtDNA. A ltrA-like protein related to splicing and mobility and a LAGLIDADG homing endonuclease have been identified in intronic regions, suggesting particular mechanisms of genome maintenance. RNA-seq analysis identified properties with putative diverse and important roles in genome regulation and evolution: 1) 672kbp (78%) of the mtDNA is covered by full-length reads; 2) most of the 243kbp intergenic regions exhibit transcripts; and 3) at least 69 novel RNA editing sites in the protein-coding genes. Additional genomic features are hypothetical ORFs (48%), chloroplast insertions, including truncated plastid genes that have been lost from the chloroplast DNA (5%), repeats (5%), relics of transposable elements mostly related to LTR retrotransposons (5%), and truncated mitovirus sequences (0.4%). Phylogenetic analysis based on 32 different Lamiales mitochondrial genomes corroborate that Lentibulariaceae is a monophyletic group. In summary, the U. reniformis mtDNA represents the eighth largest plant mtDNA described to date, shedding light on the genomic trends and evolutionary characteristics and phylogenetic history of the family Lentibulariaceae.

De novo transcriptome assembly of sugarcane leaves submitted to prolonged water-deficit stress

Sugarcane production is strongly influenced by drought, which is a limiting factor for agricultural productivity in the world. In this study, the gene expression profiles obtained by de novo assembly of the leaf transcriptome of two sugarcane cultivars that differ in their physiological response to water deficit were evaluated by the RNA-Seq method: drought-tolerant cultivar (SP81-3250) and drought-sensitive cultivar (RB855453). For this purpose, plants were grown in a greenhouse for 60 days and were then submitted to three treatments: control (-0.01 to -0.015 MPa), moderate water deficit (-0.05 to -0.055 MPa), and severe water deficit (-0.075 to -0.08 MPa). The plants were evaluated 30, 60, and 90 days after the beginning of treatment. Sequencing on an Illumina platform (RNA-Seq) generated more than one billion sequences, resulting in 177,509 and 185,153 transcripts for the tolerant and sensitive cultivar, respectively. These transcripts were aligned with sequences from Saccharum spp, Sorghum bicolor, Miscanthus giganteus, and Arabidopsis thaliana available in public databases. The differentially expressed genes detected during the prolonged period of water deficit permit to increase our understanding of the molecular patterns involved in the physiological response of the two cultivars. The tolerant cultivar differentially expressed a larger number of genes at 90 days, while in the sensitive cultivar the number of differentially expressed genes was higher in 30 days. Both cultivars perceived the lack of water, but the tolerant cultivar responded more slowly than the sensitive cultivar. The latter requires rapid activation of different water-deficit stress response mechanisms for its survival. This rapid activation of metabolic pathways in response to water stress does not appear to be the key mechanism of drought tolerance in sugarcane. There is still much to clarify on the molecular and physiological pattern of plants in response to drought.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

Transcriptional profiling reveals intrinsic mRNA alterations in multipotent mesenchymal stromal cells isolated from bone marrow of newly-diagnosed type 1 diabetes patients

Background

Bone marrow multipotent mesenchymal stromal cells (MSCs) are a diverse subset of precursors that contribute to the homeostasis of the hematopoietic niche. MSCs can be isolated and expanded in vitro and have unique immunomodulatory and regenerative properties that make them attractive for the treatment of autoimmune diseases, including type 1 diabetes (T1D). Whether autologous or allogeneic MSCs are more suitable for therapeutic purposes has not yet been established. While autologous MSCs may present abnormal function, allogeneic cells may be recognized and rejected by the host immune system. Thus, studies that investigate biological characteristics of MSCs isolated from T1D patients are essential to guide future clinical applications.

Methods

Bone marrow-derived MSCs from recently diagnosed type 1 diabetes patients (T1D-MSCs) were compared with those from healthy individuals (C-MSCs) for morphological and immunophenotypic characteristics and for differentiation potential. Bioinformatics approaches allowed us to match absolute and differential gene expression of several adhesion molecules, immune mediators, growth factors, and their receptors involved with hematopoietic support and immunomodulatory properties of MSCs. Finally, the differentially expressed genes were collated for functional pathway enrichment analysis.

Results

T1D-MSCs and C-MSCs were similar for morphology, immunophenotype, and differentiation potential. Our absolute gene expression results supported previous literature reports, while also detecting new potential molecules related to bone marrow-derived MSC functions. T1D-MSCs showed intrinsic abnormalities in mRNA expression, including the immunomodulatory molecules VCAM-1, CXCL12, HGF, and CCL2. Pathway analyses revealed activation of sympathetic nervous system and JAK STAT signaling in T1D-MSCs.

Conclusions

Collectively, our results indicate that MSCs isolated from T1D patients present intrinsic transcriptional alterations that may affect their therapeutic potential. However, the implications of these abnormalities in T1D development as well as in the therapeutic efficacy of autologous MSCs require further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13287-016-0351-y) contains supplementary material, which is available to authorized users.

MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.)

Queen and worker honeybees differ profoundly in reproductive capacity. The queen of this complex society, with 200 highly active ovarioles in each ovary, is the fertile caste, whereas the workers have approximately 20 ovarioles as a result of receiving a different diet during larval development. In a regular queenright colony, the workers have inactive ovaries and do not reproduce. However, if the queen is sensed to be absent, some of the workers activate their ovaries, producing viable haploid eggs that develop into males. Here, a deep-sequenced ovary transcriptome library of reproductive workers was used as supporting data to assess the dynamic expression of the regulatory molecules and microRNAs (miRNAs) of reproductive and nonreproductive honeybee females. In this library, most of the differentially expressed miRNAs are related to ovary physiology or oogenesis. When we quantified the dynamic expression of 19 miRNAs in the active and inactive worker ovaries and compared their expression in the ovaries of virgin and mated queens, we noted that some miRNAs (miR-1, miR-31a, miR-13b, miR-125, let-7 RNA, miR-100, miR-276, miR-12, miR-263a, miR-306, miR-317, miR-92a and miR-9a) could be used to identify reproductive and nonreproductive statuses independent of caste. Furthermore, integrative gene networks suggested that some candidate miRNAs function in the process of ovary activation in worker bees.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

Bone Marrow Mesenchymal Stromal Cells Isolated from Multiple Sclerosis Patients have Distinct Gene Expression Profile and Decreased Suppressive Function Compared with Healthy Counterparts

Multiple sclerosis (MS) is a chronic inflammatory autoimmune disease of the central nervous system, due to an immune reaction against myelin proteins. Multipotent mesenchymal stromal cells (MSCs) present immunosuppressive effects and have been used for the treatment of autoimmune diseases. In our study, gene expression profile and in vitro immunomodulatory function tests were used to compare bone marrow-derived MSCs obtained from MS patients, at pre- and postautologous hematopoietic stem cell transplantation (AHSCT) with those from healthy donors. Patient MSCs comparatively exhibited i) senescence in culture; ii) similar osteogenic and adipogenic differentiation potential; iii) decreased expression of CD105, CD73, CD44, and HLA-A/B/C molecules; iv) distinct transcription at pre-AHSCT compared with control MSCs, yielding 618 differentially expressed genes, including the downregulation of TGFB1 and HGF genes and modulation of the FGF and HGF signaling pathways; v) reduced antiproliferative effects when pre-AHSCT MSCs were cocultured with allogeneic T-lymphocytes; vi) decreased secretion of IL-10 and TGF-β in supernatants of both cocultures (pre- and post-AHSCT MSCs); and vii) similar percentages of regulatory cells recovered after MSC cocultures. The transcriptional profile of patient MSCs isolated 6 months posttransplantation was closer to pre-AHSCT samples than from healthy MSCs. Considering that patient MSCs exhibited phenotypic changes, distinct transcriptional profile and functional defects implicated in MSC immunomodulatory and immunosuppressive activity, we suggest that further MS clinical studies should be conducted using allogeneic bone marrow MSCs derived from healthy donors. We also demonstrated that treatment of MS patients with AHSCT does not reverse the transcriptional and functional alterations observed in patient MSCs.

Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera)

Major developmental transitions in multicellular organisms are driven by steroid hormones. In insects, these, together with juvenile hormone (JH), control development, metamorphosis, reproduction and aging, and are also suggested to play an important role in caste differentiation of social insects. Here, we aimed to determine how EcR transcription and ecdysteroid titers are related during honeybee postembryonic development and what may actually be the role of EcR in caste development of this social insect. In addition, we expected that knocking-down EcR gene expression would give us information on the participation of the respective protein in regulating downstream targets of EcR. We found that in Apis mellifera females, EcR-A is the predominantly expressed variant in postembryonic development, while EcR-B transcript levels are higher in embryos, indicating an early developmental switch in EcR function. During larval and pupal stages, EcR-B expression levels are very low, while EcR-A transcripts are more variable and abundant in workers compared to queens. Strikingly, these transcript levels are opposite to the ecdysteroid titer profile. 20-hydroxyecdysone (20E) application experiments revealed that low 20E levels induce EcR expression during development, whereas high ecdysteroid titers seem to be repressive. By means of RNAi-mediated knockdown (KD) of both EcR transcript variants we detected the differential expression of 234 poly-A⁺ transcripts encoding genes such as CYPs, MRJPs and certain hormone response genes (Kr-h1 and ftz-f1). EcR-KD also promoted the differential expression of 70 miRNAs, including highly conserved ones (e.g., miR-133 and miR-375), as well honeybee-specific ones (e.g., miR-3745 and miR-3761). Our results put in evidence a broad spectrum of EcR-controlled gene expression during postembryonic development of honeybees, revealing new facets of EcR biology in this social insect.

Altered expression of degranulation-related genes in CD8+ T cells in human T lymphotropic virus type I infection

Human T lymphotropic virus type I (HTLV-1) is the etiological agent of HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). CD8+ T cells may contribute to the protection or development of HAM/TSP. In this study we used SAGE methodology to screen for differentially expressed genes in CD8+ T cells isolated from HTLV-1 asymptomatic carriers (HAC) and from HAM/TSP patients to identify genes involved in HAM/TSP development. SAGE analysis was conducted by pooling samples according to clinical status. The comparison of gene expression profiles between HAC and HAM/TSP libraries identified 285 differentially expressed tags. We focus on cytotoxicity and cytokine-related genes due to their potential biological role in HTLV-1 infection. Our results showed that patients with HAM/TSP have high expression levels of degranulation-related genes, namely GZMH and PRF1, and of the cytoskeletal adaptor PXN. We found that GZMB and ZAP70 were overexpressed in HTLV-infected patients compared to the noninfected group. We also detected that CCL5 was higher in the HAM/TSP group compared to the HAC and CT groups. Our findings showed that CD8+ T cells of HAM/TSP patients have an inflammatory and active profile. PXN and ZAP70 overexpression in HTLV-1-infected patients was described for the first time here and reinforces this concept. However, although active and abundant, CD8+ T cells are not able to completely eliminate infected cells and prevent the development of HAM/TSP and, moreover, these cells might contribute to the pathogenesis of the disease by migrating to the central nervous system (CNS). These results should be further tested with biological functional assays to increase our understanding on the role of these molecules in the development of HTLV-1-related diseases.

SAGE analysis highlights the putative role of underexpression of ribosomal proteins in GH-secreting pituitary adenomas

BACKGROUND

Although the molecular pathogenesis of pituitary adenomas has been assessed by several different techniques, it still remains partially unclear. Ribosomal proteins (RPs) have been recently related to human tumorigenesis, but they have not yet been evaluated in pituitary tumorigenesis.

OBJECTIVE

The aim of this study was to introduce serial analysis of gene expression (SAGE), a high-throughput method, in pituitary research in order to compare differential gene expression.

METHODS

Two SAGE cDNA libraries were constructed, one using a pool of mRNA obtained from five GH-secreting pituitary tumors and another from three normal pituitaries. Genes differentially expressed between the libraries were further validated by real-time PCR in 22 GH-secreting pituitary tumors and in 15 normal pituitaries.

RESULTS

Computer-generated genomic analysis tools identified 13,722 and 14,993 exclusive genes in normal and adenoma libraries respectively. Both shared 6497 genes, 2188 were underexpressed and 4309 overexpressed in tumoral library. In adenoma library, 33 genes encoding RPs were underexpressed. Among these, RPSA, RPS3, RPS14, and RPS29 were validated by real-time PCR.

CONCLUSION

We report the first SAGE library from normal pituitary tissue and GH-secreting pituitary tumor, which provide quantitative assessment of cellular transcriptome. We also validated some downregulated genes encoding RPs. Altogether, the present data suggest that the underexpression of the studied RP genes possibly collaborates directly or indirectly with other genes to modify cell cycle arrest, DNA repair, and apoptosis, leading to an environment that might have a putative role in the tumorigenesis, introducing new perspectives for further studies on molecular genesis of somatotrophinomas.

High-throughput sequencing of black pepper root transcriptome

BACKGROUND

Black pepper (Piper nigrum L.) is one of the most popular spices in the world. It is used in cooking and the preservation of food and even has medicinal properties. Losses in production from disease are a major limitation in the culture of this crop. The major diseases are root rot and foot rot, which are results of root infection by Fusarium solani and Phytophtora capsici, respectively. Understanding the molecular interaction between the pathogens and the host's root region is important for obtaining resistant cultivars by biotechnological breeding. Genetic and molecular data for this species, though, are limited. In this paper, RNA-Seq technology has been employed, for the first time, to describe the root transcriptome of black pepper.

RESULTS

The root transcriptome of black pepper was sequenced by the NGS SOLiD platform and assembled using the multiple-k method. Blast2Go and orthoMCL methods were used to annotate 10338 unigenes. The 4472 predicted proteins showed about 52% homology with the Arabidopsis proteome. Two root proteomes identified 615 proteins, which seem to define the plant's root pattern. Simple-sequence repeats were identified that may be useful in studies of genetic diversity and may have applications in biotechnology and ecology.

CONCLUSIONS

This dataset of 10338 unigenes is crucially important for the biotechnological breeding of black pepper and the ecogenomics of the Magnoliids, a major group of basal angiosperms.

Distinct patterns of somatic alterations in a lymphoblastoid and a tumor genome derived from the same individual

Although patterns of somatic alterations have been reported for tumor genomes, little is known on how they compare with alterations present in non-tumor genomes. A comparison of the two would be crucial to better characterize the genetic alterations driving tumorigenesis. We sequenced the genomes of a lymphoblastoid (HCC1954BL) and a breast tumor (HCC1954) cell line derived from the same patient and compared the somatic alterations present in both. The lymphoblastoid genome presents a comparable number and similar spectrum of nucleotide substitutions to that found in the tumor genome. However, a significant difference in the ratio of non-synonymous to synonymous substitutions was observed between both genomes (P = 0.031). Protein-protein interaction analysis revealed that mutations in the tumor genome preferentially affect hub-genes (P = 0.0017) and are co-selected to present synergistic functions (P < 0.0001). KEGG analysis showed that in the tumor genome most mutated genes were organized into signaling pathways related to tumorigenesis. No such organization or synergy was observed in the lymphoblastoid genome. Our results indicate that endogenous mutagens and replication errors can generate the overall number of mutations required to drive tumorigenesis and that it is the combination rather than the frequency of mutations that is crucial to complete tumorigenic transformation.

Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach

Background

While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia.

Methodology/Principal Findings

A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue.

Conclusions/Significance

This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.

Genomics and proteomics approaches to the study of cancer-stroma interactions

Background

The development and progression of cancer depend on its genetic characteristics as well as on the interactions with its microenvironment. Understanding these interactions may contribute to diagnostic and prognostic evaluations and to the development of new cancer therapies. Aiming to investigate potential mechanisms by which the tumor microenvironment might contribute to a cancer phenotype, we evaluated soluble paracrine factors produced by stromal and neoplastic cells which may influence proliferation and gene and protein expression.

Methods

The study was carried out on the epithelial cancer cell line (Hep-2) and fibroblasts isolated from a primary oral cancer. We combined a conditioned-medium technique with subtraction hybridization approach, quantitative PCR and proteomics, in order to evaluate gene and protein expression influenced by soluble paracrine factors produced by stromal and neoplastic cells.

Results

We observed that conditioned medium from fibroblast cultures (FCM) inhibited proliferation and induced apoptosis in Hep-2 cells. In neoplastic cells, 41 genes and 5 proteins exhibited changes in expression levels in response to FCM and, in fibroblasts, 17 genes and 2 proteins showed down-regulation in response to conditioned medium from Hep-2 cells (HCM). Nine genes were selected and the expression results of 6 down-regulated genes (ARID4A, CALR, GNB2L1, RNF10, SQSTM1, USP9X) were validated by real time PCR.

Conclusions

A significant and common denominator in the results was the potential induction of signaling changes associated with immune or inflammatory response in the absence of a specific protein.

A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background

High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis.

Results

This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system.

Conclusion

These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at . S3T source code and datasets can also be downloaded from the aforementioned website.

miRNApath: a database of miRNAs, target genes and metabolic pathways

MicroRNAs (miRNAs) are small non-coding RNAs that regulate target gene expression and hence play important roles in metabolic pathways. Recent studies have evidenced the interrelation of miRNAs with cell proliferation, differentiation, development, and diseases. Since they are involved in gene regulation, they are intrinsically related to metabolic pathways. This leads to questions that are particularly interesting for investigating medical and laboratorial applications. We developed an miRNApath online database that uses miRNA target genes to link miRNAs to metabolic pathways. Currently, databases about miRNA target genes (DIANA miRGen), genomic maps (miRNAMap) and sequences (miRBase) do not provide such correlations. Additionally, miRNApath offers five search services and a download area. For each search, there is a specific type of input, which can be a list of target genes, miRNAs, or metabolic pathways, which results in different views, depending upon the input data, concerning relationships between the target genes, miRNAs and metabolic pathways. There are also internal links that lead to a deeper analysis and cross-links to other databases with more detailed information. miRNApath is being continually updated and is available at http://lgmb.fmrp.usp.br/mirnapath.

Proteomic and SAGE profiling of murine melanoma progression indicates the reduction of proteins responsible for ROS degradation

Using 2-DE of total cell protein extracts, we compared soluble proteins from murine melanoma lines Tm1 and Tm5 with proteins from the nontumoral cell melan-a from which they were derived. Seventy-one of the 452 spots (average) detected with CBB were differentially accumulated, i.e., increased or decreased twofold. Forty-four spots were identified by PMF/MALDI-TOF, 15 with increased and 29 with decreased protein levels. SAGE showed that 17/34 (50%) of the differentially accumulated proteins, pI range 4-7, presented similar differences at the mRNA level. Major reductions in protein were observed in tumor cells of proteins that degrade reactive oxygen species (ROS). Decreases of > or = twofold in GST, superoxide dismutase, aldehyde dehydrogenase, thioredoxin, peroxiredoxin 2, and peroxiredoxin 6 protein were observed. SAGE indicated the reduction of other proteins involved in ROS degradation. As expected, the accumulation of exogenous peroxides was significantly higher in the tumor cells while the levels of glutathionylation were two times lower in the tumor cells compared to melan-a. The differential accumulation of proteins involved in oncogene/tumor suppressor pathways was observed. Melanoma cells can favor survival pathways activated by ROS by inhibiting p53 pathways and activation of Ras and c-myc pathways.

Identification and complete sequencing of novel human transcripts through the use of mouse orthologs and testis cDNA sequences

The correct identification of all human genes, and their derived transcripts, has not yet been achieved, and it remains one of the major aims of the worldwide genomics community. Computational programs suggest the existence of 30,000 to 40,000 human genes. However, definitive gene identification can only be achieved by experimental approaches. We used two distinct methodologies, one based on the alignment of mouse orthologous sequences to the human genome, and another based on the construction of a high-quality human testis cDNA library, in an attempt to identify new human transcripts within the human genome sequence. We generated 47 complete human transcript sequences, comprising 27 unannotated and 20 annotated sequences. Eight of these transcripts are variants of previously known genes. These transcripts were characterized according to size, number of exons, and chromosomal localization, and a search for protein domains was undertaken based on their putative open reading frames. In silico expression analysis suggests that some of these transcripts are expressed at low levels and in a restricted set of tissues.

A transcript finishing initiative for closing gaps in the human transcriptome

We report the results of a transcript finishing initiative, undertaken for the purpose of identifying and characterizing novel human transcripts, in which RT-PCR was used to bridge gaps between paired EST clusters, mapped against the genomic sequence. Each pair of EST clusters selected for experimental validation was designated a transcript finishing unit (TFU). A total of 489 TFUs were selected for validation, and an overall efficiency of 43.1% was achieved. We generated a total of 59,975 bp of transcribed sequences organized into 432 exons, contributing to the definition of the structure of 211 human transcripts. The structure of several transcripts reported here was confirmed during the course of this project, through the generation of their corresponding full-length cDNA sequences. Nevertheless, for 21% of the validated TFUs, a full-length cDNA sequence is not yet available in public databases, and the structure of 69.2% of these TFUs was not correctly predicted by computer programs. The TF strategy provides a significant contribution to the definition of the complete catalog of human genes and transcripts, because it appears to be particularly useful for identification of low abundance transcripts expressed in a restricted set of tissues as well as for the delineation of gene boundaries and alternatively spliced isoforms.