This is just a somewhat unusual generalised linear model with a Gaussian response $Y\sim N(\mu,\sigma^2)$ where the mean of the response $\mu$ is linked to the linear predictor $\eta$ via
$$
g(\mu)=\underbrace{\beta_0 + \beta_1 x}_\eta.
$$
where $g$ is the inverse link function $g(\mu)=1/\mu$,
The linear predictor $\eta$ is clearly linear in $x$ and, as implied by its name, in the regression coefficients $\beta_0$ and $\beta_1$. However, because of the non-linear link function, $\mu=EY$ is non-linear in the model parameters and in $x$.
To fit this in R do glm(y ~ x, gaussian(link = "inverse"))
.