Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

Ascl1 coordinately regulates gene expression and the chromatin landscape during neurogenesis

The proneural transcription factor Ascl1 coordinates gene expression in both proliferating and differentiating progenitors along the neuronal lineage. Here, we used a cellular model of neurogenesis to investigate how Ascl1 interacts with the chromatin landscape to regulate gene expression when promoting neuronal differentiation. We find that Ascl1 binding occurs mostly at distal enhancers and is associated with activation of gene transcription. Surprisingly, the accessibility of Ascl1 to its binding sites in neural stem/progenitor cells remains largely unchanged throughout their differentiation, as Ascl1 targets regions of both readily accessible and closed chromatin in proliferating cells. Moreover, binding of Ascl1 often precedes an increase in chromatin accessibility and the appearance of new regions of open chromatin, associated with de novo gene expression during differentiation. Our results reveal a function of Ascl1 in promoting chromatin accessibility during neurogenesis, linking the chromatin landscape at Ascl1 target regions with the temporal progression of its transcriptional program.

Relationships Among Murine CD11chigh Dendritic Cell Subsets as Revealed by Baseline Gene Expression Patterns

The functional relationships and properties of different subtypes of dendritic cells (DC) remain largely undefined. To better characterize these cells, we used global gene analysis to determine gene expression patterns among murine CD11c(high) DC subsets. CD4(+), CD8alpha(+), and CD8alpha(-) CD4(-) (double negative (DN)) DC were purified from spleens of normal C57/BL6 mice and analyzed using Affymetrix microarrays. The CD4(+) and CD8alpha(+) DC subsets showed distinct basal expression profiles differing by >200 individual genes. These included known DC subset markers as well as previously unrecognized, differentially expressed CD Ags such as CD1d, CD5, CD22, and CD72. Flow cytometric analysis confirmed differential expression in nine of nine cases, thereby validating the microarray analysis. Interestingly, the microarray expression profiles for DN cells strongly resembled those of CD4(+) DC, differing from them by <25 genes. This suggests that CD4(+) and DN DC are closely related phylogenetically, whereas CD8alpha(+) DC represent a more distant lineage, supporting the historical distinction between CD8alpha(+) and CD8alpha(-) DC. However, staining patterns revealed that in contrast to CD4(+) DC, the DN subset is heterogeneous and comprises at least two subpopulations. Gene Ontology and literature mining analyses of genes expressed differentially among DC subsets indicated strong associations with immune response parameters as well as cell differentiation and signaling. Such associations offer clues to possible unique functions of the CD11c(high) DC subsets that to date have been difficult to define as rigid distinctions.

Bmp4 regulates chick Ebf2 and Ebf3 gene expression in somite development

The chick Early B-cell Factor-2 and 3 (cEbf2 and cEbf3) genes are members of EBF family of helix loop helix transcription factors. The expression, regulation and importance of these genes have been extensively studied in lymphatic, nervous and muscular tissues. Recently, a new role for some members of EBF in bone development has been investigated. However, the expression profile and regulation in the axial skeleton precursor, the somite, have yet to be elucidated. Therefore, this study was aimed to investigate the expression and regulation of cEbf2 and cEbf3 genes in the developing chick embryo somite from HH4 to HH28. The spatiotemporal expression study revealed predominant localization of cEbf2 and cEbf3 in the lateral sclerotomal domains and later around vertebral cartilage anlagen of the arch and the proximal rib. Subsequently, microsurgeries, ectopic gene expression experiments were performed to analyze which tissues and factors regulate cEbf2 and cEbf3 expression. Lateral barriers experiments indicated the necessity for lateral signal(s) in the regulation of cEbf2 and cEbf3 genes. Results from tissue manipulations and ectopic gene expression experiments indicate that lateral plate-derived Bmp4 signals are necessary for the initiation and maintenance of cEbf2 and cEbf3 genes in somites. In conclusion, cEbf2 and cEbf3 genes are considered as lateral sclerotome markers which their expression is regulated by Bmp4 signals from the lateral plate mesoderm.

Changes in gene expression induced by H(2)O(2) in cardiac myocytes.

Oxidative stress induces cardiac myocyte apoptosis. At least some effects are probably mediated through changes in gene expression. Using Affymetrix arrays, we examined the changes in gene expression induced by H(2)O(2) (0.04, 0.1, and 0.2mM; 2 and 4h) in rat neonatal ventricular myocytes. Changes in selected upregulated genes were confirmed by ratiometric RT-PCR. p21(Cip1/Waf1) was one of the only two genes upregulated in all conditions studied. Of the heat shock proteins, only Hsp70/70.1 was induced by H(2)O(2) with no change in the expression of Hsp25, Hsp60 or Hsp90. Heme oxygenase 1 was also potently upregulated, but not heme oxygenases 2 or 3. Of the intercellular adhesion proteins, syndecan-1 was significantly upregulated in response to H(2)O(2), with little change in the expression of other syndecans and no change in expression of any of the integrins studied. Thus, oxidative stress, exemplified by H(2)O(2), selectively promotes the expression of specific gene family members.

Using U0126 to dissect the role of the extracellular signal-regulated kinase 1/2 (ERK1/2) cascade in the regulation of gene expression by endothelin-1 in cardiac myocytes

The hypertrophic agonist endothelin-1 rapidly but transiently activates the extracellular signal-regulated kinase 1/2 (ERK1/2) cascade (and other signalling pathways) in cardiac myocytes, but the events linking this to hypertrophy are not understood. Using Affymetrix rat U34A microarrays, we identified the short-term (2-4 h) changes in gene expression induced in neonatal myocytes by endothelin-1 alone or in combination with the ERK1/2 cascade inhibitor, U0126. Expression of 15 genes was significantly changed by U0126 alone, and expression of an additional 78 genes was significantly changed by endothelin-1. Of the genes upregulated by U0126, four are classically induced through the aryl hydrocarbon receptor (AhR) by dioxins suggesting that U0126 activates the xenobiotic response element in cardiac myocytes potentially independently of effects on ERK1/2 signalling. The 78 genes showing altered expression with endothelin-1 formed five clusters: (i) three clusters showing upregulation by endothelin-1 according to time course (4 h > 2 h; 2 h > 4 h; 2 h approximately 4 h) with at least partial inhibition by U0126; (ii) a cluster of 11 genes upregulated by endothelin-1 but unaffected by U0126 suggesting regulation through signalling pathways other than ERK1/2; (iii) a cluster of six genes downregulated by endothelin-1 with attenuation by U0126. Thus, U0126 apparently activates the AhR in cardiac myocytes (which must be taken into account in protracted studies), but careful analysis allows identification of genes potentially regulated acutely via the ERK1/2 cascade. Our data suggest that the majority of changes in gene expression induced by endothelin-1 are mediated by the ERK1/2 cascade.