patients with advanced fecd ( modified krachmer grade 5 or 6 ) requiring corneal transplantation and control participants without guttae ( grade 0 ) were enrolled in a mayo clinic institutional review board - approved hereditary eye disease study .
fecd grade was established by slit lamp biomicroscopy using specular reflection techniques by one of the authors ( khb , ljm , or svp ) . in control participants ,
patients enrolled in the study agreed to a blood draw and use of their approximately 8-mm - diameter central ce / descemet membrane specimen obtained at endothelial keratoplasty for fecd .
dna was isolated from peripheral blood leukocytes , and rna was isolated from ce / descemet membrane specimens following storage in rnalater ice ( thermo fisher scientific , waltham , ma , usa ) . endothelial tissue from control subjects was obtained at the time of keratoplasty for non - fecd disease or from eyes with normal anterior segments at the time of enucleation .
total rna was isolated independently from 16 tissue samples ( 12 fecd and 4 controls ) by homogenization in qiazol lysis reagent , chloroform extraction and rneasy mini qiacube kit ( qiagen , valencia , ca , usa ) .
rna libraries were prepared for each tissue sample , using the truseq rna sample prep kit version 2 ( illumina , san diego , ca , usa ) .
all samples had rna integrity number ( rin ) values of 6.0 . for truseq stranded total rnaseq ,
ribosomal transcripts were depleted from total rna , using ribo - zero gold rna removal kit followed by replacement of deoxythymidine triphosphate ( dttp ) with deoxyuridine triphosphate ( dutp ) during reverse transcription in the second strand synthesis , using truseq stranded total library preparation kit .
the resulting libraries were minimally amplified to enrich for fragments using adapters on both ends and then quantified for sequencing at three samples / lane by using a hiseq4000 ( illumina ) sequencer .
for that set of samples , poly(a ) mrna was isolated before library preparation using oligo(dt ) magnetic beads .
whole transcriptomic sequencing data from each tissue sample was analyzed using a comprehensive computational program called the mayo analysis pipeline for rna sequencing ( map - rseq ) ( mayo clinic , rochester , mn , usa ) to align , assess , and deliver multiple genomic features .
map - rseq uses a variety of freely available bioinformatics tools along with in - house developed methods to obtain in - depth quality control data , transcriptome read alignment , abundance of gene expression , exon expression , and other transcriptomic features .
the binary alignment map files from map - rseq were analyzed using mixture of isoforms ( miso ) software ( massachusetts institute of technology , cambridge , ma , usa ) that quantitates the expression level of alternatively spliced genes between groups .
for pairwise comparisons , miso calculates bayesian probabilities that there is a genuine difference in the splicing of a given exon between two samples . to establish a differential gene database , we used stringent filtering criteria within miso ( fig .
1 ) to perform genome - wide pairwise comparisons among 11 fecd samples with tnr expansions and 4 controls ( 44 total comparisons ) using  reads to support 
isoform 1 ( inclusive exon ) > 2 ;  reads to support  isoform 2 ( exclusive exon ) > 2 ; sum of inclusive and exclusive reads of 25 ; percentage of spliced in ( psi ) changes > 0.2 ; and bayes factor > 50 as the user - defined criteria for the comparisons . for psi
, miso calculates a value for every differential splicing event , providing a range from 0 to 1 , with 0 being completely excluded and 1 being uniformly included in the splicing products .
visualization of alternatively spliced genes was generated using sashimi plots from the integrated genomics viewer . the single fecd sample without repeat expansion was not included in this pairwise comparison .
the differential gene database was compared with a similar database generated from a pilot study of 8 ce specimens from a study group consisting of 4 fecd patients with ctg tnr expansions in the tcf4 gene , 1 fecd patient without a ctg tnr expansion , and 3 control patients , which allowed twelve pairwise comparisons between fecd with tnr expansion and controls .
initial lists of candidate differential splicing events were produced by miso using stringent criteria ( see methods ) .
these lists were filtered to retain only those events which occurred in more than half of the pairwise comparisons between fecd samples and controls , yielding a combined total of 101 events from the two data sets .
twenty - four validated events present in both data sets are detailed in table 2 .
demographic , clinical and ctg trinucleotide repeat length for study participants summary of differential gene splicing events present in at least 50% of both pilot and pairwise comparisons preparation of cdna by reverse transcription ( rt - pcr ) and analysis by agarose gel electrophoresis was described previously .
specific primers for each validation assay and pcr conditions are provided in supplementary table s1 .
briefly , leukocyte - derived dna was extracted using autogen flexigene ( qiagen ) and resuspended in 1 tris - edta ( te ) buffer for a final concentration of 250 ng/l .
tnr repeats from each sample were amplified by pcr , using an i - cycler ( bio - rad , hercules , ca , usa ) by placing 100 ng of genomic dna with 10 pmol of oligonucleotide primers specific for tcf4 ( 5-tcf - fuchs ' : 5-cagatgagtttggtgtaagatg-3 ; and 3-tcf - fuchs ' 1 : 5-acaagcagaaagggggctgcaa-3 ) in the presence of platinum pcr super mix high fidelity ( invitrogen , carlsbad , ca , usa ) . the pcr program consisted of 100 ng of genomic dna and 10 pmol of each primer in 50 l of hot start ( invitrogen ) at 95c for 6 minute for 1 cycle , then 95c for 1 minute for 1 cycle ,
then 62c for 1 minute , followed by 68c for 3 minute for 35 cycles , and then 68c for 1 cycle for 7 minute , followed by a 4c hold . for short tandem repeat analysis , a 5fam primer ( 5-fam - tcf - fuchs ' : 5-cagatgagtttggtgtaagatg-3 ) was used in place of 5-tcf - fuchs ' , and pcr
after pcr amplification , 2 l of dna was mixed with 12 l of diluted map marker 1000 ( bio ventures , inc . ,
gene scan was carried out using 3730xl dna analyzer ( abi , foster city , ca , usa ) .
gene sets identified by filtering miso results were analyzed for overrepresentation of genes in specific panther system ( in the public domain , http://www.pantherdb.org/ ) families and pathways , using the default settings of the panther web portal . specifically , panther system overrepresentation test ( release date 7/15/2016 ) , annotation version , gene ontology ( go ) database ( released 08/22/2016 ; in the public domain , http://www.geneontology.org ) , and homo sapiens reference lists were used .
bonferroni correction for multiple comparison testing was performed to determine significance ( p < 0.05 ) .
patients with advanced fecd ( modified krachmer grade 5 or 6 ) requiring corneal transplantation and control participants without guttae ( grade 0 ) were enrolled in a mayo clinic institutional review board - approved hereditary eye disease study .
fecd grade was established by slit lamp biomicroscopy using specular reflection techniques by one of the authors ( khb , ljm , or svp ) . in control participants ,
patients enrolled in the study agreed to a blood draw and use of their approximately 8-mm - diameter central ce / descemet membrane specimen obtained at endothelial keratoplasty for fecd .
dna was isolated from peripheral blood leukocytes , and rna was isolated from ce / descemet membrane specimens following storage in rnalater ice ( thermo fisher scientific , waltham , ma , usa ) . endothelial tissue from control subjects was obtained at the time of keratoplasty for non - fecd disease or from eyes with normal anterior segments at the time of enucleation .
total rna was isolated independently from 16 tissue samples ( 12 fecd and 4 controls ) by homogenization in qiazol lysis reagent , chloroform extraction and rneasy mini qiacube kit ( qiagen , valencia , ca , usa ) .
rna libraries were prepared for each tissue sample , using the truseq rna sample prep kit version 2 ( illumina , san diego , ca , usa ) .
all samples had rna integrity number ( rin ) values of 6.0 . for truseq stranded total rnaseq ,
ribosomal transcripts were depleted from total rna , using ribo - zero gold rna removal kit followed by replacement of deoxythymidine triphosphate ( dttp ) with deoxyuridine triphosphate ( dutp ) during reverse transcription in the second strand synthesis , using truseq stranded total library preparation kit .
the resulting libraries were minimally amplified to enrich for fragments using adapters on both ends and then quantified for sequencing at three samples / lane by using a hiseq4000 ( illumina ) sequencer .
for that set of samples , poly(a ) mrna was isolated before library preparation using oligo(dt ) magnetic beads .
whole transcriptomic sequencing data from each tissue sample was analyzed using a comprehensive computational program called the mayo analysis pipeline for rna sequencing ( map - rseq ) ( mayo clinic , rochester , mn , usa ) to align , assess , and deliver multiple genomic features .
map - rseq uses a variety of freely available bioinformatics tools along with in - house developed methods to obtain in - depth quality control data , transcriptome read alignment , abundance of gene expression , exon expression , and other transcriptomic features .
the binary alignment map files from map - rseq were analyzed using mixture of isoforms ( miso ) software ( massachusetts institute of technology , cambridge , ma , usa ) that quantitates the expression level of alternatively spliced genes between groups .
for pairwise comparisons , miso calculates bayesian probabilities that there is a genuine difference in the splicing of a given exon between two samples . to establish a differential gene database , we used stringent filtering criteria within miso ( fig .
1 ) to perform genome - wide pairwise comparisons among 11 fecd samples with tnr expansions and 4 controls ( 44 total comparisons ) using  reads to support 
isoform 1 ( inclusive exon ) > 2 ;  reads to support  isoform 2 ( exclusive exon ) > 2 ; sum of inclusive and exclusive reads of 25 ; percentage of spliced in ( psi ) changes > 0.2 ; and bayes factor > 50 as the user - defined criteria for the comparisons . for psi
, miso calculates a value for every differential splicing event , providing a range from 0 to 1 , with 0 being completely excluded and 1 being uniformly included in the splicing products .
visualization of alternatively spliced genes was generated using sashimi plots from the integrated genomics viewer . the single fecd sample without repeat expansion was not included in this pairwise comparison .
the differential gene database was compared with a similar database generated from a pilot study of 8 ce specimens from a study group consisting of 4 fecd patients with ctg tnr expansions in the tcf4 gene , 1 fecd patient without a ctg tnr expansion , and 3 control patients , which allowed twelve pairwise comparisons between fecd with tnr expansion and controls .
initial lists of candidate differential splicing events were produced by miso using stringent criteria ( see methods ) .
these lists were filtered to retain only those events which occurred in more than half of the pairwise comparisons between fecd samples and controls , yielding a combined total of 101 events from the two data sets .
twenty - four validated events present in both data sets are detailed in table 2 .
the full set of 101 events is provided in supplementary table s2 . demographic , clinical and ctg trinucleotide repeat length for study participants summary of differential gene splicing events present in at least 50% of both pilot and pairwise comparisons
preparation of cdna by reverse transcription ( rt - pcr ) and analysis by agarose gel electrophoresis was described previously .
specific primers for each validation assay and pcr conditions are provided in supplementary table s1 .
briefly , leukocyte - derived dna was extracted using autogen flexigene ( qiagen ) and resuspended in 1 tris - edta ( te ) buffer for a final concentration of 250 ng/l .
tnr repeats from each sample were amplified by pcr , using an i - cycler ( bio - rad , hercules , ca , usa ) by placing 100 ng of genomic dna with 10 pmol of oligonucleotide primers specific for tcf4 ( 5-tcf - fuchs ' : 5-cagatgagtttggtgtaagatg-3 ; and 3-tcf - fuchs ' 1 : 5-acaagcagaaagggggctgcaa-3 ) in the presence of platinum pcr super mix high fidelity ( invitrogen , carlsbad , ca , usa ) . the pcr program consisted of 100 ng of genomic dna and 10 pmol of each primer in 50 l of hot start ( invitrogen ) at 95c for 6 minute for 1 cycle , then 95c for 1 minute for 1 cycle , then
62c for 1 minute , followed by 68c for 3 minute for 35 cycles , and then 68c for 1 cycle for 7 minute , followed by a 4c hold . for short tandem
repeat analysis , a 5fam primer ( 5-fam - tcf - fuchs ' : 5-cagatgagtttggtgtaagatg-3 ) was used in place of 5-tcf - fuchs ' , and pcr was performed as described above .
after pcr amplification , 2 l of dna was mixed with 12 l of diluted map marker 1000 ( bio ventures , inc . ,
gene scan was carried out using 3730xl dna analyzer ( abi , foster city , ca , usa ) .
gene sets identified by filtering miso results were analyzed for overrepresentation of genes in specific panther system ( in the public domain , http://www.pantherdb.org/ ) families and pathways , using the default settings of the panther web portal . specifically , panther system overrepresentation test ( release date 7/15/2016 ) , annotation version , gene ontology ( go ) database ( released 08/22/2016 ; in the public domain , http://www.geneontology.org ) , and homo sapiens reference lists were used .
bonferroni correction for multiple comparison testing was performed to determine significance ( p < 0.05 ) .
the demographic and clinical characteristics of the subjects recruited for this validation study are presented in table 1 .
leukocyte tnr lengths for 11 of the 12 fecd patients ranged from 12 to 32 repeats for the normal allele and 48 to 84 repeats for the expanded allele .
one fecd patient had repeat sizes of 15 and 23 repeats which were in the normal range .
control subjects had tnr sizes for both alleles ranging from 12 to 30 repeats which is within the normal range .
sixty - one differential splicing events in 58 genes were identified in more than half of the 44 pairwise comparisons ( 11 fecd compared to 4 controls ) from this validation sample set ( fig . 1 ; supplementary table s2 ) .
twenty - four of these events had previously met the same filtering criteria in the pilot study , validating a broad signature of differential splicing events in fecd associated with ctg tnr expansion in the tcf4 gene .
differential splicing of numa1 and ppfibp1 was identified in every fecd sample compared with normal in both the original pilot and the current validation study ( table 2 ) .
for both the numa1 and the ppfibp1 events , the median psi values for the fecd samples were much lower than for the controls , indicating preferential exclusion of the target exons in fecd .
these patterns are shown in figures 2a and 3a , which present  sashimi plot  views of the sequence data .
primers that spanned the relevant splice junctions were used for rt - pcr validation of the miso - calculated differential splicing events for numa1 and ppfibp1 ( figs .
note that the pattern of products in figures 2b and 3b , lanes 13 ( controls ) favors the larger product , whereas the pattern of products in figures 2b and 3b , lanes 5 to 7 ( fecd with tnr expansion ) favors the smaller splice product , in agreement with the psi calculations from miso .
interestingly , the same missplicing events in both numa1 and ppfibp1 were also identified in muscle cells from a mouse model of dm1 ( table 2 ) . also of interest , numa1 and ppfibp splicing products produced in the ce from a fecd patient that did not have a tnr expansion ( figs .
2b , 3b , lanes 13 ) rather than that of fecd samples from patients that have repeat expansions ( figs .
( a ) sashimi plot showing the number of reads spanning individual splice junctions next to lines connecting the exons involved .
( b ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from numa1 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
sashimi plot showing the number of reads spanning individual splice junctions are shown next to lines connecting the exons involved .
( b ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from ppfibp1 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
differential splicing of mbnl1 and mbnl2 is a well - characterized consequence of tnr expansion in dm1 and is routinely found in fecd as well .
in addition to numa1 and ppfibp1 , mbnl1 and mbnl2 transcripts were also detected as differentially spliced in all 44 comparisons in the validation data set ( table 2 ) .
experimental rt - pcr validation of the mbnl1 splicing changes in fecd was presented previously .
as noted for numa1 and ppfibp1 , the pattern of rt - pcr products in fecd patients with tnr expansion ( fig .
4a , lanes 57 ) is markedly different than the products from control ce specimens ( fig .
4a , lanes 13 ) and from an fecd patient without a tnr expansion ( fig .
the main mbnl2 product from fecd patients with a tnr expansion was 544 bp , in contrast to control samples and the fecd sample without expansion , which had several products , all less than 500 bp .
( a ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from mbnl2 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
( b ) the structure of three splice products from the mbnl2 gene is shown .
exons of the mbnl2 gene are shown as rectangles , and the length of the exon sequence included in each product is given .
the structure of the three products shown in bold ( 544 bp , 490bp , and 395 bp ) was verified by sanger sequencing . ( c ) the predicted protein sequence of the carboxyl terminus of two of the mbnl2 splice variants ( 544 and 395 bp rt - pcr variants ) is shown , with regions that are common to both products in grey .
the pattern of rt - pcr products from mbnl2 was complex , so we sought further verification of the identity of the three major bands .
sanger sequencing of the 544- , 490- , and 395-bp bands confirmed they were from mbnl2 and also identified the presence of a 54-bp insert in the 544 band which was absent from the 490- and 395-bp products ( fig .
when the 544 band sequence was compared to the 395 band sequence , the 544 band contained not only the 54-bp insert but also included a 95-bp insert near the 3 end ( fig .
comparison of the amino acid sequences showed that the 544-bp segment had an additional 18 amino acids ( aa ) in frame with the comparable region in the 395-bp sequence , but the inclusion of the 95-bp insert at the 3 end of 544 changed the reading frame and the aa sequence of the c terminus of the protein ( fig .
one of the interesting differential splicing events listed in table 2 is the splicing of vegfa ( fig .
the largest product from the vegfa gene produced in the ce is 658 bp in length and includes 7 exons .
this product was confirmed by sequencing to encode vegf165 ( enst00000523950.5 , refers to the human protein - coding splice variant ) and is the preferential splicing pattern in the ce from patients with tnr expansions ( fig .
this isoform does contain the c - terminal heparin binding domain , which helps to anchor the protein in the extracellular matrix and promotes bioavailability . in control samples ( fig .
5b , lanes 13 ) and in an fecd sample from a patient that did not have a tnr expansion ( fig .
5b , lane 8) , the product that lacks the 132-bp exon 7 is the preferential splicing pattern .
this removes 44 amino acids from the vegfa protein and creates the isoform vegf121 that lacks the heparin binding domain encoded by exons 6 and 7 ( fig .
sanger sequencing confirmed that this product is identical to enst00000372077.8 and codes for vegf121 after removal of the 26 amino acid signal peptide ( fig .
( a ) sashimi plot showing the number of reads spanning individual splice junctions next to lines connecting the exons involved .
( b ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from vegfa in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ;  45 repeats ) .
( c ) the structure of two splice products from the vegfa gene is shown .
exons of the vegfa gene are shown as rectangles , and the length of the exon sequence included in each product is given .
the structure of the two products shown in bold ( 658 bp and 526 bp ) was verified by sanger sequencing .
thus , novel splicing events that may occur in a particular tissue or disease state were not included in the analysis . in assessing the transcripts produced from the fgfr2 gene ,
a novel product was identified that used a unique 5 splice junction in the eighth exon .
sequencing of the rt - pcr product from an fecd sample confirmed the use of a canonical gt 5 splice site for this isoform , removing 51 bp from the 3 end of exon 8 ( fig 6b ) .
this novel splice event removed 17 amino acids from this receptor ( fig 6c ) .
the use of this 5 splice junction is more common in ce from fecd samples than in controls , and it is also seen to a lesser extent in the fecd sample from a patient that did not have a repeat expansion novel splicing event in fgfr2 transcripts in the ce ( a ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from fgfr2 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
( b ) sequence traces from the isolated upper ( reference ) and lower ( alternative splice form ) bands from figure 6a are shown .
using panther overrepresentation analysis , assessment of the fecd signature differential splicing gene list within the go molecular function category showed an enrichment for genes involved in cytoskeletal protein binding ( table 3 ) .
this enrichment is statistically significant ( p value , 1.9  10 ) even after bonferroni correction for multiple comparisons .
additionally , a 14.2-fold overrepresentation of genes encoding products that bind cell adhesion molecules was also identified ( p value 1  10 ) .
notably , both of these categories were also overrepresented in the entire group of 101 genes identified in the combined data from the pilot and validation sets ( supplementary table s2 ) . in that larger group , there was a 7.4-fold enrichment for cell adhesion molecule binding proteins ( p value 5.4  10 ) and a 4.2-fold enrichment for cytoskeletal binding proteins ( supplementary table s3 ) .
web - based panther system biochemical pathway analysis of the top 24 differential splicing events within the go cellular component category , the genes listed in table 2 are significantly enriched for products found in the cell cortex ( table 3 ; p value 2.6  10 ) and adherens junctions ( p value 7.4  10 ) .
these enrichments are also noted in the full list of 101 genes ( p value 2.2  10 and 9.1  10 , respectively ) ( supplementary tables s2 and s3 ) .
finally , within the go biological process category , the gene list from table 2 was enriched for golgi organization ( p value 1.7  10 ) and positive regulation of epithelial cell migration ( p value 4.4  10 ) ( table 3 ) .
neither of these enrichments was statistically significant in the full set of 101 genes ( supplementary table s2 ) .
however , the larger gene set did have a notable enrichment for positive regulation of gtpase activity ( p value 1.4  10 ) ( supplementary table s3 ) .
sixty - one differential splicing events in 58 genes were identified in more than half of the 44 pairwise comparisons ( 11 fecd compared to 4 controls ) from this validation sample set ( fig . 1 ; supplementary table s2 ) .
twenty - four of these events had previously met the same filtering criteria in the pilot study , validating a broad signature of differential splicing events in fecd associated with ctg tnr expansion in the tcf4 gene .
differential splicing of numa1 and ppfibp1 was identified in every fecd sample compared with normal in both the original pilot and the current validation study ( table 2 ) .
for both the numa1 and the ppfibp1 events , the median psi values for the fecd samples were much lower than for the controls , indicating preferential exclusion of the target exons in fecd .
these patterns are shown in figures 2a and 3a , which present  sashimi plot  views of the sequence data .
primers that spanned the relevant splice junctions were used for rt - pcr validation of the miso - calculated differential splicing events for numa1 and ppfibp1 ( figs .
note that the pattern of products in figures 2b and 3b , lanes 13 ( controls ) favors the larger product , whereas the pattern of products in figures 2b and 3b , lanes 5 to 7 ( fecd with tnr expansion ) favors the smaller splice product , in agreement with the psi calculations from miso .
interestingly , the same missplicing events in both numa1 and ppfibp1 were also identified in muscle cells from a mouse model of dm1 ( table 2 ) . also of interest , numa1 and ppfibp splicing products produced in the ce from a fecd patient that did not have a tnr expansion ( figs .
2b , 3b , lanes 13 ) rather than that of fecd samples from patients that have repeat expansions ( figs .
( a ) sashimi plot showing the number of reads spanning individual splice junctions next to lines connecting the exons involved .
( b ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from numa1 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
sashimi plot showing the number of reads spanning individual splice junctions are shown next to lines connecting the exons involved .
( b ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from ppfibp1 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
differential splicing of mbnl1 and mbnl2 is a well - characterized consequence of tnr expansion in dm1 and is routinely found in fecd as well .
in addition to numa1 and ppfibp1 , mbnl1 and mbnl2 transcripts were also detected as differentially spliced in all 44 comparisons in the validation data set ( table 2 ) .
experimental rt - pcr validation of the mbnl1 splicing changes in fecd was presented previously .
validation of the differential splicing of mbnl2 transcripts is shown in figure 4a . as noted for numa1 and ppfibp1 , the pattern of rt - pcr products in fecd patients with tnr expansion ( fig .
4a , lanes 57 ) is markedly different than the products from control ce specimens ( fig .
4a , lanes 13 ) and from an fecd patient without a tnr expansion ( fig .
4a , lane 8) . the main mbnl2 product from fecd patients with a tnr expansion was 544 bp , in contrast to control samples and the fecd sample without expansion , which had several products , all less than 500 bp .
( a ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from mbnl2 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
( b ) the structure of three splice products from the mbnl2 gene is shown .
exons of the mbnl2 gene are shown as rectangles , and the length of the exon sequence included in each product is given .
the structure of the three products shown in bold ( 544 bp , 490bp , and 395 bp ) was verified by sanger sequencing .
( c ) the predicted protein sequence of the carboxyl terminus of two of the mbnl2 splice variants ( 544 and 395 bp rt - pcr variants ) is shown , with regions that are common to both products in grey .
the pattern of rt - pcr products from mbnl2 was complex , so we sought further verification of the identity of the three major bands .
sanger sequencing of the 544- , 490- , and 395-bp bands confirmed they were from mbnl2 and also identified the presence of a 54-bp insert in the 544 band which was absent from the 490- and 395-bp products ( fig .
4a ) . when the 544 band sequence was compared to the 395 band sequence ,
the 544 band contained not only the 54-bp insert but also included a 95-bp insert near the 3 end ( fig .
comparison of the amino acid sequences showed that the 544-bp segment had an additional 18 amino acids ( aa ) in frame with the comparable region in the 395-bp sequence , but the inclusion of the 95-bp insert at the 3 end of 544 changed the reading frame and the aa sequence of the c terminus of the protein ( fig .
one of the interesting differential splicing events listed in table 2 is the splicing of vegfa ( fig . 5 ) .
the largest product from the vegfa gene produced in the ce is 658 bp in length and includes 7 exons .
this product was confirmed by sequencing to encode vegf165 ( enst00000523950.5 , refers to the human protein - coding splice variant ) and is the preferential splicing pattern in the ce from patients with tnr expansions ( fig .
this isoform does contain the c - terminal heparin binding domain , which helps to anchor the protein in the extracellular matrix and promotes bioavailability . in control samples ( fig .
5b , lanes 13 ) and in an fecd sample from a patient that did not have a tnr expansion ( fig .
5b , lane 8) , the product that lacks the 132-bp exon 7 is the preferential splicing pattern .
this removes 44 amino acids from the vegfa protein and creates the isoform vegf121 that lacks the heparin binding domain encoded by exons 6 and 7 ( fig .
sanger sequencing confirmed that this product is identical to enst00000372077.8 and codes for vegf121 after removal of the 26 amino acid signal peptide ( fig .
( a ) sashimi plot showing the number of reads spanning individual splice junctions next to lines connecting the exons involved .
( b ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from vegfa in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ;  45 repeats ) .
( c ) the structure of two splice products from the vegfa gene is shown .
exons of the vegfa gene are shown as rectangles , and the length of the exon sequence included in each product is given .
the structure of the two products shown in bold ( 658 bp and 526 bp ) was verified by sanger sequencing .
thus , novel splicing events that may occur in a particular tissue or disease state were not included in the analysis . in assessing the transcripts produced from the fgfr2 gene ,
a novel product was identified that used a unique 5 splice junction in the eighth exon .
sequencing of the rt - pcr product from an fecd sample confirmed the use of a canonical gt 5 splice site for this isoform , removing 51 bp from the 3 end of exon 8 ( fig 6b ) .
this novel splice event removed 17 amino acids from this receptor ( fig 6c ) .
the use of this 5 splice junction is more common in ce from fecd samples than in controls , and it is also seen to a lesser extent in the fecd sample from a patient that did not have a repeat expansion novel splicing event in fgfr2 transcripts in the ce ( a ) rt - pcr using primers that flank selected exons was used to assess exon inclusion from fgfr2 in corneal endothelial rna samples from controls ( lanes 13 ) and fecd patients ( lanes 58 ) .
the repeat expansion status of each sample is shown ( + , > 45 repeats ; 45 repeats ) .
( b ) sequence traces from the isolated upper ( reference ) and lower ( alternative splice form ) bands from figure 6a are shown .
using panther overrepresentation analysis , assessment of the fecd signature differential splicing gene list within the go molecular function category showed an enrichment for genes involved in cytoskeletal protein binding ( table 3 ) .
this enrichment is statistically significant ( p value , 1.9  10 ) even after bonferroni correction for multiple comparisons . additionally , a 14.2-fold overrepresentation of genes encoding products that bind cell adhesion molecules was also identified ( p value 1  10 ) .
notably , both of these categories were also overrepresented in the entire group of 101 genes identified in the combined data from the pilot and validation sets ( supplementary table s2 ) . in that larger group , there was a 7.4-fold enrichment for cell adhesion molecule binding proteins ( p value 5.4  10 ) and a 4.2-fold enrichment for cytoskeletal binding proteins ( supplementary table s3 ) .
web - based panther system biochemical pathway analysis of the top 24 differential splicing events within the go cellular component category , the genes listed in table 2 are significantly enriched for products found in the cell cortex ( table 3 ; p value 2.6  10 ) and adherens junctions ( p value 7.4  10 ) .
these enrichments are also noted in the full list of 101 genes ( p value 2.2  10 and 9.1  10 , respectively ) ( supplementary tables s2 and s3 ) .
finally , within the go biological process category , the gene list from table 2 was enriched for golgi organization ( p value 1.7  10 ) and positive regulation of epithelial cell migration ( p value 4.4  10 ) ( table 3 ) .
neither of these enrichments was statistically significant in the full set of 101 genes ( supplementary table s2 ) .
however , the larger gene set did have a notable enrichment for positive regulation of gtpase activity ( p value 1.4  10 ) ( supplementary table s3 ) .
several groups have confirmed a strong association between ctg tnr expansions in the tcf4 gene and fecd .
we have also detected the colocalization of the splicing factor mbnl1 with transcribed cug sequences in rna foci in the nuclei of ce cells from patients with ctg tnr expansions in the tcf4 gene .
these observations closely parallel those found in affected muscle tissue from patients with dm1 . in dm1 ,
more than a decade of research has established that sequestration of mbnl1 in rna foci of the effected muscle cells leads to widespread changes in rna splicing . in the present study
, we identified a signature of missplicing events for fecd ce harboring ctg repeat expansion in tcf4 . among the list of 101 missplicing events identified in both the pilot and the validation studies , 24 missplicing events ( in 23 different genes ) survived a higher level of screening by being present in at least 50% of comparisons in two independent data sets prepared with different library preparation techniques .
our results confirm that there are widespread changes in pre - mrna splicing in the ce from fecd patients who have tnr expansions .
the pathogenesis of dm1 is believed to be due to a ctg tnr expansion in the dystrophia myotonica protein kinase ( dmpk ) gene that causes rna toxicity due to the sequestration of mbnl1 in rna foci , leading to changes in rna splicing .
it is notable that at least 16 of the 24 genes listed in table 2 are also differentially spliced in a mouse model of dm1 . given the differences in species and tissue , this is a remarkable degree of concordance and is consistent with the possibility of overlapping pathogenic mechanisms between fecd and dm1 .
another two events ( in sos1 and nfix ) have also been noted in studies of patients with dm1 .
thus , at least 18 of 24 ( 75% ) of the top splicing changes noted in fecd have also been linked to those in dm1 , emphasizing the similarities between these two ctg tnr expansion diseases and suggesting a common pathogenesis , albeit in different tissues . within the validation set of samples discussed here ,
differential splicing of the transcripts of mbnl1 , mbnl2 , numa1 , and ppfibp1 was detected in every fecd sample compared to controls , and differential splicing events were also noted in 43 of 44 comparisons for syne1 and add3 .
the inclusion of mbnl1 and mbnl2 in this list is expected , given their affinity for cug - rich sequences , which are abundant in the transcripts from the expanded repeats and their known roles in mrna splicing .
it is interesting that the other four genes are also related to cellular events that have been linked to fecd .
both syne1 and numa1 produce nuclear proteins implicated in apoptosis ; ppfibp1 is thought to regulate the disassembly of focal adhesions , and adducins ( including the add3 gene product ) are important in the stabilization of endothelial junctions . nesprin-1
( the product of the syne1 gene ) has also been implicated in the regulation of cell adhesions .
overall , the finding that more than 60% ( 14 of 23 ) of the most robust differential splicing events occur in genes coding for proteins that are involved in cytoskeletal protein binding or cell adhesion is in agreement with previous proteomic studies of fecd .
quantitative proteomic analysis of ce / descemet 's membrane from fecd patients led poulsen et al . to conclude that proteins associated with cell anchoring and extracellular matrix organization were altered in fecd .
our results suggest that some of the proteins that regulate the extracellular matrix and cytoskeleton in the ce might be qualitatively different as well , by virtue of widespread differential splicing events .
the possible contribution of these changes to the pathogenesis of fecd will require further work , but the promising results in treating fecd with rho kinase ( rock ) inhibitors , which target rocks and influence cellular events involved with the cytoskeleton , cell adhesion and proliferation , suggests possible avenues for future experiments . in that regard ,
it is notable that , within the larger set of differential splicing events presented in supplementary tables s2 and s3 , the focus on cell adhesion and cytoskeletal protein binding is extended and over a 7-fold enrichment for proteins that bind to gtpases emerges .
the change in the splicing pattern of vegfa in patients with a tnr expansion is also of interest .
previous work comparing the effect of vegf165 versus that of vegf121 on the barrier function of retinal vascular endothelial cells in culture suggests that vegf165 is more disruptive to barrier function than vegf121 through an effect on the tight junction protein claudin-1 , which is expressed in ce .
the tnr - associated increase in vegf165 in the ce could therefore impair the barrier function of the ce , which has been shown to be an early defect leading to fecd .
the true functional significance of the novel splice form of fgfr2 that we have characterized is currently unknown .
the consequence of the differential splicing is the exclusion of 17 amino acids in the extracellular third immunoglobulin ( ig)-like domain .
the first 3 of the 17 amino acids that are lost in this isoform ( g345 , n346 , and s347 ) have been shown to interact with fgf2 .
it is known that other alternative splicing events in this ig - like domain do affect ligand specificity .
thus , it is reasonable to hypothesize that the function of the receptor could be changed by this splicing event .
the additional significance of this fgfr2 event is that it emphasizes the fact that the list of differential splicing events identified by miso and presented here is not exhaustive .
the software for identifying differential splicing events is not perfect , and other biologically significant splicing alterations in fecd likely exist .
we recognize that the studies presented here are based on a limited sample size ( total of 15 fecd rna samples from patients with tnr expansions and 7 independent controls ) and that the rna expression of the control samples may also be influenced by their diverse underlying diseases .
the identification of splicing changes within mbnl1 , mbnl2 , numa1 , and ppfibp1 , which were identified in all 44 comparisons between fecd with expansion and control samples plus the 12 comparisons performed in the pilot group , minimizes this concern .
however , it will be important to experimentally validate the effect of these splicing changes on the structure and function of individual proteins before the true consequences of these splicing changes for fecd can be known .
although the presence of widespread splicing abnormalities do not confer proof that missplicing is the primary or sole mechanism of disease in fecd , the parallels with well - established pathogenesis of dm1 suggest that splicing changes secondary to sequestration of mbnl1 and mbnl2 are likely to be important in fecd as well . shared missplicing events in dm1 and fecd should bring into question the possibility of shared phenotypic features .
a separate study of 52 patients with myotonic dystrophy ( unspecified whether type 1 or 2 ) found that corneal thickness was greater than normal in dm , but no other endothelial abnormalities were noted .
garcia filho et al . were unable to confirm a difference in corneal thickness in 12 dm1 subjects , but corneal endothelial parameters were not measured .
although it has not been systematically studied , the presence of neuromuscular disease in patients with fecd has not been reported , and this holds true in our own clinical experience . while larger sampling with thorough methodology is needed to confirm an association such as that described by gattey et al .
, one must remember that the repeat expansion diseases described to date are tissue - specific , with distinct clinical phenotypes despite similar genetic bases .
for example , both dm1 and spinocerebellar ataxia type 8 occur as a result of ctg expansion in the noncoding portion of their respective genes , yet the diseases are clinically distinct and involve different parts of the nervous system .
genes , including those harboring repeat expansions , exhibit tissue - specific differential expression . additionally , mosaicism is common in repeat - associated diseases , with much longer repeat lengths present in certain tissues , such as the nervous system .
genes with expansion may also use different promoters in different tissues , resulting in varied expression of the repeat sequence .
the tcf4 gene has 48 different isoforms produced from 9 different promoters , with only some of the variants having the repeat expansion sequence .
nevertheless , there are clearly many variables involved in determining the tissue specific pathogenicity of repeat expansion sequences , so shared clinical features between fecd and other repeat expansion diseases is not necessarily expected .
the linkage of the differential splicing changes reported here to tnr expansions is emphasized by the absence of mis - splicing in this gene set in the fecd sample from a patient that lacked a tnr expansion .
because we do not yet have splicing data for a meaningful number of fecd patients that lack repeat expansions , these findings must be viewed cautiously .
still , these results do emphasize the reality that there are divergent genetic variants that lead to fecd , and no unifying pathogenic mechanism has been identified .
nevertheless , the current mechanism involving missplicing of the candidate genes seems to hold true for fecd patients with tnr expansions .
therefore changes in gene function through alternate splicing induced by tnr expansion can be a valid pathogenic mechanism in most fecd patients .
furthermore , the overrepresentation of cytoskeletal and cell adhesion molecule - binding proteins in the differentially spliced 
group is further proof that rna splicing changes could contribute to the pathogenesis of fecd in the subset of patients with tnr expansions .
 purposeto identify rna missplicing events in human corneal endothelial tissue isolated from fuchs ' endothelial corneal dystrophy ( fecd).methodstotal rna was isolated and sequenced from corneal endothelial tissue obtained during keratoplasty from 12 patients with fecd and 4 patients undergoing keratoplasty or enucleation for other indications . 
 the length of the trinucleotide repeat ( tnr ) ctg in the transcription factor 4 ( tcf4 ) gene was determined using leukocyte - derived dna analyzed by a combination of southern blotting and genescan analysis . 
 commercial statistical software was used to quantify expression of alternatively spliced genes . 
 validation of selected alternative splicing events was performed by using rt - pcr . 
 gene sets identified were analyzed for overrepresentation using web - based analysis system.resultscorneal endothelial tissue from fecd patients containing a ctg tnr expansion sequence in the tcf4 gene revealed widespread changes in mrna splicing , including a novel splicing event involving fgfr2 . 
 differential splicing of numa1 , ppfibp1 , mbnl1 , and mbnl2 transcripts were identified in all fecd samples containing a tnr expansion . 
 the differentially spliced genes were enriched for products that localize to the cell cortex and bind cytoskeletal and cell adhesion proteins.conclusionscorneal endothelium from fecd patients harbors a unique signature of mis - splicing events due to ctg tnr expansion in the tcf4 gene , consistent with the hypothesis that rna toxicity contributes to the pathogenesis of fecd . 
 changes to the endothelial barrier function , a known event in the development of fecd , was identified as a key biological process influenced by the missplicing events . 
Methods
Isolation of Corneal Tissue
RNA Isolation and Sequencing
Analysis of Differentially Spliced Genes
Validation of Differential Splicing Events by RT-PCR
DNA Isolation and TNR Characterization
Pathway Analysis Using Panther system
Results
Qualitative Changes in CE Gene Expression in FECD
Differential Splicing of NUMA1 and PPFIBP1
Differential Splicing of MBNL1 and MBNL2
Differential Splicing of VEGFA
A Novel Splicing Event in FGFR2
Overrepresentation Analysis of Differential Splicing
Discussion
Supplementary Material