Next‐generation sequencing of cell‐free circulating DNA (cfDNA) has emerged as promising technique for identifying minimally invasive genomic profiling of tumor cells recently. However, it remains relatively unknown in LAM disease. In our study, paired cfDNA and genomic DNA (gDNA) in blood samples were obtained from 23 LAM patients and seven healthy controls to explore mutations profiles of targeted 70 cancer‐related genes. As results, log2‐based allele frequencies of mutations in cfDNA were significantly different from those of gDNA. By comparing the mutual mutations identified both in cfDNA and gDNA, a significant correlation was also observed. After removing mutations in gDNA, distinct somatic mutation profiles of cfDNA were observed in LAM patients. Forty of 70 targeted genes had recurrent mutations, of which ATM, BRCA2 and APC showed the highest frequency. Based on the mutation, correlation network constructed of 40 mutated genes, 11 hub genes bearing intensive interactions were highlighted, including BRCA1, BRCA2, RAD50, RB1, NF1, APC, MLH3, ATM, PDGFRA, PALB2 and BLM. Expression of the hub genes showed significant clusters between LAM patients and controls and that RAD50 and BRCA2 had the strongest associations with subject phenotypes. Myogenesis and estrogen response were confirmed to be positively regulated in LAM patients. Collectively, our study provided a landscape of genomic alterations in LAM and discovered several potential driver genes, that is, BRCA2 and RAD50, which shed a substantial light on the clinical application of key molecular markers and potential therapy targets for precision diagnosis and treatment in the future.