Low relative humidity triggers RNA-directed de novo DNA methylation and suppression of genes controlling stomatal development

Environmental cues influence the development of stomata on the leaf epidermis, and allow plants to exert plasticity in leaf stomatal abundance in response to the prevailing growing conditions. It is reported that Arabidopsis thaliana ‘Landsberg erecta’ plants grown under low relative humidity have a reduced stomatal index and that two genes in the stomatal development pathway, SPEECHLESS and FAMA, become de novo cytosine methylated and transcriptionally repressed. These environmentally-induced epigenetic responses were abolished in mutants lacking the capacity for de novo DNA methylation, for the maintenance of CG methylation, and in mutants for the production of short-interfering non-coding RNAs (siRNAs) in the RNA-directed DNA methylation pathway. Induction of methylation was quantitatively related to the induction of local siRNAs under low relative humidity. Our results indicate the involvement of both transcriptional and post-transcriptional gene suppression at these loci in response to environmental stress. Thus, in a physiologically important pathway, a targeted epigenetic response to a specific environmental stress is reported and several of its molecular, mechanistic components are described, providing a tractable platform for future epigenetics experiments. Our findings suggest epigenetic regulation of stomatal development that allows for anatomical and phenotypic plasticity, and may help to explain at least some of the plant’s resilience to fluctuating relative humidity.

Cell walls of developing wheat starchy endosperm: comparison of composition and RNA-Seq transcriptome

The transcriptome of the developing starchy endosperm of hexaploid wheat (Triticum aestivum) was determined using RNA-Seq isolated at five stages during grain fill. This resource represents an excellent way to identify candidate genes responsible for the starchy endosperm cell wall, which is dominated by arabinoxylan (AX), accounting for 70% of the cell wall polysaccharides, with 20% (1,3; 1,4)-beta-D-glucan, 7% glucomannan, and 4% cellulose. A complete inventory of transcripts of 124 glycosyltransferase (GT) and 72 glycosylhydrolase (GH) genes associated with cell walls is presented. The most highly expressed GT transcript (excluding those known to be involved in starch synthesis) was a GT47 family transcript similar to Arabidopsis (Arabidopsis thaliana) IRX10 involved in xylan extension, and the second most abundant was a GT61. Profiles for GT43 IRX9 and IRX14 putative orthologs were consistent with roles in AX synthesis. Low abundances were found for transcripts from genes in the acyl-coA transferase BAHD family, for which a role in AX feruloylation has been postulated. The relative expression of these was much greater in whole grain compared with starchy endosperm, correlating with the levels of bound ferulate. Transcripts associated with callose (GSL), cellulose (CESA), pectin (GAUT), and glucomannan (CSLA) synthesis were also abundant in starchy endosperm, while the corresponding cell wall polysaccharides were confirmed as low abundance (glucomannan and callose) or undetectable (pectin) in these samples. Abundant transcripts from GH families associated with the hydrolysis of these polysaccharides were also present, suggesting that they may be rapidly turned over. Abundant transcripts in the GT31 family may be responsible for the addition of Gal residues to arabinogalactan peptide.

Lexical access, story re-telling and sequencing skills in adults who clutter and those who do not

Cluttering is a rate-based disorder of fluency, the scope of whose diagnostic criteria currently remains unclear. This paper reports preliminary findings from a larger study which aims to determine whether cluttering can be associated with language disturbances as well as motor and rate based ones. Subtests from the Mt Wilga High Level Language Test (MWHLLT) were used to determine whether people who clutter (PWC) have word finding difficulties, and use significantly more maze behaviours compared to controls, during story re-telling and simple sequencing tasks. Independent t tests showed that PWC were significantly slower than control participants in lexical access and sentence completion tasks, but returned mixed findings when PWCs were required to name items within a semantic category. PWC produced significantly more maze behaviour than controls in a task where participants were required to explain how to undertake commonly performed actions, but no difference in use of maze behaviour was found between the two groups when retelling a story from memory. The implications of these findings are discussed

Partial identification of spirochaetes from two dairy cows with digital dermatitis by polymerase chain reaction analysis of the 16S ribosomal RNA gene

Specimens taken postmortem from typical lesions of digital dermatitis in two dairy cows were tested by the polymerase chain reaction (PCR) for the presence of a spirochaetal 16S rRNA gene. Seven different assays detected the gene in the samples from both cows. Two of the PCR products were sequenced and a comparison of the nucleotide sequences revealed that the spirochaete belonged to the genus Treponema and was closely related to Treponema denticola. A PCR specific for the detection of the digital dermatitis-associated treponeme was developed.

A Brassica exon array for whole-transcript gene expression profiling

Affymetrix GeneChip (R) arrays are used widely to study transcriptional changes in response to developmental and environmental stimuli. GeneChip (R) arrays comprise multiple 25-mer oligonucleotide probes per gene and retain certain advantages over direct sequencing. For plants, there are several public GeneChip (R) arrays whose probes are localised primarily in 39 exons. Plant whole-transcript (WT) GeneChip (R) arrays are not yet publicly available, although WT resolution is needed to study complex crop genomes such as Brassica, which are typified by segmental duplications containing paralogous genes and/or allopolyploidy. Available sequence data were sampled from the Brassica A and C genomes, and 142,997 gene models identified. The assembled gene models were then used to establish a comprehensive public WT exon array for transcriptomics studies. The Affymetrix GeneChip (R) Brassica Exon 1.0 ST Array is a 5 mu M feature size array, containing 2.4 million 25-base oligonucleotide probes representing 135,201 gene models, with 15 probes per gene distributed among exons. Discrimination of the gene models was based on an E-value cut-off of 1E(-5), with <= 98 sequence identity. The 135 k Brassica Exon Array was validated by quantifying transcriptome differences between leaf and root tissue from a reference Brassica rapa line (R-o-18), and categorisation by Gene Ontologies (GO) based on gene orthology with Arabidopsis thaliana. Technical validation involved comparison of the exon array with a 60-mer array platform using the same starting RNA samples. The 135 k Brassica Exon Array is a robust platform. All data relating to the array design and probe identities are available in the public domain and are curated within the BrassEnsembl genome viewer at http://www.brassica.info/BrassEnsembl/index.html.

RNA interference suppression of genes in glycosyl transferase families 43 and 47 in wheat starchy endosperm causes large decreases in arabinoxylan content

The cell walls of wheat (Triticum aestivum) starchy endosperm are dominated by arabinoxylan (AX), accounting for 65% to 70% of the polysaccharide content. Genes within two glycosyl transferase (GT) families, GT43 (IRREGULAR XYLEM9 [IRX9] and IRX14) and GT47 (IRX10), have previously been shown to be involved in the synthesis of the xylan backbone in Arabidopsis, and close homologs of these have been implicated in the synthesis of xylan in other species. Here, homologs of IRX10 TaGT47_2 and IRX9 TaGT43_2, which are highly expressed in wheat starchy endosperm cells, were suppressed by RNA interference (RNAi) constructs driven by a starchy endosperm-specific promoter. The total amount of AX was decreased by 40% to 50% and the degree of arabinosylation was increased by 25% to 30% in transgenic lines carrying either of the transgenes. The cell walls of starchy endosperm in sections of grain from TaGT43_2 and TaGT47_2 RNAi transgenics showed decreased immunolabeling for xylan and arabinoxylan epitopes and approximately 50% decreased cell wall thickness compared with controls. The proportion of AX that was water soluble was not significantly affected, but average AX polymer chain length was decreased in both TaGT43_2 and TaGT47_2 RNAi transgenics. However, the long AX chains seen in controls were absent in TaGT43_2 RNAi transgenics but still present in TaGT47_2 RNAi transgenics. The results support an emerging picture of IRX9-like and IRX10-like proteins acting as key components in the xylan synthesis machinery in both dicots and grasses. Since AX is the main component of dietary fiber in wheat foods, the TaGT43_2 and TaGT47_2 genes are of major importance to human nutrition.

The mysterious presence of a 5-methylcytosine oxidase in the Drosophila genome: Possible explanations

5-methylcytosine is an important epigenetic modification involved in gene control in vertebrates and many other complex living organisms. Its presence in Drosophila has been a matter of debate and recent bisulfite sequencing studies of early-stage fly embryos have concluded that the genome of Drosophila is essentially unmethylated. However, as we outline here, the Drosophila genome harbors a well-conserved homolog of the TET protein family. The mammalian orthologs TET1/2/3 are known to convert 5-methylcytosine into 5-hydroxymethylcytosine. We discuss several possible explanations for these seemingly contradictory findings. One possibility is that the 2 modified cytosine bases are generated in Drosophila only at certain developmental stages and in a cell type-specific manner during neurogenesis. Alternatively, Drosophila Tet and its mammalian homologs may carry out catalytic activity-independent functions, and the possibility that these proteins may oxidize 5-methylcytosine in RNA created by the methyltransferase Dnmt2 should also be strongly considered.

How the double spherules of infectious bronchitis virus impact our understanding of RNA virus replicative organelles

Powered by advances in electron tomography, recent studies have extended our understanding of how viruses construct "replication factories" inside infected cells. Their function, however, remains an area of speculation with important implications for human health. It is clear from these studies that whatever their purpose, organelle structure is dynamic (M. Ulasli, M. H. Verheije, C. A. de Haan, and F. Reggiori, Cell. Microbiol. 12:844-861, 2010) and intricate (K. Knoops, M. Kikkert, S. H. Worm, J. C. Zevenhoven-Dobbe, Y. van der Meer, et al., PLOS Biol. 6:e226, 2008). But by concentrating on medically important viruses, these studies have failed to take advantage of the genetic variation inherent in a family of viruses that is as diverse as the archaea, bacteria, and eukaryotes combined (C. Lauber, J. J. Goeman, M. del Carmen Parquet, P. T. Nga, E. J. Snijder, et al., PLOS Pathog. 9:e1003500, 2013). In this climate, Maier et al. (H. J. Maier, P. C. Hawes, E. M. Cottam, J. Mantell, P. Verkade, et al., mBio 4:e00801-13, 2013) explored the replicative structures formed by an avian coronavirus that appears to have diverged at an early point in coronavirus evolution and shed light on controversial aspects of viral biology.

Prediction of Staphylococcus aureus antimicrobial resistance by whole-genome sequencing

Whole-genome sequencing (WGS) could potentially provide a single platform for extracting all the information required to predict an organism’s phenotype. However, its ability to provide accurate predictions has not yet been demonstrated in large independent studies of specific organisms. In this study, we aimed to develop a genotypic prediction method for antimicrobial susceptibilities. The whole genomes of 501 unrelated Staphylococcus aureus isolates were sequenced, and the assembled genomes were interrogated using BLASTn for a panel of known resistance determinants (chromosomal mutations and genes carried on plasmids). Results were compared with phenotypic susceptibility testing for 12 commonly used antimicrobial agents (penicillin, methicillin, erythromycin, clindamycin, tetracycline, ciprofloxacin, vancomycin, trimethoprim, gentamicin, fusidic acid, rifampin, and mupirocin) performed by the routine clinical laboratory. We investigated discrepancies by repeat susceptibility testing and manual inspection of the sequences and used this information to optimize the resistance determinant panel and BLASTn algorithm. We then tested performance of the optimized tool in an independent validation set of 491 unrelated isolates, with phenotypic results obtained in duplicate by automated broth dilution (BD Phoenix) and disc diffusion. In the validation set, the overall sensitivity and specificity of the genomic prediction method were 0.97 (95% confidence interval [95% CI], 0.95 to 0.98) and 0.99 (95% CI, 0.99 to 1), respectively, compared to standard susceptibility testing methods. The very major error rate was 0.5%, and the major error rate was 0.7%. WGS was as sensitive and specific as routine antimicrobial susceptibility testing methods. WGS is a promising alternative to culture methods for resistance prediction in S. aureus and ultimately other major bacterial pathogens.

Molecular archaeology of Flaviviridae untranslated regions: duplicated RNA structures in the replication enhancer of flaviviruses and pestiviruses emerged via convergent evolution.

RNA secondary structures in the 3'untranslated regions (3'UTR) of the viruses of the family Flaviviridae, previously identified as essential (promoters) or beneficial (enhancers) for replication, have been analysed. Duplicated enhancer elements are revealed as a global feature in the evolution of the 3'UTR of distantly related viruses within the genera Flavivirus and Pestivirus. For the flaviviruses, duplicated structures occur in the 3'UTR of all four distantly related ecological virus subgroups (tick-borne, mosquito-borne, no known vector and insect-specific flaviviruses (ISFV). RNA structural differences distinguish tick-borne flaviviruses with discrete pathogenetic characteristics. For Aedes- and Culex-associated ISFV, secondary RNA structures with different conformations display numerous short ssRNA direct repeats, exposed as loops and bulges. Long quadruplicate regions comprise almost the entire 3'UTR of Culex-associated ISFV. Extended duplicated sequence and associated RNA structures were also discovered in the 3'UTR of pestiviruses. In both the Flavivirus and Pestivirus genera, duplicated RNA structures were localized to the enhancer regions of the 3'UTR suggesting an adaptive role predominantly in wild-type viruses. We propose sequence reiteration might act as a scaffold for dimerization of proteins involved in assembly of viral replicase complexes. Numerous nucleotide repeats exposed as loops/bulges might also interfere with host immune responses acting as a molecular sponge to sequester key host proteins or microRNAs.

An intronic G run within HIV-1 intron 2 is critical for splicing regulation of vif mRNA

Within target T lymphocytes, human immunodeficiency virus type I (HIV-1) encounters the retroviral restriction factor APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; A3G), which is counteracted by the HIV-1 accessory protein Vif. Vif is encoded by intron-containing viral RNAs that are generated by splicing at 3' splice site (3'ss) A1 but lack splicing at 5'ss D2, which results in the retention of a large downstream intron. Hence, the extents of activation of 3'ss A1 and repression of D2, respectively, determine the levels of vif mRNA and thus the ability to evade A3G-mediated antiviral effects. The use of 3'ss A1 can be enhanced or repressed by splicing regulatory elements that control the recognition of downstream 5'ss D2. Here we show that an intronic G run (G(I2)-1) represses the use of a second 5'ss, termed D2b, that is embedded within intron 2 and, as determined by RNA deep-sequencing analysis, is normally inefficiently used. Mutations of G(I2)-1 and activation of D2b led to the generation of transcripts coding for Gp41 and Rev protein isoforms but primarily led to considerable upregulation of vif mRNA expression. We further demonstrate, however, that higher levels of Vif protein are actually detrimental to viral replication in A3G-expressing T cell lines but not in A3G-deficient cells. These observations suggest that an appropriate ratio of Vif-to-A3G protein levels is required for optimal virus replication and that part of Vif level regulation is effected by the novel G run identified here.

Superinfection exclusion and the long-term survival of honey bees in Varroa-infested colonies

Over the past 50 years, many millions of European honey bee (Apis mellifera) colonies have died as the ectoparasitic mite, Varroa destructor, has spread around the world. Subsequent studies have indicated that the mite's association with a group of RNA viral pathogens (Deformed Wing Virus, DWV) correlates with colony death. Here, we propose a phenomenon known as superinfection exclusion that provides an explanation of how certain A. mellifera populations have survived, despite Varroa infestation and high DWV loads. Next-generation sequencing has shown that a non-lethal DWV variant 'type B' has become established in these colonies and that the lethal 'type A' DWV variant fails to persist in the bee population. We propose that this novel stable host-pathogen relationship prevents the accumulation of lethal variants, suggesting that this interaction could be exploited for the development of an effective treatment that minimises colony losses in the future.The ISME Journal advance online publication, 27 October 2015; doi:10.1038/ismej.2015.186.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

Soil bacterial phoD gene abundance and expression in response to applied phosphorus and long-term management

Bacterial transformation of phosphorus (P) compounds in soil is largely dependent on soil microbial community function, and is therefore sensitive to anthropogenic disturbances such as fertilization or cropping systems. However, the effect of soil management on the transcription of bacterial genes that encode phosphatases, such as phoD, is largely unknown. This greenhouse study examined the effect of long-term management and P amendment on potential alkaline phosphatase (ALP) activity and phoD gene (DNA) and transcript (RNA) abundance. Soil samples (0–15 cm) were collected from the Glenlea Long-term Rotation near Winnipeg, Manitoba, to compare organic, conventional and prairie management systems. In the greenhouse, pots of soil from each management system were amended with P as either soluble mineral fertilizer or cattle manure and then planted with Italian ryegrass (Lolium multiforum). Soils from each pot were sampled for analysis immediately and after 30 and 106 days. Significant differences among the soil/P treatments were detected for inorganic P, but not the organic P in NaHCO3-extracts. At day 0, ALP activity was similar among the soil/P treatments, but was higher after 30 days for all P amendments in soil from organically managed plots. In contrast, ALP activity in soils under conventional and prairie management responded to increasing rates of manure only, with significant effects from medium and high manure application rates at 30 and 106 days. Differences in ALP activity at 30 days corresponded to the abundance of bacterial phoD genes, which were also significantly higher in soils under organic management. However, this correlation was not significant for transcript abundance. Next-generation sequencing allowed the identification of 199 unique phoD operational taxonomic units (OTUs) from the metagenome (soil DNA) and 35 unique OTUs from the metatranscriptome (soil RNA), indicating that a subset of phoD genes was being transcribed in all soils.