Asian soybean rust (ASR), caused by the obligate biotrophic fungus, Phakopsora pachyrhizi , was first reported in the continental United States of America (USA) in 2004 and over the years has been of concern to soybean production in the USA. The prevailing hypothesis is that P. pachyrhizi spores were introduced into the USA via hurricanes originating from South America, particularly Hurricane Ivan. To investigate the genetic diversity and global population structure of P. pachyrhizi , we employed exome-capture based sequencing on 84 field isolates collected from different geographic regions worldwide. We compared the gene-encoding regions from all these field isolates and found that four major haplotypes are prevalent worldwide. Here, we provide genetic evidence supporting multiple incursions that have led to the currently established P. pachyrhizi population of the USA. Phylogenetic analysis of mitochondrial genes further supports this hypothesis. Notably, we observed limited genetic diversity in P. pachyrhizi populations in Brazil, suggesting a clonal population structure in that country that contrasts to populations from the USA and Africa. This study provides the first comprehensive characterization of P. pachyrhizi population structures defined by genetic evidence from populations across major soybean growing regions. ### Competing Interest Statement The authors have declared no competing interest.

Asian soybean rust (ASR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, was first reported in the continental United States of America (USA) in 2004 and over the years has been of concern to soybean production in the United States. The prevailing hypothesis is that P. pachyrhizi spores were introduced into the United States via hurricanes originating from South America, particularly hurricane Ivan. To investigate the genetic diversity and global population structure of P. pachyrhizi, we employed exome-capture based sequencing on 84 field isolates collected from different geographic regions worldwide. We compared the gene-encoding regions from all these field isolates and found that four major mitochondrial haplotypes are prevalent worldwide. Here, we provide genetic evidence supporting multiple incursions that have led to the currently established P. pachyrhizi population of the United States. Phylogenetic analysis of mitochondrial genes further supports this hypothesis. We observed limited genetic diversity in P. pachyrhizi populations across different geographic regions, suggesting a clonal population structure. Additionally, this study is the first to report the F129L mutation in the Cytb gene outside South America, which is associated with strobilurin tolerance. This study provides the first comprehensive characterisation of P. pachyrhizi population structures defined by genetic evidence from populations across major soybean-growing regions.

The biological mechanisms responsible for correct shape acquisition at the apex of the female reproductive organ, the gynoecium, remain poorly understood, despite its fundamental importance for successful plant reproduction and seed production. This process involves a rare bilateral-to-radial symmetry transition in Arabidopsis thaliana, orchestrated in part by the transcription factor SPATULA (SPT). Here, we show that SPT negatively regulates cell cycle regulators CYCLIN-P3;1 (CYCP3;1) and CYCP3;2, which are antagonistically promoted by the hormone cytokinin (CK), to control the radial style morphology by orchestrating a coherent feed-forward loop that converges on the spatial regulation of CYCP3;1 and CYCP3;2 expression. Overexpression of CYCP3s disrupts style radial symmetry, causing the split-style phenotype and hypersensitivity to CK observed in the spt mutant. Finally, we demonstrate a genetic link connecting the machinery of cell division orientation, controlled by auxin, with the cell proliferation input induced by CK, which adds robustness to the apical fusion of the carpels.

We examined the molecular basis of triazole resistance in Blumeria graminis f. sp. tritici (wheat mildew, Bgt), a model organism among powdery mildews. Four genetic models for responses to triazole fungicides were identified among US and UK isolates, involving multiple genetic mechanisms. Firstly, only two amino acid substitutions in CYP51B lanosterol demethylase, the target of triazoles, were associated with resistance, Y136F and S509T (homologous to Y137F and S524T in the reference fungus Zymoseptoria tritici). As sequence variation did not explain the wide range of resistance, we also investigated Cyp51B copy number and expression, the latter using both reverse transcription-quantitative PCR and RNA-seq. The second model for resistance involved higher copy number and expression in isolates with a resistance allele; thirdly, however, moderate resistance was associated with higher copy number of wild-type Cyp51B in some US isolates. A fourth mechanism was heteroallelism with multiple alleles of Cyp51B. UK isolates, with significantly higher mean resistance than their US counterparts, had higher mean copy number, a high frequency of the S509T substitution, which was absent from the United States, and in the most resistant isolates, heteroallelism involving both sensitivity residues Y136+S509 and resistance residues F136+T509. Some US isolates were heteroallelic for Y136+S509 and F136+S509, but this was not associated with higher resistance. The obligate biotrophy of Bgt may constrain the tertiary structure and thus the sequence of CYP51B, so other variation that increases resistance may have a selective advantage. We describe a process by which heteroallelism may be adaptive when Bgt is intermittently exposed to triazoles.

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

ABSTRACTPea,Pisum sativum, is an excellent model system through which Gregor Mendel established the foundational principles of inheritance. Surprisingly, till today, the molecular nature of the genetic differences underlying the seven pairs of contrasting traits that Mendel studied in detail remains partially understood. Here, we present a genomic and phenotypic variation map, coupled with haplotype-phenotype association analyses across a wide range of traits in a globalPisumdiversity panel. We focus on a genomics-enabled genetic dissection of each of the seven traits Mendel studied, revealing many previously undescribed alleles for the four characterized genes,R,Le,IandA, and elucidating the gene identities and mutations for the remaining three uncharacterized traits. Notably, we identify: (1) a ca. 100kb deletion upstream of theChlorophyll synthase(ChlG) gene, which generates aberrant transcripts and confers the yellow pod phenotype ofgpmutants; (2) an in-frame premature stop codon mutation in a Dodeca-CLE41/44 signalling peptide which explains the parchmentless mutant phenotype corresponding top; and (3) a 5bp in-frame deletion in aCIK-likereceptor kinase gene corresponding to the fasciated stem phenotypefa, which Mendel described in terms of flower position, and we postulate the existence of aModifier of fa(Mfa) locus that masks this meristem defect. Mendel noted the pleiotropy of theamutation, including inhibition of axil ring anthocyanin pigmentation, a trait we found to be controlled by allelic variants of the geneDwithin anR2R3-MYBgene cluster. Furthermore, we characterize and validate natural variation of a quantitative genetic locus governing both pod width and seed weight, characters that Mendel deemed were not sufficiently demarcated for his analyses. This study establishes a cornerstone for fundamental research, education in biology and genetics, and pea breeding practices.

Grasspea (Lathyrus sativus L.) is an underutilised but promising legume crop with tolerance to a wide range of abiotic and biotic stress factors, and potential for climate-resilient agriculture. Despite a long history and wide geographical distribution of cultivation, only limited breeding resources are available. This paper reports a 5.96 Gbp genome assembly of grasspea genotype LS007, of which 5.03 Gbp is scaffolded into 7 pseudo-chromosomes. The assembly has a BUSCO completeness score of 99.1% and is annotated with 31719 gene models and repeat elements. This represents the most contiguous and accurate assembly of the grasspea genome to date.

Polycomb repressive complex 2 (PRC2) mediates epigenetic silencing of target genes in animals and plants. In Arabidopsis, PRC2 is required for the cold-induced epigenetic silencing of the FLC floral repressor locus to align flowering with spring. During this process, PRC2 relies on VEL accessory factors, including the constitutively expressed VRN5 and the cold-induced VIN3. The VEL proteins are physically associated with PRC2, but their individual functions remain unclear. Here, we show an intimate association between recombinant VRN5 and multiple components within a reconstituted PRC2, dependent on a compact conformation of VRN5 central domains. Key residues mediating this compact conformation are conserved among VRN5 orthologs across the plant kingdom. In contrast, VIN3 interacts with VAL1, a transcriptional repressor that binds directly to FLC These associations differentially affect their role in H3K27me deposition: Both proteins are required for H3K27me3, but only VRN5 is necessary for H3K27me2. Although originally defined as vernalization regulators, VIN3 and VRN5 coassociate with many targets in the Arabidopsis genome that are modified with H3K27me3. Our work therefore reveals the distinct accessory roles for VEL proteins in conferring cold-induced silencing on FLC, with broad relevance for PRC2 targets generally.

Grass pea (Lathyrus sativus L.) is a rich source of protein cultivated as an insurance crop in Ethiopia, Eritrea, India, Bangladesh, and Nepal. Its resilience to both drought and flooding makes it a promising crop for ensuring food security in a changing climate. The lack of genetic resources and the crop's association with the disease neurolathyrism have limited the cultivation of grass pea. Here, we present an annotated, long read-based assembly of the 6.5 Gbp L. sativus genome. Using this genome sequence, we have elucidated the biosynthetic pathway leading to the formation of the neurotoxin, β-L-oxalyl-2,3-diaminopropionic acid (β-L-ODAP). The final reaction of the pathway depends on an interaction between L. sativus acyl-activating enzyme 3 (LsAAE3) and a BAHD-acyltransferase (LsBOS) that form a metabolon activated by CoA to produce β-L-ODAP. This provides valuable insight into the best approaches for developing varieties which produce substantially less toxin.

TenA thiamin-degrading enzymes are commonly found in prokaryotes, plants, fungi and algae and are involved in the thiamin salvage pathway. The gut symbiont Bacteroides thetaiotaomicron (Bt) produces a TenA protein (BtTenA) which is packaged into its extracellular vesicles. An alignment of BtTenA protein sequence with proteins from different databases using the basic local alignment search tool (BLAST) and the generation of a phylogenetic tree revealed that BtTenA is related to TenA-like proteins not only found in a small number of intestinal bacterial species but also in some aquatic bacteria, aquatic invertebrates, and freshwater fish. This is, to our knowledge, the first report describing the presence of TenA-encoding genes in the genome of members of the animal kingdom. By searching metagenomic databases of diverse host-associated microbial communities, we found that BtTenA homologues were mostly represented in biofilms present on the surface of macroalgae found in Australian coral reefs. We also confirmed the ability of a recombinant BtTenA to degrade thiamin. Our study shows that Bt tenA -like genes which encode a novel sub-class of TenA proteins are sparingly distributed across two kingdoms of life, a feature of accessory genes known for their ability to spread between species through horizontal gene transfer.

Cellular RNAs are heterogeneous with respect to their alternative processing and secondary structures, but the functional importance of this complexity is still poorly understood. A set of alternatively processed antisense non-coding transcripts, which are collectively called COOLAIR , are generated at the Arabidopsis floral-repressor locus FLOWERING LOCUS C ( FLC ) 1 . Different isoforms of COOLAIR influence FLC transcriptional output in warm and cold conditions 2 – 7 . Here, to further investigate the function of COOLAIR , we developed an RNA structure-profiling method to determine the in vivo structure of single RNA molecules rather than the RNA population average. This revealed that individual isoforms of the COOLAIR transcript adopt multiple structures with different conformational dynamics. The major distally polyadenylated COOLAIR isoform in warm conditions adopts three predominant structural conformations, the proportions and conformations of which change after cold exposure. An alternatively spliced, strongly cold-upregulated distal COOLAIR isoform 6 shows high structural diversity, in contrast to proximally polyadenylated COOLAIR . A hyper-variable COOLAIR structural element was identified that was complementary to the FLC transcription start site. Mutations altering the structure of this region changed FLC expression and flowering time, consistent with an important regulatory role of the COOLAIR structure in FLC transcription. Our work demonstrates that isoforms of non-coding RNA transcripts adopt multiple distinct and functionally relevant structural conformations, which change in abundance and shape in response to external conditions.

Nucleotide composition is suggested to infer gene functionality and ecological adaptation of species to distinct environments. However, the underlying biological function of nucleotide composition dictating environmental adaptations is largely unknown. Here, we systematically analyze the nucleotide composition of transcriptomes across 1000 plants (1KP) and their corresponding habitats. Intriguingly, we find that plants growing in cold climates have guanine (G)-enriched transcriptomes, which are prone to forming RNA G-quadruplex structures. Both immunofluorescence detection and in vivo structure profiling reveal that RNA G-quadruplex formation in plants is globally enhanced in response to cold. Cold-responsive RNA G-quadruplexes strongly enhanced mRNA stability, rather than affecting translation. Disruption of individual RNA G-quadruplex promotes mRNA decay in the cold, leading to impaired plant cold response. Therefore, we propose that plants adopted RNA G-quadruplex structure as a molecular signature to facilitate their adaptation to the cold during evolution.

Biparental recombinant inbred line (RIL) populations are sets of genetically stable lines and have a simple population structure that facilitates the dissection of the genetics of interesting traits. On the other hand, populations derived from multiparent intercrosses combine both greater diversity and higher numbers of recombination events than RILs. Here, we describe a simple population structure: a three-way recombinant inbred population combination. This structure was easy to produce and was a compromise between biparental and multiparent populations. We show that this structure had advantages when analyzing cultivar crosses, and could achieve a mapping resolution of a few genes.

Aegilops tauschii , the diploid wild progenitor of the D subgenome of bread wheat, is a reservoir of genetic diversity for improving bread wheat performance and environmental resilience. Here we sequenced 242 Ae. tauschii accessions and compared them to the wheat D subgenome to characterize genomic diversity. We found that a rare lineage of Ae. tauschii geographically restricted to present-day Georgia contributed to the wheat D subgenome in the independent hybridizations that gave rise to modern bread wheat. Through k -mer-based association mapping, we identified discrete genomic regions with candidate genes for disease and pest resistance and demonstrated their functional transfer into wheat by transgenesis and wide crossing, including the generation of a library of hexaploids incorporating diverse Ae. tauschii genomes. Exploiting the genomic diversity of the Ae. tauschii ancestral diploid genome permits rapid trait discovery and functional genetic validation in a hexaploid background amenable to breeding.

Agricultural soil harbors a diverse microbiome that can form beneficial relationships with plants, including the inhibition of plant pathogens. Pseudomonas spp. are one of the most abundant bacterial genera in the soil and rhizosphere and play important roles in promoting plant health. However, the genetic determinants of this beneficial activity are only partially understood. Here, we genetically and phenotypically characterize the Pseudomonas fluorescens population in a commercial potato field, where we identify strong correlations between specialized metabolite biosynthesis and antagonism of the potato pathogens Streptomyces scabies and Phytophthora infestans . Genetic and chemical analyses identified hydrogen cyanide and cyclic lipopeptides as key specialized metabolites associated with S. scabies inhibition, which was supported by in planta biocontrol experiments. We show that a single potato field contains a hugely diverse and dynamic population of Pseudomonas bacteria, whose capacity to produce specialized metabolites is shaped both by plant colonization and defined environmental inputs.

BACKGROUND:Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry.RESULTS:Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat.CONCLUSIONS:We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Background mRNA processing is critical for gene expression. A challenge in regulating mRNA processing is how to recognize the actual mRNA processing sites, such as splice and polyadenylation sites, when the sequence content is insufficient for this purpose. Previous studies suggested that RNA structure affects mRNA processing. However, the regulatory role of RNA structure in mRNA processing remains unclear. Results Here, we perform in vivo selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) chemical profiling on Arabidopsis and generate the in vivo nuclear RNA structure landscape. We find that nuclear mRNAs fold differently from cytosolic mRNAs across translation start and stop sites. Notably, we discover a two-nucleotide single-stranded RNA structure feature upstream of 5′ splice sites that is strongly associated with splicing and the selection of alternative 5′ splice sites. The regulatory role of this RNA structure feature is further confirmed by experimental validation. Moreover, we find the single-strandedness of branch sites is also associated with 3′ splice site recognition. We also identify an RNA structure feature comprising two close-by single-stranded regions that is specifically associated with both polyadenylation and alternative polyadenylation events. Conclusions We successfully identify pre-mRNA structure features associated with splicing and polyadenylation at whole-genome scale and validate an RNA structure feature which can regulate splicing. Our study unveils a new RNA structure regulatory mechanism for mRNA processing.

Aegilops tauschii, the diploid wild progenitor of the D-subgenome of bread wheat, constitutes a reservoir of genetic diversity for improving bread wheat performance and environmental resilience. To better define and understand this diversity, we sequenced 242 Ae. tauschii accessions and compared them to the wheat D-subgenome. We characterized a rare, geographically-restricted lineage of Ae. tauschii and discovered that it contributed to the wheat D-subgenome, thereby elucidating the origin of bread wheat from at least two independent hybridizations. We then used k -mer-based association mapping to identify discrete genomic regions with candidate genes for disease and pest resistance and demonstrated their functional transfer into wheat by transgenesis and wide crossing, including the generation of a library of ‘synthetic’ hexaploids incorporating diverse Ae. tauschii genomes. This pipeline permits rapid trait discovery in the diploid ancestor through to functional genetic validation in a hexaploid background amenable to breeding.

Rice blast is a devastating disease caused by the fungal pathogen Magnaporthe oryzae that threatens rice production around the world. The fungus produces a specialized infection cell, called the appressorium, that enables penetration through the plant cell wall in response to surface signals from the rice leaf. The underlying biology of plant infection, including the regulation of appressorium formation, is not completely understood. Here we report the identification of a network of temporally coregulated transcription factors that act downstream of the Pmk1 mitogen-activated protein kinase pathway to regulate gene expression during appressorium-mediated plant infection. We show that this tiered regulatory mechanism involves Pmk1-dependent phosphorylation of the Hox7 homeobox transcription factor, which regulates genes associated with induction of major physiological changes required for appressorium development—including cell-cycle control, autophagic cell death, turgor generation and melanin biosynthesis—as well as controlling a additional set of virulence-associated transcription factor–encoding genes. Pmk1-dependent phosphorylation of Mst12 then regulates gene functions involved in septin-dependent cytoskeletal re-organization, polarized exocytosis and effector gene expression, which are necessary for plant tissue invasion. Identification of this regulatory cascade provides new potential targets for disease intervention. Comparative RNA-seq, ChIP–seq and quantitative phosphoproteomics reveal how the blast fungus uses the Pmk1 MAP kinase to regulate a network of transcription factors that orchestrate the complex transcriptome changes necessary for infecting rice plants.