Integration of Metabolomes and Transcriptomes Provides Insights into Morphogenesis and Maturation in Morchella sextelata.

True morels (Morchella, Pezizales) are a popular edible and medicinal fungus with great nutritional and economic value. The dynamics and regulatory mechanisms during the morphogenesis and maturation of morels are poorly understood. In this study, the metabolomes and transcriptomes of the mycelium (MY), primordium differentiation (PR), young fruiting body (YFB), and mature fruiting body (MFB) were comprehensively analyzed to reveal the mechanism of the morphogenesis and maturation of Morchella sextelata. A total of 748 differentially expressed metabolites (DEMs) and 5342 differentially expressed genes (DEGs) were detected, mainly enriched in the carbohydrate, amino acid, and lipid metabolism pathways, with the transition from the mycelium to the primordium being the most drastic stage at both the metabolic and transcriptional levels. The integrated metabolomics and transcriptomics highlighted significant correlations between the DEMs and DEGs, and specific amino acid and nucleotide metabolic pathways were significantly co-enriched, which may play key roles in morphological development and ascocarp maturation. A conceptual model of transcriptional and metabolic regulation was proposed during morphogenesis and maturation in M. sextelata for the first time, in which environmental factors activate the regulation of transcription factors, which then promote metabolic and transcriptional regulation from vegetative to reproductive growth. These results provide insights into the metabolic dynamics and transcriptional regulation during the morphogenesis and maturation of morels and valuable resources for future breeding enhancement and sustainable artificial cultivation.

Dynamics of soil microbiome throughout the cultivation life cycle of morel (Morchella sextelata).

Although Morchella sextelata (morel) is a well-known, edible, and medicinal fungus widely cultivated in China, the dynamics and roles of its soil microbiome during cultivation are unclear. Using rhizosphere soil samples collected throughout the M. sextelata cultivation life cycle, we conducted a high-throughput metagenomic sequencing analysis, with an emphasis on variations in soil microbial composition, characteristic biomarkers, and ecological functions. We found that microbial relative abundance, alpha diversity, and structure varied significantly among fungal growth stages. A total of 47 stage-associated biomarkers were identified through a linear discriminant analysis of effect size. In addition, horizontal comparison of soil microbiomes exhibiting successful and failed primordium formation further confirmed primordium-associated microbes with possible key roles in primordium formation. A microbial function analysis revealed that nutrient metabolism-related pathways were enriched during mycelium and fruiting body stages, whereas the signal transduction pathway was enriched during the primordium stage. This result indicates that diverse microbes are required at different growth stages of M. sextelata. Our research has revealed the dynamic scenario of the soil microbiome throughout the cultivation life cycle of M. sextelata. The high-resolution microbial profiles uncovered in the present study provide novel insights that should contribute to the improvement of morel cultivation using microbial inoculants.

The Elite Alleles of OsSPL4 Regulate Grain Size and Increase Grain Yield in Rice.

Grain weight and grain number, the two important yield traits, are mainly determined by grain size and panicle architecture in rice. Herein, we report the identification and functional analysis of OsSPL4 in panicle and grain development of rice. Using CRISPR/Cas9 system, two elite alleles of OsSPL4 were obtained, which exhibited an increasing number of grains per panicle and grain size, resulting in increase of rice yield. Cytological analysis showed that OsSPL4 could regulate spikelet development by promoting cell division. The results of RNA-seq and qRT-PCR validations also demonstrated that several MADS-box and cell-cycle genes were up-regulated in the mutation lines. Co-expression network revealed that many yield-related genes were involved in the regulation network of OsSPL4. In addition, OsSPL4 could be cleaved by the osa-miR156 in vivo, and the OsmiR156-OsSPL4 module might regulate the grain size in rice. Further analysis indicated that the large-grain allele of OsSPL4 in indica rice might introgress from aus varieties under artificial selection. Taken together, our findings suggested that OsSPL4 could be as a key regulator of grain size by acting on cell division control and provided a strategy for panicle architecture and grain size modification for yield improvement in rice.

Genome Assembly and Analyses of the Macrofungus Macrocybe gigantea.

Macrocybe gigantea (M. gigantea) is a macrofungus genus that contains a big number of fairly fleshy gilled mushrooms with white spores. This macrofungus produces diverse bioactive compounds, antioxidants, and water-soluble polysaccharides. However, the genomic resources of this species remain unknown. Here, we assembled the genome of M. gigantea (41.23 Mb) into 336 scaffolds with a N50 size of 374,455 bp and compared it with the genomes of eleven other macrofungi. Comparative genomics study confirmed that M. gigantea belonged to the Macrocybe genus, a stand-alone genus different from the Tricholoma genus. In addition, we found that glycosyl hydrolase family 28 (GH28) in M. gigantea shared conserved motifs that were significantly different from their counterparts in Tricholoma. The genomic resource uncovered by this study will enhance our understanding of fungi biology, especially the differences in their growth rates and energy metabolism.

Identification of Key Genes for the Ultrahigh Yield of Rice Using Dynamic Cross-tissue Network Analysis.

Significantly increasing crop yield is a major and worldwide challenge for food supply and security. It is well-known that rice cultivated at Taoyuan in Yunnan of China can produce the highest yield worldwide. Yet, the gene regulatory mechanism underpinning this ultrahigh yield has been a mystery. Here, we systematically collected the transcriptome data for seven key tissues at different developmental stages using rice cultivated both at Taoyuan as the case group and at another regular rice planting place Jinghong as the control group. We identified the top 24 candidate high-yield genes with their network modules from these well-designed datasets by developing a novel computational systems biology method, i.e., dynamic cross-tissue (DCT) network analysis. We used one of the candidate genes, OsSPL4, whose function was previously unknown, for gene editing experimental validation of the high yield, and confirmed that OsSPL4 significantly affects panicle branching and increases the rice yield. This study, which included extensive field phenotyping, cross-tissue systems biology analyses, and functional validation, uncovered the key genes and gene regulatory networks underpinning the ultrahigh yield of rice. The DCT method could be applied to other plant or animal systems if different phenotypes under various environments with the common genome sequences of the examined sample. DCT can be downloaded from https://github.com/ztpub/DCT.

Evolution of Plant Architecture in Oryza Driven by the PROG1 Locus.

The genetic control of plant architecture in crops is critical for agriculture and understanding morphological evolution. This study showed that an open reading frame (ORF) of the rice domestication gene PROG1 appeared 3.4-3.9 million years ago (Mya). Subsequently, it acquired a novel protein-coding gene function in the genome of O. rufipogon (~0.3-0.4 Mya). This extremely young gene and its paralogous C2H2 genes located nearby define the prostrate architecture of O. rufipogon and, thus, are of adaptive significance for wild rice in swamp and water areas. However, selection for dense planting and high yield during rice domestication silenced the PROG1 gene and caused the loss of the RPAD locus containing functional C2H2 paralogs; hence, domesticated lines exhibit an erect plant architecture. Analysis of the stepwise origination process of PROG1 and its evolutionary genetics revealed that this zinc-finger coding gene may have rapidly evolved under positive selection and promoted the transition from non- or semi-prostrate growth to prostrate growth. A transgenic assay showed that PROG1 from O. rufipogon exerts a stronger function compared with PROG1 sequences from other Oryza species. However, the analysis of the expression levels of PROG1 in different Oryza species suggests that the transcriptional regulation of PROG1 has played an important role in its evolution. This study provides the first strong case showing how a fundamental morphological trait evolved in Oryza species driven by a gene locus.

Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa.

Artificially improving traits of cultivated alfalfa (Medicago sativa L.), one of the most important forage crops, is challenging due to the lack of a reference genome and an efficient genome editing protocol, which mainly result from its autotetraploidy and self-incompatibility. Here, we generate an allele-aware chromosome-level genome assembly for the cultivated alfalfa consisting of 32 allelic chromosomes by integrating high-fidelity single-molecule sequencing and Hi-C data. We further establish an efficient CRISPR/Cas9-based genome editing protocol on the basis of this genome assembly and precisely introduce tetra-allelic mutations into null mutants that display obvious phenotype changes. The mutated alleles and phenotypes of null mutants can be stably inherited in generations in a transgene-free manner by cross pollination, which may help in bypassing the debate about transgenic plants. The presented genome and CRISPR/Cas9-based transgene-free genome editing protocol provide key foundations for accelerating research and molecular breeding of this important forage crop.

Neo-functionalization of a Teosinte branched 1 homologue mediates adaptations of upland rice.

The rice orthologue of maize domestication gene Teosinte branched 1 (Tb1) affects tillering. But, unlike maize Tb1 gene, it was not selected during domestication. Here, we report that an OsTb1 duplicate gene (OsTb2) has been artificially selected during upland rice adaptation and that natural variation in OsTb2 is associated with tiller number. Interestingly, transgenic rice overexpressing this gene shows increased rather than decreased tillering, suggesting that OsTb2 gains a regulatory effect opposite to that of OsTb1 following duplication. Functional analyses suggest that the OsTb2 protein positively regulates tillering by interacting with the homologous OsTb1 protein and counteracts the inhibitory effect of OsTb1 on tillering. We further characterize two functional variations within OsTb2 that regulate protein function and gene expression, respectively. These results not only present an example of neo-functionalization that generates an opposite function following duplication but also suggest that the Tb1 homologue has been selected in upland rice.

Application of protoplast technology to CRISPR/Cas9 mutagenesis: from single-cell mutation detection to mutant plant regeneration.

Plant protoplasts are useful for assessing the efficiency of clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) mutagenesis. We improved the process of protoplast isolation and transfection of several plant species. We also developed a method to isolate and regenerate single mutagenized Nicotianna tabacum protoplasts into mature plants. Following transfection of protoplasts with constructs encoding Cas9 and sgRNAs, target gene DNA could be amplified for further analysis to determine mutagenesis efficiency. We investigated N. tabacum protoplasts and derived regenerated plants for targeted mutagenesis of the phytoene desaturase (NtPDS) gene. Genotyping of albino regenerants indicated that all four NtPDS alleles were mutated in amphidiploid tobacco, and no Cas9 DNA could be detected in most regenerated plants.

Corrigendum: Building a Genetic Manipulation Tool Box for Orchid Biology: Identification of Constitutive Promoters and Application of CRISPR/Cas9 in the Orchid, Dendrobium officinale.

[This corrects the article on p. 2036 in vol. 7, PMID: 28127299.].

Genetic Variation of Goat Interferon Regulatory Factor 3 Gene and Its Implication in Goat Evolution.

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats' populations. Fu and Li tests were significantly positive but Tajima's D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model dN/dS ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω0 (dN/dS) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.

Building a Genetic Manipulation Tool Box for Orchid Biology: Identification of Constitutive Promoters and Application of CRISPR/Cas9 in the Orchid, Dendrobium officinale.

Orchidaceae is the second largest family of flowering plants, which is highly valued for its ornamental purposes and medicinal uses. Dendrobium officinale is a special orchid species that can grow without seed vernalization. Because the whole-genome sequence of D. officinale is publicly available, this species is poised to become a convenient research model for the evolutionary, developmental, and genetic studies of Orchidaceae. Despite these advantages, the methods of genetic manipulation are poorly developed in D. officinale. In this study, based on the previously developed Agrobacterium-mediated gene transformation system, we identified several highly efficient promoters for exogenous gene expression and successfully applied the CRISPR/Cas9 system for editing endogenous genes in the genome of D. officinale. These two basic techniques contribute to the genetic manipulation toolbox of Orchidaceae. The pCambia-1301-35SN vector containing the CaMV 35S promoter and the ß-glucuronidase (GUS) and Superfolder green fluorescence protein (SG) as reporter genes were introduced into the plant tissues by the Agrobacterium-mediated transformation system. Fluorescence emission from the transformed plants confirmed the successful transcription and translation of SG genes into functional proteins. We compared the GUS activity under different promoters including four commonly used promoters (MtHP, CVMV, MMV and PCISV) with CaMV 35S promoter and found that MMV, CVMV, and PCISV were as effective as the 35S promoter. Furthermore, we applied the CRISPR/Cas9-mediated genome editing system successfully in D. officinale. By selecting five target genes (C3H, C4H, 4CL, CCR, and IRX) in the lignocellulose biosynthesis pathway, we showed that, for a given target, this system can generate edits (insertions, deletions, or substitutions) at a rate of 10 to 100%. These results showed that our two genetic manipulation tools can efficiently express exogenous genes and edit endogenous genes in D. officinale. These efficient research tools will not only help create novel D. officinale varieties, but will also facilitate the molecular genetic investigation of orchid biology.

Analysis of elite variety tag SNPs reveals an important allele in upland rice.

Elite crop varieties usually fix alleles that occur at low frequencies within non-elite gene pools. Dissecting these alleles for desirable agronomic traits can be accomplished by comparing the genomes of elite varieties with those from non-elite populations. Here we deep-sequence six elite rice varieties and use two large control panels to identify elite variety tag single-nucleotide polymorphism alleles (ETASs). Guided by this preliminary analysis, we comprehensively characterize one protein-altering ETAS in the 9-cis-epoxycarotenoid dioxygenase gene of the IRAT104 upland rice variety. This allele displays a drastic frequency difference between upland and irrigated rice, and a selective sweep is observed around this allele. Functional analysis indicates that in upland rice, this allele is associated with significantly higher abscisic acid levels and denser lateral roots, suggesting its association with upland rice suitability. This report provides a potential strategy to mine rare, agronomically important alleles.