In kidney transplantation, molecular diagnostics may be a valuable approach to improve the precision of the diagnosis. Using next-generation sequencing (NGS), we aimed to identify clinically relevant archetypes. We conducted an Illumina bulk RNA sequencing on 770 kidney biopsies (540 kidney recipients) collected between 2006 and 2021 from 11 European centers. Differentially expressed genes were determined for 11 Banff lesions. An ElasticNet model was used for feature selection, and 4 machine learning classifiers were trained to predict the probability of presence of the lesions. NGS-based classifiers were used in an unsupervised archetypal analysis to different archetypes. The association of the archetypes with allograft survival was assessed using the iBox risk prediction score. The ElasticNet feature selection reduced the number of the genes from a range of 859-10 830 to a range of 52-867 genes. NGS-based classifiers demonstrated robust performances (precision-recall area under the curves 0.708-0.980) in predicting the Banff lesions. Archetypal analysis revealed 8 distinct phenotypes, each characterized by distinct clinical, immunological, and histological features. Although the archetypes confirmed the well-defined Banff rejection phenotypes for T cell-mediated rejection and antibody-mediated rejection, equivocal histologic antibody-mediated rejection, and borderline diagnoses were reclassified into different archetypes based on their molecular signatures. The 8 NGS-based archetypes displayed distinct allograft survival profiles with incremental graft loss rates between archetypes, ranging from 90% to 56% rates 7 y after evaluation (P < 0.0001). Using molecular phenotyping, 8 archetypes were identified. These NGS-based archetypes might improve disease characterization, reclassify ambiguous Banff diagnoses, and enable patient-specific risk stratification.