Functional compensation of primary and secondary metabolites by duplicate genes in Arabidopsis thaliana

Expression profiles in knock-down transgenic plants of high and low diversified duplicate genes in Arabidopsis thaliana

Duplicated genes show various degrees of functional diversification in plants. We previously identified 1,052 pairs of high diversified duplicates (HDDs) and 600 pairs of low diversified duplicates (LDDs) in Arabidopsis thaliana. Single knock-down of HDDs induced abnormal phenotypic changes because the other gene copy could not compensate for the knock-down effect, while single knock-down of LDDs did not induce abnormal phenotypic changes because of functional compensation by the copy gene. Here, focusing on one pair each of HDDs and LDDs, we performed transcriptome analyses in single-knock-down transgenic plants. The numbers of differentially expressed genes in single-knock-down transgenic plants were not different between HDDs and LDDs. Thus, functional compensation inferred by transcriptomics was similar between HDDs and LDDs. However, the trend of differentially expressed genes was similar in the pair of LDDs, while expression profiles were dissimilar in the pair of HDDs. This result indicates that a pair of LDDs tends to share similar functions but a pair of HDDs tends to have undergone functional divergence. Taking these findings together, as the reason for no phenotypic changes in single knock-down of LDDs but phenotypic changes in double knock-down of LDDs, we concluded that phenotypic changes of LDDs were induced by decreasing gene dosage.

Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

AGAMOUS regulates various target genes via cell cycle-coupled H3K27me3 dilution in floral meristems and stamens

The MADS domain transcription factor AGAMOUS (AG) regulates floral meristem termination by preventing maintenance of the histone modification lysine 27 of histone H3 (H3K27me3) along the KNUCKLES (KNU) coding sequence. At 2 d after AG binding, cell division has diluted the repressive mark H3K27me3, allowing activation of KNU transcription prior to floral meristem termination. However, how many other downstream genes are temporally regulated by this intrinsic epigenetic timer and what their functions are remain unknown. Here, we identify direct AG targets regulated through cell cycle-coupled H3K27me3 dilution in Arabidopsis thaliana. Expression of the targets KNU, AT HOOK MOTIF NUCLEAR LOCALIZED PROTEIN18 (AHL18), and PLATZ10 occurred later in plants with longer H3K27me3-marked regions. We established a mathematical model to predict timing of gene expression and manipulated temporal gene expression using the H3K27me3-marked del region from the KNU coding sequence. Increasing the number of del copies delayed and reduced KNU expression in a polycomb repressive complex 2- and cell cycle-dependent manner. Furthermore, AHL18 was specifically expressed in stamens and caused developmental defects when misexpressed. Finally, AHL18 bound to genes important for stamen growth. Our results suggest that AG controls the timing of expression of various target genes via cell cycle-coupled dilution of H3K27me3 for proper floral meristem termination and stamen development.

The heavy metal-associated isoprenylated plant protein (HIPP) gene family plays a crucial role in cadmium resistance and accumulation in the tea plant (Camellia sinensis L.)

Heavy metal-associated isoprenylated plant proteins (HIPPs) are only distributed in vascular plants, and are essential for the detoxification and vascular transport of heavy metals in plants. However, the HIPP gene family has not been thoroughly explored in the tea plant (Camellia sinensis). In this study, we systematically identified 56C. sinensis CsHIPP genes from five groups and characterized their phylogeny, structures, and the features of the encoded proteins. The expression patterns of CsHIPP genes in various tissues of C. sinensis were investigated based on a previous RNA-seq data analysis. The expression patterns of CsHIPP genes were explored in cadmium (Cd)-treated C. sinensis roots using our RNA-seq data. Three CsHIPP genes (CsHIPP22, CsHIPP24, and CsHIPP36) with high expression levels in Cd-treated C. sinensis roots were selected as candidate genes associated with Cd tolerance. Overexpression of CsHIPP22, CsHIPP24, and CsHIPP36 in a yeast mutant (ycf1) rescued Cd-sensitive ycf1 yeast and increased the yeast resistance to Cd stress, implying that these three CsHIPPs might be involved in Cd tolerance. These findings will enable the roles of HIPPs in Cd absorption and detoxification to be better understood as well as improving our understanding of the Cd-resistance and Cd-accumulation mechanisms in tea plant.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

Combining GS-assisted GWAS and transcriptome analysis to mine candidate genes for nitrogen utilization efficiency in Populus cathayana

BACKGROUND

Forest trees such as poplar, shrub willow, et al. are essential natural resources for sustainable and renewable energy production, and their wood can reduce dependence on fossil fuels and reduce environmental pollution. However, the productivity of forest trees is often limited by the availability of nitrogen (N), improving nitrogen use efficiency (NUE) is an important way to address it. Currently, NUE genetic resources are scarce in forest tree research, and more genetic resources are urgently needed.

RESULTS

Here, we performed genome-wide association studies (GWAS) using the mixed linear model (MLM) to identify genetic loci regulating growth traits in Populus cathayana at two N levels, and attempted to enhance the signal strength of single nucleotide polymorphism (SNP) detection by performing genome selection (GS) assistance GWAS. The results of the two GWAS analyses identified 55 and 40 SNPs that were respectively associated with plant height (PH) and ground diameter (GD), and 92 and 69 candidate genes, including 30 overlapping genes. The prediction accuracy of the GS model (rrBLUP) for phenotype exceeds 0.9. Transcriptome analysis of 13 genotypes under two N levels showed that genes related to carbon and N metabolism, amino acid metabolism, energy metabolism, and signal transduction were differentially expressed in the xylem of P. cathayana under N treatment. Furthermore, we observed strong regional patterns in gene expression levels of P. cathayana, with significant differences between different regions. Among them, P. cathayana in Longquan region exhibited the highest response to N. Finally, through weighted gene co-expression network analysis (WGCNA), we identified a module closely related to the N metabolic process and eight hub genes.

CONCLUSIONS

Integrating the GWAS, RNA-seq and WGCNA data, we ultimately identified four key regulatory genes (PtrNAC123, PtrNAC025, Potri.002G233100, and Potri.006G236200) involved in the wood formation process, and they may affect P. cathayana growth and wood formation by regulating nitrogen metabolism. This study will provide strong evidence for N regulation mechanisms, and reliable genetic resources for growth and NUE genetic improvement in poplar.

Genetic Analyses of Seed Longevity in Capsicum annuum L. in Cold Storage Conditions

Seed longevity is the most important trait in the genebank management system. No seed can remain infinitely viable. There are 1241 accessions of Capsicum annuum L. available at the German Federal ex situ genebank at IPK Gatersleben. C. annuum (Capsicum) is the most economically important species of the genus Capsicum. So far, there is no report that has addressed the genetic basis of seed longevity in Capsicum. Here, we convened a total of 1152 Capsicum accessions that were deposited in Gatersleben over forty years (from 1976 to 2017) and assessed their longevity by analyzing the standard germination percentage after 5-40 years of storage at -15/-18 °C. These data were used to determine the genetic causes of seed longevity, along with 23,462 single nucleotide polymorphism (SNP) markers covering all of the 12 Capsicum chromosomes. Using the association-mapping approach, we identified a total of 224 marker trait associations (MTAs) (34, 25, 31, 35, 39, 7, 21 and 32 MTAs after 5-, 10-, 15-, 20-, 25-, 30-, 35- and 40-year storage intervals) on all the Capsicum chromosomes. Several candidate genes were identified using the blast analysis of SNPs, and these candidate genes are discussed.

Global transcriptome and coexpression network analyses reveal cultivar-specific molecular signatures associated with different rooting depth responses to drought stress in potato

Potato is one of the most important vegetable crops worldwide. Its growth, development and ultimately yield is hindered by drought stress condition. Breeding and selection of deep-rooted and drought-tolerant potato varieties has become a prime approach for improving the yield and quality of potato (Solanum tuberosum L.) in arid and semiarid areas. A comprehensive understanding of root development-related genes has enabled scientists to formulate strategies to incorporate them into breeding to improve complex agronomic traits and provide opportunities for the development of stress tolerant germplasm. Root response to drought stress is an intricate process regulated through complex transcriptional regulatory network. To understand the rooting depth and molecular mechanism, regulating root response to drought stress in potato, transcriptome dynamics of roots at different stages of drought stress were analyzed in deep (C119) and shallow-rooted (C16) cultivars. Stage-specific expression was observed for a significant proportion of genes in each cultivar and it was inferred that as compared to C16 (shallow-rooted), approximately half of the genes were differentially expressed in deep-rooted cultivar (C119). In C16 and C119, 11 and 14 coexpressed gene modules, respectively, were significantly associated with physiological traits under drought stress. In a comparative analysis, some modules were different between the two cultivars and were associated with differential response to specific drought stress stage. Transcriptional regulatory networks were constructed, and key components determining rooting depth were identified. Through the results, we found that rooting depth (shallow vs deep) was largely determined by plant-type, cell wall organization or biogenesis, hemicellulose metabolic process, and polysaccharide metabolic process. In addition, candidate genes responding to drought stress were identified in deep (C119) and shallow (C16) rooted potato varieties. The results of this study will be a valuable source for further investigations on the role of candidate gene(s) that affect rooting depth and drought tolerance mechanisms in potato.

Functional compensation of mouse duplicates by their paralogs expressed in the same tissues

Analyses in a number of organisms have shown that duplicated genes are less likely to be essential than singletons. This implies that genes can often compensate for the loss of their paralogs. However, it is unclear why the loss of some duplicates can be compensated by their paralogs, whereas the loss of other duplicates cannot. Surprisingly, initial analyses in mice did not detect differences in the essentiality of duplicates and singletons. Only subsequent analyses, using larger gene knockout datasets and controlling for a number of confounding factors, did detect significant differences. Previous studies have not taken into account the tissues in which duplicates are expressed. We hypothesized that in complex organisms, in order for a gene's loss to be compensated by one or more of its paralogs, such paralogs need to be expressed in at least the same set of tissues as the lost gene. To test our hypothesis, we classified mouse duplicates into two categories based on the expression patterns of their paralogs: "compensable duplicates" (those with paralogs expressed in all the tissues in which the gene is expressed) and "non-compensable duplicates" (those whose paralogs are not expressed in all the tissues where the gene is expressed). In agreement with our hypothesis, the essentiality of non-compensable duplicates is similar to that of singletons, whereas compensable duplicates exhibit a substantially lower essentiality. Our results imply that duplicates can often compensate for the loss of their paralogs, but only if they are expressed in the same tissues. Indeed, the compensation ability is more dependent on expression patterns than on protein sequence similarity. The existence of these two kinds of duplicates with different essentialities, which has been overlooked by prior studies, may have hindered the detection of differences between singletons and duplicates.

A Simple Evolutionary Model of Genetic Robustness After Gene Duplication

When a dispensable gene is duplicated (referred to the ancestral dispensability denoted by O⁺), genetic buffering and duplicate compensation together maintain the duplicate redundancy, whereas duplicate compensation is the only mechanism when an essential gene is duplicated (referred to the ancestral essentiality denoted by O^-). To investigate these evolutionary scenarios of genetic robustness, I formulated a simple mixture model for analyzing duplicate pairs with one of the following states: double dispensable (DD), semi-dispensable (one dispensable one essential, DE), or double essential (EE). This model was applied to the yeast duplicate pairs from a whole-genome duplication (WGD) occurred about 100 million years ago (mya), and the mouse duplicate pairs from a WGD occurred about more than 500 mya. Both case studies revealed that the proportion of essentiality for those duplicates with ancestral essentiality [P_E(O^-)] was much higher than that for those with ancestral dispensability [P_E(O⁺)]. While it was negligible in the yeast duplicate pairs, P_E(O⁺) (about 20%) was shown statistically significant in the mouse duplicate pairs. These findings, together, support the hypothesis that both sub-functionalization and neo-functionalization may play some roles after gene duplication, though the former may be much faster than the later.

Dynamic evolution of small signalling peptide compensation in plant stem cell control

Gene duplications are a hallmark of plant genome evolution and a foundation for genetic interactions that shape phenotypic diversity^1-5. Compensation is a major form of paralogue interaction^6-8 but how compensation relationships change as allelic variation accumulates is unknown. Here we leveraged genomics and genome editing across the Solanaceae family to capture the evolution of compensating paralogues. Mutations in the stem cell regulator CLV3 cause floral organs to overproliferate in many plants^9-11. In tomato, this phenotype is partially suppressed by transcriptional upregulation of a closely related paralogue¹². Tobacco lost this paralogue, resulting in no compensation and extreme clv3 phenotypes. Strikingly, the paralogues of petunia and groundcherry nearly completely suppress clv3, indicating a potent ancestral state of compensation. Cross-species transgenic complementation analyses show that this potent compensation partially degenerated in tomato due to a single amino acid change in the paralogue and cis-regulatory variation that limits its transcriptional upregulation. Our findings show how genetic interactions are remodelled following duplications and suggest that dynamic paralogue evolution is widespread over short time scales and impacts phenotypic variation from natural and engineered mutations.

A collaborative classroom investigation of the evolution of SABATH methyltransferase substrate preference shifts over 120 million years of flowering plant history

Next-generation sequencing has resulted in an explosion of available data, much of which remains unstudied in terms of biochemical function; yet, experimental characterization of these sequences has the potential to provide unprecedented insight into the evolution of enzyme activity. One way to make inroads into the experimental study of the voluminous data available is to engage students by integrating teaching and research in a college classroom such that eventually hundreds or thousands of enzymes may be characterized. In this study, we capitalize on this potential to focus on SABATH methyltransferase enzymes that have been shown to methylate the important plant hormone, salicylic acid (SA), to form methyl salicylate. We analyze data from 76 enzymes of flowering plant species in 23 orders and 41 families to investigate how widely conserved substrate preference is for SA methyltransferase orthologs. We find a high degree of conservation of substrate preference for SA over the structurally similar metabolite, benzoic acid, with recent switches that appear to be associated with gene duplication and at least three cases of functional compensation by paralogous enzymes. The presence of Met in active site position 150 is a useful predictor of SA methylation preference in SABATH methyltransferases but enzymes with other residues in the homologous position show the same substrate preference. Although our dense and systematic sampling of SABATH enzymes across angiosperms has revealed novel insights, this is merely the “tip of the iceberg” since thousands of sequences remain uncharacterized in this enzyme family alone.

Plant Metabolic Network 15: A resource of genome-wide metabolism databases for 126 plants and algae

To understand and engineer plant metabolism, we need a comprehensive and accurate annotation of all metabolic information across plant species. As a step towards this goal, we generated genome-scale metabolic pathway databases of 126 algal and plant genomes, ranging from model organisms to crops to medicinal plants (https://plantcyc.org). Of these, 104 have not been reported before. We systematically evaluated the quality of the databases, which revealed that our semi-automated validation pipeline dramatically improves the quality. We then compared the metabolic content across the 126 organisms using multiple correspondence analysis and found that Brassicaceae, Poaceae, and Chlorophyta appeared as metabolically distinct groups. To demonstrate the utility of this resource, we used recently published sorghum transcriptomics data to discover previously unreported trends of metabolism underlying drought tolerance. We also used single-cell transcriptomics data from the Arabidopsis root to infer cell type-specific metabolic pathways. This work shows the quality and quantity of our resource and demonstrates its wide-ranging utility in integrating metabolism with other areas of plant biology.

Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.

MicroRNAs Roles in Plants Secondary Metabolism

Plant growth and development is dependent on the regulation of classes of microRNAs (miRNAs) that have emerged as important gene regulators. These miRNAs can regulate plant gene expression to function. They play an important roles in biological homeostasis and environmental response controls. A wide range of plant biological and metabolic processes, including developmental timing, tissues specific development, and differentiation, depends on miRNAs. They perpetually regulate secondary metabolite functions in different plant family lines. Mapping of molecular phylogenies shows the distribution of secondary metabolism in the plant territory. More importantly, a lot of information related to miRNA regulatory processes in plants is revealed, but the role of miRNAs in secondary metabolism regulation and functions of the metabolites are still unclear. In this review, we pinnacle some potential miRNAs regulating the secondary metabolite biosynthesis activities in plants. This will provide an alternative knowledge for functional studies of secondary metabolism.

Genome-wide search and structural and functional analyses for late embryogenesis-abundant (LEA) gene family in poplar

BACKGROUND

The Late Embryogenesis-Abundant (LEA) gene families, which play significant roles in regulation of tolerance to abiotic stresses, widely exist in higher plants. Poplar is a tree species that has important ecological and economic values. But systematic studies on the gene family have not been reported yet in poplar.

RESULTS

On the basis of genome-wide search, we identified 88 LEA genes from Populus trichocarpa and renamed them as PtrLEA. The PtrLEA genes have fewer introns, and their promoters contain more cis-regulatory elements related to abiotic stress tolerance. Our results from comparative genomics indicated that the PtrLEA genes are conserved and homologous to related genes in other species, such as Eucalyptus robusta, Solanum lycopersicum and Arabidopsis. Using RNA-Seq data collected from poplar under two conditions (with and without salt treatment), we detected 24, 22 and 19 differentially expressed genes (DEGs) in roots, stems and leaves, respectively. Then we performed spatiotemporal expression analysis of the four up-regulated DEGs shared by the tissues, constructed gene co-expression-based networks, and investigated gene function annotations.

CONCLUSION

Lines of evidence indicated that the PtrLEA genes play significant roles in poplar growth and development, as well as in responses to salt stress.

Characterization of STP4 promoter in Indian mustard Brassica juncea for use as an aphid responsive promoter

OBJECTIVE

Brassica juncea, a major oilseed crop, suffers substantial yield losses due to infestation by mustard aphids (Lipaphis erysimi). Unavailability of resistance genes within the accessible gene pool underpins significance of the transgenic strategy in developing aphid resistance. In this study, we aimed for the identification of an aphid-responsive promoter from B. juncea, based on the available genomic resources.

RESULTS

A monosaccharide transporter gene, STP4 in B. juncea was activated by aphids and sustained increased expression as the aphids colonized the plants. We cloned the upstream intergenic region of STP4 and validated its stand-alone aphid-responsive promoter activity. Further, deletion analysis identified the putative cis-elements important for the aphid responsive promoter activity.

CONCLUSION

The identified STP4 promoter can potentially be used for driving high level aphid-inducible expression of transgenes in plants. Use of aphid-responsive promoter instead of constitutive promoters can potentially reduce the metabolic burden of transgene-expression on the host plant.

Influence of low O2 and high CO2 environment on changes in metabolite concentrations in harvested vegetable soybeans

Influence of atmosphere and storage period on the physicochemical and biological properties of harvested vegetable soybeans stored for 10 d at 25 °C was investigated. Storing vegetable soybeans under modified atmosphere (low O2 and high CO2), was more effective in maintaining its green color and mass than storing them under normoxia. Principal component 1 (PC1; contribution rate: 25%) was related to the atmospheres, whereas PC2 (contribution rate: 19%) was related to storage period. Cluster analysis showed that some types of sugars decreased, whereas some types of organic and amino acids increased with deterioration. Alanine, an indicator of low O2 stress, was maintained for 3 d under modified atmospheres, whereas alanine significantly decreased under normoxia. The concentrations of inositol and niacinamide (functional ingredients) under the modified atmospheres were significantly higher than those under normoxia. Thus, storage under modified atmospheres was effective in maintaining freshness and increasing the nutritional content of vegetable soybeans.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

Rgma-Induced Neo1 Proteolysis Promotes Neural Tube Morphogenesis

Neuroepithelial cell (NEC) elongation is one of several key cell behaviors that mediate the tissue-level morphogenetic movements that shape the neural tube (NT), the precursor of the brain and spinal cord. However, the upstream signals that promote NEC elongation have been difficult to tease apart from those regulating apico-basal polarity and hingepoint formation, due to their confounding interdependence. The Repulsive Guidance Molecule a (Rgma)/Neogenin 1 (Neo1) signaling pathway plays a conserved role in NT formation (neurulation) and is reported to regulate both NEC elongation and apico-basal polarity, through signal transduction events that have not been identified. We examine here the role of Rgma/Neo1 signaling in zebrafish (sex unknown), an organism that does not use hingepoints to shape its hindbrain, thereby enabling a direct assessment of the role of this pathway in NEC elongation. We confirm that Rgma/Neo1 signaling is required for microtubule-mediated NEC elongation, and demonstrate via cell transplantation that Neo1 functions cell autonomously to promote elongation. However, in contrast to previous findings, our data do not support a role for this pathway in establishing apical junctional complexes. Last, we provide evidence that Rgma promotes Neo1 glycosylation and intramembrane proteolysis, resulting in the production of a transient, nuclear intracellular fragment (NeoICD). Partial rescue of Neo1a and Rgma knockdown embryos by overexpressing neoICD suggests that this proteolytic cleavage is essential for neurulation. Based on these observations, we propose that RGMA-induced NEO1 proteolysis orchestrates NT morphogenesis by promoting NEC elongation independently of the establishment of apical junctional complexes.SIGNIFICANCE STATEMENT The neural tube, the CNS precursor, is shaped during neurulation. Neural tube defects occur frequently, yet underlying genetic risk factors are poorly understood. Neuroepithelial cell (NEC) elongation is essential for proper completion of neurulation. Thus, connecting NEC elongation with the molecular pathways that control this process is expected to reveal novel neural tube defect risk factors and increase our understanding of NT development. Effectors of cell elongation include microtubules and microtubule-associated proteins; however, upstream regulators remain controversial due to the confounding interdependence of cell elongation and establishment of apico-basal polarity. Here, we reveal that Rgma-Neo1 signaling controls NEC elongation independently of the establishment of apical junctional complexes and identify Rgma-induced Neo1 proteolytic cleavage as a key upstream signaling event.

Variations of Secondary Metabolites among Natural Populations of Sub-Antarctic Ranunculus Species Suggest Functional Redundancy and Versatility

Plants produce a high diversity of metabolites which help them sustain environmental stresses and are involved in local adaptation. However, shaped by both the genome and the environment, the patterns of variation of the metabolome in nature are difficult to decipher. Few studies have explored the relative parts of geographical region versus environment or phenotype in metabolomic variability within species and none have discussed a possible effect of the region on the correlations between metabolites and environments or phenotypes. In three sub-Antarctic Ranunculus species, we examined the role of region in metabolite differences and in the relationship between individual compounds and environmental conditions or phenotypic traits. Populations of three Ranunculus species were sampled across similar environmental gradients in two distinct geographical regions in îles Kerguelen. Two metabolite classes were studied, amines (quantified by high-performance liquid chromatography and fluorescence spectrophotometry) and flavonols (quantified by ultra-high-performance liquid chromatography with triple quadrupole mass spectrometry). Depending on regions, the same environment or the same trait may be related to different metabolites, suggesting metabolite redundancy within species. In several cases, a given metabolite showed different or even opposite relations with the same environmental condition or the same trait across the two regions, suggesting metabolite versatility within species. Our results suggest that metabolites may be functionally redundant and versatile within species, both in their response to environments and in their relation with the phenotype. These findings open new perspectives for understanding evolutionary responses of plants to environmental changes.

Diurnal changes in concerted plant protein phosphorylation and acetylation in Arabidopsis organs and seedlings

Protein phosphorylation and acetylation are the two most abundant post-translational modifications (PTMs) that regulate protein functions in eukaryotes. In plants, these PTMs have been investigated individually; however, their co-occurrence and dynamics on proteins is currently unknown. Using Arabidopsis thaliana, we quantified changes in protein phosphorylation, acetylation and protein abundance in leaf rosettes, roots, flowers, siliques and seedlings at the end of day (ED) and at the end of night (EN). This identified 2549 phosphorylated and 909 acetylated proteins, of which 1724 phosphorylated and 536 acetylated proteins were also quantified for changes in PTM abundance between ED and EN. Using a sequential dual-PTM workflow, we identified significant PTM changes and intersections in these organs and plant developmental stages. In particular, cellular process-, pathway- and protein-level analyses reveal that the phosphoproteome and acetylome predominantly intersect at the pathway- and cellular process-level at ED versus EN. We found 134 proteins involved in core plant cell processes, such as light harvesting and photosynthesis, translation, metabolism and cellular transport, that were both phosphorylated and acetylated. Our results establish connections between PTM motifs, PTM catalyzing enzymes and putative substrate networks. We also identified PTM motifs for further characterization of the regulatory mechanisms that control cellular processes during the diurnal cycle in different Arabidopsis organs and seedlings. The sequential dual-PTM analysis expands our understanding of diurnal plant cell regulation by PTMs and provides a useful resource for future analyses, while emphasizing the importance of analyzing multiple PTMs simultaneously to elucidate when, where and how they are involved in plant cell regulation.

A Genome-Wide Association Study Identifies New Loci Involved in Wound-Induced Lateral Root Formation in Arabidopsis thaliana

Root systems can display variable architectures that contribute to nutrient foraging or to increase the tolerance of abiotic stress conditions. Root tip excision promotes the developmental progression of previously specified lateral root (LR) founder cells, which allows to easily measuring the branching capacity of a given root as regards its genotype and/or growth conditions. Here, we describe the natural variation among 120 Arabidopsis thaliana accessions in root system architecture (RSA) after root tip excision. Wound-induced changes in RSA were associated with 19 genomic loci using genome-wide association mapping. Three candidate loci associated with wound-induced LR formation were investigated. Sequence variation in the hypothetical protein encoded by the At4g01090 gene affected wound-induced LR development and its loss-of-function mutants displayed a reduced number of LRs after root tip excision. Changes in a histidine phosphotransfer protein putatively involved in cytokinin signaling were significantly associated with LR number variation after root tip excision. Our results provide a better understanding of some of the genetic components involved in LR capacity variation among accessions.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

BACKGROUND

The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention.

RESULTS

We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database.

CONCLUSIONS

We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

Effect of the storage atmosphere on metabolomics of harvested tomatoes (Solanum lycopersicum L.)

Harvested tomatoes were stored under atmospheres that were normoxic, anoxic, or modified (altered O2 and CO2 concentrations). Each atmosphere was created by storing the tomatoes at 25°C for up to 8 days in different kinds of pouches. During storage, metabolites of the tomatoes were measured using metabolomics. We obtained score plots of the metabolites on eighth day of storage by principal component analysis. There was a tendency for groups to be divided on the basis of score plot according to the composition of each gas. PC1 and PC2 seemed to correspond to the influence of O2 and CO2 concentrations, respectively, and the total contribution rate of the two axes was 72%, so that we concluded that the metabolites were affected mainly by O2 and CO2 concentrations. The results indicate that metabolomics may be an effective tool to reveal the relationship between metabolic state of harvested fruits and the atmosphere.

A Highly Specific Genome-Wide Association Study Integrated with Transcriptome Data Reveals the Contribution of Copy Number Variations to Specialized Metabolites in Arabidopsis thaliana Accessions

Lineage-specific gene duplications contribute to a large variation in specialized metabolites among different plant species. There is also considerable variability in the specialized metabolites within a single plant species. However, it is unclear whether copy number variations (CNVs) derived from gene duplication events contribute to the diversity of specialized metabolites within species. We identified metabolome quantitative trait genes (mQTGs) associated with quantitative metabolite variations and examined the relationship between mQTGs and CNVs. We obtained 1,335 specialized metabolite signals from 53 worldwide A. thaliana accessions using liquid chromatography-quadrupole time-of-flight mass spectrometry. In this study, genes associated with specialized metabolites were inferred by either a generally authorized genome-wide association study (GWAS) approach or a novel analysis of the association between gene expression and metabolite accumulation. Genes qualified by both analyses are defined to be mQTGs. The integrated method enabled us to detect mQTGs with a low false positive rate (=5.71 × 10-4). We also identified 5,654 genes associated with 1,335 specialized metabolites. Of these genes, 4.4% were affected by CNVs, which was more than expected (χ2 test: P < 0.01). This result suggests that CNVs contribute to variations in specialized metabolites within a species. To assess the contribution of CNVs to adaptive evolution in A. thaliana, we examined the selective sweeps around the mQTGs. We observed that the mQTGs with CNVs tended to undergo selective sweeps. These observations imply that variations in specialized metabolites caused by CNVs contribute to the adaptive evolution of A. thaliana.

Role of duplicate genes in determining the tissue-selectivity of hereditary diseases

A longstanding puzzle in human genetics is what limits the clinical manifestation of hundreds of hereditary diseases to certain tissues, while their causal genes are expressed throughout the human body. A general conception is that tissue-selective disease phenotypes emerge when masking factors operate in unaffected tissues, but are specifically absent or insufficient in disease-manifesting tissues. Although this conception has critical impact on the understanding of disease manifestation, it was never challenged in a systematic manner across a variety of hereditary diseases and affected tissues. Here, we address this gap in our understanding via rigorous analysis of the susceptibility of over 30 tissues to 112 tissue-selective hereditary diseases. We focused on the roles of paralogs of causal genes, which are presumably capable of compensating for their aberration. We show for the first time at large-scale via quantitative analysis of omics datasets that, preferentially in the disease-manifesting tissues, paralogs are under-expressed relative to causal genes in more than half of the diseases. This was observed for several susceptible tissues and for causal genes with varying number of paralogs, suggesting that imbalanced expression of paralogs increases tissue susceptibility. While for many diseases this imbalance stemmed from up-regulation of the causal gene in the disease-manifesting tissue relative to other tissues, it was often combined with down-regulation of its paralog. Notably in roughly 20% of the cases, this imbalance stemmed only from significant down-regulation of the paralog. Thus, dosage relationships between paralogs appear as important, yet currently under-appreciated, modifiers of disease manifestation.

Functional divergence of duplicate genes several million years after gene duplication in Arabidopsis

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was ∼10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.

Genetic compensation triggered by actin mutation prevents the muscle damage caused by loss of actin protein

The lack of a mutant phenotype in homozygous mutant individuals’ due to compensatory gene expression triggered upstream of protein function has been identified as genetic compensation. Whilst this intriguing process has been recognized in zebrafish, the presence of homozygous loss of function mutations in healthy human individuals suggests that compensation may not be restricted to this model. Loss of skeletal α-actin results in nemaline myopathy and we have previously shown that the pathological symptoms of the disease and reduction in muscle performance are recapitulated in a zebrafish antisense morpholino knockdown model. Here we reveal that a genetic actc1b mutant exhibits mild muscle defects and is unaffected by injection of the actc1b targeting morpholino. We further show that the milder phenotype results from a compensatory transcriptional upregulation of an actin paralogue providing a novel approach to be explored for the treatment of actin myopathy. Our findings provide further evidence that genetic compensation may influence the penetrance of disease-causing mutations.

Many healthy individuals carry loss of function mutations in essential genes that would normally be deleterious for survival. Intriguingly, it may be the presence of the genomic lesion itself in these individuals that triggers the compensatory pathways. It is not known how widespread this phenomenon is in vertebrate populations and how genetic compensation is activated. We have shown that knockdown of actin causes nemaline myopathy as indicated by the formation of nemaline bodies within the skeletal muscle and reduced muscle function which, remarkably, we did not observe in an actin genetic mutant. We have identified that protection from the disease phenotype results from transcriptional upregulation of an actin paralogue restoring actin protein in the skeletal muscle. This study demonstrates that genetic compensation may be more prevalent than previously anticipated and highlights phenotypic differences resulting from genetic mutations versus antisense knockdown approaches. Furthermore, we suggest that activating compensatory pathways may be exploited as a potential novel therapeutic approach for human disorders caused by loss of function mutations.

Canalization of Tomato Fruit Metabolism

To explore the genetic robustness (canalization) of metabolism, we examined the levels of fruit metabolites in multiple harvests of a tomato introgression line (IL) population. The IL partitions the whole genome of the wild species Solanum pennellii in the background of the cultivated tomato (Solanum lycopersicum). We identified several metabolite quantitative trait loci that reduce variability for both primary and secondary metabolites, which we named canalization metabolite quantitative trait loci (cmQTL). We validated nine cmQTL using an independent population of backcross inbred lines, derived from the same parents, which allows increased resolution in mapping the QTL previously identified in the ILs. These cmQTL showed little overlap with QTL for the metabolite levels themselves. Moreover, the intervals they mapped to harbored few metabolism-associated genes, suggesting that the canalization of metabolism is largely controlled by regulatory genes.

High-density ddRAD linkage and yield-related QTL mapping delimits a chromosomal region responsible for oil content in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is one of the most important oil crops almost all over the world. Seed-related traits, including oil content (OC), silique length (SL), seeds per silique (SS), and seed weight (SW), are primary targets for oil yield improvement. To dissect the genetic basis of these traits, 192 recombinant inbred lines (RILs) were derived from two parents with distinct oil content and silique length. High-density linkage map with a total length of 1610.4 cM were constructed using 1,329 double-digestion restriction site associated DNA (ddRAD) markers, 107 insertion/deletions (INDELs), and 90 well-distributed simple sequence repeats (SSRs) markers. A total of 37 consensus quantitative trait loci (QTLs) were detected for the four traits, with individual QTL explained 3.1-12.8% of the phenotypic variations. Interestingly, one OC consensus QTL (cqOCA10b) on chromosome A10 was consistently detected in all three environments, and explained 9.8% to 12.8% of the OC variation. The locus was further delimited into an approximately 614 kb genomic region, in which the flanking markers could be further evaluated for marker-assisted selection in rapeseed OC improvement and the candidate genes targeted for map-based cloning and genetic manipulation.

GC-MS Metabolomic Analysis to Reveal the Metabolites and Biological Pathways Involved in the Developmental Stages and Tissue Response of Panax ginseng

Ginsenosides, the major compounds present in ginseng, are known to have numerous physiological and pharmacological effects. The physiological processes, enzymes and genes involved in ginsenoside synthesis in P. ginseng have been well characterized. However, relatively little information is known about the dynamic metabolic changes that occur during ginsenoside accumulation in ginseng. To explore this topic, we isolated metabolites from different tissues at different growth stages, and identified and characterized them by using gas chromatography coupled with mass spectrometry (GC-MS). The results showed that a total of 30, 16, 20, 36 and 31 metabolites were identified and involved in different developmental stages in leaf, stem, petiole, lateral root and main root, respectively. To investigate the contribution of tissue to the biosynthesis of ginsenosides, we examined the metabolic changes of leaf, stem, petiole, lateral root and main root during five development stages: 1-, 2-, 3-, 4- and 5-years. The score plots of partial least squares-discriminate analysis (PLS-DA) showed clear discrimination between growth stages and tissue samples. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in the same tissue at different growth stages indicated profound biochemical changes in several pathways, including carbohydrate metabolism and pentose phosphate metabolism, in addition, the tissues displayed significant variations in amino acid metabolism, sugar metabolism and energy metabolism. These results should facilitate further dissection of the metabolic flux regulation of ginsenoside accumulation in different developmental stages or different tissues of ginseng.

Gene Mining for Proline Based Signaling Proteins in Cell Wall of Arabidopsis thaliana

The cell wall (CW) as a first line of defense against biotic and abiotic stresses is of primary importance in plant biology. The proteins associated with cell walls play a significant role in determining a plant's sustainability to adverse environmental conditions. In this work, the genes encoding cell wall proteins (CWPs) in Arabidopsis were identified and functionally classified using geneMANIA and GENEVESTIGATOR with published microarrays data. This yielded 1605 genes, out of which 58 genes encoded proline-rich proteins (PRPs) and glycine-rich proteins (GRPs). Here, we have focused on the cellular compartmentalization, biological processes, and molecular functioning of proline-rich CWPs along with their expression at different plant developmental stages. The mined genes were categorized into five classes on the basis of the type of PRPs encoded in the cell wall of Arabidopsis thaliana. We review the domain structure and function of each class of protein, many with respect to the developmental stages of the plant. We have then used networks, hierarchical clustering and correlations to analyze co-expression, co-localization, genetic, and physical interactions and shared protein domains of these PRPs. This has given us further insight into these functionally important CWPs and identified a number of potentially new cell-wall related proteins in A. thaliana.

MLPA-Based Analysis of Copy Number Variation in Plant Populations

Copy number variants (CNVs) are intraspecies duplications/deletions of large DNA segments (>1 kb). A growing number of reports highlight the functional and evolutionary impact of CNV in plants, increasing the need for appropriate tools that enable locus-specific CNV genotyping on a population scale. Multiplex ligation-dependent probe amplification (MLPA) is considered a gold standard in genotyping CNV in humans. Consequently, numerous commercial MLPA assays for CNV-related human diseases have been created. We routinely genotype complex multiallelic CNVs in human and plant genomes using the modified MLPA procedure based on fully synthesized oligonucleotide probes (90-200 nt), which greatly simplifies the design process and allows for the development of custom assays. Here, we present a step-by-step protocol for gene-specific MLPA probe design, multiplexed assay setup and data analysis in a copy number genotyping experiment in plants. As a case study, we present the results of a custom assay designed to genotype the copy number status of 12 protein coding genes in a population of 80 Arabidopsis accessions. The genes were pre-selected based on whole genome sequencing data and are localized in the genomic regions that display different levels of population-scale variation (non-variable, biallelic, or multiallelic, as well as CNVs overlapping whole genes or their fragments). The presented approach is suitable for population-scale validation of the CNV regions inferred from whole genome sequencing data analysis and for focused analysis of selected genes of interest. It can also be very easily adopted for any plant species, following optimization of the template amount and design of the appropriate control probes, according to the general guidelines presented in this paper.

Two Novel Flavin-Containing Monooxygenases Involved in Biosynthesis of Aliphatic Glucosinolates

Glucosinolates, a class of secondary metabolites from cruciferous plants, are derived from amino acids and have diverse biological activities, such as in biotic defense, depending on their side chain modification. The first structural modification step in the synthesis of aliphatic (methionine-derived) glucosinolates-S-oxygenation of methylthioalkyl glucosinolates to methylsulfinylalkyl glucosinolates-was found to be catalyzed by five flavin-containing monooxygenases (FMOs), FMOGS-OX1-5. Here, we report two additional FMOGS-OX enzymes, FMOGS-OX6, and FMOGS-OX7, encoded by At1g12130 and At1g12160, respectively. The overexpression of both FMOGS-OX6 and FMOGS-OX7 decreased the ratio of methylthioalkyl glucosinolates to the sum of methylthioalkyl and methylsulfinylalkyl glucosinolates, suggesting that the introduction of the two genes converted methylthioalkyl glucosinolates into methylsulfinylalkyl glucosinolates. Analysis of expression pattern revealed that the spatial expression of the two genes is quite similar and partially overlapped with the other FMOGS-OX genes, which are primarily expressed in vascular tissue. We further analyzed the responsive expression pattern of all the seven FMOGS-OX genes to exogenous treatment with abscisic acid, 1-aminocyclopropane-1-carboxylic acid (ACC), jasmonic acid (JA), salicylic acid, indole-3-acetic acid (IAA), and low and high temperatures. Although these genes showed same tendency toward the changing stimulus, the sensitivity of each gene was quite different. The variety in spatial expression among the FMOGS-OX genes while responding to environmental stimulus indicated a complex and finely tuned regulation of glucosinolates modifications. Identification of these two novel FMOGS-OX enzymes will enhance the understanding of glucosinolates modifications and the importance of evolution of these duplicated genes.

Molecular mechanisms governing differential robustness of development and environmental responses in plants

BACKGROUND

Robustness to genetic and environmental perturbation is a salient feature of multicellular organisms. Loss of developmental robustness can lead to severe phenotypic defects and fitness loss. However, perfect robustness, i.e. no variation at all, is evolutionarily unfit as organisms must be able to change phenotype to properly respond to changing environments and biotic challenges. Plasticity is the ability to adjust phenotypes predictably in response to specific environmental stimuli, which can be considered a transient shift allowing an organism to move from one robust phenotypic state to another. Plants, as sessile organisms that undergo continuous development, are particularly dependent on an exquisite fine-tuning of the processes that balance robustness and plasticity to maximize fitness.

SCOPE AND CONCLUSIONS

This paper reviews recently identified mechanisms, both systems-level and molecular, that modulate robustness, and discusses their implications for the optimization of plant fitness. Robustness in living systems arises from the structure of genetic networks, the specific molecular functions of the underlying genes, and their interactions. This very same network responsible for the robustness of specific developmental states also has to be built such that it enables plastic yet robust shifts in response to environmental changes. In plants, the interactions and functions of signal transduction pathways activated by phytohormones and the tendency for plants to tolerate whole-genome duplications, tandem gene duplication and hybridization are emerging as major regulators of robustness in development. Despite their obvious implications for plant evolution and plant breeding, the mechanistic underpinnings by which plants modulate precise levels of robustness, plasticity and evolvability in networks controlling different phenotypes are under-studied.

Evolutionary Rate Heterogeneity of Primary and Secondary Metabolic Pathway Genes in Arabidopsis thaliana

Primary metabolism is essential to plants for growth and development, and secondary metabolism helps plants to interact with the environment. Many plant metabolites are industrially important. These metabolites are produced by plants through complex metabolic pathways. Lack of knowledge about these pathways is hindering the successful breeding practices for these metabolites. For a better knowledge of the metabolism in plants as a whole, evolutionary rate variation of primary and secondary metabolic pathway genes is a prerequisite. In this study, evolutionary rate variation of primary and secondary metabolic pathway genes has been analyzed in the model plant Arabidopsis thaliana. Primary metabolic pathway genes were found to be more conserved than secondary metabolic pathway genes. Several factors such as gene structure, expression level, tissue specificity, multifunctionality, and domain number are the key factors behind this evolutionary rate variation. This study will help to better understand the evolutionary dynamics of plant metabolism.

The Arabidopsis gene DIG6 encodes a large 60S subunit nuclear export GTPase 1 that is involved in ribosome biogenesis and affects multiple auxin-regulated development processes

The circularly permuted GTPase large subunit GTPase 1 (LSG1) is involved in the maturation step of the 60S ribosome and is essential for cell viability in yeast. Here, an Arabidopsis mutant dig6 (drought inhibited growth of lateral roots) was isolated. The mutant exhibited multiple auxin-related phenotypes, which included reduced lateral root number, altered leaf veins, and shorter roots. Genetic mapping combined with next-generation DNA sequencing identified that the mutation occurred in AtLSG1-2. This gene was highly expressed in regions of auxin accumulation. Ribosome profiling revealed that a loss of function of AtLSG1-2 led to decreased levels of monosomes, further demonstrating its role in ribosome biogenesis. Quantitative proteomics showed that the expression of certain proteins involved in ribosome biogenesis was differentially regulated, indicating that ribosome biogenesis processes were impaired in the mutant. Further investigations showed that an AtLSG1-2 deficiency caused the alteration of auxin distribution, response, and transport in plants. It is concluded that AtLSG1-2 is integral to ribosome biogenesis, consequently affecting auxin homeostasis and plant development.

Scaling nitrogen and carbon interactions: what are the consequences of biological buffering?

Understanding the consequences of elevated CO₂ (eCO₂; 800 ppm) on terrestrial ecosystems is a central theme in global change biology, but relatively little is known about how altered plant C and N metabolism influences higher levels of biological organization. Here, we investigate the consequences of C and N interactions by genetically modifying the N-assimilation pathway in Arabidopsis and initiating growth chamber and mesocosm competition studies at current CO₂ (cCO₂; 400 ppm) and eCO₂ over multiple generations. Using a suite of ecological, physiological, and molecular genomic tools, we show that a single-gene mutant of a key enzyme (nia2) elicited a highly orchestrated buffering response starting with a fivefold increase in the expression of a gene paralog (nia1) and a 63% increase in the expression of gene network module enriched for N-assimilation genes. The genetic perturbation reduced amino acids, protein, and TCA-cycle intermediate concentrations in the nia2 mutant compared to the wild-type, while eCO₂ mainly increased carbohydrate concentrations. The mutant had reduced net photosynthetic rates due to a 27% decrease in carboxylation capacity and an 18% decrease in electron transport rates. The expression of these buffering mechanisms resulted in a penalty that negatively correlated with fitness and population dynamics yet showed only minor alterations in our estimates of population function, including total per unit area biomass, ground cover, and leaf area index. This study provides insight into the consequences of buffering mechanisms that occur post-genetic perturbations in the N pathway and the associated outcomes these buffering systems have on plant populations relative to eCO₂.

Patterns of metabolite changes identified from large-scale gene perturbations in Arabidopsis using a genome-scale metabolic network

Metabolomics enables quantitative evaluation of metabolic changes caused by genetic or environmental perturbations. However, little is known about how perturbing a single gene changes the metabolic system as a whole and which network and functional properties are involved in this response. To answer this question, we investigated the metabolite profiles from 136 mutants with single gene perturbations of functionally diverse Arabidopsis (Arabidopsis thaliana) genes. Fewer than 10 metabolites were changed significantly relative to the wild type in most of the mutants, indicating that the metabolic network was robust to perturbations of single metabolic genes. These changed metabolites were closer to each other in a genome-scale metabolic network than expected by chance, supporting the notion that the genetic perturbations changed the network more locally than globally. Surprisingly, the changed metabolites were close to the perturbed reactions in only 30% of the mutants of the well-characterized genes. To determine the factors that contributed to the distance between the observed metabolic changes and the perturbation site in the network, we examined nine network and functional properties of the perturbed genes. Only the isozyme number affected the distance between the perturbed reactions and changed metabolites. This study revealed patterns of metabolic changes from large-scale gene perturbations and relationships between characteristics of the perturbed genes and metabolic changes.

Transcriptional profiling in dioecious plant Populus cathayana reveals potential and sex-related molecular adaptations to solar UV-B radiation

There is increasing evidence confirming that dioecious plants have evolved a series of sex-related physiological mechanisms allowing acclimation to environmental stresses, but few studies have attempted to thoroughly uncover the potential molecular regulation underlying these responses. Recently, we found in the dioecious plant, Populus cathayana, that males possess a greater tolerance to enhanced UV-B radiation than do females. To continue this work further, comparative transcriptome analyses were carried out. A total of 362 and 296 differentially expressed genes (DEGs) were detected in males and females, respectively. Similarly, as shown in previous studies, these DEGs were involved in a set of conserved functions and pathways related to UV-B stress, regardless of the sex. On the other hand, sex-related responses via transcriptome remodeling were detected, as changes in sex-related gene expression occurred in some pathways. For example, many DEGs involved in amino acid metabolism were mainly upregulated in males but downregulated in females. In some pathways, DEGs were expressed predominantly or exclusively in one sex. This may directly contribute to sex-related physiological responses. Taken together, our results not only provide strong evidence that males and females have adopted sex-related molecular strategies in response to solar UV-B radiation but also suggest many new insights, in particular into molecular regulation, via sex-biased gene expression. We conclude that these results can further enhance the sustainable exploitation and management of P. cathayana plants.

Unraveling the mechanism underlying the glycosylation and methylation of anthocyanins in peach

Modification of anthocyanin plays an important role in increasing its stability in plants. Here, six anthocyanins were identified in peach (Prunus persica), and their structural diversity is attributed to glycosylation and methylation. Interestingly, peach is quite similar to the wild species Prunus ferganensis but differs from both Prunus davidiana and Prunus kansueasis in terms of anthocyanin composition in flowers. This indicates that peach is probably domesticated from P. ferganensis. Subsequently, genes responsible for both methylation and glycosylation of anthocyanins were identified, and their spatiotemporal expression results in different patterns of anthocyanin accumulation in flowers, leaves, and fruits. Two tandem-duplicated genes encoding flavonoid 3-O-glycosyltransferase (F3GT) in peach, PpUGT78A1 and PpUGT78A2, showed different activity toward anthocyanin, providing an example of divergent evolution of F3GT genes in plants. Two genes encoding anthocyanin O-methyltransferase (AOMT), PpAOMT1 and PpAOMT2, are expressed in leaves and flowers, but only PpAOMT2 is responsible for the O-methylation of anthocyanins at the 3' position in peach. In addition, our study reveals a novel branch of UGT78 genes in plants that lack the highly conserved intron 2 of the UGT gene family, with a great variation of the amino acid residue at position 22 of the plant secondary product glycosyltransferase box. Our results not only provide insights into the mechanisms underlying anthocyanin glycosylation and methylation in peach but will also aid in future attempts to manipulate flavonoid biosynthesis in peach as well as in other plants.

A plant-specific HUA2-LIKE (HULK) gene family in Arabidopsis thaliana is essential for development

In Arabidopsis thaliana, the HUA2 gene is required for proper expression of FLOWERING LOCUS C (FLC) and AGAMOUS, key regulators of flowering time and reproductive development, respectively. Although HUA2 is broadly expressed, plants lacking HUA2 function have only moderately reduced plant stature, leaf initiation rate and flowering time. To better understand HUA2 activity, and to test whether redundancy with similar genes underlies the absence of strong phenotypes in HUA2 mutant plants, we identified and subsequently characterized three additional HUA2-LIKE (HULK) genes in Arabidopsis. These genes form two clades (HUA2/HULK1 and HULK2/HULK3), with members broadly conserved in both vascular and non-vascular plants, but not present outside the plant kingdom. Plants with progressively reduced HULK activity had increasingly severe developmental defects, and plants homozygous for loss-of-function mutations in all four HULK genes were not recovered. Multiple mutants displayed reproductive, embryonic and post-embryonic abnormalities, and provide detailed insights into the overlapping and unique functions of individual HULK genes. With regard to flowering time, opposing influences were apparent: hua2 hulk1 plants were early-flowering, while hulk2 hulk3 mutants were late-flowering, and hua2 acted epistatically to cause early flowering in all combinations. Genome-wide expression profiling of mutant combinations using RNA-Seq revealed complex transcriptional changes in seedlings, with FLC, a known target of HUA2, among the most affected. Our studies, which include characterization of HULK expression patterns and subcellular localization, suggest that the HULK genes encode conserved nuclear factors with partially redundant but essential functions associated with diverse genetic pathways in plants.

Network types and their application in natural variation studies in plants

We are in the age of data-driven biology. Not even a decade after the invention of high-throughput sequencing technologies, there are methods that accurately monitor DNA polymorphisms, transcription profiles, methylation states, transcription factor binding sites, chromatin compactness, nucleosome positions, dynamic histone marks, and so on. We are starting to generate comparable amounts of protein or metabolite data. A key issue is how are we going to make sense of all this information. Network analysis is the most promising method to integrate, query and display large amounts of data for human interpretation. This review shortly summarizes the basic types of networks, their properties and limitations. In addition, I introduce the application of networks to the study of the molecular mechanisms behind natural phenotypic variation.

The transcription factor titration effect dictates level of gene expression

Models of transcription are often built around a picture of RNA polymerase and transcription factors (TFs) acting on a single copy of a promoter. However, most TFs are shared between multiple genes with varying binding affinities. Beyond that, genes often exist at high copy number-in multiple identical copies on the chromosome or on plasmids or viral vectors with copy numbers in the hundreds. Using a thermodynamic model, we characterize the interplay between TF copy number and the demand for that TF. We demonstrate the parameter-free predictive power of this model as a function of the copy number of the TF and the number and affinities of the available specific binding sites; such predictive control is important for the understanding of transcription and the desire to quantitatively design the output of genetic circuits. Finally, we use these experiments to dynamically measure plasmid copy number through the cell cycle.

Statistical mechanical model of coupled transcription from multiple promoters due to transcription factor titration

Transcription factors (TFs) with regulatory action at multiple promoter targets is the rule rather than the exception, with examples ranging from the cAMP receptor protein (CRP) in E. coli that regulates hundreds of different genes simultaneously to situations involving multiple copies of the same gene, such as plasmids, retrotransposons, or highly replicated viral DNA. When the number of TFs heavily exceeds the number of binding sites, TF binding to each promoter can be regarded as independent. However, when the number of TF molecules is comparable to the number of binding sites, TF titration will result in correlation ("promoter entanglement") between transcription of different genes. We develop a statistical mechanical model which takes the TF titration effect into account and use it to predict both the level of gene expression for a general set of promoters and the resulting correlation in transcription rates of different genes. Our results show that the TF titration effect could be important for understanding gene expression in many regulatory settings.

Molecular mechanisms of paralogous compensation and the robustness of cellular networks

Robustness is the ability of a system to maintain its function despite environmental or genetic perturbation. Genetic robustness is a key emerging property of living systems and is achieved notably by the presence of partially redundant parts that result from gene duplication. Functional overlap between paralogs allows them to compensate for each other's loss, as commonly revealed by aggravating genetic interactions. However, the molecular mechanisms linking the genotype (loss of function of a gene) to the phenotype (genetic buffering by a paralog) are still poorly understood and the molecular aspects of this compensation are rarely addressed in studies of gene duplicates. Here, we review molecular mechanisms of functional compensation between paralogous genes, many of which from studies that were not meant to study this phenomenon. We propose a standardized terminology and, depending on whether or not the molecular behavior of the intact gene is modified in response to the deletion of its paralog, we classify mechanisms of compensation into passive and active events. We further describe three non-exclusive mechanisms of active paralogous compensation for which there is evidence in the literature: changes in abundance, in localization, and in protein interactions. This review will serve as a framework for the genetic and molecular analysis of paralogous compensation, one of the universal features of genetic systems.

Integrated LC-MS/MS system for plant metabolomics

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is highly sensitive, selective, and enables extensive detection of metabolites within a sample. The result allows us to characterize comprehensive metabolite accumulation patterns without dependence on authentic standard compounds and isolation of the individual metabolites. A reference database search is essential for the structural assignment process of un-targeted MS and MS/MS data. Moreover, the characterization of unknown metabolites is challenging, since these cannot be assigned a candidate structure by using a reference database. In this case study, integrated LC-MS/MS based plant metabolomics allows us to detect several hundred metabolites in a sample; and integrated omics analyses, e.g., large-scale reverse genetics, linkage mapping, and association mapping, provides a powerful tool for candidate structure selection or rejection. We also examine emerging technology and applications for LC-MS/MS-based un-targeted plant metabolomics. These activities promote the characterization of massive extended detectable metabolites.

Systems analysis of plant functional, transcriptional, physical interaction, and metabolic networks

Physiological responses, developmental programs, and cellular functions rely on complex networks of interactions at different levels and scales. Systems biology brings together high-throughput biochemical, genetic, and molecular approaches to generate omics data that can be analyzed and used in mathematical and computational models toward uncovering these networks on a global scale. Various approaches, including transcriptomics, proteomics, interactomics, and metabolomics, have been employed to obtain these data on the cellular, tissue, organ, and whole-plant level. We summarize progress on gene regulatory, cofunction, protein interaction, and metabolic networks. We also illustrate the main approaches that have been used to obtain these networks, with specific examples from Arabidopsis thaliana, and describe the pros and cons of each approach.

Gene duplicability-connectivity-complexity across organisms and a neutral evolutionary explanation

Gene duplication has long been acknowledged by biologists as a major evolutionary force shaping genomic architectures and characteristics across the Tree of Life. Major research has been conducting on elucidating the fate of duplicated genes in a variety of organisms, as well as factors that affect a gene’s duplicability–that is, the tendency of certain genes to retain more duplicates than others. In particular, two studies have looked at the correlation between gene duplicability and its degree in a protein-protein interaction network in yeast, mouse, and human, and another has looked at the correlation between gene duplicability and its complexity (length, number of domains, etc.) in yeast. In this paper, we extend these studies to six species, and two trends emerge. There is an increase in the duplicability-connectivity correlation that agrees with the increase in the genome size as well as the phylogenetic relationship of the species. Further, the duplicability-complexity correlation seems to be constant across the species. We argue that the observed correlations can be explained by neutral evolutionary forces acting on the genomic regions containing the genes. For the duplicability-connectivity correlation, we show through simulations that an increasing trend can be obtained by adjusting parameters to approximate genomic characteristics of the respective species. Our results call for more research into factors, adaptive and non-adaptive alike, that determine a gene’s duplicability.