In any gauge theory in any spacetime dimension, it is possible to construct a class of operators called the Wilson operators. For every loop $C$ in spacetime (i.e. for every embedding of a circle in spacetime), it is possible to define a Wilson operator as follows. Let us fix some classical configuration $A$ of the gauge field. Consider the holonomy $Hol(A, C)$ around $C$: if we fix a point $c$ of $C$, it is the element of the gauge group $G$ which defines the action of the parallel transport starting from $c$, following $C$ and returning to $c$. Formally, we have
$Hol(A,C) = Texp(\int_C A)$
where Texp denotes the path-ordered exponential.
The element $Hol(A,C)$ of $G$ depends of the choice of the base point $c$ but not its conjugation class. So, in order to obtain a well-defined (gauge invariant) quantity, a way to do is to take the trace $Tr_R Hol(A,C)$ in some representation R of the gauge group. For a fixed loop $C$ and a fixed representation $R$, $A \mapsto Tr_R Hol(A,C)$ is a gauge invariant function on the space of gauge fields. In particular, it is possible to take its vacuum expectation value in the quantum theory (say in Euclidean signature);
$W(C,R) = <Tr_R Hol(A,C)> = \int DA e^{-S(A)} Tr_R Hol(A,C)$
where $S(A)$ is the action of the gauge theory, $\int DA$ is the functional integral on the space of gauge fields.
Physically, the Wilson operators $Tr_R Hol(A,C)$ are of central importance to understand the dynamics of gauge theories and are kind of "test objects" to detect various phenomena such as confinement.
If one is interested by knots, one can try to apply what preceeds by taking a 3-dimensional spacetime $M$ (such as the 3-sphere) and then $C$ is simply a knot in $M$ (the possibility for a loop to be knotted is special to the dimension 3). For a general gauge theory, a quantity such as $W(C,R)$ will depend on the precise embedding of $C$ on $M$ and not only of the "topology" (more precisely the isotopy class) of the embedding. A knot theorist is interested by knot invariants, i.e. quantity associated to C depending only of the topology of the embedding. To construct such quantities using the preceding lines, we have to take a topological gauge theory. The simplest topological field theory in 3 dimensions is Chern-Simons theory, of action
$S(A,k) = k \int_M Tr(AdA + \frac{2}{3}A^3)$
where the coupling constant $k$ must be an integer (up to some universal normalization) in order to obtain a gauge invariant quantum field theory.
So, given a knot in the 3-sphere, a choice of gauge group $G$, a choice of representation $R$, and a choice of integer $k$, it is possible to define some knot invariant $W(C,R)$. For $G=SU(2)$, $R$ the fundamental representation, and $k$ varying, we obtain the famous Jones polynomial. The starting point of this story is the seminal paper of Witten http://projecteuclid.org/euclid.cmp/1104178138
From a purely mathematical point of view, the situation is not very satisfying because we don't know how to rigorously define the path integral defining the quantum Chern-Simons theory. Still, the Witten's paper had some enormous influence in knot theory and many things have been done: the Chern-Simons invariants can be rigorously defined combinatorially, it is possible to study the perturbative expansion of Chern-Simons theory...
Furthermore, the relation between knot theory and physics is still currently very active: by Witten, Chern-Simons theory is the same thing as the A-model of open topological string on the cotangent bundle of the 3-sphere, then by applying the large N transition of Gopakumar-Vafa, we obtain a relation with closed string amplitudes on the conifold... Some random names on these questions: Dijkgraaf, Aganagic, Marino, Vafa, Gukov... For another direction, more related to gauge theory, but this time in 4 or 5 dimensions, see the work of Witten on Khovanov homology http://arxiv.org/abs/1108.3103