Use of animal models to understand titin physiology and pathology

In recent years, increasing attention has been paid to titin (TTN) and its mutations. Heterozygous TTN truncating variants (TTNtv) increase the risk of a cardiomyopathy. At the same time, TTNtv and few missense variants have been identified in patients with mainly recessive skeletal muscle diseases. The pathogenic mechanisms underlying titin‐related diseases are still partly unknown. Similarly, the titin mechanical and functional role in the muscle contraction are far from being exhaustively clarified. In the last few years, several animal models carrying variants in the titin gene have been developed and characterized to study the structural and mechanical properties of specific titin domains or to mimic patients'' mutations. This review describes the main animal models so far characterized, including eight mice models and three fish models (Medaka and Zebrafish) and discusses the useful insights provided by a thorough characterization of the cell‐, tissue‐ and organism‐phenotypes in these models.     

OBSCN Mutations Associated with Dilated Cardiomyopathy and Haploinsufficiency

Background

Studies of the functional consequences of DCM-causing mutations have been limited to a few cases where patients with known mutations had heart transplants. To increase the number of potential tissue samples for direct investigation we performed whole exon sequencing of explanted heart muscle samples from 30 patients that had a diagnosis of familial dilated cardiomyopathy and screened for potentially disease-causing mutations in 58 HCM or DCM-related genes.

Results

We identified 5 potentially disease-causing OBSCN mutations in 4 samples; one sample had two OBSCN mutations and one mutation was judged to be not disease-related. Also identified were 6 truncating mutations in TTN, 3 mutations in MYH7, 2 in DSP and one each in TNNC1, TNNI3, MYOM1, VCL, GLA, PLB, TCAP, PKP2 and LAMA4. The mean level of obscurin mRNA was significantly greater and more variable in healthy donor samples than the DCM samples but did not correlate with OBSCN mutations. A single obscurin protein band was observed in human heart myofibrils with apparent mass 960 ± 60 kDa. The three samples with OBSCN mutations had significantly lower levels of obscurin immunoreactive material than DCM samples without OBSCN mutations (45±7, 48±3, and 72±6% of control level).Obscurin levels in DCM controls, donor heart and myectomy samples were the same.

Conclusions

OBSCN mutations may result in the development of a DCM phenotype via haploinsufficiency. Mutations in the obscurin gene should be considered as a significant causal factor of DCM, alone or in concert with other mutations.     

Investigating the Molecular Basis of Retinal Degeneration in a Familial Cohort of Pakistani Decent by Exome Sequencing

Purpose

To define the molecular basis of retinal degeneration in consanguineous Pakistani pedigrees with early onset retinal degeneration.

Methods

A cohort of 277 individuals representing 26 pedigrees from the Punjab province of Pakistan was analyzed. Exomes were captured with commercial kits and sequenced on an Illumina HiSeq 2500. Candidate variants were identified using standard tools and analyzed using exomeSuite to detect all potentially pathogenic changes in genes implicated in retinal degeneration. Segregation analysis was performed by dideoxy sequencing and novel variants were additionally investigated for their presence in ethnicity-matched controls.

Results

We identified a total of nine causal mutations, including six novel variants in RPE65, LCA5, USH2A, CNGB1, FAM161A, CERKL and GUCY2D as the underlying cause of inherited retinal degenerations in 13 of 26 pedigrees. In addition to the causal variants, a total of 200 variants each observed in five or more unrelated pedigrees investigated in this study that were absent from the dbSNP, HapMap, 1000 Genomes, NHLBI ESP6500, and ExAC databases were identified, suggesting that they are common in, and unique to the Pakistani population.

Conclusions

We identified causal mutations associated with retinal degeneration in nearly half of the pedigrees investigated in this study through next generation whole exome sequencing. All novel variants detected in this study through exome sequencing have been cataloged providing a reference database of variants common in, and unique to the Pakistani population.     

Increased Missense Mutation Burden of Fatty Acid Metabolism Related Genes in Nunavik Inuit Population

Background

Nunavik Inuit (northern Quebec, Canada) reside along the arctic coastline where for generations their daily energy intake has mainly been derived from animal fat. Given this particular diet it has been hypothesized that natural selection would lead to population specific allele frequency differences and unique variants in genes related to fatty acid metabolism. A group of genes, namely CPT1A, CPT1B, CPT1C, CPT2, CRAT and CROT, encode for three carnitine acyltransferases that are important for the oxidation of fatty acids, a critical step in their metabolism.

Methods

Exome sequencing and SNP array genotyping were used to examine the genetic variations in the six genes encoding for the carnitine acyltransferases in 113 Nunavik Inuit individuals.

Results

Altogether ten missense variants were found in genes CPT1A, CPT1B, CPT1C, CPT2 and CRAT, including three novel variants and one Inuit specific variant CPT1A p.P479L (rs80356779). The latter has the highest frequency (0.955) compared to other Inuit populations. We found that by comparison to Asians or Europeans, the Nunavik Inuit have an increased mutation burden in CPT1A, CPT2 and CRAT; there is also a high level of population differentiation based on carnitine acyltransferase gene variations between Nunavik Inuit and Asians.

Conclusion

The increased number and frequency of deleterious variants in these fatty acid metabolism genes in Nunavik Inuit may be the result of genetic adaptation to their diet and/or the extremely cold climate. In addition, the identification of these variants may help to understand some of the specific health risks of Nunavik Inuit.     

Whole Exome Sequencing Identifies a Troponin T Mutation Hot Spot in Familial Dilated Cardiomyopathy

Dilated cardiomyopathy (DCM) commonly causes heart failure and shows extensive genetic heterogeneity that may be amenable to newly developed next-generation DNA sequencing of the exome. In this study we report the successful use of exome sequencing to identify a pathogenic variant in the TNNT2 gene using segregation analysis in a large DCM family. Exome sequencing was performed on three distant relatives from a large family with a clear DCM phenotype. Missense, nonsense, and splice variants were analyzed for segregation among the three affected family members and confirmed in other relatives by direct sequencing. A c.517T C>T, Arg173Trp TNNT2 variant segregated with all affected family members and was also detected in one additional DCM family in our registry. The inclusion of segregation analysis using distant family members markedly improved the bioinformatics filtering process by removing from consideration variants that were not shared by all affected subjects. Haplotype analysis confirmed that the variant found in both DCM families was located on two distinct haplotypes, supporting the notion of independent mutational events in each family. In conclusion, an exome sequencing strategy that includes segregation analysis using distant affected relatives within a family represents a viable diagnostic strategy in a genetically heterogeneous disease like DCM.     

Genetic Modifiers of Duchenne Muscular Dystrophy and Dilated Cardiomyopathy

Objective

Dilated cardiomyopathy (DCM) is a major complication and leading cause of death in Duchenne muscular dystrophy (DMD). DCM onset is variable, suggesting modifier effects of genetic or environmental factors. We aimed to determine if polymorphisms previously associated with age at loss of independent ambulation (LoA) in DMD (rs28357094 in the SPP1 promoter, rs10880 and the VTTT/IAAM haplotype in LTBP4) also modify DCM onset.

Methods

A multicentric cohort of 178 DMD patients was genotyped by TaqMan assays. We performed a time-to-event analysis of DCM onset, with age as time variable, and finding of left ventricular ejection fraction < 50% and/or end diastolic volume > 70 mL/m² as event (confirmed by a previous normal exam < 12 months prior); DCM-free patients were censored at the age of last echocardiographic follow-up.

Results

Patients were followed up to an average age of 15.9 ± 6.7 years. Seventy-one/178 patients developed DCM, and median age at onset was 20.0 years. Glucocorticoid corticosteroid treatment (n = 88 untreated; n = 75 treated; n = 15 unknown) did not have a significant independent effect on DCM onset. Cardiological medications were not administered before DCM onset in this population. We observed trends towards a protective effect of the dominant G allele at SPP1 rs28357094 and recessive T allele at LTBP4 rs10880, which was statistically significant in steroid-treated patients for LTBP4 rs10880 (< 50% T/T patients developing DCM during follow-up [n = 13]; median DCM onset 17.6 years for C/C-C/T, log-rank p = 0.027).

Conclusions

We report a putative protective effect of DMD genetic modifiers on the development of cardiac complications, that might aid in risk stratification if confirmed in independent cohorts.     

Identification of Genetic Alterations,as Causative Genetic Defects in Long QT Syndrome,Using Next Generation Sequencing Technology

Background

Long QT Syndrome is an inherited channelopathy leading to sudden cardiac death due to ventricular arrhythmias. Despite that several genes have been associated with the disease, nearly 20% of cases remain without an identified genetic cause. Other genetic alterations such as copy number variations have been recently related to Long QT Syndrome. Our aim was to take advantage of current genetic technologies in a family affected by Long QT Syndrome in order to identify the cause of the disease.

Methods

Complete clinical evaluation was performed in all family members. In the index case, a Next Generation Sequencing custom-built panel, including 55 sudden cardiac death-related genes, was used both for detection of sequence and copy number variants. Next Generation Sequencing variants were confirmed by Sanger method. Copy number variations variants were confirmed by Multiplex Ligation dependent Probe Amplification method and at the mRNA level. Confirmed variants and copy number variations identified in the index case were also analyzed in relatives.

Results

In the index case, Next Generation Sequencing revealed a novel variant in TTN and a large deletion in KCNQ1, involving exons 7 and 8. Both variants were confirmed by alternative techniques. The mother and the brother of the index case were also affected by Long QT Syndrome, and family cosegregation was observed for the KCNQ1 deletion, but not for the TTN variant.

Conclusions

Next Generation Sequencing technology allows a comprehensive genetic analysis of arrhythmogenic diseases. We report a copy number variation identified using Next Generation Sequencing analysis in Long QT Syndrome. Clinical and familiar correlation is crucial to elucidate the role of genetic variants identified to distinguish the pathogenic ones from genetic noise.     

Genome-wide studies of copy number variation and exome sequencing identify rare variants in BAG3 as a cause of dilated cardiomyopathy

Dilated cardiomyopathy commonly causes heart failure and is the most frequent precipitating cause of heart transplantation. Familial dilated cardiomyopathy has been shown to be caused by rare variant mutations in more than 30 genes but only ~35% of its genetic cause has been identified, principally by using linkage-based or candidate gene discovery approaches. In a multigenerational family with autosomal dominant transmission, we employed whole-exome sequencing in a proband and three of his affected family members, and genome-wide copy number variation in the proband and his affected father and unaffected mother. Exome sequencing identified 428 single point variants resulting in missense, nonsense, or splice site changes. Genome-wide copy number analysis identified 51 insertion deletions and 440 copy number variants > 1 kb. Of these, a 8733 bp deletion, encompassing exon 4 of the heat shock protein cochaperone BCL2-associated athanogene 3 (BAG3), was found in seven affected family members and was absent in 355 controls. To establish the relevance of variants in this protein class in genetic DCM, we sequenced the coding exons in BAG3 in 311 other unrelated DCM probands and identified one frameshift, two nonsense, and four missense rare variants absent in 355 control DNAs, four of which were familial and segregated with disease. Knockdown of bag3 in a zebrafish model recapitulated DCM and heart failure. We conclude that new comprehensive genomic approaches have identified rare variants in BAG3 as causative of DCM.     

De Novo Occurrence of a Variant in ARL3 and Apparent Autosomal Dominant Transmission of Retinitis Pigmentosa

Background

Retinitis pigmentosa is a phenotype with diverse genetic causes. Due to this genetic heterogeneity, genome-wide identification and analysis of protein-altering DNA variants by exome sequencing is a powerful tool for novel variant and disease gene discovery. In this study, exome sequencing analysis was used to search for potentially causal DNA variants in a two-generation pedigree with apparent dominant retinitis pigmentosa.

Methods

Variant identification and analysis of three affected members (mother and two affected offspring) was performed via exome sequencing. Parental samples of the index case were used to establish inheritance. Follow-up testing of 94 additional retinitis pigmentosa pedigrees was performed via retrospective analysis or Sanger sequencing.

Results and Conclusions

A total of 136 high quality coding variants in 123 genes were identified which are consistent with autosomal dominant disease. Of these, one of the strongest genetic and functional candidates is a c.269A>G (p.Tyr90Cys) variant in ARL3. Follow-up testing established that this variant occurred de novo in the index case. No additional putative causal variants in ARL3 were identified in the follow-up cohort, suggesting that if ARL3 variants can cause adRP it is an extremely rare phenomenon.     

Evaluation of the Genetic Basis of Familial Aggregation of Pacemaker Implantation by a Large Next Generation Sequencing Panel

Background

The etiology of conduction disturbances necessitating permanent pacemaker (PPM) implantation is often unknown, although familial aggregation of PPM (faPPM) suggests a possible genetic basis. We developed a pan-cardiovascular next generation sequencing (NGS) panel to genetically characterize a selected cohort of faPPM.

Materials and Methods

We designed and validated a custom NGS panel targeting the coding and splicing regions of 246 genes with involvement in cardiac pathogenicity. We enrolled 112 PPM patients and selected nine (8%) with faPPM to be analyzed by NGS.

Results

Our NGS panel covers 95% of the intended target with an average of 229x read depth at a minimum of 15-fold depth, reaching a SNP true positive rate of 98%. The faPPM patients presented with isolated cardiac conduction disease (ICCD) or sick sinus syndrome (SSS) without overt structural heart disease or identifiable secondary etiology. Three patients (33.3%) had heterozygous deleterious variants previously reported in autosomal dominant cardiac diseases including CCD: LDB3 (p.D117N) and TRPM4 (p.G844D) variants in patient 4; TRPM4 (p.G844D) and ABCC9 (p.V734I) variants in patient 6; and SCN5A (p.T220I) and APOB (p.R3527Q) variants in patient 7.

Conclusion

FaPPM occurred in 8% of our PPM clinic population. The employment of massive parallel sequencing for a large selected panel of cardiovascular genes identified a high percentage (33.3%) of the faPPM patients with deleterious variants previously reported in autosomal dominant cardiac diseases, suggesting that genetic variants may play a role in faPPM.     

Low-Frequency IL23R Coding Variant Associated with Crohn’s Disease Susceptibility in Japanese Subjects Identified by Personal Genomics Analysis

Background

The common disease-common variant hypothesis is insufficient to explain the complexities of Crohn’s disease (CD) genetics; therefore, rare variants are expected to be important in the disease. We explored rare variants associated with susceptibility to CD in Japanese individuals by personal genomic analysis.

Methods

Two-step analyses were performed. The first step was a trio analysis with whole-exome sequence (WES) analysis and the second was a follow-up case-control association study. The WES analysis pipeline comprised Burrows-Wheeler Aligner, Picard, Genome Analysis Toolkit, and SAMTOOLS. Single nucleotide variants (SNVs)/indels were annotated and filtered by using programs implemented in ANNOVAR in combination with identity-by-descent (IBD), subsequently were subjected to the linkage based, and de novo based strategies. Finally, we conducted an association study that included 176 unrelated subjects with CD and 358 healthy control subjects.

Results

In family members, 234,067–297,523 SNVs/indels were detected and they were educed to 106–146 by annotation based filtering. Fifty-four CD variants common to both individuals of the affected sib pair were identified. The linkage based strategy detected five candidate variants whereas the de novo based strategy identified no variants. Consequently, five candidates were analyzed in the case-control association study. CD showed a significant association with one variant in exon 4 of IL23R, G149R [rs76418789, P = 3.9E-5, odds ratio (OR) 0.21, 95% confidence interval (CI) 0.09–0.47 for the dominant model (AA + AG versus GG), and P = 7.3E-5, OR 0.21, 95% CI 0.10–0.48 for AG versus GG, and P = 7.2E-5, OR 0.23, 95% CI 0.10–0.50 for the allele model].

Conclusions

The present study, using personal genomics analysis of a small CD pedigree, is the first to show that the low-frequency non-synonymous variant of IL23R, rs76418789, protects against CD development in Japanese subjects.     

Ubiquitin ligases of the N-end rule pathway: assessment of mutations in UBR1 that cause the Johanson-Blizzard syndrome

Background

Johanson-Blizzard syndrome (JBS; OMIM 243800) is an autosomal recessive disorder that includes congenital exocrine pancreatic insufficiency, facial dysmorphism with the characteristic nasal wing hypoplasia, multiple malformations, and frequent mental retardation. Our previous work has shown that JBS is caused by mutations in human UBR1, which encodes one of the E3 ubiquitin ligases of the N-end rule pathway. The N-end rule relates the regulation of the in vivo half-life of a protein to the identity of its N-terminal residue. One class of degradation signals (degrons) recognized by UBR1 are destabilizing N-terminal residues of protein substrates.

Methodology/Principal Findings

Most JBS-causing alterations of UBR1 are nonsense, frameshift or splice-site mutations that abolish UBR1 activity. We report here missense mutations of human UBR1 in patients with milder variants of JBS. These single-residue changes, including a previously reported missense mutation, involve positions in the RING-H2 and UBR domains of UBR1 that are conserved among eukaryotes. Taking advantage of this conservation, we constructed alleles of the yeast Saccharomyces cerevisiae UBR1 that were counterparts of missense JBS-UBR1 alleles. Among these yeast Ubr1 mutants, one of them (H160R) was inactive in yeast-based activity assays, the other one (Q1224E) had a detectable but weak activity, and the third one (V146L) exhibited a decreased but significant activity, in agreement with manifestations of JBS in the corresponding JBS patients.

Conclusions/Significance

These results, made possible by modeling defects of a human ubiquitin ligase in its yeast counterpart, verified and confirmed the relevance of specific missense UBR1 alleles to JBS, and suggested that a residual activity of a missense allele is causally associated with milder variants of JBS.     

Analysis of Genes Involved in Body Weight Regulation by Targeted Re-Sequencing

Introduction

Genes involved in body weight regulation that were previously investigated in genome-wide association studies (GWAS) and in animal models were target-enriched followed by massive parallel next generation sequencing.

Methods

We enriched and re-sequenced continuous genomic regions comprising FTO, MC4R, TMEM18, SDCCAG8, TKNS, MSRA and TBC1D1 in a screening sample of 196 extremely obese children and adolescents with age and sex specific body mass index (BMI) ≥ 99^th percentile and 176 lean adults (BMI ≤ 15^th percentile). 22 variants were confirmed by Sanger sequencing. Genotyping was performed in up to 705 independent obesity trios (extremely obese child and both parents), 243 extremely obese cases and 261 lean adults.

Results and Conclusion

We detected 20 different non-synonymous variants, one frame shift and one nonsense mutation in the 7 continuous genomic regions in study groups of different weight extremes. For SNP Arg695Cys (rs58983546) in TBC1D1 we detected nominal association with obesity (p_TDT = 0.03 in 705 trios). Eleven of the variants were rare, thus were only detected heterozygously in up to ten individual(s) of the complete screening sample of 372 individuals. Two of them (in FTO and MSRA) were found in lean individuals, nine in extremely obese. In silico analyses of the 11 variants did not reveal functional implications for the mutations. Concordant with our hypothesis we detected a rare variant that potentially leads to loss of FTO function in a lean individual. For TBC1D1, in contrary to our hypothesis, the loss of function variant (Arg443Stop) was found in an obese individual. Functional in vitro studies are warranted.     

Variants in the ASB10 Gene Are Associated with Primary Open Angle Glaucoma

Background

Recently nonsynonymous coding variants in the ankyrin repeats and suppressor of cytokine signaling box-containing protein 10 (ASB10) gene were found to be associated with primary open angle glaucoma (POAG) in cohorts from Oregon and Germany, but this finding was not confirmed in an independent cohort from Iowa. The aim of the current study was to assess the role of ASB10 gene variants in Pakistani glaucoma patients.

Methods

Sanger sequencing of the coding exons and splice junctions of the ASB10 gene was performed in 30 probands of multiplex POAG families, 208 sporadic POAG patients and 151 healthy controls from Pakistan. Genotypic associations of individual variants with POAG were analyzed with the Fisher’s exact or Chi-square test.

Results

In total 24 variants were identified in POAG probands and sporadic patients, including 11 novel variants and 13 known variants. 13 of the variants were nonsynonymous, 6 were synonymous, and 5 were intronic. Three nonsynonymous variants (p.Arg49Cys, p.Arg237Gly, p.Arg453Cys) identified in the probands were not segregating in the respective families. This is not surprising since glaucoma is a multifactorial disease, and multiple factors are likely to be involved in the disease manifestation in these families. However a nonsynonymous variant, p.Arg453Cys (rs3800791), was found in 6 sporadic POAG patients but not in controls, suggesting that it infers increased risk for the disease. In addition, one synonymous variant was found to be associated with sporadic POAG: p.Ala290Ala and the association of the variant with POAG remained significant after correction for multiple testing (uncorrected p-value 0.002, corrected p-value 0.047). The cumulative burden of rare, nonsynonymous variants was significantly higher in sporadic POAG patients compared to control individuals (p-value 0.000006).

Conclusions

Variants in ASB10 were found to be significantly associated with sporadic POAG in the Pakistani population. This supports previous findings that sequence variants in the ASB10 gene may act as a risk factor for glaucoma.     

Prevalence and Genetic Characteristics of Geographic Atrophy among Elderly Japanese with Age-Related Macular Degeneration

Objective

To investigate the prevalence and genetic characteristics of geographic atrophy (GA) among elderly Japanese with advanced age-related macular degeneration (AMD) in a clinic-based study.

Methods

Two-hundred and ninety consecutive patients with advanced AMD were classified into typical neovascular AMD, polypoidal choroidal vasculopathy (PCV), retinal angiomatous proliferation (RAP) or geographic atrophy (GA). Genetic variants of ARMS2 A69S (rs10490924) and CFH I62V (rs800292) were genotyped using TaqMan Genotyping Assays. The clinical and genetic characteristics were compared between patients with and without GA.

Results

The number of patients diagnosed as having typical neovascular AMD, PCV, RAP and GA were 98 (33.8%), 151 (52.1%), 22 (7.5%) and 19 (6.6%), respectively. Of 19 patients with GA, 13 patients (68.4%) had unilateral GA with exudative AMD in the contralateral eye. Patients with GA were significantly older, with a higher prevalence of reticular pseudodrusen, bilateral involvement of advanced AMD and T-allele frequency of ARMS2 A69S compared with those with typical AMD and PCV; although there were no differences in the genetic and clinical characteristics among patients with GA and RAP.

Conclusions

The prevalence of GA was 6.6% among elderly Japanese with AMD. Patients with GA and RAP exhibited genetic and clinical similarities.     

Applicability and Efficiency of NGS in Routine Diagnosis: In-Depth Performance Analysis of a Complete Workflow for CFTR Mutation Analysis

Background

Actually, about 2000 sequence variations have been documented in the CFTR gene requiring extensive and multi-step genetic testing in the diagnosis of cystic fibrosis and CFTR-related disorders. We present a two phases study, with validation and performance monitoring, of a single experiment methodology based on multiplex PCR and high throughput sequencing that allows detection of all variants, including large rearrangements, affecting the coding regions plus three deep intronic loci.

Methods

A total of 340 samples, including 257 patients and 83 previously characterized control samples, were sequenced in 17 MiSeq runs and analyzed with two bioinformatic pipelines in routine diagnostic conditions. We obtained 100% coverage for all the target regions in every tested sample.

Results

We correctly identified all the 87 known variants in the control samples and successfully confirmed the 62 variants identified among the patients without observing false positive results. Large rearrangements were identified in 18/18 control samples. Only 17 patient samples showed false positive signals (6.6%), 12 of which showed a borderline result for a single amplicon. We also demonstrated the ability of the assay to detect allele specific dropout of amplicons when a sequence variation occurs at a primer binding site thus limiting the risk for false negative results.

Conclusions

We described here the first NGS workflow for CFTR routine analysis that demonstrated equivalent diagnostic performances compared to Sanger sequencing and multiplex ligation-dependent probe amplification. This study illustrates the advantages of NGS in term of scalability, workload reduction and cost-effectiveness in combination with an improvement of the overall data quality due to the simultaneous detection of SNVs and large rearrangements.     

Analysis of Low Frequency Protein Truncating Stop-Codon Variants and Fasting Concentration of Growth Hormone

Background

The genetic background of Growth Hormone (GH) secretion is not well understood. Mutations giving rise to a stop codon have a high likelihood of affecting protein function.

Objectives

To analyze likely functional stop codon mutations that are associated with fasting plasma concentration of Growth Hormone.

Methods

We analyzed stop codon mutations in 5451 individuals in the Malmö Diet and Cancer study by genotyping the Illumina Exome Chip. To enrich for stop codon mutations with likely functional effects on protein function, we focused on those disrupting >80% of the predicted amino acid sequence, which were carried by ≥10 individuals. Such mutations were related to GH concentration, measured with a high sensitivity assay (hs-GH) and, if nominally significant, to GH related phenotypes, using linear regression analysis.

Results

Two stop codon mutations were associated with the fasting concentration of hs-GH. rs121909305 (NP_005370.1:p.R93*) [Minor Allele Frequency (MAF) = 0.8%] in the Myosin 1A gene (MYO1A) was associated with a 0.36 (95%CI, 0.04 to 0.54; p=0.02) increment of the standardized value of the natural logarithm of hs-GH per 1 minor allele and rs35699176 (NP_067040.1:p.Q100*) in the Zink Finger protein 77 gene (ZNF77) (MAF = 4.8%) was associated with a 0.12 (95%CI, 0.02 to 0.22; p = 0.02) increase of hs-GH. The mutated high hs-GH associated allele of MYO1A was related to lower BMI (β-coefficient, -0.22; p = 0.05), waist (β-coefficient, -0.22; p = 0.04), body fat percentage (β-coefficient, -0.23; p = 0.03) and with higher HDL (β-coefficient, 0.23; p = 0.04). The ZNF77 stop codon was associated with height (β-coefficient, 0.11; p = 0.02) but not with cardiometabolic risk factors.

Conclusion

We here suggest that a stop codon of MYO1A, disrupting 91% of the predicted amino acid sequence, is associated with higher hs-GH and GH-related traits suggesting that MYO1A is involved in GH metabolism and possibly body fat distribution. However, our results are preliminary and need replication in independent populations.     

Genome-wide next-generation DNA and RNA sequencing reveals a mutation that perturbs splicing of the phosphatidylinositol glycan anchor biosynthesis class H gene (PIGH) and causes arthrogryposis in Belgian Blue cattle

Background

Cattle populations are characterized by regular outburst of genetic defects as a result of the extensive use of elite sires. The causative genes and mutations can nowadays be rapidly identified by means of genome-wide association studies combined with next generation DNA sequencing, provided that the causative mutations are conventional loss-of-function variants. We show in this work how the combined use of next generation DNA and RNA sequencing allows for the rapid identification of otherwise difficult to identify splice-site variants.

Results

We report the use of haplotype-based association mapping to identify a locus on bovine chromosome 10 that underlies autosomal recessive arthrogryposis in Belgian Blue Cattle. We identify 31 candidate mutations by resequencing the genome of four cases and 15 controls at ~10-fold depth. By analyzing RNA-Seq data from a carrier fetus, we observe skipping of the second exon of the PIGH gene, which we confirm by RT-PCR to be fully penetrant in tissues from affected calves. We identify - amongst the 31 candidate variants - a C-to-G transversion in the first intron of the PIGH gene (c211-10C > G) that is predicted to affect its acceptor splice-site. The resulting PIGH protein is likely to be non-functional as it lacks essential domains, and hence to cause arthrogryposis.

Conclusions

This work illustrates how the growing arsenal of genome exploration tools continues to accelerate the identification of an even broader range of disease causing mutations, therefore improving the management and control of genetic defects in livestock.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1528-y) contains supplementary material, which is available to authorized users.     

TLR1, 2, 4, 6 and 9 Variants Associated with Tuberculosis Susceptibility: A Systematic Review and Meta-Analysis

Background

Studies investigating the influence of toll-like receptor (TLR) polymorphisms and tuberculosis susceptibility have yielded varying and often contradictory results in different ethnic groups. A meta-analysis was conducted to investigate the relationship between TLR variants and susceptibility to tuberculosis, both across and within specific ethnic groups.

Methods

An extensive database search was performed for studies investigating the relationship between TLR and tuberculosis (TB) susceptibility. Data was subsequently extracted from included studies and statistically analysed.

Results

32 articles involving 18907 individuals were included in this meta-analysis, and data was extracted for 14 TLR polymorphisms. Various genetic models were employed. An increased risk of TB was found for individuals with the TLR2 rs3804100 CC and the TLR9 rs352139 GA and GG genotypes, while decreased risk was identified for those with the AG genotype of TLR1 rs4833095. The T allele of TLR6 rs5743810 conferred protection across all ethnic groups. TLR2 rs5743708 subgroup analysis identified the A allele to increase susceptibility to TB in the Asian ethnic group, while conferring protection in the Hispanic group. The T allele of TLR4 rs4986791 was also found to increase the risk of TB in the Asian subgroup. All other TLR gene variants investigated were not found to be associated with TB in this meta-analysis.

Discussion

Although general associations were identified, most TLR variants showed no significant association with TB, indicating that additional studies investigating a wider range of pattern recognition receptors is required to gain a better understanding of this complex disease.