So, after proving that $\mathfrak{B}(\mathbb{R})\subset \mathfrak{L}(\mathbb{R})$, I asked myself, and now asking you, is there a set $\mathfrak{S}(\mathbb{R})$, which satisfies:
$$\mathfrak{B}(\mathbb{R} )\subset \mathfrak{S}(\mathbb{R})\subset \mathfrak{L}(\mathbb{R})$$
($\mathfrak{B}(\mathbb{R})$ is the Borel's set on $\mathbb{R}$ ; $\mathfrak{L}(\mathbb{R})$ is family of sets which are Lebesgue's measurable - which is a $\sigma$ algebra.)
In ZFC the cardinality of $\frak B(\Bbb R)$ is $2^{\aleph_0}$ while the cardinality of $\frak L(\Bbb R)$ is $2^{2^{\aleph_0}}$. By that virtue alone there are plenty of other sets in between.
If you wish to take $\frak G(\Bbb R)$ to be some sort of a $\sigma$-algebra, and not just any set, let $\cal A\subseteq\frak L(\Bbb R)$ such that $|{\cal A}|<2^{2^{\aleph_0}}$, and consider $\frak G(\Bbb R)$ to be the $\sigma$-algebra generated by $(\frak B(\Bbb R)\cup\cal A)$. If we took $\cal A$ such that $\cal A\nsubseteq\frak B(\Bbb R)$ this would be a $\sigma$-algebra which strictly contains the Borel sets and is strictly contained in the Lebesgue measurable sets.
Edit:
Some time after Byron's comment below I realized that indeed this may not be accurate. I suddenly realized that I cannot be certain that $\frak L(\Bbb R)$ is not generated by a set of less than $2^{2^{\aleph_0}}$ many elements. It's not that bad, though. One can still bound it with some certainty:
We require that $|{\cal A}|^{\aleph_0}<2^{2^{\aleph_0}}$. If $|{\cal A}|\leq\frak c$ then this is indeed true, however one can come up with models where this need not be true for any subset of $\cal P(\Bbb R)$.
For what it's worth, I asked a question on MathOverflow some time ago, but did not receive any answer regarding this question yet. I still believe this is true, though.
There is a class of sets called analytic sets which are defined to be the continuous image of Borel sets (in fact $G_\delta$ sets are enough, but it turns out to be the same thing). It is a theorem that analytic sets properly contain the Borel sets and they are Lebesgue measurable.
Therefore the complement of analytic (co-analytic) sets are also Lebesgue measurable sets (they also contain all the Borel sets, and an amazing result is that a set is Borel if and only if it is both analytic and co-analytic.)
If you consider the $\sigma$-algebra generated by the union of all analytic and co-analytic sets, you will find yourself still inside the Lebesgue measurable universe, but with a strictly larger family of sets.
Further reading on how $\sigma$-algebras are born:
Beginning in the mid 1920s several (mostly) Russian mathematicians studied quite a number of $\aleph_{1}$-length hierarchies of sigma-algebras that fit between the analytic sets and the lowest projective set levels. I researched this literature fairly extensively about 5 years ago, but I've forgotten the details and don't have my findings available with me right now. It also seems I never wrote a sci.math essay about this topic, at least I can't find one. However, the following suggestions should allow you (and others interested) to begin looking into this topic, and when I get the chance (it may be a while), I'll add a lot more detail.
A useful survey paper is
Vladimir G Kanovei, Kolmogorov's ideas in the theory of operations on sets, Russian Mathematical Surveys 43 #6 (1988), 93-128.
Unfortunately, this paper does not seem to be freely available on the internet. However, 3 useful papers published in 1983 in Fundamenta Mathematica, by John P. Burgess and titled Classical hierarchies from a modern standpoint (Part I, Part II, Part III), are all freely available on the internet.
Finally, also worth searching for are the C sets of E. Selivanovskii (search also using the spelling "Selivanovskij"), sets that were first introduced and studied by Selivanovskii in a paper (written in Russian) published in pp. 379-413 of Matematiceskij Sbornik 35 (1928).
The class of Borel sets has the same size as the set of reals. The class of Lebesgue measurable sets has the same size as the power set of the reals.
To see this: Any subset of the Cantor set is Lebesgue measurable, with measure 0. Since the Cantor set has the same size as the reals, this shows that there are at least $|{\mathcal P}({\mathbb R})|$ many Lebesgue measurable sets, but of course there can be no more.
On the other hand, the Borel sets can be seen as the product of a transfinite construction of length $\aleph_1$, the first uncountable cardinal. Note that $\aleph_1\le|{\mathbb R}|$. In this construction, we start with the open sets (there are $|{\mathbb R}|$ many such sets), and at each stage add complements, and countable unions of sets from the previous stage. At limits, we take the union of the stages so far considered. None of these stages increases the cardinality of the class of sets collected so far, since $|{\mathbb R}|=|{\mathbb R}|^{\aleph_0}$.
Suppose now that $\Sigma$ is a $\sigma$-algebra of sets of reals, and $|\Sigma|^{\aleph_0}=|\Sigma|$. Let $A$ be any set of reals not in $\Sigma$, and let $\Lambda$ be the $\sigma$-algebra generated by $\{A\}\cup\Sigma$. Then $|\Lambda|=|\Sigma|$.
This shows that there are many $\sigma$-algebras in between the Borel sets and the Lebesgue measurable sets; in fact, you can create a very long chain of $\sigma$-algebras in the Lebesgue sets, $$\Sigma_0\subsetneq \Sigma_1\subsetneq\dots\subsetneq\Sigma_\alpha\subsetneq\dots,$$ with $\Sigma_0$ the algebra of Borel sets: For $\alpha<|{\mathbb R}|^+$, the algebras all have the same size as the reals. You may go longer or stop when you reach $|{\mathbb R}|^+$, depending on the cardinality of the power set of the reals.
A less sophisticated, but potentially useful answer: You know that the intersection of $\sigma$-algebras is again a $\sigma$-algebra. If you start with $\mathfrak{B}(\mathbb{R})$ you can put some probability measure on it and take the completion under this probability measure. Doing this for every probability measure, you get a large family of $\sigma$-algebras extending $\mathfrak{B}(\mathbb{R})$. Their intersection, which you might call $\mathfrak{U}(\mathbb{R})$ is the $\sigma$-algebra of universally measurable sets. The analytic sets mentioned above are all universally measurable. This is an extremely useful fact.
In some applications such as statistics, you want t work with a whole family of measures on a $\sigma$-algebra. Now by construction, you can extend every probability measure on $\mathfrak{B}(\mathbb{R})$ in a unique way to all of $\mathfrak{U}(\mathbb{R})$. If $A\subseteq\mathfrak{B}(\mathbb{R}^2)$, the projection on the first coordinate will in general not be in $\mathfrak{B}(\mathbb{R})$, but it will be analytic and hence in $\mathfrak{U}(\mathbb{R})$. It is possible to rewrite optimization problems in terms of such projections. So let $f:A\to\mathbb{R}$ be a measurable function. Define a function $m:\mathbb{R}\to\mathbb{R}$ by $m(r)=\inf\{t:(r,t)\in A\}$. The infimum is understood to be in the extended real numbers so that $m(r)=\infty$ if there i no $t$ such that $(r,t)\in A$. In general, $m$ will not be $\mathfrak{B}(\mathbb{R})$-measurable, but $\mathfrak{U}(\mathbb{R})$-measurable. To see the latter, we note that we can rewrite $\{r:m(r)<\alpha\}$ as the projection of $A\cap\mathbb{R}\times(-\infty,\alpha)$ onto the first coordinate, a set that is analytic and hence in $\mathfrak{U}(\mathbb{R})$. This idea has very important applications. In stochastic dynamic programming, it can be that the value function is not $\mathfrak{B}(\mathbb{R})$-measurable, but $\mathfrak{U}(\mathbb{R})$-measurable for essentially the reason I just gave.
Bumping an old question, I think this is a situation where the language of elementary submodels is convenient.
It's easy to construct a set $M$ with the following properties (e.g. via repeated Lowenheim-Skolem to build an elementary chain of length $\omega_1$):
$M\prec V$ (or $M\prec H_\kappa$ for "large enough" $\kappa$ if one prefers),
$M$ is closed under $\omega$-sequences, and
$\vert M\vert = \mathfrak{c}$ (of course $\vert M\vert<2^\mathfrak{c}$ is really all we need).
These properties ensure that $M$ contains every real and every Borel set, does not contain every Lebesgue measurable set, and satisfies enough of the true theory of the universe for us to relativize basic concepts and facts to $M$. In particular, we automatically get the following:
$\mathfrak{L}(\mathbb{R})^M=\mathfrak{L}(\mathbb{R})\cap M$ is a $\sigma$-algebra strictly between $\mathfrak{B}(\mathbb{R})$ and $\mathfrak{L}(\mathbb{R})$.
Note that it would not be enough to replace condition $(2)$ by merely "$\mathfrak{B}(\mathbb{R})\subseteq M$" - we also used $\omega$-closedness to ensure that anything $M$ thinks is a $\sigma$-algebra really is a $\sigma$-algebra.
Admittedly, this isn't really that different from the existing answers since we're hiding the "close under appropriate operations" bit in the construction of $M$ itself. However, I do think it makes things simpler to think about, and yields intermediate $\sigma$-algebras of a somewhat different flavor.
