The Genomics Education Partnership: First findings on genomics research in community colleges.

The Genomics Education Partnership (GEP), a consortium of diverse colleges/universities, provides support for integrating genomics research into undergraduate curricula. To increase research opportunities for underrepresented students, GEP is expanding to more community colleges (CC). Genomics research, requiring only a computer with internet access, may be particularly accessible for 2-year institutions with limited research capacity and significant budget constraints. To understand how GEP supports student research at CCs, we analyzed student knowledge and self-reported outcomes. We found that CC student gains are comparable to non-CC student gains, with improvements in attitudes toward science and thriving in science. Our early findings suggest that the GEP model of centralized support with flexible CURE implementation benefits CC students and may help mitigate barriers to implementing research at CCs.

High-Quality Draft Genome Sequence of Actinobacterium Kibdelosporangium sp. MJ126-NF4, Producer of Type II Polyketide Azicemicins, Using Illumina and PacBio Technologies.

Here, we report the high-quality draft genome sequence of actinobacterium Kibdelosporangium sp. MJ126-NF4, producer of the type II polyketide azicemicins, obtained using Illumina and PacBio sequencing technologies. The 11.75-Mbp genome contains >11,000 genes and 22 polyketide and nonribosomal peptide natural product gene clusters.

Genomic analysis of Saccharomyces cerevisiae isolates that grow optimally with glucose as the sole carbon source.

A population of Saccharomyces cerevisiae was cultured for approximately 450 generations in the presence of high glucose to select for genetic variants that grew optimally under these conditions. Using the parental strain BY4741 as the starting population, an evolved culture was obtained after aerobic growth in a high glucose medium for approximately 450 generations. After the evolution period, three single colony isolates were selected for analysis. Next-generation Ion Torrent sequencing was used to evaluate genetic changes. Greater than 100 deletion/insertion changes were found with approximately half of these effecting genes. Additionally, over 180 SNPs were identified with more than one-quarter of these resulting in a nonsynonymous mutation. Affymetrix DNA microarrays and RNseq analysis were used to determine differences in gene expression in the evolved strains compared to the parental strain. It was established that approximately 900 genes demonstrated significantly altered expression in the evolved strains relative to the parental strain. Many of these genes showed similar alterations in their expression in all three evolved strains. Interestingly, genes with altered expression in the three evolved strains included genes with a role in oxidative metabolism. Overall these results are consistent with the physiological observations of optimal growth with glucose as the carbon source. Namely, the decreased ethanol production suggest that the underlying metabolism switched from fermentation to respiration during the selection for optimal growth on glucose.

Dual primer emulsion PCR for next- generation DNA sequencing.

We have developed a highly sensitive single-molecule clonal amplification method called dual primer emulsion PCR (DPePCR) for next-generation DNA sequencing. The approach is similar in concept to standard emulsion PCR; however, in DPePCR both primers are attached to the beads, therefore following PCR amplification, both strands of the PCR products are attached to the beads. The ends of each strand can be freed for analysis by restriction digestion of the bridged PCR fragments, which allows efficient paired-end sequencing of fragment libraries.

Differential transcriptomic responses of Biomphalaria glabrata (Gastropoda, Mollusca) to bacteria and metazoan parasites, Schistosoma mansoni and Echinostoma paraensei (Digenea, Platyhelminthes).

A 70-mer-oligonucleotide-based microarray (1152 features) that emphasizes stress and immune responses factors was constructed to study transcriptomic responses of the snail Biomphalaria glabrata to different immune challenges. In addition to sequences with relevant putative ID and Gene Ontology (GO) annotation, the array features non-immune factors and unknown B. glabrata ESTs for functional gene discovery. The transcription profiles of B. glabrata (3 biological replicates, each a pool of 5 snails) were recorded at 12h post-wounding, exposure to Gram negative or Gram positive bacteria (Escherichia coli and Micrococcus luteus, respectively), or infection with compatible trematode parasites (Schistosoma mansoni or Echinostoma paraensei, 20 miracidia/snail), relative to controls, using universal reference RNA. The data were subjected to Significance Analysis for Microarrays (SAM), with a false positive rate (FPR)

Towards an understanding of the mechanism of action of praziquantel.

Although praziquantel (PZQ) has been used to treat schistosomiasis for over 20 years its mechanism of action remains unknown. We have developed an assay based on the transcriptional response of Schistosoma mansoni PR-1 to heat shock to confirm that while 6-week post-infection (p.i.) schistosomes are sensitive to PZQ, 4-week p.i. schistosomes are not. Further, we have used this assay to demonstrate that in mice this sensitivity develops between days 37 and 40 p.i. When PZQ is linked to the fluorophore BODIPY to aid microscopic visualization, it appears to enter the cells of intact 4 and 6-week p.i. schistosomes as well as mammalian NIH 3T3 cells with ease suggesting that the differential effects of PZQ is not based on cell exclusion. A transcriptomal analysis of gene expression between 4 and 6 weeks p.i. revealed 607 up-regulated candidate genes whose products are potential PZQ targets. A comparison of this gene list with that of genes expressed by PZQ sensitive miracidia reduced this target list to 247 genes, including a number involved in aerobic metabolism and cytosolic calcium regulation. Finally, we also report the effect of an in vitro sub-lethal exposure of PZQ on the transcriptome of S. mansoni PR-1. Annotation of genes differentially regulated by PZQ exposure suggests that schistosomes may undergo a transcriptomic response similar to that observed during oxidative stress.

Microarray based analysis of temperature and oxidative stress induced messenger RNA in Schistosoma mansoni.

The body's defense against schistosome infection can take many forms. For example, upon developing acute schistosomiasis, patients often have fever coinciding with larval maturation, migration and early oviposition. As the infection becomes established, the parasite comes under oxidative stress generated by the host immune system. The most common treatment for schistosomiasis is the anti-helminthic drug praziquantel. Its effectiveness, however, is limited due to its inability to kill schistosomes 2-4 weeks post-infection. Clearly there is a need for new anti-schistosomal drugs. We hypothesize that gene products expressed as part of a protective response against heat and/or oxidative stress are potential therapeutic targets for future drug development. Using a 12,166 element oligonucleotide microarray to characterize Schistosoma mansoni genes induced by heat and oxidative stress we found that 1878 S. mansoni elements were significantly induced by heat stress. These included previously reported heat-shock genes expressing homologs of HSP40, HSP70 and HSP86. One thousand and one elements were induced by oxidative stress including those expressing homologs of superoxide dismutase, glutathione peroxidase and aldehyde dehydrogenase. Seventy-two elements were common to both stressors and could potentially be exploited in the development of novel anti-schistosomal therapeutics.

Characterization of differentiated quiescent and nonquiescent cells in yeast stationary-phase cultures.

Cells in glucose-limited Saccharomyces cerevisiae cultures differentiate into quiescent (Q) and nonquiescent (NQ) fractions before entering stationary phase. To understand this differentiation, Q and NQ cells from 101 deletion-mutant strains were tested for viability and reproductive capacity. Eleven mutants that affected one or both phenotypes in Q or NQ fractions were identified. NQ fractions exhibit a high level of petite colonies, and nine mutants affecting this phenotype were identified. Microarray analysis revealed >1300 mRNAs distinguished Q from NQ fractions. Q cell-specific mRNAs encode proteins involved in membrane maintenance, oxidative stress response, and signal transduction. NQ-cell mRNAs, consistent with apoptosis in these cells, encode proteins involved in Ty-element transposition and DNA recombination. More than 2000 protease-released mRNAs were identified only in Q cells, consistent with these cells being physiologically poised to respond to environmental changes. Our results indicate that Q and NQ cells differentiate significantly, with Q cells providing genomic stability and NQ cells providing nutrients to Q cells and a regular source of genetic diversity through mutation and transposition. These studies are relevant to chronological aging, cell cycle, and genome evolution, and they provide insight into complex responses that even simple organisms have to starvation.

Isolation of quiescent and nonquiescent cells from yeast stationary-phase cultures.

Quiescence is the most common and, arguably, most poorly understood cell cycle state. This is in part because pure populations of quiescent cells are typically difficult to isolate. We report the isolation and characterization of quiescent and nonquiescent cells from stationary-phase (SP) yeast cultures by density-gradient centrifugation. Quiescent cells are dense, unbudded daughter cells formed after glucose exhaustion. They synchronously reenter the mitotic cell cycle, suggesting that they are in a G(0) state. Nonquiescent cells are less dense, heterogeneous, and composed of replicatively older, asynchronous cells that rapidly lose the ability to reproduce. Microscopic and flow cytometric analysis revealed that nonquiescent cells accumulate more reactive oxygen species than quiescent cells, and over 21 d, about half exhibit signs of apoptosis and necrosis. The ability to isolate both quiescent and nonquiescent yeast cells from SP cultures provides a novel, tractable experimental system for studies of quiescence, chronological and replicative aging, apoptosis, and the cell cycle.

Release of extraction-resistant mRNA in stationary phase Saccharomyces cerevisiae produces a massive increase in transcript abundance in response to stress.

BACKGROUND: As carbon sources are exhausted, Saccharomyces cerevisiae cells exhibit reduced metabolic activity and cultures enter the stationary phase. We asked whether cells in stationary phase cultures respond to additional stress at the level of transcript abundance. RESULTS: Microarrays were used to quantify changes in transcript abundance in cells from stationary phase cultures in response to stress. More than 800 mRNAs increased in abundance by one minute after oxidative stress. A significant number of these mRNAs encode proteins involved in stress responses. We tested whether mRNA increases were due to new transcription, rapid poly-adenylation of message (which would not be detected by microarrays), or potential release of mature mRNA present in the cell but resistant to extraction during RNA isolation. Examination of the response to oxidative stress in an RNA polymerase II mutant, rpb1-1, suggested that new transcription was not required. Quantitative RT-PCR analysis of a subset of these transcripts further suggested that the transcripts present in isolated total RNA from stationary phase cultures were polyadenylated. In contrast, over 2,000 transcripts increased after protease treatment of cell-free lysates from stationary phase but not exponentially growing cultures. Different subsets of transcripts were released by oxidative stress and temperature upshift, suggesting that mRNA release is stress-specific. CONCLUSIONS: Cells in stationary phase cultures contain a large number of extraction-resistant mRNAs in a protease-labile, rapidly releasable form. The transcript release appears to be stress-specific. We hypothesize that these transcripts are associated with P-bodies.

An automated, pressure-driven sampling device for harvesting from liquid cultures for genomic and biochemical analyses.

Here we describe an automated, pressure-driven, sampling device for harvesting 10 to 30 ml samples, in replicate, with intervals as short as 10 s. Correlation between biological replicate time courses measured by microarrays was extremely high. The sampler enables sampling at intervals within the range of many important biological processes.

A hidden-state Markov model for cell population deconvolution.

Microarrays measure gene expression typically from a mixture of cell populations during different stages of a biological process. However, the specific effects of the distinct or pure populations on measured gene expression are difficult or impossible to determine. The ability to deconvolve measured gene expression into the contributions from pure populations is critical to maximizing the potential of microarray analysis for investigating complex biological processes. In this paper, we describe a novel approach called the multinomial hidden Markov model (MHMM) that produces: (i) a maximum a posteriori estimate of the fraction represented by each pure population and (ii) gene expression values for each pure population. Our method uses an unsupervised, probabilistic approach for handling missing data points and clusters genes based on expression in pure populations. MHMM, used with several yeast datasets, identified statistically significant temporal dynamics. This method, unlike the linear decomposition models used previously for deconvolution, can extract information from different types of data, does not require a priori identification of pure gene expression, exploits the temporal nature of time series data, and is less affected by missing data.

Hyperspectral microarray scanning: impact on the accuracy and reliability of gene expression data.

BACKGROUND: Commercial microarray scanners and software cannot distinguish between spectrally overlapping emission sources, and hence cannot accurately identify or correct for emissions not originating from the labeled cDNA. We employed our hyperspectral microarray scanner coupled with multivariate data analysis algorithms that independently identify and quantitate emissions from all sources to investigate three artifacts that reduce the accuracy and reliability of microarray data: skew toward the green channel, dye separation, and variable background emissions. RESULTS: Here we demonstrate that several common microarray artifacts resulted from the presence of emission sources other than the labeled cDNA that can dramatically alter the accuracy and reliability of the array data. The microarrays utilized in this study were representative of a wide cross-section of the microarrays currently employed in genomic research. These findings reinforce the need for careful attention to detail to recognize and subsequently eliminate or quantify the presence of extraneous emissions in microarray images. CONCLUSION: Hyperspectral scanning together with multivariate analysis offers a unique and detailed understanding of the sources of microarray emissions after hybridization. This opportunity to simultaneously identify and quantitate contaminant and background emissions in microarrays markedly improves the reliability and accuracy of the data and permits a level of quality control of microarray emissions previously unachievable. Using these tools, we can not only quantify the extent and contribution of extraneous emission sources to the signal, but also determine the consequences of failing to account for them and gain the insight necessary to adjust preparation protocols to prevent such problems from occurring.

Genomic analysis of stationary-phase and exit in Saccharomyces cerevisiae: gene expression and identification of novel essential genes.

Most cells on earth exist in a quiescent state. In yeast, quiescence is induced by carbon starvation, and exit occurs when a carbon source becomes available. To understand how cells survive in, and exit from this state, mRNA abundance was examined using oligonucleotide-based microarrays and quantitative reverse transcription-polymerase chain reaction. Cells in stationary-phase cultures exhibited a coordinated response within 5-10 min of refeeding. Levels of >1800 mRNAs increased dramatically (>or=64-fold), and a smaller group of stationary-phase mRNAs decreased in abundance. Motif analysis of sequences upstream of genes clustered by VxInsight identified an overrepresentation of Rap1p and BUF (RPA) binding sites in genes whose mRNA levels rapidly increased during exit. Examination of 95 strains carrying deletions in stationary-phase genes induced identified 32 genes essential for survival in stationary-phase at 37 degrees C. Analysis of these genes suggests that mitochondrial function is critical for entry into stationary-phase and that posttranslational modifications and protection from oxidative stress become important later. The phylogenetic conservation of stationary-phase genes, and our findings that two-thirds of the essential stationary-phase genes have human homologues and of these, many have human homologues that are disease related, demonstrate that yeast is a bona fide model system for studying the quiescent state of eukaryotic cells.

Identification and removal of contaminating fluorescence from commercial and in-house printed DNA microarrays.

Microarray analysis is a critically important technology for genome-enabled biology, therefore it is essential that the data obtained be reliable. Current software and normalization techniques for microarray analysis rely on the assumption that fluorescent background within spots is essentially the same throughout the glass slide and can be measured by fluorescence surrounding the spots. This assumption is not valid if background fluorescence is spot-localized. Inaccurate estimates of background fluorescence under the spot create a source of error, especially for low expressed genes. We have identified spot-localized, contaminating fluorescence in the Cy3 channel on several commercial and in-house printed microarray slides. We determined through mock hybridizations (without labeled target) that pre-hybridization scans could not be used to predict the contribution of this contaminating fluorescence after hybridization because the change in spot-to-spot fluorescence after hybridization was too variable. Two solutions to this problem were identified. First, allowing 4 h of exposure to air prior to printing on to Corning UltraGAPS slides significantly reduced contaminating fluorescence intensities to approximately the value of the surrounding glass. Alternatively, application of a novel, hyperspectral imaging scanner and multivariate curve resolution algorithms, allowed the spectral contributions of Cy3 signal, glass, and contaminating fluorescence to be distinguished and quantified after hybridization.