One of the most pressing questions in modified gravity is how deviations from general relativity can manifest in upcoming galaxy surveys. This is especially relevant for theories exhibiting Vainshtein screening, where such deviations are efficiently suppressed within a (typically large) Vainshtein radius. However, Vainshtein screening is known to be shape dependent: it is most effective around spherical sources, weaker around cylindrical objects and completely absent for planar sources. The Cosmic Web therefore offers a testing ground, as it displays many shapes in the form of clusters, filaments and walls. In this work, we explicitly derive the signature of the shape dependence of Vainshtein screening on the matter bispectrum, by considering a cubic Galileon model with a conformal coupling to matter and a cosmological constant. We perform a second order perturbative analysis, deriving analytic, integral expressions for the bispectrum, which we integrate using hi_class. We find that the shape dependence of Vainshtein screening enters the bispectrum with a unique scale-factor dependence of $\propto a^{3/2}$. The magnitude of the effect today is up to 2 % for a model whose linear growth rate deviates up to 5 % from $\Lambda$CDM.

We present a new framework for testing the isotropy of the Universe using cosmic microwave background data, building on the nested-sampling ANICOSMO code. Uniquely, we are able to constrain the scalar, vector and tensor degrees of freedom alike; previous studies only considered the vector mode (linked to vorticity). We employ Bianchi type VII$_h$ cosmologies to model the anisotropic Universe, from which other types may be obtained by taking suitable limits. In a separate development, we improve the statistical analysis by including the effect of Bianchi power in the high-$\ell$, as well as the low-$\ell$, likelihood. To understand the effect of all these changes, we apply our new techniques to WMAP data. We find no evidence for anisotropy, constraining shear in the vector mode to $(\sigma_V/H)_0 < 1.7 \times 10^{-10}$ (95% CL). For the first time, we place limits on the tensor mode; unlike other modes, the tensor shear can grow from a near-isotropic early Universe. The limit on this type of shear is $(\sigma_{T,\rm reg}/H)_0 < 2.4 \times 10^{-7}$ (95% CL).

A fundamental assumption in the standard model of cosmology is that the Universe is isotropic on large scales. Breaking this assumption leads to a set of solutions to Einstein's field equations, known as Bianchi cosmologies, only a subset of which have ever been tested against data. For the first time, we consider all degrees of freedom in these solutions to conduct a general test of isotropy using cosmic microwave background temperature and polarization data from Planck. For the vector mode (associated with vorticity), we obtain a limit on the anisotropic expansion of $(\sigma_V/H)_0 < 4.7 \times 10^{-11}$ (95% CI), which is an order of magnitude tighter than previous Planck results that used CMB temperature only. We also place upper limits on other modes of anisotropic expansion, with the weakest limit arising from the regular tensor mode, $(\sigma_{T,\rm reg}/H)_0<1.0 \times 10^{-6}$ (95% CI). Including all degrees of freedom simultaneously for the first time, anisotropic expansion of the Universe is strongly disfavoured, with odds of 121,000:1 against.