leading log2 fold change (FC) for the most variable genes between each pair of samples (Ritchie et al. 2015). Comparisons were adjusted for sex by implementing a batch factor according to the edgeR manual. One sample, a HAP male, was identified as an outlier and was found to have a much smaller library size with significantly more reads aligning to mitochondrial sequences when compared to other samples. Therefore, this sample was removed from the analyses. Full processed datasets and analyses may be found at https://github.com/dlhagger/AUD-Multi-Omics.