Remarkably High Internal Transcribed Spacer Haplotype Diversity of the Fungal Select Agent Coniothyrium glycines Discovered Throughout Its Range in Sub-Saharan Africa

Red leaf blotch of soybean, caused by the fungus Coniothyrium glycines, is a foliar disease characterized by blotching, necrosis, and defoliation that has only been reported from Africa. The species is listed as a Select Agent by the Federal Select Agent Program due to its potentially devastating impacts to soybean production should it spread to the United States. Despite its potential import, very few isolates are available for study. Herein, we obtained 96 new C. glycines isolates from six soybean-producing countries throughout sub-Saharan Africa. Along with 12 previously collected ones, we sequenced each at the internal transcribed spacer (ITS) region. Between all isolates, we identified a total of 28 single-nucleotide polymorphisms and 23 haplotypes. One hypothesis to explain the tremendous diversity uncovered at the ITS-which is generally conserved within a species-is that our current species concept of C. glycines is too broad and that there may be multiple species that cause red leaf blotch. Zambia contained the highest haplotype diversity, a significant fraction of which remains unsampled. Most haplotypes were specific to a single country, except for two, which were found in Zambia and either neighboring Mozambique or Zimbabwe. This geographic specificity indicates that the ITS region may be useful for identifying source populations or routes of transmission should this pathogen spread beyond Africa. The observed geographic partitioning of this pathogen is likely the result of millions of years of replication on little-studied native hosts, given that soybean has only been cultivated in Africa since the early 1900s.

Tracking flight activity of potato leafhopper (Hemiptera: Cicadellidae) with the Midwest Suction Trap Network

Potato leafhopper (PLH), Empoasca fabae Harris (Hemiptera: Cicadellidae), is an economic pest of a variety of crops that migrates between overwintering sites in the southern United States and northern breeding grounds. Since 2005, the Midwest Suction Trap Network (STN) has monitored the magnitude and timing of aerially dispersing aphids' activity, but the potential of the network to monitor other taxa is only beginning to be explored. Here, we use the Midwest STN to examine how the magnitude and timing of PLH activity vary with weather, cropland cover, and time of year. We found that weekly PLH activity increased early in the season (May-June) with increasing degree day accumulation and decreased mid-season (July-August) with increasing occurrence of rain. The first detections occurred earlier in southern latitudes, while the last detections occurred sooner, when there was more surrounding potato land cover, and later over time between 2018 and 2021 and in southern latitudes. PLH activity was thus longer in duration in southern latitudes and has continued to extend later into the year overall. Resolving uncertainty about how well the Midwest STN captures migratory activity and how closely suction trap detections reflect local population densities in crop fields remain important research priorities before the potential of the Midwest STN for PLH monitoring can be realized. Still, observed patterns suggest that PLH could increase in economic importance as insects disperse over larger portions of the growing season in the warming, agriculturally productive US Midwest and that the STN can become a useful tool to monitor these changes.

The Role of a Soybean 14-3-3 Gene (Glyma05g29080) on White Mold Resistance and Nodulation Investigations Using CRISPR-Cas9 Editing and RNA Silencing

The role of a soybean 14-3-3 gene (Glyma05g29080) in defense against white mold and in nodulation was investigated by loss-of-gene-function with CRISPR-Cas9 editing and silencing of RNA interference (RNAi). Particle bombardment was used to introduce the CRISPR expression cassette to target the soybean 14-3-3 gene and an RNAi construct to silence gene transcription. Transmission of the edited 14-3-3 gene and the RNAi construct was confirmed in their respective progeny. The recovered transgenic plants and their progeny were significantly more susceptible to Sclerotinia sclerotiorum infection and showed a significant reduction in nodulation, thus confirming the role of the 14-3-3 gene (Glyma05g29080) in both nodulation and defense.

An NBS-LRR protein in the Rpp1 locus negates the dominance of Rpp1-mediated resistance against Phakopsora pachyrhizi in soybean

The soybean Rpp1 locus confers resistance to Phakopsora pachyrhizi, causal agent of rust, and resistance is usually dominant over susceptibility. However, dominance of Rpp1-mediated resistance is lost when a resistant genotype (Rpp1 or Rpp1b) is crossed with susceptible line TMG06_0011, and the mechanism of this dominant susceptibility (DS) is unknown. Sequencing the Rpp1 region reveals that the TMG06_0011 Rpp1 locus has a single nucleotide-binding site leucine-rich repeat (NBS-LRR) gene (DS-R), whereas resistant PI 594760B (Rpp1b) is similar to PI 200492 (Rpp1) and has three NBS-LRR resistance gene candidates. Evidence that DS-R is the cause of DS was reflected in virus-induced gene silencing of DS-R in Rpp1b/DS-R or Rpp1/DS-R heterozygous plants with resistance partially restored. In heterozygous Rpp1b/DS-R plants, expression of Rpp1b candidate genes was not significantly altered, indicating no effect of DS-R on transcription. Physical interaction of the DS-R protein with candidate Rpp1b resistance proteins was supported by yeast two-hybrid studies and in silico modeling. Thus, we conclude that suppression of resistance most likely does not occur at the transcript level, but instead probably at the protein level, possibly with Rpp1 function inhibited by binding to the DS-R protein. The DS-R gene was found in other soybean lines, with an estimated allele frequency of 6% in a diverse population, and also found in wild soybean (Glycine soja). The identification of a dominant susceptible NBS-LRR gene provides insight into the behavior of NBS-LRR proteins and serves as a reminder to breeders that the dominance of an R gene can be influenced by a susceptibility allele.

Characterization and analysis of the promoter region of monodehydroascorbate reductase 4 (CpMDAR4) in papaya

Differential spatial and temporal expression patterns due to regulatory cis-elements and two different isoforms are detected among CpMDAR4 alleles in papaya. The aim of this research was to study the effects of cis-element differences between the X, Y and Y^h alleles on the expression of CpMDAR4, a potential candidate gene for sex differentiation in papaya, using a transcriptional reporter system in a model species Arabidopsis thaliana. Possible effects of a retrotransposon insertion in the Y and Y^h alleles on the transcription and expression of CpMDAR4 alleles in papaya flowers were also examined. When comparing promoters and cis-regulatory elements among genes in the non-recombining region of the sex chromosomes, paired genes exhibited differences. Our results showed that differences in the promoter sequences of the CpMDAR4 alleles drove the expression of a reporter gene to different flower tissues in Arabidopsis. β-glucuronidase staining analysis of T₂ and T₃ lines for constructs containing 5' deletions of native Y and Y^h allele promoters showed the loss of specific expression of the reporter gene in the anthers, confirming the existence and location of cis-regulatory element POLLEN1LELAT52. The expression analysis of CpMDAR4 alleles in papaya flowers also showed that all alleles are actively expressed in different flower tissues, with the existence of a shorter truncated isoform, with unknown function, for the Y and Y^h alleles due to an LTR-RT insertion in the Y and Y^h chromosomes. The observed expression patterns in Arabidopsis thaliana flowers and the expression patterns of CpMDAR4 alleles in papaya flowers suggest that MDAR4 might have a role on development of reproductive organs in papaya, and that it constitutes an important candidate for sex differentiation.

Transcriptome Profiling Reveals Molecular Players in Early Soybean-Sclerotinia sclerotiorum Interaction

Sclerotinia sclerotiorum causes Sclerotinia stem rot on soybean. Using RNA sequencing, the transcriptomes of the soybean host and the S. sclerotiorum pathogen were simultaneously determined at 4 and 8 h postinoculation (hpi). Two soybean genotypes were involved: a resistant oxalate oxidase (OxO)-transgenic line and its susceptible parent, AC Colibri (AC). Of the 594 genes that were significantly induced by S. sclerotiorum, both hosts expressed genes related to jasmonic acid, ethylene, oxidative burst, and phenylpropanoids. In all, 36% of the differentially expressed genes encoded genes associated with transcription factors, ubiquitination, or general signaling transduction such as receptor-like kinases, mitogen-activated protein kinase kinases, and hormones. No significant differentially expressed genes were identified between genotypes, suggesting that oxalic acid (OA) did not play a differential role in early disease development or primary lesion formation under the conditions used. Looking at pathogen behavior through its gene expression during infection, thousands of genes in S. sclerotiorum were induced at 8 hpi, compared with expression in culture. Many plant cell-wall-degrading enzymes (PCWDEs), sugar transport genes, and genes involved in secondary metabolism were upregulated and could contribute to early pathogenesis. When infecting the OxO plants, there was a higher induction of genes encoding OA, botcinic acid, PCWDEs, proteases, and potential effectors, revealing the wealth of virulence factors available to this pathogen as it attempts to colonize a host. Data presented identify hundreds of genes associated with the very early stages of infection for both the host and pathogen.

Characterization and Comparison of Intestinal Bacterial Microbiomes of Euschistus heros and Piezodorus guildinii Collected in Brazil and the United States

Background: Herbivorous insects are one of the main biological threats to crops. One such group of insects, stink bugs, do not eat large amounts of tissue when feeding on soybean, but are damaging to the quality of the seed yield as they feed on green developing seeds leading to poorly marketable harvests. In addition to causing physical damage during sucking-feeding activities, the insects can also transmit microbial pathogens, leading to even greater yield loss. Conducting surveys of the insect intestinal microbiome can help identify possible pathogens, as well as detail what healthy stink bug digestive systems have in common. Methods: We used the conserved V4 region of the 16S rRNA gene to characterize the bacterial microbiome of the red-banded stink bug Piezodorus guildinii collected in Brazil and the United States, as well as the neotropical brown stink bug Euschistus heros collected in Brazil. Results: After quality filtering of the data, 192 samples were kept for analyses: 117 samples from P. guildinii covering three sites in Brazil and four sites in the United States, and 75 samples for E. heros covering 10 sites in Brazil. The most interesting observations were that the diversity and abundance of some bacterial families were different in the different ecoregions of Brazil and the United States. Conclusion: Some families, such as Acetobacteraceae, Bacillaceae, Moraxellaceae, Enterobacteriaceae, and Rhodocyclaceae, may be related to the better adaptation in some localities in providing nutrients, break down cellulose, detoxify phytochemicals, and degrade organic compounds, which makes it difficult to control these species.

A first draft genome of the Sugarcane borer, Diatraea saccharalis

Background: The sugarcane borer (Diatraea saccharalis), a widely distributed moth throughout the Americas, is a pest that affects economically important crops such as sugarcane, sorghum, wheat, maize and rice. Given its significant impact on yield reduction, whole-genome information of the species is needed. Here, we report the first draft assembly of the D. saccharalis genome. Methods: The genomic sequences were obtained using the Illumina HiSeq 2500 whole-genome sequencing of a single adult male specimen. We assembled the short-reads using the SPAdes software and predicted protein-coding genes using MAKER. Genome assembly completeness was assessed through BUSCO and the repetitive content by RepeatMasker. Results: The 453 Mb assembled sequences contain 1,445 BUSCO gene orthologs and 1,161 predicted gene models identified based on homology evidence to the domestic silk moth, Bombyx mori. The repeat content composes 41.18% of the genomic sequences which is in the range of other lepidopteran species. Conclusions: Functional annotation reveals that predicted gene models are involved in important cellular mechanisms such as metabolic pathways and protein synthesis. Thus, the data generated in this study expands our knowledge on the genomic characteristics of this devastating pest and provides essential resources for future genetic studies of the species.

Genome-wide expression of low temperature response genes in Rosa hybrida L

Plants respond to low temperature stress during cold acclimation, a complex process involving changes in physiological and biochemical modifications. The rose serves as a good model to investigate low temperature responses in perennial ornamentals. In this study, a heterologous apple microarray is used to investigate genome-wide expression profiles in Rosa hybrida subjected to low temperature dark treatment. Transcriptome profiles are determined in floral buds at 0h, 2h, and 12h of low temperature treatment (4 °C). It is observed that a total of 134 transcripts are up-regulated and 169 transcripts are down-regulated in response to low temperature. Interestingly, a total of eight up-regulated genes, including those coding for two cytochrome P450 proteins, two ankyrin repeat family proteins, two metal ion binding proteins, and two zinc finger protein-related transcription factors, along with a single down-regulated gene, coding for a dynamin-like protein, are detected. Transcript profiles of 12 genes known to be involved in cold stress response are also validated using qRT-PCR. Furthermore, expression patterns of the AP2/ERF gene family of transcription factors are investigated in both floral buds and leaves. Overall, AP2/ERFs genes are more rapidly induced in leaves than in floral buds. Moreover, differential expression of several AP2/ERF genes are detected earlier in vegetative rather than in reproductive tissues. These findings highlight important roles of various low temperature response genes in mediating cold acclimation, thereby allowing roses to adapt to low temperatures, but without adversely affecting flower bud development and subsequent flowering, while vegetative tissues undergo early adaptation to low temperatures.

Population Genomics of the Neotropical Brown Stink Bug, Euschistus heros: The Most Important Emerging Insect Pest to Soybean in Brazil

Recent changes in soybean management like the adoption of transgenic crops and no-till farming, in addition to the expansion of cultivated areas into new virgin frontiers, are some of the hypotheses that can explain the rise of secondary pests, such as the Neotropical brown stink bug, Euschistus heros, in Brazil. To better access the risk of increased pests like E. heros and to determine probabilities for insecticide resistance spreading, it is necessary first to access the levels of the genetic diversity, how the genetic diversity is distributed, and how natural selection is acting upon the natural variation. Using the genotyping by sequencing (GBS) technique, we generated ~60,000 single-nucleotide polymorphisms (SNPs) distributed across the E. heros genome to answer some of those questions. The SNP data was used to investigate the pattern of genetic structure, hybridization and natural selection of this emerging pest. We found that E. heros populations presented similar levels of genetic diversity with slightly higher values at several central locations in Brazil. Our results also showed strong genetic structure separating northern and southern Brazilian regions (F_ST = 0.22; p-value = 0.000) with a very distinct hybrid zone at the central region. The analyses also suggest the possibility that GABA channels and odorant receptors might play a role in the process of natural selection. At least one marker was associated with soybean and beans crops, but no association between allele frequency and cotton was found. We discuss the implications of these findings in the management of emerging pests in agriculture, particularly in the context of large areas of monoculture such as soybean and cotton.

Patterns of Genome-Wide Variation, Population Differentiation and SNP Discovery of the Red Banded Stink Bug (Piezodorus guildinii)

Unravelling the details of range expansion and ecological dominance shifts of insect pests has been challenging due to the lack of basic knowledge about population structure, gene flow, and most importantly, how natural selection is affecting the adaptive process. Piezodous guildinii is an emerging pest of soybean in the southern region of the United States, and increasingly important in Brazil in recent years. However, the reasons P. guildinii is gradually becoming more of a problem are questions still mostly unanswered. Here, we have genotyped P. guildinii samples and discovered 1,337 loci containing 4,083 variant sites SNPs that were used to estimate genetic structure and to identify gene candidates under natural selection. Our results revealed the existence of a significant genetic structure separating populations according to their broad geographic origin, i.e., U.S. and Brazil, supported by AMOVA (FGT = 0.26), STRUCTURE, PCA, and FST analyses. High levels of gene flow or coancestry within groups (i.e., within countries) can be inferred from the data, and no spatial pattern was apparent at the finer scale in Brazil. Samples from different seasons show more heterogeneous compositions suggesting mixed ancestry and a more complex dynamic. Lastly, we were able to detect and successfully annotated 123 GBS loci (10.5%) under positive selection. The gene ontology (GO) analysis implicated candidate genes under selection with genome reorganization, neuropeptides, and energy mobilization. We discuss how these findings could be related to recent outbreaks and suggest how new efforts directed to better understand P. guildinii population dynamics.

Diatraea saccharalis history of colonization in the Americas. The case for human-mediated dispersal

The sugarcane borer moth, Diatraea saccharalis, is one of the most important pests of sugarcane and maize crops in the Western Hemisphere. The pest is widespread throughout South and Central America, the Caribbean region and the southern United States. One of the most intriguing features of D. saccharalis population dynamics is the high rate of range expansion reported in recent years. To shed light on the history of colonization of D. saccharalis, we investigated the genetic structure and diversity in American populations using single nucleotide polymorphism (SNPs) markers throughout the genome and sequences of the mitochondrial gene cytochrome oxidase (COI). Our primary goal was to propose possible dispersal routes from the putative center of origin that can explain the spatial pattern of genetic diversity. Our findings showed a clear correspondence between genetic structure and the geographical distributions of this pest insect on the American continents. The clustering analyses indicated three distinct groups: one composed of Brazilian populations, a second group composed of populations from El Salvador, Mexico, Texas and Louisiana and a third group composed of the Florida population. The predicted time of divergence predates the agriculture expansion period, but the pattern of distribution of haplotype diversity suggests that human-mediated movement was most likely the factor responsible for the widespread distribution in the Americas. The study of the early history of D. saccharalis promotes a better understanding of range expansion, the history of invasion, and demographic patterns of pest populations in the Americas.

Morphological and molecular characterization of Brazilian populations of Diatraea saccharalis (Fabricius, 1794) (Lepidoptera: Crambidae) and the evolutionary relationship among species of Diatraea Guilding

The sugarcane borer or corn stalk borer, Diatraea Guilding is polyphagous insect pest of many important crops such as corn, sorghum and sugarcane. Losses arising from the attack of Diatraea species have been a serious problem, which may cause loss in sugarcane production around 0.25% in sugar, 0.20% in alcohol and 0.77% of body weight for every 1% infestation and up to 21% in corn production fields. In Brazil, the most commonly reported species are Diatraea saccharalis (Fabricius, 1794) and Diatraea impersonatella (Walker, 1863) (= D. flavipennella). However, multiple other species of Diatraea have been identified in Brazil according to the literature. Currently, little information exists on the presence of the other species causing injury to sugarcane and corn. The objectives of this study were to improve the accuracy of species assignment, evaluate the population genetic structure, and address many of the outstanding questions of systematics and evolution of Brazilian populations of D. saccharalis. To address these main questions, classical taxonomic methods were used, focused on morphological characterization of the reproductive organs, especially the male genitalia. In addition, genetic studies were performed using simple sequence repeats (SSR) and a fragment of cytochrome C oxidase subunit I (COI) gene. The data and findings from this research will contribute to the understanding of evolutionary aspects of insect pests in order to develop more effective and sustainable population management practices.

Genome-wide association mapping of resistance to a Brazilian isolate of Sclerotinia sclerotiorum in soybean genotypes mostly from Brazil

BACKGROUND

Sclerotinia Stem Rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is ubiquitous in cooler climates where soybean crops are grown. Breeding for resistance to SSR remains challenging in crops like soybean, where no single gene provides strong resistance, but instead, multiple genes work together to provide partial resistance. In this study, a genome-wide association study (GWAS) was performed to dissect the complex genetic architecture of soybean quantitative resistance to SSR and to provide effective molecular markers that could be used in breeding programs. A collection of 420 soybean genotypes were selected based on either reports of resistance, or from one of three different breeding programs in Brazil, two commercial, one public. Plant genotype sensitivity to SSR was evaluated by the cut stem inoculation method, and lesion lengths were measured at 4 days post inoculation.

RESULTS

Genotyping-by-sequencing was conducted to genotype the 420 soybean lines. The TASSEL 5 GBSv2 pipeline was used to call SNPs under optimized parameters, and with the extra step of trimming adapter sequences. After filtering missing data, heterozygosity, and minor allele frequency, a total of 11,811 SNPs and 275 soybean genotypes were obtained for association analyses. Using a threshold of FDR-adjusted p-values <0.1, the Compressed Mixed Linear Model (CMLM) with Genome Association and Prediction Integrated Tool (GAPIT), and the Fixed and Random Model Circulating Probability Unification (FarmCPU) methods, both approaches identified SNPs with significant association to disease response on chromosomes 1, 11, and 18. The CMLM also found significance on chromosome 19, whereas FarmCPU also identified significance on chromosomes 4, 9, and 16.

CONCLUSIONS

These similar and yet different results show that the computational methods used can impact SNP associations in soybean, a plant with a high degree of linkage disequilibrium, and in SSR resistance, a trait that has a complex genetic basis. A total of 125 genes were located within linkage disequilibrium of the three loci shared between the two models. Their annotations and gene expressions in previous studies of soybean infected with S. sclerotiorum were examined to narrow down the candidates.

Resistance of soybean genotypes to Sclerotinia sclerotiorum isolates in different incubation environments

Sclerotinia sclerotiorum is an important soybean pathogen. The objectives of this study were to evaluate levels of resistance of soybean genotypes to the fungus, and to determine the effects of different incubation environments on host resistance and pathogen aggressiveness. Two experiments were conducted using 103 genotypes from the seed collection of Laboratório de Desenvolvimento de Germoplasma, Universidade Federal de Uberlândia (LAGER-UFU). The first experiment was conducted in a greenhouse, and the second in a growth chamber. Inoculations were performed by the straw test method using two Brazilian isolates of the fungus, one from Uberaba, Minas Gerais, and the other from Jataí, Goiás. The average stem-lesion length (cm) at 5 days post-inoculation was used to determine levels of resistance. Overall, the most resistant genotype was EMGOPA-316, and the most susceptible genotype was LAGER-29. Incubation in a growth chamber and use of the Jataí isolate generated the most reliable data, and multivariate analysis indicated that the genotypes were divergent under the growth chamber conditions. Therefore, when studying host resistance of soybean genotypes to S. sclerotiorum, it is important to use environmental conditions favorable to the fungus and aggressive isolates.

Investigating the Control of Chlorophyll Degradation by Genomic Correlation Mining

Chlorophyll degradation is an intricate process that is critical in a variety of plant tissues at different times during the plant life cycle. Many of the photoactive chlorophyll degradation intermediates are exceptionally cytotoxic necessitating that the pathway be carefully coordinated and regulated. The primary regulatory step in the chlorophyll degradation pathway involves the enzyme pheophorbide a oxygenase (PAO), which oxidizes the chlorophyll intermediate pheophorbide a, that is eventually converted to non-fluorescent chlorophyll catabolites. There is evidence that PAO is differentially regulated across different environmental and developmental conditions with both transcriptional and post-transcriptional components, but the involved regulatory elements are uncertain or unknown. We hypothesized that transcription factors modulate PAO expression across different environmental conditions, such as cold and drought, as well as during developmental transitions to leaf senescence and maturation of green seeds. To test these hypotheses, several sets of Arabidopsis genomic and bioinformatic experiments were investigated and re-analyzed using computational approaches. PAO expression was compared across varied environmental conditions in the three separate datasets using regression modeling and correlation mining to identify gene elements co-expressed with PAO. Their functions were investigated as candidate upstream transcription factors or other regulatory elements that may regulate PAO expression. PAO transcript expression was found to be significantly up-regulated in warm conditions, during leaf senescence, and in drought conditions, and in all three conditions significantly positively correlated with expression of transcription factor Arabidopsis thaliana activating factor 1 (ATAF1), suggesting that ATAF1 is triggered in the plant response to these processes or abiotic stresses and in result up-regulates PAO expression. The proposed regulatory network includes the freezing, senescence, and drought stresses modulating factor ATAF1 and various other transcription factors and pathways, which in turn act to regulate chlorophyll degradation by up-regulating PAO expression.

The role of invertases in plant compensatory responses to simulated herbivory

BACKGROUND

The ability of a plant to overcome animal-induced damage is referred to as compensation or tolerance and ranges from undercompensation (decreased fitness when damaged) to overcompensation (increased fitness when damaged). Although it is clear that genetic variation for compensation exists among plants, little is known about the specific genetic underpinnings leading to enhanced fitness. Our previous study identified the enzyme GLUCOSE-6-PHOSPHATE DEHYDROGENASE 1 (G6PD1) as a key regulator contributing to the phenomenon of overcompensation via its role in the oxidative pentose phosphate pathway (OPPP). Apart from G6PD1 we also identified an invertase gene which was up-regulated following damage and that potentially integrates with the OPPP. The invertase family of enzymes hydrolyze sucrose to glucose and fructose, whereby the glucose produced is shunted into the OPPP and presumably supports plant regrowth, development, and ultimately compensation. In the current study, we measured the relative expression of 12 invertase genes over the course of plant development in the Arabidopsis thaliana genotypes Columbia-4 and Landsberg erecta, which typically overcompensate and undercompensate, respectively, when damaged. We also compared the compensatory performances of a set of invertase knockout mutants to the Columbia-4 wild type.

RESULTS

We report that Columbia-4 significantly up-regulated 9 of 12 invertase genes when damaged relative to when undamaged, and ultimately overcompensated for fruit production. Landsberg erecta, in contrast, down-regulated two invertase genes following damage and suffered reduced fitness. Knockout mutants of two invertase genes both exhibited significant undercompensation for fruit production, exhibiting a complete reversal of the wild type Col-4's overcompensation.

CONCLUSION

Collectively, these results confirm that invertases are essential for not only normal plant growth and development, but also plants' abilities to regrow and ultimately compensate for fitness following apical damage.

Transcriptome analyses suggest a disturbance of iron homeostasis in soybean leaves during white mould disease establishment

Sclerotinia sclerotiorum is a serious pathogen of numerous crops around the world. The major virulence factor of this pathogen is oxalic acid (OA). Mutants that cannot produce OA do not cause disease, and plants that express enzymes that degrade OA, such as oxalate oxidase (OxO), are very resistant to S. sclerotiorum. To examine the effect of OA on plants, we infiltrated soybean leaves with 5 mm OA and examined the gene expression changes at 2 h post-infiltration. By comparing the gene expression levels between leaves of a transgenic soybean carrying an OxO gene (OxO) and its parent AC Colibri (AC) infiltrated with OA (pH 2.4) or water (pH 2.4 or 5.5), we were able to compare the effects of OA dependent or independent of its pH. Gene expression by microarray analysis identified 2390 genes that showed changes in expression, as determined using an overall F-test P-value cut-off of 0.001. The additional requirement that at least one pairwise t-test false discovery rate (FDR)-corrected P value should be less than 0.001 reduced the list of the most highly significant differentially expressed genes to 1054. Independent of pH, OA altered the expression levels of 78 genes, with ferritin showing the strongest induction by OA. The combination of OA plus its low pH caused 1045 genes (99% of all significant genes) to be differentially expressed, with many of the up-regulated genes being related to basal defence, such as genes of the phenylpropanoid pathway and various cytochrome P450s. RNA-seq was also conducted on four samples: OxO and AC genotypes infiltrated with either OA pH 2.4 or water pH 2.4. The RNA-seq analysis also identified ferritin paralogues as being strongly induced by OA. As the expression of ferritin, a gene that encodes for an iron storage protein, is induced by free iron, these results suggest that S. sclerotiorum benefits from the ability of OA to free iron from plant proteins, as this induces host cell death, and also allows the uptake and assimilation of the iron for its own metabolic needs.

Genomic evaluation of oxalate-degrading transgenic soybean in response to Sclerotinia sclerotiorum infection

Oxalate oxidases (OxO) catalyse the degradation of oxalic acid (OA). Highly resistant transgenic soybean carrying an OxO gene and its susceptible parent soybean line, AC Colibri, were tested for genome-wide gene expression in response to the necrotrophic, OA-producing pathogen Sclerotinia sclerotiorum using soybean cDNA microarrays. The genes with changed expression at statistically significant levels (overall F-test P-value cut-off of 0.0001) were classified into functional categories and pathways, and were analysed to evaluate the differences in transcriptome profiles. Although many genes and pathways were found to be similarly activated or repressed in both genotypes after inoculation with S. sclerotiorum, the OxO genotype displayed a measurably faster induction of basal defence responses, as observed by the differential changes in defence-related and secondary metabolite genes compared with its susceptible parent AC Colibri. In addition, the experiment presented provides data on several other transcripts that support the hypothesis that S. sclerotiorum at least partially elicits the hypersensitive response, induces lignin synthesis (cinnamoyl CoA reductase) and elicits as yet unstudied signalling pathways (G-protein-coupled receptor and related). Of the nine genes showing the most extreme opposite directions of expression between genotypes, eight were related to photosynthesis and/or oxidation, highlighting the importance of redox in the control of this pathogen.

Effect of Fusarium virguliforme phytotoxin on soybean gene expression suggests a role in multidimensional defence

Sudden death syndrome (SDS), caused by Fusarium virguliforme, is an important yield-limiting disease of soybean. This soil-borne fungus colonizes soybean roots causing root rot, and also releases a phytotoxin that is translocated to leaves causing interveinal chlorosis and necrosis leading to defoliation and early maturation. The objective of our study was to compare gene expression profiles during the early response of soybean leaves exposed to sterile culture filtrates of F. virguliforme in soybean genotypes with different levels of resistance to SDS. The analysis identified SDS-related defence genes that were induced in the most resistant genotype, but not in the other genotypes. Further functional annotations based on sequence homology suggested that some of the induced genes probably encode proteins involved in cell wall modification, detoxification, defence responses, primary metabolism and membrane transport. Quantitative real-time reverse-transcribed polymerase chain reaction confirmed the differential transcript accumulation of a subset of these genes. In addition, in silico mapping of differentially expressed genes to SDS-resistant quantitative trait loci allowed for the identification of new potential defence genes that could be genetically mapped to the soybean genome, and could be used further in a marker-assisted selection programme. A comparison of the response of soybean to F. virguliforme phytotoxin (Fv toxin) relative to other biotic and abiotic stresses revealed that the resistance response to Fv toxin is quite similar to the response to inoculation with an incompatible Pseudomonas syringae pv. glycinea strain, suggesting that Fv toxin might induce hypersensitive response pathways in soybean leaf tissues in the absence of pathogen in these tissues.

14-3-3 proteins SGF14c and SGF14l play critical roles during soybean nodulation

The soybean (Glycine max) genome contains 18 members of the 14-3-3 protein family, but little is known about their association with specific phenotypes. Here, we report that the Glyma0529080 Soybean G-box Factor 14-3-3c (SGF14c) and Glyma08g12220 (SGF14l) genes, encoding 14-3-3 proteins, appear to play essential roles in soybean nodulation. Quantitative reverse transcription-polymerase chain reaction and western-immunoblot analyses showed that SGF14c mRNA and protein levels were specifically increased in abundance in nodulated soybean roots 10, 12, 16, and 20 d after inoculation with Bradyrhizobium japonicum. To investigate the role of SGF14c during soybean nodulation, RNA interference was employed to silence SGF14c expression in soybean roots using Agrobacterium rhizogenes-mediated root transformation. Due to the paleopolyploid nature of soybean, designing a specific RNA interference sequence that exclusively targeted SGF14c was not possible. Therefore, two highly similar paralogs (SGF14c and SGF14l) that have been shown to function as dimers were silenced. Transcriptomic and proteomic analyses showed that mRNA and protein levels were significantly reduced in the SGF14c/SGF14l-silenced roots, and these roots exhibited reduced numbers of mature nodules. In addition, SGF14c/SGF14l-silenced roots contained large numbers of arrested nodule primordia following B. japonicum inoculation. Transmission electron microscopy further revealed that the host cytoplasm and membranes, except the symbiosome membrane, were severely degraded in the failed nodules. Altogether, transcriptomic, proteomic, and cytological data suggest a critical role of one or both of these 14-3-3 proteins in early development stages of soybean nodules.

Divergent patterns of endogenous small RNA populations from seed and vegetative tissues of Glycine max

BACKGROUND

Small non-coding RNAs (smRNAs) are known to have major roles in gene regulation in eukaryotes. In plants, knowledge of the biogenesis and mechanisms of action of smRNA classes including microRNAs (miRNAs), short interfering RNAs (siRNAs), and trans-acting siRNAs (tasiRNAs) has been gained mostly through studies with Arabidopsis. In recent years, high throughput sequencing of smRNA populations has enabled extension of knowledge from model systems to plants with larger, more complex genomes. Soybean (Glycine max) now has many genomics resources available including a complete genome sequence and predicted gene models. Relatively little is known, however, about the full complement of its endogenous smRNAs populations and the silenced genes.

RESULTS

Using Illumina sequencing and computational analysis, we characterized eight smRNA populations from multiple tissues and organs of soybean including developing seed and vegetative tissues. A total of 41 million raw sequence reads collapsed into 135,055 unique reads were mapped to the soybean genome and its predicted cDNA gene models. Bioinformatic analyses were used to distinguish miRNAs and siRNAs and to determine their genomic origins and potential target genes. In addition, we identified two soybean TAS3 gene homologs, the miRNAs that putatively guide cleavage of their transcripts, and the derived tasiRNAs that could target soybean genes annotated as auxin response factors. Tissue-differential expression based on the flux of normalized miRNA and siRNA abundances in the eight smRNA libraries was evident, some of which was confirmed by smRNA blotting. Our global view of these smRNA populations also revealed that the size classes of smRNAs varied amongst different tissues, with the developing seed and seed coat having greater numbers of unique smRNAs of the 24-nt class compared to the vegetative tissues of germinating seedlings. The 24-nt class is known to be derived from repetitive elements including transposons. Detailed analysis of the size classes associated with ribosomal RNAs and transposable element families showed greater diversity of smRNAs in the 22- and 24-nt size classes.

CONCLUSIONS

The flux of endogenous smRNAs within multiple stages and tissues of seed development was contrasted with vegetative tissues of soybean, one of the dominant sources of protein and oil in world markets. The smRNAs varied in size class, complexity of origins, and possible targets. Sequencing revealed tissue-preferential expression for certain smRNAs and expression differences among closely related miRNA family members.

The orphan gene ybjN conveys pleiotropic effects on multicellular behavior and survival of Escherichia coli

YbjN, encoding an enterobacteria-specific protein, is a multicopy suppressor of temperature sensitivity in the ts9 mutant strain of Escherichia coli. In this study, we further explored the role(s) of ybjN. First, we demonstrated that the ybjN transcript was about 10-fold lower in the ts9 strain compared to that of E. coli strain BW25113 (BW). Introduction of multiple copies of ybjN in the ts9 strain resulted in over-expression of ybjN by about 10-fold as compared to that of BW. These results suggested that temperature sensitivity of the ts9 mutant of E. coli may be related to expression levels of ybjN. Characterization of E. coli ybjN mutant revealed that ybjN mutation resulted in pleiotropic phenotypes, including increased motility, fimbriation (auto-aggregation), exopolysaccharide production, and biofilm formation. In contrast, over-expression of ybjN (in terms of multiple copies) resulted in reduced motility, fimbriation, exopolysaccharide production, biofilm formation and acid resistance. In addition, our results indicate that a ybjN-homolog gene from Erwinia amylovora, a plant enterobacterial pathogen, is functionally conserved with that of E. coli, suggesting similar evolution of the YbjN family proteins in enterobacteria. A microarray study revealed that the expression level of ybjN was inversely correlated with the expression of flagellar, fimbrial and acid resistance genes. Over-expression of ybjN significantly down-regulated genes involved in citric acid cycle, glycolysis, the glyoxylate shunt, oxidative phosphorylation, amino acid and nucleotide metabolism. Furthermore, over-expression of ybjN up-regulated toxin-antitoxin modules, the SOS response pathway, cold shock and starvation induced transporter genes. Collectively, these results suggest that YbjN may play important roles in regulating bacterial multicellular behavior, metabolism, and survival under stress conditions in E. coli. These results also suggest that ybjN over-expression-related temperature rescue of the ts9 mutant may be due to down-regulation of metabolic activity and activation of stress response genes in the ts9 mutant.

Transcriptional analysis of soybean root response to Fusarium virguliforme, the causal agent of sudden death syndrome

Sudden death syndrome (SDS) of soybean can be caused by any of four distinct Fusarium species, with F. virguliforme and F. tucumaniae being the main casual agents in North and South America, respectively. Although the fungal tissue is largely confined to the roots, the fungus releases a toxin that is translocated to leaf tissues, in which it causes interveinal chlorosis and necrosis leading to scorching symptoms and possible defoliation. In this study, we report on an Affymetrix analysis measuring transcript abundances in resistant (PI 567.374) and susceptible (Essex) roots upon infection by F. virguliforme, 5 and 7 days postinoculation. Many of the genes with increased expression were common between resistant and susceptible plants (including genes related to programmed cell death, the phenylpropanoid pathway, defense, signal transduction, and transcription factors), but some genotype-specific expression was noted. Changes in small (sm)RNA levels between inoculated and mock-treated samples were also studied and implicate a role for these molecules in this interaction. In total, 2,467 genes were significantly changing in the experiment, with 1,694 changing in response to the pathogen; 93 smRNA and 42 microRNA that have putative soybean gene targets were identified from infected tissue. Comparing genotypes, 247 genes were uniquely modulating in the resistant host, whereas 378 genes were uniquely modulating in the susceptible host. Comparing locations of differentially expressed genes to known resistant quantitative trait loci as well as identifying smRNA that increased while their putative targets decreased (or vice versa) allowed for the narrowing of candidate SDS defense-associated genes.

Biotic stress globally downregulates photosynthesis genes

To determine if damage to foliage by biotic agents, including arthropods, fungi, bacteria and viral pathogens, universally downregulates the expression of genes involved in photosynthesis, we compared transcriptome data from microarray experiments after twenty two different forms of biotic damage on eight different plant species. Transcript levels of photosynthesis light reaction, carbon reduction cycle and pigment synthesis genes decreased regardless of the type of biotic attack. The corresponding upregulation of genes coding for the synthesis of jasmonic acid and those involved in the responses to salicylic acid and ethylene suggest that the downregulation of photosynthesis-related genes was part of a defence response. Analysis of the sub-cellular targeting of co-expressed gene clusters revealed that the transcript levels of 84% of the genes that carry a chloroplast targeting peptide sequence decreased. The majority of these downregulated genes shared common regulatory elements, such as G-box (CACGTG), T-box (ACTTTG) and SORLIP (GCCAC) motifs. Strong convergence in the response of transcription suggests that the universal downregulation of photosynthesis-related gene expression is an adaptive response to biotic attack. We hypothesize that slow turnover of many photosynthetic proteins allows plants to invest resources in immediate defence needs without debilitating near term losses in photosynthetic capacity.

Differential Expression of Virulence Genes and Motility in Ralstonia (Pseudomonas) solanacearum during Exponential Growth

A complex network regulates virulence in Ralstonia solanacearum (formerly Pseudomonas solanacearum); central to this system is PhcA, a LysR-type transcriptional regulator. We report here that two PhcA-regulated virulence factors, endoglucanase (Egl) and acidic exopolysaccharide I (EPS I), and motility are expressed differentially during exponential growth in batch cultures. Tests with strains carrying lacZ fusions in a wild-type genetic background revealed that expression (on a per-cell basis) of phcA was constant but expression of egl and epsB increased 20- to 50-fold during multiplication from 1 x 10(sup7) to 5 x 10(sup8) CFU/ml. Expression of xpsR, an intermediate regulator downstream of PhcA in the regulatory cascade for eps expression, was similar to that of epsB and egl. Motility track photography revealed that all strains were essentially nonmotile at 10(sup6) CFU/ml. As cell density increased, 30 to 50% of wild-type cells were motile between 10(sup7) and 10(sup8) CFU/ml, but this population was again nonmotile at 10(sup9) CFU/ml. In contrast, about 60% of the cells of phcB and phcA mutants remained motile at 10(sup9) CFU/ml. Expression of phcB, which is not positively regulated by PhcA, was the inverse of epsB, egl, and xpsR (i.e., it decreased 20-fold at high cell density). PhcB is essential for production of an extracellular factor, tentatively identified as 3-hydroxypalmitic acid methyl ester (3-OH PAME), that might act as an exponential-phase signal to activate motility or expression of virulence genes. However, growth of the lacZ fusion strains in medium containing excess 3-OH PAME did not result in motility or expression of virulence genes at dramatically lower cell densities, suggesting that 3-OH PAME is not the only factor controlling these traits.

Complete transcriptome of the soybean root hair cell, a single-cell model, and its alteration in response to Bradyrhizobium japonicum infection

Nodulation is the result of a mutualistic interaction between legumes and symbiotic soil bacteria (e.g. soybean [Glycine max] and Bradyrhizobium japonicum) initiated by the infection of plant root hair cells by the symbiont. Fewer than 20 plant genes involved in the nodulation process have been functionally characterized. Considering the complexity of the symbiosis, significantly more genes are likely involved. To identify genes involved in root hair cell infection, we performed a large-scale transcriptome analysis of B. japonicum-inoculated and mock-inoculated soybean root hairs using three different technologies: microarray hybridization, Illumina sequencing, and quantitative real-time reverse transcription-polymerase chain reaction. Together, a total of 1,973 soybean genes were differentially expressed with high significance during root hair infection, including orthologs of previously characterized root hair infection-related genes such as NFR5 and NIN. The regulation of 60 genes was confirmed by quantitative real-time reverse transcription-polymerase chain reaction. Our analysis also highlighted changes in the expression pattern of some homeologous and tandemly duplicated soybean genes, supporting their rapid specialization.

A robust plant RNA isolation method suitable for Affymetrix GeneChip analysis and quantitative real-time RT-PCR

Microarray analysis and quantitative real-time RT-PCR are the major high-throughput techniques that are used to study transcript profiles. One of the major limitations in these technologies is the isolation of large quantities of highly pure RNA from plant tissues rich in complex polysaccharides, polyphenolics and waxes. Any contamination of the isolated RNA affects the downstream applications and requires extra cleaning procedures that result in a reduced RNA yield, especially the low molecular weight molecules. The protocol presented here is suitable for isolating high yield and clean total RNA from field-grown plants. Unlike current methods, such as LiCl and TRIZOL, with this new method, the isolated RNA can be used directly for Affymetrix GeneChip labeling or real-time RT-PCR without further purification. This fast and simple protocol provides ready-to-use RNA within 4-5 h after sampling. Additionally, the protocol described here maintains the isolation of small RNA molecules, making it an ideal choice for plant RNA preparation prior to high-throughput sequencing methods to study gene expression.

Elevated ozone alters soybean-virus interaction

Increasing concentrations of ozone (O(3)) in the troposphere affect many organisms and their interactions with each other. To analyze the changes in a plant-pathogen interaction, soybean plants were infected with Soybean mosaic virus (SMV) while they were fumigated with O(3). In otherwise natural field conditions, elevated O(3) treatment slowed systemic infection and disease development by inducing a nonspecific resistance against SMV for a period of 3 weeks. During this period, the negative effect of virus infection on light-saturated carbon assimilation rate was prevented by elevated O(3) exposure. To identify the molecular basis of a soybean nonspecific defense response, high-throughput gene expression analysis was performed in a controlled environment. Transcripts of fungal, bacterial, and viral defense-related genes, including PR-1, PR-5, PR-10, and EDS1, as well as genes of the flavonoid biosynthesis pathways (and concentrations of their end products, quercetin and kaempherol derivatives) increased in response to elevated O(3). The drastic changes in soybean basal defense response under altered atmospheric conditions suggest that one of the elements of global change may alter the ecological consequences and, eventually, coevolutionary relationship of plant-pathogen interactions in the future.

Transcriptome response to glyphosate in sensitive and resistant soybean

The majority of soybeans planted in the United States are resistant to glyphosate due to introduction of a gene encoding for a glyphosate-insensitive 5-enolypyruvylshikimate-3-phosphate synthase. Gene expression profiling was conducted using cDNA microarrays to address questions related to potential secondary effects of glyphosate. When glyphosate-sensitive plants were treated with glyphosate, 3, 170, and 311 genes were identified as having different transcript levels at 1, 4, and 24 h post-treatment (hpt), respectively. Differentially expressed genes were classified into functional categories, and their possible roles in response to glyphosate are briefly discussed. Gene expression profiling of glyphosate-resistant plants treated with glyphosate indicated that the plants were marginally affected at 1 hpt and then quickly adjusted to glyphosate treatment. Ten, four, and four genes were identified as differentially expressed at 1, 4, and 24 hpt. When gene expression profiles of cotyledons from developing seed were compared between the near-isogenic resistant and sensitive lines, two genes were identified as significantly differentially expressed out of 27000, which was less than the empirical false-discovery rate determined from a control experiment. Quantitative real-time reverse-transcribed Polymerase Chain Reaction was conducted on selected genes and yielded results consistent with those from the microarrays. Collectively, these data indicate that there are no major transcriptomic changes associated with currently used glyphosate-resistant soybean.

Transcription profiling of soybean nodulation by Bradyrhizobium japonicum

Legumes interact with nodulating bacteria that convert atmospheric nitrogen into ammonia for plant use. This nitrogen fixation takes place within root nodules that form after infection of root hairs by compatible rhizobia. Using cDNA microarrays, we monitored gene expression in soybean (Glycine max) inoculated with the nodulating bacterium Bradyrhizobium japonicum 4, 8, and 16 days after inoculation, timepoints that coincide with nodule development and the onset of nitrogen fixation. This experiment identified several thousand genes that were differentially expressed in response to B. japonicum inoculation. Expression of 27 genes was analyzed by quantitative reverse transcriptase-polymerase chain reaction, and their expression patterns mimicked the microarray results, confirming integrity of analyses. The microarray results suggest that B. japonicum reduces plant defense responses during nodule development. In addition, the data revealed a high level of regulatory complexity (transcriptional, post-transcriptional, translational, post-translational) that is likely essential for development of the symbiosis and adjustment to an altered nutritional status.

Soybean defense responses to the soybean aphid

Transcript profiles in aphid (Aphis glycines)-resistant (cv. Dowling) and -susceptible (cv. Williams 82) soybean (Glycine max) cultivars using soybean cDNA microarrays were investigated. Large-scale soybean cDNA microarrays representing approx. 18 000 genes or c. 30% of the soybean genome were compared at 6 and 12 h post-application of aphids. In a separate experiment utilizing clip cages, expression of three defense-related genes were examined at 6, 12, 24, 48, and 72 h in both cultivars by quantitative real-time PCR. One hundred and forty genes showed specific responses for resistance; these included genes related to cell wall, defense, DNA/RNA, secondary metabolism, signaling and other processes. When an extended time period of sampling was investigated, earlier and greater induction of three defense-related genes was observed in the resistant cultivar; however, the induction declined after 24 or 48 h in the resistant cultivar but continued to increase in the susceptible cultivar after 24 h. Aphid-challenged resistant plants showed rapid differential gene expression patterns similar to the incompatible response induced by avirulent Pseudomonas syringae. Five genes were identified as differentially expressed between the two genotypes in the absence of aphids.

Transcriptome changes in the phenylpropanoid pathway of Glycine max in response to Pseudomonas syringae infection

BACKGROUND

Reports of plant molecular responses to pathogenic infections have pinpointed increases in activity of several genes of the phenylpropanoid pathway leading to the synthesis of lignin and flavonoids. The majority of those findings were derived from single gene studies and more recently from several global gene expression analyses. We undertook a global transcriptional analysis focused on the response of genes of the multiple branches of the phenylpropanoid pathway to infection by the Pseudomonas syringae pv. glycinea with or without the avirulence gene avrB to characterize more broadly the contribution of the multiple branches of the pathway to the resistance response in soybean. Transcript abundance in leaves was determined from analysis of soybean cDNA microarray data and hybridizations to RNA blots with specific gene probes.

RESULTS

The majority of the genes surveyed presented patterns of increased transcript accumulation. Some increased rapidly, 2 and 4 hours after inoculation, while others started to accumulate slowly by 8-12 hours. In contrast, transcripts of a few genes decreased in abundance 2 hours post inoculation. Most interestingly was the opposite temporal fluctuation in transcript abundance between early responsive genes in defense (CHS and IFS1) and F3H, the gene encoding a pivotal enzyme in the synthesis of anthocyanins, proanthocyanidins and flavonols. F3H transcripts decreased rapidly 2 hours post inoculation and increased during periods when CHS and IFS transcripts decreased. It was also determined that all but one (CHS4) family member genes (CHS1, CHS2, CHS3, CHS5, CHS6 and CHS7/8) accumulated higher transcript levels during the defense response provoked by the avirulent pathogen challenge.

CONCLUSION

Based on the mRNA profiles, these results show the strong bias that soybean has towards increasing the synthesis of isoflavonoid phytoalexins concomitant with the down regulation of genes required for the synthesis of anthocyanins and proanthocyanins. Although proanthocyanins are known to be toxic compounds, the cells in the soybean leaves seem to be programmed to prioritize the synthesis and accumulation of isoflavonoid and pterocarpan phytoalexins during the resistance response. It was known that CHS transcripts accumulate in great abundance rapidly after inoculation of the soybean plants but our results have demonstrated that all but one (CHS4) member of the gene family member genes accumulated higher transcript levels during the defense response.

Expression profiling soybean response to Pseudomonas syringae reveals new defense-related genes and rapid HR-specific downregulation of photosynthesis

Transcript profiling during susceptible (S) and hypersensitive response-associated resistance (R) interactions was determined in soybean (Glycine max). Pseudomonas syringae pv. glycinea carrying or lacking the avirulence gene avrB, was infiltrated into cultivar Williams 82. Leaf RNA was sampled at 2, 8, and 24 h postinoculation (hpi). Significant changes in transcript abundance were observed for 3,897 genes during the experiment at P < or = 0.000005. Many of the genes showed a similar direction of increase or decrease in abundance in both the S and R responses, but the R response generally showed a significantly greater degree of differential expression. More than 25% of these responsive genes had not been previously reported as being associated with pathogen interactions, as 704 had no functional annotation and 378 had no homolog in National Center for Biotechnology Information databases. The highest number of transcriptional changes was noted at 8 hpi, including the downregulation of 94 chloroplast-associated genes specific to the R response. Photosynthetic measurements were consistent with an R-specific reduction in photosystem II operating efficiency (phiPSII) that was apparent at 8 hpi for the R response with little effect in the S or control treatments. Imaging analyses suggest that the decreased phiPSII was a result of physical damage to PSII reaction centers.

Features of a 103-kb gene-rich region in soybean include an inverted perfect repeat cluster of CHS genes comprising the I locus

The I locus in soybean (Glycine max) corresponds to a region of chalcone synthase (CHS) gene duplications affecting seed pigmentation. We sequenced and annotated BAC clone 104J7, which harbors a dominant i(i) allele from Glycine max 'Williams 82', to gain insight into the genetic structure of this multigenic region in addition to examining its flanking regions. The 103-kb BAC encompasses a gene-rich region with 11 putatively expressed genes. In addition to six copies of CHS, these genes include: a geranylgeranyltransferase type II beta subunit (E.C.2.5.1.60), a beta-galactosidase, a putative spermine and (or) spermidine synthase (E.C.2.5.1.16), and an unknown expressed gene. Strikingly, sequencing data revealed that the 10.91-kb CHS1, CHS3, CHS4 cluster is present as a perfect inverted repeat separated by 5.87 kb. Contiguous arrangement of CHS paralogs could lead to folding into multiple secondary structures, hypothesized to induce deletions that have previously been shown to effect CHS expression. BAC104J7 also contains several gene fragments representing a cation/hydrogen exchanger, a 40S ribosomal protein, a CBL-interacting protein kinase, and the amino terminus of a subtilisin. Chimeric ESTs were identified that may represent read-through transcription from a flanking truncated gene into a CHS cluster, generating aberrant CHS RNA molecules that could play a role in CHS gene silencing.

Microarrays for global expression constructed with a low redundancy set of 27,500 sequenced cDNAs representing an array of developmental stages and physiological conditions of the soybean plant

BACKGROUND

Microarrays are an important tool with which to examine coordinated gene expression. Soybean (Glycine max) is one of the most economically valuable crop species in the world food supply. In order to accelerate both gene discovery as well as hypothesis-driven research in soybean, global expression resources needed to be developed. The applications of microarray for determining patterns of expression in different tissues or during conditional treatments by dual labeling of the mRNAs are unlimited. In addition, discovery of the molecular basis of traits through examination of naturally occurring variation in hundreds of mutant lines could be enhanced by the construction and use of soybean cDNA microarrays.

RESULTS

We report the construction and analysis of a low redundancy 'unigene' set of 27,513 clones that represent a variety of soybean cDNA libraries made from a wide array of source tissue and organ systems, developmental stages, and stress or pathogen-challenged plants. The set was assembled from the 5' sequence data of the cDNA clones using cluster analysis programs. The selected clones were then physically reracked and sequenced at the 3' end. In order to increase gene discovery from immature cotyledon libraries that contain abundant mRNAs representing storage protein gene families, we utilized a high density filter normalization approach to preferentially select more weakly expressed cDNAs. All 27,513 cDNA inserts were amplified by polymerase chain reaction. The amplified products, along with some repetitively spotted control or 'choice' clones, were used to produce three 9,728-element microarrays that have been used to examine tissue specific gene expression and global expression in mutant isolines.

CONCLUSIONS

Global expression studies will be greatly aided by the availability of the sequence-validated and low redundancy cDNA sets described in this report. These cDNAs and ESTs represent a wide array of developmental stages and physiological conditions of the soybean plant. We also demonstrate that the quality of the data from the soybean cDNA microarrays is sufficiently reliable to examine isogenic lines that differ with respect to a mutant phenotype and thereby to define a small list of candidate genes potentially encoding or modulated by the mutant phenotype.

Floral dip: agrobacterium-mediated germ line transformation

Many researchers use the flowering plant Arabidopsis thaliana to study gene function and basic plant biology. This easy-to-grow, small plant is ideal for genetic studies as it has a relatively simple genome compared to crop plants and its genetic material has been recently sequenced. Another very useful feature of Arabidopsis is that it is extremely simple to transform genetically. The ability to insert genes of interest stably into a given plant is essential to understand and verify gene function. Transformation is also a means of introducing specific traits that are difficult or impossible to introduce by conventional breeding techniques. This chapter provides detailed explanations on the floral dip protocol, a simple method to transform Arabidopsis by inoculating immature flowers with Agrobacterium tumefaciens.

Tissue-specific gene silencing mediated by a naturally occurring chalcone synthase gene cluster in Glycine max

Chalcone synthase, a key regulatory enzyme in the flavonoid pathway, constitutes an eight-member gene family in Glycine max (soybean). Three of the chalcone synthase (CHS) gene family members are arranged as inverted repeats in a 10-kb region, corresponding to the I locus (inhibitor). Spontaneous mutations of a dominant allele (I or i(i)) to a recessive allele (i) have been shown to delete promoter sequences, paradoxically increasing total CHS transcript levels and resulting in black seed coats. However, it is not known which of the gene family members contribute toward pigmentation and how this locus affects CHS expression in other tissues. We investigated the unusual nature of the I locus using four pairs of isogenic lines differing with respect to alleles of the I locus. RNA gel blots using a generic open reading frame CHS probe detected similar CHS transcript levels in stems, roots, leaves, young pods, and cotyledons of the yellow and black isolines but not in the seed coats, which is consistent with the dominant I and i(i) alleles mediating CHS gene silencing in a tissue-specific manner. Using real-time RT-PCR, a variable pattern of expression of CHS genes in different tissues was demonstrated. However, increase in pigmentation in the black seed coats was associated with release of the silencing effect specifically on CHS7/CHS8, which occurred at all stages of seed coat development. These expression changes were linked to structural changes taking place at the I locus, shown to encompass a much wider region of at least 27 kb, comprising two identical 10.91-kb stretches of CHS gene duplications. The suppressive effect of this 27-kb I locus in a specific tissue of the G. max plant represents a unique endogenous gene silencing mechanism.

A second T-region of the soybean-supervirulent chrysopine-type Ti plasmid pTiChry5, and construction of a fully disarmed vir helper plasmid

Agrobacterium tumefaciens Chry5, which is particularly virulent on soybeans, induces tumors that produce a family of Amadori-type opines that includes deoxyfructosyl glutamine (Dfg) and its lactone, chrysopine (Chy). Cosmid clones mapping to the right of the known oncogenic T-region of pTiChry5 conferred Amadori opine production on tumors induced by the nopaline strain C58. Sequence analysis of DNA held in common among these cosmids identified two 25-bp, direct repeats flanking an 8.5-kb segment of pTiChry5. These probable border sequences are closely related to those of other known T-regions and define a second T-region of pTiChry5, called T-right (TR), that confers production of the Amadoriopines. The oncogenic T-left region (TL) was located precisely by identifying and sequencing the likely border repeats defining this segment. The two T-regions are separated by approximately 15 kb of plasmid DNA. Based on these results, we predicted that pKYRT1, a vir helper plasmid derived from pTiChry5, still contains all of TR and the leftmost 9 kb of TL. Consistent with this hypothesis, transgenic Arabidopsis thaliana plants selected for with a marker encoded by a binary plasmid following transformation with KYRT1 co-inherited production of the Amadori opines at high frequency. All opine-positive transgenic plants also contained TR-DNA, while those plants that lacked TR-DNA failed to produce the opines. Moreover, A. thaliana infected with KYRT1 in which an nptII gene driven by the 35S promoter of Cauliflower mosaic virus was inserted directly into the vir helper plasmid yielded kanamycin-resistant transformants at a low but detectable frequency. These results demonstrate that pKYRT1 is not disarmed, and can transfer Ti plasmid DNA to plants. A new vir helper plasmid was constructed from pTiChry5 by two rounds of sacB-mediated selection for deletion events. This plasmid, called pKPSF2, lacks both of the known T-regions and their borders. pKPSF2 failed to transfer Ti plasmid DNA to plants, but mobilized the T-region of a binary plasmid at an efficiency indistinguishable from those of pKYRT1 and the nopaline-type vir helper plasmid pMP90.

The Arabidopsis dnd1 "defense, no death" gene encodes a mutated cyclic nucleotide-gated ion channel

Gene-for-gene disease resistance typically includes a programmed cell death response known as the hypersensitive response (HR). The Arabidopsis thaliana dnd1 mutant was previously isolated as a line that failed to produce the HR in response to avirulent Pseudomonas syringae pathogens; plants homozygous for the recessive dnd1-1 mutation still carry out effective gene-for-gene resistance. The dnd1-1 mutation also causes constitutive systemic resistance and elevated levels of salicylic acid. In the present study, a positional cloning approach was used to isolate DND1. DND1 encodes the same protein as AtCNGC2, a cyclic nucleotide-gated ion channel of previously unknown organismal function that can allow passage of Ca(2+), K(+) and other cations [Leng, Q., Mercier, R. W., Yao, W. & Berkowitz, G. A. (1999) Plant Physiol. 121, 753-761]. By using a nahG transgene, we found that salicylic acid is required for the elevated resistance caused by the dnd1 mutation but that removal of salicylic acid did not completely eliminate the dwarf and loss-of-HR phenotypes of mutant dnd1 plants. A stop codon that would severely truncate the DND1 gene product was identified in the dnd1-1 allele. This demonstrates that broad-spectrum disease resistance and inhibition of the HR can be activated in plants by disruption of a cyclic nucleotide-gated ion channel.

Female reproductive tissues are the primary target of Agrobacterium-mediated transformation by the Arabidopsis floral-dip method

The floral-dip method for Agrobacterium-mediated transformation of Arabidopsis allows efficient plant transformation without need for tissue culture. To facilitate use with other plant species, we investigated the mechanisms that underlie this method. In manual outcrossing experiments, application of Agrobacterium tumefaciens to pollen donor plants did not produce any transformed progeny, whereas application of Agrobacterium to pollen recipient plants yielded transformants at a rate of 0.48%. Agrobacterium strains with T-DNA carrying gusA (encoding beta-glucuronidase [GUS]) under the control of 35S, LAT52, or ACT11 promoters revealed delivery of GUS activity to developing ovules, whereas no GUS staining of pollen or pollen tubes was observed. Transformants derived from the same seed pod contained independent T-DNA integration events. In Arabidopsis flowers, the gynoecium develops as an open, vase-like structure that fuses to form closed locules roughly 3 d prior to anthesis. In correlation with this fact, we found that the timing of Agrobacterium infection was critical. Transformants were obtained and GUS staining of ovules and embryo sacs was observed only if the Agrobacterium were applied 5 d or more prior to anthesis. A 6-fold higher rate of transformation was obtained with a CRABS-CLAW mutant that maintains an open gynoecium. Our results suggest that ovules are the site of productive transformation in the floral-dip method, and further suggest that Agrobacterium must be delivered to the interior of the developing gynoecium prior to locule closure if efficient transformation is to be achieved.

Identification of Arabidopsis mutants exhibiting an altered hypersensitive response in gene-for-gene disease resistance

A mutational study was carried out to isolate Arabidopsis thaliana plants that exhibit full or partial disruption of the RPS2-mediated hypersensitive response (HR) to Pseudomonas syringae that express avrRpt2. Five classes of mutants were identified including mutations at RPS2, dnd mutations causing a "defense, no death" loss-of-HR phenotype, a lesion-mimic mutant that also exhibited an HR- phenotype, and a number of intermediate or partial-loss-of-HR mutants. Surprisingly, many of these mutants displayed elevated resistance to virulent P. syringae and, in some cases, to Peronospora parasitica. Constitutively elevated levels of pathogenesis-related (PR) gene expression and salicylic acid were also observed. In the lesion-mimic mutant, appearance of elevated resistance was temporally correlated with appearance of lesions. For one of the intermediate lines, resistance was shown to be dependent on elevated levels of salicylic acid. A new locus was identified and named IHR1, after the mutant phenotype of "intermediate HR." Genetic analysis of the intermediate-HR plant lines was difficult due to uncertainties in distinguishing the partial/intermediate mutant phenotypes from wild type. Despite this difficulty, the intermediate-HR mutants remain of interest because they reveal potential new defense-related loci and because many of these lines exhibit partially elevated disease resistance without dwarfing or other apparent growth defects.

Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana

The Agrobacterium vacuum infiltration method has made it possible to transform Arabidopsis thaliana without plant tissue culture or regeneration. In the present study, this method was evaluated and a substantially modified transformation method was developed. The labor-intensive vacuum infiltration process was eliminated in favor of simple dipping of developing floral tissues into a solution containing Agrobacterium tumefaciens, 5% sucrose and 500 microliters per litre of surfactant Silwet L-77. Sucrose and surfactant were critical to the success of the floral dip method. Plants inoculated when numerous immature floral buds and few siliques were present produced transformed progeny at the highest rate. Plant tissue culture media, the hormone benzylamino purine and pH adjustment were unnecessary, and Agrobacterium could be applied to plants at a range of cell densities. Repeated application of Agrobacterium improved transformation rates and overall yield of transformants approximately twofold. Covering plants for 1 day to retain humidity after inoculation also raised transformation rates twofold. Multiple ecotypes were transformable by this method. The modified method should facilitate high-throughput transformation of Arabidopsis for efforts such as T-DNA gene tagging, positional cloning, or attempts at targeted gene replacement.

Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana

The Agrobacterium vacuum infiltration method has made it possible to transform Arabidopsis thaliana without plant tissue culture or regeneration. In the present study, this method was evaluated and a substantially modified transformation method was developed. The labor-intensive vacuum infiltration process was eliminated in favor of simple dipping of developing floral tissues into a solution containing Agrobacterium tumefaciens, 5% sucrose and 500 microliters per litre of surfactant Silwet L-77. Sucrose and surfactant were critical to the success of the floral dip method. Plants inoculated when numerous immature floral buds and few siliques were present produced transformed progeny at the highest rate. Plant tissue culture media, the hormone benzylamino purine and pH adjustment were unnecessary, and Agrobacterium could be applied to plants at a range of cell densities. Repeated application of Agrobacterium improved transformation rates and overall yield of transformants approximately twofold. Covering plants for 1 day to retain humidity after inoculation also raised transformation rates twofold. Multiple ecotypes were transformable by this method. The modified method should facilitate high-throughput transformation of Arabidopsis for efforts such as T-DNA gene tagging, positional cloning, or attempts at targeted gene replacement.

Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana

The Agrobacterium vacuum infiltration method has made it possible to transform Arabidopsis thaliana without plant tissue culture or regeneration. In the present study, this method was evaluated and a substantially modified transformation method was developed. The labor-intensive vacuum infiltration process was eliminated in favor of simple dipping of developing floral tissues into a solution containing Agrobacterium tumefaciens, 5% sucrose and 500 microliters per litre of surfactant Silwet L-77. Sucrose and surfactant were critical to the success of the floral dip method. Plants inoculated when numerous immature floral buds and few siliques were present produced transformed progeny at the highest rate. Plant tissue culture media, the hormone benzylamino purine and pH adjustment were unnecessary, and Agrobacterium could be applied to plants at a range of cell densities. Repeated application of Agrobacterium improved transformation rates and overall yield of transformants approximately twofold. Covering plants for 1 day to retain humidity after inoculation also raised transformation rates twofold. Multiple ecotypes were transformable by this method. The modified method should facilitate high-throughput transformation of Arabidopsis for efforts such as T-DNA gene tagging, positional cloning, or attempts at targeted gene replacement.

Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana

The Agrobacterium vacuum infiltration method has made it possible to transform Arabidopsis thaliana without plant tissue culture or regeneration. In the present study, this method was evaluated and a substantially modified transformation method was developed. The labor-intensive vacuum infiltration process was eliminated in favor of simple dipping of developing floral tissues into a solution containing Agrobacterium tumefaciens, 5% sucrose and 500 microliters per litre of surfactant Silwet L-77. Sucrose and surfactant were critical to the success of the floral dip method. Plants inoculated when numerous immature floral buds and few siliques were present produced transformed progeny at the highest rate. Plant tissue culture media, the hormone benzylamino purine and pH adjustment were unnecessary, and Agrobacterium could be applied to plants at a range of cell densities. Repeated application of Agrobacterium improved transformation rates and overall yield of transformants approximately twofold. Covering plants for 1 day to retain humidity after inoculation also raised transformation rates twofold. Multiple ecotypes were transformable by this method. The modified method should facilitate high-throughput transformation of Arabidopsis for efforts such as T-DNA gene tagging, positional cloning, or attempts at targeted gene replacement.

Identification of 3-hydroxypalmitic acid methyl ester as a novel autoregulator controlling virulence in Ralstonia solanacearum

Expression of virulence genes in Ralstonia solanacearum, a phytopathogenic bacterium, is controlled by a complex regulatory network that integrates multiple signal inputs. Production of several virulence determinants is coordinately reduced by inactivation of phcB, but is restored by growth in the presence of a volatile extracellular factor (VEF) produced by wild-type strains of R. solanacearum. The VEF was purified from spent culture broth by distillation, solvent extraction, and liquid chromatography. Gas chromatography and mass spectroscopy identified 3-hydroxypalmitic acid methyl ester (3-OH PAME) as the major component in the single peak of VEF activity. Authentic 3-OH PAME and the purified VEF were active at < or =1 nM, and had nearly equivalent specific activities for stimulating the expression of eps (the biosynthetic locus for extracellular polysaccharide) in a phcB mutant. Authentic 3-OH PAME also increased the production of three virulence factors by a phcB mutant over 20-fold to wild-type levels, restored normal cell density-associated expression of eps and increased expression of eps when delivered via the vapour phase. Reanalysis of the PhcB amino acid sequence suggested that it is a small-molecule S-adenosylmethionine-dependent methyltransferase, which might catalyse synthesis of 3-OH PAME from a naturally occurring fatty acid. Biologically active concentrations of extracellular 3-OH PAME were detected before the onset of eps expression, suggesting that it is an intercellular signal that autoregulates virulence gene expression in wild-type R. solanacearum. Other than acyl-homoserine lactones, 3-OH PAME is the only endogenous fatty acid derivative shown to be an autoregulator and may be the first example of a new family of compounds that can mediate long-distance intercellular communication.

A two-component system in Ralstonia (Pseudomonas) solanacearum modulates production of PhcA-regulated virulence factors in response to 3-hydroxypalmitic acid methyl ester

Expression of virulence factors in Ralstonia solanacearum is controlled by a complex regulatory network, at the center of which is PhcA, a LysR family transcriptional regulator. We report here that expression of phcA and production of PhcA-regulated virulence factors are affected by products of the putative operon phcBSR(Q). phcB is required for production of an extracellular factor (EF), tentatively identified as the fatty acid derivative 3-hydroxypalmitic acid methyl ester (3-OH PAME), but a biochemical function for PhcB could not be deduced from DNA sequence analysis. The other genes in the putative operon are predicted to encode proteins homologous to members of two-component signal transduction systems: PhcS has amino acid similarity to histidine kinase sensors, whereas PhcR and OrfQ are similar to response regulators. PhcR is quite unusual because its putative output domain strongly resembles the histidine kinase domain of a sensor protein. Production of the PhcA-regulated factors exopolysaccharide I, endoglucanase, and pectin methyl esterase was reduced 10- to 100-fold only in mutants with a nonpolar insertion in phcB [which express phcSR(Q) in the absence of the EF]; simultaneously, expression of phcA was reduced fivefold. Both a wild-type phenotype and phcA expression were restored by addition of 3-OH PAME to growing cultures. Mutants with polar insertions in phcB or lacking the entire phcBSR(Q) region produced wild-type levels of PhcA-regulated virulence factors. The genetic data suggest that PhcS and PhcR function together to regulate expression of phcA, but the biochemical mechanism for this is unclear. At low levels of the EF, it is likely that PhcS phosphorylates PhcR, and then PhcR interacts either with PhcA (which is required for full expression of phcA) or an unknown component of the signal cascade to inhibit expression of phcA. When the EF reaches a threshold concentration, we suggest that it reduces the ability of PhcS to phosphorylate PhcR, resulting in increased expression of phcA and production of PhcA-regulated factors.