Improving DNA array data quality by minimising 'neighbourhood' effects

Gene expression studies using cDNA arrays require robust and sensitive detection methods. Being extremely sensitive, radioactive detection suffers from the influence of signals positioned in each other's vicinity, the 'neighbourhood' effect. This limits the gene density of arrays and the quality of the results obtained. We have investigated the quantitative influence of different parameters on the 'neighbourhood' effect. By using a model experimental system, we could show that the effect is linear and depends only on the intensity of the hybridisation signal. We identified a common factor that can describe the influence of the neighbour spots based on their intensities. This factor is <1%, but it has to be taken into account if a high dynamic range of gene expression is to be detected. We could also derive the factor, although with less precision, from comparison of duplicate spots on arrays of 4565 different clones and replication of the hybridisation experiments. The calculated coefficient applied to our actual experimental results not only revealed previously undetected tissue or cell-specific expression differences, but also increased the dynamic range of detection. It thus provides a relatively simple way of improving DNA array data quality with few experimental modifications.     

Functional Similarities between Pigeon ‘Milk’ and Mammalian Milk: Induction of Immune Gene Expression and Modification of the Microbiota

Pigeon ‘milk’ and mammalian milk have functional similarities in terms of nutritional benefit and delivery of immunoglobulins to the young. Mammalian milk has been clearly shown to aid in the development of the immune system and microbiota of the young, but similar effects have not yet been attributed to pigeon ‘milk’. Therefore, using a chicken model, we investigated the effect of pigeon ‘milk’ on immune gene expression in the Gut Associated Lymphoid Tissue (GALT) and on the composition of the caecal microbiota. Chickens fed pigeon ‘milk’ had a faster rate of growth and a better feed conversion ratio than control chickens. There was significantly enhanced expression of immune-related gene pathways and interferon-stimulated genes in the GALT of pigeon ‘milk’-fed chickens. These pathways include the innate immune response, regulation of cytokine production and regulation of B cell activation and proliferation. The caecal microbiota of pigeon ‘milk’-fed chickens was significantly more diverse than control chickens, and appears to be affected by prebiotics in pigeon ‘milk’, as well as being directly seeded by bacteria present in pigeon ‘milk’. Our results demonstrate that pigeon ‘milk’ has further modes of action which make it functionally similar to mammalian milk. We hypothesise that pigeon ‘lactation’ and mammalian lactation evolved independently but resulted in similarly functional products.     

An engineered mammalian band-pass network

Gene expression circuitries, which enable cells to detect precise levels within a morphogen concentration gradient, have a pivotal impact on biological processes such as embryonic pattern formation, paracrine and autocrine signalling, and cellular migration. We present the rational synthesis of a synthetic genetic circuit exhibiting band-pass detection characteristics. The components, involving multiply linked mammalian trans-activator and -repressor control systems, were selected and fine-tuned to enable the detection of ‘low-threshold’ morphogen (tetracycline) concentrations, in which target gene expression was triggered, and a ‘high-threshold’ concentration, in which expression was muted. In silico predictions and supporting experimental findings indicated that the key criterion for functional band-pass detection was the matching of componentry that enabled sufficient separation of the low and high threshold points. Using the circuitry together with a fluorescence-encoded target gene, mammalian cells were genetically engineered to be capable of forming a band-like pattern of differentiation in response to a tetracycline chemical gradient. Synthetic gene networks designed to emulate naturally occurring gene behaviours provide not only insight into biological processes, but may also foster progress in future tissue engineering, gene therapy and biosensing applications.     

Plastid Proteomic Analysis in Tomato Fruit Development

To better understand the mechanism of plastid differentiation from chloroplast to chromoplast, we examined proteome and plastid changes over four distinct developmental stages of ‘Micro-Tom’ fruit. Additionally, to discover more about the relationship between fruit color and plastid differentiation, we also analyzed and compared ‘Micro-Tom’ results with those from two other varieties, ‘Black’ and ‘White Beauty’. We confirmed that proteins related to photosynthesis remain through the orange maturity stage of ‘Micro-Tom’, and also learned that thylakoids no longer exist at this stage. These results suggest that at a minimum there are changes in plastid morphology occurring before all related proteins change. We also compared ‘Micro-Tom’ fruits with ‘Black’ and ‘White Beauty’ using two-dimensional gel electrophoresis. We found a decrease of CHRC (plastid-lipid-associated protein) and HrBP1 (harpin binding protein-1) in the ‘Black’ and ‘White Beauty’ varieties. CHRC is involved in carotenoid accumulation and stabilization. HrBP1 in Arabidopsis has a sequence similar to proteins in the PAP/fibrillin family. These proteins have characteristics and functions similar to lipocalin, an example of which is the transport of hydrophobic molecules. We detected spots of TIL (temperature-induced lipocalin) in 2D-PAGE results, however the number of spots and their isoelectric points differed between ‘Micro-Tom’ and ‘Black’/‘White Beauty’. Lipocalin has various functions including those related to environmental stress response, apoptosis induction, membrane formation and fixation, regulation of immune response, cell growth, and metabolism adjustment. Lipocalin related proteins such as TIL and HrBP1 could be related to the accumulation of carotenoids, fruit color and the differentiation of chromoplast.     

The Expression of Pre- and Postcopulatory Sexually Selected Traits Reflects Levels of Dietary Stress in Guppies

Environmental and ecological conditions can shape the evolution of life history traits in many animals. Among such factors, food or nutrition availability can play an important evolutionary role in moderating an animal''s life history traits, particularly sexually selected traits. Here, we test whether diet quantity and/or composition in the form of omega-3 long chain polyunsaturated fatty acids (here termed ‘n3LC’) influence the expression of pre- and postcopulatory traits in the guppy (Poecilia reticulata), a livebearing poeciliid fish. We assigned males haphazardly to one of two experimental diets supplemented with n3LC, and each of these diet treatments was further divided into two diet ‘quantity’ treatments. Our experimental design therefore explored the main and interacting effects of two factors (n3LC content and diet quantity) on the expression of precopulatory (sexual behaviour and sexual ornamentation, including the size, number and spectral properties of colour spots) and postcopulatory (the velocity, viability, number and length of sperm) sexually selected traits. Our study revealed that diet quantity had significant effects on most of the pre- and postcopulatory traits, while n3LC manipulation had a significant effect on sperm traits and in particular on sperm viability. Our analyses also revealed interacting effects of diet quantity and n3LC levels on courtship displays, and the area of orange and iridescent colour spots in the males’ colour patterns. We also confirmed that our dietary manipulations of n3LC resulted in the differential uptake of n3LC in body and testes tissues in the different n3LC groups. This study reveals the effects of diet quantity and n3LC on behavioural, ornamental and ejaculate traits in P. reticulata and underscores the likely role that diet plays in maintaining the high variability in these condition-dependent sexual traits.     

Short- and long-range cis interactions between integrated HPV genomes and cellular chromatin dysregulate host gene expression in early cervical carcinogenesis

Development of cervical cancer is directly associated with integration of human papillomavirus (HPV) genomes into host chromosomes and subsequent modulation of HPV oncogene expression, which correlates with multi-layered epigenetic changes at the integrated HPV genomes. However, the process of integration itself and dysregulation of host gene expression at sites of integration in our model of HPV16 integrant clone natural selection has remained enigmatic. We now show, using a state-of-the-art ‘HPV integrated site capture’ (HISC) technique, that integration likely occurs through microhomology-mediated repair (MHMR) mechanisms via either a direct process, resulting in host sequence deletion (in our case, partially homozygously) or via a ‘looping’ mechanism by which flanking host regions become amplified. Furthermore, using our ‘HPV16-specific Region Capture Hi-C’ technique, we have determined that chromatin interactions between the integrated virus genome and host chromosomes, both at short- (<500 kbp) and long-range (>500 kbp), appear to drive local host gene dysregulation through the disruption of host:host interactions within (but not exceeding) host structures known as topologically associating domains (TADs). This mechanism of HPV-induced host gene expression modulation indicates that integration of virus genomes near to or within a ‘cancer-causing gene’ is not essential to influence their expression and that these modifications to genome interactions could have a major role in selection of HPV integrants at the early stage of cervical neoplastic progression.     

A Linear Mixed Model Spline Framework for Analysing Time Course ‘Omics’ Data

Time course ‘omics’ experiments are becoming increasingly important to study system-wide dynamic regulation. Despite their high information content, analysis remains challenging. ‘Omics’ technologies capture quantitative measurements on tens of thousands of molecules. Therefore, in a time course ‘omics’ experiment molecules are measured for multiple subjects over multiple time points. This results in a large, high-dimensional dataset, which requires computationally efficient approaches for statistical analysis. Moreover, methods need to be able to handle missing values and various levels of noise. We present a novel, robust and powerful framework to analyze time course ‘omics’ data that consists of three stages: quality assessment and filtering, profile modelling, and analysis. The first step consists of removing molecules for which expression or abundance is highly variable over time. The second step models each molecular expression profile in a linear mixed model framework which takes into account subject-specific variability. The best model is selected through a serial model selection approach and results in dimension reduction of the time course data. The final step includes two types of analysis of the modelled trajectories, namely, clustering analysis to identify groups of correlated profiles over time, and differential expression analysis to identify profiles which differ over time and/or between treatment groups. Through simulation studies we demonstrate the high sensitivity and specificity of our approach for differential expression analysis. We then illustrate how our framework can bring novel insights on two time course ‘omics’ studies in breast cancer and kidney rejection. The methods are publicly available, implemented in the R CRAN package lmms.     

DNA probes on beads arrayed in a capillary, ‘Bead-array’, exhibited high hybridization performance

A DNA analysis platform called ‘Bead-array’ is presented and its features when used in hybridization detection are shown. In ‘Bead-array’, beads of 100-µm diameter are lined in a determined order in a capillary. Each bead is conjugated with DNA probes, and can be identified by its order in the capillary. This probe array is easily produced by just arraying beads conjugated with probes into the capillary in a fixed order. The hybridization is also easily completed by introducing samples (1–300 µl) into the capillary with reciprocal flow. For hybridization detection, as little as 1 amol of fluorescent-labeled oligo DNA was detected. The hybridization reaction was completed in 1 min irrespective of the amount of target DNA. When the number of target molecules was smaller than that of probe molecules on the bead, 10 fmol, almost all targets were captured on the bead. ‘Bead-array’ enables reliable and reproducible measurement of the target quantity. This rapid and sensitive platform seems very promising for various genetic testing tasks.     

Differential Responses of CO2 Assimilation,Carbohydrate Allocation and Gene Expression to NaCl Stress in Perennial Ryegrass with Different Salt Tolerance

Little is known about the effects of NaCl stress on perennial ryegrass (Lolium perenne L.) photosynthesis and carbohydrate flux. The objective of this study was to understand the carbohydrate metabolism and identify the gene expression affected by salinity stress. Seventy-four days old seedlings of two perennial ryegrass accessions (salt-sensitive ‘PI 538976’ and salt-tolerant ‘Overdrive’) were subjected to three levels of salinity stress for 5 days. Turf quality in all tissues (leaves, stems and roots) of both grass accessions negatively and significantly correlated with GFS (Glu+Fru+Suc) content, except for ‘Overdrive’ stems. Relative growth rate (RGR) in leaves negatively and significantly correlated with GFS content in ‘Overdrive’ (P<0.01) and ‘PI 538976’ (P<0.05) under salt stress. ‘Overdrive’ had higher CO₂ assimilation and F_v/F_m than ‘PI 538976’. Intercellular CO₂ concentration, however, was higher in ‘PI 538976’ treated with 400 mM NaCl relative to that with 200 mM NaCl. GFS content negatively and significantly correlated with RGR in ‘Overdrive’ and ‘PI 538976’ leaves and in ‘PI 538976’ stems and roots under salt stress. In leaves, carbohydrate allocation negatively and significantly correlated with RGR (r² = 0.83, P<0.01) and turf quality (r² = 0.88, P<0.01) in salt-tolerant ‘Overdrive’, however, the opposite trend for salt-sensitive ‘PI 538976’ (r² = 0.71, P<0.05 for RGR; r² = 0.62, P>0.05 for turf quality). A greater up-regulation in the expression of SPS, SS, SI, 6-SFT gene was observed in ‘Overdrive’ than ‘PI 538976’. A higher level of SPS and SS expression in leaves was found in ‘PI 538976’ relative to ‘Overdrive’. Accumulation of hexoses in roots, stems and leaves can induce a feedback repression to photosynthesis in salt-stressed perennial ryegrass and the salt tolerance may be changed with the carbohydrate allocation in leaves and stems.     

Human and mouse oligonucleotide-based array CGH

Array-based comparative genomic hybridization is a high resolution method for measuring chromosomal copy number changes. Here we present a validated protocol using in-house spotted oligonucleotide libraries for array comparative genomic hybridization (CGH). This oligo array CGH platform yields reproducible results and is capable of detecting single copy gains, multi-copy amplifications as well as homozygous and heterozygous deletions as small as 100 kb with high resolution. A human oligonucleotide library was printed on amine binding slides. Arrays were hybridized using a hybstation and analysed using BlueFuse feature extraction software, with >95% of spots passing quality control. The protocol allows as little as 300 ng of input DNA and a 90% reduction of Cot-1 DNA without compromising quality. High quality results have also been obtained with DNA from archival tissue. Finally, in addition to human oligo arrays, we have applied the protocol successfully to mouse oligo arrays. We believe that this oligo-based platform using ‘off-the-shelf’ oligo libraries provides an easy accessible alternative to BAC arrays for CGH, which is cost-effective, available at high resolution and easily implemented for any sequenced organism without compromising the quality of the results.     

Hypervariable genes—experimental error or hidden dynamics

In a homogeneous group of samples, not all genes of high variability stem from experimental errors in microarray experiments. These expression variations can be attributed to many factors including natural biological oscillations or metabolic processes. The behavior of these genes can tease out important clues about naturally occurring dynamic processes in the organism or experimental system under study. We developed a statistical procedure for the selection of genes with high variability denoted hypervariable (HV) genes. After the exclusion of low expressed genes and a stabilizing log-transformation, the majority of genes have comparable residual variability. Based on an F-test, HV genes are selected as having a statistically significant difference from the majority of variability stabilized genes measured by the ‘reference group’. A novel F-test clustering technique, further noted as ‘F-means clustering’, groups HV genes with similar variability patterns, presumably from their participation in a common dynamic biological process. F-means clustering establishes, for the first time, groups of co-expressed HV genes and is illustrated with microarray data from patients with juvenile rheumatoid arthritis and healthy controls.     

Breaking an Epigenetic Chromatin Switch: Curious Features of Hysteresis in Saccharomyces cerevisiae Telomeric Silencing

In addition to gene network switches, local epigenetic modifications to DNA and histones play an important role in all-or-none cellular decision-making. Here, we study the dynamical design of a well-characterized epigenetic chromatin switch: the yeast SIR system, in order to understand the origin of the stability of epigenetic states. We study hysteresis in this system by perturbing it with a histone deacetylase inhibitor. We find that SIR silencing has many characteristics of a non-linear bistable system, as observed in conventional genetic switches, which are based on activities of a few promoters affecting each other through the abundance of their gene products. Quite remarkably, our experiments in yeast telomeric silencing show a very distinctive pattern when it comes to the transition from bistability to monostability. In particular, the loss of the stable silenced state, upon increasing the inhibitor concentration, does not seem to show the expected saddle node behavior, instead looking like a supercritical pitchfork bifurcation. In other words, the ‘off’ state merges with the ‘on’ state at a threshold concentration leading to a single state, as opposed to the two states remaining distinct up to the threshold and exhibiting a discontinuous jump from the ‘off’ to the ‘on’ state. We argue that this is an inevitable consequence of silenced and active regions coexisting with dynamic domain boundaries. The experimental observations in our study therefore have broad implications for the understanding of chromatin silencing in yeast and beyond.     

Tethering-facilitated DNA ‘opening’ and complementary roles of β-hairpin motifs in the Rad4/XPC DNA damage sensor protein

XPC/Rad4 initiates eukaryotic nucleotide excision repair on structurally diverse helix-destabilizing/distorting DNA lesions by selectively ‘opening’ these sites while rapidly diffusing along undamaged DNA. Previous structural studies showed that Rad4, when tethered to DNA, could also open undamaged DNA, suggesting a ‘kinetic gating’ mechanism whereby lesion discrimination relied on efficient opening versus diffusion. However, solution studies in support of such a mechanism were lacking and how ‘opening’ is brought about remained unclear. Here, we present crystal structures and fluorescence-based conformational analyses on tethered complexes, showing that Rad4 can indeed ‘open’ undamaged DNA in solution and that such ‘opening’ can largely occur without one or the other of the β-hairpin motifs in the BHD2 or BHD3 domains. Notably, the Rad4-bound ‘open’ DNA adopts multiple conformations in solution notwithstanding the DNA’s original structure or the β-hairpins. Molecular dynamics simulations reveal compensatory roles of the β-hairpins, which may render robustness in dealing with and opening diverse lesions. Our study showcases how fluorescence-based studies can be used to obtain information complementary to ensemble structural studies. The tethering-facilitated DNA ‘opening’ of undamaged sites and the dynamic nature of ‘open’ DNA may shed light on how the protein functions within and beyond nucleotide excision repair in cells.     

From DNA biosensors to gene chips

Wide-scale DNA testing requires the development of small, fast and easy-to-use devices. This article describes the preparation, operation and applications of biosensors and gene chips, which provide fast, sensitive and selective detection of DNA hybridization. Various new strategies for DNA biosensors and gene chips are examined, along with recent trends and future directions. The integration of hybridization detection schemes with the sample preparation process in a ‘Lab-on-a-Chip’ format is also covered. While the use of DNA biosensors and gene chips is at an early stage, such devices are expected to have an enormous effect on future DNA diagnostics.