Suppose $X^1, \cdots, X^n \sim N_k(\mu, \Sigma)$, where $\mu$ is an unknown mean vector and $\Sigma$ is an unknown positive definite matrix.
Let $\hat \Sigma = \sum_{i=1}^{n}{(X^i-\bar X)( X^i-\bar X)^{\mathrm{T}}}/n$ be the MLE of $\Sigma$, where $\bar X = \sum_{i=1}^{n}{X^i}/n$.
How to prove it: $\mathbb{E}(\hat \Sigma) = (1-1/n)\Sigma$.