Flavonoids Modulate the Accumulation of Toxins From Aspergillus flavus in Maize Kernels

Aspergillus flavus is an opportunistic fungal pathogen capable of producing aflatoxins, potent carcinogenic toxins that accumulate in maize kernels after infection. To better understand the molecular mechanisms of maize resistance to A. flavus growth and aflatoxin accumulation, we performed a high-throughput transcriptomic study in situ using maize kernels infected with A. flavus strain 3357. Three maize lines were evaluated: aflatoxin-contamination resistant line TZAR102, semi-resistant MI82, and susceptible line Va35. A modified genotype-environment association method (GEA) used to detect loci under selection via redundancy analysis (RDA) was used with the transcriptomic data to detect genes significantly influenced by maize line, fungal treatment, and duration of infection. Gene ontology enrichment analysis of genes highly expressed in infected kernels identified molecular pathways associated with defense responses to fungi and other microbes such as production of pathogenesis-related (PR) proteins and lipid bilayer formation. To further identify novel genes of interest, we incorporated genomic and phenotypic field data from a genome wide association analysis with gene expression data, allowing us to detect significantly expressed quantitative trait loci (eQTL). These results identified significant association between flavonoid biosynthetic pathway genes and infection by A. flavus. In planta fungal infections showed that the resistant line, TZAR102, has a higher fold increase of the metabolites naringenin and luteolin than the susceptible line, Va35, when comparing untreated and fungal infected plants. These results suggest flavonoids contribute to plant resistance mechanisms against aflatoxin contamination through modulation of toxin accumulation in maize kernels.

Functional Analysis of A Soybean Ferredoxin-NADP Reductase (FNR) Gene in Response to Soybean Mosaic Virus

The Ferredoxin-NADP reductase (FNR) gene plays a significant role in NADPH production, carbon assimilation, antioxidation, and cross-talking between chloroplasts and mitochondria in plants. This study aims to know the functional response of the soybean FNR gene (GmFNR) during a soybean mosaic virus (SMV) infection. For this purpose, we developed the bean pod mottle virus (BPMV)-based gene construct (BPMV-GmFNR) and used it to silence the GmFNR gene in resistant and susceptible lines. The results showed that GmFNR expression decreased to 50% in the susceptible line, compared to 40% in the resistant line. The silencing of GmFNR reduces the photosynthetic capacity and CAT activity of both lines compared to their respective controls. In addition, the H2O2 content increased significantly in the susceptible line, whereas the resistant line did not exhibit any change. Further, an SMV infection in the silencing plants of the susceptible line resulted in serious morphological changes and increased the SMV NIa-protease transcript accumulation compared to its control plants. However, the same impact was not observed in the resistant line. The yeast two-hybrid system, BIFC assay, and quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed that the GmFNR was interacting with EF1A and coincided with the increased SMV accumulation. The results obtained in this study improve the understanding of the soybean FNR gene response during SMV infection and provide a novel insight into the SMV resistance mechanism.

Transcriptome Analysis of Tolerant and Susceptible Maize Genotypes Reveals Novel Insights about the Molecular Mechanisms Underlying Drought Responses in Leaves

Maize (Zea mays L.) is the most essential food crop in the world. However, maize is highly susceptible to drought stress, especially at the seedling stage, and the molecular mechanisms underlying drought tolerance remain elusive. In this study, we conducted comparative transcriptome and physiological analyses of drought-tolerant (CML69) and susceptible (LX9801) inbred lines subjected to drought treatment at the seedling stage for three and five days. The tolerant line had significantly higher relative water content in the leaves, as well as lower electrolyte leakage and malondialdehyde levels, than the susceptible line. Using an RNA-seq-based approach, we identified 10,084 differentially expressed genes (DEGs) with 6906 and 3178 DEGs been annotated and unannotated, respectively. Two critical sets of drought-responsive DEGs, including 4687 genotype-specific and 2219 common drought-responsive genes, were mined out of the annotated DEGs. The tolerant-line DEGs were predominantly associated with the cytoskeleton, cell wall modification, glycolysis/gluconeogenesis, transport, osmotic regulation, drought avoidance, ROS scavengers, defense, and transcriptional factors. For the susceptible line, the DEGs were highly enriched in the photosynthesis, histone, and carbon fixation pathways. The unannotated DEGs were implicated in lncRNAs, including 428 previously reported and 22% putative TE-lncRNAs. There was consensus on both the physiological response and RNA-seq outcomes. Collectively, our findings will provide a comprehensive basis of the molecular networks mediating drought stress tolerance of maize at the seedling stage.

Identification of long noncoding RNAs involved in resistance to downy mildew in Chinese cabbage

Brassica downy mildew, a severe disease caused by Hyaloperonospora brassicae, can cause enormous economic losses in Chinese cabbage (Brassica rapa L. ssp. pekinensis) production. Although some research has been reported recently concerning the underlying resistance to this disease, no studies have identified or characterized long noncoding RNAs involved in this defense response. In this study, using high-throughput RNA sequencing, we analyzed the disease-responding mRNAs and long noncoding RNAs in two resistant lines (T12–19 and 12–85) and one susceptible line (91–112). Clustering and Gene Ontology analysis of differentially expressed genes (DEGs) showed that more DEGs were involved in the defense response in the two resistant lines than in the susceptible line. Different expression patterns and proposed functions of differentially expressed long noncoding RNAs among T12–19, 12–85, and 91–112 indicated that each has a distinct disease response mechanism. There were significantly more cis- and trans-functional long noncoding RNAs in the resistant lines than in the susceptible line, and the genes regulated by these RNAs mostly participated in the disease defense response. Furthermore, we identified a candidate resistance-related long noncoding RNA, MSTRG.19915, which is a long noncoding natural antisense transcript of a MAPK gene, BrMAPK15. Via an agroinfiltration-mediated transient overexpression system and virus-induced gene silencing technology, BrMAPK15 was indicated to have a greater ability to defend against pathogens. MSTRG.19915-silenced seedlings showed enhanced resistance to downy mildew, probably because of the upregulated expression of BrMAPK15. This research identified and characterized long noncoding RNAs involved in resistance to downy mildew, laying a foundation for future in-depth studies of disease resistance mechanisms in Chinese cabbage.

Glucosinolate profile and Myrosinase gene expression are modulated upon Plasmodiophora brassicae infection in cabbage

Clubroot is a devastating disease of Brassicaceae caused by the biotrophic protist Plasmodiophora brassicae. The progression of clubroot disease is modulated by the glucosinolate (GSL) profile of the host plant. GSL is hydrolysed by the enzyme myrosinase upon cell disruption and gives rise to metabolites like isothiocyanate, nitriles, thiocyanates, epithionitriles and oxazolidines. Some of these metabolites play important roles in the plant’s defence mechanism. We identified 13 Myrosinase (Myro) and 28 Myrosinase-Binding Protein-like (MBP) genes from Brassica oleracea L. using a comparative genomics approach and characterised them through in silico analyses. We compared the expression patterns of these genes in a clubroot-susceptible line and a resistant line following inoculation with P. brassicae. Two BolMyro and 12 BolMBP genes were highly expressed in the susceptible line, whereas only one BolMyro and five BolMBP genes were highly expressed in the resistant line. Principal component analysis confirmed that specific GSL profiles and gene expression were modulated due to pathogen infection. Plants with higher levels of neoglucobrassicin, glucobrassicin and methooxyglucobrassicin produced disease symptoms and formed galls, whereas, plants with higher levels of sinigrin, hydroxyglucobrassicin and progoitrin produced less symptoms with almost no galls. Our results provide insights into the roles of Myro and MBP genes in GSL hydrolysis during P. brassicae infection, which will help for developing clubroot resistant cabbage lines.

Global mRNA and microRNA expression dynamics in response to anthracnose infection in sorghum

Background Anthracnose is a damaging disease of sorghum caused by the fungal pathogen Colletotrichum sublineolum. Genome-wide mRNA and microRNA (miRNA) profiles of resistant and susceptible sorghum genotypes were studied to understand components of immune responses, and fungal induced miRNA and target gene networks. Results A total of 18 mRNA and 12 miRNA libraries from resistant and susceptible sorghum lines were sequenced prior to and after inoculation with C. sublineolum. Significant differences in transcriptomes of the susceptible and resistant genotypes were observed with dispersion distance and hierarchical cluster tree analyses. Of the total 33,032 genes predicted in the sorghum genome, 19,593 were induced by C. sublineolum, and 15,512 were differentially expressed (DEGs) between the two genotypes. The resistant line was marked by significant reprogramming of the transcriptome at 24 h post inoculation (hpi), and a decrease at 48 hpi, whereas the susceptible line displayed continued changes in gene expression concordant with elevated fungal growth in the susceptible genotype. DEGs encode proteins implicated in diverse functions including photosynthesis, synthesis of tetrapyrrole, carbohydrate and secondary metabolism, immune signaling, and chitin binding. Genes encoding immune receptors, MAPKs, pentatricopeptide repeat proteins, and WRKY transcription factors were induced in the resistant genotype. In a parallel miRNA profiling, the susceptible line displayed greater number of differentially expressed miRNAs than the resistant line indicative of a widespread suppression of gene expression. Interestingly, we found 75 miRNAs, including 36 novel miRNAs, which were differentially expressed in response to fungal inoculation. The expression of 50 miRNAs was significantly different between resistant and susceptible lines. Subsequently, for 35 differentially expressed miRNAs, the corresponding 149 target genes were identified. Expression of 56 target genes were significantly altered after inoculation, showing inverse expression with the corresponding miRNAs. Conclusions We provide insights into genome wide dynamics of mRNA and miRNA profiles, biological and cellular processes underlying host responses to fungal infection in sorghum. Resistance is correlated with early transcriptional reprogramming of genes in various pathways. Fungal induced genes, miRNAs and their targets with a potential function in host responses to anthracnose were identified, opening avenues for genetic dissection of resistance mechanisms.

The pod shattering resistance of soybean lines based on the shattering incidence and severity

The study is aimed at evaluating the pod shattering resistance of F₈ soybean lines based on the shattering incidence and shattering severity. The materials consist of fourteen F₈ soybean lines and two check cultivars. The pod shattering incidence was examined by using the oven-dry method, meanwhile, the shattering severity was evaluated based on the severity of the pod opening. The pod shattering resistance based on the shattering incidence resulted in five resistant lines (7–10% shattering), seven moderate lines (13–23% shattering), one susceptible line (53% shattering), and one very susceptible line (100% shattering). The pod shattering resistance based on the shattering severity showed that the pod opening on the ventral side differed between the lines and between the shattering degree, and it tends to form sigmoid curves with a different peak position for each shattering degree. The shattering severity of the resistant, moderate, and susceptible lines reached a peak at 60 °C, 50 °C, and 40 °C, respectively. A longer pod length indicated by the length of the dorsal (r = 0.827**) and ventral (r = 0.880**) sides of the pod, a higher total pod weight (0.827**), and a larger seed size (0.794**) will increase the degree of susceptibility to pod shattering. Those characteristics were considered to be the ones that should be used as the selection criteria in the breeding programme for pod shattering resistance in soybeans.

Transcriptome-based Analysis of Tomato Genotypes Resistant to Bacterial spot (Xanthomonas perforans) Race T4

Bacterial spot (BS) is one of the most devastating foliar bacterial diseases of tomato caused by multiple species of Xanthomonas. We performed the RNA-Seq analysis of three tomato lines with different level of resistance to Xanthomonas perforans race T4 to study the differentially expressed genes (DEGs) and transcript-based sequence variations.Analysis between inoculated and control samples revealed that resistant line PI 270443 had more DEGs (834), followed by susceptible line NC 714 (373), and intermediate line NC 1CELBR (154). Gene functional analysis based on Gene Ontology (GO) terms revealed that more GO terms (51) were enriched for up-regulated DEGs in the resistant line PI 270443, and more down-regulated DEGs (67) were enriched in the susceptible line NC 714. The specific analysis for DEGs in biotic stress pathway using MapMan software showed more up-regulated biotic stress pathway DEGs (67) for PI 270443 compared to more down-regulated DEGs (125) for susceptible NC 714 line. One interesting feature was that resistant PI 270443 has three up-regulated DEGs for PR-protein, and susceptible line NC 714 has one down-regulated R gene, which is disease-related.Analysis of sequence variations called from RNA-Seq reads against the reference genome of susceptible Heinz 1706 showed that chr11 which has multiple reported resistance QTLs to BS race T4 is identical between two resistant lines, PI 270443 and NC 1CELBR, suggesting that these two lines share the same resistance QTLs on this chromosome. Several loci for PR-resistance proteins with sequence variation between the resistant and susceptible tomato lines were identified near the known Rx4 resistance gene on chr11. These findings may be useful for further molecular breeding of tomato.

Marek’s Disease Virus Infection Induced Mitochondria Changes in Chickens

Mitochondria are crucial cellular organelles in eukaryotes and participate in many cell processes including immune response, growth development, and tumorigenesis. Marek’s disease (MD), caused by an avian alpha-herpesvirus Marek’s disease virus (MDV), is characterized with lymphomas and immunosuppression. In this research, we hypothesize that mitochondria may play roles in response to MDV infection. To test it, mitochondrial DNA (mtDNA) abundance and gene expression in immune organs were examined in two well-defined and highly inbred lines of chickens, the MD-susceptible line 72 and the MD-resistant line 63. We found that mitochondrial DNA contents decreased significantly at the transformation phase in spleen of the MD-susceptible line 72 birds in contrast to the MD-resistant line 63. The mtDNA-genes and the nucleus-genes relevant to mtDNA maintenance and transcription, however, were significantly up-regulated. Interestingly, we found that POLG2 might play a potential role that led to the imbalance of mtDNA copy number and gene expression alteration. MDV infection induced imbalance of mitochondrial contents and gene expression, demonstrating the indispensability of mitochondria in virus-induced cell transformation and subsequent lymphoma formation, such as MD development in chicken. This is the first report on relationship between virus infection and mitochondria in chicken, which provides important insights into the understanding on pathogenesis and tumorigenesis due to viral infection.

Variation in tolerance mechanisms to fluazifop-P-butyl among selected zoysiagrass lines

Breeding herbicide tolerance into new cultivars can improve safety and weed control in turfgrass systems. The sensitivity to fluazifop-P-butyl of 27 zoysiagrass (Zoysiaspp.) lines was screened under greenhouse conditions to identify potential tolerant germplasm for breeding programs. The herbicide rate that caused 50% biomass reduction (GR50) and the rate that caused 50% injury (ID50) were calculated to select the three most-tolerant and the five most-susceptible lines for studying the physiological mechanisms responsible for fluazifop-P-butyl tolerance. The differences in GR50and ID50between susceptible and tolerant lines ranged from 4-fold to more than 10-fold. Cytochrome P450–mediated metabolism was not detected in fluazifop-P-butyl–tolerant lines. Sequencing of theACCasegene confirmed that none of the seven previously reported mutations conferring resistance to acetyl-CoA carboxylase (ACCase)-inhibiting herbicides in other species were present in any of the tolerant or susceptible zoysiagrass lines studied. An Ala-2073-Thr substitution was identified in two tolerant lines, but this mutation did not completely explain the tolerant phenotype. No clear differences in absorption and translocation rates of14C-radiolabeled fluazifop-P-butyl were observed among most lines, with the exception of a susceptible line that exhibited greater translocation than two of the tolerant lines. Metabolite profiles did not differ between tolerant and susceptible lines. Our results suggest that the diversity in tolerance to fluazifop-P-butyl in zoysiagrass germplasm is most likely the result of a combination of different, minor, additive non–target site mechanisms such as translocation rate and compartmentation after absorption.