Next-generation sequencing to genetically diagnose a diverse range of inherited eye disorders in 15 consanguineous families from Pakistan.

Inherited retinal dystrophies (IRDs) are characterized by photoreceptor dysfunction or degeneration. Clinical and phenotypic overlap between IRDs makes the genetic diagnosis very challenging and comprehensive genomic approaches for accurate diagnosis are frequently required. While there are previous studies on IRDs in Pakistan, causative genes and variants are still unknown for a significant portion of patients. Therefore, there is a need to expand the knowledge of the genetic spectrum of IRDs in Pakistan. Here, we recruited 52 affected and 53 normal individuals from 15 consanguineous Pakistani families presenting non-syndromic and syndromic forms of IRDs. We employed single molecule Molecular Inversion Probes (smMIPs) based panel sequencing and whole genome sequencing to identify the probable disease-causing variants in these families. Using this approach, we obtained a 93% genetic solve rate and identified 16 (likely) causative variants in 14 families, of which seven novel variants were identified in ATOH7, COL18A1, MERTK, NDP, PROM1, PRPF8 and USH2A while nine recurrent variants were identified in CNGA3, CNGB1, HGSNAT, NMNAT1, SIX6 and TULP1. The novel MERTK variant and one recurrent TULP1 variant explained the intra-familial locus heterogeneity in one of the screened families while two recurrent CNGA3 variants explained compound heterozygosity in another family. The identification of variants in known disease-associated genes emphasizes the utilization of time and cost-effective screening approaches for rapid diagnosis. The timely genetic diagnosis will not only identify any associated systemic issues in case of syndromic IRDs, but will also aid in the acceleration of personalized medicine for patients affected with IRDs.

Structural Variants Create New Topological-Associated Domains and Ectopic Retinal Enhancer-Gene Contact in Dominant Retinitis Pigmentosa.

The cause of autosomal-dominant retinitis pigmentosa (adRP), which leads to loss of vision and blindness, was investigated in families lacking a molecular diagnosis. A refined locus for adRP on Chr17q22 (RP17) was delineated through genotyping and genome sequencing, leading to the identification of structural variants (SVs) that segregate with disease. Eight different complex SVs were characterized in 22 adRP-affected families with >300 affected individuals. All RP17 SVs had breakpoints within a genomic region spanning YPEL2 to LINC01476. To investigate the mechanism of disease, we reprogrammed fibroblasts from affected individuals and controls into induced pluripotent stem cells (iPSCs) and differentiated them into photoreceptor precursor cells (PPCs) or retinal organoids (ROs). Hi-C was performed on ROs, and differential expression of regional genes and a retinal enhancer RNA at this locus was assessed by qPCR. The epigenetic landscape of the region, and Hi-C RO data, showed that YPEL2 sits within its own topologically associating domain (TAD), rich in enhancers with binding sites for retinal transcription factors. The Hi-C map of RP17 ROs revealed creation of a neo-TAD with ectopic contacts between GDPD1 and retinal enhancers, and modeling of all RP17 SVs was consistent with neo-TADs leading to ectopic retinal-specific enhancer-GDPD1 accessibility. qPCR confirmed increased expression of GDPD1 and increased expression of the retinal enhancer that enters the neo-TAD. Altered TAD structure resulting in increased retinal expression of GDPD1 is the likely convergent mechanism of disease, consistent with a dominant gain of function. Our study highlights the importance of SVs as a genomic mechanism in unsolved Mendelian diseases.

Optical genome mapping and revisiting short-read genome sequencing data reveal previously overlooked structural variants disrupting retinal disease-associated genes.

PURPOSE: Structural variants (SVs) play an important role in inherited retinal diseases (IRD). Although the identification of SVs significantly improved upon the availability of genome sequencing, it is expected that involvement of SVs in IRDs is higher than anticipated. We revisited short-read genome sequencing data to enhance the identification of gene-disruptive SVs. METHODS: Optical genome mapping was performed to improve SV detection in short-read genome sequencing-negative cases. In addition, reanalysis of short-read genome sequencing data was performed to improve the interpretation of SVs and to re-establish SV prioritization criteria. RESULTS: In a monoallelic USH2A case, optical genome mapping identified a pericentric inversion (173 megabase), with 1 breakpoint disrupting USH2A. Retrospectively, the variant could be observed in genome sequencing data but was previously deemed false positive. Reanalysis of short-read genome sequencing data (427 IRD cases) was performed which yielded 30 pathogenic SVs affecting, among other genes, USH2A (n = 15), PRPF31 (n = 3), and EYS (n = 2). Eight of these (>25%) were overlooked during previous analyses. CONCLUSION: Critical evaluation of our findings allowed us to re-establish and improve our SV prioritization and interpretation guidelines, which will prevent missing pathogenic events in future analyses. Our data suggest that more attention should be paid to SV interpretation and the current contribution of SVs in IRDs is still underestimated.

BBS1 branchpoint variant is associated with non-syndromic retinitis pigmentosa.

BACKGROUND: Inherited retinal diseases (IRDs) can be caused by variants in >270 genes. The Bardet-Biedl syndrome 1 (BBS1) gene is one of these genes and may be associated with syndromic and non-syndromic autosomal recessive retinitis pigmentosa (RP). Here, we identified a branchpoint variant in BBS1 and assessed its pathogenicity by in vitro functional analysis. METHODS: Whole genome sequencing was performed for three unrelated monoallelic BBS1 cases with non-syndromic RP. A fourth case received MGCM 105 gene panel analysis. Functional analysis using a midigene splice assay was performed for the putative pathogenic branchpoint variant in BBS1. After confirmation of its pathogenicity, patients were clinically re-evaluated, including assessment of non-ocular features of Bardet-Biedl syndrome. RESULTS: Clinical assessments of probands showed that all individuals displayed non-syndromic RP with macular involvement. Through detailed variant analysis and prioritisation, two pathogenic variants in BBS1, the most common missense variant, c.1169T>G (p.(Met390Arg)), and a branchpoint variant, c.592-21A>T, were identified. Segregation analysis confirmed that in all families, probands were compound heterozygous for c.1169T>G and c.592-21A>T. Functional analysis of the branchpoint variant revealed a complex splicing defect including exon 8 and exon 7/8 skipping, and partial in-frame deletion of exon 8. CONCLUSION: A putative severe branchpoint variant in BBS1, together with a mild missense variant, underlies non-syndromic RP in four unrelated individuals. To our knowledge, this is the first report of a pathogenic branchpoint variant in IRDs that results in a complex splice defect. In addition, this research highlights the importance of the analysis of non-coding regions in order to provide a conclusive molecular diagnosis.

Clinical exome sequencing-Mistakes and caveats.

Massive parallel sequencing technology has become the predominant technique for genetic diagnostics and research. Many genetic laboratories have wrestled with the challenges of setting up genetic testing workflows based on a completely new technology. The learning curve we went through as a laboratory was accompanied by growing pains while we gained new knowledge and expertise. Here we discuss some important mistakes that have been made in our laboratory through 10 years of clinical exome sequencing but that have given us important new insights on how to adapt our working methods. We provide these examples and the lessons that we learned to help other laboratories avoid to make the same mistakes.

Family-based exome sequencing identifies rare coding variants in age-related macular degeneration.

Genome-wide association studies (GWAS) have identified 52 independent variants at 34 genetic loci that are associated with age-related macular degeneration (AMD), the most common cause of incurable vision loss in the elderly worldwide. However, causal genes at the majority of these loci remain unknown. In this study, we performed whole exome sequencing of 264 individuals from 63 multiplex families with AMD and analyzed the data for rare protein-altering variants in candidate target genes at AMD-associated loci. Rare coding variants were identified in the CFH, PUS7, RXFP2, PHF12 and TACC2 genes in three or more families. In addition, we detected rare coding variants in the C9, SPEF2 and BCAR1 genes, which were previously suggested as likely causative genes at respective AMD susceptibility loci. Identification of rare variants in the CFH and C9 genes in our study validated previous reports of rare variants in complement pathway genes in AMD. We then extended our exome-wide analysis and identified rare protein-altering variants in 13 genes outside the AMD-GWAS loci in three or more families. Two of these genes, SCN10A and KIR2DL4, are of interest because variants in these genes also showed association with AMD in case-control cohorts, albeit not at the level of genome-wide significance. Our study presents the first large-scale, exome-wide analysis of rare variants in AMD. Further independent replications and molecular investigation of candidate target genes, reported here, would assist in gaining novel insights into mechanisms underlying AMD pathogenesis.

Systematic analysis of short tandem repeats in 38,095 exomes provides an additional diagnostic yield.

PURPOSE: Expansions of a subset of short tandem repeats (STRs) have been implicated in approximately 30 different human genetic disorders. Despite extensive application of exome sequencing (ES) in routine diagnostic genetic testing, STRs are not routinely identified from these data. METHODS: We assessed diagnostic utility of STR analysis in exome sequencing by applying ExpansionHunter to 2,867 exomes from movement disorder patients and 35,228 other clinical exomes. RESULTS: We identified 38 movement disorder patients with a possible aberrant STR length. Validation by polymerase chain reaction (PCR) and/or repeat-primed PCR technologies confirmed the presence of aberrant expansion alleles for 13 (34%). For seven of these patients the genotype was compatible with the phenotypic description, resulting in a molecular diagnosis. We subsequently tested the remainder of our diagnostic ES cohort, including over 30 clinically and genetically heterogeneous disorders. Optimized manual curation yielded 167 samples with a likely aberrant STR length. Validations confirmed 93/167 (56%) aberrant expansion alleles, of which 48 were in the pathogenic range and 45 in the premutation range. CONCLUSION: Our work provides guidance for the implementation of STR analysis in clinical ES. Our results show that systematic STR evaluation may increase diagnostic ES yield by 0.2%, and recommend making STR evaluation a routine part of ES interpretation in genetic testing laboratories.

Identification of ADHD risk genes in extended pedigrees by combining linkage analysis and whole-exome sequencing.

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (Ntotal = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.

Sudden death in epilepsy and ectopic neurohypophysis in Joubert syndrome 23 diagnosed using SNVs/indels and structural variants pipelines on WGS data: a case report.

BACKGROUND: Joubert syndrome (JBTS) is a genetically heterogeneous group of neurodevelopmental syndromes caused by primary cilia dysfunction. Usually the neurological presentation starts with abnormal neonatal breathing followed by muscular hypotonia, psychomotor delay, and cerebellar ataxia. Cerebral MRI shows mid- and hindbrain anomalies including the molar tooth sign. We report a male patient with atypical presentation of Joubert syndrome type 23, thus expanding the phenotype. CASE PRESENTATION: Clinical features were consistent with JBTS already from infancy, yet the syndrome was not suspected before cerebral MRI later in childhood showed the characteristic molar tooth sign and ectopic neurohypophysis. From age 11 years seizures developed and after few years became increasingly difficult to treat, also related to inadequate compliance to therapy. He died at 23 years of sudden unexpected death in epilepsy (SUDEP). The genetic diagnosis remained elusive for many years, despite extensive genetic testing. We reached the genetic diagnosis by performing whole genome sequencing of the family trio and analyzing the data with the combination of one analysis pipeline for single nucleotide variants (SNVs)/indels and one for structural variants (SVs). This lead to the identification of the most common variant detected in patients with JBTS23 (OMIM# 616490), rs534542684, in compound heterozygosity with a 8.3 kb deletion in KIAA0586, not previously reported. CONCLUSIONS: We describe for the first time ectopic neurohypophysis and SUDEP in JBTS23, expanding the phenotype of this condition and raising the attention on the possible severity of the epilepsy in this disease. We also highlight the diagnostic power of WGS, which efficiently detects SNVs/indels and in addition allows the identification of SVs.

Genome-Wide Association Study Reveals Variants in CFH and CFHR4 Associated with Systemic Complement Activation: Implications in Age-Related Macular Degeneration.

PURPOSE: To identify genetic variants associated with complement activation, which may help to select age-related macular degeneration (AMD) patients for complement-inhibiting therapies. DESIGN: Genome-wide association study (GWAS) followed by replication and meta-analysis. PARTICIPANTS: AMD patients and controls (n = 2245). METHODS: A GWAS on serum C3d-to-C3 ratio was performed in 1548 AMD patients and controls. For replication and meta-analysis, 697 additional individuals were genotyped. A model for complement activation including genetic and non-genetic factors was built, and the variance explained was estimated. Haplotype analysis was performed for 8 SNPs across the CFH/CFHR locus. Association with AMD was performed for the variants and haplotypes found to influence complement activation. MAIN OUTCOME MEASURES: Normalized C3d/C3 ratio as a measure of systemic complement activation. RESULTS: Complement activation was associated independently with rs3753396 located in CFH (Pdiscovery = 1.09 × 10-15; Pmeta = 3.66 × 10-21; ß = 0.141; standard error [SE] = 0.015) and rs6685931 located in CFHR4 (Pdiscovery = 8.18 × 10-7; Pmeta = 6.32 × 10-8; ß = 0.054; SE = 0.010). A model including age, AMD disease status, body mass index, triglycerides, rs3753396, rs6685931, and previously identified SNPs explained 18.7% of the variability in complement activation. Haplotype analysis revealed 3 haplotypes (H1-2 and H6 containing rs6685931 and H3 containing rs3753396) associated with complement activation. Haplotypes H3 and H6 conferred stronger effects on complement activation compared with the single variants (P = 2.53 × 10-14; ß = 0.183; SE = 0.024; and P = 4.28 × 10-4; ß = 0.144; SE = 0.041; respectively). Association analyses with AMD revealed that SNP rs6685931 and haplotype H1-2 containing rs6685931 were associated with a risk for AMD development, whereas SNP rs3753396 and haplotypes H3 and H6 were not. CONCLUSIONS: The SNP rs3753396 in CFH and SNP rs6685931 in CFHR4 are associated with systemic complement activation levels. The SNP rs6685931 in CFHR4 and its linked haplotype H1-2 also conferred a risk for AMD development, and therefore could be used to identify AMD patients who would benefit most from complement-inhibiting therapies.

Whole-Exome Sequencing in Age-Related Macular Degeneration Identifies Rare Variants in COL8A1, a Component of Bruch's Membrane.

PURPOSE: Genome-wide association studies and targeted sequencing studies of candidate genes have identified common and rare variants that are associated with age-related macular degeneration (AMD). Whole-exome sequencing (WES) studies allow a more comprehensive analysis of rare coding variants across all genes of the genome and will contribute to a better understanding of the underlying disease mechanisms. To date, the number of WES studies in AMD case-control cohorts remains scarce and sample sizes are limited. To scrutinize the role of rare protein-altering variants in AMD cause, we performed the largest WES study in AMD to date in a large European cohort consisting of 1125 AMD patients and 1361 control participants. DESIGN: Genome-wide case-control association study of WES data. PARTICIPANTS: One thousand one hundred twenty-five AMD patients and 1361 control participants. METHODS: A single variant association test of WES data was performed to detect variants that are associated individually with AMD. The cumulative effect of multiple rare variants with 1 gene was analyzed using a gene-based CMC burden test. Immunohistochemistry was performed to determine the localization of the Col8a1 protein in mouse eyes. MAIN OUTCOME MEASURES: Genetic variants associated with AMD. RESULTS: We detected significantly more rare protein-altering variants in the COL8A1 gene in patients (22/2250 alleles [1.0%]) than in control participants (11/2722 alleles [0.4%]; P = 7.07×10-5). The association of rare variants in the COL8A1 gene is independent of the common intergenic variant (rs140647181) near the COL8A1 gene previously associated with AMD. We demonstrated that the Col8a1 protein localizes at Bruch's membrane. CONCLUSIONS: This study supported a role for protein-altering variants in the COL8A1 gene in AMD pathogenesis. We demonstrated the presence of Col8a1 in Bruch's membrane, further supporting the role of COL8A1 variants in AMD pathogenesis. Protein-altering variants in COL8A1 may alter the integrity of Bruch's membrane, contributing to the accumulation of drusen and the development of AMD.

Genetic screening for macular dystrophies in patients clinically diagnosed with dry age-related macular degeneration.

It can be clinically challenging to distinguish dry age-related macular degeneration (AMD) from AMD-mimicking dystrophies, and sometimes misdiagnosis occurs. With upcoming therapies for dry AMD it is important to exclude patients with a different retinal disease from clinical trials. In this study we evaluated the occurrence of AMD-mimicking dystrophies in an AMD cohort. Whole-exome sequencing (WES) was performed in 218 patients with intermediate AMD or geographic atrophy secondary to AMD and 133 control individuals. WES data was analyzed for rare variants in 19 genes associated with autosomal dominant and recessive macular dystrophies mimicking AMD. In three (1.4%) of 218 cases we identified a pathogenic heterozygous variant (PRPH2 c.424C > T; p.R142W) causal for autosomal dominant central areolar choroidal dystrophy (CACD). Phenotypically, these patients all presented with geographic atrophy. In 12 (5.5%) of 218 cases we identified a heterozygous variant of unknown clinical significance, but predicted to be highly deleterious, in genes previously associated with autosomal dominant macular dystrophies. The distinction between AMD and AMD-mimicking dystrophies, such as CACD, can be challenging based on fundus examination alone. Genetic screening for genes associated with macular dystrophies, especially PRPH2, can be beneficial to help identify AMD-mimicking dystrophies.

Epigenetic regulation of the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs.

BACKGROUND: In previous studies on an Iberian x Landrace cross, we have provided evidence that supported the porcine ELOVL6 gene as the major causative gene of the QTL on pig chromosome 8 for palmitic and palmitoleic acid contents in muscle and backfat. The single nucleotide polymorphism (SNP) ELOVL6:c.-533C > T located in the promoter region of ELOVL6 was found to be highly associated with ELOVL6 expression and, accordingly, with the percentages of palmitic and palmitoleic acids in longissimus dorsi and adipose tissue. The main goal of the current work was to further study the role of ELOVL6 on these traits by analyzing the regulation of the expression of ELOVL6 and the implication of ELOVL6 polymorphisms on meat quality traits in pigs. RESULTS: High-throughput sequencing of BAC clones that contain the porcine ELOVL6 gene coupled to RNAseq data re-analysis showed that two isoforms of this gene are expressed in liver and adipose tissue and that they differ in number of exons and 3'UTR length. Although several SNPs in the 3'UTR of ELOVL6 were associated with palmitic and palmitoleic acid contents, this association was lower than that previously observed with SNP ELOVL6:c.-533C > T. This SNP is in full linkage disequilibrium with SNP ELOVL6:c.-394G > A that was identified in the binding site for estrogen receptor alpha (ERα). Interestingly, the ELOVL6:c.-394G allele is associated with an increase in methylation levels of the ELOVL6 promoter and with a decrease of ELOVL6 expression. Therefore, ERα is clearly a good candidate to explain the regulation of ELOVL6 expression through dynamic epigenetic changes in the binding site of known regulators of ELOVL6 gene, such as SREBF1 and SP1. CONCLUSIONS: Our results strongly suggest the ELOVL6:c.-394G > A polymorphism as the causal mutation for the QTL on pig chromosome 8 that affects fatty acid composition in pigs.

New insight into the SSC8 genetic determination of fatty acid composition in pigs.

BACKGROUND: Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F2 intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C < T), which is associated with percentage of palmitic and palmitoleic acids in muscle and adipose tissues. Here, a combination of single-marker association and the haplotype-based approach was used to analyze backfat fatty acid composition in 470 animals of an Iberian x Landrace F2 intercross genotyped with 144 SNPs (single nucleotide polymorphisms) distributed along SSC8. RESULTS: Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat. CONCLUSIONS: Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region.

Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition.

BACKGROUND: In pigs, adipose tissue is one of the principal organs involved in the regulation of lipid metabolism. It is particularly involved in the overall fatty acid synthesis with consequences in other lipid-target organs such as muscles and the liver. With this in mind, we have used massive, parallel high-throughput sequencing technologies to characterize the porcine adipose tissue transcriptome architecture in six Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition (three per group). RESULTS: High-throughput RNA sequencing was used to generate a whole characterization of adipose tissue (backfat) transcriptome. A total of 4,130 putative unannotated protein-coding sequences were identified in the 20% of reads which mapped in intergenic regions. Furthermore, 36% of the unmapped reads were represented by interspersed repeats, SINEs being the most abundant elements. Differential expression analyses identified 396 candidate genes among divergent animals for intramuscular fatty acid composition. Sixty-two percent of these genes (247/396) presented higher expression in the group of pigs with higher content of intramuscular SFA and MUFA, while the remaining 149 showed higher expression in the group with higher content of PUFA. Pathway analysis related these genes to biological functions and canonical pathways controlling lipid and fatty acid metabolisms. In concordance with the phenotypic classification of animals, the major metabolic pathway differentially modulated between groups was de novo lipogenesis, the group with more PUFA being the one that showed lower expression of lipogenic genes. CONCLUSIONS: These results will help in the identification of genetic variants at loci that affect fatty acid composition traits. The implications of these results range from the improvement of porcine meat quality traits to the application of the pig as an animal model of human metabolic diseases.

Nutrition and Hydration for High-Altitude Alpinism: A Narrative Review.

This report aims to summarise the scientific knowledge around hydration, nutrition, and metabolism at high altitudes and to transfer it into the practical context of extreme altitude alpinism, which, as far as we know, has never been considered before in the literature. Maintaining energy balance during alpine expeditions is difficult for several reasons and requires a deep understanding of human physiology and the biological basis for altitude acclimation. However, in these harsh conditions it is difficult to reconcile our current scientific knowledge in sports nutrition or even for mountaineering to high-altitude alpinism: extreme hypoxia, cold, and the logistical difficulties intrinsic to these kinds of expeditions are not considered in the current literature. Requirements for the different stages of an expedition vary dramatically with increasing altitude, so recommendations must differentiate whether the alpinist is at base camp, at high-altitude camps, or attempting the summit. This paper highlights nutritional recommendations regarding prioritising carbohydrates as a source of energy and trying to maintain a protein balance with a practical contextualisation in the extreme altitude environment in the different stages of an alpine expedition. More research is needed regarding specific macro and micronutrient requirements as well as the adequacy of nutritional supplementations at high altitudes.

Combined Single Gene Testing and Genome Sequencing as an Effective Diagnostic Approach for Anophthalmia and Microphthalmia Patients.

Anophthalmia and microphthalmia (A/M) are among the most severe congenital developmental eye disorders. Despite the advancements in genome screening technologies, more than half of A/M patients do not receive a molecular diagnosis. We included seven consanguineous families affected with A/M from Pakistani cohort and an unknown molecular basis. Single gene testing of FOXE3 was performed, followed by genome sequencing for unsolved probands in order to establish a genetic diagnosis for these families. All seven families were provided with a genetic diagnosis. The identified variants were all homozygous, classified as (likely) pathogenic and present in an A/M-associated gene. Targeted FOXE3 sequencing revealed two previously reported pathogenic FOXE3 variants in four families. In the remaining families, genome sequencing revealed a known pathogenic PXDN variant, a novel 13bp deletion in VSX2, and one novel deep intronic splice variant in PXDN. An in vitro splice assay was performed for the PXDN splice variant which revealed a severe splicing defect. Our study confirmed the utility of genome sequencing as a diagnostic tool for A/M-affected individuals. Furthermore, the identification of a novel deep intronic pathogenic variant in PXDN highlights the role of non-coding variants in A/M-disorders and the value of genome sequencing for the identification of this type of variants.

ABCA4 c.6480-35A>G, a novel branchpoint variant associated with Stargardt disease.

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 (ABCA4) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone-rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants. Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro. The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48. Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699_768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans. For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases. Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis.

Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition.

BACKGROUND: New advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq. RESULTS: The liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism. CONCLUSIONS: In the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.

Diagnostic yield of next-generation sequencing in 87 families with neurodevelopmental disorders.

BACKGROUND: Neurodevelopmental disorders (NDDs) are a group of heterogeneous conditions, which include mainly intellectual disability, developmental delay (DD) and autism spectrum disorder (ASD), among others. These diseases are highly heterogeneous and both genetic and environmental factors play an important role in many of them. The introduction of next generation sequencing (NGS) has lead to the detection of genetic variants in several genetic diseases. The main aim of this report is to discuss the impact and advantages of the implementation of NGS in the diagnosis of NDDs. Herein, we report diagnostic yields of applying whole exome sequencing in 87 families affected by NDDs and additional data of whole genome sequencing (WGS) from 12 of these families. RESULTS: The use of NGS technologies allowed identifying the causative gene alteration in approximately 36% (31/87) of the families. Among them, de novo mutation represented the most common cause of genetic alteration found in 48% (15/31) of the patients with diagnostic mutations. The majority of variants were located in known neurodevelopmental disorders genes. Nevertheless, some of the diagnoses were made after the use of GeneMatcher tools which allow the identification of additional patients carrying mutations in THOC2, SETD1B and CHD9 genes. Finally the use of WGS only allowed the identification of disease causing variants in 8% (1/12) of the patients in which previous WES failed to identify a genetic aetiology. CONCLUSION: NGS is more powerful in identifying causative pathogenic variant than conventional algorithms based on chromosomal microarray as first-tier test. Our results reinforce the implementation of NGS as a first-test in genetic diagnosis of NDDs.