Salmeterol enhances the cardiac response to gene therapy in Pompe disease.

Enzyme replacement therapy (ERT) with recombinant human (rh) acid α-glucosidase (GAA) has prolonged the survival of patients. However, the paucity of cation-independent mannose-6-phosphate receptor (CI-MPR) in skeletal muscle, where it is needed to take up rhGAA, correlated with a poor response to ERT by muscle in Pompe disease. Clenbuterol, a selective β2 receptor agonist, enhanced the CI-MPR expression in striated muscle through Igf-1 mediated muscle hypertrophy, which correlated with increased CI-MPR (also the Igf-2 receptor) expression. In this study we have evaluated 4 new drugs in GAA knockout (KO) mice in combination with an adeno-associated virus (AAV) vector encoding human GAA, 3 alternative β2 agonists and dehydroepiandrosterone (DHEA). Mice were injected with AAV2/9-CBhGAA (1E+11 vector particles) at a dose that was not effective at clearing glycogen storage from the heart. Heart GAA activity was significantly increased by either salmeterol (p<0.01) or DHEA (p<0.05), in comparison with untreated mice. Furthermore, glycogen content was reduced in the heart by treatment with DHEA (p<0.001), salmeterol (p<0.05), formoterol (p<0.01), or clenbuterol (p<0.01) in combination with the AAV vector, in comparison with untreated GAA-KO mice. Wirehang testing revealed that salmeterol and the AAV vector significantly increased performance, in comparison with the AAV vector alone (p<0.001). Similarly, salmeterol with the vector increased performance significantly more than any of the other drugs. The most effective individual drugs had no significant effect in absence of vector, in comparison with untreated mice. Thus, salmeterol should be further developed as adjunctive therapy in combination with either ERT or gene therapy for Pompe disease.

The Role of TET Proteins in the Epigenetic Regulation of Neural Gene Expression and Behavior

Understanding how genes affect behavior is critical to develop precise therapies for human behavioral disorders. The ability to investigate the relationship between genes and behavior has been greatly advanced over the last few decades due to progress in gene-targeting technology. Recently, the Tet gene family was discovered and implicated in epigenetic modification of DNA methylation by converting 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). 5hmC and its catalysts, the TET proteins, are highly abundant in the postnatal brain but with unclear functions. To investigate their neural functions, we generated new lines of Tet1 and Tet3 mutant mice using a gene targeting approach. We designed both mutations to cause a frameshift by deleting the largest coding exon of Tet1 (Tet1Δe4) and the catalytic domain of Tet3 (Tet3Δe7-9). As Tet1 is also highly expressed in embryonic stem cells (ESCs), we generated Tet1 homozygous deleted ESCs through sequential targeting to compare the function of Tet1 in the brain to its role in ESCs. To test our hypothesis that TET proteins epigenetically regulate transcription of key neural genes important for normal brain function, we examined transcriptional and epigenetic differences in the Tet1Δe4 mouse brain. The oxytocin receptor (OXTR), a neural gene implicated in social behaviors, is suggested to be epigenetically regulated by an unknown mechanism. Interestingly, several human studies have found associations between OXTR DNA hypermethylation and a wide spectrum of behavioral traits and neuropsychiatric disorders including autism spectrum disorders. Here we report the first evidence for an epigenetic mechanism of Oxtr transcription as expression of Oxtr is reduced in the brains of Tet1Δe4-/- mice. Likewise, the CpG island overlapping the promoter of Oxtr is hypermethylated during early embryonic development and persists into adulthood. We also discovered altered histone modifications at the hypermethylated regions, indicating the loss of TET1 has broad effects on the chromatin structure at Oxtr. Unexpectedly, we discovered an array of novel mRNA isoforms of Oxtr that are selectively reduced in Tet1Δe4-/- mice. Additionally, Tet1Δe4-/- mice display increased agonistic behaviors and impaired maternal care and short-term memory. Our findings support a novel role for TET1 in regulating Oxtr expression by preventing DNA hypermethylation and implicate TET1 in social behaviors, offering novel insight into Oxtr epigenetic regulation and its role in neuropsychiatric disorders.

The Influence of Gene Environment Interaction on the Risk of Cognitive Impairment: Reducing Sexual Risk Behaviors and Alcohol Use in HIV-infected Adults

Memory deficits and executive dysfunction are highly prevalent among HIV-infected adults. These conditions can affect their quality of life, antiretroviral adherence, and HIV risk behaviors. Several factors have been suggested including the role of genetics in relation to HIV disease progression. This dissertation aimed to determine whether genetic differences in HIV-infected individuals were correlated with impaired memory, cognitive flexibility and executive function and whether cognitive decline moderated alcohol use and sexual transmission risk behaviors among HIV-infected alcohol abusers participating in an NIH-funded clinical trial comparing the efficacy of the adapted Holistic Health Recovery Program (HHRP-A) intervention to a Health Promotion Control (HPC) condition in reducing risk behaviors. A total of 267 individuals were genotyped for polymorphisms in the dopamine and serotonin gene systems. Results yielded significant associations for TPH2, GALM, DRD2 and DRD4 genetic variants with impaired executive function, cognitive flexibility and memory. SNPs TPH2 rs4570625 and DRD2 rs6277 showed a risk association with executive function (odds ratio = 2.5, p = .02; 3.6, p = .001). GALM rs6741892 was associated with impaired memory (odds ratio = 1.9, p = .006). At the six-month follow-up, HHRP-A participants were less likely to report trading sex for food, drugs and money (20.0%) and unprotected insertive or receptive oral (11.6%) or vaginal and/or anal sex (3.2%) than HPC participants (49.4%, p

Impacts of seawater acidification on mantle gene expression patterns of the Baltic Sea blue mussel: implications for shell formation and energy metabolism, link to supplementary material

Marine organisms have to cope with increasing CO2 partial pressures and decreasing pH in the oceans. We elucidated the impacts of an 8-week acclimation period to four seawater pCO2 treatments (39, 113, 243 and 405 Pa/385, 1,120, 2,400 and 4,000 µatm) on mantle gene expression patterns in the blue mussel Mytilus edulis from the Baltic Sea. Based on the M. edulis mantle tissue transcriptome, the expression of several genes involved in metabolism, calcification and stress responses was assessed in the outer (marginal and pallial zone) and the inner mantle tissues (central zone) using quantitative real-time PCR. The expression of genes involved in energy and protein metabolism (F-ATPase, hexokinase and elongation factor alpha) was strongly affected by acclimation to moderately elevated CO2 partial pressures. Expression of a chitinase, potentially important for the calcification process, was strongly depressed (maximum ninefold), correlating with a linear decrease in shell growth observed in the experimental animals. Interestingly, shell matrix protein candidate genes were less affected by CO2 in both tissues. A compensatory process toward enhanced shell protection is indicated by a massive increase in the expression of tyrosinase, a gene involved in periostracum formation (maximum 220-fold). Using correlation matrices and a force-directed layout network graph, we were able to uncover possible underlying regulatory networks and the connections between different pathways, thereby providing a molecular basis of observed changes in animal physiology in response to ocean acidification.

Enhancement of a novel gene therapy approach for Sandhoff disease through complimentary drug therapy

GM2 gangliosidoses is a family of severe, neurodegenerative disorders resulting from a deficiency in the β-hexosaminidase A (Hex A) enzyme. This disorder is typically caused by a mutation to either the HEXA gene, causing Tay Sachs disease, or a mutation to the HEXB gene, causing Sandhoff disease. The HEXA and HEXB genes are required to produce the α and β subunits of the Hex A enzyme respectively. Using a Sandhoff disease (SD) mouse model (Hexb-/-) we tested the potential of a low dose of systemically delivered single stranded adeno-associated virus 9 (ssAAV9) expressing human HEXB and human HEXA cDNA under the control of a single promoter through the use of a bicistronic vector design with a P2A linker to correct the neurological phenotype. Neonatal mice were injected with either this ssAAV9-HexB-P2A-HexA vector (HexB-HexA) or a vehicle solution via the superficial temporal vein. HexB-HexA treatment alone conferred an increase in survival of 56% compared to vehicle-injected controls and biochemical analysis of the brain tissue and serum revealed an increase in HexA activity and a decrease in brain GM2 ganglioside buildup. Additionally, treatments with the non-steroidal anti-inflammatory drug indomethacin (Indo), the histone deactylase inhibitor ITF2357 (ITF) and the pharmacological chaperone pyrimethamine (Pyr) were tested. The anti-inflammatory treatments of Indo and ITF conferred an increase in survival of 12% and 8% respectively while causing no alteration in the HexA activity or GM2 ganglioside buildup. Pyr had no observable effect on disease progression. Lastly HexB-HexA treatment was tested in conjunction with Indo, ITF and Pyr individually. Additive increases in survival and behavioural testing results were observed with Indo and ITF treatments while no additional benefit to HexA activity or GM2 ganglioside levels in the brain tissue was observed. This indicates the two treatments slowed the progression of the disease through a different mechanism than the reduction of the GM2 ganglioside substrate. Pyr treatment was shown to have no effect when combined with HexB-HexA treatment. This study demonstrates the potential amelioration of SD with a novel AAV9 gene therapy approach as well as helped to identify the additive potential of anti-inflammatory treatments in gene therapy of GM2 gangliosidoses.

Over 30% of patients with splenic marginal zone lymphoma express the same immunoglobulin heavy variable gene: ontogenetic implications.

We performed an immunogenetic analysis of 345 IGHV-IGHD-IGHJ rearrangements from 337 cases with primary splenic small B-cell lymphomas of marginal-zone origin. Three immunoglobulin (IG) heavy variable (IGHV) genes accounted for 45.8% of the cases (IGHV1-2, 24.9%; IGHV4-34, 12.8%; IGHV3-23, 8.1%). Particularly for the IGHV1-2 gene, strong biases were evident regarding utilization of different alleles, with 79/86 rearrangements (92%) using allele (*)04. Among cases more stringently classified as splenic marginal-zone lymphoma (SMZL) thanks to the availability of splenic histopathological specimens, the frequency of IGHV1-2(*)04 peaked at 31%. The IGHV1-2(*)04 rearrangements carried significantly longer complementarity-determining region-3 (CDR3) than all other cases and showed biased IGHD gene usage, leading to CDR3s with common motifs. The great majority of analyzed rearrangements (299/345, 86.7%) carried IGHV genes with some impact of somatic hypermutation, from minimal to pronounced. Noticeably, 75/79 (95%) IGHV1-2(*)04 rearrangements were mutated; however, they mostly (56/75 cases; 74.6%) carried few mutations (97-99.9% germline identity) of conservative nature and restricted distribution. These distinctive features of the IG receptors indicate selection by (super)antigenic element(s) in the pathogenesis of SMZL. Furthermore, they raise the possibility that certain SMZL subtypes could derive from progenitor populations adapted to particular antigenic challenges through selection of VH domain specificities, in particular the IGHV1-2(*)04 allele.

Specific JAK2 mutation (JAK2R683) and multiple gene deletions in Down syndrome acute lymphoblastic leukemia

Children with Down syndrome (DS) have a greatly increased risk of acute megakaryoblastic leukemia (AMKL) and acute lymphoblastic leukemia (ALL). Both DS-AMKL and the related transient myeloproliferative disorder (TMD) have GATA1 mutations as obligatory, early events. To identify mutations contributing to leukemogenesis in DS-ALL, we undertook sequencing of candidate genes, including FLT3, RAS, PTPN11, BRAF, and JAK2. Sequencing of the JAK2 pseudokinase domain identified a specific, acquired mutation, JAK2R683, in 12 (28%) of 42 DS-ALL cases. Functional studies of the common JAK2R683G mutation in murine Ba/F3 cells showed growth factor independence and constitutive activation of the JAK/STAT signaling pathway. High-resolution SNP array analysis of 9 DS-ALL cases identified additional submicroscopic deletions in key genes, including ETV6, CDKN2A, and PAX5. These results infer a complex molecular pathogenesis for DS-ALL leukemogenesis, with trisomy 21 as an initiating or first hit and with chromosome aneuploidy, gene deletions, and activating JAK2 mutations as complementary genetic events. (Blood. 2009; 113: 646-648)

Gene mapping and expression analysis of 16q loss of heterozygosity identifies WWOX and CYLD as being important in determining clinical outcome in multiple myeloma.

We performed fluorescent in situ hybridization (FISH) for 16q23 abnormalities in 861 patients with newly diagnosed multiple myeloma and identified deletion of 16q [del(16q)] in 19.5%. In 467 cases in which demographic and survival data were available, del(16q) was associated with a worse overall survival (OS). It was an independent prognostic marker and conferred additional adverse survival impact in cases with the known poor-risk cytogenetic factors t(4;14) and del(17p). Gene expression profiling and gene mapping using 500K single-nucleotide polymorphism (SNP) mapping arrays revealed loss of heterozygosity (LOH) involving 3 regions: the whole of 16q, a region centered on 16q12 (the location of CYLD), and a region centered on 16q23 (the location of the WW domain-containing oxidoreductase gene WWOX). CYLD is a negative regulator of the NF-kappaB pathway, and cases with low expression of CYLD were used to define a "low-CYLD signature." Cases with 16q LOH or t(14;16) had significantly reduced WWOX expression. WWOX, the site of the translocation breakpoint in t(14;16) cases, is a known tumor suppressor gene involved in apoptosis, and we were able to generate a "low-WWOX signature" defined by WWOX expression. These 2 genes and their corresponding pathways provide an important insight into the potential mechanisms by which 16q LOH confers poor prognosis.

Immunoglobulin gene rearrangements and the pathogenesis of multiple myeloma.

The ability to rearrange the germ-line DNA to generate antibody diversity is an essential prerequisite for the production of a functional repertoire. While this is essential to prevent infections, it also represents the "Achilles heel" of the B-cell lineage, occasionally leading to malignant transformation of these cells by translocation of protooncogenes into the immunoglobulin (Ig) loci. However, in evolutionary terms this is a small price to pay for a functional immune system. The study of the configuration and rearrangements of the Ig gene loci has contributed extensively to our understanding of the natural history of development of myeloma. In addition to this, the analysis of Ig gene rearrangements in B-cell neoplasms provides information about the clonal origin of the disease, prognosis, as well as providing a clinical useful tool for clonality detection and minimal residual disease monitoring. Herein, we review the data currently available on both Ig gene rearrangements and protein patterns seen in myeloma with the aim of illustrating how this knowledge has contributed to our understanding of the pathobiology of myeloma.

Molecular characterization of complete and incomplete immunoglobulin heavy chain gene rearrangements in hairy cell leukemia.

PURPOSE: We analyzed patients with hairy cell leukemia (HCL) to achieve a better understanding of the differentiation stage reached by HCL cells and to define the key role of the diversification of cell surface makers, especially CD25 expression. PATIENTS AND METHODS: We analyzed 38 previously untreated patients with HCL to characterize their complete (VDJ(H)) and incomplete (DJ(H)) immunoglobulin (Ig) heavy chain (IgH) rearrangements, including somatic hypermutation pattern and gene segment use. RESULTS: A correlation between immunophenotypic profile and molecular data was seen. All 38 cases showed monoclonal amplifications: VDJ(H) in 97%, DJ(H) in 42%, and both in 39%. Segments from the D(H)3 family were used more in complete compared with incomplete rearrangements (45% vs. 12%; P <.005). Furthermore, comparison between molecular and immunophenotypic characteristics disclosed differences in the expression of CD25 antigen; CD25(-) cases, a phenotype associated with HCL variant, showed complete homology to the germline in 3 of 5 cases (60%), whereas this characteristic was never observed in CD25(+) cases (P <.005). Moreover, V(H)4-34, V(H)1-08, and J(H)3 segments appeared in 2, 1, and 2 CD25(-) cases, respectively, whereas they were absent in all CD25(+) cases. CONCLUSION: These results support that HCL is a heterogeneous entity including subgroups with different molecular characteristics, which reinforces the need for additional studies with a larger number of patients to clarify the real role of gene rearrangements in HCL.

The Calreticulin gene and myeloproliferative neoplasms

The Philadelphia negative myeloproliferative neoplasms include polycythaemia vera (PV), essential thrombocytopenia (ET) and primary myelofibrosis (PMF). Patients with these conditions were mainly thought to harbour JAK2V617F mutations or an Myeloproliferative leukaemia (MPL) substitution. In 2013, two revolutionary studies identified recurrent mutations in a gene that encodes the protein calreticulin (CALR). This mutation was detected in patients with PMF and ET with non-mutated JAK2 or MPL but was absent in patients with PV. The CALR gene encodes the calreticulin protein, which is a multifactorial protein, mainly located in the endoplasmic reticulum in chromosome 19 and regulates calcium homeostasis, chaperones and has also been implicated in multiple cellular processes including cell signalling, regulation of gene expression, cell adhesion, autoimmunity and apoptosis. Somatic 52 bp deletions and recurrent 52 bp insertion mutations in CALR were detected and all resulted in frameshift and clusters in exon 9 of the gene. This review will summarise the current knowledge on the CALR gene and mutation of the gene in pathological conditions and patient phenotypes.

Delivery of nucleic acids for cancer gene therapy: overcoming extra- and intra-cellular barriers

The therapeutic potential of cancer gene therapy has been limited by the difficulty of delivering genetic material to target sites. Various biological and molecular barriers exist which need to be overcome before effective nonviral delivery systems can be applied successfully in oncology. Herein, various barriers are described and strategies to circumvent such obstacles are discussed, considering both the extracellular and intracellular setting. Development of multifunctional delivery systems holds much promise for the progression of gene delivery, and a growing body of evidence supports this approach involving rational design of vectors, with a unique molecular architecture. In addition, the potential application of composite gene delivery platforms is highlighted which may provide an alternative delivery strategy to traditional systemic administration.

Analyzing Gene Centres with the Help of the Checklist Method – the Case of Syria

The present survey of species diversity of cultivated plants is the first for Syria. Some cultivated species will be added in the future, because due to the civil war in Syria, it was not possible to visit the country in the frame of the present work, as initially planned. Checklists proved to be a useful tool for overviewing the cultivated plants of selected areas and allow a characterization of the state of plant genetic resources of Syria. Syria has experienced several civilizations. Man settled in this productive land since ancient times and used its resources. However, such use has led to changes in vegetation and decline of wildlife through the country, in seashore areas, interior, mountains, and grassland. Plant domestication and growing started more than 10,000 years ago in West Asia. Since then, plentiful of economic plant species were present and used by man and his domesticated animals. Forming a part of the Fertile Crescent, where many of the world’s agricultural plants have evolved, Syria is extremely rich in agrobiodiversity. Wild progenitors of wheat and barley and wild relatives of many fruit trees such as almonds and pistachio as well as forage species are still found in marginal lands and less disturbed areas. These are threatened by a wide range of human activities, notably modern, extensive agriculture, overgrazing, overcutting and urban expansion. Syria is also considered as part of one of the main centres of origin, according to Vavilov, who had collected in Syria in 1926. The first expeditions to crop fields showed the exclusive nature of cultivated plants in Syria with a high number of endemic forms. Furthermore, Syria is a part of a biodiversity hotspot. Several studies have been performed to study agrobiodiversity in different parts of Syria, but usually on wild species. Many collections have been carried out; however, they focussed preferably on cereals and pulses, and particularly on wheat, like Vavilov’s expedition. Only 30 crops make up the major part of the conserved Syrian crop plant material in the genebank, indicating that most of the remaining 7,000 species of cultivated plants and many other valuable genetic resources species have only been included on a limited scale in the genebank collections. Although a small country (185,180 km2), Syria accommodates numerous ecosystems that allow for a large diversity of plant genetic resources for agriculture ranging from cold-requiring to subtropical crops to live and thrive. Only few references are available in this respect. The aim of the present study was to complete a checklist of Syria’s cultivated plants of agriculture and horticulture excluding plants only grown as ornamental or for forestry. Furthermore, plants taken for reforestation have not been included, if they do not have also agricultural or horticultural uses. Therefore, the inclusion of plants into the checklist follows the same principles as “Mansfeld’s Encyclopedia”. Main sources of information were published literature, floras of Syria, Lebanon and the Mediterranean, as well as Syrian printed sources in Arabic and/or English, reports from FAO on agricultural statistics in Syria, and data from ICARDA and Bioversity International. In addition, personal observations gathered during professional work in the General Commission for Scientific Agricultural Research (GCSAR) in Syria (since 1989) and participation in projects were taken into account. These were: (1) A project on “Conservation and Sustainable Use of Dry Land Agrobiodiversity in the Near East” with participation of Jordan, Lebanon, Syria, and the Palestinian Authority, focussing on landraces and wild relatives of barley, wheat, lentil, alliums, feed legumes, and fruit trees (1999–2005). (2) A project for vegetable landraces (1993–1995) in collaboration with the former International Plant Genetic Resources Institute and the UN Development Programme, in which 380 local vegetable accessions were evaluated. For medicinal plants and fruit trees I was in personal contact with departments of GCSAR and the Ministry of Agriculture and Agrarian Reform, as well as with private organizations. The resulting checklist was compared with the catalogues of crop plants of Italy and a checklist of cultivated plants of Iraq. The cultivated plant species are presented in alphabetical order according to their accepted scientific names. Each entry consists of a nomenclatural part, folk names, details of plant uses, the distribution in Syria (by provinces), a textual description, and references to literature. In total, 262 species belonging to 146 genera and 57 families were identified. Within-species (intraspecific) diversity is a significant measure of the biodiversity. Intraspecific diversity for wild plants has been and remains to be well studied, but for crop plants there are only few results. Mansfeld’s method is an actual logical contribution to such studies. Among the families, the following have the highest number of crop species: Leguminosae (34 spp.), Rosaceae (24), Gramineae (18), Labiatae (18), Compositae (14), Cruciferae (14), Cucurbitaceae (11), Rutaceae (10), Malvaceae (9), Alliaceae (7), and Anacardiaceae (7). The establishment of an effective programme for the maintenance of plant genetic resources in Syria started in the mid-1970s. This programme considered ex situ and in situ collection of the genetic resources of various field crops, fruit trees and vegetables. From a plant genetic resources viewpoint, it is clear that the homegarden is an important location for the cultivation of so-called neglected and underutilized species (neglected from a research side and underutilized from a larger economic side). Such species have so far not received much care from ecologists, botanists and agronomists, and they are considerably under-represented in genebanks.

Fructose Alters Adiponectin, Haptoglobin and Angiotensinogen Gene Expression in 3T3-L1 Adipocytes

Fructose- or sucrose-rich diets can cause insulin resistance and increase the risk of cardiovascular disease. Adipokines are correlated with the development of these diseases in obesity. We hypothesize that fructose and sucrose induce insulin resistance via effects on adipokine gene expression in adipocytes. This study analyzed the effect of fructose or glucose on adiponectin, haptoglobin, and angiotensinogen gene expression in 3T3-L1 adipocytes. Ten days after differentiation, the cells were pretreated with serum- and glucose-free medium. Twenty-four hours later, fructose or glucose (0, 5, 10, or 20 mmol) was added into the medium, and the cells were collected after a further 24 hours. Adiponectin, haptoglobin, and angiotensinogen gene expression were determined. Adiponectin gene expression increased when 10 or 20 mmol glucose was added compared with that observed for the non–hexose-treated cells. A similar effect occurred when 5 mmol fructose was added. Glucose (10 mmol) and fructose (20 mmol) stimulated haptoglobin gene expression in 3T3-L1 adipocytes compared with 0 mmol, with glucose producing a more pronounced effect. Although 20 mmol fructose caused an increase in angiotensinogen gene expression, glucose did not. In conclusion, in this study of 2 hexoses revealed an increase in adiponectin gene expression, suggesting that the effect of a glucose-rich diet on the development of insulin resistance is not related to the effect of these hexoses on adipocyte adiponectin gene expression. However, insulin resistance and cardiovascular disease promoted by fructose-rich diets could be partially related to the effect of fructose on adiponectin and angiotensinogen gene expression.

Data supporting the co-expression of PDHA1 gene and of its paralogue PDHA2 in somatic cells of a family

This article presents a dataset proving the simultaneous presence of a 5′UTR-truncated PDHA1 mRNA and a full-length PDHA2 mRNA in the somatic cells of a PDC-deficient female patient and all members of her immediate family (parents and brother). We have designed a large set of primer pairs in order to perform detailed RT-PCR assays allowing the clear identification of both PDHA1 and PDHA2 mRNA species in somatic cells. In addition, two different experimental approaches were used to elucidate the copy number of PDHA1 gene in the patient and her mother. The interpretation and discussion of these data, along with further extensive experiments concerning the origin of this altered gene expression and its potential therapeutic consequences, can be found in “Complex genetic findings in a female patient with pyruvate dehydrogenase complex deficiency: null mutations in the PDHX gene associated with unusual expression of the testis-specific PDHA2 gene in her somatic cells” (A. Pinheiro, M.J. Silva, C. Florindo, et al., 2016).