Development of a Common Oligonucleotide Reference Standard for Microarray Data Normalization and Comparison across Different Microbial Communities

High-density functional gene arrays have become a powerful tool for environmental microbial detection and characterization. However, microarray data normalization and comparison for this type of microarray remain a challenge in environmental microbiology studies because some commonly used normalization methods (e.g., genomic DNA) for the study of pure cultures are not applicable. In this study, we developed a common oligonucleotide reference standard (CORS) method to address this problem. A unique 50-mer reference oligonucleotide probe was selected to co-spot with gene probes for each array feature. The complementary sequence was synthesized and labeled for use as the reference target, which was then spiked and cohybridized with each sample. The signal intensity of this reference target was used for microarray data normalization and comparison. The optimal amount or concentration were determined to be ca. 0.5 to 2.5% of a gene probe for the reference probe and ca. 0.25 to 1.25 fmol/μl for the reference target based on our evaluation with a pilot array. The CORS method was then compared to dye swap and genomic DNA normalization methods using the Desulfovibrio vulgaris whole-genome microarray, and significant linear correlations were observed. This method was then applied to a functional gene array to analyze soil microbial communities, and the results demonstrated that the variation of signal intensities among replicates based on the CORS method was significantly lower than the total intensity normalization method. The developed CORS provides a useful approach for microarray data normalization and comparison for studies of complex microbial communities.Microarray-based technology has become a robust genomic tool to detect, track, and profile hundreds to thousands of different microbial populations simultaneously in complex environments such as soils and sediments. For example, GeoChip, a comprehensive functional gene array, has been developed for investigating biogeochemical, ecological, and environmental processes (12, 18, 23, 27, 29, 32). Although a massive amount of microarray data can be generated rapidly, one of the bottlenecks in using microarrays for environmental microbial community studies is the lack of an appropriate standard for data comparison and normalization (6). Currently, it is difficult to compare microarray data across different sites, experiments, laboratories, and/or time periods (10). This limits the power of the technology to address ecological and environmental questions.In pure culture-based functional genomics studies, genomic DNAs (gDNAs) have been used as a common reference for hybridizations in which the same amount of gDNAs are used to cohybridize with each target cDNA sample and then to normalize different target cDNAs based on the gDNA standard (4, 5, 8, 9, 19, 21, 23). Several normalization methods such as scale normalization, quantile normalization, and Lowess normalization have been used for gene expression studies (2). Using the gDNA standard method can minimize or eliminate differences in target cDNA quantity, spot morphology, uneven hybridization, labeling, and sequence-specific hybridization behaviors (5), and this allows the comparison of microarray data across different sites, laboratories, experiments, and/or times. The main rationale for gDNA as a common reference is that it provides complete coverage for all genes represented on the array because the DNA composition from a particular organism should be identical across different treatment samples even though RNA expression is different (8). However, this approach is not applicable to microbial community studies because not all communities have identical DNA compositions. Pooling of equal amounts of gDNA or RNA from every target sample to make a common sample could be used as an alternative reference for cohybridization (1, 22). However, the disadvantage of the sample pooling approach is that samples do not provide large amounts of DNA or RNA in a reliable and reproducible way. For example, groundwater samples usually have a very low biomass and thus would not provide enough DNA for pooling. In addition, the sample pool itself is uncharacterized, and gene abundance may be diluted out so that insufficient DNA is present to result in a positive signal some array features, especially for those genes in low abundance. Moreover, a new sample pool would be required for every new experiment, making comparison across experiments difficult. Thus, other approaches need to be developed for microbial community studies.Dudley et al. (7) used a 25-mer oligonucleotide that matched a small portion of the parental EST clone vector contained in every PCR product printed on the array for normalization of pure culture RNA expression. Although the oligonucleotide generated a stable hybridization signal on every array feature, this method requires a universal sequence tag as a “capture” sequence, limiting its general use in microbial community studies. Thus, in the present study, we developed a common oligonucleotide reference standard (CORS) approach by co-spotting a common oligonucleotide with each array feature to improve the accuracy and comparability of microarray data for microbial community studies. This method was evaluated by using a pilot array, a whole-genome array, and a functional gene array, and all results demonstrate that the developed CORS is a reliable and reproducible method for microarray data normalization and comparison for microbial community studies.     

An Invasive Whitefly Feeding on a Virus-Infected Plant Increased Its Egg Production and Realized Fecundity

Background

Plant-pathogenic begomoviruses have a complex association with their insect vectors. The interactions of begomoviruses and reproduction of their vectors are poorly understood. Bemisia tabaci is known to transmit many begomoviruses, and the spread of B. tabaci, especially the B and Q ‘biotypes’, has been accompanied by the epidemics of begomoviruses. One of these identified disease-causing agents was Tomato yellow leaf curl China virus (TYLCCNV).

Methodology/Principal Findings

In this study, we compared the egg production and realized fecundity of two ‘biotypes’ or putative species of the whitefly B. tabaci, including the alien invasive B and the indigenous ZHJ1 from Zhejiang, China, feeding on either healthy or TYLCCNV-infected tobacco plants. The ovary of the whitefly was composed of 12–22 telotrophic ovarioles. According to the morphology of the oocytes and level of yolk content, oocytes in ovarioles were divided into four developmental phases (I-IV). Significantly higher proportion of immature oocytes (phase II, III) and mature oocytes (phase IV) was observed in ovary of females that fed on TYLCCNV-infected tobacco compared to that on healthy plants. Moreover, there was significant increase of eggs laid of B whitefly that fed on TYLCCNV-infected tobacco plants during the early developmental stages. In contrast, the proportion of oocytes of different developmental phases and eggs laid had no significant differences between ZHJ1 whiteflies feeding on TYLCCNV-infected and non-infected host plants.

Conclusions/Significance

The invasive B whitefly benefits from feeding on a begomovirus-infected plant through increased egg production and realized fecundity.     

Prevalence of Antibiotic Resistance in Drinking Water Treatment and Distribution Systems

The occurrence and spread of antibiotic-resistant bacteria (ARB) are pressing public health problems worldwide, and aquatic ecosystems are a recognized reservoir for ARB. We used culture-dependent methods and quantitative molecular techniques to detect and quantify ARB and antibiotic resistance genes (ARGs) in source waters, drinking water treatment plants, and tap water from several cities in Michigan and Ohio. We found ARGs and heterotrophic ARB in all finished water and tap water tested, although the amounts were small. The quantities of most ARGs were greater in tap water than in finished water and source water. In general, the levels of bacteria were higher in source water than in tap water, and the levels of ARB were higher in tap water than in finished water, indicating that there was regrowth of bacteria in drinking water distribution systems. Elevated resistance to some antibiotics was observed during water treatment and in tap water. Water treatment might increase the antibiotic resistance of surviving bacteria, and water distribution systems may serve as an important reservoir for the spread of antibiotic resistance to opportunistic pathogens.The occurrence and spread of antibiotic-resistant bacteria (ARB) are pressing public health problems worldwide, and aquatic ecosystems are a recognized reservoir for ARB and antibiotic resistance genes (ARGs) (4, 6, 8, 11, 12, 15, 39). Naturally occurring ARB and ARGs in the aquatic environment are selected for and enriched for by antibiotics found in sewage and agricultural runoff, which result from the widespread and increased use of antibiotics (4, 11, 12, 15, 38). Historically, concerns about the microbial quality of drinking water have focused on the occurrence of pathogens in drinking water distribution systems (5, 34). However, the presence of trace levels of antibiotics and ARB in source water and finished drinking water may also greatly affect public health and is an emerging issue for the general public and the drinking water industry (3, 30). Although several studies have detected ARB in drinking water systems (2, 3, 20, 30, 38), most previous studies focused on cultivable bacteria and/or indicator organisms. Little is known about the fate of ARGs in drinking water systems, and it was recently proposed that ARGs are emerging contaminants (24).We used culture-dependent methods and molecular techniques to investigate the prevalence and dynamics of heterotrophic ARB and ARGs in a drinking water source (source RW-P) and treated drinking water (source DW-P) (see Materials and Methods in the supplemental material). We tested water from a drinking water plant located in Michigan and tap water from several small cities located in Michigan and Ohio (sources TW-1, TW-2, TW-3, and TW-4). Two independent samples were collected each time at each collection site at three different times, and we used four replicates from each sample for tests. We tested bacterial resistance to the following antibiotics: amoxicillin (amoxicilline), chloramphenicol, ciprofloxacin, gentamicin, rifampin (rifampicin), sulfisoxazole, and tetracycline. We also examined the presence of eight ARGs, including beta-lactam resistance genes (bla_TEM and bla_SHV), chloramphenicol resistance genes (cat and cmr), sulfonamide resistance genes (sulI and sulII), and tetracycline resistance genes (tetO and tetW).Total heterotrophic plate counts (HPC) were determined using R₂A agar without added antibiotics. The water treatment process reduced the total HPC from 9.9 × 10⁶ CFU/100 ml in source water to 68 CFU/100 ml in treated drinking water, indicating that there was efficient removal and/or deactivation of total HPC (Table ​(Table1).1). In contrast, the total 16S rRNA gene copy number decreased from 3.4 × 10⁷ copies/100 ml in source water to 1.6 × 10⁶ copies/100 ml in treated drinking water (Fig. ​(Fig.1).1). The discrepancy between the reduction in the HPC and the reduction in the total 16S rRNA gene copy number suggests that the final disinfection step effectively inactivated bacteria but most of the dead or damaged cells were still present in finished drinking water. The number of HPC in tap water ranged from 3.44 × 10² to 6.1 × 10⁴ CFU/100 ml water, values that are lower than those for source water but significantly higher than those for treated drinking water, indicating that there is regrowth of bacteria in drinking water distribution systems. The copy numbers of total 16S rRNA genes in tap water ranged from 2.45 × 10⁵ to 1.02 × 10⁷ copies/100 ml water. The higher levels suggested by the 16S rRNA data are consistent with results of previous studies demonstrating that only 5 to 10% and 1% of bacteria in wastewater and soil, respectively, can be cultivated or identified by culture-based methods (9, 37). A significant correlation (P < 0.05, R² = 0.78) was found between the 16S rRNA gene copy number and the total HPC if treated drinking water (DW-P) data were not included (Fig. ​(Fig.1).1). This suggests that cultivable bacteria in drinking water represent only a small portion of the total bacterial biomass. Including treated drinking water (DW-P) data resulted in a distorted correlation, suggesting that a large proportion of the 16S rRNA genes present came from dead and/or damaged cells. The levels of total heterotrophic bacteria were significantly higher in tap water (TW-1) than in treated drinking water (DW-P), indicating that there was bacterial regrowth in the water distribution system.Open in a separate windowFIG. 1.Heterotrophic bacteria and the 16S rRNA gene in different water samples. (A) Copy numbers of the 16S rRNA gene and numbers of heterotrophic bacteria (CFU) in 100 ml water. (B) Correlation (P < 0.05, R² = 0.78) between the copy number of the 16S rRNA gene and the number of heterotrophic bacteria in different water samples (without the data for DW-P). RW-P, source water from the drinking water treatment plant; DW-P, finished drinking water from the drinking water treatment plant; TW-1, tap water from the city where the drinking water treatment plant is located; TW-2, TW-3, and TW-4, tap water from three towns in Michigan and Ohio close to the city where the TW-1 drinking water treatment plant is located. The statistical analysis was done using six samples for each type of water sample. Lg, log₁₀.

TABLE 1.

Prevalence of ARB HPC in source water, finished drinking water, and tap water from four towns^a

Airborne hyperspectral imaging for estimating acorn yield based on the PLS B-matrix calibration technique

Alternate bearing of acorn is a well-marked yield variability phenomenon in forest production. In Japan, this phenomenon is also related to wildlife management (e.g. of animals such as wild pigs, that rely on acorn as their major feed source). Effective management of animals dependent on acorn will require accurate estimation of acorn yield at an early stage. In this paper, we proposed a way to estimate acorn yield from the canopy reflectance values of individual trees. Using an Airborne Imaging Spectrometer for Application (AISA) Eagle System, hyperspectral images in 72 visible and near-infrared wavelengths (407–898 nm) were acquired over an acorn forest in Japan 10 times over three consecutive years (2003–2005) during the early acorn growing season. The canopy spectral reflectance values for individual trees at each wavelength were extracted from the images, and important wavelengths were determined as estimating factors by the B-matrix technique based on partial least squares (PLS) analysis. Yield-estimating models were then developed by multiple linear regression (MLR). Three models obtained from images acquired on June 27 in 2003, July 13 in 2004 and June 21 in 2005 estimated acorn yield well in comparison with ground truth, indicating that the procedure has considerable potential. The study also demonstrated the B-matrix technique based on PLS analysis to be reliable and efficient in identifying important wavelengths for determining suitable estimating factors that best contribute to the estimation model.     

The effect of mutation at M307 in <Emphasis Type="Italic">FUT1</Emphasis> gene on susceptibility of <Emphasis Type="Italic">Escherichia coli</Emphasis> F18 and gene expression in Sutai piglets

Escherichia coli F18 (ECF18) is a common porcine enteric pathogen. The pathogenicity of ECF18 bacteria depends on the existence of ECF18 receptor in the brush border membranes of piglet’s small intestinal mucosa. Alpha (1) fucosyltransferase gene (FUT1) has been identified as the candidate gene controlling the adhesion to ECF18 receptor. The genetic variations in the position of M307 nucleotide in open reading frame of FUT1 have been proposed as a marker for selecting resistant pigs. The piglets were divided into three groups, AA, AG and GG, according to the genotypes present at M307 of FUT1. Small intestinal epithelium cells of piglets with AA, AG and GG genotypes were selected to test the adhesion capability of the wild type E.coli expressing F18ab fimbriae, the recombinant E. coli expressing F18ac fimbriae or the recombinant E. coli secreting and surface-displaying the FedF subunit of F18ab fimbriae, respectively. Here, we examined the distribution and expression of porcine FUT1 mRNA in different tissues in Sutai pigs using real-time PCR. The results showed that piglets with AA genotype show resistance, whereas piglets with GG or AG genotypes are sensitive to the pathogenic E. coli F18 in Sutai piglets. FUT1 was expressed in all the tissues that were examined, and the gene’s expression was highest in the lungs. There was no significant difference in expression level among the three genotypes in the liver, lung, stomach and duodenum, where the gene expression was relatively high. The present analysis suggested that mutation at M307 in FUT1 gene determines susceptibility of small intestinal epithelium to E. coli F18 adhesion in Sutai piglet and the expression of FUT1 gene may be regulated by other factors or the mutation was likely to be in linkage disequilibrium with some cis-regulatory variants.     

RyR2 and calpain-10 delineate a novel apoptosis pathway in pancreatic islets

Cells are programmed to die when critical signaling and metabolic pathways are disrupted. Inhibiting the type 2 ryanodine receptor (RyR2) in human and mouse pancreatic beta-cells markedly increased apoptosis. This mode of programmed cell death was not associated with robust caspase-3 activation prompting a search for an alternative mechanism. Increased calpain activity and calpain gene expression suggested a role for a calpain-dependent death pathway. Using a combination of pharmacological and genetic approaches, we demonstrated that the calpain-10 isoform mediated ryanodine-induced apoptosis. Apoptosis induced by the fatty acid palmitate and by low glucose also required calpain-10. Ryanodine-induced calpain activation and apoptosis were reversed by glucagon-like peptide or short-term exposure to high glucose. Thus RyR2 activity seems to play an essential role in beta-cell survival in vitro by suppressing a death pathway mediated by calpain-10, a type 2 diabetes susceptibility gene with previously unknown function.     

PKCα translocation from mitochondria to nucleus is closely related to induction of apoptosis in gastric cancer cells

PKCs have been implicated in the regulation of cellular differentiation, proliferation, apoptosis and signal transduction. It was demonstrated in this study that PKCα was located both at mitochondria and in cytosol in gastric cancer cell line BGC-823. Treatment of cells with 12-O-tetradecanoylphorbol-13-acetate (TPA) resulted in the translocation of PKCα from both mitochondria and cytosol to nucleus as clearly shown by laserscanningconfocal microscopy, while the protein level of PKCα was not changed by TPA treatment as detected by Western blot. The results also revealed that TPA-induced translocation of PKCα was in close association with apoptosis induction, and such association was further affirmed by other experiments where various apoptotic stimuli and specific inhibitors of PKC were used. Taken together, these findings indicate that translocation of PKCα from both mitochondria and cytosol to nucleus in gastric cancer cell is accompanied by induction of apoptosis, and may imply a new mechanism of the potential linking between cell apoptosis and PKCα translocation.     

Specific inhibition effects of N-pentafluorobenzyl-1-deoxynojirimycin on human CD4+ T cells

We first synthesized N-pentafluorobenzyl-1-deoxynojirimycin (5F-DNM), one new derivative of 1-deoxynojirimycin (DNM). Effects on human peripheral blood mononuclear cells (PMBC) and secretion of cytokines from human PBMC by 5F-DNM were investigated. It was first found that 5F-DNM remarkably inhibited the secretion of interleukin-4 (IL-4) and had a specific inhibition on the expression of CD4 molecules. 5F-DNM, much less toxic than cyclosporin A, might be used as a new candidate of immunosuppressant for specifically treating Th2-mediated immune diseases.