BackgroundDrought stress is one of the major limiting factors for maize production. With the availability of maize B73 reference genome and whole-genome resequencing of 15 maize inbreds, common variants (CV) and clustering analyses were applied to identify non-synonymous SNPs (nsSNPs) and corresponding candidate genes for drought tolerance.ResultsA total of 524 nsSNPs that were associated with 271 candidate genes involved in plant hormone regulation, carbohydrate and sugar metabolism, signaling molecules regulation, redox reaction and acclimation of photosynthesis to environment were detected by CV and cluster analyses. Most of the nsSNPs identified were clustered in bin 1.07 region that harbored six previously reported QTL with relatively high phenotypic variation explained for drought tolerance. Genes Ontology (GO) analysis of candidate genes revealed that there were 35 GO terms related to biotic stimulus and membrane-bounded organelle, showing significant differences between the candidate genes and the reference B73 background. Changes of expression level in these candidate genes for drought tolerance were detected using RNA sequencing for fertilized ovary, basal leaf meristem tissue and roots collected under drought stressed and well-watered conditions. The results indicated that 70% of candidate genes showed significantly expression changes under two water treatments and our strategies for mining candidate genes are feasible and relatively efficient.ConclusionsOur results successfully revealed candidate nsSNPs and associated genes for drought tolerance by comparative sequence analysis of 16 maize inbred lines. Both methods we applied were proved to be efficient for identifying candidate genes for complex traits through the next-generation sequencing technologies (NGS). These selected genes will not only facilitate understanding of genetic basis of drought stress response, but also accelerate genetic improvement through marker-assisted selection in maize.

Heterochromatin is a tightly packed form of chromatin that is associated with DNA methylation and histone 3 lysine 9 methylation (H3K9me). Here, we identify an H3K9me2-binding protein, Agenet domain (AGD)-containing p1 (AGDP1), in Arabidopsis thaliana. Here we find that AGDP1 can specifically recognize the H3K9me2 mark by its three pairs of tandem AGDs. We determine the crystal structure of the Agenet domain 1 and 2 cassette (AGD12) of Raphanus sativus AGDP1 in complex with an H3K9me2 peptide. In the complex, the histone peptide adopts a unique helical conformation. AGD12 specifically recognizes the H3K4me0 and H3K9me2 marks by hydrogen bonding and hydrophobic interactions. In addition, we find that AGDP1 is required for transcriptional silencing, non-CG DNA methylation, and H3K9 dimethylation at some loci. ChIP-seq data show that AGDP1 preferentially occupies long transposons and is associated with heterochromatin marks. Our findings suggest that, as a heterochromatin-binding protein, AGDP1 links H3K9me2 to DNA methylation in heterochromatin regions.

Simple sequence repeats (SSRs) have been widely used in maize genetics and breeding, because they are co-dominant, easy to score, and highly abundant. In this study, we used whole-genome sequences from 16 maize inbreds and 1 wild relative to determine SSR abundance and to develop a set of high-density polymorphic SSR markers. A total of 264 658 SSRs were identified across the 17 genomes, with an average of 135 693 SSRs per genome. Marker density was one SSR every of 15.48 kb. (C/G)n, (AT)n, (CAG/CTG)n, and (AAAT/ATTT)n were the most frequent motifs for mono, di-, tri-, and tetra-nucleotide SSRs, respectively. SSRs were most abundant in intergenic region and least frequent in untranslated regions, as revealed by comparing SSR distributions of three representative resequenced genomes. Comparing SSR sequences and e-polymerase chain reaction analysis among the 17 tested genomes created a new database, including 111 887 SSRs, that could be develop as polymorphic markers in silico. Among these markers, 58.00, 26.09, 7.20, 3.00, 3.93, and 1.78% of them had mono, di-, tri-, tetra-, penta-, and hexa-nucleotide motifs, respectively. Polymorphic information content for 35 573 polymorphic SSRs out of 111 887 loci varied from 0.05 to 0.83, with an average of 0.31 in the 17 tested genomes. Experimental validation of polymorphic SSR markers showed that over 70% of the primer pairs could generate the target bands with length polymorphism, and these markers would be very powerful when they are used for genetic populations derived from various types of maize germplasms that were sampled for this study.

Although the root system is indispensable for absorption of nutrients and water, it is poorly studied in maize owing to the difficulties of direct measurement of roots. Here, 103 maize lines were used to compare root architectures under well-watered and water-stressed conditions. Significant genetic variation, with medium to high heritability and significant correlations, was observed for root traits. Total root length (TRL) and total root surface area (TSA) had high phenotypical diversity, and TRL was positively correlated with TSA, root volume, and root forks. The first two principal components explained 94.01% and 91.15% of total root variation in well-watered and water-stressed conditions, respectively. Thus, TRL and TSA, major contributors to root variation, can be used as favorable selection criteria at the seedling stage. We found that stiff stalk and non-stiff stalk groups (temperate backgrounds) showed relatively higher mean values for root morphological diversity than the TST group (tropical/subtropical background). Of the tested lines, 7, 42, 45, and 9 were classified as drought sensitive, moderately sensitive, moderately drought tolerant, and highly drought tolerant, respectively. Seven of the 9 extremely drought tolerant lines were from the TST group, suggesting that TST germplasms harbor valuable genetic resources for drought tolerance that could be used in breeding to improve abiotic stress tolerance in maize.