One conclusion is that an impenetrable solenoid with a given flux can be replaced by a space R³ with a line R removed, i.e. R³ \ R on which a non-vanishing A-field exists. This A-field is "pure gauge" locally, i.e. A ~ U* dU with dA = 0, but not globally b/c the line R prevents us from constructing a global gauge function U' to gauge away this A → A' = 0 globally.
This can be extended to more than one flux line R, R', R'', ..., i.e. several lines removed: R³ \ R \R' \ R'' ...
Then there is a non-local flux observable ∫A which corresponds to a topological quantity, the cohomology of the one-form A of the fibre bundle of the U(1) gauge group over the base manifold R³ \ R. This ∫A depends on A and on the the path C via an integer winding number w[C] of the path C w.r.t. the line R.
This can be expressed in terms of the gauge function U = exp(iχ) from which A can be constructed as U* dU, therefore already this U carries the full information on the topological structure.
Note that in principle it's not ∫A which is observed directly but only the fractional part of the flux mod 2π/q (where q is the value of the electric charge of the particle). This is due to the fact that we study the compact gauge group U(1) with U ~ exp(iχ), not only the gauge function χ. This is related to the problem of the logarithm χ of U on R³ \ R. There are continuous single-valued gauge functions U where χ jumps by 2πw along a path C for winding number w[C], the simplest one is χ ~ θ where θ is the cylindrical angle coordinate w.r.t. the z-axis.
So this results in a quite complex structure of the gauge group defined over R³ \ R; one may see this by applying the following gauge fixings:
1) gauging away A° = 0
2) gauging away A³ = 0 (we assume that the line R corresponds to the z-axis) with a time-indep. gauge function respecting A° = 0
(up to now this is standard)
3) now there is a residual gauge symmetry with t- and z-indep. gauge transformations U respecting A° = 0 and A³ = 0 such that the integral part of the flux is still unobservable, i.e. related to a global residual gauge symmetry
So the Aharonov-Bohm effect is related to a physical observable
\phi_C = \oint_C A \;\text{mod}\; 2\pi
which is the "fractional part of the flux".
This is defined in terms of a pure gauge vector potential
A = U^\dagger i\nabla U
which is sensitive to the topology of the base manifold, i.e. we have
\phi_C = \phi_C
So neither the base manifold nor the gauge group alone are sufficient top describe the Aharonov-Bohm effect; what really matters is the topology of the fibre bundle F which locally (i.e. for each point of the base manifold M) looks like a direct product F ~ M * U(1); but a non-trivial structure of M, i.e. R³ \ R induces a non-trivial structure of F, and we finally arrive at the conclusion that the physical observable is sensitive to
\phi_C = \phi_C[F]