We compute Euclidean Wilson loops for BMN Plane Wave Matrix Models by probing 1/2 BPS Type IIA geometries with F strings and D2 branes. These geometries fall under the general category of Lin and Maldacena (LM). As a special class of solutions within the LM category, we construct Wilson operators corresponding to the non-Abelian T duals (NATDs) of AdS 5 × S 5 . Considering various open string embeddings within Type IIA bulk, we show the existence Euclidean Wilson loops in the theory some of which further break the N = 2 SUSY of the BMN matrix model and are in fact turn out to be 1 8 BPS. Moreover, considering the 11d SUGRA uplift of the NATDs of AdS 5 × S 5 , we compute Wilson loops using M2 branes and also comment on the supersymmetry preserved by these loops. Finally, we set an argument for Wilson loop calculations within the framework of BMN matrix model.