Polymorphisms in the ASMT and ADAMTS1 gene may increase litter size in goats

Abstract Prolificacy of most local goat breeds in China is low. Jining Grey goat is one of the most prolific goat breeds in China, it is an important goat breed for the rural economy. ASMT (acetylserotonin O‐methyltransferase) and ADAMTS1 (ADAM metallopeptidase with thrombospondin type 1 motif) are essential for animal reproduction. Single nucleotide polymorphisms (SNPs) of ASMT and ADAMTS1 genes in the highly prolific breed (Jining Grey goats), medium prolific breed (Boer goats and Guizhou White goats) and low prolific breeds (Angora goats, Liaoning Cashmere goats and Inner Mongolia Cashmere goats) were detected by polymerase chain reaction‐restriction fragment length polymorphism and sequencing. Two SNPs (g.158122T>C, g.158700G>A) of ASMT gene and two SNPs (g.7979798A>G, g.7979477C>T) of ADAMTS1 gene were identified. For g.158122T>C of ASMT gene, further analysis revealed that genotype TC or CC had 0.66 (p < 0.05) or 0.75 (p < 0.05) kids more than those with genotype TT in Jining Grey goats. No significant difference (p > 0.05) was found in litter size between TC and CC genotypes. The SNP (g.158122T>C) caused a p.Tyr298His change and this SNP mutation resulted in changes in protein binding sites and macromolecule‐binding sites. The improvement in reproductive performance may be due to changes in the structure of ASMT protein. For g.7979477C>T of ADAMTS1 gene, Jining Grey does with genotype CT or TT had 0.82 (p < 0.05) or 0.86 (p < 0.05) more kids than those with genotype CC. No significant difference (p > 0.05) was found in litter size between CT or TT genotypes. These results preliminarily indicated that C allele (g.158122T>C) of ASMT gene and T allele (g.7979477C>T) of ADAMTS1 gene are potential molecular markers which could improve litter size of Jining Grey goats and be used in goat breeding.

The biosynthesis of MLT from serotonin needs two enzymatic steps. First, serotonin N-acetyltransferase catalyses serotonin to N-acetylserotonin (NAS), and then ASMT catalyses NAS to O-methylated NAS, also known as MLT (Byeon & Back, 2016).
Melatonin is an important hormone synthesized and secreted by the pineal gland and retina in dark. It plays important roles in physiological functions including regulating the biological clock, circadian rhythms, animal reproduction and other activities (Dollins, Zhdanova, Wurtman, Lynch, & Deng, 1994). One of the important roles of MLT is to regulate the reproductive system, which can directly act on the gonads to control the secretion of reproductive hormones and also can indirectly influence the gonad to regulate the secretion of reproductive hormones through acting on the anterior pituitary firstly (Recio, Mediavilla, Cardinali, & Sanchez-Barcelo, 1994). Animals with seasonal reproductive characteristics can be divided into short-day and long-day reproductive animals.
The role of MLT on the two kinds of animals is the opposite, which can inhibit sexual activities of the long-day animals but stimulate the short-day ones.
Up to now, ASMT gene of goat had been cloned and sequenced, which is 28,511 bp (Capra hircus, NW_017189541.1). ASMT gene was expressed in bovine cumulus oocyte complexes (COCs). Melatonin receptors presented in COCs, and MLT could significantly enhance oocyte nuclear maturation and cumulus cell expansion, it indicated the potentially important roles of MLT in regulating bovine oocyte maturation (El-Raey et al., 2011). ASMT is important for animal reproduction. It was expressed in human placenta throughout pregnancy (Soliman et al., 2015). There is however paucity of information on ASMT gene in goats.
Insufficiency of ADAMTS1 expression in Sertoli cells may be related to male infertility, and it could be a potential diagnostic marker in male azoospermia (Aydos, Yukselten, Ozkavukcu, Sunguroglu, & Aydos, 2019).
The expression of ADAMTS1 gene in ovaries of biparous Mongolian sheep was about 30 fold higher than that of monotocous Mongolian sheep (He et al., 2012). ADAMTS1 mRNA was expressed in endometria, conceptus, and placentomes (Dunlap et al., 2010).
Both ASMT and ADAMTS1 are important for animal reproduction (Brown & Russell, 2014;El-Raey et al., 2011;Mishra et al., 2013). Literature on polymorphisms of caprine ASMT and ADAMTS1 genes and their association with reproductive traits is rare. Goats reared in P.R. China display different litter sizes, which provide the materials to analyse the association of ASMT and ADAMTS1 genes with prolificacy in goats. Most of the local goat breeds are low prolific in China, Jining Grey goat is one of the most prolific caprine breeds in P.R. China, and Jining Grey goat displays significant characteristics of year-round oestrus, and sexual precocity, it is an important goat breed for the rural economy. Sexual maturity of Jining Grey goats is at 3-4 months, the first mating age of that is 5-7 months. Sexual maturity of Guizhou White goats and Boer goats is 4-6 months, the first mating age of that is 6-8 months. Sexual maturity and the first mating age of Angora goats, Liaoning Cashmere goats, and Inner Mongolia Cashmere goats are the latest, about 6-8 months and 18 months respectively. The meanlitter sizes of Jining Grey goats, Guizhou White goats, Boer goats, Angora goats, Liaoning Cashmere goats, and Inner Mongolia Cashmere goats were reported to be 2.83, 2.13, 2.10, 1.31, 1.15 and 1.05 respectively (Malan, 2000;Roberts & Reeves, 1988;Tu, 1989). Single nucleotide polymorphisms (SNPs) of ASMT and ADAMTS1 genes were identified and the association of polymorphisms with litter size in goat were investigated in this study. So as to acquire molecular markers related to prolificacy for marker-assisted selection. LTD.) and then dissolved in TE buffer and stored at −20°C.

| MATERIAL S AND ME
The 296 Jining Grey does were selected at random and were the progeny of five goat bucks (n = 55, 57, 60, 61, 63). No selection on litter size or other fertility traits was conducted in the flock over previous years. Kidding seasons consisted of 3-month groups starting with March to May as season 1 (spring, n = 76), June to August as season 2 (summer, n = 68), September to November as season 3 (autumn, n = 88) and December to February as season 4 (winter, n = 64). Genomic DNA isolation method was the same with AMST gene.

| Genomic DNA isolation for ADAMTS1 gene
The 243 Jining Grey does were randomly selected from the progeny of five goat bucks (n = 42, 46, 50, 51, 54). There was no selection on litter size or other fertility traits in this population over the years.

| Cloning and sequencing
PCR products were recovered using Geneclean Ⅱ kit (Promega), and then ligated into the pGEM-T Easy vector (Promega) at 16°C overnight according to the manufacturer's instructions. After ligation, DNA was transformed into the competent cell (Escherichia coli DH5α). Positive clones were identified by the restriction enzyme and then sequenced by Sangon Biotech (Shanghai) Co. Ltd.

| Restriction fragment length polymorphism analysis
After sequence alignment, polymorphisms of ASMT and ADAMTS1 gene were screened between Liaoning Cashmere goat and Jining Grey goat. Primers P3 and P7 were used for polymerase chain reactionrestriction fragment length polymorphism (PCR-RFLP) to detect the polymorphisms for five goat breeds respectively. The mixture for PCR-RFLP of ASMT gene was: 5 U of the restriction enzyme Afa I or Msp I (Takara), 5 μl of PCR products and 1 μl 10× reaction buffer. The mixture for PCR-RFLP of ADAMTS1 gene was: 5 U of the restriction enzyme Nco Ⅰ or Hha I (Takara), 5 μl of PCR products, and 1 μl 10× reaction buffer. The mixtures were incubated at 37°C for 4 hr, and then separated on a 12% polyacrylamide gel at 120 V. After electrophoresis, the DNA fragments in the gels were visualized by silver nitrate staining, photographed and analysed using an AlphaImager ™ 2,200 and 1,220 Documentation and Analysis Systems (Alpha Innotech Corporation).

| Statistical analysis
Analysis of litter size in Jining Grey goat was performed using the following fixed effects model. Least squares mean was used for multiple comparisons in litter size among different genotypes.
where Y ijklm is the phenotypic value of litter size; µ is the population mean; S i is the fixed effect of the ith sire (i = 1, 2, 3, 4, 5); KS j is the fixed effect of the jth kidding season (j = 1, 2, 3, 4); P k is the fixed effect of the kth parity (k = 1, 2, 3); G l is the fixed effect of the lth genotype (l = three different genotypes); and e ijklm is the random residual effect of each observation. The general linear model and mean separation procedures of SAS (Ver 8.1; SAS Institute Inc.) were used to analyse the least significant differences.

| PCR amplicons of goat ASMT and ADAMTS1 genes
Using PCR with the primers P1-P3, the ASMT gene was successfully amplified (Figure 1a). ADAMTS1 gene was also successfully amplified by primers P1-P7 ( Figure 1b). The 2% agarose gels were used to separate the PCR products. It showed that the sizes of amplified fragments and the target ones were consistent and the specificity of amplification results was apparent. It could be directly analysed by RFLP and sequencing.

| RFLP for ASMT and ADAMTS1 genes
The PCR products of ASMT gene amplified by primer P2 were di-

| SNPs identified by sequencing
The PCR products with different genotypes were sequenced to confirm the mutations. The sequences of different genotypes are shown in Figure 4. One SNP g.158122T>C was found in the sequences amplified with primer P2 of the ASMT gene, which was located in exon 8 ( Figure 4a). One SNP g.158700G>A was found in the sequences amplified with primer P3 of the ASMT gene, which was located in 3′-regulatory region ( Figure 4b) For primer P3 of ADAMTS1 gene, the sequencing results of three genotypes (AA, AG and GG) revealed one mutation (g.7979798A>G; Figure 4c). For primer P7 of ADAMTS1 gene, three genotypes (CC, CT and TT) were sequenced, and one mutation (g.7979477C>T) was revealed ( Figure 4d).

| Allele and genotype frequencies of ASMT and ADAMTS1 genes in five goat breeds
Allele and genotype frequencies of the ASMT gene in five goat breeds are presented in Table 3. The results indicated that at g.158122T>C,

| Influence of different genotypes on litter size in Jining Grey goats
The least squares means and standard error for litter size of different ASMT genotypes in Jining Grey goats are presented in Table 5.

| Association between ASMT and reproductive performance
ASMT is the last step key enzyme for catalysis of 5-serotonin to In early research, ASMT mRNA was detected in the pineal gland, retina and ovary in rat (Gauer & Craft, 1996;Itoh et al., 1997), as well as oocyte, blastoderm and ovarian follicles in Japanese quail (Olszanska, Majewski, Lewczuk, & Stepinska, 2007). ASMT gene in the testes of rat expressed with 24-hr rhythms and reached the maximal values during the dark phase (Coelho et al., 2019). But the mRNA levels and enzyme activities of the ASMT and endogenous MLT displayed no significant variation throughout the oestrous cycle of rat (Skorupa et al., 2003) From sequencing data, it showed ASMT mRNA has high expression level in the ovary, testis, adrenal and mammary gland in mouse (Yue et al., 2014). In human, ASMT mRNA has the highest expression level in the ovary (Fagerberg et al., 2014). The enzyme activity of ASMT and MLT were higher in gestation-matched normotensive controls than the preeclamptic pregnancy placentas, despite insignificant expression difference for both transcript and protein of ASMT in placentas between the affected and controls (Lanoix, Guerin, & Vaillancourt, 2012). In sheep, it has a high expression level in skin, kidney, testes, and hypothalamus (Jiang et al., 2014). A higher concentration of MLT was found in ram seminal plasma than that in blood, gene expression of ASMT was high in the testis, and ASMT protein in the testis was found in the Leydig cells, spermatocytes, and spermatids (Gonzalez-Arto et al., 2016).. Furthermore, our previous study showed that ASMT mRNA expressed dramatically distinctive between Jining Grey goat and non-prolific Liaoning Cashmere goat (Huang., 2012).  Etain et al., 2012;Pagan et al., 2011). A large number of mutations of human ASMT gene have been screened to determine candidate alleles with some mental diseases, such as intellectual disability , attention-deficit/hyperactivity disorder (Chaste et al., 2011), bipolar disorder (Etain et al., 2012), autism (Jonsson et al., 2010;Wang et al., 2013). A polymorphism (rs4446909) of the promoter of the ASMT gene associated with bipolar disorder influences sleep and circadian rhythms, and it associated with lower ASMT transcription level and weaker activity in lymphoblastoid cell lines (Geoffroy et al., 2014). More than ten nonsynonymous variants of ASMT identified through the 1,000 genomes project, stopped or reduced ASMT activity in patients with neuropsychiatric disorders, and one of these variants (N17K, rs17149149) is a relatively frequent polymorphism in the Han Chinese population (Botros et al., 2013).

| Polymorphisms of ASMT gene
A restriction enzyme BglⅡ site in intron 8 of the chicken ASMT gene was found (Grechez-Cassiau, Bernard, Ladjali, Rodriguez, & F I G U R E 5 Secondary structure and tertiary structure of ASMT protein before and after the mutation at g.158122T>C based on its predicted amino acid sequence. (a) Protein secondary structure before the mutation (T Allele). (b) Secondary protein structure after the mutation (C Allele). (c) CDS sequence alignment before and after mutation. (d) The 3-dimension structure prediction before mutation (T Allele). (e) The 3-dimension structure prediction after mutation (C Allele) Voisin, 1998). Two mutations in exon 5 (G606A) and exon 6 (A737C) of sheep ASMT gene were identified, which was unassociated with seasonal reproduction and litter size (Sun., 2013). However, polymorphism of the goat ASMT gene was rare reported, the association of polymorphism with reproduction in goats was unclear.
The present study found two SNPs and one (g.158122T>C) was located in the coding region and the other (g.158700G>A) was in 3′-regulatory region of the goat ASMT gene. For g.158122T>C, the genotype distribution was different between prolific Jining Grey goat and the other four non-prolific goat breeds, and C allele frequency in Jining Grey goat was higher than other breeds with less litter size. g.158122T>C locus was moderately polymorphic (0.25 < PIC < 0.5) in Jining Grey goat, it indicated that the locus had strong selection potential in these sheep populations. And the Jining Grey goats with genotype CC and TC had more litter size than those with TT for g.158122T>C. And SNP g.158122T>C caused an amino acid change from weakly polar aromatic Tyr to polar basic His, in which amino acid residue His may participate in the enzyme activity centre. And this mutation would cause a protein binding site and five macromolecule-binding sites lost, and also obtain a new protein binding site and three new macromolecule-binding sites.
The improvement in reproductive performance may be due to changes in the structure of the ASMT protein.
So we speculated that the C allele of g.158122T>C locus may be one of the causal candidates for prolificacy in Jining Grey goat. The selection of individuals with CC would be more profitable. Further study should validate the association study, given that TT individuals were fewer than TC or CC.

| Association between ADAMTS1 and reproductive performance
As a multifunctional protease, ADAMTS1 is capable of cleaving matrix proteoglycans such as aggrecan, versican and brevican. During Drosophila oogenesis, ADAMTS1 motif A (AdamTS-A) was identified as a novel target of Janus kinase/signal transducer and activator of transcription (JAK/STAT) in epithelial follicle cells which regulates egg chamber shape by remodelling the basement membrane (Wittes & Schupbach, 2019). To the female mice which lack ADAMTS1, mature oocytes trapped in the follicles, resulted in impaired ovulation and subfertility (Mittaz et al., 2004). ADAMTS1 null ovaries had some unusual atretic follicles (Shozu et al., 2005). The morphological assessment revealed peri-ovulatory ovaries had abnormal morphogenesis (Brown et al., 2010). In ovulation, ADAMTS1 could cleave versican in the mature COC matrix. However, the ADAMTS1 null mice lost this function partially (Russell, Doyle, Ochsner, Sandy, & Richards, 2003). Ovulation rate reduced 77%, and the fertilization rate of oocytes reduced about 63% in ADAMTS1 null mice, it caused reduced litter size and littered number. Shindo et al. (2000) also found ADAMTS1 null mice had a significantly low number of pups and delivery rate, it suggested ADAMTS1 null female mice were infertile. It also reported that the number of implantation sites was decreased in ADAMTS1 null mice (Mittaz et al., 2004;Shindo et al., 2000), and loss of the mature form of ADAMTS1 caused the developmental arrest of early secondary follicles (Meng et al., 2017). The newest research also indicated ADAMTS family members play important roles in follicle rupture in cattle (Willis et al., 2017).

| Polymorphisms of ADAMTS1gene
Current studies have found that the polymorphism of the ADAMTS1 gene is associated with many diseases. Polymorphisms of the ADAMTS1 gene (rs416905 and rs402007) may be associated with ischemic stroke caused by arge artery atherosclerosis (Lyu et al., 2015). Two SNPs (rs2738, rs229038) of ADAMTS1 were significantly associated with mandibular prognathism (Guan et al., 2015;Liu et al., 2017). One SNP of ADAMTS1 (rs12140)  with genotype TT had 0.04 (p > 0.05) more kids than those with genotype CT. So we speculated that T allele of g.7979477C>T locus may be one of the causal candidates for prolificacy in Jining Grey goat. The selection of individuals with TT would be more profitable. In goats, the T allele at the g.7979477C>T locus of ADAMTS1 gene could be a potential marker for improving litter size of goat.

| CON CLUS ION
In this study, two SNPs in goat ASMT gene and two SNPs in goat ADAMTS1 gene were identified. The results indicated that C allele of the g.158122T>C locus of ASMT gene and the T allele at the g.7979477C>T locus of ADAMTS1 gene were potential molecular markers which could improve litter size of Jining Grey goats and be used in goat breeding. This novel mutation provides further evidence that ASMT and ADAMTS1 gene may play key roles in reproductive function.

| E THI C S APPROVAL AND CON S ENT TO PARTI CIPATE
All procedures involving animals were authorized and approved by the Animal Ethics Committee of the Institute of Animal Science, Chinese Academy of Agricultural Sciences with the following number: IASCAAS-AE-03.

CO N FLI C T O F I NTE R E S T
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

DATA AVA I L A B I L I T Y S TAT E M E N T
All public data generated or analyzed during this study are included in this article. Data sharing is not applicable to this article as no new data were created or analyzed in this study.