An improved genome release (version Mt4.0) for the model legume Medicago truncatula

The Chrysosplenium sinicum genome provides insights into adaptive evolution of shade plants

Chrysosplenium sinicum, a traditional Tibetan medicinal plant, can successfully thrive in low-light environments for long periods of time. To investigate the adaptive evolution of shade plants in low-light environments, we generated a chromosome-scale genome assembly (~320 Mb) for C. sinicum by combining PacBio sequencing and Hi-C technologies. Based on our results, gene families related to photosynthesis and cell respiration greatly expanded and evolved in C. sinicum genome due to intracellular DNA transfer from organelle genome to nuclear genome. Under positive selective pressure, adaptive evolution of light-harvesting complex II (LHCII) component protein CsLhcb1s resulted in the expansion of threonine residues at the phosphorylation site of STN7 kinase, potentially establishing a crucial genomic foundation for enhancing C. sinicum's adaptability in low-light environments. Through transcriptome and metabolome analysis, we identified chrysosplenol and chrysosplenoside as predominant flavonoid metabolites of C. sinicum and predicted their synthesis pathways. In addition, analysis of alternative splicing (AS) revealed that AS events help regulate state transition and flavonoid biosynthesis. The present study provides new insights into the genomes of shade plants exposed to low-light conditions and adaptive evolution of these genomes; in addition, the results improve our current knowledge on the biosynthetic and regulatory processes of chrysosplenol and chrysosplenoside.

Species-specific microsymbiont discrimination mediated by a Medicago receptor kinase

Host range specificity is a prominent feature of the legume-rhizobial symbiosis. Sinorhizobium meliloti and Sinorhizobium medicae are two closely related species that engage in root nodule symbiosis with legume plants of the Medicago genus, but certain Medicago species exhibit selectivity in their interactions with the two rhizobial species. We have identified a Medicago receptor-like kinase, which can discriminate between the two bacterial species, acting as a genetic barrier against infection by most S. medicae strains. Activation of this receptor-mediated nodulation restriction requires a bacterial gene that encodes a glycine-rich octapeptide repeat protein with distinct variants capable of distinguishing S. medicae from S. meliloti. This study sheds light on the coevolution of host plants and rhizobia, shaping symbiotic selectivity in their respective ecological niches.

afila, the origin and nature of a major innovation in the history of pea breeding

The afila (af) mutation causes the replacement of leaflets by a branched mass of tendrils in the compound leaves of pea - Pisum sativum L. This mutation was first described in 1953, and several reports of spontaneous af mutations and induced mutants with a similar phenotype exist. Despite widespread introgression into breeding material, the nature of af and the origin of the alleles used remain unknown. Here, we combine comparative genomics with reverse genetic approaches to elucidate the genetic determinants of af. We also investigate haplotype diversity using a set of AfAf and afaf cultivars and breeding lines and molecular markers linked to seven consecutive genes. Our results show that deletion of two tandemly arranged genes encoding Q-type Cys(2)His(2) zinc finger transcription factors, PsPALM1a and PsPALM1b, is responsible for the af phenotype in pea. Eight haplotypes were identified in the af-harbouring genomic region on chromosome 2. These haplotypes differ in the size of the deletion, covering more or less genes. Diversity at the af locus is valuable for crop improvement and sheds light on the history of pea breeding for improved standing ability. The results will be used to understand the function of PsPALM1a/b and to transfer the knowledge for innovation in related crops.

Haplotype-resolved genome assembly for tetraploid Chinese cherry (Prunus pseudocerasus) offers insights into fruit firmness

Chinese cherry (Prunus pseudocerasus) holds considerable importance as one of the primary stone fruit crops in China. However, artificially improving its traits and genetic analysis are challenging due to lack of high-quality genomic resources, which mainly result from difficulties associated with resolving its tetraploid and highly heterozygous genome. Herein, we assembled a chromosome-level, haplotype-resolved genome of the cultivar 'Zhuji Duanbing', comprising 993.69 Mb assembled into 32 pseudochromosomes using PacBio HiFi, Oxford Nanopore, and Hi-C. Intra-haplotype comparative analyses revealed extensive intra-genomic sequence and expression consistency. Phylogenetic and comparative genomic analyses demonstrated that P. pseudocerasus was a stable autotetraploid species, closely related to wild P. pusilliflora, with the two diverging ~18.34 million years ago. Similar to other Prunus species, P. pseudocerasus underwent a common whole-genome duplication event that occurred ~139.96 million years ago. Because of its low fruit firmness, P. pseudocerasus is unsuitable for long-distance transportation, thereby restricting its rapid development throughout China. At the ripe fruit stage, P. pseudocerasus cv. 'Zhuji Duanbing' was significantly less firm than P. avium cv. 'Heizhenzhu'. The difference in firmness is attributed to the degree of alteration in pectin, cellulose, and hemicellulose contents. In addition, comparative transcriptomic analyses identified GalAK-like and Stv1, two genes involved in pectin biosynthesis, which potentially caused the difference in firmness between 'Zhuji Duanbing' and 'Heizhenzhu'. Transient transformations of PpsGalAK-like and PpsStv1 increase protopectin content and thereby enhance fruit firmness. Our study lays a solid foundation for functional genomic studies and the enhancement of important horticultural traits in Chinese cherries.

Evolutionary assembly of the plant terrestrialization toolkit from protein domains

Land plants (embryophytes) came about in a momentous evolutionary singularity: plant terrestrialization. This event marks not only the conquest of land by plants but also the massive radiation of embryophytes into a diverse array of novel forms and functions. The unique suite of traits present in the earliest land plants is thought to have been ushered in by a burst in genomic novelty. Here, we asked the question of how these bursts were possible. For this, we explored: (i) the initial emergence and (ii) the reshuffling of domains to give rise to hallmark environmental response genes of land plants. We pinpoint that a quarter of the embryophytic genes for stress physiology are specific to the lineage, yet a significant portion of this novelty arises not de novo but from reshuffling and recombining of pre-existing domains. Our data suggest that novel combinations of old genomic substrate shaped the plant terrestrialization toolkit, including hallmark processes in signalling, biotic interactions and specialized metabolism.

Age-related compositional and functional changes in the adult and breastfed buffalo rumen microbiome

Introduction

The buffalo is an important domestic animal globally, providing milk, meat, and labor to more than 2 billion people in 67 countries. The rumen microorganisms of buffaloes play an indispensable role in enabling the healthy functionality and digestive function of buffalo organisms. Currently, there is a lack of clarity regarding the differences in the composition and function of rumen microorganisms among buffaloes at different growth stages.

Methods

In this study, metagenomics sequencing technology was applied to examine the compositional and functional differences of rumen microorganisms in adult and breastfed buffaloes.

Results

The results revealed that the rumen of adult buffaloes had significantly higher levels of the following dominant genera: Prevotella, UBA1711, RF16, Saccharofermentans, F23-D06, UBA1777, RUG472, and Methanobrevibacter_A. Interestingly, the dominant genera specific to the rumen of adult buffaloes showed a significant positive correlation (correlation>0.5, p-value<0.05) with both lignocellulose degradation-related carbohydrate-active enzymes (CAZymes) and immune signaling pathways activated by antigenic stimulation. The rumen of breastfed buffaloes had significantly higher levels of the following dominant genera: UBA629, CAG- 791, Selenomonas_C, Treponema_D, Succinivibrio, and RC9. Simultaneously, the rumen-dominant genera specific to breastfed buffaloes were significantly positively correlated (correlation>0.5, p-value<0.05) with CAZymes associated with lactose degradation, amino acid synthesis pathways, and antibiotic-producing pathways.

Discussion

This indicates that rumen microorganisms in adult buffaloes are more engaged in lignocellulose degradation, whereas rumen microorganisms in breastfed buffaloes are more involved in lactose and amino acid degradation, as well as antibiotic production. In conclusion, these findings suggest a close relationship between differences in rumen microbes and the survival needs of buffaloes at different growth stages.

Phaseolus vulgaris STP13.1 is an H+-coupled monosaccharide transporter, present in source leaves and seed coats, with higher substrate affinity at depolarized potentials

Sugar transport proteins (STPs) are high-affinity H+-coupled hexose symporters. Recently, the contribution of STP13 to bacterial and fungal pathogen resistance across multiple plant species has garnered significant interest. Quantitative PCR analysis of source leaves, developing embryos, and seed coats of Phaseolus vulgaris L. (common bean) revealed that PvSTP13.1 was expressed in source leaves and seed coats throughout seed development. In contrast, PvSTP13.1 transcripts were detected at exceedingly low levels in developing embryos. To characterize the transport mechanism, PvSTP13.1 was expressed in Xenopus laevis oocytes, and inward-directed currents were analyzed using two-electrode voltage clamping. PvSTP13.1 was shown to function as an H+-coupled monosaccharide symporter exhibiting a unique high affinity for hexoses and aldopentoses at depolarized membrane potentials. Specifically, of the 31 assessed substrates, which included aldohexoses, deoxyhexoses, fructose, 3-O-methyl-D-glucose, aldopentoses, polyols, glycosides, disaccharides, trisaccharides, and glucuronic acid, PvSTP13.1 displayed the highest affinity (K 0.5) for glucose (43 μM), mannose (92 μM), galactose (145 μM), fructose (224 μM), xylose (1.0 mM), and fucose (3.7 mM) at pH 5.6 at a depolarized membrane potential of -40 mV. The results presented here suggest PvSTP13.1 contributes to retrieval of hexoses from the apoplasmic space in source leaves and coats of developing seeds.

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

Haplotype-resolved genome assembly of Phanera championii reveals molecular mechanisms of flavonoid synthesis and adaptive evolution

Phanera championii is a medicinal liana plant that has successfully adapted to hostile karst habitats. Despite extensive research on its medicinal components and pharmacological effects, the molecular mechanisms underlying the biosynthesis of critical flavonoids and its adaptation to karst habitats remain elusive. In this study, we performed high-coverage PacBio and Hi-C sequencing of P. championii, which revealed its high heterozygosity and phased the genome into two haplotypes: Hap1 (384.60 Mb) and Hap2 (383.70 Mb), encompassing a total of 58 612 annotated genes. Comparative genomes analysis revealed that P. championii experienced two whole-genome duplications (WGDs), with approximately 59.59% of genes originating from WGD events, thereby providing a valuable genetic resource for P. championii. Moreover, we identified a total of 112 genes that were strongly positively selected. Additionally, about 81.60 Mb of structural variations between the two haplotypes. The allele-specific expression patterns suggested that the dominant effect of P. championii was the elimination of deleterious mutations and the promotion of beneficial mutations to enhance fitness. Moreover, our transcriptome and metabolome analysis revealed alleles in different tissues or different haplotypes collectively regulate the synthesis of flavonoid metabolites. In summary, our comprehensive study highlights the significance of genomic and morphological adaptation in the successful adaptation of P. championii to karst habitats. The high-quality phased genomes obtained in this study serve as invaluable genomic resources for various applications, including germplasm conservation, breeding, evolutionary studies, and elucidation of pathways governing key biological traits of P. championii.

Chromosome-level genome of the transformable northern wattle, Acacia crassicarpa

The genus Acacia is a large group of woody legumes containing an enormous amount of morphological diversity in leaf shape. This diversity is at least in part the result of an innovation in leaf development where many Acacia species are capable of developing leaves of both bifacial and unifacial morphologies. While not unique in the plant kingdom, unifaciality is most commonly associated with monocots, and its developmental genetic mechanisms have yet to be explored beyond this group. In this study, we identify an accession of Acacia crassicarpa with high regeneration rates and isolate a clone for genome sequencing. We generate a chromosome-level assembly of this readily transformable clone, and using comparative analyses, confirm a whole-genome duplication unique to Caesalpinoid legumes. This resource will be important for future work examining genome evolution in legumes and the unique developmental genetic mechanisms underlying unifacial morphogenesis in Acacia.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

Physical seed dormancy in pea is genetically separable from seed coat thickness and roughness

Introduction

The seeds of wild pea (Pisum) exhibit marked physical dormancy due to impermeability of the seed coat to water, and the loss of this dormancy is thought to have been critical for domestication. Wild pea seed coats are also notably thick and rough, traits that have also reduced during domestication and are anecdotally linked to increased permeability. However, how these traits specifically interact with permeability is unclear.

Methods

To investigate this, we examined the genetic control of differences in seed coat characteristics between wild P. sativum ssp. humile and a non-dormant domesticated P. s. sativum accession in a recombinant inbred population. QTL effects were confirmed and their locations refined in segregating F4/5 populations.

Results

In this population we found a moderate correlation between testa thickness and permeability, and identified loci that affect them independently, suggesting no close functional association. However, the major loci affecting both testa thickness and permeability collocated closely with Mendel's pigmentation locus A, suggesting flavonoid compounds under its control might contribute significantly to both traits. We also show that seed coat roughness is oligogenic in this population, with the major locus independent of both testa thickness and permeability, suggesting selection for smooth seed was unlikely to be due to effects on either of these traits.

Discussion

Results indicate loss of seed coat dormancy during domestication was not primarily driven by reduced testa thickness or smooth seededness. The close association between major permeability and thickness QTL and Mendel's 'A' warrant further study, particularly regarding the role of flavonoids.

Light-sensitive short hypocotyl genes confer symbiotic nodule identity in the legume Medicago truncatula

Legumes produce specialized root nodules that are distinct from lateral roots in morphology and function, with nodules intracellularly hosting nitrogen-fixing bacteria. We have previously shown that a lateral root program underpins nodule initiation, but there must be additional developmental regulators that confer nodule identity. Here, we show two members of the LIGHT-SENSITIVE SHORT HYPOCOTYL (LSH) transcription factor family, predominantly known to define shoot meristem complexity and organ boundaries, function as regulators of nodule organ identity. In parallel to the root initiation program, LSH1/LSH2 recruit a program into the root cortex that mediates the divergence into nodules, in particular with cell divisions in the mid-cortex. This includes regulation of auxin and cytokinin, promotion of NODULE ROOT1/2 and Nuclear Factor YA1, and suppression of the lateral root program. A principal outcome of LSH1/LSH2 function is the production of cells able to accommodate nitrogen-fixing bacteria, a key feature unique to nodules.

Gene-edited Mtsoc1 triple mutant Medicago plants do not flower

Optimized flowering time is an important trait that ensures successful plant adaptation and crop productivity. SOC1-like genes encode MADS transcription factors, which are known to play important roles in flowering control in many plants. This includes the best-characterized eudicot model Arabidopsis thaliana (Arabidopsis), where SOC1 promotes flowering and functions as a floral integrator gene integrating signals from different flowering-time regulatory pathways. Medicago truncatula (Medicago) is a temperate reference legume with strong genomic and genetic resources used to study flowering pathways in legumes. Interestingly, despite responding to similar floral-inductive cues of extended cold (vernalization) followed by warm long days (VLD), such as in winter annual Arabidopsis, Medicago lacks FLC and CO which are key regulators of flowering in Arabidopsis. Unlike Arabidopsis with one SOC1 gene, multiple gene duplication events have given rise to three MtSOC1 paralogs within the Medicago genus in legumes: one Fabaceae group A SOC1 gene, MtSOC1a, and two tandemly repeated Fabaceae group B SOC1 genes, MtSOC1b and MtSOC1c. Previously, we showed that MtSOC1a has unique functions in floral promotion in Medicago. The Mtsoc1a Tnt1 retroelement insertion single mutant showed moderately delayed flowering in long- and short-day photoperiods, with and without prior vernalization, compared to the wild-type. In contrast, Mtsoc1b Tnt1 single mutants did not have altered flowering time or flower development, indicating that it was redundant in an otherwise wild-type background. Here, we describe the generation of Mtsoc1a Mtsoc1b Mtsoc1c triple mutant lines using CRISPR-Cas9 gene editing. We studied two independent triple mutant lines that segregated plants that did not flower and were bushy under floral inductive VLD. Genotyping indicated that these non-flowering plants were homozygous for the predicted strong mutant alleles of the three MtSOC1 genes. Gene expression analyses using RNA-seq and RT-qPCR indicated that these plants remained vegetative. Overall, the non-flowering triple mutants were dramatically different from the single Mtsoc1a mutant and the Arabidopsis soc1 mutant; implicating multiple MtSOC1 genes in critical overlapping roles in the transition to flowering in Medicago.

Genome assembly of Medicago truncatula accession SA27063 provides insight into spring black stem and leaf spot disease resistance

Medicago truncatula, model legume and alfalfa relative, has served as an essential resource for advancing our understanding of legume physiology, functional genetics, and crop improvement traits. Necrotrophic fungus, Ascochyta medicaginicola, the causal agent of spring black stem (SBS) and leaf spot is a devasting foliar disease of alfalfa affecting stand survival, yield, and forage quality. Host resistance to SBS disease is poorly understood, and control methods rely on cultural practices. Resistance has been observed in M. truncatula accession SA27063 (HM078) with two recessively inherited quantitative-trait loci (QTL), rnpm1 and rnpm2, previously reported. To shed light on host resistance, we carried out a de novo genome assembly of HM078. The genome, referred to as MtHM078 v1.0, is comprised of 23 contigs totaling 481.19 Mbp. Notably, this assembly contains a substantial amount of novel centromere-related repeat sequences due to deep long-read sequencing. Genome annotation resulted in 98.4% of BUSCO fabales proteins being complete. The assembly enabled sequence-level analysis of rnpm1 and rnpm2 for gene content, synteny, and structural variation between SBS-resistant accession SA27063 (HM078) and SBS-susceptible accession A17 (HM101). Fourteen candidate genes were identified, and some have been implicated in resistance to necrotrophic fungi. Especially interesting candidates include loss-of-function events in HM078 because they fit the inverse gene-for-gene model, where resistance is recessively inherited. In rnpm1, these include a loss-of-function in a disease resistance gene due to a premature stop codon, and a 10.85 kbp retrotransposon-like insertion disrupting a ubiquitin conjugating E2. In rnpm2, we identified a frameshift mutation causing a loss-of-function in a glycosidase, as well as a missense and frameshift mutation altering an F-box family protein. This study generated a high-quality genome of HM078 and has identified promising candidates, that once validated, could be further studied in alfalfa to enhance disease resistance.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

CRISPR/Cas9-Mediated Generation of Mutant Lines in Medicago truncatula Indicates a Symbiotic Role of MtLYK10 during Nodule Formation

CRISPR/Cas9 systems are commonly used for plant genome editing; however, the generation of homozygous mutant lines in Medicago truncatula remains challenging. Here, we present a CRISPR/Cas9-based protocol that allows the efficient generation of M. truncatula mutants. Gene editing was performed for the LysM receptor kinase gene MtLYK10 and two major facilitator superfamily transporter genes. The functionality of CRISPR/Cas9 vectors was tested in Nicotiana benthamiana leaves by editing a co-transformed GUSPlus gene. Transformed M. truncatula leaf explants were regenerated to whole plants at high efficiency (80%). An editing efficiency (frequency of mutations at a given target site) of up to 70% was reached in the regenerated plants. Plants with MtLYK10 knockout mutations were propagated, and three independent homozygous mutant lines were further characterized. No off-target mutations were identified in these lyk10 mutants. Finally, the lyk10 mutants and wild-type plants were compared with respect to the formation of root nodules induced by nitrogen-fixing Sinorhizobium meliloti bacteria. Nodule formation was considerably delayed in the three lyk10 mutant lines. Surprisingly, the size of the rare nodules in mutant plants was higher than in wild-type plants. In conclusion, the symbiotic characterization of lyk10 mutants generated with the developed CRISPR/Cas9 protocol indicated a role of MtLYK10 in nodule formation.

The WOX Genes from the Intermediate Clade: Influence on the Somatic Embryogenesis in Medicago truncatula

Transcription factors from the WOX family are well-known regulators of cell proliferation and differentiation in plants. Herein, we focused on several WOX genes from the intermediate clade and checked their impact on somatic embryogenesis using the model legume object Medicago truncatula. As a result, we show that MtWOX9-1 overexpression not only stimulates somatic embryogenesis in the embryogenic M. truncatula line, as it was shown previously, but can also induce somatic embryogenesis in the non-embryogenic line. Other intermediate clade WOX, including the close paralog of MtWOX9-1, as well as WOX11 homologs, did not have any significant impact on somatic embryogenesis in our in vitro cultivation system. Together, our results give new information about the diversity of the WOX family proteins and their specific functions. These data can be used for the search of new regeneration stimulators.

An emerging role of heterotrimeric G-proteins in nodulation and nitrogen sensing

MAIN CONCLUSION

A comprehensive understanding of nitrogen signaling cascades involving heterotrimeric G-proteins and their putative receptors can assist in the production of nitrogen-efficient plants. Plants are immobile in nature, so they must endure abiotic stresses including nutrient stress. Plant development and agricultural productivity are frequently constrained by the restricted availability of nitrogen in the soil. Non-legume plants acquire nitrogen from the soil through root membrane-bound transporters. In depleted soil nitrogen conditions, legumes are naturally conditioned to fix atmospheric nitrogen with the aid of nodulation elicited by nitrogen-fixing bacteria. Moreover, apart from the symbiotic nitrogen fixation process, nitrogen uptake from the soil can also be a significant secondary source to satisfy the nitrogen requirements of legumes. Heterotrimeric G-proteins function as molecular switches to help plant cells relay diverse stimuli emanating from external stress conditions. They are comprised of Gα, Gβ and Gγ subunits, which cooperate with several downstream effectors to regulate multiple plant signaling events. In the present review, we concentrate on signaling mechanisms that regulate plant nitrogen nutrition. Our review highlights the potential of heterotrimeric G-proteins, together with their putative receptors, to assist the legume root nodule symbiosis (RNS) cascade, particularly during calcium spiking and nodulation. Additionally, the functions of heterotrimeric G-proteins in nitrogen acquisition by plant roots as well as in improving nitrogen use efficiency (NUE) have also been discussed. Future research oriented towards heterotrimeric G-proteins through genome editing tools can be a game changer in the enhancement of the nitrogen fixation process. This will foster the precise manipulation and production of plants to ensure global food security in an era of climate change by enhancing crop productivity and minimizing reliance on external inputs.

The multi-kingdom microbiome of the goat gastrointestinal tract

BACKGROUND

Goat is an important livestock worldwide, which plays an indispensable role in human life by providing meat, milk, fiber, and pelts. Despite recent significant advances in microbiome studies, a comprehensive survey on the goat microbiomes covering gastrointestinal tract (GIT) sites, developmental stages, feeding styles, and geographical factors is still unavailable. Here, we surveyed its multi-kingdom microbial communities using 497 samples from ten sites along the goat GIT.

RESULTS

We reconstructed a goat multi-kingdom microbiome catalog (GMMC) including 4004 bacterial, 71 archaeal, and 7204 viral genomes and annotated over 4,817,256 non-redundant protein-coding genes. We revealed patterns of feeding-driven microbial community dynamics along the goat GIT sites which were likely associated with gastrointestinal food digestion and absorption capabilities and disease risks, and identified an abundance of large intestine-enriched genera involved in plant fiber digestion. We quantified the effects of various factors affecting the distribution and abundance of methane-producing microbes including the GIT site, age, feeding style, and geography, and identified 68 virulent viruses targeting the methane producers via a comprehensive virus-bacterium/archaea interaction network.

CONCLUSIONS

Together, our GMMC catalog provides functional insights of the goat GIT microbiota through microbiome-host interactions and paves the way to microbial interventions for better goat and eco-environmental qualities. Video Abstract.

Biological and genomic analysis of a symbiotic nitrogen fixation defective mutant in Medicago truncatula

Medicago truncatula has been selected as one of the model legume species for gene functional studies. To elucidate the functions of the very large number of genes present in plant genomes, genetic mutant resources are very useful and necessary tools. Fast Neutron (FN) mutagenesis is effective in inducing deletion mutations in genomes of diverse species. Through this method, we have generated a large mutant resource in M. truncatula. This mutant resources have been used to screen for different mutant using a forward genetics methods. We have isolated and identified a large amount of symbiotic nitrogen fixation (SNF) deficiency mutants. Here, we describe the detail procedures that are being used to characterize symbiotic mutants in M. truncatula. In recent years, whole genome sequencing has been used to speed up and scale up the deletion identification in the mutant. Using this method, we have successfully isolated a SNF defective mutant FN007 and identified that it has a large segment deletion on chromosome 3. The causal deletion in the mutant was confirmed by tail PCR amplication and sequencing. Our results illustrate the utility of whole genome sequencing analysis in the characterization of FN induced deletion mutants for gene discovery and functional studies in the M. truncatula. It is expected to improve our understanding of molecular mechanisms underlying symbiotic nitrogen fixation in legume plants to a great extent.

Balancing nitrate acquisition strategies in symbiotic legumes

MAIN CONCLUSION

Legumes manage both symbiotic (indirect) and non-symbiotic (direct) nitrogen acquisition pathways. Understanding and optimising the direct pathway for nitrate uptake will support greater legume growth and seed yields. Legumes have multiple pathways to acquire reduced nitrogen to grow and set seed. Apart from the symbiotic N₂-fixation pathway involving soil-borne rhizobia bacteria, the acquisition of nitrate and ammonia from the soil can also be an important secondary nitrogen source to meet plant N demand. The balance in N delivery between symbiotic N (indirect) and inorganic N uptake (direct) remains less clear over the growing cycle and with the type of legume under cultivation. In fertile, pH balanced agricultural soils, NO₃^- is often the predominant form of reduced N available to crop plants and will be a major contributor to whole plant N supply if provided at sufficient levels. The transport processes for NO₃^- uptake into legume root cells and its transport between root and shoot tissues involves both high and low-affinity transport systems called HATS and LATS, respectively. These proteins are regulated by external NO₃^- availability and by the N status of the cell. Other proteins also play a role in NO₃^- transport, including the voltage dependent chloride/nitrate channel family (CLC) and the S-type anion channels of the SLAC/SLAH family. CLC's are linked to NO₃^- transport across the tonoplast of vacuoles and the SLAC/SLAH's with NO₃^- efflux across the plasma membrane and out of the cell. An important step in managing the N requirements of a plant are the mechanisms involved in root N uptake and the subsequent cellular distribution within the plant. In this review, we will present the current knowledge of these proteins and what is understood on how they function in key model legumes (Lotus japonicus, Medicago truncatula and Glycine sp.). The review will examine their regulation and role in N signalling, discuss how post-translational modification affects NO₃^- transport in roots and aerial tissues and its translocation to vegetative tissues and storage/remobilization in reproductive tissues. Lastly, we will present how NO₃^-influences the autoregulation of nodulation and nitrogen fixation and its role in mitigating salt and other abiotic stresses.

Genome-Wide Identification and Comparative Analysis of RALF Gene Family in Legume and Non-Legume Species

Rapid alkalinization factor (RALF) are small secreted peptide hormones that can induce rapid alkalinization in a medium. They act as signaling molecules in plants, playing a critical role in plant development and growth, especially in plant immunity. Although the function of RALF peptides has been comprehensively analyzed, the evolutionary mechanism of RALFs in symbiosis has not been studied. In this study, 41, 24, 17 and 12 RALFs were identified in Arabidopsis, soybean, Lotus and Medicago, respectively. A comparative analysis including the molecular characteristics and conserved motifs suggested that the RALF pre-peptides in soybean represented a higher value of isoelectric point and more conservative motifs/residues composition than other species. All 94 RALFs were divided into two clades according to the phylogenetic analysis. Chromosome distribution and synteny analysis suggested that the expansion of the RALF gene family in Arabidopsis mainly depended on tandem duplication, while segment duplication played a dominant role in legume species. The expression levels of most RALFs in soybean were significantly affected by the treatment of rhizobia. Seven GmRALFs are potentially involved in the release of rhizobia in the cortex cells. Overall, our research provides novel insights into the understanding of the role of the RALF gene family in nodule symbiosis.

The function of the Medicago truncatula ZIP transporter MtZIP14 is linked to arbuscular mycorrhizal fungal colonization

Soil micronutrient availability, including zinc (Zn), is a limiting factor for crop yield. Arbuscular mycorrhizal (AM) fungi can improve host plant growth and nutrition through the mycorrhizal pathway of nutrient uptake. Although the physiology of Zn uptake through the mycorrhizal pathway is well established, the identity of the related molecular components are unknown. Here, RNA-seq analysis was used to identify genes differentially-regulated by AM colonization and soil Zn concentration in roots of Medicago truncatula. The putative Zn transporter gene MtZIP14 was markedly up-regulated in M. truncatula roots when colonized by Rhizophagus irregularis. MtZIP14 restored yeast growth under low Zn availability. Loss-of-function mutant plants (mtzip14) had reduced shoot biomass compared to the wild-type when colonized by AM fungi and grown under low and sufficient soil Zn concentration; at high soil Zn concentration, there were no genotypic differences in shoot biomass. The vesicular and arbuscular colonization of roots was lower in the mtzip14 plants regardless of soil Zn concentration. We propose that MtZIP14 is linked to AM colonization in M. truncatula plants, with the possibility that MtZIP14 function with AM colonization is linked to plant Zn nutrition.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Variability in Responses to Phoma medicaginis Infection in a Tunisian Collection of Three Annual Medicago Species

Spring black stem and leaf spot, caused by Phoma medicaginis, is an issue in annual Medicago species. Therefore, in this study, we analyzed the response to P. medicaginis infection in a collection of 46 lines of three annual Medicago species (M. truncatula, M. ciliaris, and M. polymorpha) showing different geographic distribution in Tunisia. The reaction in the host to the disease is explained by the effects based on plant species, lines nested within species, treatment, the interaction of species × treatment, and the interaction of lines nested within species × treatment. Medicago ciliaris was the least affected for aerial growth under infection. Furthermore, the largest variation within species was found for M. truncatula under both conditions. Principal component analysis and hierarchical classification showed that M. ciliaris lines formed a separate group under control treatment and P. medicaginis infection and they are the most vigorous in growth. These results indicate that M. ciliaris is the least susceptible in response to P. medicaginis infection among the three Medicago species investigated here, which can be used as a good candidate in crop rotation to reduce disease pressure in the field and as a source of P. medicaginis resistance for the improvement of forage legumes.

Genome-Wide Mining of the Tandem Duplicated Type III Polyketide Synthases and Their Expression, Structure Analysis of Senna tora

Senna tora is one of the homologous crops used as a medicinal food containing an abundance of anthraquinones. Type III polyketide synthases (PKSs) are key enzymes that catalyze polyketide formation; in particular, the chalcone synthase-like (CHS-L) genes are involved in anthraquinone production. Tandem duplication is a fundamental mechanism for gene family expansion. However, the analysis of the tandem duplicated genes (TDGs) and the identification and characterization of PKSs have not been reported for S. tora. Herein, we identified 3087 TDGs in the S. tora genome; the synonymous substitution rates (Ks) analysis indicated that the TDGs had recently undergone duplication. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the type III PKSs were the most enriched TDGs involved in the biosynthesis of the secondary metabolite pathways, as evidenced by 14 tandem duplicated CHS-L genes. Subsequently, we identified 30 type III PKSs with complete sequences in the S. tora genome. Based on the phylogenetic analysis, the type III PKSs were classified into three groups. The protein conserved motifs and key active residues showed similar patterns in the same group. The transcriptome analysis showed that the chalcone synthase (CHS) genes were more highly expressed in the leaves than in the seeds in S. tora. The transcriptome and qRT-PCR analysis showed that the CHS-L genes had a higher expression in the seeds than in other tissues, particularly seven tandem duplicated CHS-L2/3/5/6/9/10/13 genes. The key active-site residues and three-dimensional models of the CHS-L2/3/5/6/9/10/13 proteins showed slight variation. These results indicated that the rich anthraquinones in S. tora seeds might be ascribed to the PKSs' expansion from tandem duplication, and the seven key CHS-L2/3/5/6/9/10/13 genes provide candidate genes for further research. Our study provides an important basis for further research on the regulation of anthraquinones' biosynthesis in S. tora.

Medicago truncatula resources to study legume biology and symbiotic nitrogen fixation

Medicago truncatula is a chosen model for legumes towards deciphering fundamental legume biology, especially symbiotic nitrogen fixation. Current genomic resources for M. truncatula include a completed whole genome sequence information for R108 and Jemalong A17 accessions along with the sparse draft genome sequences for other 226 M. truncatula accessions. These genomic resources are complemented by the availability of mutant resources such as retrotransposon (Tnt1) insertion mutants in R108 and fast neutron bombardment (FNB) mutants in A17. In addition, several M. truncatula databases such as small secreted peptides (SSPs) database, transporter protein database, gene expression atlas, proteomic atlas, and metabolite atlas are available to the research community. This review describes these resources and provide information regarding how to access these resources.

Long-lasting impact of chitooligosaccharide application on strigolactone biosynthesis and fungal accommodation promotes arbuscular mycorrhiza in Medicago truncatula

The establishment of arbuscular mycorrhiza (AM) between plants and Glomeromycotina fungi is preceded by the exchange of chemical signals: fungal released Myc-factors, including chitooligosaccharides (CO) and lipo-chitooligosaccharides (LCO), activate plant symbiotic responses, while root-exuded strigolactones stimulate hyphal branching and boost CO release. Furthermore, fungal signaling reinforcement through CO application was shown to promote AM development in Medicago truncatula, but the cellular and molecular bases of this effect remained unclear. Here, we focused on long-term M. truncatula responses to CO treatment, demonstrating its impact on the transcriptome of both mycorrhizal and nonmycorrhizal roots over several weeks and providing an insight into the mechanistic bases of the CO-dependent promotion of AM colonization. CO treatment caused the long-lasting regulation of strigolactone biosynthesis and fungal accommodation-related genes. This was mirrored by an increase in root didehydro-orobanchol content, and the promotion of accommodation responses to AM fungi in root epidermal cells. Lastly, an advanced downregulation of AM symbiosis marker genes was observed at the latest time point in CO-treated plants, in line with an increased number of senescent arbuscules. Overall, CO treatment triggered molecular, metabolic, and cellular responses underpinning a protracted acceleration of AM development.

Genome-Wide Identification and Expression Analysis of the Ammonium Transporter Family Genes in Soybean

Ammonium transporters (AMTs) are responsible for ammonium absorption and utilization in plants. As a high-nitrogen-demand crop and a legume, soybean can also obtain ammonium from symbiotic root nodules in which nitrogen-fixing rhizobia convert atmospheric nitrogen (N₂) into ammonium. Although increasing evidence implicates vital roles of ammonium transport in soybean, no systematic analyses of AMTs in soybean (named GmAMTs) or functional analyses of GmAMTs are available. In this study, we aimed to identify all GmAMT family genes and gain a better understanding of the characteristics of GmAMT genes in soybean. Here, due to the improved genome assembly and annotation of soybean, we tried to generate a phylogenetic tree of 16 GmAMTs based on new information. Consistent with reported data, GmAMT family members can be divided into two subfamilies of GmAMT1 (6 genes) and GmAMT2 (10 genes). Interestingly, unlike Arabidopsis, which has only one AMT2, soybean has substantially increased the number of GmAMT2s, suggesting enhanced demand for ammonium transport. These genes were distributed on nine chromosomes, of which GmAMT1.3, GmAMT1.4, and GmAMT1.5 were three tandem repeat genes. The gene structures and conserved protein motifs of the GmAMT1 and GmAMT2 subfamilies were different. All the GmAMTs were membrane proteins with varying numbers of transmembrane domains ranging from 4 to 11. Promoter analysis found that these GmAMT genes have phytohormone-, circadian control-, and organ expression-related cis-elements in their promoters, and notably, there were nodulation-specific and nitrogen-responsive elements in the promoters of the GmAMT1 and GmAMT2 genes. Further expression data showed that these GmAMT family genes exhibited different spatiotemporal expression patterns across tissues and organs. In addition, GmAMT1.1, GmAMT1.2, GmAMT2.2, and GmAMT2.3 were responsive to nitrogen treatment, while GmAMT1.2, GmAMT1.3, GmAMT1.4, GmAMT1.5, GmAMT1.6, GmAMT2.1, GmAMT2.2, GmAMT2.3, GmAMT3.1, and GmAMT4.6 showed circadian rhythms in transcription. RT-qPCR validated the expression patterns of GmAMTs in response to different forms of nitrogen and exogenous ABA treatments. Gene expression analysis also confirmed that GmAMTs are regulated by key nodulation gene GmNINa, indicating a role of GmAMTs in symbiosis. Together, these data indicate that GmAMTs may differentially and/or redundantly regulate ammonium transport during plant development and in response to environmental factors. These findings provide a basis for future research on the functions of GmAMTs and the mechanisms through which GmAMTs regulate ammonium metabolism and nodulation in soybean.

A Germin-Like Protein GLP1 of Legumes Mediates Symbiotic Nodulation by Interacting with an Outer Membrane Protein of Rhizobia

Rhizobia can infect legumes and induce the coordinated expression of symbiosis and defense genes for the establishment of mutualistic symbiosis. Numerous studies have elucidated the molecular interactions between rhizobia and host plants, which are associated with Nod factor, exopolysaccharide, and T3SS effector proteins. However, there have been relatively few reports about how the host plant recognizes the outer membrane proteins (OMPs) of rhizobia to mediate symbiotic nodulation. In our previous work, a gene (Mhopa22) encoding an OMP was identified in Mesorhizobium huakuii 7653R, whose homologous genes are widely distributed in Rhizobiales. In this study, a germin-like protein GLP1 interacting with Mhopa22 was identified in Astragalus sinicus. RNA interference of AsGLP1 resulted in a decrease in nodule number, whereas overexpression of AsGLP1 increased the number of nodules in the hairy roots of A. sinicus. Consistent symbiotic phenotypes were identified in Medicago truncatula with MtGLPx (refer to medtr7g111240.1, the isogeny of AsGLP1) overexpression or Tnt1 mutant (glpx-1) in symbiosis with Sinorhizobium meliloti 1021. The glpx-1 mutant displayed hyperinfection and the formation of more infection threads but a decrease in root nodules. RNA sequencing analysis showed that many differentially expressed genes were involved in hormone signaling and symbiosis. Taken together, AsGLP1 and its homology play an essential role in mediating the early symbiotic process through interacting with the OMPs of rhizobia. IMPORTANCE This study is the first report to characterize a legume host plant protein to sense and interact with an outer membrane protein (OMP) of rhizobia. It can be speculated that GLP1 plays an essential role to mediate early symbiotic process through interacting with OMPs of rhizobia. The results provide deeper understanding and novel insights into the molecular interactive mechanism of a legume symbiosis signaling pathway in recognition with rhizobial OMPs. Our findings may also provide a new perspective to improve the symbiotic compatibility and nodulation of legume.

Genome-wide identification of the TGA genes in common bean (Phaseolus vulgaris) and revealing their functions in response to Fusarium oxysporum f. sp. phaseoli infection

Fusarium wilt, which affects common bean all across the world, is caused by Fusarium oxysporum f. sp. Phaseoli (Fop). It is necessary to have functional genes in response to Fop infection because they might be used to manage disease. As a crucial regulator, TGA-binding transcription factor (TGA) is engaged in the defense mechanism of plants against pathogens. The role of TGA regulators in common bean in response to Fop infection, however, has not been documented. Hence, we performed genome-wide identified and characterized eight TGA genes in common bean. In this study, eight PvTGA genes were distributed on six chromosomes and classified into four subgroups. The PvTGA genes have the same conserved bZIP and DOG1 domains, but there are specific sequence structures in different PvTGAs. Phylogenetic and synteny analysis explained that PvTGA gene has a close genetic relationship with legume TGAs and that PvTGA03 and PvTGA05 may play an important role in evolution. Transcriptome data explained that expression levels of PvTGA genes showed diversity in different tissues. After Fop inoculation, the expression levels of PvTGA03 and PvTGA07 were significantly different between resistant and susceptible genotypes. Under SA treatment, the expression levels of PvTGA03, PvTGA04, PvTGA06, PvTGA07 and PvTGA08 were significantly different. These results imply that PvTGA03 and PvTGA07 play key roles in SA-mediated resistance to Fusarium wilt. Together, these findings advance knowledge of the PvTGA gene family in common bean and will help future studies aimed at reducing Fusarium wilt.

MtCLE08, MtCLE16, and MtCLE18 Transcription Patterns and Their Possible Functions in the Embryogenic Calli of Medicago truncatula

CLE peptides are well-known hormonal regulators of plant development, but their role in somatic embryogenesis remains undetermined. CLE genes are often regulated by WOX transcription factors and, in their turn, regulate the expression level of WOX genes. In this study, we used in vitro cultivation, as well as qPCR and transcriptomic analysis, to find CLE peptides which could regulate the MtWOX9-1 gene, stimulating somatic embryogenesis in Medicago truncatula. Three CLE peptides were found which could probably be such regulators, but none of them was found to influence MtWOX9-1 expression in the embryogenic calli. Nevertheless, overexpression of one of CLE genes under study, MtCLE16, decreased somatic embryogenesis intensity. Additionally, overexpression of MtCLE08 was found to suppress expression of MtWOX13a, a supposed antagonist of somatic embryo development. Our findings can be helpful for the search for new regeneration regulators which could be used for plant transformation.

Genome-wide identification, evolution, and expression analysis of carbonic anhydrases genes in soybean (Glycine max)

Carbonic anhydrases (CAs), as zinc metalloenzymes, are ubiquitous in nature and play essential roles in diverse biological processes. Although CAs have been broadly explored and studied, comprehensive characteristics of CA gene family members in the soybean (Glycine max) are still lacking. A total of 35 CA genes (GmCAs) were identified; they distributed on sixteen chromosomes of the soybean genome and can be divided into three subfamilies (α-type, β-type, and γ-type). Bioinformatics analysis showed that the specific GmCA gene subfamily or clade exhibited similar characteristics and that segmental duplications took the major role in generating new GmCAs. Furthermore, the synteny and evolutionary constraints analyses of CAs among soybean and distinct species provided more detailed evidence for GmCA gene family evolution. Cis-element analysis of promoter indicated that GmCAs may be responsive to abiotic stress and regulate photosynthesis. Moreover, the expression patterns of GmCAs varied in different tissues at diverse developmental stages in soybean. Additionally, we found that eight representative GmCAs may be involved in the response of soybean to low phosphorus stress. The systematic investigation of the GmCA gene family in this study will provide a valuable basis for further functional research on soybean CA genes.

Methodology for Constructing a Knowledgebase for Plant Gene Regulation Information

The amount of biological data is growing at a rapid pace as many high-throughput omics technologies and data pipelines are developed. This is resulting in the growth of databases for DNA and protein sequences, gene expression, protein accumulation, structural, and localization information. The diversity and multi-omics nature of such bioinformatic data requires well-designed databases for flexible organization and presentation. Besides general-purpose online bioinformatic databases, users need narrowly focused online databases to quickly access a meaningful collection of related data for their research. Here, we describe the methodology used to implement a plant gene regulatory knowledgebase, with data, query, and tool features, as well as the ability to expand to accommodate future datasets. We exemplify this methodology for the GRASSIUS knowledgebase, but it is applicable to developing and updating similar plant gene regulatory knowledgebases. GRASSIUS organizes and presents gene regulatory data from grass species with a central focus on maize (Zea mays). The main class of data presented include not only the families of transcription factors (TFs) and co-regulators (CRs) but also protein-DNA interaction data, where available.

Medicago truncatula quantitative resistance to a new strain of Verticillium alfalfae from Iran revealed by a genome-wide association study

Verticillium wilt is a major threat to many crops, among them alfalfa (Medicago sativa). The model plant Medicago truncatula, a close relative of alfalfa was used to study the genetic control of resistance towards a new Verticillium alfalfae isolate. The accidental introduction of pathogen strains through global trade is a threat to crop production and such new strains might also be better adapted to global warming. Isolates of V. alfalfae were obtained from alfalfa fields in Iran and characterized. The Iranian isolate AF1 was used in a genome-wide association study (GWAS) involving 242 accessions from the Mediterranean region. Root inoculations were performed with conidia at 25°C and symptoms were scored regularly. Maximum Symptom Score and Area under Disease Progess Curve were computed as phenotypic traits to be used in GWAS and for comparison to a previous study with French isolate V31.2 at 20°C. This comparison showed high correlation with a shift to higher susceptibility, and similar geographical distribution of resistant and susceptible accessions to AF1 at 25°C, with resistant accessions mainly in the western part. GWAS revealed 30 significant SNPs linked to resistance towards isolate AF1. None of them were common to the previous study with isolate V31.2 at 20°C. To confirm these loci, the expression of nine underlying genes was studied. All genes were induced in roots following inoculation, in susceptible and resistant plants. However, in resistant plants induction was higher and lasted longer. Taken together, the use of a new pathogen strain and a shift in temperature revealed a completely different genetic control compared to a previous study that demonstrated the existence of two major QTLs. These results can be useful for Medicago breeding programs to obtain varieties better adapted to future conditions.

Alfalfa genomic selection for different stress-prone growing regions

Alfalfa (Medicago sativa L.) selection for stress-prone regions has high priority for sustainable crop-livestock systems. This study assessed the genomic selection (GS) ability to predict alfalfa breeding values for drought-prone agricultural sites of Algeria, Morocco, and Argentina; managed-stress (MS) environments of Italy featuring moderate or intense drought; and one Tunisian site irrigated with moderately saline water. Additional aims were to investigate genotype × environment interaction (GEI) patterns and the effect on GS predictions of three single-nucleotide polymorphism (SNP) calling procedures, 12 statistical models that exclude or incorporate GEI, and allele dosage information. Our study included 127 genotypes from a Mediterranean reference population originated from three geographically contrasting populations, genotyped via genotyping-by-sequencing and phenotyped based on multi-year biomass dry matter yield of their dense-planted half-sib progenies. The GEI was very large, as shown by 27-fold greater additive genetic variance × environment interaction relative to the additive genetic variance and low genetic correlation for progeny yield responses across environments. The predictive ability of GS (using at least 37,969 SNP markers) exceeded 0.20 for moderate MS (representing Italian stress-prone sites) and the sites of Algeria and Argentina while being quite low for the Tunisian site and intense MS. Predictions of GS were complicated by rapid linkage disequilibrium decay. The weighted GBLUP model, GEI incorporation into GS models, and SNP calling based on a mock reference genome exhibited a predictive ability advantage for some environments. Our results support the specific breeding for each target region and suggest a positive role for GS in most regions when considering the challenges associated with phenotypic selection.

Sequence analysis and crystal structure of a glycosylated protease from Euphorbia resinifera latex for its proteolytic activity aspect

The investigation of a plant glycosylated serine protease (EuRP-61) isolated from Euphorbia resinifera latex for potential antiplatelet and anticoagulation activities has been previously reported. In the present study, the protein sequence and native crystal structure of EuRP-61 were characterized. The structure was identified using single-wavelength anomalous diffraction with a refinement resolution of 1.7 Å (PDB ID: 7EOX). The main structural components of EuRP-61 were composed of three domains: catalytic, protease-associated (PA), and fibronectin type III (Fn3)-like domains. The crystal structure revealed that some loops in the PA and catalytic domains of EuRP-61 were different from the other subtilisin-like proteases (cucumisin and SBT3). These different loops might be involved in the general monomer formation of EuRP-61, substrate specificity, and maintenance of the catalytic domain. The Fn3-like domain may provide flexibility to the enzyme to bind with various substrates and cell receptors. Additionally, the active site of EuRP-61 consisted of the catalytic triad of Ser434, His106, and Asp32, similar to other serine proteases. The present study provides additional information and insight into the protease and antithrombotic activities of EuRP-61, which could contribute to further development of this enzyme for biomedical treatment.

Phenotypic variation and quantitative trait loci for resistance to southern anthracnose and clover rot in red clover

High variability for and candidate loci associated with resistance to southern anthracnose and clover rot in a worldwide collection of red clover provide a first basis for genomics-assisted breeding. Red clover (Trifolium pratense L.) is an important forage legume of temperate regions, particularly valued for its high yield potential and its high forage quality. Despite substantial breeding progress during the last decades, continuous improvement of cultivars is crucial to ensure yield stability in view of newly emerging diseases or changing climatic conditions. The high amount of genetic diversity present in red clover ecotypes, landraces, and cultivars provides an invaluable, but often unexploited resource for the improvement of key traits such as yield, quality, and resistance to biotic and abiotic stresses. A collection of 397 red clover accessions was genotyped using a pooled genotyping-by-sequencing approach with 200 plants per accession. Resistance to the two most pertinent diseases in red clover production, southern anthracnose caused by Colletotrichum trifolii, and clover rot caused by Sclerotinia trifoliorum, was assessed using spray inoculation. The mean survival rate for southern anthracnose was 22.9% and the mean resistance index for clover rot was 34.0%. Genome-wide association analysis revealed several loci significantly associated with resistance to southern anthracnose and clover rot. Most of these loci are in coding regions. One quantitative trait locus (QTL) on chromosome 1 explained 16.8% of the variation in resistance to southern anthracnose. For clover rot resistance we found eight QTL, explaining together 80.2% of the total phenotypic variation. The SNPs associated with these QTL provide a promising resource for marker-assisted selection in existing breeding programs, facilitating the development of novel cultivars with increased resistance against two devastating fungal diseases of red clover.

Validation of a QTL on Chromosome 1DS Showing a Major Effect on Salt Tolerance in Winter Wheat

Salt stress is one the most destructive abiotic stressors, causing yield losses in wheat worldwide. A prerequisite for improving salt tolerance is the identification of traits for screening genotypes and uncovering causative genes. Two populations of F₃ lines developed from crosses between sensitive and tolerant parents were tested for salt tolerance at the seedling stage. Based on their response, the offspring were classified as salt sensitive and tolerant. Under saline conditions, tolerant genotypes showed lower Na⁺ and proline content but higher K⁺, higher chlorophyll content, higher K⁺/Na⁺ ratio, higher PSII activity levels, and higher photochemical efficiency, and were selected for further molecular analysis. Five stress responsive QTL identified in a previous study were validated in the populations. A QTL on the short arm of chromosome 1D showed large allelic effects in several salt tolerant related traits. An expression analysis of associated candidate genes showed that TraesCS1D02G052200 and TraesCS5B02G368800 had the highest expression in most tissues. Furthermore, qRT-PCR expression analysis revealed that ZIP-7 had higher differential expressions under saline conditions compared to KefC, AtABC8 and 6-SFT. This study provides information on the genetic and molecular basis of salt tolerance that could be useful in development of salt-tolerant wheat varieties.

MtING2 encodes an ING domain PHD finger protein which affects Medicago growth, flowering, global patterns of H3K4me3, and gene expression

Flowering of the reference legume Medicago truncatula is promoted by winter cold (vernalization) followed by long-day photoperiods (VLD) similar to winter annual Arabidopsis. However, Medicago lacks FLC and CO, key regulators of Arabidopsis VLD flowering. Most plants have two INHIBITOR OF GROWTH (ING) genes (ING1 and ING2), encoding proteins with an ING domain with two anti-parallel alpha-helices and a plant homeodomain (PHD) finger, but their genetic role has not been previously described. In Medicago, Mting1 gene-edited mutants developed and flowered normally, but an Mting2-1 Tnt1 insertion mutant and gene-edited Mting2 mutants had developmental abnormalities including delayed flowering particularly in VLD, compact architecture, abnormal leaves with extra leaflets but no trichomes, and smaller seeds and barrels. Mting2 mutants had reduced expression of activators of flowering, including the FT-like gene MtFTa1, and increased expression of the candidate repressor MtTFL1c, consistent with the delayed flowering of the mutant. MtING2 overexpression complemented Mting2-1, but did not accelerate flowering in wild type. The MtING2 PHD finger bound H3K4me2/3 peptides weakly in vitro, but analysis of gene-edited mutants indicated that it was dispensable to MtING2 function in wild-type plants. RNA sequencing experiments indicated that >7000 genes are mis-expressed in the Mting2-1 mutant, consistent with its strong mutant phenotypes. Interestingly, ChIP-seq analysis identified >5000 novel H3K4me3 locations in the genome of Mting2-1 mutants compared to wild type R108. Overall, our mutant study has uncovered an important physiological role of a plant ING2 gene in development, flowering, and gene expression, which likely involves an epigenetic mechanism.

Nutrient regulation of lipochitooligosaccharide recognition in plants via NSP1 and NSP2

Many plants associate with arbuscular mycorrhizal fungi for nutrient acquisition, while legumes also associate with nitrogen-fixing rhizobial bacteria. Both associations rely on symbiosis signaling and here we show that cereals can perceive lipochitooligosaccharides (LCOs) for activation of symbiosis signaling, surprisingly including Nod factors produced by nitrogen-fixing bacteria. However, legumes show stringent perception of specifically decorated LCOs, that is absent in cereals. LCO perception in plants is activated by nutrient starvation, through transcriptional regulation of Nodulation Signaling Pathway (NSP)1 and NSP2. These transcription factors induce expression of an LCO receptor and act through the control of strigolactone biosynthesis and the karrikin-like receptor DWARF14-LIKE. We conclude that LCO production and perception is coordinately regulated by nutrient starvation to promote engagement with mycorrhizal fungi. Our work has implications for the use of both mycorrhizal and rhizobial associations for sustainable productivity in cereals.

α-CAs from Photosynthetic Organisms

Carbonic anhydrases (CAs) are ubiquitous enzymes that catalyze the reversible carbon dioxide hydration reaction. Among the eight different CA classes existing in nature, the α-class is the largest one being present in animals, bacteria, protozoa, fungi, and photosynthetic organisms. Although many studies have been reported on these enzymes, few functional, biochemical, and structural data are currently available on α-CAs isolated from photosynthetic organisms. Here, we give an overview of the most recent literature on the topic. In higher plants, these enzymes are engaged in both supplying CO₂ at the Rubisco and determining proton concentration in PSII membranes, while in algae and cyanobacteria they are involved in carbon-concentrating mechanism (CCM), photosynthetic reactions and in detecting or signaling changes in the CO₂ level in the environment. Crystal structures are only available for three algal α-CAs, thus not allowing to associate specific structural features to cellular localizations or physiological roles. Therefore, further studies on α-CAs from photosynthetic organisms are strongly needed to provide insights into their structure–function relationship.

Improved pea reference genome and pan-genome highlight genomic features and evolutionary characteristics

Complete and accurate reference genomes and annotations provide fundamental resources for functional genomics and crop breeding. Here we report a de novo assembly and annotation of a pea cultivar ZW6 with contig N50 of 8.98 Mb, which features a 243-fold increase in contig length and evident improvements in the continuity and quality of sequence in complex repeat regions compared with the existing one. Genome diversity of 118 cultivated and wild pea demonstrated that Pisum abyssinicum is a separate species different from P. fulvum and P. sativum within Pisum. Quantitative trait locus analyses uncovered two known Mendel's genes related to stem length (Le/le) and seed shape (R/r) as well as some candidate genes for pod form studied by Mendel. A pan-genome of 116 pea accessions was constructed, and pan-genes preferred in P. abyssinicum and P. fulvum showed distinct functional enrichment, indicating the potential value of them as pea breeding resources in the future.

Interference between ER stress-related bZIP-type and jasmonate-inducible bHLH-type transcription factors in the regulation of triterpene saponin biosynthesis in Medicago truncatula

Triterpene saponins (TS) are a structurally diverse group of metabolites that are widely distributed in plants. They primarily serve as defense compounds and their production is often triggered by biotic stresses through signaling cascades that are modulated by phytohormones such as the jasmonates (JA). Two JA-modulated basic helix-loop-helix (bHLH) transcription factors (TFs), triterpene saponin biosynthesis activating regulator 1 (TSAR1) and TSAR2, have previously been identified as direct activators of TS biosynthesis in the model legume Medicago truncatula. Here, we report on the involvement of the core endoplasmic reticulum (ER) stress-related basic leucine zipper (bZIP) TFs bZIP17 and bZIP60 in the regulation of TS biosynthesis. Expression and processing of M. truncatula bZIP17 and bZIP60 proteins were altered in roots with perturbed TS biosynthesis or treated with JA. Accordingly, such roots displayed an altered ER network structure. M. truncatula bZIP17 and bZIP60 proteins were shown to localize in the nucleus and appeared to be capable of interfering with the TSAR-mediated transactivation of TS biosynthesis genes. Furthermore, interference between ER stress-related bZIP and JA-modulated bHLH TFs in the regulation of JA-dependent terpene biosynthetic pathways may be widespread in the plant kingdom, as we demonstrate that it also occurs in the regulation of monoterpene indole alkaloid biosynthesis in the medicinal plant Catharanthus roseus.

MtNF-YC6 and MtNF-YC11 are involved in regulating the transcriptional program of arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal fungi are obligate symbionts that transfer mineral nutrients to host plants through arbuscules, a fungal structure specialized for exchange for photosynthetic products. MtNF-YC6 and MtNF-YC11, which encode the C subunits of nuclear factor Y (NF-Y) family in Medicago truncatula are induced specifically by arbuscular mycorrhizal symbiosis (AMS). A previous study showed that MtNF-YC6 and MtNF-YC11 are activated in cortical cells of mycorrhizal roots, but the gene functions were unknown. Herein, we identified both MtNF-YB17 and MtNF-YB12 as the interacting partners of MtNF-YC6 and MtNF-YC11 in yeast and plants. MtNF-YB17 was highly induced by AMS and activated in cortical cells only in mycorrhizal roots but MtNF-YB12 was not affected. The formation of B/C heterodimers led the protein complexes to transfer from the cytoplasm to the nucleus. Silencing MtNF-YC6 and C11 by RNA interference (RNAi) resulted in decreased colonization efficiency and arbuscule richness. Coincidently, genes associated with arbuscule development and degeneration in RNAi roots were also downregulated. In silico analysis showed CCAAT-binding motifs in the promoter regions of downregulated genes, further supporting the involvement of NF-Y complexes in transcriptional regulation of symbiosis. Taken together, this study identifies MtNF-YC6- or MtNF-YC11-containing protein complexes as novel transcriptional regulators of symbiotic program and provides a list of potential downstream target genes. These data will help to further dissect the AMS regulatory network.

The importance of the urea cycle and its relationships to polyamine metabolism during ammonium stress in Medicago truncatula

The ornithine-urea cycle (urea cycle) makes a significant contribution to the metabolic responses of lower photosynthetic eukaryotes to episodes of high nitrogen availability. In this study, we compared the role of the plant urea cycle and its relationships to polyamine metabolism in ammonium-fed and nitrate-fed Medicago truncatula plants. High ammonium resulted in the accumulation of ammonium and pathway intermediates, particularly glutamine, arginine, ornithine, and putrescine. Arginine decarboxylase activity was decreased in roots, suggesting that the ornithine decarboxylase-dependent production of putrescine was important in situations of ammonium stress. The activity of copper amine oxidase, which releases ammonium from putrescine, was significantly decreased in both shoots and roots. In addition, physiological concentrations of ammonium inhibited copper amine oxidase activity in in vitro assays, supporting the conclusion that high ammonium accumulation favors putrescine synthesis. Moreover, early supplementation of plants with putrescine avoided ammonium toxicity. The levels of transcripts encoding urea-cycle-related proteins were increased and transcripts involved in polyamine catabolism were decreased under high ammonium concentrations. We conclude that the urea cycle and associated polyamine metabolism function as important protective mechanisms limiting ammonium toxicity in M. truncatula. These findings demonstrate the relevance of the urea cycle to polyamine metabolism in higher plants.

Tolerant mechanism of model legume plant Medicago truncatula to drought, salt, and cold stresses

Legume plants produce one-third of the total yield of primary crops and are important food sources for both humans and animals worldwide. Frequent exposure to abiotic stresses, such as drought, salt, and cold, greatly limits the production of legume crops. Several morphological, physiological, and molecular studies have been conducted to characterize the response and adaptation mechanism to abiotic stresses. The tolerant mechanisms of the model legume plant Medicago truncatula to abiotic stresses have been extensively studied. Although many potential genes and integrated networks underlying the M. truncatula in responding to abiotic stresses have been identified and described, a comprehensive summary of the tolerant mechanism is lacking. In this review, we provide a comprehensive summary of the adaptive mechanism by which M. truncatula responds to drought, salt, and cold stress. We also discuss future research that need to be explored to improve the abiotic tolerance of legume plants.

Flavonoid Accumulation Varies in Medicago truncatula in Response to Mercury Stress

Mercury (Hg) contamination is increasing worldwide in both wild ecosystems and agricultural soils due to natural processes, but mostly to anthropic activities. The molecular mechanisms involved in Hg toxicity and tolerance in plants have been extensively studied; however, the role of flavonoids in response to Hg stress remains to be investigated. We conducted a metabolomic study to analyze the changes induced at the secondary metabolite level in three Hg-tolerant and one Hg-sensitive Medicago truncatula cultivars. A total of 46 flavonoid compounds, classified into five different flavonoid families: anthocyanidins, flavones, isoflavones, pterocarpan flavonoids, and flavanones, along with their respective glycoconjugate derivatives, were identified in leaf and root tissues. The synthesis of free isoflavones, followed by monoglycosylation and further malonylation was shown to be characteristic of root samples, whereas higher glycosylation, followed by further acylation with coumaric and ferulic acid was characteristic of leaf tissues. While minor changes were observed in leaves, significant quantitative changes could be observed in roots upon Hg treatment. Some flavonoids were strongly upregulated in roots, including malonylglucosides of biochanin A, formononetin and medicarpin, and aglycones biochanin, daidzein, and irisolidone. Hg tolerance appeared to be mainly associated to the accumulation of formononetin MalGlc, tricin GlcAGlcA, and afrormosin Glc II in leaves, whereas aglycone accumulation was associated with tolerance to Hg stress in roots. The results evidence the alteration of the flavonoid metabolic profile and their glycosylation processes in response to Hg stress. However, notable differences existed between varieties, both in the basal metabolic profile and in the response to treatment with Hg. Overall, we observed an increase in flavonoid production in response to Hg stress, and Hg tolerance appeared to be associated to a characteristic glycosylation pattern in roots, associated with the accumulation of aglycones and monoglycosylated flavonoids. The findings are discussed in the context of the flavonoid biosynthetic pathway to provide a better understanding of the role of these secondary metabolites in the response and tolerance to Hg stress in M. truncatula.

Enabling Medicago truncatula forward genetics: identification of genetic crossing partner for R108 and development of mapping resources for Tnt1 mutants

Though Medicago truncatula Tnt1 mutants are widely used by researchers in the legume community, they are mainly used for reverse genetics because of the availability of the BLAST-searchable large-scale flanking sequence tags database. However, these mutants should have also been used extensively for forward genetic screens, an effort that has been hindered due to the lack of a compatible genetic crossing partner for the M. truncatula genotype R108, from which Tnt1 mutants were generated. In this study, we selected three Medicago HapMap lines (HM017, HM018 and HM022) and performed reciprocal genetic crosses with R108. After phenotypic analyses in F1 and F2 progenies, HM017 was identified as a compatible crossing partner with R108. By comparing the assembled genomic sequences of HM017 and R108, we developed and confirmed 318 Indel markers evenly distributed across the eight chromosomes of the M. truncatula genome. To validate the effectiveness of these markers, by employing the map-based cloning approach, we cloned the causative gene in the dwarf mutant crs isolated from the Tnt1 mutant population, identifying it as gibberellin 3-β-dioxygenase 1, using some of the confirmed Indel markers. The primer sequences and the size difference of each marker were made available for users in the web-based database. The identification of the crossing partner for R108 and the generation of Indel markers will enhance the forward genetics and the overall usage of the Tnt1 mutants.