RAPD markers and sequences of chloroplast DNA (cpDNA) atp B-rbc L intergenic spacers were used to characterize the pattern of genetic variation and the phylogenetic relationships of the relict populations of Alsophila spinulosa located in Jian Feng Ling (JFL) and Diao Luo Shan (DLS), Hainan, and Tang Lang Shan (TLS), Ding Hu Shan (DHS), and Da Xi Shan (DXS), Guangdong, of southern China. 28 random primers generated 118 bands, out of which 26 (22.03%) were polymorphic loci, distinguishing 17 different RAPD phenotypes. Percentage of polymorphic loci, Shannon phenotypic diversity and Nei's gene diversity comprehensively indicated that JFL possessed the highest diversity, TLS and DHS in intermediate and DLS or DXS the least; the corresponding values of the population appeared correlated with the population size. Differentiation was detected among populations of A. spinulosa (1 - Hpop/Hsp = 0.7453, G(ST) = 0.7763, and phi(st) = 0.8145). AMOVA showed that 47.44% of the variance was partitioned among regions (Hainan and Guangdong), 34.01% attributed among populations within regions, whereas only 18.55% occurring within populations. Low level of intra-specific diversity was maintained in A. spinulosa with Shannon diversity and gene diversity merely 0.0560 and 0.0590, repectively. Sequence length of atp B-rbc L intergenic spacer varied from 724 bp to 730 bp. Base composition was with A+T content between 63.17% and 63.70%. 13 haplotypes of atp B-rbc L noncoding spacers were identified. UPGMA dendrogram of RAPD phenotypes, principal components analysis based on RAPD patterns, minimum spanning network and neighbour-joining (NJ) tree established on atp B-rbc L haplotypes consistently suggested the geographical subdivision of populations of A. spinulosa between Hainan and Guangdong. Breeding system and conservation strategy of A. spinulosa was discussed based on the information of population genetic structure and variation.