A Deep Neural Network for Predicting and Engineering Alternative Polyadenylation.

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3'UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3' end processing, and integrates these features into a comprehensive, interpretable, cis-regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease.

A common single nucleotide variant in the cytokine receptor-like factor-3 (CRLF3) gene causes neuronal deficits in human and mouse cells.

Single nucleotide variants in the general population are common genomic alterations, where the majority are presumed to be silent polymorphisms without known clinical significance. Using human induced pluripotent stem cell (hiPSC) cerebral organoid modeling of the 1.4 megabase Neurofibromatosis type 1 (NF1) deletion syndrome, we previously discovered that the cytokine receptor-like factor-3 (CRLF3) gene, which is co-deleted with the NF1 gene, functions as a major regulator of neuronal maturation. Moreover, children with NF1 and the CRLF3L389P variant have greater autism burden, suggesting that this gene might be important for neurologic function. To explore the functional consequences of this variant, we generated CRLF3L389P-mutant hiPSC lines and Crlf3L389P-mutant genetically engineered mice. While this variant does not impair protein expression, brain structure, or mouse behavior, CRLF3L389P-mutant human cerebral organoids and mouse brains exhibit impaired neuronal maturation and dendrite formation. In addition, Crlf3L389P-mutant mouse neurons have reduced dendrite lengths and branching, without any axonal deficits. Moreover, Crlf3L389P-mutant mouse hippocampal neurons have decreased firing rates and synaptic current amplitudes relative to wild type controls. Taken together, these findings establish the CRLF3L389P variant as functionally deleterious and suggest that it may be a neurodevelopmental disease modifier.

Improving somatic exome sequencing performance by biological replicates.

BACKGROUND: Next-generation sequencing (NGS) technologies offer fast and inexpensive identification of DNA sequences. Somatic sequencing is among the primary applications of NGS, where acquired (non-inherited) variants are based on comparing diseased and healthy tissues from the same individual. Somatic mutations in genetic diseases such as cancer are tightly associated with genomic instability. Genomic instability increases heterogenity, complicating sequencing efforts further, a task already challenged by the presence of short reads and repetitions in human DNA. This leads to low concordance among studies and limits reproducibility. This limitation is a significant problem since identified mutations in somatic sequencing are major biomarkers for diagnosis and the primary input of targeted therapies. Benchmarking studies were conducted to assess the error rates and increase reproducibility. Unfortunately, the number of somatic benchmarking sets is very limited due to difficulties in validating true somatic variants. Moreover, most NGS benchmarking studies are based on relatively simpler germline (inherited) sequencing. Recently, a comprehensive somatic sequencing benchmarking set was published by Sequencing Quality Control Phase 2 (SEQC2). We chose this dataset for our experiments because it is a well-validated, cancer-focused dataset that includes many tumor/normal biological replicates. Our study has two primary goals. First goal is to determine how replicate-based consensus approaches can improve the accuracy of somatic variant detection systems. Second goal is to develop highly predictive machine learning (ML) models by employing replicate-based consensus variants as labels during the training phase. RESULTS: Ensemble approaches that combine alternative algorithms are relatively common; here, as an alternative, we study the performance enhancement potential of biological replicates. We first developed replicate-based consensus approaches that utilize the biological replicates available in this study to improve variant calling performance. Subsequently, we trained ML models using these biological replicates and achieved performance comparable to optimal ML models, those trained using high-confidence variants identified in advance. CONCLUSIONS: Our replicate-based consensus approach can be used to improve variant calling performance and develop efficient ML models. Given the relative ease of obtaining biological replicates, this strategy allows for the development of efficient ML models tailored to specific datasets or scenarios.

Identification of novel single nucleotide variants in the drug resistance mechanism of Mycobacterium tuberculosis isolates by whole-genome analysis.

BACKGROUND: Tuberculosis (TB) represents a major global health challenge. Drug resistance in Mycobacterium tuberculosis (MTB) poses a substantial obstacle to effective TB treatment. Identifying genomic mutations in MTB isolates holds promise for unraveling the underlying mechanisms of drug resistance in this bacterium. METHODS: In this study, we investigated the roles of single nucleotide variants (SNVs) in MTB isolates resistant to four antibiotics (moxifloxacin, ofloxacin, amikacin, and capreomycin) through whole-genome analysis. We identified the drug-resistance-associated SNVs by comparing the genomes of MTB isolates with reference genomes using the MuMmer4 tool. RESULTS: We observed a strikingly high proportion (94.2%) of MTB isolates resistant to ofloxacin, underscoring the current prevalence of drug resistance in MTB. An average of 3529 SNVs were detected in a single ofloxacin-resistant isolate, indicating a mutation rate of approximately 0.08% under the selective pressure of ofloxacin exposure. We identified a set of 60 SNVs associated with extensively drug-resistant tuberculosis (XDR-TB), among which 42 SNVs were non-synonymous mutations located in the coding regions of nine key genes (ctpI, desA3, mce1R, moeB1, ndhA, PE_PGRS4, PPE18, rpsA, secF). Protein structure modeling revealed that SNVs of three genes (PE_PGRS4, desA3, secF) are close to the critical catalytic active sites in the three-dimensional structure of the coding proteins. CONCLUSION: This comprehensive study elucidates novel resistance mechanisms in MTB against antibiotics, paving the way for future design and development of anti-tuberculosis drugs.

MiRiQ Database: A Platform for In Silico Rice Mutant Screening.

Genetic studies using mutant resources have significantly contributed to elucidating plant gene function. Massive mutant libraries sequenced by next-generation sequencing technology facilitate mutant identification and functional analysis of genes of interest. Here, we report the creation and release of an open-access database (https://miriq.agr.kyushu-u.ac.jp/index.php), called Mutation-induced Rice in Kyushu University (MiRiQ), designed for in silico mutant screening based on a whole-genome-sequenced mutant library. This database allows any user to easily find mutants of interest without laborious efforts such as large-scale screening by PCR. The initial version of the MiRiQ database (version 1.0) harbors a total of 1.6 million single-nucleotide variants (SNVs) and InDels of 721 M1 plants that were mutagenized by N-methyl-N-nitrosourea treatment of the rice cultivar Nipponbare (Oryza sativa ssp. japonica). The SNVs were distributed among 87% of all 35,630 annotated protein-coding genes of the Nipponbare genome and were predicted to induce missense and nonsense mutations. The MiRiQ database provides built-in tools, such as a search tool by keywords and JBrowse for mutation searches. Users can request mutant seeds in the M2 or M3 generations from a request form linked to this database. We believe that the availability of a wide range of gene mutations in this database will benefit the plant science community and breeders worldwide by accelerating functional genomic research and crop improvement.

New approaches to predict the effect of co-occurring variants on protein characteristics.

Predicting the effect of a mutated gene before the onset of symptoms of genetic diseases would greatly facilitate diagnosis and potentiate early intervention. There have been myriad attempts to predict the effects of single-nucleotide variants. However, the applicability of these efforts does not scale to co-occurring variants. Furthermore, an increasing number of protein therapeutics contain co-occurring nucleotide variations, adding uncertainty during development to the safety and efficiency of these drugs. Co-occurring nucleotide variants may often have synergistic, additive, or antagonistic effects on protein attributes, further complicating the task of outcome prediction. We tested four models based on the cooperative and antagonistic effects of co-occurring variants to predict pathogenicity and effectiveness of protein therapeutics. A total of 30 attributes, including amino acid and nucleotide features, as well as existing single-variant effect prediction tools, were considered on the basis of previous studies on single-nucleotide variants. Importantly, the effects of synonymous variants, often seen in protein therapeutics, were also included in our models. We used 12 datasets of people with monogenic diseases and controls with co-occurring genetic variants to evaluate the accuracy of our models, accomplishing a degree of accuracy comparable to that of prediction tools for single-nucleotide variants. More importantly, our framework is generalizable to new, well-curated datasets of monogenic diseases and new variant scoring tools. This approach successfully assists in addressing the challenging task of predicting the effect of co-occurring variants on pathogenicity and protein effectiveness and is applicable for a wide range of protein therapeutics and genetic diseases.

Comprehensive genetic analysis of facioscapulohumeral muscular dystrophy by Nanopore long-read whole-genome sequencing.

BACKGROUND: Facioscapulohumeral muscular dystrophy (FSHD) is a high-prevalence autosomal dominant neuromuscular disease characterized by significant clinical and genetic heterogeneity. Genetic diagnosis of FSHD remains a challenge because it cannot be detected by standard sequencing methods and requires a complex diagnosis workflow. METHODS: We developed a comprehensive genetic FSHD detection method based on Oxford Nanopore Technologies (ONT) whole-genome sequencing. Using a case-control design, we applied this procedure to 29 samples and compared the results with those from optical genome mapping (OGM), bisulfite sequencing (BSS), and whole-exome sequencing (WES). RESULTS: Using our ONT-based method, we identified 59 haplotypes (35 4qA and 24 4qB) among the 29 samples (including a mosaic sample), as well as the number of D4Z4 repeat units (RUs). The pathogenetic D4Z4 RU contraction identified by our ONT-based method showed 100% concordance with OGM results. The methylation levels of the most distal D4Z4 RU and the double homeobox 4 gene (DUX4) detected by ONT sequencing are highly consistent with the BSS results and showed excellent diagnostic efficiency. Additionally, our ONT-based method provided an independent methylation profile analysis of two permissive 4qA alleles, reflecting a more accurate scenario than traditional BSS. The ONT-based method detected 17 variations in three FSHD2-related genes from nine samples, showing 100% concordance with WES. CONCLUSIONS: Our ONT-based FSHD detection method is a comprehensive method for identifying pathogenetic D4Z4 RU contractions, methylation level alterations, allele-specific methylation of two 4qA haplotypes, and variations in FSHD2-related genes, which will all greatly improve genetic testing for FSHD.

Association of rare single nucleotide variant MUC5B rs35705950 with interstitial lung disease in Japanese rheumatoid arthritis.

OBJECTIVES: Rheumatoid arthritis (RA) is sometimes complicated by interstitial lung disease (ILD) with a poor prognosis. A single nucleotide variant (SNV) in MUC5B was associated with ILD in European RA patients. However, associations of this SNV were not found in Japanese RA patients, because its frequency in Japanese populations is very low. We investigated the associations of candidate SNVs including the MUC5B variant with ILD in Japanese RA. METHODS: Genotyping of MUC5B rs35705950, MUC2 rs7934606, MAD1L1 rs12699415, and PPFIBP2 rs6578890 in Japanese RA patients was conducted for association analyses. RESULTS: MUC5B rs35705950 was associated with usual interstitial pneumonia (UIP) (p = 0.0039, Pc = 0.0156, odds ratio [OR] 10.66, 95% confidence interval [CI] 2.05-55.37) or ILD (p = 0.0071, Pc = 0.0284, OR 7.33, 95%CI 1.52-35.44) in Japanese RA under the allele model. MUC2 rs7934606 was associated with UIP (p = 0.0072, Pc = 0.0288, OR 29.55, 95%CI 1.52-574.57) or ILD (p = 0.0037, Pc = 0.0148, OR 22.95, 95%CI 1.27-416.13) in RA. Haplotype analyses suggested the primary association of MUC5B rs35705950 with UIP in Japanese RA. No significant association of MAD1L1 rs12699415 or PPFIBP2 rs6578890 with UIP, nonspecific interstitial pneumonia, or ILD in RA was observed. CONCLUSIONS: MUC5B rs35705950 is associated with, and might be involved in the pathogenesis of ILD, especially UIP, in Japanese RA.

TNFRSF1B Gene Variants in Clinicopathological Aspects and Prognosis of Patients with Cutaneous Melanoma.

Regulatory T lymphocytes play a critical role in immune regulation and are involved in the aberrant cell elimination by facilitating tumor necrosis factor connection to the TNFR2 receptor, encoded by the TNFRSF1B polymorphic gene. We aimed to examine the effects of single nucleotide variants TNFRSF1B c.587T>G, c.*188A>G, c.*215C>T, and c.*922C>T on the clinicopathological characteristics and survival of cutaneous melanoma (CM) patients. Patients were genotyped using RT-PCR. TNFRSF1B levels were measured using qPCR. Luciferase reporter assay evaluated the interaction of miR-96 and miR-1271 with the 3'-UTR of TNFRSF1B. The c.587TT genotype was more common in patients younger than 54 years old than in older patients. Patients with c.*922CT or TT, c.587TG or GG + c.*922CT or TT genotypes, as well as those with the haplotype TATT, presented a higher risk of tumor progression and death due to the disease effects. Individuals with the c.*922TT genotype had a higher TNFRSF1B expression than those with the CC genotype. miR-1271 had less efficient binding with the 3'-UTR of the T allele when compared with the C allele of the SNV c.*922C>T. Our findings, for the first time, demonstrate that TNFRSF1B c.587T>G and c.*922C>T variants can serve as independent prognostic factors in CM patients.

Single-Base Gene Variants in MIR-146A and SCN1A Genes Related to the Epileptogenic Process in Drug-Responsive and Drug-Resistant Temporal Lobe Epilepsy-A Preliminary Study in a Brazilian Cohort Sample.

The drug-resistant temporal lobe epilepsy (TLE) has recently been associated with single nucleotide variants (SNVs) in microRNA(miR)-146a (MIR-146A) (rs2910164) and Sodium Voltage-Gated Channel Alpha Subunit 1 (SCN1A) (rs2298771 and rs3812718) genes. Moreover, no studies have shown an association between these SNVs and susceptibility to drug-resistant and drug-responsive TLE in Brazil. Thus, deoxyribonucleic acid (DNA) samples from 120 patients with TLE (55 drug-responsive and 65 drug-resistant) were evaluated by real-time polymerase chain reaction (RT-PCR). A total of 1171 healthy blood donor individuals from the Online Archive of Brazilian Mutations (ABraOM, from Portuguese Arquivo Brasileiro On-line de Mutações), a repository containing genomic variants of the Brazilian population, were added as a control population for the studied SNVs. MIR-146A and SCN1A relative expression was performed by quantitative RT-PCR (qRT-PCR). The statistical analysis protocol was performed using an alpha error of 0.05. TLE patient samples and ABraOM control samples were in Hardy-Weinberg equilibrium for all studied SNVs. For rs2910164, the frequencies of the homozygous genotype (CC) (15.00% vs. 9.65%) and C allele (37.80% vs. 29.97%) were superior in patients with TLE compared to controls with a higher risk for TLE disease [odds ratio (OR) = 1.89 (95% confidence interval (95%CI) = 1.06-3.37); OR = 1.38 (95%CI = 1.04-1.82), respectively]. Drug-responsive patients also presented higher frequencies of the CC genotype [21.81% vs. 9.65%; OR = 2.58 (95%CI = 1.25-5.30)] and C allele [39.09% vs. 29.97%; OR = 1.50 (95%CI = 1.01-2.22)] compared to controls. For rs2298771, the frequency of the heterozygous genotype (AG) (51.67% vs. 40.40%) was superior in patients with TLE compared to controls with a higher risk for TLE disease [OR = 2.42 (95%CI = 1.08-5.41)]. Drug-resistant patients presented a higher AG frequency [56.92% vs. 40.40%; OR = 3.36 (95%CI = 1.04-17.30)] compared to the control group. For rs3812718, the prevalence of genotypes and alleles were similar in both studied groups. The MIR-146A relative expression level was lower in drug-resistant compared to drug-responsive patients for GC (1.6 vs. 0.1, p-value = 0.049) and CC (1.8 vs. 0.6, p-value = 0.039). Also, the SCN1A relative expression levels in samples from TLE patients were significantly higher in AG [2.09 vs. 1.10, p-value = 0.038] and GG (3.19 vs. 1.10, p-value < 0.001) compared to the AA genotype. In conclusion, the rs2910164-CC and rs2298771-AG genotypes are exerting significant risk influence, respectively, on responsive disease and resistant disease, probably due to an upregulated nuclear factor kappa B (NF-kB) and SCN1A loss of function.

Association between mtDNA haplogroups and skeletal fluorosis in Han population residing in drinking water endemic fluorosis area of northern China.

To investigate the association between mtDNA genetic information and the risk of SF, individuals were conducted in the drinking water endemic fluorosis area in northern China, sequenced the whole genome of mtDNA, identified the SNPs and SNVs, analyzed the haplogroups, and diagnosed SF, and then, the effect of mtDNA genetic information on the risk of SF was evaluated. We find that, D5 haplogroup and its specific SNPs reduced the risk, while the D4 haplogroup and its specific SNPs increased the risk of SF. The number of SNVs in coding regions of mitochondrial respiratory chain (MRC) is different between the controls and cases. This suggests that D5 haplogroup may play a protective role in the risk of SF, while the opposite is observed for the D4 haplogroup, this may relate to their specific SNPs. And SNVs that encode the MRC complex may also be associated with the risk of SF.

Identification of Novel Genetic Risk Factors for Focal Segmental Glomerulosclerosis in Children: Results From the Chronic Kidney Disease in Children (CKiD) Cohort.

RATIONALE & OBJECTIVE: Focal segmental glomerulosclerosis (FSGS) is a major cause of pediatric nephrotic syndrome, and African Americans exhibit an increased risk for developing FSGS compared with other populations. Predisposing genetic factors have previously been described in adults. Here we performed genomic screening of primary FSGS in a pediatric African American population. STUDY DESIGN: Prospective cohort with case-control genetic association study design. SETTING & PARTICIPANTS: 140 African American children with chronic kidney disease from the Chronic Kidney Disease in Children (CKiD) cohort, including 32 cases with FSGS. PREDICTORS: Over 680,000 common single-nucleotide polymorphisms (SNPs) were tested for association. We also ran a pathway enrichment analysis and a human leucocyte antigen (HLA)-focused association study. OUTCOME: Primary biopsy-proven pediatric FSGS. ANALYTICAL APPROACH: Multivariate logistic regression models. RESULTS: The genome-wide association study revealed 169 SNPs from 14 independent loci significantly associated with FSGS (false discovery rate [FDR]<5%). We observed notable signals for genetic variants within the APOL1 (P=8.6×10-7; OR, 25.8 [95% CI, 7.1-94.0]), ALMS1 (P=1.3×10-7; 13.0% in FSGS cases vs 0% in controls), and FGFR4 (P=4.3×10-6; OR, 24.8 [95% CI, 6.3-97.7]) genes, all of which had previously been associated with adult FSGS, kidney function, or chronic kidney disease. We also highlighted novel, functionally relevant genes, including GRB2 (which encodes a slit diaphragm protein promoting podocyte structure through actin polymerization) and ITGB1 (which is linked to renal injuries). Our results suggest a major role for immune responses and antigen presentation in pediatric FSGS through (1) associations with SNPs in PTPRJ (or CD148, P=3.5×10-7), which plays a role in T-cell receptor signaling, (2) HLA-DRB1∗11:01 association (P=6.1×10-3; OR, 4.5 [95% CI, 1.5-13.0]), and (3) signaling pathway enrichment (P=1.3×10-6). LIMITATIONS: Sample size and no independent replication cohort with genomic data readily available. CONCLUSIONS: Our genetic study has identified functionally relevant risk factors and the importance of immune regulation for pediatric primary FSGS, which contributes to a better description of its molecular pathophysiological mechanisms. PLAIN-LANGUAGE SUMMARY: We assessed the genetic risk factors for primary focal segmental glomerulosclerosis (FSGS) by simultaneously testing over 680,000 genetic markers spread across the genome in 140 children, including 32 with FSGS lesions. Fourteen independent genetic regions were significantly associated with pediatric FSGS, including APOL1 and ALMS1-NAT8, which were previously found to be associated with FSGS and chronic kidney diseases in adults. Novel genes with relevant biological functions were also highlighted, such as GRB2 and FGFR4, which play a role in the kidney filtration barrier and in kidney cell differentiation, respectively. Finally, we revealed the importance of immune regulation in pediatric FSGS through associations involving cell surface proteins presenting antigens to the immune system and interacting with T-cell receptors.

Whole exome sequencing in unexplained recurrent miscarriage families identified novel pathogenic genetic causes of euploid miscarriage.

STUDY QUESTION: Can whole exome sequencing (WES) followed by trio bioinformatics analysis identify novel pathogenic genetic causes of first trimester euploid miscarriage? SUMMARY ANSWER: We identified genetic variants in six candidate genes that indicated plausible underlying causes of first-trimester euploid miscarriage. WHAT IS KNOWN ALREADY: Previous studies have identified several monogenic causes of Mendelian inheritance in euploid miscarriages. However, most of these studies are without trio analyses and lack cellular and animal models to validate the functional effect of putative pathogenic variants. STUDY DESIGN, SIZE, DURATION: Eight unexplained recurrent miscarriage (URM) couples and corresponding euploid miscarriages were included in our study for whole genome sequencing (WGS) and WES followed by trio bioinformatics analysis. Knock-in mice with Rry2 and Plxnb2 variants and immortalized human trophoblasts were utilized for functional study. Additional 113 unexplained miscarriages were included to identify the mutation prevalence of specific genes by multiplex PCR. PARTICIPANTS/MATERIALS, SETTING, METHODS: Whole blood from URM couples and their <13 weeks gestation miscarriage products were both collected for WES, and all variants in selected genes were verified by Sanger sequencing. Different stage C57BL/6J wild-type mouse embryos were collected for immunofluorescence. Ryr2N1552S/+, Ryr2R137W/+, Plxnb2D1577E/+, and Plxnb2R465Q/+ point mutation mice were generated and backcrossed. Matrigel-coated transwell invasion assays and wound-healing assays were performed using HTR-8/SVneo cells transfected with PLXNB2 small-interfering RNA and negative control. Multiplex PCR was performed focusing on RYR2 and PLXNB2. MAIN RESULTS AND THE ROLE OF CHANCE: Six novel candidate genes, including ATP2A2, NAP1L1, RYR2, NRK, PLXNB2, and SSPO, were identified. Immunofluorescence staining showed that ATP2A2, NAP1L1, RyR2, and PLXNB2 were widely expressed from the zygote to the blastocyst stage in mouse embryos. Although compound heterozygous mice with Rry2 and Plxnb2 variants did not show embryonic lethality, the number of pups per litter was significantly reduced when backcrossing Ryr2N1552S/+ ♂ with Ryr2R137W/+ ♀ or Plxnb2D1577E/+ ♂ with Plxnb2R465Q/+ ♀ (P < 0.05), which were in accordance with the sequencing results of Family 2 and Family 3, and the proportion of Ryr2N1552S/+ offspring was significantly lower when Ryr2N1552S/+ female mice were backcrossed with Ryr2R137W/+ male mice (P < 0.05). Moreover, siRNA-mediated PLXNB2 knockdown inhibited the migratory and invasive abilities of immortalized human trophoblasts. Besides, additional 10 variants of RYR2 and PLXNB2 were detected in 113 unexplained euploid miscarriages by multiplex PCR. LIMITATIONS, REASONS FOR CAUTION: The relatively small number of samples is a limitation of our study which may result in the identification of variants in unique candidate genes with no definitive although plausible causal effect. Larger cohorts are needed to replicate these findings and additional functional research is needed to confirm the pathogenic effects of these variants. Moreover, the sequencing coverage restricted the detection of low-level parental mosaic variants. WIDER IMPLICATIONS OF THE FINDINGS: For first-trimester euploid miscarriage, variants in unique genes may be underlying genetic etiologies and WES on trio could be an ideal model to identify potential genetic causes, which could facilitate individualized precise diagnostic and therapeutic regimens in the future. STUDY FUNDING/COMPETING INTERESTS: This study was supported by grants from the National Key Research and Development Program of China (2021YFC2700604), National Natural Science Foundation of China (31900492, 82101784, 82171648), Basic Science Center Program of the National Natural Science Foundation of China (31988101), Key Research and Development Program of Shandong Province (2021LCZX02), Natural Science Foundation of Shandong Province (ZR2020QH051), Natural Science Foundation of Jiangsu Province (BK20200223), Taishan Scholars Program for Young Experts of Shandong Province (tsqn201812154) and Young Scholars Program of Shandong University. The authors declare no conflicts of interest. TRIAL REGISTRATION NUMBER: N/A.

Diminished TMEM100 Expression in a Newborn With Acinar Dysplasia and a Novel TBX4 Variant: A Case Report.

Acinar dysplasia (AcDys) of the lung is a rare lethal developmental disorder in neonates characterized by severe respiratory failure and pulmonary arterial hypertension refractory to treatment. Recently, abnormalities of TBX4-FGF10-FGFR2-TMEM100 signaling regulating lung development have been reported in patients with AcDys due to heterozygous single-nucleotide variants or copy-number variant deletions involving TBX4, FGF10, or FGFR2. Here, we describe a female neonate who died at 4 hours of life due to severe respiratory distress related to AcDys diagnosed by postmortem histopathologic evaluation. Genomic analyses revealed a novel deleterious heterozygous missense variant c.728A>C (p.Asn243Thr) in TBX4 that arose de novo on paternal chromosome 17. We also identified 6 candidate hypomorphic rare variants in the TBX4 enhancer in trans to TBX4 coding variant. Gene expression analyses of proband's lung tissue showed a significant reduction of TMEM100 expression with near absence of TMEM100 within the endothelium of arteries and capillaries by immunohistochemistry. These results support the pathogenicity of the detected TBX4 variant and provide further evidence that disrupted signaling between TBX4 and TMEM100 may contribute to severe lung phenotypes in humans, including AcDys.

Complex Analysis of Single-Cell RNA Sequencing Data.

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

Association of Locus Xq26.1-27.3 with Erosive Tooth Wear Phenotypes in a Group of Adolescents.

Erosive tooth wear is a multifactorial condition of an increasing prevalence. There is a need for discovering individual genetic predisposition for the development of this condition. Considering that the chromosome X locus was previously shown to be associated with dental caries, the aim of the present study was to look for the association between this locus and erosive tooth wear when dietary habits are considered as a co-factor. Saliva samples, erosive wear experience data, and dietary information from 16- to 18-year-old dental patients (n = 705) were used. Genotyping analyses were performed, and thereafter, analyses considering diet and oral hygiene data, using logistic regression, with the assumption that erosive tooth wear is a complex gene-environment model. Genotypic analyses revealed an association between chromosome X marker rs1324156 and erosive tooth wear phenotype. Logistic regression analysis showed that, in the presence of less common allele of rs12687601 and rs1324156, erosive tooth wear more likely develops when associated with numerous dietary variables from the questionnaire. These results indicate that erosive tooth wear may be the result of gene-environment interactions.

Impact of single-nucleotide variants and nutritional status on population pharmacokinetics of Doxorubicin, and its effect on cardiotoxicity in children with leukemia.

PURPOSE:  Doxorubicin is an important antineoplastic agent with wide interindividual variability in response to treatment and in its cardiotoxic effects. To determine the effect of genotypic status of three single-nucleotide variants in ABCC1, NCF4, and CBR3 genes and nutritional status assessed by body mass index, on the population pharmacokinetics of Doxorubicin and its cardiotoxic effects in pediatric patients with leukemia. PATIENTS AND METHODS: Seventy pediatric patients treated with Doxorubicin were studied, in which 189 biological samples were obtained to determine Doxorubicin concentrations (1 to 3 samples per patient) at different times, for 20 h. RESULTS: Low body mass index and age ≤ 7 years were associated with decreased clearance of Doxorubicin, and female gender was associated with increased clearance of Doxorubicin. Low BMI and low height were associated with a decrease and increase, respectively, in the intercompartmental clearance (Q) of Doxorubicin. TT homozygosity of the single-nucleotide variant rs3743527 of the ABCC1 gene was associated with an increase in clearance and decreased area under the curve, AA homozygosity of the single-nucleotide variant rs1883112 of the NCF4 gene was associated with a decrease in the volume of distribution in the peripheral compartment (V2), and GG homozygosity of CBR3 rs1056892 with increasing area under the curve. CONCLUSION: Some covariates studied are directly related to the increase or decrease of the pharmacokinetic parameters of Doxorubicin. Decreased clearance, V2, and increased area under the curve were associated with systolic dysfunction, and decreased Q and V2 were associated with diastolic dysfunction. These results may contribute to the effective and safe use of Doxorubicin in pediatric patients with leukemia.

The ACE rs1799752 Variant Is Associated with COVID-19 Severity but Is Independent of Serum ACE Activity in Hospitalized and Recovered Patients.

This paper assesses the association of the insertion/deletion ACE (angiotensin-converting enzyme) variant (rs1799752 I/D) and the serum ACE activity with the severity of COVID-19 as well as its impact on post-COVID-19, and we compare these associations with those for patients with non-COVID-19 respiratory disorders. We studied 1252 patients with COVID-19, 104 subjects recovered from COVID-19, and 74 patients hospitalized with a respiratory disease different from COVID-19. The rs1799752 ACE variant was assessed using TaqMan® Assays. The serum ACE activity was determined using a colorimetric assay. The DD genotype was related to risk for invasive mechanical ventilation (IMV) requirement as an indicator of COVID-19 severity when compared to the frequencies of II + ID genotypes (p = 0.025, OR = 1.428, 95% CI = 1.046-1.949). In addition, this genotype was significantly higher in COVID-19 and post-COVID-19 groups than in the non-COVID-19 subjects. The serum ACE activity levels were lower in the COVID-19 group (22.30 U/L (13.84-32.23 U/L)), which was followed by the non-COVID-19 (27.94 U/L (20.32-53.36 U/L)) and post-COVID-19 subjects (50.00 U/L (42.16-62.25 U/L)). The DD genotype of the rs1799752 ACE variant was associated with the IMV requirement in patients with COVID-19, and low serum ACE activity levels could be related to patients with severe disease.

The rs16969968 Tobacco Smoking-Related Single-Nucleotide Variant Is Associated with Clinical Markers in Patients with Severe COVID-19.

Tobacco smoking is the leading risk factor for many respiratory diseases. Several genes are associated with nicotine addiction, such as CHRNA5 and ADAM33. This research aims to evaluate the association of the polymorphisms rs16969968 (CHRNA5) and rs3918396 (ADAM33) in patients who developed severe COVID-19. We included 917 COVID-19 patients hospitalized with critical disease and oxygenation impairment. They were divided into two groups, tobacco-smoking (n = 257) and non-smoker (n = 660) patients. The genotype and allele frequencies of two single nucleotide variants, the rs16969968 (CHRNA5) and rs3918396 (ADAM33), were evaluated. The rs3918396 in ADAM33 does not show a significative association. We analyzed the study population according to the rs16969968 genotype (GA + AA, n = 180, and GG, n = 737). The erythrocyte sedimentation rate (ESR) shows statistical differences; the GA + AA group had higher values than the GG group (p = 0.038, 32 vs. 26 mm/h, respectively). The smoking patients and GA or AA genotype carriers had a high positive correlation (p < 0.001, rho = 0.753) between fibrinogen and C-reactive protein. COVID-19 patients and smokers carriers of one or two copies of the risk allele (rs16969968/A) have high ESR and a positive correlation between fibrinogen and C-reactive protein.

Predicting the Effect of Single Mutations on Protein Stability and Binding with Respect to Types of Mutations.

The development of methods and algorithms to predict the effect of mutations on protein stability, protein-protein interaction, and protein-DNA/RNA binding is necessitated by the needs of protein engineering and for understanding the molecular mechanism of disease-causing variants. The vast majority of the leading methods require a database of experimentally measured folding and binding free energy changes for training. These databases are collections of experimental data taken from scientific investigations typically aimed at probing the role of particular residues on the above-mentioned thermodynamic characteristics, i.e., the mutations are not introduced at random and do not necessarily represent mutations originating from single nucleotide variants (SNV). Thus, the reported performance of the leading algorithms assessed on these databases or other limited cases may not be applicable for predicting the effect of SNVs seen in the human population. Indeed, we demonstrate that the SNVs and non-SNVs are not equally presented in the corresponding databases, and the distribution of the free energy changes is not the same. It is shown that the Pearson correlation coefficients (PCCs) of folding and binding free energy changes obtained in cases involving SNVs are smaller than for non-SNVs, indicating that caution should be used in applying them to reveal the effect of human SNVs. Furthermore, it is demonstrated that some methods are sensitive to the chemical nature of the mutations, resulting in PCCs that differ by a factor of four across chemically different mutations. All methods are found to underestimate the energy changes by roughly a factor of 2.