linear mixed model for gene expression data in R

Question

I have gene expression data for multiple genotypes as TPM matrix:

gene    Sample1 Sample2 Sample3
gene A  12.6    3.5     6.7
gene B  7.8     9.7     4.3
<