Whole-genome Tiling Microarrays Research Articles

Transcriptome profile in Drosophila Kc and S2 embryonic cell lines.

Drosophila melanogaster cell lines are an important resource for a range of studies spanning genomics, molecular genetics, and cell biology. Amongst these valuable lines are Kc167 (Kc) and Schneider 2 (S2) cells, which were originally isolated in the late 1960s from embryonic sources and have been used extensively to investigate a broad spectrum of biological activities including cell-cell signaling and immune system function. Whole-genome tiling microarray analysis of total RNA from these two cell types was performed as part of the modENCODE project over a decade ago and revealed that they share a number of gene expression features. Here, we expand on these earlier studies by using deep-coverage RNA-sequencing approaches to investigate the transcriptional profile in Kc and S2 cells in detail. Comparison of the transcriptomes reveals that ∼75% of the 13,919 annotated genes are expressed at a detectable level in at least one of the cell lines, with the majority of these genes expressed at high levels in both cell lines. Despite the overall similarity of the transcriptional landscape in the two cell types, 2,588 differentially expressed genes are identified. Many of the genes with the largest fold change are known only by their "CG" designations, indicating that the molecular control of Kc and S2 cell identity may be regulated in part by a cohort of relatively uncharacterized genes. Our data also indicate that both cell lines have distinct hemocyte-like identities, but share active signaling pathways and express a number of genes in the network responsible for dorsal-ventral patterning of the early embryo.

Thermodynamically optimal whole-genome tiling microarray design and validation.

BackgroundMicroarray is an efficient apparatus to interrogate the whole transcriptome of species. Microarray can be designed according to annotated gene sets, but the resulted microarrays cannot be used to identify novel transcripts and this design method is not applicable to unannotated species. Alternatively, a whole-genome tiling microarray can be designed using only genomic sequences without gene annotations, and it can be used to detect novel RNA transcripts as well as known genes. The difficulty with tiling microarray design lies in the tradeoff between probe-specificity and coverage of the genome. Sequence comparison methods based on BLAST or similar software are commonly employed in microarray design, but they cannot precisely determine the subtle thermodynamic competition between probe targets and partially matched probe nontargets during hybridizations.FindingsUsing the whole-genome thermodynamic analysis software PICKY to design tiling microarrays, we can achieve maximum whole-genome coverage allowable under the thermodynamic constraints of each target genome. The resulted tiling microarrays are thermodynamically optimal in the sense that all selected probes share the same melting temperature separation range between their targets and closest nontargets, and no additional probes can be added without violating the specificity of the microarray to the target genome.ConclusionsThis new design method was used to create two whole-genome tiling microarrays for Escherichia coli MG1655 and Agrobacterium tumefaciens C58 and the experiment results validated the design.Electronic supplementary materialThe online version of this article (doi:10.1186/s13104-016-2113-4) contains supplementary material, which is available to authorized users.

Chromatin remodeling in the aging genome of Drosophila

Chromatin structure affects the accessibility of DNA to transcription, repair, and replication. Changes in chromatin structure occur during development, but less is known about changes during aging. We examined the state of chromatin structure and its effect on gene expression during aging in Drosophila at the whole genome and cellular level using whole-genome tiling microarrays of activation and repressive chromatin marks, whole-genome transcriptional microarrays and single-cell immunohistochemistry. We found dramatic reorganization of chromosomal regions with age. Mapping of H3K9me3 and HP1 signals to fly chromosomes reveals in young flies the expected high enrichment in the pericentric regions, the 4th chromosome, and islands of facultative heterochromatin dispersed throughout the genome. With age, there is a striking reduction in this enrichment resulting in a nearly equivalent level of H3K9me3 and HP1 in the pericentric regions, the 4th chromosome, facultative heterochromatin, and euchromatin. These extensive changes in repressive chromatin marks are associated with alterations in age-related gene expression. Large-scale changes in repressive marks with age are further substantiated by single-cell immunohistochemistry that shows changes in nuclear distribution of H3K9me3 and HP1 marks with age. Such epigenetic changes are expected to directly or indirectly impinge upon important cellular functions such as gene expression, DNA repair, and DNA replication. The combination of genome-wide approaches such as whole-genome chromatin immunoprecipitation and transcriptional studies in conjunction with single-cell immunohistochemistry as shown here provide a first step toward defining how changes in chromatin may contribute to the process of aging in metazoans.

Comprehensive Identification of Salmonella enterica Serovar Typhimurium Genes Required for Infection of BALB/c Mice

Genes required for infection of mice by Salmonella Typhimurium can be identified by the interrogation of random transposon mutant libraries for mutants that cannot survive in vivo. Inactivation of such genes produces attenuated S. Typhimurium strains that have potential for use as live attenuated vaccines. A quantitative screen, Transposon Mediated Differential Hybridisation (TMDH), has been developed that identifies those members of a large library of transposon mutants that are attenuated. TMDH employs custom transposons with outward-facing T7 and SP6 promoters. Fluorescently-labelled transcripts from the promoters are hybridised to whole-genome tiling microarrays, to allow the position of the transposon insertions to be determined. Comparison of microarray data from the mutant library grown in vitro (input) with equivalent data produced after passage of the library through mice (output) enables an attenuation score to be determined for each transposon mutant. These scores are significantly correlated with bacterial counts obtained during infection of mice using mutants with individual defined deletions of the same genes. Defined deletion mutants of several novel targets identified in the TMDH screen are effective live vaccines.

Evidence for Antisense Transcription Associated with MicroRNA Target mRNAs in Arabidopsis

Antisense transcription is a pervasive phenomenon, but its source and functional significance is largely unknown. We took an expression-based approach to explore microRNA (miRNA)-related antisense transcription by computational analyses of published whole-genome tiling microarray transcriptome and deep sequencing small RNA (smRNA) data. Statistical support for greater abundance of antisense transcription signatures and smRNAs was observed for miRNA targets than for paralogous genes with no miRNA cleavage site. Antisense smRNAs were also found associated with MIRNA genes. This suggests that miRNA-associated “transitivity” (production of small interfering RNAs through antisense transcription) is more common than previously reported. High-resolution (3 nt) custom tiling microarray transcriptome analysis was performed with probes 400 bp 5′ upstream and 3′ downstream of the miRNA cleavage sites (direction relative to the mRNA) for 22 select miRNA target genes. We hybridized RNAs labeled from the smRNA pathway mutants, including hen1-1, dcl1-7, hyl1-2, rdr6-15, and sgs3-14. Results showed that antisense transcripts associated with miRNA targets were mainly elevated in hen1-1 and sgs3-14 to a lesser extent, and somewhat reduced in dcl11-7, hyl11-2, or rdr6-15 mutants. This was corroborated by semi-quantitative reverse transcription PCR; however, a direct correlation of antisense transcript abundance in MIR164 gene knockouts was not observed. Our overall analysis reveals a more widespread role for miRNA-associated transitivity with implications for functions of antisense transcription in gene regulation. HEN1 and SGS3 may be links for miRNA target entry into different RNA processing pathways.

Genome-wide identification of binding sites for Kaposi's sarcoma-associated herpesvirus lytic switch protein, RTA

Kaposi's sarcoma-associated herpesvirus (KSHV) replication and transcription activator (RTA) encoded by ORF50 is a lytic switch protein for viral reactivation from latency. The expression of RTA activates the expression of downstream viral genes, and is necessary for triggering the full viral lytic program. Using chromatin immunoprecipitation assay coupled with a KSHV whole-genome tiling microarray (ChIP-on-chip) approach, we identified a set of 19 RTA binding sites in the KSHV genome in a KSHV-infected cell line BCBL-1. These binding sites are located in the regions of promoters, introns, or exons of KSHV genes including ORF8, ORFK4.1, ORFK5, PAN, ORF16, ORF29, ORF45, ORF50, ORFK8, ORFK10.1, ORF59, ORFK12, ORF71/72, ORFK14/ORF74, and ORFK15, the two origins of lytic replication OriLyt-L and OriLyt-R, and the microRNA cluster. We confirmed these RTA binding sites by ChIP and quantitative real-time PCR. We further mapped the RTA binding site in the first intron of the ORFK15 gene, and determined that it is RTA-responsive. The ORFK15 RTA binding sequence TTCCAGGAA TTCCTGGAA consists of a palindromic structure of two tandem repeats, of which each itself is also an imperfect inverted repeat. Reporter assay and electrophoretic mobility shift assay confirmed the binding of the RTA protein to this sequence in vitro. Sequence alignment with other RTA binding sites identified the RTA consensus binding motif as TTCCAGGAT(N) 0–16TTCCTGGGA. Interestingly, most of the identified RTA binding sites contain only half or part of this RTA binding motif. These results suggest the complexity of RTA binding in vivo, and the involvement of other cellular or viral transcription factors during RTA transactivation of target genes.

Comprehensive polymorphism survey elucidates population structure of Saccharomyces cerevisiae

Comprehensive identification of polymorphisms among individuals within a species is essential both for studying the genetic basis of phenotypic differences and for elucidating the evolutionary history of the species. Large-scale polymorphism surveys have recently been reported for human, mouse and Arabidopsis thaliana. Here we report a nucleotide-level survey of genomic variation in a diverse collection of 63 Saccharomyces cerevisiae strains sampled from different ecological niches (beer, bread, vineyards, immunocompromised individuals, various fermentations and nature) and from locations on different continents. We hybridized genomic DNA from each strain to whole-genome tiling microarrays and detected 1.89 million single nucleotide polymorphisms, which were grouped into 101,343 distinct segregating sites. We also identified 3,985 deletion events of length >200 base pairs among the surveyed strains. We analysed the genome-wide patterns of nucleotide polymorphism and deletion variants, and measured the extent of linkage disequilibrium in S. cerevisiae. These results and the polymorphism resource we have generated lay the foundation for genome-wide association studies in yeast. We also examined the population structure of S. cerevisiae, providing support for multiple domestication events as well as insight into the origins of pathogenic strains.

Genome-wide analysis of mono-, di- and trimethylation of histone H3 lysine 4 in Arabidopsis thaliana

Analysis of the genome-wide distribution patterns of histone H3 lysine4 methylation in Arabidopsis thaliana seedlings shows that it has widespread roles in regulating gene expression.

Utilizing tiling microarrays for whole‐genome analysis in plants

The recent explosion in available genome sequence data has ushered in an era in which analysis of a whole genome can be performed in a single experiment. While DNA microarrays have long been the established technology for measuring gene expression levels, standard expression arrays use relatively few probes for each gene and are typically biased toward known and predicted gene structures. Recently, with the availability of complete genome sequences for many organisms, very-high-density oligonucleotide-based microarrays that span the entire genome have emerged as the preferred platform for genomic analysis. Whole-genome tiling microarrays can be employed for a myriad of purposes, including empirical annotation of the transcriptome, chromatin immunoprecipitation-chip studies, analysis of alternative RNA splicing, characterization of the methylation state of cytosine bases throughout a genome (methylome), and DNA polymorphism discovery. Here, we review several applications of whole-genome technology to obtain a variety of genomic-scale information in plants.

Development and evaluation of a high-throughput, low-cost genotyping platform based on oligonucleotide microarrays in rice

BackgroundWe report the development of a microarray platform for rapid and cost-effective genetic mapping, and its evaluation using rice as a model. In contrast to methods employing whole-genome tiling microarrays for genotyping, our method is based on low-cost spotted microarray production, focusing only on known polymorphic features.ResultsWe have produced a genotyping microarray for rice, comprising 880 single feature polymorphism (SFP) elements derived from insertions/deletions identified by aligning genomic sequences of the japonica cultivar Nipponbare and the indica cultivar 93-11. The SFPs were experimentally verified by hybridization with labeled genomic DNA prepared from the two cultivars. Using the genotyping microarrays, we found high levels of polymorphism across diverse rice accessions, and were able to classify all five subpopulations of rice with high bootstrap support. The microarrays were used for mapping of a gene conferring resistance to Magnaporthe grisea, the causative organism of rice blast disease, by quantitative genotyping of samples from a recombinant inbred line population pooled by phenotype.ConclusionWe anticipate this microarray-based genotyping platform, based on its low cost-per-sample, to be particularly useful in applications requiring whole-genome molecular marker coverage across large numbers of individuals.

Mapping the C. elegans noncoding transcriptome with a whole-genome tiling microarray

The number of annotated protein coding genes in the genome of Caenorhabditis elegans is similar to that of other animals, but the extent of its non-protein-coding transcriptome remains unknown. Expression profiling on whole-genome tiling microarrays applied to a mixed-stage C. elegans population verified the expression of 71% of all annotated exons. Only a small fraction (11%) of the polyadenylated transcription is non-annotated and appears to consist of approximately 3200 missed or alternative exons and 7800 small transcripts of unknown function (TUFs). Almost half (44%) of the detected transcriptional output is non-polyadenylated and probably not protein coding, and of this, 70% overlaps the boundaries of protein-coding genes in a complex manner. Specific analysis of small non-polyadenylated transcripts verified 97% of all annotated small ncRNAs and suggested that the transcriptome contains approximately 1200 small (<500 nt) unannotated noncoding loci. After combining overlapping transcripts, we estimate that at least 70% of the total C. elegans genome is transcribed.

Mapping the genome landscape using tiling array technology

With the availability of complete genome sequences for a growing number of organisms, high-throughput methods for gene annotation and analysis of genome dynamics are needed. The application of whole-genome tiling microarrays for studies of global gene expression is providing a more unbiased view of the transcriptional activity within genomes. For example, this approach has led to the identification and isolation of many novel non-protein-coding RNAs (ncRNAs), which have been suggested to comprise a major component of the transcriptome that have novel functions involved in epigenetic regulation of the genome. Additionally, tiling arrays have been recently applied to the study of histone modifications and methylation of cytosine bases (DNA methylation). Surprisingly, recent studies combining the analysis of gene expression (transcriptome) and DNA methylation (methylome) using whole-genome tiling arrays revealed that DNA methylation regulates the expression levels of many ncRNAs. Further capture and integration of additional types of genome-wide data sets will help to illuminate additional hidden features of the dynamic genomic landscape that are regulated by both genetic and epigenetic pathways in plants.

Whole-Genome Analysis of Histone H3 Lysine 27 Trimethylation in Arabidopsis

Trimethylation of histone H3 lysine 27 (H3K27me3) plays critical roles in regulating animal development, and in several cases, H3K27me3 is also required for the proper expression of developmentally important genes in plants. However, the extent to which H3K27me3 regulates plant genes on a genome-wide scale remains unknown. In addition, it is not clear whether the establishment and spreading of H3K27me3 occur through the same mechanisms in plants and animals. We identified regions containing H3K27me3 in the genome of the flowering plant Arabidopsis thaliana using a high-density whole-genome tiling microarray. The results suggest that H3K27me3 is a major silencing mechanism in plants that regulates an unexpectedly large number of genes in Arabidopsis (~4,400), and that the maintenance of H3K27me3 is largely independent of other epigenetic pathways, such as DNA methylation or RNA interference. Unlike in animals, where H3K27m3 occupies large genomic regions, in Arabidopsis, we found that H3K27m3 domains were largely restricted to the transcribed regions of single genes. Furthermore, unlike in animals systems, H3K27m3 domains were not preferentially associated with low–nucleosome density regions. The results suggest that different mechanisms may underlie the establishment and spreading of H3K27me3 in plants and animals.