Targeted sequencing identifies a novel SH2D1A pathogenic variant in a Chinese family: Carrier screening and prenatal genetic testing

X-linked lymphoproliferative disease type 1 (XLP1) is a rare primary immunodeficiency characterized by a clinical triad consisting of severe EBV-induced hemophagocytic lymphohistiocytosis, B-cell lymphoma, and dysgammaglobulinemia. Mutations in SH2D1A gene have been revealed as the cause of XLP1. In this study, a pregnant woman with recurrence history of birthing immunodeficiency was screened for pathogenic variant because the proband sample was unavailable. We aimed to clarify the genetic diagnosis and provide prenatal testing for the family. Next-generation sequencing (NGS)-based multigene panel was used in carrier screening of the pregnant woman. Variants of immunodeficiency related genes were analyzed and prioritized. Candidate variant was verified by using Sanger sequencing. The possible influence of the identified variant was evaluated through RNA assay. Amniocentesis, karyotyping, and Sanger sequencing were performed for prenatal testing. We identified a novel de novo frameshift SH2D1A pathogenic variant (c.251_255delTTTCA) in the pregnant carrier. Peripheral blood RNA assay indicated that the mutant transcript could escape nonsense-mediated mRNA decay (NMD) and might encode a C-terminal truncated protein. Information of the variant led to success prenatal diagnosis of the fetus. In conclusion, our study clarified the genetic diagnosis and altered disease prevention for a pregnant carrier of XLP1.

The SH2D1A is extremely conserved among species and found to be highly nonpolymorphic [4,8]. To date, over 100 SH2D1A mutations have been included in HGMD (the Human Gene Mutation Database) [9]. Although it has been appreciated for nearly two decades that mutation of SH1D1A result in XLP1, the management of XLP1 is still difficult and death usually occurs within 2 months from patients presenting with EBV-induced HLH [10]. Hematopoietic stem cell transplantation has been considered as the only treatment against XLP1 [11], just like many primary immunodeficiency diseases. The mean age at death reported for individuals with SH2D1A pathogenic variant is 11 years (ranging from 2 years to 69 years) [12].
The use of next-generation sequencing (NGS) technology to move from testing small panels of genes to large multi-gene panels made its clinical application possible. Performing a larger panel but then restricting analysis to a disease-associated set of genes based on the subject's clinical phenotypes has been recommend to be more efficient in clinical diagnosis [13].
In this study, we reported the identification of a novel frameshift SH2D1A pathogenic variant in a pregnant carrier of a Chinese family. The consequence of the variant was analyzed through RNA assay. Moreover, the potential transmitting of the variant was prenatal tested by use of fetal DNA derived from amniocytes.

Subjects and ethics statement
All the subjects of this study were recruited from the outpatient department of the International Peace Maternity & Child Health Hospital (IPMCH), Shanghai Jiao Tong University School of Medicine. Peripheral blood samples were collected from all members of the pedigree (Fig 1) as well as 192 ethnically matched unrelated healthy female controls. Genomic DNA was extracted according to standard techniques. This study was prospectively reviewed and approved by the Ethics Review Committee of IPMCH of Shanghai Jiao Tong University School of Medicine, and conducted according to the Declaration of Helsinki Principles. Written informed consents were obtained from all participants, including 192 healthy female control individuals.
Targeted sequence capture and next-generation sequencing Sequencing library preparation and target capture had been described in detail previously [14]. Briefly, DNA was sheared to create fragments of 200 to 250 bp. We captured target regions of 2181 genes associated with 561 Mendelian diseases. The total size of target regions of the capture array was 6.19 Mb. Paired-end sequencing (2×90 bp) was performed with Illumina HiSeq2500 Analyzers (Illumina, San Diego, CA, USA) to generate average sequence coverage of more than 100× with more than 97% of target bases having at least 10× coverage.

Data filtering, mapping and variant detection
Sequencing output data were converted from bcl files into FastQ files by using the Illumina consensus assessment of sequence and variation (CASAVA) software. Alignment of sequence reads to the human reference genome (hg19) was performed by using Burrows-Wheeler Aligner (BWA; v.0.7.12) [15]. The Genome Analysis Toolkit (GATK; v.3.5) was applied for variants calling. SNVs and indels were annotated with ANNOVAR [16] and prioritized with inhouse developed scripts (MultiOmics One) restricted to immunodeficiency related genes

RT-PCR analysis of SH2D1A mRNA in peripheral blood mononuclear cells
Total RNA from peripheral blood mononuclear cells (PBMCs) were extracted and reversely transcribed from the carrier (II-1) and a healthy control. The RT-PCR were performed with the primers SH2D1A-RT-F: 5'-CCAGGCGTGTACTGCCTATG-3' and SH2D1A-RT-R: AGCTGAGGACTTCTTCTCAACTG. The amplified DNA fragments were directly subcloned into the TA cloning vector pGM-T (TIANGEN Biotech, Beijing, China) and sequenced by using the universal sequencing T7 primer.

Identity testing
A short tandem repeat (STR) typing human personal identification detection kit (R1004T; GENESKY, Shanghai, China) was used for identity testing in de novo variant confirmation and potential maternity contamination analysis in prenatal testing, according to the instruction manuals. The amplicons were analyzed by using Applied Biosystems 3500Dx sequencer.

Prenatal testing
Amniocentesis was performed in 18 weeks of gestation using a 22-gauge spinal needle and continuous ultrasonographic guidance [19]. An amniotic fluid sample of 30 mL was obtained. Chromosomal karyotyping and sanger sequencing were performed in parallel.

Clinical characteristics of a Chinese immunodeficiency family
A healthy 38-year-old pregnant woman of 8 weeks' gestation, with twice abnormal reproduction history, was referred to reproductive genetics clinic for counseling regarding her risk for having a child with immunodeficiency. The pedigree comprised 9 individuals spanning three generations (Fig 1). The proband boy (III-2) and his elder brother (III-1) were both diagnosed with immunodeficiency, and died with rapid disease progression from EBV infection, bronchopneumonia, severe anemia, and multiple organ dysfunction at around 2 years of age. Unfortunately, no DNA analysis had been conducted on the proband and his brother to determine which type of pathogenic gene they had. No other family member had this disorder. Family histories of immunodeficiency and abnormal reproduction were denied. The clinical features of the family supported a provisional diagnosis of immunodeficiency disease.

Identification of a novel SH2D1A pathogenic variant by NGS-based carrier screening
No proband sample was preserved, direct proband genetic diagnosis was prevented, consequently. On the basis of the clinical features of the family, NGS-based carrier screening of 2181 genes including 34 immunodeficiency-related genes was carried out with peripheral blood of the pregnant woman (II-2). Totally, 2.4 Gb raw data was obtained, with a mean coverage of 122.99× and 97.7% of bases in the defined target reached at least 10× coverage. Date priority filtering revealed that the pregnant woman was heterozygous for a novel deletion variant NM_002351.4(SH2D1A):c.
Subsequent Sanger sequencing of the variant-containing PCR amplicons was performed. Consistent with the NGS results, the variant was detected in the pregnant woman ( Fig 2B). Notably, the variant was not detected in other alive family members, including the parents of the pregnant women. Identity testing of both parental samples using short tandem repeat (STR) loci and amelogenin confirmed the de novo characteristics (Table 1), which was a strong evidence of pathogenicity according to ACMG guideline [20].

RT-PCR analysis of SH2D1A mRNA
To evaluate the influence of c.251_255delTTTCA on SH2D1A mRNA, RT-PCR analysis of total RNA extracted from peripheral blood mononuclear cells (PBMCs) of the pregnant woman (II-1) and a healthy female control was performed with the primers of SH2D1A-RT-F and SH2D1A-RT-R. The same splicing pattern was found between the wild-type control and the c.251_255delTTTCA heterozygous pregnant woman (Fig 3A). Further cDNA sequencing revealed that SH2D1A transcripts of the pregnant carrier were heterozygous for c.251_255delTTTCA (Fig 3B). RT-PCR products were then subcloned into pGM-T vector. Totally, 34 clones were sequenced using T7 universal primer. Sequence analysis showed an approximately 1:1 proportion of wild-type and c.251_255delTTTCA clones (16:18) (Fig 3C). These results suggested that the c.251_255delTTTCA mutant transcripts could escape the NMD and might encode a C-terminal truncated protein.

Variant classifying and prenatal testing
The variant can be classified into pathogenic variant according to ACMG guideline for sequence variants interpretation [20]. The family was therefore counseled that the risk for affected sons was 50%, the risk for carrier daughters was 50%, and the risk for affected daughters was low but dependent on the degree of skewing of X inactivation. Note: The analysis was performed using a human identification kit that amplified 16 short tandem repeat loci and the amelogenin locus (for sex typing). Values shown are the numbers of short tandem repeats detected at each allele of each locus. doi:10.1371/journal.pone.0172173.t001 After determining the pathogenicity of the variant, amniocentesis was performed in 18 weeks of gestation of the pregnant carrier. Both variant analysis and chromosomal karyotyping were performed with amniotic fluid fetal cells (Fig 1, III-3). STR analysis ruled out the possibility that maternal blood would contaminate the sample on testing (Table 1). Sanger sequencing and karyotyping results demonstrated that the fetus (III-3) was unaffected (NM_002351.4 (SH2D1A):c.
[251_255delTTTCA];[0]) with normal male karyotype (Fig 4). Pregnancy with unaffected fetus continued after genetic counseling. The neonate was born at 40 weeks' gestation after natural delivery with the Apgar score of 10 points.

Discussion
The success in carrier screening and prenatal diagnosis of Tay-Sachs disease (TSD; OMIM 272800) in Ashkenazi Jewish population has promoted the application of the screening panel strategy in disease prevention [25]. Prenatal genetic testing can be performed by targeted genetics testing of fetal cells, thus providing the opportunity for pregnancy termination if the pathogenic variant in the affected family has been identified.
XLP is a rare but life-threaten disorder with unknown ethnic or racial predilection. Only a few cases of XLP were reported in Chinese population [26,27]. Both SH2D1A and XIAP pathogenic variants were associated with XLP. Approximately, 83-97% XLP patients were attributed to pathogenic variants of SH2D1A gene and subclassified as XLP1; while 12% attributed to XIAP gene and subclassified as XLP2 (OMIM 300635) [28][29][30]. XLP1 is inherited in an Xlinked recessive manner. SH2D1A pathogenic variant carrier females have a half chance of transmitting the pathogenic allele in each pregnancy. In this report, the lack of affected male proband sample brought about the necessity of female carrier screening for further disease prevention. The prenatal genetic diagnostic testing usually cannot be performed until the pathogenic variant has been identified, although NGS-based prenatal screening has been reported [31,32]. Affected male proband is the most suitable for initial gene testing. For lethal or severe recessive diseases, however, proband usually dies in prenatal or neonatal periods. Thus, phenotype-based carrier screening prior to conception, likely in the form of exome or multigene panel sequencing, could give the family an opportunity for genetic diagnosis and birth defects prevention.
Apart from the risk for recurrence of XLP1 mentioned above, lymphoma was also been reported in SH2D1A female carrier due to nonrandom X-inactivation [33], suggesting that carrier females in an XLP kindred might be at the risk of developing XLP phenotype and should be monitored for early symptoms. Thus, screening SH2D1A may provide early molecular diagnosis for both XLP1 patients and the potential XLP1 carriers.
The de novo pathogenic variant of SH2D1A identified in this report has not been previously described. RNA assay indicated that the variant could escape NMD and might introduce a truncated protein. Mutations in SH2D1A are primarily loss-of-function result in diminished or absent protein expression [34]. There is no apparent genotype-phenotype correlation between SH2D1A mutations and severity of the disease. The most frequently reported mutation involved the Arginine at codon 55 (exon 2) in previous publications worldwide, as well as in China [10,26].
In summary, we reported a family that the proband died neonatally with immunodeficiency features. We identified a novel de novo frameshift SH2D1A pathogenic variant (c.251_255delTTTCA) in the pregnant carrier by using multigene panel carrier screening. Peripheral blood RNA assay indicated that the mutant transcript could escape NMD and might encode an C-terminal truncated protein. Further, information of the variant led to successful prenatal diagnosis of the fetus. Our study clarified the clinical diagnosis and altered disease prevention for the family. The novel pathogenic variant identified in this report can enrich the mutation database of SH2D1A. Moreover, our approach of multigene panel carrier screening and prenatal genetic testing would provide valuable insights in genetic disease families with inheritance risk.