Universal behavior related to continuous symmetry breaking in nematic liquid crystals is studied using Brownian molecular dynamics. A three-dimensional lattice system of rod-like objects interacting via the Lebwohl-Lasher interaction is considered. We test the applicability of predictions originally derived in cosmology and magnetism. In the first part we focus on coarsening dynamics following the temperature driven isotropic-nematic phase transition for different quench rates. The behavior in the early coarsening regime supports predictions made originally by Kibble in cosmology. For fast enough quenches, symmetry breaking and causality give rise to a dense tangle of defects. When the degree of orientational ordering is large enough, well defined protodomains characterized by a single average domain length are formed. With time subcritical domains gradually vanish and supercritical domains grow with time, exhibiting a universal scaling law. In the second part of the paper we study the impact of random-field-type disorder on a range of ordering in the (symmetry broken) nematic phase. We demonstrate that short-range order is observed even for a minute concentration of impurities, giving rise to disorder in line with the Imry-Ma theorem prediction only for the appropriate history of systems.