In principle, we could directly check for significant differences in the demographic groups and use Bonferroni correction to control the family-wise error (i.e., the probability of obtaining a false positive).
given the large number of treelets, the correction for multiple comparisons would underpower our analyses and potentially cause us to miss many significant differences.