A versatile transformation system for thraustochytrids, a promising producer for polyunsaturated fatty acids and fatty acid-derived fuels, was established. G418, hygromycin B, blasticidin, and zeocin inhibited the growth of thraustochytrids, indicating that multiple selectable marker genes could be used in the transformation system. A neomycin resistance gene (neo(r)), driven with an ubiquitin or an EF-1α promoter-terminator from Thraustochytrium aureum ATCC 34304, was introduced into representatives of two thraustochytrid genera, Aurantiochytrium and Thraustochytrium. The neo(r) marker was integrated into the chromosomal DNA by random recombination and then functionally translated into neo(r) mRNA. Additionally, we confirmed that another two genera, Parietichytrium and Schizochytrium, could be transformed by the same method. By this method, the enhanced green fluorescent protein was functionally expressed in thraustochytrids. Meanwhile, T. aureum ATCC 34304 could be transformed by two 18S ribosomal DNA-targeting vectors, designed to cause single- or double-crossover homologous recombination. Finally, the fatty acid Δ5 desaturase gene was disrupted by double-crossover homologous recombination in T. aureum ATCC 34304, resulting in an increase of dihomo-γ-linolenic acid (C(20:3n-6)) and eicosatetraenoic acid (C(20:4n-3)), substrates for Δ5 desaturase, and a decrease of arachidonic acid (C(20:4n-6)) and eicosapentaenoic acid (C(20:5n-3)), products for the enzyme. These results clearly indicate that a versatile transformation system which could be applicable to both multiple transgene expression and gene targeting was established for thraustochytrids.