I am using the deseq2 function plotPCA to visualize the principal components of my count data. I would like to extract the list of geneIDs that are contributing most to each component. I can get the value of PC1 and PC2 for each sample using returnData=TRUE, but I would like to extract the top and bottom genes from each component. Any ideas for me? Thank you!