Genomic clusters of immune genes, including those encoding nucleotide-binding leucine-rich repeat (NLR) proteins, are a model for exploring the dynamics of genomic regions in flux. Rapid sequence evolution of immune genes, including NLRs, and variation in their gene content, may enable long-lived plants, which lack adaptive immune systems, to keep pace with the fast evolution of pathogens. To explore the patterns and processes shaping the evolution of NLR gene content in a genus of long-lived tree species, we unified the annotation of NLR genes across 19 accessions (or 28 haplotypes) from the genus Citrus and its relatives, including three new diploid genome assemblies. A majority of NLRs were arranged in genomic clusters composed of paralogous genes, typically from a single gene family. Even larger clusters, with 10 or more NLRs, were limited to genes derived from one or a few gene families. These patterns suggested that genomic clustering of NLRs arose through local expansion of phylogenetically related NLRs, but the mechanistic processes driving these patterns are not clear. Local gene duplication can be mediated by multiple processes, including transposon-mediated gene capture and subsequent proliferation, and non-allelic repair of double-stranded breaks, including unequal recombination. Examples of retrotransposon-mediated duplication of NLRs were identified, but these were not sufficient to explain massive regional expansions. Signatures of unequal recombination are challenging to identify. Focusing on recent lineage-specific sequence duplications, at least one case of unequal recombination was identified, supporting a role for unequal recombination in shaping genomic variation in these regions.