Insight into the composition and differentiation of endophytic microbial communities in kernels via 368 maize transcriptomes.

INTRODUCTION: Kernels are important reproductive organs in maize, yet there is a lack of systematic investigation on the differences in the composition of endophytic microorganisms in plants from a population perspective. OBJECTIVES: We aimed to elucidate the composition of endophytic microorganisms in developing maize kernels, emphasizing differences among various inbred lines. METHODS: The transcriptomic data of 368 maize inbred lines were used to explore the composition and diversity of endophytic microorganisms. RESULTS: The findings revealed a higher abundance of fungi than bacteria in developing maize kernels, followed by protozoa, while viruses were less abundant. There were significant differences in the composition and relative abundance of endophytic microorganisms among different maize lines. Diversity analysis revealed overall similarity in the community composition structure between tropical/subtropical (TST) and temperate (NSS) maize germplasm with apparent variations in community richness and abundance. The endophytic microorganisms network in the kernels from TST genotypes exhibited higher connectivity and stability compared to NSS kernels. Bacteria dominated the highly connected species in the networks, and different core species showed microbial phylum specificity. Some low-abundance species acted as core species, contributing to network stability. Beneficial bacteria were predominant in the core species of networks in TST kernels, while pathogenic bacteria were more abundant in the core species of networks in NSS kernels. CONCLUSION: Tropical maize germplasm may have advantages in resisting the invasion of pathogenic microorganisms, providing excellent genetic resources for disease-resistant breeding.

HacA, a key transcription factor for the unfolded protein response, is required for fungal development, aflatoxin biosynthesis and pathogenicity of Aspergillus flavus.

Aspergillus flavus is a fungus notorious for contaminating food and feed with aflatoxins. As a saprophytic fungus, it secretes large amounts of enzymes to access nutrients, making endoplasmic reticulum (ER) homeostasis important for protein folding and secretion. The role of HacA, a key transcription factor in the unfolded protein response pathway, remains poorly understood in A. flavus. In this study, the hacA gene in A. flavus was knockout. Results showed that the absence of hacA led to a decreased pathogenicity of the strain, as it failed to colonize intact maize kernels. This may be due to retarded vegetable growth, especially the abnormal development of swollen tips and shorter hyphal septa. Deletion of hacA also hindered conidiogenesis and sclerotial development. Notably, the mutant strain failed to produce aflatoxin B1. Moreover, compared to the wild type, the mutant strain showed increased sensitivity to ER stress inducer such as Dithiothreitol (DTT), and heat stress. It also displayed heightened sensitivity to other environmental stresses, including cell wall, osmotic, and pH stresses. Further transcriptomic analysis revealed the involvement of the hacA in numerous biological processes, including filamentous growth, asexual reproduction, mycotoxin biosynthetic process, signal transduction, budding cell apical bud growth, invasive filamentous growth, response to stimulus, and so on. Taken together, HacA plays a vital role in fungal development, pathogenicity and aflatoxins biosynthesis. This highlights the potential of targeting hacA as a novel approach for early prevention of A. flavus contamination.

Metagenomic insights into the diversity of 2,4-dichlorophenol degraders and the cooperation patterns in a bacterial consortium.

2,4-Dichlorophenol, which is largely employed in herbicides and industrial production, is frequently detected in ecosystems and poses risks to human health and environmental safety. Microbial communities are thought to perform better than individual strains in the complete degradation of organic contaminants. However, the synergistic degradation mechanisms of the microbial consortia involved in 2,4-dichlorophenol degradation are still not widely understood. In this study, a bacterial consortium named DCP-2 that is capable of degrading 2,4-dichlorophenol was obtained. Metagenomic analysis, cultivation-dependent functional verification, and co-occurrence network analysis were combined to reveal the primary 2,4-dichlorophenol degraders and the cooperation patterns in the consortium DCP-2. Metagenomic analysis showed that Pseudomonas, Achromobacter, and Pigmentiphaga were the primary degraders for the complete degradation of 2,4-dichlorophenol. Thirty-nine phylogenetically diverse bacterial genera, such as Brucella, Acinetobacter, Aeromonas, Allochromatium and Bosea, were identified as keystone taxa for 2,4-dichlorophenol degradation by keystone taxa analysis of the co-occurrence networks. In addition, a stable synthetic consortium of isolates from DCP-2 was constructed, consisting of Pseudomonas sp. DD-13 and Brucella sp. FZ-1; this synthetic consortium showed superior degradation capability for 2,4-dichlorophenol in both mineral salt medium and wastewater compared with monoculture. The findings provide valuable insights into the practical bioremediation of 2,4-dichlorophenol-contaminated sites.

Biological control of wheat powdery mildew disease by the termite-associated fungus Aspergillus chevalieri BYST01 and potential role of secondary metabolites.

BACKGROUND: Wheat powdery mildew, caused by the biotrophic pathogen Blumeria graminis f. sp. tritici (Bgt) is a serious fungal disease. Natural metabolites produced by microorganisms are beneficial biological control agents to inhibit Bgt. In the present study, we investigated the effects of Aspergillus chevalieri BYST01 on wheat powdery mildew. RESULTS: A strain isolated from the termite was identified as A. chevalieri BYST01 by morphological characteristics and phylogenetic analysis. The fermentation broth of BYST01 showed good biocontrol effect on the Bgt in vivo with the control efficiencies of 81.59% and 71.34% under the protective and therapeutic tests, respectively. Four known metabolites, including the main compound physcion (30 mg/L), were isolated from the fermentation broth of BYST01 extracted with ethyl acetate. Importantly, under a concentration of 0.1 mM, physcion repressed conidial germination of Bgt with an inhibition rate of 77.04% in vitro and showed important control efficiencies of 80.36% and 74.64% in vivo under the protective and therapeutic tests, respectively. Hence, the BYST01 showed important potential as a microbial cell factory for the high yield of the green natural fungicide physcion. Finally, the biosynthetic gene clusters responsible for physicon production in BYST01 was predicted by analyzing a chromosome-scale genome obtained using a combination of Illumina, PacBio, and Hi-C sequencing technologies. CONCLUSION: Aspergillus chevalieri BYST01 and its main metabolite physcion had a significant control effect on wheat powdery mildew. The biosynthesis pathway of physcion in BYST01 was predicted. © 2023 Society of Chemical Industry.

Identification and Functional Characterization of ZmSCYL2 Involved in Phytosterol Accumulation in Plants.

Phytosterols are natural active substances widely found in plants and play an important role in hypolipidemia, antioxidants, antitumor, immunomodulation, plant growth, and development. In this study, phytosterols were extracted and identified from the seed embryos of 244 maize inbred lines. Based on this, a genome-wide association study (GWAS) was used to predict the possible candidate genes responsible for phytosterol content; 9 SNPs and 32 candidate genes were detected, and ZmSCYL2 was identified to be associated with phytosterol accumulation. We initially confirmed its functions in transgenic Arabidopsis and found that mutation of ZmSCYL2 resulted in slow plant growth and a significant reduction in sterol content, while overexpression of ZmSCYL2 accelerated plant growth and significantly increased sterol content. These results were further confirmed in transgenic tobacco and suggest that ZmSCYL2 was closely related to plant growth; overexpression of ZmSCYL2 not only facilitated plant growth and development but also promoted the accumulation of phytosterols.

Establishment and Validation of a New Analysis Strategy for the Study of Plant Endophytic Microorganisms.

Amplicon sequencing of bacterial or fungal marker sequences is currently the main method for the study of endophytic microorganisms in plants. However, it cannot obtain all types of microorganisms, including bacteria, fungi, protozoa, etc., in samples, nor compare the relative content between endophytic microorganisms and plants and between different types of endophytes. Therefore, it is necessary to develop a better analysis strategy for endophytic microorganism investigation. In this study, a new analysis strategy was developed to obtain endophytic microbiome information from plant transcriptome data. Results showed that the new strategy can obtain the composition of microbial communities and the relative content between plants and endophytic microorganisms, and between different types of endophytic microorganisms from the plant transcriptome data. Compared with the amplicon sequencing method, more endophytic microorganisms and relative content information can be obtained with the new strategy, which can greatly broaden the research scope and save the experimental cost. Furthermore, the advantages and effectiveness of the new strategy were verified with different analysis of the microbial composition, correlation analysis, inoculant content test, and repeatability test.

Influence of Oxygen Management during the Post-Fermentation Stage on Acetaldehyde, Color, and Phenolics of Vitis vinifera L. Cv. Cabernet Sauvignon Wine.

Oxygen exposure is unavoidable and the impact of its management during the post-fermentation stage (PFS) on dry red wine is poorly investigated. This study was dedicated to the variation of acetaldehyde, color and phenolics of Cabernet Sauvignon dry red wine during five discontinuous oxidation cycles of four levels of controlled oxygen supply, which were carried out to simulate probable oxidation during the PFS. Free SO2 disappeared after the first, second and third oxidation cycles in wines with high, medium and low levels of oxygen exposure severally, but subsequent oxygen exposure below or equal to 2 mg O2/L per cycle had little effect while 3-3.9 mg O2/L per cycle dramatically facilitated acetaldehyde accumulation, which was accompanied by an enormous variation in color and pigments, especially when total oxygen consumption was above 10 mg/L. The utilization of clustered heatmap and partial least square regression demonstrated the feasibility of characterization of wine oxidation degree using the chemical parameters measured by UV-spectrophotometry. Oxygen exposure during the PFS should be emphatically controlled, and chemical indexes determined by the UV-spectrophotometric method can be used for a scientific and effective description of wine oxidation degree.

The chromosome-scale genome provides insights into pigmentation in Acer rubrum.

Acer rubrum L. is one of the most prevalent ornamental species of the genus Acer, due to its straight and tall stems and beautiful leaf colors. For this study, the Oxford Nanopore platform and Hi-C technology were employed to obtain a chromosome-scale genome for A. rubrum. The genome size of A. rubrum was 1.69 Gb with an N50 of 549.44 Kb, and a total of 39 pseudochromosomes were generated with a 99.61% genome. The A. rubrum genome was predicted to have 64644 genes, of which 97.34% were functionally annotated. Genome annotation identified 67.14% as the transposable element (TE) repeat sequence, with long terminal repeats (LTR) being the richest (55.68%). Genome evolution analysis indicated that A. rubrum diverged from A. yangbiense ∼6.34 million years ago. We identified 13 genes related to pigment synthesis in A. rubrum leaves, where the expressions of four ArF3'H genes were consistent with the synthesis of cyanidin (a key pigment) in red leaves. Correlation analysis verified that the pigmentation of A. rubrum leaves was under the coordinated regulation of non-structural carbohydrates and hormones. The genomic sequence of A. rubrum will facilitate genomic breeding research for this species, while providing the valuable utilization of Aceraceae resources.

Linkage mapping combined with GWAS revealed the genetic structural relationship and candidate genes of maize flowering time-related traits.

BACKGROUND: Flowering time is an important agronomic trait of crops and significantly affects plant adaptation and seed production. Flowering time varies greatly among maize (Zea mays) inbred lines, but the genetic basis of this variation is not well understood. Here, we report the comprehensive genetic architecture of six flowering time-related traits using a recombinant inbred line (RIL) population obtained from a cross between two maize genotypes, B73 and Abe2, and combined with genome-wide association studies to identify candidate genes that affect flowering time. RESULTS: Our results indicate that these six traits showed extensive phenotypic variation and high heritability in the RIL population. The flowering time of this RIL population showed little correlation with the leaf number under different environmental conditions. A genetic linkage map was constructed by 10,114 polymorphic markers covering the whole maize genome, which was applied to QTL mapping for these traits, and identified a total of 82 QTLs that contain 13 flowering genes. Furthermore, a combined genome-wide association study and linkage mapping analysis revealed 17 new candidate genes associated with flowering time. CONCLUSIONS: In the present study, by using genetic mapping and GWAS approaches with the RIL population, we revealed a list of genomic regions and candidate genes that were significantly associated with flowering time. This work provides an important resource for the breeding of flowering time traits in maize.

Impact of the Acetaldehyde-Mediated Condensation on the Phenolic Composition and Antioxidant Activity of Vitis vinifera L. Cv. Merlot Wine.

Acetaldehyde is a critical reactant on modifying the phenolic profile during red wine aging, suggesting that the acetaldehyde-mediated condensation can be responsible for the variation of antioxidant activity during the aging of this beverage. The present study employs exogenous acetaldehyde at six levels of treatment (7.86 ± 0.10-259.02 ± 4.95 mg/L) before the bottle aging of Merlot wines to encourage phenolic modification. Acetaldehyde and antioxidant activity of wine were evaluated at 0, 15, 30, 45, 60 and 75 days of storage, while monomeric and polymeric phenolics were analyzed at 0, 30 and 75 days of storage. The loss of acetaldehyde was fitted to a first-order reaction model, the rate constant (k) demonstrated that different chemical reaction happened in wines containing a different initial acetaldehyde. The disappearance of monomeric phenolics and the formation of polymeric phenolics induced by acetaldehyde could be divided into two phases, the antioxidant activity of wine did not alter significantly in the first phase, although most monomeric phenolics vanished, but the second phase would dramatically reduce the antioxidant activity of wine. Furthermore, a higher level of acetaldehyde could shorten the reaction time of the first phase. These results indicate that careful vinification handling aiming at controlling the acetaldehyde allows one to maintain prolonged biological activity during wine aging.

Systematic Characterization of bZIP Transcription Factors Required for Development and Aflatoxin Generation by High-Throughput Gene Knockout in Aspergillus flavus.

The basic leucine zipper (bZIP) is an important transcription factor required for fungal development, nutrient utilization, biosynthesis of secondary metabolites, and defense against various stresses. Aspergillus flavus is a major producer of aflatoxin and an opportunistic fungus on a wide range of hosts. However, little is known about the role of most bZIP genes in A. flavus. In this study, we developed a high-throughput gene knockout method based on an Agrobacterium-mediated transformation system. Gene knockout construction by yeast recombinational cloning and screening of the null mutants by double fluorescence provides an efficient way to construct gene-deleted mutants for this multinucleate fungus. We deleted 15 bZIP genes in A. flavus. Twelve of these genes were identified and characterized in this strain for the first time. The phenotypic analysis of these mutants showed that the 15 bZIP genes play a diverse role in mycelial growth (eight genes), conidiation (13 genes), aflatoxin biosynthesis (10 genes), oxidative stress response (11 genes), cell wall stress (five genes), osmotic stress (three genes), acid and alkali stress (four genes), and virulence to kernels (nine genes). Impressively, all 15 genes were involved in the development of sclerotia, and the respective deletion mutants of five of them did not produce sclerotia. Moreover, MetR was involved in this biological process. In addition, HapX and MetR play important roles in the adaptation to excessive iron and sulfur metabolism, respectively. These studies provide comprehensive insights into the role of bZIP transcription factors in this aflatoxigenic fungus of global significance.

Identification and characterization of heat-responsive lncRNAs in maize inbred line CM1.

BACKGROUND: Frequent occurrence of extreme high temperature is a major threat to crop production. Increasing evidence demonstrates that long non-coding RNAs (lncRNAs) have important biological functions in the regulation of the response to heat stress. However, the regulatory mechanism of lncRNAs involved in heat response requires further exploration and the regulatory network remains poorly understood in maize. RESULTS: In this research, high-throughput sequencing was adopted to systematically identify lncRNAs in maize inbred line CM1. In total, 53,249 lncRNAs (259 known lncRNAs and 52,990 novel lncRNAs) were detected, of which 993 lncRNAs showed significantly differential expression (DElncRNAs) under heat stress. By predicting the target genes, 953 common targets shared by cis- and trans-regulation of the DElncRNAs were identified, which exhibited differential expression between the control and the heat stress treatments. Functional annotation indicated that a number of important biological processes and pathways, including photosynthesis, metabolism, translation, stress response, hormone signal transduction, and spliceosome, were enriched for the common targets, suggesting that they play important roles in heat response. A lncRNA-mediated regulatory network was constructed to visualize the molecular response mechanism in response to heat stress, which represented the direct regulatory relationships of DElncRNAs, differentially expressed miRNAs, target genes, and functional annotations. CONCLUSIONS: This study lays a foundation for further elucidation of the regulatory mechanism for the response to heat stress in the maize inbred line CM1. The findings provide important information for identification of heat-responsive genes, which will be beneficial for the molecular breeding in the cultivation of heat-tolerant maize germplasm.

Evaluation of Anthocyanin Profile and Color in Sweet Cherry Wine: Effect of Sinapic Acid and Grape Tannins during Aging.

Cherries are rich in bioactive phenolic compounds and are often fermented into cherry wines. The degradation of anthocyanins during storage will cause color deterioration. The study aimed to utilize sinapic acid and grape tannins in cherry wine to maintain a high fraction in the colored forms of anthocyanins, in order to maximize the color intensity, the latter being associated with good product quality. The effects on the anthocyanin profile and on color parameters of copigments, utilizing spectral measurement combined with UPLC-MS quantitative analysis, have been evaluated in sweet cherry wines. The copigmentation effect of sinapic acid and grape tannin was accompanied by the bathochromic shift and the hyperchromic effect, which lead to an increase in color intensity (lower L*, higher a* and b*). During the aging process, sinapic and grape tannin increased the content of pyranoanthocyanins in cherry wine, especially the addition of sinapic acid makes the cherry wine generate 10-syringyl-pyranocyanidin-3-rutinoside. These results demonstrate that sinapic acid is suitable for adding before alcohol fermentation, while grape tannins can be added before aging.

The chromosome-scale reference genome of Rubus chingii Hu provides insight into the biosynthetic pathway of hydrolyzable tannins.

Rubus chingii Hu (Fu-Pen-Zi), a perennial woody plant in the Rosaceae family, is a characteristic traditional Chinese medicinal plant because of its unique pharmacological effects. There are abundant hydrolyzable tannin (HT) components in R. chingii that provide health benefits. Here, an R. chingii chromosome-scale genome and related functional analysis provide insights into the biosynthetic pathway of HTs. In total, sequence data of 231.21 Mb (155 scaffolds with an N50 of 8.2 Mb) were assembled into seven chromosomes with an average length of 31.4 Mb, and 33 130 protein-coding genes were predicted, 89.28% of which were functionally annotated. Evolutionary analysis showed that R. chingii was most closely related to Rubus occidentalis, from which it was predicted to have diverged 22.46 million years ago (Table S8). Comparative genomic analysis showed that there was a tandem gene cluster of UGT, carboxylesterase (CXE) and SCPL genes on chromosome 02 of R. chingii, including 11 CXE, eight UGT, and six SCPL genes, which may be critical for the synthesis of HTs. In vitro enzyme assays indicated that the proteins encoded by the CXE (LG02.4273) and UGT (LG02.4102) genes have tannin hydrolase and gallic acid glycosyltransferase functions, respectively. The genomic sequence of R. chingii will be a valuable resource for comparative genomic analysis within the Rosaceae family and will be useful for understanding the biosynthesis of HTs.

A predicted protein functional network aids in novel gene mining for characteristic secondary metabolites in tea plant (Camellia sinensis).

Modeling a protein functional network in concerned species is an efficient approach for identifying novel genes in certain biological pathways. Tea plant (Camellia sinensis) is an important commercial crop abundant in numerous characteristic secondary metabolites (e.g., polyphenols, alkaloids, alkaloids) that confer tea quality and health benefits. Decoding novel genes responsible for tea characteristic components is an important basis for applied genetic improvement and metabolic engineering. Herein, a high-quality protein functional network for tea plant (TeaPoN) was predicted using cross-species protein functional associations transferring and integration combined with a stringent biological network criterion control. TeaPoN contained 31,273 nonredundant functional interactions among 6,634 tea proteins (or genes), with general network topological properties such as scale-free and small-world. We revealed the modular organization of genes related to the major three tea characteristic components (theanine, caffeine, catechin) in TeaPoN, which served as strong evidence for the utility of TeaPoN in novel gene mining. Importantly, several case studies regarding gene identification for tea characteristic components were presented. To aid in the use of TeaPoN, a concise web interface for data deposit and novel gene screening was developed (http://teapon.wchoda.com). We believe that TeaPoN will serve as a useful platform for functional genomics studies associated with characteristic secondary metabolites in tea plant.

Insights into the mechanism of cyanobacteria removal by the algicidal fungi Bjerkandera adusta and Trametes versicolor.

Fungal mycelia can eliminate almost all cocultured cyanobacterial cells within a short time. However, molecular mechanisms of algicidal fungi are poorly understood. In this study, a time-course transcriptomic analysis of algicidal fungus Bjerkandera adusta T1 was applied to investigate gene expression and regulation. A total of 132, 300, 422, and 823 differentially expressed genes (DEGs) were identified at 6, 12, 24, and 48 hr, respectively. Most DEGs exhibited high endopeptidase activity, cellulose catabolic process, and transmembrane transporter activity by using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Many decomposition genes encoding endopeptidases were induced a little later in B. adusta T1 when compared with previously investigated algicidal fungus Trametes versicolor F21a. Besides, the accumulated expression of Polysaccharide lyases8 (PL8) gene with peptidoglycan and alginate decomposition abilities was greatly delayed in B. adusta T1 relative to T. versicolor F21a. It was implied that endopeptidases and enzymes of PL8 might be responsible for the strong algicidal ability of B. adusta T1 as well as T. versicolor F21a.

Genome-wide association study leads to novel genetic insights into resistance to Aspergillus flavus in maize kernels.

BACKGROUND: Fungus infection in staple grains affects the food storage and threatens food security. The Aspergillus flavus is known to infect multiple grains and produce mycotoxin Aflatoxin B1, which is mutagenic, teratogenic and causes immunosuppression in animals. However, the molecular mechanism of maize resistance to A. flavus is largely unknown. RESULTS: Here we used corn kernels to investigate resistance genes to A. flavus using genome-wide association study (GWAS) of 313 inbred lines. We characterized the resistance levels of kernels after inoculating with A. flavus. The GWAS with 558,529 SNPs identified four associated loci involving 29 candidate genes that were linked to seed development, resistance or infection, and involved in signal pathways, seed development, germination, dormancy, epigenetic modification, and antimicrobial activity. In addition, a few candidate genes were also associated with several G-protein signaling and phytohormones that might involve in synergistic work conferring different resistance during seed development. Expression of 16 genes out of 29 during kernel development was also associated with resistance levels. CONCLUSIONS: We characterized the resistance levels of 313 maize kernels after inoculating with A. flavus, and found four associated loci and 16 candidate maize genes. The expressed 16 genes involved in kernel structure and kernel composition most likely contribute to mature maize kernels' resistance to A. flavus, and in particular, in the development of pericarp. The linked candidate genes could be experimentally transformed to validate and manipulate fungal resistance. Thus this result adds value to maize kernels in breeding programs.

Exceptional subgenome stability and functional divergence in the allotetraploid Ethiopian cereal teff.

Teff (Eragrostis tef) is a cornerstone of food security in the Horn of Africa, where it is prized for stress resilience, grain nutrition, and market value. Here, we report a chromosome-scale assembly of allotetraploid teff (variety Dabbi) and patterns of subgenome dynamics. The teff genome contains two complete sets of homoeologous chromosomes, with most genes maintaining as syntenic gene pairs. TE analysis allows us to estimate that the teff polyploidy event occurred ~1.1 million years ago (mya) and that the two subgenomes diverged ~5.0 mya. Despite this divergence, we detect no large-scale structural rearrangements, homoeologous exchanges, or biased gene loss, in contrast to many other allopolyploids. The two teff subgenomes have partitioned their ancestral functions based on divergent expression across a diverse expression atlas. Together, these genomic resources will be useful for accelerating breeding of this underutilized grain crop and for fundamental insights into polyploid genome evolution.

A Novel Site-Specific Integration System for Genetic Modification of Aspergillus flavus.

Aspergillus flavus is a fungus that produces aflatoxin B1, one of the most carcinogenic secondary metabolites. Understanding the regulation mechanism of aflatoxin biosynthesis in this fungus requires precise methods for genomic integration of mutant alleles. To avoid the disadvantage of DNA integration into the genome by non-homologous or ectopic recombination, we developed a novel strategy for site-specific integration of foreign DNA by using a carboxin-resistant sdh2R allele (His 249 Leu). Our results demonstrated that the transformants were generated with a high efficiency (>96%) of correct integration into the sdh2-lcus of the genome of A. flavus NRRL 3357. The advantage of this method is that introduction of the eGFP expression cassette into the sdh2-locus had little effect on fungal growth and virulence while also being rapid and efficient. This system will be a valuable tool for genetic manipulation in A. flavus To the best of our knowledge, this is the first report on the efficient site-specific integration at the sdh2-locus in the genome of Aspergillus.

Identification of Long Non-Coding RNAs and the Regulatory Network Responsive to Arbuscular Mycorrhizal Fungi Colonization in Maize Roots.

Recently, long noncoding RNAs (lncRNAs) have emerged as vital regulators of many biological processes in animals and plants. However, to our knowledge no investigations on plant lncRNAs which respond to arbuscular mycorrhizal (AM) fungi have been reported thus far. In this study, maize roots colonized with AM fungus were analyzed by strand-specific RNA-Seq to identify AM fungi-responsive lncRNAs and construct an associated regulatory network. A total of 1837 differentially expressed protein coding genes (DEGs) were identified from maize roots with Rhizophagus irregularis inoculation. Many AM fungi-responsive genes were homologs to MtPt4, STR, STR2, MtFatM, and enriched pathways such as fatty acid biosynthesis, response to phosphate starvation, and nitrogen metabolism are consistent with previous studies. In total, 5941 lncRNAs were identified, of which more than 3000 were new. Of those, 63 lncRNAs were differentially expressed. The putative target genes of differentially expressed lncRNAs (DELs) were mainly related to phosphate ion transmembrane transport, cellular response to potassium ion starvation, and lipid catabolic processes. Regulatory network analysis showed that DELs might be involved in the regulation of bidirectional nutrient exchange between plant and AM fungi as mimicry of microRNA targets. The results of this study can broaden our knowledge on the interaction between plant and AM fungi.