Combinatorial optimization of gene expression through recombinase-mediated promoter and terminator shuffling in yeast.

Microbes are increasingly employed as cell factories to produce biomolecules. This often involves the expression of complex heterologous biosynthesis pathways in host strains. Achieving maximal product yields and avoiding build-up of (toxic) intermediates requires balanced expression of every pathway gene. However, despite progress in metabolic modeling, the optimization of gene expression still heavily relies on trial-and-error. Here, we report an approach for in vivo, multiplexed Gene Expression Modification by LoxPsym-Cre Recombination (GEMbLeR). GEMbLeR exploits orthogonal LoxPsym sites to independently shuffle promoter and terminator modules at distinct genomic loci. This approach facilitates creation of large strain libraries, in which expression of every pathway gene ranges over 120-fold and each strain harbors a unique expression profile. When applied to the biosynthetic pathway of astaxanthin, an industrially relevant antioxidant, a single round of GEMbLeR improved pathway flux and doubled production titers. Together, this shows that GEMbLeR allows rapid and efficient gene expression optimization in heterologous biosynthetic pathways, offering possibilities for enhancing the performance of microbial cell factories.

Mitochondrial GpC and CpG DNA Hypermethylation Cause Metabolic Stress-Induced Mitophagy and Cholestophagy.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is characterized by a constant accumulation of lipids in the liver. This hepatic lipotoxicity is associated with a dysregulation of the first step in lipid catabolism, known as beta oxidation, which occurs in the mitochondrial matrix. Eventually, this dysregulation will lead to mitochondrial dysfunction. To evaluate the possible involvement of mitochondrial DNA methylation in this lipid metabolic dysfunction, we investigated the functional metabolic effects of mitochondrial overexpression of CpG (MSssI) and GpC (MCviPI) DNA methyltransferases in relation to gene expression and (mito)epigenetic signatures. Overall, the results show that mitochondrial GpC and, to a lesser extent, CpG methylation increase bile acid metabolic gene expression, inducing the onset of cholestasis through mito-nuclear epigenetic reprogramming. Moreover, both increase the expression of metabolic nuclear receptors and thereby induce basal overactivation of mitochondrial respiration. The latter promotes mitochondrial swelling, favoring lipid accumulation and metabolic-stress-induced mitophagy and autophagy stress responses. In conclusion, both mitochondrial GpC and CpG methylation create a metabolically challenging environment that induces mitochondrial dysfunction, which may contribute to the progression of MASLD.

Mitochondrial DNA methylation in metabolic associated fatty liver disease.

Introduction: Hepatic lipid accumulation and mitochondrial dysfunction are hallmarks of metabolic associated fatty liver disease (MAFLD), yet molecular parameters underlying MAFLD progression are not well understood. Differential methylation within the mitochondrial DNA (mtDNA) has been suggested to be associated with dysfunctional mitochondria, also during progression to Metabolic Steatohepatitis (MeSH). This study further investigates whether mtDNA methylation is associated with hepatic lipid accumulation and MAFLD. Methods: HepG2 cells were constructed to stably express mitochondria-targeted viral and prokaryotic cytosine DNA methyltransferases (mtM.CviPI or mtM.SssI for GpC or CpG methylation, respectively). A catalytically inactive variant (mtM.CviPI-Mut) was constructed as a control. Mouse and human patients' samples were also investigated. mtDNA methylation was assessed by pyro- or nanopore sequencing. Results and discussion: Differentially induced mtDNA hypermethylation impaired mitochondrial gene expression and metabolic activity in HepG2-mtM.CviPI and HepG2-mtM.SssI cells and was associated with increased lipid accumulation, when compared to the controls. To test whether lipid accumulation causes mtDNA methylation, HepG2 cells were subjected to 1 or 2 weeks of fatty acid treatment, but no clear differences in mtDNA methylation were detected. In contrast, hepatic Nd6 mitochondrial gene body cytosine methylation and Nd6 gene expression were increased in mice fed a high-fat high cholesterol diet (HFC for 6 or 20 weeks), when compared to controls, while mtDNA content was unchanged. For patients with simple steatosis, a higher ND6 methylation was confirmed using Methylation Specific PCR, but no additional distinctive cytosines could be identified using pyrosequencing. This study warrants further investigation into a role for mtDNA methylation in promoting mitochondrial dysfunction and impaired lipid metabolism in MAFLD.

Whole-exome rare-variant analysis of Alzheimer's disease and related biomarker traits.

INTRODUCTION: Despite increasing evidence of a role of rare genetic variation in the risk of Alzheimer's disease (AD), limited attention has been paid to its contribution to AD-related biomarker traits indicative of AD-relevant pathophysiological processes. METHODS: We performed whole-exome gene-based rare-variant association studies (RVASs) of 17 AD-related traits on whole-exome sequencing (WES) data generated in the European Medical Information Framework for Alzheimer's Disease Multimodal Biomarker Discovery (EMIF-AD MBD) study (n = 450) and whole-genome sequencing (WGS) data from ADNI (n = 808). RESULTS: Mutation screening revealed a novel probably pathogenic mutation (PSEN1 p.Leu232Phe). Gene-based RVAS revealed the exome-wide significant contribution of rare coding variation in RBKS and OR7A10 to cognitive performance and protection against left hippocampal atrophy, respectively. DISCUSSION: The identification of these novel gene-trait associations offers new perspectives into the role of rare coding variation in the distinct pathophysiological processes culminating in AD, which may lead to identification of novel therapeutic and diagnostic targets.

Rare variants in IFFO1, DTNB, NLRC3 and SLC22A10 associate with Alzheimer's disease CSF profile of neuronal injury and inflammation.

Alzheimer's disease (AD) biomarkers represent several neurodegenerative processes, such as synaptic dysfunction, neuronal inflammation and injury, as well as amyloid pathology. We performed an exome-wide rare variant analysis of six AD biomarkers (ß-amyloid, total/phosphorylated tau, NfL, YKL-40, and Neurogranin) to discover genes associated with these markers. Genetic and biomarker information was available for 480 participants from two studies: EMIF-AD and ADNI. We applied a principal component (PC) analysis to derive biomarkers combinations, which represent statistically independent biological processes. We then tested whether rare variants in 9576 protein-coding genes associate with these PCs using a Meta-SKAT test. We also tested whether the PCs are intermediary to gene effects on AD symptoms with a SMUT test. One PC loaded on NfL and YKL-40, indicators of neuronal injury and inflammation. Four genes were associated with this PC: IFFO1, DTNB, NLRC3, and SLC22A10. Mediation tests suggest, that these genes also affect dementia symptoms via inflammation/injury. We also observed an association between a PC loading on Neurogranin, a marker for synaptic functioning, with GABBR2 and CASZ1, but no mediation effects. The results suggest that rare variants in IFFO1, DTNB, NLRC3, and SLC22A10 heighten susceptibility to neuronal injury and inflammation, potentially by altering cytoskeleton structure and immune activity disinhibition, resulting in an elevated dementia risk. GABBR2 and CASZ1 were associated with synaptic functioning, but mediation analyses suggest that the effect of these two genes on synaptic functioning is not consequential for AD development.

Dissection of TAF1 neuronal splicing and implications for neurodegeneration in X-linked dystonia-parkinsonism.

X-linked dystonia-parkinsonism (XDP) is a monogenic neurodegenerative disorder of the basal ganglia, which presents as a combination of hyperkinetic movements and parkinsonian features. The underlying genetic mechanism involves the insertion of a SINE-VNTR-Alu retrotransposon within the TAF1 gene. Interestingly, alterations of TAF1 have been involved in multiple neurological diseases. In XDP, the SINE-VNTR-Alu insertion in TAF1 has been proposed to result in alternative splicing defects, including the decreased incorporation of a neuron-specific microexon annotated as 34'. This mechanism has become controversial as recent studies failed to provide support. In order to resolve this conundrum, we examined the alternative splicing patterns of TAF1 mRNAs in XDP and control brains. The impact of the disease-associated SINE-VNTR-Alu on alternative splicing of microexon 34' was further investigated in cellular assays. Subsequently, microexon 34' incorporation was explored by RT-PCR and Nanopore long-read sequencing of TAF1 mRNAs from XDP and control brains tissues. Using cell-based splicing assays, we demonstrate that presence of the disease-associated SINE-VNTR-Alu does not affect the inclusion of microexon 34'. In addition, we show that (1) microexon 34'-containing TAF1 mRNAs are detected at similar levels in XDP as in controls and that (2) the architecture of TAF1 transcripts is remarkably similar between XDP and controls brains. These results indicate that microexon 34' incorporation into TAF1 mRNA is not affected in XDP brains. Our findings shift the current paradigm of XDP by discounting alternative splicing of TAF1 microexon 34' as the molecular basis for this disease.

Long-Read Sequencing to Unravel Complex Structural Variants of CEP78 Leading to Cone-Rod Dystrophy and Hearing Loss.

Inactivating variants as well as a missense variant in the centrosomal CEP78 gene have been identified in autosomal recessive cone-rod dystrophy with hearing loss (CRDHL), a rare syndromic inherited retinal disease distinct from Usher syndrome. Apart from this, a complex structural variant (SV) implicating CEP78 has been reported in CRDHL. Here we aimed to expand the genetic architecture of typical CRDHL by the identification of complex SVs of the CEP78 region and characterization of their underlying mechanisms. Approaches used for the identification of the SVs are shallow whole-genome sequencing (sWGS) combined with quantitative polymerase chain reaction (PCR) and long-range PCR, or ExomeDepth analysis on whole-exome sequencing (WES) data. Targeted or whole-genome nanopore long-read sequencing (LRS) was used to delineate breakpoint junctions at the nucleotide level. For all SVs cases, the effect of the SVs on CEP78 expression was assessed using quantitative PCR on patient-derived RNA. Apart from two novel canonical CEP78 splice variants and a frameshifting single-nucleotide variant (SNV), two SVs affecting CEP78 were identified in three unrelated individuals with CRDHL: a heterozygous total gene deletion of 235 kb and a partial gene deletion of 15 kb in a heterozygous and homozygous state, respectively. Assessment of the molecular consequences of the SVs on patient's materials displayed a loss-of-function effect. Delineation and characterization of the 15-kb deletion using targeted LRS revealed the previously described complex CEP78 SV, suggestive of a recurrent genomic rearrangement. A founder haplotype was demonstrated for the latter SV in cases of Belgian and British origin, respectively. The novel 235-kb deletion was delineated using whole-genome LRS. Breakpoint analysis showed microhomology and pointed to a replication-based underlying mechanism. Moreover, data mining of bulk and single-cell human and mouse transcriptional datasets, together with CEP78 immunostaining on human retina, linked the CEP78 expression domain with its phenotypic manifestations. Overall, this study supports that the CEP78 locus is prone to distinct SVs and that SV analysis should be considered in a genetic workup of CRDHL. Finally, it demonstrated the power of sWGS and both targeted and whole-genome LRS in identifying and characterizing complex SVs in patients with ocular diseases.

Critical length in long-read resequencing.

Long-read sequencing has substantial advantages for structural variant discovery and phasing of variants compared to short-read technologies, but the required and optimal read length has not been assessed. In this work, we used long reads simulated from human genomes and evaluated structural variant discovery and variant phasing using current best practice bioinformatics methods. We determined that optimal discovery of structural variants from human genomes can be obtained with reads of minimally 20 kb. Haplotyping variants across genes only reaches its optimum from reads of 100 kb. These findings are important for the design of future long-read sequencing projects.

NanoSatellite: accurate characterization of expanded tandem repeat length and sequence through whole genome long-read sequencing on PromethION.

Technological limitations have hindered the large-scale genetic investigation of tandem repeats in disease. We show that long-read sequencing with a single Oxford Nanopore Technologies PromethION flow cell per individual achieves 30× human genome coverage and enables accurate assessment of tandem repeats including the 10,000-bp Alzheimer's disease-associated ABCA7 VNTR. The Guppy "flip-flop" base caller and tandem-genotypes tandem repeat caller are efficient for large-scale tandem repeat assessment, but base calling and alignment challenges persist. We present NanoSatellite, which analyzes tandem repeats directly on electric current data and improves calling of GC-rich tandem repeats, expanded alleles, and motif interruptions.

Structural variants identified by Oxford Nanopore PromethION sequencing of the human genome.

We sequenced the genome of the Yoruban reference individual NA19240 on the long-read sequencing platform Oxford Nanopore PromethION for evaluation and benchmarking of recently published aligners and germline structural variant calling tools, as well as a comparison with the performance of structural variant calling from short-read sequencing data. The structural variant caller Sniffles after NGMLR or minimap2 alignment provides the most accurate results, but additional confidence or sensitivity can be obtained by a combination of multiple variant callers. Sensitive and fast results can be obtained by minimap2 for alignment and a combination of Sniffles and SVIM for variant identification. We describe a scalable workflow for identification, annotation, and characterization of tens of thousands of structural variants from long-read genome sequencing of an individual or population. By discussing the results of this well-characterized reference individual, we provide an approximation of what can be expected in future long-read sequencing studies aiming for structural variant identification.

De novo single-nucleotide and copy number variation in discordant monozygotic twins reveals disease-related genes.

Recent studies have demonstrated genetic differences between monozygotic (MZ) twins. To test the hypothesis that early post-twinning mutational events associate with phenotypic discordance, we investigated a cohort of 13 twin pairs (n = 26) discordant for various clinical phenotypes using whole-exome sequencing and screened for copy number variation (CNV). We identified a de novo variant in PLCB1, a gene involved in the hydrolysis of lipid phosphorus in milk from dairy cows, associated with lactase non-persistence, and a variant in the mitochondrial complex I gene MT-ND5 associated with amyotrophic lateral sclerosis (ALS). We also found somatic variants in multiple genes (TMEM225B, KBTBD3, TUBGCP4, TFIP11) in another MZ twin pair discordant for ALS. Based on the assumption that discordance between twins could be explained by a common variant with variable penetrance or expressivity, we screened the twin samples for known pathogenic variants that are shared and identified a rare deletion overlapping ARHGAP11B, in the twin pair manifesting with either schizotypal personality disorder or schizophrenia. Parent-offspring trio analysis was implemented for two twin pairs to assess potential association of variants of parental origin with susceptibility to disease. We identified a de novo variant in RASD2 shared by 8-year-old male twins with a suspected diagnosis of autism spectrum disorder (ASD) manifesting as different traits. A de novo CNV duplication was also identified in these twins overlapping CD38, a gene previously implicated in ASD. In twins discordant for Tourette's syndrome, a paternally inherited stop loss variant was detected in AADAC, a known candidate gene for the disorder.

Loss of DPP6 in neurodegenerative dementia: a genetic player in the dysfunction of neuronal excitability.

Emerging evidence suggested a converging mechanism in neurodegenerative brain diseases (NBD) involving early neuronal network dysfunctions and alterations in the homeostasis of neuronal firing as culprits of neurodegeneration. In this study, we used paired-end short-read and direct long-read whole genome sequencing to investigate an unresolved autosomal dominant dementia family significantly linked to 7q36. We identified and validated a chromosomal inversion of ca. 4 Mb, segregating on the disease haplotype and disrupting the coding sequence of dipeptidyl-peptidase 6 gene (DPP6). DPP6 resequencing identified significantly more rare variants-nonsense, frameshift, and missense-in early-onset Alzheimer's disease (EOAD, p value = 0.03, OR = 2.21 95% CI 1.05-4.82) and frontotemporal dementia (FTD, p = 0.006, OR = 2.59, 95% CI 1.28-5.49) patient cohorts. DPP6 is a type II transmembrane protein with a highly structured extracellular domain and is mainly expressed in brain, where it binds to the potassium channel Kv4.2 enhancing its expression, regulating its gating properties and controlling the dendritic excitability of hippocampal neurons. Using in vitro modeling, we showed that the missense variants found in patients destabilize DPP6 and reduce its membrane expression (p < 0.001 and p < 0.0001) leading to a loss of protein. Reduced DPP6 and/or Kv4.2 expression was also detected in brain tissue of missense variant carriers. Loss of DPP6 is known to cause neuronal hyperexcitability and behavioral alterations in Dpp6-KO mice. Taken together, the results of our genomic, genetic, expression and modeling analyses, provided direct evidence supporting the involvement of DPP6 loss in dementia. We propose that loss of function variants have a higher penetrance and disease impact, whereas the missense variants have a variable risk contribution to disease that can vary from high to low penetrance. Our findings of DPP6, as novel gene in dementia, strengthen the involvement of neuronal hyperexcitability and alteration in the homeostasis of neuronal firing as a disease mechanism to further investigate.

miRVaS: a tool to predict the impact of genetic variants on miRNAs.

Genetic variants in or near miRNA genes can have profound effects on miRNA expression and targeting. As user-friendly software for the impact prediction of miRNA variants on a large scale is still lacking, we created a tool called miRVaS. miRVaS automates this prediction by annotating the location of the variant relative to functional regions within the miRNA hairpin (seed, mature, loop, hairpin arm, flanks) and by annotating all predicted structural changes within the miRNA due to the variant. In addition, the tool defines the most important region that is predicted to have structural changes and calculates a conservation score that is indicative of the reliability of the structure prediction. The output is presented in a tab-separated file, which enables fast screening, and in an html file, which allows visual comparison between wild-type and variant structures. All separate images are provided for downstream use. Finally, we tested two different approaches on a small test set of published functionally validated genetic variants for their capacity to predict the impact of variants on miRNA expression.

Novel mutations in genes causing hereditary spastic paraplegia and Charcot-Marie-Tooth neuropathy identified by an optimized protocol for homozygosity mapping based on whole-exome sequencing.

PURPOSE: Homozygosity mapping is an effective approach for detecting molecular defects in consanguineous families by delineating stretches of genomic DNA that are identical by descent. Constant developments in next-generation sequencing created possibilities to combine whole-exome sequencing (WES) and homozygosity mapping in a single step. METHODS: Basic optimization of homozygosity mapping parameters was performed in a group of families with autosomal-recessive (AR) mutations for which both single-nucleotide polymorphism (SNP) array and WES data were available. We varied the criteria for SNP extraction and PLINK thresholds to estimate their effect on the accuracy of homozygosity mapping based on WES. RESULTS: Our protocol showed high specificity and sensitivity for homozygosity detection and facilitated the identification of novel mutations in GAN, GBA2, and ZFYVE26 in four families affected by hereditary spastic paraplegia or Charcot-Marie-Tooth disease. Filtering and mapping with optimized parameters was integrated into the HOMWES (homozygosity mapping based on WES analysis) tool in the GenomeComb package for genomic data analysis. CONCLUSION: We present recommendations for detection of homozygous regions based on WES data and a bioinformatics tool for their identification, which can be widely applied for studying AR disorders.Genet Med 18 6, 600-607.

Schizophrenia-Associated MIR204 Regulates Noncoding RNAs and Affects Neurotransmitter and Ion Channel Gene Sets.

As regulators of gene expression, microRNAs (miRNAs) are likely to play an important role in the development of disease. In this study we present a large-scale strategy to identify miRNAs with a role in the regulation of neuronal processes. Thereby we found variant rs7861254 located near the MIR204 gene to be significantly associated with schizophrenia. This variant resulted in reduced expression of miR-204 in neuronal-like SH-SY5Y cells. Analysis of the consequences of the altered miR-204 expression on the transcriptome of these cells uncovered a new mode of action for miR-204, being the regulation of noncoding RNAs (ncRNAs), including several miRNAs, such as MIR296. Furthermore, pathway analysis showed downstream effects of miR-204 on neurotransmitter and ion channel related gene sets, potentially mediated by miRNAs regulated through miR-204.

Genetic variants in microRNA genes: impact on microRNA expression, function, and disease.

MicroRNAs (miRNAs) are important regulators of gene expression and like any other gene, their coding sequences are subject to genetic variation. Variants in miRNA genes can have profound effects on miRNA functionality at all levels, including miRNA transcription, maturation, and target specificity, and as such they can also contribute to disease. The impact of variants in miRNA genes is the focus of the present review. To put these effects into context, we first discuss the requirements of miRNA transcripts for maturation. In the last part an overview of available databases and tools and experimental approaches to investigate miRNA variants related to human disease is presented.

Copy Number Variations in DISC1 and DISC1-Interacting Partners in Major Mental Illness.

Robust statistical, genetic and functional evidence supports a role for DISC1 in the aetiology of major mental illness. Furthermore, many of its protein-binding partners show evidence for involvement in the pathophysiology of a range of neurodevelopmental and psychiatric disorders. Copy number variants (CNVs) are suspected to play an important causal role in these disorders. In this study, CNV analysis of DISC1 and its binding partners PAFAH1B1, NDE1, NDEL1, FEZ1, MAP1A, CIT and PDE4B in Scottish and Northern Swedish population-based samples was carried out using multiplex amplicon quantification. Here, we report the finding of rare CNVs in DISC1, NDE1 (together with adjacent genes within the 16p13.11 duplication), NDEL1 (including the overlapping MYH10 gene) and CIT. Our findings provide further evidence for involvement of DISC1 and its interaction partners in neuropsychiatric disorders and also for a role of structural variants in the aetiology of these devastating diseases.

CMT-associated mutations in glycyl- and tyrosyl-tRNA synthetases exhibit similar pattern of toxicity and share common genetic modifiers in Drosophila.

Aminoacyl-tRNA synthetases are ubiquitously expressed proteins that charge tRNAs with their cognate amino acids. By ensuring the fidelity of protein synthesis, these enzymes are essential for the viability of every cell. Yet, mutations in six tRNA synthetases specifically affect the peripheral nerves and cause Charcot-Marie-Tooth (CMT) disease. The CMT-causing mutations in tyrosyl- and glycyl-tRNA synthetases (YARS and GARS, respectively) alter the activity of the proteins in a range of ways (some mutations do not impact charging function, while others abrogate it), making a loss of function in tRNA charging unlikely to be the cause of disease pathology. It is currently unknown which cellular mechanisms are triggered by the mutant enzymes and how this leads to neurodegeneration. Here, by expressing two pathogenic mutations (G240R, P234KY) in Drosophila, we generated a model for GARS-associated neuropathy. We observed compromised viability, and behavioral, electrophysiological and morphological impairment in flies expressing the cytoplasmic isoform of mutant GARS. Their features recapitulated several hallmarks of CMT pathophysiology and were similar to the phenotypes identified in our previously described Drosophila model of YARS-associated neuropathy. Furthermore, CG8316 and CG15599 - genes identified in a retinal degeneration screen to modify mutant YARS, also modified the mutant GARS phenotypes. Our study presents genetic evidence for common mutant-specific interactions between two CMT-associated aminoacyl-tRNA synthetases, lending support for a shared mechanism responsible for the synthetase-induced peripheral neuropathies.

A qPCR assay to detect and quantify Shiga toxin-producing E. coli (STEC) in cattle and on farms: a potential predictive tool for STEC culture-positive farms.

Shiga toxin-producing E. coli (STEC), of various serogroups harboring the intimin gene, form a serious threat to human health. They are asymptomatically carried by cattle. In this study, a quantitative real-time PCR (qPCR) method was developed as a molecular method to detect and quantify Shiga toxin genes stx1 and stx2 and the intimin gene eae. Subsequently, 59 fecal samples from six farms were tested using qPCR and a culture method as a reference. Three farms had contaminated animals as demonstrated by the culture method. Culture-positive farms showed moderate significantly higher stx prevalences than culture-negative farms (p = 0.05). This is the first study which showed preliminary results that qPCR can predict STEC farm contamination, with a specificity of 77% and a sensitivity of 83%, as compared with the culture method. Furthermore, the presence or quantity of stx genes in feces was not correlated to the isolation of STEC from the individual animal. Quantitative data thus did not add value to the results. Finally, the detection of both stx and eae genes within the same fecal sample or farm using qPCR was not correlated with the isolation of an eae-harboring STEC strain from the respective sample or farm using the culture method.

Identification of rare copy number variants in high burden schizophrenia families.

Over the last years, genome-wide studies consistently showed an increased burden of rare copy number variants (CNVs) in schizophrenia patients, supporting the "common disease, rare variant" hypothesis in at least a subset of patients. We hypothesize that in families with a high burden of disease, and thus probably a high genetic load influencing disease susceptibility, rare CNVs might be involved in the etiology of schizophrenia. We performed a genome-wide CNV analysis in the index patients of eight families with multiple schizophrenia affected members, and consecutively performed a detailed family analysis for the most relevant CNVs. One index patient showed a DRD5 containing duplication. A second index patient presented with an NRXN1 containing deletion and two adjacent duplications containing MYT1L and SNTG2. Detailed analysis in the subsequent families showed segregation of the identified CNVs. With this study we show the importance of screening high burden families for rare CNVs, which will not only broaden our knowledge concerning the molecular genetic mechanisms involved in schizophrenia but also allow the use of the obtained genetic data to provide better clinical care to these families in general and to non-symptomatic causal CNV carriers in particular.