Integrating GWAS and TWAS to elucidate the genetic architecture of maize leaf cuticular conductance.

The cuticle, a hydrophobic layer of cutin and waxes synthesized by plant epidermal cells, is the major barrier to water loss when stomata are closed. Dissecting the genetic architecture of natural variation for maize (Zea mays L.) leaf cuticular conductance (gc) is important for identifying genes relevant to improving crop productivity in drought-prone environments. To this end, we performed an integrated genome- and transcriptome-wide association studies (GWAS and TWAS) to identify candidate genes putatively regulating variation in leaf gc. Of the 22 plausible candidate genes identified, 4 were predicted to be involved in cuticle precursor biosynthesis and export, 2 in cell wall modification, 9 in intracellular membrane trafficking, and 7 in the regulation of cuticle development. A gene encoding an INCREASED SALT TOLERANCE1-LIKE1 (ISTL1) protein putatively involved in intracellular protein and membrane trafficking was identified in GWAS and TWAS as the strongest candidate causal gene. A set of maize nested near-isogenic lines that harbor the ISTL1 genomic region from eight donor parents were evaluated for gc, confirming the association between gc and ISTL1 in a haplotype-based association analysis. The findings of this study provide insights into the role of regulatory variation in the development of the maize leaf cuticle and will ultimately assist breeders to develop drought-tolerant maize for target environments.

Structure-function analysis of the maize bulliform cell cuticle and its potential role in dehydration and leaf rolling.

The hydrophobic cuticle of plant shoots serves as an important interaction interface with the environment. It consists of the lipid polymer cutin, embedded with and covered by waxes, and provides protection against stresses including desiccation, UV radiation, and pathogen attack. Bulliform cells form in longitudinal strips on the adaxial leaf surface, and have been implicated in the leaf rolling response observed in drought-stressed grass leaves. In this study, we show that bulliform cells of the adult maize leaf epidermis have a specialized cuticle, and we investigate its function along with that of bulliform cells themselves. Bulliform cells displayed increased shrinkage compared to other epidermal cell types during dehydration of the leaf, providing a potential mechanism to facilitate leaf rolling. Analysis of natural variation was used to relate bulliform strip patterning to leaf rolling rate, providing further evidence of a role for bulliform cells in leaf rolling. Bulliform cell cuticles showed a distinct ultrastructure with increased cuticle thickness compared to other leaf epidermal cells. Comparisons of cuticular conductance between adaxial and abaxial leaf surfaces, and between bulliform-enriched mutants versus wild-type siblings, showed a correlation between elevated water loss rates and presence or increased density of bulliform cells, suggesting that bulliform cuticles are more water-permeable. Biochemical analysis revealed altered cutin composition and increased cutin monomer content in bulliform-enriched tissues. In particular, our findings suggest that an increase in 9,10-epoxy-18-hydroxyoctadecanoic acid content, and a lower proportion of ferulate, are characteristics of bulliform cuticles. We hypothesize that elevated water permeability of the bulliform cell cuticle contributes to the differential shrinkage of these cells during leaf dehydration, thereby facilitating the function of bulliform cells in stress-induced leaf rolling observed in grasses.

Transcriptomic network analyses shed light on the regulation of cuticle development in maize leaves.

Plant cuticles are composed of wax and cutin and evolved in the land plants as a hydrophobic boundary that reduces water loss from the plant epidermis. The expanding maize adult leaf displays a dynamic, proximodistal gradient of cuticle development, from the leaf base to the tip. Laser microdissection RNA Sequencing (LM-RNAseq) was performed along this proximodistal gradient, and complementary network analyses identified potential regulators of cuticle biosynthesis and deposition. A weighted gene coexpression network (WGCN) analysis suggested a previously undescribed function for PHYTOCHROME-mediated light signaling during the regulation of cuticular wax deposition. Genetic analyses reveal that phyB1 phyB2 double mutants of maize exhibit abnormal cuticle composition, supporting the predictions of our coexpression analysis. Reverse genetic analyses also show that phy mutants of the moss Physcomitrella patens exhibit abnormal cuticle composition, suggesting an ancestral role for PHYTOCHROME-mediated, light-stimulated regulation of cuticle development during plant evolution.

Genome-Wide Association Study for Maize Leaf Cuticular Conductance Identifies Candidate Genes Involved in the Regulation of Cuticle Development.

The cuticle, a hydrophobic layer of cutin and waxes synthesized by plant epidermal cells, is the major barrier to water loss when stomata are closed at night and under water-limited conditions. Elucidating the genetic architecture of natural variation for leaf cuticular conductance (gc) is important for identifying genes relevant to improving crop productivity in drought-prone environments. To this end, we conducted a genome-wide association study of gc of adult leaves in a maize inbred association panel that was evaluated in four environments (Maricopa, AZ, and San Diego, CA, in 2016 and 2017). Five genomic regions significantly associated with gc were resolved to seven plausible candidate genes (ISTL1, two SEC14 homologs, cyclase-associated protein, a CER7 homolog, GDSL lipase, and ß-D-XYLOSIDASE 4). These candidates are potentially involved in cuticle biosynthesis, trafficking and deposition of cuticle lipids, cutin polymerization, and cell wall modification. Laser microdissection RNA sequencing revealed that all these candidate genes, with the exception of the CER7 homolog, were expressed in the zone of the expanding adult maize leaf where cuticle maturation occurs. With direct application to genetic improvement, moderately high average predictive abilities were observed for whole-genome prediction of gc in locations (0.46 and 0.45) and across all environments (0.52). The findings of this study provide novel insights into the genetic control of gc and have the potential to help breeders more effectively develop drought-tolerant maize for target environments.

Constructing functional cuticles: analysis of relationships between cuticle lipid composition, ultrastructure and water barrier function in developing adult maize leaves.

BACKGROUND AND AIMS: Prior work has examined cuticle function, composition and ultrastructure in many plant species, but much remains to be learned about how these features are related. This study aims to elucidate relationships between these features via analysis of cuticle development in adult maize (Zea mays L.) leaves, while also providing the most comprehensive investigation to date of the composition and ultrastructure of adult leaf cuticles in this important crop plant. METHODS: We examined water permeability, wax and cutin composition via gas chromatography, and ultrastructure via transmission electron microscopy, along the developmental gradient of partially expanded adult maize leaves, and analysed the relationships between these features. KEY RESULTS: The water barrier property of the adult maize leaf cuticle is acquired at the cessation of cell expansion. Wax types and chain lengths accumulate asynchronously over the course of development, while overall wax load does not vary. Cutin begins to accumulate prior to establishment of the water barrier and continues thereafter. Ultrastructurally, pavement cell cuticles consist of an epicuticular layer, and a thin cuticle proper that acquires an inner, osmiophilic layer during development. CONCLUSIONS: Cuticular waxes of the adult maize leaf are dominated by alkanes and alkyl esters. Unexpectedly, these are localized mainly in the epicuticular layer. Establishment of the water barrier during development coincides with a switch from alkanes to esters as the major wax type, and the emergence of an osmiophilic (likely cutin-rich) layer of the cuticle proper. Thus, alkyl esters and the deposition of the cutin polyester are implicated as key components of the water barrier property of adult maize leaf cuticles.

Machine Learning Enables High-Throughput Phenotyping for Analyses of the Genetic Architecture of Bulliform Cell Patterning in Maize.

Bulliform cells comprise specialized cell types that develop on the adaxial (upper) surface of grass leaves, and are patterned to form linear rows along the proximodistal axis of the adult leaf blade. Bulliform cell patterning affects leaf angle and is presumed to function during leaf rolling, thereby reducing water loss during temperature extremes and drought. In this study, epidermal leaf impressions were collected from a genetically and anatomically diverse population of maize inbred lines. Subsequently, convolutional neural networks were employed to measure microscopic, bulliform cell-patterning phenotypes in high-throughput. A genome-wide association study, combined with RNAseq analyses of the bulliform cell ontogenic zone, identified candidate regulatory genes affecting bulliform cell column number and cell width. This study is the first to combine machine learning approaches, transcriptomics, and genomics to study bulliform cell patterning, and the first to utilize natural variation to investigate the genetic architecture of this microscopic trait. In addition, this study provides insight toward the improvement of macroscopic traits such as drought resistance and plant architecture in an agronomically important crop plant.

Dominant, Heritable Resistance to Stewart's Wilt in Maize Is Associated with an Enhanced Vascular Defense Response to Infection with Pantoea stewartii.

Vascular wilt bacteria such as Pantoea stewartii, the causal agent of Stewart's bacterial wilt of maize (SW), are destructive pathogens that are difficult to control. These bacteria colonize the xylem, where they form biofilms that block sap flow leading to characteristic wilting symptoms. Heritable forms of SW resistance exist and are used in maize breeding programs but the underlying genes and mechanisms are mostly unknown. Here, we show that seedlings of maize inbred lines with pan1 mutations are highly resistant to SW. However, current evidence suggests that other genes introgressed along with pan1 are responsible for resistance. Genomic analyses of pan1 lines were used to identify candidate resistance genes. In-depth comparison of P. stewartii interaction with susceptible and resistant maize lines revealed an enhanced vascular defense response in pan1 lines characterized by accumulation of electron-dense materials in xylem conduits visible by electron microscopy. We propose that this vascular defense response restricts P. stewartii spread through the vasculature, reducing both systemic bacterial colonization of the xylem network and consequent wilting. Though apparently unrelated to the resistance phenotype of pan1 lines, we also demonstrate that the effector WtsE is essential for P. stewartii xylem dissemination, show evidence for a nutritional immunity response to P. stewartii that alters xylem sap composition, and present the first analysis of maize transcriptional responses to P. stewartii infection.

Integration of omic networks in a developmental atlas of maize.

Coexpression networks and gene regulatory networks (GRNs) are emerging as important tools for predicting functional roles of individual genes at a system-wide scale. To enable network reconstructions, we built a large-scale gene expression atlas composed of 62,547 messenger RNAs (mRNAs), 17,862 nonmodified proteins, and 6227 phosphoproteins harboring 31,595 phosphorylation sites quantified across maize development. Networks in which nodes are genes connected on the basis of highly correlated expression patterns of mRNAs were very different from networks that were based on coexpression of proteins. Roughly 85% of highly interconnected hubs were not conserved in expression between RNA and protein networks. However, networks from either data type were enriched in similar ontological categories and were effective in predicting known regulatory relationships. Integration of mRNA, protein, and phosphoprotein data sets greatly improved the predictive power of GRNs.

A high-resolution tissue-specific proteome and phosphoproteome atlas of maize primary roots reveals functional gradients along the root axes.

A high-resolution proteome and phosphoproteome atlas of four maize (Zea mays) primary root tissues, the cortex, stele, meristematic zone, and elongation zone, was generated. High-performance liquid chromatography coupled with tandem mass spectrometry identified 11,552 distinct nonmodified and 2,852 phosphorylated proteins across the four root tissues. Two gradients reflecting the abundance of functional protein classes along the longitudinal root axis were observed. While the classes RNA, DNA, and protein peaked in the meristematic zone, cell wall, lipid metabolism, stress, transport, and secondary metabolism culminated in the differentiation zone. Functional specialization of tissues is underscored by six of 10 cortex-specific proteins involved in flavonoid biosynthesis. Comparison of this data set with high-resolution seed and leaf proteome studies revealed 13% (1,504/11,552) root-specific proteins. While only 23% of the 1,504 root-specific proteins accumulated in all four root tissues, 61% of all 11,552 identified proteins accumulated in all four root tissues. This suggests a much higher degree of tissue-specific functionalization of root-specific proteins. In summary, these data illustrate the remarkable plasticity of the proteomic landscape of maize primary roots and thus provide a starting point for gaining a better understanding of their tissue-specific functions.

The SCAR/WAVE complex polarizes PAN receptors and promotes division asymmetry in maize.

Pre-mitotic establishment of polarity is a key event in the preparation of mother cells for asymmetric cell divisions that produce daughters of distinct fates, and ensures correct cellular patterning of tissues and eventual organ function. Previous work has shown that two receptor-like kinases, PANGLOSS2 (PAN2) and PAN1, and the small GTPase RHO GTPASE OF PLANTS (ROP) promote mother cell polarity and subsequent division asymmetry in developing maize stomata. PAN proteins become polarized prior to asymmetric cell division, however, the mechanism of this polarization is unknown. Here we show that the SCAR/WAVE regulatory complex, which activates the actin-nucleating ARP2/3 complex, is the first known marker of polarity in this asymmetric division model and is required for PAN polarization. These findings implicate actin, and specifically branched actin networks, in PAN polarization and asymmetric cell division.

Unraveling genomic complexity at a quantitative disease resistance locus in maize.

Multiple disease resistance has important implications for plant fitness, given the selection pressure that many pathogens exert directly on natural plant populations and indirectly via crop improvement programs. Evidence of a locus conditioning resistance to multiple pathogens was found in bin 1.06 of the maize genome with the allele from inbred line "Tx303" conditioning quantitative resistance to northern leaf blight (NLB) and qualitative resistance to Stewart's wilt. To dissect the genetic basis of resistance in this region and to refine candidate gene hypotheses, we mapped resistance to the two diseases. Both resistance phenotypes were localized to overlapping regions, with the Stewart's wilt interval refined to a 95.9-kb segment containing three genes and the NLB interval to a 3.60-Mb segment containing 117 genes. Regions of the introgression showed little to no recombination, suggesting structural differences between the inbred lines Tx303 and "B73," the parents of the fine-mapping population. We examined copy number variation across the region using next-generation sequencing data, and found large variation in read depth in Tx303 across the region relative to the reference genome of B73. In the fine-mapping region, association mapping for NLB implicated candidate genes, including a putative zinc finger and pan1. We tested mutant alleles and found that pan1 is a susceptibility gene for NLB and Stewart's wilt. Our data strongly suggest that structural variation plays an important role in resistance conditioned by this region, and pan1, a gene conditioning susceptibility for NLB, may underlie the QTL.

Divergent roles for maize PAN1 and PAN2 receptor-like proteins in cytokinesis and cell morphogenesis.

Pangloss1 (PAN1) and PAN2 are leucine-rich repeat receptor-like proteins that function cooperatively to polarize the divisions of subsidiary mother cells (SMCs) during stomatal development in maize (Zea mays). PANs colocalize in SMCs, and both PAN1 and PAN2 promote polarization of the actin cytoskeleton and nuclei in these cells. Here, we show that PAN1 and PAN2 have additional functions that are unequal or divergent. PAN1, but not PAN2, is localized to cell plates in all classes of dividing cells examined. pan1 mutants exhibited no defects in cell plate formation or in the recruitment or removal of a variety of cell plate components; thus, they did not demonstrate a function for PAN1 in cytokinesis. PAN2, in turn, plays a greater role than PAN1 in directing patterns of postmitotic cell expansion that determine the shapes of mature stomatal subsidiary cells and interstomatal cells. Localization studies indicate that PAN2 impacts subsidiary cell shape indirectly by stimulating localized cortical actin accumulation and polarized growth in interstomatal cells. Localization of PAN1, Rho of Plants2, and PIN1a suggests that PAN2-dependent cell shape changes do not involve any of these proteins, indicating that PAN2 function is linked to actin polymerization by a different mechanism in interstomatal cells compared with SMCs. Together, these results demonstrate that PAN1 and PAN2 are not dedicated to SMC polarization but instead play broader roles in plant development. We speculate that PANs may function in all contexts to regulate polarized membrane trafficking either directly or indirectly via their influence on actin polymerization.

Reconstruction of protein networks from an atlas of maize seed proteotypes.

A comprehensive knowledge of proteomic states is essential for understanding biological systems. Using mass spectrometry, we mapped an atlas of developing maize seed proteotypes comprising 14,165 proteins and 18,405 phosphopeptides (from 4,511 proteins), quantified across eight tissues. We found that many of the most abundant proteins are not associated with detectable levels of their mRNAs, and we provide evidence for three potential explanations: transport of proteins between tissues; diurnal, out-of-phase accumulation of mRNAs and cognate proteins; and differential lifetimes of mRNAs compared with proteins. Likewise, many of the most abundant mRNAs were not associated with detectable levels of their proteins. Across the entire dataset, protein abundance was poorly correlated with mRNA levels and was largely independent of phosphorylation status. Comparisons between proteotypes revealed the quantitative contribution of specific proteins and phosphorylation events to the spatially and temporally regulated starch and oil biosynthetic pathways. Reconstruction of signaling networks established associations of proteins and phosphoproteins with distinct biological processes acting during seed development. Additionally, a protein kinase substrate network was reconstructed, enabling the identification of 762 potential substrates of specific protein kinases. Finally, examination of 694 transcription factors revealed remarkable constraints on patterns of expression and phosphorylation within transcription factor families. These results provide a resource for understanding seed development in a crop that is the foundation of modern agriculture.

Parallel proteomic and phosphoproteomic analyses of successive stages of maize leaf development.

We performed large-scale, quantitative analyses of the maize (Zea mays) leaf proteome and phosphoproteome at four developmental stages. Exploiting the developmental gradient of maize leaves, we analyzed protein and phosphoprotein abundance as maize leaves transition from proliferative cell division to differentiation to cell expansion and compared these developing zones to one another and the mature leaf blade. Comparison of the proteomes and phosphoproteomes suggests a key role for posttranslational regulation in developmental transitions. Analysis of proteins with cell wall- and hormone-related functions illustrates the utility of the data set and provides further insight into maize leaf development. We compare phosphorylation sites identified here to those previously identified in Arabidopsis thaliana. We also discuss instances where comparison of phosphorylated and unmodified peptides from a particular protein indicates tissue-specific phosphorylation. For example, comparison of unmodified and phosphorylated forms of PINFORMED1 (PIN1) suggests a tissue-specific difference in phosphorylation, which correlates with changes in PIN1 polarization in epidermal cells during development. Together, our data provide insights into regulatory processes underlying maize leaf development and provide a community resource cataloging the abundance and phosphorylation status of thousands of maize proteins at four leaf developmental stages.

Identification of PAN2 by quantitative proteomics as a leucine-rich repeat-receptor-like kinase acting upstream of PAN1 to polarize cell division in maize.

Mechanisms governing the polarization of plant cell division are poorly understood. Previously, we identified pangloss1 (PAN1) as a leucine-rich repeat-receptor-like kinase (LRR-RLK) that promotes the polarization of subsidiary mother cell (SMC) divisions toward the adjacent guard mother cell (GMC) during stomatal development in maize (Zea mays). Here, we identify pangloss2 (PAN2) as a second LRR-RLK promoting SMC polarization. Quantitative proteomic analysis identified a PAN2 candidate by its depletion from membranes of pan2 single and pan1;pan2 double mutants. Genetic mapping and sequencing of mutant alleles confirmed the identity of this protein as PAN2. Like PAN1, PAN2 has a catalytically inactive kinase domain and accumulates in SMCs at sites of GMC contact before nuclear polarization. The timing of polarized PAN1 and PAN2 localization is very similar, but PAN2 acts upstream because it is required for polarized accumulation of PAN1 but is independent of PAN1 for its own localization. We find no evidence that PAN2 recruits PAN1 to the GMC contact site via a direct or indirect physical interaction, but PAN2 interacts with itself. Together, these results place PAN2 at the top of a cascade of events promoting the polarization of SMC divisions, potentially functioning to perceive or amplify GMC-derived polarizing cues.

Division polarity in developing stomata.

Stomata are generated via asymmetric cell division in both dicots and monocots. Intrinsic or extrinsic polarity cues are perceived and acted upon to generate mother cell polarity and determine asymmetric division planes. Arabidopsis employs both intrinsic and extrinsic cues to orient a variable series of asymmetric stomatal divisions, using novel proteins such as BASL and POLAR to generate polarity. In contrast, maize appears to employ only extrinsic cues to orient the polarities of divisions occurring in an invariant sequence to generate stomatal complexes. Although both plants use receptor-like kinases to generate or orient division polarity in developing stomata, there are few similarities in the proteins and pathway identified to date as regulators of these processes.

ROP GTPases act with the receptor-like protein PAN1 to polarize asymmetric cell division in maize.

Plant Rho family GTPases (ROPs) have been investigated primarily for their functions in polarized cell growth. We previously showed that the maize (Zea mays) Leu-rich repeat receptor-like protein PANGLOSS1 (PAN1) promotes the polarization of asymmetric subsidiary mother cell (SMC) divisions during stomatal development. Here, we show that maize Type I ROPs 2 and 9 function together with PAN1 in this process. Partial loss of ROP2/9 function causes a weak SMC division polarity phenotype and strongly enhances this phenotype in pan1 mutants. Like PAN1, ROPs accumulate in an asymmetric manner in SMCs. Overexpression of yellow fluorescent protein-ROP2 is associated with its delocalization in SMCs and with aberrantly oriented SMC divisions. Polarized localization of ROPs depends on PAN1, but PAN1 localization is insensitive to depletion and depolarization of ROP. Membrane-associated Type I ROPs display increased nonionic detergent solubility in pan1 mutants, suggesting a role for PAN1 in membrane partitioning of ROPs. Finally, endogenous PAN1 and ROP proteins are physically associated with each other in maize tissue extracts, as demonstrated by reciprocal coimmunoprecipitation experiments. This study demonstrates that ROPs play a key role in polarization of plant cell division and cell growth and reveals a role for a receptor-like protein in spatial localization of ROPs.

Determination of symmetric and asymmetric division planes in plant cells.

The cellular organization of plant tissues is determined by patterns of cell division and growth coupled with cellular differentiation. Cells proliferate mainly via symmetric division, whereas asymmetric divisions are associated with initiation of new developmental patterns and cell types. Division planes in both symmetrically and asymmetrically dividing cells are established through the action of a cortical preprophase band (PPB) of cytoskeletal filaments, which is disassembled upon transition to metaphase, leaving behind a cortical division site (CDS) to which the cytokinetic phragmoplast is later guided to position the cell plate. Recent progress has been made in understanding PPB formation and function as well as the nature and function of the CDS. In asymmetrically dividing cells, division plane establishment is governed by cell polarity. Recent work is beginning to shed light on polarization mechanisms in asymmetrically dividing cells, with receptor-like proteins and potential downstream effectors emerging as important players in this process.

Tangled localization at the cortical division site of plant cells occurs by several mechanisms.

TANGLED (TAN) is the founding member of a family of plant-specific proteins required for correct orientation of the division plane. Arabidopsis thaliana TAN is localized before prophase until the end of cytokinesis at the cortical division site (CDS), where it appears to help guide the cytokinetic apparatus towards the cortex. We show that TAN is actively recruited to the CDS by distinct mechanisms before and after preprophase band (PPB) disassembly. Colocalization with the PPB is mediated by one region of TAN, whereas another region mediates its recruitment to the CDS during cytokinesis. This second region binds directly to POK1, a kinesin that is required for TAN localization. Although this region of TAN is recruited to the CDS during cytokinesis without first colocalizing with the PPB, pharmacological evidence indicates that the PPB is nevertheless required for both early and late localization of TAN at the CDS. Finally, we show that phosphatase activity is required for maintenance of early but not late TAN localization at the CDS. We propose a new model in which TAN is actively recruited to the CDS by several mechanisms, indicating that the CDS is dynamically modified from prophase through to the completion of cytokinesis.

Visualization of F-actin localization and dynamics with live cell markers in Neurospora crassa.

Filamentous actin (F-actin) plays essential roles in filamentous fungi, as in all other eukaryotes, in a wide variety of cellular processes including cell growth, intracellular motility, and cytokinesis. We visualized F-actin organization and dynamics in living Neurospora crassa cells via confocal microscopy of growing hyphae expressing GFP fusions with homologues of the actin-binding proteins fimbrin (FIM) and tropomyosin (TPM-1), a subunit of the Arp2/3 complex (ARP-3) and a recently developed live cell F-actin marker, Lifeact (ABP140 of Saccharomyces cerevisiae). FIM-GFP, ARP-3-GFP, and Lifeact-GFP associated with small patches in the cortical cytoplasm that were concentrated in a subapical ring, which appeared similar for all three markers but was broadest in hyphae expressing Lifeact-GFP. These cortical patches were short-lived, and a subset was mobile throughout the hypha, exhibiting both anterograde and retrograde motility. TPM-1-GFP and Lifeact-GFP co-localized within the Spitzenkörper (Spk) core at the hyphal apex, and were also observed in actin cables throughout the hypha. All GFP fusion proteins studied were also transiently localized at septa: Lifeact-GFP first appeared as a broad ring during early stages of contractile ring formation and later coalesced into a sharper ring, TPM-1-GFP was observed in maturing septa, and FIM-GFP/ARP3-GFP-labeled cortical patches formed a double ring flanking the septa. Our observations suggest that each of the N. crassa F-actin-binding proteins analyzed associates with a different subset of F-actin structures, presumably reflecting distinct roles in F-actin organization and dynamics. Moreover, Lifeact-GFP marked the broadest spectrum of F-actin structures; it may serve as a global live cell marker for F-actin in filamentous fungi.