Novel cis-acting regulatory elements in wild Oryza species impart improved rice bran quality by lowering the expression of phospholipase D alpha1 enzyme (OsPLDα1).

Rice bran oil is good quality edible oil, rich in antioxidants and comprised typically of oleic-linoleic type fatty acids. However, presence of a highly lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) increases free fatty acid content in the oil which further leads to stale flavor and rancidity of the oil, making it unfit for human consumption. In this study, we compared the upstream regions of OsPLDα1 orthologs across 34 accessions representing 5 wild Oryza species and 8 cultivars, to uncover sequence variations and identify cis-elements involved in differential transcription of orthologs. Alignment of the upstream sequences to the Nipponbare OsPLDα1 reference sequence revealed the presence of 39 SNPs. Phylogenetic analysis showed that all the selected cultivars and wild species accessions are closely related to the reference except for three accessions of O. rufipogon (IRGC89224, IRGC104425, and IRGC105902). Furthermore, using exon-specific qRT-PCR, OsPLDα1 expression patterns in immature grains indicated significant differences in transcript abundance between the wild species accessions. In comparison to the control, lowest gene expression was observed in IRGC89224 accession (0.20-fold) followed by IRGC105902 (0.26-fold) and IRGC104425 (0.41-fold) accessions. In-silico analysis of the OsPLDα1 promoter revealed that the copy number variations of CGCGBOXAT, GT1CONSENSUS, IBOXCORE, NODCON2GM, OSE2ROOTNODULE, SURECOREATSULTR11, and SORLIP1AT cis-elements play an important role in the transcriptional activities of orthologous genes. Owing to the presence of ARFAT and SEBF elements only in the IRGC89224 accession, which had the lowest gene expression as well, these putative upstream regulatory sequences have been identified as novel cis-elements which may act as repressors in regulating the OsPLDα1 gene expression. The accessions identified with low OsPLDα1 expressions could be further deployed as potential donors of ideal OsPLDα1 allele for transfer of the desired trait into elite rice cultivars.

Implications of miR166 and miR159 induction to the basal response mechanisms of an andigena potato (Solanum tuberosum subsp. andigena) to salinity stress, predicted from network models in Arabidopsis.

MicroRNA (miRNA) mediated changes in gene expression by post-transcriptional modulation of major regulatory transcription factors is a potent mechanism for integrating growth and stress-related responses. Exotic plants including many traditional varieties of Andean potatoes (Solanum tuberosum subsp. andigena) are known for better adaptation to marginal environments. Stress physiological studies confirmed earlier reports on the salinity tolerance potentials of certain andigena cultivars. Guided by the hypothesis that certain miRNAs play important roles in growth modulation under suboptimal conditions, we identified and characterized salinity stress-responsive miRNA-target gene pairs in the andigena cultivar Sullu by parallel analysis of noncoding and coding RNA transcriptomes. Inverse relationships were established by the reverse co-expression between two salinity stress-regulated miRNAs (miR166, miR159) and their target transcriptional regulators HD-ZIP-Phabulosa/Phavulota and Myb101, respectively. Based on heterologous models in Arabidopsis, the miR166-HD-ZIP-Phabulosa/Phavulota network appears to be involved in modulating growth perhaps by mediating vegetative dormancy, with linkages to defense-related pathways. The miR159-Myb101 network may be important for the modulation of vegetative growth while also controlling stress-induced premature transition to reproductive phase. We postulate that the induction of miR166 and miR159 under salinity stress represents important network hubs for balancing gene expression required for basal growth adjustments.

Metagenomic investigation reveals bacteriophage-mediated horizontal transfer of antibiotic resistance genes in microbial communities of an organic agricultural ecosystem.

Agricultural microbiomes are major reservoirs of antibiotic resistance genes (ARGs), posing continuous risks to human health. To understand the role of bacteriophages as vehicles for the horizontal transfer of ARGs in the agricultural microbiome, we investigated the diversity of bacterial and viral microbiota from fecal and environmental samples on an organic farm. The profiles of the microbiome indicated the highest abundance of Bacteroidetes, Firmicutes, and Proteobacteria phyla in animal feces, with varying Actinobacteria and Spirochaetes abundance across farm animals. The most predominant composition in environmental samples was the phylum Proteobacteria. Compared to the microbiome profiles, the trends in virome indicated much broader diversity with more specific signatures between the fecal and environmental samples. Overall, viruses belonging to the order Caudovirales were the most prevalent across the agricultural samples. Additionally, the similarities within and between fecal and environmental components of the agricultural environment based on ARG-associated bacteria alone were much lower than those of total microbiome composition. However, there were significant similarities in the profiles of ARG-associated viruses across the fecal and environmental components. Moreover, the predictive models of phage-bacterial interactions on bipartite ARG transfer networks indicated that phages belonging to the order Caudovirales, particularly in the Siphoviridae family, contained diverse ARG types in different samples. Their interaction with various bacterial hosts further implied the important role of bacteriophages in ARG transmission across bacterial populations. Our findings provided a novel insight into the potential mechanisms of phage-mediated ARG transmission and their correlation with resistome evolution in natural agricultural environments. IMPORTANCE Antibiotic resistance has become a serious health concern worldwide. The potential impact of viruses, bacteriophages in particular, on spreading antibiotic resistance genes is still controversial due to the complexity of bacteriophage-bacterial interactions within diverse environments. In this study, we determined the microbiome profiles and the potential antibiotic resistance gene (ARG) transfer between bacterial and viral populations in different agricultural samples using a high-resolution analysis of the metagenomes. The results of this study provide compelling genetic evidence for ARG transfer through bacteriophage-bacteria interactions, revealing the inherent risks associated with bacteriophage-mediated ARG transfer across the agricultural microbiome.

Cis-regulatory signatures of orthologous stress-associated bZIP transcription factors from rice, sorghum and Arabidopsis based on phylogenetic footprints.

BACKGROUND: The potential contribution of upstream sequence variation to the unique features of orthologous genes is just beginning to be unraveled. A core subset of stress-associated bZIP transcription factors from rice (Oryza sativa) formed ten clusters of orthologous groups (COG) with genes from the monocot sorghum (Sorghum bicolor) and dicot Arabidopsis (Arabidopsis thaliana). The total cis-regulatory information content of each stress-associated COG was examined by phylogenetic footprinting to reveal ortholog-specific, lineage-specific and species-specific conservation patterns. RESULTS: The most apparent pattern observed was the occurrence of spatially conserved 'core modules' among the COGs but not among paralogs. These core modules are comprised of various combinations of two to four putative transcription factor binding site (TFBS) classes associated with either developmental or stress-related functions. Outside the core modules are specific stress (ABA, oxidative, abiotic, biotic) or organ-associated signals, which may be functioning as 'regulatory fine-tuners' and further define lineage-specific and species-specific cis-regulatory signatures. Orthologous monocot and dicot promoters have distinct TFBS classes involved in disease and oxidative-regulated expression, while the orthologous rice and sorghum promoters have distinct combinations of root-specific signals, a pattern that is not particularly conserved in Arabidopsis. CONCLUSIONS: Patterns of cis-regulatory conservation imply that each ortholog has distinct signatures, further suggesting that they are potentially unique in a regulatory context despite the presumed conservation of broad biological function during speciation. Based on the observed patterns of conservation, we postulate that core modules are likely primary determinants of basal developmental programming, which may be integrated with and further elaborated by additional intrinsic or extrinsic signals in conjunction with lineage-specific or species-specific regulatory fine-tuners. This synergy may be critical for finer-scale spatio-temporal regulation, hence unique expression profiles of homologous transcription factors from different species with distinct zones of ecological adaptation such as rice, sorghum and Arabidopsis. The patterns revealed from these comparisons set the stage for further empirical validation by functional genomics.

DECUSSATE network with flowering genes explains the variable effects of qDTY12.1 to rice yield under drought across genetic backgrounds.

The impact of qDTY12.1 in maintaining yield under drought has not been consistent across genetic backgrounds. We hypothesized that synergism or antagonism with additive-effect peripheral genes across the background genome either enhances or undermines its full potential. By modeling the transcriptional networks across sibling qDTY12.1-introgression lines with contrasting yield under drought (LPB = low-yield penalty; HPB = high-yield penalty), the qDTY12.1-encoded DECUSSATE gene (OsDEC) was revealed as the core of a synergy with other genes in the genetic background. OsDEC is expressed in flag leaves and induced by progressive drought at booting stage in LPB but not in HPB. The unique OsDEC signature in LPB is coordinated with 35 upstream and downstream peripheral genes involved in floral development through the cytokinin signaling pathway. Results support the differential network rewiring effects through genetic coupling-uncoupling between qDTY12.1 and other upstream and downstream peripheral genes across the distinct genetic backgrounds of LPB and HPB. The functional DEC-network in LPB defines a mechanism for early flowering as a means for avoiding the drought-induced depletion of photosynthate needed for reproductive growth. Its impact is likely through the timely establishment of stronger source-sink dynamics that sustains a robust reproductive transition under drought.

Responses of Escherichia coli and Listeria monocytogenes to ozone treatment on non-host tomato: Efficacy of intervention and evidence of induced acclimation.

Because of the continuous rise of foodborne illnesses caused by the consumption of raw fruits and vegetables, effective post-harvest anti-microbial strategies are necessary. The aim of this study was to evaluate the anti-microbial efficacy of ozone (O3) against two common causes of fresh produce contamination, the Gram-negative Escherichia coli O157:H7 and Gram-positive Listeria monocytogenes, and to relate its effects to potential mechanisms of xenobiosis by transcriptional network modeling. The study on non-host tomato environment correlated the dose × time aspects of xenobiosis by examining the correlation between bacterial survival in terms of log-reduction and defense responses at the level of gene expression. In E. coli, low (1 µg O3/g of fruit) and moderate (2 µg O3/g of fruit) doses caused insignificant reduction in survival, while high dose (3 µg/g of fruit) caused significant reduction in survival in a time-dependent manner. In L. monocytogenes, moderate dose caused significant reduction even with short-duration exposure. Distinct responses to O3 xenobiosis between E. coli and L. monocytogenes are likely related to differences in membrane and cytoplasmic structure and components. Transcriptome profiling by RNA-Seq showed that primary defenses in E. coli were attenuated after exposure to a low dose, while the responses at moderate dose were characterized by massive upregulation of pathogenesis and stress-related genes, which implied the activation of defense responses. More genes were downregulated during the first hour at high dose, with a large number of such genes getting significantly upregulated after 2 hr and 3 hr. This trend suggests that prolonged exposure led to potential adaptation. In contrast, massive downregulation of genes was observed in L. monocytogenes regardless of dose and exposure duration, implying a mechanism of defense distinct from that of E. coli. The nature of bacterial responses revealed by this study should guide the selection of xenobiotic agents for eliminating bacterial contamination on fresh produce without overlooking the potential risks of adaptation.

Novel and Transgressive Salinity Tolerance in Recombinant Inbred Lines of Rice Created by Physiological Coupling-Uncoupling and Network Rewiring Effects.

The phenomenon of transgressive segregation, where a small minority of recombinants are outliers relative to the range of parental phenotypes, is commonly observed in plant breeding populations. While this phenomenon has been attributed to complementation and epistatic effects, the physiological and developmental synergism involved have not been fully illuminated by the QTL mapping approach alone, especially for stress-adaptive traits involving highly complex interactions. By systems-level profiling of the IR29 × Pokkali recombinant inbred population of rice, we addressed the hypothesis that novel salinity tolerance phenotypes are created by reconfigured physiological networks due to positive or negative coupling-uncoupling of developmental and physiological attributes of each parent. Real-time growth and hyperspectral profiling distinguished the transgressive individuals in terms of stress penalty to growth. Non-parental network signatures that led to either optimal or non-optimal integration of developmental with stress-related mechanisms were evident at the macro-physiological, biochemical, metabolic, and transcriptomic levels. Large positive net gain in super-tolerant progeny was due to ideal complementation of beneficial traits while shedding antagonistic traits. Super-sensitivity was explained by the stacking of multiple antagonistic traits and loss of major beneficial traits. The synergism uncovered by the phenomics approach in this study supports the modern views of the Omnigenic Theory, emphasizing the synergy or lack thereof between core and peripheral components. This study also supports a breeding paradigm rooted on genomic modeling from multi-dimensional genetic, physiological, and phenotypic profiles to create novel adaptive traits for new crop varieties of the 21st century.

Contributions of Adaptive Plant Architecture to Transgressive Salinity Tolerance in Recombinant Inbred Lines of Rice: Molecular Mechanisms Based on Transcriptional Networks.

Genetic novelties are important nucleators of adaptive speciation. Transgressive segregation is a major mechanism that creates genetic novelties with morphological and developmental attributes that confer adaptive advantages in certain environments. This study examined the morpho-developmental and physiological profiles of recombinant inbred lines (RILs) from the salt-sensitive IR29 and salt-tolerant Pokkali rice, representing the total range of salt tolerance including the outliers at both ends of the spectrum. Morpho-developmental and physiological profiles were integrated with a hypothesis-driven interrogation of mRNA and miRNA transcriptomes to uncover the critical genetic networks that have been rewired for novel adaptive architecture. The transgressive super-tolerant FL510 had a characteristic small tiller angle and wider, more erect, sturdier, and darker green leaves. This unique morphology resulted in lower transpiration rate, which also conferred a special ability to retain water more efficiently for osmotic avoidance. The unique ability for water retention conferred by such adaptive morphology appeared to enhance the efficacy of defenses mediated by Na+ exclusion mechanism (SalTol-effects) inherited from Pokkali. The super-tolerant FL510 and super-sensitive FL499 had the smallest proportions of differentially expressed genes with little overlaps. Genes that were steadily upregulated in FL510 comprised a putative cytokinin-regulated genetic network that appeared to maintain robust growth under salt stress through well-orchestrated cell wall biogenesis and cell expansion, likely through major regulatory (OsRR23, OsHK5) and biosynthetic (OsIPT9) genes in the cytokinin signaling pathway. Meanwhile, a constitutively expressed cluster in FL510 prominently featured two transcription factors (OsIBH1, TAC3) that control tiller angle and growth habit through the brassinosteroid signaling pathway. Both the putative cytokinin-mediated and brassinosteroid-mediated clusters appeared to function as highly coordinated network synergies in FL510. In contrast, both networks appeared to be sub-optimal and inferior in the other RILs and parents as they were disjointed and highly fragmented. Transgressively expressed miRNAs (miR169, miR397, miR827) were also identified as prominent signatures of FL510, with functional implications to mechanisms that support robust growth, homeostasis, and osmotic stress avoidance. Results of this study demonstrate how genetic recombination creates novel morphology that complements inducible defenses hence transgressive adaptive phenotypes.

Novel allelic variation in the Phospholipase D alpha1 gene (OsPLDα1) of wild Oryza species implies to its low expression in rice bran.

Rice bran, a by-product after milling, is a rich source of phytonutrients like oryzanols, tocopherols, tocotrienols, phytosterols, and dietary fibers. Moreover, exceptional properties of the rice bran oil make it unparalleled to other vegetable oils. However, a lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) causes rancidity and 'stale flavor' in the oil, and thus limits the rice bran usage for human consumption. To improve the rice bran quality, sequence based allele mining at OsPLDα1 locus (3.6 Kb) was performed across 48 accessions representing 11 wild Oryza species, 8 accessions of African cultivated rice, and 7 Oryza sativa cultivars. From comparative sequence analysis, 216 SNPs and 30 InDels were detected at the OsPLDα1 locus. Phylogenetic analysis revealed 20 OsPLDα1 cDNA variants which further translated into 12 protein variants. The O. officinalis protein variant, when compared to Nipponbare, showed maximum variability comprising 22 amino acid substitutions and absence of two peptides and two ß-sheets. Further, expression profiling indicated significant differences in transcript abundance within as well as between the OsPLDα1 variants. Also, a new OsPLDα1 transcript variant having third exon missing in it, Os01t0172400-06, has been revealed. An O. officinalis accession (IRGC101152) had lowest gene expression which suggests the presence of novel allele, named as OsPLDα1-1a (GenBank accession no. MF966931). The identified novel allele could be further deployed in the breeding programs to overcome rice bran rancidity in elite cultivars.

Double-Locking Mechanism of Self-Compatibility in Arabidopsis thaliana: The Synergistic Effect of Transcriptional Depression and Disruption of Coding Region in the Male Specificity Gene.

Self-compatibility in Arabidopsis thaliana represents the relatively recent disruption of ancestral obligate cross pollination, recognized as one of the prevalent evolutionary pathways in flowering plants, as noted by Darwin. Our previous study found that inversion of the male specificity gene (SP11/SCR) disrupted self-incompatibility, which was restored by overexpressing the SCR with the reversed inversion. However, SCR in A. thaliana has other mutations aside from the pivotal inversion, in both promoter and coding regions, with probable effects on transcriptional regulation. To examine the functional consequences of these mutations, we conducted reciprocal introductions of native promoters and downstream sequences from orthologous loci of self-compatible A. thaliana and self-incompatible A. halleri. Use of this inter-species pair enabled us to expand the scope of the analysis to transcriptional regulation and deletion in the intron, in addition to inversion in the native genomic background. Initial analysis revealed that A. thaliana has a significantly lower basal expression level of SCR transcripts in the critical reproductive stage compared to that of A. halleri, suggesting that the promoter was attenuated in inducing transcription in A. thaliana. However, in reciprocal transgenic experiments, this A. thaliana promoter was able to restore partial function if coupled with the functional A. halleri coding sequence, despite extensive alterations due to the self-compatible mode of reproduction in A. thaliana. This represents a synergistic effect of the promoter and the inversion resulting in fixation of self-compatibility, primarily enforced by disruption of SCR. Our findings elucidate the functional and evolutionary context of the historical transition in A. thaliana thus contributing to the understanding of the molecular events leading to development of self-compatibility.

Xenobiotic Effects of Chlorine Dioxide to Escherichia coli O157:H7 on Non-host Tomato Environment Revealed by Transcriptional Network Modeling: Implications to Adaptation and Selection.

Escherichia coli serotype O157:H7 is one of the major agents of pathogen outbreaks associated with fresh fruits and vegetables. Gaseous chlorine dioxide (ClO2) has been reported to be an effective intervention to eliminate bacterial contamination on fresh produce. Although remarkable positive effects of low doses of ClO2 have been reported, the genetic regulatory machinery coordinating the mechanisms of xenobiotic effects and the potential bacterial adaptation remained unclear. This study examined the temporal transcriptome profiles of E. coli O157:H7 during exposure to different doses of ClO2 in order to elucidate the genetic mechanisms underlying bacterial survival under such harsh conditions. Dosages of 1 µg, 5 µg, and 10 µg ClO2 per gram of tomato fruits cause different effects with dose-by-time dynamics. The first hour of exposure to 1 µg and 5 µg ClO2 caused only partial killing with significant growth reduction starting at the second hour, and without further significant reduction at the third hour. However, 10 µg ClO2 exposure led to massive bacterial cell death at 1 h with further increase in cell death at 2 and 3 h. The first hour exposure to 1 µg ClO2 caused activation of primary defense and survival mechanisms. However, the defense response was attenuated during the second and third hours. Upon treatment with 5 µg ClO2, the transcriptional networks showed massive downregulation of pathogenesis and stress response genes at the first hour of exposure, with decreasing number of differentially expressed genes at the second and third hours. In contrast, more genes were further downregulated with exposure to 10 µg ClO2 at the first hour, with the number of both upregulated and downregulated genes significantly decreasing at the second hour. A total of 810 genes were uniquely upregulated at the third hour at 10 µg ClO2, suggesting that the potency of xenobiotic effects had led to potential adaptation. This study provides important knowledge on the possible selection of target molecules for eliminating bacterial contamination on fresh produce without overlooking potential risks of adaptation.

Potential of Oryza officinalis to augment the cold tolerance genetic mechanisms of Oryza sativa by network complementation.

Oryza officinalis is an accessible alien donor for genetic improvement of rice. Comparison across a representative panel of Oryza species showed that the wild O. officinalis and cultivated O. sativa ssp. japonica have similar cold tolerance potentials. The possibility that either distinct or similar genetic mechanisms are involved in the low temperature responses of each species was addressed by comparing their transcriptional networks. General similarities were supported by shared transcriptomic signatures indicative of equivalent metabolic, hormonal, and defense status. However, O. officinalis has maintained an elaborate cold-responsive brassinosteroid-regulated BES1-network that appeared to have been fragmented in O. sativa. BES1-network is potentially important for integrating growth-related responses with physiological adjustments and defenses through the protection of photosynthetic machinery and maintenance of stomatal aperture, oxidative defenses, and osmotic adjustment. Equivalent physiological processes are functional in O. sativa but their genetic mechanisms are under the direct control of ABA-dependent, DREB-dependent and/or oxidative-mediated networks uncoupled to BES1. While O. officinalis and O. sativa represent long periods of speciation and domestication, their comparable cold tolerance potentials involve equivalent physiological processes but distinct genetic networks. BES1-network represents a novel attribute of O. officinalis with potential applications in diversifying or complementing other mechanisms in the cultivated germplasm.

Heat-induced mortality and expression of heat shock proteins in Colorado potato beetles treated with imidacloprid.

The Colorado potato beetle is an important pest of solanaceous plants in the Northern Hemisphere. Better understanding of its physiological responses to temperature stress and their interactions with still-prevalent chemical control has important implications for the management of this insect. We measured mortality and expression of the Hsp70 heat shock proteins in the Colorado potato beetle larvae exposed to sublethal concentration of the commonly used insecticide imidacloprid, and to supraoptimal temperatures. Both turned out to be significant stress factors, although induction of Hsp70 by imidacloprid observed in the present study was low compared to its induction by the heat. The two factors also interacted with each other. At an extreme temperature of 43 °C, exposure to a sublethal dose of imidacloprid resulted in a significant rise in larval mortality, which was not observed at an optimal temperature of 25 °C. Heat-stressed larvae also failed to respond to imidacloprid by producing more Hsp70. These findings suggest that when field rates of insecticides become insufficient for killing the exposed beetles under optimal temperature conditions due to the evolution of resistance in beetle populations, they may still reduce the probability of resistant beetles surviving the heat shock created by using propane flamers as a rescue treatment.

Identification of candidate network hubs involved in metabolic adjustments of rice under drought stress by integrating transcriptome data and genome-scale metabolic network.

In this study, we have integrated a rice genome-scale metabolic network and the transcriptome of a drought-tolerant rice line, DK151, to identify the major transcriptional regulators involved in metabolic adjustments necessary for adaptation to drought. This was achieved by examining the differential expressions of transcription factors and metabolic genes in leaf, root and young panicle of rice plants subjected to drought stress during tillering, booting and panicle elongation stages. Critical transcription factors such as AP2/ERF, bZIP, MYB and NAC that control the important nodes in the gene regulatory pathway were identified through correlative analysis of the patterns of spatio-temporal expression and cis-element enrichment. We showed that many of the candidate transcription factors involved in metabolic adjustments were previously linked to phenotypic variation for drought tolerance. This approach represents the first attempt to integrate models of transcriptional regulation and metabolic pathways for the identification of candidate regulatory genes for targeted selection in rice breeding.