Tianyu Wang

High-resolution genetic mapping of maize pan-genome sequence anchors

In addition to single-nucleotide polymorphisms, structural variation is abundant in many plant genomes. The structural variation across a species can be represented by a ‘pan-genome, which is essential to fully understand the genetic control of phenotypes. However, the pan-genomes complexity hinders its accurate assembly via sequence alignment. Here we demonstrate an approach to facilitate pan-genome construction in maize. By performing 18 trillion association tests we map 26 million tags generated by reduced representation sequencing of 14,129 maize inbred lines. Using machine-learning models we select 4.4 million accurately mapped tags as sequence anchors, 1.1 million of which are presence/absence variations. Structural variations exhibit enriched association with phenotypic traits, indicating that it is a significant source of adaptive variation in maize. The ability to efficiently map ultrahigh-density pan-genome sequence anchors enables fine characterization of structural variation and will advance both genetic research and breeding in many crops.

Genetic Control of the Leaf Angle and Leaf Orientation Value as Revealed by Ultra-High Density Maps in Three Connected Maize Populations

Plant architecture is a key factor for high productivity maize because ideal plant architecture with an erect leaf angle and optimum leaf orientation value allow for more efficient light capture during photosynthesis and better wind circulation under dense planting conditions. To extend our understanding of the genetic mechanisms involved in leaf-related traits, three connected recombination inbred line (RIL) populations including 538 RILs were genotyped by genotyping-by-sequencing (GBS) method and phenotyped for the leaf angle and related traits in six environments. We conducted single population quantitative trait locus (QTL) mapping and joint linkage analysis based on high-density recombination bin maps constructed from GBS genotype data. A total of 45 QTLs with phenotypic effects ranging from 1.2% to 29.2% were detected for four leaf architecture traits by using joint linkage mapping across the three populations. All the QTLs identified for each trait could explain approximately 60% of the phenotypic variance. Four QTLs were located on small genomic regions where candidate genes were found. Genomic predictions from a genomic best linear unbiased prediction (GBLUP) model explained 45±9% to 68±8% of the variation in the remaining RILs for the four traits. These results extend our understanding of the genetics of leaf traits and can be used in genomic prediction to accelerate plant architecture improvement.

Joint‐linkage mapping and GWAS reveal extensive genetic loci that regulate male inflorescence size in maize

Summary Both insufficient and excessive male inflorescence size leads to a reduction in maize yield. Knowledge of the genetic architecture of male inflorescence is essential to achieve the optimum inflorescence size for maize breeding. In this study, we used approximately eight thousand inbreds, including both linkage populations and association populations, to dissect the genetic architecture of male inflorescence. The linkage populations include 25 families developed in the U.S. and 11 families developed in China. Each family contains approximately 200 recombinant inbred lines (RILs). The association populations include approximately 1000 diverse lines from the U.S. and China. All inbreds were genotyped by either sequencing or microarray. Inflorescence size was measured as the tassel primary branch number (TBN) and tassel length (TL). A total of 125 quantitative trait loci (QTLs) were identified (63 for TBN, 62 for TL) through linkage analyses. In addition, 965 quantitative trait nucleotides (QTNs) were identified through genomewide study (GWAS) at a bootstrap posterior probability (BPP) above a 5% threshold. These QTLs/QTNs include 24 known genes that were cloned using mutants, for example Ramosa3 (ra3), Thick tassel dwarf1 (td1), tasselseed2 (ts2), liguleless2 (lg2), ramosa1 (ra1), barren stalk1 (ba1), branch silkless1 (bd1) and tasselseed6 (ts6). The newly identified genes encode a zinc transporter (e.g. GRMZM5G838098 and GRMZM2G047762), the adapt in terminal region protein (e.g. GRMZM5G885628), O‐methyl‐transferase (e.g. GRMZM2G147491), helix‐loop‐helix (HLH) DNA‐binding proteins (e.g. GRMZM2G414252 and GRMZM2G042895) and an SBP‐box protein (e.g. GRMZM2G058588). These results provide extensive genetic information to dissect the genetic architecture of inflorescence size for the improvement of maize yield.

Construction of high-quality recombination maps with low-coverage genomic sequencing for joint linkage analysis in maize

BackgroundA genome-wide association study (GWAS) is the foremost strategy used for finding genes that control human diseases and agriculturally important traits, but it often reports false positives. In contrast, its complementary method, linkage analysis, provides direct genetic confirmation, but with limited resolution. A joint approach, using multiple linkage populations, dramatically improves resolution and statistical power. For example, this approach has been used to confirm that many complex traits, such as flowering time controlling adaptation in maize, are controlled by multiple genes with small effects. In addition, genotyping by sequencing (GBS) at low coverage not only produces genotyping errors, but also results in large datasets, making the use of high-throughput sequencing technologies computationally inefficient or unfeasible.ResultsIn this study, we converted raw SNPs into effective recombination bins. The reduced bins not only retain the original information, but also correct sequencing errors from low-coverage genomic sequencing. To further increase the statistical power and resolution, we merged a new temperate maize nested association mapping (NAM) population derived in China (CN-NAM) with the existing maize NAM population developed in the US (US-NAM). Together, the two populations contain 36 families and 7,000 recombinant inbred lines (RILs). One million SNPs were generated for all the RILs with GBS at low coverage. We developed high-quality recombination maps for each NAM population to correct genotyping errors and improve the computational efficiency of the joint linkage analysis. The original one million SNPs were reduced to 4,932 and 5,296 recombination bins with average interval distances of 0.34xa0cM and 0.28xa0cM for CN-NAM and US-NAM, respectively. The quantitative trait locus (QTL) mapping for flowering time (days to tasseling) indicated that the high-density, recombination bin map improved resolution of QTL mapping by 50xa0% compared with that using a medium-density map. We also demonstrated that combining the CN-NAM and US-NAM populations improves the power to detect QTL by 50xa0% compared to single NAM population mapping. Among the QTLs mapped by joint usage of the US-NAM and CN-NAM maps, 25xa0% of the QTLs overlapped with known flowering-time genes in maize.ConclusionThis study provides directions and resources for the research community, especially maize researchers, for future studies using the recombination bin strategy for joint linkage analysis. Available resources include efficient usage of low-coverage genomic sequencing, detailed positions for genes controlling maize flowering, and recombination bin maps and flowering- time data for both CN and US NAMs. Maize researchers even have the opportunity to grow both CN and US NAM populations to study the traits of their interest, as the seeds of both NAM populations are available from the seed repository in China and the US.

Identification of genetic variants associated with maize flowering time using an extremely large multi-genetic background population.

Flowering time is one of the major adaptive traits in domestication of maize and an important selection criterion in breeding. To detect more maize flowering time variants we evaluated flowering time traits using an extremely large multi- genetic background population that contained more than 8000 lines under multiple Sino-United States environments. The population included two nested association mapping (NAM) panels and a natural association panel. Nearly 1 million single-nucleotide polymorphisms (SNPs) were used in the analyses. Through the parallel linkage analysis of the two NAM panels, both common and unique flowering time regions were detected. Genome wide, a total of 90 flowering time regions were identified. One-third of these regions were connected to traits associated with the environmental sensitivity of maize flowering time. The genome-wide association study of the three panels identified nearly 1000 flowering time-associated SNPs, mainly distributed around 220 candidate genes (within a distance of 1xa0Mb). Interestingly, two types of regions were significantly enriched for these associated SNPs - one was the candidate gene regions and the other was the approximately 5xa0kb regions away from the candidate genes. Moreover, the associated SNPs exhibited high accuracy for predicting flowering time.

Analysis of recombination QTLs, segregation distortion, and epistasis for fitness in maize multiple populations using ultra-high-density markers

Key messageUsing two nested association mapping populations and high-density markers, some important genomic regions controlling recombination frequency and segregation distortion were detected.AbstractUnderstanding the maize genomic features would be useful for the study of genetic diversity and evolution and for maize breeding. Here, we used two maize nested association mapping (NAM) populations separately derived in China (CN-NAM) and the US (US-NAM) to explore the maize genomic features. The two populations containing 36 families and about 7000 recombinant inbred lines were evaluated with genotyping-by-sequencing. Through the comparison between the two NAMs, we revealed that segregation distortion is little, whereas epistasis for fitness is present in the two maize NAM populations. When conducting quantitative trait loci (QTL) mapping for the total number of recombination events, we detected 14 QTLs controlling recombination. Using high-density markers to identify segregation distortion regions (SDRs), a total of 445 SDRs were detected within the 36 families, among which 15 common SDRs were found in at least ten families. About 80xa0% of the known maize gametophytic factors (ga) genes controlling segregation distortion were overlapped with highly significant SDRs. In addition, we also found that the regions with high recombination rate and high gene density usually tended to have little segregation distortion. This study will facilitate population genetic studies and gene cloning affecting recombination variation and segregation distortion in maize, which can improve plant breeding progress.

Transcriptome Sequencing Identified Genes and Gene Ontologies Associated with Early Freezing Tolerance in Maize

Originating in a tropical climate, maize has faced great challenges as cultivation has expanded to the majority of the worlds temperate zones. In these zones, frost and cold temperatures are major factors that prevent maize from reaching its full yield potential. Among 30 elite maize inbred lines adapted to northern China, we identified two lines of extreme, but opposite, freezing tolerance levels—highly tolerant and highly sensitive. During the seedling stage of these two lines, we used RNA-seq to measure changes in maize whole genome transcriptome before and after freezing treatment. In total, 19,794 genes were expressed, of which 4550 exhibited differential expression due to either treatment (before or after freezing) or line type (tolerant or sensitive). Of the 4550 differently expressed genes, 948 exhibited differential expression due to treatment within line or lines under freezing condition. Analysis of gene ontology found that these 948 genes were significantly enriched for binding functions (DNA binding, ATP binding, and metal ion binding), protein kinase activity, and peptidase activity. Based on their enrichment, literature support, and significant levels of differential expression, 30 of these 948 genes were selected for quantitative real-time PCR (qRT-PCR) validation. The validation confirmed our RNA-Seq-based findings, with squared correlation coefficients of 80% and 50% in the tolerance and sensitive lines, respectively. This study provided valuable resources for further studies to enhance understanding of the molecular mechanisms underlying maize early freezing response and enable targeted breeding strategies for developing varieties with superior frost resistance to achieve yield potential.

Numerous genetic loci identified for drought tolerance in the maize nested association mapping populations

BackgroundMaize requires more water than most other crops; therefore, the water use efficiency of this crop must be improved for maize production under undesirable land and changing environmental conditions.ResultsTo elucidate the genetic control of drought in maize, we evaluated approximately 5000 inbred lines from 30 linkage-association joint mapping populations under two contrasting water regimes for seven drought-related traits, including yield and anthesis-silking interval (ASI). The joint linkage analysis was conducted to identify 220 quantitative trait loci (QTLs) under well-watered conditions and 169 QTLs under water-stressed conditions. The genome-wide association analysis identified 365 single nucleotide polymorphisms (SNPs) associated with drought-related traits, and these SNPs were located in 354 candidate genes. Fifty-two of these genes showed significant differential expression in the inbred line B73 under the well-watered and water-stressed conditions. In addition, genomic predictions suggested that the moderate-density SNPs obtained through genotyping-by-sequencing were able to make accurate predictions in the nested association mapping population for drought-related traits with moderate-to-high heritability under the water-stressed conditions.ConclusionsThe results of the present study provide important information that can be used to understand the genetic basis of drought stress responses and facilitate the use of beneficial alleles for the improvement of drought tolerance in maize.

Fine-mapping of qGW4.05, a major QTL for kernel weight and size in maize

BackgroundKernel weight and size are important components of grain yield in cereals. Although some information is available concerning the map positions of quantitative trait loci (QTL) for kernel weight and size in maize, little is known about the molecular mechanisms of these QTLs. qGW4.05 is a major QTL that is associated with kernel weight and size in maize. We combined linkage analysis and association mapping to fine-map and identify candidate gene(s) at qGW4.05.ResultsQTL qGW4.05 was fine-mapped to a 279.6-kb interval in a segregating population derived from a cross of Huangzaosi with LV28. By combining the results of regional association mapping and linkage analysis, we identified GRMZM2G039934 as a candidate gene responsible for qGW4.05. Candidate gene-based association mapping was conducted using a panel of 184 inbred lines with variable kernel weights and kernel sizes. Six polymorphic sites in the gene GRMZM2G039934 were significantly associated with kernel weight and kernel size.ConclusionThe results of linkage analysis and association mapping revealed that GRMZM2G039934 is the most likely candidate gene for qGW4.05. These results will improve our understanding of the genetic architecture and molecular mechanisms underlying kernel development in maize.

Analysis of genetic differentiation and genomic variation to reveal potential regions of importance during maize improvement

BackgroundExploring genetic differentiation and genomic variation is important for both the utilization of heterosis and the dissection of the genetic bases of complex traits.MethodsWe integrated 1857 diverse maize accessions from America, Africa, Europe and Asia to investigatetheir genetic differentiation, genomic variation using 43,252 high-quality single-nucleotide polymorphisms(SNPs),combing GWAS and linkage analysis strategy to exploring the function of relevant genetic segments.ResultsWe uncovered many more subpopulations that recently or historically formed during the breeding process. These patterns are represented by the following lines: Mo17, GB, E28, Ye8112, HZS, Shen137, PHG39, B73, 207, A634, Oh43, Reid Yellow Dent, and the Tropical/subtropical (TS) germplasm. A total of 85 highly differentiated regions with a DEST of more than 0.2 were identified between the TS and temperate subpopulations. These regions comprised 79xa0% of the genetic variation, and most were significantly associated with adaptive traits. For example, the region containing the SNP tag PZE.108075114 was highly differentiated, and this region was significantly associated with flowering time (FT)-related traits, as supported by a genome-wide association study (GWAS) within the interval of FT-related quantitative trait loci (QTL). This region was also closely linked to zcn8 and vgt1, which were shown to be involved in maize adaptation. Most importantly, 197 highly differentiated regions between different subpopulation pairs were located within an FT- or plant architecture-related QTL.ConclusionsHere we reported that 700–1000 SNPs were necessary needed to robustly estimate the genetic differentiation of a naturally diverse panel. In addition, 13 subpopulations were observed in maize germplasm, 85 genetic regions with higher differentiation between TS and temperate maize germplasm, 197 highly differentiated regions between different subpopulation pairs, which contained some FT- related QTNs/QTLs/genes supported by GWAS and linkage analysis, and these regions were expected to play important roles in maize adaptation.