We present a novel method for the mechanical simulation of slender, elastic, spatial rods and rod structures subject to large deformation and rotation. We develop an isogeometric collocation method for the geometrically exact, nonlinear Cosserat rod theory. The rod centerlines are represented as spatial NURBS curves and cross-section orientations are parameterized in terms of unit quaternions as 4-dimensional NURBS curves. Within the isogeometric framework, the strong forms of the equilibrium equations of forces and moments of the discretized Cosserat model are collocated, leading to an efficient method for higher-order discretizations. For rod structures consisting of multiple, connected rods we introduce a formulation with rigid, quasi-G1-coupling. It is based on the strong enforcement of continuity of displacement and change of cross-section orientation at interfaces. We also develop a mixed isogeometric formulation, which is based on an independent discretization of internal forces and moments and alleviates shear locking for thin rods. The novel rod simulation methods are verified by numerical convergence studies. Further computational examples include realistic applications with large deformations and rotations, as well as a large-scale rod structure with several hundreds of coupled rods and complex buckling behavior.