Genetic aetiologies for childhood speech disorder: novel pathways co-expressed during brain development.

Childhood apraxia of speech (CAS), the prototypic severe childhood speech disorder, is characterized by motor programming and planning deficits. Genetic factors make substantive contributions to CAS aetiology, with a monogenic pathogenic variant identified in a third of cases, implicating around 20 single genes to date. Here we aimed to identify molecular causation in 70 unrelated probands ascertained with CAS. We performed trio genome sequencing. Our bioinformatic analysis examined single nucleotide, indel, copy number, structural and short tandem repeat variants. We prioritised appropriate variants arising de novo or inherited that were expected to be damaging based on in silico predictions. We identified high confidence variants in 18/70 (26%) probands, almost doubling the current number of candidate genes for CAS. Three of the 18 variants affected SETBP1, SETD1A and DDX3X, thus confirming their roles in CAS, while the remaining 15 occurred in genes not previously associated with this disorder. Fifteen variants arose de novo and three were inherited. We provide further novel insights into the biology of child speech disorder, highlighting the roles of chromatin organization and gene regulation in CAS, and confirm that genes involved in CAS are co-expressed during brain development. Our findings confirm a diagnostic yield comparable to, or even higher, than other neurodevelopmental disorders with substantial de novo variant burden. Data also support the increasingly recognised overlaps between genes conferring risk for a range of neurodevelopmental disorders. Understanding the aetiological basis of CAS is critical to end the diagnostic odyssey and ensure affected individuals are poised for precision medicine trials.

Sex-Specific Control of Human Heart Maturation by the Progesterone Receptor.

BACKGROUND: Despite in-depth knowledge of the molecular mechanisms controlling embryonic heart development, little is known about the signals governing postnatal maturation of the human heart. METHODS: Single-nucleus RNA sequencing of 54 140 nuclei from 9 human donors was used to profile transcriptional changes in diverse cardiac cell types during maturation from fetal stages to adulthood. Bulk RNA sequencing and the Assay for Transposase-Accessible Chromatin using sequencing were used to further validate transcriptional changes and to profile alterations in the chromatin accessibility landscape in purified cardiomyocyte nuclei from 21 human donors. Functional validation studies of sex steroids implicated in cardiac maturation were performed in human pluripotent stem cell-derived cardiac organoids and mice. RESULTS: Our data identify the progesterone receptor as a key mediator of sex-dependent transcriptional programs during cardiomyocyte maturation. Functional validation studies in human cardiac organoids and mice demonstrate that the progesterone receptor drives sex-specific metabolic programs and maturation of cardiac contractile properties. CONCLUSIONS: These data provide a blueprint for understanding human heart maturation in both sexes and reveal an important role for the progesterone receptor in human heart development.

Hepatitis C virus leaves an epigenetic signature post cure of infection by direct-acting antivirals.

The increasing worldwide prevalence of Hepatocellular carcinoma (HCC), characterized by resistance to conventional chemotherapy, poor prognosis and eventually mortality, place it as a prime target for new modes of prevention and treatment. Hepatitis C Virus (HCV) is the predominant risk factor for HCC in the US and Europe. Multiple epidemiological studies showed that sustained virological responses (SVR) following treatment with the powerful direct acting antivirals (DAAs), which have replaced interferon-based regimes, do not eliminate tumor development. We aimed to identify an HCV-specific pathogenic mechanism that persists post SVR following DAAs treatment. We demonstrate that HCV infection induces genome-wide epigenetic changes by performing chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) for histone post-translational modifications that are epigenetic markers for active and repressed chromatin. The changes in histone modifications correlate with reprogramed host gene expression and alter signaling pathways known to be associated with HCV life cycle and HCC. These epigenetic alterations require the presence of HCV RNA or/and expression of the viral proteins in the cells. Importantly, the epigenetic changes induced following infection persist as an "epigenetic signature" after virus eradication by DAAs treatment, as detected using in vitro HCV infection models. These observations led to the identification of an 8 gene signature that is associated with HCC development and demonstrate persistent epigenetic alterations in HCV infected and post SVR liver biopsy samples. The epigenetic signature was reverted in vitro by drugs that inhibit epigenetic modifying enzyme and by the EGFR inhibitor, Erlotinib. This epigenetic "scarring" of the genome, persisting following HCV eradication, suggest a novel mechanism for the persistent pathogenesis of HCV after its eradication by DAAs. Our study offers new avenues for prevention of the persistent oncogenic effects of chronic hepatitis infections using specific drugs to revert the epigenetic changes to the genome.

Infanticide vs. inherited cardiac arrhythmias.

AIMS: In 2003, an Australian woman was convicted by a jury of smothering and killing her four children over a 10-year period. Each child died suddenly and unexpectedly during a sleep period, at ages ranging from 19 days to 18 months. In 2019 we were asked to investigate if a genetic cause could explain the children's deaths as part of an inquiry into the mother's convictions. METHODS AND RESULTS: Whole genomes or exomes of the mother and her four children were sequenced. Functional analysis of a novel CALM2 variant was performed by measuring Ca2+-binding affinity, interaction with calcium channels and channel function. We found two children had a novel calmodulin variant (CALM2 G114R) that was inherited maternally. Three genes (CALM1-3) encode identical calmodulin proteins. A variant in the corresponding residue of CALM3 (G114W) was recently reported in a child who died suddenly at age 4 and a sibling who suffered a cardiac arrest at age 5. We show that CALM2 G114R impairs calmodulin's ability to bind calcium and regulate two pivotal calcium channels (CaV1.2 and RyR2) involved in cardiac excitation contraction coupling. The deleterious effects of G114R are similar to those produced by G114W and N98S, which are considered arrhythmogenic and cause sudden cardiac death in children. CONCLUSION: A novel functional calmodulin variant (G114R) predicted to cause idiopathic ventricular fibrillation, catecholaminergic polymorphic ventricular tachycardia, or mild long QT syndrome was present in two children. A fatal arrhythmic event may have been triggered by their intercurrent infections. Thus, calmodulinopathy emerges as a reasonable explanation for a natural cause of their deaths.

Dysregulation of the cohesin subunit RAD21 by Hepatitis C virus mediates host-virus interactions.

Hepatitis C virus (HCV) infection is the leading cause of chronic hepatitis, which often results in liver fibrosis, cirrhosis and hepatocellular carcinoma (HCC). HCV possesses an RNA genome and its replication is confined to the cytoplasm. Yet, infection with HCV leads to global changes in gene expression, and chromosomal instability (CIN) in the host cell. The mechanisms by which the cytoplasmic virus affects these nuclear processes are elusive. Here, we show that HCV modulates the function of the Structural Maintenance of Chromosome (SMC) protein complex, cohesin, which tethers remote regions of chromatin. We demonstrate that infection of hepatoma cells with HCV leads to up regulation of the expression of the RAD21 cohesin subunit and changes cohesin residency on the chromatin. These changes regulate the expression of genes associated with virus-induced pathways. Furthermore, siRNA downregulation of viral-induced RAD21 reduces HCV infection. During mitosis, HCV infection induces hypercondensation of chromosomes and the appearance of multi-centrosomes. We provide evidence that the underlying mechanism involves the viral NS3/4 protease and the cohesin regulator, WAPL. Altogether, our results provide the first evidence that HCV induces changes in gene expression and chromosome structure of infected cells by modulating cohesin.

mitch: multi-contrast pathway enrichment for multi-omics and single-cell profiling data.

BACKGROUND: Inference of biological pathway activity via gene set enrichment analysis is frequently used in the interpretation of clinical and other omics data. With the proliferation of new omics profiling approaches and ever-growing size of data sets generated, there is a lack of tools available to perform and visualise gene set enrichments in analyses involving multiple contrasts. RESULTS: To address this, we developed mitch, an R package for multi-contrast gene set enrichment analysis. It uses a rank-MANOVA statistical approach to identify sets of genes that exhibit joint enrichment across multiple contrasts. Its unique visualisation features enable the exploration of enrichments in up to 20 contrasts. We demonstrate the utility of mitch with case studies spanning multi-contrast RNA expression profiling, integrative multi-omics, tool benchmarking and single-cell RNA sequencing. Using simulated data we show that mitch has similar accuracy to state of the art tools for single-contrast enrichment analysis, and superior accuracy in identifying multi-contrast enrichments. CONCLUSION: mitch is a versatile tool for rapidly and accurately identifying and visualising gene set enrichments in multi-contrast omics data. Mitch is available from Bioconductor ( https://bioconductor.org/packages/mitch ).

Genomic DNA methylation distinguishes subtypes of human focal cortical dysplasia.

OBJECTIVES: Focal cortical dysplasia (FCD) is a major cause of drug-resistant focal epilepsy in children, and the clinicopathological classification remains a challenging issue in daily practice. With the recent progress in DNA methylation-based classification of human brain tumors we examined whether genomic DNA methylation and gene expression analysis can be used to also distinguish human FCD subtypes. METHODS: DNA methylomes and transcriptomes were generated from massive parallel sequencing in 15 surgical FCD specimens, matched with 5 epilepsy and 6 nonepilepsy controls. RESULTS: Differential hierarchical cluster analysis of DNA methylation distinguished major FCD subtypes (ie, Ia, IIa, and IIb) from patients with temporal lobe epilepsy patients and nonepileptic controls. Targeted panel sequencing identified a novel likely pathogenic variant in DEPDC5 in a patient with FCD type IIa. However, no enrichment of differential DNA methylation or gene expression was observed in mechanistic target of rapamycin (mTOR) pathway-related genes. SIGNIFICANCE: Our studies extend the evidence for disease-specific methylation signatures toward focal epilepsies in favor of an integrated clinicopathologic and molecular classification system of FCD subtypes incorporating genomic methylation.

Lipoxins Regulate the Early Growth Response-1 Network and Reverse Diabetic Kidney Disease.

Background The failure of spontaneous resolution underlies chronic inflammatory conditions, including microvascular complications of diabetes such as diabetic kidney disease. The identification of endogenously generated molecules that promote the physiologic resolution of inflammation suggests that these bioactions may have therapeutic potential in the context of chronic inflammation. Lipoxins (LXs) are lipid mediators that promote the resolution of inflammation.Methods We investigated the potential of LXA4 and a synthetic LX analog (Benzo-LXA4) as therapeutics in a murine model of diabetic kidney disease, ApoE-/- mice treated with streptozotocin.Results Intraperitoneal injection of LXs attenuated the development of diabetes-induced albuminuria, mesangial expansion, and collagen deposition. Notably, LXs administered 10 weeks after disease onset also attenuated established kidney disease, with evidence of preserved kidney function. Kidney transcriptome profiling defined a diabetic signature (725 genes; false discovery rate P≤0.05). Comparison of this murine gene signature with that of human diabetic kidney disease identified shared renal proinflammatory/profibrotic signals (TNF-α, IL-1ß, NF-κB). In diabetic mice, we identified 20 and 51 transcripts regulated by LXA4 and Benzo-LXA4, respectively, and pathway analysis identified established (TGF-ß1, PDGF, TNF-α, NF-κB) and novel (early growth response-1 [EGR-1]) networks activated in diabetes and regulated by LXs. In cultured human renal epithelial cells, treatment with LXs attenuated TNF-α-driven Egr-1 activation, and Egr-1 depletion prevented cellular responses to TGF-ß1 and TNF-αConclusions These data demonstrate that LXs can reverse established diabetic complications and support a therapeutic paradigm to promote the resolution of inflammation.

Multicellular Transcriptional Analysis of Mammalian Heart Regeneration.

BACKGROUND: The inability of the adult mammalian heart to regenerate following injury represents a major barrier in cardiovascular medicine. In contrast, the neonatal mammalian heart retains a transient capacity for regeneration, which is lost shortly after birth. Defining the molecular mechanisms that govern regenerative capacity in the neonatal period remains a central goal in cardiac biology. Here, we assemble a transcriptomic framework of multiple cardiac cell populations during postnatal development and following injury, which enables comparative analyses of the regenerative (neonatal) versus nonregenerative (adult) state for the first time. METHODS: Cardiomyocytes, fibroblasts, leukocytes, and endothelial cells from infarcted and noninfarcted neonatal (P1) and adult (P56) mouse hearts were isolated by enzymatic dissociation and fluorescence-activated cell sorting at day 3 following surgery. RNA sequencing was performed on these cell populations to generate the transcriptome of the major cardiac cell populations during cardiac development, repair, and regeneration. To complement our transcriptomic data, we also surveyed the epigenetic landscape of cardiomyocytes during postnatal maturation by performing deep sequencing of accessible chromatin regions by using the Assay for Transposase-Accessible Chromatin from purified mouse cardiomyocyte nuclei (P1, P14, and P56). RESULTS: Profiling of cardiomyocyte and nonmyocyte transcriptional programs uncovered several injury-responsive genes across regenerative and nonregenerative time points. However, the majority of transcriptional changes in all cardiac cell types resulted from developmental maturation from neonatal stages to adulthood rather than activation of a distinct regeneration-specific gene program. Furthermore, adult leukocytes and fibroblasts were characterized by the expression of a proliferative gene expression network following infarction, which mirrored the neonatal state. In contrast, cardiomyocytes failed to reactivate the neonatal proliferative network following infarction, which was associated with loss of chromatin accessibility around cell cycle genes during postnatal maturation. CONCLUSIONS: This work provides a comprehensive framework and transcriptional resource of multiple cardiac cell populations during cardiac development, repair, and regeneration. Our findings define a regulatory program underpinning the neonatal regenerative state and identify alterations in the chromatin landscape that could limit reinduction of the regenerative program in adult cardiomyocytes.

Evaluation of microRNA alignment techniques.

Genomic alignment of small RNA (smRNA) sequences such as microRNAs poses considerable challenges due to their short length (∼21 nucleotides [nt]) as well as the large size and complexity of plant and animal genomes. While several tools have been developed for high-throughput mapping of longer mRNA-seq reads (>30 nt), there are few that are specifically designed for mapping of smRNA reads including microRNAs. The accuracy of these mappers has not been systematically determined in the case of smRNA-seq. In addition, it is unknown whether these aligners accurately map smRNA reads containing sequence errors and polymorphisms. By using simulated read sets, we determine the alignment sensitivity and accuracy of 16 short-read mappers and quantify their robustness to mismatches, indels, and nontemplated nucleotide additions. These were explored in the context of a plant genome (Oryza sativa, ∼500 Mbp) and a mammalian genome (Homo sapiens, ∼3.1 Gbp). Analysis of simulated and real smRNA-seq data demonstrates that mapper selection impacts differential expression results and interpretation. These results will inform on best practice for smRNA mapping and enable more accurate smRNA detection and quantification of expression and RNA editing.

Vascular histone deacetylation by pharmacological HDAC inhibition.

HDAC inhibitors can regulate gene expression by post-translational modification of histone as well as nonhistone proteins. Often studied at single loci, increased histone acetylation is the paradigmatic mechanism of action. However, little is known of the extent of genome-wide changes in cells stimulated by the hydroxamic acids, TSA and SAHA. In this article, we map vascular chromatin modifications including histone H3 acetylation of lysine 9 and 14 (H3K9/14ac) using chromatin immunoprecipitation (ChIP) coupled with massive parallel sequencing (ChIP-seq). Since acetylation-mediated gene expression is often associated with modification of other lysine residues, we also examined H3K4me3 and H3K9me3 as well as changes in CpG methylation (CpG-seq). RNA sequencing indicates the differential expression of ∼30% of genes, with almost equal numbers being up- and down-regulated. We observed broad deacetylation and gene expression changes conferred by TSA and SAHA mediated by the loss of EP300/CREBBP binding at multiple gene promoters. This study provides an important framework for HDAC inhibitor function in vascular biology and a comprehensive description of genome-wide deacetylation by pharmacological HDAC inhibition.

Dynamic changes in the cardiac methylome during postnatal development.

Relatively little is known about the epigenetic control mechanisms that guide postnatal organ maturation. The goal of this study was to determine whether DNA methylation plays an important role in guiding transcriptional changes during the first 2 wk of mouse heart development, which is an important period for cardiomyocyte maturation, loss of proliferative capacity and loss of regenerative potential. Gene expression profiling (RNA-seq) and genome-wide sequencing of methylated DNA (MBD-seq) identified dynamic changes in the cardiac methylome during postnatal development [2545 differentially methylated regions (DMRs) from P1 to P14 in the mouse]. The vast majority (~80%) of DMRs were hypermethylated between P1 and P14, and these hypermethylated regions were associated with transcriptional shut down of important developmental signaling pathways, including Hedgehog, bone morphogenetic protein, TGF-ß, fibroblast growth factor, and Wnt/ß-catenin signaling. Postnatal inhibition of DNA methylation with 5-aza-2'-deoxycytidine induced a marked increase (~3-fold) in cardiomyocyte proliferation and ~50% reduction in the percentage of binucleated cardiomyocytes compared with saline-treated controls. This study provides novel evidence for widespread alterations in DNA methylation during postnatal heart maturation and suggests that cardiomyocyte cell cycle arrest during the neonatal period is subject to regulation by DNA methylation.

Genome-wide analysis distinguishes hyperglycemia regulated epigenetic signatures of primary vascular cells.

Emerging evidence suggests that poor glycemic control mediates post-translational modifications to the H3 histone tail. We are only beginning to understand the dynamic role of some of the diverse epigenetic changes mediated by hyperglycemia at single loci, yet elevated glucose levels are thought to regulate genome-wide changes, and this still remains poorly understood. In this article we describe genome-wide histone H3K9/K14 hyperacetylation and DNA methylation maps conferred by hyperglycemia in primary human vascular cells. Chromatin immunoprecipitation (ChIP) as well as CpG methylation (CpG) assays, followed by massive parallel sequencing (ChIP-seq and CpG-seq) identified unique hyperacetylation and CpG methylation signatures with proximal and distal patterns of regionalization associative with gene expression. Ingenuity knowledge-based pathway and gene ontology analyses indicate that hyperglycemia significantly affects human vascular chromatin with the transcriptional up-regulation of genes involved in metabolic and cardiovascular disease. We have generated the first installment of a reference collection of hyperglycemia-induced chromatin modifications using robust and reproducible platforms that allow parallel sequencing-by-synthesis of immunopurified content. We uncover that hyperglycemia-mediated induction of genes and pathways associated with endothelial dysfunction occur through modulation of acetylated H3K9/K14 inversely correlated with methyl-CpG content.

Epigenetic regulation of cell type-specific expression patterns in the human mammary epithelium.

Differentiation is an epigenetic program that involves the gradual loss of pluripotency and acquisition of cell type-specific features. Understanding these processes requires genome-wide analysis of epigenetic and gene expression profiles, which have been challenging in primary tissue samples due to limited numbers of cells available. Here we describe the application of high-throughput sequencing technology for profiling histone and DNA methylation, as well as gene expression patterns of normal human mammary progenitor-enriched and luminal lineage-committed cells. We observed significant differences in histone H3 lysine 27 tri-methylation (H3K27me3) enrichment and DNA methylation of genes expressed in a cell type-specific manner, suggesting their regulation by epigenetic mechanisms and a dynamic interplay between the two processes that together define developmental potential. The technologies we developed and the epigenetically regulated genes we identified will accelerate the characterization of primary cell epigenomes and the dissection of human mammary epithelial lineage-commitment and luminal differentiation.

Inherited PURA Pathogenic Variant Associated With a Mild Neurodevelopmental Disorder.

Objectives: Purine-rich element-binding protein alpha (PURA) regulates gene expression and is ubiquitously expressed with an enrichment in neural tissues. Pathogenic variants in PURA cause the neurodevelopmental disorder PURA syndrome that has a variable phenotype but typically comprises moderate-to-severe global developmental delay, intellectual disability, early-onset hypotonia and hypothermia, epilepsy, feeding difficulties, movement disorders, and subtle facial dysmorphism. Speech is reportedly absent in most, but the specific linguistic phenotype is not well described. Methods: We used genome sequencing to identify a pathogenic gene variant as part of a study of children ascertained for severe primary speech disorder in the absence of moderate or severe ID. Results: The novel PURA c.296G>T (p.Arg99Leu) pathogenic missense variant segregated in the female proband and her affected mother. The proband had dysarthria; phonological disorder; and severe receptive and expressive language impairment, borderline intellect, attention difficulties, oropharyngeal dysmotility, and dysmorphic facial features. Her mother had dysarthria, moderate receptive language impairment, and borderline intellect. Both the proband and her mother completed mainstream schooling with classroom support. Discussion: This is the first inherited PURA pathogenic germline variant in over 600 unrelated families documented on ClinVar or reported in the literature. PURA testing should be considered in families with primary speech disorder and borderline intellectual disability, given the specific genetic counseling implications.

Creating reusable tools from scripts: the Galaxy Tool Factory.

MOTIVATION: Galaxy is a software application supporting high-throughput biology analyses and work flows, available as a free on-line service or as source code for local deployment. New tools can be written to extend Galaxy, and these can be shared using public Galaxy Tool Shed (GTS) repositories, but converting even simple scripts into tools requires effort from a skilled developer. RESULTS: The Tool Factory is a novel Galaxy tool that automates the generation of all code needed to execute user-supplied scripts, and wraps them into new Galaxy tools for upload to a GTS, ready for review and installation through the Galaxy administrative interface. AVAILABILITY AND IMPLEMENTATION: The Galaxy administrative interface supports automated installation from the main GTS. Source code and support are available at the project website, https://bitbucket.org/fubar/galaxytoolfactory. The Tool Factory is implemented as an installable Galaxy tool. CONTACT: ross.lazarus@channing.harvard.edu.

Deep sequencing reveals increased DNA methylation in chronic rat epilepsy.

Epilepsy is a frequent neurological disorder, although onset and progression of seizures remain difficult to predict in affected patients, irrespective of their epileptogenic condition. Previous studies in animal models as well as human epileptic brain tissue revealed a remarkably diverse pattern of gene expression implicating epigenetic changes to contribute to disease progression. Here we mapped for the first time global DNA methylation patterns in chronic epileptic rats and controls. Using methyl-CpG capture associated with massive parallel sequencing (Methyl-Seq) we report the genomic methylation signature of the chronic epileptic state. We observed a predominant increase, rather than loss of DNA methylation in chronic rat epilepsy. Aberrant methylation patterns were inversely correlated with gene expression changes using mRNA sequencing from same animals and tissue specimens. Administration of a ketogenic, high-fat, low-carbohydrate diet attenuated seizure progression and ameliorated DNA methylation mediated changes in gene expression. This is the first report of unsupervised clustering of an epigenetic mark being used in epilepsy research to separate epileptic from non-epileptic animals as well as from animals receiving anti-convulsive dietary treatment. We further discuss the potential impact of epigenetic changes as a pathogenic mechanism of epileptogenesis.

Reduced methylation correlates with diabetic nephropathy risk in type 1 diabetes.

Diabetic nephropathy (DN) is a polygenic disorder with few risk variants showing robust replication in large-scale genome-wide association studies. To understand the role of DNA methylation, it is important to have the prevailing genomic view to distinguish key sequence elements that influence gene expression. This is particularly challenging for DN because genome-wide methylation patterns are poorly defined. While methylation is known to alter gene expression, the importance of this causal relationship is obscured by array-based technologies since coverage outside promoter regions is low. To overcome these challenges, we performed methylation sequencing using leukocytes derived from participants of the Finnish Diabetic Nephropathy (FinnDiane) type 1 diabetes (T1D) study (n = 39) that was subsequently replicated in a larger validation cohort (n = 296). Gene body-related regions made up more than 60% of the methylation differences and emphasized the importance of methylation sequencing. We observed differentially methylated genes associated with DN in 3 independent T1D registries originating from Denmark (n = 445), Hong Kong (n = 107), and Thailand (n = 130). Reduced DNA methylation at CTCF and Pol2B sites was tightly connected with DN pathways that include insulin signaling, lipid metabolism, and fibrosis. To define the pathophysiological significance of these population findings, methylation indices were assessed in human renal cells such as podocytes and proximal convoluted tubule cells. The expression of core genes was associated with reduced methylation, elevated CTCF and Pol2B binding, and the activation of insulin-signaling phosphoproteins in hyperglycemic cells. These experimental observations also closely parallel methylation-mediated regulation in human macrophages and vascular endothelial cells.

Combining target enrichment with barcode multiplexing for high throughput SNP discovery.

BACKGROUND: The primary goal of genetic linkage analysis is to identify genes affecting a phenotypic trait. After localisation of the linkage region, efficient genetic dissection of the disease linked loci requires that functional variants are identified across the loci. These functional variations are difficult to detect due to extent of genetic diversity and, to date, incomplete cataloguing of the large number of variants present both within and between populations. Massively parallel sequencing platforms offer unprecedented capacity for variant discovery, however the number of samples analysed are still limited by cost per sample. Some progress has been made in reducing the cost of resequencing using either multiplexing methodologies or through the utilisation of targeted enrichment technologies which provide the ability to resequence genomic areas of interest rather that full genome sequencing. RESULTS: We developed a method that combines current multiplexing methodologies with a solution-based target enrichment method to further reduce the cost of resequencing where region-specific sequencing is required. Our multiplex/enrichment strategy produced high quality data with nominal reduction of sequencing depth. We undertook a genotyping study and were successful in the discovery of novel SNP alleles in all samples at uniplex, duplex and pentaplex levels. CONCLUSION: Our work describes the successful combination of a targeted enrichment method and index barcode multiplexing to reduce costs, time and labour associated with processing large sample sets. Furthermore, we have shown that the sequencing depth obtained is adequate for credible SNP genotyping analysis at uniplex, duplex and pentaplex levels.