Search for: All records

Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene’s expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.

 

Dry beans (Phaseolus vulgarisL.) 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.

 

Tepary bean (Phaseolus acutifoliusA. Gray), indigenous to the arid climates of northern Mexico and the Southwest United States, diverged from common bean (Phaseolus vulgarisL.), approximately 2 million years ago and exhibits a wide range of resistance to biotic stressors. The tepary genome is highly syntenic to the common bean genome providing a foundation for discovery and breeding of agronomic traits between these two crop species. Although a limited number of adaptive traits from tepary bean have been introgressed into common bean, hybridization barriers between these two species required the development of bridging lines to alleviate this barrier. Thus, to fully utilize the extant tepary bean germplasm as both a crop and as a donor of adaptive traits, we developed a diversity panel of 422 cultivated, weedy, and wild tepary bean accessions which were then genotyped and phenotyped to enable population genetic analyses and genome‐wide association studies for their response to a range of biotic stressors. Population structure analyses of the panel revealed eight subpopulations and the differentiation of botanical varieties withinP. acutifolius. Genome‐wide association studies revealed loci and candidate genes underlying biotic stress resistance including quantitative trait loci for resistance to weevils, common bacterial blight, Fusarium wilt, and bean common mosaic necrosis virus that can be harnessed not only for tepary bean but also common bean improvement.

 

Common bean (Phaseolus vulgarisL.) is a nutrient-rich food, but its long cooking times hinder its wider utilization. The Yellow Bean Collection (YBC) was assembled with 295 genotypes from global sources to assess the genetic and phenotypic diversity for end-use quality traits in yellow beans. The panel was genotyped with over 2,000 SNPs identified via Genotyping-By-Sequencing (GBS). Through population structure analyses with the GBS markers, the YBC was determined to be 69% Andean, 26% Middle American, and 5% admixture. The YBC was grown in two major bean production regions in the U.S., Michigan (MI) and Nebraska (NE) over two years. The genotypes exhibited a wide diversity in days to flower, seed weight, water uptake, and cooking time. The cooking times of the YBC ranged from 17–123 min. The cooking time were longer and varied more widely in NE with many more genotypes exhibiting hardshell than in MI. Fast-cooking genotypes were identified with various yellow colors; 20 genotypes cooked within 20 min in MI, and eight genotypes cooked within 31 min in NE. Water uptake and cooking time were significantly affected by the environment, which included both the growing and cooking environment, and notably in relation to cooking, NE is higher elevation than MI. SNPs associated with cooking time were identified with genome-wide association analyses and a polygalacturonase gene on Pv04 was considered to be a candidate gene. The genotypic and phenotypic variability, fast-cooking genotypes, and the associated SNPs of the YBC will lay the foundation for utilizing yellow beans for breeding and genetic analyses.

 

Leaf morphology is dynamic, continuously deforming during leaf expansion and among leaves within a shoot. Here, we measured the leaf morphology of more than 200 grapevines (Vitisspp.) over four years and modeled changes in leaf shape along the shoot to determine whether a composite leaf shape comprising all the leaves from a single shoot can better capture the variation and predict species identity compared with individual leaves.

Using homologous universal landmarks found in grapevine leaves, we modeled various morphological features as polynomial functions of leaf nodes. The resulting functions were used to reconstruct modeled leaf shapes across the shoots, generating composite leaves that comprehensively capture the spectrum of leaf morphologies present.

We found that composite leaves are better predictors of species identity than individual leaves from the same plant. We were able to use composite leaves to predict the species identity of previously unassigned grapevines, which were verified with genotyping.

Observations of individual leaf shape fail to capture the true diversity between species. Composite leaf shape—an assemblage of modeled leaf snapshots across the shoot—is a better representation of the dynamic and essential shapes of leaves, in addition to serving as a better predictor of species identity than individual leaves.

 

The regulation of gene expression is central to many biological processes. Gene regulatory networks (GRNs) link transcription factors (TFs) to their target genes and represent maps of potential transcriptional regulation. Here, we analyzed a large number of publically available maize (Zea mays) transcriptome data sets including >6000 RNA sequencing samples to generate 45 coexpression-based GRNs that represent potential regulatory relationships between TFs and other genes in different populations of samples (cross-tissue, cross-genotype, and tissue-and-genotype samples). While these networks are all enriched for biologically relevant interactions, different networks capture distinct TF-target associations and biological processes. By examining the power of our coexpression-based GRNs to accurately predict covarying TF-target relationships in natural variation data sets, we found that presence/absence changes rather than quantitative changes in TF gene expression are more likely associated with changes in target gene expression. Integrating information from our TF-target predictions and previous expression quantitative trait loci (eQTL) mapping results provided support for 68 TFs underlying 74 previously identified trans-eQTL hotspots spanning a variety of metabolic pathways. This study highlights the utility of developing multiple GRNs within a species to detect putative regulators of important plant pathways and provides potential targets for breeding or biotechnological applications.

 

Abstract The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris , Plectranthus barbatus , and Leonotis leonurus . We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites. 

Abstract Motivation The accurate prediction of complex phenotypes such as metabolic fluxes in living systems is a grand challenge for systems biology and central to efficiently identifying biotechnological interventions that can address pressing industrial needs. The application of gene expression data to improve the accuracy of metabolic flux predictions using mechanistic modeling methods such as flux balance analysis (FBA) has not been previously demonstrated in multi-tissue systems, despite their biotechnological importance. We hypothesized that a method for generating metabolic flux predictions informed by relative expression levels between tissues would improve prediction accuracy. Results Relative gene expression levels derived from multiple transcriptomic and proteomic datasets were integrated into FBA predictions of a multi-tissue, diel model of Arabidopsis thaliana’s central metabolism. This integration dramatically improved the agreement of flux predictions with experimentally based flux maps from 13C metabolic flux analysis compared with a standard parsimonious FBA approach. Disagreement between FBA predictions and MFA flux maps was measured using weighted averaged percent error values, and for parsimonious FBA this was169%–180% for high light conditions and 94%–103% for low light conditions, depending on the gene expression dataset used. This fell to 10%-13% and 9%-11% upon incorporating expression data into the modeling process, which also substantially altered the predicted carbon and energy economy of the plant. Availability and implementation Code and data generated as part of this study are available from https://github.com/Gibberella/ArabidopsisGeneExpressionWeights. 

Tar spot is a devasting corn disease caused by the obligate fungal pathogen Phyllachora maydis. Since its initial identification in the United States in 2015, P. maydis has become an increasing threat to corn production. Despite this, P. maydis has remained largely understudied at the molecular level due to difficulties surrounding its obligate lifestyle. Here, we generated a significantly improved P. maydis nuclear and mitochondrial genome using a combination of long- and short-read technologies and also provide the first transcriptomic analysis of primary tar spot lesions. Our results show that P. maydis is deficient in inorganic nitrogen utilization, is likely heterothallic, and encodes for significantly more protein coding genes, including secreted enzymes and effectors, than previous determined. Furthermore, our expression analysis suggests that following primary tar spot lesion formation, P. maydis might reroute carbon flux away from DNA replication and cell division pathways and towards pathways previously implicated in having significant roles in pathogenicity, such as autophagy and secretion. Together, our results identified several highly expressed unique secreted factors that likely contribute to host recognition and subsequent infection, greatly increasing our knowledge of the biological capacity of P. maydis, which have much broader implications for mitigating tar spot of corn.