Blood Phytosterol Concentration and Genetic Variant Associations in a Sample Population

The main objective of this study was to determine plasma levels of PS and to study SNVs rs41360247, rs4245791, rs4148217, and rs11887534 of ABCG8 and the r657152 SNV at the ABO blood group locus in a sample of a population treated at our hospital, and to determine whether these SNVs are related to plasma PS concentrations. The secondary objective was to establish the variables associated with plasma PS concentrations in adults. Participants completed a dietary habit questionnaire and a blood sample was collected to obtain the following variables: campesterol, sitosterol, sitostanol, lanosterol, stigmasterol, biochemical parameters, and the SNVs. In addition, biometric and demographic variables were also recorded. In the generalized linear model, cholesterol and age were positively associated with total PS levels, while BMI was negatively related. For rs4245791, homozygous T allele individuals showed a significantly lower campesterol concentration compared with C homozygotes, and the GG alleles of rs657152 had the lowest levels of campesterol compared with the other alleles of the SNV. Conclusions: The screening of certain SNVs could help prevent the increase in plasma PS and maybe PNALD in some patients. However, further studies on the determinants of plasma phytosterol concentrations are needed.

Genetic characteristics involved in COVID-19 severity. The CARGENCORS case-control study and meta-analysis

Pre-existing coronary artery disease (CAD), and thrombotic, inflammatory, or virus infectivity response phenomena have been associated with COVID-19 disease severity. However, the association of candidate single nucleotide variants (SNVs) related to mechanisms of COVID-19 complications has been seldom analysed. Our aim was to test and validate the effect of candidate SNVs on COVID-19 severity. CARGENCORS (CARdiovascular GENetic risk score for Risk Stratification of patients positive for SARS-CoV-2 [COVID-19] virus) is an age- and sex-matched case-control study with 818 COVID-19 cases hospitalized with hypoxemia, and 1636 controls with COVID-19 treated at home. The association between severity and SNVs related to CAD (n = 32), inflammation (n = 19), thrombosis (n = 14), virus infectivity (n = 11), and two published to be related to COVID-19 severity was tested with adjusted logistic regression models. Two external independent cohorts were used for meta-analysis (SCOURGE and UK Biobank). After adjustment for potential confounders, 14 new SNVs were associated with COVID-19 severity in the CARGENCORS Study. These SNVs were related to CAD (n = 10), thrombosis (n = 2), and inflammation (n = 2). We also confirmed eight SNVs previously related to severe COVID-19 and virus infectivity. The meta-analysis showed five SNVs associated with severe COVID-19 in adjusted analyses (rs11385942, rs1561198, rs6632704, rs6629110, and rs12329760). We identified 14 novel SNVs and confirmed eight previously related to COVID-19 severity in the CARGENCORS data. In the meta-analysis, five SNVs were significantly associated to COVID-19 severity, one of them previously related to CAD.

Selected cachaça yeast strains share a genomic profile related to traits relevant to industrial fermentation processes

The isolation and selection of yeast strains to improve the quality of the cachaça-Brazilian Spirit-have been studied in our research group. Our strategy considers Saccharomyces cerevisiae as the predominant species involved in sugarcane juice fermentation and the presence of different stressors (osmolarity, temperature, ethanol content, and competition with other microorganisms). It also considers producing balanced concentrations of volatile compounds (higher alcohols and acetate and/or ethyl esters), flocculation capacity, and ethanol production. Since the genetic bases behind these traits of interest are not fully established, the whole genome sequencing of 11 different Saccharomyces cerevisiae strains isolated and selected from different places was analyzed to identify the presence of a specific genetic variation common to cachaça yeast strains. We have identified 20,128 single-nucleotide variants shared by all genomes. Of these shared variants, 37 were new variants (being six missenses), and 4,451 were identified as missenses. We performed a detailed functional annotation (using enrichment analysis, protein-protein interaction network analysis, and database and in-depth literature searches) of these new and missense variants. Many genes carrying these variations were involved in the phenotypes of flocculation, tolerance to fermentative stresses, and production of volatile compounds and ethanol. These results demonstrate the existence of a genetic profile shared by the 11 strains under study that could be associated with the applied selective strategy. Thus, this study points out genes and variants that may be used as molecular markers for selecting strains well suited to the fermentation process, including genetic improvement by genome editing, ultimately producing high-quality beverages and adding value.IMPORTANCEThis work demonstrates the existence of new genetic markers related to different phenotypes used to select yeast strains and mutations in genes directly involved in producing flavoring compounds and ethanol, and others related to flocculation and stress resistance.

MRAS in coronary artery disease-Unchartered territory

Genome-wide association studies (GWAS) have identified coronary artery disease (CAD) susceptibility locus on chromosome 3q22.3. This locus contains a cluster of several genes that includes muscle rat sarcoma virus (MRAS). Common MRAS variants are also associated with CAD causing risk factors such as hypertension, dyslipidemia, obesity, and type II diabetes. The MRAS gene is an oncogene that encodes a membrane-bound small GTPase. It is involved in a variety of signaling pathways, regulating cell differentiation and cell survival (mitogen-activated protein kinase [MAPK]/extracellular signal-regulated kinase and phosphatidylinositol 3-kinase) as well as acute phase response signaling (tumor necrosis factor [TNF] and interleukin 6 [IL6] signaling). In this review, we will summarize the role of genetic MRAS variants in the etiology of CAD and its comorbidities with the focus on tissue distribution of MRAS isoforms, cell type/tissue specificity, and mode of action of single nucleotide variants in MRAS associated complex traits. Finally, we postulate that CAD risk variants in the MRAS locus are specific to smooth muscle cells and lead to higher levels of MRAS, particularly in arterial and cardiac tissue, resulting in MAPK-dependent tissue hypertrophy or hyperplasia.

Association Study between Antioxidant Nutrient Intake and Low Bone Mineral Density with Oxidative Stress-Single Nucleotide Variants: GPX1 (rs1050450 and rs17650792), SOD2 (rs4880) and CAT (rs769217) in Mexican Women

Oxidative stress is essential in developing multiple bone metabolism diseases, including osteoporosis. Single-nucleotide variants (SNVs) have been associated with oxidative stress, promoting an imbalance between the production of reactive oxygen species and the ability to neutralize them, and it has been reported that antioxidant nutrient intake can influence bone mineral density (BMD). This work reports the association between oxidative stress-related SNVs (GPX1-rs1050450, rs17650792, SOD2-rs4880, and CAT-rs769217), BMD, and antioxidant nutrient intake. The study included 1269 Mexican women from the Health Workers Cohort Study. Genotyping was performed using predesigned TaqMan assays. Dietary data were collected using a 116-item semi-quantitative food frequency questionnaire. A dietary antioxidant quality score (DAQS) was used to estimate antioxidant-nutrient intake. Association analysis was estimated via linear, logistic, or quantile regression models. The results showed an association of the rs1050450-A and rs17650792-A alleles with femoral neck BMD (p = 0.038 and p = 0.017, respectively) and the SNV rs4880-A allele with total hip BMD (p = 0.026) in respondents aged 45 years or older. In addition, antioxidant-nutrient intake was associated with the rs4880-GG genotype, being significant for fiber (p = 0.007), riboflavin (p = 0.005), vitamin B6 (p = 0.034), and vitamin D (p = 0.002). The study showed an association between oxidative stress-related SNVs, BMD, and antioxidant-nutrient intake in Mexican women. Therefore, treatments for low BMD could be developed based on antioxidant supplementation.

Genomic architecture of autism spectrum disorder in Qatar: The BARAKA-Qatar Study

BACKGROUND

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by impaired social and communication skills, restricted interests, and repetitive behaviors. The prevalence of ASD among children in Qatar was recently estimated to be 1.1%, though the genetic architecture underlying ASD both in Qatar and the greater Middle East has been largely unexplored. Here, we describe the first genomic data release from the BARAKA-Qatar Study-a nationwide program building a broadly consented biorepository of individuals with ASD and their families available for sample and data sharing and multi-omics research.

METHODS

In this first release, we present a comprehensive analysis of whole-genome sequencing (WGS) data of the first 100 families (372 individuals), investigating the genetic architecture, including single-nucleotide variants (SNVs), copy number variants (CNVs), tandem repeat expansions (TREs), as well as mitochondrial DNA variants (mtDNA) segregating with ASD in local families.

RESULTS

Overall, we identify potentially pathogenic variants in known genes or regions in 27 out of 100 families (27%), of which 11 variants (40.7%) were classified as pathogenic or likely-pathogenic based on American College of Medical Genetics (ACMG) guidelines. Dominant variants, including de novo and inherited, contributed to 15 (55.6%) of these families, consisting of SNVs/indels (66.7%), CNVs (13.3%), TREs (13.3%), and mtDNA variants (6.7%). Moreover, homozygous variants were found in 7 families (25.9%), with a sixfold increase in homozygous burden in consanguineous versus non-consanguineous families (13.6% and 1.8%, respectively). Furthermore, 28 novel ASD candidate genes were identified in 20 families, 23 of which had recurrent hits in MSSNG and SSC cohorts.

CONCLUSIONS

This study illustrates the value of ASD studies in under-represented populations and the importance of WGS as a comprehensive tool for establishing a molecular diagnosis for families with ASD. Moreover, it uncovers a significant role for recessive variation in ASD architecture in consanguineous settings and provides a unique resource of Middle Eastern genomes for future research to the global ASD community.

Detection of Cancer-Associated Gene Mutations in Urinary Cell-Free DNA among Prostate Cancer Patients in South Africa

Prostate cancer (PCa) is the most common cause of cancer death among African men. The presence of tumor-specific variations in cell-free DNA (cfDNA), such as mutations, microsatellite instability, and DNA methylation, has been explored as a source of biomarkers for cancer diagnosis. In this study, we investigated the diagnostic role of cfDNA among South African PCa patients. We performed whole exome sequencing (WES) of urinary cfDNA. We identified a novel panel of 31 significantly deregulated somatic mutated genes between PCa and benign prostatic hyperplasia (BPH). Additionally, we performed whole-genome sequencing (WGS) on matching PCa and normal prostate tissue in an independent PCa cohort from South Africa. Our results suggest that the mutations are of germline origin as they were also found in the normal prostate tissue. In conclusion, our study contributes to the knowledge of cfDNA as a biomarker for diagnosing PCa in the South African population.

An in-silico study to determine susceptibility to cancer by evaluating the coding and non-coding non-synonymous single nucleotide variants in the SOCS3 gene

Single Nucleotide Variant (SNVs) affect gene expression as well as protein structure and activity, leading to reduced signaling capabilities and ultimately, increasing cancer risk. SOCS3 (suppressor of cytokine signaling 3), a critical tumor suppressor providing a substantial part in the feedback loop of the JAK/STAT pathway, is abnormally suppressed in various cancer. This study aims to screen non-coding and potentially deleterious coding SNVs in the SOCS3 gene. We performed six programs: PredictSNP1.0 (predicting Deleterious nsSNVs), ConSurf (analyzing sequence conservation), ModPred (analyzing SNVS in PTMs sites), I-Mutant and MUpro (to analyze SNVs effecting protein stability), and molecular docking and molecular dynamics (MD) (to assess the consequences of SOCS3 genetic variations on JAK interactions) for coding regions and three programs (UTRSite, SNP2TFBS, miRNA SNP) (to analyze SNVs effecting the gene expression) in non-coding regions, respectively. After screening 2786 SOCS3 SNVs, we found 10 SNVs, as well as 49 SNPs that change the function of non-coding areas. Out of 10 selected nsSNVs, 3 SNVs (W48R, R71C, N198S) predicted to be the most damaging by all the software programs, as well as one nsSNV (R194W) could be highly deleterious from Molecular Docking analysis combined with MD Simulations. Our findings propose a procedure for studying the structure-related consequences of SNVs on protein function in the future.Communicated by Ramaswamy H. Sarma.

CRISPR-Mediated Base Editing: Promises and Challenges for a Viable Oncotherapy Strategy

Base editing technology, developed from the CRISPR/Cas9 system, is able to efficiently implement single-base substitutions at specific DNA or RNA sites without generating double-strand breaks with precision and efficiency. Point mutations account for 58% of disease-causing genetic mutations in humans, and single nucleotide variants are an important cause of tumorigenesis, and the advent of base editors offers new hope for the study or treatment of such diseases. Although it has some limitations, base editors have been continuously improved in terms of editing efficiency, specificity, and product purity since their development. In this review, we examine the main base editing technologies and discuss their applications and prospects in tumor research and therapy, as well as elaborate on their mode of delivery.

Could exposure to spaceflight cause mutations in genes that affect male fertility?

Recently, a study reported that upon analyzing blood samples from 14 astronauts that flew Space Transportation System missions between 1998 and 2001, 34 somatic nonsynonymous single nucleotide variants were detected in 17 CH-driver genes. Of interest is that the cohort consisted of relatively young astronauts, 85% of which were males of reproductive age. Having investigated the genes with nonsynonymous substitutes from the literature, it was found that twelve of these 17 genes appear to play essential roles in male reproduction. Changes in telomere length and gene regulation were also reported in another study conducted on an astronaut during a long duration stay on the International Space Station. Realizing the impact of spaceflight on gene sequence with potential influence on male fertility, it is important that more studies are conducted in this field. Specifically, in light of ultimately colonizing space, multi-generational survival is crucial and strategies to mitigate or counteract such effects should be explored.

Genetic analysis and allele-specific expression of SMAD7 3'UTR variants in human colorectal cancer reveal a novel somatic variant exhibiting allelic imbalance

BACKGROUND

Considering the impact of SMAD7 deregulation in colorectal cancer (CRC) progression and the significance of single nucleotide variant (SNV)-mediated disruptions of microRNA (miRNA)-dependent regulation for cancer susceptibility, our study aimed to analyze genetic variation in the SMAD7 3' untranslated region ( 3'UTR) in CRC, measure differences in allelic mRNA expression, and evaluate its interference with miRNA-mediated post-transcriptional regulation.

PATIENTS AND METHODS

This study included 80 patients with different CRC stages and six human colon cancer cell lines of various histological origins. SMAD7 3'UTR was analyzed by direct sequencing, followed by the relative quantification of differential allelic expression of detected variants by allele-specific qRT-PCR. In silico tools were employed for predictions of regulatory consequences of detected variants.

RESULTS

A total of four different SNVs in one cell line and nine patients were found, among which were a novel somatic point variant and three already known germline variants (rs16950113, rs1050799536, and rs1043778717). All evaluated SNVs exhibited variable extents of allelic imbalance in expression. In silico analysis predicted significant effects of SNVs on miRNA binding efficiency, with each SNV disrupting existing and creating new target sites for one or more miRNAs.

CONCLUSION

Imbalance observed in the expression of SNV alleles altering miRNA binding suggests that all investigated SNVs are potential contributing factors impacting SMAD7 expression regulation in CRC that further studies should investigate.

Two single nucleotide variants in the miR-23a promoter affect granulosa cell apoptosis

microRNAs (miRNAs) are well known to be important in mammalian female fertility. However, the genetic regulation of miRNAs associated with female fertility remains largely unknown. Here, we report that two single-nucleotide variants (SNVs) in the miR-23a promoter strongly influence miR-23a transcription and function in granulosa cell (GC) apoptosis. Two novel SNVs, g.-283G>C and g.-271C>T, were detected in the porcine miR-23a promoter by pooled-DNA sequencing. Furthermore, SNVs in the promoter region influenced miR-23a transcription in porcine GCs by altering its promoter activity. Functionally, SNVs in the promoter strongly influenced miR-23a regulation of early apoptosis in porcine GCs cultured in vitro. In addition, a preliminary association analysis showed that the combined genotypes of the two SNVs, rather than a single SNV, were tentatively associated with sow fertility traits in a Large White population. Overall, our findings suggest that the SNVs g.-283G>C and g.-271C>T in the miR-23a promoter are causal variants affecting GC apoptosis and miR-23a may be a potential small-molecule nonhormonal drug for regulating female fertility.

Influence of Genetic Variations in miRNA and Genes Encoding Proteins in the miRNA Synthesis Complex on Toxicity of the Treatment of Pediatric B-Cell ALL in the Brazilian Amazon

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer in the world. Single nucleotide variants (SNVs) in miRNA and genes encoding proteins of the miRNA synthesis complex (SC) may affect the processing of drugs used in the treatment of ALL, resulting in treatment-related toxicities (TRTs). We investigated the role of 25 SNVs in microRNA genes and genes encoding proteins of the miRNA SC, in 77 patients treated for ALL-B from the Brazilian Amazon. The 25 SNVs were investigated using the TaqMan^® OpenArray™ Genotyping System. SNVs rs2292832 (MIR149), rs2043556 (MIR605), and rs10505168 (MIR2053) were associated with an increased risk of developing Neurological Toxicity, while rs2505901 (MIR938) was associated with protection from this toxicity. MIR2053 (rs10505168) and MIR323B (rs56103835) were associated with protection from gastrointestinal toxicity, while DROSHA (rs639174) increased the risk of development. The rs2043556 (MIR605) variant was related to protection from infectious toxicity. SNVs rs12904 (MIR200C), rs3746444 (MIR499A), and rs10739971 (MIRLET7A1) were associated with a lower risk for severe hematologic toxicity during ALL treatment. These findings reveal the potential for the use of these genetic variants to understand the development of toxicities related to the treatment of ALL in patients from the Brazilian Amazon region.

Analysis of Doberman Pinscher and Toy Poodle samples with targeted next-generation sequencing

Next-generation sequencing (NGS) technologies have enabled the identification of many causal variants of genetic disorders, the development of parentage tests and the analysis of multiple traits in domestic animals. In this study, we evaluated the performance of a Canine Targeted Genotyping-by-Sequencing (GBS) custom panel (Thermo Fisher Scientific, Waltham, Ma, USA) in a cohort of 95 dog DNA samples, comprising 76 Doberman Pinschers and 19 Toy Poodles from Argentina. The used panel included 383 targets (228 parentage SNVs, 137 genetic disorder markers and 18 trait markers). While paternity analysis showed correct duo (97.4%; LOD > 2.98E⁺¹³) and trio (100%; LOD > 2.20E⁺¹⁵) parentage assignment, the panel resulted still insufficient for excluding close relatives in inbred populations. In this sense, close relatives were wrongly assigned as parents in 12.6% of duos and 0.3% of trios. We detected 17 polymorphic markers (genetic disorders, n = 4; hair type, n = 3; coat color, n = 10) and estimated their allele frequencies in the studied breeds. The accuracy of targeted GBS results were evaluated for three markers that were associated with Progressive rod-cone degeneration, von Willebrand disease type 1 and dilated cardiomyopathy by pyrosequencing and Sanger sequencing genotyping, showing 94-100% concordance among assays. The targeted GBS custom panel resulted cost-effective strategy to study the prevalence of genetic disorders and traits in a large number of samples and to analyze genetic interactions between previously reported variants. Once assays based on AgriSeq technology were standardized, their uses are a good strategy for large-scale routine genetic evaluation of animal populations.

CVID-Associated B Cell Activating Factor Receptor Variants Change Receptor Oligomerization, Ligand Binding, and Signaling Responses

PURPOSE

Binding of the B cell activating factor (BAFF) to its receptor (BAFFR) activates in mature B cells many essential pro-survival functions. Null mutations in the BAFFR gene result in complete BAFFR deficiency and cause a block in B cell development at the transition from immature to mature B cells leading therefore to B lymphopenia and hypogammaglobulinemia. In addition to complete BAFFR deficiency, single nucleotide variants encoding BAFFR missense mutations were found in patients suffering from common variable immunodeficiency (CVID), autoimmunity, or B cell lymphomas. As it remained unclear to which extent such variants disturb the activity of BAFFR, we performed genetic association studies and developed a cellular system that allows the unbiased analysis of BAFFR variants regarding oligomerization, signaling, and ectodomain shedding.

METHODS

In addition to genetic association studies, the BAFFR variants P21R, A52T, G64V, DUP92-95, P146S, and H159Y were expressed by lentiviral gene transfer in DG-75 Burkitt's lymphoma cells and analyzed for their impacts on BAFFR function.

RESULTS

Binding of BAFF to BAFFR was affected by P21R and A52T. Spontaneous oligomerization of BAFFR was disturbed by P21R, A52T, G64V, and P146S. BAFF-dependent activation of NF-κB2 was reduced by P21R and P146S, while interactions between BAFFR and the B cell antigen receptor component CD79B and AKT phosphorylation were impaired by P21R, A52T, G64V, and DUP92-95. P21R, G64V, and DUP92-95 interfered with phosphorylation of ERK1/2, while BAFF-induced shedding of the BAFFR ectodomain was only impaired by P21R.

CONCLUSION

Although all variants change BAFFR function and have the potential to contribute as modifiers to the development of primary antibody deficiencies, autoimmunity, and lymphoma, P21R is the only variant that was found to correlate positively with CVID.

In Silico Analysis: HLA-DRB1 Gene's Variants and Their Clinical Impact

The HLA-DRB1 gene encodes a protein that is essential for the immune system. This gene is important in organ transplant rejection and acceptance, as well as multiple sclerosis, systemic lupus erythematosus, Addison's disease, rheumatoid arthritis, caries susceptibility, and Aspirin-exacerbated respiratory disease. The following Homo sapiens variants were investigated: single-nucleotide variants (SNVs), multi-nucleotide variants (MNVs), and small insertions-deletions (Indels) in the HLA-DRB1 gene via coding and untranslated regions. The current study sought to identify functional variants that could affect gene expression and protein product function/structure. ALL target variants available until April 14, 2022, were obtained from the Single Nucleotide Polymorphism database (dbSNP). Out of all the variants in the coding region, 91 nsSNVs were considered highly deleterious by seven prediction tools and instability index; 25 of them are evolutionary conserved and located in domain regions. Furthermore, 31 indels were predicted as harmful, potentially affecting a few amino acids or even the entire protein. Last, within the coding sequence (CDS), 23 stop-gain variants (SNVs/indels) were predicted as high impact. High impact refers to the assumption that the variant will have a significant (disruptive) effect on the protein, likely leading to protein truncation or loss of function. For untranslated regions, functional 55 single-nucleotide polymorphisms (SNPs), and 16 indels located within microRNA binding sites, furthermore, 10 functionally verified SNPs were predicted at transcription factor-binding sites. The findings demonstrate that employing in silico methods in biomedical research is extremely successful and has a major influence on the capacity to identify the source of genetic variation in diverse disorders. In conclusion, these previously functional identified variants could lead to gene alteration, which may directly or indirectly contribute to the occurrence of many diseases. The study's results could be an important guide in the research of potential diagnostic and therapeutic interventions that require experimental mutational validation and large-scale clinical trials.

Association between PPARγ, PPARGC1A, and PPARGC1B genetic variants and susceptibility of gastric cancer in an Eastern Chinese population

PURPOSE

Previous studies showed that peroxisome proliferator-activated receptor gamma (PPARγ) and PPARγ coactivator1 family (PPARGC1A and PPARGC1B) gene single nucleotide variants (SNVs)were strongly associated with cancer susceptibility. The purpose of this study was to investigate the association of PPARγ, PPARGC1A, and PPARGC1B variants with the risk of gastric cancer (GC).

PATIENTS AND METHODS

We performed a case-control study of 490 GC cases and 1,476 healthy controls from eastern China. PPARγ rs1801282 C > G, rs3856806 C > T, PPARGC1A rs2970847 C > T, rs8192678 C > T and PPARGC1B rs7732671 G > C, rs17572019 G > A SNVs were selected to investigate the association between these SNVs and GC susceptibility. Genotypes of the SNVs were assessed by multiplex fluorescent PCR using a custom-by-design 48-Plex SNPscan^tm Kit.

RESULTS

The PPARγ rs1801282 SNV was associated with a decreased risk for GC (GC vs. CC: odds ratio (OR) = 0.62, 95% confidence interval (95%CI) = 0.42-0.93, adjusted P = 0.019; GC + GG vs. GG: OR = 0.63 95%CI = 0.42-0.93, adjusted P = 0.019; respectively). In addition, stratified analysis revealed that the PPARγ rs1801282 SNV was correlated with the risk of GC in subgroups of age ≥ 61, no smoking, and no alcohol consuming. We also confirmed that the PPARγ rs3856806 C > T SNV promoted the risk of GC in women. The PPARGC1A rs8192678 TT genotype decreased the susceptibility of GC in men. The PPARGC1A rs2970847 C > T SNV decreased the susceptibility of GC in the subgroup of BMI ≥ 24 kg/m². The PPARGC1B rs7732671 G > C and rs17572019 G > A SNVs promoted the risk of GC in the subgroup of BMI ≥ 24 kg/m².

CONCLUSION

This study indicates that the PPARγ, PPARGC1A, and PPARGC1B SNVs may be associated with the susceptibility of GC in eastern Chinese population. Future studies with larger populations, detailed H. pylori infection status for subgroup analysis, and functional study are needed to further clarify the relationship between these SNVs and GC risk.

Genome stability during serial subculturing in hyperepidemic multidrug-resistant Klebsiella pneumoniae and Escherichia coli

OBJECTIVES

Core-genome single nucleotide variant (cgSNV) analysis represents a powerful tool for epidemiological investigations of multidrug-resistant (MDR) bacteria. However, cgSNV thresholds to confirm whether isolates are the same clone are not formally defined.

METHODS

We implemented hybrid whole-genome sequencing to study the genomic changes of four MDR isolates belonging to hyperepidemic sequence types (STs) during 20 propagation steps (T20) on MacConkey and CHROMID^(R) ESBL plates. The following strains were analyzed: Klebsiella pneumoniae AE-2247421 (OXA-48/NDM-1-producing, ST101), K. pneumoniae MCL-2017-2 (CTX-M-15-producing, ST307), Escherichia coli Ec-042 (OXA-181-producing, ST410), and E. coli Ec-050 (NDM-5-producing, ST167). The genome assembly at T5 and T20 was compared to that at time point zero (T0) and to two reference genomes.

RESULTS

At T20, AE-2247421 lost the IncL bla_OXA-48-carrying plasmid when grown on CHROMID^(R) ESBL plates, while a large fragment encompassing bla_NDM-1 was lost from its IncC plasmid when grown on both plates. In contrast, no structural changes were noted for the other three strains. Regarding the cgSNVs, the following results were obtained at T5 and T20 (ranges considering the different agar plates and reference genomes): AE-2247421 (1-8 and 2-12 cgSNVs), MCL-2017-2 (both 1-2 cgSNVs), Ec-042 (both 0 cgSNVs), and Ec-050 (0-6 and 0-9 cgSNVs).

CONCLUSION

We showed that structural changes and accumulation of cgSNVs can occur in few propagation steps under laboratory conditions. These changes might also arise in the clinical context in a short time, especially under antibiotics treatment. This phenomenon should be carefully considered because it might affect the final interpretation of epidemiological genomic analyses.

Whole-exome sequencing reveals novel variants associated with abnormal uterine bleeding caused by copper intrauterine device

Aim: This study aimed to explore the genetic risk factors and validate variants of abnormal uterine bleeding after copper intrauterine device insertion. Methods: Whole-exome sequencing was performed and several variants were validated by Sequenom MassARRAY. Results: Eight variants showed potential clinical damage according to American College of Medical Genetics and Genomics criteria. By combined analysis of screening and validation, NFASC RS2802808 C>G p.Ile971Met (P_allele = 0.009 and P_genotype = 0.027) and PIGR RS2275531 C>T p.Gly365Ser (P_allele = 0.009 and P_genotype = 0.013) variants were identified as significantly associated with abnormal uterine bleeding with a false discovery rate <0.05. NFASC and PIGR may play a role in abnormal uterine bleeding by regulating coagulation fibrinolysis and endometrial epithelium inflammation functions. Conclusion: These findings provide a genetic basis for clinical individualization and precision of intrauterine device implantation.

Landscapes of synchronous multiple primary cancers detected by next-generation sequencing

An increase in the detection rate of multiple primary cancers has been accompanied with declining cancer death rates over the past few decades. However, synchronous multiple primary tumors have gradually increased, and the molecular mechanisms involved in the synchronous occurrence of multiple primary cancers of different origins are unclear. To investigate these mechanisms, we sequenced cancer tissues by FoundationOne CDx. Data were annotated with annovar, and we then performed pathway enrichment analysis. A total of 109 genes that were mutated in all samples were clustered into different diseases, biological processes, and molecular functions. GO and KEGG analyses indicated that the P53 and PKB signaling pathways may be relevant to the occurrence of synchronous multiple primary cancers. In summary, patients with a concordance of mutations in pathogenetic genes may have a higher risk of developing a second cancer. Our research may provide a basis for the development of individualized treatments for synchronous multiple primary cancers.

A comprehensive analysis of the mutational landscape of the newly emerging Omicron (B.1.1.529) variant and comparison of mutations with VOCs and VOIs

The Omicron variant is spreading rapidly throughout several countries. Thus, we comprehensively analyzed Omicron's mutational landscape and compared mutations with VOC/VOI. We analyzed SNVs throughout the genome, and AA variants (NSP and SP) in VOC/VOI, including Omicron. We generated heat maps to illustrate the AA variants with high mutation prevalence (> 75% frequency) of Omicron, which demonstrated eight mutations with > 90% prevalence in ORF1a and 29 mutations with > 75% prevalence in S-glycoprotein. A scatter plot for Omicron and VOC/VOI's cluster evaluation was computed. We performed a risk analysis of the antibody-binding risk among four mutations (L452, F490, P681, D614) and observed three mutations (L452R, F490S, D614G) destabilized antibody interactions. Our comparative study evaluated the properties of 28 emerging mutations of the S-glycoprotein of Omicron, and the ΔΔG values. Our results showed K417N with minimum and Q954H with maximum ΔΔG value. Furthermore, six important RBD mutations (G339D, S371L, N440K, G446S, T478K, Q498R) were chosen for comprehensive analysis for stabilizing/destabilizing properties and molecular flexibility. The G339D, S371L, N440K, and T478K were noted as stable mutations with 0.019 kcal/mol, 0.127 kcal/mol, 0.064 kcal/mol, and 1.009 kcal/mol. While, G446S and Q498R mutations showed destabilizing results. Simultaneously, among six RBD mutations, G339D, G446S, and Q498R mutations increased the molecular flexibility of S-glycoprotein. This study depicts the comparative mutational pattern of Omicron and other VOC/VOI, which will help researchers to design and deploy novel vaccines and therapeutic antibodies to fight against VOC/VOI, including Omicron.

Human ACE2 Genetic Polymorphism Affecting SARS-CoV and SARS-CoV-2 Entry into Cells

Severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 have a single envelope glycoprotein (S protein) that binds to human angiotensin-converting enzyme 2 (ACE2) on the host cell membrane. Previous mutational scanning studies have suggested that some substitutions corresponding to single nucleotide variants (SNVs) in human ACE2 affect the binding affinity to the receptor binding domain (RBD) of the SARS-CoV-2 S protein. However, the importance of these substitutions in actual virus infection is still unclear. In this study, we investigated the effects of the reported ACE2 SNV substitutions on the entry of SARS-CoV and SARS-CoV-2 into cells, using vesicular stomatitis Indiana virus (VSIV) pseudotyped with S proteins of these coronaviruses (CoVs). HEK293T cells transfected with plasmids expressing ACE2 having each SNV substitution were infected with the pseudotyped VSIVs and relative infectivities were determined compared to the cells expressing wild-type ACE2. We found that some of the SNV substitutions positively or negatively affected the infectivities of the pseudotyped viruses. Particularly, the H505R substitution significantly enhanced the infection with the pseudotyped VSIVs, including those having the substitutions found in the S protein RBD of SARS-CoV-2 variants of concern. Our findings suggest that human ACE2 SNVs may potentially affect cell susceptibilities to SARS-CoV and SARS-CoV-2. IMPORTANCE SARS-CoV and SARS-CoV-2 are known to cause severe pneumonia in humans. The S protein of these CoVs binds to the ACE2 molecule on the plasma membrane and mediates virus entry into cells. The interaction between the S protein and ACE2 is thought to be important for host susceptibility to these CoVs. Although previous studies suggested that some SNV substitutions in ACE2 might affect the binding to the S protein, it remains elusive whether these SNV substitutions actually alter the efficiency of the entry of SARS CoVs into cells. We analyzed the impact of the ACE2 SNVs on the cellular entry of SARS CoVs using pseudotyped VSIVs having the S protein on the viral surface. We found that some of the SNV substitutions positively or negatively affected the infectivities of the viruses. Our data support the notion that genetic polymorphisms of ACE2 may potentially influence cell susceptibilities to SARS CoVs.

Novel Screening System of Virulent Strains for the Establishment of a Mycobacterium avium Complex Lung Disease Mouse Model Using Whole-Genome Sequencing

The establishment of animal models reflecting human Mycobacterium avium complex (MAC) lung disease (LD) pathology has the potential to expand our understanding of the disease pathophysiology. However, inducing sustained infection in immunocompetent mice is difficult since MAC generally shows less virulence and higher genetic variability than M. tuberculosis. To overcome this hurdle, we developed a screening system for identifying virulent MAC strains using whole-genome sequencing (WGS). We obtained nine clinical strains from Mycobacterium avium complex lung disease (MAC-LD) patients and divided them into two groups to make the mixed strain inocula for infection. Intranasal infection with the strain mixture of both groups in BALB/c mice resulted in progressive infection and extensive granuloma formation in the lungs, suggesting the existence of highly pathogenic strains in each group. We hypothesized that the change in the abundance of strain-specific single-nucleotide variants (SNVs) reflects the change in bacterial number of each strain in infected lungs. Based on this hypothesis, we quantified individual strain-specific SNVs in bacterial DNA from infected lungs. Specific SNVs for four strains were detected, suggesting the pathogenicity of these four strains. Consistent with these results, individual infection with these four strains induced a high lung bacterial burden, forming extensive peribronchial granuloma, while the other strains showed a decreased lung bacterial burden. The current method combining mixed infection and WGS accurately identified virulent strains that induced sustained infection in mice. This method will contribute to the establishment of mouse models that reflect human MAC-LD and lead to antimycobacterial drug testing. IMPORTANCE To promote research on Mycobacterium avium complex (MAC) pathogenicity, animal models reflecting human progressive MAC lung disease (MAC-LD) are needed. Because there is high genetic and virulence diversity among clinical MAC strains, choosing a suitable strain is an important process for developing a mouse model. In this study, we developed a screening system for virulent strains in mice by combining mixed infection and whole-genome sequencing analysis. This approach is designed on the hypothesis that in vivo virulence of MAC strains can be examined simultaneously by comparing changes in the abundance of strain-specific single-nucleotide variants in the mouse lungs after infection with mixed strains. The identified strains were shown to induce high bacterial burdens and cause extensive peribronchial granuloma resembling the pulmonary pathology of human MAC-LD. The current method will help researchers develop mouse models that reflect human MAC-LD and will lead to further investigation of MAC pathogenicity.

Genomic Profile in a Non-Seminoma Testicular Germ-Cell Tumor Cohort Reveals a Potential Biomarker of Sensitivity to Platinum-Based Therapy

Despite having a favorable response to platinum-based chemotherapies, ~15% of Testicular Germ-Cell Tumor (TGCT) patients are platinum-resistant. Mortality rates among Latin American countries have remained constant over time, which makes the study of this population of particular interest. To gain insight into this phenomenon, we conducted whole-exome sequencing, microarray-based comparative genomic hybridization, and copy number analysis of 32 tumors from a Mexican cohort, of which 18 were platinum-sensitive and 14 were platinum-resistant. We incorporated analyses of mutational burden, driver mutations, and SNV and CNV signatures. DNA breakpoints in genes were also investigated and might represent an interesting research opportunity. We observed that sensitivity to chemotherapy does not seem to be explained by any of the mutations detected. Instead, we uncovered CNVs, particularly amplifications on segment 2q11.1 as a novel variant with chemosensitivity biomarker potential. Our data shed light into understanding platinum resistance in a Latin-origin population.

A Comprehensive Survey of Genomic Mutations in Breast Cancer Reveals Recurrent Neoantigens as Potential Therapeutic Targets

Neoantigens are mutated antigens specifically generated by cancer cells but absent in normal cells. With high specificity and immunogenicity, neoantigens are considered as an ideal target for immunotherapy. This study was aimed to investigate the signature of neoantigens in breast cancer. Somatic mutations, including SNVs and indels, were obtained from cBioPortal of 5991 breast cancer patients. 738 non-silent somatic variants present in at least 3 patients for neoantigen prediction were selected. PIK3CA (38%), the highly mutated gene in breast cancer, could produce the highest number of neoantigens per gene. Some pan-cancer hotspot mutations, such as PIK3CA E545K (6.93%), could be recognized by at least one HLA molecule. Since there are more SNVs than indels in breast cancer, SNVs are the major source of neoantigens. Patients with hormone receptor-positive or HER2 negative are more competent to produce neoantigens. Age, but not the clinical stage, is a significant contributory factor of neoantigen production. We believe a detailed description of breast cancer neoantigen signatures could contribute to neoantigen-based immunotherapy development.

TERT Genotype Polymorphism: A Glance of Change Egyptian MDS Outcomes

Background:

Myelodysplastic Syndromes (MDS)are clonal hematologic disorders characterized by genetic instability and ineffective hematopoiesis associated with telomere dysfunction. We aimed at investigating the association between the rs2242652 single nucleotide variant of the TERT gene and susceptibility for MDS, as well as its prognostic impact and relation to disease phenotype.

Methods:

Genotyping analysis was carried on 100 MDS patients recruited at Mansoura Oncology center, in addition to 100 healthy subjects for detection of rs2242652 variant of TERT gene on chromosome 5 by real time PCR following the protocol of Custom TaqMan^® SNP Genotyping.

Results:

The rs2242652 TERT genetic polymorphism was associated with an increased risk of MDS (odds ratios 2.6 for genotype GA, 6.4 for genotype AA). The majority of AA homozygous mutant variant were associated pancytopenia (88%), poor risk cytogenetics (92%) and High/very high IPSS-R score (88%). At the end of follow-up (median 30 months), 14% of the cases transformed to secondary AML. The rate of leukemic transformation was significantly associated with the mutant AA genotype (93% of transformed cases, 52% of AA genotype cases; P< 0.0001). Survival outcome was inferior in AA mutant genotype (median 14 months, 95% CI: 12-16 months) to the GA genotype (median 30 months, 95% CI: 26-33 months) and those of the GG genotype (median not reached), P<0.001.

Conclusion:

Our study shows an intriguing and previously unrecognized association between rs2242652 TERT mutation and MDS risk. The presence of rs2242652 mutation defines a subgroup of patients with aggressive disease phenotype and dismal outcome. Further research is recommended to elucidate underlying pathologic mechanisms and to define an efficient therapeutic target.

Novel MYT1 variants expose the complexity of oculo-auriculo-vertebral spectrum genetic mechanisms

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder characterized by anomalies mainly involving the structures derived from the first and second pharyngeal arches. The spectrum presents with heterogeneous clinical features and complex etiology with genetic factors not yet completely understood. To date, MYT1 is the most important gene unambiguously associated with the spectrum and with functional data confirmation. In this work, we aimed to identify new single nucleotide variants (SNVs) affecting MYT1 in a cohort of 73 Brazilian patients diagnosed with OAVS. In addition, we investigated copy number variations (CNVs) encompassing this gene or its cis-regulatory elements and compared the frequency of these events in patients versus a cohort of 455 Brazilian control individuals. A new SNV, predicted as likely deleterious, was identified in five unrelated patients with OAVS. All five patients presented hearing impairment and orbital asymmetry suggesting an association with the variant. CNVs near MYT1, located in its neighboring topologically associating domain (TAD), were found to be enriched in patients when compared to controls, indicating a possible involvement of this region with OAVS pathogenicity. Our findings highlight the genetic complexity of the spectrum that seems to involve more than one variant type and inheritance patterns.

Single nucleotide variants in immune-response genes and the tumor microenvironment composition predict progression of mantle cell lymphoma

Background

There is evidence to consider that the tumor microenvironment (TME) composition associates with antitumor immune response, and may predict the outcome of various non-Hodgkin lymphoma subtypes. However, in the case of mantle cell lymphoma (MCL), a rare and aggressive disease, there is lacking a detailed study of the TME components, as well as an integrative approach among them in patients’ samples. Also, from the genetic point of view, it is known that single nucleotide variants (SNVs) in immune-response genes are among important regulators of immunity. At present, it is uncertain whether SNVs in candidate immune-response genes and the TME composition are able to alter the prognosis in MCL.

Methods

We assessed a detailed TME composition in 88 MCL biopsies using immunohistochemistry, which was automatically analyzed by pixel counting (Aperio system). We also genotyped SNVs located in candidate immune-response genes (IL12A, IL2, IL10, TGFB1, TGFBR1, TGFBR2, IL17A, IL17F) in 95 MCL patients. We tested whether the SNVs could modulate the respective protein expression and TME composition in the tumor compartment. Finally, we proposed survival models in rituximab-treated patients, considering immunohistochemical and SNV models.

Results

High FOXP3/CD3 ratios (p = 0.001), high IL17A levels (p = 0.003) and low IL2 levels (p = 0.03) were individual immunohistochemical predictors of poorer survival. A principal component, comprising high quantities of macrophages and high Ki-67 index, also worsened outcome (p = 0.02). In the SNV model, the CC haplotype of IL10 (p < 0.01), the GG genotype of IL2 rs2069762 (p = 0.02) and the AA+AG genotypes of TGFBR2 rs3087465 (p < 0.01) were independent predictors of outcome. Finally, the GG genotype of TGFB1 rs6957 associated with lower tumor TGFβ levels (p = 0.03) and less CD163+ macrophages (p = 0.01), but did not modulate patients’ survival.

Conclusions

Our results indicate that the TME composition has relevant biological roles in MCL. In this setting, immunohistochemical detection of T-reg cells, IL17A and IL2, coupled with SNV genotyping in IL10, TGFBR2 and IL2, may represent novel prognostic factors in this disease, following future validations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-07891-9.

Rescue of Senecavirus A to uncover mutation profiles of its progenies during 80 serial passages in vitro

Senecavirus A (SVA), also known as Seneca Valley virus, belongs to the genus Senecavirus in the family Picornaviridae. In this study, a China SVA isolate (CH-LX-01-2016) was rescued from its cDNA clone, and then identified by RT-PCR, indirect immunofluorescence assay and mass spectrometry. The rescued SVA could separately induce typical plaque formations and cytopathic effects in cell monolayers. In order to uncover its evolutionary dynamics, the SVA was subjected to eighty serial passages in vitro. Its progenies per ten passages were analyzed by next-generation sequencing (NGS). The NGS analyses showed that neither sequence-deleting nor -inserting phenotype was detectable in eight progenies, within which a total of forty-one intra-host single-nucleotide variations (SNVs) arose with passaging. Almost all SNVs were identified as the single-nucleotide polymorphism with mixture of two nucleotides. SNVs led to eighteen nonsynonymous mutations, out of which sixteen could directly reflect their own frequencies of amino acid mutation, due to only one SNV occurring in their individual codons. Compared with its parental virus without passaging, the passage-80 SVA progeny had formed a viral quasispecies, as evidenced by a total of twenty-eight SNVs identified in it.

microRNAs deregulation in head and neck squamous cell carcinoma

Head and neck (HN) squamous cell carcinoma (SCC) is the eighth most common human cancer worldwide. Besides tobacco and alcohol consumption, genetic and epigenetic alterations play an important role in HNSCC occurrence and progression. microRNAs (miRNAs) are small noncoding RNAs that regulate cell cycle, proliferation, development, differentiation, and apoptosis by interfering in gene expression. Expression profiling of miRNAs showed that some miRNAs are upregulated or downregulated in tumor cells when compared with the normal cells. The present review focuses on the role of miRNAs deregulations in HNSCC, enrolled in risk, development, outcome, and therapy sensitivity. Moreover, the influence of single nucleotide variants in miRNAs target sites, miRNAs seed sites, and miRNAs-processing genes in HNSCC was also revised. Due to its potential for cancer diagnosis, progression, and as a therapeutic target, miRNAs may bring new perspectives in HNSCC understanding and therapy, especially for those patients with no or insufficient treatment options.

Identification of single nucleotide variants in the Moroccan population by whole-genome sequencing

BACKGROUND

Large-scale human sequencing projects have described around a hundred-million single nucleotide variants (SNVs). These studies have predominately involved individuals with European ancestry despite the fact that genetic diversity is expected to be highest in Africa where Homo sapiens evolved and has maintained a large population for the longest time. The African Genome Variation Project examined several African populations but these were all located south of the Sahara. Morocco is on the northwest coast of Africa and mostly lies north of the Sahara, which makes it very attractive for studying genetic diversity. The ancestry of present-day Moroccans is unknown and may be substantially different from Africans found South of the Sahara desert, Recent genomic data of Taforalt individuals in Eastern Morocco revealed 15,000-year-old modern humans and suggested that North African individuals may be genetically distinct from previously studied African populations.

RESULTS

We present SNVs discovered by whole genome sequencing (WGS) of three Moroccans. From a total of 5.9 million SNVs detected, over 200,000 were not identified by 1000G and were not in the extensive gnomAD database. We summarise the SNVs by genomic position, type of sequence gene context and effect on proteins encoded by the sequence. Analysis of the overall genomic information of the Moroccan individuals to individuals from 1000G supports the Moroccan population being distinct from both sub-Saharan African and European populations.

CONCLUSIONS

We conclude that Moroccan samples are genetically distinct and lie in the middle of the previously observed cline between populations of European and African ancestry. WGS of Moroccan individuals can identify a large number of novel SNVs and aid in functional characterisation of the genome.

Cryptorchidism and Testicular Tumor: Comprehensive Analysis of Common Clinical Features and Search of SNVs in the KIT and AR Genes

Allelic variants in genes implicated in the development of testicular germ cell tumor (TGCT) could be present in patients with cryptorchidism (CO). Currently; the mechanisms explaining this relationship are still unknown. In this study the common clinical features in patients with CO and TGCT and 6 variants of KIT and AR genes associated to TGCT were analyzed. Population analyzed included 328 individuals: 91 patients with CO; 79 with TGCT, 13 of them with previous CO diagnosis, and 158 healthy males. Of the 13 patients with TGCT and history of CO, one patient (7.7%) presented the heterozygous form of the variant rs121913507 and two patients (15.4%) presented homozygote genotype for the variant rs121913506 in KIT gene. Interestingly, the heterozygous form for the variant rs121913506 of KIT gene was identifying in all of 13 patients. The rs201934623, rs774171864, and rs12014709 variants of the AR gene did not show any clinical association. Our results strongly support that genetic component in CO could be conditioning for the development of TGCT. Notably, KIT gene variants might be determinants in the pathological association between TGCT and CO.

Genome-Wide DNA Alterations in X-Irradiated Human Gingiva Fibroblasts

While ionizing radiation (IR) is a powerful tool in medical diagnostics, nuclear medicine, and radiology, it also is a serious threat to the integrity of genetic material. Mutagenic effects of IR to the human genome have long been the subject of research, yet still comparatively little is known about the genome-wide effects of IR exposure on the DNA-sequence level. In this study, we employed high throughput sequencing technologies to investigate IR-induced DNA alterations in human gingiva fibroblasts (HGF) that were acutely exposed to 0.5, 2, and 10 Gy of 240 kV X-radiation followed by repair times of 16 h or 7 days before whole-genome sequencing (WGS). Our analysis of the obtained WGS datasets revealed patterns of IR-induced variant (SNV and InDel) accumulation across the genome, within chromosomes as well as around the borders of topologically associating domains (TADs). Chromosome 19 consistently accumulated the highest SNVs and InDels events. Translocations showed variable patterns but with recurrent chromosomes of origin (e.g., Chr7 and Chr16). IR-induced InDels showed a relative increase in number relative to SNVs and a characteristic signature with respect to the frequency of triplet deletions in areas without repetitive or microhomology features. Overall experimental conditions and datasets the majority of SNVs per genome had no or little predicted functional impact with a maximum of 62, showing damaging potential. A dose-dependent effect of IR was surprisingly not apparent. We also observed a significant reduction in transition/transversion (Ti/Tv) ratios for IR-dependent SNVs, which could point to a contribution of the mismatch repair (MMR) system that strongly favors the repair of transitions over transversions, to the IR-induced DNA-damage response in human cells. Taken together, our results show the presence of distinguishable characteristic patterns of IR-induced DNA-alterations on a genome-wide level and implicate DNA-repair mechanisms in the formation of these signatures.

Genetic polymorphisms in vitamin D pathway influence 25(OH)D levels and are associated with atopy and asthma

Background

Vitamin D deficiency or insufficiency, has been associated with atopy and lack of asthma control. Our objective was to investigate associations between variants in genes of vitamin D pathway with serum levels of 25-hydroxyvitamin D (25(OH)D), atopy, asthma and asthma severity in teenagers from Northeast Brazil.

Methods

This is a cross sectional study nested in a cohort population of asthma. 25(OH)D was quantified from 968 of 11–17 years old individuals by ELISA. Asthma diagnosis was obtained by using the ISAAC Phase III questionnaire. Specific IgE was determined by ImmunoCAP; genotyping was performed using the 2.5 HumanOmni Biochip from Illumina. Statistical analyses were performed in PLINK 1.07 and SPSS 22.1.

Results

After quality control, 104 Single Nucleotides Variants (SNVs) in vitamin D pathway genes, typed in 792 individuals, were included in the analysis. The allele A of rs10875694 on VDR was positively associated with atopy (OR = 1.35; 95% CI 1.01–1.81). The allele C of rs9279 on VDR, was negatively associated with asthma risk (OR = 0.66; 95% CI 0.45–0.97), vitamin D insufficiency (OR = 0.78; 95% CI 0.70–0.96) and higher VDR expression. Two variants in VDR were associated with asthma severity, the allele A of rs2189480 (OR = 0.34; 95% CI 0.13–0.89) and the allele G of rs4328262 (OR = 3.18; 95% CI 1.09–9.28). The combination of variants in CYP2R1 and CYP24A1 (GAC, to rs10500804, rs12794714 and rs3886163, respectively) was negatively associated with vitamin D production (β = − 1.24; 95% CI − 2.42 to − 0.06).

Conclusions

Genetic variants in the vitamin D pathway affect vitamin D serum levels and, thus, atopy and asthma.

Spectrum of gene mutations identified by targeted next-generation sequencing in Chinese leukemia patients

BACKGROUND

Despite targeted sequencing have identified several mutations for leukemia, there is still a limit of mutation screening for Chinese leukemia. Here, we used targeted next-generation sequencing for testing the mutation patterns of Chinese leukemia patients.

METHODS

We performed targeted sequencing of 504 tumor-related genes in 109 leukemia samples to identify single-nucleotide variants (SNVs) and insertions and deletions (INDELs). Pathogenic variants were assessed based on the American College of Medical Genetics and Genomics (ACMG) guidelines. The functional impact of pathogenic genes was explored through gene ontology (GO), pathway analysis, and protein-protein interaction network in silico.

RESULTS

We identified a total of 4,655 SNVs and 614 INDELs in 419 genes, in which PDE4DIP, NOTCH2, FANCA, BCR, and ROS1 emerged as the highly mutated genes. Of note, we were the first to demonstrate an association of PDE4DIP mutation and leukemia. Based on ACMG guidelines, 39 pathogenic and likely pathogenic mutations in 27 genes were found. GO annotation showed that the biological process including gland development, leukocyte differentiation, respiratory system development, myeloid leukocyte differentiation, mesenchymal to epithelial transition, and so on were involved.

CONCLUSION

Our study provided a map of gene mutations in Chinese patients with leukemia and gave insights into the molecular pathogenesis of leukemia.

Origins and characterization of variants shared between databases of somatic and germline human mutations

Background

Mutations arise in the human genome in two major settings: the germline and the soma. These settings involve different inheritance patterns, time scales, chromatin structures, and environmental exposures, all of which impact the resulting distribution of substitutions. Nonetheless, many of the same single nucleotide variants (SNVs) are shared between germline and somatic mutation databases, such as between the gnomAD database of 120,000 germline exomes and the TCGA database of 10,000 somatic exomes. Here, we sought to explain this overlap.

Results

After strict filtering to exclude common germline polymorphisms and sites with poor coverage or mappability, we found 336,987 variants shared between the somatic and germline databases. A uniform statistical model explains 34% of these shared variants; a model that incorporates the varying mutation rates of the basic mutation types explains another 50% of shared variants; and a model that includes extended nucleotide contexts (e.g. surrounding 3 bases on either side) explains an additional 4% of shared variants. Analysis of read depth finds mixed evidence that up to 4% of the shared variants may represent germline variants leaked into somatic call sets. 9% of the shared variants are not explained by any model. Sequencing errors and convergent evolution did not account for these. We surveyed other factors as well: Cancers driven by endogenous mutational processes share a greater fraction of variants with the germline, and recently derived germline variants were more likely to be somatically shared than were ancient germline ones.

Conclusions

Overall, we find that shared variants largely represent bona fide biological occurrences of the same variant in the germline and somatic setting and arise primarily because DNA has some of the same basic chemical vulnerabilities in either setting. Moreover, we find mixed evidence that somatic call-sets leak appreciable numbers of germline variants, which is relevant to genomic privacy regulations. In future studies, the similar chemical vulnerability of DNA between the somatic and germline settings might be used to help identify disease-related genes by guiding the development of background-mutation models that are informed by both somatic and germline patterns of variation.

Insights on early mutational events in SARS-CoV-2 virus reveal founder effects across geographical regions

Here we aim to describe early mutational events across samples from publicly available SARS-CoV-2 sequences from the sequence read archive and GenBank repositories. Up until 27 March 2020, we downloaded 50 illumina datasets, mostly from China, USA (WA State) and Australia (VIC). A total of 30 datasets (60%) contain at least a single founder mutation and most of the variants are missense (over 63%). Five-point mutations with clonal (founder) effect were found in USA next-generation sequencing samples. Sequencing samples from North America in GenBank (22 April 2020) present this signature with up to 39% allele frequencies among samples (n = 1,359). Australian variant signatures were more diverse than USA samples, but still, clonal events were found in these samples. Mutations in the helicase, encoded by the ORF1ab gene in SARS-CoV-2 were predominant, among others, suggesting that these regions are actively evolving. Finally, we firmly urge that primer sets for diagnosis be carefully designed, since rapidly occurring variants would affect the performance of the reverse transcribed quantitative PCR (RT-qPCR) based viral testing.

Rare variants in dynein heavy chain genes in two individuals with situs inversus and developmental dyslexia: a case report

Background

Developmental dyslexia (DD) is a neurodevelopmental learning disorder with high heritability. A number of candidate susceptibility genes have been identified, some of which are linked to the function of the cilium, an organelle regulating left-right asymmetry development in the embryo. Furthermore, it has been suggested that disrupted left-right asymmetry of the brain may play a role in neurodevelopmental disorders such as DD. However, it is unknown whether there is a common genetic cause to DD and laterality defects or ciliopathies.

Case presentation

Here, we studied two individuals with co-occurring situs inversus (SI) and DD using whole genome sequencing to identify genetic variants of importance for DD and SI. Individual 1 had primary ciliary dyskinesia (PCD), a rare, autosomal recessive disorder with oto-sino-pulmonary phenotype and SI. We identified two rare nonsynonymous variants in the dynein axonemal heavy chain 5 gene (DNAH5): a previously reported variant c.7502G > C; p.(R2501P), and a novel variant c.12043 T > G; p.(Y4015D). Both variants are predicted to be damaging. Ultrastructural analysis of the cilia revealed a lack of outer dynein arms and normal inner dynein arms. MRI of the brain revealed no significant abnormalities. Individual 2 had non-syndromic SI and DD. In individual 2, one rare variant (c.9110A > G;p.(H3037R)) in the dynein axonemal heavy chain 11 gene (DNAH11), coding for another component of the outer dynein arm, was identified.

Conclusions

We identified the likely genetic cause of SI and PCD in one individual, and a possibly significant heterozygosity in the other, both involving dynein genes. Given the present evidence, it is unclear if the identified variants also predispose to DD and further studies into the association between laterality, ciliopathies and DD are needed.

Meta-analysis of massively parallel reporter assays enables prediction of regulatory function across cell types

Deciphering the potential of noncoding loci to influence gene regulation has been the subject of intense research, with important implications in understanding genetic underpinnings of human diseases. Massively parallel reporter assays (MPRAs) can measure regulatory activity of thousands of DNA sequences and their variants in a single experiment. With increasing number of publically available MPRA data sets, one can now develop data-driven models which, given a DNA sequence, predict its regulatory activity. Here, we performed a comprehensive meta-analysis of several MPRA data sets in a variety of cellular contexts. We first applied an ensemble of methods to predict MPRA output in each context and observed that the most predictive features are consistent across data sets. We then demonstrate that predictive models trained in one cellular context can be used to predict MPRA output in another, with loss of accuracy attributed to cell-type-specific features. Finally, we show that our approach achieves top performance in the Fifth Critical Assessment of Genome Interpretation "Regulation Saturation" Challenge for predicting effects of single-nucleotide variants. Overall, our analysis provides insights into how MPRA data can be leveraged to highlight functional regulatory regions throughout the genome and can guide effective design of future experiments by better prioritizing regions of interest.

Rare single-nucleotide variants in oculo-auriculo-vertebral spectrum (OAVS)

Background

Oculo‐auriculo‐vertebral spectrum (OAVS) is a craniofacial developmental disorder that affects structures derived from the first and second pharyngeal arches. The clinically heterogeneous phenotype involves mandibular, oral, and ear development anomalies. Etiology is complex and poorly understood. Genetic factors have been associated, evidenced by chromosomal abnormalities affecting different genomic regions and genes. However, known pathogenic single‐nucleotide variants (SNVs) have only been identified in MYT1 in a restricted number of patients. Therefore, investigations of SNVs on candidate genes may reveal other pathogenic mechanisms.

Methods

In a cohort of 73 patients, coding and untranslated regions (UTR) of 10 candidate genes (CRKL, YPEL1, MAPK1, NKX3‐2, HMX1, MYT1, OTX2, GSC, PUF60, HOXA2) were sequenced. Rare SNVs were selected and in silico predictions were performed to ascertain pathogenicity. Likely pathogenic variants were validated by Sanger sequencing and heritability was assessed when possible.

Results

Four likely pathogenic variants in heterozygous state were identified in different patients. Two SNVs were located in the 5’UTR of YPEL1; one in the 3’UTR of CRKL and one in the 3’UTR of OTX2.

Conclusion

Our work described variants in candidate genes for OAVS and supported the genetic heterogeneity of the spectrum.

Deciphering the structural basis for glucocorticoid resistance caused by missense mutations in the ligand binding domain of glucocorticoid receptor

The glucocorticoid resistance hereditary condition may emerge from the occurrence of point mutations in the glucocorticoid receptor (GR), which could impair its functionality. Because the main feature of such pathology is the resistance of the hypothalamic-pituitary-adrenal axis to the hormone cortisol, we used the GR ligand binding domain three-dimensional structure to perform computational analysis for eight variants known to cause this clinical condition (I559 N, V571A, D641V, G679S, F737L, I747 M, L753F and L773P), aiming to understand, on the atom scale, how they cause glucocorticoid resistance. We observed that the mutations generated a reduced affinity to cortisol and they alter some loop conformations, which could be a consequence from changes in protein motion, which in turn could result from the reduced stability of mutant GR structures. Therefore, the analyzed mutations compromise the GR ligand binding domain structure and cortisol binding, which could characterize the glucocorticoid resistance phenotype.

Beyond sequencing: re-visiting annotations for PJL as a test case

Objectives

Current developments in sequencing techniques have enabled rapid and high-throughput generation of sequence data. However, there is a growing gap between the generation of raw sequencing data and the extraction of meaningful biological information. Variant annotation is a crucial step in the analysis of genome sequencing data. Incorrect or incomplete annotations can cause researchers to dilute interesting variants in a pool of false positives. We require consistent, accurate and reliable annotation of variants for making diagnostic and treatment decisions. Current annotation depends on the set of transcripts, and software used can be managed, with sufficient care, in the research context. Careful thought needs to be given to the choice of transcript sets and software packages for variant annotation in sequencing studies. In this project, the main objective is to analyze the genetic variants observed in Pakistani population data within the 1000 genomes project (1KGP).

Results

We characterized only SNVs and InDels types of genetic variations, in total ~ 1.4 million variants. Besides this, we also annotated the genetic variants with multiple annotations tools, ANNOVAR and SnpEff and compared the differential results. Our population-specific catalogue will enhance future studies on the functional impact at protein level.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4508-5) contains supplementary material, which is available to authorized users.

Mitochondrial DNA D-loop sequencing reveals obesity variants in an Arab population

Background: The association of mitochondrial DNA (mtDNA) variations with obesity has been investigated in diverse populations across the world. However, such obesity-associated mtDNA examinations are rarely conducted in Arab populations. Materials and methods: We re-sequenced mtDNA displacement loop (D-loop) region of 395 Arab individuals of Kuwait. We categorized the individuals based on their BMI scores as obese (n=232; BMI ≥30 kg/m²), overweight (n=110; BMI ≥25 kg/m² and <30 kg/m²), and lean (n=53; BMI <25 kg/m²). We performed all the statistical tests by combining obese and overweight individuals in one group. Association analyses were conducted applying Fisher's exact test and logistic regression model. Results: We identified that the mtDNA variations m.73A>G, and m.523delAC were positively correlated with obesity, while m.310T>C, and m.16318A>T were negatively associated. All these variants, except m.16318A>T, remain statistically significant after adjusting for age and gender. We found that the variant m.73A>G increases the likelihood of being obese by 6-fold, whereas haplogroup H decreases the probability of being obese in Arab individuals of Kuwait. Haplotype analysis revealed that a haplotype, A263G-C309CT-T310C, defining the H2a clade of H haplogroup, reduces the probability of being obese. Conclusion: Our study reports, for the first time, the obesity-related mtDNA variants in Arabs of Kuwait. Based on the mtDNA D-loop region variations, we detected particular variants and haplogroup that are related with increased and decreased probability of being obese in the Kuwait Arab population.

Comparative genome-wide survey of single nucleotide variation uncovers the genetic diversity and potential biomedical applications among six Macaca species

Macaca is of great importance in evolutionary and biomedical research. Aiming at elucidating genetic diversity patterns and potential biomedical applications of macaques, we characterized single nucleotide variations (SNVs) of six Macaca species based on the reference genome of Macaca mulatta. Using eight whole-genome sequences, representing the most comprehensive genomic SNV study in Macaca to date, we focused on discovery and comparison of nonsynonymous SNVs (nsSNVs) with bioinformatic tools. We observed that SNV distribution patterns were generally congruent among the eight individuals. Outlier tests of nsSNV distribution patterns detected 319 bins with significantly distinct genetic divergence among macaques, including differences in genes associated with taste transduction, homologous recombination, and fat and protein digestion. Genes with specific nsSNVs in various macaques were differentially enriched for metabolism pathways, such as glycolysis, protein digestion and absorption. On average, 24.95% and 11.67% specific nsSNVs were putatively deleterious according to PolyPhen2 and SIFT4G, respectively, among which the shared deleterious SNVs were located in 564–1981 genes. These genes displayed enrichment signals in the ‘obesity-related traits’ disease category for all surveyed macaques, confirming that they were suitable models for obesity related studies. Additional enriched disease categories were observed in some macaques, exhibiting promising potential for biomedical application. Positively selected genes identified by PAML in most tested Macaca species played roles in immune and nervous system, growth and development, and fat metabolism. We propose that metabolism and body size play important roles in the evolutionary adaptation of macaques.

Lightning-fast genome variant detection with GROM

Current human whole genome sequencing projects produce massive amounts of data, often creating significant computational challenges. Different approaches have been developed for each type of genome variant and method of its detection, necessitating users to run multiple algorithms to find variants. We present Genome Rearrangement OmniMapper (GROM), a novel comprehensive variant detection algorithm accepting aligned read files as input and finding SNVs, indels, structural variants (SVs), and copy number variants (CNVs). We show that GROM outperforms state-of-the-art methods on 7 validated benchmarks using 2 whole genome sequencing (WGS) data sets. Additionally, GROM boasts lightning-fast run times, analyzing a 50× WGS human data set (NA12878) on commonly available computer hardware in 11 minutes, more than an order of magnitude (up to 72 times) faster than tools detecting a similar range of variants. Addressing the needs of big data analysis, GROM combines in 1 algorithm SNV, indel, SV, and CNV detection, providing superior speed, sensitivity, and precision. GROM is also able to detect CNVs, SNVs, and indels in non-paired-read WGS libraries, as well as SNVs and indels in whole exome or RNA sequencing data sets.

Application of next-generation sequencing to characterize novel mutations in clarithromycin-susceptible Helicobacter pylori strains with A2143G of 23S rRNA gene

Background

Clarithromycin (CLR) resistance has become a predominant factor for treatment failure of Helicobacter pylori eradication. Although the molecular mechanism of CLR resistance has been clearly understood in H. pylori, it is lack of evidence of other genes involved in drug resistance. Furthermore, the molecular mechanism of phenotype susceptible to CLR while genotype of 23S rRNA is mutant with A2143G is unclear. Here, we characterized the mutations of CLR-resistant and -susceptible H. pylori strains to explore bacterial resistance.

Methods

In the present study, the whole genomes of twelve clinical isolated H. pylori strains were sequenced, including two CLR-susceptible strains with mutation of A2143G. Single nucleotide variants (SNVs) were extracted and analyzed from multidrug efflux transporter genes.

Results

We did not find mutations associated with known CLR-resistant sites except for controversial T2182C outside of A2143G in the 23S rRNA gene. Although total SNVs of multidrug efflux transporter gene and the SNVs of HP0605 were significant differences (P ≤ 0.05) between phenotype resistant and susceptible strains. There is no significant difference in SNVs of RND or MFS (HP1181) family. However, the number of mutations in the RND family was significantly higher in the mutant strain (A2143G) than in the wild type. In addition, three special variations from two membrane proteins of mtrC and hefD were identified in both CLR-susceptible strains with A2143G.

Conclusions

Next-generation sequencing is a practical strategy for analyzing genomic variation associated with antibiotic resistance in H. pylori. The variations of membrane proteins of the RND family may be able to participate in the regulation of clinical isolated H. pylori susceptibility profiles.

Tracing the origin of Treponema pallidum in China using next-generation sequencing

Syphilis is a systemic sexually transmitted disease caused by Treponema pallidum ssp. pallidum (TPA). The origin and genetic background of Chinese TPA strains remain unclear. We identified a total of 329 single-nucleotide variants (SNVs) in eight Chinese TPA strains using next-generation sequencing. All of the TPA strains were clustered into three lineages, and Chinese TPA strains were grouped in Lineage 2 based on phylogenetic analysis. The phylogeographical data showed that TPA strains originated earlier than did T. pallidum ssp. pertenue (TPE) and T. pallidum ssp. endemicum (TPN) strains and that Chinese TPA strains might be derived from recombination between Lineage 1 and Lineage 3. Moreover, we found through a homology modeling analysis that a nonsynonymous substitution (I415F) in the PBP3 protein might affect the structural flexibility of PBP3 and the binding constant for substrates based on its possible association with penicillin resistance in T. pallidum. Our findings provide new insight into the molecular foundation of the evolutionary origin of TPA and support the development of novel diagnostic/therapeutic technology for syphilis.

Using RNA-Seq Data for the Detection of a Panel of Clinically Relevant Mutations

Somatic single nucleotide variants (SNVs) are genomic events with increasing implications in cancer treatment. The clinical standard for SNVs detection is whole genome/exome sequencing (WGS/WES) in matched tumor-normal samples. Yet, this is a very costly approach both economically and biologically and very often only tumor samples are sequenced. On the other hand, RNA sequencing (RNA-Seq) is the most popular technology to study gene expression, and has also the potential for a cost-effective identification of SNVs as an alternative to tumor-only WES. Here we present a method for the identification of SNVs in tumor-only RNA-Seq data putting a special focus on a small panel of clinically relevant SNVs. For evaluation purposeswe analyzed matched tumor-normal WEStumor-only RNA-Seq data from 14 cancer patients. We compared SNVs detected in i) RNA-Seq by our method, ii) WES tumor-only by Mutect2 and iii) WES matched tumor-normal by Mutect2. We did a detailed evaluation for a reduced panel of clinically relevant SNVs and reliably identified in RNA-Seq data a subset of mutations for which we had pathological annotation. Hence, RNA-Seq rises as a cost-effective option to detect in parallel gene expression as well as a small panel of clinically relevant SNVs in research.

Genetics of Schizophrenia: Overview of Methods, Findings and Limitations

Genetics constitute a crucial risk factor to schizophrenia. In the last decade, molecular genetic research has produced novel findings, infusing optimism about discovering the biological roots of schizophrenia. However, the complexity of the object of inquiry makes it almost impossible for non-specialists in genetics (e.g., many clinicians and researchers) to get a proper understanding and appreciation of the genetic findings and their limitations. This study aims at facilitating such an understanding by providing a brief overview of some of the central methods and findings in schizophrenia genetics, from its historical origins to its current status, and also by addressing some limitations and challenges that confront this field of research. In short, the genetic architecture of schizophrenia has proven to be highly complex, heterogeneous and polygenic. The disease risk is constituted by numerous common genetic variants of only very small individual effect and by rare, highly penetrant genetic variants of larger effects. In spite of recent advances in molecular genetics, our knowledge of the etiopathogenesis of schizophrenia and the genotype-environment interactions remain limited.

Gene structure of the pregnancy-associated glycoprotein-like (PAG-L) in the Eurasian beaver (Castor fiber L.)

The pregnancy-associated glycoprotein-like family (PAG-L) is a large group of chorionic products, expressed in the pre-placental trophoblast and later in the post-implantational chorionic epithelium, and are involved in proper placenta development and embryo-maternal interaction in eutherians. This study describes identification of the PAG-L family in the genome of the Eurasian beaver (Castor fiber L.), named CfPAG-L. We identified 7657 bp of the CfPAG-L gDNA sequence (Acc. No. KX377932), encompassing nine exons (1-9) and eight introns (A-H). The length of the CfPAG-L exons (59-200 bp) was equivalently similar to the only known counterparts of bPAG1, bPAG2, and pPAG2. The length of the CfPAG-L introns ranged 288-1937 bp and was completely different from previously known PAG introns. The exonic CfPAG-L regions revealed 50.3-72.9% homology with equivalent segments of bPAG1 and pPAG2 structure. The intronic CfPAG-L regions alignments revealed a lack of homology. Within the entire CfPAG-L gene, 31 potential single nucleotide variants (SNV: 7 transversions and 24 transitions) were predicted. The identified exonic polymorphic loci did not affect the amino acid sequence of the CfPAG-L polypeptide precursor. This is the first report describing the CfPAG-L gene sequence, structural organization, and SNVs in the Eurasian beaver, one of the largest rodents.