Genome-wide analysis of transcription factors involved in maize embryonic callus formation

GWAS determined genetic loci associated with callus induction in oil palm tissue culture

KEY MESSAGE

GWAS identified six loci at 25 kb downstream of WAK2, a crucial gene for cell wall and callus formation, enabling development of a SNP marker for enhanced callus induction potential. Efficient callus induction is vital for successful oil palm tissue culture, yet identifying genomic loci and markers for early detection of genotypes with high potential of callus induction remains unclear. In this study, immature male inflorescences from 198 oil palm accessions (dura, tenera and pisifera) were used as explants for tissue culture. Callus induction rates were collected at one-, two- and three-months after inoculation (C1, C2 and C3) as phenotypes. Resequencing generated 11,475,258 high quality single nucleotide polymorphisms (SNPs) as genotypes. GWAS was then performed, and correlation analysis revealed a positive association of C1 with both C2 (R = 0.81) and C3 (R = 0.50), indicating that C1 could be used as the major phenotype for callus induction rate. Therefore, only significant SNPs (P ≤ 0.05) in C1 were identified to develop markers for screening individuals with high potential of callus induction. Among 21 significant SNPs in C1, LD block analysis revealed six SNPs on chromosome 12 (Chr12) potentially linked to callus formation. Subsequently, 13 SNP markers were identified from these loci and electrophoresis results showed that marker C-12 at locus Chr12_12704856 can be used effectively to distinguish the GG allele, which showed the highest probability (69%) of callus induction. Furthermore, a rapid SNP variant detection method without electrophoresis was established via qPCR-based melting curve analysis. Our findings facilitated marker-assisted selection for specific palms with high potential of callus induction using immature male inflorescence as explant, aiding ortet palm selection in oil palm tissue culture.

Genome-Wide Association Studies of Embryogenic Callus Induction Rate in Peanut (Arachis hypogaea L.)

The capability of embryogenic callus induction is a prerequisite for in vitro plant regeneration. However, embryogenic callus induction is strongly genotype-dependent, thus hindering the development of in vitro plant genetic engineering technology. In this study, to examine the genetic variation in embryogenic callus induction rate (CIR) in peanut (Arachis hypogaea L.) at the seventh, eighth, and ninth subcultures (T7, T8, and T9, respectively), we performed genome-wide association studies (GWAS) for CIR in a population of 353 peanut accessions. The coefficient of variation of CIR among the genotypes was high in the T7, T8, and T9 subcultures (33.06%, 34.18%, and 35.54%, respectively), and the average CIR ranged from 1.58 to 1.66. A total of 53 significant single-nucleotide polymorphisms (SNPs) were detected (based on the threshold value -log10(p) = 4.5). Among these SNPs, SNPB03-83801701 showed high phenotypic variance and neared a gene that encodes a peroxisomal ABC transporter 1. SNPA05-94095749, representing a nonsynonymous mutation, was located in the Arahy.MIX90M locus (encoding an auxin response factor 19 protein) at T8, which was associated with callus formation. These results provide guidance for future elucidation of the regulatory mechanism of embryogenic callus induction in peanut.

Comprehensive Transcriptomic and Metabolic Profiling of Agrobacterium-tumefaciens-Infected Immature Wheat Embryos

The transformation efficiency (TE) was improved by a series of special chemical and physical methods using immature embryos from the cultivar Fielder, with the PureWheat technique. To analyze the reaction of immature embryos infected, which seemed to provide the necessary by Agrobacterium tumefaciens in PureWheat, a combination of scanning electron microscopy (SEM), complete transcriptome analysis, and metabolome analysis was conducted to understand the progress. The results of the SEM analysis revealed that Agrobacterium tumefaciens were deposited under the damaged cortex of immature embryos as a result of pretreatment and contacted the receptor cells to improve the TE. Transcriptome analysis indicated that the differentially expressed genes were mainly enriched in phenylpropanoid biosynthesis, starch and sucrose metabolism, plant-pathogen interaction, plant hormone signal transduction, and the MAPK (Mitogen-activated protein kinase) signaling pathway. By analyzing the correlation between differentially expressed genes and metabolites, the expression of many genes and the accumulation of metabolites were changed in glucose metabolism and the TCA cycle (Citrate cycle), as well as the amino acid metabolism; this suggests that the infection of wheat embryos with Agrobacterium is an energy-demanding process. The shikimate pathway may act as a hub between glucose metabolism and phenylpropanoid metabolism during Agrobacterium infection. The downregulation of the F5H gene and upregulation of the CCR gene led to the accumulation of lignin precursors through phenylpropanoid metabolism. In addition, several metabolic pathways and oxidases were found to be involved in the infection treatment, including melatonin biosynthesis, benzoxazinoid biosynthesis, betaine biosynthesis, superoxide dismutase, and peroxidase, suggesting that wheat embryos may be under the stress of Agrobacterium and, thus, undergo an oxidative stress response. These findings explore the physiological and molecular changes of immature embryos during the co-culture stage of the PureWheat technique and provide insights for Agrobacterium-mediated transgenic wheat experiments.

GWAS across multiple environments and WGCNA suggest the involvement of ZmARF23 in embryonic callus induction from immature maize embryos

Combined GWAS, WGCNA, and gene-based association studies identified the co-expression network and hub genes for maize EC induction. ZmARF23 bound to ZmSAUR15 promoter and regulated its expression, affecting EC induction. Embryonic callus (EC) induction in immature maize embryos shows high genotype dependence, which limits the application of genetic transformation in transgenic breeding and gene function elucidation in maize. Herein, we conducted a genome-wide association mapping (GWAS) for four EC induction-related traits, namely rate of embryonic callus induction (REC), increased callus diameter (ICD), ratio of shoot formation (RSF), and length of shoot (LS) across different environments. A total of 77 SNPs were significantly associated these traits under three environments and using the averages (across environments). Among these significant SNPs, five were simultaneously detected under multiple environments and 11 had respective phenotypic variation explained > 10%. A total of 257 genes were located in the linkage disequilibrium decay of these REC- and ICD-associated SNPs, of which 178 were responsive to EC induction. According to the expression values of the 178 genes, we performed a weighted gene co-expression network analysis (WGCNA) and revealed an EC induction-associated module and five hub genes. Hub gene-based association studies uncovered that the intragenic variations in GRMZM2G105473 and ZmARF23 influenced EC induction efficiency among different maize lines. Dual-luciferase reporter assay indicated that ZmARF23 bound to the promoter of a known causal gene (ZmSAUR15) for EC induction and positively regulated its expression on the transcription level. Our study will deepen the understanding of genetic and molecular mechanisms underlying EC induction and contribute to the use of genetic transformation in maize.

Combined QTL Mapping across Multiple Environments and Co-Expression Network Analysis Identified Key Genes for Embryogenic Callus Induction from Immature Maize Embryos

The ability of immature embryos to induce embryogenic callus (EC) is crucial for genetic transformation in maize, which is highly genotype-dependent. To dissect the genetic basis of maize EC induction, we conducted QTL mapping for four EC induction-related traits, the rate of embryogenic callus induction (REC), rate of shoot formation (RSF), length of shoot (LS), and diameter of callus (DC) under three environments by using an IBM Syn10 DH population derived from a cross of B73 and Mo17. These EC induction traits showed high broad-sense heritability (>80%), and significantly negative correlations were observed between REC and each of the other traits across multiple environments. A total of 41 QTLs for EC induction were identified, among which 13, 12, 10, and 6 QTLs were responsible for DC, RSF, LS, and REC, respectively. Among them, three major QTLs accounted for >10% of the phenotypic variation, including qLS1-1 (11.54%), qLS1-3 (10.68%), and qREC4-1 (11.45%). Based on the expression data of the 215 candidate genes located in these QTL intervals, we performed a weighted gene co-expression network analysis (WGCNA). A combined use of KEGG pathway enrichment and eigengene-based connectivity (KME) values identified the EC induction-associated module and four hub genes (Zm00001d028477, Zm00001d047896, Zm00001d034388, and Zm00001d022542). Gene-based association analyses validated that the variations in Zm00001d028477 and Zm00001d034388, which were involved in tryptophan biosynthesis and metabolism, respectively, significantly affected EC induction ability among different inbred lines. Our study brings novel insights into the genetic and molecular mechanisms of EC induction and helps to promote marker-assisted selection of high-REC varieties in maize.

New Insights Into Tissue Culture Plant-Regeneration Mechanisms

Plant regeneration occurs when plants repair or replace damaged structures based on the totipotency and pluripotency of their cells. Tissue culture is one of the most widely used regenerative technologies. Recently, a series of breakthroughs were made in the study of plant regeneration. This review summarizes two regenerative pathways in tissue culture: somatic embryogenesis and de novo organogenesis. Furthermore, we review the environmental factors influencing plant regeneration from explant sources, basal culture medium, plant growth regulators, and light/dark treatment. Additionally, we analyse the molecular mechanisms underlying two pathways. This knowledge will promote an understanding of the fundamental principles of plant regeneration from precursor cells and lay a solid foundation for applying plant micropropagation and genetic modification.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Effect of aneuploidy of a non-essential chromosome on gene expression in maize

The non-essential supernumerary maize (Zea mays) B chromosome (B) has recently been shown to contain active genes and to be capable of impacting gene expression of the A chromosomes. However, the effect of the B chromosome on gene expression is still unclear. In addition, it is unknown whether the accumulation of the B chromosome has a cumulative effect on gene expression. To examine these questions, the global expression of genes, microRNAs (miRNAs), and transposable elements (TEs) of leaf tissue of maize W22 plants with 0-7 copies of the B chromosome was studied. All experimental genotypes with B chromosomes displayed a trend of upregulated gene expression for a subset of A-located genes compared to the control. Over 3000 A-located genes are significantly differentially expressed in all experimental genotypes with the B chromosome relative to the control. Modulations of these genes are largely determined by the presence rather than the copy number of the B chromosome. By contrast, the expression of most B-located genes is positively correlated with B copy number, showing a proportional gene dosage effect. The B chromosome also causes increased expression of A-located miRNAs. Differentially expressed miRNAs potentially regulate their targets in a cascade of effects. Furthermore, the varied copy number of the B chromosome leads to the differential expression of A-located and B-located TEs. The findings provide novel insights into the function and properties of the B chromosome.

Genetic variations in ZmSAUR15 contribute to the formation of immature embryo-derived embryonic calluses in maize

The ability of immature maize (Zea mays) embryos to form embryonic calluses (ECs) is highly genotype dependent, which limits transgenic breeding development in maize. Here, we report the association map-based cloning of ZmSAUR15 using an association panel (AP) consisting of 309 inbred lines with diverse formation abilities for ECs. We demonstrated that ZmSAUR15, which encodes a small auxin-upregulated RNA, acts as a negative effector in maize EC induction. Polymorphisms in the ZmSAUR15 promoter that influence the expression of ZmSAUR15 transcripts modulate the EC induction capacity in maize. ZmSAUR15 is involved in indole-3-acetic acid biosynthesis and cell division in immature embryo-derived callus. The ability of immature embryos to induce EC formation can be improved by the knockout of ZmSAUR15, which consequently increases the callus regeneration efficiency. Our study provides new insights into overcoming the genotypic limitations associated with EC formation and improving genetic transformation in maize.

Transcriptomic analysis of the maize inbred line Chang7-2 and a large-grain mutant tc19

Backgrounds

Grain size is a key factor in crop yield that gradually develops after pollination. However, few studies have reported gene expression patterns in maize grain development using large-grain mutants. To investigate the developmental mechanisms of grain size, we analyzed a large-grain mutant, named tc19, at the morphological and transcriptome level at five stages corresponding to days after pollination (DAP).

Results

After maturation, the grain length, width, and thickness in tc19 were greater than that in Chang7-2 (control) and increased by 3.57, 8.80, and 3.88%, respectively. Further analysis showed that grain width and 100-kernel weight in tc19 was lower than in Chang7-2 at 14 and 21 DAP, but greater than that in Chang7-2 at 28 DAP, indicating that 21 to 28 DAP was the critical stage for kernel width and weight development. For all five stages, the concentrations of auxin and brassinosteroids were significantly higher in tc19 than in Chang7-2. Gibberellin was higher at 7, 14, and 21 DAP, and cytokinin was higher at 21 and 35 DAP, in tc19 than in Chang7-2. Through transcriptome analysis at 14, 21, and 28 DAP, we identified 2987, 2647 and 3209 differentially expressed genes (DEGs) between tc19 and Chang7-2. By using KEGG analysis, 556, 500 and 633 DEGs at 14, 21 and 28 DAP were pathway annotated, respectively, 77 of them are related to plant hormone signal transduction pathway. ARF3, AO2, DWF4 and XTH are higher expressed in tc19 than that in Chang7-2.

Conclusions

We found some DEGs in maize grain development by using Chang7-2 and a large-grain mutant tc19. These DEGs have potential application value in improving maize performance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08230-9.

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.

miRNA expression profiling and zeatin dynamic changes in a new model system of in vivo indirect regeneration of tomato

Callus formation and adventitious shoot differentiation could be observed on the cut surface of completely decapitated tomato plants. We propose that this process can be used as a model system to investigate the mechanisms that regulate indirect regeneration of higher plants without the addition of exogenous hormones. This study analyzed the patterns of trans-zeatin and miRNA expression during in vivo regeneration of tomato. Analysis of trans-zeatin revealed that the hormone cytokinin played an important role in in vivo regeneration of tomato. Among 183 miRNAs and 1168 predicted target genes sequences identified, 93 miRNAs and 505 potential targets were selected based on differential expression levels for further characterization. Expression patterns of six miRNAs, including sly-miR166, sly-miR167, sly-miR396, sly-miR397, novel 156, and novel 128, were further validated by qRT-PCR. We speculate that sly-miR156, sly-miR160, sly-miR166, and sly-miR397 play major roles in callus formation of tomato during in vivo regeneration by regulating cytokinin, IAA, and laccase levels. Overall, our microRNA sequence and target analyses of callus formation during in vivo regeneration of tomato provide novel insights into the regulation of regeneration in higher plants.

Integrated transcriptome, small RNA, and degradome analysis reveals the complex network regulating starch biosynthesis in maize

Background

Starch biosynthesis in endosperm is a key process influencing grain yield and quality in maize. Although a number of starch biosynthetic genes have been well characterized, the mechanisms by which the expression of these genes is regulated, especially in regard to microRNAs (miRNAs), remain largely unclear.

Results

Sequence data for small RNAs, degradome, and transcriptome of maize endosperm at 15 and 25 d after pollination (DAP) from inbred lines Mo17 and Ji419, which exhibit distinct starch content and starch granule structure, revealed the mediation of starch biosynthetic pathways by miRNAs. Transcriptome analysis of these two lines indicated that 33 of 40 starch biosynthetic genes were differentially expressed, of which 12 were up-regulated in Ji419 at 15 DAP, one was up-regulated in Ji419 at 25 DAP, 14 were up-regulated in Ji419 at both 15 and 25 DAP, one was down-regulated in Ji419 at 15 DAP, two were down-regulated in Ji419 at 25 DAP, and three were up-regulated in Ji419 at 15 DAP and down-regulated in Ji419 at 25 DAP, compared with Mo17. Through combined analyses of small RNA and degradome sequences, 22 differentially expressed miRNAs were identified, including 14 known and eight previously unknown miRNAs that could target 35 genes. Furthermore, a complex co-expression regulatory network was constructed, in which 19 miRNAs could modulate starch biosynthesis in endosperm by tuning the expression of 19 target genes. Moreover, the potential operation of four miRNA-mediated pathways involving transcription factors, miR169a-NF-YA1-GBSSI/SSIIIa and miR169o-GATA9-SSIIIa/SBEIIb, was validated via analyses of expression pattern, transient transformation assays, and transactivation assays.

Conclusion

Our results suggest that miRNAs play a critical role in starch biosynthesis in endosperm, and that miRNA-mediated networks could modulate starch biosynthesis in this tissue. These results have provided important insights into the molecular mechanism of starch biosynthesis in developing maize endosperm.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5945-1) contains supplementary material, which is available to authorized users.

Transcriptome sequencing analysis of maize embryonic callus during early redifferentiation

BACKGROUND

Maize is one of the primary crops of genetic manipulation, which provides an excellent means of promoting stress resistance and increasing yield. However, the differences in induction and regeneration capacity of embryonic callus (EC) among various genotypes result in genotypic dependence in genetic transformation.

RESULTS

In this study, embryonic calli of two maize inbred lines with strong redifferentiation capacity and two lines with weak redifferentiation capability were separately subjected to transcriptome sequencing analysis during the early redifferentiation stages (stage I, 1-3 d; stage II, 4-6 d; stage III, 7-9 d) along with their corresponding controls. A total of ~ 654.72 million cDNA clean reads were yielded, and 62.64%~ 69.21% clean reads were mapped to the reference genome for each library. In comparison with the control, the numbers of differentially expressed genes (DEGs) for the four inbred lines identified in the three stages ranged from 1694 to 7193. By analyzing the common and specific DEGs of the four materials, we found that there were 321 upregulated genes and 386 downregulated genes identified in the high-regeneration lines (141 and DH40), whereas 611 upregulated genes and 500 downregulated genes were specifically expressed in the low-regeneration lines (ZYDH381-1 and DH3732). Analysis of the DEG expression patterns indicated a sharp change at stage I in both the high- and low-regeneration lines, which suggested that stage I constitutes a crucial period for EC regeneration. Notably, the specific common DEGs of 141 and DH40 were mainly associated with photosynthesis, porphyrin and chlorophyll metabolism, ribosomes, and plant hormone signal transduction. In contrast, the DEGs in ZYDH381-1 and DH3732 were mainly related to taurine and hypotaurine metabolism, nitrogen metabolism, fatty acid elongation, starch and sucrose metabolism, phenylpropanoid biosynthesis, and plant circadian rhythm. More importantly, WOX genes, which have an ancestral role in embryo development in seed plants and promote the regeneration of transformed calli, were specifically upregulated in the two high-regeneration lines.

CONCLUSIONS

Our research contributes to the elucidation of molecular regulation during early redifferentiation in the maize embryonic callus.

Genome-wide characterization and analysis of bHLH transcription factors in Panax ginseng

Ginseng (Panax ginseng C.A. Meyer) is one of the best-selling herbal medicines, with ginsenosides as its main pharmacologically active constituents. Although extensive chemical and pharmaceutical studies of these compounds have been performed, genome-wide studies of the basic helix-loop-helix (bHLH) transcription factors of ginseng are still limited. The bHLH transcription factor family is one of the largest transcription factor families found in eukaryotic organisms, and these proteins are involved in a myriad of regulatory processes. In our study, 169 bHLH transcription factor genes were identified in the genome of P. ginseng, and phylogenetic analysis indicated that these PGbHLHs could be classified into 24 subfamilies. A total of 21 RNA-seq data sets, including two sequencing libraries for jasmonate (JA)-responsive and 19 reported libraries for organ-specific expression analyses were constructed. Through a combination of gene-specific expression patterns and chemical contents, 6 PGbHLH genes from 4 subfamilies were revealed to be potentially involved in the regulation of ginsenoside biosynthesis. These 6 PGbHLHs, which had distinct target genes, were further divided into two groups depending on the absence of MYC-N structure. Our results would provide a foundation for understanding the molecular basis and regulatory mechanisms of bHLH transcription factor action in P. ginseng.

Effect of plant growth regulators on indirect shoot organogenesis of Ficus religiosa through seedling derived petiole segments

Ficus religiosa is known as a long-lived multipurpose forest tree. The tree plays an important role for religious, medicinal, and ornamental purposes. However, the propagation rate of Ficus religiosa is low in natural habitat so the plant tissue culture techniques are an applicable method for multiplication of this valuable medicinal plants. Thus, the aim of this study is to understand the effect of different auxin/cytokinin ratios on indirect shoot organogenesis of this plant. According to our results, the maximum callus induction frequency (100%) was obtained on Murashige and Skoog (MS) medium supplemented with 0.5 mg/l 2,4-dichlorophenoxyacetic acid (2,4-D) plus 0.05 mg/l 6-benzylaminopurine (BAP) from petiole segments. For shoot induction purpose, the yellow-brownish, friable, organogenic calli were inoculated on shoot induction medium. On MS medium supplemented with 1.5 mg/l BAP and 0.15 mg/l Indole-3-butyric acid (IBA), 96.66% of the petiole-derived calli responded with an average number of 3.56 shoots per culture. The highest root formation frequency (96.66%), root number (5.5), and root length (4.83 cm) were achieved on MS medium containing 2.0 mg/l IBA plus 0.1 mg/l Naphthaleneacetic acid (NAA). The rooted shoots were successfully transferred to field condition and the substrate with the mixture of cocopeat and perlite (1:1) had the highest survival rate (96.66%). This is the first report of an effective in vitro organogenesis protocol for F. religiosa by indirect shoot organogenesis through axenic seedling derived petiole explants, which can be efficiently employed for conservation of this important medicinal plant species as well as the utilization of active biomolecules.

Genetic Dissection of Maize Embryonic Callus Regenerative Capacity Using Multi-Locus Genome-Wide Association Studies

The regenerative capacity of the embryonic callus, a complex quantitative trait, is one of the main limiting factors for maize transformation. This trait was decomposed into five traits, namely, green callus rate (GCR), callus differentiating rate (CDR), callus plantlet number (CPN), callus rooting rate (CRR), and callus browning rate (CBR). To dissect the genetic foundation of maize transformation, in this study multi-locus genome-wide association studies (GWAS) for the five traits were performed in a population of 144 inbred lines genotyped with 43,427 SNPs. Using the phenotypic values in three environments and best linear unbiased prediction (BLUP) values, as a result, a total of 127, 56, 160, and 130 significant quantitative trait nucleotides (QTNs) were identified by mrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB, respectively. Of these QTNs, 63 QTNs were commonly detected, including 15 across multiple environments and 58 across multiple methods. Allele distribution analysis showed that the proportion of superior alleles for 36 QTNs was <50% in 31 elite inbred lines. Meanwhile, these superior alleles had obviously additive effect on the regenerative capacity. This indicates that the regenerative capacity-related traits can be improved by proper integration of the superior alleles using marker-assisted selection. Moreover, a total of 40 candidate genes were found based on these common QTNs. Some annotated genes were previously reported to relate with auxin transport, cell fate, seed germination, or embryo development, especially, GRMZM2G108933 (WOX2) was found to promote maize transgenic embryonic callus regeneration. These identified candidate genes will contribute to a further understanding of the genetic foundation of maize embryonic callus regeneration.

Small RNA differential expression and regulation in Tuxpeño maize embryogenic callus induction and establishment

Somatic embryogenesis represents an alternative developmental process used to achieve genetic transformation and to approach key questions in maize development. It is known that embryogenic callus induction and plant regeneration are accompanied by microRNA expression changes. However, small RNA (sRNA) populations have not been explored during the proliferative callus subculture establishment and their impact on maintaining the dedifferentiated status and embryogenic potential is far from being completely understood. Here we globally tested the sRNA populations in explants (immature embryos), induced and established maize embryogenic callus from the Mexican cultivar VS-535, Tuxpeño landrace. We detected readjustments in 24 nt and 21-22 nt sRNAs during the embryogenic callus (EC) establishment and maintenance. A follow up on specific microRNAs (miRNAs) indicated that miRNAs related to stress response substantially increase upon the callus proliferation establishment, correlating with a reduction in some of their target levels. On the other hand, while 24 nt-long heterochromatic small interfering RNAs (hc-siRNAs) derived from transposable retroelements transiently decreased in abundance during the EC establishment, a population of 22 nt-hc-siRNAs increased. This was accompanied by reduction in transposon expression in the established callus subcultures. We conclude that stress- and development-related miRNAs are highly expressed upon maize EC callus induction and during maintenance of the subcultures, while miRNAs involved in hormone response only transiently increase during induction. In addition, the establishment of a proliferative status in embryogenic callus is accompanied by important readjustments in hc-siRNAs mapping to long tandem repeat (LTR) retrotransposons, and their expression regulation.

Transcription Factors Responding to Pb Stress in Maize

Pb can damage the physiological function of human organs by entering the human body via food-chain enrichment. Revealing the mechanisms of maize tolerance to Pb is critical for preventing this. In this study, a Pb-tolerant maize inbred line, 178, was used to analyse transcription factors (TFs) expressed under Pb stress based on RNA sequencing data. A total of 464 genes expressed in control check (CK) or Pb treatment samples were annotated as TFs. Among them, 262 differentially expressed transcription factors (DETs) were identified that responded to Pb treatment. Furthermore, the DETs were classified into 4 classes according to their expression patterns, and 17, 12 and 2 DETs were significantly annotated to plant hormone signal transduction, basal transcription factors and base excision repair, respectively. Seventeen DETs were found to participate in the plant hormone signal transduction pathway, where basic leucine zippers (bZIPs) were the most significantly enriched TFs, with 12 members involved. We further obtained 5 Arabidopsis transfer DNA (T-DNA) mutants for 6 of the maize bZIPs, among which the mutants atbzip20 and atbzip47, representing ZmbZIP54 and ZmbZIP107, showed obviously inhibited growth of roots and above-ground parts, compared with wild type. Five highly Pb-tolerant and 5 highly Pb-sensitive in maize lines were subjected to DNA polymorphism and expression level analysis of ZmbZIP54 and ZmbZIP107. The results suggested that differences in bZIPs expression partially accounted for the differences in Pb-tolerance among the maize lines. Our results contribute to the understanding of the molecular regulation mechanisms of TFs in maize under Pb stress.

Identification and characterization of microRNAs in maize endosperm response to exogenous sucrose using small RNA sequencing

Sucrose acts as a signaling molecule for genes critical to starch biosynthesis in maize endosperm. Previously, we showed that sucrose could regulate starch biosynthesis in maize via transcription factors. To better understand the complex regulation of starch biosynthesis, the 10days after pollination endosperm from Zea mays L. B73 inbred line was collected and treated with sucrose for small RNA sequencing. The sequencing results revealed that 24 known miRNAs and 190 novel miRNAs were significantly differentially expressed in response to sucrose. In addition, most of target mRNAs were characterized as transcription factors, mainly including, MYB, ARF, NAC, AP2/ERF, WRKY, and GRAS, which play important roles in starch biosynthesis and seed development in maize endosperm. The expression profiles of miR398a/b and miR159b/j/k followed opposite expression trends to their target genes when analyzed by qPCR. In conclusion, these results show that sucrose regulates the expression of starch synthetic genes through miRNAs.