Goat RFRP gene was cloned and its mutations were detected in thirteen goat breeds whose reproductive seasonality and litter size were different. Then sequence characteristics were analyzed and association analyses were performed to reveal the relationships between mutations of RFRP gene and average daily sunshine duration, reproductive seasonality as well as litter size in goats. A 4,862 bp DNA fragment of goat RFRP gene was obtained and the complete CDS of 591 bp encodes 196 amino acids, having high homology with that of other mammals. The protein was predicted to be a secreted protein with a signal peptide of 21 amino acids. Moreover, two mutations (A712G, T1493C) in 5' regulatory region and one mutation (A3438T) in exon 2 were detected. The test of genotype distribution in six selective goat breeds showed that there was no uniform significant association between the three polymorphisms and seasonal reproduction. The association just existed in some goat breeds for each locus. Interestingly, however, there was a strong positive correlation (r = 0.830, P = 0.003) between the G allele frequency of the A712G locus and average daily sunshine duration in ten local goat breeds, suggesting that RFRP gene has undergone a selective pressure in sunshine duration and may have indirect relationship with reproductive seasonality in goats. Additionally, no significant difference was found in litter size between genotypes in prolific Jining Grey goats.