Genome Transcriptional Profiling Research Articles

Genetic variability and phylogeny of Torque teno sus virus 1 (TTSuV1) and 2 (TTSuV2) based on complete genomes

Thirteen genomes of Torque teno sus virus 1 (TTSuV1) and 2 (TTSuV2) were obtained to examine the diversity and evolution of swine TTVs. Despite the low nucleotide identity reported, the genomic organization and transcriptional profiles of TTVs are similar. The nucleotide diversity for TTSuV1 was higher than TTSuV2, and the pattern of mutation among the ORFs was also different. Phylogenetic and genetic analyses support the proposed division of TTV into two species. TTSuV1 showed high levels of variability (>30%), with three different types (the third described for the first time) that may display a geographical structure. In contrast, TTSuV2 showed lower levels of variability (<15%), and no different types could be described. Larger values for the ratios of synonymous (dS) to non-synonymous (dN) base substitutions (dS/dN) were reported for the ORFs pointing to a certain level of selective constraint in TTV genomes.

Comparative transcriptional profiling identifies takeout as a gene that regulates life span

A major challenge in translating the positive effects of dietary restriction (DR) for the improvement of human health is the development of therapeutic mimics. One approach to finding DR mimics is based upon identification of the proximal effectors of DR life span extension. Whole genome profiling of DR in Drosophila shows a large number of changes in gene expression, making it difficult to establish which changes are involved in life span determination as opposed to other unrelated physiological changes. We used comparative whole genome expression profiling to discover genes whose change in expression is shared between DR and two molecular genetic life span extending interventions related to DR, increased dSir2 and decreased Dmp53 activity. We find twenty-one genes shared among the three related life span extending interventions. One of these genes, takeout, thought to be involved in circadian rhythms, feeding behavior and juvenile hormone binding is also increased in four other life span extending conditions: Rpd3, Indy, chico and methuselah. We demonstrate takeout is involved in longevity determination by specifically increasing adult takeout expression and extending life span. These studies demonstrate the power of comparative whole genome transcriptional profiling for identifying specific downstream elements of the DR life span extending pathway.

Oxidative stress protection and the repair response to hydrogen peroxide in the hyperthermophilic archaeon Pyrococcus furiosus and in related species

Pyrococcus furiosus is a shallow marine, anaerobic archaeon that grows optimally at 100 degrees C. Addition of H(2)O(2) (0.5 mM) to a growing culture resulted in the cessation of growth with a 2-h lag before normal growth resumed. Whole genome transcriptional profiling revealed that the main response occurs within 30 min of peroxide addition, with the up-regulation of 62 open reading frames (ORFs), 36 of which are part of 10 potential operons. More than half of the up-regulated ORFs are of unknown function, while some others encode proteins that are involved potentially in sequestering iron and sulfide, in DNA repair and in generating NADPH. This response is thought to involve primarily damage repair rather than protection, since cultures exposed to sub-toxic levels of H(2)O(2) were not more resistant to the subsequent addition of H(2)O(2) (0.5-5.0 mM). Consequently, there is little if any induced protective response to peroxide. The organism maintains a constitutive protective mechanism involving high levels of oxidoreductase-type enzymes such as superoxide reductase, rubrerythrin, and alkyl hydroperoxide reductase. Related hyperthermophiles contain homologs of the proteins involved in the constitutive protective mechanism but these organisms were more sensitive to peroxide than P. furiosus and lack several of its peroxide-responsive ORFs.

Correcting for intra-experiment variation in Illumina BeadChip data is necessary to generate robust gene-expression profiles

BackgroundMicroarray technology is a popular means of producing whole genome transcriptional profiles, however high cost and scarcity of mRNA has led many studies to be conducted based on the analysis of single samples. We exploit the design of the Illumina platform, specifically multiple arrays on each chip, to evaluate intra-experiment technical variation using repeated hybridisations of universal human reference RNA (UHRR) and duplicate hybridisations of primary breast tumour samples from a clinical study.ResultsA clear batch-specific bias was detected in the measured expressions of both the UHRR and clinical samples. This bias was found to persist following standard microarray normalisation techniques. However, when mean-centering or empirical Bayes batch-correction methods (ComBat) were applied to the data, inter-batch variation in the UHRR and clinical samples were greatly reduced. Correlation between replicate UHRR samples improved by two orders of magnitude following batch-correction using ComBat (ranging from 0.9833-0.9991 to 0.9997-0.9999) and increased the consistency of the gene-lists from the duplicate clinical samples, from 11.6% in quantile normalised data to 66.4% in batch-corrected data. The use of UHRR as an inter-batch calibrator provided a small additional benefit when used in conjunction with ComBat, further increasing the agreement between the two gene-lists, up to 74.1%.ConclusionIn the interests of practicalities and cost, these results suggest that single samples can generate reliable data, but only after careful compensation for technical bias in the experiment. We recommend that investigators appreciate the propensity for such variation in the design stages of a microarray experiment and that the use of suitable correction methods become routine during the statistical analysis of the data.

Whole genome transcriptional profiling identifies novel differentiation regulated genes in keratinocytes

Keratinocyte differentiation plays a pivotal role in the epidermal barrier. Single keratinocyte differentiation genes have already been studied, but many important constituents of this process may have been missed so far. Gene expression profiling by microarray was carried out in cultured normal human epidermal keratinocytes undergoing confluence-induced differentiation to find novel differentiation genes. Candidate gene lists were established and genes of potential dermatological interest were validated by quantitative reverse transcription polymerase chain reaction and immunohistochemical analysis. Some of these points lead to the identification of counter-regulation of heme oxygenase and biliverdin reductase as well as glutaredoxin and glutathione reductase indicative of potential novel redox signaling in differentiating human keratinocytes. Others indicate a strong concert down-regulation of interleukin-1 signaling at previously unidentified levels during keratinocyte differentiation. We believe that identified genes contribute to a more comprehensive understanding of the complicated epidermal differentiation process and lead to better understanding of dermatological diseases.

Two Regulators of Vibrio parahaemolyticus Play Important Roles in Enterotoxicity by Controlling the Expression of Genes in the Vp-PAI Region

Vibrio parahaemolyticus is an important pathogen causing food-borne disease worldwide. An 80-kb pathogenicity island (Vp-PAI), which contains two tdh (thermostable direct hemolysin) genes and a set of genes for the type III secretion system (T3SS2), is closely related to the pathogenicity of this bacterium. However, the regulatory mechanisms of Vp-PAI's gene expression are poorly understood. Here we report that two novel ToxR-like transcriptional regulatory proteins (VtrA and VtrB) regulate the expression of the genes encoded within the Vp-PAI region, including those for TDH and T3SS2-related proteins. Expression of vtrB was under control of the VtrA, as vector-expressed vtrB was able to recover a functional protein secretory capacity for T3SS2, independent of VtrA. Moreover, these regulatory proteins were essential for T3SS2-dependent biological activities, such as in vitro cytotoxicity and in vivo enterotoxicity. Enterotoxic activities of vtrA and/or vtrB deletion strains derived from the wild-type strain were almost absent, showing fluid accumulation similar to non-infected control. Whole genome transcriptional profiling of vtrA or vtrB deletion strains revealed that the expression levels of over 60 genes were downregulated significantly in these deletion mutant strains and that such genes were almost exclusively located in the Vp-PAI region. These results strongly suggest that VtrA and VtrB are master regulators for virulence gene expression in the Vp-PAI and play critical roles in the pathogenicity of this bacterium.

Transcriptional profiling reveals divergent roles of PPARα and PPARβ/δ in regulation of gene expression in mouse liver

Little is known about the role of the transcription factor peroxisome proliferator-activated receptor (PPAR) beta/delta in liver. Here we set out to better elucidate the function of PPARbeta/delta in liver by comparing the effect of PPARalpha and PPARbeta/delta deletion using whole genome transcriptional profiling and analysis of plasma and liver metabolites. In fed state, the number of genes altered by PPARalpha and PPARbeta/delta deletion was similar, whereas in fasted state the effect of PPARalpha deletion was much more pronounced, consistent with the pattern of gene expression of PPARalpha and PPARbeta/delta. Minor overlap was found between PPARalpha- and PPARbeta/delta-dependent gene regulation in liver. Pathways upregulated by PPARbeta/delta deletion were connected to innate immunity and inflammation. Pathways downregulated by PPARbeta/delta deletion included lipoprotein metabolism and various pathways related to glucose utilization, which correlated with elevated plasma glucose and triglycerides and reduced plasma cholesterol in PPARbeta/delta-/- mice. Downregulated genes that may underlie these metabolic alterations included Pklr, Fbp1, Apoa4, Vldlr, Lipg, and Pcsk9, which may represent novel PPARbeta/delta target genes. In contrast to PPARalpha-/- mice, no changes in plasma free fatty acid, plasma beta-hydroxybutyrate, liver triglycerides, and liver glycogen were observed in PPARbeta/delta-/- mice. Our data indicate that PPARbeta/delta governs glucose utilization and lipoprotein metabolism and has an important anti-inflammatory role in liver. Overall, our analysis reveals divergent roles of PPARalpha and PPARbeta/delta in regulation of gene expression in mouse liver.

In vitromanipulation of gene expression in larvalSchistosoma: a model for postgenomic approaches in Trematoda

With rapid developments in DNA and protein sequencing technologies, combined with powerful bioinformatics tools, a continued acceleration of gene identification in parasitic helminths is predicted, potentially leading to discovery of new drug and vaccine targets, enhanced diagnostics and insights into the complex biology underlying host-parasite interactions. For the schistosome blood flukes, with the recent completion of genome sequencing and comprehensive transcriptomic datasets, there has accumulated massive amounts of gene sequence data, for which, in the vast majority of cases, little is known about actual functions within the intact organism. In this review we attempt to bring together traditional in vitro cultivation approaches and recent emergent technologies of molecular genomics, transcriptomics and genetic manipulation to illustrate the considerable progress made in our understanding of trematode gene expression and function during development of the intramolluscan larval stages. Using several prominent trematode families (Schistosomatidae, Fasciolidae, Echinostomatidae), we have focused on the current status of in vitro larval isolation/cultivation as a source of valuable raw material supporting gene discovery efforts in model digeneans that include whole genome sequencing, transcript and protein expression profiling during larval development, and progress made in the in vitro manipulation of genes and their expression in larval trematodes using transgenic and RNA interference (RNAi) approaches.

Whole genome transcript profiling from fingerstick blood samples: a comparison and feasibility study

BackgroundWhole genome gene expression profiling has revolutionized research in the past decade especially with the advent of microarrays. Recently, there have been significant improvements in whole blood RNA isolation techniques which, through stabilization of RNA at the time of sample collection, avoid bias and artifacts introduced during sample handling. Despite these improvements, current human whole blood RNA stabilization/isolation kits are limited by the requirement of a venous blood sample of at least 2.5 mL. While fingerstick blood collection has been used for many different assays, there has yet to be a kit developed to isolate high quality RNA for use in gene expression studies from such small human samples. The clinical and field testing advantages of obtaining reliable and reproducible gene expression data from a fingerstick are many; it is less invasive, time saving, more mobile, and eliminates the need of a trained phlebotomist. Furthermore, this method could also be employed in small animal studies, i.e. mice, where larger sample collections often require sacrificing the animal. In this study, we offer a rapid and simple method to extract sufficient amounts of high quality total RNA from approximately 70 μl of whole blood collected via a fingerstick using a modified protocol of the commercially available Qiagen PAXgene RNA Blood Kit.ResultsFrom two sets of fingerstick collections, about 70 uL whole blood collected via finger lancet and capillary tube, we recovered an average of 252.6 ng total RNA with an average RIN of 9.3. The post-amplification yields for 50 ng of total RNA averaged at 7.0 ug cDNA. The cDNA hybridized to Affymetrix HG-U133 Plus 2.0 GeneChips had an average % Present call of 52.5%. Both fingerstick collections were highly correlated with r2 values ranging from 0.94 to 0.97. Similarly both fingerstick collections were highly correlated to the venous collection with r2 values ranging from 0.88 to 0.96 for fingerstick collection 1 and 0.94 to 0.96 for fingerstick collection 2.ConclusionsOur comparisons of RNA quality and gene expression data of the fingerstick method with traditionally processed sample workflows demonstrate excellent RNA quality from the capillary collection as well as very high correlations of gene expression data.

Abstract C36: Using the Rb/E2f3 mutant mouse model to reveal markers for the progression to metastasis

Abstract The retinoblastoma tumor suppressor protein (pRb) was the first tumor suppressor to be identified. The E2F family of transcription factors is considered the best-studied class of interaction partners of pRB. Specifically, E2F1-3 are all capable to control re-entry and progression through the S-phase of the cell cycle. The simplified view is that hypo-phosphorylated pRb is bound by E2F and recruits chromatin remodeling complexes to represses the transcription of E2F target genes. Upon cell cycle progression pRb is phosphorylated leading to a release of E2Fs and the subsequent progression through the G1/S-transition. In order to understand the consequences of deregulation of the pRb/E2F-pathway, a number of groups have established mouse knock-out models of various members of this pathway. Using a combined knockout mouse model it was demonstrated that E2F3 contributes to both the inappropriate proliferation and p53-dependent apoptosis that arises in pRb-deficient embryos. Furthermore, it was shown that absence of E2F3 results on one hand in the attenuation of the growth of pituitary adenomas, on the other hand it promotes the development of medullary thyroid carcinomas (MTCs) yielding metastases at an increased frequency. Since this is one of the few existing mouse models to study metastasis, we wished to use this model to identify novel markers and mechanisms that regulate this process. We established the whole genomic transcriptional profile of mouse MTCs harvested from Rb/E2f3 mutant mice. The comparison of metastatic MTCs with non-metastatic tumors defined a data set of 123 putative metastatic markers. Among this data set, we identified a cluster of 22 genes that are directly involved in the formation of the kinetochore, establishment of the centromere or play a distinct role in the spindle assembly checkpoint. We validated the expression of these novel metastatic markers using mouse RNA, as well as material from human metastatic MTC biopsies. Importantly, we were able to establish that the same genes were only upregulated in metastatic tumors utilizing immunhistochemistry on sections of human metastatic MTCs but not in non-metastatic counterparts. Moreover, chromatin immunoprecipitation assays support that these genes are direct targets of the pRB/E2F-pathway. In summary we identified novel markers of metastasis. Moreover, these markers are also putative E2F target genes, the fact of which emphasizes the important contribution of E2Fs to tumor progression. Citation Information: Cancer Res 2009;69(23 Suppl):C36.

Increased nitrogen‐use efficiency in transgenic rice plants over‐expressing a nitrogen‐responsive early nodulin gene identified from rice expression profiling

Development of genetic varieties with improved nitrogen-use efficiency (NUE) is essential for sustainable agriculture. In this study, we developed a growth system for rice wherein N was the growth-limiting factor, and identified N-responsive genes by a whole genome transcriptional profiling approach. Some genes were selected to test their functionality in NUE by a transgenic approach. One such example with positive effects on NUE is an early nodulin gene OsENOD93-1. This OsENOD93-1 gene responded significantly to both N induction and N reduction. Transgenic rice plants over-expressing the OsENOD93-1 gene had increased shoot dry biomass and seed yield. This OsENOD93-1 gene was expressed at high levels in roots of wild-type (WT) plants, and its protein product was localized in mitochondria. Transgenic plants accumulated higher concentrations of total amino acids and total N in roots. A higher concentration of amino acids in xylem sap was detected in transgenic plants, especially under N stress. In situ hybridization revealed that OsENOD93-1 is expressed in vascular bundles, as well as in epidermis and endodermis. This work demonstrates that transcriptional profiling, coupled with a transgenic validation approach, is an effective strategy for gene discovery. The knowledge gained from this study could be applied to other important crops.

Lifestyle of Lactobacillus plantarum in the mouse caecum

Lactobacillus plantarum is a common inhabitant of mammalian gastrointestinal tracts. Strains of L. plantarum are also marketed as probiotics intended to confer beneficial health effects upon delivery to the human gut. To understand how L. plantarum adapts to its gut habitat, we used whole genome transcriptional profiling to characterize the transcriptome of strain WCFS1 during colonization of the caeca of adult germ-free C57Bl/6 J mice fed a standard low-fat rodent chow diet rich in complex plant polysaccharides or a prototypic Western diet high in simple sugars and fat. Lactobacillus plantarum colonized the digestive tracts of these animals to high levels, although L. plantarum was found in 10-fold higher amounts in the caeca of mice fed the standard chow. Metabolic reconstructions based on the transcriptional data sets revealed that genes involved in carbohydrate transport and metabolism form the principal functional group that is upregulated in vivo compared with exponential phase cells grown in three different culture media, and that a Western diet provides a more nutritionally restricted, growth limiting milieu for the microbe in the distal gut. A set of bacterial genes encoding cell surface-related functions were differentially regulated in both groups of mice. This set included downregulated genes required for the d-alanylation of lipoteichoic acids, extracellular structures of L. plantarum that mediate interactions with the host immune system. These results, obtained in a reductionist gnotobiotic mouse model of the gut ecosystem, provide insights about the niches (professions) of this lactic acid bacterium, and a context for systematically testing features that affect epithelial and immune cell responses to this organism in the digestive tract.

Transcriptional Profiling of the Dose Response: A More Powerful Approach for Characterizing Drug Activities

The dose response curve is the gold standard for measuring the effect of a drug treatment, but is rarely used in genomic scale transcriptional profiling due to perceived obstacles of cost and analysis. One barrier to examining transcriptional dose responses is that existing methods for microarray data analysis can identify patterns, but provide no quantitative pharmacological information. We developed analytical methods that identify transcripts responsive to dose, calculate classical pharmacological parameters such as the EC50, and enable an in-depth analysis of coordinated dose-dependent treatment effects. The approach was applied to a transcriptional profiling study that evaluated four kinase inhibitors (imatinib, nilotinib, dasatinib and PD0325901) across a six-logarithm dose range, using 12 arrays per compound. The transcript responses proved a powerful means to characterize and compare the compounds: the distribution of EC50 values for the transcriptome was linked to specific targets, dose-dependent effects on cellular processes were identified using automated pathway analysis, and a connection was seen between EC50s in standard cellular assays and transcriptional EC50s. Our approach greatly enriches the information that can be obtained from standard transcriptional profiling technology. Moreover, these methods are automated, robust to non-optimized assays, and could be applied to other sources of quantitative data.

Comparison of the Histology, Gene Expression Profile, and Phenotype of Cultured Human Limbal Epithelial Cells from Different Limbal Regions

To investigate whether human limbal epithelial cells (HLECs) derived from various regions of the limbus exhibit differences in gene expression and epithelial characteristics. HLECs were derived from explants taken from the superior, nasal, inferior, and temporal limbus and cultured for 21 days. Whole genome transcript profiling was performed with a gene microarray. The microarray results were validated by using RT-PCR. Epithelial morphology was studied with light microscopy and transmission electron microscopy, and phenotype was evaluated by immunohistochemistry. Epithelial outgrowth was present in most cultures of superior origin (88%) in contrast to cultures of temporal origin (38%). The epithelial thickness and number of cell layers were significantly greater in cultures of superior origin than in cultures from inferior and temporal areas. TRIM36, OSR2, and RHOU, which are involved in morphogenesis, were significantly differentially expressed in the superior region, compared with the other regions. Proposed limbal stem cell, progenitor, and differentiation markers were not differentially expressed. The uniform gene expression of ocular surface markers correlated with homogeneous immunostaining of corresponding protein markers in HLEC cultures from all regions, demonstrating an undifferentiated phenotype (p63(+), DeltaNp63alpha(+), ABCG2(+), K19(+), vimentin(+), integrin beta1(+), nestin(-), K3(-), K5(+), and E-cadherin(+)). No major transcriptional or phenotypic differences were observed in cultured HLECs derived from different regions of the limbus. However, explants of superior origin demonstrated the highest outgrowth success rate and generated epithelia with greater epithelial thickness and number of cell layers, which may prove useful for transplantation purposes.

Characterization of a Naturally Occurring Breast Cancer Subset Enriched in Epithelial-to-Mesenchymal Transition and Stem Cell Characteristics

Metaplastic breast cancers (MBC) are aggressive, chemoresistant tumors characterized by lineage plasticity. To advance understanding of their pathogenesis and relatedness to other breast cancer subtypes, 28 MBCs were compared with common breast cancers using comparative genomic hybridization, transcriptional profiling, and reverse-phase protein arrays and by sequencing for common breast cancer mutations. MBCs showed unique DNA copy number aberrations compared with common breast cancers. PIK3CA mutations were detected in 9 of 19 MBCs (47.4%) versus 80 of 232 hormone receptor-positive cancers (34.5%; P = 0.32), 17 of 75 HER-2-positive samples (22.7%; P = 0.04), 20 of 240 basal-like cancers (8.3%; P < 0.0001), and 0 of 14 claudin-low tumors (P = 0.004). Of 7 phosphatidylinositol 3-kinase/AKT pathway phosphorylation sites, 6 were more highly phosphorylated in MBCs than in other breast tumor subtypes. The majority of MBCs displayed mRNA profiles different from those of the most common, including basal-like cancers. By transcriptional profiling, MBCs and the recently identified claudin-low breast cancer subset constitute related receptor-negative subgroups characterized by low expression of GATA3-regulated genes and of genes responsible for cell-cell adhesion with enrichment for markers linked to stem cell function and epithelial-to-mesenchymal transition (EMT). In contrast to other breast cancers, claudin-low tumors and most MBCs showed a significant similarity to a "tumorigenic" signature defined using CD44(+)/CD24(-) breast tumor-initiating stem cell-like cells. MBCs and claudin-low tumors are thus enriched in EMT and stem cell-like features, and may arise from an earlier, more chemoresistant breast epithelial precursor than basal-like or luminal cancers. PIK3CA mutations, EMT, and stem cell-like characteristics likely contribute to the poor outcomes of MBC and suggest novel therapeutic targets.

Temporal genomic transcriptional profiling of DSS induced colitis by microarray

We have previously reported that the dextran sodium sulfate (DSS) induced colitis model involves a complex relationship between inflammation and angiogenic activity in that inflammation and angiogenesis serve to perpetuate chronic inflammation. Here we performed genome transcription profile analysis of 3% DSS induced colitis using the Affymetrix mouse genome 430 2.0 Array from colon tissues taken at day 0 (D0), 2 (D2), 4 (D4), and 6 (D6). Several inflammatory genes were greatly up‐regulated, such as Interleukin‐6 (0.9, 8.6 and 80.7 fold increase at D2, D4, and D6), Mmp 10 (1, 3, 22 fold increase at D2, D4, and D6), and Mmp 13 (2 and 21 fold increase at D4 and D6). Angiogenesis related genes were also significantly up‐regulated, such as HGF (0.4 and 2.8 fold increase at D4 and D6), Map3k8 (3 fold increase at D6), and Epiregulin (3.7 fold increase at D6). Increased inflammatory and angiogenic gene expression was observed in a proportional manner with earlier changes in inflammatory gene expression. Several genes were also significantly down regulated such as Trpm6 (38%, 43% and 71% at D2, D4 and D6), Repin1 (16%, 54%, and 76% at D2, D4 and D6), and Pnpla6 (41%, 59% and 56% at D2, D4 and D6). Our study demonstrates that inflammatory and angiogenic activities are linked during experimental colitis and delineate distinct time points of genome transcriptional change which may suggest new targets for the treatment of colitis.

DNA microarray technology and its application in fish biology and aquaculture

Fishery is an important industry in China as well as in the rest of the world, and it provides a human food resource containing high-quality protein. Best practice in aquaculture requires a full understanding of the genomic controls and transcriptional profiles of cultured fish species. Improvements in aquaculture can be made by regulation of the expression of functional genes. Microarray technology is a powerful tool for rapid screening of genes or transcriptional profiles in a particular fish or for a particular economic character; for example, genes that are related to growth and disease control in the fish. This review provides a brief introduction to microarray technology and its methods and applications, together with a discussion of the achievements in fish biology that have resulted from this technology.

Glucose and Auxin Signaling Interaction in Controlling Arabidopsis thaliana Seedlings Root Growth and Development

BackgroundPlant root growth and development is highly plastic and can adapt to many environmental conditions. Sugar signaling has been shown to affect root growth and development by interacting with phytohormones such as gibberellins, cytokinin and abscisic acid. Auxin signaling and transport has been earlier shown to be controlling plant root length, number of lateral roots, root hair and root growth direction.Principal FindingsIncreasing concentration of glucose not only controls root length, root hair and number of lateral roots but can also modulate root growth direction. Since root growth and development is also controlled by auxin, whole genome transcript profiling was done to find out the extent of interaction between glucose and auxin response pathways. Glucose alone could transcriptionally regulate 376 (62%) genes out of 604 genes affected by IAA. Presence of glucose could also modulate the extent of regulation 2 fold or more of almost 63% genes induced or repressed by IAA. Interestingly, glucose could affect induction or repression of IAA affected genes (35%) even if glucose alone had no significant effect on the transcription of these genes itself. Glucose could affect auxin biosynthetic YUCCA genes family members, auxin transporter PIN proteins, receptor TIR1 and members of a number of gene families including AUX/IAA, GH3 and SAUR involved in auxin signaling. Arabidopsis auxin receptor tir1 and response mutants, axr2, axr3 and slr1 not only display a defect in glucose induced change in root length, root hair elongation and lateral root production but also accentuate glucose induced increase in root growth randomization from vertical suggesting glucose effects on plant root growth and development are mediated by auxin signaling components.ConclusionOur findings implicate an important role of the glucose interacting with auxin signaling and transport machinery to control seedling root growth and development in changing nutrient conditions.

Transcript length bias in RNA-seq data confounds systems biology

BackgroundSeveral recent studies have demonstrated the effectiveness of deep sequencing for transcriptome analysis (RNA-seq) in mammals. As RNA-seq becomes more affordable, whole genome transcriptional profiling is likely to become the platform of choice for species with good genomic sequences. As yet, a rigorous analysis methodology has not been developed and we are still in the stages of exploring the features of the data.ResultsWe investigated the effect of transcript length bias in RNA-seq data using three different published data sets. For standard analyses using aggregated tag counts for each gene, the ability to call differentially expressed genes between samples is strongly associated with the length of the transcript.ConclusionTranscript length bias for calling differentially expressed genes is a general feature of current protocols for RNA-seq technology. This has implications for the ranking of differentially expressed genes, and in particular may introduce bias in gene set testing for pathway analysis and other multi-gene systems biology analyses.ReviewersThis article was reviewed by Rohan Williams (nominated by Gavin Huttley), Nicole Cloonan (nominated by Mark Ragan) and James Bullard (nominated by Sandrine Dudoit).

Genome-wide innate immune responses in HIV-1-infected macrophages are preserved despite attenuation of the NF-kappa B activation pathway.

Macrophages contribute to HIV-1 infection at many levels. They provide permissive cells at the site of inoculation, augment virus transfer to T cells, generate long-lived viral reservoirs, and cause bystander cell apoptosis. A body of evidence suggests that the role of macrophages in cellular host defense is also compromised by HIV-1 infection. In this respect, macrophages are potent cells of the innate immune system that initiate and regulate wide-ranging immunological responses. This study focuses on the effect of HIV-1 infection on innate immune responses by macrophages at the level of signal transduction, whole genome transcriptional profiling, and cytokine secretion. We show that in an ex vivo model, M-CSF-differentiated monocyte-derived macrophages uniformly infected with replicating CCR5-tropic HIV-1, without cytopathic effect, exhibit selective attenuation of the NF-kappaB activation pathway in response to TLR4 and TLR2 stimulation. However, functional annotation clustering analysis of genome-wide transcriptional responses to LPS stimulation suggests substantial preservation of gene expression changes at the systems level, with modest attenuation of a subset of up-regulated LPS-responsive genes, and no effect on a selection of inflammatory cytokine responses at the protein level. These results extend existing reports of inhibitory interactions between HIV-1 accessory proteins and NF-kappaB signaling pathways, and whole genome expression profiling provides comprehensive assessment of the consequent effects on immune response gene expression. Unexpectedly, our data suggest innate immune responses are broadly preserved with limited exceptions, and pave the way for further study of the complex relationship between HIV-1 and immunological pathways within macrophages.