Patients
This study was conducted in strict accordance with the principles of the Helsinki Declaration. The patients and their families signed consent for using their data, and the work was approved by the ethical review committee of our hospital.
The proband (II-1) was a 2-month-old female infant, who was admitted to the hospital with impaired consciousness, paroxysmal shortness of breath, and vomiting (occurred 2 times within 12 h of observation). Neonatal screening results showed 6.98 µmol/L of Met (reference value 7.5–45 µmol/L), 0.22 µmol/L of Met/Cys (reference value 0.5-4 µmol/L), and 0.03 µmol/L of Met/Leu (reference value 0.06–0.24 µmol/L). The chest radiograph (CR) showed bronchitis. Computed tomography of her head revealed an enlarged bilateral frontotemporal sub-arachnoid space, slightly enlarged lateral ventricles, left lateral ventricle subependymal cyst, and bilateral maxillary sinusitis. Upper gastrointestinal angiography indicated gastric volvulus. The symptoms did not improve after 1 month of treatment, as evidenced by the examination of the head magnetic resonance imaging (MRI) revealing traffic hydrocephalus, bilateral frontotemporal apical space enlargement, and a left lateral ventricle subependymal cyst. After the following treatment for another month, the examination of the head MRI revealed a left ventricle subependymal hemorrhage and cerebral atrophy. The patient eventually died as a consequence of multiple symptoms, possibly including bronchopneumonia-induced type I respiratory failure and atelectasis, communicating hydrocephalus, and cerebral atrophy.
The second child (II-2) was a 2-month-old male infant who was admitted to the hospital after choking on milk for 1 week and showed shortness of breath during the day. CR revealed a patchy-like high-density shadow in both lungs. Head MRI T2WI showed increased white matter signal, widening and deepening of the cerebral cortex, thinning of the cortex, and widening of subdural space, indicating cerebral dysplasia (Fig. 1B). Inherited metabolic diseases determination report showed 7.93 µmol/L of Met (reference value 8–35 µmol/L), 0.18 µmol/L of Met/Phe (reference value 0.2–0.6), and 242.58 µmol/L of Glu (reference value 45–200 µmol/L). Chromosomal analysis revealed a normal karyotype (46, XY). After one month, the patient died of multiple symptoms, including severe bronchopneumonia, cerebral dysplasia, bronchopulmonary dysplasia, and moderate anemia.
The third child (II-3) was a fetus at 19 weeks of gestation, of whom MRI examination showed bilateral lateral ventricle trigonometry enlargement. The parents of the fetus decided to choose a medical abortion to avoid birth defects.
The fourth male child (II-4) was delivered by cesarean section at 40 weeks of gestation. The newborn weighed 4060 g and was diagnosed with macrosomia. There was no perinatal abnormality (such as hypoxia and asphyxia) at birth. The level of homocysteine in the postnatal peripheral blood was 10.6 µmol/L (reference value 0–15 µmol/L). Four months later, the homocysteine level was found to be elevated to 16.4 µmol/L (reference value 0–15µmol/L) but the levels of complete blood count, vitamin B12, folic acid, and ferritin were all normal. The father of the proband had a slight increase in peripheral blood homocysteine (Hcy: 27.6 µmol/L) and no other obvious clinical manifestations. It was improved after folic acid and vitamin B12 supplementation but returned to a higher Hcy level (Hcy: 52.2 µmol/L) after withdrawal. The proband’s mother is in good health condition. The pedigree chart of the family was presented in Fig. 1 A.
Library preparation
Firstly, genomic DNA was extracted from 200µL peripheral blood, using a Qiagen DNA Blood Midi/Mini kit (Qiagen GmbH, Hilden, Germany). About 50 ng of genomic DNA was interrupted to around 200 bp fragments by the enzyme. The DNA fragments were end-repaired by adding one A base at the 3’end.
Secondly, the DNA fragments were ligated with barcoded sequencing adaptors, and ligated fragments (about 320 bp) were collected by XP beads. After PCR amplification, the DNA fragments were hybridized and captured by NanoWES(Berry Genomics, China) according to the manufacturer’s Protocol. The hybrid products were eluted and collected, and then subjected to PCR amplification and purification. Next, the libraries were quantified by qPCR, and size distribution was determined using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Finally, the Novaseq 6000 platform (Illumina, San Diego, USA) with 150 bp pair-end sequencing mode was used for sequencing the genomic DNA samples of the family. Raw image files were processed using CASAVA v1.82 for base calling and generating raw data [8].
Data analysis
The sample genomes were aligned to the human reference genome (hg19/GRCh37) using the Burrows-Wheeler Aligner tool and PCR duplicates were removed by using Picard v1.57 (http://picard.sourceforge.net/). Verita Trekker® Variants Detection System by Berry Genomics and the third-party software GATK (https://software.broadinstitute.org/gatk/) were employed for variant calling. Variant annotation and interpretation were performed with ANNOVAR (Wang, et al., 2010) and the Enliven® Variants Annotation Interpretation System authorized by Berry Genomics.
Annotation databases mainly included:
-
(i)
Human population databases, including gnomAD (http://gnomad.broadinstitute.org/), the 1000 Genome Project (http://browser.1000genomes.org), Berrybig data population database, dbSNP (http://www.ncbi.nlm.nih.gov/snp) etc.;
-
(ii)
Databases of in silico prediction algorithms, including SIFT (http://sift.jcvi.org), FAT-HMM (http://fathmm.biocompute.org.uk), MutationAssessor (http://mutationassessor.o-rg), CADD (http://cadd.gs.washington.edu), SPIDEX (Xiong et al., Science 2015);
-
(iii)
Disease and phenotype databases, including OMIM (http://www.omim.org), ClinVar (http://www.ncbi.nlm.nih.gov/clinvar), HGMD (http://www.hgmd.orgO (http://hpo.jax.org/app/). The variants were classified into five categories: pathogenic, likely pathogenic, uncertain significance, likely benign, and benign, according to the American College of Medical Genetics and Genomics (ACMG) guidelines for the interpretation of genetic variants. Variants with minor allele frequencies (MAF) < 1% in the exonic region or with splicing impact were taken for deep interpretation based on the ACMG category, evidence of pathogenicity, clinical synopsis, and inheritance model of the associated disease. Sanger sequencing was further used to identify the presence of the pathogenic variants in the family members.