A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

Regional reemergence of a SARS-CoV-2 Delta lineage amid an Omicron wave detected by wastewater sequencing

The implementation and integration of wastewater-based epidemiology constitutes a valuable addition to existing pathogen surveillance systems, such as clinical surveillance for SARS-CoV-2. In the Netherlands, SARS-CoV-2 variant circulation is monitored by performing whole-genome sequencing on wastewater samples. In this manuscript, we describe the detection of an AY.43 lineage (Delta variant) amid a period of BA.5 (Omicron variant) dominance in wastewater samples from two wastewater treatment plants (WWTPs) during the months of August and September of 2022. Our results describe a temporary emergence, which was absent in samples from other WWTPs, and which coincided with peaks in viral load. We show how these lineage estimates can be traced back to lineage-specific substitution patterns. The absence of this variant from reported clinical data, but high associated viral loads suggest cryptic transmission. Our findings highlight the additional value of wastewater surveillance for generating insights into circulating pathogens.

Identification of a Novel FOXP1 Variant in a Patient with Hypotonia, Intellectual Disability, and Severe Speech Impairment

The FOXP subfamily includes four different transcription factors: FOXP1, FOXP2, FOXP3, and FOXP4, all with important roles in regulating gene expression from early development through adulthood. Haploinsufficiency of FOXP1, due to deleterious variants (point mutations, copy number variants) disrupting the gene, leads to an emerging disorder known as "FOXP1 syndrome", mainly characterized by intellectual disability, language impairment, dysmorphic features, and multiple congenital abnormalities with or without autistic features in some affected individuals (MIM 613670). Here we describe a 10-year-old female patient, born to unrelated parents, showing hypotonia, intellectual disability, and severe language delay. Targeted resequencing analysis allowed us to identify a heterozygous de novo FOXP1 variant c.1030C>T, p.(Gln344Ter) classified as likely pathogenetic according to the American College of Medical Genetics and Genomics guidelines. To the best of our knowledge, our patient is the first to date to report carrying this stop mutation, which is, for this reason, useful for broadening the molecular spectrum of FOXP1 clinically relevant variants. In addition, our results highlight the utility of next-generation sequencing in establishing an etiological basis for heterogeneous conditions such as neurodevelopmental disorders and providing additional insight into the phenotypic features of FOXP1-related syndrome.

capTEs enables locus-specific dissection of transcriptional outputs from reference and nonreference transposable elements

Transposable elements (TEs) serve as both insertional mutagens and regulatory elements in cells, and their aberrant activity is increasingly being revealed to contribute to diseases and cancers. However, measuring the transcriptional consequences of nonreference and young TEs at individual loci remains challenging with current methods, primarily due to technical limitations, including short read lengths generated and insufficient coverage in target regions. Here, we introduce a long-read targeted RNA sequencing method, Cas9-assisted profiling TE expression sequencing (capTEs), for quantitative analysis of transcriptional outputs for individual TEs, including transcribed nonreference insertions, noncanonical transcripts from various transcription patterns and their correlations with expression changes in related genes. This method selectively identified TE-containing transcripts and outputted data with up to 90% TE reads, maintaining a comparable data yield to whole-transcriptome sequencing. We applied capTEs to human cancer cells and found that internal and inserted Alu elements may employ distinct regulatory mechanisms to upregulate gene expression. We expect that capTEs will be a critical tool for advancing our understanding of the biological functions of individual TEs at the locus level, revealing their roles as both mutagens and regulators in biological and pathogenic processes.

Targeted Resequencing of Otosclerosis Patients from Different Populations Replicates Results from a Previous Genome-Wide Association Study

Otosclerosis is one of the most common causes of hearing loss in young adults. It has a prevalence of 0.3-0.4% in the European population. Clinical symptoms usually occur between the second and fifth decade of life. Different studies have been performed to unravel the genetic architecture of the disease. Recently, a genome-wide association study (GWAS) identified 15 novel risk loci and replicated the regions of three previously reported candidate genes. In this study, seven candidate genes from the GWAS were resequenced using single molecule molecular inversion probes (smMIPs). smMIPs were used to capture the exonic regions and the 3' and 5' untranslated regions (UTR). Discovered variants were tested for association with the disease using single variant and gene-based association analysis. The single variant results showed that 13 significant variants were associated with otosclerosis. Associated variants were found in five of the seven genes studied here, including AHSG, LINC01482, MARK3, SUPT3H and RELN. Conversely, burden testing did not show a major role of rare variants in the disease. In conclusion, this study was able to replicate five out of seven candidate genes reported in the previous GWAS. This association is likely mainly driven by common variants.

Association of soluble epoxide hydrolase 2 gene with the risk of non-syndromic cleft lip with or without cleft palate in western Han Chinese population

OBJECTIVES

This study aimed to explore the associations between soluble epoxide hydrolase 2 gene (EPHX2) variants and non-syndromic cleft lip with or without cleft palate (NSCL/P) in Chinese Han population.

METHODS

We recruited 159 NSCL/P cases from Chinese Han population and carried out targeted resequencing using the whole genome sequencing data of 542 healthy Chinese individuals from Novegene internal database as controls. We classified EPHX2 variants as common or rare according to their minor allele frequency and performed an association analysis for common variations and a burden analysis for rare variations.

RESULTS

The lowest P-value in NSCL/P was observed at rs57699806 (P=0.000 13, OR=2.849 and 95% CI: 1.691-4.800), followed by rs4732723 (P=0.006 50, OR=0.662 and 95%CI: 0.491-0.892), rs7829267 (P=0.009 20, OR=1.496 and 95%CI: 1.117-2.005), rs721619 (P=0.011 00, OR=1.474 and 95%CI: 1.098-1.980), and rs7816586 (P=0.040 00, OR=1.310 and 95%CI: 1.015-1.691). The odds ratios suggested the C allele at rs4732723 as a protective factor for NSCL/P and the reference alleles at other single nucleotide polymorphisms (SNPs) as the risk factors for NSCL/P. Burden analysis showed no statistical significance (P>0.05).

CONCLUSIONS

Through targeted resequencing, this study identified five SNPs named rs57699806, rs4732723, rs7829267, rs721619, and rs7816586 around the region of EPHX2 gene associated with NSCL/P in Chinese Han population. Four SNPs of rs57699806, rs4732723, rs7829267, and rs7816586 were first identified.

Deep Resequencing of 9 Candidate Genes Identifies a Role for ARAP1 and IGF2BP2 in Modulating Insulin Secretion Adjusted for Insulin Resistance in Obese Southern Europeans

Type 2 diabetes is characterized by impairment in insulin secretion, with an established genetic contribution. We aimed to evaluate common and low-frequency (1–5%) variants in nine genes strongly associated with insulin secretion by targeted sequencing in subjects selected from the extremes of insulin release measured by the disposition index. Collapsing data by gene and/or function, the association between disposition index and nonsense variants were significant, also after adjustment for confounding factors (OR = 0.25, 95% CI = 0.11–0.59, p = 0.001). Evaluating variants individually, three novel variants in ARAP1, IGF2BP2 and GCK, out of eight reaching significance singularly, remained associated after adjustment. Constructing a genetic risk model combining the effects of the three variants, only carriers of the ARAP1 and IGF2BP2 variants were significantly associated with a reduced probability to be in the lower, worst, extreme of insulin secretion (OR = 0.223, 95% CI = 0.105–0.473, p < 0.001). Observing a high number of normal glucose tolerance between carriers, a regression posthoc analysis was performed. Carriers of genetic risk model variants had higher probability to be normoglycemic, also after adjustment (OR = 2.411, 95% CI = 1.136–5.116, p = 0.022). Thus, in our southern European cohort, nonsense variants in all nine candidate genes showed association with better insulin secretion adjusted for insulin resistance, and we established the role of ARAP1 and IGF2BP2 in modulating insulin secretion.

Sequence Variant in the TRIM39-RPP21 Gene Readthrough is Shared Across a Cohort of Arabian Foals Diagnosed with Juvenile Idiopathic Epilepsy

Juvenile idiopathic epilepsy (JIE) is a self-limiting neurological disorder with a suspected genetic predisposition affecting young Arabian foals of the Egyptian lineage. The condition is characterized by tonic-clonic seizures with intermittent post-ictal blindness, in which most incidents are sporadic and unrecognized. This study aimed to identify genetic components shared across a local cohort of Arabian foals diagnosed with JIE via a combined whole genome and targeted resequencing approach: Initial whole genome comparisons between a small cohort of nine diagnosed foals (cases) and 27 controls from other horse breeds identified variants uniquely shared amongst the case cohort. Further validation via targeted resequencing of these variants, that pertain to non-intergenic regions, on additional eleven case individuals revealed a single 19bp deletion coupled with a triple-C insertion (Δ19InsCCC) within the TRIM39-RPP21 gene readthrough that was uniquely shared across all case individuals, and absent from three additional Arabian controls. Furthermore, we have confirmed recent findings refuting potential linkage between JIE and other inherited diseases in the Arabian lineage, and refuted the potential linkage between JIE and genes predisposing a similar disorder in human newborns. This is the first study to report a genetic variant to be shared in a sub-population cohort of Arabian foals diagnosed with JIE. Further evaluation of the sensitivity and specificity of the Δ19InsCCC allele within additional cohorts of the Arabian horse is warranted in order to validate its credibility as a marker for JIE, and to ascertain whether it has been introduced into other horse breeds by Arabian ancestry.

Multi-omics mapping of human papillomavirus integration sites illuminates novel cervical cancer target genes

BACKGROUND

Integration of human papillomavirus (HPV) into the host genome is a dominant feature of invasive cervical cancer (ICC), yet the tumorigenicity of cis genomic changes at integration sites remains largely understudied.

METHODS

Combining multi-omics data from The Cancer Genome Atlas with patient-matched long-read sequencing of HPV integration sites, we developed a strategy for using HPV integration events to identify and prioritise novel candidate ICC target genes (integration-detected genes (IDGs)). Four IDGs were then chosen for in vitro functional studies employing small interfering RNA-mediated knockdown in cell migration, proliferation and colony formation assays.

RESULTS

PacBio data revealed 267 unique human-HPV breakpoints comprising 87 total integration events in eight tumours. Candidate IDGs were filtered based on the following criteria: (1) proximity to integration site, (2) clonal representation of integration event, (3) tumour-specific expression (Z-score) and (4) association with ICC survival. Four candidates prioritised based on their unknown function in ICC (BNC1, RSBN1, USP36 and TAOK3) exhibited oncogenic properties in cervical cancer cell lines. Further, annotation of integration events provided clues regarding potential mechanisms underlying altered IDG expression in both integrated and non-integrated ICC tumours.

CONCLUSIONS

HPV integration events can guide the identification of novel IDGs for further study in cervical carcinogenesis and as putative therapeutic targets.

Contribution of rare variant associations to neurodegenerative disease presentation

Genetic factors contribute to neurodegenerative diseases, with high heritability estimates across diagnoses; however, a large portion of the genetic influence remains poorly understood. Many previous studies have attempted to fill the gaps by performing linkage analyses and association studies in individual disease cohorts, but have failed to consider the clinical and pathological overlap observed across neurodegenerative diseases and the potential for genetic overlap between the phenotypes. Here, we leveraged rare variant association analyses (RVAAs) to elucidate the genetic overlap among multiple neurodegenerative diagnoses, including Alzheimer's disease, amyotrophic lateral sclerosis, frontotemporal dementia (FTD), mild cognitive impairment, and Parkinson's disease (PD), as well as cerebrovascular disease, using the data generated with a custom-designed neurodegenerative disease gene panel in the Ontario Neurodegenerative Disease Research Initiative (ONDRI). As expected, only ~3% of ONDRI participants harboured a monogenic variant likely driving their disease presentation. Yet, when genes were binned based on previous disease associations, we observed an enrichment of putative loss of function variants in PD genes across all ONDRI cohorts. Further, individual gene-based RVAA identified significant enrichment of rare, nonsynonymous variants in PARK2 in the FTD cohort, and in NOTCH3 in the PD cohort. The results indicate that there may be greater heterogeneity in the genetic factors contributing to neurodegeneration than previously appreciated. Although the mechanisms by which these genes contribute to disease presentation must be further explored, we hypothesize they may be a result of rare variants of moderate phenotypic effect contributing to overlapping pathology and clinical features observed across neurodegenerative diagnoses.

Analysis of HLA gene polymorphisms in East Africans reveals evidence of gene flow in two Semitic populations from Sudan

Sudan, a northeastern African country, is characterized by high levels of cultural, linguistic, and genetic diversity, which is believed to be affected by continuous migration from neighboring countries. Consistent with such demographic effect, genome-wide SNP data revealed a shared ancestral component among Sudanese Afro-Asiatic speaking groups and non-African populations, mainly from West Asia. Although this component is shared among all Afro-Asiatic speaking groups, the extent of this sharing in Semitic groups, such as Sudanese Arab, is still unknown. Using genotypes of six polymorphic human leukocyte antigen (HLA) genes (i.e., HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1), we examined the genetic structure of eight East African ethnic groups with origins in Sudan, South Sudan, and Ethiopia. We identified informative HLA alleles using principal component analysis, which revealed that the two Semitic groups (Gaalien and Shokrya) constituted a distinct cluster from the other Afro-Asiatic speaking groups in this study. The HLA alleles that distinguished Semitic Arabs co-exist in the same extended HLA haplotype, and those alleles are in strong linkage disequilibrium. Interestingly, we find the four-locus haplotype "C*12:02-B*52:01-DRB1*15:02-DQB1*06:01" exclusively in non-African populations and it is widely spread across Asia. The identification of this haplotype suggests a gene flow from Asia, and likely these haplotypes were brought to Africa through back migration from the Near East. These findings will be of interest to biomedical and anthropological studies that examine the demographic history of northeast Africa.

Interspecific introgression and natural selection in the evolution of Japanese apricot (Prunus mume)

Domestication and population differentiation in crops involve considerable phenotypic changes. The logs of these evolutionary paths, including natural/artificial selection, can be found in the genomes of the current populations. However, these profiles have been little studied in tree crops, which have specific characters, such as long generation time and clonal propagation, maintaining high levels of heterozygosity. We conducted exon-targeted resequencing of 129 genomes in the genus Prunus, mainly Japanese apricot (Prunus mume), and apricot (Prunus armeniaca), plum (Prunus salicina), and peach (Prunus persica). Based on their genome-wide single-nucleotide polymorphisms merged with published resequencing data of 79 Chinese P. mume cultivars, we inferred complete and ongoing population differentiation in P. mume. Sliding window characterization of the indexes for genetic differentiation identified interspecific fragment introgressions between P. mume and related species (plum and apricot). These regions often exhibited strong selective sweeps formed in the paths of establishment or formation of substructures of P. mume, suggesting that P. mume has frequently imported advantageous genes from other species in the subgenus Prunus as adaptive evolution. These findings shed light on the complicated nature of adaptive evolution in a tree crop that has undergone interspecific exchange of genome fragments with natural/artificial selections.

Use of Targeted Amplicon Sequencing in Peanut to Generate Allele Information on Allotetraploid Sub-Genomes

The use of molecular markers in plant breeding has become a routine practice, but the cost per accession can be a hindrance to the routine use of Quantitative Trait Loci (QTL) identification in breeding programs. In this study, we demonstrate the use of targeted re-sequencing as a proof of concept of a cost-effective approach to retrieve highly informative allele information, as well as develop a bioinformatics strategy to capture the genome-specific information of a polyploid species. SNPs were identified from alignment of raw transcriptome reads (2 × 50 bp) to a synthetic tetraploid genome using BWA followed by a GATK pipeline. Regions containing high polymorphic SNPs in both A genome and B genomes were selected as targets for the resequencing study. Targets were amplified using multiplex PCR followed by sequencing on an Illumina HiSeq. Eighty-one percent of the SNP calls in diploids and 68% of the SNP calls in tetraploids were confirmed. These results were also confirmed by KASP validation. Based on this study, we find that targeted resequencing technologies have potential for obtaining maximum allele information in allopolyploids at reduced cost.

Circulating tumor cells and breast cancer-specific mutations in primary breast cancer

Circulating tumor cells (CTCs) play a pivotal role in tumor dissemination and progression, and are considered to be a critical part of the metastatic cascade. The aim of the present research article was to examine breast cancer-specific mutations in primary breast cancer (PBC) using targeted resequencing. A total of 78 patients with PBC were enrolled into this translational study. Reverse transcription-quantitative PCR assay for the expression of epithelial markers (CK19) or epithelial-to-mesenchymal transition (EMT)-related genes (TWIST1, SNAIL1, SLUG and ZEB1) was applied for identification of CTCs prior to surgery. Total DNA was isolated from fresh frozen primary tumors. Sequencing was performed by Agilent SureSelect target enrichment and Illumina paired-end sequencing on the MiSeq platform. The most commonly affected genes were TP53 (mutated in 21 tumors; 26.9%), followed by PIK3CA (mutated in 16 tumors; 20.5%) and BRCA1/2 (mutated in 7 tumors, BRCA1 n=2 and BRCA2 n=5; 9.0%). In our cohort, a significantly higher proportion of patients with epithelial CTCs harbored mutations in the BRCA1/2 genes in the tumor tissue. There were no mutations in specific genes associated with CTCs with the EMT phenotype. To the best of our knowledge, this study is the first to report a correlation between the presence of epithelial CTCs in the peripheral blood and mutations of the BRCA1/2 genes in primary tumor tissue.

Resequencing Study Confirms That Host Defense and Cell Senescence Gene Variants Contribute to the Risk of Idiopathic Pulmonary Fibrosis

Rationale: Several common and rare genetic variants have been associated with idiopathic pulmonary fibrosis, a progressive fibrotic condition that is localized to the lung. Objectives: To develop an integrated understanding of the rare and common variants located in multiple loci that have been reported to contribute to the risk of disease. Methods: We performed deep targeted resequencing (3.69 Mb of DNA) in cases (n = 3,624) and control subjects (n = 4,442) across genes and regions previously associated with disease. We tested for associations between disease and 1) individual common variants via logistic regression and 2) groups of rare variants via sequence kernel association tests. Measurements and Main Results: Statistically significant common variant association signals occurred in all 10 of the regions chosen based on genome-wide association studies. The strongest risk variant is the MUC5B promoter variant rs35705950, with an odds ratio of 5.45 (95% confidence interval, 4.91-6.06) for one copy of the risk allele and 18.68 (95% confidence interval, 13.34-26.17) for two copies of the risk allele (P = 9.60 × 10-295). In addition to identifying for the first time that rare variation in FAM13A is associated with disease, we confirmed the role of rare variation in the TERT and RTEL1 gene regions in the risk of IPF, and found that the FAM13A and TERT regions have independent common and rare variant signals. Conclusions: A limited number of common and rare variants contribute to the risk of idiopathic pulmonary fibrosis in each of the resequencing regions, and these genetic variants focus on biological mechanisms of host defense and cell senescence.

Targeted resequencing of 358 candidate genes for autism spectrum disorder in a Chinese cohort reveals diagnostic potential and genotype-phenotype correlations

Autism spectrum disorder (ASD) is a childhood neuropsychiatric disorder with a complex genetic architecture. The diagnostic potential of a targeted panel of ASD genes has only been evaluated in small cohorts to date and is especially understudied in the Chinese population. Here, we designed a capture panel with 358 genes (111 syndromic and 247 nonsyndromic) for ASD and sequenced a Chinese cohort of 539 cases evaluated with the Autism Diagnostic Interview‐Revised (ADI‐R) and the Autism Diagnostic Observation Schedule (ADOS) as well as 512 controls. ASD cases were found to carry significantly more ultra‐rare functional variants than controls. A subset of 78 syndromic and 54 nonsyndromic genes was the most significantly associated and should be given high priority in the future screening of ASD patients. Pathogenic and likely pathogenic variants were detected in 9.5% of cases. Variants in SHANK3 and SHANK2 were the most frequent, especially in females, and occurred in 1.2% of cases. Duplications of 15q11–13 were detected in 0.8% of cases. Variants in CNTNAP2 and MEF2C were correlated with epilepsy/tics in cases. Our findings reveal the diagnostic potential of ASD genetic panel testing and new insights regarding the variant spectrum. Genotype–phenotype correlations may facilitate the diagnosis and management of ASD.

Overcoming challenges in variant calling: exploring sequence diversity in candidate genes for plant development in perennial ryegrass (Lolium perenne)

Revealing DNA sequence variation within the Lolium perenne genepool is important for genetic analysis and development of breeding applications. We reviewed current literature on plant development to select candidate genes in pathways that control agronomic traits, and identified 503 orthologues in L. perenne. Using targeted resequencing, we constructed a comprehensive catalogue of genomic variation for a L. perenne germplasm collection of 736 genotypes derived from current cultivars, breeding material and wild accessions. To overcome challenges of variant calling in heterogeneous outbreeding species, we used two complementary strategies to explore sequence diversity. First, four variant calling pipelines were integrated with the VariantMetaCaller to reach maximal sensitivity. Additional multiplex amplicon sequencing was used to empirically estimate an appropriate precision threshold. Second, a de novo assembly strategy was used to reconstruct divergent alleles for each gene. The advantage of this approach was illustrated by discovery of 28 novel alleles of LpSDUF247, a polymorphic gene co-segregating with the S-locus of the grass self-incompatibility system. Our approach is applicable to other genetically diverse outbreeding species. The resulting collection of functionally annotated variants can be mined for variants causing phenotypic variation, either through genetic association studies, or by selecting carriers of rare defective alleles for physiological analyses.

Discovery of rare variants implicated in schizophrenia using next-generation sequencing

Schizophrenia is a highly heritable psychiatric disorder that affects 1% of the population. Genome-wide association studies have identified common variants in candidate genes associated with schizophrenia, but the genetics mechanisms of this disorder have not yet been elucidated. The discovery of rare genetic variants that contribute to schizophrenia symptoms promises to help explain the missing heritability of the disease. Next generation sequencing techniques are revolutionizing the field of psychiatric genetics. Various statistical approaches have been developed for rare variant association testing in case-control and family studies. Targeted resequencing, whole exome sequencing and whole genome sequencing combined with these computational tools are used for the discovery of rare genetic variations in schizophrenia. The findings provide useful information for characterizing the rare mutations and elucidating the genetic mechanisms by which the variants cause schizophrenia.

Novel SYNGAP1 variant in a patient with intellectual disability and distinctive dysmorphisms

We describe a novel de novo heterozygous variant in SYNGAP1 (c.1741C>T, p.R581W), identified through targeted resequencing in an 8-year-old boy with intellectual disability, autism spectrum disorder, distinctive dysmorphic features, and no seizures. Our data strongly suggest that the SYNGAP1 variant is causative of intellectual disability in this patient.

Molecular Genetic Characterization of Patients With Focal Epilepsy Using a Customized Targeted Resequencing Gene Panel

Objective: Focal epilepsy is the most common subtype of epilepsies in which the influence of underlying genetic factors is emerging but remains largely uncharacterized. The purpose of this study is to determine the contribution of currently known disease-causing genes in a large cohort (n = 593) of common focal non-lesional epilepsy patients.

Methods: The customized focal epilepsy gene panel (21 genes) was based on multiplex polymerase chain reaction (PCR) and sequenced by Illumina MiSeq platform.

Results: Eleven variants (1.85%) were considered as pathogenic or likely pathogenic, including seven novel mutations. There were three SCN1A (p.Leu890Pro, p.Arg1636Ter, and p.Met1714Val), three PRRT2 (two p.Arg217Profs^*8 and p.Leu298Pro), two CHRNA4 (p.Ser284Leu, p.Ile321Asn), one DEPDC5 (p.Val516Ter), one PCDH19 (p.Asp233Asn), and one SLC2A1 (p.Ser414Ter) variants. Additionally, 16 other rare variants were classified as unknown significance due to inconsistent phenotype or lack of segregation data.

Conclusion: Currently known focal epilepsy genes only explained a very small subset of focal epilepsy patients. This indicates that the underlying genetic architecture of focal epilepsies is very heterogeneous and more novel genes are likely to be discovered. Our study highlights the usefulness, challenges and limitations of using the multi-gene panel as a diagnostic test in routine clinical practice in patients with focal epilepsy.

GSTP1 rs1695 is associated with both hematological toxicity and prognosis of ovarian cancer treated with paclitaxel plus carboplatin combination chemotherapy: a comprehensive analysis using targeted resequencing of 100 pharmacogenes

Purpose

To find genetic variants that predicted toxicity and/or efficacy of paclitaxel plus carboplatin combination therapy (TC therapy).

Patients and methods

In a retrospective case-control study, we analyzed 320 patients who had received TC therapy for gynecological cancers (ovarian, fallopian tube, peritoneal, uterine, and cervical cancers) and collected their germline DNA. We performed a comprehensive pharmacogenomic analysis using a targeted resequencing panel of 100 pharmacogenes. For 1,013 variants passing QC, case-control association studies and survival analyses were conducted.

Results

GSTP1 rs1695 showed the smallest p value for hematotoxicity association, and the ¹⁰⁵Ile wild type allele had a significantly higher risk of severe hematotoxicity (neutropenia G4, thrombocytopenia ≥ G3 and anemia ≥ G3) than the ¹⁰⁵Val allele (p=0.00034, odds ratio=5.71 (95% confidence interval:1.77-18.44)). Next, we assessed 5-year progression-free survival (PFS) and overall survival (OS) in 56 advanced ovarian cancer patients who received tri-weekly TC as a first-line chemotherapy. Patients with the ¹⁰⁵Ile/¹⁰⁵Ile genotype showed significantly better PFS (p=0.00070) and OS (p=0.0012) than those with the ¹⁰⁵Ile/¹⁰⁵Val or ¹⁰⁵Val/¹⁰⁵Val genotype.

Conclusion

Our study indicates that the GSTP1 rs1695 ¹⁰⁵Ile/¹⁰⁵Ile genotype is associated with both severe hematotoxicity and high efficacy of TC therapy, identifying a possible prognostic indicator for patients with TC therapy.

Targeted resequencing reveals genomic signatures of barley domestication

Barley (Hordeum vulgare) is an established model to study domestication of the Fertile Crescent cereals. Recent molecular data suggested that domesticated barley genomes consist of the ancestral blocks descending from multiple wild barley populations. However, the relationship between the mosaic ancestry patterns and the process of domestication itself remained unclear. To address this knowledge gap, we identified candidate domestication genes using selection scans based on targeted resequencing of 433 wild and domesticated barley accessions. We conducted phylogenetic, population structure, and ancestry analyses to investigate the origin of the domesticated barley haplotypes separately at the neutral and candidate domestication loci. We discovered multiple selective sweeps that occurred on all barley chromosomes during domestication in the background of several ancestral wild populations. The ancestry analyses demonstrated that, although the ancestral blocks of the domesticated barley genomes were descended from all over the Fertile Crescent, the candidate domestication loci originated specifically in its eastern and western parts. These findings provided the first molecular evidence implicating multiple wild or protodomesticated lineages in the process of barley domestication initiated in the Levantine and Zagros clusters of the origin of agriculture.

Discovery of MYH14 as an important and unique deafness gene causing prelingually severe autosomal dominant nonsyndromic hearing loss

BACKGROUND

Pathogenic variants of MYH14 are known to be associated (in either a syndromic or nonsyndromic manner) with hearing loss. Interestingly, all reported cases to date of MYH14-related nonsyndromic hearing loss with detailed phenotypes have demonstrated mild-to-moderate progressive hearing loss with postlingual onset.

METHODS

In the present study, targeted resequencing (TRS) of known deafness genes was performed to identify the causative variant in two multiplex families segregating autosomal dominant (AD) inherited hearing loss.

RESULTS

TRS uncovered two novel variants of MYH14 (c.A572G: p.Asp191Gly in the myosin head domain and c.C73T:p.Gln25* in exon 2) from two multiplex deafness Korean families. Notably, both probands showed phenotypes of congenital or prelingual severe hearing loss. It is remarkably uncommon to encounter such a severe-to-profound, prelingual, AD hearing loss. Given that the first variant, p. Asp191Gly, was the first documented missense allele discovered in the myosin head domain of this gene related to either congenital or prelingual severe nonsyndromic hearing loss, and also that the second variant, p. Gln25*, lead to a null allele, more severe phenotypes from our probands may have been the result of either genotype-phenotype correlation or genetic backgrounds, or both.

CONCLUSIONS

In the present study, we report that MYH14 can manifest as nonsyndromic prelingual severe sensorineural hearing loss in an AD fashion in Koreans. The results of the present study suggest that further genetic studies of similar patients should consider MYH14 as a causative gene, and cochlear implantation during infant or early childhood should be indicated for those patients with certain MYH14 pathogenic variants.

An investigation of obesity susceptibility genes in Northern Han Chinese by targeted resequencing

Our earlier genome-wide linkage study of body mass index (BMI) showed strong signals from 7q36.3 and 8q21.13. This case-control study set to investigate 2 genomic regions which may harbor variants contributed to development of obesity.We employed targeted resequencing technology to detect single nucleotide polymorphisms (SNPs) in 7q36.3 and 8q21.13 from 16 individuals with obesity. These were compared with 504 East Asians in the 1000 Genomes Project as a reference panel. Linkage disequilibrium (LD) block analysis was performed for the significant SNPs located near the same gene. Genes involved in statistically significant loci were then subject to gene set enrichment analysis (GSEA).The 16 individuals aged between 30 and 60 years with BMI = 33.25 ± 2.22 kg/m. A total of 12,131 genetic variants across all of samples were found. After correcting for multiple testing, 65 SNPs from 25 nearest genes (INSIG1, FABP5, PTPRN2, VIPR2, WDR60, SHH, UBE3C, LMBR1, PAG1, IMPA1, CHMP4, SNX16, BLACE, EN2, CNPY1, LOC100506302, RBM33, LOC389602, LOC285889, LINC01006, NOM1, DNAJB6, LOC101927914, ESYT2, LINC00689) were associated with obesity at significant level q-value ≤ 0.05. LD block analysis showed there were 10 pairs of loci with D' ≥ 0.8 and r ≥ 0.8. GSEA further identified 2 major related gene sets, involving lipid raft and lipid metabolic process, with FDR values <0.12 and <0.4, respectively.Our data are the first documentation of genetic variants in 7q36.3 and 8q21.13 associated with obesity using target capture sequencing and Northern Han Chinese samples. Additional replication and functional studies are merited to validate our findings.

Targeted Resequencing of Deafness Genes Reveals a Founder MYO15A Variant in Northeastern Brazil

Identifying the genetic etiology in a person with hearing loss (HL) is challenging due to the extreme genetic heterogeneity in HL and the population-specific variability. In this study, after excluding GJB2 variants, targeted resequencing of 180 deafness-related genes revealed the causative variants in 11 of 19 (58%) Brazilian probands with autosomal recessive HL. Identified pathogenic variants were in MYO15A (10 families) and CLDN14 (one family). Remarkably, the MYO15A p.(Val1400Met) variant was identified in eight families from the city of Monte Santo in the northeast region of Brazil. Haplotype analysis of this variant was consistent with a single founder. No other cases with this variant were detected among 105 simplex cases from other cities of northeastern Brazil, suggesting that this variant is confined to a geographical region. This study suggests that it is feasible to develop population-specific screening for deafness variants once causative variants are identified in different geographical groups.

Targeted re-sequencing confirms the importance of chemosensory genes in aphid host race differentiation

Host‐associated races of phytophagous insects provide a model for understanding how adaptation to a new environment can lead to reproductive isolation and speciation, ultimately enabling us to connect barriers to gene flow to adaptive causes of divergence. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on legume species and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. As host choice produces assortative mating, understanding the underlying mechanisms of choice will contribute directly to understanding of speciation. As host choice in the pea aphid is likely mediated by smell and taste, we use capture sequencing and SNP genotyping to test for the role of chemosensory genes in the divergence between eight host plant species across the continuum of differentiation and sampled at multiple locations across western Europe. We show high differentiation of chemosensory loci relative to control loci in a broad set of pea aphid races and localities, using a model‐free approach based on principal component analysis. Olfactory and gustatory receptors form the majority of highly differentiated genes and include loci that were already identified as outliers in a previous study focusing on the three most closely related host races. Consistent indications that chemosensory genes may be good candidates for local adaptation and barriers to gene flow in the pea aphid open the way to further investigations aiming to understand their impact on gene flow and to determine their precise functions in response to host plant metabolites.

Mutations in the IRBIT domain of ITPR1 are a frequent cause of autosomal dominant nonprogressive congenital ataxia

Congenital ataxias are nonprogressive neurological disorders characterized by neonatal hypotonia, developmental delay and ataxia, variably associated with intellectual disability and other neurological or extraneurological features. We performed trio-based whole-exome sequencing of 12 families with congenital cerebellar and/or vermis atrophy in parallel with targeted next-generation sequencing of known ataxia genes (CACNA1A, ITPR1, KCNC3, ATP2B3 and GRM1) in 12 additional patients with a similar phenotype. Novel pathological mutations of ITPR1 (inositol 1,4,5-trisphosphate receptor, type 1) were found in seven patients from four families (4/24, ∼16.8%) all localized in the IRBIT (inositol triphosphate receptor binding protein) domain which plays an essential role in the regulation of neuronal plasticity and development. Our study expands the mutational spectrum of ITPR1-related congenital ataxia and indicates that ITPR1 gene screening should be implemented in this subgroup of ataxias.

DNA Diagnostics of Hereditary Hearing Loss: A Targeted Resequencing Approach Combined with a Mutation Classification System

Although there are nearly 100 different causative genes identified for nonsyndromic hearing loss (NSHL), Sanger sequencing-based DNA diagnostics usually only analyses three, namely, GJB2, SLC26A4, and OTOF. As this is seen as inadequate, there is a need for high-throughput diagnostic methods to detect disease-causing variations, including single-nucleotide variations (SNVs), insertions/deletions (Indels), and copy-number variations (CNVs). In this study, a targeted resequencing panel for hearing loss was developed including 79 genes for NSHL and selected forms of syndromic hearing loss. One-hundred thirty one presumed autosomal-recessive NSHL (arNSHL) patients of Western-European ethnicity were analyzed for SNVs, Indels, and CNVs. In addition, we established a straightforward variant classification system to deal with the large number of variants encountered. We estimate that combining prescreening of GJB2 with our panel leads to a diagnosis in 25%-30% of patients. Our data show that after GJB2, the most commonly mutated genes in a Western-European population are TMC1, MYO15A, and MYO7A (3.1%). CNV analysis resulted in the identification of causative variants in two patients in OTOA and STRC. One of the major challenges for diagnostic gene panels is assigning pathogenicity for variants. A collaborative database collecting all identified variants from multiple centers could be a valuable resource for hearing loss diagnostics.

Reduced levels of protein recoding by A-to-I RNA editing in Alzheimer's disease

Adenosine to inosine (A-to-I) RNA editing, catalyzed by the ADAR enzyme family, acts on dsRNA structures within pre-mRNA molecules. Editing of the coding part of the mRNA may lead to recoding, amino acid substitution in the resulting protein, possibly modifying its biochemical and biophysical properties. Altered RNA editing patterns have been observed in various neurological pathologies. Here, we present a comprehensive study of recoding by RNA editing in Alzheimer's disease (AD), the most common cause of irreversible dementia. We have used a targeted resequencing approach supplemented by a microfluidic-based high-throughput PCR coupled with next-generation sequencing to accurately quantify A-to-I RNA editing levels in a preselected set of target sites, mostly located within the coding sequence of synaptic genes. Overall, editing levels decreased in AD patients' brain tissues, mainly in the hippocampus and to a lesser degree in the temporal and frontal lobes. Differential RNA editing levels were observed in 35 target sites within 22 genes. These results may shed light on a possible association between the neurodegenerative processes typical for AD and deficient RNA editing.

Flexible, scalable, and efficient targeted resequencing on a benchtop sequencer for variant detection in clinical practice

The release of benchtop next-generation sequencing (NGS) instruments has paved the way to implement the technology in clinical setting. The need for flexible, qualitative, and cost-efficient workflows is high. We used singleplex-PCR for highly efficient target enrichment, allowing us to reach the quality standards set in Sanger sequencing-based diagnostics. For the library preparation, a modified NexteraXT protocol was used, followed by sequencing on a MiSeq instrument. With an innovative pooling strategy, high flexibility, scalability, and cost-efficiency were obtained, independent of the availability of commercial kits. The approach was validated for ∼250 genes associated with monogenic disorders. An overall sensitivity (>99%) similar to Sanger sequencing was observed in combination with a positive predictive value of >98%. The distribution of coverage was highly uniform, guaranteeing a minimal number of gaps to be filled with alternative methods. ISO15189-accreditation was obtained for the workflow. A major asset of the singleplex PCR-based enrichment is that new targets can be easily implemented. Diagnostic laboratories have validated assays available ensuring that the proposed workflow can easily be adopted. Although our platform was optimized for constitutional variant detection of monogenic disease genes, it is now also used as a model for somatic mutation detection in acquired diseases.

The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

Dynamics of copy number variation in host races of the pea aphid

Copy number variation (CNV) makes a major contribution to overall genetic variation and is suspected to play an important role in adaptation. However, aside from a few model species, the extent of CNV in natural populations has seldom been investigated. Here, we report on CNV in the pea aphid Acyrthosiphon pisum, a powerful system for studying the genetic architecture of host-plant adaptation and speciation thanks to multiple host races forming a continuum of genetic divergence. Recent studies have highlighted the potential importance of chemosensory genes, including the gustatory and olfactory receptor gene families (Gr and Or, respectively), in the process of host race formation. We used targeted resequencing to achieve a very high depth of coverage, and thereby revealed the extent of CNV of 434 genes, including 150 chemosensory genes, in 104 individuals distributed across eight host races of the pea aphid. We found that CNV was widespread in our global sample, with a significantly higher occurrence in multigene families, especially in Ors. We also observed a decrease in the gene probability of being completely duplicated or deleted (CDD) with increase in coding sequence length. Genes with CDD variants were usually more polymorphic for copy number, especially in the P450 gene family where toxin resistance may be related to gene dosage. We found that Gr were overrepresented among genes discriminating host races, as were CDD genes and pseudogenes. Our observations shed new light on CNV dynamics and are consistent with CNV playing a role in both local adaptation and speciation.

Targeted next-generation sequencing on Hirschsprung disease: a pilot study exploits DNA pooling

To adopt an efficient approach of identifying rare variants possibly related to Hirschsprung disease (HSCR), a pilot study was set up to evaluate the performance of a newly designed protocol for next generation targeted resquencing. In total, 20 Chinese HSCR patients and 20 Chinese sex-matched individuals with no HSCR were included, for which coding sequences (CDS) of 62 genes known to be in signaling pathways relevant to enteric nervous system development were selected for capture and sequencing. Blood DNAs from eight pools of five cases or controls were enriched by PCR-based RainDance technology (RDT) and then sequenced on a 454 FLX platform. As technical validation, five patients from case Pool-3 were also independently enriched by RDT, indexed with barcode and sequenced with sufficient coverage. Assessment for CDS single nucleotide variants showed DNA pooling performed well (specificity/sensitivity at 98.4%/83.7%) at the common variant level; but relatively worse (specificity/sensitivity at 65.5%/61.3%) at the rare variant level. Further Sanger sequencing only validated five out of 12 rare damaging variants likely involved in HSCR. Hence more improvement at variant detection and sequencing technology is needed to realize the potential of DNA pooling for large-scale resequencing projects.

Assessment of a targeted resequencing assay as a support tool in the diagnosis of lysosomal storage disorders

BACKGROUND

With over 50 different disorders and a combined incidence of up to 1/3000 births, lysosomal storage diseases (LSDs) constitute a major public health problem and place an enormous burden on affected individuals and their families. Many factors make LSD diagnosis difficult, including phenotype and penetrance variability, shared signs and symptoms, and problems inherent to biochemical diagnosis. Developing a powerful diagnostic tool could mitigate the protracted diagnostic process for these families, lead to better outcomes for current and proposed therapies, and provide the basis for more appropriate genetic counseling.

METHODS

We have designed a targeted resequencing assay for the simultaneous testing of 57 lysosomal genes, using in-solution capture as the enrichment method and two different sequencing platforms. A total of 84 patients with high to moderate-or low suspicion index for LSD were enrolled in different centers in Spain and Portugal, including 18 positive controls.

RESULTS

We correctly diagnosed 18 positive blinded controls, provided genetic diagnosis to 25 potential LSD patients, and ended with 18 diagnostic odysseys.

CONCLUSION

We report the assessment of a next-generation-sequencing-based approach as an accessory tool in the diagnosis of LSDs, a group of disorders which have overlapping clinical profiles and genetic heterogeneity. We have also identified and quantified the strengths and limitations of next generation sequencing (NGS) technology applied to diagnosis.

Accurate molecular diagnosis of phenylketonuria and tetrahydrobiopterin-deficient hyperphenylalaninemias using high-throughput targeted sequencing

Genetic diagnostics of phenylketonuria (PKU) and tetrahydrobiopterin (BH4) deficient hyperphenylalaninemia (BH4DH) rely on methods that scan for known mutations or on laborious molecular tools that use Sanger sequencing. We have implemented a novel and much more efficient strategy based on high-throughput multiplex-targeted resequencing of four genes (PAH, GCH1, PTS, and QDPR) that, when affected by loss-of-function mutations, cause PKU and BH4DH. We have validated this approach in a cohort of 95 samples with the previously known PAH, GCH1, PTS, and QDPR mutations and one control sample. Pooled barcoded DNA libraries were enriched using a custom NimbleGen SeqCap EZ Choice array and sequenced using a HiSeq2000 sequencer. The combination of several robust bioinformatics tools allowed us to detect all known pathogenic mutations (point mutations, short insertions/deletions, and large genomic rearrangements) in the 95 samples, without detecting spurious calls in these genes in the control sample. We then used the same capture assay in a discovery cohort of 11 uncharacterized HPA patients using a MiSeq sequencer. In addition, we report the precise characterization of the breakpoints of four genomic rearrangements in PAH, including a novel deletion of 899 bp in intron 3. Our study is a proof-of-principle that high-throughput-targeted resequencing is ready to substitute classical molecular methods to perform differential genetic diagnosis of hyperphenylalaninemias, allowing the establishment of specifically tailored treatments a few days after birth.

Targeted deep resequencing identifies MID2 mutation for X-linked intellectual disability with varied disease severity in a large kindred from India

We report a novel missense mutation (c.1040G>A, p.Arg347Gln) in MID2, which encodes ubiquitin ligase E3, as the likely cause of X-linked mental retardation in a large kindred. The mutation was observed in all affected and obligate carriers but not in any unaffected males of the family or in population controls (n = 200). When transiently expressed in HEK293T cell line, the mutation was found to abolish the function of the COS domain in the protein. The GFP-tagged mutant protein accumulated in the cytoplasm instead of binding to the cytoskeleton resulting in its altered subcellular localization. Screening of coding exons of this gene in additional 480 unrelated individuals with idiopathic intellectual disability identified another novel variation p.Asn343Ser. This study highlights the growing role of the ubiquitin pathway in intellectual disability and also, the difference in MID2 determined phenotype observed in this study compared with that of its paralogue MID1 reported in literature.

Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond

Advanced resources for genome-assisted research in barley (Hordeum vulgare) including a whole-genome shotgun assembly and an integrated physical map have recently become available. These have made possible studies that aim to assess genetic diversity or to isolate single genes by whole-genome resequencing and in silico variant detection. However such an approach remains expensive given the 5 Gb size of the barley genome. Targeted sequencing of the mRNA-coding exome reduces barley genomic complexity more than 50-fold, thus dramatically reducing this heavy sequencing and analysis load. We have developed and employed an in-solution hybridization-based sequence capture platform to selectively enrich for a 61.6 megabase coding sequence target that includes predicted genes from the genome assembly of the cultivar Morex as well as publicly available full-length cDNAs and de novo assembled RNA-Seq consensus sequence contigs. The platform provides a highly specific capture with substantial and reproducible enrichment of targeted exons, both for cultivated barley and related species. We show that this exome capture platform provides a clear path towards a broader and deeper understanding of the natural variation residing in the mRNA-coding part of the barley genome and will thus constitute a valuable resource for applications such as mapping-by-sequencing and genetic diversity analyzes.

Using Ascertainment for Targeted Resequencing to Increase Power to Identify Causal Variants

Researchers continue to use genome-wide association studies (GWAS) to find the genetic markers associated with disease. Recent studies have added to the typical two-stage analysis a third stage that uses targeted resequencing on a randomly selected subset of the cases to detect the causal single-nucleotide polymorphism (SNP). We propose a design for targeted resequencing that increases the power to detect the causal variant. The design features an ascertainment scheme wherein only those cases with the presence of a risk allele are selected for targeted resequencing. We simulated a disease with a single causal SNP to evaluate our method versus a targeted resequencing design using randomly selected individuals. The simulation studies showed that ascertaining individuals for the targeted resequencing can substantially increase the power to detect a causal SNP, without increasing the false-positive rate.