Presumed COL4A3/COL4A4 Missense/Synonymous Variants Induce Aberrant Splicing

Background The incorrect interpretation of missense and synonymous variants can lead to improper molecular diagnosis and subsequent faulty genetic counselling. The aim of this study was to evaluate the pathogenicity of presumed COL4A3/COL4A4 missense and synonymous variants detected by next-generation sequencing to provide evidence for diagnosis and genetic counselling. Methods Patients' clinical findings and genetic data were analysed retrospectively. An in vitro minigene assay was conducted to assess the effect of presumed COL4A3/COL4A4 missense and synonymous variants on RNA splicing. Results Five unclassified COL4A3/COL4A4 variants, which were detected in five of 343 patients with hereditary kidney diseases, were analysed. All of them were predicted to affect splicing by Human Splicing Finder. The presumed COL4A3 missense variant c.4793T > G [p. (Leu1598Arg)] resulted in a loss of alternative full-length transcript during the splicing process. The COL4A3 transcript carried synonymous variant c.765G > A [p. (Thr255Thr)], led to an in-frame deletion of exon 13. Nevertheless, variants c.3566G > A [p. (Gly1189Glu)] in COL4A3 and c.3990G > A [p. (Pro1330Pro)], c.4766C > T [p. (Pro1589Leu)] in COL4A4 exhibited no deleterious effect on splicing. Among the five patients harbouring the abovementioned COL4A3/COL4A4 variants, three patients were genetically diagnosed with autosomal recessive Alport syndrome, one patient was highly suspected of having thin basement membrane nephropathy, and the other patient was clinically diagnosed with Alport syndrome. Conclusions COL4A3 presumed missense variant p. (Leu1598Arg) and synonymous variant p. (Thr255Thr) affect RNA splicing, which highlights the prime importance of transcript analysis of unclassified exonic sequence variants for better molecular diagnosis and genetic counselling. Meanwhile, the reliability of splicing predictions by predictive tools for exonic substitutions needs to be improved.

(Pro1589Leu)] in COL4A4 exhibited no deleterious effect on splicing. Among the five patients harbouring the abovementioned COL4A3/COL4A4 variants, three patients were genetically diagnosed with autosomal recessive Alport syndrome, one patient was highly suspected of having thin basement membrane nephropathy, and the other patient was clinically diagnosed with Alport syndrome.

INTRODUCTION
The widespread application of next-generation sequencing (NGS) in clinical settings is leading to incredible progress in the molecular diagnostics of human inherited disorders. However, a large number of unknown significance or unclassified variants in disease genes pose a challenge for genetic diagnosis and subsequent genetic counselling. For example, COL4A3/COL4A4 are the causative genes of both autosomal inherited Alport syndrome and thin basement membrane nephropathy (1), which are chronic haematuria nephropathies with significantly different prognoses. Currently, 22% of the reported COL4A3/COL4A4 variants are classified as uncertain significance, and most of them are missense variants (refer to ClinVar database, https://www.ncbi.nlm.nih.gov/clinvar). It is difficult to assess such variants as pathogenic or neutral based solely on genomic analysis.
A study showed (2) that 5-10% of human genes contain at least one region where synonymous variants could be harmful. Some missense and synonymous variants could interfere with protein expression, conformation or function by affecting precursor mRNA (pre-mRNA) splicing, although they may only change one amino acid or not at genomic DNA level (3,4). Our previous study demonstrated that 6 presumed missense variants caused splicing defects by analysing COL4A5 mRNA in skin fibroblasts (5). Therefore, transcript analysis to demonstrate the effect of unknown significance variants on pre-mRNA splicing is crucial during the assessment of their pathogenicity.
In recent years, a variety of bioinformatics tools have been developed and applied to predict their effects on pre-mRNA splicing. Nevertheless, the predicted in silico results may not necessarily match the transcription process in vivo (6), and functional experiments are still needed. Our previous study increased the detection rate of COL4A5 splicing mutations to 25% by analysing the cDNA of skin fibroblasts (17.6% on average) (7). However, since COL4A3/COL4A4-encoding proteins are expressed in the kidney, cochlea, and eyes, it is difficult to obtain appropriate material from patients for transcript analysis. A minigene assay in vitro, an attractive alternative method to assess the impact of unclassified variants on splicing, is a simple and useful solution to overcome the unavailability of patientderived specimens.
In this study, we used a minigene assay to evaluate the pathogenicity of five presumed missense and synonymous variants in COL4A3/COL4A4 detected by NGS to provide evidence for diagnosis and genetic counselling.

MATERIALS AND METHODS
The clinical data and gene testing results of all patients in this study were retrospectively collected from the online hereditary renal diseases registry database established in 2012 by Peking University First Hospital, the procedures were approved by the Ethical Committee of Peking University First Hospital (approval number: 2016[1179]), and informed consent was obtained from the participants or their parents.

Patients
Patients who were registered in the online registry of hereditary kidney diseases, admitted to our hospital from February 2015 to September 2018 and met the inclusion and exclusion criteria were screened in our study. The inclusion criteria were as follows: (1) haematuria; (2) missense or synonymous variants in COL4A3/COL4A4 detected by NGS; (3) extremely rare variants with minor allele frequency <0.01% (8); (4) variants in a highly conserved amino acid of the protein; (5) harmful variants predicted by bioinformatics tools (Mutation Taster/SIFT/Polyphen-2/GERP++_RS/CADD_phred); and (6) Human Splicing Finder predicted that the variants affected splicing. We excluded patients detected with COL4A3/COL4A4 pathogenic/likely pathogenic variants evaluated by American College of Medical Genetics and Genomics (ACMG) guidelines or those with unavailable medical records.

Clinical Data and Genotypes
Patient information, including sex, age, height, weight, family history, clinical findings (age at onset of disease, renal and extrarenal features), skin and renal biopsy results, and genotypes, was reviewed.
All COL4A3/COL4A4 variants were verified by Sanger sequencing, and direct sequencing of corresponding variants was performed if parental genomic DNA samples were available. Meanwhile, the effect on RNA splicing of these variants was assessed using another three available splicing prediction tools including NNSPLICE (http://www.fruitfly. org/seq_tools/splice.html), NetGene2 (http://www.cbs.dtu.dk/ services/NetGene2/), and MaxEntScan (http://hollywood.mit. edu/burgelab/maxent/Xmaxentscan_scoreseq.html). NNSPLICE and NetGene2 present scores of 0-1 for the predicted site; the higher the score the more likely a variant is a splicing site.

Nomenclature
All genetic sequences involved in this study were referred to the UCSC Genome Browser and human GRCh37/hg19 (http:// genome.ucsc.edu/). The DNA sequence numbering was based on the cDNA sequence for COL4A3 (NM_000091.5)/COL4A4 (NM_000092.5), following the recommendations of the Human Genome Variation Society (HGVS, c.1 represented the first position of the translation initiation codon).

Minigene Assay Generation of Minigene Constructs
The assay relied on the use of a minigene vector, which contained a fragment of COL4A3/COL4A4. The genomic segment encompassing the variant exon sequence and at least one exon nearby it, along with 150 bp of the corresponding flanking intronic sequence, were inserted into the minigene expression plasmid pcDNA3.1 (ViGene Biosciences, Shandong, China). Two minigene constructs (wild-type and mutant-type) were built, and the sequence was verified by direct sequencing after synthesis.

RNA Isolation and Reverse Transcription-Polymerase Chain Reaction (RT-PCR) Analysis
Total RNA was extracted from the transfected cells using TRIzol reagent (Invitrogen, California, USA) according to the manufacturer's instructions. The RT experiment was performed using TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China). The RT-PCR thermal profile was as follows: 25 • C for 10 min, 42 • C for 30 min, and 85 • C for 5 s.

PCR Amplification of cDNA
Five pairs of primers were designed to amplify the resulting cDNA by PCR ( Table 1). The PCR parameters were optimised as follows: 94 • C for 5 min, 35 cycles of denaturation (94 • C for 30 s), annealing (58 • C for 30 s), and elongation (72 • C for 45 s), and a final elongation at 72 • C for 7 min. The amplification products were separated by 2% agarose gel electrophoresis and sequenced using a 3730XL automatic sequencer (Applied Biosystems). Splicing defects were determined according to agarose gel image patterns of the amplification products containing the wild-type and variant alleles. The criteria were based on those proposed by Aoto et al. (9).

RESULTS
According to the inclusion and exclusion criteria, our attention was caught by five COL4A3/COL4A4 variants detected in five patients from 343 patients with hereditary kidney diseases. These variants were evaluated as uncertain significance by ACMG guidelines. COL4A3 gene variants c.4793T>G [p. (Leu1598Arg)] (109) and p. (Thr255Thr) (10) have been previously reported, and COL4A4 variant c.3990G > A [p. (Pro1330Pro)] has been reported in ClinVar, whereas the effect of these three variants on pre-mRNA splicing has not been reported. The splicing effects predicted by in silico splice tools are shown in Table 2. Aberrant RNA splicing effects for these five variants were predicted by Human Splicing Finder, and 4 out of the 5 variants were predicted as deleterious by at least another tool.

Characterisation of Phenotypes and Genotypes
Patient 1 (III1 of family 1) was a 1.4-year-old boy who visited the hospital due to a darkening urine colour. The laboratory investigation showed haematuria [urinary erythrocytes were full of view under urinary sediment microscopy] and proteinuria [protein-to-creatinine ratio (PCR): 1.55 g/gcr; albumin-to-creatinine ratio (ACR): 295.74 mg/g]. His nonconsanguineous parents both had haematuria ( Figure 1A). NGS and Sanger sequencing revealed that Patient 1 had compound heterozygous COL4A3 missense variants [exon 51: c.4793T > G, p. (Leu1598Arg), maternal; exon 37: c.3143G > A, p. (Gly1048Asp), paternal] ( Figure 1A). These two variants were both evaluated as uncertain significance by ACMG guidelines. Among them, the allele frequency of c.4793T > G was 0 in the 1000 Genomes Project, Exome Sequencing Project (ESP) database, and the gnomAD South Asian, African/African-American, Lation/Admixed American, Ashkenazi Jewish, European and other database, 0.0007678 in the gnomAD East Asian database, and 0.00005341 in the gnomAD database. Amino acid Leucine (Leu)-1598 is highly conserved in vertebrates. In silico predictors such as MutationTaster, SIFT, PolyPhen-2 and CADD_phred all predicted it to be deleterious. Only Human Splicing Finder indicated that it may have an effect on splicing.
Patient 2 (III1 of family 2) presented with oedema, haematuria [20-30 red blood cells per high-power field (RBCs/HPF)], and proteinuria (1+, detected by the dipstick test) when she was 7 years old. Prednisone at a standard dose and mycophenolate mofetil treatment did not improve the proteinuria. Subsequently, renal biopsy performed in the local hospital indicated focal segmental glomerular sclerosis (FSGS) under a light microscope and segmental thickening and a lamellated glomerular basement membrane (GBM) under an electron microscope. The uncle of patient 2 had a history of haematuria for 4-5 years, but his urine protein was negative. Renal biopsy indicated thin basement membrane nephropathy (TBMN). Neither ocular nor hearing abnormalities were found ( Figure 1B). Two  was 3 years old. His renal function was normal, whereas at the age of 23 years old, he developed gout, and an elevated level of serum creatinine (150 µmol/L) was detected. Renal biopsy revealed mild mesangial proliferative glomerulonephritis under a light microscope and segmental thickening and thinning and lamellated GBM under an electron microscope. He had no ocular or hearing abnormalities. Laboratory findings showed his serum creatinine was 163 µmol/L, urinary protein was 0.76 g/24 h and estimated glomerular filtration rate was 49.6 ml/min/1.73 m 2 when he was 25 years old. His multiple maternal members (including his mother, aunt, uncle, and grandfather) only presented with haematuria ( Figure 1C). 3-year-old boy. From the age of 3, multiple laboratory tests showed haematuria, no proteinuria, and normal renal function. Renal biopsy was performed on his grandmother at the local hospital due to positive urine occult blood when she was 53 years old. Her GBM was observed to have irregular thickening under electron microscopy. Her urinary protein was negative, and her renal function was normal. At the age of 50, she was found to have hearing loss and needed to wear a hearing aid. Haematuria without proteinuria was detected in the grandmother's sister when she was 55 years old. She had normal renal function without hearing loss or ocular abnormalities. Patient 4's mother had increased numbers of urinary red blood cells without proteinuria or hearing loss (Figure 1D) Figure 1D) in a compound heterozygous state were detected in patient 4's peripheral blood leukocytes genomic DNA. They were evaluated as uncertain significance and pathogenic by ACMG guidelines. Among them, the allele frequency of c.3990G > A was 0 in the 1000 Genomes Project, ESP databases and the gnomAD African/African-American, Lation/Admixed American, Ashkenazi Jewish, European and other database, 0.0006783 in the gnomAD East Asian database, 0.00009902 in the gnomAD South Asian database, and 0.00006097 in the gnomAD database.
In silico MutationTaster predicted it to be deleterious. Human Splicing Finder and NetGene2 indicated that it may have an effect on splicing.
Patient 5 (II2 of family 5) was a 34-year-old woman who manifested haematuria and proteinuria for seven years. The Splicing Effect of Two Variants in COL4A3: c.4793T > G and c.765G > A Agarose gel electrophoresis showed that the wild-type minigene transcription product had two bands, suggesting alternative splicing, while the variant transcription product (c.4793T > G) exhibited only one band after PCR amplification (Figure 2A). Direct sequencing of the PCR products showed that there were two wild-type transcripts (679 bp and 506 bp long) and a variant transcript (506 bp). The longer wild-type transcript contained exons 49, 50, 51 and 52, while the shorter transcript skipped exon 51 (173 bp); the variant transcript also skipped exon 51 during the splicing process ( Figure 2B). These results indicated that the variant c.4793T > G could affect splicing at the mRNA level, which was consistent with the prediction of Human Splicing Finder.
The variant (c.4793T > G) was reclassified according to the ACMG guidelines and the minigene assay results. It could result in a deletion of the normal transcript, which belonged to the null variant (PVS1), and was re-evaluated as "pathogenic" based on the existing evidence (PVS1, PM2, PP3, and PP5). The variant in trans c.3143G > A was reclassified as "likely pathogenic" (in line with PM2, PM3, PP3 and PP4). Autosomal recessive Alport syndrome was diagnosed in patient 1 according to his clinical manifestation, family history of haematuria, and genetic testing results. The recurrence risk of his siblings was 25%.
If primers were located at the inserted exons, the variant transcript (c.765G > A) would be <100 bp and could not be sequenced, so we set the forward primer on the plasmid. The wild-type minigene transcription products were all longer than those of the variants (c.765G > A) after PCR amplification and agarose gel electrophoresis ( Figure 2C). Direct sequencing of the PCR products showed that the two variant transcripts all skipped exon 13 during the splicing process (Figure 2D), which is predicted to cause an in-frame deletion of 26 amino acids during translation. The results indicated that it affected splicing at the mRNA level, consistent with the prediction of four splice site prediction tools.
The variant c.765G > A was reclassified as "pathogenic" (PVS1, PS3, PM2, PP3) according to the ACMG guidelines and minigene assay results. The variant c.4664C > T in trans was classified as likely pathogenic (PM2, PM3, PP3, PP4). Autosomal recessive Alport syndrome could be Variants c.3566G > A in COL4A3 and c.3990G > A and c.4766C > T in COL4A4 Exhibited no Deleterious Effect on Splicing Figure 3A shows that the length of the wild-type minigene transcription product was equal to that of the variant (c.3566G > A in COL4A3) after PCR amplification and agarose gel electrophoresis. Direct sequencing indicated that the two transcripts were all 320 bp and made up of partial exon 41, exon 42 and partial exon 43 (Figure 3B). Only the first nucleotide of exon 42 was discrepant between the wild-type and mutant-type transcripts, which suggested that the variant was missense, not compatible with the prediction of Human Splicing Finder and MaxEntScan. The variant c.4793T > G was the same as patient 1 and reclassified as "pathogenic"; thus, the variant c.3566G > A in trans was re-evaluated as "likely pathogenic" (PM2, PM3, PP3 and PP4). Autosomal recessive Alport syndrome was diagnosed in patient 2 according to the clinical manifestations, renal biopsy results, family history of haematuria, and genetic testing results. The recurrence risk of her siblings was 25%.
The length of the wild-type minigene transcription product was equal to that of the variant (c.3990G > A in COL4A4) after PCR amplification and agarose gel electrophoresis ( Figure 3C). Direct sequencing indicated that the two transcripts were all 229 bp and contained part of exons 41 and 42 ( Figure 3D). Only the mutant-type transcript harboured the variant c.3990G > A, which suggested that the variant was synonymous and not compatible with the prediction of Human Splicing Finder and NetGene2. Considering the clinical manifestations and minigene assay results, it was believed that c.3990G > A was not pathogenic, and TBMN was highly suspected in patient 4, but long-term follow-up was needed to confirm the diagnosis.
The length of the wild-type minigene transcription product was equal to that of the variant (c.4766C > T in COL4A4) after PCR amplification and agarose gel electrophoresis (Figure 3E). Direct sequencing showed that the two transcripts were all 592 bp and contained partial exon 46, exon 47 and partial 48 (Figure 3F). Only the mutant-type transcript harboured the variant c.4766C > T, which suggested that the variant was missense and not compatible with the prediction of Human Splicing Finder and MaxEntScan. The pathogenicity re-evaluation had not been changed. Considering the clinical manifestations, renal biopsy and minigene assay results, Alport syndrome could be diagnosed in patient 5. However, the inheritance pattern could not be determined since her family members had no haematuria and we lacked parental DNA samples.

DISCUSSION
In this study, we assessed the pathogenicity of five COL4A3/COL4A4 unclassified presumed missense and synonymous variants that may affect splicing as predicted by the in silico tool Human Splicing Finder. It was confirmed by the minigene assay that the COL4A3 missense variant p. (Leu1598Arg) and synonymous variant p. (Thr255Thr) theoretically speculated based on the nucleotide substitutions at the genomic DNA level may be splicing variants. Similar phenomena had been observed in COL4A5 gene using in vitro minigene analysis combined with or without analysis of patient-derived RNA from urine sediments and peripheral blood leukocytes (9,11). These pathogenic variants cannot be accurately evaluated based on the NGS results of genomic DNA alone if mRNA analysis is not performed. In addition, it has not yet been reported that reclassification of single base variants in the COL4A3 coding region based on functional experimental results.
Splicing is catalysed by a dynamic spliceosome, which is comprised of five different small nuclear ribonucleoproteins (snRNPs). snRNPs are mainly synthesised by five small nuclear RNAs (snRNAs) named U1, U2, U4, U5, and U6 due to their uracil richness. They can accurately recognise5 ′ splice sites, 3 ′ splice sites and branch points, thereby catalysing the premRNA splicing process (12). The cis-regulatory elements comprise exon splicing enhancers (ESEs), intron splicing enhancers (ISEs), exon splicing silencers (ESSs) and intron splicing silencers (ISSs) located in either exons or introns. These cis-splicing regulatory elements work by recruiting trans-acting factors to activate or suppress splice site recognition or spliceosome assembly (13,14). Therefore, variants in the genome could affect the splicing process by repressing 5 ′ or 3 ′ splice sites, thus generating new splice sites, or having an impact on recruiting trans-acting factors. At present, 619 COL4A3/COL4A4 variants have been reported, of which splicing variants account for ∼13.6% (refer to the human gene mutation database, HGMD), and most of these variants are located at typical splice sites, namely highly conserved 5 ′ donor or 3 ′ acceptor sites (↓GT or AG↓). Studies have shown that 15-60% of human genetic diseases are caused by splicing variants (15), and it has been speculated that the true proportion of splicing variants in COL4A3/COL4A4 may be underestimated. Daga et al. found that a deep intronic variant in COL4A4 could alter splicing products by analysing RNA extracted from Alport patient's urine-derived podocyte lineage cells (16). In our previous study, atypical intronic variants were detected in 10.7% (3/28) of patients with X-linked Alport syndrome by analysing the cDNA of skin fibroblasts (7), and deep intronic splicing variants in Alport syndrome causative genes could be detected by the analysis of patient-derived RNA (17 Many bioinformatics tools have been used to predict splicing variants. Among them, Human Splicing Finder has been widely used for evaluating variant effects on 5 ′ splice sites, 3 ′ splice sites, enhancers, silencers, and branch points (19). However, no splice site prediction tool can be 100% accurate. It was suggested to using multiple in silico tools to pinpoint unclassified variants worthy of transcript analyses (20). Nonetheless, except COL4A3 variant c.765G > A [p. (Thr255Thr)], no splicing abnormalities were observed for the variants (c.3566G > A in COL4A3 and c.3990G > A and c.4766C > T in COL4A4) with a splicing defect predicted using at least two splice site prediction tools. Additionally, in vitro minigene assay analysis of DYSF missense and intronic variants showed that 41% of the tested variants were correctly predicted by Human Splicing Finder (6). In the present study, a false-positive prediction was observed in three of five variants, in accordance with the previous report (6). Therefore, in a molecular diagnostic setting, it is essential to compare in silico splicing predictions and transcript analysis results. The minigene assay is a popular approach to examine whether a splicing variant affects splicing of the neighbouring exon (21,22), and it provides good agreement with experiments using patient-derived RNA (23,24). Given the unavailability of specimens from our patients, such as urine or kidney tissue, and our previous study showed that COL4A3/COL4A4 variants detected in mRNA isolated from Epstein Barr virus-transfected lymphocytes were not always confirmed at the genomic DNA level (data not reported), we used a minigene assay to analyse the consequences of the selected variants. Our study showed that the corresponding transcripts were obtained successfully, and as described above, two predictions were confirmed, which indicated that it was possible and feasible to investigate the effect of COL4A3/COL4A4 variants on splicing by the minigene assay.
It is worth discussing the interpretation of the results of the minigene assay. The alternative transcripts, which are revealed by using different vectors and transfecting cell lines, are slightly different from those found by analysing RNA extracted from the peripheral blood (25,26). This phenomenon may be due to differences in the expression of splicing factors or splicing regulatory factors in different cells. In addition, it may be difficult to evaluate genotype-phenotype correlations when the splicing variants change the relative proportion of different transcripts, and the minigene assay cannot precisely simulate the tissue/cell environment in which the gene is normally expressed. For example, minigene assay analysis demonstrated that the variant c.388 + 5G > A in haemophilia A causative gene F8 leads to two transcripts: one to an abnormal transcript with skipping of the associated exon and one corresponding to the wild-type transcript, which was inconsistent with the patient's severe phenotype (18). In this study, we showed that COL4A3 variant c. 4793T > G led to the production of a transcript lacking exon 51, whereas the wild-type minigene recombinant clone caused the expression of two transcripts: one containing exon 51 and another without exon 51. It has been previously shown that there are three distinct COL4A3 transcripts in foetal (12 weeks) and adult kidney tissues: alternative transcript I containing exons 48-52; alternative transcript II without exon 49; and alternative transcript III without exon 51 (27). Variations could be seen in the expression of these isoforms: alternative transcript I was dominant and accounted for more than 70%, and the relative expression of alternative transcript III in foetal kidney (5.5%) was higher than that in the adult kidney (4.5%), which suggested a slight change in alternative splicing along with organ development. Therefore, due to the loss of the full-length transcript, the COL4A3 variant c.4793T > G was considered a splicing variant. Three patients harbouring this heterozygous variant plus another heterozygous COL4A3 nonsense variant p. (Gly709Ter), intronic splicing variant c. 3752-511_3955 + 576 del, or missense variant p. (Gly1155Asp) demonstrated the typical ultrastructural basket-weave change in the GBM or a thin GBM and an absence of type IV collagen α5 chain in the GBM and only normal immunostaining of Bowman's capsule (typical immunofluorescence staining pattern for autosomal recessive Alport syndrome) (28), suggesting these COL4A3 variants affect type IV collagen α3 chain production.
There were some limitations in this study. First, only unclassified missense and synonymous variants in COL4A3/COL4A4 were screened for splicing analysis. Some frameshifts (29) or nonsense variants (30) may also interfere with splicing and need to be identified. Second, only HEK293T cells were employed to assess the impact on splicing of selected exonic substitutions detected in the patients without immunostaining for type IV collagen α5 chain in renal tissues, and the samples from our patients such as urine or kidney specimen which specific expresses type IV collagen α3 and α4 chains are unavailable, to avoid possible false-negative results, at least one other human cell line, such as HeLa cells, should be used to verify the results from HEK293T cells. Third, Horinouchi et al. showed different nucleotide changes affecting the same codon may or may not affect splicing of an exon using minigenes (11), whereas there were no previously published pathogenic variants at the same codons for the variants in this study, we could not compare the effect on RNA splicing.

CONCLUSION
For the first time, we demonstrated that the presumed missense variant p. Leu1598Arg and synonymous variant p. Thr255Thr in COL4A3 affected pre-mRNA splicing, which extends this gene mutational spectrum and highlights the importance of transcript analysis of exonic sequence variants in the gene to assist with molecular diagnosis and genetic counselling. Meanwhile, we showed that the accuracy of in silico predictive tools for splicing variant prediction was limited.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih. gov/genbank/, OL998453-OL998462.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Ethical Committee of Peking University First Hospital. Written informed consent from the participants' legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

AUTHOR CONTRIBUTIONS
HD and FW designed the research study and wrote the manuscript. HD performed the research. HD and YZ analysed the data. JD provided help and advice on designing the research and writing the manuscript. All authors contributed to editorial changes in the manuscript. All authors read and approved the final manuscript.