Research Article | Open Access
Common Variable Immunodeficiency with Genetic Defects Identified by Whole Exome Sequencing
Common variable immunodeficiency (CVID) belongs to the primary immunodeficiency disorders (PIDs), presenting a profound heterogeneity in phenotype and genotype, with monogenic or complex causes. Recurrent respiratory infections are the most common clinical manifestations. CVID patients can also develop various autoimmune and lymphoproliferative complications. Genetic testing such as whole exome sequencing (WES) can be utilized to investigate likely genetic defects, helping for better clinical management. We described the clinical phenotypes of three sporadic cases of CVID, who developed recurrent respiratory infections with different autoimmune and lymphoproliferative complications. WES was applied to screen disease-causing or disease-associated mutations. Two patients were identified to have monogenic disorders, with compound heterozygous mutations in LRBA for one patient and a frameshift insertion in NFKB1 for another. The third patient was identified to be a complex form of CVID. Two novel mutations were identified, respectively, in LRBA and NFKB1. A combination of clinical and genetic diagnosis can be more extensively utilized in the clinical practice due to the complexity and heterogeneity of CVID.
Common variable immunodeficiency (CVID) is the most common primary immunodeficiency disorder (PID) with an estimated prevalence of 1:50,000 to 1:25,000 . CVID is characterized by remarkable hypogammaglobulinemia resulting from B-lymphocyte dysfunction, which also involves T-lymphocyte abnormalities . CVID is usually complicated by recurrent infections, autoimmune diseases, malignancies, and lymphoproliferative diseases. Pulmonary infections are most commonly developed in CVID patients, which may further lead to bronchiectasis. Noninfectious pulmonary complications exhibit as granulomatous-lymphocytic interstitial lung disease (GLILD). Multiple systems can be involved as a result of immune dysregulation, such as cytopenia, splenomegaly, enteropathy, and various autoimmune diseases.
CVID shows a considerable heterogeneity in phenotype and genotype. Most cases of CVID occur sporadically, with approximately 5%~25% of patients having a familial tendency however . Most cases have an unknown genetic cause, with monogenic diseases accounting for only 2%~10% with both autosomal recessive and dominant inherited mutations . Monogenic defects that have been implicated include recessively inherited mutations in ICOS, CD19, CD20, CD21, CD27, CD81, IL21, IL21R, LRBA, PRKCD, and RAC2, and dominantly inherited mutations in TNFSF12, CTLA4, PLCG2, NFKB1, NFKB2, PIK3CD, PIK3R1, VAV1, BLK, IKZF1, and IRF2BP2, as well as monoallelic or biallelic mutations in TNFRSF13B and TNFRSF13C . A subgroup of patients demonstrates a complex disorder, rather than a Mendelian inheritance, with possible polygenic, epigenetic, and/or environmental factors participating in CVID pathogenesis.
Genetic testing is applied to unravel the heterogeneity of CVID, including genetic linkage studies in families, genomewide association study (GWAS), whole genome sequencing (WGS), whole exome sequencing (WES), RNA sequencing, and epigenetic studies [5, 6]. CVID patients may benefit from genetic testing to identify harmful variations and allow genetic counseling. WES is a potential effective tool to help discover genetic defects for monogenic cases, which targets protein-coding sequences accounting for 1% of the whole genome, but reported to harbor 85% of disease-causing variants . Here we describe the phenotypes of three sporadic CVID cases admitted to the Department of Respiratory and Critical Care Medicine, Peking University People’s Hospital, and further investigate the genetic disorders using the approach of WES.
2. Materials and Methods
2.1. Patient Phenotypes
ase 1, female, 46 years old, presented with a history of recurrent lower respiratory tract infections at the age of 38. She developed suppurative meningitis at the age of 43, with a sequela of blindness. Other clinical conditions included chronic sinusitis, pancytopenia, splenomegaly, and sensorineural hearing loss. Immunological findings showed decreased levels of IgG, IgA, and IgM in serum and a low proportion of B cells (Table 1). GLILD was suspected according to the chest CT without a pathologic confirmation (Figure 1(a)). Pulmonary function tests (PFTs) demonstrated a mild restrictive ventilatory defect and a diffusion impairment. Abdominal contrast-enhanced CT revealed multiple hypodense lesions in the spleen, which mimicked splenic infarction (Figure 1(b)).
ase 2, female, 54 years old, presented with a history of recurrent lower respiratory tract infections at the age of 39. She suffered from tuberculous pleuritis at the age of 46. Other clinical conditions included chronic sinusitis, intermittent gastrointestinal infections, neutropenia, sensorineural hearing loss, and splenomegaly. The patient had extremely low levels of IgG, IgA, and IgM, a low proportion of B cells, and an inverted CD4+/CD8+ ratio (Table 1). Chest CT showed bilateral bronchiectasis with multiple infiltrates (Figure 1(c)). PFTs demonstrated a severe obstructive ventilatory defect and a diffusion impairment.
ase 3, female, 34 years old, presented as recurrent lower respiratory tract infections with an onset age of seven. Other clinical conditions included autoimmune hemolytic anemia, splenomegaly, hypothyroidism, and nephrotic syndrome. Decreased levels of IgG, IgA, and IgM and an inverted CD4+/CD8+ ratio were also detected (Table 1). Chest CT demonstrated diffused nodules, bronchiectasis, and mediastinal lymphadenopathy (Figure 1(d)), with a severe restrictive ventilatory defect and a diffusion impairment confirmed by PFTs. Wedge resection of the right middle lobe and right lower lobe was performed; the diagnosis of GLILD was confirmed by pathology subsequently (Figures 1(e) and 1(f)).
The three cases met the criteria for CVID established by European Society for Immunodeficiencies/Pan-American Group for Immunodeficiency . All the cases received antibiotics and immunoglobulin replacement therapy and survived to date. ase 3 received extra corticosteroids treatment for GLILD with clinical improvement. Normal IgG, IgA, and IgM levels were detected in the offspring of the probands, including the daughter and the son of ase 1, the son of ase 2, and the son of ase 3.
2.2. Whole Exome Sequencing
Genomic DNA of the patients and offspring was isolated from peripheral blood mononuclear cells. The protocol was approved by the Ethics Committee of Peking University People’s Hospital, and informed consent was obtained from the subjects. The exome sequences were efficiently enriched from 1.0 μg genomic DNA using the Agilent SureSelect Human All Exon V6 kit (Agilent Technologies). Qualified DNA was randomly fragmented to an average size of 180~280 bp; then DNA fragment was end-repaired and phosphorylated, followed by A-tailing and ligation at the 3’ ends with paired-end adaptors (Illumina). At last, DNA library was sequenced on a Hiseq 4000 (Illumina) for paired-end 150 bp reads. Valid sequencing data were mapped to the reference genome (UCSC hg19) by Burrows-Wheeler Aligner software (BWA) . SAMtools  and Picard Tools (http://broadinstitute.github.io/picard/) were utilized to sort the results and mark duplicate reads, respectively. Single nucleotide variants (SNVs) and insertions and deletions (Indels) detected were annotated with the ANNOVAR software (http://annovar.openbioinformatics.org/en/latest/) .
2.3. Candidate Gene Screening
Patients’ exomes were filtered for mutations associated with immunodeficiency . Synonymous SNVs were discarded. SNVs obtained with the minor allele frequency (MAF)<0.01 in the general population according to the Exome Aggregation Consortium database (ExAC, Broad Institute) were supposed to be novel or rare, which were potentially significant. Deleterious variations were predicted utilizing SIFT, PolyPhen-2, MutationTaster, and CADD [13–16]. Afterwards, candidate variants were screened based on the phenotypes and the inheritance patterns of the patients. Deleterious indels associated with the phenotypes were also screened. As the three cases were all sporadic without a familial inheritance tendency, we firstly hypothesized that the patients had monogenic disorders with an autosomal recessive pattern caused by a homozygous or compound heterozygous inheritance or with a dominant pattern with an incomplete penetrance. If no causative mutations were found, we considered the case as a complex form of CVID rather than a Mendelian disease. Top likely disease-associated mutations were confirmed by Sanger sequencing.
The filtering results for candidate SNVs were shown in Figure 2. Three SNVs in case 1, one SNV in case 2, one SNV, and one insertion in case 3 were most likely to be disease-associated (Table 2). All the candidate mutations have been confirmed by Sanger sequencing (supplementary figure (available here)).
sing SIFT, PolyPhen, MutationTaster, and CADD to predict deleterious SNVs. SIFT (T, tolerated; D, deleterious); PolyPhen (D, probably damaging; P, possibly damaging; B, benign); MutationTaster (D, disease-causing); CADD (score>15 implied deleterious variations).
In case 1, heterozygous variants in LRBA (c.8436G>C and c.4089A>T) and in TNFRSF13B (c.226G>A) were identified to be disease-associated. LRBA deficiency was reported to be recessively inherited accounting for 26.74% of monogenic causes of CVID [4, 12, 17]. Further sequencing of her offspring showed biallelic variants in LRBA of the proband were separately inherited by her unaffected daughter (c.8436G>C) and her unaffected son (c.4089A>T), demonstrating the LRBA mutations were probably compound heterozygous which led to CVID. The variant of c.4089A>T is novel and has not been reported before. Biallelic and monoallelic TNFRSF13B variants are both reported to be disease-associated [4, 18]. However, the TNFRSF13B variant in this patient alone is insufficient to cause a CVID phenotype, which may play a role in CVID development but was not crucial. We also found that the TNFRSF13B variant was inherited by her unaffected daughter. Therefore, we identified compound heterozygous mutations in LABA most likely to be disease-causing in this patient.
In ase 2, a heterozygous LRBA variant (c.3764G>C) was filtered out. However, a monoallelic LRBA mutation cannot explain the phenotype of the patient. As a result, we assumed ase 2 was a complex form of CVID rather than a monogenic disease, although the monoallelic variant in LRBA may play a minor role. Further investigations found that the LRBA variant was not transmitted to her unaffected son.
In case 3, a heterozygous LRBA variant (c.5084T>C) was identified, which was inherited by her unaffected son. Likewise, this monoallelic LRBA mutation was insufficient to explain a disease-causing effect. We also found a monoallelic insertion in NFKB1 (c.666dupG), resulting in a frameshift mutation. Considering mutations in NFKB1 were reported to be inherited with an autosomal dominant trait , the insertion in NFKB1 was considered to be causative in case 3. This NFKB1 insertion is also novel and has not been reported before. Since mutations in NFKB1 were inherited with an incomplete penetrance , her son also inherited the mutation without a clinical phenotype of CVID.
Overall, two patients with monogenic causes and one patient with a complex cause of CVID were identified using WES. Two novel mutations were found, respectively, in LRBA (c.4089A>T) and in NFKB1 (c.666dupG).
The three cases of CVID we discussed had different clinical phenotypes with different genetic defects. Respiratory tract infections and various noninfectious complications were noted, with GLILD suspected or confirmed in two cases. The radiological abnormalities in CVID patients with GLILD mainly manifest as diffuse reticulation or nodules, while isolated bronchiectasis can also be seen . Lymphadenopathy is sometimes accompanied by interstitial changes, and lymphoma should be excluded, with a higher risk in CVID patients . The type of impaired lung function was consistent with the radiological findings in our cases, with interstitial lung disease (ILD) outweighing the impact of bronchiectasis in case 3. GLILD exhibits both granulomatous and lymphoproliferative patterns in histologic examinations, consisting of lymphocytic interstitial pneumonia, follicular bronchiolitis, and lymphoid hyperplasia, with a potential role of human herpesvirus 8 in the pathogenesis . At present, GLILD is defined as a distinct clinico-radio-pathological ILD in CVID patients, with corticosteroids as the first-line treatment . Splenomegaly in our patients was also a manifestation of the lymphoproliferative and granulomatous disease. Hypodense splenic lesions were considered nonspecific with a nonneoplastic etiology, since spontaneous resolution was reported for similar lesions . Recognition of the noninfectious complications, which is not completely understood currently, is important for including an accurate evaluation for the disease, as noninfectious complications reduce the overall survival in CVID patients [24, 25].
As an increased number of genetic defects was found, genetic diagnosis is more and more important in PID classification, because of the broad overlapping in clinical and immunological features. With the help of genetic testing, the diagnosis of CVID could be more precise, such as LRBA deficiency or NFκB1 deficiency. WES was reported to identify 30% of disease-causing mutations in CVID patients with severe phenotypes .
Deleterious LRBA mutations were found in all three patients. Lipopolysaccharide-responsive beige-like anchor protein (LRBA) is a cytosolic protein expressed by immune effector cells. LRBA participates in the CTLA-4 pathways, which negatively regulates immune responses . LRBA deficiency results in a loss of CTLA-4 protein. Hence, LRBA deficiency usually results in immune dysregulation and autoimmunity in CVID patients. LRBA mutations were also associated with inflammatory bowel disease-like disorder, and immune dysregulation, polyendocrinopathy, enteropathy, and X-linked-like disease [28, 29]. Homozygous mutations in LRBA were reported to result in loss of function in multiple consanguineous families , while compound heterozygous mutations can also cause a CVID phenotype as we identified in case 1 . Monoallelic LRBA mutations are insufficient to be disease-causing and reported to be disease-associated, causing recurrent pulmonary infections, organomegaly, and autoimmune cytopenia . An explanation is that a monoallelic mutation may also influence the protein stability of LRBA .
TNFRSF13B belongs to the tumor necrosis factor receptor (TNFR) superfamily, which encodes transmembrane activator and calcium-modulator and cyclophilin ligand interactor (TACI) and plays a vital role in the maturation and survival of peripheral B cells . Heterozygous variations in TNFRSF13B are disease-modifying mutations rather than disease-causing mutations, which may increase the risk for developing CVID and are also found in healthy individuals [4, 18, 32]. The TNFRSF13B mutations can participate in the pathogenesis of CVID through the epistatic interactions with mutations of other genes . In addition, the TNFRSF13B mutations were reported to have an incomplete penetrance, which may explain the different phenotypes between family members harboring the same variants [4, 31, 34, 35].
NF-κB signaling pathways participate in the process of B cell differentiation and function, playing a pivotal role in the pathogenesis of multiple diseases of immune dysfunction [36, 37]. NFKB1 encodes the mature p52 subunit and its precursor p105 subunit, which belongs to the NF-κB transcription factor family, reported to be associated with CVID in multiple consanguine families or sporadic cases [26, 37]. NFKB1 mutations have an autosomal dominant inheritance with an incomplete penetrance, since the variants were also found in unaffected family members .
Our study has some limitations. Firstly, the functional effects of the found mutations are not verified at the expression level but have provided additional information for clinical management. Secondly, since the parents of the probands were unavailable for genetic testing, the interpretation for the mode of inheritance and gene functions is limited.
In conclusion, we described the clinical phenotypes of three sporadic cases of CVID and tried to identify genetic defects in these patients using WES. Since CVID forms a heterogeneous group of phenotypes and genotypes, genetic testing promotes the diagnosis of CVID to the genetic level, as well as profoundly improves our understanding for CVID. A combination of clinical and genetic diagnosis can be more extensively utilized in the clinical practice of CVID. However, because of the complexity of the disease, genetic investigation is still a great challenge.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Ran Li and Yali Zheng contributed equally to this work.
Supplementary figure. Results of Sanger sequencing for mutations likely associated with the disease. A. LRBA variant (c.8436G>C) in case 1; B. LRBA variant (c.4089A>T) in case 1; C. TNFRSF13B variant (c.226G>A) in case 1; D. LRBA variant (c.3764G>C) in case 2; E. LRBA variant (c.5084T>C) in case 3; F. NFKB1 variant in case 3 (c.666dupG). (Supplementary Materials)
- B. Gathmann, N. Mahlaoui, L. Gérard et al., “Clinical picture and treatment of 2212 patients with common variable immunodeficiency,” The Journal of Allergy and Clinical Immunology, vol. 134, pp. 116–126, 2014.
- A. J. Stagg, M. Funauchi, S. C. Knight, A. D. B. Webster, and J. Farrant, “Failure in antigen responses by T cells from patients with common variable immunodeficiency (CVID),” Clinical & Experimental Immunology, vol. 96, no. 1, pp. 48–53, 1994.
- F. A. Bonilla, I. Barlan, H. Chapel et al., “International Consensus Document (ICON): common variable immunodeficiency disorders,” Journal of Allergy and Clinical Immunology: In Practice, vol. 4, no. 1, pp. 38–59, 2016.
- D. J. A. Bogaert, M. Dullaers, B. N. Lambrecht, K. Y. Vermaelen, E. De Baere, and F. Haerynck, “Genes associated with common variable immunodeficiency: one diagnosis to rule them all?” Journal of Medical Genetics, vol. 53, no. 9, pp. 575–590, 2016.
- A. Kienzler, C. E. Hargreaves, and S. Y. Patel, “The role of genomics in common variable immunodeficiency disorders,” Clinical & Experimental Immunology, vol. 188, no. 3, pp. 326–332, 2017.
- G. de Valles-Ibáñez, A. Esteve-Solé, M. Piquer et al., “Evaluating the genetics of common variable immunodeficiency: monogenetic model and beyond,” Frontiers in Immunology, vol. 9, article 636, 2018.
- E. Hodges, Z. Xuan, V. Balija et al., “Genome-wide in situ exon capture for selective resequencing,” Nature Genetics, vol. 39, no. 12, pp. 1522–1527, 2007.
- M. E. Conley, L. D. Notarangelo, and A. Etzioni, “Diagnostic criteria for primary immunodeficiencies. Representing PAGID (Pan-American Group for Immunodeficiency) and ESID (European Society for Immunodeficiencies),” Clinical Immunology, vol. 93, no. 3, pp. 190–197, 1999.
- H. Li and R. Durbin, “Fast and accurate long-read alignment with Burrows-Wheeler transform,” Bioinformatics, vol. 26, no. 5, pp. 589–595, 2010.
- H. Li, B. Handsaker, A. Wysoker et al., “The sequence alignment/map format and SAMtools,” Bioinformatics, vol. 25, no. 16, pp. 2078-2079, 2009.
- K. Wang, M. Li, and H. Hakonarson, “ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data,” Nucleic Acids Research, vol. 38, no. 16, article e164, 2010.
- C. Picard, W. Al-Herz, A. Bousfiha et al., “Primary immunodeficiency diseases: an update on the classification from the international union of Immunological societies expert committee for primary immunodeficiency 2015,” Journal of Clinical Immunology, vol. 35, no. 8, pp. 696–726, 2015.
- P. C. Ng and S. Henikoff, “SIFT: predicting amino acid changes that affect protein function,” Nucleic Acids Research, vol. 31, no. 13, pp. 3812–3814, 2003.
- I. Adzhubei, D. M. Jordan, and S. R. Sunyaev, “Predicting functional effect of human missense mutations using PolyPhen-2,” Current Protocols in Human Genetics, no. 76, article no. 7.20, 2013.
- J. M. Schwarz, C. Rödelsperger, M. Schuelke, and D. Seelow, “MutationTaster evaluates disease-causing potential of sequence alterations,” Nature Methods, vol. 7, no. 8, pp. 575-576, 2010.
- M. Kircher, D. M. Witten, P. Jain, B. J. O'Roak, G. M. Cooper, and J. Shendure, “A general framework for estimating the relative pathogenicity of human genetic variants,” Nature Genetics, vol. 46, no. 3, pp. 310–315, 2014.
- G. Lopez-Herrera, G. Tampella, Q. Pan-Hammarström et al., “Deleterious mutations in LRBA are associated with a syndrome of immune deficiency and autoimmunity,” American Journal of Human Genetics, vol. 90, no. 6, pp. 986–1001, 2012.
- U. Salzer, C. Bacchelli, S. Buckridge et al., “Relevance of biallelic versus monoallelic TNFRSF13B mutations in distinguishing disease-causing from risk-increasing TNFRSF13B variants in antibody deficiency syndromes,” Blood, vol. 113, no. 9, pp. 1967–1976, 2009.
- J. E. S. Park, I. Beal, J. P. Dilworth, V. Tormey, and J. Haddock, “The HRCT appearances of granulomatous pulmonary disease in common variable immune deficiency,” European Journal of Radiology, vol. 54, no. 3, pp. 359–364, 2005.
- L. J. Kinlen, A. D. B. Webster, A. G. Bird et al., “Prospective study of cancer in patients with hypogammaglobulinaemia,” The Lancet, vol. 325, no. 8423, pp. 263–266, 1985.
- J. H. Park and A. I. Levinson, “Granulomatous-lymphocytic interstitial lung disease (GLILD) in common variable immunodeficiency (CVID),” Clinical Immunology, vol. 134, no. 2, pp. 97–103, 2010.
- J. R. Hurst, N. Verma, D. Lowe et al., “British lung foundation/united kingdom primary immunodeficiency network consensus statement on the definition, diagnosis, and management of granulomatous-lymphocytic interstitial lung disease in common variable immunodeficiency disorders,” The Journal of Allergy and Clinical Immunology: In Practice, vol. 5, no. 4, pp. 938–945, 2017.
- D. A. Torigian, D. F. LaRosa, A. I. Levinson, L. A. Litzky, and W. T. Miller Jr., “Granulomatous-lymphocytic interstitial lung disease associated with common variable immunodeficiency: CT findings,” Journal of Thoracic Imaging, vol. 23, no. 3, pp. 162–169, 2008.
- P. J. Maglione, “Autoimmune and lymphoproliferative complications of common variable immunodeficiency,” Current Allergy and Asthma Reports, vol. 16, article 19, 2016.
- E. S. Resnick, E. L. Moshier, J. H. Godbold, and C. Cunningham-Rundles, “Morbidity and mortality in common variable immune deficiency over 4 decades,” Blood, vol. 119, no. 7, pp. 1650–1657, 2012.
- P. Maffucci, C. A. Filion, B. Boisson et al., “Genetic diagnosis using whole exome sequencing in common variable immunodeficiency,” Frontiers in Immunology, vol. 7, article 220, 2016.
- B. Lo, K. Zhang, W. Lu et al., “AUTOIMMUNE DISEASE, Patients with LRBA deficiency show CTLA4 loss and immune dysregulation responsive to abatacept therapy,” Science, vol. 349, pp. 436–440, 2015.
- A. Alangari, A. Alsultan, N. Adly et al., “LPS-responsive beige-like anchor (LRBA) gene mutation in a family with inflammatory bowel disease and combined immunodeficiency,” The Journal of Allergy and Clinical Immunology, vol. 130, no. 2, pp. 481–488, 2012.
- L.-M. Charbonnier, E. Janssen, J. Chou et al., “Regulatory T-cell deficiency and immune dysregulation, polyendocrinopathy, enteropathy, X-linked-like disorder caused by loss-of-function mutations in LRBA,” The Journal of Allergy and Clinical Immunology, vol. 135, no. 1, pp. 217–227, 2015.
- P. A. van Schouwenburg, E. E. Davenport, A.-K. Kienzler et al., “Application of whole genome and RNA sequencing to investigate the genomic landscape of common variable immunodeficiency disorders,” Clinical Immunology, vol. 160, no. 2, pp. 301–314, 2015.
- M. L. Waldrep, Y. Zhuang, and H. W. Schroeder Jr., “Analysis of TACI mutations in CVID & RESPI patients who have inherited HLA B⁎44 or HLA⁎B8,” BMC Medical Genetics, vol. 10, p. 100, 2009.
- H. D. Sathkumara, N. R. De Silva, S. Handunnetti, and A. D. De Silva, “Genetics of common variable immunodeficiency: role of transmembrane activator and calcium modulator and cyclophilin ligand interactor,” International Journal of Immunogenetics, vol. 42, no. 4, pp. 239–253, 2015.
- R. Ameratunga, W. Koopmans, S. Woon et al., “Epistatic interactions between mutations of TACI (TNFRSF13B) and TCF3 result in a severe primary immunodeficiency disorder and systemic lupus erythematosus,” Clinical & Translational Immunology, vol. 6, no. 10, p. e159, 2017.
- U. Salzer, H. M. Chapel, A. D. B. Webster et al., “Mutations in TNFRSF13B encoding TACI are associated with common variable immunodeficiency in humans,” Nature Genetics, vol. 37, no. 8, pp. 820–828, 2005.
- E. Castigli, S. A. Wilson, L. Garibyan et al., “TACI is mutant in common variable immunodeficiency and IgA deficiency,” Nature Genetics, vol. 37, no. 8, pp. 829–834, 2005.
- M. S. Hayden and S. Ghosh, “NF-κB in immunobiology,” Cell Research, vol. 21, no. 2, pp. 223–244, 2011.
- M. Fliegauf, V. L. Bryant, and N. Frede, “Haploinsufficiency of the NF-κB1 Subunit p50 in Common Variable Immunodeficiency,” American Journal of Human Genetics, vol. 97, pp. 389–403, 2015.
Copyright © 2018 Ran Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.