Gray matter (GM) forms the ‘computational engine’ of our brain and plays the key role in brain function. Measures derived from MRI (e.g., structural MRI (sMRI) and diffusion MRI (dMRI)) provide a unique opportunity to non-invasively study GM structure in-vivo and thus can be used to probe GM pathology in development, aging and neuropsychiatric disorders. Investigation of the influence of various factors on MRI measures in GM is critical to facilitate their use for future non-invasive studies in healthy and diseased populations. In this dissertation, GM structure was studied with MRI to understand how it is influenced by genetic and environmental factors. Validation of dMRI- derived measures was conducted by comparing them with histological data from monkeys to better understand the cytoarchitectural features that influence GM measures. First, the influence of genetic and environmental factors was quantified on gray matter macrostructure and microstructure measures using phenotypic modelling of structural and diffusion MRI data obtained from a large twin and sibling population (N = 840). Results of this study showed that in GM, while sMRI measures like cortical thickness and GM volume are mainly affected by genetic factors, advanced dMRI measures of mean squared displacement (MSD) and return to origin probability (RTOP) derived from advanced biexponential model can tap into regionally specific patterns of both genetic and environmental influence in cortical and subcortical GM. Our results thus highlight the potential of these advanced dMRI measures for use in future studies that aim to investigate and follow in healthy and clinical population changes in GM microstructure linked with both genes and environment. Second, using data from a large healthy population (n=550), we investigated changes in sMRI tissue contrast at the gray-white matter boundary with biological development during adolescence to assess how this affects estimation of the developmental trajectory of cortical thickness. Results of this study suggest that increased myelination during brain development contributes to age-related changes in gray-white boundary contrast in sMRI scans causing an apparent shift of the estimated gray-white boundary towards the cortical surface, in turn reducing estimations of cortical thickness and its developmental trajectory. Based on these results, we emphasize the importance of accounting for the effects of myelination on T1 gray-white matter boundary contrast to enable more precise estimation of cortical thickness during neurodevelopment. Finally, we conducted histological validation of dMRI measures in gray matter by comparing dMRI measures derived from two models, conventional Diffusion Tensor Imaging (DTI) model and an advanced biexponential model with histology acquired from the same 4 rhesus monkeys. Results demonstrate differences in the ability of distinct dMRI measures including DTI-derived measures of fractional anisotropy (FA), Trace and advanced Biexponential model-derived measures of MSD and RTOP to capture the biological features of underlying cytoarchitecture and identify the dMRI measures that best reflect underlying gray matter cytoarchitectural properties. Investigation of the contribution of underlying cytoarchitecture (cellular organization) to dMRI measures in gray matter provides validation of dMRI measures of average and regional heterogeneity in MSD & Trace as markers of cytoarchitecture as measured by regional average and heterogeneity in cell area density. This postmortem validation of these dMRI measures makes their use possible for treatment monitoring of various GM pathologies. These studies and their results together demonstrate the utility of imaging measures to investigate the complex relationships between GM cellular organization, brain development, environment and genes.