We consider model order reduction of bilinear descriptor systems using an interpolatory projection framework. Such nonlinear descriptor systems can be represented by a series of generalized linear descriptor systems (also called subsystems) by utilizing the Volterra-Wiener approach [22]. Standard projection techniques for bilinear systems utilize the generalized transfer function of these subsystems to construct an interpolating approximation. However, the resulting reduced-order system may not match the polynomial part of the generalized transfer functions. This may result in an unbounded error in terms of ℋ₂ or ℋ∞ norms. In this paper, we derive an explicit expression for the polynomial part of each subsystem by assuming a special structure of the bilinear system which reduces to an index-1 linear DAE if the bilinear term is zero. This allows us to propose an interpolatory technique for bilinear DAEs which not only achieves interpolation but also retains the polynomial part of the bilinear system. The approach extends the interpolatory technique for index-1 linear DAEs [18] to bilinear DAEs. Numerical examples are used to illustrate the theoretical results.