ZmHMA3, a Member of the Heavy-Metal-Transporting ATPase Family, Regulates Cd and Zn Tolerance in Maize.

The pollution of heavy metals is extremely serious in China, including zinc (Zn), copper (Cu), lead (Pb), and cadmium (Cd). Heavy-metal-transporting ATPase (HMA) belongs to a subfamily of the P-ATPase family, which absorbs and transports Zn, Cu, Pb, and Cd in plants. Here, we describe a ZmHMA-encoding HMA family protein that positively regulates Cd and Zn tolerance. The real-time fluorescence quantification (RT-PCR) results revealed that ZmHMA3 had a high expression in B73, and the expression of ZmHMA3 was sensitive to Cd in yeast cells, which was related to Cd accumulation in yeast. Additionally, the Arabidopsis thaliana homologous mutants of AtHMA2 showed Cd sensitivity compared with WT. The overexpressing ZmHMA3 plants showed higher tolerance under Cd and Zn stresses than the wild type. The overexpression of ZmHMA3 led to higher Cd and Zn accumulation in tissues based on the subcellular distribution analysis. We propose that ZmHMA3 improves maize tolerance to Cd and Zn stresses by absorbing and transporting Cd and Zn ions. This study elucidates the gene function of the ZmHMA3 response to Cd and Zn stress and provides a reference for improving the characteristics of heavy metals enrichment in existing maize varieties and the plant remediation technology of heavy-metal-contaminated soil.

Identification and Fine Mapping of the Recessive Gene BK-5, Which Affects Cell Wall Biosynthesis and Plant Brittleness in Maize.

The cellulose of the plant cell wall indirectly affects the cell shape and straw stiffness of the plant. Here, the novel brittleness mutant brittle stalk-5 (bk-5) of the maize inbred line RP125 was characterized. We found that the mutant displayed brittleness of the stalk and even the whole plant, and that the brittleness phenotype existed during the whole growth period from germination to senescence. The compressive strength was reduced, the cell wall was thinner, and the cellulose content was decreased compared to that of the wild type. Genetic analysis and map-based cloning indicated that bk-5 was controlled by a single recessive nuclear gene and that it was located in a 90.2-Kb region on chromosome 3 that covers three open reading frames (ORFs). Sequence analysis revealed a single non-synonymous missense mutation, T-to-A, in the last exon of Zm00001d043477 (B73: version 4, named BK-5) that caused the 951th amino acid to go from leucine to histidine. BK-5 encodes a cellulose synthase catalytic subunit (CesA), which is involved with cellulose synthesis. We found that BK-5 was constitutively expressed in all tissues of the germinating stage and silking stage, and highly expressed in the leaf, auricula, and root of the silking stage and the 2-cm root and bud of the germinating stage. We found that BK-5 mainly localized to the Golgi apparatus, suggesting that the protein might move to the plasma membrane with the aid of Golgi in maize. According to RNA-seq data, bk-5 had more downregulated genes than upregulated genes, and many of the downregulated genes were enzymes and transcription factors related to cellulose, hemicellulose, and lignin biosynthesis of the secondary cell wall. The other differentially expressed genes were related to metabolic and cellular processes, and were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and the plant-pathogen interaction pathway. Taken together, we propose that the mutation of gene BK-5 causes the brittle stalk phenotype and provides important insights into the regulatory mechanism of cellulose biosynthesis and cell wall development in maize.

Genome assembly of the Chinese maize elite inbred line RP125 and its EMS mutant collection provide new resources for maize genetics research and crop improvement.

Maize is an important crop worldwide, as well as a valuable model with vast genetic diversity. Accurate genome and annotation information for a wide range of inbred lines would provide valuable resources for crop improvement and pan-genome characterization. In this study, we generated a high-quality de novo genome assembly (contig N50 of 15.43 Mb) of the Chinese elite inbred line RP125 using Nanopore long-read sequencing and Hi-C scaffolding, which yield highly contiguous, chromosome-length scaffolds. Global comparison of the RP125 genome with those of B73, W22, and Mo17 revealed a large number of structural variations. To create new germplasm for maize research and crop improvement, we carried out an EMS mutagenesis screen on RP125. In total, we obtained 5818 independent M2 families, with 946 mutants showing heritable phenotypes. Taking advantage of the high-quality RP125 genome, we successfully cloned 10 mutants from the EMS library, including the novel kernel mutant qk1 (quekou: "missing a small part" in Chinese), which exhibited partial loss of endosperm and a starch accumulation defect. QK1 encodes a predicted metal tolerance protein, which is specifically required for Fe transport. Increased accumulation of Fe and reactive oxygen species as well as ferroptosis-like cell death were detected in qk1 endosperm. Our study provides the community with a high-quality genome sequence and a large collection of mutant germplasm.

Identification and characterization of the glutathione S-Transferase (GST) family in radish reveals a likely role in anthocyanin biosynthesis and heavy metal stress tolerance.

Glutathione S-transferases (GSTs) are a large complex family of enzymes (EC 2.5.1.18) that play vital roles in flavonoid metabolism and plant growth and development and are responsive to heavy metal stress. However, knowledge about GST genes in radish (a vegetable crop with an extraordinary capacity to adapt to heavy metal stresses) is limited. Therefore, it is critical to identify putative candidate GST genes responsible for heavy metal stress tolerance and anthocyanin biosynthesis. In this study, we first identified 82 R. sativus GST (RsGST) genes using various bioinformatic approaches, and their expression profiles were characterized from RNAseq data. These RsGST genes could be grouped into 7 major subclasses: tau (43 members), phi (21 members), tetrachlorohydroquinone dehalogenase (7 members), dehydroascorbat reductase (5 members), zeta (3 members), lambda (2 members) and theta (1 member). In addition, most of the RsGST genes showed organ-specific expression in our study. Moreover, the transcripts of RsGSTF12-1 and RsGSTF12-2, belonging to the phi class, might be candidates encoding anthocyanin transporters in carmine radish, whereas the tau class, consisting of RsGSTU13-1, RsGSTU19, RsGSTU24-1, and RsGSTU3, and theta class, consisting of RsGSTT1-1, might be defend radish against adverse heavy metal stresses. These results will aid in understanding the functions of the GST family related to heavy metal stress and anthocyanin biosynthesis, thereby potentially improving radish breeding programs for high-pigment-content material as well as HM-tolerant material.

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping.

Kernel size-related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10-6 ) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10-3 ) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable-effect SNPs and the co-localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co-localized SNPs, of which zma-miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma-miR164e resulted in the down-regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker-assisted selection (MAS) for high-yield breeding in maize.

Combined GWAS and QTL analysis for dissecting the genetic architecture of kernel test weight in maize.

Kernel weight in a unit volume is referred to as kernel test weight (KTW) that directly reflects maize (Zea mays L.) grain quality. In this study, an inter-mated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population and an association panel were used to identify loci responsible for KTW of maize across multiple environments. A total of 18 significant KTW-related single-nucleotide polymorphisms (SNPs) were identified using genome-wide association study (GWAS); they were closely linked to 12 candidate genes. In the IBM Syn10 DH population, linkage analysis detected 19 common quantitative trait loci (QTL), five of which were repeatedly detected among multiple environments. Several verified genes that regulate maize seed development were found in the confidence intervals of the mapped QTL and the LD regions of GWAS, such as ZmYUC1, BAP2, ZmTCRR-1, dek36 and ZmSWEET4c. Combined QTL mapping and GWAS identified one significant SNP that was co-identified in the both populations. Based on the co-localized SNP across the both populations, 17 candidate genes were identified. Of them, Zm00001d044075, Zm00001d044086, and Zm00001d044081 were further identified by candidate gene association study for KTW. Zm00001d044081 encodes homeobox-leucine zipper protein ATHB-4, which has been demonstrated to control apical embryo development in Arabidopsis. Our findings provided insights into the mechanism underlying maize KTW and contributed to the application of molecular-assisted selection of high KTW breeding in maize.

Genome-wide identification of ZmHMAs and association of natural variation in ZmHMA2 and ZmHMA3 with leaf cadmium accumulation in maize.

P1B-type ATPases, known as heavy metal ATPases (HMAs), play an important role in the control of cadmium (Cd) accumulation in plants. In this study, a total of 12 ZmHMA genes were identified in the maize genome and particularly classified into six clusters based on their phylogenetic relationship and motif compositions. Furthermore, the expression patterns of different ZmHMA genes varied with developmental stages, and were tissue specific under normal conditions. ZmHMA2 and ZmHMA3 genes exhibited significant up-regulation under Cd treatment. Eventually, the association analysis between 103 inbred lines and alleles in ZmHMA2 and ZmHMA3 revealed that one insertion-deletion (InDel) in the intron from ZmHMA2 was associated with leaf Cd concentration under low Cd condition at the seedling stage. Twenty polymorphisms in ZmHMA3 were significantly associated with leaf Cd concentration under various Cd levels at seedling and maturing stages. Five single nucleotide polymorphisms (SNPs) and two InDels of these significantly associated polymorphic loci from ZmHMA3 caused the amino acid substitutions and insertion or deletion events. Importantly, the proteins encoded by ZmHMA2 and ZmHMA3 genes were located in the plasma membrane. This comprehensive analysis will provide an important theoretical basis for future functional verification of ZmHMA genes to unravel the mechanisms of Cd accumulation in leaves of maize. Additionally, the favorable alleles in ZmHMA3 will lay a foundation for the marker-assisted selection of low Cd accumulation in maize.

Genetic dissection of stalk lodging-related traits using an IBM Syn10 DH population in maize across three environments (Zea mays L.).

Stalk lodging severely limits the grain yield of maize (Zea mays L.). Mechanical stalk strength can be reflected by the traits of stalk diameter (SD), stalk bending strength (SBS), and lodging rind penetrometer resistance (RPR). To determine the genetic basis of maize stalk lodging, quantitative trait loci (QTLs) were mapped for these three traits using the IBM Syn10 DH population in three environments. The results indicated that there were strong genetic correlations among the three traits, and the analyses of phenotypic variations for SD, SBS, and RPR across the three environments showed high broad-sense heritability (0.6843, 0.5175, and 0.7379, respectively). In total, 44 significant QTLs were identified control the above traits across the 3 environments. A total of 14, 14, and 16 QTLs were identified for SD, SBS, and RPR across single-environment mapping, respectively. Notably, ten QTLs were stably expressed across multiple-environments, including two QTLs for SD, three for SBS, and five for RPR. Three major QTLs each accounting for over 10% of the phenotypic variation were qSD6-2 (10.03%), qSD8-2 (13.73%), and qSBS1-2 (11.89%). Comprehensive analysis of all QTLs in this study revealed that 5 QTL clusters including 12 QTLs were located on chromosomes 1, 3, 7, and 8, respectively. Among these 44 QTLs, 9 harbored 13 stalk lodging-associated SNPs that were detected by our recently published work, with 1 SNP successfully validated in the IBM Syn10 DH population. These chromosomal regions will be useful for marker-assisted selection and fine mapping of stalk lodging-related traits in maize.

Metabolite profiling and genome-wide association studies reveal response mechanisms of phosphorus deficiency in maize seedling.

Inorganic phosphorus (Pi) is an essential element in numerous metabolic reactions and signaling pathways, but the molecular details of these pathways remain largely unknown. In this study, metabolite profiles of maize (Zea mays L.) leaves and roots were compared between six low-Pi-sensitive lines and six low-Pi-tolerant lines under Pi-sufficient and Pi-deficient conditions to identify pathways and genes associated with the low-Pi stress response. Results showed that under Pi deprivation the concentrations of nucleic acids, organic acids and sugars were increased, but that the concentrations of phosphorylated metabolites, certain amino acids, lipid metabolites and nitrogenous compounds were decreased. The levels of secondary metabolites involved in plant immune reactions, including benzoxazinoids and flavonoids, were significantly different in plants grown under Pi-deficient conditions. Among them, the 11 most stable metabolites showed significant differences under low- and normal-Pi conditions based on the coefficient of variation (CV). Isoleucine and alanine were the most stable metabolites for the identification of Pi-sensitive and Pi-resistant maize inbred lines. With the significant correlation between morphological traits and metabolites, five low-Pi-responding consensus genes associated with morphological traits and simultaneously involved in metabolic pathways were mined by combining metabolites profiles and genome-wide association study (GWAS). The consensus genes induced by Pi deficiency in maize seedlings were also validated by reverse-transcription quantitative polymerase chain reaction (RT-qPCR). Moreover, these genes were further validated in a recombinant inbred line (RIL) population, in which the glucose-6-phosphate-1-epimerase encoding gene mediated yield and correlated traits to phosphorus availability. Together, our results provide a framework for understanding the metabolic processes underlying Pi-deficient responses and give multiple insights into improving the efficiency of Pi use in maize.

Genome-wide association analysis and QTL mapping reveal the genetic control of cadmium accumulation in maize leaf.

BACKGROUND: Accumulation of cadmium (Cd) in maize (Zea mays L.) poses a significant risk to human health as it is ingested via the food chain. A genome-wide association study (GWAS) was conducted in a population of 269 maize accessions with 43,737 single nucleotide polymorphisms (SNPs) to identify candidate genes and favorable alleles for controlling Cd accumulation in maize. RESULTS: When grown in contaminated soil, accessions varied significantly in leaf Cd concentration at both the seeding and maturing stages with phenotypic variation and the coefficient of variation all above 48%. The co-localized region between SYN27837 (147,034,650 bp) and SYN36598 (168,551,327 bp) on chromosome 2 was associated with leaf Cd under three soil conditions varying in Cd content in 2015 and 2016. The significant SNP (SYN25051) at position 161,275,547 could explained 27.1% of the phenotype variation. Through QTL mapping using the IBMSyn10 double haploid (DH) population, we validated the existence of a major QTL identified by GWAS; qLCd2 could explain the 39.8% average phenotype variation across the experiments. Expression of GRMZM2G175576 encoding a cadmium/zinc-transporting ATPase underlying the QTL was significantly increased in roots, stems and leaves of B73, a low Cd accumulation line in response to Cd stress. CONCLUSIONS: Our findings provide new insights into the genetic control of Cd accumulation and could aid rapid development of maize genotypes with low-Cd accumulation by manipulation of the favorable alleles.

Genome-wide association analysis of lead accumulation in maize.

Large phenotypic variations in the lead (Pb) concentration were observed in grains and leaves of maize plants. A further understanding of inheritance of Pb accumulation may facilitate improvement of low-Pb-accumulating cultivars in maize. A genome-wide association study was conducted in a population of 269 maize accessions with 43,737 single-nucleotide polymorphisms (SNPs). The Pb concentrations in leaves and kernels of 269 accessions were collected in pot-culture and field experiments in years of 2015 and 2016. Significant differences in Pb accumulation were found among individuals under different environments. Using the structure and kinship model, a total of 21 SNPs significantly associated with the Pb accumulation were identified with P < 2.28 × 10-5 and FDR < 0.05 in the pot-culture and field experiments across 2 years. Three SNPs on chromosome 4 had significant associations simultaneously with the Pb concentrations of kernels and leaves and were co-localized with the previously detected quantitative trait loci. Through ridge regression best linear unbiased prediction Pb accumulation in the association population, the prediction accuracies by cross validation were 0.18-0.59 and 0.17-0.64, depending on the k-fold and the size of the training population. The results are helpful for genetic improvement and genomic prediction of Pb accumulation in maize.

Integrative analysis of DNA methylation, mRNAs, and small RNAs during maize embryo dedifferentiation.

BACKGROUND: Maize (Zea mays) is an important model crop for transgenic studies. However, genetic transformation of maize requires embryonic calli derived from immature embryo, and the impact of utilizing tissue culture methods on the maize epigenome is poorly understood. Here, we generated whole-genome MeDIP-seq data examining DNA methylation in dedifferentiated and normal immature maize embryos. RESULTS: We observed that most of the dedifferentiated embryos exhibited a methylation increase compared to normal embryos. Increased methylation at promoters was associated with down-regulated protein-coding gene expression; however, the correlation was not strong. Analysis of the callus and immature embryos indicated that the methylation increase was induced during induction of embryonic callus, suggesting phenotypic consequences may be caused by perturbations in genomic DNA methylation levels. The correlation between the 21-24nt small RNAs and DNA methylation regions were investigated but only a statistically significant correlation for 24nt small RNAs was observed. CONCLUSIONS: These data extend the significance of epigenetic changes during maize embryo callus formation, and the methylation changes might explain some of the observed embryonic callus variation in callus formation.

Molecular Evolution and Association of Natural Variation in ZmARF31 with Low Phosphorus Tolerance in Maize.

Low-phosphorus (P) stress is one of the major factors constraining plant growth and yield. Improving plant tolerance to P starvation through molecular breeding is an efficient alternative to increase grain production. In the study, 331 diverse maize inbreds were used to detect nucleotide diversity and favorable alleles of ZmARF31, which plays a key role in low P responses and root architecture regulation. Significant phenotypic variation was found in each of 11 tested traits under both P and no-P treatments, and 30 single nucleotide polymorphisms (SNPs) and 14 insertion-deletions (InDels) were detected in ZmARF31 among the 331 maize inbreds. The 5'-untranslated region (UTR) of ZmARF31 showed a small linkage disequilibrium (LD) block under significant purifying selection, whereas the 3'-UTR showed the most abundant diversity and a larger LD block. Thirty, fourteen, and nine natural variations were identified in ZmARF31 that were associated with P-deficiency-tolerance traits (P ≤ 0.01) by using the general linear model (GLM), GLM incorporated with population structure, and mixed linear model, respectively. Four SNPs were significantly associated with the total dry weight (TDW) in the three models, of which SNPs S410 and S462 were located in a complete LD block. A further verification conducted in a recombinant inbred line population revealed that favorable allele G/G of non-synonymous mutation S410 and favorable allele with a 38 bp insertion of InDel S1442 exhibited positive genetic effects on the TDW and total root tips, respectively. Expression analysis further confirmed that ZmARF31 was highly expressed in the roots of low-P-tolerant inbred 178. The protein encoded by ZmARF31 was located both in the nucleus and cytoplasm. Haplotypes carrying more favorable alleles showed a greater effect on phenotypic variation than single loci. Such haplotypes should be helpful to develop valuable genetic markers and perform maize molecular breeding.

Correction: The Impact of Genetic Relationship and Linkage Disequilibrium on Genomic Selection.

[This corrects the article DOI: 10.1371/journal.pone.0132379.].

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

In this study, a maize inbred line with a strong capacity to induce embryonic callus, 18-599R, was used to analyze the transcription factors expressed during embryonic callus formation. A total of 1180 transcription factors were found to be expressed during three key stages of callus induction. Of these, compared with control, 361, 346 and 328 transcription factors were significantly downregulated during stages I, II and III, respectively. In contrast, 355, 372 and 401 transcription factors (TFs) were upregulated during the respective stages. We constructed a transcription factor-mediated regulatory network and found that plant hormone signal transduction was the pathway most significantly enriched among TFs. This pathway includes 48 TFs regulating cell enlargement, cell differentiation, cell division and cell dedifferentiation via the response to plant hormones. Through real-time polymerase chain reaction (PCR) and degradome sequencing, we identified 23 transcription factors that are regulated by miRNA. Through further analysis, ZmMYB138, a member of the MYB transcription factor family localized in the nucleus, was verified to promote embryonic callus formation in the maize embryo through GA signal transduction.

Comparative RNA-Seq Analysis Reveals That Regulatory Network of Maize Root Development Controls the Expression of Genes in Response to N Stress.

Nitrogen (N) is an essential nutrient for plants, and it directly affects grain yield and protein content in cereal crops. Plant root systems are not only critical for anchorage in the soil, but also for N acquisition. Therefore, genes controlling root development might also affect N uptake by plants. In this study, the responses of nitrogen on root architecture of mutant rtcs and wild-type of maize were investigated by morphological and physiological analysis. Subsequently, we performed a comparative RNA-Seq analysis to compare gene expression profiles between mutant rtcs roots and wild-type roots under different N conditions. We identified 786 co-modulated differentially expressed genes (DEGs) related to root development. These genes participated in various metabolic processes. A co-expression cluster analysis and a cis-regulatory motifs analysis revealed the importance of the AP2-EREBP transcription factor family in the rtcs-dependent regulatory network. Some genotype-specific DEGs contained at least one LBD motif in their promoter region. Further analyses of the differences in gene transcript levels between rtcs and wild-type under different N conditions revealed 403 co-modulated DEGs with distinct functions. A comparative analysis revealed that the regulatory network controlling root development also controlled gene expression in response to N-deficiency. Several AP2-EREBP family members involved in multiple hormone signaling pathways were among the DEGs. These transcription factors might play important roles in the rtcs-dependent regulatory network related to root development and the N-deficiency response. Genes encoding the nitrate transporters NRT2-1, NAR2.1, NAR2.2, and NAR2.3 showed much higher transcript levels in rtcs than in wild-type under normal-N conditions. This result indicated that the LBD gene family mainly functions as transcriptional repressors, as noted in other studies. In summary, using a comparative RNA-Seq-based approach, we identified DEGs related to root development that also participated in the N-deficiency response in maize. These findings will increase our understanding of the molecular regulatory networks controlling root development and N-stress responses.

An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize.

BACKGROUND: To safeguard the food supply for the growing human population, it is important to understand and exploit the genetic basis of quantitative traits. Next-generation sequencing technology performs advantageously and effectively in genetic mapping and genome analysis of diverse genetic resources. Hence, we combined re-sequencing technology and a bin map strategy to construct an ultra-high-density bin map with thousands of bin markers to precisely map a quantitative trait locus. RESULTS: In this study, we generated a linkage map containing 1,151,856 high quality SNPs between Mo17 and B73, which were verified in the maize intermated B73 × Mo17 (IBM) Syn10 population. This resource is an excellent complement to existing maize genetic maps available in an online database  (iPlant, http://data.maizecode.org/maize/qtl/syn10/ ). Moreover, in this population combined with the IBM Syn4 RIL population, we detected 135 QTLs for flowering time and plant height traits across the two populations. Eighteen known functional genes and twenty-five candidate genes for flowering time and plant height trait were fine-mapped into a 2.21-4.96 Mb interval. Map expansion and segregation distortion were also analyzed, and evidence for inadvertent selection of early flowering time in the process of mapping population development was observed. Furthermore, an updated integrated map with 1,151,856 high-quality SNPs, 2,916 traditional markers and 6,618 bin markers was constructed. The data were deposited into the iPlant Discovery Environment (DE), which provides a fundamental resource of genetic data for the maize genetic research community. CONCLUSIONS: Our findings provide basic essential genetic data for the maize genetic research community. An updated IBM Syn10 population and a reliable, verified high-quality SNP set between Mo17 and B73 will aid in future molecular breeding efforts.

The Impact of Genetic Relationship and Linkage Disequilibrium on Genomic Selection.

Genomic selection is a promising research area due to its practical application in breeding. In this study, impact of realized genetic relationship and linkage disequilibrium (LD) on marker density and training population size required was investigated and their impact on practical application was further discussed. This study is based on experimental data of two populations derived from the same two founder lines (B73, Mo17). Two populations were genotyped with different marker sets at different density: IBM Syn4 and IBM Syn10. A high-density marker set in Syn10 was imputed into the Syn4 population with low marker density. Seven different prediction scenarios were carried out with a random regression best linear unbiased prediction (RR-BLUP) model. The result showed that the closer the real genetic relationship between training and validation population, the fewer markers were required to reach a good prediction accuracy. Taken the short-term cost for consideration, relationship information is more valuable than LD information. Meanwhile, the result indicated that accuracies based on high LD between QTL and markers were more stable over generations, thus LD information would provide more robust prediction capacity in practical applications.

Transcriptomic changes during maize roots development responsive to Cadmium (Cd) pollution using comparative RNAseq-based approach.

The heavy metal cadmium (Cd), acts as a widespread environmental contaminant, which has shown to adversely affect human health, food safety and ecosystem safety in recent years. However, research on how plant respond to various kinds of heavy metal stress is scarcely reported, especially for understanding of complex molecular regulatory mechanisms and elucidating the gene networks of plant respond to Cd stress. Here, transcriptomic changes during Mo17 and B73 seedlings development responsive to Cd pollution were investigated and comparative RNAseq-based approach in both genotypes were performed. 115 differential expression genes (DEGs) with significant alteration in expression were found co-modulated in both genotypes during the maize seedling development; of those, most of DGEs were found comprised of stress and defense responses proteins, transporters, as well as transcription factors, such as thaumatin-like protein, ZmOPR2 and ZmOPR5. More interestingly, genotype-specific transcriptional factors changes induced by Cd stress were found contributed to the regulatory mechanism of Cd sensitivity in both different genotypes. Moreover, 12 co-expression modules associated with specific biological processes or pathways (M1 to M12) were identified by consensus co-expression network. These results will expand our understanding of complex molecular mechanism of response and defense to Cd exposure in maize seedling roots.

Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development.

The present study profiled and analyzed gene expression of the maize ear at four key developmental stages. Based on genome-wide profile analysis, we detected differential mRNA of maize genes. Some of the differentially expressed genes (DEGs) were predicted to be potential candidates of maize ear development. Several well-known genes were found with reported mutant analyses, such as, compact plant2 (ct2), zea AGAMOUS homolog1 (zag1), bearded ear (bde), and silky1 (si1). MicroRNAs such as microRNA156 were predicted to target genes involved in maize ear development. Antisense transcripts were widespread throughout all the four stages, and are suspected to play important roles in maize ear development. Thus, identification and characterization of important genes and regulators at all the four developmental stages will contribute to an improved understanding of the molecular mechanisms responsible for maize ear development.