Insights into corn genes derived from large-scale cDNA sequencing

Biochemical and Epigenetic Regulation of Glutamate Metabolism in Maize (Zea mays L.) Leaves under Salt Stress

The effect of salt stress (150 mM NaCl) on the expression of genes, methylation of their promoters, and enzymatic activity of glutamate dehydrogenase (GDH), glutamate decarboxylase (GAD), and the 2-oxoglutarate (2-OG)-dehydrogenase (2-OGDH) complex was studied in maize (Zea mays L.). GDH activity increased continuously under salt stress, being 3-fold higher after 24 h. This was accompanied by the appearance of a second isoform with lower electrophoretic mobility. The expression of the Gdh1 gene strongly increased after 6-12 h of incubation, which corresponded to the demethylation of its promoter, while Gdh2 gene expression slightly increased after 2-6 h and then decreased. GAD activity gradually increased in the first 12 h, and then returned to the control level. This corresponded to the increase of Gad expression and its demethylation. Salt stress led to a 2-fold increase in the activity of 2-OGDH during the first 6 h of NaCl treatment, then the activity returned to the control level. Expression of the genes Ogdh1 and Ogdh3 peaked after 1-2 h of incubation. After 6-8 h with NaCl, the expression of these genes declined below the control levels, which correlated with the higher methylation of their promoters. We conclude that salt stress causes a redirection of the 2-OG flux to the γ-aminobutyric acid shunt via its amination to glutamate, by altering the expression of the Gdh1 and Gdh2 genes, which likely promotes the assembly of the native GDH molecule having a different subunit composition and greater affinity for 2-OG.

A teosinte-derived allele of ZmSC improves salt tolerance in maize

Maize, a salt-sensitive crop, frequently suffers severe yield losses due to soil salinization. Enhancing salt tolerance in maize is crucial for maintaining yield stability. To address this, we developed an introgression line (IL76) through introgressive hybridization between maize wild relatives Zea perennis, Tripsacum dactyloides, and inbred Zheng58, utilizing the tri-species hybrid MTP as a genetic bridge. Previously, genetic variation analysis identified a polymorphic marker on Zm00001eb244520 (designated as ZmSC), which encodes a vesicle-sorting protein described as a salt-tolerant protein in the NCBI database. To characterize the identified polymorphic marker, we employed gene cloning and homologous cloning techniques. Gene cloning analysis revealed a non-synonymous mutation at the 1847th base of ZmSCIL76 , where a guanine-to-cytosine substitution resulted in the mutation of serine to threonine at the 119th amino acid sequence (using ZmSCZ58 as the reference sequence). Moreover, homologous cloning demonstrated that the variation site derived from Z. perennis. Functional analyses showed that transgenic Arabidopsis lines overexpressing ZmSCZ58 exhibited significant reductions in leaf number, root length, and pod number, alongside suppression of the expression of genes in the SOS and CDPK pathways associated with Ca2+ signaling. Similarly, fission yeast strains expressing ZmSCZ58 displayed inhibited growth. In contrast, the ZmSCIL76 allele from Z. perennis alleviated these negative effects in both Arabidopsis and yeast, with the lines overexpressing ZmSCIL76 exhibiting significantly higher abscisic acid (ABA) content compared to those overexpressing ZmSCZ58 . Our findings suggest that ZmSC negatively regulates salt tolerance in maize by suppressing downstream gene expression associated with Ca2+ signaling in the CDPK and SOS pathways. The ZmSCIL76 allele from Z. perennis, however, can mitigate this negative regulatory effect. These results provide valuable insights and genetic resources for future maize salt tolerance breeding programs.

A reference-free algorithm discovers regulation in the plant transcriptome

Most plant genomes and their regulation remain unknown. We used SPLASH - a new, reference-genome free sequence variation detection algorithm - to analyze transcriptional and post-transcriptional regulation from RNA-seq data. We discovered differential homolog expression during maize pollen development, and imbibition-dependent cryptic splicing in Arabidopsis seeds. SPLASH enables discovery of novel regulatory mechanisms, including differential regulation of genes from hybrid parental haplotypes, without the use of alignment to a reference genome.

Transcriptome differential expression analysis of defoliation in different lemon varieties under drought treatment

'Allen Eureka' is a bud variety of Eureka lemon with excellent fruiting traits, but severe winter defoliation affects the following year's yield, and the response mechanism of lemon defoliation is currently unknown. Two lemon cultivars ('Allen Eureka' and 'Yunning No. 1') with different defoliation traits were used as materials to investigate the molecular regulatory mechanisms of different leaf abscission periods in lemons. The petiole abscission zone was collected at three different defoliation stages, namely, the predefoliation stage (k15), the middefoliation stage (k30), and the postdefoliation stage (k45). Transcriptome sequencing was performed to analyze the gene expression differences between these two cultivars. A total of 1141, 2695, and 1433 differentially expressed genes (DEGs) were obtained in k15, k30, and k45, respectively, and the number of DEGs in k30 was the largest. GO analysis revealed that the DEGs between the two cultivars were mainly enriched in processes related to hydrolase activity, chitinase activity, oxidoreductase activity, and transcription regulator activity in the defoliation stages. KEGG analysis showed that the DEGs were concentrated in k30, which involved plant hormone signal transduction, phenylpropanoid biosynthesis, and biosynthesis of amino acids. The expression trends of some DEGs suggested their roles in regulating defoliation in Lemon. Seven genes were obtained by WGCNA, including sorbitol dehydrogenase (CL9G068822012_alt, CL9G068820012_alt, CL9G068818012_alt), abscisic acid 8'-hydroxylase (CL8G064053012_alt, CL8G064054012_alt), and asparagine synthetase (CL8G065162012_alt, CL8G065151012_alt), suggesting that these genes may be involved in the regulation of lemon leaf abscission.

Studies on Candidate Genes Related to Flowering Time in a Multiparent Population of Maize Derived from Tropical and Temperate Germplasm

A comprehensive study on maize flowering traits, focusing on the regulation of flowering time and the elucidation of molecular mechanisms underlying the genes controlling flowering, holds the potential to significantly enhance our understanding of the associated regulatory gene network. In this study, three tropical maize inbreds, CML384, CML171, and CML444, were used, along with a temperate maize variety, Shen137, as parental lines to cross with Ye107. The resulting F1s underwent seven consecutive generations of self-pollination through the single-seed descent (SSD) method to develop a multiparent population. To investigate the regulation of maize flowering time-related traits and to identify loci and candidate genes, a genome-wide association study (GWAS) was conducted. GWAS analysis identified 556 SNPs and 12 candidate genes that were significantly associated with flowering time-related traits. Additionally, an analysis of the effect of the estimated breeding values of the subpopulations on flowering time was conducted to further validate the findings of the present study. Collectively, this study offers valuable insights into novel candidate genes, contributing to an improved understanding of maize flowering time-related traits. This information holds practical significance for future maize breeding programs aimed at developing high-yielding hybrids.

Plant root associated chitinases: structures and functions

Chitinases degrade chitin, a linear homopolymer of β-1,4-linked N-acetyl-D-glucosamine (GlcNAc) residues found in the cell walls of fungi and the exoskeletons of arthropods. They are secreted by the roots into the rhizosphere, a complex and dynamic environment where intense nutrient exchange occurs between plants and microbes. Here we modeled, expressed, purified, and characterized Zea mays and Oryza sativa root chitinases, and the chitinase of a symbiotic bacterium, Chitinophaga oryzae 1303 for their activities with chitin, di-, tri-, and tetra-saccharides and Aspergillus niger, with the goal of determining their role(s) in the rhizosphere and better understanding the molecular mechanisms underlying plant-microbe interactions. We show that Zea mays basic endochitinase (ZmChi19A) and Oryza sativa chitinase (OsChi19A) are from the GH19 chitinase family. The Chitinophaga oryzae 1303 chitinase (CspCh18A) belongs to the GH18 family. The three enzymes have similar apparent K M values of (20-40 µM) for the substrate 4-MU-GlcNAc3. They vary in their pH and temperature optima with OsChi19A activity optimal between pH 5-7 and 30-40°C while ZmChi19A and CspCh18A activities were optimal at pH 7-9 and 50-60°C. Modeling and site-directed mutation of ZmChi19A identified the catalytic cleft and the active residues E147 and E169 strategically positioned at ~8.6Å from each other in the folded protein. Cleavage of 4-MU-GlcNAc3 was unaffected by the absence of the CBD but diminished in the absence of the flexible C-terminal domain. However, unlike for the soluble substrate, the CBD and the newly identified flexible C-terminal domain were vital for inhibiting Aspergillus niger growth. The results are consistent with the involvement of the plant chitinases in defense against pathogens like fungi that have chitin exoskeletons. In summary, we have characterized the functional features and structural domains necessary for the activity of two plant root chitinases that are believed to be involved in plant defense and a bacterial chitinase that, along with the plant chitinases, may participate in nutrient recycling in the rhizosphere.

Unravelling the genetic framework associated with grain quality and yield-related traits in maize (Zea mays L.)

Maize serves as a crucial nutrient reservoir for a significant portion of the global population. However, to effectively address the growing world population's hidden hunger, it is essential to focus on two key aspects: biofortification of maize and improving its yield potential through advanced breeding techniques. Moreover, the coordination of multiple targets within a single breeding program poses a complex challenge. This study compiled mapping studies conducted over the past decade, identifying quantitative trait loci associated with grain quality and yield related traits in maize. Meta-QTL analysis of 2,974 QTLs for 169 component traits (associated with quality and yield related traits) revealed 68 MQTLs across different genetic backgrounds and environments. Most of these MQTLs were further validated using the data from genome-wide association studies (GWAS). Further, ten MQTLs, referred to as breeding-friendly MQTLs (BF-MQTLs), with a significant phenotypic variation explained over 10% and confidence interval less than 2 Mb, were shortlisted. BF-MQTLs were further used to identify potential candidate genes, including 59 genes encoding important proteins/products involved in essential metabolic pathways. Five BF-MQTLs associated with both quality and yield traits were also recommended to be utilized in future breeding programs. Synteny analysis with wheat and rice genomes revealed conserved regions across the genomes, indicating these hotspot regions as validated targets for developing biofortified, high-yielding maize varieties in future breeding programs. After validation, the identified candidate genes can also be utilized to effectively model the plant architecture and enhance desirable quality traits through various approaches such as marker-assisted breeding, genetic engineering, and genome editing.

A gene regulatory network inference model based on pseudo-siamese network

MOTIVATION

Gene regulatory networks (GRNs) arise from the intricate interactions between transcription factors (TFs) and their target genes during the growth and development of organisms. The inference of GRNs can unveil the underlying gene interactions in living systems and facilitate the investigation of the relationship between gene expression patterns and phenotypic traits. Although several machine-learning models have been proposed for inferring GRNs from single-cell RNA sequencing (scRNA-seq) data, some of these models, such as Boolean and tree-based networks, suffer from sensitivity to noise and may encounter difficulties in handling the high noise and dimensionality of actual scRNA-seq data, as well as the sparse nature of gene regulation relationships. Thus, inferring large-scale information from GRNs remains a formidable challenge.

RESULTS

This study proposes a multilevel, multi-structure framework called a pseudo-Siamese GRN (PSGRN) for inferring large-scale GRNs from time-series expression datasets. Based on the pseudo-Siamese network, we applied a gated recurrent unit to capture the time features of each TF and target matrix and learn the spatial features of the matrices after merging by applying the DenseNet framework. Finally, we applied a sigmoid function to evaluate interactions. We constructed two maize sub-datasets, including gene expression levels and GRNs, using existing open-source maize multi-omics data and compared them to other GRN inference methods, including GENIE3, GRNBoost2, nonlinear ordinary differential equations, CNNC, and DGRNS. Our results show that PSGRN outperforms state-of-the-art methods. This study proposed a new framework: a PSGRN that allows GRNs to be inferred from scRNA-seq data, elucidating the temporal and spatial features of TFs and their target genes. The results show the model's robustness and generalization, laying a theoretical foundation for maize genotype-phenotype associations with implications for breeding work.

Analysis of co-expression and gene regulatory networks associated with sterile lemma development in rice

BACKGROUND

The sterile lemma is a unique organ of the rice (Oryza sativa L.) spikelet. However, the characteristics and origin of the rice sterile lemma have not been determined unequivocally, so it is important to elucidate the molecular mechanism of the development of the sterile lemma.

RESULTS

In the paper, we outline the regulatory mechanism of sterile lemma development by LONG STERILE LEMMA1 (G1), which has been identified as the gene controlling sterile lemma development. Based on the comprehensive analyses of transcriptome dynamics during sterile lemma development with G1 alleles between wild-type (WT) and mutant (MT) in rice, we obtained co-expression data and regulatory networks related to sterile lemma development. Co-transfection assays of rice protoplasts confirmed that G1 affects the expression of various phytohormone-related genes by regulating a number of critical transcription factors, such as OsLBD37 and OSH1. The hormone levels in sterile lemmas from WT and MT of rice supports the hypotheses that lower auxin, lower gibberellin, and higher cytokinin concentrations are required to maintain a normal phenotype of sterile lemmas.

CONCLUSION

The regulatory networks have considerable reference value, and some of the regulatory relationships exhibiting strong correlations are worthy of further study. Taken together, these work provided a detailed guide for further studies into the molecular mechanism of sterile lemma development.

Identifying Signal-Crosstalk Mechanism in Maize Plants during Combined Salinity and Boron Stress Using Integrative Systems Biology Approaches

Combined stress has been seen as a major threat to world agriculture production. Maize is one of the leading cereal crops of the world due to its wide spectrum of growth conditions and is moderately sensitive to salt stress. A saline soil environment is a major factor that hinders its growth and overall yield and causes an increase in the concentration of micronutrients like boron, leading to excess over the requirement of the plant. Boron toxicity combined with salinity has been reported to be a serious threat to the yield and quality of maize. The response signatures of the maize plants to the combined effect of salinity and boron stress have not been studied well. We carried out an integrative systems-level analysis of the publicly available transcriptomic data generated on tolerant maize (Lluteño maize from the Atacama Desert, Chile) landrace under combined salt and boron stress. We identified significant biological processes that are differentially regulated in combined salt and boron stress in the leaves and roots of maize, respectively. Protein-protein interaction network analysis identified important roles of aldehyde dehydrogenase (ALDH), galactinol synthase 2 (GOLS2) proteins of leaf and proteolipid membrane potential regulator (pmpm4), metallothionein lea protein group 3 (mlg3), and cold regulated 410 (COR410) proteins of root in salt tolerance and regulating boron toxicity in maize. Identification of transcription factors coupled with regulatory network analysis using machine learning approach identified a few heat shock factors (HSFs) and NAC (NAM (no apical meristem, Petunia), ATAF1-2 (Arabidopsis thaliana activating factor), and CUC2 (cup-shaped cotyledon, Arabidopsis)) family transcription factors (TFs) to play crucial roles in salt tolerance, maintaining reactive oxygen species (ROS) levels and minimizing oxidative damage to the cells. These findings will provide new ways to design targeted functional validation experiments for developing multistress-resistant maize crops.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Critical assessment of computational tools for prokaryotic and eukaryotic promoter prediction

Promoters are crucial regulatory DNA regions for gene transcriptional activation. Rapid advances in next-generation sequencing technologies have accelerated the accumulation of genome sequences, providing increased training data to inform computational approaches for both prokaryotic and eukaryotic promoter prediction. However, it remains a significant challenge to accurately identify species-specific promoter sequences using computational approaches. To advance computational support for promoter prediction, in this study, we curated 58 comprehensive, up-to-date, benchmark datasets for 7 different species (i.e. Escherichia coli, Bacillus subtilis, Homo sapiens, Mus musculus, Arabidopsis thaliana, Zea mays and Drosophila melanogaster) to assist the research community to assess the relative functionality of alternative approaches and support future research on both prokaryotic and eukaryotic promoters. We revisited 106 predictors published since 2000 for promoter identification (40 for prokaryotic promoter, 61 for eukaryotic promoter, and 5 for both). We systematically evaluated their training datasets, computational methodologies, calculated features, performance and software usability. On the basis of these benchmark datasets, we benchmarked 19 predictors with functioning webservers/local tools and assessed their prediction performance. We found that deep learning and traditional machine learning-based approaches generally outperformed scoring function-based approaches. Taken together, the curated benchmark dataset repository and the benchmarking analysis in this study serve to inform the design and implementation of computational approaches for promoter prediction and facilitate more rigorous comparison of new techniques in the future.

Genome-Wide Identification, Characterization and Function Analysis of Lineage-Specific Genes in the Tea Plant Camellia sinensis

Genes that have no homologous sequences with other species are called lineage-specific genes (LSGs), are common in living organisms, and have an important role in the generation of new functions, adaptive evolution and phenotypic alteration of species. Camellia sinensis var. sinensis (CSS) is one of the most widely distributed cultivars for quality green tea production. The rich catechins in tea have antioxidant, free radical elimination, fat loss and cancer prevention potential. To further understand the evolution and utilize the function of LSGs in tea, we performed a comparative genomics approach to identify Camellia-specific genes (CSGs). Our result reveals that 1701 CSGs were identified specific to CSS, accounting for 3.37% of all protein-coding genes. The majority of CSGs (57.08%) were generated by gene duplication, and the time of duplication occurrence coincide with the time of two genome-wide replication (WGD) events that happened in CSS genome. Gene structure analysis revealed that CSGs have shorter gene lengths, fewer exons, higher GC content and higher isoelectric point. Gene expression analysis showed that CSG had more tissue-specific expression compared to evolutionary conserved genes (ECs). Weighted gene co-expression network analysis (WGCNA) showed that 18 CSGs are mainly associated with catechin synthesis-related pathways, including phenylalanine biosynthesis, biosynthesis of amino acids, pentose phosphate pathway, photosynthesis and carbon metabolism. Besides, we found that the expression of three CSGs (CSS0030246, CSS0002298, and CSS0030939) was significantly down-regulated in response to both types of stresses (salt and drought). Our study first systematically identified LSGs in CSS, and comprehensively analyzed the features and potential functions of CSGs. We also identified key candidate genes, which will provide valuable assistance for further studies on catechin synthesis and provide a molecular basis for the excavation of excellent germplasm resources.

Expression Elements Derived From Plant Sequences Provide Effective Gene Expression Regulation and New Opportunities for Plant Biotechnology Traits

Plant biotechnology traits provide a means to increase crop yields, manage weeds and pests, and sustainably contribute to addressing the needs of a growing population. One of the key challenges in developing new traits for plant biotechnology is the availability of expression elements for efficacious and predictable transgene regulation. Recent advances in genomics, transcriptomics, and computational tools have enabled the generation of new expression elements in a variety of model organisms. In this study, new expression element sequences were computationally generated for use in crops, starting from native Arabidopsis and maize sequences. These elements include promoters, 5' untranslated regions (5' UTRs), introns, and 3' UTRs. The expression elements were demonstrated to drive effective transgene expression in stably transformed soybean plants across multiple tissues types and developmental stages. The expressed transcripts were characterized to demonstrate the molecular function of these expression elements. The data show that the promoters precisely initiate transcripts, the introns are effectively spliced, and the 3' UTRs enable predictable processing of transcript 3' ends. Overall, our results indicate that these new expression elements can recapitulate key functional properties of natural sequences and provide opportunities for optimizing the expression of genes in future plant biotechnology traits.

Cloning and expression analysis of mevalonate kinase and phosphomevalonate kinase genes associated with the MVA pathway in Santalum album

Sandalwood (Santalum album L.) is highly valued for its fragrant heartwood and extracted oil. Santalols, which are the main components of that oil, are terpenoids, and these are biosynthesized via the mevalonic acid (MVA) pathway. Mevalonate kinase (MK) and phosphomevalonate kinase (PMK) are key enzymes in the MVA pathway. Little is known about the genes that encode MK and PMK in S. album or the mechanism that regulates their expression. To isolate and identify the functional genes involved in santalol biosynthesis in S. album, an MK gene designated as SaMK, and a PMK gene designated as SaPMK, were cloned from S. album. The sequences of these genes were analyzed. A bioinformatics analysis was conducted to assess the homology of SaMK and SaPMK with MK and PMK genes from other plants. The subcellular localization of SaMK and SaPMK proteins was also investigated, as was the functional complementation of SaMK and SaPMK in yeast. Our results show that the full-length cDNA sequences of SaMK and SaPMK were 1409 bp and 1679 bp long, respectively. SaMK contained a 1381 bp open reading frame (ORF) encoding a polypeptide of 460 amino acids and SaPMK contained a 1527 bp ORF encoding a polypeptide of 508 amino acids. SaMK and SaPMK showed high homology with MK and PMK genes of other plant species. Functional complementation of SaMK in a MK-deficient mutant yeast strain YMR208W and SaPMK in a PMK-deficient mutant yeast strain YMR220W confirmed that cloned SaMK and SaPMK cDNA encode a functional MK and PMK, respectively, mediating MVA biosynthesis in yeast. An analysis of tissue expression patterns revealed that SaMK and SaPMK were constitutively expressed in all the tested tissues. SaMK was highly expressed in young leaves but weakly expressed in sapwood. SaPMK was highly expressed in roots and mature leaves, but weakly expressed in young leaves. Induction experiments with several elicitors showed that SaMK and SaPMK expression was upregulated by methyl jasmonate. These results will help to further study the role of MK and PMK genes during santalol biosynthesis in S. album.

MicroRNA Zma-miR528 Versatile Regulation on Target mRNAs during Maize Somatic Embryogenesis

MicroRNAs (miRNAs) are small non-coding RNAs that regulate the accumulation and translation of their target mRNAs through sequence complementarity. miRNAs have emerged as crucial regulators during maize somatic embryogenesis (SE) and plant regeneration. A monocot-specific miRNA, mainly accumulated during maize SE, is zma-miR528. While several targets have been described for this miRNA, the regulation has not been experimentally confirmed for the SE process. Here, we explored the accumulation of zma-miR528 and several predicted targets during embryogenic callus induction, proliferation, and plantlet regeneration using the maize cultivar VS-535. We confirmed the cleavage site for all tested zma-miR528 targets; however, PLC1 showed very low levels of processing. The abundance of zma-miR528 slightly decreased in one month-induced callus compared to the immature embryo (IE) explant tissue. However, it displayed a significant increase in four-month sub-cultured callus, coincident with proliferation establishment. In callus-regenerated plantlets, zma-miR528 greatly decreased to levels below those observed in the initial explant. Three of the target transcripts (MATE, bHLH, and SOD1a) showed an inverse correlation with the miRNA abundance in total RNA samples at all stages. Using polysome fractionation, zma-miR528 was detected in the polysome fraction and exhibited an inverse distribution with the PLC1 target, which was not observed at total RNA. Accordingly, we conclude that zma-miR528 regulates multiple target mRNAs during the SE process by promoting their degradation, translation inhibition or both.

Perspective: 50 years of plant chromosome biology

The past 50 years has been the greatest era of plant science discovery, and most of the discoveries have emerged from or been facilitated by our knowledge of plant chromosomes. At last we have descriptive and mechanistic outlines of the information in chromosomes that programs plant life. We had almost no such information 50 years ago when few had isolated DNA from any plant species. The important features of genes have been revealed through whole genome comparative genomics and testing of variants using transgenesis. Progress has been enabled by the development of technologies that had to be invented and then become widely available. Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) have played extraordinary roles as model species. Unexpected evolutionary dramas were uncovered when learning that chromosomes have to manage constantly the vast numbers of potentially mutagenic families of transposons and other repeated sequences. The chromatin-based transcriptional and epigenetic mechanisms that co-evolved to manage the evolutionary drama as well as gene expression and 3-D nuclear architecture have been elucidated these past 20 years. This perspective traces some of the major developments with which I have become particularly familiar while seeking ways to improve crop plants. I draw some conclusions from this look-back over 50 years during which the scientific community has (i) exposed how chromosomes guard, readout, control, recombine, and transmit information that programs plant species, large and small, weed and crop, and (ii) modified the information in chromosomes for the purposes of genetic, physiological, and developmental analyses and plant improvement.

Prediction of Rice Transcription Start Sites Using TransPrise: A Novel Machine Learning Approach

As the interest in genetic resequencing increases, so does the need for effective mathematical, computational, and statistical approaches. One of the difficult problems in genome annotation is determination of precise positions of transcription start sites. In this paper, we present TransPrise-an efficient deep learning tool for predicting positions of eukaryotic transcription start sites. TransPrise offers significant improvement over existing promoter-prediction methods. To illustrate this, we compared predictions of TransPrise with the TSSPlant approach for well-annotated genome of Oryza sativa. Using a computer with a graphics processing unit, the run time of TransPrise is 250 min on a genome of 374 Mb long.We provide the full basis for the comparison and encourage users to freely access a set of our computational tools to facilitate and streamline their own analyses. The ready-to-use Docker image with all the necessary packages, models, and code as well as the source code of the TransPrise algorithm are available at http://compubioverne.group/ . The source code is ready to use and to be customized to predict TSS in any eukaryotic organism.

Variation in Maize Chlorophyll Biosynthesis Alters Plant Architecture

Chlorophyll is a tetrapyrrole metabolite essential for photosynthesis in plants. The first committed step of chlorophyll biosynthesis is catalyzed by a multimeric enzyme, magnesium chelatase, the subunit I of which is encoded by the oil yellow1 (oy1) gene in maize (Zea mays). A range of chlorophyll contents and net CO2 assimilation rates can be achieved in maize by combining a semidominant mutant allele of oy1 (Oy1-N1989) and a cis-regulatory modifier named very oil yellow1 (vey1) that varies between different inbred lines. We previously demonstrated that these allelic interactions can delay reproductive maturity. In this study, we demonstrate that multiple gross morphological traits respond to a reduction in chlorophyll. We found that stalk width, number of lateral branches (tillers), and branching of the inflorescence decline with a decrease in chlorophyll level. Chlorophyll deficit suppressed tillering in multiple maize mutants, including teosinte branched1, Tillering1, and grassy tillers1 In contrast to these traits, plant height showed a nonlinear response to chlorophyll levels. Weak suppression of Oy1-N1989 by vey1 B73 resulted in a significant increase in mutant plant height. By contrast, enhancement of the severity of the Oy1-N1989 phenotype by the vey1 Mo17 allele resulted in reduced plant height. We demonstrate that the effects of reduced chlorophyll contents on plant growth and development are complex and depend on the trait being measured. We propose that the lack of chlorophyll exerts growth control via energy balance sensing, which is upstream of the known genetic networks for branching and architecture.

Genome-wide identification, characterization and expression analysis of lineage-specific genes within Hanseniaspora yeasts

Lineage-specific genes (LSGs) are defined as genes with sequences that are not significantly similar to those in any other lineage. LSGs have been proposed, and sometimes shown, to have significant effects in the evolution of biological function. In this study, two sets of Hanseniaspora spp. LSGs were identified by comparing the sequences of the Kloeckera apiculata genome and of 80 other yeast genomes. This study identified 344 Hanseniaspora-specific genes (HSGs) and 109 genes ('orphan genes') specific to K. apiculata. Three thousand three hundred thirty-one K. apiculata genes that showed significant similarity to at least one sequence outside the Hanseniaspora were classified into evolutionarily conserved genes. We analyzed their sequence features, functional categories, gene origin, gene structure and gene expression. We also investigated the predicted cellular roles and Gene Ontology categories of the LSGs using functional inference. The patterns of the functions of LSGs do not deviate significantly from genome-wide average. The results showed that a few LSGs were formed by gene duplication, followed by rapid sequence divergence. Many of the HSGs and orphan genes exhibited altered expression in response to abiotic stress. Studying these LSGs might be helpful for understanding the molecular mechanism of yeast adaption.

Identification of Structural Variants in Two Novel Genomes of Maize Inbred Lines Possibly Related to Glyphosate Tolerance

To study genetic variations between genomes of plants that are naturally tolerant and sensitive to glyphosate, we used two Zea mays L. lines traditionally bred in Poland. To overcome the complexity of the maize genome, two sequencing technologies were employed: Illumina and Single Molecule Real-Time (SMRT) PacBio. Eleven thousand structural variants, 4 million SNPs and approximately 800 thousand indels differentiating the two genomes were identified. Detailed analyses allowed to identify 20 variations within the EPSPS gene, but all of them were predicted to have moderate or unknown effects on gene expression. Other genes of the shikimate pathway encoding bifunctional 3-dehydroquinate dehydratase/shikimate dehydrogenase and chorismate synthase were altered by variants predicted to have a high impact on gene expression. Additionally, high-impact variants located within the genes involved in the active transport of glyphosate through the cell membrane encoding phosphate transporters as well as multidrug and toxic compound extrusion have been identified.

Functional Aspects of Early Light-Induced Protein (ELIP) Genes from the Desiccation-Tolerant Moss Syntrichia caninervis

The early light-induced proteins (ELIPs) are postulated to act as transient pigment-binding proteins that protect the chloroplast from photodamage caused by excessive light energy. Desert mosses such as Syntrichia caninervis, that are desiccation-tolerant and homoiochlorophyllous, are often exposed to high-light conditions when both hydrated and dry ELIP transcripts are accumulated in response to dehydration. To gain further insights into ELIP gene function in the moss S. caninervis, two ELIP cDNAs cloned from S. caninervis, ScELIP1 and ScELIP2 and both sequences were used as the basis of a transcript abundance assessment in plants exposed to high-light, UV-A, UV-B, red-light, and blue-light. ScELIPs were expressed separately in an ArabidopsisELIP mutant Atelip. Transcript abundance for ScELIPs in gametophytes respond to each of the light treatments, in similar but not in identical ways. Ectopic expression of either ScELIPs protected PSII against photoinhibition and stabilized leaf chlorophyll content and thus partially complementing the loss of AtELIP2. Ectopic expression of ScELIPs also complements the germination phenotype of the mutant and improves protection of the photosynthetic apparatus of transgenic Arabidopsis from high-light stress. Our study extends knowledge of bryophyte photoprotection and provides further insight into the molecular mechanisms related to the function of ELIPs.

Assembly of chromosome-scale contigs by efficiently resolving repetitive sequences with long reads

The abundant repetitive sequences in complex eukaryotic genomes cause fragmented assemblies, which lose value as reference genomes, often due to incomplete gene sequences and unanchored or mispositioned contigs on chromosomes. Here we report a genome assembly method HERA, which resolves repeats efficiently by constructing a connection graph from an overlap graph. We test HERA on the genomes of rice, maize, human, and Tartary buckwheat with single-molecule sequencing and mapping data. HERA correctly assembles most of the previously unassembled regions, resulting in dramatically improved, highly contiguous genome assemblies with newly assembled gene sequences. For example, the maize contig N50 size reaches 61.2 Mb and the Tartary buckwheat genome comprises only 20 contigs. HERA can also be used to fill gaps and fix errors in reference genomes. The application of HERA will greatly improve the quality of new or existing assemblies of complex genomes.

TransPrise: a novel machine learning approach for eukaryotic promoter prediction

As interest in genetic resequencing increases, so does the need for effective mathematical, computational, and statistical approaches. One of the difficult problems in genome annotation is determination of precise positions of transcription start sites. In this paper we present TransPrise-an efficient deep learning tool for prediction of positions of eukaryotic transcription start sites. Our pipeline consists of two parts: the binary classifier operates the first, and if a sequence is classified as TSS-containing the regression step follows, where the precise location of TSS is being identified. TransPrise offers significant improvement over existing promoter-prediction methods. To illustrate this, we compared predictions of TransPrise classification and regression models with the TSSPlant approach for the well annotated genome of Oryza sativa. Using a computer equipped with a graphics processing unit, the run time of TransPrise is 250 minutes on a genome of 374 Mb long. The Matthews correlation coefficient value for TransPrise is 0.79, more than two times larger than the 0.31 for TSSPlant classification models. This represents a high level of prediction accuracy. Additionally, the mean absolute error for the regression model is 29.19 nt, allowing for accurate prediction of TSS location. TransPrise was also tested in Homo sapiens, where mean absolute error of the regression model was 47.986 nt. We provide the full basis for the comparison and encourage users to freely access a set of our computational tools to facilitate and streamline their own analyses. The ready-to-use Docker image with all necessary packages, models, code as well as the source code of the TransPrise algorithm are available at (http://compubioverne.group/). The source code is ready to use and customizable to predict TSS in any eukaryotic organism.

A 6&1-FEH Encodes an Enzyme for Fructan Degradation and Interact with Invertase Inhibitor Protein in Maize (Zea mays L.)

About 15% of higher plants have acquired the ability to convert sucrose into fructans. Fructan degradation is catalyzed by fructan exohydrolases (FEHs), which are structurally related to cell wall invertases (CWI). However, the biological function(s) of FEH enzymes in non-fructan species have remained largely enigmatic. In the present study, one maize CWI-related enzyme named Zm-6&1-FEH1, displaying FEH activity, was explored with respect to its substrate specificities, its expression during plant development, and its possible interaction with CWI inhibitor protein. Following heterologous expression in Pichia pastoris and in N. benthamiana leaves, recombinant Zm-6&1-FEH1 revealed substrate specificities of levan and inulin, and also displayed partially invertase activity. Expression of Zm-6&1-FEH1 as monitored by qPCR was strongly dependent on plant development and was further modulated by abiotic stress. To explore whether maize FEH can interact with invertase inhibitor protein, Zm-6&1-FEH1 and maize invertase inhibitor Zm-INVINH1 were co-expressed in N. benthamiana leaves. Bimolecular fluorescence complementation (BiFC) analysis and in vitro enzyme inhibition assays indicated productive complex formation. In summary, the results provide support to the hypothesis that in non-fructan species FEH enzymes may modulate the regulation of CWIs.

Physical Mapping of Peroxidase Genes and Development of Functional Markers for TaPod-D1 on Bread Wheat Chromosome 7D

Peroxidase (POD) activity in wheat (Triticum aestivum L.) grain influences natural carotenoid pigment content and is associated with the color of flour, and processing and product quality. Here, we report the molecular characterization and physical mapping of POD genes in bread wheat. The complete genomic DNA (gDNA) sequence of two POD genes (TaPod-A2 and TaPod-D1), and the partial gDNA sequence of two additional POD genes (TaPod-A3 and TaPod-B1) from wheat were characterized using in silico cloning and validated through laboratory experiments. Using a set of 21 nullisomic-tetrasomic (NT) lines, six group-7 ditelosomic (Dt) lines, and 38 group-7 deletion (Del) lines of Chinese Spring (CS), TaPod-A2 and TaPod-D1 were found to be physically located on 0.73-0.83 and on the most distal 0.39 fraction arm length (FL) of 7AS and 7DS in cv. CS, respectively; whereas, TaPod-A3 and TaPod-B1 were assigned to the 0.40-0.49 and 0.40-0.48 FL of 7AL and 7BL, respectively. Based on single nucleotide polymorphisms (SNPs) of two alleles at the TaPod-D1 locus, two functional markers POD-7D1 and POD-7D6 were developed, amplifying 540- and 640-bp, fragments in varieties with higher and lower POD activities, respectively. A total of 224 wheat varieties were analyzed and showed a significant association between the polymorphic fragments and POD activity using POD-7D1 and POD-7D6 markers. The analysis of variance (ANOVA) indicated the average POD activities of 115 varieties with TaPod-D1a were significantly lower than 109 varieties with TaPod-D1b (P < 0.01). This study provides useful information of the POD genes in bread wheat, insight into wheat genome synteny and structure, gene-specific markers, and contributes a valuable resource for quality improvement in wheat breeding programs.

Salt-Tolerant Antifungal and Antibacterial Activities of the Corn Defensin ZmD32

Pathogenic microbes are developing resistance to established antibiotics, making the development of novel antimicrobial molecules paramount. One major resource for discovery of antimicrobials is the arsenal of innate immunity molecules that are part of the first line of pathogen defense in many organisms. Gene encoded cationic antimicrobial peptides are a major constituent of innate immune arsenals. Many of these peptides exhibit potent antimicrobial activity in vitro. However, a major hurdle that has impeded their development for use in the clinic is the loss of activity at physiological salt concentrations, attributed to weakening of the electrostatic interactions between the cationic peptide and anionic surfaces of the microbial cells in the presence of salt. Using plant defensins we have investigated the relationship between the charge of an antimicrobial peptide and its activity in media with elevated salt concentrations. Plant defensins are a large class of antifungal peptides that have remarkable stability at extremes of pH and temperature as well as resistance to protease digestion. A search of a database of over 1200 plant defensins identified ZmD32, a defensin from Zea mays, with a predicted charge of +10.1 at pH 7, the highest of any defensin in the database. Recombinant ZmD32 retained activity against a range of fungal species in media containing elevated concentrations of salt. In addition, ZmD32 was active against Candida albicans biofilms as well as both Gram negative and Gram-positive bacteria. This broad spectrum antimicrobial activity, combined with a low toxicity on human cells make ZmD32 an attractive lead for development of future antimicrobial molecules.

A k-mer grammar analysis to uncover maize regulatory architecture

BACKGROUND

Only a small percentage of the genome sequence is involved in regulation of gene expression, but to biochemically identify this portion is expensive and laborious. In species like maize, with diverse intergenic regions and lots of repetitive elements, this is an especially challenging problem that limits the use of the data from one line to the other. While regulatory regions are rare, they do have characteristic chromatin contexts and sequence organization (the grammar) with which they can be identified.

RESULTS

We developed a computational framework to exploit this sequence arrangement. The models learn to classify regulatory regions based on sequence features - k-mers. To do this, we borrowed two approaches from the field of natural language processing: (1) "bag-of-words" which is commonly used for differentially weighting key words in tasks like sentiment analyses, and (2) a vector-space model using word2vec (vector-k-mers), that captures semantic and linguistic relationships between words. We built "bag-of-k-mers" and "vector-k-mers" models that distinguish between regulatory and non-regulatory regions with an average accuracy above 90%. Our "bag-of-k-mers" achieved higher overall accuracy, while the "vector-k-mers" models were more useful in highlighting key groups of sequences within the regulatory regions.

CONCLUSIONS

These models now provide powerful tools to annotate regulatory regions in other maize lines beyond the reference, at low cost and with high accuracy.

A rice Serine/Threonine receptor-like kinase regulates arbuscular mycorrhizal symbiosis at the peri-arbuscular membrane

In terrestrial ecosystems most plant species live in mutualistic symbioses with nutrient-delivering arbuscular mycorrhizal (AM) fungi. Establishment of AM symbioses includes transient, intracellular formation of fungal feeding structures, the arbuscules. A plant-derived peri-arbuscular membrane (PAM) surrounds the arbuscules, mediating reciprocal nutrient exchange. Signaling at the PAM must be well coordinated to achieve this dynamic cellular intimacy. Here, we identify the PAM-specific Arbuscular Receptor-like Kinase 1 (ARK1) from maize and rice to condition sustained AM symbiosis. Mutation of rice ARK1 causes a significant reduction in vesicles, the fungal storage structures, and a concomitant reduction in overall root colonization by the AM fungus Rhizophagus irregularis. Arbuscules, although less frequent in the ark1 mutant, are morphologically normal. Co-cultivation with wild-type plants restores vesicle and spore formation, suggesting ARK1 function is required for the completion of the fungal life-cycle, thereby defining a functional stage, post arbuscule development.

The peri-arbuscular membrane (PAM) mediates mutually-beneficial nutrient exchange between plants and arbuscular mycorrhizal (AM) fungi. Here the authors identify ARK1, a PAM-specific receptor-like kinase from rice that sustains AM symbiosis post-arbuscule development.

Endogenous tassel-specific small RNAs-mediated RNA interference enables a novel glyphosate-inducible male sterility system for commercial production of hybrid seed in Zea mays L

Hybrid crops produce higher yields than their inbred parents due to heterosis. For high purity of hybrid seeds, it is critical to eliminate self-pollination. Manual or mechanical removal of male parts (such as detasseling in maize) is labor-intensive, fuel and time-consuming, and can cause physical damage to female plants, resulting in significant seed yield reductions. Many male-sterility systems either require a maintainer for male-sterile line propagation or are often affected by environmental factors. Roundup® Hybridization System (RHS) utilizes glyphosate to induce male sterility, which effectively eliminates the need for maintainer lines and removal of male parts for commercial hybrid seed production. The first-generation RHS (RHS1) is based on low expression of a glyphosate-insensitive 5-enolpyruvylshikimate-3-phosphate synthase (CP4 EPSPS) in pollen. This report presents the second-generation RHS (RHS2) technology built on RNA interference (RNAi) combined with CP4 EPSPS. It utilizes maize endogenous male tissue-specific small interfering RNAs (mts-siRNAs) to trigger cleavage of the CP4 EPSPS mRNA specifically in tassels, resulting in glyphosate-sensitive male cells due to lack of the CP4 EPSPS protein. Male sterility is then induced by glyphosate application at the stages critical for pollen development, and the male-sterile plants are used as the female parent to produce hybrid seed. The endogenous mts-siRNAs are conserved across maize germplasms, and the inducible male sterility was replicated in representative germplasms through introgression of a CP4 EPSPS transgene containing the mts-siRNA target sequence. This technology combines the relative simplicity and convenience of a systemic herbicide spray methodology with targeted protein expression to create an inducible male sterility system for industrial production of row crop hybrid seeds in an environmentally-independent manner.

Where are the drought tolerant crops? An assessment of more than two decades of plant biotechnology effort in crop improvement

Since the dawn of modern biotechnology public and private enterprise have pursued the development of a new breed of drought tolerant crop products. After more than 20 years of research and investment only a few such products have reached the market. This is due to several technical and market constraints. The technical challenges include the difficulty in defining tractable single-gene trait development strategies, the logistics of moving traits from initial to commercial genetic backgrounds, and the disconnect between conditions in farmer's fields and controlled environments. Market constraints include the significant difficulty, and associated costs, in obtaining access to markets around the world. Advances in the biology of plant water management, including response to water deficit reveal new opportunities to improve crop response to water deficit and new genome-based tools promise to usher in the next era of crop improvement. As biotechnology looks to improve crop productivity under drought conditions, the environmental and food security advantages will influence public perception and shift the debate toward benefits rather than risks.

Combined analysis of genome-wide expression profiling of maize (Zea mays L.) leaves infected with Ustilago maydis

Although many gene expression profiling studies of maize leaves infected with Ustilago maydis have been published, heterogeneity of the results, caused by various data processing methods and pathogenic strains in different data sets, remains strong. Hence, we conducted a combined analysis of six genome-wide expression data sets of maize leaves infected with five different U. maydis strains by using the same pre-processing and quality control procedures. Six data sets were regrouped into five groups according to pathogenic strain used. Subsequently, each group of data set was processed by Multi-array Average for pre-processing and by pair-wise Pearson correlation for quality control. The differentially expressed genes were calculated by a standard linear mixed-effect model and then validated by various sensitivity analysis and multiple evidences. Finally, 44 unique differentially expressed genes were identified. Pathway enrichment analysis indicated that these genes related to response to fungus, oxidation-reduction, transferase activity, and several carbohydrate metabolic and catabolic processes. In addition, the hub genes within protein-protein interaction networks showed high relevance with the basic pathogenesis. We report a highly credible differentially expressed list, and the genes with multiple validations may denote a common signature of U. maydis in maize, which provides a new window for disease-resistant protection of maize plants.

Comparison of relative efficiency of genomic SSR and EST-SSR markers in estimating genetic diversity in sugarcane

Twenty-five primer pairs developed from genomic simple sequence repeats (SSR) were compared with 25 expressed sequence tags (EST) SSRs to evaluate the efficiency of these two sets of primers using 59 sugarcane genetic stocks. The mean polymorphism information content (PIC) of genomic SSR was higher (0.72) compared to the PIC value recorded by EST-SSR marker (0.62). The relatively low level of polymorphism in EST-SSR markers may be due to the location of these markers in more conserved and expressed sequences compared to genomic sequences which are spread throughout the genome. Dendrogram based on the genomic SSR and EST-SSR marker data showed differences in grouping of genotypes. A total of 59 sugarcane accessions were grouped into 6 and 4 clusters using genomic SSR and EST-SSR, respectively. The highly efficient genomic SSR could subcluster the genotypes of some of the clusters formed by EST-SSR markers. The difference in dendrogram observed was probably due to the variation in number of markers produced by genomic SSR and EST-SSR and different portion of genome amplified by both the markers. The combined dendrogram (genomic SSR and EST-SSR) more clearly showed the genetic relationship among the sugarcane genotypes by forming four clusters. The mean genetic similarity (GS) value obtained using EST-SSR among 59 sugarcane accessions was 0.70, whereas the mean GS obtained using genomic SSR was 0.63. Although relatively lower level of polymorphism was displayed by the EST-SSR markers, genetic diversity shown by the EST-SSR was found to be promising as they were functional marker. High level of PIC and low genetic similarity values of genomic SSR may be more useful in DNA fingerprinting, selection of true hybrids, identification of variety specific markers and genetic diversity analysis. Identification of diverse parents based on cluster analysis can be effectively done with EST-SSR as the genetic similarity estimates are based on functional attributes related to morphological/agronomical traits.

Identification of genes from the ICE-CBF-COR pathway under cold stress in Aegilops-Triticum composite group and the evolution analysis with those from Triticeae

Adverse environmental conditions limit various aspects of plant growth, productivity, and ecological distribution. To get more insights into the signaling pathways under low temperature, we identified 10 C-repeat binding factors (CBFs), 9 inducer of CBF expression (ICEs) and 10 cold-responsive (CORs) genes from Aegilops-Triticum composite group under cold stress. Conserved amino acids analysis revealed that all CBF, ICE, COR contained specific and typical functional domains. Phylogenetic analysis of CBF proteins from Triticeae showed that these CBF homologs were divided into 11 groups. CBFs from Triticum were found in every group, which shows that these CBFs generated prior to the divergence of the subfamilies of Triticeae. The evolutionary relationship among the ICE and COR proteins in Poaceae were divided into four groups with high multispecies specificity, respectively. Moreover, expression analysis revealed that mRNA accumulation was altered by cold treatment and the genes of three types involved in the ICE-CBF-COR signaling pathway were induced by cold stress. Together, the results make CBF, ICE, COR genes family in Triticeae more abundant, and provide a starting point for future studies on transcriptional regulatory network for improvement of chilling tolerance in crop.

Identification and Expression Analysis of Wheat TaGF14 Genes

The 14-3-3 gene family members play key roles in various cellular processes. However, little is known about the numbers and roles of 14-3-3 genes in wheat. The aims of this study were to identify TaGF14 numbers in wheat by searching its whole genome through blast, to study the phylogenetic relationships with other plant species and to discuss the functions of TaGF14s. The results showed that common wheat harbored 20 TaGF14 genes, located on wheat chromosome groups 2, 3, 4, and 7. Out of them, eighteen TaGF14s are non-ε proteins, and two wheat TaGF14 genes, TaGF14i and TaGF14f, are ε proteins. Phylogenetic analysis indicated that these genes were divided into six clusters: cluster 1 (TaGF14d, TaGF14g, TaGF14j, TaGF14h, TaGF14c, and TaGF14n); cluster 2 (TaGF14k); cluster 3 (TaGF14b, TaGF14l, TaGF14m, and TaGF14s); cluster 4 (TaGF14a, TaGF14e, and TaGF14r); cluster 5 (TaGF14i and TaGF14f); and cluster 6 (TaGF14o, TaGF14p, TaGF14q, and TaGF14t). Tissue-specific gene expressions suggested that all TaGF14s were likely constitutively expressed, except two genes, i.e., TaGF14p and TaGF14f. And the highest amount of TaGF14 transcripts were observed in developing grains at 20 days post anthesis (DPA), especially for TaGF14j and TaGF14l. After drought stress, five genes, i.e., TaGF14c, TaGF14d, TaGF14g, TaGF14h, and TaGF14j, were up-regulated expression under drought stress for both 1 and 6 h, suggesting these genes played vital role in combating against drought stress. However, all the TaGF14s were down-regulated expression under heat stress for both 1 and 6 h, indicating TaGF14s may be negatively associated with heat stress by reducing the expression to combat heat stress or through other pathways. These results suggested that cluster 1, e.g., TaGF14j, may participate in the whole wheat developing stages, e.g., grain-filling (starch biosynthesis) and may also participate in combating against drought stress. Subsequently, a homolog of TaGF14j, TaGF14-JM22, were cloned by RACE and used to validate its function. Immunoblotting results showed that TaGF14-JM22 protein, closely related to TaGF14d, TaGF14g, and TaGF14j, can interact with AGP-L, SSI, SSII, SBEIIa, and SBEIIb in developing grains, suggesting that TaGF14s located on group 4 may be involved in starch biosynthesis. Therefore, it is possible to develop starch-rich wheat cultivars by modifying TaGF14s.

Nucleotide patterns aiding in prediction of eukaryotic promoters

Computational analysis of promoters is hindered by the complexity of their architecture. In less studied genomes with complex organization, false positive promoter predictions are common. Accurate identification of transcription start sites and core promoter regions remains an unsolved problem. In this paper, we present a comprehensive analysis of genomic features associated with promoters and show that probabilistic integrative algorithms-driven models allow accurate classification of DNA sequence into “promoters” and “non-promoters” even in absence of the full-length cDNA sequences. These models may be built upon the maps of the distributions of sequence polymorphisms, RNA sequencing reads on genomic DNA, methylated nucleotides, transcription factor binding sites, as well as relative frequencies of nucleotides and their combinations. Positional clustering of binding sites shows that the cells of Oryza sativa utilize three distinct classes of transcription factors: those that bind preferentially to the [-500,0] region (188 “promoter-specific” transcription factors), those that bind preferentially to the [0,500] region (282 “5′ UTR-specific” TFs), and 207 of the “promiscuous” transcription factors with little or no location preference with respect to TSS. For the most informative motifs, their positional preferences are conserved between dicots and monocots.

Purification and biochemical characterization of a novel transglutaminase from Mythimna separata larvae (Noctuidae, Lepidoptera)

A novel transglutaminase (MsTGase) from Mythimna separata larvae was separated and purified; its biochemical property and enzymatic catalytic activities were investigated. MsTGase was obtained chromatographically by the precipitation of Sephadex G-100 gel and DEAE-Cellulose-52 ion-exchange column with 48-fold purification and a reproducible yield of approximately 12%. Molecular weight of the MsTGase was 63.5 KDa and its N-terminal amino acid sequence was GKIEEG-LVI. Michaelis constant of the MsTGase for the substrate N-CBZ-Gln-Gly was 12.83mM with a V_max of 7.99U/mL. Optimum conditions for MsTGase activity were at 42°C and pH7.5. The enzyme didn't possess metal ion at its catalytic active site; its activity could be significantly inhibited by Mg²⁺, but activated by Ca²⁺. Chlorpyrifos and spinosad showed a strong potential to increase MsTGase activity, supporting the view that MsTGase was a novel target. Moreover, the formation of intermolecular cross-links of casein and bovine serum albumin polymerized by MsTGase in the presence of DTT was observed. These findings pave the way for future studies on the physiological role of MsTGase and the potential impact of its regulation on MsTGase-associated pest management.

Evidence-based gene models for structural and functional annotations of the oil palm genome

Background

Oil palm is an important source of edible oil. The importance of the crop, as well as its long breeding cycle (10-12 years) has led to the sequencing of its genome in 2013 to pave the way for genomics-guided breeding. Nevertheless, the first set of gene predictions, although useful, had many fragmented genes. Classification and characterization of genes associated with traits of interest, such as those for fatty acid biosynthesis and disease resistance, were also limited. Lipid-, especially fatty acid (FA)-related genes are of particular interest for the oil palm as they specify oil yields and quality. This paper presents the characterization of the oil palm genome using different gene prediction methods and comparative genomics analysis, identification of FA biosynthesis and disease resistance genes, and the development of an annotation database and bioinformatics tools.

Results

Using two independent gene-prediction pipelines, Fgenesh++ and Seqping, 26,059 oil palm genes with transcriptome and RefSeq support were identified from the oil palm genome. These coding regions of the genome have a characteristic broad distribution of GC₃ (fraction of cytosine and guanine in the third position of a codon) with over half the GC₃-rich genes (GC₃ ≥ 0.75286) being intronless. In comparison, only one-seventh of the oil palm genes identified are intronless. Using comparative genomics analysis, characterization of conserved domains and active sites, and expression analysis, 42 key genes involved in FA biosynthesis in oil palm were identified. For three of them, namely EgFABF, EgFABH and EgFAD3, segmental duplication events were detected. Our analysis also identified 210 candidate resistance genes in six classes, grouped by their protein domain structures.

Conclusions

We present an accurate and comprehensive annotation of the oil palm genome, focusing on analysis of important categories of genes (GC₃-rich and intronless), as well as those associated with important functions, such as FA biosynthesis and disease resistance. The study demonstrated the advantages of having an integrated approach to gene prediction and developed a computational framework for combining multiple genome annotations. These results, available in the oil palm annotation database (http://palmxplore.mpob.gov.my), will provide important resources for studies on the genomes of oil palm and related crops.

Reviewers

This article was reviewed by Alexander Kel, Igor Rogozin, and Vladimir A. Kuznetsov.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-017-0191-4) contains supplementary material, which is available to authorized users.

A novel Sugarcane bacilliform virus promoter confers gene expression preferentially in the vascular bundle and storage parenchyma of the sugarcane culm

BACKGROUND

Saccharum species such as sugarcane and energy cane are key players in the expanding bioeconomy for sugars, bioenergy, and production of high-value proteins. Genomic tools such as culm-regulated promoters would be of great value in terms of improving biomass characteristics through enhanced carbon metabolism for sugar accumulation and/or fiber content for biofuel feedstock. Unlike the situation in dicots, monocot promoters currently used are limited and mostly derived from highly expressed constitutive plant genes and viruses. In this study, a novel promoter region of Sugarcane bacilliform virus (SCBV; genus Badnavirus, family Caulimoviridae), SCBV21 was cloned and mapped by deletion analysis and functionally characterized transiently in monocot and dicot species and stably in sugarcane.

RESULTS

In silico analysis of SCBV21 [1816 base pair (bp)] identified two putative promoter regions (PPR1 and PPR2) with transcription start sites (TSS1 and TSS2) and two TATA-boxes (TATAAAT and ATATAA), and several vascular-specific and regulatory elements. Deletion analysis revealed that the 710 bp region spanning PPR2 (with TSS2 and ATATAA) at the 3' end of SCBV21 retained the full promoter activity in both dicots and monocots, as shown by transient expression of the enhanced yellow fluorescent protein (EYFP) gene. In sugarcane young leaf segments, SCBV21 directed a 1.8- and 2.4-fold higher transient EYFP expression than the common maize ubiquitin 1 (Ubi1) and Cauliflower mosaic virus 35S promoters, respectively. In transgenic sugarcane, SCBV21 conferred a preferential expression of the β-glucuronidase (GUS) gene in leaves and culms and specifically in the culm storage parenchyma surrounding the vascular bundle and in vascular phloem cells. Among the transgenic events and tissues characterized in this study, the SCBV21 promoter frequently produced higher GUS activity than the Ubi1 or 35S promoters in a manner that was not obviously correlated with the transgene copy number.

CONCLUSIONS

The newly developed plant viral SCBV21 promoter is distinct from the few existing SCBV promoters in its sequence and expression pattern. The potential of SCBV21 as a tissue-regulated promoter with a strong activity in the culm vascular bundle and its storage parenchyma makes it useful in sugarcane engineering for improved carbon metabolism, increased bioenergy production, and enhanced stress tolerance.

Expression of four phosphate transporter genes from Finger millet (Eleusine coracana L.) in response to mycorrhizal colonization and Pi stress

Phosphorus (P) is a vital nutrient for plant growth and development, and is absorbed in cells with the help of membrane-spanning inorganic phosphate transporter (Pht) protein. Symbiosis with arbuscular mycorrhiza (AM) also helps in transporting P from the soil to plant and Pht proteins play an important role in it. To understand this phenomenon in Finger Mille plant, we have cloned four Pht genes from Finger millet, which shares the homology with Pht1 protein family of cereals. Expression pattern analysis during the AM infection indicated that EcPT4 gene was AM specific, and its expression was higher in roots where AM colonization percentage was high. The expression level of EcPT1-4 gene under the phosphorous (Pi) stress in seedlings was found to be consistent with its role in acquisition of phosphorus. Homology study of the EcPt proteins with Pht proteins of cereals shows close relationship. The findings of the study indicate that Pht1 family genes from finger millet can serve to be an important resource for the better understanding of phosphorus use efficiency.

Assessment of Natural Variability of Maize Lipid Transfer Protein Using a Validated Sandwich ELISA

Lipid transfer protein (LTP) is the main causative agent for rare food allergic reactions to maize. This paper describes a new, validated ELISA that accurately measures maize LTP concentrations from 0.2 to 6.4 ng/mL. The levels of LTP ranged from 171 to 865 μg/g of grain, a 5.1-fold difference, across a set of 49 samples of maize B73 hybrids derived from the Nested Association Mapping (NAM) founder lines and a diverse collection of landrace accessions from North and South America. A second set of 107 unique samples from 18 commercial hybrids grown over two years across 10 U.S. states showed a comparable range of LTP level (212-751 μg/g of grain). Statistical analysis showed that genetic and environmental factors contributed 63 and 6%, respectively, to the variance in LTP levels. Therefore, the natural variation of maize LTP is up to 5-fold different across a diverse collection of varieties that have a history of safe cultivation and consumption.

Seqping: gene prediction pipeline for plant genomes using self-training gene models and transcriptomic data

BACKGROUND

Gene prediction is one of the most important steps in the genome annotation process. A large number of software tools and pipelines developed by various computing techniques are available for gene prediction. However, these systems have yet to accurately predict all or even most of the protein-coding regions. Furthermore, none of the currently available gene-finders has a universal Hidden Markov Model (HMM) that can perform gene prediction for all organisms equally well in an automatic fashion.

RESULTS

We present an automated gene prediction pipeline, Seqping that uses self-training HMM models and transcriptomic data. The pipeline processes the genome and transcriptome sequences of the target species using GlimmerHMM, SNAP, and AUGUSTUS pipelines, followed by MAKER2 program to combine predictions from the three tools in association with the transcriptomic evidence. Seqping generates species-specific HMMs that are able to offer unbiased gene predictions. The pipeline was evaluated using the Oryza sativa and Arabidopsis thaliana genomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the pipeline was able to identify at least 95% of BUSCO's plantae dataset. Our evaluation shows that Seqping was able to generate better gene predictions compared to three HMM-based programs (MAKER2, GlimmerHMM and AUGUSTUS) using their respective available HMMs. Seqping had the highest accuracy in rice (0.5648 for CDS, 0.4468 for exon, and 0.6695 nucleotide structure) and A. thaliana (0.5808 for CDS, 0.5955 for exon, and 0.8839 nucleotide structure).

CONCLUSIONS

Seqping provides researchers a seamless pipeline to train species-specific HMMs and predict genes in newly sequenced or less-studied genomes. We conclude that the Seqping pipeline predictions are more accurate than gene predictions using the other three approaches with the default or available HMMs.

Main regulatory pathways, key genes and microRNAs involved in flower formation and development of moso bamboo (Phyllostachys edulis)

Moso bamboo is characterized by infrequent sexual reproduction and erratic flowering habit; however, the molecular biology of flower formation and development is not well studied in this species. We studied the molecular regulation mechanisms of moso bamboo development and flowering by selecting three key regulatory pathways: plant-pathogen interaction, plant hormone signal transduction and protein processing in endoplasmic reticulum at different stages of flowering in moso bamboo. We selected PheDof1, PheMADS14 and six microRNAs involved in the three pathways through KEGG pathway and cluster analysis. Subcellular localization, transcriptional activation, Western blotting, in situ hybridization and qRT-PCR were used to further investigate the expression patterns and regulatory roles of pivotal genes at different flower development stages. Differential expression patterns showed that PheDof1, PheMADS14 and six miRNAs may play vital regulatory roles in flower development and floral transition in moso bamboo. Our research paves way for further studies on metabolic regulatory networks and provides insight into the molecular regulation mechanisms of moso bamboo flowering and senescence.

Nucleotide diversity analysis highlights functionally important genomic regions

We analyzed functionality and relative distribution of genetic variants across the complete Oryza sativa genome, using the 40 million single nucleotide polymorphisms (SNPs) dataset from the 3,000 Rice Genomes Project (http://snp-seek.irri.org), the largest and highest density SNP collection for any higher plant. We have shown that the DNA-binding transcription factors (TFs) are the most conserved group of genes, whereas kinases and membrane-localized transporters are the most variable ones. TFs may be conserved because they belong to some of the most connected regulatory hubs that modulate transcription of vast downstream gene networks, whereas signaling kinases and transporters need to adapt rapidly to changing environmental conditions. In general, the observed profound patterns of nucleotide variability reveal functionally important genomic regions. As expected, nucleotide diversity is much higher in intergenic regions than within gene bodies (regions spanning gene models), and protein-coding sequences are more conserved than untranslated gene regions. We have observed a sharp decline in nucleotide diversity that begins at about 250 nucleotides upstream of the transcription start and reaches minimal diversity exactly at the transcription start. We found the transcription termination sites to have remarkably symmetrical patterns of SNP density, implying presence of functional sites near transcription termination. Also, nucleotide diversity was significantly lower near 3′ UTRs, the area rich with regulatory regions.

Genome-Wide Identification and Characterization of MicroRNAs and Target Genes in Lonicera japonica

MiRNAs function in post-transcriptional regulation of gene expression and play very important roles in plant development. Lonicera japonica is one of the important medicinal plants in China. However, few studies on the discovery of conserved and novel miRNAs from L. japonica were reported. In this study, we employed deep sequencing technology to identify miRNAs in leaf and flower tissues of L. japonica. A total of 22.97 million clean reads from flower and leaf tissues were obtained, which generated 146 conserved miRNAs distributed in 20 families and 110 novel miRNAs. Accordingly, 72 differentially expressed miRNAs (P≤0.001) between leaves and flowers and their potential target genes were identified and validated. The qRT-PCR validation showed that majority of the differentially expressed miRNAs showed significant tissue-specific expression in L. japonica. Furthermore, the miRNA-mRNA and mRNA-mRNA regulatory networks were constructed using Cytoscape software. Taken together, this study identified a large number of miRNAs and target genes in L. japonica, which not only provides the first global miRNA expression profiles, but also sheds light on functional genomics research on L. japonica in the future.