Regarding how to fix the issue, it depends on the level of precission you want. From the nVenn algorithm, I authored the nVennR package to create quasi-proportional Euler diagrams. With the caveats mentioned in the link, you can represent larger numbers of sets and show the relative size of each region. In your example,
library(nVennR)
myV <- createVennObj(nSets = 4, sNames = c('A', 'B', 'C', 'D'), sSizes = c(0, 26, 53, 7, 22, 5, 16, 3, 54, 10, 29, 4, 20, 5, 14, 3))
myV <- plotVenn(nVennObj = myV)
And the result would be:

Depending on your requirements, this may not be satisfactory. The proportionality is in the area of the circles, not the regions (you can see that the region 1, 2, 3, 4 - A&B&C&D - has empty space. However, this strategy overcomes the limitations of regular shapes in these representations mentioned by Johan Larsson. If you are interested, there are more details in the vignette.
The reason why some areas are left out is simple: the diagram is inexact and is missing some areas. There is no place to put the label for B&C so that's why B and C are missing 2 units. There likely isn't any way (or at least eulerr cannot find it) to perfectly represent your combination with an Euler diagram using ellipses. You either have to accept that it is inexact or try another solution.
Similarly, the residual plot cannot show the missing residuals graphically because there is no area representing them. I am, by the way, the author of this package and I do have something better in mind for the residual plot which would display missing areas as well, but I haven't had time to implement it yet.

