Tingzhao Rong

Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.

Joint linkage-linkage disequilibrium mapping is a powerful approach to detecting quantitative trait loci underlying drought tolerance in maize.

This paper describes two joint linkage–linkage disequilibrium (LD) mapping approaches: parallel mapping (independent linkage and LD analysis) and integrated mapping (datasets analyzed in combination). These approaches were achieved using 2,052 single nucleotide polymorphism (SNP) markers, including 659 SNPs developed from drought-response candidate genes, screened across three recombinant inbred line (RIL) populations and 305 diverse inbred lines, with anthesis-silking interval (ASI), an important trait for maize drought tolerance, as the target trait. Mapping efficiency was improved significantly due to increased population size and allele diversity and balanced allele frequencies. Integrated mapping identified 18 additional quantitative trait loci (QTL) not detected by parallel mapping. The use of haplotypes improved mapping efficiency, with the sum of phenotypic variation explained (PVE) increasing from 5.4% to 23.3% for single SNP-based analysis. Integrated mapping with haplotype further improved the mapping efficiency, and the most significant QTL had a PVE of up to 34.7%. Normal allele frequencies for 113 of 277 (40.8%) SNPs with minor allele frequency (<5%) in 305 lines were recovered in three RIL populations, three of which were significantly associated with ASI. The candidate genes identified by two significant haplotype loci included one for a SET domain protein involved in the control of flowering time and the other encoding aldo/keto reductase associated with detoxification pathways that contribute to cellular damage due to environmental stress. Joint linkage–LD mapping is a powerful approach for detecting QTL underlying complex traits, including drought tolerance.

Identification of candidate genes for drought tolerance by whole-genome resequencing in maize

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.

Comparative SNP and Haplotype Analysis Reveals a Higher Genetic Diversity and Rapider LD Decay in Tropical than Temperate Germplasm in Maize

Understanding of genetic diversity and linkage disequilibrium (LD) decay in diverse maize germplasm is fundamentally important for maize improvement. A total of 287 tropical and 160 temperate inbred lines were genotyped with 1943 single nucleotide polymorphism (SNP) markers of high quality and compared for genetic diversity and LD decay using the SNPs and their haplotypes developed from genic and intergenic regions. Intronic SNPs revealed a substantial higher variation than exonic SNPs. The big window size haplotypes (3-SNP slide-window covering 2160 kb on average) revealed much higher genetic diversity than the 10 kb-window and gene-window haplotypes. The polymorphic information content values revealed by the haplotypes (0.436–0.566) were generally much higher than individual SNPs (0.247–0.259). Cluster analysis classified the 447 maize lines into two major groups, corresponding to temperate and tropical types. The level of genetic diversity and subpopulation structure were associated with the germplasm origin and post-domestication selection. Compared to temperate lines, the tropical lines had a much higher level of genetic diversity with no significant subpopulation structure identified. Significant variation in LD decay distance (2–100 kb) was found across the genome, chromosomal regions and germplasm groups. The average of LD decay distance (10–100 kb) in the temperate germplasm was two to ten times larger than that in the tropical germplasm (5–10 kb). In conclusion, tropical maize not only host high genetic diversity that can be exploited for future plant breeding, but also show rapid LD decay that provides more opportunity for selection.

Comparative LD mapping using single SNPs and haplotypes identifies QTL for plant height and biomass as secondary traits of drought tolerance in maize

Drought often delays developmental events so that plant height and above-ground biomass are reduced, resulting in yield loss due to inadequate photosynthate. In this study, plant height and biomass measured by the Normalized Difference Vegetation Index (NDVI) were used as criteria for drought tolerance. A total of 305 lines representing temperate, tropical and subtropical maize germplasm were genotyped using two single nucleotide polymorphism (SNP) chips each containing 1536 markers, from which 2052 informative SNPs and 386 haplotypes each constructed with two or more SNPs were used for linkage disequilibrium (LD) or association mapping. Single SNP- and haplotype-based LD mapping identified two significant SNPs and three haplotype loci [a total of four quantitative trait loci (QTL)] for plant height under well-watered and water-stressed conditions. For biomass, 32 SNPs and 12 haplotype loci (30 QTL) were identified using NDVIs measured at seven stages under the two water regimes. Some significant SNP and haplotype loci for NDVI were shared by different stages. Comparing significant loci identified by single SNP- and haplotype-based LD mapping, we found that six out of the 14 chromosomal regions defined by haplotype loci each included at least one significant SNP for the same trait. Significant SNP haplotype loci explained much higher phenotypic variation than individual SNPs. Moreover, we found that two significant SNPs (two QTL) and one haplotype locus were shared by plant height and NDVI. The results indicate the power of comparative LD mapping using single SNPs and SNP haplotypes with QTL shared by plant height and biomass as secondary traits for drought tolerance in maize.

Validation of Potential Reference Genes for qPCR in Maize across Abiotic Stresses, Hormone Treatments, and Tissue Types

The reverse transcription quantitative polymerase chain reaction (RT-qPCR) is a powerful and widely used technique for the measurement of gene expression. Reference genes, which serve as endogenous controls ensure that the results are accurate and reproducible, are vital for data normalization. To bolster the literature on reference gene selection in maize, ten candidate reference genes, including eight traditionally used internal control genes and two potential candidate genes from our microarray datasets, were evaluated for expression level in maize across abiotic stresses (cold, heat, salinity, and PEG), phytohormone treatments (abscisic acid, salicylic acid, jasmonic acid, ethylene, and gibberellins), and different tissue types. Three analytical software packages, geNorm, NormFinder, and Bestkeeper, were used to assess the stability of reference gene expression. The results revealed that elongation factor 1 alpha (EF1α), tubulin beta (β-TUB), cyclophilin (CYP), and eukaryotic initiation factor 4A (EIF4A) were the most reliable reference genes for overall gene expression normalization in maize, while GRP (Glycine-rich RNA-binding protein), GLU1(beta-glucosidase), and UBQ9 (ubiquitin 9) were the least stable and most unsuitable genes. In addition, the suitability of EF1α, β-TUB, and their combination as reference genes was confirmed by validating the expression of WRKY50 in various samples. The current study indicates the appropriate reference genes for the urgent requirement of gene expression normalization in maize across certain abiotic stresses, hormones, and tissue types.

Development and Characterization of Simple Sequence Repeat Markers Providing Genome-Wide Coverage and High Resolution in Maize

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.

Cytoplasmic male sterility-regulated novel microRNAs from maize

In higher plants, microRNA (miRNA) is involved in regulation of developmental processes, including sexual organ development. Seven novel miRNA families with one known miRNA were isolated by constructing a small RNA library from a mixture of anther from a cytoplasmic male sterile line and its maintainer. Two miRNAs are conserved in plant species. A total of 18 potential targets were identified for the eight miRNA families, including 15 proteins annotated with function and three unknown proteins. The known proteins include several proteins relevant to cell structure and stress response, transcription factors, and enzymes associated with metabolic and signaling pathways, playing important roles in microspore development. Quantitative real-time PCR assay revealed different expression patterns of the miRNAs between the cytoplasmic male sterile line and its maintainer. Each of the miRNAs tended to be down-regulated after the tetrad stage in a fertile line. However, most of the miRNAs in the cytoplasmic male sterile line were shown to be up-regulated from the tetrad to mononuclear stage, displaying special expression patterns differing from the ones in fertile line. We conclude that additional inactive miRNA pathways are essential during pollen development for a fertile line to ensure male fertility. Contrarily, miRNAs are up-regulated during the period from the tetrad to mononuclear stage, which contributes to pollen abortion for a cytoplasmic male sterile line.

Combined small RNA and degradome sequencing reveals microRNA regulation during immature maize embryo dedifferentiation

Genetic transformation of maize is highly dependent on the development of embryonic calli from the dedifferentiated immature embryo. To better understand the regulatory mechanism of immature embryo dedifferentiation, we generated four small RNA and degradome libraries from samples representing the major stages of dedifferentiation. More than 186 million raw reads of small RNA and degradome sequence data were generated. We detected 102 known miRNAs belonging to 23 miRNA families. In total, we identified 51, 70 and 63 differentially expressed miRNAs (DEMs) in the stage I, II, III samples, respectively, compared to the control. However, only 6 miRNAs were continually up-regulated by more than fivefold throughout the process of dedifferentiation. A total of 87 genes were identified as the targets of 21 DEM families. This group of targets was enriched in members of four significant pathways including plant hormone signal transduction, antigen processing and presentation, ECM-receptor interaction, and alpha-linolenic acid metabolism. The hormone signal transduction pathway appeared to be particularly significant, involving 21 of the targets. While the targets of the most significant DEMs have been proved to play essential roles in cell dedifferentiation. Our results provide important information regarding the regulatory networks that control immature embryo dedifferentiation in maize.

The genetic architecture of flowering time and photoperiod sensitivity in maize as revealed by QTL review and meta analysis.

The control of flowering is not only important for reproduction, but also plays a key role in the processes of domestication and adaptation. To reveal the genetic architecture for flowering time and photoperiod sensitivity, a comprehensive evaluation of the relevant literature was performed and followed by meta analysis. A total of 25 synthetic consensus quantitative trait loci (QTL) and four hot-spot genomic regions were identified for photoperiod sensitivity including 11 genes related to photoperiod response or flower morphogenesis and development. Besides, a comparative analysis of the QTL for flowering time and photoperiod sensitivity highlighted the regions containing shared and unique QTL for the two traits. Candidate genes associated with maize flowering were identified through integrated analysis of the homologous genes for flowering time in plants and the consensus QTL regions for photoperiod sensitivity in maize (Zea mays L.). Our results suggest that the combination of literature review, meta-analysis and homologous blast is an efficient approach to identify new candidate genes and create a global view of the genetic architecture for maize photoperiodic flowering. Sequences of candidate genes can be used to develop molecular markers for various models of marker-assisted selection, such as marker-assisted recurrent selection and genomic selection that can contribute significantly to crop environmental adaptation.