The reciprocal basis and dual basis are essentially the same. Here is why.
To explicate the "abstract, symmetric definition of duality" mentioned in the comments, the idea here is that an inner product space is canonically isomorphic to its dual under the map $v \mapsto \langle v, \_\rangle$. Since we are in $\mathbb R^n$, there is a canonical inner product given by the dot product (in general, this comes from the Rietz Representation Theorem).
Thus, in an inner product space $V$ with inner product $\langle \_, \_ \rangle$ and orthonormal basis $e_i$ (meaning $\langle e_i, e_j \rangle = \delta_i^j$, i.e. take a set of basis vectors, make them orthogonal via Gram-Schmidt, then normalize) we get a dual basis $e^i$ in $V^*$ given by $e^i = \langle e_i, \_ \rangle$. Indeed, $e^i(e_j) = \langle e_i, e_j \rangle = \delta_i^j$ as desired.
Thus the dot product gives us an isomorphism from $\mathbb R^n$ to its dual, meaning $e^i$ should themselves be a basis for $\mathbb R^n$. As such, we can calculate $e^i$ using a change of basis matrix.
While this was a bit abstract, let us now think about how to represent dual vectors a bit more concretely in $\mathbb R^n$. To start, we usually start by writing a vector in a finite dimensional space as a column vector, e.g. in $\mathbb R^3$ we write $e_2 = [0,1,0]^T$. We then write dual vectors as row vectors, since it is forced by compatibility with the identification of linear maps with matrices. (Indeed, this is essentially forced from the fact we want to write $Ax$ for a matrix $A$ and vector $x$ instead of $xA$.) Recall the standard representation of a linear map from $\mathbb R^n$ to $\mathbb R^m$, or more generally from an $n$ dimensional vector space to an $m$ dimensional vector space, is given by an $m \times n$ matrix (whose columns are the images of the basis vectors). Since a dual vector is a linear map from $\mathbb R^n$ to $\mathbb R$, it must be represented by a $1 \times n$ matrix, i.e. a row vector.
Next, note the way matrix multiplication works means a row vector acting a column vector is basically just a dot product. Indeed, $[a_1, \ldots, a_i][b_1, \ldots, b_i]^T = a_1 b_1 + a_2 b_2 + \cdots + a_i b_i = (a_1, \ldots, a_i)\cdot (b_1, \ldots, b_i)$. We can thus think of vectors as the things on the "right" of a dot product, and dual vectors as the things on the "left". Thinking back to our earlier isomorphism $e^i = \langle e_i, \_ \rangle$, since our inner product is here given by the dot product, the dual vector $e^i$ is simply given by $e^i(v) = e_i \cdot v$
The reason this is called the "reciprocal basis" comes from how we actually calculate the dual basis from the change of basis matrix. Indeed, given some basis $\{b_i\}$ of a vector space $V$ of dimension $n$ we can form the matrix $P = [b_1, \ldots, b_n]$ and the matrix $B = [b_1^*, \ldots, b_n^*]$, where $\{b_i^*\}$ is the dual basis (to $\{b_i\}$) which we seek (the choice of the letter $P$ comes from it being a change of basis matrix, explained below). We now consider the transpose matrix $B^T$ whose rows are the vectors $b_i^*$. Since matrix multiplication takes rows on the left and dot products them with columns on the right and $b_i^* \cdot b_j = \delta^i_j$ per our earlier discussion, we find:
$$B^T P = [b_1^*, \ldots, b_n^*]^T [b_1, \ldots, b_n] = \begin{bmatrix} b_1^* \cdot b_1 & b_1^* \cdot b_2 & \cdots & b_1^* \cdot b_n \\ b_2^* \cdot b_1 & b_2^* \cdot b_2 & \cdots & b_2^* \cdot b_n \\ \vdots & \vdots & \ddots & \vdots \\ b_n^* \cdot b_1 & b_n^* \cdot b_2 & \cdots & b_n^* \cdot b_n \end{bmatrix} = \begin{bmatrix} 1 & 0 & \cdots & 0 \\ 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & 1 \end{bmatrix}$$
i.e. $B^T P = I$ so $B = (P^{-1})^T$. In other words, the dual basis we seek is given by the columns of $(P^{-1})^T$, or equivalently (and more in line with our view of dual vectors as row vectors) is given by the rows of $P^{-1}$. Since taking the inverse of a real number is often called taking its "reciprocal", by analogy these dual vectors are sometimes called the "reciprocal basis" since they come from taking the inverse.
Finally, note this is the same thing we would get from the change of basis matrix idea mentioned earlier, and to get there we need a few facts. First, recall that if $T \colon A \to B$ is a linear map with matrix $M$ then the tranpose of $T$ is the map $T^* \colon B^* \to A^*$ given by $f \mapsto f \circ T$ whose matrix representation in the dual basis is given by $M^T$ (note the transpose matrix "flips" the vector spaces along with dualizing them). Second, recall that if $T$ is invertible then $M^{-1}$ is the matrix associated to $T^{-1} \colon B \to A$ (note the inverse also "flips" the spaces). Finally, between any two bases $B_1$ and $B_2$ on a vector space $V$, there exists a change of basis matrix $P$ whose columns are vectors $B_1$ expressed in terms of $B_2$. Multiplication on the left by $P$ takes a vector expressed in $B_1$ coordinates and returns the same vector expressed in $B_2$ coordinates, so $P$ is referred to as the identity map from $(V, B_1)$ to $(V, B_2)$
Proceeding, let $\scr E$ be the canonical basis $\{e_i\}$ on $V$ with dimension $n$ and $\scr B$ be the given basis $\{b_i\}$ at hand. Since we seek the coordinates for the dual basis $\{b_i^*\}$ in terms of the canonical dual basis $\{e_i^*\}$, what we need to calculate is the identity map from $(V^*,\scr B^*)$ to $(V^*, \scr E^*)$. To this end, note the the identity map from $(V,\scr B)$ to $(V, \scr E)$. is given by $P := [b_1, \ldots, b_n]$, i.e. the columns are the vectors $b_i$. We would like the dualize immediately via the transpose, but this would flip the spaces; to fix this, we flip them an extra time via the inverse. Indeed, $P^{-1}$ is the identity map from $(V,\scr E)$ to $(V, \scr B)$ so $(P^{-1})^T$ is the identity map from $(V^*,\scr B^*)$ to $(V^*, \scr E^*)$. We once again conclude the columns of $(P^{-1})^T$ are the dual basis to $\{b_i\}$, i.e. that the rows of $P^{-1}$ are the dual basis to $\{b_i\}$.
As an example, note the claim made on Wikipedia that
$$b^1 = \left(\frac{b_2 \times b_3}{V}\right)^T,\;\;
b^2 = \left(\frac{b_3 \times b_1}{V}\right)^T,\;\;
b^3 = \left(\frac{b_1 \times b_2}{V}\right)^T.$$
(where $V$ is the volume of the parallelepiped spanned by the $b_i$) simply follows from the following formula for the inverse of a $3 \times 3$ matrix:
$$P^{-1} = [b_1, b_2, b_3]^{-1} = \frac{1}{\det P} [b_2 \times b_3, b_3 \times b_1, b_1 \times b_2]^T$$
(Note $V = \det P$ by the standard geometric interpretation of the determinant.) For a proof of this formula you can use the standard expression for the cross product as a determinant, which breaks into a sum of cofactors. When one writes out the formula for the the cofactor matrix, the columns are immediately seen to be the cross products, and the inverse matrix is classically given by the tranpose of the cofactor matrix scaled by the multiplicative inverse of the determinant. Alternatively, a more geometric proof can be found here.
Tl;dr:
Q) How should I reconcile these superficially similar descriptions of what seem to be fundamentally different things?
A) They are the same thing, depending on what you mean by "same" (equivalent up to some dualities). It is thus expected they have the same description.
Q) Are they linked by an isomorphism between the dual basis in the dual space and the reciprocal basis in original vector space?
A) Yes.
Q) Does this relate to the canonical isomorphism (independent of basis) between a vector space and the dual of its dual (something I am still trying to get my head around)?
A) Not in any meaningful way.