Optimizing freeform lenses for extended sources with algorithmic differentiable ray tracing and truncated hierarchical B-splines.

We propose a method for optimizing the geometry of a freeform lens to redirect the light emitted from an extended source into a desired irradiance distribution. We utilize a gradient-based optimization approach with MITSUBA 3, an algorithmic differentiable non-sequential ray tracer that allows us to obtain the gradients of the freeform surface parameters with respect to the produced irradiance distribution. To prevent the optimizer from getting trapped in local minima, we gradually increase the number of degrees of freedom of the surface by using Truncated Hierarchical B-splines (THB-splines) during optimization. The refinement locations are determined by analyzing the gradients of the surface vertices. We first design a freeform using a collimated beam (zero-etendue source) for a complex target distribution to demonstrate the method's effectiveness. Then, we demonstrate the ability of this approach to create a freeform that can project the light of an extended Lambertian source into a prescribed target distribution.

Clustering of Cardiac Transcriptome Profiles Reveals Unique: Subgroups of Dilated Cardiomyopathy Patients.

Dilated cardiomyopathy is a heterogeneous disease characterized by multiple genetic and environmental etiologies. The majority of patients are treated the same despite these differences. The cardiac transcriptome provides information on the patient's pathophysiology, which allows targeted therapy. Using clustering techniques on data from the genotype, phenotype, and cardiac transcriptome of patients with early- and end-stage dilated cardiomyopathy, more homogeneous patient subgroups are identified based on shared underlying pathophysiology. Distinct patient subgroups are identified based on differences in protein quality control, cardiac metabolism, cardiomyocyte function, and inflammatory pathways. The identified pathways have the potential to guide future treatment and individualize patient care.

Cost-effective sequence analysis of 113 genes in 1,192 probands with retinitis pigmentosa and Leber congenital amaurosis.

Introduction: Retinitis pigmentosa (RP) and Leber congenital amaurosis (LCA) are two groups of inherited retinal diseases (IRDs) where the rod photoreceptors degenerate followed by the cone photoreceptors of the retina. A genetic diagnosis for IRDs is challenging since >280 genes are associated with these conditions. While whole exome sequencing (WES) is commonly used by diagnostic facilities, the costs and required infrastructure prevent its global applicability. Previous studies have shown the cost-effectiveness of sequence analysis using single molecule Molecular Inversion Probes (smMIPs) in a cohort of patients diagnosed with Stargardt disease and other maculopathies. Methods: Here, we introduce a smMIPs panel that targets the exons and splice sites of all currently known genes associated with RP and LCA, the entire RPE65 gene, known causative deep-intronic variants leading to pseudo-exons, and part of the RP17 region associated with autosomal dominant RP, by using a total of 16,812 smMIPs. The RP-LCA smMIPs panel was used to screen 1,192 probands from an international cohort of predominantly RP and LCA cases. Results and discussion: After genetic analysis, a diagnostic yield of 56% was obtained which is on par with results from WES analysis. The effectiveness and the reduced costs compared to WES renders the RP-LCA smMIPs panel a competitive approach to provide IRD patients with a genetic diagnosis, especially in countries with restricted access to genetic testing.

Correction: Evaluation of molecular inversion probe versus TruSeq® custom methods for targeted next-generation sequencing.

[This corrects the article DOI: 10.1371/journal.pone.0238467.].

Diagnostic exome-based preconception carrier testing in consanguineous couples: results from the first 100 couples in clinical practice.

PURPOSE: Consanguineous couples are at increased risk of being heterozygous for the same autosomal recessive (AR) disorder(s), with a 25% risk of affected offspring as a consequence. Until recently, comprehensive preconception carrier testing (PCT) for AR disorders was unavailable in routine diagnostics. Here we developed and implemented such a test in routine clinical care. METHODS: We performed exome sequencing (ES) for 100 consanguineous couples. For each couple, rare variants that could give rise to biallelic variants in offspring were selected. These variants were subsequently filtered against a gene panel consisting of ~2,000 genes associated with known AR disorders (OMIM-based). Remaining variants were classified according to American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines, after which only likely pathogenic and pathogenic (class IV/V) variants, present in both partners, were reported. RESULTS: In 28 of 100 tested consanguineous couples (28%), likely pathogenic and pathogenic variants not previously known in the couple or their family were reported conferring 25% risk of affected offspring. CONCLUSION: ES-based PCT provides a powerful diagnostic tool to identify AR disease carrier status in consanguineous couples. Outcomes provided significant reproductive choices for a higher proportion of these couples than previous tests.

Dutch genome diagnostic laboratories accelerated and improved variant interpretation and increased accuracy by sharing data.

Each year diagnostic laboratories in the Netherlands profile thousands of individuals for heritable disease using next-generation sequencing (NGS). This requires pathogenicity classification of millions of DNA variants on the standard 5-tier scale. To reduce time spent on data interpretation and increase data quality and reliability, the nine Dutch labs decided to publicly share their classifications. Variant classifications of nearly 100,000 unique variants were catalogued and compared in a centralized MOLGENIS database. Variants classified by more than one center were labeled as "consensus" when classifications agreed, and shared internationally with LOVD and ClinVar. When classifications opposed (LB/B vs. LP/P), they were labeled "conflicting", while other nonconsensus observations were labeled "no consensus". We assessed our classifications using the InterVar software to compare to ACMG 2015 guidelines, showing 99.7% overall consistency with only 0.3% discrepancies. Differences in classifications between Dutch labs or between Dutch labs and ACMG were mainly present in genes with low penetrance or for late onset disorders and highlight limitations of the current 5-tier classification system. The data sharing boosted the quality of DNA diagnostics in Dutch labs, an initiative we hope will be followed internationally. Recently, a positive match with a case from outside our consortium resulted in a more definite disease diagnosis.

Whole Exome Sequencing Is the Preferred Strategy to Identify the Genetic Defect in Patients With a Probable or Possible Mitochondrial Cause.

Mitochondrial disorders, characterized by clinical symptoms and/or OXPHOS deficiencies, are caused by pathogenic variants in mitochondrial genes. However, pathogenic variants in some of these genes can lead to clinical manifestations which overlap with other neuromuscular diseases, which can be caused by pathogenic variants in non-mitochondrial genes as well. Mitochondrial pathogenic variants can be found in the mitochondrial DNA (mtDNA) or in any of the 1,500 nuclear genes with a mitochondrial function. We have performed a two-step next-generation sequencing approach in a cohort of 117 patients, mostly children, in whom a mitochondrial disease-cause could likely or possibly explain the phenotype. A total of 86 patients had a mitochondrial disorder, according to established clinical and biochemical criteria. The other 31 patients had neuromuscular symptoms, where in a minority a mitochondrial genetic cause is present, but a non-mitochondrial genetic cause is more likely. All patients were screened for pathogenic variants in the mtDNA and, if excluded, analyzed by whole exome sequencing (WES). Variants were filtered for being pathogenic and compatible with an autosomal or X-linked recessive mode of inheritance in families with multiple affected siblings and/or consanguineous parents. Non-consanguineous families with a single patient were additionally screened for autosomal and X-linked dominant mutations in a predefined gene-set. We identified causative pathogenic variants in the mtDNA in 20% of the patient-cohort, and in nuclear genes in 49%, implying an overall yield of 68%. We identified pathogenic variants in mitochondrial and non-mitochondrial genes in both groups with, obviously, a higher number of mitochondrial genes affected in mitochondrial disease patients. Furthermore, we show that 31% of the disease-causing genes in the mitochondrial patient group were not included in the MitoCarta database, and therefore would have been missed with MitoCarta based gene-panels. We conclude that WES is preferable to panel-based approaches for both groups of patients, as the mitochondrial gene-list is not complete and mitochondrial symptoms can be secondary. Also, clinically and genetically heterogeneous disorders would require sequential use of multiple different gene panels. We conclude that WES is a comprehensive and unbiased approach to establish a genetic diagnosis in these patients, able to resolve multi-genic disease-causes.

Genetic defects in mtDNA-encoded protein translation cause pediatric, mitochondrial cardiomyopathy with early-onset brain disease.

This study aims to identify gene defects in pediatric cardiomyopathy and early-onset brain disease with oxidative phosphorylation (OXPHOS) deficiencies. We applied whole-exome sequencing in three patients with pediatric cardiomyopathy and early-onset brain disease with OXPHOS deficiencies. The brain pathology was studied by MRI analysis. In consanguineous patient 1, we identified a homozygous intronic variant (c.850-3A > G) in the QRSL1 gene, which was predicted to cause abnormal splicing. The variant segregated with the disease and affected the protein function, which was confirmed by complementation studies, restoring OXPHOS function only with wild-type QRSL1. Patient 2 was compound heterozygous for two novel affected and disease-causing variants (c.[253G > A];[938G > A]) in the MTO1 gene. In patient 3, we detected one unknown affected and disease-causing variants (c.2872C > T) and one known disease-causing variant (c.1774C > T) in the AARS2 gene. The c.1774C > T variant was present in the paternal copy of the AARS2 gene, the c.2872C > T in the maternal copy. All genes were involved in translation of mtDNA-encoded proteins. Defects in mtDNA-encoded protein translation lead to severe pediatric cardiomyopathy and brain disease with OXPHOS abnormalities. This suggests that the heart and brain are particularly sensitive to defects in mitochondrial protein synthesis during late embryonic or early postnatal development, probably due to the massive mitochondrial biogenesis occurring at that stage. If both the heart and brain are involved, the prognosis is poor with a likely fatal outcome at young age.

A comprehensive strategy for exome-based preconception carrier screening.

PURPOSE: Whole-exome sequencing (WES) provides the possibility of genome-wide preconception carrier screening (PCS). Here, we propose a filter strategy to rapidly identify the majority of relevant pathogenic mutations. METHODS: Our strategy was developed using WES data from eight consanguineous and five fictive nonconsanguineous couples and was subsequently applied to 20 other fictive nonconsanguineous couples. Presumably pathogenic variants based on frequency and database annotations or generic characteristics and mutation type were selected in genes shared by the couple and in the female's X-chromosome. Unclassified variants were not included. RESULTS: This yielded an average of 29 (19-51) variants in genes shared by the consanguineous couples and 15 (6-30) shared by the nonconsanguineous couples. For X-linked variants, the numbers per female were 3 (1-5) and 1 (0-3), respectively. Remaining variants were verified manually. The majority were able to be quickly discarded, effectively leaving true pathogenic variants. CONCLUSION: We conclude that WES is applicable for PCS, both for consanguineous and nonconsanguineous couples, with the remaining number of variants being manageable in a clinical setting. The addition of gene panels for filtering was not favorable because it resulted in missing pathogenic variants. It is important to develop and continuously curate databases with pathogenic mutations to further increase the sensitivity of WES-based PCS.Genet Med advance online publication 27 October 2016.

Rapid Resolution of Blended or Composite Multigenic Disease in Infants by Whole-Exome Sequencing.

Whole-exome sequencing identified multiple genetic causes in 2 infants with heterogeneous disease. Three gene defects in the first patient explained all symptoms, but manifestations were overlapping (blended phenotype). Two gene defects in the second patient explained nonoverlapping symptoms (composite phenotype). Whole-exome sequencing rapidly and comprehensively resolves heterogeneous genetic disease.

Specific MRI Abnormalities Reveal Severe Perrault Syndrome due to CLPP Defects.

In establishing a genetic diagnosis in heterogeneous neurological disease, clinical characterization and whole exome sequencing (WES) go hand-in-hand. Clinical data are essential, not only to guide WES variant selection and define the clinical severity of a genetic defect but also to identify other patients with defects in the same gene. In an infant patient with sensorineural hearing loss, psychomotor retardation, and epilepsy, WES resulted in identification of a novel homozygous CLPP frameshift mutation (c.21delA). Based on the gene defect and clinical symptoms, the diagnosis Perrault syndrome type 3 (PRLTS3) was established. The patient's brain-MRI revealed specific abnormalities of the subcortical and deep cerebral white matter and the middle blade of the corpus callosum, which was used to identify similar patients in the Amsterdam brain-MRI database, containing over 3000 unclassified leukoencephalopathy cases. In three unrelated patients with similar MRI abnormalities the CLPP gene was sequenced, and in two of them novel missense mutations were identified together with a large deletion that covered part of the CLPP gene on the other allele. The severe neurological and MRI abnormalities in these young patients were due to the drastic impact of the CLPP mutations, correlating with the variation in clinical manifestations among previously reported patients. Our data show that similarity in brain-MRI patterns can be used to identify novel PRLTS3 patients, especially during early disease stages, when only part of the disease manifestations are present. This seems especially applicable to the severely affected cases in which CLPP function is drastically affected and MRI abnormalities are pronounced.

Pathogenic CWF19L1 variants as a novel cause of autosomal recessive cerebellar ataxia and atrophy.

Autosomal recessive cerebellar ataxia (ARCA) is a group of neurological disorders characterized by degeneration or abnormal development of the cerebellum and spinal cord. ARCA is clinically and genetically highly heterogeneous, with over 20 genes involved. Exome sequencing of a girl with ARCA from non-consanguineous Dutch parents revealed two pathogenic variants c.37G>C; p.D13H and c.946A>T; p.K316* in CWF19L1, a gene with an unknown function, recently reported to cause ARCA in a Turkish family. Sanger sequencing showed that the c.37G>C variant was inherited from the father and the c.946A>T variant from the mother. Pathogenicity was based on the damaging effect on protein function as the c.37G>C variant changed the highly conserved, negatively charged aspartic acid to the positively charged histidine and the c.946A>T variant introduced a premature stop codon. In addition, 27 patients with ARCA were tested for pathogenic variants in CWF19L1, however, no pathogenic variants were identified. Our data confirm CWF19L1 as a novel but rare gene causing ARCA.

Germ-line variants identified by next generation sequencing in a panel of estrogen and cancer associated genes correlate with poor clinical outcome in Lynch syndrome patients.

BACKGROUND: The risk to develop colorectal and endometrial cancers among subjects testing positive for a pathogenic Lynch syndrome mutation varies, making the risk prediction difficult. Genetic risk modifiers alter the risk conferred by inherited Lynch syndrome mutations, and their identification can improve genetic counseling. We aimed at identifying rare genetic modifiers of the risk of Lynch syndrome endometrial cancer. METHODS: A family based approach was used to assess the presence of genetic risk modifiers among 35 Lynch syndrome mutation carriers having either a poor clinical phenotype (early age of endometrial cancer diagnosis or multiple cancers) or a neutral clinical phenotype. Putative genetic risk modifiers were identified by Next Generation Sequencing among a panel of 154 genes involved in endometrial physiology and carcinogenesis. RESULTS: A simple pipeline, based on an allele frequency lower than 0.001 and on predicted non-conservative amino-acid substitutions returned 54 variants that were considered putative risk modifiers. The presence of two or more risk modifying variants in women carrying a pathogenic Lynch syndrome mutation was associated with a poor clinical phenotype. CONCLUSION: A gene-panel is proposed that comprehends genes that can carry variants with putative modifying effects on the risk of Lynch syndrome endometrial cancer. Validation in further studies is warranted before considering the possible use of this tool in genetic counseling.

Next-generation sequencing-based genome diagnostics across clinical genetics centers: implementation choices and their effects.

Implementation of next-generation DNA sequencing (NGS) technology into routine diagnostic genome care requires strategic choices. Instead of theoretical discussions on the consequences of such choices, we compared NGS-based diagnostic practices in eight clinical genetic centers in the Netherlands, based on genetic testing of nine pre-selected patients with cardiomyopathy. We highlight critical implementation choices, including the specific contributions of laboratory and medical specialists, bioinformaticians and researchers to diagnostic genome care, and how these affect interpretation and reporting of variants. Reported pathogenic mutations were consistent for all but one patient. Of the two centers that were inconsistent in their diagnosis, one reported to have found 'no causal variant', thereby underdiagnosing this patient. The other provided an alternative diagnosis, identifying another variant as causal than the other centers. Ethical and legal analysis showed that informed consent procedures in all centers were generally adequate for diagnostic NGS applications that target a limited set of genes, but not for exome- and genome-based diagnosis. We propose changes to further improve and align these procedures, taking into account the blurring boundary between diagnostics and research, and specific counseling options for exome- and genome-based diagnostics. We conclude that alternative diagnoses may infer a certain level of 'greediness' to come to a positive diagnosis in interpreting sequencing results. Moreover, there is an increasing interdependence of clinic, diagnostics and research departments for comprehensive diagnostic genome care. Therefore, we invite clinical geneticists, physicians, researchers, bioinformatics experts and patients to reconsider their role and position in future diagnostic genome care.

Archaeal MBF1 binds to 30S and 70S ribosomes via its helix-turn-helix domain.

MBF1 (multi-protein bridging factor 1) is a protein containing a conserved HTH (helix-turn-helix) domain in both eukaryotes and archaea. Eukaryotic MBF1 has been reported to function as a transcriptional co-activator that physically bridges transcription regulators with the core transcription initiation machinery of RNA polymerase II. In addition, MBF1 has been found to be associated with polyadenylated mRNA in yeast as well as in mammalian cells. aMBF1 (archaeal MBF1) is very well conserved among most archaeal lineages; however, its function has so far remained elusive. To address this, we have conducted a molecular characterization of this aMBF1. Affinity purification of interacting proteins indicates that aMBF1 binds to ribosomal subunits. On sucrose density gradients, aMBF1 co-fractionates with free 30S ribosomal subunits as well as with 70S ribosomes engaged in translation. Binding of aMBF1 to ribosomes does not inhibit translation. Using NMR spectroscopy, we show that aMBF1 contains a long intrinsically disordered linker connecting the predicted N-terminal zinc-ribbon domain with the C-terminal HTH domain. The HTH domain, which is conserved in all archaeal and eukaryotic MBF1 homologues, is directly involved in the association of aMBF1 with ribosomes. The disordered linker of the ribosome-bound aMBF1 provides the N-terminal domain with high flexibility in the aMBF1-ribosome complex. Overall, our findings suggest a role for aMBF1 in the archaeal translation process.

Exome sequencing reveals a novel Moroccan founder mutation in SLC19A3 as a new cause of early-childhood fatal Leigh syndrome.

Leigh syndrome is an early onset, often fatal progressive neurodegenerative disorder caused by mutations in the mitochondrial or nuclear DNA. Until now, mutations in more than 35 genes have been reported to cause Leigh syndrome, indicating an extreme genetic heterogeneity for this disorder, but still only explaining part of the cases. The possibility of whole exome sequencing enables not only mutation detection in known candidate genes, but also the identification of new genes associated with Leigh syndrome in small families and isolated cases. Exome sequencing was combined with homozygosity mapping to identify the genetic defect in a Moroccan family with fatal Leigh syndrome in early childhood and specific magnetic resonance imaging abnormalities in the brain. We detected a homozygous nonsense mutation (c.20C>A; p.Ser7Ter) in the thiamine transporter SLC19A3. In vivo overexpression of wild-type SLC19A3 showed an increased thiamine uptake, whereas overexpression of mutant SLC19A3 did not, confirming that the mutation results in an absent or non-functional protein. Seventeen additional patients with Leigh syndrome were screened for mutations in SLC19A3 using conventional Sanger sequencing. Two unrelated patients, both from Moroccan origin and one from consanguineous parents, were homozygous for the same p.Ser7Ter mutation. One of these patients showed the same MRI abnormalities as the patients from the first family. Strikingly, patients receiving thiamine had an improved life-expectancy. One patient in the third family deteriorated upon interruption of the thiamine treatment and recovered after reinitiating. Although unrelated, all patients came from the province Al Hoceima in Northern Morocco. Based on the recombination events the mutation was estimated to have occurred 1250-1750 years ago. Our data shows that SLC19A3 is a new candidate for mutation screening in patients with Leigh syndrome, who might benefit from high doses of thiamine and/or biotin. Especially, Moroccan patients with Leigh syndrome should be tested for the c.20C>A founder mutation in SLC19A3.

Fidelity in archaeal information processing.

A key element during the flow of genetic information in living systems is fidelity. The accuracy of DNA replication influences the genome size as well as the rate of genome evolution. The large amount of energy invested in gene expression implies that fidelity plays a major role in fitness. On the other hand, an increase in fidelity generally coincides with a decrease in velocity. Hence, an important determinant of the evolution of life has been the establishment of a delicate balance between fidelity and variability. This paper reviews the current knowledge on quality control in archaeal information processing. While the majority of these processes are homologous in Archaea, Bacteria, and Eukaryotes, examples are provided of nonorthologous factors and processes operating in the archaeal domain. In some instances, evidence for the existence of certain fidelity mechanisms has been provided, but the factors involved still remain to be identified.

Role of multiprotein bridging factor 1 in archaea: bridging the domains?

MBF1 (multiprotein bridging factor 1) is a highly conserved protein in archaea and eukaryotes. It was originally identified as a mediator of the eukaryotic transcription regulator BmFTZ-F1 (Bombyx mori regulator of fushi tarazu). MBF1 was demonstrated to enhance transcription by forming a bridge between distinct regulatory DNA-binding proteins and the TATA-box-binding protein. MBF1 consists of two parts: a C-terminal part that contains a highly conserved helix-turn-helix, and an N-terminal part that shows a clear divergence: in eukaryotes, it is a weakly conserved flexible domain, whereas, in archaea, it is a conserved zinc-ribbon domain. Although its function in archaea remains elusive, its function as a transcriptional co-activator has been deduced from thorough studies of several eukaryotic proteins, often indicating a role in stress response. In addition, MBF1 was found to influence translation fidelity in yeast. Genome context analysis of mbf1 in archaea revealed conserved clustering in the crenarchaeal branch together with genes generally involved in gene expression. It points to a role of MBF1 in transcription and/or translation. Experimental data are required to allow comparison of the archaeal MBF1 with its eukaryotic counterpart.

Relative quantification of proteins across the species boundary through the use of shared peptides.

We show that shared peptides of proteins that are encoded in different species are suitable for cross-species relative protein quantification. A 14N-containing proteome from the thermoacidophilic archaeon Sulfolobus tokodaii was mixed with a 15N-labeled proteome from Sulfolobus solfataricus. Using three shared peptides per protein, the relative abundance of six orthologous proteins was calculated. Observed standard deviations were approximately 10%, indicating that the trypsin accessibility to cleavage sites was not altered in the orthologs. The abundance ratios of the and subunits of the Thermosome were 0.64 and 1.24 in Sulfolobus tokodaii compared to Sulfolobus solfataricus, suggesting a different stoichiometry of the complex in both species. In addition, an in silico study was performed on the occurrence of shared peptides. Inter- and intra-species peptide redundancy was investigated in the model organisms Homo sapiens, Mus musculus, Escherichia coli K12, Escherichia coli O157:H7, S. solfataricus, and S. tokodaii. M. musculus and H. sapiens share 30-50% of all peptides (6-15 residues). Moreover, approximately one-third of all proteins shared > or = 40% of their peptides with at least one other protein in the related species, thus offering strong potential for cross-species relative protein quantification. Conversely, approximately 40% of all peptides (6-15 residues) encoded in H. sapiens are encoded multiple times and therefore complicate identification and quantification.