Student or educator? Sign up with your school email and get 6 months of Pro free.
Cited by 1,210

The DOI matches its Crossref record; 51 of the 51 references it cites were found in a registry or an index.

Improved maize reference genome with single-molecule technologies

Jiao, Yinping, Peluso, Paul, Shi, Jinghua, Liang, Tiffany, Stitzer, Michelle C., Wang, Bo, Campbell, Michael S., Stein, Joshua C., Wei, Xuehong, Chin, Chen-Shan, Guill, Katherine, Regulski, Michael, Kumari, Sunita, Olson, Andrew, Gent, Jonathan, Schneider, Kevin L., Wolfgruber, Thomas K., May, Michael R., Springer, Nathan M., Antoniou, Eric, McCombie, W. Richard, Presting, Gernot G., McMullen, Michael, Ross-Ibarra, Jeffrey, Dawe, R. Kelly, Hastie, Alex, Rank, David R., Ware, Doreen·2017·Nature

Published version of a preprint

Abstract

ABSTRACT Complete and accurate reference genomes and annotations provide fundamental tools for characterization of genetic and functional variation. These resources facilitate elucidation of biological processes and support translation of research findings into improved and sustainable agricultural technologies. Many reference genomes for crop plants have been generated over the past decade, but these genomes are often fragmented and missing complex repeat regions. Here, we report the assembly and annotation of maize, a genetic and agricultural model species, using Single Molecule Real-Time (SMRT) sequencing and high-resolution optical mapping. Relative to the previous reference genome, our assembly features a 52-fold increase in contig length and significant improvements in the assembly of intergenic spaces and centromeres. Characterization of the repetitive portion of the genome revealed over 130,000 intact transposable elements (TEs), allowing us to identify TE lineage expansions unique to maize. Gene annotations were updated using 111,000 full-length transcripts obtained by SMRT sequencing. In addition, comparative optical mapping of two other inbreds revealed a prevalence of deletions in the low gene density region and maize lineage-specific genes.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/nature22971

Cites

51 of 51 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Genome assembly layout. a, Workflow for genome construction. b, Ideograms of maize B73 version 4 reference pseudomolecules. The top track shows positions of 2,522 gaps in the pseudomolecules, including 1,115 gaps in which the lengths were estimated using optical genome maps (orange), whereas the remainder (purple) have undetermined lengths. More than half of the assembly is constituted of contigs longer than 1 Mb, which are shown as light grey bars in the bottom track. PowerPoint slide Source data
Figure 1 Genome assembly layout. a, Workflow for genome construction. b, Ideograms of maize B73 version 4 reference pseudomolecules. The top track shows positions of 2,522 gaps in the pseudomolecules, including 1,115 gaps in which the lengths were estimated using optical genome maps (orange), whereas the remainder (purple) have undetermined lengths. More than half of the assembly is constituted of contigs longer than 1 Mb, which are shown as light grey bars in the bottom track. PowerPoint slide Source dataCC BY · PubMed Central
Figure 2 Phylogeny of maize and sorghum LTR retrotransposon families. a, b, Both Ty3/Gypsy (a) and Ty1/Copia (b) superfamilies are present at higher copy number in maize (red) than in sorghum (blue). Bars (log10-scaled) depict family copy numbers. PowerPoint slide Source data
Figure 2 Phylogeny of maize and sorghum LTR retrotransposon families. a, b, Both Ty3/Gypsy (a) and Ty1/Copia (b) superfamilies are present at higher copy number in maize (red) than in sorghum (blue). Bars (log10-scaled) depict family copy numbers. PowerPoint slide Source dataCC BY · PubMed Central
Figure 3 Structural variation from Ki11 and W22. a, Alignment and structural variation called from Ki11 and W22 optical maps on chromosome 10. b, Size distribution of the insertion and deletions in Ki11 and W22. c, Example of using short-read alignment to verify a missing region mapped in Ki11. PowerPoint slide Source data
Figure 3 Structural variation from Ki11 and W22. a, Alignment and structural variation called from Ki11 and W22 optical maps on chromosome 10. b, Size distribution of the insertion and deletions in Ki11 and W22. c, Example of using short-read alignment to verify a missing region mapped in Ki11. PowerPoint slide Source dataCC BY · PubMed Central
Extended Data Figure 1 Summary of data generated for genome construction. a, Size distribution of single molecules for the optical maps. A total of 150 Gb (~60-fold coverage) of single-molecule raw data from BioNano chips was collected for map construction. The N50 of the single molecules was ~261 kb, and the label density was 11.6 per 100 kb. After assembly, the total size of the map reached 2.12 Gb with an N50 of 2.47 Mb. b, Length distribution of SMRT sequencing reads. Sequencing of 212 P6-C4 SMRT cells on the PacBio platform generated ~65-fold depth-of-coverage of the nuclear genome. Read lengths averaged 11.7 kb, with reads above 10 kb providing 53-fold depth-of-coverage. c, The accuracy of SMRT sequencing from a representative run. The sequencing error rate was estimated at 10% from the alignment with the maize B73 RefGen_v3 by BLASR. d, Plot of the fraction of alignable data per run (alignable bases/total bases per chip) versus total raw bases (per chip) for each B73 sequencing run. As the plot shows, the trend in the data suggests that as the overall per run raw base yield increases, the fraction of alignable bases decreases. This is owing to the fact that in all runs, a subset of the zero-mode waveguide (ZMWs) will initially have more than one active sequencing enzyme in the observation field at the start of the sequencing run. A ZMW with more than one active polymerase will create unalignable bases while the two polymerases are simultaneously synthesizing DNA and yield a ‘merged sequencing signal from two independent polymerases’. As the loading of a chips increases (yield of bases), the probability of having two or more polymerases in a single ZMW increases.
Extended Data Figure 1 Summary of data generated for genome construction. a, Size distribution of single molecules for the optical maps. A total of 150 Gb (~60-fold coverage) of single-molecule raw data from BioNano chips was collected for map construction. The N50 of the single molecules was ~261 kb, and the label density was 11.6 per 100 kb. After assembly, the total size of the map reached 2.12 Gb with an N50 of 2.47 Mb. b, Length distribution of SMRT sequencing reads. Sequencing of 212 P6-C4 SMRT cells on the PacBio platform generated ~65-fold depth-of-coverage of the nuclear genome. Read lengths averaged 11.7 kb, with reads above 10 kb providing 53-fold depth-of-coverage. c, The accuracy of SMRT sequencing from a representative run. The sequencing error rate was estimated at 10% from the alignment with the maize B73 RefGen_v3 by BLASR. d, Plot of the fraction of alignable data per run (alignable bases/total bases per chip) versus total raw bases (per chip) for each B73 sequencing run. As the plot shows, the trend in the data suggests that as the overall per run raw base yield increases, the fraction of alignable bases decreases. This is owing to the fact that in all runs, a subset of the zero-mode waveguide (ZMWs) will initially have more than one active sequencing enzyme in the observation field at the start of the sequencing run. A ZMW with more than one active polymerase will create unalignable bases while the two polymerases are simultaneously synthesizing DNA and yield a ‘merged sequencing signal from two independent polymerases’. As the loading of a chips increases (yield of bases), the probability of having two or more polymerases in a single ZMW increases.CC BY · PubMed Central
Extended Data Figure 2 Construction of pseudomolecules. a, Summary of the three assembly sets. b, How the scaffolds were ordered according to the order of the BACs. c, Size distribution of gaps in the pseudomolecules estimated using the optical map.
Extended Data Figure 2 Construction of pseudomolecules. a, Summary of the three assembly sets. b, How the scaffolds were ordered according to the order of the BACs. c, Size distribution of gaps in the pseudomolecules estimated using the optical map.CC BY · PubMed Central
Extended Data Figure 3 Quality assessment and comparison of the assembly in centromere and telomere regions in maize B73 RefGen_v3 and v4. a, Quality assessment of centromere and telomere using optical genome map. b, Locations of centromeres on pseudomolecules defined by ChIP–seq in the B73 RefGen_v3 and v4. c, Telomere repeats found in the B73 RefGen_v4 pseudomolecules.
Extended Data Figure 3 Quality assessment and comparison of the assembly in centromere and telomere regions in maize B73 RefGen_v3 and v4. a, Quality assessment of centromere and telomere using optical genome map. b, Locations of centromeres on pseudomolecules defined by ChIP–seq in the B73 RefGen_v3 and v4. c, Telomere repeats found in the B73 RefGen_v4 pseudomolecules.CC BY · PubMed Central
Extended Data Figure 4 Details of the gene annotation of maize B73 RefGen_v4. a, The pipeline used to characterize high confidence gene models. b, Summary of B73 RefGen_v4 protein-coding gene annotation, and comparison with RefGen_v3 annotation.
Extended Data Figure 4 Details of the gene annotation of maize B73 RefGen_v4. a, The pipeline used to characterize high confidence gene models. b, Summary of B73 RefGen_v4 protein-coding gene annotation, and comparison with RefGen_v3 annotation.CC BY · PubMed Central
Extended Data Figure 5 Improvement of the annotation of alternative splicing and completeness of regulatory regions of maize RefGen_v4 genes. a, Number of transcripts of each gene in v3 and v4 annotation. b, Percentages of genes with gaps in flanking regions in the v3 and v4 annotations.
Extended Data Figure 5 Improvement of the annotation of alternative splicing and completeness of regulatory regions of maize RefGen_v4 genes. a, Number of transcripts of each gene in v3 and v4 annotation. b, Percentages of genes with gaps in flanking regions in the v3 and v4 annotations.CC BY · PubMed Central
Extended Data Figure 6 Comparative analysis of the maize B73 RefGen_v4 genes with other grasses. a, Species-membership in orthologue sets, giving counts and percentage of orthologue sets of which each species is a member. Numbers in parentheses give the percentage of orthologue sets with membership of all species and versions within the clade. na, not applicable. b, Venn diagram showing overlap of 6,539 orthologue sets rooted in the Poaceae (true grasses) that are deficient in gene membership among five species.
Extended Data Figure 6 Comparative analysis of the maize B73 RefGen_v4 genes with other grasses. a, Species-membership in orthologue sets, giving counts and percentage of orthologue sets of which each species is a member. Numbers in parentheses give the percentage of orthologue sets with membership of all species and versions within the clade. na, not applicable. b, Venn diagram showing overlap of 6,539 orthologue sets rooted in the Poaceae (true grasses) that are deficient in gene membership among five species.CC BY · PubMed Central
Extended Data Figure 7 Structural variation characterized from the Ki11 and W22 optical maps.
Extended Data Figure 7 Structural variation characterized from the Ki11 and W22 optical maps.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Cited by

See all 1,210

Metadata

Authors
Jiao, Y., Peluso, P., Shi, J., Liang, T., Stitzer, M. C., Wang, B., Campbell, M. S., Stein, J. C., Wei, X., Chin, C.-S., Guill, K., Regulski, M., Kumari, S., Olson, A., Gent, J., Schneider, K. L., Wolfgruber, T. K., May, M. R., Springer, N. M., Antoniou, E., McCombie, W. R., Presting, G. G., McMullen, M., Ross-Ibarra, J., Dawe, R. K., Hastie, A., Rank, D. R., & Ware, D.
Venue
Nature
Published
2017-06
Volume
546
Issue
7659
Pages
524-527
ISSN
  • 0028-0836
  • 1476-4687
Cited by
1,210
References
51
Language
en
Licence
creativecommons.org/licenses/by/4.0 (tdm)

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

Data sources

Bibliographic data: Semantic Scholar, OpenAlex, Crossref.

Contains information from the Semantic Scholar Academic Graph Datasets, made available under the ODC Attribution License (ODC-BY 1.0).

Semantic Scholar’s data is described in: Kinney, R., et al. (2023). The Semantic Scholar Open Data Platform. arXiv:2301.10140.