by Pavlos Alexandros Dimitriou, Matteo D’Alessandro, Brian L. Chang, Valentinos Silvestros, Elisavet Constantinou, Costas Pitris, Panayiotis Kolios, Piet Van MieghemUsing contact tracing data provided by the Cyprus Ministry of Health, infection trees for the first four waves of the COVID-19 epidemic are constructed. In these trees, nodes represent infected individuals, while links indicate the direction of transmission between them. For each infection tree of N nodes, the hopcount distribution from the root node to all other nodes is calculated. The empirical distribution is then compared to the hopcount distribution of infection trees generated by a non-Markovian SI process on a complete graph, with Weibull infection times characterized by a shape parameter α. We compute the values of the shape parameter α that best fit the empirical distribution and find that only values of α>1 are obtained, while the Markovian case is characterized by α=1. A Weibull distributed infection time with shape parameter α>1 is characterized by a unimodal density function with a peak at finite time, consistent with previous findings in the literature. Our analysis therefore suggests that the spreading process is most likely governed by non-Markovian dynamics, and that non-Markovianity can be detected solely from the topology of the infection trees. Finally, we analyze the evolution of the empirical distribution of the number of secondary infections caused by each node in the infection trees across different time windows to estimate the effective reproduction number. In practice, the average number of secondary infections seems to often provide a lower bound of the reproduction number computed by the Cyprus Ministry of Health. When the last level of the trees, composed predominantly of terminal nodes that do not generate further infections, is excluded, the estimate reflects more accurately the dynamics of the epidemic.