Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing

Genome-wide somatic mutation analysis of sinonasal adenocarcinoma with and without wood dust exposure

BACKGROUND

Sinonasal adenocarcinoma is a rare cancer, encompassing two different entities, the intestinal-type sinonasal adenocarcinoma (ITAC) and the non-intestinal-type sinonasal adenocarcinoma (non-ITAC). Occurrence of ITAC is strongly associated with exposure to hardwood dusts. In countries with predominant exposure to softwood dust the occurrence of sinonasal adenocarcinomas is lower and the relative amount of non-ITACs to ITACs is higher. The molecular mechanisms behind the tumorigenic effects of wood dust remain largely unknown.

METHODS

We carried out whole-genome sequencing of formalin-fixed paraffin-embedded (FFPE) samples of sinonasal adenocarcinomas from ten wood dust-exposed and six non-exposed individuals, with partial tobacco exposure data. Sequences were analyzed for the presence of mutational signatures matching COSMIC database signatures. Driver mutations and CN variant regions were characterized.

RESULTS

Mutation burden was higher in samples of wood dust-exposed patients (p = 0.016). Reactive oxygen species (ROS) damage-related mutational signatures were almost exclusively identified in ITAC subtype samples (p = 0.00055). Tobacco smoke mutational signatures were observed in samples of patients with tobacco exposure or missing information, but not in samples from non-exposed patients. A tetraploidy copy number (CN) signature was enriched in ITAC subtype (p = 0.042). CN variation included recurrent gains in COSMIC Cancer Gene Census genes TERT, SDHA, RAC1, ETV1, PCM1, and MYC. Pathogenic variants were observed most frequently in TP53, NF1, CHD2, BRAF, APC, and LRP1B. Driver mutations and copy number gains did not segregate by subtype.

CONCLUSIONS

Our analysis identified distinct mutational characteristics in ITAC and non-ITAC. Mutational signature analysis may eventually become useful for documentation of occupation-related cancer, while the exact mechanisms behind wood dust-driven carcinogenesis remain elusive. The presence of homologous recombination deficiency signatures implies a novel opportunity for treatment, but further studies are needed.

Identification and classification of the genomes of novel microviruses in poultry slaughterhouse

Microviridae is a family of phages with circular ssDNA genomes and they are widely found in various environments and organisms. In this study, virome techniques were employed to explore potential members of Microviridae in a poultry slaughterhouse, leading to the identification of 98 novel and complete microvirus genomes. Using a similarity clustering network classification approach, these viruses were found to belong to at least 6 new subfamilies within Microviridae and 3 higher-level taxonomic units. Genome size, GC content and genome structure of these new taxa showed evident regularities, validating the rationality of our classification method. Our method can divide microviruses into about 45 additional detailed clusters, which may serve as a new standard for classifying Microviridae members. Furthermore, by addressing the scarcity of host information for microviruses, the current study significantly broadened their host range and discovered over 20 possible new hosts, including important pathogenic bacteria such as Helicobacter pylori and Vibrio cholerae, as well as different taxa demonstrated different host specificities. The findings of this study effectively expand the diversity of the Microviridae family, providing new insights for their classification and identification. Additionally, it offers a novel perspective for monitoring and controlling pathogenic microorganisms in poultry slaughterhouse environments.

Identification and Characterization of a Plant Endophytic Fungus Paraphaosphaeria sp. JRF11 and Its Growth-Promoting Effects

Endophytic fungi establish mutualistic relationships with host plants and can promote the growth and development of plants. In this study, the endophytic fungus JRF11 was isolated from Carya illinoinensis. Sequence analysis of the internal transcribed spacer (ITS) region and 18S rRNA gene combined with colonial and conidial morphology identified JRF11 as a Paraphaosphaeria strain. Plant-fungus interaction assays revealed that JRF11 showed significant growth-promoting effects on plants. In particular, JRF11 significantly increased the root biomass and soluble sugar content of plants. Furthermore, transcriptome analysis demonstrated that JRF11 treatment reprogrammed a variety of genes involved in plant mitogen-activated protein kinase (MAPK) signaling and starch and sucrose metabolism pathways through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Our research indicates that beneficial endophytic fungi are able to interact with plants and exhibit outstanding plant growth-promoting activities.

Comparison of the DNBSEQ platform and Illumina HiSeq 2000 for bacterial genome assembly

The Illumina HiSeq platform has been a commonly used option for bacterial genome sequencing. Now the BGI DNA nanoball (DNB) nanoarrays platform may provide an alternative platform for sequencing of bacterial genomes. To explore the impact of sequencing platforms on bacterial genome assembly, quality assessment, sequence alignment, functional annotation, mutation detection, and metagenome mapping, we compared genome assemblies based on sequencing of cultured bacterial species using the HiSeq 2000 and BGISEQ-500 platforms. In addition, simulated reads were used to evaluate the impact of insert size on genome assembly. Genome assemblies based on BGISEQ-500 sequencing exhibited higher completeness and fewer N bases in high GC genomes, whereas HiSeq 2000 assemblies exhibited higher N50. The majority of assembly assessment parameters, sequences of 16S rRNA genes and genomes, numbers of single nucleotide variants (SNV), and mapping to metagenome data did not differ significantly between platforms. More insertions were detected in HiSeq 2000 genome assemblies, whereas more deletions were detected in BGISEQ-500 genome assemblies. Insert size had no significant impact on genome assembly. Taken together, our results suggest that DNBSEQ platforms would be a valid substitute for HiSeq 2000 for bacterial genome sequencing.

Comparative transcriptomic analysis of Illumina and MGI next-generation sequencing platforms using RUNX3- and ZBTB46-instructed embryonic stem cells

Introduction: We have previously observed phenotypic and developmental changes upon the ectopic expression of the RUNX3 or the ZBTB46 transcription factors in mouse embryonic stem cell (ESC) derived progenitors. In this study, we evaluated the gene expression profiles of the RUNX3- and the ZBTB46-instructed murine ESCs with RNA-seq testing two next-generation sequencing technologies. Methods: We compared the DNA nanoball-based DNBSEQ G400 sequencer (MGI) with the bridge-PCR-based NextSeq 500 instrument (Illumina) for RNA sequencing. Moreover, we also compared two types of MGI sequencing reagents (Standard versus Hot-massive parallel sequencing (MPS)) with the DNBSEQ G400. Results: We observed that both sequencing platforms showed comparable levels of quality, sequencing uniformity, and gene expression profiles. For example, highly overlapping RUNX3- and ZBTB46-regulated gene lists were obtained from both sequencing datasets. Moreover, we observed that the Standard and the Hot-MPS-derived RUNX3- and ZBTB46-regulated gene lists were also considerably overlapped. This transcriptome analysis also helped us to identify differently expressed genes in the presence of the transgenic RUNX3 or ZBTB46. For example, we found that Gzmb, Gzmd, Gzme, Gdf6, and Ccr7 genes were robustly upregulated upon the forced expression of Runx3; on the other hand, Gpx2, Tdpoz4, and Arg2 were induced alongside the ectopic expression of Zbtb46. Discussion: Similar gene expression profile and greatly overlapping RUNX3- and ZBTB46-regulated gene sets were detected with both DNA sequencing platforms. Our analyses demonstrate that both sequencing technologies are suitable for transcriptome profiling and target gene selection. These findings suggest that DNBSEQ G400 represents a cost-effective alternative sequencing platform for gene expression monitoring. Moreover, this analysis provides a resource for exploration of the RUNX3- and ZBTB46-dependent gene regulatory networks.

A Network Landscape of HPVOPC Reveals Methylation Alterations as Significant Drivers of Gene Expression via an Immune-Mediated GPCR Signal

HPV-associated oropharynx carcinoma (HPVOPC) tumors have a relatively low mutational burden. Elucidating the relative contributions of other tumor alterations, such as DNA methylation alterations, alternative splicing events (ASE), and copy number variation (CNV), could provide a deeper understanding of carcinogenesis drivers in this disease. We applied network propagation analysis to multiple classes of tumor alterations in a discovery cohort of 46 primary HPVOPC tumors and 25 cancer-unaffected controls and validated our findings with TCGA data. We identified significant overlap between differential gene expression networks and all alteration classes, and this association was highest for methylation and lowest for CNV. Significant overlap was seen for gene clusters of G protein-coupled receptor (GPCR) pathways. HPV16-human protein interaction analysis identified an enriched cluster defined by an immune-mediated GPCR signal, including CXCR3 cytokines CXCL9, CXCL10, and CXCL11. CXCR3 was found to be expressed in primary HPVOPC, and scRNA-seq analysis demonstrated CXCR3 ligands to be highly expressed in M2 macrophages. In vivo models demonstrated decreased tumor growth with antagonism of the CXCR3 receptor in immunodeficient but not immunocompetent mice, suggesting that the CXCR3 axis can drive tumor proliferation in an autocrine fashion, but the effect is tempered by an intact immune system. In conclusion, methylation, ASE, and SNV alterations are highly associated with network gene expression changes in HPVOPC, suggesting that ASE and methylation alterations have an important role in driving the oncogenic phenotype. Network analysis identifies GPCR networks, specifically the CXCR3 chemokine axis, as modulators of tumor-immune interactions that may have proliferative effects on primary tumors as well as a role for immunosurveillance; however, CXCR3 inhibition should be used with caution, as these agents may both inhibit and stimulate tumor growth considering the competing effects of this cytokine axis. Further investigation is needed to explore opportunities for targeted therapy in this setting.

Cathepsin F genetic mutation is associated with familial papillary thyroid cancer

BACKGROUND

Thyroid cancer is one of the most common cancers in the world. Genetic factors are important in the occurrence and development of thyroid cancer, and genetic diagnosis has become an important basis for the prognosis of benign and malignant nodules. We identify a family of six siblings with inherited thyroid cancer susceptibility. All six members of this generation have been definitely diagnosed with papillary thyroid carcinoma. This work aims at confirming the relevant causative genes for thyroid cancer in this pedigree.

METHODS

We extract DNA from the peripheral blood of six individuals and perform whole genome sequencing. Sanger sequencing and immunohistochemistry further testify the cathepsin F (CTSF) mutation and expression.

RESULTS

We identify 57 single nucleotide variations (SNVs) out of at least 4 affected family members via certain filter criteria. The CTSF gene found in five of the six family members is here considered the most promising candidate gene mutation for familial thyroid cancer. Besides, our research also proves several known genes including CTSB, TEKT4, ESR1, MSH6, DIRC3, GNAS, and BANCR that act as probable oncogenic drivers in this family. The Sanger sequencing identifies the existence and veracity of CTSF somatic mutations. The CTSF immunohistochemistry of thyroid cancer tissue specimens displays that higher CTSF expression in mutated patients than that in wild-type patient as well as pericarcinomatous tissue.

CONCLUSIONS

We conclude that the evaluation of CTSF gene mutations of patients in thyroid cancer families may be predictive and valuable for the familial heredity of thyroid cancer.

Draft Genome Sequences of the Kocuria subflava Type Strain KCTC 39547 and Kocuria sp. Strain JC486, a Newly Isolated Strain from a Wild Ass Sanctuary in Gujarat, India

Here, we report a 2.86-Mbp genome sequence of Kocuria sp. strain JC486, which was isolated from a salt marsh, and a 3.03-Mbp sequence of the type strain Kocuria subflava KCTC 39547. Prediction from their genomes indicates that both strains are nonpathogenic.

Transcriptome Analysis of Tryptophan-Induced Resistance against Potato Common Scab

Potato common scab (CS) is a worldwide soil-borne disease that severely reduces tuber quality and market value. We observed that foliar application of tryptophan (Trp) could induce resistance against CS. However, the mechanism of Trp as an inducer to trigger host immune responses is still unclear. To facilitate dissecting the molecular mechanisms, the transcriptome of foliar application of Trp and water (control, C) was compared under Streptomyces scabies (S) inoculation and uninoculation. Results showed that 4867 differentially expressed genes (DEGs) were identified under S. scabies uninoculation (C-vs-Trp) and 2069 DEGs were identified under S. scabies inoculation (S-vs-S+Trp). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that Trp induced resistance related to the metabolic process, response to stimulus, and biological regulation. As phytohormone metabolic pathways related to inducing resistance, the expression patterns of candidate genes involved in salicylic acid (SA) and jasmonic acid/ethylene (JA/ET) pathways were analyzed using qRT-PCR. Their expression patterns showed that the systemic acquired resistance (SAR) and induced systemic resistance (ISR) pathways could be co-induced by Trp under S. scabies uninoculation. However, the SAR pathway was induced by Trp under S. scabies inoculation. This study will provide insights into Trp-induced resistance mechanisms of potato for controlling CS, and extend the application methods of Trp as a plant resistance inducer in a way that is cheap, safe, and environmentally friendly.

Formation of Blood Neutrophil Extracellular Traps Increases the Mastitis Risk of Dairy Cows During the Transition Period

The current study was conducted to analyze the functions of blood neutrophils in transition cows and their association with postpartum mastitis risk as indicated by somatic cell counts (SCCs) in milk. Seventy-six healthy Holstein dairy cows were monitored from Week 4 prepartum to Week 4 postpartum. Five dairy cows with low SCCs (38 ± 6.0 × 10³/mL) and five with high SCCs (3,753 ± 570.0 × 10³/mL) were selected based on milk SCCs during the first three weeks of lactation. At Week 1 pre- and postpartum, serum samples were obtained from each cow to measure neutrophil extracellular trap (NET)-related variables, and blood neutrophils were collected for transcriptome analysis by RNA sequencing. The serum concentration of NETs was significantly higher (P < 0.05) in cows with high SCCs than in cows with low SCCs (36.5 ± 2.92 vs. 18.4 ± 1.73 ng/mL). The transcriptomic analysis revealed that the transcriptome differences in neutrophils between high- and low-SCC cows were mainly in cell cycle-related pathways (42.6%), including the cell cycle, DNA damage, and chromosomal conformation, at Week 1 prepartum. The hub genes of these pathways were mainly involved in both the cell cycle and NETosis. These results indicated that the formation of NETs in the blood of transition dairy cows was different between cows with low and high SCCs, which may be used as a potential indicator for the prognosis of postpartum mastitis risk and management strategies of perinatal dairy cows.

The mutationathon highlights the importance of reaching standardization in estimates of pedigree-based germline mutation rates

In the past decade, several studies have estimated the human per-generation germline mutation rate using large pedigrees. More recently, estimates for various nonhuman species have been published. However, methodological differences among studies in detecting germline mutations and estimating mutation rates make direct comparisons difficult. Here, we describe the many different steps involved in estimating pedigree-based mutation rates, including sampling, sequencing, mapping, variant calling, filtering, and appropriately accounting for false-positive and false-negative rates. For each step, we review the different methods and parameter choices that have been used in the recent literature. Additionally, we present the results from a ‘Mutationathon,’ a competition organized among five research labs to compare germline mutation rate estimates for a single pedigree of rhesus macaques. We report almost a twofold variation in the final estimated rate among groups using different post-alignment processing, calling, and filtering criteria, and provide details into the sources of variation across studies. Though the difference among estimates is not statistically significant, this discrepancy emphasizes the need for standardized methods in mutation rate estimations and the difficulty in comparing rates from different studies. Finally, this work aims to provide guidelines for computational and statistical benchmarks for future studies interested in identifying germline mutations from pedigrees.

Canine Oral Melanoma Genomic and Transcriptomic Study Defines Two Molecular Subgroups with Different Therapeutical Targets

Simple Summary

In humans, mucosal melanoma (MM) is a rare and aggressive cancer. The canine model is frequently and spontaneously affected by MM, thus facilitating the collection of samples and the study of its genetic bases. Thanks to an integrative genomic and transcriptomic analysis of 32 canine MM samples, we identified two molecular subgroups of MM with a different microenvironment and structural variant (SV) content. We demonstrated that SVs are associated with recurrently amplified regions, and identified new candidate oncogenes (TRPM7, GABPB1, and SPPL2A) for MM. Our findings suggest the existence of two MM molecular subgroups that could benefit from dedicated therapies, such as immune checkpoint inhibitors or targeted therapies, for both human and veterinary medicine.

Abstract

Mucosal melanoma (MM) is a rare, aggressive clinical cancer. Despite recent advances in genetics and treatment, the prognosis of MM remains poor. Canine MM offers a relevant spontaneous and immunocompetent model to decipher the genetic bases and explore treatments for MM. We performed an integrative genomic and transcriptomic analysis of 32 canine MM samples, which identified two molecular subgroups with a different microenvironment and structural variant (SV) content. The overexpression of genes related to the microenvironment and T-cell response was associated with tumors harboring a lower content of SVs, whereas the overexpression of pigmentation-related pathways and oncogenes, such as TERT, was associated with a high SV burden. Using whole-genome sequencing, we showed that focal amplifications characterized complex chromosomal rearrangements targeting oncogenes, such as MDM2 or CDK4, and a recurrently amplified region on canine chromosome 30. We also demonstrated that the genes TRPM7, GABPB1, and SPPL2A, located in this CFA30 region, play a role in cell proliferation, and thus, may be considered as new candidate oncogenes for human MM. Our findings suggest the existence of two MM molecular subgroups that may benefit from dedicated therapies, such as immune checkpoint inhibitors or targeted therapies, for both human and veterinary medicine.

Assessing reproducibility of inherited variants detected with short-read whole genome sequencing

BACKGROUND

Reproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS.

RESULTS

To dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×.

CONCLUSIONS

Our findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

Evaluation of the MGISEQ-2000 Sequencing Platform for Illumina Target Capture Sequencing Libraries

Illumina is the leading sequencing platform in the next-generation sequencing (NGS) market globally. In recent years, MGI Tech has presented a series of new sequencers, including DNBSEQ-T7, MGISEQ-2000 and MGISEQ-200. As a complex application of NGS, cancer-detecting panels pose increasing demands for the high accuracy and sensitivity of sequencing and data analysis. In this study, we used the same capture DNA libraries constructed based on the Illumina protocol to evaluate the performance of the Illumina Nextseq500 and MGISEQ-2000 sequencing platforms. We found that the two platforms had high consistency in the results of hotspot mutation analysis; more importantly, we found that there was a significant loss of fragments in the 101-133 bp size range on the MGISEQ-2000 sequencing platform for Illumina libraries, but not for the capture DNA libraries prepared based on the MGISEQ protocol. This phenomenon may indicate fragment selection or low fragment ligation efficiency during the DNA circularization step, which is a unique step of the MGISEQ-2000 sequence platform. In conclusion, these different sequencing libraries and corresponding sequencing platforms are compatible with each other, but protocol and platform selection need to be carefully evaluated in combination with research purpose.

Extrachromosomal DNA in HPV-Mediated Oropharyngeal Cancer Drives Diverse Oncogene Transcription

PURPOSE

Human papillomavirus (HPV) plays a major role in oncogenesis and circular extrachromosomal DNA (ecDNA) is found in many cancers. However, the relationship between HPV and circular ecDNA in human cancer is not understood.

EXPERIMENTAL DESIGN

Forty-four primary tumor tissue samples were obtained from a cohort of patients with HPV-positive oropharynx squamous cell carcinoma (OPSCC). Twenty-eight additional HPV oropharyngeal cancer (HPVOPC) tumors from The Cancer Genome Atlas (TCGA) project were analyzed as a separate validation cohort. Genomic, transcriptomic, proteomic, computational, and functional analyses of HPVOPC were applied to these datasets.

RESULTS

Our analysis revealed circular, oncogenic DNA in nearly all HPVOPC, with circular human and human-viral hybrid ecDNA present in over a third of HPVOPC and viral circular DNA in remaining tumors. Hybrid ecDNA highly express fusion transcripts from HPV promoters and HPV oncogenes linked to downstream human transcripts that drive oncogenic transformation and immune evasion, and splice multiple, diverse human acceptors to a canonical SA880 viral donor site. HPVOPC have high E6*I expression with specific viral oncogene expression pattern related to viral or hybrid ecDNA composition.

CONCLUSIONS

Nonchromosomal circular oncogenic DNA is a dominant feature of HPVOPC, revealing an unanticipated link between HPV and ecDNA that leverages the power of extrachromosomal inheritance to drive HPV and somatic oncogene expression.

Comparison between MGI and Illumina sequencing platforms for whole genome sequencing

BACKGROUND

Illumina next generation sequencing (NGS) systems are the major sequencing platform in worldwide next-generation sequencing market. On the other hand, MGI Tech launched a series of new NGS equipment that promises to deliver high-quality sequencing data faster and at lower prices than Illumina's sequencing instruments.

OBJECTIVE

In this study, we compared the performance of the two platform's major sequencing instruments-Illumina's NovaSeq 6000 and MGI's MGISEQ-2000 and DNBSEQ-T7-to test whether the MGISEQ-2000 and DNBSEQ-T7 sequencing instruments are also suitable for whole genome sequencing.

METHODS

We sequenced two pairs of normal and tumor tissues from Korean lung cancer patients using the three platforms. Then, we called single nucleotide variants (SNVs) and insertion and deletion (indels) for somatic and germline variants to compare the performance among the three platforms.

RESULTS

In quality control analysis, all of the three platforms showed high-quality scores and deep coverages. Comparison among the three platforms revealed that MGISEQ-2000 is most concordant with NovaSeq 6000 for germline SNVs and indels, and DNBSEQ-T7 is most concordant with NovaSeq 6000 for somatic SNVs and indels.

CONCLUSIONS

These results suggest that the performances of the MGISEQ-2000 and DNBSEQ-T7 platforms are comparable to that of the Illumina NovaSeq 6000 platform and support the potential applicability of the MGISEQ-2000 and DNBSEQ-T7 platforms in actual genome analysis fields.

Comparative analysis of 7 short-read sequencing platforms using the Korean Reference Genome: MGI and Illumina sequencing benchmark for whole-genome sequencing

BACKGROUND

DNBSEQ-T7 is a new whole-genome sequencer developed by Complete Genomics and MGI using DNA nanoball and combinatorial probe anchor synthesis technologies to generate short reads at a very large scale-up to 60 human genomes per day. However, it has not been objectively and systematically compared against Illumina short-read sequencers.

FINDINGS

By using the same KOREF sample, the Korean Reference Genome, we have compared 7 sequencing platforms including BGISEQ-500, DNBSEQ-T7, HiSeq2000, HiSeq2500, HiSeq4000, HiSeqX10, and NovaSeq6000. We measured sequencing quality by comparing sequencing statistics (base quality, duplication rate, and random error rate), mapping statistics (mapping rate, depth distribution, and percent GC coverage), and variant statistics (transition/transversion ratio, dbSNP annotation rate, and concordance rate with single-nucleotide polymorphism [SNP] genotyping chip) across the 7 sequencing platforms. We found that MGI platforms showed a higher concordance rate for SNP genotyping than HiSeq2000 and HiSeq4000. The similarity matrix of variant calls confirmed that the 2 MGI platforms have the most similar characteristics to the HiSeq2500 platform.

CONCLUSIONS

Overall, MGI and Illumina sequencing platforms showed comparable levels of sequencing quality, uniformity of coverage, percent GC coverage, and variant accuracy; thus we conclude that the MGI platforms can be used for a wide range of genomics research fields at a lower cost than the Illumina platforms.

ReSeq simulates realistic Illumina high-throughput sequencing data

In high-throughput sequencing data, performance comparisons between computational tools are essential for making informed decisions at each step of a project. Simulations are a critical part of method comparisons, but for standard Illumina sequencing of genomic DNA, they are often oversimplified, which leads to optimistic results for most tools. ReSeq improves the authenticity of synthetic data by extracting and reproducing key components from real data. Major advancements are the inclusion of systematic errors, a fragment-based coverage model and sampling-matrix estimates based on two-dimensional margins. These improvements lead to more faithful performance evaluations. ReSeq is available at https://github.com/schmeing/ReSeq.

Abundancy of polymorphic CGG repeats in the human genome suggest a broad involvement in neurological disease

Expanded CGG-repeats have been linked to neurodevelopmental and neurodegenerative disorders, including the fragile X syndrome and fragile X-associated tremor/ataxia syndrome (FXTAS). We hypothesized that as of yet uncharacterised CGG-repeat expansions within the genome contribute to human disease. To catalogue the CGG-repeats, 544 human whole genomes were analyzed. In total, 6101 unique CGG-repeats were detected of which more than 93% were highly variable in repeat length. Repeats with a median size of 12 repeat units or more were always polymorphic but shorter repeats were often polymorphic, suggesting a potential intergenerational instability of the CGG region even for repeats units with a median length of four or less. 410 of the CGG repeats were associated with known neurodevelopmental disease genes or with strong candidate genes. Based on their frequency and genomic location, CGG repeats may thus be a currently overlooked cause of human disease.

Family-Based Whole Genome Sequencing Identified Novel Variants in ABCA5 Gene in a Patient with Idiopathic Ventricular Tachycardia

Idiopathic ventricular tachycardia (IVT) is the major cause of sudden cardiac death. Patients with IVT were usually manifested without structural heart disease. In this present study, we performed family-based whole genome sequencing (WGS) and Sanger sequencing for a 5-year-old Chinese boy with IVT and all the unaffected family members in order to identify the candidate gene and disease-causing mutation underlying the disease phenotype. Results showed that a novel heterozygous single-nucleotide duplication (c.128dup) and a novel heterozygous missense (c.3328A > G) variant in ABCA5 gene were identified in the proband. The single-nucleotide duplication (c.128dupT), inherited from his father and patrilineal grandfather, leads to a frameshift which results into the formation of a truncated ABCA5 protein of 50 (p.Leu43Phefs*8) amino acids. Hence, it is a loss-of-function mutation. The missense (c.3328A > G) variant, inherited from his mother, leads to the replacement of isoleucine by valine at the position of 1110 (p.Ile1110Val) of the ABCA5 protein. Multiple sequence alignment showed that p.Ile1110 is evolutionarily conserved among several species indicating both the structural and functional significance of the p.Ile1110 residue in the wild-type ABCA5 protein. Quantitative RT-PCR showed that the ABCA5 mRNA expression levels were decreased in the proband. These two novel variants of ABCA5 gene were co-segregated well among all the members of this family. Our present study also strongly supports the importance of using family-based whole genome sequencing for identifying novel candidate genes associated with IVT.

Cytogenetic and molecular genetic methods for chromosomal translocations detection with reference to the KMT2A/MLL gene

Acute leukemias (ALs) are often associated with chromosomal translocations, in particular, KMT2A/MLL gene rearrangements. Identification or confirmation of these translocations is carried out by a number of genetic and molecular methods, some of which are routinely used in clinical practice, while others are primarily used for research purposes. In the clinic, these methods serve to clarify diagnoses and monitor the course of disease and therapy. On the other hand, the identification of new translocations and the confirmation of known translocations are of key importance in the study of disease mechanisms and further molecular classification. There are multiple methods for the detection of rearrangements that differ in their principle of operation, the type of problem being solved, and the cost-result ratio. This review is intended to help researchers and clinicians studying AL and related chromosomal translocations to navigate this variety of methods. All methods considered in the review are grouped by their principle of action and include karyotyping, fluorescence in situ hybridization (FISH) with probes for whole chromosomes or individual loci, PCR and reverse transcription-based methods, and high-throughput sequencing. Another characteristic of the described methods is the type of problem being solved. This can be the discovery of new rearrangements, the determination of unknown partner genes participating in the rearrangement, or the confirmation of the proposed rearrangement between the two genes. We consider the specifics of the application, the basic principle of each method, and its pros and cons. To illustrate the application, examples of studying the rearrangements of the KMT2A/MLL gene, one of the genes that are often rearranged in AL, are mentioned.

Advantages and Limitations of 16S rRNA Next-Generation Sequencing for Pathogen Identification in the Diagnostic Microbiology Laboratory: Perspectives from a Middle-Income Country

Bacterial culture and biochemical testing (CBtest) have been the cornerstone of pathogen identification in the diagnostic microbiology laboratory. With the advent of Sanger sequencing and later, next-generation sequencing, 16S rRNA next-generation sequencing (16SNGS) has been proposed to be a plausible platform for this purpose. Nevertheless, usage of the 16SNGS platform has both advantages and limitations. In addition, transition from the traditional methods of CBtest to 16SNGS requires procurement of costly equipment, timely and sustainable maintenance of these platforms, specific facility infrastructure and technical expertise. All these factors pose a challenge for middle-income countries, more so for countries in the lower middle-income range. In this review, we describe the basis for CBtest and 16SNGS, and discuss the limitations, challenges, advantages and future potential of using 16SNGS for bacterial pathogen identification in diagnostic microbiology laboratories of middle-income countries.

Novel GLI3 Mutations in Chinese Patients with Non-syndromic Post-axial Polydactyly

BACKGROUND

Polydactyly, characterized by supernumerary digits in the upper or lower extremities, is the most common congenital digital abnormalities. It derives from the defective patterning of anteroposterior axis of the developing limb, with various etiology and clinical heterogeneity. The patients with post-axial polydactyly type A (PAPA) have the typical symptom of a well-formed supernumerary digit outside the fifth digit.

OBJECTIVE

The aim of present study was to identify the causative mutations of two unrelated Han Chinese patients with non-syndromic PAPA.

METHODS

Two unrelated Han Chinese patients and 100 ethnicity-matched, unrelated normal controls were recruited for this study. BGISEQ-500 exome sequencing was performed in the two patients, followed by validation in the patients and 100 controls by using Sanger sequencing.

RESULTS

Two mutations in the GLI family zinc finger 3 gene (GLI3), including a frameshift mutation c.3437_3453delTCGAGCAGCCCTGCCCC (p.L1146RfsX95) and a nonsense mutation c.3997C>T (p.Q1333X), were identified in two patients but were absent in the 100 healthy controls.

CONCLUSION

The two GLI3 mutations, p.L1146RfsX95 and p.Q1333X, may account for non-syndromic PAPA in the two patients, respectively. The findings of this study may expand the mutational spectrum of GLI3-PAPA and provide novel insights into the genetic basis of polydactyly.

Estimating the costs of genomic sequencing in cancer control

BACKGROUND

Despite the rapid uptake of genomic technologies within cancer care, few studies provide detailed information on the costs of sequencing across different applications. The objective of the study was to examine and categorise the complete costs involved in genomic sequencing for a range of applications within cancer settings.

METHODS

We performed a cost-analysis using gross and micro-costing approaches for genomic sequencing performed during 2017/2018 across different settings in Brisbane, Australia. Sequencing was undertaken for patients with lung, breast, oesophageal cancers, melanoma or mesothelioma. Aggregated resource data were captured for a total of 1433 patients and point estimates of per patient costs were generated. Deterministic sensitivity analyses addressed the uncertainty in the estimates. Estimated costs to the public health system for resources were categorised into seven distinct activities in the sequencing process: sampling, extraction, library preparation, sequencing, analysis, data storage and clinical reporting. Costs were also aggregated according to labour, consumables, testing, equipment and 'other' categories.

RESULTS

The per person costs were AU$347-429 (2018 US$240-297) for targeted panels, AU$871-$2788 (2018 US$604-1932) for exome sequencing, and AU$2895-4830 (2018 US$2006-3347) for whole genome sequencing. Cost proportions were highest for library preparation/sequencing materials (average 76.8% of total costs), sample extraction (8.1%), data analysis (9.2%) and data storage (2.6%). Capital costs for the sequencers were an additional AU$34-197 (2018 US$24-67) per person.

CONCLUSIONS

Total costs were most sensitive to consumables and sequencing activities driven by commercial prices. Per person sequencing costs for cancer are high when tumour/blood pairs require testing. Using the natural steps involved in sequencing and categorising resources accordingly, future evaluations of costs or cost-effectiveness of clinical genomics across cancer projects could be more standardised and facilitate easier comparison of cost drivers.

Genome sequencing of human in vitro fertilisation embryos for pathogenic variation screening

Whole-genome sequencing of preimplantation human embryos to detect and screen for genetic diseases is a technically challenging extension to preconception screening. Combining preconception genetic screening with preimplantation testing of human embryos facilitates the detection of de novo mutations and self-validates transmitted variant detection in both the reproductive couple and the embryo’s samples. Here we describe a trio testing workflow that involves whole-genome sequencing of amplified DNA from biopsied embryo trophectoderm cells and genomic DNA from both parents. Variant prediction software and annotation databases were used to assess variants of unknown significance and previously not described de novo variants in five single-gene preimplantation genetic testing couples and eleven of their embryos. Pathogenic variation, tandem repeat, copy number and structural variations were examined against variant calls for compound heterozygosity and predicted disease status was ascertained. Multiple trio testing showed complete concordance with known variants ascertained by single-nucleotide polymorphism array and uncovered de novo and transmitted pathogenic variants. This pilot study describes a method of whole-genome sequencing and analysis for embryo selection in high-risk couples to prevent early life fatal genetic conditions that adversely affect the quality of life of the individual and families.

Sub-Exome Target Sequencing in a Family With Syndactyly Type IV Due to a Novel Partial Duplication of the LMBR1 Gene: First Case Report in Fujian Province of China

Syndactyly is one of the most frequent hereditary limb malformations with clinical and genetical complexity. Autosomal dominant syndactyly type IV (SD4) is a rare form of syndactyly, caused by heterozygous mutations in a sonic hedgehog (SHH) regulatory element (ZRS) which resides in intron 5 of the LMBR1 gene on chromosome 7q36.3. SD4 is characterized by complete cutaneous syndactyly of the fingers, accompanied by cup-shaped hands due to flexion of the fingers and polydactyly. Here, for the first time, we reported a large Chinese family from Fujian province, manifesting cup-shaped hands consistent with SD4 and intrafamilial heterogeneity in clinical phenotype of tibial and fibulal shortening, triphalangeal thumb-polysyndactyly syndrome (TPTPS). We identified a novel duplication of ∼222 kb covering exons 2–17 of the LMBR1 gene in this family by sub-exome target sequencing. This case expands our new clinical understanding of SD4 phenotype and again confirms the feasibility to detect copy number variation by sub-exome target sequencing.

Molecular digitization of a botanical garden: high-depth whole-genome sequencing of 689 vascular plant species from the Ruili Botanical Garden

Background

Genome sequencing has been widely used in plant research to construct reference genomes and provide evolutionary insights. However, few plant species have had their whole genome sequenced, thus restraining the utility of these data. We collected 1,093 samples of vascular plant species growing in the Ruili Botanical Garden, located in southwest China. Of these, we sequenced 761 samples and collected voucher specimens stored in the Herbarium of China National GeneBank.

Results

The 761 sequenced samples represented 689 vascular plant species from 137 families belonging to 49 orders. Of these, 257 samples were identified to the species level and 504 to the family level, using specimen and chloroplast sequences. In total, we generated 54 Tb of sequencing data, with an average sequencing depth of 60X per species, as estimated from genome sizes. A reference phylogeny was reconstructed with 78 chloroplast genes for molecular identification and other possible applications.

Conclusions

The large dataset of vascular plant genomes generated in this study, which includes both high-depth whole-genome sequencing data and associated voucher specimens, is valuable for plant genome research and other applications. This project also provides insight into the feasibility and technical requirements for “planetary-scale” projects such as the 10,000 Plant Genomes Project and the Earth BioGenome Project.

SEQdata-BEACON: a comprehensive database of sequencing performance and statistical tools for performance evaluation and yield simulation in BGISEQ-500

BACKGROUND

The sequencing platform BGISEQ-500 is based on DNBSEQ technology and provides high throughput with low costs. This sequencer has been widely used in various areas of scientific and clinical research. A better understanding of the sequencing process and performance of this system is essential for stabilizing the sequencing process, accurately interpreting sequencing results and efficiently solving sequencing problems. To address these concerns, a comprehensive database, SEQdata-BEACON, was constructed to accumulate the run performance data in BGISEQ-500.

RESULTS

A total of 60 BGISEQ-500 instruments in the BGI-Wuhan lab were used to collect sequencing performance data. Lanes in paired-end 100 (PE100) sequencing using 10 bp barcode were chosen, and each lane was assigned a unique entry number as its identification number (ID). From November 2018 to April 2019, 2236 entries were recorded in the database containing 65 metrics about sample, yield, quality, machine state and supplies information. Using a correlation matrix, 52 numerical metrics were clustered into three groups signifying yield-quality, machine state and sequencing calibration. The distributions of the metrics also delivered information about patterns and rendered clues for further explanation or analysis of the sequencing process. Using the data of a total of 200 cycles, a linear regression model well simulated the final outputs. Moreover, the predicted final yield could be provided in the 15th cycle of the early stage of sequencing, and the corresponding R2 of the 200th and 15th cycle models were 0.97 and 0.81, respectively. The model was run with the test sets obtained from May 2019 to predict the yield, which resulted in an R2 of 0.96. These results indicate that our simulation model was reliable and effective.

CONCLUSIONS

Data sources, statistical findings and application tools provide a constantly updated reference for BGISEQ-500 users to comprehensively understand DNBSEQ technology, solve sequencing problems and optimize run performance. These resources are available on our website http://seqBEACON.genomics.cn:443/home.html.

Comparison of the MGISEQ-2000 and Illumina HiSeq 4000 sequencing platforms for RNA sequencing

Currently, Illumina sequencers are the globally leading sequencing platform in the next-generation sequencing market. Recently, MGI Tech launched a series of new sequencers, including the MGISEQ-2000, which promise to deliver high-quality sequencing data faster and at lower prices than Illumina's sequencers. In this study, we compared the performance of two major sequencers (MGISEQ-2000 and HiSeq 4000) to test whether the MGISEQ-2000 sequencer delivers high-quality sequence data as suggested. We performed RNA sequencing of four human colon cancer samples with the two platforms, and compared the sequencing quality and expression values. The data produced from the MGISEQ-2000 and HiSeq 4000 showed high concordance, with Pearson correlation coefficients ranging from 0.98 to 0.99. Various quality control (QC) analyses showed that the MGISEQ-2000 data fulfilled the required QC measures. Our study suggests that the performance of the MGISEQ-2000 is comparable to that of the HiSeq 4000 and that the MGISEQ-2000 can be a useful platform for sequencing.

Impact of sequencing depth and technology on de novo RNA-Seq assembly

BACKGROUND

RNA-Seq data is inherently nonuniform for different transcripts because of differences in gene expression. This makes it challenging to decide how much data should be generated from each sample. How much should one spend to recover the less expressed transcripts? The sequencing technology used is another consideration, as there are inevitably always biases against certain sequences. To investigate these effects, we first looked at high-depth libraries from a set of well-annotated organisms to ascertain the impact of sequencing depth on de novo assembly. We then looked at libraries sequenced from the Universal Human Reference RNA (UHRR) to compare the performance of Illumina HiSeq and MGI DNBseq™ technologies.

RESULTS

On the issue of sequencing depth, the amount of exomic sequence assembled plateaued using data sets of approximately 2 to 8 Gbp. However, the amount of genomic sequence assembled did not plateau for many of the analyzed organisms. Most of the unannotated genomic sequences are single-exon transcripts whose biological significance will be questionable for some users. On the issue of sequencing technology, both of the analyzed platforms recovered a similar number of full-length transcripts. The missing "gap" regions in the HiSeq assemblies were often attributed to higher GC contents, but this may be an artefact of library preparation and not of sequencing technology.

CONCLUSIONS

Increasing sequencing depth beyond modest data sets of less than 10 Gbp recovers a plethora of single-exon transcripts undocumented in genome annotations. DNBseq™ is a viable alternative to HiSeq for de novo RNA-Seq assembly.

Systematic comparison of germline variant calling pipelines cross multiple next-generation sequencers

The development and innovation of next generation sequencing (NGS) and the subsequent analysis tools have gain popularity in scientific researches and clinical diagnostic applications. Hence, a systematic comparison of the sequencing platforms and variant calling pipelines could provide significant guidance to NGS-based scientific and clinical genomics. In this study, we compared the performance, concordance and operating efficiency of 27 combinations of sequencing platforms and variant calling pipelines, testing three variant calling pipelines—Genome Analysis Tool Kit HaplotypeCaller, Strelka2 and Samtools-Varscan2 for nine data sets for the NA12878 genome sequenced by different platforms including BGISEQ500, MGISEQ2000, HiSeq4000, NovaSeq and HiSeq Xten. For the variants calling performance of 12 combinations in WES datasets, all combinations displayed good performance in calling SNPs, with their F-scores entirely higher than 0.96, and their performance in calling INDELs varies from 0.75 to 0.91. And all 15 combinations in WGS datasets also manifested good performance, with F-scores in calling SNPs were entirely higher than 0.975 and their performance in calling INDELs varies from 0.71 to 0.93. All of these combinations manifested high concordance in variant identification, while the divergence of variants identification in WGS datasets were larger than that in WES datasets. We also down-sampled the original WES and WGS datasets at a series of gradient coverage across multiple platforms, then the variants calling period consumed by the three pipelines at each coverage were counted, respectively. For the GIAB datasets on both BGI and Illumina platforms, Strelka2 manifested its ultra-performance in detecting accuracy and processing efficiency compared with other two pipelines on each sequencing platform, which was recommended in the further promotion and application of next generation sequencing technology. The results of our researches will provide useful and comprehensive guidelines for personal or organizational researchers in reliable and consistent variants identification.

Novel and Recurring NOTCH3 Mutations in Two Chinese Patients with CADASIL

BACKGROUND

Cerebral autosomal-dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is an autosomal-dominant, inherited, systemic, vascular disorder primarily involving the small arteries. It is characterized by migraine, recurrent ischemic strokes, cognitive decline, and dementia. Mutations in the Notch receptor 3 gene (NOTCH3) and the HtrA serine peptidase 1 gene (HTRA1) are 2 genetic causes for CADASIL. The NOTCH3 gene, located on chromosome 19p13.12, is the most common disease-causing gene in CADASIL.

OBJECTIVE

To investigate genetic causes in 2 unrelated Han-Chinese patients with presentations strongly suggestive of CADASIL.

METHODS

Exome sequencing was performed on both patients and potential pathogenic mutations were validated by Sanger sequencing.

RESULTS

This study reports on 2 unrelated Han-Chinese patients with presentations strongly suggestive of CADASIL, identifying that NOTCH3 mutations were the genetic cause. A common mutation, c.268C>T (p.Arg90Cys), and a novel mutation, c.331G>T (p.Gly111Cys) in the NOTCH3 gene, were detected and confirmed in the patients, respectively, and were predicted to be deleterious based on bioinformation analyses.

CONCLUSIONS

We identified 2 NOTCH3 mutations as likely genetic causes for CADASIL in these 2 patients. Our findings broaden the mutational spectrum of the NOTCH3 gene accountable for CADASIL. Clinical manifestations supplemented with molecular genetic analyses are critical for accurate diagnosis, the provision of genetic counseling, and the development of therapies for CADASIL.

Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan

The primary goal of cattle genomics is the identification of genome-wide polymorphism associated with economically important traits. The bovine genome sequencing project was completed in 2009. Since then, using massively parallel sequencing technologies, a large number of Bos taurus cattle breeds have been resequenced and scanned for genome-wide polymorphisms. As a result, a substantial number of single nucleotide polymorphisms (SNPs) have been discovered across European Bos taurus genomes, whereas extremely less number of SNPs are cataloged for Bos indicus breeds. In this study, we performed whole-genome resequencing, reference-based mapping, functional annotation and gene enrichment analysis of 20 sires representing eleven important Bos indicus (indicine) breeds of Pakistan. The breeds sequenced here include: Sahiwal, Red Sindhi, Tharparkar and Cholistani (tropically adapted dairy and dual purpose breeds), Achai, Bhagnari, Dajal and Lohani (high altitude adapted dual and drought purpose breeds); Dhanni, Hisar Haryana and Gabrali (dairy and light drought purpose breeds). A total of 17.4 billion QC passed reads were produced using BGISEQ-500 next generation sequencing platform to generate 9 to 27-fold genome coverage (average ~16×) for each of the 20 sequenced sires. A total of 67,303,469 SNPs were identified, of which 3,850,365 were found novel and 1,083,842 insertions-deletions (InDels) were detected across the whole sequenced genomes (491,247 novel). Comparative analysis using coding region SNPs revealed a close relationship between the best milking indicine breeds; Red Sindhi and Sahiwal. On the other hand, Bhagnari and Tharparkar being popular for their adaptation to dry and extremely hot climates were found to share the highest number of SNPs. Functional annotation identified a total of 3,194 high-impact (disruptive) SNPs and 745 disruptive InDels (in 275 genes) that may possibly affect economically important dairy and beef traits. Functional enrichment analysis was performed and revealed that high or moderate impact variants in wingless-related integration site (Wnt) and vascular smooth muscle contraction (VSMC) signaling pathways were significantly over-represented in tropically adapted heat tolerant Pakistani-indicine breeds. On the other hand, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 (HIF-1) signaling pathways were found over-represented in highland adapted Pakistani-indicine breeds. Similarly, the ECM-receptor interaction and Jak-STAT signaling pathway were significantly enriched in dairy and beef purpose Pakistani-indicine cattle breeds. The Toll-like receptor signaling pathway was significantly enriched in most of the Pakistani-indicine cattle. Therefore, this study provides baseline data for further research to investigate the molecular mechanisms of major traits and to develop potential genomic markers associated with economically important breeding traits, particularly in indicine cattle.

Mutation of chromatin regulators and focal hotspot alterations characterize human papillomavirus-positive oropharyngeal squamous cell carcinoma

BACKGROUND

Human papillomavirus (HPV)-associated oropharyngeal cancer is a disease clinically and biologically distinct from smoking-related head and neck squamous cell carcinoma (HNSCC). Despite its rapidly increasing incidence, the mutational landscape of HPV+ oropharyngeal squamous cell carcinoma (OPSCC) remains understudied.

METHODS

This article presents the first mutational analysis of the 46 HPV+ OPSCC tumors within the newly expanded cohort of 530 HNSCC tumors from The Cancer Genome Atlas. A separate exome sequencing analysis was also performed for 46 HPV+ OPSCCs matched to their normal lymphocyte controls from the Johns Hopkins University cohort.

RESULTS

There was a strikingly high 33% frequency of mutations within genes associated with chromatin regulation, including mutations in lysine methyltransferase 2C (KMT2C), lysine methyltransferase 2D (KMT2D), nuclear receptor binding SET domain protein 1 (NSD1), CREB binding protein (CREBBP), E1A-associated protein p300 (EP300), and CCCTC-binding factor (CTCF). In addition, the commonly altered genes phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit α (PIK3CA) and fibroblast growth factor receptor 3 (FGFR3) showed distinct domain-specific hotspot mutations in comparison with their HPV- counterparts. PIK3CA showed a uniquely high rate of mutations within the helicase domain, and FGFR3 contained a predominance of hotspot S249C alterations that were not found in HPV- HNSCC.

CONCLUSIONS

This analysis represents one of the largest studies to date of HPV+ OPSCC and lends novel insight into the genetic landscape of this biologically distinct disease, including a high rate of mutations in histone- and chromatin-modifying genes, which may offer novel therapeutic targets.

Actualités du séquençage d’ADN

The planned acquisition of Pacific Biosciences by Illumina could signal a big change in the DNA sequencing landscape, and provides the incentive for an update on a field that still has a number of players and could be disrupted by innovative technologies.

Expanding the Clinical Spectrum of Osteogenesis Imperfecta Type V: 13 Additional Patients and Review

Osteogenesis imperfecta (OI) is an inherited connective tissue disorder characterized by bone fragility and is characterized by clinical and genetic heterogeneity. Previous studies showed that the same mutation (c.-14C> T) of the IFITM5 gene is responsible for autosomal dominant OI type V. However, the mutation has a variable expressivity. Clinical heterogeneity has been recognized in OI type V. In this study, we investigated 13 individuals with molecularly confirmed OI type V from seven Chinese families and explored the genotype-phenotype relationship. Increased callus formation is not observed in all individuals, and several novel clinical features were described: joint contractures (three individuals) and unexplained hip arthritis (six individuals). Significant clinical variability was observed even within families. Specific facial features were observed in six individuals from two families consistent with the facial features associated with OI type V reported so far in the literature. Interestingly, we report the process of hypertrophic callus formation in detail for the first time, and in five individuals with hyperplastic callus, increased erythrocyte sedimentation rate (ESR) and levels of C-reactive protein (C-RP) were measured, suggestive of inflammatory activation.

The application of NIPT using combinatorial probe-anchor synthesis to identify sex chromosomal aneuploidies (SCAs) in a cohort of 570 pregnancies

Background

Non-invasive prenatal testing (NIPT) as alternative screening method had been proven to have very high sensitivity and specificity for detecting common aneuploidies such as T21, T18, and T13, with low false positive and false negative rates. Unfortunately, recent studies suggested that the NIPT achieved lower accuracy in sex chromosomal aneuploidies (SCAs) detection than autosomal aneuploidies detection. BGISEQ-500 powered by Combinatorial Probe-Anchor Synthesis (CPAS) and DNA Nanoballs (DNBs) technology that combined linear amplification and rolling circle replication to reduce the error rate while enhancing the signal. Therefore, NIPT based on CPAS might be a good method for SCAs screening in routine clinical practice. In the study, we intended to evaluate the clinical utility of NIPT based on CPAS on screening for fetal SCAs.

Results

A total of 570 pregnant women were included in the retrospective study. Maternal blood samples were collected for NIPT; amniocentesis was performed on all pregnant women. NIPT was carried out by BGISEQ-500 sequencing platform based on CPAS. Karyotype analysis of amniotic cells was performed by standard G-banding techniques. 43 out of the total 570 pregnant women tested by NIPT showed fetal SCAs (19 of 45,X, 12 of 47,XXY, 10 of 47,XXX, and 2 of 47,XYY). The following amniocentesis confirmed that 26 cases were true positive (7 of true positive 45,X, 9 of true positive 47,XXY, 9 of true positive 47,XXX as well as 1 of 47,XYY) and the positive predictive value (PPV) for fetal SCAs was 60.47%. In addition, the PPV of advanced maternal age group (67.74%) was higher than the other indications group (45.45%) or serological screening high-risk /critical-risk group (0%).

Conclusions

NIPT based on CPAS could be a potential method for SCAs screening. However, it still had high false positive rates, especially for 45,X. The pregnant women with fetal SCAs detected by NIPT, especially those with non-age-related prenatal diagnostic indications, should be advised to accept invasive prenatal karyotype analysis.

Genomics and Epigenetics of Malignant Mesothelioma

Malignant mesothelioma is an aggressive and lethal asbestos-related disease. Diagnosis of malignant mesothelioma is particularly challenging and is further complicated by the lack of disease subtype-specific markers. As a result, it is especially difficult to distinguish malignant mesothelioma from benign reactive mesothelial proliferations or reactive fibrosis. Additionally, mesothelioma diagnoses can be confounded by other anatomically related tumors that can invade the pleural or peritoneal cavities, collectively resulting in delayed diagnoses and greatly affecting patient management. High-throughput analyses have uncovered key genomic and epigenomic alterations driving malignant mesothelioma. These molecular features have the potential to better our understanding of malignant mesothelioma biology as well as to improve disease diagnosis and patient prognosis. Genomic approaches have been instrumental in identifying molecular events frequently occurring in mesothelioma. As such, we review the discoveries made using high-throughput technologies, including novel insights obtained from the analysis of the non-coding transcriptome, and the clinical potential of these genetic and epigenetic findings in mesothelioma. Furthermore, we aim to highlight the potential of these technologies in the future clinical applications of the novel molecular features in malignant mesothelioma.