study include the availability of GWAS data in a large sample well-characterized in terms of psychiatric diagnoses, the use of different sub-phenotyping strategies to identify MDD subtypes and GPRS based on large international consortia. An important limitation is that the sample size was still largely underpowered to perform bivariate-GREML analyses to calculate the SNP correlation (average genome-wide relationship)(41) between the subtypes.