A pan-grass transcriptome reveals patterns of cellular divergence in crops.

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

From buds to shoots: insights into grapevine development from the Witch's Broom bud sport.

BACKGROUND: Bud sports occur spontaneously in plants when new growth exhibits a distinct phenotype from the rest of the parent plant. The Witch's Broom bud sport occurs occasionally in various grapevine (Vitis vinifera) varieties and displays a suite of developmental defects, including dwarf features and reduced fertility. While it is highly detrimental for grapevine growers, it also serves as a useful tool for studying grapevine development. We used the Witch's Broom bud sport in grapevine to understand the developmental trajectories of the bud sports, as well as the potential genetic basis. We analyzed the phenotypes of two independent cases of the Witch's Broom bud sport, in the Dakapo and Merlot varieties of grapevine, alongside wild type counterparts. To do so, we quantified various shoot traits, performed 3D X-ray Computed Tomography on dormant buds, and landmarked leaves from the samples. We also performed Illumina and Oxford Nanopore sequencing on the samples and called genetic variants using these sequencing datasets. RESULTS: The Dakapo and Merlot cases of Witch's Broom displayed severe developmental defects, with no fruit/clusters formed and dwarf vegetative features. However, the Dakapo and Merlot cases of Witch's Broom studied were also phenotypically different from one another, with distinct differences in bud and leaf development. We identified 968-974 unique genetic mutations in our two Witch's Broom cases that are potential causal variants of the bud sports. Examining gene function and validating these genetic candidates through PCR and Sanger-sequencing revealed one strong candidate mutation in Merlot Witch's Broom impacting the gene GSVIVG01008260001. CONCLUSIONS: The Witch's Broom bud sports in both varieties studied had dwarf phenotypes, but the two instances studied were also vastly different from one another and likely have distinct genetic bases. Future work on Witch's Broom bud sports in grapevine could provide more insight into development and the genetic pathways involved in grapevine.

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 epigenomic 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 a characteristic of more recent single-gene duplicates. Core angiosperm gene families were 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.

Rhythmic lipid and gene expression responses to chilling in panicoid grasses.

Chilling stress threatens plant growth and development, particularly affecting membrane fluidity and cellular integrity. Understanding plant membrane responses to chilling stress is important for unraveling the molecular mechanisms of stress tolerance. Whereas core transcriptional responses to chilling stress and stress tolerance are conserved across species, the associated changes in membrane lipids appear to be less conserved, as which lipids are affected by chilling stress varies by species. Here, we investigated changes in gene expression and membrane lipids in response to chilling stress during one 24 hour cycle in chilling-tolerant foxtail millet (Setaria italica), and chilling-sensitive sorghum (Sorghum bicolor), and Urochloa (browntop signal grass, Urochloa fusca, lipids only), leveraging their evolutionary relatedness and differing levels of chilling-stress tolerance. We show that most chilling-induced lipid changes are conserved across the three species, while we observed distinct, time-specific responses in chilling-tolerant foxtail millet, indicating the presence of a finely orchestrated adaptive mechanism. We detected rhythmicity in lipid responses to chilling stress in the three grasses, which were also present in Arabidopsis (Arabidopsis thaliana), suggesting the conservation of rhythmic patterns across species and highlighting the importance of accounting for time of day. When integrating lipid datasets with gene expression profiles, we identified potential candidate genes that showed corresponding transcriptional changes in response to chilling stress, providing insights into the differences in regulatory mechanisms between chilling-sensitive sorghum and chilling-tolerant foxtail millet.

Characterization of huntingtin interactomes and their dynamic responses in living cells by proximity proteomics.

Huntingtin (Htt) is a large protein without clearly defined molecular functions. Mutation in this protein causes Huntington's disease (HD), a fatal inherited neurodegenerative disorder. Identification of Htt-interacting proteins by the traditional approaches including yeast two-hybrid systems and affinity purifications has greatly facilitated the understanding of Htt function. However, these methods eliminated the intracellular spatial information of the Htt interactome during sample preparations. Moreover, the temporal changes of the Htt interactome in response to acute cellular stresses cannot be easily resolved with these approaches. Ascorbate peroxidase (APEX2)-based proximity labeling has been used to spatiotemporally investigate protein-protein interactions in living cells. In this study, we generated stable human SH-SY5Y cell lines expressing full-length Htt23Q and Htt145Q with N-terminus tagged Flag-APEX2 to quantitatively map the spatiotemporal changes of Htt interactome to a mild acute proteotoxic stress. Our data revealed that normal and mutant Htt (muHtt) are associated with distinct intracellular microenvironments. Specifically, mutant Htt is preferentially associated with intermediate filaments and myosin complexes. Furthermore, the dynamic changes of Htt interactomes in response to stress are different between normal and mutant Htt. Vimentin is identified as one of the most significant proteins that preferentially interacts with muHtt in situ. Further functional studies demonstrated that mutant Htt affects the vimentin's function of regulating proteostasis in healthy and HD human neural stem cells. Taken together, our data offer important insights into the molecular functions of normal and mutant Htt by providing a list of Htt-interacting proteins in their natural cellular context for further studies in different HD models.

Treatment of Low-grade Intermediate-risk Nonmuscle-invasive Bladder Cancer With UGN-102 ± Transurethral Resection of Bladder Tumor Compared to Transurethral Resection of Bladder Tumor Monotherapy: A Randomized, Controlled, Phase 3 Trial (ATLAS).

PURPOSE: Low-grade intermediate-risk nonmuscle-invasive bladder cancer is a chronic illness commonly treated by repetitive transurethral resection of bladder tumor. We compared the efficacy and safety of intravesical chemoablation with UGN-102 (a reverse thermal gel containing mitomycin), with or without subsequent transurethral resection of bladder tumor, to transurethral resection of bladder tumor alone in patients with low-grade intermediate-risk nonmuscle-invasive bladder cancer. MATERIALS AND METHODS: This prospective, randomized, phase 3 trial recruited patients with new or recurrent low-grade intermediate-risk nonmuscle-invasive bladder cancer to receive initial treatment with either UGN-102 once weekly for 6 weeks or transurethral resection of bladder tumor. Patients were followed quarterly by endoscopy, cytology, and for-cause biopsy. The primary end point was disease-free survival. All patients were followed for adverse events. RESULTS: Trial enrollment was halted by the sponsor to pursue an alternative development strategy after 282 of a planned 632 patients were randomized to UGN-102 ± subsequent transurethral resection of bladder tumor (n=142) or transurethral resection of bladder tumor monotherapy (n=140), rendering the trial underpowered to perform hypothesis testing. Patients were predominantly male and ≥65 years of age. Tumor-free complete response 3 months after initial treatment was achieved by 92 patients (65%) who received UGN-102 and 89 patients (64%) treated by transurethral resection of bladder tumor. The estimated probability of disease-free survival 15 months after randomization was 72% for UGN-102 ± transurethral resection of bladder tumor and 50% for transurethral resection of bladder tumor (hazard ratio 0.45). The most common adverse events (incidence ≥10%) in the UGN-102 group were dysuria, micturition urgency, nocturia, and pollakiuria. CONCLUSIONS: Primary, nonsurgical chemoablation with UGN-102 for the management of low-grade intermediate-risk nonmuscle-invasive bladder cancer offers a potential therapeutic alternative to immediate transurethral resection of bladder tumor monotherapy and warrants further investigation.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs.

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Primary Chemoablation of Low-Grade Intermediate-Risk Nonmuscle-Invasive Bladder Cancer Using UGN-102, a Mitomycin-Containing Reverse Thermal Gel (Optima II): A Phase 2b, Open-Label, Single-Arm Trial.

PURPOSE: Low-grade intermediate-risk nonmuscle-invasive bladder cancer (LG IR NMIBC) is a recurrent disease, thus requiring repeated transurethral resection of bladder tumor under general anesthesia. We evaluated the efficacy and safety of UGN-102, a mitomycin-containing reverse thermal gel, as a primary chemoablative therapeutic alternative to transurethral resection of bladder tumor for patients with LG IR NMIBC. MATERIALS AND METHODS: This prospective, phase 2b, open-label, single-arm trial recruited patients with biopsy-proven LG IR NMIBC to receive 6 once-weekly instillations of UGN-102. The primary end point was complete response (CR) rate, defined as the proportion of patients with negative endoscopic examination, negative cytology and negative for-cause biopsy 3 months after treatment initiation. Patients with CR were followed quarterly up to 12 months to assess durability of treatment effect. Safety and adverse events were monitored throughout the trial. RESULTS: A total of 63 patients (38 males and 25 females 33-96 years old) enrolled and received ≥1 instillation of UGN-102. Among the patients 41 (65%) achieved CR at 3 months, of whom 39 (95%), 30 (73%) and 25 (61%) remained disease-free at 6, 9 and 12 months after treatment initiation, respectively. A total of 13 patients had documented recurrences. The probability of durable response 9 months after CR (12 months after treatment initiation) was estimated to be 73% by Kaplan-Meier analysis. Common adverse events (incidence ≥10%) included dysuria, urinary frequency, hematuria, micturition urgency, urinary tract infection and fatigue. CONCLUSIONS: Nonsurgical primary chemoablation of LG IR NMIBC using UGN-102 resulted in significant treatment response with sustained durability. UGN-102 may provide an alternative to repetitive surgery for patients with LG IR NMIBC.

Parallels between natural selection in the cold-adapted crop-wild relative Tripsacum dactyloides and artificial selection in temperate adapted maize.

Artificial selection has produced varieties of domesticated maize that thrive in temperate climates around the world. However, the direct progenitor of maize, teosinte, is indigenous only to a relatively small range of tropical and subtropical latitudes and grows poorly or not at all outside of this region. Tripsacum, a sister genus to maize and teosinte, is naturally endemic to the majority of areas in the western hemisphere where maize is cultivated. A full-length reference transcriptome for Tripsacum dactyloides generated using long-read Iso-Seq data was used to characterize independent adaptation to temperate climates in this clade. Genes related to phospholipid biosynthesis, a critical component of cold acclimation in other cold-adapted plant lineages, were enriched among those genes experiencing more rapid rates of protein sequence evolution in T. dactyloides. In contrast with previous studies of parallel selection, we find that there is a significant overlap between the genes that were targets of artificial selection during the adaptation of maize to temperate climates and those that were targets of natural selection in temperate-adapted T. dactyloides. Genes related to growth, development, response to stimulus, signaling, and organelles were enriched in the set of genes identified as both targets of natural and artificial selection.

An epigenetic breeding system in soybean for increased yield and stability.

Epigenetic variation has been associated with a wide range of adaptive phenotypes in plants, but there exist few direct means for exploiting this variation. RNAi suppression of the plant-specific gene, MutS HOMOLOG1 (MSH1), in multiple plant species produces a range of developmental changes accompanied by modulation of defence, phytohormone and abiotic stress response pathways along with methylome repatterning. This msh1-conditioned developmental reprogramming is retained independent of transgene segregation, giving rise to transgene-null 'memory' effects. An isogenic memory line crossed to wild type produces progeny families displaying increased variation in adaptive traits that respond to selection. This study investigates amenability of the MSH1 system for inducing agronomically valuable epigenetic variation in soybean. We developed MSH1 epi-populations by crossing with msh1-acquired soybean memory lines. Derived soybean epi-lines showed increase in variance for multiple yield-related traits including pods per plant, seed weight and maturity time in both glasshouse and field trials. Selected epi-F2:4 and epi-F2:5 lines showed an increase in seed yield over wild type. By epi-F2:6, we observed a return of MSH1-derived enhanced growth back to wild-type levels. Epi-populations also showed evidence of reduced epitype-by-environment (e × E) interaction, indicating higher yield stability. Transcript profiling of epi-lines identified putative signatures of enhanced growth behaviour across generations. Genes related to cell cycle, abscisic acid biosynthesis and auxin response, particularly SMALL AUXIN UP RNAs (SAURs), were differentially expressed in epi-F2:4 lines that showed increased yield when compared to epi-F2:6 . These data support the potential of MSH1-derived epigenetic variation in plant breeding for enhanced yield and yield stability.

Stress-responsive pathways and small RNA changes distinguish variable developmental phenotypes caused by MSH1 loss.

BACKGROUND: Proper regulation of nuclear-encoded, organelle-targeted genes is crucial for plastid and mitochondrial function. Among these genes, MutS Homolog 1 (MSH1) is notable for generating an assortment of mutant phenotypes with varying degrees of penetrance and pleiotropy. Stronger phenotypes have been connected to stress tolerance and epigenetic changes, and in Arabidopsis T-DNA mutants, two generations of homozygosity with the msh1 insertion are required before severe phenotypes begin to emerge. These observations prompted us to examine how msh1 mutants contrast according to generation and phenotype by profiling their respective transcriptomes and small RNA populations. RESULTS: Using RNA-seq, we analyze pathways that are associated with MSH1 loss, including abiotic stresses such as cold response, pathogen defense and immune response, salicylic acid, MAPK signaling, and circadian rhythm. Subtle redox and environment-responsive changes also begin in the first generation, in the absence of strong phenotypes. Using small RNA-seq we further identify miRNA changes, and uncover siRNA trends that indicate modifications at the chromatin organization level. In all cases, the magnitude of changes among protein-coding genes, transposable elements, and small RNAs increases according to generation and phenotypic severity. CONCLUSION: Loss of MSH1 is sufficient to cause large-scale regulatory changes in pathways that have been individually linked to one another, but rarely described all together within a single mutant background. This study enforces the recognition of organelles as critical integrators of both internal and external cues, and highlights the relationship between organelle and nuclear regulation in fundamental aspects of plant development and stress signaling. Our findings also encourage further investigation into potential connections between organelle state and genome regulation vis-á-vis small RNA feedback.

Regulatory dynamics distinguishing desiccation tolerance strategies within resurrection grasses.

Desiccation tolerance has evolved recurrently in grasses using two unique strategies of either protecting or dismantling the photosynthetic apparatus to minimize photooxidative damage under life without water (anhydrobiosis). Here, we surveyed chromatin architecture and gene expression during desiccation in two closely related grasses with distinguishing desiccation tolerance strategies to identify regulatory dynamics underlying these unique adaptations. In both grasses, we observed a strong association between nearby chromatin accessibility and gene expression in desiccated tissues compared to well-watered, reflecting an unusual chromatin stability under anhydrobiosis. Integration of chromatin accessibility (ATACseq) and expression data (RNAseq) revealed a core desiccation response across these two grasses. This includes many genes with binding sites for the core seed development transcription factor ABI5, supporting the long-standing hypothesis that vegetative desiccation tolerance evolved from rewiring seed pathways. Oropetium thomaeum has a unique set of desiccation induced genes and regulatory elements associated with photoprotection, pigment biosynthesis, and response to high light, reflecting its adaptation of protecting the photosynthetic apparatus under desiccation (homoiochlorophyly). By contrast, Eragrostis nindensis has unique accessible and expressed genes related to chlorophyll catabolism, scavenging of amino acids, and hypoxia, highlighting its poikilochlorophyllous adaptations of dismantling the photosynthetic apparatus and degrading chlorophyll under desiccation. Together, our results highlight the complex regulatory and expression dynamics underlying desiccation tolerance in grasses.