Three-dimensional direct numerical simulation (DNS) of turbulent combustion under moderate and intense low-oxygen dilution (MILD) conditions has been carried out inside a cuboid with inflow and outflow boundaries on the upstream and downstreamfaces respectively. The initial and inflowing mixture and turbulence fields are constructed carefully to be representative of MILD conditions involving partially mixed pockets of unburnt and burnt gases. The combustion kinetics is modelled using a skeletal mechanism for methane-air combustion, including non-unity Lewis numbers for species and temperature dependent transport properties. The DNS data is analysed to study theMILD reaction zone structure and its behaviour. The results show that the instantaneous reaction zones are convoluted and the degree of convolution increases with dilution and turbulence levels. Interactions of reaction zones occur frequently and are spread out in a large portion of the computational domain due to the mixture non-uniformity and high turbulence level. These interactions lead to local thickening of reaction zones yielding an appearance of distributed combustion despite the presence of local thin reaction zones. A canonical MILD flame element, called as MIFE, is proposed which represents the averaged mass fraction variation for major species reasonably well, although a fully representative canonical element needs to include the effect of reaction zone interactions and associated thickening effects on the mean reaction rate.