How do I keep all genes after integration workflow?

[f6v]

This has probably already been asked before, but I can't find the right answer. I'm running the Seurat v5 integration workflow and would like to visualize some genes after the integration.

Seurat object before the integration contains CD79A gene:

'CD79A' %in% rownames(seurat_obj)
TRUE

I ran the integration:

seurat_obj[["RNA"]] <- split(seurat_obj[["RNA"]], f = seurat_obj$sample_id)
seurat_obj

seurat_obj <- SCTransform(seurat_obj, return.only.var.genes = T, verbose = F)
seurat_obj <- RunPCA(seurat_obj, npcs = 30, verbose = F)
seurat_obj <- IntegrateLayers(
  object = seurat_obj,
  method = RPCAIntegration,
  normalization.method = "SCT",
  verbose = F,
  k.weight = 40
)
seurat_obj <- FindNeighbors(seurat_obj, dims = 1:30, reduction = "integrated.dr", verbose = F)
seurat_obj <- FindClusters(seurat_obj, resolution = 0.3, verbose = F)
seurat_obj <- RunUMAP(seurat_obj, dims = 1:30, reduction = "integrated.dr", verbose = F)

The resulting object:

An object of class Seurat 
28856 features across 4792 samples within 2 assays 
Active assay: SCT (11996 features, 3000 variable features)
 3 layers present: counts, data, scale.data
 1 other assay present: RNA
 3 dimensional reductions calculated: pca, integrated.dr, umap

Now I try to plot CD79A expression:

FeaturePlot(seurat_obj, features = c('CD79A'))

And get the error:

Warning message:
“No layers found matching search pattern provided”
Warning message:
“The following requested variables were not found: CD79A”
Error in `FeaturePlot()`:
! None of the requested features were found: CD79A in slot  data
Traceback:

1. abort(message = paste("None of the requested features were found:", 
 .     paste(features, collapse = ", "), "in slot ", slot))
2. signal_abort(cnd, .file)
3. signalCondition(cnd)

I tied setting return.only.var.genes = F to return all genes, but it doesn't help. What am I missing?