diff --git a/src/sage/libs/symmetrica/symmetrica.pxi b/src/sage/libs/symmetrica/symmetrica.pxi
index 7244b7e9655..06330bffbd4 100644
--- a/src/sage/libs/symmetrica/symmetrica.pxi
+++ b/src/sage/libs/symmetrica/symmetrica.pxi
@@ -401,7 +401,7 @@ cdef void late_import():
            SymmetricFunctions, \
            sqrt, \
            builtinlist, \
-           MPolynomialRing_base, is_MPolynomial,\
+           MPolynomialRing_base, MPolynomial,\
            SchubertPolynomialRing, SchubertPolynomial_class,\
            two, fifteen, thirty, zero, sage_maxint
 
@@ -454,7 +454,7 @@ cdef void late_import():
     import sage.rings.polynomial.multi_polynomial_ring
     MPolynomialRing_base = sage.rings.polynomial.multi_polynomial_ring.MPolynomialRing_base
     import sage.rings.polynomial.multi_polynomial_element
-    is_MPolynomial = sage.rings.polynomial.multi_polynomial_element.is_MPolynomial
+    MPolynomial = sage.structure.element.MPolynomial
 
     import sage.combinat.schubert_polynomial
     SchubertPolynomialRing = sage.combinat.schubert_polynomial.SchubertPolynomialRing
diff --git a/src/sage/rings/polynomial/multi_polynomial.pyx b/src/sage/rings/polynomial/multi_polynomial.pyx
index 9ce51a354b3..f4bf728d739 100644
--- a/src/sage/rings/polynomial/multi_polynomial.pyx
+++ b/src/sage/rings/polynomial/multi_polynomial.pyx
@@ -18,6 +18,9 @@ from sage.misc.derivative import multi_derivative
 from sage.misc.misc_c import prod
 
 def is_MPolynomial(x):
+    from sage.misc.superseded import deprecation
+    deprecation(32709, "the function is_MPolynomial is deprecated; use isinstance(x, sage.structure.element.MPolynomial) instead")
+
     return isinstance(x, MPolynomial)
 
 from sage.rings.polynomial.polynomial_ring_constructor import PolynomialRing
diff --git a/src/sage/rings/polynomial/polynomial_element.pyx b/src/sage/rings/polynomial/polynomial_element.pyx
index f73b13bebc5..852cfb13fe7 100644
--- a/src/sage/rings/polynomial/polynomial_element.pyx
+++ b/src/sage/rings/polynomial/polynomial_element.pyx
@@ -142,20 +142,24 @@ cpdef is_Polynomial(f):
     """
     Return True if f is of type univariate polynomial.
 
+    This function is deprecated.
+
     INPUT:
 
     -  ``f`` -- an object
 
     EXAMPLES::
 
-        sage: from sage.structure.element import Polynomial
+        sage: from sage.rings.polynomial.polynomial_element import is_Polynomial
         sage: R.<x> = ZZ[]
-        sage: isinstance(x^3 + x + 1, Polynomial)
+        sage: is_Polynomial(x^3 + x + 1)
+        doctest:...: DeprecationWarning: the function is_Polynomial is deprecated; use isinstance(x, sage.structure.element.Polynomial) instead
+        See https://trac.sagemath.org/32709 for details.
         True
         sage: S.<y> = R[]
         sage: f = y^3 + x*y -3*x; f
         y^3 + x*y - 3*x
-        sage: isinstance(f, Polynomial)
+        sage: is_Polynomial(f)
         True
 
     However this function does not return True for genuine multivariate
@@ -165,15 +169,18 @@ cpdef is_Polynomial(f):
         sage: R.<x,y> = QQ[]
         sage: f = y^3 + x*y -3*x; f
         y^3 + x*y - 3*x
-        sage: isinstance(f, Polynomial)
+        sage: is_Polynomial(f)
         False
         sage: var('x,y')
         (x, y)
         sage: f = y^3 + x*y -3*x; f
         y^3 + x*y - 3*x
-        sage: isinstance(f, Polynomial)
+        sage: is_Polynomial(f)
         False
     """
+    from sage.misc.superseded import deprecation
+    deprecation(32709, "the function is_Polynomial is deprecated; use isinstance(x, sage.structure.element.Polynomial) instead")
+
     return isinstance(f, Polynomial)
 
 from .polynomial_compiled cimport CompiledPolynomialFunction