Second order β-decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double-β decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double-β decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double-β decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double-β decay. LQCD calculations of the hadronic contributions to the neutrinoful process nn → ppe − e −ν eνe and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double-β decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double-β decay processes themselves or using inputs from LQCD. Future prospects for direct calculations in larger nuclei are also discussed. many-body methods and effective field theory (EFT) analysis of the transition operators, where lattice QCD can provide key input. To improve upon this situation, recent efforts have advocated an "endto-end" EFT analysis of 0νββ to link the scale Λ of LNV to nuclear scales. This multi-prong approach includes various steps:1. The use of the Standard Model EFT to link the scale Λ of LNV to the hadronic scale Λ χ ∼ O(1) GeV, where non-perturbative QCD effects arise. This step is by now mature: the operator basis (to which any underlying model can be matched) is known up to dimension-nine and the renormalization group evolution of these operators under strong interactions is known. Light new degrees of freedom (such as sterile neutrinos) can also be included in this framework.2. The matching of the quark-gluon level EFT to hadronic EFTs such as Chiral Perturbation Theory (χPT) in the meson and single nucleon sector, and chiral EFT and pionless EFT in the multinucleon sector. This step can be performed consistently in the strong and weak sectors of the theory, which in the case of interest here involves ∆L = 2 transition operators. The form of the transition operators is known to...