Is it reasonable to use natural orbitals of an unrestricted HF calculation as a basis to start a CASSCF calculation for a radical
The best starting orbitals for a CASSCF calculation are optimised orbitals of another CASSCF (or RASSCF) calculation. That sounds a bit ridiculous, but this is probably the best way to figure out the active space.

At first you are probably choosing something quite small, like CAS(2,2) to CAS(4,4). For these calculations it barely matters what kind of starting orbitals you are choosing. The change will be negligible and probably depending heavily on the system you are using. 

The advantage of calculating the cations by dismissing the radicals is often a good choice since you can run a very fast RHF calculation and the orbitals are symmetric, which is what you need for CAS.
Natural orbitals of an UHF calculation are also a quite sensible choice, because they already allow for fractional occupation and should also be symmetric. The only downside I can think of is, that this kind of start orbital generation takes a lot longer.
Similarly you can choose to run a ROHF calculation, if you are uncomfortable with too much cationic charge.
As you have said you can choose also an UHF calculation by throwing away one set. I would advise caution here. Sometimes - probably especially when CAS is necessary - spin contamination in UHF is significant. Therefore one set might not give you anything close and you could end up choosing the wrong CAS.

I am not aware of any major pitfalls except for the just mentioned choosing of a completely wrong active space. 

In any case it is a good idea to give the CASSCF calculation a few cycles and see how the active space is developing. Sometimes you need to adjust the active space long before the calculation converges and then you will have already a much better guess than with your initial starting orbitals.

I personally use the following routine 


R(O)HF [x+] -> CAS(x, &lceil;&frac12;x&rceil; +1) 
            -> RAS(x+2n, &lceil;&frac12;x&rceil; +n +1) 
            -> CAS(x+2n, &lceil;&frac12;x&rceil; +n +1) 
            -> RAS(...) 
            -> CAS(...) -> ...


to narrow down the active space and include all orbitals I want and need. In most cases I am not converging the intermediate calculations, except for the minimal CAS.
As you can see the procedure that follows after your start orbitals is much&#8482; more time consuming, so it really does not matter.