Myofibrillar myopathy in the genomic context

Myofibrillar myopathy (MFM) is a group of inherited muscular disorders characterized by myofibril dissolution and abnormal accumulation of degradation products. The diagnosis of muscular disorders based on clinical presentation is difficult due to phenotypic heterogeneity and overlapping symptoms. In addition, precise diagnosis does not always explain the disease etiopathology or the highly variable clinical course even among patients diagnosed with the same type of myopathy. The advent of high-throughput next-generation sequencing (NGS) has provided a successful and cost-effective strategy for identification of novel causative genes in myopathies, including MFM. So far, pathogenic mutations associated with MFM phenotype, including atypical MFM-like cases, have been identified in 17 genes: DES, CRYAB, MYOT, ZASP, FLNC, BAG3, FHL1, TTN, DNAJB6, PLEC, LMNA, ACTA1, HSPB8, KY, PYROXD1, and SQSTM + TIA1 (digenic). Most of these genes are also associated with other forms of muscle diseases. In addition, in many MFM patients, numerous genomic variants in muscle-related genes have been identified. The various myopathies and muscular dystrophies seem to form a single disease continuum; therefore, gene identification in one disease impacts the genetic etiology of the others. In this review, we describe the heterogeneity of the MFM genetic background focusing on the role of rare variants, the importance of whole genome sequencing in the identification of novel disease-associated mutations, and the emerging concept of variant load as the basis of the phenotypic heterogeneity.     

Identification novel LQT syndrome-associated variants in Polish population and genotype-phenotype correlations in eight families

Congenital long QT syndrome (LQTS) is a primary cardiac channelopathy. Genetic testing has not only diagnostic but also prognostic and therapeutic implications. At present, 15 genes have been associated with the disease, with most mutations located in 3 major LQTS-susceptibility genes. During a routine genetic screening for KCNQ1, KCNH2 and SCN5A genes in index cases with LQTS, seven novel variants in KCNH2 and SCN5A genes were found. Genotype-phenotype correlations were analysed in these patients and their families. An open reading frame and splice site analysis of the exons was conducted using next-generation sequencing. In novel variants, phenotypes of carriers and their affected relatives were analysed. In 39 unrelated patients, 40 pathogenic/putative pathogenic mutations were found. Thirty-three of them, predominantly missense, were reported previously: 11 were in the KCNQ, 17 in the KCNH2 and 5 in the SCN5A gene. Seven novel missense variants were found in eight families. Among them, four variants were in typical for LQTS location. Two variants in the KCNH2 gene (p.D803Y and p.D46F) and one in the SCN5A gene (G1391R) were in amino acid (AA) position which up to present has not been reported in LQTS. Phenotype analysis showed the life-threatening course of the disease in index cases with a history of sudden cardiac death in six families. Mutation carriers presented with ECG abnormalities and some of them received beta-blocker therapy. We report three novel variants (KCNQ1 p.46, KCNH2 p.D803Y, SCN5A p.G1391R) which have never been reported for this AA location in LQTS; the phenotype-genotype correlation suggests their pathogenicity.     

Evolutionary divergence of the rye <Emphasis Type="Italic">Pm17</Emphasis> and <Emphasis Type="Italic">Pm8</Emphasis> resistance genes reveals ancient diversity

Key message

We have isolated a novel powdery mildew resistance gene in wheat that was originally introgressed from rye. Further analysis revealed evolutionary divergent history of wheat and rye orthologous resistance genes.

Abstract

Wheat production is under constant threat from a number of fungal pathogens, among them is wheat powdery mildew (Blumeria graminis f. sp. tritici). Deployment of resistance genes is the most economical and sustainable method for mildew control. However, domestication and selective breeding have narrowed genetic diversity of modern wheat germplasm, and breeders have relied on wheat relatives for enriching its gene pool through introgression. Translocations where the 1RS chromosome arm was introgressed from rye to wheat have improved yield and resistance against various pathogens. Here, we isolated the Pm17 mildew resistance gene located on the 1RS introgression in wheat cultivar ‘Amigo’ and found that it is an allele or a close paralog of the Pm8 gene isolated earlier from ‘Petkus’ rye. Functional validation using transient and stable transformation confirmed the identity of Pm17. Analysis of Pm17 and Pm8 coding regions revealed an overall identity of 82.9% at the protein level, with the LRR domains being most divergent. Our analysis also showed that the two rye genes are much more diverse compared to the variants encoded by the Pm3 gene in wheat, which is orthologous to Pm17/Pm8 as concluded from highly conserved upstream sequences in all these genes. Thus, the evolutionary history of these orthologous loci differs in the cereal species rye and wheat and demonstrates that orthologous resistance genes can take different routes towards functionally active genes. These findings suggest that the isolation of Pm3/Pm8/Pm17 orthologs from other grass species, additional alleles from the rye germplasm as well as possibly synthetic variants will result in novel resistance genes useful in wheat breeding.

     

Identification and Phylogenetic Classification of <Emphasis Type="Italic">Pennisetum</Emphasis> (Poaceae) Ornamental Grasses Based on SSR Locus Polymorphisms

Pennisetum species are widely used as ornamental grasses and may be a valuable genetic resource for their breeding to broaden its genetic basis. At present, new ornamental Pennisetum cultivars are primarily bred via somaclone, which increases the number of variants. It is difficult to estimate whether the suspected variants are authentic in genetic features by morphological traits because of their many limitations. Moreover, although the phylogenetic classification of the Pennisetum genus has been approved in some morphological and cytological studies, genetic evidence is lacking. In the present study, we developed 15 specific simple sequence repeat (SSR) markers with a large amount of polymorphisms and strong distinguishing abilities for Pennisetum ornamental grasses using magnetic bead enrichment. These markers, together with the other 11 reported polymorphic SSRs, were further used for the identification of a broad collection of 55 Pennisetum samples, including nine original taxa and 46 suspected variants. After comparing the genetic characteristics between each variant and its corresponding original taxon, we verified 20 suspected variants that possess the potential to become new, commercially desirable cultivars. The nine original taxa and the 20 verified variants were identified based on the polymorphisms of six core loci, and unique molecular identities with 15 denary digits for each taxon were further established. The rationality of the traditional phylogenetic classification system of the Pennisetum genus was further verified using 147 polymorphic alleles. The present study promotes the protection, registration, breeding, and international communication of Pennisetum ornamental grasses.     

LncRNA <Emphasis Type="Italic">ANRIL</Emphasis> Expression and <Emphasis Type="Italic">ANRIL</Emphasis> Gene Polymorphisms Contribute to the Risk of Ischemic Stroke in the Chinese Han Population

The aim of the present study was to explore the role of lncRNA ANRIL in the pathogenesis of ischemic stroke (IS) and coronary artery disease (CAD) and to determine the association between ANRIL variants and the genetic susceptibility of IS and CAD in the Chinese Han population. A genetic association study including 550 IS patients, 550 CAD patients, and 550 healthy controls was conducted. The expression levels of lncRNA ANRIL, CDKN2A, and CDKN2B were detected using qRT-PCR. Genotyping was performed by Sequenom MassARRAY on an Agena platform. Our study showed that IS patients had an increased lncRNA ANRIL expression (P?=?0.002) and a decreased CDKN2A expression (P?<?0.001) compared with normal controls. A significant difference with regard to the genotype distribution of rs2383207 was found between male IS patients and controls (P?=?0.011). The minor allele of rs2383207 significantly increased the IS risk under a recessive model (OR?=?1.52, 95% CI?=?1.05–2.21, P?=?0.027). The minor allele of rs1333049 was significantly associated with the risk of IS among the male patients under a recessive model (OR?=?1.56, 95% CI?=?1.04–2.35, P?=?0.031). However, no significant association was found between the ANRIL variants and the risk of CAD (all P?>?0.050). In addition, we found a decreased lncRNA ANRIL expression in IS patients who carried the GG genotype of rs1333049 compared with IS patients who carried the CC or CG genotype (P?=?0.041). In summary, we found that IS patients had an increased lncRNA ANRIL expression and a decreased CDKN2A expression compared with the controls, which might play an impellent role in pathological processes of IS. The ANRIL variants rs2383207 and rs1333049 were significantly associated with the risk of IS among males but not females in the Chinese Han population.     

Modifier locus of the skeletal muscle involvement in Emery–Dreifuss muscular dystrophy

Autosomal dominant Emery–Dreifuss muscular dystrophy is caused by mutations in LMNA gene encoding lamins A and C. The disease is characterized by early onset joint contractures during childhood associated with humero-peroneal muscular wasting and weakness, and by the development of a cardiac disease in adulthood. Important intra-familial variability characterized by a wide range of age at onset of myopathic symptoms (AOMS) has been recurrently reported, suggesting the contribution of a modifier gene. Our objective was to identify a modifier locus of AOMS in relation with the LMNA mutation. To map the modifier locus, we genotyped 291 microsatellite markers in 59 individuals of a large French family, where 19 patients carrying the same LMNA mutation, exhibited wide range of AOMS. We performed Bayesian Markov Chain Monte Carlo-based joint segregation and linkage methods implemented in the Loki^© software, and detected a strong linkage signal on chromosome 2 between markers D2S143 and D2S2244 (211 cM) with a Bayes factor of 28.7 (empirical p value = 0.0032). The linked region harbours two main candidate genes, DES and MYL1 encoding desmin and light chain of myosin. Importantly, the impact of the genotype on the phenotype for this locus showed an overdominant effect with AOMS 2 years earlier for the homozygotes of the rare allele and 37 years earlier for the heterozygotes than the homozygotes for the common allele. These results provide important highlights for the natural history and for the physiopathology of Emery–Dreifuss muscular dystrophy.     

Three Novel <Emphasis Type="Italic">NF1</Emphasis> Gene Mutations in a Cohort of Bulgarian Neurofibromatoses Patients

Neurofibromatosis (NF) is a clinically heterogeneous autosomal dominant disorder. Three distinct forms have been identified: neurofibromatosis type 1 (NF1), type 2 (NF2) and schwannomatosis. In the present study, we report clinical and genetic findings in the NF1 and NF2 genes in a cohort of 27 Bulgarian patients, with 18 cases (67%) genetically verified. Both NF1 and NF2 genes were screened by Sanger sequencing on DNA samples. The Sanger negative samples were screened by Multiplex Ligation-dependent Probe Amplification (MLPA) for deletions and duplications. The results from genetic testing revealed three novel mutations and fifteen previously reported ones (13 in the NF1 gene and 2 in the NF2 gene). The novel variants in the NF1 gene are a splice site mutation c.4725-1G>A, a small deletion of five bases c.823delATCTT, p.Leu275ValfsTer14, and a single base duplication c.6547dupC, p.Arg2183ProfsTer11. The novel splice site mutation is manifested by multiple “café au lait” macules and neurofibromas. Both novel out of frame mutations were found in patients with multiple “café au lait” spots and focal epilepsy. A segmental neurofibromatosis (SNF1) is restricted to one or more body segments. Here we present a case with SNF1 caused by a somatic deletion of exons 1 to 12 of the NF1 gene which is manifested by multiple neurofibromas in the right hand. Two nonsense mutations are found in the NF2 gene. Our study adds three novel mutations to the NF1 mutation spectra and contributes to the clinical-genetic NF1-characterization. Here we report strikingly different phenotypic spectra caused by the same mutation in a single family. Our findings contribute to the genotype- phenotype correlations which are difficult to establish, due to the extremely complex NF phenotype being a combination of clinical features.     

Prevalence and Phylogenetic Analysis of <Emphasis Type="Italic">Bartonella</Emphasis> Species of Wild Carnivores and Their Fleas in Northwestern Mexico

The host–parasite–vector relationship of Bartonella spp. system in wild carnivores and their fleas from northwestern Mexico was investigated. Sixty-six carnivores belonging to eight species were sampled, and 285 fleas belonging to three species were collected during spring (April–May) and fall (October–November) seasons. We detected Bartonella species in 7 carnivores (10.6%) and 27 fleas (9.5%) through either blood culture or PCR. Of the 27 Bartonella-positive fleas, twenty-two were Pulex simulans, three were Pulex irritans and one was Echidnophaga gallinacea. The gltA gene and ITS region sequences alignment revealed six and eight genetic variants of Bartonella spp., respectively. These variants were clustered into Bartonella rochalimae, Bartonella vinsonii subsp. berkhoffii and another genotype, which likely represents a novel species of Bartonella spp. Although experimental infection studies are required to prove the vector role of P. simulans, our results suggest that this flea may play an important role in the Bartonella transmission. The results indicated possible host-specific relationships between Bartonella genotypes and the families of the carnivores, but further studies are needed to verify this finding. The presence of zoonotic species of Bartonella spp. in wild carnivores raises the issue of their potential risk for humans in fragmented ecosystems.     

<Emphasis Type="Italic">Agrobacterium tumefaciens</Emphasis>-mediated transformation: an efficient tool for targeted gene disruption in <Emphasis Type="Italic">Talaromyces marneffei</Emphasis>

Talaromyces marneffei causes life-threatening infections in immunocompromised hosts. An efficient tool for genetic manipulation of T. marneffei will allow for increased understanding of this thermally dimorphic fungus. Agrobacterium tumefaciens-mediated transformation (ATMT) was optimized for targeted gene disruption in T. marneffei using the plasmid pDHt/acuD::pyrG. Molecular analyses of transformants were performed by PCR, Southern blot and semi-quantitative RT-PCR. A. tumefaciens strain EHA105 was more efficient at transformation than strain AGL-1 in ATMT via solid co-cultivation. An A. tumefaciens:T. marneffei ratio of 1000:1 in an ATMT liquid co-cultivation led to a relatively high transformation efficiency of 90 transformants per 10⁶ yeast cells. Using ATMT-mediated knockout mutagenesis, we successfully deleted the acuD gene in T. marneffei. PCR and Southern blot analysis confirmed that acuD was disrupted and that the foreign pyrG gene was integrated into T. marneffei. Semi-quantitative RT-PCR analysis further confirmed that pyrG was expressed normally. These results suggest that ATMT can be a potential platform for targeted gene disruption in T. marneffei and that liquid co-cultivation may provide new opportunities to develop clinical treatments.     

Spatial genetic structure in seed stands of <Emphasis Type="Italic">Pinus lumholtzii</Emphasis> B.L. Rob. &amp; Fernald in Durango,Mexico

Sad pine is one of the most prominent pine species in Mexico due to its conspicuous pendulous foliage and extreme habitat. However, scientific studies of the species are scarce, and genetic information on sad pine populations is lacking. This endemic tree species occurs naturally on the Sierra Madre Occidental where it covers a total area of about 1,600,000 ha. It typically grows with several species of Quercus and Pinus or in pure stands in uneven-aged forests. Pinus lumholtzii is naturally spatially fragmented, and genetic research on seed and pollen dispersal patterns and spatial genetic structure (SGS)—and the possible implications of these in terms of evolution, conservation and breeding management—is particularly important. Given the fragmented occurrence of the preferred soil type, the goal of this research was to use amplified fragment length polymorphism (AFLP) markers to identify potential differences in spatial genetic structure within and between five P. lumholtzii seed stands at fine and large scales. At the fine scale, we almost always observed non-significant autocorrelation, suggesting that the genetic variants of P. lumholtzii are randomly distributed in space within each sampled seed stand. At the larger scale, our findings provide strong support for the theory of isolation by distance that predicts the expected pattern of SGS at drift–dispersal equilibrium. We recommend a network of P. lumholtzii seed stands of maximum distances of 100 km among stands to prevent greater loss of local genetic variants and use the seeds for reforestations in a radius of maximal 50 km from their proveniences.     

Biosynthesis of <Emphasis Type="Italic">Veratrum californicum</Emphasis> specialty chemicals in <Emphasis Type="Italic">Camelina sativa</Emphasis> seed

Economically feasible systems for heterologous production of complex secondary metabolites originating from difficult to cultivate species are in demand since Escherichia coli and Saccharomyces cerevisiae are not always suitable for expression of plant and animal genes. An emerging oilseed crop, Camelina sativa, has recently been engineered to produce novel oil profiles, jet fuel precursors, and small molecules of industrial interest. To establish C. sativa as a system for the production of medicinally relevant compounds, we introduced four genes from Veratrum californicum involved in steroid alkaloid biosynthesis. Together, these four genes produce verazine, the hypothesized precursor to cyclopamine, a medicinally relevant steroid alkaloid whose analogs are currently being tested for cancer therapy in clinical trials. The future supply of this potential cancer treatment is uncertain as V. californicum is slow-growing and not amendable to cultivation. Moreover, the complex stereochemistry of cyclopamine results in low-yield syntheses. Herein, we successfully engineered C. sativa to synthesize verazine, as well as other V. californicum secondary metabolites, in seed. In addition, we have clarified the stereochemistry of verazine and related V. californicum metabolites.     

Metabolic engineering to produce γ-linolenic acid in <Emphasis Type="Italic">Brassica napus</Emphasis> using a Δ6-desaturase from pike eel

We generated γ-linolenic acid (GLA, C18:3^Δ6,9,12)-producing transgenic Brassica napus transformed with McD6DES, the Δ6-desaturase gene identified from pike eel (Muraenesox cinereus) under the control of the seed-specific vicilin promoter. Seed-specific expression of McD6DES in B. napus produced up to 8.4% of GLA by creating a double bond at the sixth position from the carboxyl end of linoleic acid (LA, C18:2^Δ9,12) in seeds. These results demonstrate that McD6DES expression enables to reconstitute in polyunsaturated fatty acid biosynthetic pathways, highlighting the potential of GLA biosynthesis as a target for metabolic engineering of oilseed crops.     

<Emphasis Type="Italic">Agrobacterium-medlated</Emphasis> high-efficiency transformation of creeping bentgrass with herbicide resistance

We performedAgrobacterium-mediated genetic transformation of creeping bentgrass(Agrostis stolonifera L.) and produced herbicide-resistant transformants from commercial cultivars Crenshaw and Penncross. Seed-derived embryogenie calli were infected withA. tumefaciens EHA105 harboring pCAMBIA 3301, which includes an intron-containinggus reporter and abar selection marker. To establish a stable system, we examined various factors that could potentially influence transformation efficiency during the pre-culture, infection, and co-cultivation steps. The addition of kinetin to the callus pre-culture media increased efficiency about three-fold. Once the optimum infection and co-cultivation conditions were identified, this protocol was used successfully to bulk-produce herbicide-resistant transgenic plants whose herbicide resistance was confirmed using the BASTA® resistance test. Southern blot analysis demonstrated integration and low copy numbers of the integrated transgenes, and northern blot analysis verified their expression. Thus, we have established an efficient genetic transformation system for creeping bentgrass and confirmed a high frequency of single-copy transgene integration and functional gene expression.     

Genetic analysis of <Emphasis Type="Italic">COL11A2</Emphasis> in Korean patients with autosomal dominant non-syndromic hearing loss

The collagen type XI alpha 2 gene (COL11A2) is associated with autosomal dominant non-syndromic hearing loss (ADNSHL), and all mutations of this gene in ADNSHL are missense mutations. To evaluate its potential as a major causative gene of ADNSHL in the Korean population, we performed genetic analysis of COL11A2 in 75 unrelated Korean patients with ADNSHL. Consequently, 5 non-synonymous variants, 7 synonymous variants, and 6 intronic variants were identified in COL11A2. Among them, a novel variant, p.G829R (c.2485G>C) was found in a patient as a heterozygote. However, pedigree analysis showed this variation was not co-segregated with hearing loss. Previously reported variants p.G230W (c.688G>T) and p.P1422L (c.4265C>T) were discovered in Korean patients. However, these variants were also detected in normal individuals. These results suggest that COL11A2 is not a major causative gene of ADNSHL in the Korean population.     

Phenol Oxidase Activity of <Emphasis Type="Italic">Azospirillum brasilense</Emphasis> Sp245 Mutants with Modified Motility and <Emphasis Type="Italic">Azospirillum brasilense</Emphasis> Sp7 Phase Variants with Different Plasmid Composition

Herein, we reveal the alteration in phenol oxidase enzymes complex production from Azospirillum brasilense Sp245 omegon mutants with polar and lateral flagella dysfunction and from A. brasilense Sp7 phase variants with different plasmid composition. The enzymatic activities for various laccases, tyrosinases, Mnperoxidases, and lignin peroxidases as well as the isomorphic composition of intracellular laccases and tyrosinases were estimated for the studied variants and the parent strains. It was noted that various genetic events correlating with phenotypic heterogeneity in A. brasilense populations affect their phenol oxidase activity level.     

The response of transgenic <Emphasis Type="Italic">Brassica</Emphasis> species to salt stress: a review

Salt stress is considered one of the main abiotic factors to limit crop growth and productivity by affecting morpho-physiological and biochemical processes. Genetically, a number of salt tolerant Brassica varieties have been developed and introduced, but breeding of such varieties is time consuming. Therefore, current focus is on transgenic technology, which plays an important role in the development of salt tolerant varieties. Various salt tolerant genes have been characterized and incorporated into Brassica. Therefore, such genetic transformation of Brassica species is a significant step for improvement of crops, as well as conferring salt stress resistance qualities to Brassica species. Complete genome sequencing has made the task of genetically transforming Brassica species easier, by identifying desired candidate genes. The present review discusses relevant information about the principles which should be employed to develop transgenic Brassica species, and also will recommend tools for improved tolerance to salinity.     

A novel set of 223 EST-SSR markers in <Emphasis Type="Italic">Casuarina</Emphasis> L. ex Adans.: polymorphisms,cross-species transferability,and utility for commercial clone genotyping

Simple sequence repeat (SSR) markers are very useful for genetic applications in plants, but SSR resource for the important tree genus Casuarina L. ex Adans. is still limited. In this study, we report a novel set of 223 SSR markers in Casuarina developed from expressed sequence tag (EST) resource of GenBank. The 223 EST-SSR markers were polymorphic among 10 unrelated individuals of C. equisetifolia L. Johnson, with the number of alleles per locus (N_a), observed heterozygosity (H_o), expected heterozygosity (H_e), and polymorphic information content (PIC) averaging at 5.5, 0.72, 0.86, and 0.63, respectively. The rates of cross-species transferability ranged from 96.9% (C. glauca Sieber ex Sprengel) through 97.8% (C. cunninghamiana Miquel) to 99.1% (C. junghuhniana Miquel). Fifty-five C. equisetifolia clones widely planted in China were successfully genotyped with a subset of 20 EST-SSRs. These newly developed markers will have a great potential for genetic and breeding applications in Casuarina species and related taxa.     

Genetic structure of cultivated <Emphasis Type="Italic">Zanthoxylum</Emphasis> species investigated with SSR markers

Zanthoxylum is an economically and ecologically important genus of the Rutaceae family, of which Z. bungeanum and Z. armatum have a long history of cultivation in China. However, how the natural processes such as selection and drift and agriculture practices have influenced the genetic variation of cultivated Zanthoxylum species during long-term domestication remains elusive. Herein, we determined the population genetic structure of current widely cultivated Zanthoxylum species, Z. bungeanum and Z. armatum. Microsatellite markers revealed a high level of genetic variation and significant genetic differentiation for both species despite Z. bungeanum showed higher genetic diversity than Z. armatum. AMOVA indicated that most of the genetic variation exists within individuals rather than among provenances for both species. Population structure analyses generated three distinct groups within the entire accessions. All Z. bungeanum accessions were distinguished into two major geographic groups, north and south groups, with Qinling Mountains as the main geographic barrier to gene flow while a significant genetic differentiation was observed between cultivated and wild Z. armatum accessions. Mantel test of Z. bungeanum displayed a significant correlation between genetic and geographic distances within each inferred group but no correlation between genetic and geographic distance was observed when comparing genetic and geographic distances focusing only on pairwise of north vs. south provenances, ruling out the hypothesis that gene flow between north and south provenances followed an isolation-by-distance model. Our research provided a fundamental genetic profile that will improve the conservation and responsible exploitation of the extant germplasm of Zanthoxylum.     

Haplotype-specific <Emphasis Type="Italic">MAPT</Emphasis> exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

Background

Genome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson’s disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.

Methods

In this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.

Results

Using whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.

Conclusions

We identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.