Tiling Array Analysis Research Articles

Sodium Fluoride Exposure Leads to ATP Depletion and Altered RNA Decay in Escherichia coli under Anaerobic Conditions.

Although fluoride-containing compounds are widely used to inhibit bacterial growth, the reprogramming of gene expression underlying cellular responses to fluoride, especially under anaerobic conditions, is still poorly understood. Here, we compare the genome-wide transcriptomic profiles of E. coli grown in the absence (control) or presence (20 and 70 mM) of sodium fluoride (NaF) under anaerobic conditions and assess the impact of fluoride-dependent ATP depletion on RNA turnover. Tiling array analysis revealed transcripts displaying altered abundance in response to NaF treatments. Quantile-based K-means clustering uncovered a subset of genes that were highly upregulated and then downregulated in response to increased and subsequently decreased fluoride concentrations, many of which (~40%) contained repetitive extragenic palindromic (REP) sequences. Northern blot analysis of some of these highly upregulated REP-containing transcripts (i.e., osmC, proP, efeO and yghA) confirmed their considerably enhanced abundance in response to NaF treatment. An mRNA stability analysis of osmC and yghA transcripts demonstrated that fluoride treatment slows down RNA degradation, thereby enhancing RNA stability and steady-state mRNA levels. Moreover, we demonstrate that turnover of these transcripts depends on RNase E activity and RNA degradosome. Thus, we show that NaF exerts significant effects at the whole-transcriptome level under hypoxic growth (i.e., mimicking the host environment), and fluoride can impact gene expression posttranscriptionally by slowing down ATP-dependent degradation of structured RNAs. IMPORTANCE Gram-negative Escherichia coli is a rod-shaped facultative anaerobic bacterium commonly found in microaerobic/anaerobic environments, including the dental plaques of warm-blooded organisms. These latter can be treated efficiently with fluoride-rich compounds that act as anticaries agents to prevent tooth decay. Although fluoride inhibits microbial growth by affecting metabolic pathways, the molecular mechanisms underlying its activity under anaerobic conditions remain poorly defined. Here, using genome-wide transcriptomics, we explore the impact of fluoride treatments on E. coli gene expression under anaerobic conditions. We reveal key gene clusters associated with cellular responses to fluoride and define its ATP-dependent stabilizing effects on transcripts containing repetitive extragenic palindromic sequences. We demonstrate the mechanisms controlling the RNA stability of these REP-containing mRNAs. Thus, fluoride can affect gene expression posttranscriptionally by stabilizing structured RNAs.

Detecting SARS-CoV-2 and its variant strains with a full genome tiling array.

Coronavirus disease 2019 pandemic is the most damaging pandemic in recent human history. Rapid detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and variant strains is paramount for recovery from this pandemic. Conventional SARS-CoV-2 tests interrogate only limited regions of the whole SARS-CoV-2 genome, which are subjected to low specificity and miss the opportunity of detecting variant strains. In this work, we developed the first SARS-CoV-2 tiling array that captures the entire SARS-CoV-2 genome at single nucleotide resolution and offers the opportunity to detect point mutations. A thorough bioinformatics protocol of two base calling methods has been developed to accompany this array. To demonstrate the effectiveness of the tiling array, we genotyped all genomic positions of eight SARS-CoV-2 samples. Using high-throughput sequencing as the benchmark, we show that the tiling array had a genome-wide accuracy of at least 99.5%. From the tiling array analysis results, we identified the D614G mutation in the spike protein in four of the eight samples, suggesting the widespread distribution of this variant at the early stage of the outbreak in the United States. Two additional nonsynonymous mutations were identified in one sample in the nucleocapsid protein (P13L and S197L), which may complicate future vaccine development. With around $5 per array, supreme accuracy, and an ultrafast bioinformatics protocol, the SARS-CoV-2 tiling array makes an invaluable toolkit for combating current and future pandemics. Our SARS-CoV-2 tilting array is currently utilized by Molecular Vision, a CLIA-certified lab for SARS-CoV-2 diagnosis.

Genome-Wide Tiling Array Analysis of HPV-Induced Warts Reveals Aberrant Methylation of Protein-Coding and Non-Coding Regions

The human papillomaviruses (HPV) are a group of double-stranded DNA viruses that exhibit an exclusive tropism for squamous epithelia. HPV can either be low- or high-risk depending on its ability to cause benign lesions or cancer, respectively. Unsurprisingly, the majority of epigenetic research has focused on the high-risk HPV types, neglecting the low-risk types in the process. Therefore, the main objective of this study is to better understand the epigenetics of wart formation by investigating the differences in methylation between HPV-induced cutaneous warts and normal skin. A number of clear and very significant differences in methylation patterns were found between cutaneous warts and normal skin. Around 55% of the top-ranking 100 differentially methylated genes in warts were protein coding, including the EXOC4, KCNU, RTN1, LGI1, IRF2, and NRG1 genes. Additionally, non-coding RNA genes, such as the AZIN1-AS1, LINC02008, and MGC27382 genes, constituted 11% of the top-ranking 100 differentially methylated genes. Warts exhibited a unique pattern of methylation that is a possible explanation for their transient nature. Since the genetics of cutaneous wart formation are not completely known, the findings of the present study could contribute to a better understanding of how HPV infection modulates host methylation to give rise to warts in the skin.

Transcriptional Regulation of Ovarian Steroidogenic Genes: Recent Findings Obtained from Stem Cell-Derived Steroidogenic Cells.

Ovaries represent one of the primary steroidogenic organs, producing estrogen and progesterone under the regulation of gonadotropins during the estrous cycle. Gonadotropins fluctuate the expression of various steroidogenesis-related genes, such as those encoding steroidogenic enzymes, cholesterol deliverer, and electronic transporter. Steroidogenic factor-1 (SF-1)/adrenal 4-binding protein (Ad4BP)/NR5A1 and liver receptor homolog-1 (LRH-1) play important roles in these phenomena via transcriptional regulation. With the aid of cAMP, SF-1/Ad4BP and LRH-1 can induce the differentiation of stem cells into steroidogenic cells. This model is a useful tool for studying the molecular mechanisms of steroidogenesis. In this article, we will provide insight into the transcriptional regulation of steroidogenesis-related genes in ovaries that are revealed from stem cell-derived steroidogenic cells. Using the cells derived from the model, novel SF-1/Ad4BP- and LRH-1-regulated genes were identified by combined DNA microarray and promoter tiling array analyses. The interaction of SF-1/Ad4BP and LRH-1 with transcriptional regulators in the regulation of ovarian steroidogenesis was also revealed.

Genome-wide screening of potential RNase Y-processed mRNAs in the M49 serotype Streptococcus pyogenes NZ131.

RNase Y is a major endoribonuclease in Group A streptococcus (GAS) and other Gram‐positive bacteria. Our previous study showed that RNase Y was involved in mRNA degradation and processing in GAS. We hypothesized that mRNA processing regulated the expression of important GAS virulence factors via altering their mRNA stabilities and that RNase Y mediated at least some of the mRNA‐processing events. The aims of this study were to (1) identify mRNAs that were processed by RNase Y and (2) confirm the mRNA‐processing events. The transcriptomes of Streptococcus pyogenes NZ131 wild type and its RNase Y mutant (Δrny) were examined with RNA‐seq. The data were further analyzed to define GAS operons. The mRNA stabilities of the wild type and Δrny at subgene level were determined with tiling array analysis. Operons displaying segmental stability in the wild type but not in the Δrny were predicted to be RNase Y processed. Overall 865 operons were defined and their boundaries predicted. Further analysis narrowed down 15 mRNAs potentially processed by RNase Y. A selection of four candidates including folC1 (folylpolyglutamate synthetase), prtF (fibronectin‐binding protein), speG (streptococcal exotoxin G), ropB (transcriptional regulator of speB), and ypaA (riboflavin transporter) mRNAs was examined with Northern blot analysis. However, only folC1 was confirmed to be processed, but it is unlikely that RNase Y is responsible. We conclude that GAS use RNase Y to selectively process mRNA, but the overall impact is confined to selected virulence factors.

Histone acetyltransferase 7 (KAT7)-dependent intragenic histone acetylation regulates endothelial cell gene regulation

Although the functional role of chromatin marks at promoters in mediating cell-restricted gene expression has been well characterized, the role of intragenic chromatin marks is not well understood, especially in endothelial cell (EC) gene expression. Here, we characterized the histone H3 and H4 acetylation profiles of 19 genes with EC-enriched expression via locus-wide chromatin immunoprecipitation followed by ultra-high-resolution (5 bp) tiling array analysis in ECs versus non-ECs throughout their genomic loci. Importantly, these genes exhibit differential EC enrichment of H3 and H4 acetylation in their promoter in ECs versus non-ECs. Interestingly, VEGFR-2 and VEGFR-1 show EC-enriched acetylation across broad intragenic regions and are up-regulated in non-ECs by histone deacetylase inhibition. It is unclear which histone acetyltransferases (KATs) are key to EC physiology. Depletion of KAT7 reduced VEGFR-2 expression and disrupted angiogenic potential. Microarray analysis of KAT7-depleted ECs identified 263 differentially regulated genes, many of which are key for growth and angiogenic potential. KAT7 inhibition in zebrafish embryos disrupted vessel formation and caused loss of circulatory integrity, especially hemorrhage, all of which were rescued with human KAT7. Notably, perturbed EC-enriched gene expression, especially the VEGFR-2 homologs, contributed to these vascular defects. Mechanistically, KAT7 participates in VEGFR-2 transcription by mediating RNA polymerase II binding, H3 lysine 14, and H4 acetylation in its intragenic region. Collectively, our findings support the importance of differential histone acetylation at both promoter and intragenic regions of EC genes and reveal a previously underappreciated role of KAT7 and intragenic histone acetylation in regulating VEGFR-2 and endothelial function.

Global transcriptional landscape and promoter mapping of the gut commensal Bifidobacterium breve UCC2003

BackgroundBifidobacterium breve represents a common member of the infant gut microbiota and its presence in the gut has been associated with host well being. For this reason it is relevant to investigate and understand the molecular mechanisms underlying the establishment, persistence and activities of this gut commensal in the host environment.ResultsThe assessment of vegetative promoters in the bifidobacterial prototype Bifidobacterium breve UCC2003 was performed employing a combination of RNA tiling array analysis and cDNA sequencing. Canonical −10 (TATAAT) and −35 (TTGACA) sequences were identified upstream of transcribed genes or operons, where deviations from this consensus correspond to transcription level variations. A Random Forest analysis assigned the −10 region of B. breve promoters as the element most impacting on the level of transcription, followed by the spacer length and the 5’-UTR length of transcripts. Furthermore, our transcriptome study also identified rho-independent termination as the most common and effective termination signal of highly and moderately transcribed operons in B. breve.ConclusionThe present study allowed us to identify genes and operons that are actively transcribed in this organism during logarithmic growth, and link promoter elements with levels of transcription of essential genes in this organism. As homologs of many of our identified genes are present across the whole genus Bifidobacterium, our dataset constitutes a transcriptomic reference to be used for future investigations of gene expression in members of this genus.

Conserved two-component Hik34-Rre1 module directly activates heat-stress inducible transcription of major chaperone and other genes in Synechococcus elongatus PCC 7942.

Bacteria and other organisms, including cyanobacteria, employ two-component signal transducing modules comprising histidine kinases and response regulators to acclimate to changing environments. While the number and composition of these modules differ among cyanobacteria, two response regulators that contain DNA binding domains, RpaB and Rre1, are conserved in all sequenced cyanobacterial genomes and are essential for viability. Although RpaB negatively or positively regulates high light and other stress-responsive gene expression, little is known about the function of Rre1. Here, they investigated the direct regulatory targets of Rre1 in the cyanobacterium Synechococcus elongatus PCC 7942. Chromatin immunoprecipitation and high-density tiling array analysis were used to map Rre1 binding sites. The sites included promoter regions for chaperone genes such as dnaK2, groESL-1, groEL-2, hspA and htpG, as well as the group 2 sigma factor gene rpoD2. In vivo and in vitro analyses revealed that Rre1 phosphorylation level, DNA binding activity and adjacent gene transcription increased in response to heat stress. These responses were much diminished in a knock-out mutant of Hik34, a previously identified heat shock regulator. Based on our results, we propose Hik34-Rre1 is the heat shock-responsive signaling module that positively regulates major chaperone and other genes in cyanobacteria.

Novel Epigenetic Markers on Chromosome 21 for Noninvasive Prenatal Testing of Fetal Trisomy 21

Until now, fetal placenta-specific epigenetic markers for noninvasive prenatal testing of fetal trisomy 21 (T21) have been identified based only on differences in tissue-specific epigenetic characteristics between placenta and maternal blood, but these characteristics have not been validated in T21 placenta. We aimed to discover novel epigenetic markers on chromosome 21 that show a hypermethylated pattern in fetal placenta compared with blood, regardless of the presence of T21. We performed a high-resolution tiling array analysis of chromosome 21 using the methylated-CpG binding domain protein-based method. We identified 93 epigenetic regions that showed fetal placenta-specific differential methylation patterns; among these, three regions showed fetal placenta-specific methylation patterns in T21 placenta samples. The methylation patterns of these three regions in the array were confirmed by bisulfite direct sequencing. The three regions were detectable in first-trimester maternal plasma. Moreover, a combination of their methylation ratio achieved high diagnostic accuracy for noninvasive prenatal testing of fetal T21 by further statistical analysis. These three novel regions with fetal placenta-specific differential methylation patterns on chromosome 21 were identified irrespective of the presence of T21. Our findings suggest that epigenetic characteristics of markers according to the presence or absence of T21 should be considered in the development of noninvasive prenatal testing of fetal T21 using fetal placenta-specific epigenetic markers.

The first decade of estrogen receptor cistromics in breast cancer.

The advent of genome-wide transcription factor profiling has revolutionized the field of breast cancer research. Estrogen receptor α (ERα), the major drug target in hormone receptor-positive breast cancer, has been known as a key transcriptional regulator in tumor progression for over 30 years. Even though this function of ERα is heavily exploited and widely accepted as an Achilles heel for hormonal breast cancer, only since the last decade we have been able to understand how this transcription factor is functioning on a genome-wide scale. Initial ChIP-on-chip (chromatin immunoprecipitation coupled with tiling array) analyses have taught us that ERα is an enhancer-associated factor binding to many thousands of sites throughout the human genome and revealed the identity of a number of directly interacting transcription factors that are essential for ERα action. More recently, with the development of massive parallel sequencing technologies and refinements thereof in sample processing, a genome-wide interrogation of ERα has become feasible and affordable with unprecedented data quality and richness. These studies have revealed numerous additional biological insights into ERα behavior in cell lines and especially in clinical specimens. Therefore, what have we actually learned during this first decade of cistromics in breast cancer and where may future developments in the field take us?

Prenatal diesel exhaust exposure disrupts the DNA methylation profile in the brain of mouse offspring.

Prenatal diesel exhaust (DE) exposure is associated with detrimental health effects in offspring. Although previous reports suggest that DE exposure affects the brain of offspring in the developmental period, the molecular events associated with the health effects have largely remained unclear. We hypothesized that the DNA methylation state would be disrupted by prenatal DE exposure. In the present study, the authors examined the genome-wide DNA methylation state of the gene promoter and bioinformatically analyzed the obtained data to identify the molecular events related to disrupted DNA methylation. Pregnant C57BL/6J mice were exposed to DE (DEP; 0.1 mg/m3) in an inhalation chamber on gestational days 0-16. Brains were collected from 1-day-old and 21-day-old offspring. The genomewide DNA methylation state of the brain genome was analyzed by methylation-specific DNA immunoprecipitation and subsequent promoter tiling array analysis. The genes in which the DNA methylation level was affected by prenatal DE exposure were bioinformatically categorized using Gene Ontology (GO). Differentially methylated DNA regions were detected in all chromosomes in brains collected from both 1-day-old and 21-day-old offspring. Altered DNA methylation was observed independently of the presence of CpG island. Bioinformatic interpretation using GO terms showed that differentially methylated genes with CpG islands in their promoter were commonly enriched in neuronal differentiation and neurogenesis. The results suggest that prenatal DE exposure causes genome-wide disruption of DNA methylation in the brain. Disrupted DNA methylation would disturb neuronal development in the developmental period and may be associated with health and disease in later life.

Oncogenic Ras influences the expression of multiple lncRNAs.

Recent ultrahigh-density tiling array and large-scale transcriptome analysis have revealed that large numbers of long non-coding RNAs (lncRNAs) are transcribed in mammals. Several lncRNAs have been implicated in transcriptional regulation, organization of nuclear structure, and post-transcriptional processing. However, the regulation of expression of lncRNAs is less well understood. Here, we show that the exogenous and endogenous expression of an oncogenic form of small GTPase Ras (called oncogenic Ras) decrease the expression of lncRNA ANRIL (antisense non-coding RNA in the INK4 locus), which is involved in the regulation of cellular senescence. We also show that forced expression of oncogenic Ras increases the expression of lncRNA PANDA (p21 associated ncRNA DNA damage activated), which is involved in the regulation of apoptosis. Microarray analysis demonstrated that expression of multiple lncRNAs fluctuated by forced expression of oncogenic Ras. These findings indicate that oncogenic Ras regulates the expression of a large number of lncRNAs including functional lncRNAs, such as ANRIL and PANDA.

Potential epigenetic dysregulation of genes associated with MODY and type 2 diabetes in humans exposed to a diabetic intrauterine environment: An analysis of genome-wide DNA methylation

ObjectiveThe aim of this study is to investigate the potential role of DNA methylation in mediating the increased risk of developing type 2 diabetes in offspring of mothers who had diabetes during pregnancy. Materials and MethodsPeripheral blood leukocytes were collected from non-diabetic Pima Indians who were either offspring of diabetic mothers (ODM; n=14) or offspring of nondiabetic mothers (ONDM; n=14). The two groups were matched for age, sex, age of mother, and fraction of Pima ethnicity. Differentially methylated regions were determined using a MeDIP-chip assay on an Affymetrix Human Tiling 2.0R Array. Data were analyzed using the model based analysis of tiling arrays (MAT) algorithm, and 4883 regions overlapping with putative promoters, were identified as differentially methylated, having met an empirically derived threshold (nominal p<0.0077). The list of genes with differentially methylated promoters was subjected to KEGG pathway analysis to determine canonical metabolic pathways enriched by these genes. ResultsPathway analysis of genes with differentially methylated promoters identified the top 3 enriched pathways as maturity onset diabetes of the young (MODY), type 2 diabetes, and Notch signaling. Several genes in these pathways are known to affect pancreatic development and insulin secretion. ConclusionsThese findings support the hypothesis that epigenetic changes may increase the risk of type 2 diabetes via an effect on β-cell function in the offspring of mothers with diabetes during pregnancy.

Disease specific characteristics of fetal epigenetic markers for non-invasive prenatal testing of trisomy 21

BackgroundNon-invasive prenatal testing of trisomy 21 (T21) is being actively investigated using fetal-specific epigenetic markers (EPs) that are present in maternal plasma. Recently, 12 EPs on chromosome 21 were identified based on tissue-specific epigenetic characteristics between placenta and blood, and demonstrated excellent clinical performance in the non-invasive detection of fetal T21. However, the disease-specific epigenetic characteristics of the EPs have not been established. Therefore, we validated the disease-specific epigenetic characteristics of these EPs for use in non-invasive detection of fetal T21.MethodsWe performed a high-resolution tiling array analysis of human chromosome 21 using a methyl-CpG binding domain-based protein (MBD) method with whole blood samples from non-pregnant normal women, whole blood samples from pregnant normal women, placenta samples of normal fetuses, and placenta samples of T21 fetuses. Tiling array results were validated by bisulfite direct sequencing and qPCR.ResultsAmong 12 EPs, only four EPs were confirmed to be hypermethylated in normal placenta and hypomethylated in blood. One of these four showed a severe discrepancy in the methylation patterns of T21 placenta samples, and another was located within a region of copy number variations. Thus, two EPs were confirmed to be potential fetal-specific markers based on their disease-specific epigenetic characteristics. The array results of these EPs were consisted with the results obtained by bisulfite direct sequencing and qPCR. Moreover, the two EPs were detected in maternal plasma.ConclusionsWe validated that two EPs have the potential to be fetal-specific EPs which is consistent with their disease-specific epigenetic characteristics. The findings of this study suggest that disease-specific epigenetic characteristics should be considered in the development of fetal-specific EPs for non-invasive prenatal testing of T21.

Transcriptional regulation of human ferredoxin reductase through an intronic enhancer in steroidogenic cells

Ferredoxin reductase (FDXR, also known as adrenodoxin reductase) is a mitochondrial flavoprotein that transfers electrons from NADPH to mitochondrial cytochrome P450 enzymes, mediating the function of an iron–sulfur cluster protein, ferredoxin. FDXR functions in various metabolic processes including steroidogenesis. It is well known that multiple steroidogenic enzymes are regulated by a transcription factor steroidogenic factor-1 (SF-1, also known as Ad4BP). Previously, we have shown that SF-1 transduction causes human mesenchymal stem cell differentiation into steroidogenic cells. Genome-wide analysis of differentiated cells, using a combination of DNA microarray and promoter tiling array analyses, showed that FDXR is a novel SF-1 target gene. In this study, the transcriptional regulatory mechanism of FDXR was examined in steroidogenic cells. A chromatin immunoprecipitation assay revealed that a novel SF-1 binding region was located within intron 2 of the human FDXR gene. Luciferase reporter assays showed that FDXR transcription was activated through the novel SF-1 binding site within intron 2. Endogenous SF-1 knockdown in human adrenocortical H295R and KGN cells decreased FDXR expression. In H295R cells, strong binding of two histone markers of active enhancers, histones H3K27ac and H3K4me2, were detected near the SF-1 binding site within intron 2. Furthermore, the binding of these histone markers was decreased concurrent with SF-1 knockdown in H295R cells. These results indicated that abundant FDXR expression in these steroidogenic cells was maintained through SF-1 binding to the intronic enhancer of the FDXR gene.

The Role of Ctk1 Kinase in Termination of Small Non-Coding RNAs

Transcription termination in Saccharomyces cerevisiae can be performed by at least two distinct pathways and is influenced by the phosphorylation status of the carboxy-terminal domain (CTD) of RNA polymerase II (Pol II). Late termination of mRNAs is performed by the CPF/CF complex, the recruitment of which is dependent on CTD-Ser2 phosphorylation (Ser2P). Early termination of shorter cryptic unstable transcripts (CUTs) and small nucleolar/nuclear RNAs (sno/snRNAs) is performed by the Nrd1-Nab3-Sen1 (NNS) complex that binds phosphorylated CTD-Ser5 (Ser5P) via the CTD-interacting domain (CID) of Nrd1p. In this study, mutants of the different termination pathways were compared by genome-wide expression analysis. Surprisingly, the expression changes observed upon loss of the CTD-Ser2 kinase Ctk1p are more similar to those derived from alterations in the Ser5P-dependent NNS pathway, than from loss of CTD-Ser2P binding factors. Tiling array analysis of ctk1Δ cells reveals readthrough at snoRNAs, at many cryptic unstable transcripts (CUTs) and stable uncharacterized transcripts (SUTs), but only at some mRNAs. Despite the suggested predominant role in termination of mRNAs, we observed that a CTK1 deletion or a Pol II CTD mutant lacking all Ser2 positions does not result in a global mRNA termination defect. Rather, termination defects in these strains are widely observed at NNS-dependent genes. These results indicate that Ctk1p and Ser2 CTD phosphorylation have a wide impact in termination of small non-coding RNAs but only affect a subset of mRNA coding genes.

Mining the genome of Arabidopsis thaliana as a basis for the identification of novel bioactive peptides involved in oxidative stress tolerance

Although evidence has accumulated on the role of plant peptides in the response to external conditions, the number of peptide-encoding genes in the genome is still underestimated. Using tiling arrays, we identified 176 unannotated transcriptionally active regions (TARs) in Arabidopsis thaliana that were induced upon oxidative stress generated by the herbicide paraquat (PQ). These 176 TARs could be translated into 575 putative oxidative stress-induced peptides (OSIPs). A high-throughput functional assay was used in the eukaryotic model organism Saccharomyces cerevisiae allowing us to test for bioactive peptides that increase oxidative stress tolerance. In this way, we identified three OSIPs that, upon overexpression in yeast, resulted in a significant rise in tolerance to hydrogen peroxide (H2O2). For one of these peptides, the decapeptide OSIP108, exogenous application to H2O2-treated yeast also resulted in significantly increased survival. OSIP108 is contained within a pseudogene and is induced in A. thaliana leaves by both the reactive oxygen species-inducer PQ and the necrotrophic fungal pathogen Botrytis cinerea. Moreover, infiltration and overexpression of OSIP108 in A. thaliana leaves resulted in increased tolerance to treatment with PQ. In conclusion, the identification and characterization of OSIP108 confirms the validity of our high-throughput approach, based on tiling array analysis in A. thaliana and functional screening in yeast, to identify bioactive peptides.

Lytic Infection of Lactococcus lactis by Bacteriophages Tuc2009 and c2 Triggers Alternative Transcriptional Host Responses

Here we present an entire temporal transcriptional profile of Lactococcus lactis subsp. cremoris UC509.9 undergoing lytic infection with two distinct bacteriophages, Tuc2009 and c2. Furthermore, corresponding high-resolution whole-phage genome tiling arrays of both bacteriophages were performed throughout lytic infection. Whole-genome microarrays performed at various time points postinfection demonstrated a rather modest impact on host transcription. The majority of changes in the host transcriptome occur during late infection stages; few changes in host gene transcription occur during the immediate and early infection stages. Alterations in the L. lactis UC509.9 transcriptome during lytic infection appear to be phage specific, with relatively few differentially transcribed genes shared between cells infected with Tuc2009 and those infected with c2. Despite the apparent lack of a coordinated general phage response, three themes common to both infections were noted: alternative transcription of genes involved in catabolic flux and energy production, differential transcription of genes involved in cell wall modification, and differential transcription of genes involved in the conversion of ribonucleotides to deoxyribonucleotides. The transcriptional profiles of both bacteriophages during lytic infection generally correlated with the findings of previous studies and allowed the confirmation of previously predicted promoter sequences. In addition, the host transcriptional response to lysogenization with Tuc2009 was monitored along with tiling array analysis of Tuc2009 in the lysogenic state. Analysis identified 44 host genes with altered transcription during lysogeny, 36 of which displayed levels of transcription significantly reduced from those for uninfected cells.

Microarray-Based Capture of Novel Expressed Cell Type–Specific Transfrags (CoNECT) to Annotate Tissue-Specific Transcription inDrosophila melanogaster

Faithful annotation of tissue-specific transcript isoforms is important not only to understand how genes are organized and regulated but also to identify potential novel, unannotated exons of genes, which may be additional targets of mutation in disease states or while performing mutagenic screens. We have developed a microarray enrichment methodology followed by long-read, next-generation sequencing for identification of unannotated transcript isoforms expressed in two Drosophila tissues, the ovary and the testis. Even with limited sequencing, these studies have identified a large number of novel transcription units, including 5′ exons and extensions, 3′ exons and extensions, internal exons and exon extensions, gene fusions, and both germline-specific splicing events and promoters. Additionally, comparing our capture dataset with tiling array and traditional RNA-seq analysis, we demonstrate that our enrichment strategy is able to capture low-abundance transcripts that cannot readily be identified by the other strategies. Finally, we show that our methodology can help identify transcriptional signatures of minority cell types within the ovary that would otherwise be difficult to reveal without the CoNECT enrichment strategy. These studies introduce an efficient methodology for cataloging tissue-specific transcriptomes in which specific classes of genes or transcripts can be targeted for capture and sequence, thus reducing the significant sequencing depth normally required for accurate annotation.

Bimodal control of Hoxd gene transcription in the spinal cord defines two regulatory subclusters

The importance of Hox genes in the specification of neuronal fates in the spinal cord has long been recognized. However, the transcriptional controls underlying their collinear expression domains remain largely unknown. Here we show in mice that the correspondence between the physical order of Hoxd genes and their rostral expression boundaries, although respecting spatial collinearity, does not display a fully progressive distribution. Instead, two major anteroposterior boundaries are detected, coinciding with the functional subdivision of the spinal cord. Tiling array analyses reveal two distinct blocks of transcription, regulated independently from one another, that define the observed expression boundaries. Targeted deletions in vivo that remove the genomic fragments separating the two blocks induce ectopic expression of posterior genes. We further evaluate the independent regulatory potential and transcription profile of each gene locus by a tiling array approach using a contiguous series of transgenes combined with locus-specific deletions. Our work uncovers a bimodal type of HoxD spatial collinearity in the developing spinal cord that relies on two separate 'enhancer mini-hubs' to ensure correct Hoxd gene expression levels while maintaining their appropriate anteroposterior boundaries.