Optimising the analysis of transcript data using high density oligonucleotide arrays and genomic DNA-based probe selection

Background: Affymetrix GeneChip arrays are widely used for transcriptomic studies in a diverse range of species. Each gene is represented on a GeneChip array by a probe- set, consisting of up to 16 probe-pairs. Signal intensities across probe- pairs within a probe-set vary in part due to different physical hybridisation characteristics of individual probes with their target labelled transcripts. We have previously developed a technique to study the transcriptomes of heterologous species based on hybridising genomic DNA (gDNA) to a GeneChip array designed for a different species, and subsequently using only those probes with good homology. Results: Here we have investigated the effects of hybridising homologous species gDNA to study the transcriptomes of species for which the arrays have been designed. Genomic DNA from Arabidopsis thaliana and rice (Oryza sativa) were hybridised to the Affymetrix Arabidopsis ATH1 and Rice Genome GeneChip arrays respectively. Probe selection based on gDNA hybridisation intensity increased the number of genes identified as significantly differentially expressed in two published studies of Arabidopsis development, and optimised the analysis of technical replicates obtained from pooled samples of RNA from rice. Conclusion: This mixed physical and bioinformatics approach can be used to optimise estimates of gene expression when using GeneChip arrays.

Using genomic DNA-based probe-selection to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species

High-density oligonucleotide (oligo) arrays are a powerful tool for transcript profiling. Arrays based on GeneChip® technology are amongst the most widely used, although GeneChip® arrays are currently available for only a small number of plant and animal species. Thus, we have developed a method to improve the sensitivity of high-density oligonucleotide arrays when applied to heterologous species and tested the method by analysing the transcriptome of Brassica oleracea L., a species for which no GeneChip® array is available, using a GeneChip® array designed for Arabidopsis thaliana (L.) Heynh. Genomic DNA from B. oleracea was labelled and hybridised to the ATH1-121501 GeneChip® array. Arabidopsis thaliana probe-pairs that hybridised to the B. oleracea genomic DNA on the basis of the perfect-match (PM) probe signal were then selected for subsequent B. oleracea transcriptome analysis using a .cel file parser script to generate probe mask files. The transcriptional response of B. oleracea to a mineral nutrient (phosphorus; P) stress was quantified using probe mask files generated for a wide range of gDNA hybridisation intensity thresholds. An example probe mask file generated with a gDNA hybridisation intensity threshold of 400 removed > 68 % of the available PM probes from the analysis but retained >96 % of available A. thaliana probe-sets. Ninety-nine of these genes were then identified as significantly regulated under P stress in B. oleracea, including the homologues of P stress responsive genes in A. thaliana. Increasing the gDNA hybridisation intensity thresholds up to 500 for probe-selection increased the sensitivity of the GeneChip® array to detect regulation of gene expression in B. oleracea under P stress by up to 13-fold. Our open-source software to create probe mask files is freely available http://affymetrix.arabidopsis.info/xspecies/ webcite and may be used to facilitate transcriptomic analyses of a wide range of plant and animal species in the absence of custom arrays.

Genomic signals of migration and continuity in Britain before the Anglo-Saxons

The purported migrations that have formed the peoples of Britain have been the focus of generations of scholarly controversy. However, this has not benefited from direct analyses of ancient genomes. Here we report nine ancient genomes (~1 x) of individuals from northern Britain: seven from a Roman era York cemetery, bookended by earlier Iron-Age and later Anglo-Saxon burials. Six of the Roman genomes show affinity with modern British Celtic populations, particularly Welsh, but significantly diverge from populations from Yorkshire and other eastern English samples. They also show similarity with the earlier Iron-Age genome, suggesting population continuity, but differ from the later Anglo-Saxon genome. This pattern concords with profound impact of migrations in the Anglo-Saxon period. Strikingly, one Roman skeleton shows a clear signal of exogenous origin, with affinities pointing towards the Middle East, confirming the cosmopolitan character of the Empire, even at its northernmost fringes.

Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Early membrane estrogenic effects required for full expression of slower genomic actions in a nerve cell line.

Interpretations of steroid hormone actions as slow, nuclear, transcriptional events have frequently been seen as competing against inferences of rapid membrane actions. We have discovered conditions where membrane-limited effects potentiate later transcriptional actions in a nerve cell line. Making use of a two-pulse hormonal schedule in a transfection system, early and brief administration of conjugated, membrane-limited estradiol was necessary but not sufficient for full transcriptional potency of the second estrogen pulse. Efficacy of the first pulse depended on intact signal transduction pathways. Surprisingly, the actions of both pulses were blocked by a classical estrogen receptor (ER) antagonist. Thus, two different modes of steroid hormone action can synergize.

Integration of steroid hormone initiated membrane action to genomic function in the brain

Estrogen is a ligand for the estrogen receptor (ER), which on binding 17beta-estradiol, functions as a ligand-activated transcription factor and regulates the transcription of target genes. This is the slow genomic mode of action. However, rapid non-genomic actions of estrogen also exist at the cell membrane. Using a novel two-pulse paradigm in which the first pulse rapidly initiates non-genomic actions using a membrane-limited estrogen conjugate (E-BSA), while the second pulse promotes genomic transcription from a consensus estrogen response element (ERE), we have demonstrated that rapid actions of estrogen potentiate the slower transcriptional response from an ERE-reporter in neuroblastoma cells. Since rapid actions of estrogen activate kinases, we used selective inhibitors in the two-pulse paradigm to determine the intracellular signaling cascades important in such potentiation. Inhibition of protein kinase A (PKA), PKC, mitogen activated protein kinase (MAPK) or phosphatidylinositol 3-OH kinase (PI-3K) in the first pulse decreases potentiation of transcription. Also, our data with both dominant negative and constitutive mutants of Galpha subunits show that Galpha(q) initiates the rapid signaling cascade at the membrane in SK-N-BE(2)C neuroblastoma cells. We discuss two models of multiple kinase activation at the membrane Pulses of estrogen induce lordosis behavior in female rats. Infusion of E-BSA into the ventromedial hypothalamus followed by 17beta-estradiol in the second pulse could induce lordosis behavior, demonstrating the applicability of this paradigm in vivo. A model where non-genomic actions of estrogen couple to genomic actions unites both aspects of hormone action.

Calcium Flux in Neuroblastoma Cells Is a Coupling Mechanism between Non-Genomic and Genomic Modes of Estrogens

Estrogens have been demonstrated to rapidly modulate calcium levels in a variety of cell types. However, the significance of estrogen-mediated calcium flux in neuronal cells is largely unknown. The relative importance of intra- and extracellular sources of calcium in estrogenic effects on neurons is also not well understood. Previously, we have demonstrated that membrane-limited estrogens, such as E-BSA given before an administration of a 2-hour pulse of 17beta-estradiol (E(2)), can potentiate the transcription mediated by E(2) from a consensus estrogen response element (ERE)-driven reporter gene. Inhibitors to signal transduction cascades given along with E-BSA or E(2) demonstrated that calcium flux is important for E-BSA-mediated potentiation of transcription in a transiently transfected neuroblastoma cell line. In this report, we have used inhibitors to different voltage-gated calcium channels (VGCCs) and to intracellular store receptors along with E-BSA in the first pulse or with E(2) in the second pulse to investigate the relative importance of these channels to estrogen-mediated transcription. Neither L- nor P-type VGCCs seem to play a role in estrogen action in these cells; while N-type VGCCs are important in both the non-genomic and genomic modes of estrogen action. Specific inhibitors also showed that the ryanodine receptor and the inositol trisphosphate receptor are important to E-BSA-mediated transcriptional potentiation. This report provides evidence that while intracellular stores of calcium are required to couple non-genomic actions of estrogen initiated at the membrane to transcription in the nucleus, extracellular sources of calcium are also important in both non-genomic and genomic actions of estrogens. Copyright (c) 2005 S. Karger AG, Basel.

Estrogens and thyroid hormone: Non-genomic effects are coupled to transcription.

Estrogens and thyroid hormones are regulators of important diverse physiological processes such as reproduction, thermogenesis, neural development, neural differentiation and cardiovascular functions. Both are ligands for receptors in the nuclear receptor superfamily, which act as ligand-dependent transcription factors, regulating transcription. However, estrogens and thyroid hormones also rapidly (within minutes or seconds) activate kinase cascades and calcium increases, presumably initiated at the cell membrane. We discuss the relevance of both modes of hormone action, including the membrane estrogen receptor, to physiology, with particular reference to lordosis behavior. We first showed that estrogen restricted to the membrane can, in fact, lead to subsequent increases in transcription from a consensus estrogen response element-based reporter in the neuroblastoma cell line, SK-N-BE(2)C. Using a novel hormonal paradigm, we also showed that the activation of protein kinase A, protein kinase C, mitogen activated protein kinase and increases in calcium were important in the ability of the membrane-limited estrogen to potentiate transcription. We discuss the source of calcium important in transcriptional potentiation. Since estrogens and thyroid hormones have common effects on neuroprotection, cognition and mood, we also hypothesized that crosstalk could occur between the rapid actions of thyroid hormones and the genomic actions of estrogens. In neural cells, we showed that triiodothyronine acting rapidly via MAPK can increase transcription by the nuclear estrogen receptor ERa from a consensus estrogen response element, possibly by the phosphorylation of the ERa. Novel mechanisms that link signals initiated by hormones from the membrane to the nucleus are physiologically relevant and can achieve neuroendocrine integration

Non-genomic actions of estrogens and their interaction with genomic actions in the brain

Ligands for the nuclear receptor superfamily have at least two mechanisms of action: (a) classical transcriptional regulation of target genes (genomic mechanisms); and (b) non-genomic actions, which are initiated at the cell membrane, which could also impact transcription. Though transcriptional mechanisms are increasingly well understood, membrane-initiated actions of these ligands are incompletely understood. This has led to considerable debate over the physiological relevance of membrane-initiated actions of hormones versus genomic actions of hormones, with genomic actions predominating in the endocrine field. There is good evidence that the membrane-limited actions of hormones, particularly estrogens, involve the rapid activation of kinases and the release of calcium and that these are linked to physiologically relevant scenarios in the brain. We show evidence in this review, that membrane actions of estrogens, which activate these rapid signaling cascades, can also potentiate nuclear transcription in both the central nervous system and in non-neuronal cell lines. We present a theoretical scenario which can be used to understand this phenomenon. These signaling cascades may occur in parallel or in series but subsequently, converge at the modification of transcriptionally relevant molecules such as nuclear receptors and/or coactivators. In addition, other non-cognate hormones or neurotransmitters may also activate cascades to crosstalk with estrogen receptor-mediated transcription, though the relevance of this is less clear. The idea that coupling between membrane-initiated and genomic actions of hormones is a novel idea in neuroendocrinology and provides us with a unified view of hormone action in the central nervous system.

Membrane-initiated non-genomic signaling by estrogens in the hypothalamus: cross-talk with glucocorticoids with implications for behavior

The estrogen receptor and glucocorticoid receptor are members of the nuclear receptor superfamily that can signal using both non-genomic and genomic transcriptional modes. Though genomic modes of signaling have been well characterized and several behaviors attributed to this signaling mechanism, the physiological significance of non-genomic modes of signaling has not been well understood. This has partly been due to the controversy regarding the identity of the membrane ER (mER) or membrane GR (mGR) that may mediate rapid, non-genomic signaling and the downstream signaling cascades that may result as a consequence of steroid ligands binding the mER or the mGR. Both estrogens and glucocorticoids exert a number of actions on the hypothalamus, including feedback. This review focuses on the various candidates for the mER or mGR in the hypothalamus and the contribution of non-genomic signaling to classical hypothalamically driven behaviors and changes in neuronal morphology. It also attempts to categorize some of the possible functions of non-genomic signaling at both the cellular level and at the organismal level that are relevant for behavior, including some behaviors that are regulated by both estrogens and glucocorticoids in a potentially synergistic manner. Lastly, it attempts to show that steroid signaling via non-genomic modes may provide the organism with rapid behavioral responses to stimuli.

The use of fluorescein for labeling genomic probes in the checkerboard DNA-DNA hybridization method

Molecular methods that permit the simultaneous detection and quantification of a large number of microbial species are currently employed in the evaluation of complex ecosystems. The checkerboard DNA-DNA hybridization technique enables the simultaneous identification of distinct bacterial. species in a large number of dental samples. The original technique employed digoxigenin-labeled whole genomic DNA probes which were detected by chemiluminescence. In this study, we present an alternative protocol for labeling and detecting whole genomic DNA probes in the Checkerboard DNA-DNA hybridization method. Whole genomic DNA was extracted from five bacterial species and labeled with fluorescein. The fluorescein labeled whole genomic DNA probes were hybridized against whole genomic DNA or subgingival plaque samples in a checkerboard hybridization format, followed by chemiluminescent detection. Our results reveal that fluorescein is a viable and adequate alternative labeling reagent to be employed in the checkerboard DNA-DNA hybridization technique. (c) 2007 Elsevier GmbH. All rights reserved.

Genomic and nongenomic stimulatory effect of aldosterone on H(+)-ATPase in proximal S3 segments

Leite-Dellova DC, Malnic G, Mello-Aires M. Genomic and non-genomic stimulatory effect of aldosterone on H(+)-ATPase in proximal S3 segments. Am J Physiol Renal Physiol 300: F682-F691, 2011. First published December 29, 2010; doi:10.1152/ajprenal.00172.2010.-The genomic and nongenomic effects of aldosterone on the intracellular pH recovery rate (pHirr) via H(-)(+)ATPase and on cytosolic free calcium concentration ([Ca(2+)](i)) were investigated in isolated proximal S3 segments of rats during superfusion with an Na(+)-free solution, by using the fluorescent probes BCECF-AM and FLUO-4-AM, respectively. The pHirr, after cellular acidification with a NH(4)Cl pulse, was 0.064 +/- 0.003 pH units/min (n = 17/74) and was abolished with concanamycin. Aldosterone (10(-12), 10(-10),10(-8), or 10(-6) M with 1-h or 15- or 2-min preincubation) increased the pHirr. The baseline [Ca(2+)](i) was 103 +/- 2 nM (n = 58). After 1 min of aldosterone preincubation, there was a transient and dose-dependent increase in [Ca(2+)](i) and after 6-min preincubation there was a new increase in [Ca(2+)](i) that persisted after 1 h. Spironolactone [mineralocorticoid (MR) antagonist], actinomycin D, or cycloheximide did not affect the effects of aldosterone (15- or 2-min preincubation) on pHirr and on [Ca(2+)](i) but inhibited the effects of aldosterone (1-h preincubation) on these parameters. RU 486 [glucocorticoid (GR) antagonist] and dimethyl-BAPTA (Ca(2+) chelator) prevented the effect of aldosterone on both parameters. The data indicate a genomic (1 h, via MR) and a nongenomic action (15 or 2 min, probably via GR) on the H(+)-ATPase and on [Ca(2+)](i). The results are compatible with stimulation of the H(+)-ATPase by increases in [Ca(2+)](i) (at 10(-12)-10(-6) M aldosterone) and inhibition of the H(+)-ATPase by decreases in [Ca(2+)](i) (at 10(-12) or 10(-6) M aldosterone plus RU 486).

Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair

Recurrent submicroscopic genomic copy number changes are the result of nonallelic homologous recombination (NAHR). Nonrecurrent aberrations, however, can result from different nonexclusive recombination-repair mechanisms. We previously described small microduplications at Xq28 containing MECP2 in four male patients with a severe neurological phenotype. Here, we report on the fine-mapping and breakpoint analysis of 16 unique microduplications. The size of the overlapping copy number changes varies between 0.3 and 2.3 Mb, and FISH analysis on three patients demonstrated a tandem orientation. Although eight of the 32 breakpoint regions coincide with low-copy repeats, none of the duplications are the result of NAHR. Bioinformatics analysis of the breakpoint regions demonstrated a 2.5-fold higher frequency of Alu interspersed repeats as compared with control regions, as well as a very high GC content (53%). Unexpectedly, we obtained the junction in only one patient by long-range PCR, which revealed nonhomologous end joining as the mechanism. Breakpoint analysis in two other patients by inverse PCR and subsequent array comparative genomic hybridization analysis demonstrated the presence of a second duplicated region more telomeric at Xq28, of which one copy was inserted in between the duplicated MECP2 regions. These data suggest a two-step mechanism in which part of Xq28 is first inserted near the MECP2 locus, followed by breakage-induced replication with strand invasion of the normal sister chromatid. Our results indicate that the mechanism by which copy number changes occur in regions with a complex genomic architecture can yield complex rearrangements.

Genomic imbalances associated with mullerian aplasia

Background: Aplasia of the mullerian ducts leads to absence of the uterine corpus, uterine cervix, and upper (superior) vagina. Patients with mullerian aplasia (MA) often exhibit additional clinical features such as renal, vertebral and cardiac defects. A number of different syndromes have been associated with MA, and in most cases its aetiology remains poorly understood. Objective and methods: 14 syndromic patients with MA and 46, XX G-banded karyotype were screened for DNA copy number changes by similar to 1 Mb whole genome bacterial artificial chromosome (BAC) array based comparative genomic hybridisation (CGH). The detected alterations were validated by an independent method and further mapped by high resolution oligo-arrays. Results: Submicroscopic genomic imbalances affecting the 1q21.1, 17q12, 22q11.21, and Xq21.31 chromosome regions were detected in four probands. Presence of the alterations in the normal mother of one patient suggests incomplete penetrance and/or variable expressivity. Conclusion: 4 of the 14 patients (29%) were found to have cryptic genomic alterations. The imbalances on 22q11.21 support recent findings by us and others that alterations in this chromosome region may result in impairment of mullerian duct development. The remaining imbalances indicate involvement of previously unknown chromosome regions in MA, and point specifically to LHX1 and KLHL4 as candidate genes.

Genomic Analysis of Wild Tomato Introgressions Determining Metabolism- and Yield-Associated Traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.