ISMI-VAE: A deep learning model for classifying disease cells using gene expression and SNV data

Various studies have linked several diseases, including cancer and COVID-19, to single nucleotide variations (SNV). Although single-cell RNA sequencing (scRNA-seq) technology can provide SNV and gene expression data, few studies have integrated and analyzed these multimodal data. To address this issue, we introduce Interpretable Single-cell Multimodal Data Integration Based on Variational Autoencoder (ISMI-VAE). ISMI-VAE leverages latent variable models that utilize the characteristics of SNV and gene expression data to overcome high noise levels and uses deep learning techniques to integrate multimodal information, map them to a low-dimensional space, and classify disease cells. Moreover, ISMI-VAE introduces an attention mechanism to reflect feature importance and analyze genetic features that could potentially cause disease. Experimental results on three cancer data sets and one COVID-19 data set demonstrate that ISMI-VAE surpasses the baseline method in terms of both effectiveness and interpretability and can effectively identify disease-causing gene features.

A Predictive Noninvasive Single-Nucleotide Variation-Based Biomarker Signature for Resectable Pancreatic Cancer: Protocol for a Prospective Validation Study

BACKGROUND

Single-nucleotide variations (SNVs; formerly SNPs) are inherited genetic variants that can be easily determined in routine clinical practice using a simple blood or saliva test. SNVs have potential to serve as noninvasive biomarkers for predicting cancer-specific patient outcomes after resection of pancreatic ductal adenocarcinoma (PDAC). Two recent analyses led to the identification and validation of three SNVs in the CD44 and CHI3L2 genes (rs187115, rs353630, and rs684559), which can be used as predictive biomarkers to help select patients most likely to benefit from pancreatic resection. These variants were associated with an over 2-fold increased risk for tumor-related death in three independent PDAC study cohorts from Europe and the United States, including The Cancer Genome Atlas cohorts (reaching a P value of 1×10-8). However, these analyses were limited by the inherent biases of a retrospective study design, such as selection and publication biases, thereby limiting the clinical use of these promising biomarkers in guiding PDAC therapy.

OBJECTIVE

To overcome the limitations of previous retrospectively designed studies and translate the findings into clinical practice, we aim to validate the association of the identified SNVs with survival in a controlled setting using a prospective cohort of patients with PDAC following pancreatic resection.

METHODS

All patients with PDAC who will undergo pancreatic resection at three participating hospitals in Switzerland and fulfill the inclusion criteria will be included in the study consecutively. The SNV genotypes will be determined using standard genotyping techniques from patient blood samples. For each genotyped locus, log-rank and Cox multivariate regression tests will be performed, accounting for the relevant covariates American Joint Committee on Cancer stage and resection status. Clinical follow-up data will be collected for at least 3 years. Sample size calculation resulted in a required sample of 150 patients to sufficiently power the analysis.

RESULTS

The follow-up data collection started in August 2019 and the estimated end of data collection will be in May 2027. The study is still recruiting participants and 142 patients have been recruited as of November 2023. The DNA extraction and genotyping of the SNVs will be performed after inclusion of the last patient. Since no SNV genotypes have been determined, no data analysis has been performed to date. The results are expected to be published in 2027.

CONCLUSIONS

This is the first prospective study of the CD44 and CHI3L2 SNV-based biomarker signature in PDAC. A prospective validation of this signature would enable its clinical use as a noninvasive predictive biomarker of survival after pancreatic resection that is readily available at the time of diagnosis and can assist in guiding PDAC therapy. The results of this study may help to individualize treatment decisions and potentially improve patient outcomes.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/54042.

A Protective Role of Low Polygenic Risk Score in Healthy Individuals Carrying Attention-Deficit/Hyperactivity Disorder-Associated Copy Number Variations

BACKGROUND

Previous studies have implicated both rare copy number variations (CNVs) and common variants in liability for attention-deficit/hyperactivity disorder (ADHD). However, how common and rare genetic variants jointly contribute to individual liability requires further investigation in larger cohorts.

METHODS

This study comprises 9385 participants of European descent and 7810 participants of African American ancestry who were recruited from the greater Philadelphia area by the Children's Hospital of Philadelphia. The polygenic risk score (PRS) of each participant was estimated by linkage disequilibrium pruning and p-value thresholding (P + T) methods using PRSice-2. We investigated whether the risk of ADHD follows a polygenic liability threshold model wherein 1) the risk of ADHD requires less contribution from common variants in the presence of a rare CNV, and 2) control carriers of ADHD-associated CNVs have lower common risk allele burden than noncarriers.

RESULTS

CNVs previously reported in ADHD cases were significantly associated with ADHD risk in both the European American cohort and the African American cohort. Healthy control participants carrying those same risk CNVs had lower PRSs than those without risk CNVs in the European American cohort. This result was replicated in the African American cohort. However, PRSs were not significantly different in case participants carrying risk CNVs versus those without risk CNVs.

CONCLUSIONS

These findings provide evidence in support of interactive effects of PRS and ADHD-associated CNVs on disease risk and add novel insights into the genetic basis of ADHD by highlighting a protective role of low PRS in ADHD.

Changes in intrahost genetic diversity according to lesion severity in longitudinal HPV16 samples

Human papillomavirus type 16 (HPV16) is the most common cause of cervical cancer, but most infections are transient with lesions not progressing to cancer. There is a lack of specific biomarkers for early cancer risk stratification. This study aimed to explore the intrahost HPV16 genomic variation in longitudinal samples from HPV16-infected women with different cervical lesion severity (normal, low-grade, and high-grade). The TaME-seq deep sequencing protocol was used to generate whole genome HPV16 sequences of 102 samples collected over time from 40 individuals. Single nucleotide variants (SNVs) and intrahost SNVs (iSNVs) were identified in the viral genomes. A majority of individuals had a unique set of SNVs and these SNVs were stable over time. Overall, the number of iSNVs and APOBEC3-induced iSNVs were significantly lower in high-grade relative to normal and low-grade samples. A significant increase in the number of APOBEC3-induced iSNVs over time was observed for normal samples when compared to high-grade. Our results indicates that the lower incidence of iSNVs and APOBEC3-induced iSNVs in high-grade lesions may have implications for novel biomarkers discoveries, potentially aiding early stratification of HPV-induced cervical precancerous lesions.

ACE Gene Mutations (rs577350502) in Early-Onset and Recurrent Myocardial Infarction: A Case Report and Review

Background

Acute myocardial infarction (AMI) is a severe acute coronary syndrome, demonstrating a trend toward affecting younger individuals in recent years. The association between early-onset myocardial infarction and single nucleotide polymorphism necessitates further exploration and evaluation.

Case description

We present a case of a patient experiencing early-onset and recurrent myocardial infarction. The patient underwent stent implantation for myocardial infarction at the age of 53 and subsequently encountered two more myocardial infarctions within a span of 16 years. Following interventional therapy, genetic testing was conducted to assess the efficacy of subsequent anti-heart failure medications, with the aim to preemptively address heart failure risks. Genetic testing revealed a mutation in the angiotensin-converting enzyme (ACE) gene (rs577350502, g.63488533C>A), characterized by an intron-deletion single nucleotide variant.

Conclusion

While this variant has not been previously reported to be associated with any specific disease, we hypothesize that it may contribute to the susceptibility and risk of myocardial infarction and coronary heart disease in the patient under consideration. This observation underscores the significance of investigating the insertion/deletion polymorphisms of the ACE gene in the context of AMI and emphasizes the necessity for further validation of this variant and other genetic markers associated with AMI in related diseases.

Single nucleotide variants in the CCL2, OAS1 and DPP9 genes and their association with the severity of COVID-19 in an Ecuadorian population

COVID-19 has a broad clinical spectrum, ranging from asymptomatic-mild form to severe phenotype. The severity of COVID-19 is a complex trait influenced by various genetic and environmental factors. Ethnic differences have been observed in relation to COVID-19 severity during the pandemic. It is currently unknown whether genetic variations may contribute to the increased risk of severity observed in Latin-American individuals The aim of this study is to investigate the potential correlation between gene variants at CCL2, OAS1, and DPP9 genes and the severity of COVID-19 in a population from Quito, Ecuador. This observational case-control study was conducted at the Carrera de Biologia from the Universidad Central del Ecuador and the Hospital Quito Sur of the Instituto Ecuatoriano de Seguridad Social (Quito-SUR-IESS), Quito, Ecuador. Genotyping for gene variants at rs1024611 (A>G), rs10774671 (A>G), and rs10406145 (G>C) of CCL2, OAS1, and DPP9 genes was performed on 100 COVID-19 patients (43 with severe form and 57 asymptomatic-mild) using RFLP-PCR. The genotype distribution of all SNVs throughout the entire sample of 100 individuals showed Hardy Weinberg equilibrium (P=0.53, 0.35, and 0.4 for CCL2, OAS1, and DPP9, respectively). The HWE test did not find any statistically significant difference in genotype distribution between the study and control groups for any of the three SNVs. The multivariable logistic regression analysis showed that individuals with the GG of the CCL2 rs1024611 gene variant had an increased association with the severe COVID-19 phenotype in a recessive model (P = 0.0003, OR = 6.43, 95% CI 2.19-18.89) and for the OAS1 rs10774671 gene variant, the log-additive model showed a significant association with the severe phenotype of COVID-19 (P=0.0084, OR=3.85, 95% CI 1.33-11.12). Analysis of haplotype frequencies revealed that the coexistence of GAG at CCL2, OAS1, and DPP9 variants, respectively, in the same individual increased the presence of the severe COVID-19 phenotype (OR=2.273, 95% CI: 1.271-4.068, P=0.005305). The findings of the current study suggests that the ethnic background affects the allele and genotype frequencies of genes associated with the severity of COVID-19. The experience with COVID-19 has provided an opportunity to identify an ethnicity-based approach to recognize genetically high-risk individuals in different populations for emerging diseases.

Unexpected mutations occurred in CRISPR/Cas9 edited Drosophila analyzed by deeply whole genomic sequencing

CRISPR/Cas9 possesses the most promising prospects as a gene-editing tool in post-genomic researches. It becomes an epoch-marking technique for the features of speed and convenience of genomic modification. However, it is still unclear whether CRISPR/Cas9 gene editing can cause irreversible damage to the genome. In this study, we successfully knocked out the WHITE gene in Drosophila, which governs eye color, utilizing CRISPR/Cas9 technology. Subsequently, we conducted high-throughput sequencing to assess the impact of this editing process on the stability of the entire genomic profile. The results revealed the presence of numerous unexpected mutations in the Drosophila genome, including 630 SNVs (Single Nucleotide Variants), 525 Indels (Insertion and Deletion) and 425 MSIs (microsatellite instability). Although the KO (knockout) specifically occurred on chromosome X, the majority of mutations were observed on chromosome 3, indicating that this effect is genome-wide and associated with the spatial structure between chromosomes, rather than being solely limited to the location of the KO gene. It is worth noting that most of the mutations occurred in the intergenic and intron regions, without exerting any significant on the function or healthy of the animal. In addition, the mutations downstream of the knockout gene well beyond the upstream. This study has found that gene editing can lead to unexpected mutations in the genome, but most of these mutations are harmless. This research has deepened our understanding of CRISPR/Cas9 and broadened its application prospects.

Reviewing critical TRPM2 variants through a structure-function lens

The transient receptor potential melastatin 2 (TRPM2) channel plays a central role in connecting redox state with calcium signaling in living cells. This coupling makes TRPM2 essential for physiological functions such as pancreatic insulin secretion or cytokine production, but also allows it to contribute to pathological processes, including neuronal cell death or ischemia-reperfusion injury. Genetic deletion of the channel, albeit not lethal, alters physiological functions in mice. In humans, population genetic studies and whole-exome sequencing have identified several common and rare genetic variants associated with mental disorders and neurodegenerative diseases, including single nucleotide variants (SNVs) in exonic regions. In this review, we summarize available information on the four best-documented SNVs: one common (rs1556314) and three rare genetic variants (rs139554968, rs35288229, and rs145947009), manifested in amino acid substitutions D543E, R707C, R755C, and P1018L respectively. We discuss existing evidence supporting or refuting the associations between SNVs and disease. Furthermore, we aim to interpret the molecular impacts of these amino acid substitutions based on recently published structures of human TRPM2. Finally, we formulate testable hypotheses and suggest means to investigate them. Studying the function of proteins with rare mutations might provide insight into disease etiology and delineate new drug targets.

Factors Influencing Mortality in Children with Central Nervous System Tumors: A Cohort Study on Clinical Characteristics and Genetic Markers

Multidrug resistance (MDR) commonly leads to cancer treatment failure because cancer cells often expel chemotherapeutic drugs using ATP-binding cassette (ABC) transporters, which reduce drug levels within the cells. This study investigated the clinical characteristics and single nucleotide variant (SNV) in ABCB1, ABCC1, ABCC2, ABCC4, and ABCG2, and their association with mortality in pediatric patients with central nervous system tumors (CNST). Using TaqMan probes, a real-time polymerase chain reaction genotyped 15 SNPs in 111 samples. Patients were followed up until death or the last follow-up day using the Cox proportional hazards model. An association was found between the rs1045642 (ABCB1) in the recessive model (HR = 2.433, 95% CI 1.098-5.392, p = 0.029), and the ICE scheme in the codominant model (HR = 9.810, 95% CI 2.74-35.06, p ≤ 0.001), dominant model (HR = 6.807, 95% CI 2.87-16.103, p ≤ 0.001), and recessive model (HR = 6.903, 95% CI 2.915-16.544, p = 0.038) significantly increased mortality in this cohort of patients. An association was also observed between the variant rs3114020 (ABCG2) and mortality in the codominant model (HR = 5.35, 95% CI 1.83-15.39, p = 0.002) and the dominant model (HR = 4.421, 95% CI 1.747-11.185, p = 0.002). A significant association between the ICE treatment schedule and increased mortality risk in the codominant model (HR = 6.351, 95% CI 1.831-22.02, p = 0.004, HR = 9.571, 95% CI 2.856-32.07, p ≤ 0.001), dominant model (HR = 6.592, 95% CI 2.669-16.280, p ≤ 0.001), and recessive model (HR = 5.798, 95% CI 2.411-13.940, p ≤ 0.001). The genetic variants rs3114020 in the ABCG2 gene and rs1045642 in the ABCB1 gene and the ICE chemotherapy schedule were associated with an increased mortality risk in this cohort of pediatric patients with CNST.

Genome-Wide Meta-Analysis Identifies 11 Susceptibility Variants of Vitiligo in the Chinese Han Population

Vitiligo is an autoimmune disease involving loss of melanocytes. Although several genetic studies have confirmed that genetic factors play an important role, its pathogenesis remains incompletely characterized. In this study, a genome-wide meta-analysis was conducted to search for more susceptibility variants of vitiligo. Tang et al performed a GWAS for cohort I (1117 vitiligo cases and 1701 healthy controls) previously, and we conducted a GWAS for cohort II (3323 vitiligo cases and 7186 healthy controls) in this study, with the results subjected to a genome-wide meta-analysis and linkage disequilibrium analysis. We identify, to our knowledge, 11 previously unreported susceptibility variants, of which 6 variants are located in the intronic regions, and the remaining 5 variants are located within intergenic regions between genes. In addition, the results of polygenic risk score show that the best evaluated effect for target data is among significant SNVs of the base data. The susceptibility genes of vitiligo are mainly enriched in the immune-related functions and pathways. The susceptibility variants expand the role of genetic factors associated with vitiligo. The bioinformatics analysis for risk genes provides further insight into the pathogenesis of vitiligo.

Single Nucleotide Variants (SNVs) of the Mesocorticolimbic System Associated with Cardiovascular Diseases and Type 2 Diabetes: A Systematic Review

The mesocorticolimbic (MCL) system is crucial in developing risky health behaviors which lead to cardiovascular diseases (CVDs) and type 2 diabetes (T2D). Although there is some knowledge of the MCL system genes linked to CVDs and T2D, a comprehensive list is lacking, underscoring the significance of this review. This systematic review followed PRISMA guidelines and the Cochrane Handbook for Systematic Reviews of Interventions. The PubMed and Web of Science databases were searched intensively for articles related to the MCL system, single nucleotide variants (SNVs, formerly single nucleotide polymorphisms, SNPs), CVDs, T2D, and associated risk factors. Included studies had to involve a genotype with at least one MCL system gene (with an identified SNV) for all participants and the analysis of its link to CVDs, T2D, or associated risk factors. The quality assessment of the included studies was performed using the Q-Genie tool. The VEP and DAVID tools were used to annotate and interpret genetic variants and identify enriched pathways and gene ontology terms associated with the gene list. The review identified 77 articles that met the inclusion criteria. These articles provided information on 174 SNVs related to the MCL system that were linked to CVDs, T2D, or associated risk factors. The COMT gene was found to be significantly related to hypertension, dyslipidemia, insulin resistance, obesity, and drug abuse, with rs4680 being the most commonly reported variant. This systematic review found a strong association between the MCL system and the risk of developing CVDs and T2D, suggesting that identifying genetic variations related to this system could help with disease prevention and treatment strategies.

scSNV-seq: high-throughput phenotyping of single nucleotide variants by coupled single-cell genotyping and transcriptomics

CRISPR screens with single-cell transcriptomic readouts are a valuable tool to understand the effect of genetic perturbations including single nucleotide variants (SNVs) associated with diseases. Interpretation of these data is currently limited as genotypes cannot be accurately inferred from guide RNA identity alone. scSNV-seq overcomes this limitation by coupling single-cell genotyping and transcriptomics of the same cells enabling accurate and high-throughput screening of SNVs. Analysis of variants across the JAK1 gene with scSNV-seq demonstrates the importance of determining the precise genetic perturbation and accurately classifies clinically observed missense variants into three functional categories: benign, loss of function, and separation of function.

A Bioinformatics Toolkit for Next-Generation Sequencing in Clinical Oncology

Next-generation sequencing (NGS) has taken on major importance in clinical oncology practice. With the advent of targeted therapies capable of effectively targeting specific genomic alterations in cancer patients, the development of bioinformatics processes has become crucial. Thus, bioinformatics pipelines play an essential role not only in the detection and in identification of molecular alterations obtained from NGS data but also in the analysis and interpretation of variants, making it possible to transform raw sequencing data into meaningful and clinically useful information. In this review, we aim to examine the multiple steps of a bioinformatics pipeline as used in current clinical practice, and we also provide an updated list of the necessary bioinformatics tools. This resource is intended to assist researchers and clinicians in their genetic data analyses, improving the precision and efficiency of these processes in clinical research and patient care.

The predictive and prognostic role of single nucleotide gene variants of PD-1 and PD-L1 in patients with advanced melanoma treated with PD-1 inhibitors

Background

Despite having revolutionized the treatment paradigm for advanced melanoma, not all patients benefit from immune checkpoint inhibitor therapy. To date, there are no predictive biomarkers for response or the occurrence of immune-related adverse events (irAEs) to programmed cell death protein 1 (PD-1) inhibitors. Our aim was to investigate the predictive and prognostic role of single nucleotide variants (SNVs) of genes involved in the PD-1 axis.

Methods

We analysed, in metastatic melanoma patients treated with nivolumab or pembrolizumab, five PD-1 SNVs, namely PD1.3 G>A (rs11568821), PD1.5 C>T (rs2227981), PD1.6 G>A (rs10204525), PD1.7 T>C(rs7421861), PD1.10 C>G (rs5582977) and three programmed death-ligand 1 (PD-L1) SNVs: +8293 C>A (rs2890658), PD-L1 C>T (rs2297136) and PD-L1 G>C (rs4143815). Association of SNV genotypic frequencies with best overall response to PD-1 inhibitors and development of irAEs were estimated through a modified Poisson regression. A Cox regression modelling approach was applied to evaluate the SNV association with OS.

Results

A total of 125 patients with advanced melanoma were included in the analysis. A reduction in irAEs risk was observed in patients carrying the PD-L1 +8293 C/A genotype compared with those carrying the C/C genotype (risk ratio = 0.45; 95% CL 0.22-0.93; P = 0.031). A trend for a reduction in irAEs was also observed with the PD1.5 T allele (risk ratio = 0.70, 95% confidence limits 0.48-1.01 versus C allele). None of the SNVs was associated with response to therapy. Finally, a survival benefit was observed in patients harbouring the PD1.7 C/C genotype (hazard ratio = 0.37; 95% confidence limits 0.14-0.96; P = 0.028) in the homozygous model.

Conclusions

Our study showed that PD-1.5 and PD-L1 +8293 SNVs may play a role as a predictive biomarker of development of irAEs to PD-1 inhibitors. PD1.7 SNV may also be associated with a reduction of the risk of death, although further translational research is needed to confirm these results.

Exploring Novel Variants of the Cytochrome P450 Reductase Gene (POR) from the Genome Aggregation Database by Integrating Bioinformatic Tools and Functional Assays

Cytochrome P450 oxidoreductase (POR) is an essential redox partner for steroid and drug-metabolizing cytochromes P450 located in the endoplasmic reticulum. Mutations in POR lead to metabolic disorders, including congenital adrenal hyperplasia, and affect the metabolism of steroids, drugs, and xenobiotics. In this study, we examined approximately 450 missense variants of the POR gene listed in the Genome Aggregation Database (gnomAD) using eleven different in silico prediction tools. We found that 64 novel variants were consistently predicted to be disease-causing by most tools. To validate our findings, we conducted a population analysis and selected two variations in POR for further investigation. The human POR wild type and the R268W and L577P variants were expressed in bacteria and subjected to enzyme kinetic assays using a model substrate. We also examined the activities of several cytochrome P450 proteins in the presence of POR (WT or variants) by combining P450 and reductase proteins in liposomes. We observed a decrease in enzymatic activities (ranging from 35% to 85%) of key drug-metabolizing enzymes, supported by POR variants R288W and L577P compared to WT-POR. These results validate our approach of curating a vast amount of data from genome projects and provide an updated and reliable reference for diagnosing POR deficiency.

Diversity and Classification of Genetic Variations in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a neurodevelopmental condition with symptoms that affect the whole personality and all aspects of life. Although there is a high degree of heterogeneity in both its etiology and its characteristic behavioral patterns, the disorder is well-captured along the autistic triad. Currently, ASD status can be confirmed following an assessment of behavioral features, but there is a growing emphasis on conceptualizing autism as a spectrum, which allows for establishing a diagnosis based on the level of support need, free of discrete categories. Since ASD has a high genetic predominance, the number of genetic variations identified in the background of the condition is increasing exponentially as genetic testing methods are rapidly evolving. However, due to the huge amount of data to be analyzed, grouping the different DNA variations is still challenging. Therefore, in the present review, a multidimensional classification scheme was developed to accommodate most of the currently known genetic variants associated with autism. Genetic variations have been grouped according to six criteria (extent, time of onset, information content, frequency, number of genes involved, inheritance pattern), which are themselves not discrete categories, but form a coherent continuum in line with the autism spectrum approach.

Analysis of Nucleotide Variations in Human G-Quadruplex Forming Regions Associated with Disease States

While the role of G quadruplex (G4) structures has been identified in cancers and metabolic disorders, single nucleotide variations (SNVs) and their effect on G4s in disease contexts have not been extensively studied. The COSMIC and CLINVAR databases were used to detect SNVs present in G4s to identify sequence level changes and their effect on the alteration of the G4 secondary structure. A total of 37,515 G4 SNVs in the COSMIC database and 2378 in CLINVAR were identified. Of those, 7236 COSMIC (19.3%) and 457 (19%) of the CLINVAR variants result in G4 loss, while 2728 (COSMIC) and 129 (CLINVAR) SNVs gain a G4 structure. The remaining variants potentially affect the folding energy without affecting the presence of a G4. Analysis of mutational patterns in the G4 structure shows a higher selective pressure (3-fold) in the coding region on the template strand compared to the reverse strand. At the same time, an equal proportion of SNVs were observed among intronic, promoter, and enhancer regions across strands.

A common variant close to the "tripwire" linker region of NLRP1 contributes to severe COVID-19

OBJECTIVE AND DESIGN

The heterogeneity of response to SARS-CoV-2 infection is directly linked to the individual genetic background. Genetic variants of inflammasome-related genes have been pointed as risk factors for several inflammatory sterile and infectious disease. In the group of inflammasome receptors, NLRP1 stands out as a good novel candidate as severity factor for COVID-19 disease.

METHODS

To address this question, we performed an association study of NLRP1, DPP9, CARD8, IL1B, and IL18 single nucleotide variants (SNVs) in a cohort of 945 COVID-19 patients.

RESULTS

The NLRP1 p.Leu155His in the linker region, target of viral protease, was significantly associated to COVID-19 severity, which could contribute to the excessive cytokine release reported in severe cases.

CONCLUSION

Inflammasome genetic background contributes to individual response to SARS-CoV-2.

Deep Sequencing of Porcine Reproductive and Respiratory Syndrome Virus ORF7: A Promising Tool for Diagnostics and Epidemiologic Surveillance

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major concern worldwide. Control of PRRSV is a challenging task due to various factors, including the viral diversity and variability. In this study, we evaluated an amplicon library preparation protocol targeting the ORF7 region of both PRRSV species, Betaarterivirus suid 1 and Betaarterivirus suid 2. We designed tailed primers for a two-step PCR procedure that generates ORF7-specific amplicon libraries suitable for use on Illumina sequencers. We tested the method with serum samples containing common laboratory strains and with pooled serum samples (n = 15) collected from different pig farms during 2019-2021 in Hungary. Testing spiked serum samples showed that the newly designed method is highly sensitive and detects the viral RNA even at low copy numbers (corresponding to approx. Ct 35). The ORF7 sequences were easily assembled even from clinical samples. Two different sequence variants were identified in five samples, and the Porcilis MLV vaccine strain was identified as the minor variant in four samples. An in-depth analysis of the deep sequencing results revealed numerous polymorphic sites along the ORF7 gene in a total of eight samples, and some sites (positions 12, 165, 219, 225, 315, 345, and 351) were found to be common in several clinical specimens. We conclude that amplicon deep sequencing of a highly conserved region of the PRRSV genome could support both laboratory diagnosis and epidemiologic surveillance of the disease.

Intra-host mutation rate of acute SARS-CoV-2 infection during the initial pandemic wave

SARS-CoV-2 mutation is minimized through a proofreading function encoded by NSP-14. Most estimates of the SARS-CoV-2 mutation rate are derived from population based sequence data. Our understanding of SARS-CoV-2 evolution might be enhanced through analysis of intra-host viral mutation rates in specific populations. Viral genome analysis was performed between paired samples and mutations quantified at allele frequencies (AF) ≥ 0.25, ≥ 0.5 and ≥ 0.75. Mutation rate was determined employing F81 and JC69 evolution models and compared between isolates with (ΔNSP-14) and without (wtNSP-14) non-synonymous mutations in NSP-14 and by patient comorbidity. Forty paired samples with median interval of 13 days [IQR 8.5-20] were analyzed. The estimated mutation rate by F81 modeling was 93.6 (95%CI 90.8-96.4], 40.7 (95%CI 38.9-42.6) and 34.7 (95%CI 33.0-36.4) substitutions/genome/year at AF ≥ 0.25, ≥ 0.5, ≥ 0.75 respectively. Mutation rate in ΔNSP-14 were significantly elevated at AF ≥ 0.25 vs wtNSP-14. Patients with immune comorbidities had higher mutation rate at all allele frequencies. Intra-host SARS-CoV-2 mutation rates are substantially higher than those reported through population analysis. Virus strains with altered NSP-14 have accelerated mutation rate at low AF. Immunosuppressed patients have elevated mutation rate at all AF. Understanding intra-host virus evolution will aid in current and future pandemic modeling.

Transcription-wide impact by RESCUE-induced off-target single-nucleotide variants in mammalian cells

RNA base editing is a promising tool in precise molecular therapy. Currently, there are two widely used RNA base editors, REPAIR and RESCUE. REPAIR only facilitates A-to-I conversions, while RESCUE performs both A-to-I and C-to-U conversions. Thus, RESCUE can generate twice the number of mutations compared to REPAIR. However, transcription-wide impact due to RESCUE-induced off-target single-nucleotide variants (SNVs) is not fully appreciated. Therefore, to determine the off-target effects of RESCUE-mediated editing, we employed transcription-wide sequencing on cells edited by RESCUE. The SNVs showed different off-target effects on mRNA, circRNA, lncRNA, and miRNA expression patterns and their interacting networks. Our results illustrate the transcription-wide impact of RESCUE-induced off-target SNVs and highlight the need for careful characterization of the off-target impact by this editing platform.

Congruence Concept for Comparison of Spectra

This paper introduces an alternative, easy-to-implement spectrum comparison concept. The evaluation procedure is illustrated by artificial and attenuated total reflection Fourier transform infrared (ATR FT-IR) spectra, which it can also be extended to other spectrometries (e.g., ultraviolet-visible or UV-Vis and Raman). The evaluation for the comparison of two spectra is divided into four phases: (i) spectrum pre-treatment (e.g., smoothing and background correction), (ii) standard normal variate (SNV) transformation, (iii) regression analysis of SNV spectra, and (iv) calculation of the quantification index (FG). The FG is derived from the formula of R². It characterizes and quantifies the identity and/or similarity of the compared spectra.

Pharmacogenetics of ABCB1, CDA, DCK, GSTT1, GSTM1 and outcomes in a cohort of pediatric acute myeloid leukemia patients from Colombia

BACKGROUND AND AIM

Different studies have shown pharmacogenetic variants related to drug toxicity in acute myeloid leukemia (AML) patients. Our aim was to identify the association between ABCB1, CDA, DCK, GSTT1, and GSTM1 variants with clinical outcomes and toxicity in pediatric patients with AML.

METHODS

Fifty-one confirmed de novo AML pediatric patients were included. A SNaPshot™ assay and conventional PCR were used to evaluate ABCB1, CDA, DCK, GSTT1, and GSTM1 variants. Clinical outcomes and toxicity associations were evaluated using odds ratios and Chi-square analysis.

RESULTS

Patients carrying ABCB1 (1236C > T, rs1128503) GG genotype in had a 6.8 OR (CI 95% 1.08-42.73, p = .044) for cardiotoxicity as compared to patients carrying either AA or GA genotypes 0.14 OR (CI 95% 0.023-0.92, p = .044). For ABCB1 (1236G > A rs1128503/2677C > A/T rs2032582/3435G > A rs1045642) AA/AA/AA combined genotypes had a strong association with death after HSTC OR 13.73 (CI 95% 1.94-97.17, p = .009). Combined genotypes GG/CC/GG with CDA (79A > C, rs2072671) CA genotype or CDA (-451G > A, rs532545) CT genotype, had a 4.11 OR (CI 95% 2.32-725, p = .007) and 3.8 OR (CI 95% 2.23-6.47, p = .027) with MRD >0.1% after first chemotherapy cycle, respectively.

CONCLUSION

Our results highlight the importance of pharmacogenetic analysis in pediatric AML, particularly in populations with a high degree of admixture, and might be useful as a future tool for patient stratification for treatment.

FixItFelix: improving genomic analysis by fixing reference errors

The current version of the human reference genome, GRCh38, contains a number of errors including 1.2 Mbp of falsely duplicated and 8.04 Mbp of collapsed regions. These errors impact the variant calling of 33 protein-coding genes, including 12 with medical relevance. Here, we present FixItFelix, an efficient remapping approach, together with a modified version of the GRCh38 reference genome that improves the subsequent analysis across these genes within minutes for an existing alignment file while maintaining the same coordinates. We showcase these improvements over multi-ethnic control samples, demonstrating improvements for population variant calling as well as eQTL studies.

Leishmaniasis: Omics approaches to understand its biology from molecule to cell level

Leishmaniasis is the second deadliest vector-borne, neglected tropical zoonotic disease and is found in a variety of clinical forms based on genetic background. Its endemic type is present in tropical, sub-tropical and Mediterranean areas around the world which accounts for a lot of deaths every year. Currently, a variety of techniques are present to detect leishmaniasis with advantages and disadvantages. The advancing next-generation sequencing (NGS) techniques are employed to find out novel diagnostic markers based on single nucleotide variants. Total 274 NGS studies are available in European Nucleotide Archive (ENA) portal (https://www.ebi.ac.uk/ena/browser/home) that focused on wild-type and mutated Leishmania, differential gene expression, miRNA expression, and detection of aneuploidy mosaicism by omics approaches. These studies have provided insights into the population structure, virulence, and extensive structural variation, including known and suspected drug resistance loci, mosaic aneuploidy and hybrid formation under stressed conditions and inside the midgut of the sandfly. The complex interactions occurring within the parasite-host-vector triangle can be better understood by omics approaches. Further, advanced CRISPR technology allows researchers to delete and modify each gene individually to know the importance of genes in the virulence and survival of the disease-causing protozoa. In vitro Generation of Leishmania, hybrids are helping to understand the mechanism of disease progression in its different stages of infection. This review will give a comprehensive picture of the available omics data of various Leishmania spp. which helped to reveal the effect of climate change on the spread of its vector, the pathogen survival strategies, emerging antimicrobial resistance and its clinical importance.

CriSNPr, a single interface for the curated and de novo design of gRNAs for CRISPR diagnostics using diverse Cas systems

CRISPR-based diagnostics (CRISPRDx) have improved clinical decision-making, especially during the COVID-19 pandemic, by detecting nucleic acids and identifying variants. This has been accelerated by the discovery of new and engineered CRISPR effectors, which have expanded the portfolio of diagnostic applications to include a broad range of pathogenic and non-pathogenic conditions. However, each diagnostic CRISPR pipeline necessitates customized detection schemes based on the fundamental principles of the Cas protein used, its guide RNA (gRNA) design parameters, and the assay readout. This is especially relevant for variant detection, a low-cost alternative to sequencing-based approaches for which no in silico pipeline for the ready-to-use design of CRISPRDx currently exists. In this manuscript, we fill this lacuna using a unified web server, CriSNPr (CRISPR-based SNP recognition), which provides the user with the opportunity to de novo design gRNAs based on six CRISPRDx proteins of choice (Fn/enFnCas9, LwCas13a, LbCas12a, AaCas12b, and Cas14a) and query for ready-to-use oligonucleotide sequences for validation on relevant samples. Furthermore, we provide a database of curated pre-designed gRNAs as well as target/off-target for all human and SARS-CoV-2 variants reported thus far. CriSNPr has been validated on multiple Cas proteins, demonstrating its broad and immediate applicability across multiple detection platforms. CriSNPr can be found at http://crisnpr.igib.res.in/.

Evidence for etiologic field changes in tongue distant from tumor in patients with squamous cell carcinoma of the oral tongue

Oral cancer is a paradigm of Slaughter's concept of field cancerization, where tumors are thought to originate within an area of cells containing genetic alterations that predispose to cancer development. The field size is unclear but may represent a large area of tissue, and the origin of mutations is also unclear. Here, we analyzed whole exome and transcriptome features in contralateral tumor-distal tongue (i.e. distant from the tumor, not tumor-adjacent) and corresponding tumor tissues of 15 patients with squamous cell carcinoma of the oral tongue. The number of point mutations ranged from 41 to 237 in tumors and from one to 78 in tumor-distal samples. Tumor-distal samples showed mainly clock-like (associated with aging) or tobacco smoking mutational signatures. Tumors additionally showed mutations that associate with cytidine deaminase AID/APOBEC enzyme activities or a UV-like signature. Importantly, no point mutations were shared between a tumor and the matched tumor-distal sample in any patient. TP53 was the most frequently mutated gene in tumors (67%), whereas a TP53 mutation was detected in only one tumor-distal sample, and this mutation was not shared with the matched tumor. Arm-level copy number variation (CNV) was found in 12 tumors, with loss of chromosome (Chr) 8p or gain of 8q being the most frequent events. Two tumor-distal samples showed a gain of Chr8, which was associated with increased expression of Chr8-located genes in these samples, although gene ontology did not show a role for these genes in oncogenic processes. In situ hybridization revealed a mixed pattern of Chr8 gain and neutral copy number in both tumor cells and adjacent nontumor epithelium in one patient. We conclude that distant field cancerization exists but does not present as tumor-related mutational events. The data are compatible with etiologic field effects, rather than classical monoclonal field cancerization theory. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Genetic diversity in the metronidazole metabolism genes nitroreductases and pyruvate ferredoxin oxidoreductases in susceptible and refractory clinical samples of Giardia lamblia

The effectiveness of metronidazole against the tetraploid intestinal parasite Giardia lamblia is dependent on its activation/inactivation within the cytoplasm. There are several activating enzymes, including pyruvate ferredoxin reductase (PFOR) and nitroreductase (NR) 1 which metabolize metronidazole into toxic forms, while NR2 on the other hand inactivates it. Metronidazole treatment failures have been increasing rapidly over the last decade, indicating genetic resistance mechanisms. Analyzing genetic variation in the PFOR and NR genes in susceptible and refractory Giardia isolates may help identify potential markers of resistance. Full length PFOR1, PFOR2, NR1 and NR2 genes from clinical culturable isolates and non-cultured clinical Giardia assemblage B samples were cloned, sequenced and single nucleotide variants (SNVs) were analyzed to assess genetic diversity and alleles. A similar ratio of amino acid changing SNVs per gene length was found for the NRs; 4.2% for NR1 and 6.4% for NR2, while the PFOR1 and PFOR2 genes had less variability with a ratio of 1.1% and 1.6%, respectively. One of the samples from a refractory case had a nonsense mutation which caused a truncated NR1 gene in one out of six alleles. Further, we found three NR2 alleles with frameshift mutations, possibly causing a truncated protein in two susceptible isolates. One of these isolates was homozygous for the affected NR2 allele. Three nsSNVs with potential for affecting protein function were found in the ferredoxin domain of the PFOR2 gene. The considerable variation and discovery of mutations possibly causing dysfunctional NR proteins in clinical Giardia assemblage B isolates, reveal a potential for genetic link to metronidazole susceptibility and resistance.

Modulators of radiation-induced cardiopulmonary toxicities for non-small cell lung cancer: Integrated cytokines, single nucleotide variants, and HBP systems imaging

Radiation therapy (RT)-induced cardiopulmonary toxicities remain dose-limiting toxicities for patients receiving radiation dosages to the thorax, especially for lung cancer. Means of monitoring and predicting for those receiving RT or concurrent chemoradiation therapy before treatment begins in individual patients could benefit early intervention to prevent or minimize RT-induced side effects. Another aspect of an individual's susceptibility to the adverse effects of thoracic irradiation is the immune system as reflected by phenotypic factors (patterns of cytokine expressions), genotypic factors (single nucleotide variants SNVs; formerly single nucleotide polymorphisms [SNPs]), and aspects of quantitative cellular imaging. Levels of transcription, production, and functional activity of cytokines are often influenced by SNVs that affect coding regions in the promoter or regulatory regions of cytokine genes. SNVs can also lead to changes in the expression of the inflammatory cytokines, interferons, interleukins (IL-6, IL-17) and tumor necrosis factors (TNF-α) at the protein level. RT-induced cardiopulmonary toxicities could be quantified by the uptake of ¹⁸F-fluorodeoxyglucose (FDG), however, FDG is a sensitive but not specific biomarker in differential diagnosis between inflammation/infection and tumor recurrence. FDG is suitable for initial diagnosis of predisposed tissue injuries in non-small cell lung cancer (NSCLC). ^99mTc-ethylenedicysteine-glucosamine (^99mTc-EC-G) was able to measure tumor DNA proliferation and myocardial ischemia via hexosamine biosynthetic pathways (HBP). Thus, ^99mTc-EC-G could be an alternative to FDG in the assessment of RT doses and select patients in HBP-directed targets for optimal outcomes. This article reviewed correlative analyses of pro-inflammatory cytokines, genotype SNVs, and cellular imaging to improve the diagnosis, prognosis, monitoring, and prediction of RT-induced cardiopulmonary toxicities in NSCLC.

A Single Variant in Pri-miRNA-155 Associated with Susceptibility to Hereditary Breast Cancer Promotes Aggressiveness in Breast Cancer Cells

Variants in genes encoding for microRNAs have been associated with their deregulation in breast cancer (BC). Sequencing of microRNAs deregulated in BC was performed using DNA from Chilean patients with a strong family history and negative for mutations in BRCA1/BRCA2. Seventeen variants were identified, three of which were selected for a case-control association study: rs376491654 (miR-335), rs755634302 (miR-497), and rs190708267 (miR-155). For rs190708267 C>T, the heterozygous T allele was detected in four BC cases and absent in controls, while homozygous TT cases were not detected. Variants were modelled in silico, cloned in a plasmid, expressed in BC cell lines, and functional in vitro assays were performed. Overexpression of the miR-155-T allele increased mature miR-155-5p levels in both BC cell lines, suggesting that its presence alters pre-miR-155 processing. Moreover, BC cells overexpressing the miR-155-T allele showed increased proliferation, migration, and resistance to cisplatin-induced death compared to miR-155-C overexpressing cells. Of note, the 3′UTR of APC, GSK3β, and PPP1CA genes, all into the canonical Wnt signaling pathway, were identified as direct targets. APC and GSK3β mRNA levels decreased while PP1 levels increased. These results suggest a pathogenic role of the variant rs190708267 (miR-155) in BRCA 1/2 negative BC, conferring susceptibility and promoting traits of aggressiveness.

Transgenerational genomic analyses reveal allelic oscillation and purifying selection in a gut parasite Nosema ceranae

Standing genetic variation is the predominant source acted on by selection. Organisms with high genetic diversity generally show faster responses toward environmental change. Nosema ceranae is a microsporidian parasite of honey bees, infecting midgut epithelial cells. High genetic diversity has been found in this parasite, but the mechanism for the parasite to maintain this diversity remains unclear. This study involved continuous inoculation of N. ceranae to honey bees. We found that the parasites slowly increased genetic diversity over three continuous inoculations. The number of lost single nucleotide variants (SNVs) was balanced with novel SNVs, which were mainly embedded in coding regions. Classic allele frequency oscillation was found at the regional level along the genome, and the associated genes were enriched in apoptosis regulation and ATP binding. The ratio of synonymous and non-synonymous substitution suggests a purifying selection, and our results provide novel insights into the evolutionary dynamics in microsporidian parasites.

Random distribution of nucleotide polymorphism throughout the genome of tomato-infecting begomovirus species occurring in India: implication in PCR based diagnosis

Multiple begomovirus species are known to cause leaf curl disease in tomato in India. In order to develop specific and generic PCR based diagnostics for the tomato-infecting begomoviruses, in this study, we attempted to design primers initially based on the multiple alignment of the complete genome sequence of DNA-A component. However, the specific nucleotide stretches adequate for preparing specific primers could not be obtained. Alternatively, the online Primer-BLAST tool that offers designing of target-specific PCR primers was attempted to prepare specific primers targeting three clones (DNA-A) of tomato-infecting begomovirus species (Tomato leaf curl New Delhi virus, Tomato leaf curl Palampur virus and Tomato leaf curl Joydebpur virus) selected based on their sequence identity and phylogenetic relatedness. The primers derived from Primer-BLAST tool showed high level of cross-reaction among these begomovirus species and therefore were not able to differentiate these target begomovirus species. In order to understand the reason of cross-reactivity further sequence analysis revealed the high occurrence of single nucleotide variations (SNVs) compared to the multi-nucleotide stretches. There was no SNV hot-spot in the genome, rather the SNVs were randomly distributed throughout the genome of these begomovirus species. This pattern of nucleotide diversities among these tomato-infecting begomoviruses seriously implicated on developing specific PCR diagnostics. On the contrary, sequence analysis showed high sequence conservancy, which enabled to develop a generic PCR diagnostic for these begomoviruses. Our study, thus showed that the genome sequence diversity pattern among the tomato-infecting begomoviruses in India poses challenges in developing PCR based specific diagnostics.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-022-00785-9.

High-coverage whole-genome sequencing of the expanded 1000 Genomes Project cohort including 602 trios

The 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies.

Copy Number Analysis in a Large Cohort Suggestive of Inborn Errors of Immunity Indicates a Wide Spectrum of Relevant Chromosomal Losses and Gains

Inborn errors of immunity (IEI) are genetically driven disorders. With the advancement of sequencing technologies, a rapidly increasing number of gene defects has been identified, thereby mirroring the high heterogeneity in immunological and clinical presentations observed in patients. However, for a large majority of patients, no causative single nucleotide variant (SNV) or small indel can be identified using next-generation sequencing. First studies have shown that also copy number variants (CNVs) can cause IEI. Unfortunately, CNVs are not well examined in many routine diagnostic settings and the aim of this study was to assess the number of clinically relevant chromosomal losses and gains in a large cohort. We identified a total of 20 CNVs using whole exome sequencing data of a cohort of 191 patients with a suspected IEI. A definite molecular diagnosis could be made in five patients (2.6%), including pathogenic deletions affecting ICOS, TNFAIP3, and 22q11.2. CNVs of uncertain significance were observed in fifteen patients (7.9%), including deletions of 11q22.1q22.3 and 16p11.2 but also duplications affecting entire or parts of genes previously associated with IEI. Importantly, five patients carrying a CNV of uncertain significance also carried pathogenic or likely pathogenic SNVs (PIK3R1, NFKB1, NLRC4, DOCK2), or SNVs of unknown significance (NFKB2). This cooccurrence of SNVs and CNVs suggests modifying effects in some patients, and functional follow-up is warranted now in order to better understand phenotypic heterogeneity. In summary, the diagnostic yield of IEI can be increased substantially by evaluating CNVs, which allows an improved therapeutic management in those patients.

Characterization of the Immunological Status of Hypermutated Solid Tumors in the Cancer Genome Analysis Project HOPE

BACKGROUND

Many reports demonstrate that a high tumor mutation burden (TMB-H) is closely associated with good prognosis of cancer. However, specific studies investigating the association of various TMB statuses with overall survival in patients with solid tumors are scarce.

PATIENTS AND METHODS

In the present study, we investigated the association of TMB status with overall survival in 5,072 patients with cancer from the HOPE project and clarified the specific mechanism responsible for the good prognosis of the TMB-H group. All tumors were classified into one of four groups based on TMB: ultralow (UL), low (L), intermediate (I) and high (H).

RESULTS

The TMB-H group had a better prognosis than the TMB-I and TMB-L groups, but not than the TMB-UL group. Analyzing the expression of 293 immune response-associated genes, 17 genes were up-regulated in the TMB-H group compared to the TMB-I and TNB-L groups, and two genes [CD274 and interferon-γ (IFNG)] were identified as good prognostic factors. Analysis of immune cell populations inside tumors demonstrated that the frequencies of exhausted CD8⁺ T-cells, activated effector CD8⁺ T-cells and natural killer cells were significantly higher in the TMB-H group. The T-cell receptor repertoire numbers and the diversity evenness score (DE50) were lower in the TMB-H group than in TMB-UL group; however, no association of the DE50 value with the binding or elution affinity of epitope peptides from neoantigens was found.

CONCLUSION

One possible mechanism for the good prognosis of the TMB-UL group compared to the TMB-H group might be that the TMB-UL group features a balance between immunosuppression and immunostimulation.

Genetic Alterations in Papillary Thyroid Carcinoma With Hashimoto's Thyroiditis： ANK3, an Indolent Maintainer of Papillary Thyroid Carcinoma

Hashimoto’s thyroiditis (TH) is a risk factor for the occurrence of papillary thyroid carcinoma (PTC), which is considered to be the most common type of thyroid cancer. In recent years, the prevalence of PTC with TH has been increasing, but little is known about the genetic alteration in PTC with TH. This study analyzed the mutation spectrum and mutation signature of somatic single nucleotide variants (SNV) for 10 non-tumor and tumor pair tissues of PTC with TH using whole-exome sequencing. The ANK3 protein expression was evaluated by immunohistochemistry in PTC with TH and PTC samples. Moreover, the functional role of ANK3 in PTC cells was determined by CCK-8 proliferation assay, colony formation assays, cell cycle analysis, cell invasion and migration and in vivo study through overexpression assay. Our results showed three distinct mutational signatures and the C>T/G>A substitution was the most common type of SNV. Gene-set enrichment analysis showed that most of the significantly mutated genes were enriched in the regulation of actin cytoskeleton signaling. Moreover, NCOR2, BPTF, ANK3, and PCSK5 were identified as the significantly mutated genes in PTC with TH, most of which have not been previously characterized. Unexpectedly, it was found that ANK3 was overexpressed in cytoplasm close to the membrane of PTC cells with TH and in almost all PTC cases, suggesting its role as a diagnostic marker of PTC. Ectopic expression of ANK3 suppressed invasion and migration, increased apoptosis of B-CPAP and TPC-1 cells. Moreover, our findings revealed that enhanced ANK3 expression inhibits growth of PTC cells both in vitro and in vivo. Ectopic expression of ANK3 significantly enhanced E-cadherin protein expression and inhibited PTC progression, at least in part, by suppression of epithelial-mesenchymal transition (EMT). Our study shows that ANK3 exerts an anti-oncogenic role in the development of PTC and might be an indolent maintainer of PTC.

Multi-Thoracolumbar Variations and NR6A1 Gene Polymorphisms Potentially Associated with Body Size and Carcass Traits of Dezhou Donkey

The number of thoracolumbar vertebrae is a quantitative trait positively correlated with the economic traits of livestock. More thoracolumbar vertebrae individuals could genetically be used to improve the livestock population, as more thoracolumbar vertebrae means a longer carcass, which could bring more meat production. Nuclear receptor subfamily 6 group A member 1 (NR6A1) is considered a strong candidate gene for effecting the number of vertebrae in livestock. The purposes of this study are as follows: (a) Analyzing the effect of TLN variation on body size and carcass traits of Dezhou donkey; (b) Studying the distribution of seven single nucleotide variants (SNVs) in NR6A1 gene of Dezhou donkey; (c) Exploring the relationship between latent SNVs and TLN, the body size and carcass traits. We examined the thoracic and lumbar vertebrae number and seven SNVs in NR6A1 gene of 455 Dezhou donkeys, and analyzed the relationships between them. Five types of thoracolumbar combinations (T17L5 (individual with 17 thoracic and five lumbar vertebrae) 2.4%, T18L5 75.8%, T19L5 1.1%, T17L6 11.9%, and T18L6 8.8%) of Dezhou donkeys were detected in this study. For one thoracolumbar vertebra added, the body length of Dezhou donkey increases by 3 cm and the carcass weight increases by 6 kg. Seven SNVs (g.18093100G > T, g.18094587G > T, g.18106043G > T, g.18108764G > T, g.18110615T > G, g.18112000C > T and g.18114954T > G) of the NR6A1 gene were found to have a significant association with the TLN, body size and carcass traits of Dezhou donkey (p < 0.05), respectively. For instance, g.18114954C > T is significantly associated with lumber vertebrae number, the total number of thoracolumbar, and carcass weight, and individuals with TT genotype had significantly larger value than CC genotype (p < 0.05). Using these 7SNVs, 16 different haplotypes were estimated. Compared to Hap3Hap3, individuals homozygous for Hap2Hap2 showed significantly longer length in one thoracic spine (STL), the total thoracic vertebrae and one thoracolumbar spine. Our study will not only extend the understanding of genetic variation in the NR6A1 gene of Dezhou donkey, but also provide useful information for marker assisted selection in donkey breeding program.

A Systematic Screening of ADHD-Susceptible Variants From 25 Chinese Parents-Offspring Trios

Attention-deficit/hyperactivity disorder (ADHD) is one of the most prevalent and heritable childhood behavioral disorders. Although a number of ADHD-susceptible regions had been identified, details about the variations of genes and their related patterns involved in ADHD are still lacking. In this study, we collected 25 Chinese parents-offspring trios, each of which consisted of a child diagnosed with ADHD and his/her unaffected parents, and analyzed the variations from whole-genome sequencing data. SNVs in reported ADHD-susceptible regions and on the genes whose functions were related to dopamine were screened, and we identified a set of variants with functional annotations which were specifically detected in ADHD children, including most SNVs in the gene coding region that might impair protein functions and a few SNVs in promoter or 3' untranslated region (3'-UTR) that might affect the regulation of relative gene expression in a transcriptional or posttranscriptional level. All the information may further contribute to the understanding, prediction, prevention, and treatment of ADHD in clinical.

FAT1 and MSH2 Are Predictive Prognostic Markers for Chinese Osteosarcoma Patients Following Chemotherapeutic Treatment

Osteosarcoma is characterized by diverse genetic mutations, including single-nucleotide variants (SNVs), which can complicate clinical outcomes of the treatment. This study identified key mutations or polymorphisms in genes that correlate with osteosarcoma prognoses. A total of 110 patients with osteosarcoma were assigned to "good" or "poor" cohorts depending on their 5-year disease-free survival (DFS) after surgery and chemotherapeutic treatment. We performed next-generation sequencing analysis of tumor tissues for prognosis-associated SNVs in 315 tumorigenesis-related genes, followed by modeling of clinical outcomes for these patients using random forest classification via a support vector machine (SVM). Data from the Chinese Millionome Database were used to compare SNV frequency in osteosarcoma patients and healthy people. SVM screening identified 17 nonsynonymous SNVs located in 15 genes, of which rs17224367 and rs3733406 (located in MSH2 and FAT1, respectively) were strongly correlated with osteosarcoma prognosis. These results were verified in a 26-patient validation cohort, confirming that these SNVs could be used to predict prognosis. These results demonstrated that two SNVs located in MSH2 and FAT1 are associated with prognosis of osteosarcoma patients. © 2022 American Society for Bone and Mineral Research (ASBMR).

Homophilic Interaction Between Transmembrane-JAM-A and Soluble JAM-A Regulates Thrombo-Inflammation: Implications for Coronary Artery Disease

Genetic predisposition through F11R-single-nucleotide variation (SNV) influences circulatory soluble junctional adhesion molecule-A (sJAM-A) levels in coronary artery disease (CAD) patients. Homozygous carriers of the minor alleles (F11R-SNVs rs2774276, rs790056) show enhanced levels of thrombo-inflammatory sJAM-A. Both F11R-SNVs and sJAM-A are associated with worse prognosis for recurrent myocardial infarction in CAD patients. Platelet surface-associated JAM-A correlate with platelet activation markers in CAD patients. Activated platelets shed transmembrane-JAM-A, generating proinflammatory sJAM-A and JAM-A-bearing microparticles. Platelet transmembrane-JAM-A and sJAM-A as homophilic interaction partners exaggerate thrombotic and thrombo-inflammatory platelet monocyte interactions. Therapeutic strategies interfering with this homophilic interface may regulate thrombotic and thrombo-inflammatory platelet response in cardiovascular pathologies where circulatory sJAM-A levels are elevated.

Clinical impact of variants in non-coding regions of SHOX - Current knowledge

The short stature homeobox-containing (SHOX) is the most frequently analysed gene in patients classified as short stature patients (ISS) or diagnosed with Leri-Weill dyschondrosteosis (LWD), Langer mesomelic dysplasia (LMD), or Madelung deformity (MD). However, clinical testing of this gene focuses primarily on single nucleotide variants (SNV) in its coding sequences and copy number variants (CNV) overlapping SHOX gene. This review summarizes the clinical impact of variants in noncoding regions of SHOX. RECENT FINDINGS: CNV extending exclusively into the regulatory elements (i.e., not interrupting the coding sequence) are found more frequently in downstream regulatory elements of SHOX. Further, duplications are more frequent than deletions. Interestingly, downstream duplications are more common than deletions in patients with ISS or LWD but no such differences exist for upstream CNV. Moreover, the presence of specific CNVs in the patient population suggests the involvement of additional unknown factors. Some of its intronic variants, notably NM_000451.3(SHOX):c.-9delG and c.-65C>A in the 5'UTR, have unclear clinical roles. However, these intronic SNV may increase the probability that other CNV will arise de novo in the SHOX gene based on homologous recombination or incorrect splicing of mRNA. SUMMARY: This review highlights the clinical impact of noncoding changes in the SHOX gene and the need to apply new technologies and genotype-phenotype correlation in their analysis.

Advanced Functions Embedded in the Second Version of Database, Global Evaluation of SARS-CoV-2/hCoV-19 Sequences 2

The Global Evaluation of SARS-CoV-2/hCoV-19 Sequences 2 (GESS v2 https://shiny.ph.iu.edu/GESS_v2/) is an updated version of GESS, which has offered a handy query platform to analyze single-nucleotide variants (SNVs) on millions of high coverages and high-quality severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) complete genomes provided by the Global Initiative on Sharing Avian Influenza Data (GISAID). Including the tools in the first version, the GESS v2 is embedded with new functions, which allow users to search SNVs, given the viral nucleotide or amino acid sequence. The GESS v2 helps users to identify SNVs or SARS-CoV-2 lineages enriched in countries of user's interest and show the migration path of a selected lineage on a world map during specific time periods chosen by the users. In addition, the GESS v2 can recognize the dynamic variations of newly emerging SNVs in each month to help users monitor SNVs, which will potentially become dominant soon. More importantly, multiple sets of analyzed results about SNVs can be downloaded directly from the GESS v2 by which users can conduct their own independent research. With these significant updates, the GESS v2 will continue to serve as a public open platform for researchers to explore SARS-CoV-2 evolutionary patterns from the perspectives of the prevalence and impact of SNVs.

Single-Cell Sequencing Technologies in Precision Oncology

Single-cell sequencing technologies are revolutionizing cancer research and are poised to become the standard for translational cancer studies. Rapidly decreasing costs and increasing throughput and resolution are paving the way for the adoption of single-cell technologies in clinical settings for personalized medicine applications. In this chapter, we review the state of the art of single-cell DNA and RNA sequencing technologies, the computational tools to analyze the data, and their potential application to precision oncology. We also discuss the advantages of single-cell over bulk sequencing for the dissection of intra-tumor heterogeneity and the characterization of subclonal cell populations, the implementation of targeted drug repurposing approaches, and describe advanced methodologies for multi-omics data integration and to assess cell signaling at single-cell resolution.

Association between WWOX/MAF variants and dementia-related neuropathologic endophenotypes

The genetic locus containing the WWOX and MAF genes was implicated as a clinical Alzheimer's disease (AD) risk locus in two recent large meta-analytic genome wide association studies (GWAS). In a prior GWAS, we identified a variant in WWOX as a suggestive risk allele for hippocampal sclerosis. We hypothesized that the WWOX/MAF locus may be preferentially associated with non-plaque- and non-tau-related neuropathological changes (NC). Data from research participants with GWAS and autopsy measures from the National Alzheimer's Coordinating Center and the Religious Orders Study and the Rush Memory and Aging Project were meta-analyzed. Notably, no variants in the locus were significantly associated with ADNC. However, several WWOX/MAF variants had significant adjusted associations with limbic-predominant age-related TDP-43 encephalopathy NC (LATE-NC), HS, and brain arteriolosclerosis. These associations remained largely unchanged after adjustment for ADNC (operationalized with standard semiquantitative staging), suggesting that these associations are independent of ADNC. Thus, WWOX genetic variants were associated pathologically with LATE-NC, not ADNC.

Variant Interpretation: UCSC Genome Browser Recommended Track Sets

The UCSC Genome Browser has been an important tool for genomics and clinical genetics since the sequence of the human genome was first released in 2000. As it has grown in scope to display more types of data it has also grown more complicated. The data, which are dispersed at many locations worldwide, are collected into one view on the Browser, where the graphical interface presents the data in one location. This supports the expertise of the researcher to interpret variants in the genome. Because the analysis of Single Nucleotide Variants (SNVs) and Copy Number Variants (CNVs) require interpretation of data at very different genomic scales, different data resources are required. We present here several Recommended Track Sets designed to facilitate the interpretation of variants in the clinic, offering quick access to datasets relevant to the appropriate scale. This article is protected by copyright. All rights reserved.

Association of αENaC p. Ala663Thr Gene Polymorphism With Sudden Sensorineural Hearing Loss

Objective: The etiology of sudden sensorineural hearing loss (SSNHL) is still unknown. It has been demonstrated that normal endolymph metabolism is essential for inner ear function and that epithelial sodium channels (ENaC) may play an important role in the regulation of endolymphatic Na⁺. This study aimed to explore the potential association between αENaC p. Ala663Thr gene polymorphism and SSNHL. Methods: Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) was used to examine the genotype and allele frequency of the αENaC p. Ala663Thr polymorphism in 20 cases of low-frequency SSNHL (LF-SSNHL), 19 cases of high-frequency SSNHL (HF-SSNHL), 31 cases of all frequency SSNHL (AF-SSNHL), 42 cases of profound deafness SSNHL (PD-SSNHL), and 115 normal controls. Results: The T663 allele was found to be significantly associated with an increased risk of LF-SSNHL (p = 0.046, OR = 2.16, 95% CI = 1.01-4.62). The TT genotype and T663 allele, on the other hand, conferred a protective effect for PD-SSNHL (AA vs. TT: p = 0.012, OR = 0.25, 95% CI = 0.08-0.74; A vs. T: p = 0.001, OR = 0.36, 95% CI = 0.21-0.61). However, there was no statistically significant difference in genotype or allele frequency between the two groups (HF-SSNHL and AF-SSNHL) and the control group. Conclusion: The αENaC p. Ala663Thr gene polymorphism plays different roles in different types of SSNHL.

Computational Pipeline for Next-Generation Sequencing (NGS) Studies in Genetics of NASH

High-throughput sequencing (HTS) technologies have contributed to expand current knowledge of the biology of complex diseases, including nonalcoholic fatty liver disease (NAFLD). Genome-wide association studies, whole exome sequencing, and sequencing of entire genes are used to identify variants and/or mutations that predispose to the disease pathogenesis. Here, we present a tutorial that may guide readers to manage high volume of genetics data in the context of Next-Generation Sequencing (NGS) studies.

High-Throughput Multimodal Single-Cell Targeted DNA and Surface Protein Analysis Using the Mission Bio Tapestri Platform

An important aspect of understanding cancer biology is to connect the diverse repertoire of genotype-to-phenotype displays in individual specimens and ultimately resolve disease course outcome through informative datasets. A focus of cancer genomics has strived to provide predictive capabilities using genomic information to further inform therapeutic strategies. The advent of single-cell sequencing and analysis now provides a route to decipher high-resolution genomic diversity in individual samples and facilitate detailed understanding of clonal evolution in clinical research settings. In addition to generating high-throughput single-cell genomic SNV and CNV data, this protocol describes a new analytical ability that adds a second dimension which provides for interrogation of surface protein marker expression. The first immediate application of this technology is quite suitable to heme cancer cell studies. This multimodal approach allows for correlation of diverse genomic signatures to key phenotypic biomarkers such as immunophenotypes in leukemic diseases.

Whole-exome sequencing reveals novel candidate single nucleotide variations for preventing adverse effects of levonorgestrel implantation

Aim: To identify novel genes associated with adverse effects of levonorgestrel (LNG) implants based on comparative whole-exome sequencing. Materials & methods: A cohort comprising 104 participants, including 52 controls and 52 women with LNG-related adverse effects, was recruited. Seven cases and eight controls were selected for whole-exome sequencing. We verified 13 single nucleotide variations (SNVs) related with integrin-mediated signaling pathway and cell proliferation using the MassARRAY platform. Results: Finally, we screened 49 cases and 52 controls for analyses. Two SNVs (rs7255721 and rs1042522) were located in ADAMTS10 and TP53, respectively, and significantly different between two groups. These two SNVs lead to changes in protein structure and physicochemical parameters. Conclusion: Here, we defined two pathogenic mutations related to adverse LNG effects.