Topological data analysis reveals a core gene expression backbone that defines form and function across flowering plants.

Since they emerged approximately 125 million years ago, flowering plants have evolved to dominate the terrestrial landscape and survive in the most inhospitable environments on earth. At their core, these adaptations have been shaped by changes in numerous, interconnected pathways and genes that collectively give rise to emergent biological phenomena. Linking gene expression to morphological outcomes remains a grand challenge in biology, and new approaches are needed to begin to address this gap. Here, we implemented topological data analysis (TDA) to summarize the high dimensionality and noisiness of gene expression data using lens functions that delineate plant tissue and stress responses. Using this framework, we created a topological representation of the shape of gene expression across plant evolution, development, and environment for the phylogenetically diverse flowering plants. The TDA-based Mapper graphs form a well-defined gradient of tissues from leaves to seeds, or from healthy to stressed samples, depending on the lens function. This suggests that there are distinct and conserved expression patterns across angiosperms that delineate different tissue types or responses to biotic and abiotic stresses. Genes that correlate with the tissue lens function are enriched in central processes such as photosynthetic, growth and development, housekeeping, or stress responses. Together, our results highlight the power of TDA for analyzing complex biological data and reveal a core expression backbone that defines plant form and function.

Gene expression profiling of soaked dry beans (Phaseolus vulgaris L.) reveals cell wall modification plays a role in cooking time.

Dry beans (Phaseolus vulgaris L.) are a nutritious food, but their lengthy cooking requirements are barriers to consumption. Presoaking is one strategy to reduce cooking time. Soaking allows hydration to occur prior to cooking, and enzymatic changes to pectic polysaccharides also occur during soaking that shorten the cooking time of beans. Little is known about how gene expression during soaking influences cooking times. The objectives of this study were to (1) identify gene expression patterns that are altered by soaking and (2) compare gene expression in fast-cooking and slow-cooking bean genotypes. RNA was extracted from four bean genotypes at five soaking time points (0, 3, 6, 12, and 18 h) and expression abundances were detected using Quant-seq. Differential gene expression analysis and weighted gene coexpression network analysis were used to identify candidate genes within quantitative trait loci for water uptake and cooking time. Genes related to cell wall growth and development as well as hypoxic stress were differentially expressed between the fast- and slow-cooking beans due to soaking. Candidate genes identified in the slow-cooking beans included enzymes that increase intracellular calcium concentrations and cell wall modification enzymes. The expression of cell wall-strengthening enzymes in the slow-cooking beans may increase their cooking time and ability to resist osmotic stress by preventing cell separation and water uptake in the cotyledon.