Isogeometric Analysis (IGA) is a computational methodology for the numerical approximation of Partial Differential Equations (PDEs). IGA is based on the isogeometric concept, for which the same basis functions, usually Non-Uniform Rational B-Splines (NURBS), are used both to represent the geometry and to approximate the unknown solutions of PDEs. Compared to the standard Finite Element method, NURBS-based IGA offers several advantages: ideally a direct interface with CAD tools, exact geometrical representation, simple refinement procedures, and smooth basis functions allowing to easily solve higher order problems, including structural shell problems. In these contexts, repeatedly solving a problem for a large set of geometric parameters might lead to high and eventually prohibitive computational cost. To cope with this problem, we consider in this work the Reduced Basis (RB) method for the solution of parameter dependent PDEs, specifically for which the NURBS representation of the computational domain is parameter dependent. RB refers to a technique that enables a rapid and reliable approximation of parametrized PDEs by constructing low dimensional approximation spaces. In this work, for the construction of the reduced spaces we adopt two different strategies, namely the Proper Orthogonal Decomposition and the greedy algorithm. In this thesis we combine RB and IGA for the efficient solution of parametrized problems for all the possible cases of NURBS geometrical parametrizations, which specifically include the NURBS control points, the weights, and both the control points and weights. In particular, we first focus on the solution of second order PDEs on parametrized lower dimensional manifolds, specifically surfaces in the three dimensional space. We consider geometrical parametrizations that entail a nonaffine dependence of the variational forms on the spatial coordinates and the geometric parameters. Thus, depending on the parametrization at hand and in order to ensure a suitable Offline/Online decomposition between the reduced order model construction and solution, we resort to the Empirical Interpolation Method (EIM) or the Matrix Discrete Empirical Interpolation Method (MDEIM), by comparing their performances. As application, we solve a class of benchmark structural problems modeled by Kirchoff-Love shells for which we consider NURBS geometric parametrizations and we apply the RB method to the solution of this class of fourth order PDEs. We highlight by means of numerical tests, the performance of the RB method applied to standard IGA approximation of parametrized shell geometries.