A validated HPLC-DAD-ESI-MSn method for the analysis of non-anthocyanin flavonoids was applied to nine different tissues of twelve lotus genotypes of Nelumbo nucifera and N. lutea, together with an optimized anthocyanin extraction and separation protocol for lotus petals. A total of five anthocyanins and twenty non-anthocyanin flavonoids was identified and quantified. Flavonoid contents and compositions varied with cultivar and tissue and were used as a basis to divide tissues into three groups characterized by kaempferol and quercetin derivatives. Influences on flower petal coloration were investigated by principal components analyses. High contents of kaempferol glycosides were detected in the petals of N. nucifera while high quercetin glycoside concentrations occurred in N. lutea. Based on these results, biosynthetic pathways leading to specific compounds in lotus tissues are deduced through metabolomic analysis of different genotypes and tissues and correlations among flavonoid compounds.