Early detection of active Human CytomegaloVirus (hCMV) infection in pregnant women using data generated for noninvasive fetal aneuploidy testing.

BACKGROUND: Prenatal hCMV infections can lead to severe embryopathy and neurological sequelae in neonates. Screening during pregnancy is not recommended by global societies, as there is no effective therapy. Recently, several groups showed that maternal-fetal hCMV transmission can be strongly reduced by administering anti-viral agents early in pregnancy. This calls for a screening method to identify at risk pregnancies at an appropriate gestational age, with the possibility for large-scale enrolment. Non-Invasive Prenatal Testing (NIPT) for fetal aneuploidy screening early in pregnancy is already implemented in many countries and performed on a large-scale basis. We investigated the use of whole genome cell-free DNA (cfDNA) sequencing data, generated for the purpose of NIPT, as (pre-)screening tool to identify women with active hCMV-infections, eligible for therapy. METHODS: Coded raw sequencing NIPT data from 204,818 pregnant women from three testing laboratories were analyzed for the presence of hCMV-cfDNA. Samples were stratified by cfDNA-hCMV load. For validation and interpretation, diagnostic hCMV-qPCR and serology testing were performed on a subset of cfDNA-hCMV-positive (n = 112) and -negative (n = 127) samples. FINDINGS: In 1930 samples (0.94%) hCMV fragments were detected. Validation by hCMV-qPCR showed that samples with high cfDNA-hCMV load tested positive and cfDNA-hCMV-negative samples tested negative. In 32/112 cfDNA-hCMV-positive samples (28.6%) the serological profile suggested a recent primary infection: this was more likely in samples with high cfDNA-hCMV load (78.6%) than in samples with low cfDNA-hCMV load (11.0%). In none of the cfDNA-hCMV-negative samples serology was indicative of a recent primary infection. INTERPRETATION: Our study shows that large-scale (pre-)screening for both genetic fetal aberrations and active maternal hCMV infections during pregnancy can be combined in one cfDNA sequencing test, performed on a single blood sample, drawn in the first trimester of pregnancy. FUNDING: This work was partly funded by the Prenatal Screening Foundation Nijmegen, the Netherlands.

The cell-free DNA virome of 108,349 Dutch pregnant women.

OBJECTIVE: Viral infections during pregnancy are a major health concern to mother and fetus. By repurposing cell-free Non Invasive Prenatal Testing (NIPT) sequencing data, we investigated prevalence and abundance of viral DNA in a cohort of 108,349 pregnant women. METHOD: Cell-free DNA (cfDNA) sequencing reads that did not map to any of the human chromosomes or mitochondrial DNA of the human reference genome build GRCh38 were aligned to 224 DNA viruses selected from the NCBI refseq viral database. RESULTS: In total 443,665 reads of viral origin were detected across 42,273 samples representing 165 viral species. Several are known to be potentially harmful during pregnancy and/or childbirth, including Cytomegalovirus, Parvovirus B19 and Hepatitis B. Viral sequences were mostly detected at very low abundance. However, several cases had exceptionally high viral loads for Parvovirus B19, Hepatitis B and others. We found statistically significant associations between presence of viral DNA and gestational age, maternal age, fetal fraction, cfDNA concentration and others. CONCLUSION: We demonstrate the feasibility to detect viral DNA from typical genome-wide NIPT cfDNA sequencing and describe the main characteristics of the viral DNA in our cohort. Our dataset of detected viral sequence reads is made publicly available to guide future clinical implementations.

Clinically relevant DNA viruses in pregnancy.

Infections by DNA viruses during pregnancy are associated with increased health risks to both mother and fetus. Although not all DNA viruses are related to an increased risk of complications during pregnancy, several can directly infect the fetus and/or cause placental dysfunction. During Non-Invasive Prenatal Testing analysis, the presence of viral DNA can be detected, theoretically allowing screening early in pregnancy. Although treatment options are currently limited, this might rapidly change in the near future. It is therefore important to be aware of the potential impact of these viruses on feto-maternal health. In this manuscript we provide a brief introduction into the most commonly detected DNA viruses in human cell-free DNA sequencing experiments and their pathogenic potential during pregnancy.

Variability in Fetal Fraction Estimation: Comparing Fetal Fractions Reported by Noninvasive Prenatal Testing Providers Globally.

BACKGROUND: Fetal fraction (FF) measurement is considered important for reliable noninvasive prenatal testing (NIPT). Using minimal FF threshold as a quality parameter is under debate. We evaluated the variability in reported FFs of individual samples between providers and laboratories and within a single laboratory. METHODS: Genomic quality assessment and European Molecular Genetics Quality Network provide joint proficiency testing for NIPT. We compared reported FFs across all laboratories and stratified according to test methodologies. A single sample was sequenced repeatedly and FF estimated by 2 bioinformatics methods: Veriseq2 and SeqFF. Finally, we compared FFs by Veriseq and SeqFF in 87 351 NIPT samples. RESULTS: For each proficiency test sample we observed a large variability in reported FF, SDs and CVs ranging from 1.7 to 3.6 and 17.0 to 35.8, respectively. FF measurements reported by single nucleotide polymorphism-based methods had smaller SDs (0.5 to 2.4) compared to whole genome sequencing-based methods (1.8 to 2.9). In the internal quality assessment, SDs were similar between SeqFF (SD 1.0) and Veriseq v2 (SD 0.9), but mean FF by Veriseq v2 was higher compared to SeqFF (9.0 vs 6.4, P 0.001). In patient samples, reported FFs were on average 1.12-points higher in Veriseq than in SeqFF (P 0.001). CONCLUSIONS: Current methods do not allow for a reliable and consistent FF estimation. Our data show estimated FF should be regarded as a laboratory-specific range, rather than a precise number. Applying strict universal minimum thresholds might result in unnecessary test failures and should be used with caution.

Extreme enrichment of VNTR-associated polymorphicity in human subtelomeres: genes with most VNTRs are predominantly expressed in the brain.

The human genome harbors numerous structural variants (SVs) which, due to their repetitive nature, are currently underexplored in short-read whole-genome sequencing approaches. Using single-molecule, real-time (SMRT) long-read sequencing technology in combination with FALCON-Unzip, we generated a de novo assembly of the diploid genome of a 115-year-old Dutch cognitively healthy woman. We combined this assembly with two previously published haploid assemblies (CHM1 and CHM13) and the GRCh38 reference genome to create a compendium of SVs that occur across five independent human haplotypes using the graph-based multi-genome aligner REVEAL. Across these five haplotypes, we detected 31,680 euchromatic SVs (>50 bp). Of these, ~62% were comprised of repetitive sequences with 'variable number tandem repeats' (VNTRs), ~10% were mobile elements (Alu, L1, and SVA), while the remaining variants were inversions and indels. We observed that VNTRs with GC-content >60% and repeat patterns longer than 15 bp were 21-fold enriched in the subtelomeric regions (within 5 Mb of the ends of chromosome arms). VNTR lengths can expand to exceed a critical length which is associated with impaired gene transcription. The genes that contained most VNTRs, of which PTPRN2 and DLGAP2 are the most prominent examples, were found to be predominantly expressed in the brain and associated with a wide variety of neurological disorders. Repeat-induced variation represents a sizeable fraction of the genetic variation in human genomes and should be included in investigations of genetic factors associated with phenotypic traits, specifically those associated with neurological disorders. We make available the long and short-read sequence data of the supercentenarian genome, and a compendium of SVs as identified across 5 human haplotypes.

CHOP: haplotype-aware path indexing in population graphs.

The practical use of graph-based reference genomes depends on the ability to align reads to them. Performing substring queries to paths through these graphs lies at the core of this task. The combination of increasing pattern length and encoded variations inevitably leads to a combinatorial explosion of the search space. Instead of heuristic filtering or pruning steps to reduce the complexity, we propose CHOP, a method that constrains the search space by exploiting haplotype information, bounding the search space to the number of haplotypes so that a combinatorial explosion is prevented. We show that CHOP can be applied to large and complex datasets, by applying it on a graph-based representation of the human genome encoding all 80 million variants reported by the 1000 Genomes Project.

Genomic and molecular characterization of preterm birth.

Preterm birth (PTB) complications are the leading cause of long-term morbidity and mortality in children. By using whole blood samples, we integrated whole-genome sequencing (WGS), RNA sequencing (RNA-seq), and DNA methylation data for 270 PTB and 521 control families. We analyzed this combined dataset to identify genomic variants associated with PTB and secondary analyses to identify variants associated with very early PTB (VEPTB) as well as other subcategories of disease that may contribute to PTB. We identified differentially expressed genes (DEGs) and methylated genomic loci and performed expression and methylation quantitative trait loci analyses to link genomic variants to these expression and methylation changes. We performed enrichment tests to identify overlaps between new and known PTB candidate gene systems. We identified 160 significant genomic variants associated with PTB-related phenotypes. The most significant variants, DEGs, and differentially methylated loci were associated with VEPTB. Integration of all data types identified a set of 72 candidate biomarker genes for VEPTB, encompassing genes and those previously associated with PTB. Notably, PTB-associated genes RAB31 and RBPJ were identified by all three data types (WGS, RNA-seq, and methylation). Pathways associated with VEPTB include EGFR and prolactin signaling pathways, inflammation- and immunity-related pathways, chemokine signaling, IFN-γ signaling, and Notch1 signaling. Progress in identifying molecular components of a complex disease is aided by integrated analyses of multiple molecular data types and clinical data. With these data, and by stratifying PTB by subphenotype, we have identified associations between VEPTB and the underlying biology.