Itraconazole is currently used to treat paracoccidioidomycosis. The mechanism of action of azoles has been elucidated in some fungi, although little is known regarding its mechanism of action in Paracoccidioides spp. The present work focused on identification of regulated transcripts using representational difference analysis of Paracoccidioides spp. yeast cells treated with itraconazole for 1 and 2 h.
Paracoccidioides Pb01 genes up-regulated by itraconazole included genes involved in cellular transport, metabolism/energy, transcription, cell rescue, defense and virulence. ERG11, ERG6, ERG3, ERG5 and ERG25 were up-regulated at multiple time points. In vivo infection experiments in mice corroborated the in vitro results. Ergosterol levels and distribution were evaluated in Paracoccidioides Pb18 yeast cells, and the results demonstrate that both factors were changed in the fungus treated with itraconazole.
To our knowledge, this is the first transcriptional analysis of Paracoccidioides spp. exposed to a triazole drug. Here acetyl seems to be intensively produced from different metabolic pathways to produce ergosterol by the action of ergosterol synthesis related enzymes, which were also affected in other fungi. Among the genes affected, we identified genes in common with other fungi, as well as genes unique to Paracoccidioides Pb01. Those genes could be considered target to new drugs. Voltage-gated Ca2+ alpha subunit (CAV), Tetracycline resistance protein (TETA) and Hemolisyn-iii channel protein (HLYiii) were found only here and a probably involvement with resistence to itraconazole could be investigated in the future. However our findings do not permit inference to current clinical practice.
Paracoccidioides spp.; Transcriptional response; Itraconazole; Ergosterol
Human RUNX1 gene is one of the most frequent target for chromosomal translocations associated with acute myeloid leukemia (AML) and acute lymphoid leukemia (ALL). The highest prevalence in AML is noted with (8; 21) translocation; which represents 12 to 15% of all AML cases. Interestingly, all the breakpoints mapped to date in t(8;21) are clustered in intron 5 of the RUNX1 gene and intron 1 of the ETO gene. No homologous sequences have been found at the recombination regions; but DNase I hypersensitive sites (DHS) have been mapped to the areas of the genes involved in t(8;21). Presence of DHS sites is commonly associated with regulatory elements such as promoters, enhancers and silencers, among others.
In this study we used a combination of comparative genomics, cloning and transfection assays to evaluate potential regulatory elements located in intron 5 of the RUNX1 gene. Our genomic analysis identified nine conserved non-coding sequences that are evolutionarily conserved among rat, mouse and human. We cloned two of these regions in pGL-3 Promoter plasmid in order to analyze their transcriptional regulatory activity. Our results demonstrate that the identified regions can indeed regulate transcription of a reporter gene in a distance and position independent manner; moreover, their transcriptional effect is cell type specific.
We have identified nine conserved non coding sequence that are harbored in intron 5 of the RUNX1 gene. We have also demonstrated that two of these regions can regulate transcriptional activity in vitro. Taken together our results suggest that intron 5 of the RUNX1 gene contains multiple potential cis-regulatory elements.
RUNX1 gene; AML; (8; 21) translocation; Non coding sequences; Enhancer modules; Chromatin structure
Yellowtail (Seriola quinqueradiata) are an economically important species in Japan. However, there are currently no methods for captive breeding and early rearing for yellowtail. Thus, the commercial cultivation of this species is reliant upon the capture of wild immature fish. Given this, there is a need to develop captive breeding techniques to reduce pressure on wild stocks and facilitate the sustainable development of yellowtail aquaculture. We constructed a whole genome radiation hybrid (RH) panel for yellowtail gene mapping and developed a framework physical map using a nanofluidic dynamic array to use SNPs (single nucleotide polymorphisms) in ESTs (expressed sequence tags) for the DNA-assisted breeding of yellowtail.
Clonal RH cell lines were obtained after ionizing radiation; specifically, 78, 64, 129, 55, 42, and 53 clones were isolated after treatment with 3,000, 4,000, 5,000, 6,000, 8,000, or 10,000 rads, respectively. A total of 421 hybrid cell lines were obtained by fusion with mouse B78 cells. Ninety-four microsatellite markers used in the genetic linkage map were genotyped using the 421 hybrid cell lines. Based upon marker retention and genome coverage, we selected 93 hybrid cell lines to form an RH panel. Importantly, we performed the first genotyping of yellowtail markers in an RH panel using a nanofluidic dynamic array (Fluidigm, CA, USA). Then, 580 markers containing ESTs and SNPs were mapped in the first yellowtail RH map.
We successfully developed a yellowtail RH panel to facilitate the localization of markers. Using this, a framework RH map was constructed with 580 markers. This high-density physical map will serve as a useful tool for the identification of genes related to important breeding traits using genetic structural information, such as conserved synteny. Moreover, in a comparison of 30 sequences in the RH group 1 (SQ1), yellowtail appeared to be evolutionarily closer to medaka and the green-spotted pufferfish than to zebrafish. We suggest that synteny analysis may be potentially useful as a tool to investigate chromosomal evolution by comparison with model fish.
Yellowtail; Radiation hybrid (RH) panel; RH map; Synteny
The basidiomycete Moniliophthora roreri is the causal agent of Frosty pod rot (FPR) disease of cacao (Theobroma cacao), the source of chocolate, and FPR is one of the most destructive diseases of this important perennial crop in the Americas. This hemibiotroph infects only cacao pods and has an extended biotrophic phase lasting up to sixty days, culminating in plant necrosis and sporulation of the fungus without the formation of a basidiocarp.
We sequenced and assembled 52.3 Mb into 3,298 contigs that represent the M. roreri genome. Of the 17,920 predicted open reading frames (OFRs), 13,760 were validated by RNA-Seq. Using read count data from RNA sequencing of cacao pods at 30 and 60 days post infection, differential gene expression was estimated for the biotrophic and necrotrophic phases of this plant-pathogen interaction. The sequencing data were used to develop a genome based secretome for the infected pods. Of the 1,535 genes encoding putative secreted proteins, 1,355 were expressed in the biotrophic and necrotrophic phases. Analysis of the data revealed secretome gene expression that correlated with infection and intercellular growth in the biotrophic phase and invasive growth and plant cellular death in the necrotrophic phase.
Genome sequencing and RNA-Seq was used to determine and validate the Moniliophthora roreri genome and secretome. High sequence identity between Moniliophthora roreri genes and Moniliophthora perniciosa genes supports the taxonomic relationship with Moniliophthora perniciosa and the relatedness of this fungus to other basidiomycetes. Analysis of RNA-Seq data from infected plant tissues revealed differentially expressed genes in the biotrophic and necrotrophic phases. The secreted protein genes that were upregulated in the biotrophic phase are primarily associated with breakdown of the intercellular matrix and modification of the fungal mycelia, possibly to mask the fungus from plant defenses. Based on the transcriptome data, the upregulated secreted proteins in the necrotrophic phase are hypothesized to be actively attacking the plant cell walls and plant cellular components resulting in necrosis. These genes are being used to develop a new understanding of how this disease interaction progresses and to identify potential targets to reduce the impact of this devastating disease.
The Type VI secretion system (T6SS) has been identified in several different bacteria, including the plant pathogenPantoea ananatis. Previous in silico analyses described three different T6SS loci present in the pathogenic strain of P. ananatis LMG 20103. This initial investigation has been extended to include an additional seven sequenced strains of P. ananatis together with 39 strains from different ecological niches. Comparative and phylogenetic analyses were used to investigate the distribution, evolution, intra-strain variability and operon structure of the T6SS in the sequenced strains.
Three different T6SS loci were identified in P. ananatis strain LMG 20103 and designated PA T6SS 1-3. PA T6SS-1 was present in all sequenced strains of P. ananatis and in all 39 additional strains examined in this study. In addition, PA T6SS-1 included all 13 core T6SS genes required for synthesis of a functional T6SS. The plasmid-borne PA T6SS-2 also included all 13 core T6SS genes but was restricted to only 33% (15/46) of the strains examined. In addition, PA T6SS-2 was restricted to strains of P. ananatis isolated from symptomatic plant material. This finding raises the possibility of an association between PA T6SS-2 and either pathogenicity or host specificity. The third cluster PA T6SS-3 was present in all strains analyzed in this study but lacked 11 of the 13 core T6SS genes suggesting it may not encoded a functional T6SS. Inter-strain variability was also associated with hcp and vgrG islands, which are associated with the T6SS and encode a variable number of proteins usually of unknown function. These proteins may play a role in the fitness of different strains in a variety of ecological niches or as candidate T6SS effectors. Phylogenetic analysis indicated that PA T6SS-1 and PA T6SS-2 are evolutionarily distinct.
Our analysis indicates that the three T6SSs of P. ananatis appear to have been independently acquired and may play different roles relating to pathogenicity, host range determination and/or niche adaptation. Future work will be directed toward understanding the roles that these T6SSs play in the biology of P. ananatis.
Type VI secretion system; T6SS; Pantoea ananatis; Phytopathogen
Gene duplication is a source of evolutionary innovation and can contribute to the divergence of lineages; however, the relative importance of this process remains to be determined. The explosive divergence of the African cichlid adaptive radiations provides both a model for studying the general role of gene duplication in the divergence of lineages and also an exciting foray into the identification of genomic features that underlie the dramatic phenotypic and ecological diversification in this particular lineage. We present the first genome-wide study of gene duplication in African cichlid fishes, identifying gene duplicates in three species belonging to the Lake Malawi adaptive radiation (Metriaclima estherae, Protomelas similis, Rhamphochromis “chilingali”) and one closely related species from a non-radiated riverine lineage (Astatotilapia tweddlei).
Using Astatotilapia burtoni as reference, microarray comparative genomic hybridization analysis of 5689 genes reveals 134 duplicated genes among the four cichlid species tested. Between 51 and 55 genes were identified as duplicated in each of the three species from the Lake Malawi radiation, representing a 38%–49% increase in number of duplicated genes relative to the non-radiated lineage (37 genes). Duplicated genes include several that are involved in immune response, ATP metabolism and detoxification.
These results contribute to our understanding of the abundance and type of gene duplicates present in cichlid fish lineages. The duplicated genes identified in this study provide candidates for the analysis of functional relevance with regard to phenotype and divergence. Comparative sequence analysis of gene duplicates can address the role of positive selection and adaptive evolution by gene duplication, while further study across the phylogenetic range of cichlid radiations (and more generally in other adaptive radiations) will determine whether the patterns of gene duplication seen in this study consistently accompany rapid radiation.
Phylogenetic trees are widely used for genetic and evolutionary studies in various organisms. Advanced sequencing technology has dramatically enriched data available for constructing phylogenetic trees based on single nucleotide polymorphisms (SNPs). However, massive SNP data makes it difficult to perform reliable analysis, and there has been no ready-to-use pipeline to generate phylogenetic trees from these data.
We developed a new pipeline, SNPhylo, to construct phylogenetic trees based on large SNP datasets. The pipeline may enable users to construct a phylogenetic tree from three representative SNP data file formats. In addition, in order to increase reliability of a tree, the pipeline has steps such as removing low quality data and considering linkage disequilibrium. A maximum likelihood method for the inference of phylogeny is also adopted in generation of a tree in our pipeline.
Using SNPhylo, users can easily produce a reliable phylogenetic tree from a large SNP data file. Thus, this pipeline can help a researcher focus more on interpretation of the results of analysis of voluminous data sets, rather than manipulations necessary to accomplish the analysis.
Polymorphisms; Linkage disequilibrium; Maximum likelihood
The unfolded protein response (UPR) is a network of intracellular signaling pathways that supports the ability of the secretory pathway to maintain a balance between the load of proteins entering the endoplasmic reticulum (ER) and the protein folding capacity of the ER lumen. Current evidence indicates that several pathogenic fungi rely heavily on this pathway for virulence, but there is limited understanding of the mechanisms involved. The best known functional output of the UPR is transcriptional upregulation of mRNAs involved in ER homeostasis. However, this does not take into account mechanisms of translational regulation that involve differential loading of ribosomes onto mRNAs. In this study, a global analysis of transcript-specific translational regulation was performed in the pathogenic mold Aspergillus fumigatus to determine the nature and scope of the translational response to ER stress.
ER stress was induced by treating the fungus with dithiothreitol, tunicamycin, or a thermal up-shift. The mRNAs were then fractionated on the basis of ribosome occupancy into an under-translated pool (U) and a well-translated pool (W). The mRNAs were used to interrogate microarrays and the ratio of the hybridization signal (W/U) was used as an indicator of the relative translational efficiency of a mRNA under each condition. The largest category of translationally upregulated mRNAs during ER stress encoded proteins involved in translation. Components of the ergosterol and GPI anchor biosynthetic pathways also showed increased polysome association, suggesting an important role for translational regulation in membrane and cell wall homeostasis. ER stress induced limited remodeling of the secretory pathway translatome. However, a select group of transcription factors was translationally upregulated, providing a link to subsequent modification of the transcriptome. Finally, we provide evidence that one component of the ER stress translatome is a novel mRNA isoform from the yvc1 gene that is induced by ER stress in a UPR-dependent manner.
Together, these findings define a core set of mRNAs subject to translational control during the adaptive response to acute ER stress in A. fumigatus and reveal a remarkable breadth of functions that are needed to resolve ER stress in this organism.
Aspergillus fumigatus; UPR; Unfolded protein response; ER stress; Translational regulation; Polysome profiling; Yvc1
Synthetic lethality is an appealing technique for selectively targeting cancer cells which have acquired molecular changes that distinguish them from normal cells. High-throughput RNAi-based screens have been successfully used to identify synthetic lethal pathways with well-characterized tumor suppressors and oncogenes. The recent identification of metabolic tumor suppressors suggests that the concept of synthetic lethality can be applied to selectively target cancer metabolism as well.
Here, we perform a high-throughput RNAi screen to identify synthetic lethal genes with fumarate hydratase (FH), a metabolic tumor suppressor whose loss-of-function has been associated with hereditary leiomyomatosis and renal cell carcinoma (HLRCC). Our unbiased screen identified synthetic lethality between FH and several genes in heme metabolism, in accordance with recent findings. Furthermore, we identified an enrichment of synthetic lethality with adenylate cyclases. The effects were validated in an embryonic kidney cell line (HEK293T) and in HLRCC-patient derived cells (UOK262) via both genetic and pharmacological inhibition. The reliance on adenylate cyclases in FH-deficient cells is consistent with increased cyclic-AMP levels, which may act to regulate cellular energy metabolism.
The identified synthetic lethality of FH with adenylate cyclases suggests a new potential target for treating HLRCC patients.
HLRCC; Fumarate hydratase-deficiency; High-throughput RNAi screen; Adenylate cyclase
Litter size in pigs is a major factor affecting the profitability in the pig industry. The peri-implantation window in pigs is characterized by the coordinated interactions between the maternal uterine endometrium and the rapidly elongating conceptuses and represents a period of time during which a large percentage of the developing conceptuses are lost. However, the gene expression and regulatory networks in the endometrium contributing to the establishment of the maternal: placental interface remain poorly understood.
We characterized the endometrial gene expression profile during the peri-implantation stage of development by comparing two breeds that demonstrate very different reproductive efficiencies. We employed the porcine Affymetrix GeneChip® to assay the transcriptomic profiles of genes expressed in the uterine endometrium obtained from Meishan and Yorkshire gilts (n = 4 for each breed) on day 12 of gestation (M12 and Y12, respectively). Total of 17,076 probesets were identified as "present" in at least two arrays. A mixed model-based statistical analysis predicted a total of 2,656 (q < 0.1) transcripts as differentially expressed between Meishan and Yorkshire pigs. Eighteen differentially expressed transcripts of interest were validated by quantitative real-time PCR. Gene ontology (GO) annotation revealed that the known functions of the differentially expressed genes were involved in a series of important biological processes relevant to early pregnancy establishment in the pig.
The results identified endometrial gene expression profiles of two breeds differing in litter size and identified candidate genes that are related to known physiological pathways related to reproductive prolificacy. These findings provide a deeper understanding of molecular pathways differing between two breeds at the critical peri-implantation stage of pregnancy, which can be utilized to better understand the events contributing to pregnancy establishment in the pig.
Pig; Prolificacy; Endometrium; Differentially expressed genes; Implantation
Defining cell type-specific transcriptomes in mammals can be challenging, especially for unannotated regions of the genome. We have developed an analytical pipeline called groHMM for annotating primary transcripts using global nuclear run-on sequencing (GRO-seq) data. Herein, we use this pipeline to characterize the transcriptome of an immortalized adult human ventricular cardiomyocyte cell line (AC16) in response to signaling by tumor necrosis factor alpha (TNFα), which is controlled in part by NF-κB, a key transcriptional regulator of inflammation. A unique aspect of this work is the use of the RNA polymerase II (Pol II) inhibitor α-amanitin, which we used to define a set of RNA polymerase I and III (Pol I and Pol III) transcripts.
Using groHMM, we identified ~30,000 coding and non-coding transcribed regions in AC16 cells, which includes a set of unique Pol I and Pol III primary transcripts. Many of these transcripts have not been annotated previously, including enhancer RNAs originating from NF-κB binding sites. In addition, we observed that AC16 cells rapidly and dynamically reorganize their transcriptomes in response to TNFα stimulation in an NF-κB-dependent manner, switching from a basal state to a proinflammatory state affecting a spectrum of cardiac-associated protein-coding and non-coding genes. Moreover, we observed distinct Pol II dynamics for up- and downregulated genes, with a rapid release of Pol II into productive elongation for TNFα-stimulated genes. As expected, the TNFα-induced changes in the AC16 transcriptome resulted in corresponding changes in cognate mRNA and protein levels in a similar manner, but with delayed kinetics.
Our studies illustrate how computational genomics can be used to characterize the signal-regulated transcriptome in biologically relevant cell types, providing new information about how the human genome is organized, transcribed and regulated. In addition, they show how α-amanitin can be used to reveal the Pol I and Pol III transcriptome. Furthermore, they shed new light on the regulation of the cardiomyocyte transcriptome in response to a proinflammatory signal and help to clarify the link between inflammation and cardiomyocyte function at the transcriptional level.
Expressed sequence tag (EST)-based markers are preferred because they reflect transcribed portions of the genome. We report the development of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers derived from transcriptome sequences in cabbage, and their utility for map construction.
Transcriptome sequences were obtained from two cabbage parental lines, C1184 and C1234, which are susceptible and resistant to black rot disease, respectively, using the 454 platform. A total of 92,255 and 127,522 reads were generated and clustered into 34,688 and 40,947 unigenes, respectively. We identified 2,405 SSR motifs from the unigenes of the black rot-resistant parent C1234. Trinucleotide motifs were the most abundant (66.15%) among the repeat motifs. In addition, 1,167 SNPs were detected between the two parental lines. A total of 937 EST-based SSR and 97 SNP-based dCAPS markers were designed and used for detection of polymorphism between parents. Using an F2 population, we built a genetic map comprising 265 loci, and consisting of 98 EST-based SSRs, 21 SNP-based dCAPS, 55 IBP markers derived from B. rapa genome sequence and 91 public SSRs, distributed on nine linkage groups spanning a total of 1,331.88 cM with an average distance of 5.03 cM between adjacent loci. The parental lines used in this study are elite breeding lines with little genetic diversity; therefore, the markers that mapped in our genetic map will have broad spectrum utility.
This genetic map provides additional genetic information to the existing B. oleracea map. Moreover, the new set of EST-based SSR and dCAPS markers developed herein is a valuable resource for genetic studies and will facilitate cabbage breeding. Additionally, this study demonstrates the usefulness of NGS transcriptomes for the development of genetic maps even with little genetic diversity in the mapping population.
Cabbage; EST; Genetic linkage map; SSR; SNP; Transcriptome sequencing
Regular smoking is associated with a wide variety of syndromes with prominent inflammatory components such as cancer, obesity and type 2 diabetes. Heavy regular smoking is also associated with changes in the DNA methylation of peripheral mononuclear cells. However, in younger smokers, inflammatory epigenetic findings are largely absent which suggests the inflammatory response(s) to smoking may be dose dependent. To help understand whether peripheral mononuclear cells have a role in mediating these responses in older smokers with higher cumulative smoke exposure, we examined genome-wide DNA methylation in a group of well characterized adult African American subjects informative for smoking, as well as serum C-reactive protein (CRP) and interleukin-6 receptor (IL6R) levels. In addition, complementary bioinformatic analyses were conducted to delineate possible pathways affected by long-term smoking.
Genome-wide DNA methylation analysis with respect to smoking status yielded 910 significant loci after Benjamini-Hochberg correction. In particular, two loci from the AHRR gene (cg05575921 and cg23576855) and one locus from the GPR15 gene (cg19859270) were identified as highly significantly differentially methylated between smokers and non-smokers. The bioinformatic analyses showed that long-term chronic smoking is associated with altered promoter DNA methylation of genes coding for proteins mapping to critical sub-networks moderating inflammation, immune function, and coagulation.
We conclude that chronic regular smoking is associated with changes in peripheral mononuclear cell methylation signature which perturb inflammatory and immune function pathways and may contribute to increased vulnerability for complex illnesses with inflammatory components.
microRNAs (miRNAs) are endogenous, noncoding, small RNAs that have essential regulatory functions in plant growth, development, and stress response processes. However, limited information is available about their functions in sexual reproduction of flowering plants. Pollen development is an important process in the life cycle of a flowering plant and is a major factor that affects the yield and quality of crop seeds.
This study aims to identify miRNAs involved in pollen development. Two independent small RNA libraries were constructed from the flower buds of the male sterile line (Bcajh97-01A) and male fertile line (Bcajh97-01B) of Brassica campestris ssp. chinensis. The libraries were subjected to high-throughput sequencing by using the Illumina Solexa system. Eight novel miRNAs on the other arm of known pre-miRNAs, 54 new conserved miRNAs, and 8 novel miRNA members were identified. Twenty-five pairs of novel miRNA/miRNA* were found. Among all the identified miRNAs, 18 differentially expressed miRNAs with over two-fold change between flower buds of male sterile line (Bcajh97-01A) and male fertile line (Bcajh97-01B) were identified. qRT-PCR analysis revealed that most of the differentially expressed miRNAs were preferentially expressed in flower buds of the male fertile line (Bcajh97-01B). Degradome analysis showed that a total of 15 genes were predicted to be the targets of seven miRNAs.
Our findings provide an overview of potential miRNAs involved in pollen development and interactions between miRNAs and their corresponding targets, which may provide important clues on the function of miRNAs in pollen development.
Brassica campestris; Brassica rapa; miRNAs; Pollen development; High-throughput sequencing; Deep sequencing; Degradome analysis
Penicillium chrysogenum has been used in producing penicillin and derived β-lactam antibiotics for many years. Although the genome of the mutant strain P. chrysogenum Wisconsin 54-1255 has already been sequenced, the versatility and genetic diversity of this species still needs to be intensively studied. In this study, the genome of the wild-type P. chrysogenum strain KF-25, which has high activity against Ustilaginoidea virens, was sequenced and characterized.
The genome of KF-25 was about 29.9 Mb in size and contained 9,804 putative open reading frames (orfs). Thirteen genes were predicted to encode two-component system proteins, of which six were putatively involved in osmolarity adaption. There were 33 putative secondary metabolism pathways and numerous genes that were essential in metabolite biosynthesis. Several P. chrysogenum virus untranslated region sequences were found in the KF-25 genome, suggesting that there might be a relationship between the virus and P. chrysogenum in evolution. Comparative genome analysis showed that the genomes of KF-25 and Wisconsin 54-1255 were highly similar, except that KF-25 was 2.3 Mb smaller. Three hundred and fifty-five KF-25 specific genes were found and the biological functions of the proteins encoded by these genes were mainly unknown (232, representing 65%), except for some orfs encoding proteins with predicted functions in transport, metabolism, and signal transduction. Numerous KF-25-specific genes were found to be associated with the pathogenicity and virulence of the strains, which were identical to those of wild-type P. chrysogenum NRRL 1951.
Genome sequencing and comparative analysis are helpful in further understanding the biology, evolution, and environment adaption of P. chrysogenum, and provide a new tool for identifying further functional metabolites.
Penicillium chrysogenum; Genome; Comparative genome
The set of all mRNA molecules present in a cell constitute the transcriptome. The transcriptome varies depending on cell type as well as in response to internal and external stimuli during development. Here we present a study of the changes that occur in the transcriptome of chili pepper fruit during development and ripening.
RNA-Seq was used to obtain transcriptomes of whole Serrano-type chili pepper fruits (Capsicum annuum L.; ‘Tampiqueño 74’) collected at 10, 20, 40 and 60 days after anthesis (DAA). 15,550,468 Illumina MiSeq reads were assembled de novo into 34,066 chili genes. We classified the expression patterns of individual genes as well as genes grouped into Biological Process ontologies and Metabolic Pathway categories using statistical criteria. For the analyses of gene groups we added the weighted expression of individual genes. This method was effective in interpreting general patterns of expression changes and increased the statistical power of the analyses. We also estimated the variation in diversity and specialization of the transcriptome during chili pepper development. Approximately 17% of genes exhibited a significant change of expression in at least one of the intervals sampled. In contrast, significant differences in approximately 63% of the Biological Processes and 80% of the Metabolic Pathways studied were detected in at least one interval. Confirming previous reports, genes related to capsaicinoid and ascorbic acid biosynthesis were significantly upregulated at 20 DAA while those related to carotenoid biosynthesis were highly expressed in the last period of fruit maturation (40–60 DAA). Our RNA-Seq data was validated by examining the expression of nine genes involved in carotenoid biosynthesis by qRT-PCR.
In general, more profound changes in the chili fruit transcriptome were observed in the intervals between 10 to 20 and 40 to 60 DAA. The last interval, between 40 to 60 DAA, included 49% of all significant changes detected, and was characterized predominantly by a global decrease in gene expression. This period signals the end of maturation and the beginning of senescence of chili pepper fruit. The transcriptome at 60 DAA was the most specialized and least diverse of the four states sampled.
Capsicum; Chili pepper; Transcriptome; Fruit development; Fruit ripening; Gene expression; RNA-Seq; High-throughput sequencing; Pathway analysis; Illumina MiSeq
In teleosts such as Atlantic salmon (Salmo salar L.), segmentation and subsequent mineralisation of the notochord during embryonic stages are essential for normal vertebrae formation. However, the molecular mechanisms leading to segmentation and mineralisation of the notochord are poorly understood. The aim of this study was to identify genes/pathways acting in gradients over time and along the anterior-posterior axis during notochord segmentation and immediately prior to mineralisation of the vertebral bodies in Atlantic salmon.
Notochord samples were collected from unsegmented, pre-segmented and segmented developmental stages. In each stage, the cellular core of the notochord was cut into three pieces along the longitudinal axis (anterior, mid, posterior). RNA was sequenced (22 million pair-end 100 bp/ library) and mapped to the salmon genome. 66569 transcripts were predicted and 55775 were annotated. In order to identify possible gradients leading to segmentation of the notochord, all 71 notochord-expressed hox genes were investigated, most of them displaying a typical anterior-posterior expression pattern along the notochord axis. The clustering of hox genes revealed a pattern that could be related to notochord segmentation. We further investigated how mineralisation is initiated in the notochord, and several factors related to chondrogenic lineage were identified (sox9, sox5, sox6, tgfb3, ihhb and col2a1), suggesting a cartilage-like character of the notochord. KEGG analysis of differentially expressed genes between stages revealed down-regulation of pathways associated with ECM, cell division, metabolism and development at onset of notochord segmentation. This implies that inhibitory signals produce segmentation of the notochord. One such potential inhibitory signal was identified, col11a2, which was detected in segments of non-mineralising notochord.
An incomplete salmon genome was successfully used to analyse RNA-seq data from the cellular core of the Atlantic salmon notochord. In transcriptome we found; hox gene patterns possibly linked to segmentation; down-regulation of pathways in the notochord at onset of segmentation; segmented expression of col11a2 in non-mineralised segments of the notochord; and a chondroblast-like footprint in the notochord.
Atlantic salmon; Notochord; RNA-seq; col11a2; Hox
Even in the age of next-generation sequencing (NGS), it has been unclear whether or not cells within a single organism have systematically distinctive genomes. Resolving this question, one of the most basic biological problems associated with DNA mutation rates, can assist efforts to elucidate essential mechanisms of cancer.
Using genome profiling (GP), we detected considerable systematic variation in genome sequences among cells in individual woody plants. The degree of genome sequence difference (genomic distance) varied systematically from the bottom to the top of the plant, such that the greatest divergence was observed between leaf genomes from uppermost branches and the remainder of the tree. This systematic variation was observed within both Yoshino cherry and Japanese beech trees.
As measured by GP, the genomic distance between two cells within an individual organism was non-negligible, and was correlated with physical distance (i.e., branch-to-branch distance). This phenomenon was assumed to be the result of accumulation of mutations from each cell division, implying that the degree of divergence is proportional to the number of generations separating the two cells.
Genome sequence; Genomic distance; Mutation rate; Japanese beech; Yoshino cherry; Leaf genomes; Genome profiling (GP); Next-generation sequencing (NGS)
MicroRNAs (miRNAs) are small non-coding RNAs found to regulate several biological processes including adipogenesis. Understanding adipose tissue regulation is critical for beef cattle as fat is an important determinant of beef quality and nutrient value. This study analyzed the association between genomic context characteristics of miRNAs with their expression and function in bovine adipose tissue. Twenty-four subcutaneous adipose tissue biopsies were obtained from eight British-continental crossbred steers at 3 different time points. Total RNA was extracted and miRNAs were profiled using a miRNA microarray with expression further validated by qRT-PCR.
A total of 224 miRNAs were detected of which 155 were expressed in all steers (n = 8), and defined as the core miRNAs of bovine subcutaneous adipose tissue. Core adipose miRNAs varied in terms of genomic location (59.5% intergenic, 38.7% intronic, 1.2% exonic, and 0.6% mirtron), organization (55.5% non-clustered and 44.5% clustered), and conservation (49% highly conserved, 14% conserved and 37% poorly conserved). Clustered miRNAs and highly conserved miRNAs were more highly expressed (p < 0.05) and had more predicted targets than non-clustered or less conserved miRNAs (p < 0.001). A total of 34 miRNAs were coordinately expressed, being part of six identified relevant networks. Two intronic miRNAs (miR-33a and miR-1281) were confirmed to have coordinated expression with their host genes, transcriptional factor SREBF2 and EP300 (a transcriptional co-activator of transcriptional factor C/EBPα), respectively which are involved in lipid metabolism, suggesting these miRNAs may also play a role in regulation of bovine lipid metabolism/adipogenesis. Furthermore, a total of 17 bovine specific miRNAs were predicted to be involved in the regulation of energy balance in adipose tissue.
These findings improve our understanding on the behavior of miRNAs in the regulation of bovine adipogenesis and fat metabolism as it reveals that miRNA expression patterns and functions are associated with miRNA genomic location, organization and conservation.
Adipogenesis; Adipose tissue; Bovine; Fat metabolism; Genomic context; microRNA; Cluster; Co-expression; Species specific
Aneuploidy, a karyotype deviating from multiples of a haploid chromosome set, affects the physiology of eukaryotes. In humans, aneuploidy is linked to pathological defects such as developmental abnormalities, mental retardation or cancer, but the underlying mechanisms remain elusive. There are many different types and origins of aneuploidy, but whether there is a uniform cellular response to aneuploidy in human cells has not been addressed so far.
Here we evaluate the transcription profiles of eleven trisomic and tetrasomic cell lines and two cell lines with complex aneuploid karyotypes. We identify a characteristic aneuploidy response pattern defined by upregulation of genes linked to endoplasmic reticulum, Golgi apparatus and lysosomes, and downregulation of DNA replication, transcription as well as ribosomes. Strikingly, complex aneuploidy elicits the same transcriptional changes as trisomy. To uncover the triggers of the response, we compared the profiles with transcription changes in human cells subjected to stress conditions. Interestingly, we found an overlap only with the response to treatment with the autophagy inhibitor bafilomycin A1. Finally, we identified 23 genes whose expression is significantly altered in all aneuploids and which may thus serve as aneuploidy markers.
Our analysis shows that despite the variability in chromosome content, aneuploidy triggers uniform transcriptional response in human cells. A common response independent of the type of aneuploidy might be exploited as a novel target for cancer therapy. Moreover, the potential aneuploidy markers identified in our analysis might represent novel biomarkers to assess the malignant potential of a tumor.
Investigation of tomato genetic resources is a crucial issue for better straight evolution and genetic studies as well as tomato breeding strategies. Traditional Vesuviano and San Marzano varieties grown in Campania region (Southern Italy) are famous for their remarkable fruit quality. Owing to their economic and social importance is crucial to understand the genetic basis of their unique traits.
Here, we present the draft genome sequences of tomato Vesuviano and San Marzano genome. A 40x genome coverage was obtained from a hybrid Illumina paired-end reads assembling that combines de novo assembly with iterative mapping to the reference S. lycopersicum genome (SL2.40). Insertions, deletions and SNP variants were carefully measured. When assessed on the basis of the reference annotation, 30% of protein-coding genes are predicted to have variants in both varieties. Copy genes number and gene location were assessed by mRNA transcripts mapping, showing a closer relationship of San Marzano with reference genome. Distinctive variations in key genes and transcription/regulation factors related to fruit quality have been revealed for both cultivars.
The effort performed highlighted varieties relationships and important variants in fruit key processes useful to dissect the path from sequence variant to phenotype.
Combined assembling; Fruit quality; NGS sequencing; SNPs; Solanum lycopersicum
Saccharomyces cerevisiae strains isolated from natural settings form structured biofilm colonies that are equipped with intricate protective mechanisms. These wild strains are able to reprogram themselves with a certain frequency during cultivation in plentiful laboratory conditions. The resulting domesticated strains switch off certain protective mechanisms and form smooth colonies that resemble those of common laboratory strains.
Here, we show that domestication can be reversed when a domesticated strain is challenged by various adverse conditions; the resulting feral strain restores its ability to form structured biofilm colonies. Phenotypic, microscopic and transcriptomic analyses show that phenotypic transition is a complex process that affects various aspects of feral strain physiology; it leads to a phenotype that resembles the original wild strain in some aspects and the domesticated derivative in others. We specify the genetic determinants that are likely involved in the formation of a structured biofilm colonies. In addition to FLO11, these determinants include genes that affect the cell wall and membrane composition. We also identify changes occurring during phenotypic transitions that affect other properties of phenotypic strain-variants, such as resistance to the impact of environmental stress. Here we document the regulatory role of the histone deacetylase Hda1p in developing such a resistance.
We provide detailed analysis of transcriptomic and phenotypic modulations of three related S. cerevisiae strains that arose by phenotypic switching under diverse environmental conditions. We identify changes specifically related to a strain’s ability to create complex structured colonies; we also show that other changes, such as genome rearrangement(s), are unrelated to this ability. Finally, we identify the importance of histone deacetylase Hda1p in strain resistance to stresses.
Biofilm colony; Histone deacetylase; Phenotypic switching; Wild yeast strains
Royal jelly (RJ) is a proteinaceous secretion produced from the hypopharyngeal and mandibular glands of nurse bees. It plays vital roles in honeybee biology and in the improvement of human health. However, some proteins remain unknown in RJ, and mapping N-glycosylation modification sites on RJ proteins demands further investigation. We used two different liquid chromatography-tandem mass spectrometry techniques, complementary N-glycopeptide enrichment strategies, and bioinformatic approaches to gain a better understanding of novel and glycosylated proteins in RJ.
A total of 25 N-glycosylated proteins, carrying 53 N-glycosylation sites, were identified in RJ proteins, of which 42 N-linked glycosylation sites were mapped as novel on RJ proteins. Most of the glycosylated proteins were related to metabolic activities and health improvement. The 13 newly identified proteins were also mainly associated with metabolic processes and health improvement activities.
Our in-depth, large-scale mapping of novel glycosylation sites represents a crucial step toward systematically revealing the functionality of N-glycosylated RJ proteins, and is potentially useful for producing a protein with desirable pharmacokinetic and biological activity using a genetic engineering approach. The newly-identified proteins significantly extend the proteome coverage of RJ. These findings contribute vital and new knowledge to our understanding of the innate biochemical nature of RJ at both the proteome and glycoproteome levels.
Royal jelly; N-glycosylation; Hydrazide chemistry; Lectin affinity; Tandem mass spectrometry
Honey bees (Apis mellifera) contribute substantially to the worldwide economy and ecosystem health as pollinators. Pollen is essential to the bee’s diet, providing protein, lipids, and micronutrients. The dramatic shifts in physiology, anatomy, and behavior that accompany normal worker development are highly plastic and recent work demonstrates that development, particularly the transition from nurse to foraging roles, is greatly impacted by diet. However, the role that diet plays in the developmental transition of newly eclosed bees to nurse workers is poorly understood. To further understand honey bee nutrition and the role of diet in nurse development, we used a high-throughput screen of the transcriptome of 3 day and 8 day old worker bees fed either honey and stored pollen (rich diet) or honey alone (poor diet) within the hive. We employed a three factor (age, diet, age x diet) analysis of the transcriptome to determine whether diet affected nurse worker physiology and whether poor diet altered the developmental processes normally associated with aging.
Substantial changes in gene expression occurred due to starvation. Diet-induced changes in gene transcription occurring in younger bees were largely a subset of those occurring in older bees, but certain signatures of starvation were only evident 8 day old workers. Of the 18,542 annotated transcripts in the A. mellifera genome, 150 transcripts exhibited differential expression due to poor diet at 3d of age compared with 17,226 transcripts that differed due to poor diet at 8d of age, and poor diet caused more frequent down-regulation of gene expression in younger bees compared to older bees. In addition, the age-related physiological changes that accompanied early adult development differed due to the diet these young adult bees were fed. More frequent down-regulation of gene expression was observed in developing bees fed a poor diet compared to those fed an adequate diet. Functional analyses also suggest that the physiological and developmental processes occurring in well-fed bees are vastly different than those occurring in pollen deprived bees. Our data support the hypothesis that poor diet causes normal age-related development to go awry.
Poor nutrition has major consequences for the expression of genes underlying the physiology and age-related development of nurse worker bees. More work is certainly needed to fully understand the consequences of starvation and the complex biology of nutrition and development in this system, but the genes identified in the present study provide a starting point for understanding the consequences of poor diet and for mitigating the economic costs of colony starvation.
Transcriptome; Apis mellifera; Nutrition; Starvation; Nurse; Development; Physiology
Genome-wide expression profiles are altered during biological aging and can describe molecular regulation of tissue degeneration. Age-regulated mRNA expression trends from cross-sectional studies could describe how aging progresses. We developed a novel statistical methodology to identify age-regulated expression trends in cross-sectional datasets.
We studied six cross-sectional RNA expression profiles from different human tissues. Our methodology, capable of overcoming technical and genetic background differences, identified an age-regulation in four of the tissues. For the identification of expression trends, five regression models were compared and the quadratic model was found as the most suitable for this study. After k-means clustering of the age-associated probes, expression trends were found to change at two major age-positions in brain cortex and in Vastus lateralis muscles. The first age-position was found to occur during the fifth decade and a later one during the eighth decade. In kidney cortex, however, only one age-position was identified correlating with a late age-position. Functional mapping of genes at each age-position suggests that calcium homeostasis and lipid metabolisms are initially affected and subsequently, in elderly mitochondria, apoptosis and hormonal signaling pathways are affected.
Our results suggest that age-associated temporal changes in human tissues progress at distinct age-positions, which differ between tissues and in their molecular composition.
Human aging; Expression profiles; Quadratic regression model; Kmeans clustering