Drought 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. A 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. Our 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.

DNA methylation is important for plant growth, development, and stress response. To understand DNA methylation dynamics in maize roots under water stress (WS), we reanalyzed DNA methylation sequencing data to profile DNA methylation and the gene expression landscape of two inbred lines with different drought sensitivities, as well as two of their derived recombination inbred lines (RILs). Combined with genotyping-by-sequencing, we found that the inheritance pattern of DNA methylation between RILs and parental lines was sequence-dependent. Increased DNA methylation levels were observed under WS and the methylome of drought-tolerant inbred lines were much more stable than that of the drought-sensitive inbred lines. Distinctive differentially methylated genes were found among diverse genetic backgrounds, suggesting that inbred lines with different drought sensitivities may have responded to stress in varying ways. Gene body DNA methylation showed a negative correlation with gene expression but a positive correlation with exon splicing events. Furthermore, a positive correlation of a varying extent was observed between small interfering RNA (siRNA) and DNA methylation, which at different genic regions. The response of siRNAs under WS was consistent with the differential DNA methylation. Taken together, our data can be useful in deciphering the roles of DNA methylation in plant drought-tolerance variations and in emphasizing its function in alternative splicing.

Long noncoding RNAs (lncRNAs) are transcripts with lengths of more than 200 nt and limited protein-coding potential. They were found to play important roles in plant stress responses. In this study, the maize drought-tolerant inbred line AC7643 and drought-sensitive inbred line AC7729/TZSRW, as well as their recombinant inbred lines (RILs) were selected to identify drought-responsive lncRNAs in roots. Compared with non-responsive lncRNAs, drought-responsive lncRNAs had different sequence characteristics in length of genes and number of exons. The ratio of down-regulated lncRNAs induced by drought was significantly higher than that of coding genes; and lncRNAs were more widespread expressed in recombination sites in the RILs. Additionally, by integration of the modifications of DNA 5-methylcytidine (5mC), histones, and RNA N6-methyladenosine (m⁶A), it was found that the enrichment of histone modifications associated with transcriptional activation in the genes generated lncRNAs was lower that coding genes. The lncRNAs-mRNAs co-expression network, containing 15,340 coding genes and 953 lncRNAs, was constructed to investigate the molecular functions of lncRNAs. There are 13 modules found to be associated with survival rate under drought. We found nine SNPs located in lncRNAs among the modules associated with plant survival under drought. In conclusion, we revealed the characteristics of lncRNAs responding to drought in maize roots based on multiomics studies. These findings enrich our understanding of lncRNAs under drought and shed light on the complex regulatory networks that are orchestrated by the noncoding RNAs in response to drought stress.

Fish are rich sources of amino acids (AAs), particularly human essential amino acids (HEAAs). Exploring the regulatory mechanisms behind the changes in the combined AA content in the fillet and enhancing the content of AAs, especially HEAAs, in fillets of farmed fish is crucial for meeting human nutritional needs. After hot acidic hydrolysis of 304 common carp fillets, we quantified the contents of 17 single AAs and 5 AA groups and observed significant variations among them. Except for Pro, 16 single AAs and all AA groups showed medium-to-high heritabilities over 0.2. Through a genome-wide association study (GWAS), we identified 1974 SNPs and candidate genes associated with at least one AA content. Using transcriptome data from groups with the highest and lowest contents for each AA, 7089 candidate genes were related to the concentrations of at least two AAs. For the total HEAA content, 121 SNPs and their associated genes preferred ATPase-coupled transmembrane transporter activity, and 4727 differentially expressed genes were enriched in cytokine activity, chemokine activity, oxidoreductase activity, and ion binding. With the optimal genomic selection programs and associated SNPs, the correlation between the actual AA contents and estimated breeding values was high and positive, ranging from 0.76 to 0.90. These findings revealed the major-effect processes and regulatory mechanisms modulating the differences in fillet AA contents. The genomic selection programs will guide the future selection of common carp with high AA contents.