Nooshin Omranian

Identification and Mode of Inheritance of Quantitative Trait Loci for Secondary Metabolite Abundance in Tomato

Metabolite QTL for secondary metabolism in a Solanum pennelli introgression line population show different modes of inheritance and network properties linking the traits. A large-scale metabolic quantitative trait loci (mQTL) analysis was performed on the well-characterized Solanum pennellii introgression lines to investigate the genomic regions associated with secondary metabolism in tomato fruit pericarp. In total, 679 mQTLs were detected across the 76 introgression lines. Heritability analyses revealed that mQTLs of secondary metabolism were less affected by environment than mQTLs of primary metabolism. Network analysis allowed us to assess the interconnectivity of primary and secondary metabolism as well as to compare and contrast their respective associations with morphological traits. Additionally, we applied a recently established real-time quantitative PCR platform to gain insight into transcriptional control mechanisms of a subset of the mQTLs, including those for hydroxycinnamates, acyl-sugar, naringenin chalcone, and a range of glycoalkaloids. Intriguingly, many of these compounds displayed a dominant-negative mode of inheritance, which is contrary to the conventional wisdom that secondary metabolite contents decreased on domestication. We additionally performed an exemplary evaluation of two candidate genes for glycolalkaloid mQTLs via the use of virus-induced gene silencing. The combined data of this study were compared with previous results on primary metabolism obtained from the same material and to other studies of natural variance of secondary metabolism.

Genetic Determinants of the Network of Primary Metabolism and Their Relationships to Plant Performance in a Maize Recombinant Inbred Line Population

Elucidation of the genetic determinants of maize primary metabolism and a metabolite-metabolite-agronomic trait network will promote efficient use of metabolites in maize improvement. Deciphering the influence of genetics on primary metabolism in plants will provide insights useful for genetic improvement and enhance our fundamental understanding of plant growth and development. Although maize (Zea mays) is a major crop for food and feed worldwide, the genetic architecture of its primary metabolism is largely unknown. Here, we use high-density linkage mapping to dissect large-scale metabolic traits measured in three different tissues (leaf at seedling stage, leaf at reproductive stage, and kernel at 15 d after pollination [DAP]) of a maize recombinant inbred line population. We identify 297 quantitative trait loci (QTLs) with moderate (86.2% of the mapped QTL, R2 = 2.4 to 15%) to major effects (13.8% of the mapped QTL, R2 >15%) for 79 primary metabolites across three tissues. Pairwise epistatic interactions between these identified loci are detected for more than 25.9% metabolites explaining 6.6% of the phenotypic variance on average (ranging between 1.7 and 16.6%), which implies that epistasis may play an important role for some metabolites. Key candidate genes are highlighted and mapped to carbohydrate metabolism, the tricarboxylic acid cycle, and several important amino acid biosynthetic and catabolic pathways, with two of them being further validated using candidate gene association and expression profiling analysis. Our results reveal a metabolite-metabolite-agronomic trait network that, together with the genetic determinants of maize primary metabolism identified herein, promotes efficient utilization of metabolites in maize improvement.

Gene regulatory network inference using fused LASSO on multiple data sets

Devising computational methods to accurately reconstruct gene regulatory networks given gene expression data is key to systems biology applications. Here we propose a method for reconstructing gene regulatory networks by simultaneous consideration of data sets from different perturbation experiments and corresponding controls. The method imposes three biologically meaningful constraints: (1) expression levels of each gene should be explained by the expression levels of a small number of transcription factor coding genes, (2) networks inferred from different data sets should be similar with respect to the type and number of regulatory interactions, and (3) relationships between genes which exhibit similar differential behavior over the considered perturbations should be favored. We demonstrate that these constraints can be transformed in a fused LASSO formulation for the proposed method. The comparative analysis on transcriptomics time-series data from prokaryotic species, Escherichia coli and Mycobacterium tuberculosis, as well as a eukaryotic species, mouse, demonstrated that the proposed method has the advantages of the most recent approaches for regulatory network inference, while obtaining better performance and assigning higher scores to the true regulatory links. The study indicates that the combination of sparse regression techniques with other biologically meaningful constraints is a promising framework for gene regulatory network reconstructions.

Effect of salt stress on genes encoding translation-associated proteins in Arabidopsis thaliana.

Salinity negatively affects plant growth and disturbs chloroplast integrity. Here, we aimed at identifying salt-responsive translation-related genes in Arabidopsis thaliana with an emphasis on those encoding plastid-located proteins. We used quantitative real-time PCR to test the expression of 170 genes after short-term salt stress (up to 24 h) and identified several genes affected by the stress including: PRPL11, encoding plastid ribosomal protein L11, ATAB2, encoding a chloroplast-located RNA-binding protein presumably functioning as an activator of translation, and PDF1B, encoding a peptide deformylase involved in N-formyl group removal from nascent proteins synthesized in chloroplasts. These genes were previously shown to have important functions in chloroplast biology and may therefore represent new targets for biotechnological optimization of salinity tolerance.

Segmentation of biological multivariate time-series data.

Time-series data from multicomponent systems capture the dynamics of the ongoing processes and reflect the interactions between the components. The progression of processes in such systems usually involves check-points and events at which the relationships between the components are altered in response to stimuli. Detecting these events together with the implicated components can help understand the temporal aspects of complex biological systems. Here we propose a regularized regression-based approach for identifying breakpoints and corresponding segments from multivariate time-series data. In combination with techniques from clustering, the approach also allows estimating the significance of the determined breakpoints as well as the key components implicated in the emergence of the breakpoints. Comparative analysis with the existing alternatives demonstrates the power of the approach to identify biologically meaningful breakpoints in diverse time-resolved transcriptomics data sets from the yeast Saccharomyces cerevisiae and the diatom Thalassiosira pseudonana.

Differential metabolic and coexpression networks of plant metabolism

Recent analyses have demonstrated that plant metabolic networks do not differ in their structural properties and that genes involved in basic metabolic processes show smaller coexpression than genes involved in specialized metabolism. By contrast, our analysis reveals differences in the structure of plant metabolic networks and patterns of coexpression for genes in (non)specialized metabolism. Here we caution that conclusions concerning the organization of plant metabolism based on network-driven analyses strongly depend on the computational approaches used.

Network-based segmentation of biological multivariate time series.

Molecular phenotyping technologies (e.g., transcriptomics, proteomics, and metabolomics) offer the possibility to simultaneously obtain multivariate time series (MTS) data from different levels of information processing and metabolic conversions in biological systems. As a result, MTS data capture the dynamics of biochemical processes and components whose couplings may involve different scales and exhibit temporal changes. Therefore, it is important to develop methods for determining the time segments in MTS data, which may correspond to critical biochemical events reflected in the coupling of the system’s components. Here we provide a novel network-based formalization of the MTS segmentation problem based on temporal dependencies and the covariance structure of the data. We demonstrate that the problem of partitioning MTS data into segments to maximize a distance function, operating on polynomially computable network properties, often used in analysis of biological network, can be efficiently solved. To enable biological interpretation, we also propose a breakpoint-penalty (BP-penalty) formulation for determining MTS segmentation which combines a distance function with the number/length of segments. Our empirical analyses of synthetic benchmark data as well as time-resolved transcriptomics data from the metabolic and cell cycles of Saccharomyces cerevisiae demonstrate that the proposed method accurately infers the phases in the temporal compartmentalization of biological processes. In addition, through comparison on the same data sets, we show that the results from the proposed formalization of the MTS segmentation problem match biological knowledge and provide more rigorous statistical support in comparison to the contending state-of-the-art methods.

A network property necessary for concentration robustness

Maintenance of functionality of complex cellular networks and entire organisms exposed to environmental perturbations often depends on concentration robustness of the underlying components. Yet, the reasons and consequences of concentration robustness in large-scale cellular networks remain largely unknown. Here, we derive a necessary condition for concentration robustness based only on the structure of networks endowed with mass action kinetics. The structural condition can be used to design targeted experiments to study concentration robustness. We show that metabolites satisfying the necessary condition are present in metabolic networks from diverse species, suggesting prevalence of this property across kingdoms of life. We also demonstrate that our predictions about concentration robustness of energy-related metabolites are in line with experimental evidence from Escherichia coli. The necessary condition is applicable to mass action biological systems of arbitrary size, and will enable understanding the implications of concentration robustness in genetic engineering strategies and medical applications.

AtRsgA from Arabidopsis thaliana is important for maturation of the small subunit of the chloroplast ribosome

Plastid ribosomes are very similar in structure and function to the ribosomes of their bacterial ancestors. Since ribosome biogenesis is not thermodynamically favorable under biological conditions it requires the activity of many assembly factors. Here we have characterized a homolog of bacterial RsgA in Arabidopsis thaliana and show that it can complement the bacterial homolog. Functional characterization of a strong mutant in Arabidopsis revealed that the protein is essential for plant viability, while a weak mutant produced dwarf, chlorotic plants that incorporated immature pre-16S ribosomal RNA into translating ribosomes. Physiological analysis of the mutant plants revealed smaller, but more numerous, chloroplasts in the mesophyll cells, reduction of chlorophyll a and b, depletion of proplastids from the rib meristem and decreased photosynthetic electron transport rate and efficiency. Comparative RNA sequencing and proteomic analysis of the weak mutant and wild-type plants revealed that various biotic stress-related, transcriptional regulation and post-transcriptional modification pathways were repressed in the mutant. Intriguingly, while nuclear- and chloroplast-encoded photosynthesis-related proteins were less abundant in the mutant, the corresponding transcripts were increased, suggesting an elaborate compensatory mechanism, potentially via differentially active retrograde signaling pathways. To conclude, this study reveals a chloroplast ribosome assembly factor and outlines the transcriptomic and proteomic responses of the compensatory mechanism activated during decreased chloroplast function.

AtRsgA from Arabidopsis thaliana controls maturation of the small subunit of the chloroplast ribosome

Summary Plastid ribosomes are very similar in structure and function to ribosomes of their bacterial ancestors. Since ribosome biogenesis is not thermodynamically favourable at biological conditions, it requires activity of many assembly factors. Here, we have characterized a homolog of bacterial rsgA in Arabidopsis thaliana and show that it can complement the bacterial homolog. Functional characterization of a strong mutant in Arabidopsis revealed that the protein is essential for plant viability, while a weak mutant produced dwarf, chlorotic plants that incorporated immature pre-16S ribosomal RNA into translating ribosomes. Physiological analysis of the mutant plants revealed smaller, but more numerous chloroplasts in the mesophyll cells, reduction of chlorophyll a and b, depletion of proplastids from the rib meristem and decreased photosynthetic electron transport rate and efficiency. Comparative RNA-sequencing and proteomic analysis of the weak mutant and wild-type plants revealed that various biotic stress-related, transcriptional regulation and post-transcriptional modification pathways were repressed in the mutant. Intriguingly, while nuclear- and chloroplast-encoded photosynthesis-related proteins were less abundant in the mutant, the corresponding transcripts were upregulated, suggesting an elaborate compensatory mechanism, potentially via differentially active retrograde signalling pathways. To conclude, this study reveals a new chloroplast ribosome assembly factor and outlines the transcriptomic and proteomic responses of the compensatory mechanism activated during decreased chloroplast function. Significance statement AtRsgA is an assembly factor necessary for maturation of the small subunit of the chloroplast ribosome. Depletion of AtRsgA leads to dwarfed, chlorotic plants and smaller, but more numerous chloroplasts. Large-scale transcriptomic and proteomic analysis revealed that chloroplast-encoded and - targeted proteins were less abundant, while the corresponding transcripts were upregulated in the mutant. We analyse the transcriptional responses of several retrograde signalling pathways to suggest a mechanism underlying this compensatory response.