CNV-seq, a new method to detect copy number variation using high-throughput sequencing

UBTF haploinsufficiency associated with UBTF-related global developmental delay and distinctive facial features without neuroregression

BACKGROUND

The Upstream Binding Transcription Factor (UBTF) gene encodes two nucleolar proteins, UBTF1 and UBTF2. UBTF1 regulates rRNA transcription by RNA polymerase I, while UBTF2 regulates mRNA transcription by RNA polymerase II. A recurrent de novo dominant mutation c.628G>A (p.Glu210Lys) has been identified as a gain-of-function mutation associated with childhood onset neurodegeneration with brain atrophy (CONDBA). Evidence from large-scale population databases and Ubtf+/- mouse models indicates that UBTF haploinsufficiency is not tolerated.

METHODS

Three unrelated patients with global developmental delay and distinctive facial features were recruited for the study. Whole exome sequencing (WES) was performed to identify potential genetic abnormalities. Additionally, copy number variation analysis was conducted based on the WES data.

RESULTS

All three patients exhibited intellectual disabilities, social challenges and developmental delays in language and gross motor skills. Distinctive facial features included a wide forehead, sparse eyebrows, hypertelorism, narrow palpebral fissures, single-fold eyelids, a flat nasal bridge, anteverted nares, a long philtrum and a thin upper lip. Additionally, patient C presented with more severe language delay, recurrent hepatic dysfunction and an atrial septal defect. Patient A was found to have a nonsense variant, c.1327C>T (p.R443Ter), in the exon 13 of UBTF. Patients B and C both carried a heterozygous deletion encompassing the UBTF gene.

CONCLUSION

In this study, we analysed the detailed phenotypes associated with UBTF haploinsufficiency, which, to our knowledge, have not been previously reported. We propose that UBTF haploinsufficiency-related global developmental delay and distinctive facial features, without neuroregression, constitute a new syndrome distinct from CONDBA.

Propagation path of a flowering cherry (Cerasus × yedoensis) cultivar 'Somei-Yoshino' traced by somatic mutations

In the long history of human relations with flowering cherry trees in Japan, 'Somei-Yoshino' occupies an exceptional position among a variety of flowering trees: it is a self-incompatible interspecific hybrid but has been enthusiastically planted by grafting throughout Japan, due most likely to its flamboyant appearance upon full bloom. Thus, 'Somei-Yoshino' gives us a rare opportunity to trace and investigate the occurrence and distribution of somatic mutations within a single plant species through analysis of the genomes of the clonally propagated trees grown under a variety of geographical and artificial environments. In the studies presented here, a total of 46 samples of 'Somei-Yoshino' trees were collected and their genomes were analysed. We identified 684 single nucleotide mutations, of which 71 were present in more than two samples. Clustering analysis of the mutations indicated that the 46 samples were classified into eight groups, four of which included 36 of the 46 samples analysed. Interestingly, all the four tree samples collected in Ueno Park of Tokyo were members of the four groups mentioned above. Based on comparative analysis of their mutations, one of the four trees growing in Ueno Park was concluded to be the closest to the original ancestor. We propose that somatic mutations may be used as tracers to establish the ancestral relationship amongst clonally propagated individuals.

The cryptic complex rearrangements involving the DMD gene: etiologic clues about phenotypical differences revealed by optical genome mapping

BACKGROUND

Deletion or duplication in the DMD gene is one of the most common causes of Duchenne and Becker muscular dystrophy (DMD/BMD). However, the pathogenicity of complex rearrangements involving DMD, especially segmental duplications with unknown breakpoints, is not well understood. This study aimed to evaluate the structure, pattern, and potential impact of rearrangements involving DMD duplication.

METHODS

Two families with DMD segmental duplications exhibiting phenotypical differences were recruited. Optical genome mapping (OGM) was used to explore the cryptic pattern of the rearrangements. Breakpoints were validated using long-range polymerase chain reaction combined with next-generation sequencing and Sanger sequencing.

RESULTS

A multi-copy duplication involving exons 64-79 of DMD was identified in Family A without obvious clinical symptoms. Family B exhibited typical DMD neuromuscular manifestations and presented a duplication involving exons 10-13 of DMD. The rearrangement in Family A involved complex in-cis tandem repeats shown by OGM but retained a complete copy (reading frame) of DMD inferred from breakpoint validation. A reversed insertion with a segmental repeat was identified in Family B by OGM, which was predicted to disrupt the normal structure and reading frame of DMD after confirming the breakpoints.

CONCLUSIONS

Validating breakpoint and rearrangement pattern is crucial for the functional annotation and pathogenic classification of genomic structural variations. OGM provides valuable insights into etiological analysis of DMD/BMD and enhances our understanding for cryptic effects of complex rearrangements.

Implementing next-generation sequencing for diagnosis and management of hereditary hearing impairment: a comprehensive review

INTRODUCTION

Sensorineural hearing impairment (SNHI), a common childhood disorder with heterogeneous genetic causes, can lead to delayed language development and psychosocial problems. Next-generation sequencing (NGS) offers high-throughput screening and high-sensitivity detection of genetic etiologies of SNHI, enabling clinicians to make informed medical decisions, provide tailored treatments, and improve prognostic outcomes.

AREAS COVERED

This review covers the diverse etiologies of HHI and the utility of different NGS modalities (targeted sequencing and whole exome/genome sequencing), and includes HHI-related studies on newborn screening, genetic counseling, prognostic prediction, and personalized treatment. Challenges such as the trade-off between cost and diagnostic yield, detection of structural variants, and exploration of the non-coding genome are also highlighted.

EXPERT OPINION

In the current landscape of NGS-based diagnostics for HHI, there are both challenges (e.g. detection of structural variants and non-coding genome variants) and opportunities (e.g. the emergence of medical artificial intelligence tools). The authors advocate the use of technological advances such as long-read sequencing for structural variant detection, multi-omics analysis for non-coding variant exploration, and medical artificial intelligence for pathogenicity assessment and outcome prediction. By integrating these innovations into clinical practice, precision medicine in the diagnosis and management of HHI can be further improved.

High positive predictive value of CNVs detected by clinical exome sequencing in suspected genetic diseases

BACKGROUND

Genetic disorders often manifest as abnormal fetal or childhood development. Copy number variations (CNVs) represent a significant genetic mechanism underlying such disorders. Despite their importance, the effectiveness of clinical exome sequencing (CES) in detecting CNVs, particularly small ones, remains incompletely understood. We aimed to evaluate the detection of both large and small CNVs using CES in a substantial clinical cohort, including parent-offspring trios and proband only analysis.

METHODS

We conducted a retrospective analysis of CES data from 2428 families, collected from 2018 to 2021. Detected CNV were categorized as large or small, and various validation techniques including chromosome microarray (CMA), Multiplex ligation-dependent probe amplification assay (MLPA), and/or PCR-based methods, were employed for cross-validation.

RESULTS

Our CNV discovery pipeline identified 171 CNV events in 154 cases, resulting in an overall detection rate of 6.3%. Validation was performed on 113 CNVs from 103 cases to assess CES reliability. The overall concordance rate between CES and other validation methods was 88.49% (100/113). Specifically, CES demonstrated complete consistency in detecting large CNV. However, for small CNVs, consistency rates were 81.08% (30/37) for deletions and 73.91% (17/23) for duplications.

CONCLUSION

CES demonstrated high sensitivity and reliability in CNV detection. It emerges as an economical and dependable option for the clinical CNV detection in cases of developmental abnormalities, especially fetal structural abnormalities.

Optical Genome Mapping Identifies a Novel Unbalanced Translocation Between Chromosomes 4q and 6q Leading to Feeding Difficulties and Hypotonia in a Neonate: A Case Report

Optical Genome Mapping (OGM) technology has garnered growing interest for the identification of chromosomal structural variations (SVs), particularly complex ones that are implicated in genetic diseases in humans. In this study, we performed genetic diagnostics on a neonatal patient who presented with feeding difficulties, hypotonia, and an atrial septal defect. We utilized a combination of trio-whole exome sequencing and OGM for our analysis. The results revealed an unbalanced translocation between maternal chromosomes 4 and 6 in the proband, ogm[GRch38]t(4:6)(q35.2;q25.3), resulting in a 2.8 Mb deletion at the 4q35 terminal and a 10.2 Mb duplication at the 6q25 terminal. In summary, this study highlights how OGM, in conjunction with other genetic approaches, can unveil the genetic etiology of complex clinical syndromes. Neonatal patients often exhibit low specific phenotypes, underlining the significance of SV detection.

LoRA-TV: read depth profile-based clustering of tumor cells in single-cell sequencing

Single-cell sequencing has revolutionized our ability to dissect the heterogeneity within tumor populations. In this study, we present LoRA-TV (Low Rank Approximation with Total Variation), a novel method for clustering tumor cells based on the read depth profiles derived from single-cell sequencing data. Traditional analysis pipelines process read depth profiles of each cell individually. By aggregating shared genomic signatures distributed among individual cells using low-rank optimization and robust smoothing, the proposed method enhances clustering performance. Results from analyses of both simulated and real data demonstrate its effectiveness compared with state-of-the-art alternatives, as supported by improvements in the adjusted Rand index and computational efficiency.

Retrospective study revealed integration of CNV-seq and karyotype analysis is an effective strategy for prenatal diagnosis of chromosomal abnormalities

Fetal chromosomal abnormalities are the main cause of adverse pregnancy outcomes and are the focus of invasive prenatal diagnosis. Recent studies have demonstrated that various techniques have distinct advantages. Achieving high-resolution and effective prenatal chromosomal abnormality diagnosis requires a multi-technology integration strategy. Based on retrospective samples from a single center, we propose that integrating CNV-seq and karyotype analysis is an effective strategy for prenatal diagnosis of chromosomal abnormalities. In this study, 13.80% of the pregnant women (347/2514) were found to have likely pathogenic or pathogenic fetal chromosomal abnormalities using this integrated approach. Among these cases, 53.89% (187/347) had consistent chromosomal abnormalities detected by both CNV-seq and karyotyping analysis, while 19.02% (66/347) and 27.09% (94/347) of cases were diagnosed solely by CNV-seq or karyotyping, respectively. Fetal chromosomal abnormalities were identified in 18.39% of samples with abnormal ultrasound, which was significantly higher than the percentage found in samples with normal ultrasound (p < 0.001). Samples with multiple ultrasound abnormalities and single-indicator ultrasound abnormalities such as nasal bone dysplasia, renal dysplasia, or echogenic fetal bowel also had higher rates of chromosomal abnormalities (p < 0.05) compared to normal samples. Analyzing samples with Trio family data (N = 521) revealed that about 94% of variants of uncertain significance were inherited from parents and were non-pathogenic. Overall, integrating CNV-seq and karyotype analysis is an effective strategy for prenatal diagnosis of chromosomal abnormalities. This study provides valuable insights for correlating prenatal screening indicators with chromosomal abnormalities.

A Novel 165 Kb Duplication Involving the α-Globin Gene Cluster Is Identified by Low-Pass Whole Genome Sequencing in a Chinese Thalassemia Intermedia Patient

Copy number variations (CNVs) involving the α-globin gene cluster can lead to an imbalance in the proportion of α- and β-globin chains and consequently cause clinical symptoms of β-thalassemia. In our case, a 6-year-old boy, clinically diagnosed with β thalassemia intermedia, was admitted for further genetic diagnosis with his family. Targeted sequencing and third generation sequencing (TGS) were used to detect the possible variants of the thalassemia genes. Low-pass whole genome sequencing (lpWGS) was conducted to specify the exact location of relevant CNVs across the genome, which was then validated by multiplex ligation-dependent probe amplification.The results revealed that the patient had a heterozygous β0 mutation of Codon17 (A > T) and a full duplication of the α-globin gene cluster, inherited from his mother and father, respectively. Besides, a novel point mutation within the 5' untranslated region of β-Globin (HBB: c. -175 (G > A) was only detected in the patient. This study suggests that lpWGS seems a powerful alternative to detect large CNVs related to thalassemia with second intention for more information of the breakpoints and a simultaneous genome-scale detection of other pathogenic CNVs.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

The reckoning of chromosomal instability: past, present, future

Quantitative measures of CIN are crucial to our understanding of its role in cancer. Technological advances have changed the way CIN is quantified, offering increased accuracy and insight. Here, we review measures of CIN through its rise as a field, discuss considerations for its measurement, and look forward to future quantification of CIN.

On the core segmentation algorithms of copy number variation detection tools

Shotgun sequencing is a high-throughput method used to detect copy number variants (CNVs). Although there are numerous CNV detection tools based on shotgun sequencing, their quality varies significantly, leading to performance discrepancies. Therefore, we conducted a comprehensive analysis of next-generation sequencing-based CNV detection tools over the past decade. Our findings revealed that the majority of mainstream tools employ similar detection rationale: calculates the so-called read depth signal from aligned sequencing reads and then segments the signal by utilizing either circular binary segmentation (CBS) or hidden Markov model (HMM). Hence, we compared the performance of those two core segmentation algorithms in CNV detection, considering varying sequencing depths, segment lengths and complex types of CNVs. To ensure a fair comparison, we designed a parametrical model using mainstream statistical distributions, which allows for pre-excluding bias correction such as guanine-cytosine (GC) content during the preprocessing step. The results indicate the following key points: (1) Under ideal conditions, CBS demonstrates high precision, while HMM exhibits a high recall rate. (2) For practical conditions, HMM is advantageous at lower sequencing depths, while CBS is more competitive in detecting small variant segments compared to HMM. (3) In case involving complex CNVs resembling real sequencing, HMM demonstrates more robustness compared with CBS. (4) When facing large-scale sequencing data, HMM costs less time compared with the CBS, while their memory usage is approximately equal. This can provide an important guidance and reference for researchers to develop new tools for CNV detection.

Molecular underpinnings and environmental drivers of loss of heterozygosity in Drosophila intestinal stem cells

During development and aging, genome mutation leading to loss of heterozygosity (LOH) can uncover recessive phenotypes within tissue compartments. This phenomenon occurs in normal human tissues and is prevalent in pathological genetic conditions and cancers. While studies in yeast have defined DNA repair mechanisms that can promote LOH, the predominant pathways and environmental triggers in somatic tissues of multicellular organisms are not well understood. Here, we investigate mechanisms underlying LOH in intestinal stem cells in Drosophila. Infection with the pathogenic bacteria, Erwinia carotovora carotovora 15, but not Pseudomonas entomophila, increases LOH frequency. Using whole genome sequencing of somatic LOH events, we demonstrate that they arise primarily via mitotic recombination. Molecular features and genetic evidence argue against a break-induced replication mechanism and instead support cross-over via double Holliday junction-based repair. This study provides a mechanistic understanding of mitotic recombination, an important mediator of LOH, and its effects on stem cells in vivo.

Exome and genome sequencing to unravel the precise breakpoints of partial trisomy 6q and partial Monosomy 2q

BACKGROUND

Patients with complex phenotypes and a chromosomal translocation are particularly challenging, since several potentially pathogenic mechanisms need to be investigated.

CASE PRESENTATION

Here, we combined exome and genome sequencing techniques to identify the precise breakpoints of heterozygous microduplications in the 6q25.3-q27 region and microdeletions in the 2q37.1-q37.3 region in a proband. The 5-year-old girl exhibited a severe form of congenital cranial dysinnervation disorder (CCDD) in addition to skeletal dysmorphism anomalies and severe intellectual disability. This is the second case affecting chromosomes 2q and 6q. The individual's karyotype showed an unbalanced translocation 46,XX,del(2)t(2;6)(q37.1;q25.3), which was inherited from her unaffected father [46,XY,t(2;6)(q37.1;q25.3)]. We also obtained the precise breakpoints of a de novo heterozygous copy number deletion [del(2)(q37.1q37.3)chr2:g.232963568_24305260del] and a copy number duplication [dup(6)(q25.3q27)chr6:g.158730978_170930050dup]. The parental origin of the observed balanced translocation was not clear because the parents declined genetic testing.

CONCLUSION

Patients with a 2q37 deletion and 6q25.3 duplication may exhibit severe significant neurological and skeletal dysmorphisms, and the utilization of exome and genome sequencing techniques has the potential to unveil the entire translocation of the CNV and the precise breakpoint.

Association between the copy number variation of CCSER1 gene and growth traits in Chinese Capra hircus (goat) populations

Recently, Coiled-coil serine-rich protein 1 (CCSER1) gene is reported to be related to economic traits in livestock, and become a hotspot. In our study, we detected CCSER1 gene CNV in 693 goats from six breeds (GZB, GZW, AN, BH, HG, TH) by quantitative real-time PCR (qPCR) and the association analysis between the types of CNV and growth traits. Then, CCSER1 gene expression pattern was discovered in seven tissues from NB goats. Our results showed that the CCSER1 gene copy numbers were distributed differently in the aforementioned six breeds. The type of CCSER1 gene CNV was significantly associated with body weight and heart girth traits in GZW goat, in which individuals with deletion type were dominant in body weight trait (P < 0.05), while the normal type individuals were more advantageous in heart girth trait (P < 0.01); and there was a significant association with heart girth in TH goat (P < 0.05), which normal type was the dominant one. The expression profile revealed that CCSER1 gene has the highest level in the lung, followed by the small intestine and heart. In conclusion, our result is dedicated to an in-depth study of the novel CCSER1 gene CNV site and to provide essential information for Chinese goats molecular selective breeding in the future.

Analysis of copy number variations and possible candidate genes in spontaneous abortion by copy number variation sequencing

Introduction

Embryonic chromosomal abnormalities represent a major causative factor in early pregnancy loss, highlighting the importance of understanding their role in spontaneous abortion. This study investigates the potential correlation between chromosomal abnormalities and spontaneous abortion using copy number variation sequencing (CNV-seq), a Next-Generation Sequencing (NGS) technology.

Methods

We analyzed Copy Number Variations (CNVs) in 395 aborted fetal specimens from spontaneous abortion patients by CNV-seq. And collected correlated data, including maternal age, gestational week, and Body Mass Index (BMI), and analyzed their relationship with the CNVs.

Results

Out of the 395 cases, 67.09% of the fetuses had chromosomal abnormalities, including numerical abnormalities, structural abnormalities, and mosaicisms. Maternal age was found to be an important risk factor for fetal chromosomal abnormalities, with the proportion of autosomal trisomy in abnormal karyotypes increasing with maternal age, while polyploidy decreased. The proportion of abnormal karyotypes with mosaic decreased as gestational age increased, while the frequency of polyploidy and sex chromosome monosomy increased. Gene enrichment analysis identified potential miscarriage candidate genes and functions, as well as pathogenic genes and pathways associated with unexplained miscarriage among women aged below or over 35 years old. Based on our study, it can be inferred that there is an association between BMI values and the risk of recurrent miscarriage caused by chromosomal abnormalities.

Discussion

Overall, these findings provide important insights into the understanding of spontaneous abortion and have implications for the development of personalized interventions for patients with abnormal karyotypes.

Genomic Variations Explorer (GenVarX): a toolset for annotating promoter and CNV regions using genotypic and phenotypic differences

The rapid growth of sequencing technology and its increasing popularity in biology-related research over the years has made whole genome re-sequencing (WGRS) data become widely available. A large amount of WGRS data can unlock the knowledge gap between genomics and phenomics through gaining an understanding of the genomic variations that can lead to phenotype changes. These genomic variations are usually comprised of allele and structural changes in DNA, and these changes can affect the regulatory mechanisms causing changes in gene expression and altering the phenotypes of organisms. In this research work, we created the GenVarX toolset, that is backed by transcription factor binding sequence data in promoter regions, the copy number variations data, SNPs and Indels data, and phenotypes data which can potentially provide insights about phenotypic differences and solve compelling questions in plant research. Analytics-wise, we have developed strategies to better utilize the WGRS data and mine the data using efficient data processing scripts, libraries, tools, and frameworks to create the interactive and visualization-enhanced GenVarX toolset that encompasses both promoter regions and copy number variation analysis components. The main capabilities of the GenVarX toolset are to provide easy-to-use interfaces for users to perform queries, visualize data, and interact with the data. Based on different input windows on the user interface, users can provide inputs corresponding to each field and submit the information as a query. The data returned on the results page is usually displayed in a tabular fashion. In addition, interactive figures are also included in the toolset to facilitate the visualization of statistical results or tool outputs. Currently, the GenVarX toolset supports soybean, rice, and Arabidopsis. The researchers can access the soybean GenVarX toolset from SoyKB via https://soykb.org/SoybeanGenVarX/, rice GenVarX toolset, and Arabidopsis GenVarX toolset from KBCommons web portal with links https://kbcommons.org/system/tools/GenVarX/Osativa and https://kbcommons.org/system/tools/GenVarX/Athaliana, respectively.

Validation and depth evaluation of low-pass genome sequencing in prenatal diagnosis using 387 amniotic fluid samples

BACKGROUND

Low-pass genome sequencing (LP GS) is an alternative to chromosomal microarray analysis (CMA). However, validations of LP GS as a prenatal diagnostic test for amniotic fluid are rare. Moreover, sequencing depth of LP GS in prenatal diagnosis has not been evaluated.

OBJECTIVE

The diagnostic performance of LP GS was compared with CMA using 375 amniotic fluid samples. Then, sequencing depth was evaluated by downsampling.

RESULTS

CMA and LP GS had the same diagnostic yield (8.3%, 31/375). LP GS showed all copy number variations (CNVs) detected by CMA and six additional variant of uncertain significance CNVs (>100 kb) in samples with negative CMA results; CNV size influenced LP GS detection sensitivity. CNV detection was greatly influenced by sequencing depth when the CNV size was small or the CNV was located in the azoospermia factor c (AZFc) region of the Y chromosome. Large CNVs were less affected by sequencing depth and more stably detected. There were 155 CNVs detected by LP GS with at least a 50% reciprocal overlap with CNVs detected by CMA. With 25 M uniquely aligned high-quality reads (UAHRs), the detection sensitivity for the 155 CNVs was 99.14%. LP GS using samples with 25 M UAHRs showed the same performance as LP GS using total UAHRs. Considering the detection sensitivity, cost and interpretation workload, 25 M UAHRs are optimal for detecting most aneuploidies and microdeletions/microduplications.

CONCLUSION

LP GS is a promising, robust alternative to CMA in clinical settings. A total of 25 M UAHRs are sufficient for detecting aneuploidies and most microdeletions/microduplications.

Hub gene associated with prognosis in bladder cancer is a novel therapeutic target

Objective

Bladder cancer is a clinical and social conundrum due to its high incidence and recurrence rate. It is urgent to find new targets for the diagnosis and treatment of bladder cancer and improve the prognosis and survival rate of bladder cancer patients. We sought a prognosis-related gene, built related models of evaluated bladder cancer and identified the function of the hub gene in bladder cancer.

Methods

We downloaded the data of bladder cancer patients from the TCGA database, and used differentially expressed genes (DEGs), copy number variation (CNV) and survival analysis to scan the hub genes associated with prognosis in bladder cancer. Then, multi-factor cox regression was used to obtain the bladder cancer prognosis correlation model. Then, we analyzed the relationship between the expression of hub gene and immune microenvironment of bladder cancer. The relationship between the expression of hub gene and prognosis in bladder cancer patients was verified by immunohistochemistry. Cell proliferation assay and drug sensitivity test in vivo were used to verify the inhibition of bladder cancer by targeted inhibitors.

Results

In bladder cancer, we screened seven hub genes (ACLY, CNP, NKIRAS2, P3H4, PDIA6, VPS25 and XPO1) associated with survival. Moreover, the multifactor regression model constructed with hub gene can well distinguish the prognosis of bladder cancer. Hub gene is mostly associated with immune microenvironment. Immunohistochemical results basically confirmed the importance of XPO1 in bladder cancer. Selinexor (an inhibitor of XPO1) could effectively inhibit the proliferation of bladder cancer in the cell proliferation experiments by CCK-8 assays and it could suppress the growth of bladder cancer in mouse bladder cancer model.

Conclusions

In this study, a prognostic model with seven hub genes has provided great help for the prognosis prediction of bladder cancer patients. And XPO1 is an important target affecting the prognosis of bladder cancer, and inhibition of XPO1 can effectively inhibit bladder cancer proliferation and growth.

Sen1 and Rrm3 ensure permissive topological conditions for replication termination

Replication forks terminate at TERs and telomeres. Forks that converge or encounter transcription generate topological stress. Combining genetics, genomics, and transmission electron microscopy, we find that Rrm3^hPif1 and Sen1^hSenataxin helicases assist termination at TERs; Sen1 specifically acts at telomeres. rrm3 and sen1 genetically interact and fail to terminate replication, exhibiting fragility at termination zones (TERs) and telomeres. sen1rrm3 accumulates RNA-DNA hybrids and X-shaped gapped or reversed converging forks at TERs; sen1, but not rrm3, builds up RNA polymerase II (RNPII) at TERs and telomeres. Rrm3 and Sen1 restrain Top1 and Top2 activities, preventing toxic accumulation of positive supercoil at TERs and telomeres. We suggest that Rrm3 and Sen1 coordinate the activities of Top1 and Top2 when forks encounter transcription head on or codirectionally, respectively, thus preventing the slowing down of DNA and RNA polymerases. Hence Rrm3 and Sen1 are indispensable to generate permissive topological conditions for replication termination.

Prenatal diagnosis of hereditary diffuse gastric cancer: a case report

BACKGROUND

Hereditary diffuse gastric cancer(HDGC) is a kind of malignant gastric cancer that is difficult to find in the early stage. However, this late onset and incomplete penetrance hereditary cancer, and its prenatal diagnosis have rarely been reported previously.

CASE PRESENTATION

A 26-year-old woman was referred to genetic counseling for an ultrasonography of fetal choroid plexus cyst at 17 weeks of gestation. The ultrasonographic evaluation showed bilateral choroid plexus cysts(CPC) in the lateral ventricles, and the women showed a family history of gastric cancer and breast cancer. Trio copy number sequencing identified a pathogenic CDH1 deletion in the fetus and unaffected mother. The CDH1 deletion was found in three of the five family members tested, segregation among affected family members. The couple finally decided to terminate the pregnancy after genetic counseling by hospital geneticists due to the uncertainty of the occurrence of HDGC in the future.

CONCLUSIONS

In prenatal diagnosis, a family history of cancer should be widely concerned, and prenatal diagnosis of hereditary tumors requires extensive cooperation between the prenatal diagnosis structure and the pathology department.

Novel maternal duplication of 6p22.3-p25.3 with subtelomeric 6p25.3 deletion: new clinical findings and genotype-phenotype correlations

BACKGROUND

Copy-number variants (CNVs) drive many neurodevelopmental-related disorders. Although many neurodevelopmental-related CNVs can give rise to widespread phenotypes, it is necessary to identify the major genes contributing to phenotypic presentation. Copy-number variations in chromosome 6, such as independent 6p deletion and 6p duplication, have been reported in several live-born infants and present widespread abnormalities such as intellectual disability, growth deficiency, developmental delay, and multiple dysmorphic facial features. However, a contiguous deletion and duplication in chromosome 6p regions have been reported in only a few cases.

CASE PRESENTATION

In this study, we reported the first duplication of chromosome band 6p25.3-p22.3 with deletion of 6p25.3 in a pedigree. This is the first case reported involving CNVs in these chromosomal regions. In this pedigree, we reported a 1-year-old boy with maternal 6p25-pter duplication characterized by chromosome karyotype. Further analysis using CNV-seq revealed a 20.88-Mb duplication at 6p25.3-p22.3 associated with a contiguous 0.66-Mb 6p25.3 deletion. Whole exome sequencing confirmed the deletion/duplication and identified no pathogenic or likely pathogenic variants related with the patient´s phenotype. The proband presented abnormal growth, developmental delay, skeletal dysplasia, hearing loss, and dysmorphic facial features. Additionally, he presented recurrent infection after birth. CNV-seq using the proband´s parental samples showed that the deletion/duplication was inherited from the proband´s mother, who exhibited a similar phenotype to the proband. When compared with other cases, this proband and his mother presented a new clinical finding: forearm bone dysplasia. The major candidate genes contributing to recurrent infection, eye development, hearing loss features, neurodevelopmental development, and congenital bone dysplasia were further discussed.

CONCLUSIONS

Our results showed a new clinical finding of a contiguous deletion and duplication in chromosome 6p regions and suggested candidate genes associated with phenotypic features, such as FOXC1, SERPINB6, NRN1, TUBB2A, IRF4, and RIPK1.

Triplication of the interferon receptor locus contributes to hallmarks of Down syndrome in a mouse model

Down syndrome (DS), the genetic condition caused by trisomy 21, is characterized by variable cognitive impairment, immune dysregulation, dysmorphogenesis and increased prevalence of diverse co-occurring conditions. The mechanisms by which trisomy 21 causes these effects remain largely unknown. We demonstrate that triplication of the interferon receptor (IFNR) gene cluster on chromosome 21 is necessary for multiple phenotypes in a mouse model of DS. Whole-blood transcriptome analysis demonstrated that IFNR overexpression associates with chronic interferon hyperactivity and inflammation in people with DS. To define the contribution of this locus to DS phenotypes, we used genome editing to correct its copy number in a mouse model of DS, which normalized antiviral responses, prevented heart malformations, ameliorated developmental delays, improved cognition and attenuated craniofacial anomalies. Triplication of the Ifnr locus modulates hallmarks of DS in mice, suggesting that trisomy 21 elicits an interferonopathy potentially amenable to therapeutic intervention.

Clinical diagnosis of genetic disorders at both single-nucleotide and chromosomal levels based on BGISEQ-500 platform

Most variations in the human genome refer to single-nucleotide variation (SNV), small fragment insertions and deletions, and genomic copy number variation (CNV). Many human diseases including genetic disorders are associated with variations in the genome. These disorders are often difficult to be diagnosed because of their complex clinical conditions, therefore, an effective detection method is needed to facilitate clinical diagnosis and prevent birth defects. With the development of high-throughput sequencing technology, the method of targeted sequence capture chip has been extensively used owing to its high throughput, high accuracy, fast speed, and low cost. In this study, we designed a chip that potentially captured the coding region of 3043 genes associated with 4013 monogenic diseases, with an addition of 148 chromosomal abnormalities that can be identified by targeting specific regions. To assess the efficiency, a strategy of combining the BGISEQ500 sequencing platform with the designed chip was utilized to screen variants in 63 patients. Eventually, 67 disease-associated variants were found, 31 of which were novel. The results of the evaluation test also show that this combined strategy complies with the requirements of clinical testing and has proper clinical application value.

A comparison of algorithms for identifying copy number variants in family-based whole-exome sequencing data and its implications in inheritance pattern analysis

There remain challenges in accurately identifying constitutional or germline copy number variants (gCNVs) based on whole-exome sequencing data that have implications for genetic diagnosis for 'rare undiagnosed disease' in the clinical setting. Although multiple algorithms have been proposed, a systematic comparison of these algorithms for calling gCNVs and analyzing inherited pattern have yet to be fully conducted. Therefore, we empirically compared seven exome-based algorithms, including XHMM, CLAMMS, CODEX2, ExomeDepth, DECoN, CN.MOPS, and GATK gCNV, for calling gCNVs in 151 individuals from 44 pedigrees, together with the gold standard of genotyping-derived gCNVs in the same cohort for the performance assessment. These algorithms demonstrated varied powers in identifying gCNVs, although the distribution of gCNVs size was similar. The number of shared gCNVs across these algorithms was limited (e.g., only four gCNVs shared among seven algorithms); however, several algorithms showed varying degrees of consistency (e.g., 1,843 gCNVs shared between DECoN and ExomeDepth). CLAMMS and CODEX2 outperformed the remaining algorithms according to a relatively higher F-score (i.e., 0.145 and 0.152, respectively). In addition, these algorithms exhibited different Mendelian inconsistencies of gCNVs and significant challenges remained in inheritance pattern analysis. In conclusion, selecting good algorithms may have important implications in gCNVs-based inheritance pattern analysis for family-based studies.

Additional diagnostic value of CNV-seq over conventional karyotyping in prenatal diagnosis: A systematic review and meta-analysis

OBJECTIVE

To identify the additional diagnostic value of CNV-seq over conventional karyotyping on the part of chromosomal abnormalities in prenatal diagnosis.

METHOD

This was a systematic review conducted in accordance with PRISMA criteria. In order to clarify related research, PubMed, Web of Science databases (including Core Collection, BIOSIS Previews, MEDLINE, and so on), The Cochrane Library and Wiley Online Library were searched with the terms: "prenatal diagnosis," "CNV-seq," "karyotyping," published from January 2010 to May 2022. No language restrictions. RenMan 5.4 was used for the meta-analysis.

RESULTS

Eight studies were included in this systemic review and meta-analysis, including 11 091 pregnant women with high-risk pregnancy factors or with structurally abnormal fetus under ultrasound. CNV-seq detected a 2% (95% CI, -0% to 4%) additional chromosomal anomalies over conventional karyotyping in the six series. A 4% (95% CI, 3%-6%) pooled mean incremental yield of pathogenic CNVs by CNV-seq over karyotyping was observed, with a 1%-16% range.

CONCLUSION

CNV-seq, applied in prenatal diagnosis, may detect more chromosomal abnormalities when compared with karyotyping. With the advantages of wide coverage, high throughput, high resolution, no culture, good compatibility, and adjustable sequencing depth, CNV-seq has high application value in prenatal diagnosis.

Local augmented graph neural network for multi-omics cancer prognosis prediction and analysis

Cancer prognosis prediction and analysis can help patients understand expected life and help clinicians provide correct therapeutic guidance. Thanks to the development of sequencing technology, multi-omics data, and biological networks have been used for cancer prognosis prediction. Besides, graph neural networks can simultaneously consider multi-omics features and molecular interactions in biological networks, becoming mainstream in cancer prognosis prediction and analysis. However, the limited number of neighboring genes in biological networks restricts the accuracy of graph neural networks. To solve this problem, a local augmented graph convolutional network named LAGProg is proposed in this paper for cancer prognosis prediction and analysis. The process follows: first, given a patient's multi-omics data features and biological network, the corresponding augmented conditional variational autoencoder generates features. Then, the generated augmented features and the original features are fed into a cancer prognosis prediction model to complete the cancer prognosis prediction task. The conditional variational autoencoder consists of two parts: encoder-decoder. In the encoding phase, an encoder learns the conditional distribution of the multi-omics data. As a generative model, a decoder takes the conditional distribution and the original feature as inputs to generate the enhanced features. The cancer prognosis prediction model consists of a two-layer graph convolutional neural network and a Cox proportional risk network. The Cox proportional risk network consists of fully connected layers. Extensive experiments on 15 real-world datasets from TCGA demonstrated the effectiveness and efficiency of the proposed method in predicting cancer prognosis. LAGProg improved the C-index values by an average of 8.5% over the state-of-the-art graph neural network method. Moreover, we confirmed that the local augmentation technique could enhance the model's ability to represent multi-omics features, improve the model's robustness to missing multi-omics features, and prevent the model's over-smoothing during training. Finally, based on genes identified through differential expression analysis, we discovered 13 prognostic markers highly associated with breast cancer, among which ten genes have been proved by literature review.

An improved reference genome for Trifolium subterraneum L. provides insight into molecular diversity and intra-specific phylogeny

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

Copy Number Variation Analysis Revealed the Evolutionary Difference between Chinese Indigenous Pigs and Asian Wild Boars

Copy number variation (CNV) has been widely used to study the evolution of different species. We first discovered different CNVs in 24 Anqingliubai pigs and 6 Asian wild boars using next-generation sequencing at the whole-genome level with 10× depth to understand the relationship between genetic evolution and production traits in wild boars and domestic pigs. A total of 97,489 CNVs were identified and divided into 10,429 copy number variation regions (CNVRs), occupying 32.06% of the porcine genome. Chromosome 1 had the most CNVRs, and chromosome 18 had the least. Ninety-six CNVRs were selected using VST 1% based on the signatures of all CNVRs, and sixty-five genes were identified in the selected regions. These genes were strongly correlated with traits distinguishing groups by enrichment in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways, such as growth (CD36), reproduction (CIT, RLN), detoxification (CYP3A29), and fatty acid metabolism (ELOVL6). The QTL overlapping regions were associated with meat traits, growth, and immunity, which was consistent with CNV analysis. Our findings increase the understanding of evolved genome structural variations between wild boars and domestic pigs, and provide new molecular biomarkers to guide breeding and the efficient use of available genetic resources.

Pain-Related Gene Solute Carrier Family 24 Member 3 Is a Prognostic Biomarker and Correlated with Immune Infiltrates in Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma: A Study via Integrated Bioinformatics Analyses and Experimental Verification

The aim of this study was to explore cervical carcinoma and screen a suitable gene as the biomarker used for prognosis evaluation as well as pain therapy. Low expression levels of solute carrier family 24 member 3 (SLC24A3) was involved in the appearance and development of numerous malignancies. Nevertheless, the prognostic value of SLC24A3 expression with cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) patients remains uncertain. During the present study, SLC24A3 expression in CESC was retrieved from TCGA, GEO, and MSigDB databases. Based on TCGA and GEO profiles, we performed survival and difference analyses about SLC24A3 both in two GEO (GSE44001 and GSE63514) and TCGA-CESC cohorts (all p < 0.05), indicating that SLC24A3 was low expressed in tumors and associated with higher overall survival in CESC patients. Additionally, we programmed a series of analyses, including genomic profiling, enrichment analysis, immune infiltration analysis, and therapy-related analysis to identify the mechanism of the SLC24A3 in the process of cancer in CESC. Meanwhile, qRT-PCR was used to validate that the expression of SLC24A3 mRNA in Hela and SiHa cell lines was significantly lower than in PANC-1 and HUCEC cell lines. Our finding elucidated that the SLC24A3, a sodium-calcium regulator of cells, is an indispensable factor which can significantly influence the prognosis of patients with CESC and could provide novel clinical evidence to serve as a potential biological indicator for future diagnosis and pain therapy.

Evaluation and Analysis of Absence of Homozygosity (AOH) Using Chromosome Analysis by Medium Coverage Whole Genome Sequencing (CMA-seq) in Prenatal Diagnosis

OBJECTIVE

Absence of homozygosity (AOH) is a genetic characteristic known to cause human diseases mainly through autosomal recessive or imprinting mechanisms. The importance and necessity of accurate AOH detection has become more clinically significant in recent years. However, it remains a challenging task for sequencing-based methods thus far.

METHODS

In this study, we developed and optimized a new bioinformatic algorithm based on the assessment of minimum sequencing coverage, optimal bin size, the Z-score threshold of four types of allele count and the frequency for accurate genotyping using 28 AOH negative samples, and redefined the AOH detection cutoff value. We showed the performance of chromosome analysis by five-fold coverage whole genome sequencing (CMA-seq) for AOH identification in 27 typical prenatal/postnatal AOH positive samples, which were previously confirmed by chromosomal microarray analysis with single nucleotide polymorphism array (CMA/SNP array).

RESULTS

The blinded study indicated that for all three forms of AOH, including whole genomic AOH, single chromosomal AOH and segmental AOH, and all kinds of sample types, including chorionic villus sampling, amniotic fluid, cord blood, peripheral blood and abortive tissue, CMA-seq showed equivalent detection power to that of routine CMA/SNP arrays (750K). The subtle difference between the two methods is that CMA-seq is prone to detect small inconsecutive AOHs, while CMA/SNP array reports it as a whole.

CONCLUSION

Based on our newly developed bioinformatic algorithm, it is feasible to detect clinically significant AOH using CMA-seq in prenatal diagnosis.

Comparison of the combined use of CNV-seq and karyotyping or QF-PCR in prenatal diagnosis: a retrospective study

To elevate the accuracy of diagnostic results, CNV-seq is usually performed simultaneously with karyotyping or QF-PCR. Although several studies have investigated the performance of the combined use of CNV-seq with karyotyping or QF-PCR, there have been no reports focusing on the comparison of these 2 diagnostic strategies. In our study, 2507 pregnant women were included to investigate these 2 strategies. The detection rates of foetal genetic abnormalities and turnaround time were compared between these 2 groups. Moreover, the detection rates of foetal genetic abnormalities in different indications were analyzed. Our results unveiled that the detection rates of numerical chromosomal abnormalities were nearly the same in these 2 groups. In addition to numerical chromosomal abnormalities, 39 balanced karyotypic changes and chromosome polymorphisms were detected via the combined use of karyotyping and CNV-seq. Further investigation revealed that the vast majority of these karyotypic changes were inherited from parents. Compared with the karyotyping group, the combination of QF-PCR and CNV-seq reduced the reporting time from 31.593 ± 4.944 days to 11.460 ± 4.894 days. Meanwhile, NIPT, maternal serum screening and ultrasound scan significantly improved the detection of foetal genetic abnormalities. In conclusion, our results revealed that parental karyotyping is a useful supplementary method for CNV-seq and systematic prenatal examinations improved the detection of foetal genetic defects.

Combined use of karyotyping and copy number variation sequencing technology in prenatal diagnosis

Background

Karyotyping and genome copy number variation sequencing (CNV-seq) are two techniques frequently used in prenatal diagnosis. This study aimed to explore the diagnostic potential of using a combination of these two methods in order to provide a more accurate clinical basis for prenatal diagnosis.

Methods

We selected 822 pregnant women undergoing amniocentesis and separated them into six groups according to different risk indicators. Karyotyping and CNV-seq were performed simultaneously to compare the diagnostic performance of the two methods.

Results

Among the different amniocentesis indicators, abnormal fetal ultrasounds accounted for 39.29% of the total number of examinees and made up the largest group. The abnormal detection rate of non-invasive prenatal testing (NIPT) high risk was 37.93% and significantly higher than the other five groups (P < 0.05). The abnormal detection rate of mixed indicators was significantly higher than the history of the adverse reproductive outcomes group (P = 0.0151). The two methods combined found a total of 119 abnormal cases (14.48%). Karyotyping detected 57 cases (6.93%) of abnormal karyotypes, 30 numerical aberrations, and 27 structural aberrations. CNV-seq identified 99 cases (12.04%) with altered CNVs, 30 cases of chromosome aneuploidies, and 69 structural aberrations (28 pathogenic, eight that were likely pathogenic, and 33 microdeletion/duplication variants of uncertain significance (VUS)). Thirty-seven cases were found abnormal by both methods, 20 cases were detected abnormally by karyotyping (mainly mutual translocation and mostly balanced), and 62 cases of microdeletion/duplication were detected by CNV-seq. Steroid sulfatase gene (STS) deletion was identified at chromosome Xp22.31 in three cases. Postnatal follow-up confirmed that babies manifested skin abnormalities one week after birth. Six fetuses had Xp22.31 duplications ranging from 1.5 Kb to 1.7 Mb that were detected by CNV-seq. Follow-up showed that five babies presented no abnormalities during follow-up, except for one terminated pregnancy due to a history of adverse reproductive outcomes.

Conclusion

The combination of using CNV-seq and karyotyping significantly improved the detection rate of fetal pathogenic chromosomal abnormalities. CNV-seq is an effective complement to karyotyping and improves the accuracy of prenatal diagnosis.

Copy Number Variations in the MICALL2 and MOGAT2 Genes Are Associated with Ashidan Yak Growth Traits

Copy number variations (CNVs) are a result of genomic rearrangement affecting DNA regions over 1 kb in length, and can include inversions, translocations, deletions, and duplications. The molecule interacting with CasL-like protein 2 (MICALL2) gene is primarily associated with mitochondrial protein targeting and exhibits predicted stress fiber colocalization. The monoacylglycerol O-acyltransferase 2 (MOGAT2) gene encodes an enzyme responsible for catalyzing diacylglycerol synthesis from 2-monoacylglycerol and fatty acyl-CoA. For this study, blood samples were obtained from 315 yaks, and the body weight, body length, withers height, and chest girth of these animals were measured at 6, 12, 18, and 30 months of age. Genomic DNA was harvested from the collected blood samples, and CNVs in these samples were detected by qPCR. The resultant data were compared using ANOVAs, revealing significant associations between MICALL2 gene CNVs and body weight at 6 months of age (p < 0.05), body length and chest girth at 30 months of age (p < 0.05), and withers height at 18 months of age (p < 0.01) in Ashidan yaks. Similarly, MOGAT2 CNVs were significantly associated with body weight at 6 and 30 months of age (p < 0.05), and with withers height at 18 months of age (p < 0.01) in these Ashidan yaks. MICALL2 and MOGAT2 gene expression was further analyzed in yak tissue samples, revealing that MICALL2 was most highly expressed in the adipose tissue, whereas MOGAT2 was most highly expressed in the lung. These results thus confirmed the relationship between CNVs in the MICALL2 and MOGAT2 genes and Ashidan yak growth traits, providing a valuable gene locus that can be leveraged for future marker-assisted yak breeding efforts.

Early-onset ophthalmoplegia, cervical dyskinesia, and lower extremity weakness due to partial deletion of chromosome 16: A case report

BACKGROUND

We explored the genotype-phenotype correlation of the novel deletion 16p13.2p12.3 in an 8-year-old child with progressive total ophthalmoplegia, cervical dyskinesia, and lower limb weakness by comparing the patient’s clinical features with previously reported data on adjacent copy number variation (CNV) regions.

CASE SUMMARY

Specifically, we first performed whole-exome sequencing, CNV-sequencing, and mitochondrial genome sequencing on the patient and his parents, then applied “MitoExome” (the entire mitochondrial genome and exons of nuclear genes encoding the mitochondrial proteome) analysis to screen for genetic mitochondrial diseases. We identified a de novo 7.23 Mb deletion, covering 16p13.2p12.3, by both whole-exome sequencing and CNV sequencing. We also detected 16p13.11 in the deleted region, which is the recurrent distinct region associated with neurodevelopmental disorder. However, the patient only displayed features of progressive total ophthalmoplegia, cervical dyskinesia, and weakness in his lower limbs without neurodevelopmental disorder. The “MitoExome” sequencing was negative. Brain magnetic resonance imaging revealed non-specific sporadic changes in the occipital parietal lobe and basal ganglia.

CONCLUSION

Taken together, these results indicated that 16p13.2p12.3 deletion causes a syndrome with the phenotype of early-onset total ophthalmoplegia. The “MitoExome” analysis is powerful for the differential diagnosis of mitochondrial diseases. We report a novel copy number variant in this case, but further confirmation is required.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

Origins, genomic structure and copy number variation of snake venom myotoxins

Crotamine, myotoxin a and homologs are short peptides that often comprise major fractions of rattlesnake venoms and have been extensively studied for their bioactive properties. These toxins are thought to be important for rapidly immobilizing mammalian prey and are implicated in serious, and sometimes fatal, responses to envenomation in humans. While high quality reference genomes for multiple venomous snakes are available, the loci that encode myotoxins have not been successfully assembled in any existing genome assembly. Here, we integrate new and existing genomic and transcriptomic data from the Prairie Rattlesnake (Crotalus viridis viridis) to reconstruct, characterize, and infer the chromosomal locations of myotoxin-encoding loci. We integrate long-read transcriptomics (Pacific Bioscience's Iso-Seq) and short-read RNA-seq to infer gene sequence diversity and characterize patterns of myotoxin and paralogous β-defensin expression across multiple tissues. We also identify two long non-coding RNA sequences which both encode functional myotoxins, demonstrating a newly discovered source of venom coding sequence diversity. We also integrate long-range mate-pair chromatin contact data and linked-read sequencing to infer the structure and chromosomal locations of the three myotoxin-like loci. Further, we conclude that the venom-associated myotoxin is located on chromosome 1 and is adjacent to non-venom paralogs. Consistent with this locus contributing to venom composition, we find evidence that the promoter of this gene is selectively open in venom gland tissue and contains transcription factor binding sites implicated in broad trans-regulatory pathways that regulate snake venoms. This study provides the best genomic reconstruction of myotoxin loci to date and raises questions about the physiological roles and interplay between myotoxin and related genes, as well as the genomic origins of snake venom variation.

NCAPG2 Is a Novel Prognostic Biomarker and Promotes Cancer Stem Cell Maintenance in Low-Grade Glioma

Gliomas account for 75% of all primary malignant brain tumors in adults and are associated with high mortality. Mounting evidence has shown that NCAPG2 is differentially expressed in various cancers. However, the prognostic value and immune functions of NCAPG2 in low-grade glioma (LGG) remain unresolved. In the present study, we revealed that NCAPG2 was up-regulated in LGG, and its higher expression was associated with adverse clinical outcomes and poor clinical characteristics, including WHO grade, IDH mutation, 1p/19q codeletion, and primary therapy outcome. The results of the Cox regression analysis revealed that NCAPG2 was an independent factor for the prognosis of low-grade glioma. Meanwhile, we also established a nomogram based on NCAPG2 to predict the 1-, 3-, or 5-year survival in LGG patients. Furthermore, we found that Copy number variation (CNV) and DNA hypomethylation results in its overexpression in LGG. In addition, functional annotation confirmed that NCAPG2 was mainly involved in the immune regulation and WNT signaling pathways. Finally, we determined that increased expression of NCAPG2 was correlated with infiltration levels of various immune cells and immune checkpoint in LGG. Importantly, we found that NCAPG2 was highly expressed in glioma stem cells lines and knockdown of NCAPG2 significantly inhibited the self-renewal ability of GSC. This is the first study to identify NCAPG2 as a new potential prognostic biomarker and characterize the functional roles of NCAPG2 in the progression of LGG, and provides a novel potential diagnostic and therapeutic biomarker for LGG in the future.

Experience of Low-Pass Whole-Genome Sequencing-Based Copy Number Variant Analysis: A Survey of Chinese Tertiary Hospitals

In China, low-pass whole-genome sequencing (low-pass WGS) is emerging as an alternative diagnostic test to detect copy number variants (CNVs). This survey aimed to study the laboratory practice, service quality, and case volumes of low-pass WGS-based CNV analysis among national accredited Chinese tertiary hospitals that have routinely applied low-pass WGS for more than a year and that have been certified in next-generation sequencing (NGS) clinical applications for more than three years. The questionnaire focused on (1) the composition of patients’ referral indications for testing and annual case volumes; (2) the capacity of conducting laboratory assays, bioinformatic analyses, and reporting; (3) the sequencing platforms and parameters utilized; and (4) CNV nomenclature in reports. Participants were required to respond based on their routine laboratory practices and data audited in a 12-month period from February 2019 to January 2020. Overall, 24 participants representing 24 tertiary referral hospitals from 21 provincial administrative regions in China returned the questionnaires. Excluding three hospitals routinely applying low-pass WGS for non-invasive prenatal testing (NIPT) only, the analysis only focused on the data submitted by the rest 21 hospitals. These hospitals applied low-pass WGS-based CNV analysis for four primary applications: high-risk pregnancies, spontaneous abortions, couples with adverse pregnancy history, and children with congenital birth defects. The overall estimated annual sample volume was over 36,000 cases. The survey results showed that the most commonly reported detection limit for CNV size (resolution) was 100 kb; however, the sequencing methods utilized by the participants were variable (single-end: 61.90%, 13/21; paired-end: 28.57%, 6/21; both: 9.52%, 2/21). The diversity was also reflected in the sequencing parameters: the mean read count was 13.75 million reads/case (95% CI, 9.91–17.60) and the read-length median was 65 bp (95% CI, 75.17–104.83). To assess further the compliance of the CNV reporting nomenclature according to the 2016 edition of International System for Human Cytogenomics Nomenclature (ISCN 2016), a scoring metric was applied and yielded responses from 19 hospitals; the mean compliance score was 7.79 out of 10 points (95% CI, 6.78–8.80). Our results indicated that the low-pass WGS-based CNV analysis service is in great demand in China. From a quality control perspective, challenges remain regarding the establishment of standard criteria for low-pass WGS-based CNV analysis and data reporting formats. In summary, the low-pass WGS-based method is becoming a common diagnostic approach, transforming the possibilities for genetic diagnoses for patients in China.

OBIF: an omics-based interaction framework to reveal molecular drivers of synergy

Bioactive molecule library screening may empirically identify effective combination therapies, but molecular mechanisms underlying favorable drug–drug interactions often remain unclear, precluding further rational design. In the absence of an accepted systems theory to interrogate synergistic responses, we introduce Omics-Based Interaction Framework (OBIF) to reveal molecular drivers of synergy through integration of statistical and biological interactions in synergistic biological responses. OBIF performs full factorial analysis of feature expression data from single versus dual exposures to identify molecular clusters that reveal synergy-mediating pathways, functions and regulators. As a practical demonstration, OBIF analyzed transcriptomic and proteomic data of a dyad of immunostimulatory molecules that induces synergistic protection against influenza A and revealed unanticipated NF-κB/AP-1 cooperation that is required for antiviral protection. To demonstrate generalizability, OBIF analyzed data from a diverse array of Omics platforms and experimental conditions, successfully identifying the molecular clusters driving their synergistic responses. Hence, unlike existing synergy quantification and prediction methods, OBIF is a phenotype-driven systems model that supports multiplatform interrogation of synergy mechanisms.

Zero-inflated poisson models with measurement error in the response

Zero-inflated count data arise frequently from genomics studies. Analysis of such data is often based on a mixture model which facilitates excess zeros in combination with a Poisson distribution, and various inference methods have been proposed under such a model. Those analysis procedures, however, are challenged by the presence of measurement error in count responses. In this article, we propose a new measurement error model to describe error-contaminated count data. We show that ignoring the measurement error effects in the analysis may generally lead to invalid inference results, and meanwhile, we identify situations where ignoring measurement error can still yield consistent estimators. Furthermore, we propose a Bayesian method to address the effects of measurement error under the zero-inflated Poisson model and discuss the identifiability issues. We develop a data-augmentation algorithm that is easy to implement. Simulation studies are conducted to evaluate the performance of the proposed method. We apply our method to analyze the data arising from a prostate adenocarcinoma genomics study. This article is protected by copyright. All rights reserved.

Determining the Accuracy of Next Generation Sequencing Based Copy Number Variation Analysis in Hereditary Breast and Ovarian Cancer

BACKGROUND

Copy number variations (CNVs) are commonly associated with malignancies, including hereditary breast and ovarian cancers. Next generation sequencing (NGS) provides solutions for CNV detection in a single run. This study aimed to compare the accuracy of CNV detection by NGS analysing tool against Multiplex Ligation Dependent Probe Amplification (MLPA).

RESEARCH DESIGN AND METHODS

In total, 1276 cases were studied by targeted NGS panels and 691 cases (61 calls in 58 NGS-CNV positive and 633 NGS-CNV negative cases) were validated by MLPA.

RESULTS

Twenty-eight (46%) NGS-CNV positive calls were consistent, whereas 33 (54%) calls showed discordance with MLPA. Two cases were detected as SNV by the NGS and CNV by the MLPA analysis. In total, 2% of the cases showed an MLPA confirmed CNV region in BRCA1/2. The results of this study showed that despite the high false positive call rate of the NGS-CNV algorithm, there were no false negative calls. The cases that were determined to be negative by the NGS and positive by the MLPA were actually carrying SNVs that were located on the MLPA probe binding sites.

CONCLUSION

The diagnostic performance of NGS-CNV analysis is promising; however, the need for confirmation by different methods remains.

Progress in Methods for Copy Number Variation Profiling

Copy number variations (CNVs) are the predominant class of structural genomic variations involved in the processes of evolutionary adaptation, genomic disorders, and disease progression. Compared with single-nucleotide variants, there have been challenges associated with the detection of CNVs owing to their diverse sizes. However, the field has seen significant progress in the past 20–30 years. This has been made possible due to the rapid development of molecular diagnostic methods which ensure a more detailed view of the genome structure, further complemented by recent advances in computational methods. Here, we review the major approaches that have been used to routinely detect CNVs, ranging from cytogenetics to the latest sequencing technologies, and then cover their specific features.

Clinical efficiency of simultaneous CNV-seq and whole-exome sequencing for testing fetal structural anomalies

Background

Birth defects are responsible for approximately 7% of neonatal deaths worldwide by World Health Organization in 2004. Many methods have been utilized for examining the congenital anomalies in fetuses. This study aims to investigate the efficiency of simultaneous CNV-seq and whole-exome sequencing (WES) in the diagnosis of fetal anomaly based on a large Chinese cohort.

Methods

In this cohort study, 1800 pregnant women with singleton fetus in Hubei Province were recruited from 2018 to 2020 for prenatal ultrasonic screening. Those with fetal structural anomalies were transferred to the Maternal and Child Health Hospital of Hubei Province through a referral network in Hubei, China. After multidisciplinary consultation and decision on fetal outcome, products of conception (POC) samples were obtained. Simultaneous CNV-seq and WES was conducted to identify the fetal anomalies that can compress initial DNA and turnaround time of reports.

Results

In total, 959 couples were finally eligible for the enrollment. A total of 227 trios were identified with a causative alteration (CNV or variant), among which 191 (84.14%) were de novo. Double diagnosis of pathogenic CNVs and variants have been identified in 10 fetuses. The diagnostic yield of multisystem anomalies was significantly higher than single system anomalies (32.28% vs. 22.36%, P  = 0.0183). The diagnostic rate of fetuses with consistent intra- and extra-uterine phenotypes (172/684) was significantly higher than the rate of these with inconsistent phenotypes (17/116, P  = 0.0130).

Conclusions

Simultaneous CNV-seq and WES analysis contributed to fetal anomaly diagnosis and played a vital role in elucidating complex anomalies with compound causes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03202-9.

High resolution copy number inference in cancer using short-molecule nanopore sequencing

Genome copy number is an important source of genetic variation in health and disease. In cancer, Copy Number Alterations (CNAs) can be inferred from short-read sequencing data, enabling genomics-based precision oncology. Emerging Nanopore sequencing technologies offer the potential for broader clinical utility, for example in smaller hospitals, due to lower instrument cost, higher portability, and ease of use. Nonetheless, Nanopore sequencing devices are limited in the number of retrievable sequencing reads/molecules compared to short-read sequencing platforms, limiting CNA inference accuracy. To address this limitation, we targeted the sequencing of short-length DNA molecules loaded at optimized concentration in an effort to increase sequence read/molecule yield from a single nanopore run. We show that short-molecule nanopore sequencing reproducibly returns high read counts and allows high quality CNA inference. We demonstrate the clinical relevance of this approach by accurately inferring CNAs in acute myeloid leukemia samples. The data shows that, compared to traditional approaches such as chromosome analysis/cytogenetics, short molecule nanopore sequencing returns more sensitive, accurate copy number information in a cost effective and expeditious manner, including for multiplex samples. Our results provide a framework for short-molecule nanopore sequencing with applications in research and medicine, which includes but is not limited to, CNAs.

Current status of structural variation studies in plants

Structural variations (SVs) including gene presence/absence variations and copy number variations are a common feature of genomes in plants and, together with single nucleotide polymorphisms and epigenetic differences, are responsible for the heritable phenotypic diversity observed within and between species. Understanding the contribution of SVs to plant phenotypic variation is important for plant breeders to assist in producing improved varieties. The low resolution of early genetic technologies and inefficient methods have previously limited our understanding of SVs in plants. However, with the rapid expansion in genomic technologies, it is possible to assess SVs with an ever-greater resolution and accuracy. Here, we review the current status of SV studies in plants, examine the roles that SVs play in phenotypic traits, compare current technologies and assess future challenges for SV studies.

High resolution estimates of relative gene abundance with quantitative ratiometric regression PCR (qRR-PCR)

Quantification of the relative abundance of genetic traits has broad applications for biomarker discovery, diagnostics, and assessing gene expression in biological research. Relative quantification of genes is traditionally done with the 2-ΔΔCT method using quantitative real-time polymerase chain reaction (qPCR) data, which is often limited in resolution beyond orders of magnitude difference. The latest techniques for quantification of nucleic acids employ digital PCR or microarrays which involve lengthy sample preparation and complex instrumentation. In this work, we describe a quantitative ratiometric regression PCR (qRR-PCR) method for computing relative abundance of genetic traits in a sample with high resolution from a single duplexed real-time quantitative PCR assay. Instead of comparing the individual cycle threshold (Ct) values as is done for the 2-ΔΔCT method, our qRR-PCR algorithm leverages the innate relationship of co-amplified PCR targets to measure their relative quantities using characteristic curves derived from the normalized ratios of qPCR fluorescence curves. We demonstrate the utility of this technique for discriminating the fractional abundance of mixed alleles with resolution below 5%.

A chromosome-scale draft genome sequence of horsegram (Macrotyloma uniflorum)

Horsegram (Macrotyloma uniflorum [Lam.] Verdc.) is an underutilized warm-season diploid legume (2n = 20, 22). Because of its ability to grow under water-deficient and marginal soil conditions, horsegram is a preferred choice in the era of global climate change. In recognition of its potential as a crop species, we generated and analyzed a draft genome sequence for a horsegram variety, HPK-4. Ten chromosome-scale pseudomolecules were created by aligning Illumina scaffold sequences onto a linkage map. The total length of the ten pseudomolecules was 259.2 Mbp, covering 89% of the total length of the assembled sequences. A total of 36,105 genes were predicted on the assembled sequences. Diversity analysis of 89 horsegram accessions by dd-RAD-Seq identified 277 single nucleotide polymorphisms (SNPs), suggesting narrow genetic diversity among the horsegram accessions. This is the first attempt to generate a draft genome sequence of horsegram and will provide a reference for sequence-based analysis of horsegram germplasm.