The main producing areas of maize is mostly located on the arid or semi-arid region that relying on the rainfed farming in China. The maize production losses caused by drought is a great threaten to food security. As a cross-pollinating crop, maize is mostly sensitive to water stress during flowering time. Drought at flowering stage will lead to asynchronous development between the male and female flower and cause massive grain yield loss. Thus, mining drought resistance related genes at flowering stage is important for maize drought resistance improvement and breeding.

In the present study, the phylogenic tree of 24 ZCN genes in maize genome, which is homologs of Arabidopsis FT gene, was build. The gene expression patterns of ZCN7 were analysis using qRT-PCR and in vivo GFP fluorescence imaging. A maize natural population consisting of 118 diverse inbred lines were planted in three environments, Beijing in 2021 and 2022 and Urumqi in 2022, to identify the flowering time related traits under different water treatments. The genomic variants around ZCN7 were detected by PCR and Sanger sequencing. The candidate gene association analysis was performed based on mixed linear model and the significant associated variants with drought induced anthesis-silking interval was obtained. The gene expression level of ZCN7 in natural population at flowering time was also measured by qRT-PCR. The differences of drought resistance traits and ZCN7 expression were compared between different haplotypes of significant associated variant. The Ubi1:ZCN7 overexpression transgenic maize were obtained, and the phenotypic performance was identified under different water treatments.

The 24 ZCN genes in maize genome included 15 FT like genes, 6 TFL1 like genes and 3 MFT like genes. The protein sequence of ZCN genes varied from 111 nn to 193 nn. The ZCN7 showed close relationship with ZCN8 and the protein sequence identity was 83.3% between the two genes. ZCN7 showed highest gene expression in the leaf blade at V12 stage. And the ZCN7-promoter:GFP vector was transformed to Arabidopsis and the GFP showed enriched signal at the blade edge of mature leaf. The candidate gene association analysis revealed a SNP variant at 1001 bp upstream of ZCN7 start codon had highest association signal with drought induced anthesis-silking interval under drought. The A/A and G/G haplotypes of SNP-1001 included 78 and 27 inbred lines, respectively. The anthesis-silking interval of A/A haplotype lines were significantly lower than G/G lines. And the ZCN7 gene expression of A/A haplotype lines were significantly higher than G/G lines. In addition, the ZCN7 overexpression transgenic lines showed significantly decreased anthesis-silking interval than wild type lines. Under drought, the anthesis-silking intervals of OE1 and OE2 were 2.3 and 2.6 days shorter than wild type lines. And the grain yield per plant and kernel number per plant of transgenic lines were significantly higher than wild type lines under drought, while the hundred kernel weight, kernel length and kernel width showed no significant difference.

The maize ZCN7 played positive role in drought resistance and its overexpression improved grain yield by reducing anthesis-silking interval under drought.

Drought is one of the serious abiotic stresses influencing maize production worldwide. Understanding the molecular mechanisms underlying drought tolerance is of great importance in maize improvement.

In the present study, representative maize inbred lines were selected for field drought experiment and the drought tolerance was estimated based on leaf relative water content and anthesis-silking interval. Two inbred lines with contrasting drought tolerance was used for genome resequencing and transposable element insertions were identified. The DNA methylation level of leaf and root tissues under different water treatments of the two lines were measured using Whole Genome Bisulfite Sequencing (WGBS). And the gene expression profiles of these samples were detected by RNA sequencing. The integrative atlas of transposable element insertion/deletion variants, differentially methylated regions and differentially expressed genes in the two lines were constructed. In addition, the transposable element insertion/deletion variant mediated epigenetic regulation of ZCN7, which has been conferred the drought tolerance function in our previous study, was analyzed.

The field experiment showed inbred line H082183 showed highest drought tolerance, in which the leaf relative water content and anthesis-silking interval had no significant difference between drought and well-watered treatments. While the Lü28 displayed lowest leaf relative water content and largest anthesis-silking interval under drought. Thus, these two lines were selected for further analysis. A total of 333 754 and 333 296 transposable element insertions were identified in the genome of H082183 and Lü28, respectively. And 89 954 transposable element insertions were polymorphism between two lines. The transposable element insertions, introns and promoters showed higher CG and CHG methylation level than exons and untranslated regions. Furthermore, 41 352 differentially methylated regions were identified between H082183 and Lü28. And 60% of the differentially methylated regions were located in the transposable element insertion\deletion variants and 5 kb flanking regions. The gene expression level showed negatively correlated with CG and CHG methylation. Differentially expression analysis between H082183 and Lü28 obtained 4 196 and 3 500 differentially expressed genes in leaf and root under drought, respectively. The 19.5% and 19.7% of these differentially expressed genes were located in differentially methylated regions. Three LTR transposable element insertions were identified in the 34 kb region of ZCN7 in Lü28 but absent in the genome H082183. And the DNA methylation levels of CG and CHG in this genomic region were significantly higher in Lü28 than H082183 under both drought and well-watered environments, which conferred higher ZCN7 expression in the drought tolerant line H082183.

Our results highlight the important role of interplay of transposable element insertions, DNA methylation and gene expression under drought. And gene expression regulation mechanism of ZCN7 relied on the transposable element insertion/deletion variants mediated DNA methylation was proposed.