In particular, RVB states are considered to be an important system to study the ground state properties of the doped quantum spin1/2 ladder. It is therefore interesting to understand how quantum correlations are distributed among the constituents of these composite systems. In this regard, we formulate an analytical recursive method to generate the wave function of doped shortrange resonating valence bond (RVB) states as a tool to efficiently estimate multisite entanglement as well as other physical quantities in doped quantum spin ladders. Importantly, our results show that within a specific doping concentration and model parameter regimes, the doped RVB state essentially characterizes the trends of genuine multiparty entanglement in the exact ground states of a Hubbard model with large onsite interactions. Moreover, we consider an isotropic RVB network of spin1/2 particles with a finite fraction of defects, where the corresponding wave function of the network is rotationally invariant under the action of local unitaries. By using quantum informationtheoretic concepts like strong subadditivity of von Neumann entropy and approximate quantum telecloning, we prove analytically that in the presence of defects, caused by loss of a finite fraction of spins, the RVB network sustains genuine multisite entanglement, and at the same time may exhibit finite moderaterange bipartite entanglement, in contrast to the case with no defects.