The presence of a subduction signature with variable isotopic patterns in the melt products of mantle origin within the Aegean-Anatolian domain reflects the differential effects of slab-derived material on the chemical composition of the supraslab mantle region. Here, we use Sr-Nd-Pb isotopic compositions and major and trace element abundances to evaluate the characteristics of the source components and possible roles of fluid/melt metasomatism in the genesis of the Middle to Late Miocene mafic-intermediate lavas from the Soke Basin in western Anatolia. This volcanic suite comprises a number of individual lava flows in which the dominant rock types are potassium-rich basaltic andesites and andesites with trace element relative abundances indicative of melt derivation from subduction-modified sources. These samples have significantly heterogeneous Sr and Nd isotopic ratios, with Sr-87/Sr-86 and Nd-143/Nd-144 values ranging from 0.70475 to 0.70927 and from 0.51218 to 0.51270, respectively; they also yield a range of Pb isotope ratios, with Pb-206/Pb-204 = 18.963-19.344, Pb-207/Pb-204 = 15.672-15.733, and Pb-208/Pb-204 = 39.096-39.264. Remarkably large variations in isotopic ratios with relatively limited changes in major elements suggest the involvement of heterogeneously enriched sources with diverse isotopic signatures. The modelling of trace element and isotopic data indicates that parts of the lithospheric mantle in this region are isotopically enriched, with geochemical signatures that reflect the interaction of variably depleted mantle peridotites with melts from subducted sediments of continental crustal origin, while other parts host phlogopite-rich metasomatic components with signatures of devolatilization reactions from ancient subduction events that led to the preferential loss of Pb and Rb to create metasomatic domains with high U/Pb and low Rb/Sr values that ultimately generated highly radiogenic Pb-206/Pb-204 and relatively unradiogenic Sr-87/Sr-86 values. The observed isotopic diversity in relatively short-time- and length-scales may be interpreted to have resulted largely from the variable extent of melt contributions from a number of discrete lithological domains with contrasting long-term isotopic evolution and is consistent with the view that the continental lithosphere in this region comprises a number of disparate lithospheric fragments with different petrogenetic histories.