Identification of intermediate-sized deletions and inference of their impact on gene expression in a human population

Mapping crossover events of mouse meiotic recombination by restriction fragment ligation-based Refresh-seq

Single-cell whole-genome sequencing methods have undergone great improvements over the past decade. However, allele dropout, which means the inability to detect both alleles simultaneously in an individual diploid cell, largely restricts the application of these methods particularly for medical applications. Here, we develop a new single-cell whole-genome sequencing method based on third-generation sequencing (TGS) platform named Refresh-seq (restriction fragment ligation-based genome amplification and TGS). It is based on restriction endonuclease cutting and ligation strategy in which two alleles in an individual cell can be cut into equal fragments and tend to be amplified simultaneously. As a new single-cell long-read genome sequencing method, Refresh-seq features much lower allele dropout rate compared with SMOOTH-seq. Furthermore, we apply Refresh-seq to 688 sperm cells and 272 female haploid cells (secondary polar bodies and parthenogenetic oocytes) from F1 hybrid mice. We acquire high-resolution genetic map of mouse meiosis recombination at low sequencing depth and reveal the sexual dimorphism in meiotic crossovers. We also phase the structure variations (deletions and insertions) in sperm cells and female haploid cells with high precision. Refresh-seq shows great performance in screening aneuploid sperm cells and oocytes due to the low allele dropout rate and has great potential for medical applications such as preimplantation genetic diagnosis.

Long-read-based single sperm genome sequencing for chromosome-wide haplotype phasing of both SNPs and SVs

Although localized haploid phasing can be achieved using long read genome sequencing without parental data, reliable chromosome-scale phasing remains a great challenge. Given that sperm is a natural haploid cell, single-sperm genome sequencing can provide a chromosome-wide phase signal. Due to the limitation of read length, current short-read-based single-sperm genome sequencing methods can only achieve SNP haplotyping and come with difficulties in detecting and haplotyping structural variations (SVs) in complex genomic regions. To overcome these limitations, we developed a long-read-based single-sperm genome sequencing method and a corresponding data analysis pipeline that can accurately identify crossover events and chromosomal level aneuploidies in single sperm and efficiently detect SVs within individual sperm cells. Importantly, without parental genome information, our method can accurately conduct de novo phasing of heterozygous SVs as well as SNPs from male individuals at the whole chromosome scale. The accuracy for phasing of SVs was as high as 98.59% using 100 single sperm cells, and the accuracy for phasing of SNPs was as high as 99.95%. Additionally, our method reliably enabled deduction of the repeat expansions of haplotype-resolved STRs/VNTRs in single sperm cells. Our method provides a new opportunity for studying haplotype-related genetics in mammals.

Targeted deletion of ecto-5'-nucleotidase results in retention of inosine monophosphate content in postmortem muscle of medaka (Oryzias latipes)

Inosine monophosphate (IMP) is an important indicator of meat freshness and contributes to its umami taste. An attractive strategy for enhancing umami is to suppress the IMP-degrading activity and increase the IMP content in the skeletal muscle through genome editing technology using the CRISPR-Cas9 system. However, the molecular mechanisms underlying IMP degradation remain unclear. We cloned two ecto-5'-nucleotidase genes, designated as ecto-5'-nucleotidase-a (nt5ea) and ecto-5'-nucleotidase-b (nt5eb), from medaka (Oryzias latipes), a vertebrate model organism. Expression analysis using embryos showed that nt5ea or nt5eb overexpression remarkably upregulated IMP degradation, and that the IMP-degrading activity was higher in Nt5ea than in Nt5eb. Furthermore, we established frame-shifted or large deletion (lacking nt5ea or nt5eb locus) mutant strains and assayed the effects of gene disruptions on the amount of IMP in skeletal muscle. The nt5ea-deficient medaka showed considerable higher levels of IMP at 48 h postmortem than did the wild-type fish. The nt5eb mutants also exhibited higher IMP contents than that in the wild types, but the increase was less than that in the nt5ea mutants. Our results demonstrated that nt5e is an important regulator of IMP levels in skeletal muscle and that its loss of function was effective in maintaining IMP content.

VAMP2 Expression and Genotype Are Possible Discriminators in Different Forms of Dementia

Vascular alterations often overlap with neurodegeneration, resulting in mixed forms of dementia (MD) that are hard to differentiate from Alzheimer’s Disease (AD). The 26 bp intergenic polymorphism of VAMP2, a key component of SNARE complex, as well as its mRNA and protein levels are associated with neurological diseases. We evaluated ApoE4 and VAMP2 26 bp Ins/Del genotype distribution in 177 AD, 132 MD, 115 Mild Cognitive Impairment (MCI) and 250 individuals without cognitive decline (CT), as well as VAMP2 gene expression in a subset of 73 AD, 122 MD, 103 MCI and 140 CT. Forty-two MCI evolved to AD (22 MCI-AD) or MD (20 MCI-MD) over time. VAMP2 mRNA was higher in MD compared to AD (p = 0.0013) and CT (p = 0.0017), and in MCI-MD compared to MCI-AD (p < 0.001) after correcting for age, gender, MMSE and ApoE4 +/− covariates (p_c = 0.004). A higher VAMP2 expression was observed in subjects carrying the minor allele Del compared to those carrying the Ins/Ins genotype (p = 0.012). Finally, Del/Del genotype was more frequently carried by MD/MCI-MD compared to CT (p_c = 0.036). These results suggest that VAMP2 mRNA expression can discriminate mixed form of dementia from AD, possibly being a biomarker of AD evolution in MCI patients.

Compound genetic etiology in a patient with a syndrome including diabetes, intellectual deficiency and distichiasis

BACKGROUND

We studied a young woman with atypical diabetes associated with mild intellectual disability, lymphedema distichiasis syndrome (LDS) and polymalformative syndrome including distichiasis. We used different genetic tools to identify causative pathogenic mutations and/or copy number variations.

RESULTS

Although proband's, diabetes mellitus occurred during childhood, type 1 diabetes was unlikely due to the absence of detectable autoimmunity. DNA microarray analysis first identified a de novo, heterozygous deletion at the chr16q24.2 locus. Previously, thirty-three pathogenic or likely pathogenic deletions encompassing this locus have been reported in patients presenting with intellectual deficiency, obesity and/or lymphedema but not with diabetes. Of note, the deletion encompassed two topological association domains, whose one included FOXC2 that is known to be linked with LDS. Via whole-exome sequencing, we found a heterozygous, likely pathogenic variant in WFS1 (encoding wolframin endoplasmic reticulum [ER] transmembrane glycoprotein) which was inherited from her father who also had diabetes. WFS1 is known to be involved in monogenic diabetes. We also found a likely pathogenic variant in USP9X (encoding ubiquitin specific peptidase 9 X-linked) that is involved in X-linked intellectual disability, which was inherited from her mother who had dyscalculia and dyspraxia.

CONCLUSIONS

Our comprehensive genetic analysis suggested that the peculiar phenotypes of our patient were possibly due to the combination of multiple genetic causes including chr16q24.2 deletion, and two likely pathogenic variants in WFS1 and USP9X.

Three Novel Structural Variations at MHC and IL12B Predisposing to Psoriasis

BACKGROUND

Structural variations (SVs, defined as DNA variants ≥50 bp) have been associated with various complex human diseases. However, research to screen the whole genome for SVs predisposing to psoriasis is still lacking.

OBJECTIVES

This study aimed to investigate the association of SVs and psoriasis.

METHODS

We performed a genome-wide screen on SVs using an imputation method on 5 independent cohorts with 45,386 subjects from the Chinese Han population. Fine mapping analysis, genetic interaction analysis and RNA expression analysis were conducted to explore the mechanism of SVs.

RESULTS

We obtained 4,535 SVs in total and identified 2 novel deletions (esv3608550, OR=2.73, P<2.00×10^-308 ; esv3608542, OR=0.47, P=7.40×10^-28 ) at 6q21.33 (MHC), 1 novel Alu element insertion (esv3607339, OR=1.22, P=1.18×10^-35 ) at 5q33.3 (IL12B), and confirmed 1 previously reported deletion (esv3587563, OR=1.30, P=9.52×10^-60 ) at 1q21.2 (LCE) for psoriasis. Fine mapping analysis including SNPs and small Insertions/Deletions (InDels) revealed that esv3608550 and esv3608542 were independently associated with psoriasis, and a novel independent SNP (rs9378188, OR=1.65, P=3.46×10^-38 ) was identified at 6q21.33. By genetic interaction analysis and RNA expression analysis, we speculate that the association of 2 deletions at 6q21.33 with psoriasis might relate to their influence on the expression of HLA-C.

CONCLUSIONS

Our study constructed the most comprehensive SV map for psoriasis thus far and enriched the genetic architecture and pathogenesis of psoriasis as well as highlighted the nonnegligible impact of SVs on complex diseases.

Characterization of intermediate-sized insertions using whole-genome sequencing data and analysis of their functional impact on gene expression

Intermediate-sized insertions are one of the structural variants contributing to genome diversity. However, due to technical difficulties in identifying them, their importance in disease pathogenicity and gene expression regulation remains unclear. We used whole-genome sequencing data of 174 Japanese samples to characterize intermediate-sized insertions using a highly-accurate insertion calling method (IMSindel software and joint-call recovery) and obtained a catalogue of 4254 insertions. We constructed an imputation panel comprising of insertions and SNVs from all samples, and conducted imputation of intermediate-sized insertions for 82 publicly-available Japanese samples. Positive Predictive Value of imputation, evaluated using Nanopore long-read sequencing data, was 97%. Subsequent eQTL analysis predicted 128 (~ 3.0%) insertions as causative for gene expression level changes. Enrichment analysis of causal insertions for genome regulatory elements showed significant associations with CTCF-binding sites, super-enhancers, and promoters. Among 17 causal insertions found in the same causal set with GWAS hits, there were insertions associated with changes in expression of cancer-related genes such as BRCA1, ZNF222, and ABCB10. Analysis of insertions sequences revealed that 461 insertions were short tandem duplications frequently found in early-replicating regions of genome. Furthermore, comparison of functional importance of intermediate-sized insertions with that of intermediate-sized deletions detected in the same sample set in our previous study showed that insertions were more frequent in genic regions, and proportion of functional candidates was smaller in insertions. Here, we characterize a high-confidence set of intermediate-sized insertions and indicate their importance in gene expression regulation. Our results emphasize the importance of considering intermediate-sized insertions in trait association studies.

Whole-genome sequencing with long reads reveals complex structure and origin of structural variation in human genetic variations and somatic mutations in cancer

BACKGROUND

Identification of germline variation and somatic mutations is a major issue in human genetics. However, due to the limitations of DNA sequencing technologies and computational algorithms, our understanding of genetic variation and somatic mutations is far from complete.

METHODS

In the present study, we performed whole-genome sequencing using long-read sequencing technology (Oxford Nanopore) for 11 Japanese liver cancers and matched normal samples which were previously sequenced for the International Cancer Genome Consortium (ICGC). We constructed an analysis pipeline for the long-read data and identified germline and somatic structural variations (SVs).

RESULTS

In polymorphic germline SVs, our analysis identified 8004 insertions, 6389 deletions, 27 inversions, and 32 intra-chromosomal translocations. By comparing to the chimpanzee genome, we correctly inferred events that caused insertions and deletions and found that most insertions were caused by transposons and Alu is the most predominant source, while other types of insertions, such as tandem duplications and processed pseudogenes, are rare. We inferred mechanisms of deletion generations and found that most non-allelic homolog recombination (NAHR) events were caused by recombination errors in SINEs. Analysis of somatic mutations in liver cancers showed that long reads could detect larger numbers of SVs than a previous short-read study and that mechanisms of cancer SV generation were different from that of germline deletions.

CONCLUSIONS

Our analysis provides a comprehensive catalog of polymorphic and somatic SVs, as well as their possible causes. Our software are available at https://github.com/afujimoto/CAMPHOR and https://github.com/afujimoto/CAMPHORsomatic .