Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively.
To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements.
Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads.

Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6
This is a list of supplementary files associated with this preprint. Click to download.
Loading...
On 18 Jan, 2021
Received 12 Jan, 2021
On 04 Jan, 2021
Invitations sent on 04 Jan, 2021
On 04 Jan, 2021
On 04 Jan, 2021
On 04 Jan, 2021
Received 31 Dec, 2020
Received 27 Dec, 2020
On 23 Dec, 2020
On 21 Dec, 2020
On 21 Dec, 2020
Invitations sent on 21 Dec, 2020
On 21 Dec, 2020
On 21 Dec, 2020
Posted 10 Jun, 2020
On 20 Nov, 2020
Received 17 Nov, 2020
On 04 Nov, 2020
Received 23 Aug, 2020
On 01 Jul, 2020
Invitations sent on 29 Jun, 2020
On 29 May, 2020
On 28 May, 2020
On 28 May, 2020
On 27 May, 2020
On 18 Jan, 2021
Received 12 Jan, 2021
On 04 Jan, 2021
Invitations sent on 04 Jan, 2021
On 04 Jan, 2021
On 04 Jan, 2021
On 04 Jan, 2021
Received 31 Dec, 2020
Received 27 Dec, 2020
On 23 Dec, 2020
On 21 Dec, 2020
On 21 Dec, 2020
Invitations sent on 21 Dec, 2020
On 21 Dec, 2020
On 21 Dec, 2020
Posted 10 Jun, 2020
On 20 Nov, 2020
Received 17 Nov, 2020
On 04 Nov, 2020
Received 23 Aug, 2020
On 01 Jul, 2020
Invitations sent on 29 Jun, 2020
On 29 May, 2020
On 28 May, 2020
On 28 May, 2020
On 27 May, 2020
Availability of plant genome sequences has led to significant advances. However, with few exceptions, the great majority of existing genome assemblies are derived from short read sequencing technologies with highly uneven read coverages indicative of sequencing and assembly issues that could significantly impact any downstream analysis of plant genomes. In tomato for example, 0.6% (5.1 Mb) and 9.7% (79.6 Mb) of short-read based assembly had significantly higher and lower coverage compared to background, respectively.
To understand what the causes may be for such uneven coverage, we first established machine learning models capable of predicting genomic regions with variable coverages and found that high coverage regions tend to have higher simple sequence repeat and tandem gene densities compared to background regions. To determine if the high coverage regions were misassembled, we examined a recently available tomato long-read based assembly and found that 27.8% (1.41 Mb) of high coverage regions were potentially misassembled of duplicate sequences, compared to 1.4% in background regions. In addition, using a predictive model that can distinguish correctly and incorrectly assembled high coverage regions, we found that misassembled, high coverage regions tend to be flanked by simple sequence repeats, pseudogenes, and transposon elements.
Our study provides insights on the causes of variable coverage regions and a quantitative assessment of factors contributing to plant genome misassembly when using short reads.

Figure 1

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6
This is a list of supplementary files associated with this preprint. Click to download.
Loading...