I expect not. For a counter example I would consider a differential operator where the coefficients have differing differentiability conditions so the formal adjoint might not exist in a normal sense such as out of the definition of the adjoint associated with a ordinary operator where the coefficients have no smoothness and the formal adjoint is a "quasi differential operator" in the verbiage of Naimark. Naimark covered this topic in detail in his books Linear Differential Operators.
That Naimark's work in differential equations was not better known is probably a victim of the Cold War. My wife at the time was a Russian linguist and I didn't have to depend on translations and not all of his better works were - but those two volumes above were.
In the thesis of G. C. Rota, a student of Jacob Schwartz, on the extension theory of ordinary differential operators he develops a theory which might point the way to the counter example you seek. His results were published in the Communications on Pure and Applied Mathematics in 1958. Another place to look is in the classic Dunford and Schwartz "Linear Operators", vol II where there is an detailed development of the extension theory of ordinary operators on Banach spaces (L^p).
My gut tells me without other conditions the conjecture is false and the counter example lies in differential operators.
Thank you very much for this detailed answer, which is also very informative.
In the first place, I thought a counterexample would be fairly simple, and I then thought that I just was not aware of it. Now, I think it has to be constructed in several steps. In Reed-Simon's (Vol. 2), the authors assumed somewhere $T^2$ is densely defined even for a symmetric operator $T$, and that has made me think that the answer to my question was probably negative (I have just emailed Prof. B. Simon who is probably aware of such a counterexample).
I will have a look at Naimark's paper and your other suggestions. If anything is found, I would certainly let you know.
If the answer to my question is positive (which is rather unlikely), then it would also be an interesting question to investigate the question for a closable and densely defined $T$?
This question spurred my interest. As for as closure of the formal adjoint - see Theorem 1.1 of Part II of Article Application of exponential dichotomies to asymptotic integra...
For background and if you cannot find Naimark's work see W. T. Reid's paper
is useful. When the coefficients of a differential operator are sufficient smooth so that D^j(ak(t) f(t)) can be calculated and the domains of definition of the formal adjoint are obvious. If not - then great care needs to be taken. The same would be true when one defines the square of such an operator if the coefficients are not sufficiently smooth. For a counter example, I expect one could find it in Strum-Louisville theory.
I wanted to tell you that I have already found a densely defined (unclosable) T such that $D(T^2)=D(T^*)=D(TT^*)=D(T^*T)=\{0\}$. Perhaps this example may be used in some other way (e.g. using a matrix of operators) to obtain a symmetric operator with the very above required properties.