In this thesis we present our contributions to symbolic summation, extending WZFasenmyer methods to handle definite hypergeometric sums with nonstandard boundary conditions and to compute recurrences for multiple Mellin-Barnes integrals over hypergeometric terms. We also include concrete applications of these methods to Feynman integral calculus, as well as for proving identities involving definite integrals and special functions.First we give a short introduction to WZ-summation methods, including K. Wegschaider's approach to this method and its implementation in the package MultiSum. Inspired by work of Sister Celine Fasenmyer, these techniques were introduced by H. Wilf and D. Zeilberger to algorithmically compute recurrences for multiple sums over hypergeometric terms. Their procedure is based on finding a certificate recurrence satisfied by the hypergeometric summand and summing over this difference equation to obtain a recurrence for the nested sum. Our proofs of two nontrivial special function identities involving Gegenbauer polynomials, provide classic applications of the method.As part of the collaboration between RISC and DESY coordinated by C. Schneider, we developed an algorithmic approach to compute Feynman parameter integrals after rewriting them as multisums over hypergeometric terms to fit the input class of classic summation algorithms.Since these definite sums have nonstandard boundary conditions, the WZ-method delivers inhomogeneous recurrence relations. We designed a recursive procedure to determine the inhomogeneous parts of these recurrences and implemented it in the Mathematica package FSums which builds on the already existing packages, MultiSum and C. Schneider's Sigma.Another approach to evaluate Feynman integrals is by representating them in terms of nested Mellin-Barnes integrals. These complex contour integrals can also be viewed as sums of residues at certain poles of the integrands and they are connected to the inversion formula for the Mellin transform.In the last part, we show how WZ-methods can be used to compute recurrences for multiple Mellin-Barnes integrals over hypergeometric terms, eliminating the need to search for sum representations. We applied this new algorithmic technique to prove typical entries from the Gradshteyn-Ryzhik