Polymorphism and Allele Frequency Variation in Proximal End of Promoter and Exon 1 of Human Growth Hormone Gene in Pakistan- A Pilot Study

Madeeha Khalid and Amtul Jamil Sami*

Institute of Biochemistry and Biotechnology, Quaid-e-Azam Campus, University of the Punjab, 54590, Lahore

ABSTRACT

Growth hormone gene belongs to a highly polymorphic gene family whose members contain regions of very similar conserved sequences making this gene a challenging candidate to analyze. Almost 90% of the sequence of growth hormone gene family members is similar and prone to single nucleotide variations. The sequence variation occurring in the non-coding region of the GH-1 gene leads to altered transcriptional regulation causing various developmental defects in humans. In this report, we have analyzed a short sequence of ~200-240 bp spanning the proximal promoter and exon 1 (proximal portion) of the GH-1 gene to analyze the SNPs and allele frequencies of respective SNPs. Genomic DNA from blood was extracted, and the selected region was amplified using PCR. Sequencing of the amplified product was performed using Sanger sequencing, and the resulting sequence was analyzed. We have recorded 11 SNPs in total that occur in the VDRE region, 5’ UTR and untranslated region of exon 1. A greater number of the SNPs were observed to coincide with the reported data from various sources. Two SNPs at position -62 and -56 (relative to transcription start site) are novel and not reported elsewhere up to the best of our knowledge. The minor allele frequencies of SNPs at position -57, -50, -48, -40, -34, -1 and +58, respectively are higher than other populations. All the observed SNPs are located in the region involved in transcription factor binding. This study is the first report regarding the allele frequency data and SNPs of the GH promoter region and exon 1 in the Local population.

Article Information

Received 26 May 2018

Revised 23 July 2018

Accepted 10 August 2018

Available online 28 March 2019

Authors’ Contribution

AJS designed the experiments, analyzed the results and reviewed the manuscript. MK performed the experimental work and drafted the manuscript.

Key words

Polymorphism, Growth hormone gene family, Minor allele frequency, Vitamin D response elements.

DOI: http://dx.doi.org/10.17582/journal.pjz/2019.51.3.995.1001

* Corresponding author: [email protected]

0030-9923/2019/0003-0995 $ 9.00/0

Copyright 2019 Zoological Society of Pakistan

Introduction

The human growth hormone (GH) gene family is a family of homologous gene members. The gene family arose as result of gene duplication of a single ancestral gene followed by further duplication events. The GH gene family is organized into five highly similar genes namely growth hormone-1 (GH-1), chorionic somatomammotropin hormone-like-1 (CSHP-1), chorionic somatomammotropin hormone-1 (CSH-1); growth hormone-2 (GH-2) and chorionic somatomammotropin hormone-2 (CSH-2), spanning a region of 66-67kb on chromosome 17q23. The 5’ end of the GH gene cluster is occupied by the GH-1 gene which is the most comprehensively studied member of the gene family (Chen et al., 1989). GH-1 is expressed in somatotroph cells of pituitary whereas other paralogues are all placentally expressed.

GH-1 is the primary gene responsible for the secretion of human growth hormone for regulating all postnatal growth. Growth hormone controls the overall growth pattern, stature etc. in human beings. GH mainly acts by binding to its receptor and stimulating the production of insulin-like growth factor (IGF). GH along with IGF are known to operate in the pathways that regulate cell proliferation. GH-1 gene is strictly regulated to achieve the desired gene expression under various physiological states and during the different stages of life in humans. A locus control region is present at ~14.5 kb upstream of the gene that interacts with the promoter of the GH-1 to control transcription (Jones et al., 1995). Various studies have identified SNPs in the LCR region that act in conjunction with the promoter SNPs and affect the activity of the promoter (Horan et al., 2003) and thus makes it a difficult candidate to study (Kopchick, 2004).

The promoter of the GH-1 is highly polymorphic, harboring a cluster of SNPs in different haplotypic combinations (Krawczak et al., 1999). Several studies have linked the promoter single nucleotide variations with increased and decreased risk of neoplasia. In Polish and Chinese cohorts the promoter SNPs were linked to breast cancer susceptibility (Millar et al., 2003; Le Marchand et al., 2002; Wagner et al., 2005; Ren et al., 2004; Canzian et al., 2005). Association between one SNP in the promoter and Intron 4 was found to cause variation in the secretion of GH and IGF profiles and linked with accelerated bone loss and short stature in children (Hasegawa et al., 2000). Various other SNPs in the promoter region including SNPs in pituitary specific transcription factor (Pit-1), Vitamin D response elements (VDRE), nuclear factor (NF) etc. are known to cause idiopathic short stature, Isolated growth hormone deficiency, dwarfism and other growth-related changes (Madeira et al., 2016; Giordano et al., 2006; Bona et al., 2004).

In the current study, we have investigated the SNPs of the proximal promoter and exon-1 of the human growth hormone gene (GH-1) in a preliminary study in the local population to analyze the allelic variations of the SNP of GH-1 gene in Pakistan in comparison to reported data in other populations.

 

Materials and Methods

Human subjects

One hundred human subjects that fit the criteria of exclusion and inclusion were recruited for the study. The study group consisted of normal statured local Pakistani population (male and female). Height and weight of the subjects were recorded. Venous blood was collected from the subjects. Informed consent according to Helsinki’s guidelines was taken from the subjects. 3 ml of blood was collected from the subjects and added to K2EDTA vacutainers (BD®). The vacutainers containing blood were stored at 4ºC and used for isolation of genomic DNA.

The criteria for inclusion of subjects was normal stature, no genetic abnormality, no chromosomal non-disjunction cases, no cases of dwarfism or gigantism, no history of breast cancer and no history of isolated growth hormone deficiency up to two previous generations in the family respectively. The criteria of exclusion was growth abnormalities, stunted growth, abnormally low or high BMI and genetic susceptibility to cancer.

Isolation of genomic DNA and restriction digestion

Genomic DNA was isolated from the blood collected from all subjects. DNA extraction was performed as described previously (Farooq and Sami, 2013). Restriction of DNA was performed with EcoR1 to obtain a 2.6 kb fragment specific to GH-1 gene. Restriction was performed in a 50 µl reaction containing 2 µg genomic DNA, 4 µl of EcoR1 enzyme (Thermo Scientific), 4 µl of 10X EcoR1 Buffer and nuclease-free water. The reaction mixture was incubated at 37°C for 16 h. The restriction pattern was observed on a 1.2% agarose gel by electrophoresis.

Primer designing

Primers were designed against the gene sequence of human growth hormone-1 gene (Accession number NG_011676.1). Primers were designed to amplify exon 1 and the 5’ proximal promoter region of the gene for polymorphism analysis. The primers were validated by primer blast software.

Amplification and sequencing

DNA was amplified in a 25µl reaction using 50 ng of genomic DNA, 0.4 µM of each primer (Forward and Reverse), 12.5 µl GoTaq Green Master mix (Promega) and nuclease-free water, using PCR amplification conditions as described previously (Farooq and Sami, 2013). The resulting PCR products (240bp) were visualized on 1.2% agarose gel by electrophoresis. The amplified product was excised from the gel and cleaned using GeneJET PCR purification kit (ThermoFischer Scientific). The amplified product was subjected to DNA sequencing in a 10µl reaction by Sanger sequencing method using Applied Biosystems BigDye® Terminator Cycle Sequencing Kit v3.1 (Applied Biosystems, Foster City, CA, USA) per the manufacturer’s instructions. Sense and Antisense sequencing were performed for confirmation of single nucleotide changes. The sequences were assembled by ABI PRISM sequence analysis software. The final interpretation of sequence was carried out using Chromas software visually by reading the electropherograms.

Bioinformatics and statistical analysis

Individual sequences were aligned with reference sequence NG_011676.1 using multiple sequence alignment by Clustal Omega. Individual SNPs were analyzed, and genome coordinates were found using Ensemble Genome browser version 92 (https://asia.ensembl.org/index.html) from Human genome (Build GRCh38.p12). Existing Minor allele frequency data (1000 Genomes Project) of individual SNPs was downloaded and compared with the SNP data of this study.

 

Table I.- Primers for amplification of Exon 1 and promoter of GH-1.

	Sequence of primer (5’-3’)	Location*	Tm (°C)	GC (%)
Forward primer	AGCACAAGCCCGTCAGTGGC	-112 to -93	68.22	65
Reverse primer	GGTCGCTGCCTCTCCCCTCA	+111 to +130	69.50	70

*Relative to transcription start site.

 

The age and height data of the subjects included in the study were statistically analyzed including standard error and standard deviation. Frequencies of SNPs were calculated, and each SNP was tested for Hardy Weinberg Equilibrium.

 

Results and discussion

GH-1 along with other members of the paralogue is a set of highly polymorphic genes (Chen et al., 1989). The non-coding (especially the proximal promoter) and coding regions of the GH gene has a high frequency of occurrence of allelic variations and SNPs. The GH promoter contains various closely spaced cis elements for different trans-acting transcriptional regulators and are prone to sequence variations. In this study, we analyzed the proximal region of the promoter including the VDRE, 5’UTR, TATA box and exon 1 of the GH-1 gene.

PCR product of the GH-1 promoter and exon 1 was validated by electrophoresis. A 240 bp fragment was observed on the agarose gel (Farooq and Sami, 2013). The fragment was cleaned and sequenced using the same primers (forward and reverse).

The growth hormone-1 gene region under study is reported to contain over 20 single nucleotide polymorphisms in the ensemble genome viewer and the single nucleotide polymorphism database (bdSNP) of NCBI (Stenson et al., 2017; Zheng-Bradley et al., 2017; Day, 2010). 11 SNPs were recorded in the proximal end of the gene (region: vitamin D response element to intron 1). 9 of the observed SNPs are similar to the previously reported SNPs in the SNP database by different researchers. (Giordano et al., 1997; Wagner et al., 2005; Esteban et al., 2007; Adkins et al., 2005; Millar et al., 2003; Hasegawa et al., 2000; Horan et al., 2003). One novel SNP is still under investigation for a possible triallelic variation. All the allele frequencies were found to be in Hardy-Weinberg equilibrium (Table III). The SNPs tend to be clustered in the vitamin D response elements (VDRE) and near the transcription start site (Filion et al., 2006).

 

Table II.- Height and weight data of subjects.

Gender	Age (year)	Height (cm)	Weight (kg)
Male	37.5 ± 9.0	174.6 ± 7.0	73.3 ± 12.3
Female	32.5 ± 8.0	162.5 ± 5.0	55.5 ± 8.2

Values are Mean ± SD.

 

The sequence changes observed were compared with the already reported data, and the sequence variation is found to agree with the published data. However, the minor allele frequencies at various SNPs are observed to be higher in the group of the local population. SNP 1 (T→G) is present at the 5’ end of the VDRE site 1 (position -62) and is not reported elsewhere to the best of our knowledge; it is detected first time in the local population with a MAF of G=0.17 (Table III). SNP 2 at position -57, a change of T→C was observed with a frequency of 0.7 (Fig. 2). SNP 3 at position -56 contains a possible novel triallelic

 

Table III.- Single neucletoide polymorphism data with chromosomal coordinates and minor allele frequency data.

 

Table IV.- Comparison of single polymorphism with GH paralouge genes.

S. No	SNP ID	Chromos-omal location	Position on gene	GH1 Allele	Frequ-ency	GH Paralogues
						GH2	CSH1	CSH2	CSHP1
1	not reported	17:63918900	-62	T	0.83	T	T	T	T
2	rs2005172,	17:63918895	-57	G	0.227	A	T	T	G
	CR031008			C	0.773
3	not reported	17:63918894	-56			G	C	C	C
4	not reported	17:63918892	-54	No data		A	A	A	A
5	rs788661715	17:63918888	-50	G	0.852	C	C	C	C
				C	0.148
6	rs61762493	17:63918886	-48	G	0.67	A	A	A	A
				A	0.33
7	rs71651678	17:63918878	-40	G	0.5	G	C	C	G
				C	0.5
8	rs759606522	17:63918872	-34	A	0.886	A	A	A	A
				C	0.113
9	rs6171,	17:63918844	-6	G	0.307	A	G	G	A
	CR031006			A	0.693
10	Rs695	17:63918839	-1	A	0.625	A	T	T	C
				T	0.329
11	Rs771309275	17:63918781	58	C	0.77	C	T	C	C
				T	0.23

 

variation and could not be confirmed by conventional sequencing methods, so it is under further study and analysis. SNP 4 at position -54 is a transition SNP A→G with a frequency of ~1%. SNPs 2, 5-8, 10 and 11 are similar to previously reported SNPs, however, the frequencies of minor alleles are higher as compared to Polish, German, Japanese and Caucasian populations (Table III) (Wagner et al., 1997, 2005).

Another critical SNP at position -34 (A→C) lies in the VDRE site2 of the gene with an MAF of 0.11 for C allele. At all instances the C allele was present in the heterozygous form. The SNPs 1-8 all lie in the vitamin D response element and should be considered significant as they can cause a variation in the binding of trans-element (Vitamin D-Vitamin D receptor complex) of the VDRE thus causing a change in transcriptional regulation of GH-1 (Seoane et al., 2002). The Human GH VDRE is an imperfect direct repeat, and the Glu42 of the VDR interacts with the -57 T with higher affinity to exert a negative regulatory effect. The change at -57 from T to any other base affect the binding of cis-trans elements creating a change in the expression profile of the GH-1 gene (Shaffer and Gewrith, 2002).

The SNPs 10 and 11 at positions -6 and -1 respectively are present near the start of transcription start site. SNP 10 is at -1 position, A to T change was recorded with a frequency of 0.329 for T allele, and the rate of heterozygosity was 100%. SNP at position -6 that is present in the 5’ Untranslated Region (UTR) and has a similar MAF as reported by 1000 genomes project (Sayers et al., 2011).

The portion of the GH-1 gene analyzed for SNPs was aligned with the paralogues of the GH gene family and most of the variations sites that exist in the GH1 sequence differ from that of the paralogue genes.

The alignment of the SNPs with the paralogues indicated that 4 out of 11 SNPs (~36%) in the GH-1 gene does not correspond to any of the other four genes (Table IV) other SNPs match with one of the four paralogues, respectively. The possible mechanism of the sequence recombination and high variation rate is due to possible gene conversion or transfer where the other four paralogue genes act as donors. Our results confirm the existing data of the GH gene polymorphism and highlight the variations in allele frequency data of the local population from other studied populations including German, Polish, Japanese, Chinese, and Caucasian etc. (Esteban et al., 2007; Farooq and Sami, 2013).

SNP 1 and 3 have not been reported in any population till date, all other SNPs are published by various researchers in different populations. Some of the SNPs of the GH-1 promoter and exon 1 region are observed and reported in short stature subjects in multiple studies furthermore the haplotype map of the patients with Idiopathic short stature, and isolated growth hormone deficiency varies from that of the normal statured population (Hirschhorn et al., 2001; Dennison et al., 2004).

 

Conclusion

GH-1 promoter is susceptible to sequence variations and polymorphisms. There exists a high density of closely spaced and clustered single nucleotide polymorphisms in the proximal region of promoter including transcriptional regulator binding elements and the 5’UTR. Along with the promoter, the exon 1 (proximal portion) including the signal peptide also contains sequence changes and are frequently reported till date. The comparison of the observed SNP in the current study indicates same changes as given in the previous work. Some base changes are also comparable to the SNPs found in idiopathic short stature and isolated growth hormone deficiency patients. In future the results reported could be helpful to design mutant variants to perform functional studies and analyze the effect of SNPs on the expression and functionality of the gene.

 

Acknowledgement

Authors are grateful to Higher Education Commission for providing indigenous PhD. Fellowship to the first author to carry out this research work.

 

Statement of conflict of interest

Authors declare no conflict of interest.

 

References

Adkins, R.M., Campese, C., Vaidya, R. and Boyd, T.K., 2005. Association between fetal growth restriction and polymorphisms at sites-1 and +3 of pituitary growth hormone: A case-control study. BMC Pregn. Childb., 5: 2. https://doi.org/10.1186/1471-2393-5-2

Bona, G., Paracchini, R., Giordano, M. and Momigliano-Richiardi, P., 2004. Genetic defects in GH synthesis and secretion. Eur. J. Endocrinol., 151(Suppl 1): S3-S9. https://doi.org/10.1530/eje.0.151S003

Canzian, F., McKay, J.D., Cleveland, R.J., Dossus, L., Biessy, C., Boillot, C., Rinaldi, S., Llewellyn, M., Chajès, V., Clavel-Chapelon, F. and Téhard, B., 2005. Genetic variation in the growth hormone synthesis pathway in relation to circulating insulin-like growth factor-I, insulin-like growth factor binding protein-3, and breast cancer risk: results from the European prospective investigation into cancer and nutrition study. Cancer Epidemiol. Prevent. Biomark., 14: 2316-2325. https://doi.org/10.1158/1055-9965.EPI-04-0874

Chen, E.Y., Liao, Y.C., Smith, D.H., Barrera-Saldaña, H.A., Gelinas, R.E. and Seeburg, P.H., 1989. The human growth hormone locus: nucleotide sequence, biology, and evolution. Genomics, 4: 479-497. https://doi.org/10.1016/0888-7543(89)90271-1

Day, I.N., 2010. dbSNP in the detail and copy number complexities. Hum. Mut., 31: 2-4. https://doi.org/10.1002/humu.21149

Dennison, E.M., Syddall, H.E., Rodriguez, S., Voropanov, A., Day, I.N.M., Cooper, C. and Southampton Genetic Epidemiology Research Group, 2004. Polymorphism in the growth hormone gene, weight in infancy, and adult bone mass. J. clin. Endocrinol. Metab., 89: 4898-4903. https://doi.org/10.1210/jc.2004-0151

Esteban, C., Audí, L., Carrascosa, A., Fernández-Cancio, M., Pérez-Arroyo, A., Ulied, A., Andaluz, P., Arjona, R., Albisu, M., Clemente, M. and Gussinyé, M., 2007. Human growth hormone (GH1) gene polymorphism map in a normal-statured adult population. Clin. Endocrinol., 66: 258-268. https://doi.org/10.1111/j.1365-2265.2006.02718.x

Farooq, A. and Jamil-Sami, A., 2013. Targeting of GH gene at the proximal end for identification of markers for breast cancer among Pakistani women. Pakistan J. Zool., 45: 569-573.

Filion, G.J., Zhenilo, S., Salozhin, S., Yamada, D., Prokhortchouk, E. and Defossez, P.A., 2006. A family of human zinc finger proteins that bind methylated DNA and repress transcription. Mol. Cell. Biol., 26: 169-181. https://doi.org/10.1128/MCB.26.1.169-181.2006

Giordano, M., Marchetti, C., Chiorboli, E., Bona, G. and Richiardi, P.M., 1997. Evidence for gene conversion in the generation of extensive polymorphism in the promoter of the growth hormone gene. Hum. Genet., 100: 249-255. https://doi.org/10.1007/s004390050500

Giordano, M., Godi, M., Giacopelli, F., Lessi, M., Mellone, S., Paracchini, R., Petri, A., Bellone, J., Ravazzolo, R., Bona, G. and Momigliano-Richiardi, P., 2006. A variation in a Pit-1 site in the growth hormone gene (GH1) promoter induces a differential transcriptional activity. Mol. Cell. Endocrinol., 249: 51-57. https://doi.org/10.1016/j.mce.2006.01.006

Jones, B.K., Monks, B.R., Liebhaber, S.A. and Cooke, N.E., 1995. The human growth hormone gene is regulated by a multicomponent locus control region. Mol. Cell. Biol., 15: 7010-7021.

Hasegawa, Y., Fujii, K., Yamada, M., Igarashi, Y., Tachibana, K., Tanaka, T., Onigata, K., Nishi, Y., Kato, S. and Hasegawa, T., 2000. Identification of novel human GH-1 gene polymorphisms that are associated with growth hormone secretion and height. J. clin. Endocrinol. Metab., 85: 1290-1295. https://doi.org/10.1210/jc.85.3.1290

Hirschhorn, J.N., Lindgren, C.M., Daly, M.J., Kirby, A., Schaffner, S.F., Burtt, N.P., Altshuler, D., Parker, A., Rioux, J.D., Platko, J. and Gaudet, D., 2001. Genomewide linkage analysis of stature in multiple populations reveals several regions with evidence of linkage to adult height. Am. J. Hum. Genet., 69: 106-116. https://doi.org/10.1086/321287

Horan, M., Millar, D.S., Hedderich, J., Lewis, G., Newsway, V., Mo, N., Fryklund, L., Procter, A.M., Krawczak, M. and Cooper, D.N., 2003. Human growth hormone 1 (GH1) gene expression: Complex haplotype-dependent influence of polymorphic variation in the proximal promoter and locus control region. Hum. Mut., 21: 408-423. https://doi.org/10.1002/humu.10167

International HapMap Consortium, 2003. The international HapMap project. Nature, 426: 789. https://doi.org/10.1038/nature02168

Kopchick, J.J., 2004. History and future of growth hormone research. Horm. Res. Paediatr., 60: 103-112. https://doi.org/10.1159/000074510

Krawczak, M., Chuzhanova, N.A. and Cooper, D.N., 1999. Evolution of the proximal promoter region of the mammalian growth hormone gene. Gene, 237: 143-151. https://doi.org/10.1016/S0378-1119(99)00313-3

Le Marchand, L., Donlon, T., Seifried, A., Kaaks, R., Rinaldi, S. and Wilkens, L.R., 2002. Association of a common polymorphism in the human GH1 gene with colorectal neoplasia. J. Nat. Cancer Inst., 94: 454-460. https://doi.org/10.1093/jnci/94.6.454

Madeira, J. L., Jorge, A. A., Martin, R. M., Montenegro, L.R., Franca, M.M., Costalonga, E.F., Correa, F.A., Otto, A.P., Arnhold, I.J., Freitas, H.S. and Machado, U.F., 2016. A homozygous point mutation in the GH1 promoter (c.-223C> T) leads to reduced GH1 expression in siblings with isolated GH deficiency (IGHD). Eur. J. Endocrinol., 175: K7-K15. https://doi.org/10.1530/EJE-15-0149

Millar, D.S., Lewis, M.D., Horan, M., Newsway, V., Easter, T.E., Gregory, J.W., Fryklund, L., Norin, M., Crowne, E.C., Davies, S.J. and Edwards, P., 2003. Novel mutations of the growth hormone 1 (GH1) gene disclosed by modulation of the clinical selection criteria for individuals with short stature. Hum. Mut., 21: 424-440. https://doi.org/10.1002/humu.10168

Ren, Z., Cai, Q., Shu, X.O., Cai, H., Cheng, J.R., Wen, W. Q., Gao, Y.T. and Zheng, W., 2004. Genetic polymorphisms in the human growth hormone-1 gene (GH1) and the risk of breast carcinoma. Cancer, 101: 251-257. https://doi.org/10.1002/cncr.20375

Sayers, E.W., Barrett, T., Benson, D.A., Bolton, E., Bryant, S.H., Canese, K., Chetvernin, V., Church, D.M., DiCuccio, M., Federhen, S. and Feolo, M., 2011. Database resources of the national center for biotechnology information. Nucl. Acids Res., 39: D38. https://doi.org/10.1093/nar/gkq1172

Seoane, S., Alonso, M., Segura, C. and Pérez-Fernández, R., 2002. Localization of a negative vitamin D response sequence in the human growth hormone gene. Biochem. Biophy. Res. Comm., 292: 250-255. https://doi.org/10.1006/bbrc.2002.6641

Shaffer, P.L. and Gewirth, D.T., 2002. Structural basis of VDR–DNA interactions on direct repeat response elements. EMBO J., 21: 2242-2252. https://doi.org/10.1093/emboj/21.9.2242

Stenson, P.D., Mort, M., Ball, E.V., Evans, K., Hayden, M., Heywood, S., Hussain, M., Phillips, A.D. and Cooper, D.N., 2017. The human gene mutation database: Towards a comprehensive repository of inherited mutation data for medical research, genetic diagnosis and next-generation sequencing studies. Hum. Genet., 136: 665-677. https://doi.org/10.1007/s00439-017-1779-6

Wagner, J.K., Eblé, A., Cogan, J.D., Prince, M.A., Phillips, J. and Mullis, P.E., 1997. Allelic variations in the human growth hormone-1 gene promoter of growth hormone-deficient patients and normal controls. Eur. J. Endocrinol., 137: 474-481. https://doi.org/10.1530/eje.0.1370474

Wagner, K., Hemminki, K., Israelsson, E., Grzybowska, E., Klaes, R., Chen, B., Butkiewicz, D., Pamula, J., Pekala, W. and Försti, A., 2005. Association of polymorphisms and haplotypes in the human growth hormone 1 (GH1) gene with breast cancer. Endoc-Relat. Cancer, 12: 917-928.

Zheng-Bradley, X., Streeter, I., Fairley, S., Richardson, D., Clarke, L., Flicek, P. and 1000 Genomes Project Consortium, 2017. Alignment of 1000 genomes project reads to reference assembly GRCh38. Gigascience, 6: 1-8. https://doi.org/10.1093/gigascience/gix038