NBPF9 Gene May Be Involved in Congenital Hypopituitarism: A Whole-Genome Study of a Boy with Pituitary Stalk Interruption Syndrome and His Family

Pituitary stalk interruption syndrome (PSIS) is a rare congenital defect manifesting as various degrees of anterior pituitary hormone deficiency. Scattered familial cases have been found, revealing some genetic variants. However, most of the previous research studies involved an affected sibling, and the gene spectra of the patients' entire family have rarely been reported. We conducted a study of a family consisting of a PSIS patient with his unaffected sibling and healthy parents of Han Chinese background using whole-genome sequencing. Bioinformatic analysis was carried out, and mutations related to PSIS, single-nucleotide variants (SNVs), insertion-deletion (InDELs), and structural variations (SVs) in all the four samples were filtered. After Sanger sequencing, we confirmed the variants obtained and selected three candidate genes for functional verification. The gene variations in this boy with PSIS and his lineal relatives are reported herein; de novo sequencing revealed that the NBPF9 gene may be involved in the pathogenesis of PSIS.


Introduction
Pituitary stalk interruption syndrome (PSIS) is a rare congenital defect with an estimated incidence of 0.5/100,000 births. PSIS causes various symptoms of pituitary hormone deficiency, including growth retardation and infertility [1][2][3].
Previous studies have revealed that genetic changes may contribute to its etiology. Studies on patients born to consanguineous parents have revealed germline mutations in some genes that may be responsible for the disease. For example, Fernandez-Rodriguez et al. reported a homozygous mutation, 301-302delAG, in the PROP1 gene in two consanguineous sisters [4]. Reynaud et al. identified a novel PROKR2 variant (p.Ala51 r) in two familial cases [1]. In addition, Tatsi et al. [5] found a novel heterozygous nonsense mutation (c.799C > T, p.Q267X) in the TGIF gene in a female PSIS patient with a single central incisor; her father had the same mutation but was asymptomatic.
However, the above studies only provide partial details of the genetic variants responsible, as the samples were limited to siblings or one parent. In our study, we examined DNA from a PSIS boy, his unaffected parents, and his sister. We hope that our findings help in obtaining more information about this rare disease. Moreover, we performed functional verification of a novel gene found by sequencing to explain the pathogenesis in-depth. unaffected sister and healthy parents, were also included. e study was approved by the Ethics Committee of the General Hospital of Chinese PLA, and informed consent was obtained from all participants prior to their participation in the study. e family lineage is displayed in Figure 1.

Whole-Genome Sequencing and Analysis.
Genomic DNA was extracted from peripheral blood mononuclear cells (PBMCs) from the family members in our discovery cohort. Exome capturing was performed to collect whole exons of human genomic DNA. e exon-enriched DNA libraries were sequenced in 100 bp paired-end reads using a HiSeq2000 sequencer (Illumina, San Diego, CA); coverage of 100x was typically achieved for each sample.
Qualified DNA was randomly fragmented using a Covaris Ultrasonic Processor. DNA libraries were constructed through end repair, addition of A to tails, purification, and PCR amplification. e pair-end libraries were sequenced using an Illumina HiSeq X Ten machine.

Genome Sequencing and Analysis.
After verification, the different libraries were sequenced based on the effective concentration and required raw data by Illumina HiSeq X Ten sequencing platform after being verified. Clean reads were obtained after removing adaptor sequences and low mapping quality reads. Read alignment to the reference genome hg38 was performed with Burrows-Wheeler Aligner MEM (BWA-MEM) (0.7.15). Duplicated reads were marked and removed by Picard (2.5.0). e alignment results were sorted using SAMtools (0.1.18). Recalibration of base quality scores and local realignment around InDELs were performed using Genome Analysis Toolkit (GATK) (3.6) [6]. SNVs and InDELs for each sample were determined using GATK HaplotypeCaller. Structural variants (SVs) were evaluated via LUMPY [7]. e database annotation, analysis, and global overview of variation were performed by ANNOVAR [8] and Variant Effect Predictor [9]. e genes SOX3, PROP1, LHX3, LHX4, HESX1, OTX2, PROKR2, TGIF, GPR161, and CDON have been analyzed in gene variation sets [1,2,[10][11][12][13].

Sanger Sequencing.
Candidate variants were further confirmed by Sanger sequencing executed by Taihegene Biotechnology (Beijing). Purified PCR amplification products were obtained after amplification of target gene fragments in DNA samples using a Biometra PCR instrument. Chromas 2 software was used for manual interpretation, and sequences were obtained from NCBI GenBank (http://www. ncbi.nlm.nih.gov/sites/entrez/). e first-generation sequencing results were compared with the second-generation sequencing results, and genes were screened for functional verification. e primer sequences are shown in Table 1.

Functional Verification
2.5.1. Cell Culture and Transfection. Human embryonic stem cell (HESCS) H1 cells were cultured in mTeSR1 medium in pretreated culture plates covered with Matrigel (BD Company). e cells were passaged at approximately 70% confluence. siRNAs targeting MUC12, NBPF9, and TYW1B were designed by Gemma Company. Transfection was based on optimization methods reported in domestic and foreign literature [14][15][16]. For further experiments, RNA was extracted at 24 h after transfection and protein at 48 h after transfection.

Western
Blotting. Protein samples were boiled (5 min, 95°C), separated by SDS-PAGE, and transferred onto polyvinylidene difluoride membranes. e membranes were blocked for 1 h in 5% bovine serum albumin (BSA) in TBST (20 mM Tris-HCl, pH 7.6, 150 mM NaCl, and 0.1% Tween 20). Primary antibodies against LHX3 and PITX1 (Enogene Biotech, Nanjing, China) were diluted in TBST containing 5% BSA and then applied to the membranes overnight at 4°C at a dilution of 1 : 1,000. After being washed three times with TBST, the membranes were incubated with horseradish peroxidase-conjugated secondary antibodies for 1 h. e membranes were washed three times with TBST and developed with Chemiluminescence ECL Plus-Western Blotting detection reagents (Amersham Biosciences, Piscataway, NJ, USA).

Clinical
Features. e boy with PSIS was referred to our hospital due to growth retardation. He was born in a nonconsanguineous family, and none of his relatives had a history of growth retardation or short stature. He was born via breech delivery with a body weight of 3,200 g and a height of 50 cm. After he was 5 years old, his parents noticed growth retardation compared with other children of his age. When he was then admitted to a local children's hospital, his height was 100 cm (−2.15 SDS), and he showed a delayed bone age of 3. According to GH stimulation tests, the maximum value detected at the time was 0.69 ng/mL, along with a low IGF-1 level estimated at 24.45 ng/mL. Pituitary CT revealed a hypoplastic pituitary gland. He was diagnosed with GH deficiency, and hormone replacement therapy was initiated.
In 2013, when the boy was 8 years old, he was admitted to Chinese PLA General Hospital; his height was 106 cm (−1.07 SDS), and he had a body weight of 16.5 kg. He had a single central incisor, sharp palate, and single transverse palmar crease; he could not easily bend the fourth and fifth fingers of both hands. In addition to GH deficiency, he was diagnosed with hypothyroidism (free T4: 8.35 pmol/mL (reference interval: 12-22 pmol/mL); TSH: 6.54 mU/L (reference interval: 0.35-5.5 mU/L)) and hypoadrenalism (F 8 am: 128.6 nmol/L (reference interval: 198.7-797.5 nmol/L at 24 h UFC) at admission. His bone age was that of a 5-year-old boy. e patient was treated with GH, hydrocortisone, and levothyroxine thereafter. Records from the latest follow-up in 2015 showed a height of 131 cm (−1.06 SDS) and a bone age of approximately 8 years.
His younger sister, who was the second child in this family, exhibited normal growth with no deformities. As a 5year-old girl, her height was 111.5 cm and her weight was 18 kg, which were in the normal range for her age. e results of her hormone laboratory tests were normal. e boy's parents denied consanguineous marriage or family histories of similar manifestations. e mother confirmed that she did not receive radiation or take any drugs during pregnancy. Figure 2 shows the appearance of the PSIS boy and his healthy sister.

Variant Statistics
(1) SNVs and InDELs. We obtained 4,935,685 SNVs from the four samples as well as 9.644 InDELs; based on genetics, genomic region, and functional annotation, we filtered SNV results for each sample. 48.1% were intron variants and 37.7% were intergenic variants. In exon region, 52.3% were synonymous variants.
(2) Structural Variation (SV). LUMPY and the DGV database were used to assess all SVs in the 4 samples, SVs were detected in different regions of the genomic region, and 97.06% SVs were located in intergenic region.

Analysis of Reported Genes.
Mutations in the genes SOX3, PROP1, LHX3, LHX4, HESX1, OTX2, PROKR2, TGIF, GPR161, and CDON have been reported. e current database is not related to PSIS mutation locus information matching, according to the variation function for filtering exonic regions of the mutation loci. ere are five genes that showed regional variation, including LHX4, PROP1, CDON, TGIF1, and PROKR2, as shown in Table 2.

De Novo Mutation Analysis.
We aimed to find mutations that manifested as hybrid-type genome mutation loci in the PSIS boy but not genomic mutations in his parents. A total of 19 specific mutations were found (Table 3).

Inheritance Pattern Analysis.
According to the pedigree information of the family, the variation and variation annotation information for each individual were used for grade selection; mutation loci were filtered by grade and function (nonsynonymous mutations) to obtain 267 candidate mutation loci. e 267 variation sites are located in 207 genes. Further analysis was performed by dbSNP, ClinVar, Exome Sequencing Project, 1000 Genomes Project, and Exome Aggregation Consortium database. Table 4 presents 11 filtered recessive genetic variations.    As most of the above genes have rarely been studied, we selected NBPF9, MUC12, and TYW1B as candidate genes to further analyze their roles in nervous system development and cancer [17,18].

Sanger Sequencing.
We further confirmed the candidate variations by Sanger sequencing, and the results are depicted in Figure 3.

Protein Structure Prediction.
After protein structure prediction, it was found that the NBPF9 p.L279W point mutation would cause protein structure variation. e conserved amino acid sequence is illustrated in Figure 4.
is sequence shows good conservation.

Cell Culture and Western
Blotting. Figure 5 shows the cell culture of HESCS (H1 cells). siRNAs targeting the candidate genes NBPF9, MUC12, and TYW1B were transfected into cells, and inhibition rates of more than 50% were obtained with the NBPF9 and MUC12 groups. At the protein level, LHX3 abundance decreased after transfection with NBPF9 siRNA. For PITX1, no significant change was observed between the NBPF9 transfection group and the control group.

Discussion
As a rare and complicated defect, PSIS is thought to be caused by genetic changes. However, genetic variants are found in only 5% of cases. Sporadic cases have been collected and studied, but little information was found. rough International Journal of Endocrinology familial cases, we obtained some reported mutations, but in these studies, only the father or mother's DNA was examined. Here, we evaluate a family with a PSIS boy and a healthy sibling and parents, which may reveal more details of the etiology of PSIS. We analyzed the reported genes LHX4, PROP1, CDON, TGIF, and PROKR2 in our study. LHX4/LHX3 are members of the LIM-homeobox family and regulate the proliferation and differentiation of pituitary lineages [1,2,12]. PROP1 is involved in early patterning, proliferation, positional determination, and terminal differentiation of the pituitary gland [4,13]. CDON is a component of a cell-surface receptor complex that mediates cell-cell interactions between muscle precursor cells [11]. Tatsi reported TGIF as an active transcriptional corepressor of SMAD2 that may participate in the transmission of nuclear signals during development and in adults [5]. PROKR2 is thought to be involved in endocrine angiogenesis and neuronal migration and is expressed in the hypothalamus and pituitary gland [1]. However, none of the mutation frequencies were significant. In other words, the mutation frequencies were equal to those of the healthy group. We further performed inheritance pattern analysis and found 11 genes with related variations but that have rarely been studied; thus, little is known about their functions.
In our previous study related to sporadic cases, we proposed that compromised genetic interactions among multiple signaling pathways are likely to underlie the pathogenesis of sporadic PSIS [19]. erefore, multigenetic causes may explain the pathogenesis of PSIS. erefore, in this study, we attempted to identify novel gene variations. After bioinformatics analysis, we selected candidate variants for further research.
NBPF9 belongs to the neuroblastoma breakpoint family (NBPF) family, the members of which are characterized by binding to the repetitive DUF1220 protein domain. Gene copy number variation exists in the region of human chromosome 1q21. 1. DUF1220 domains have been associated with many diseases related to development and neurogenesis, such as small head deformity, autism and schizophrenia, mental retardation, congenital heart disease, neuroblastoma, and congenital kidney and urinary tract abnormalities, and previous studies have reported that DUF1220 is associated with lung cancer [17].
We found that the level of LHX3 decreased after transfection with NBPF9 siRNA. e LHX3 and PITX1 genes appear to play important roles during pituitary development. LHX3 is a LIM homologous domain transcription factor that is one of the earliest pituitary transcription factors expressed in development; it determines the specific differentiation of tissues and cells in vertebrates and invertebrates by regulating the expression of some genes. Currently, three types of LHX3 autosomal recessive inheritance have been identified in humans, and all patients have multiple pituitary hormone deficiencies [21][22][23]. PITX1, a member of   In our study, we describe the landscape of the gene spectrum of a PSIS family and, based on further experiments, speculate that the NBPF9 gene is likely a regulatory factor affecting the development of the pituitary. Nonetheless, more samples of entire PSIS families are needed to explain the exact pathogenesis.

Conclusions
e findings of this study provide considerable information about the gene spectrum of PSIS as a rare disease. However, there is a lack of evidence regarding all of the reported genes. More samples are needed to assess clinical applicability. e NBPF9 gene may be involved in congenital hypopituitarism [20].

Data Availability
e data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
Yi-Ming Mu and Qing-Hua Guo had designed the experiments; Cheng-Zhi Wang and Ling-Ling Guo accomplished the experiments; Cheng-Zhi Wang wrote the manuscript and revised it. Qing-Hua Guo and Yi-Ming Mu contributed equally to the work.