Genome sequencing in cytogenetics: Comparison of short-read and linked-read approaches for germline structural variant detection and characterization.

BACKGROUND: Structural variants (SVs) include copy number variants (CNVs) and apparently balanced chromosomal rearrangements (ABCRs). Genome sequencing (GS) enables SV detection at base-pair resolution, but the use of short-read sequencing is limited by repetitive sequences, and long-read approaches are not yet validated for diagnosis. Recently, 10X Genomics proposed Chromium, a technology providing linked-reads to reconstruct long DNA fragments and which could represent a good alternative. No study has compared short-read to linked-read technologies to detect SVs in a constitutional diagnostic setting yet. The aim of this work was to determine whether the 10X Genomics technology enables better detection and comprehension of SVs than short-read WGS. METHODS: We included 13 patients carrying various SVs. Whole genome analyses were performed using paired-end HiSeq X sequencing with (linked-read strategy) or without (short-read strategy) Chromium library preparation. Two different bioinformatic pipelines were used: Variants are called using BreakDancer for short-read strategy and LongRanger for long-read strategy. Variant interpretations were first blinded. RESULTS: The short-read strategy allowed diagnosis of known SV in 10/13 patients. After unblinding, the linked-read strategy identified 10/13 SVs, including one (patient 7) missed by the short-read strategy. CONCLUSION: In conclusion, regarding the results of this study, 10X Genomics solution did not improve the detection and characterization of SV.

Familial breast cancer and DNA repair genes: Insights into known and novel susceptibility genes from the GENESIS study, and implications for multigene panel testing.

Pathogenic variants in BRCA1 and BRCA2 only explain the underlying genetic cause of about 10% of hereditary breast and ovarian cancer families. Because of cost-effectiveness, multigene panel testing is often performed even if the clinical utility of testing most of the genes remains questionable. The purpose of our study was to assess the contribution of rare, deleterious-predicted variants in DNA repair genes in familial breast cancer (BC) in a well-characterized and homogeneous population. We analyzed 113 DNA repair genes selected from either an exome sequencing or a candidate gene approach in the GENESIS study, which includes familial BC cases with no BRCA1 or BRCA2 mutation and having a sister with BC (N = 1,207), and general population controls (N = 1,199). Sequencing data were filtered for rare loss-of-function variants (LoF) and likely deleterious missense variants (MV). We confirmed associations between LoF and MV in PALB2, ATM and CHEK2 and BC occurrence. We also identified for the first time associations between FANCI, MAST1, POLH and RTEL1 and BC susceptibility. Unlike other associated genes, carriers of an ATM LoF had a significantly higher risk of developing BC than carriers of an ATM MV (ORLoF = 17.4 vs. ORMV = 1.6; p Het = 0.002). Hence, our approach allowed us to specify BC relative risks associated with deleterious-predicted variants in PALB2, ATM and CHEK2 and to add MAST1, POLH, RTEL1 and FANCI to the list of DNA repair genes possibly involved in BC susceptibility. We also highlight that different types of variants within the same gene can lead to different risk estimates.

Performance comparison of three DNA extraction kits on human whole-exome data from formalin-fixed paraffin-embedded normal and tumor samples.

Next-generation sequencing (NGS) studies are becoming routinely used for the detection of novel and clinically actionable DNA variants at a pangenomic scale. Such analyses are now used in the clinical practice to enable precision medicine. Formalin-fixed paraffin-embedded (FFPE) tissues are still one of the most abundant source of cancer clinical specimen, unfortunately this method of preparation is known to degrade DNA and therefore compromise subsequent analysis. Some studies have reported that variant detection can be performed on FFPE samples sequenced with NGS techniques, but few or none have done an in-depth coverage analysis and compared the influence of different state-of-the-art FFPE DNA extraction kits on the quality of the variant calling. Here, we generated 42 human whole-exome sequencing data sets from fresh-frozen (FF) and FFPE samples. These samples include normal and tumor tissues from two different organs (liver and colon), that we extracted with three different FFPE extraction kits (QIAamp DNA FFPE Tissue kit and GeneRead DNA FFPE kit from Qiagen, Maxwell™ RSC DNA FFPE Kit from Promega). We determined the rate of concordance of called variants between matched FF and FFPE samples on all common variants (representing at least 86% of the total number of variants for SNVs). The concordance rate is very high between all matched FF / FFPE pairs, with equivalent values for the three kits we analyzed. On the other hand, when looking at the difference between the total number of variants in FF and FFPE, we find a significant variation for the three different FFPE DNA extraction kits. Coverage analysis shows that FFPE samples have less good indicators than FF samples, yet the coverage quality remains above accepted thresholds. We detect limited but statistically significant variations in coverage indicator values between the three FFPE extraction kits. Globally, the GeneRead and QIAamp kits have better variant calling and coverage indicators than the Maxwell kit on the samples used in this study, although this kit performs better on some indicators and has advantages in terms of practical usage. Taken together, our results confirm the potential of FFPE samples analysis for clinical genomic studies, but also indicate that the choice of a FFPE DNA extraction kit should be done with careful testing and analysis beforehand in order to maximize the accuracy of the results.

Genome-wide nucleosome specificity and function of chromatin remodellers in ES cells.

ATP-dependent chromatin remodellers allow access to DNA for transcription factors and the general transcription machinery, but whether mammalian chromatin remodellers target specific nucleosomes to regulate transcription is unclear. Here we present genome-wide remodeller-nucleosome interaction profiles for the chromatin remodellers Chd1, Chd2, Chd4, Chd6, Chd8, Chd9, Brg1 and Ep400 in mouse embryonic stem (ES) cells. These remodellers bind one or both full nucleosomes that flank micrococcal nuclease (MNase)-defined nucleosome-free promoter regions (NFRs), where they separate divergent transcription. Surprisingly, large CpG-rich NFRs that extend downstream of annotated transcriptional start sites are nevertheless bound by non-nucleosomal or subnucleosomal histone variants (H3.3 and H2A.Z) and marked by H3K4me3 and H3K27ac modifications. RNA polymerase II therefore navigates hundreds of base pairs of altered chromatin in the sense direction before encountering an MNase-resistant nucleosome at the 3' end of the NFR. Transcriptome analysis after remodeller depletion reveals reciprocal mechanisms of transcriptional regulation by remodellers. Whereas at active genes individual remodellers have either positive or negative roles via altering nucleosome stability, at polycomb-enriched bivalent genes the same remodellers act in an opposite manner. These findings indicate that remodellers target specific nucleosomes at the edge of NFRs, where they regulate ES cell transcriptional programs.

Deep sequencing of the Trypanosoma cruzi GP63 surface proteases reveals diversity and diversifying selection among chronic and congenital Chagas disease patients.

BACKGROUND: Chagas disease results from infection with the diploid protozoan parasite Trypanosoma cruzi. T. cruzi is highly genetically diverse, and multiclonal infections in individual hosts are common, but little studied. In this study, we explore T. cruzi infection multiclonality in the context of age, sex and clinical profile among a cohort of chronic patients, as well as paired congenital cases from Cochabamba, Bolivia and Goias, Brazil using amplicon deep sequencing technology. METHODOLOGY/ PRINCIPAL FINDINGS: A 450bp fragment of the trypomastigote TcGP63I surface protease gene was amplified and sequenced across 70 chronic and 22 congenital cases on the Illumina MiSeq platform. In addition, a second, mitochondrial target--ND5--was sequenced across the same cohort of cases. Several million reads were generated, and sequencing read depths were normalized within patient cohorts (Goias chronic, n = 43, Goias congenital n = 2, Bolivia chronic, n = 27; Bolivia congenital, n = 20), Among chronic cases, analyses of variance indicated no clear correlation between intra-host sequence diversity and age, sex or symptoms, while principal coordinate analyses showed no clustering by symptoms between patients. Between congenital pairs, we found evidence for the transmission of multiple sequence types from mother to infant, as well as widespread instances of novel genotypes in infants. Finally, non-synonymous to synonymous (dn:ds) nucleotide substitution ratios among sequences of TcGP63Ia and TcGP63Ib subfamilies within each cohort provided powerful evidence of strong diversifying selection at this locus. CONCLUSIONS/SIGNIFICANCE: Our results shed light on the diversity of parasite DTUs within each patient, as well as the extent to which parasite strains pass between mother and foetus in congenital cases. Although we were unable to find any evidence that parasite diversity accumulates with age in our study cohorts, putative diversifying selection within members of the TcGP63I gene family suggests a link between genetic diversity within this gene family and survival in the mammalian host.

Differentially expressed genes in autosomal dominant osteopetrosis type II osteoclasts reveal known and novel pathways for osteoclast biology.

Autosomal dominant osteopetrosis type II (ADO II) is a rare, heritable bone disorder characterized by a high bone mass and insufficient osteoclast activity. Mutations in the CLCN7 gene have been reported to cause ADO II. To gain novel insights into the pathways dysregulated in ADOII osteoclasts, we identified changes in gene expression in osteoclasts from patients with a heterozygous mutation of CLCN7. To do this, we carried out a transcriptomic study comparing gene expression in the osteoclasts of patients with ADO II and healthy donors. Our data show that, according to our selection criteria, 182 genes were differentially expressed in osteoclasts from patients and controls. From the 18 displaying the highest change in microarray, we confirmed differential expression for seven by qPCR. Although two of them have previously been found to be expressed in osteoclasts (ITGB5 and SERPINE2), the other five (CES1 (carboxyl esterase 1), UCHL1 (ubiquitin carboxy-terminal esterase L1, also known as ubiquitin thiolesterase), WARS (tryptophanyl-tRNA synthetase), GBP4 (guanylate-binding protein 4), and PRF1) are not yet known to have a role in this cell type. At the protein level, we confirmed elevated expression of ITGB5 and reduced expression of WARS, PRF1, and SERPINE2. Transfection of ClC-7 harboring the G215R mutation into osteoclasts resulted in an increased ITGB5 and reduced PRF1 expression of borderline significance. Finally, we observed that the ADO II patients presented a normal or increased serum level of bone formation markers, demonstrating a coupling between dysfunctional osteoclasts and osteoblasts. Sphingosine kinase 1 mRNA was expressed at the same level in ADO II and control osteoclasts. In conclusion, these data suggest that in addition to an acidification dysfunction caused by the CLCN7 mutation, a change in ITGB5, PRF1, WARS, and SERPINE2 expression could be part of the osteoclastic phenotype of ADO II.

Effects of mechanical wounding on Carica papaya cysteine endopeptidases accumulation and activity.

The mechanical wounding impact on the Carica papaya latex protein pattern was investigated by analyzing three latexes. A first one commercially available, a second harvested from unripe but fully grown fruits, both obtained from regularly tapped fruits. A third one was collected from similar fruits but wounded for the first time. The results demonstrated both quantitative and qualitative changes in the protein content and in the enzymatic activity. Repeated wounding results in either, accumulation or activation (or both of them) of papain, chymopapain and caricain. Furthermore, new cysteine protease activity was found to transiently accumulate in the latex collected from newly wounded fruits. The possible implication of this enzymatic material in the papaya cysteine endopeptidases pro-forms activation is discussed.

Large-scale pre-diagnosis study of fetal RHD genotyping by PCR on plasma DNA from RhD-negative pregnant women.

BACKGROUND: The routine prenatal determination of fetal RhD blood group would be very useful in the management of pregnancies in RhD-negative women, as up to 40% of these pregnancies bear a RhD-negative fetus. The fetal DNA present in maternal plasma offers an opportunity for risk-free prenatal diagnosis. AIM: This study focused on the feasibility and accuracy of large-scale RhD fetal diagnosis in non-immunized and anti-D immunized RhD-negative women. METHODS: Plasma DNA was extracted from 893 RhD-negative pregnant women and amplified in exons 7 and 10 of the RHD gene using conventional and real-time PCR. The results were then compared with the RHD fetal genotype determined on amniotic cells and/or the RhD phenotype of the red blood cells of the infants at birth. RESULTS: After exclusion of 42 samples from women exhibiting a nonfunctional or rearranged RHD gene, fetal RhD status was predicted with a 99.5% accuracy. A strategy is also proposed to avoid the small number of false-positive and -negative results. CONCLUSION: Fetal RHD genotyping from maternal plasma DNA in different clinical situations may be used with confidence.