Ensembl variation resources

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

The effectiveness of selection in a species affects the direction of amino acid frequency evolution

Nearly neutral theory predicts that species with higher effective population size (N e ) are better able to purge slightly deleterious mutations. We compare evolution in high-N e vs. low-N e vertebrates to reveal which amino acid frequencies are subject to subtle selective preferences. We take three complementary approaches, two measuring flux and one measuring outcomes. First, we fit non-stationary substitution models of amino acid flux using maximum likelihood, comparing the high-N e clade of rodents and lagomorphs to its low-N e sister clade of primates and colugos. Second, we compare evolutionary outcomes across a wider range of vertebrates, via correlations between amino acid frequencies and N e . Third, we dissect the details of flux in human, chimpanzee, mouse, and rat, as scored by parsimony - this also enables comparison to a historical paper. All three methods agree on which amino acids are preferred under more effective selection. Preferred amino acids tend to be smaller, less costly to synthesize, and to promote intrinsic structural disorder. Parsimony-induced bias in the historical study produces an apparent reduction in structural disorder, perhaps driven by slightly deleterious substitutions. Within highly exchangeable pairs of amino acids, arginine is strongly preferred over lysine, and valine over isoleucine, consistent with more effective selection preferring a marginally larger free energy of folding. These two preferences match differences between thermophiles and mesophilic relatives. These results reveal the biophysical consequences of mutation-selection-drift balance, and demonstrate the utility of nearly neutral theory for understanding protein evolution.

De novo structure prediction of meteorin and meteorin-like protein for identification of domains, functional receptor binding regions, and their high-risk missense variants

Meteorin (Metrn) and Meteorin-like (Metrnl) are homologous secreted proteins involved in neural development and metabolic regulation. In this study, we have performed de novo structure prediction and analysis of both Metrn and Metrnl using Alphafold2 (AF2) and RoseTTAfold (RF). Based on the domain and structural homology analysis of the predicted structures, we have identified that these proteins are composed of two functional domains, a CUB domain and an NTR domain, connected by a hinge/loop region. We have identified the receptor binding regions of Metrn and Metrnl using the machine-learning tools ScanNet and Masif. These were further validated by docking Metrnl with its reported KIT receptor, thus establishing the role of each domain in the receptor interaction. Also, we have studied the effect of non-synonymous SNPs on the structure and function of these proteins using an array of bioinformatics tools and selected 16 missense variants in Metrn and 10 in Metrnl that can affect the protein stability. This is the first study to comprehensively characterize the functional domains of Metrn and Metrnl at their structural level and identify the functional domains, and protein binding regions. This study also highlights the interaction mechanism of the KIT receptor and Metrnl. The predicted deleterious SNPs will allow further understanding of the role of these variants in modulating the plasma levels of these proteins in disease conditions such as diabetes.Communicated by Ramaswamy H. Sarma.

Network-based approach for drug repurposing against mpox

The current outbreak of mpox presents a significant threat to the global community. However, the lack of mpox-specific drugs necessitates the identification of additional candidates for clinical trials. In this study, a network medicine framework was used to investigate poxviruses-human interactions to identify potential drugs effective against the mpox virus (MPXV). The results indicated that poxviruses preferentially target hubs on the human interactome, and that these virally-targeted proteins (VTPs) tend to aggregate together within specific modules. Comorbidity analysis revealed that mpox is closely related to immune system diseases. Based on predicted drug-target interactions, 268 drugs were identified using the network proximity approach, among which 23 drugs displaying the least side-effects and significant proximity to MPXV were selected as the final candidates. Lastly, specific drugs were explored based on VTPs, differentially expressed proteins, and intermediate nodes, corresponding to different categories. These findings provide novel insights that can contribute to a deeper understanding of the pathogenesis of MPXV and development of ready-to-use treatment strategies based on drug repurposing.

Molecular evolution of transcription factors AF4/FMR2 family member (AFF) gene family and the role of lamprey AFF3 in cell proliferation

AF4/FMR2 family member (AFF) proteins are a group of transcriptional regulators that can regulate gene transcription and play an important role in cellular physiological processes such as proliferation and differentiation. The transcriptome data of the lamprey spinal cord injury were analyzed in previous research. We then identified a hub gene, Lr-AFF3, from this dataset. Phylogenetic tree analysis determined the evolutionary relationships of the AFF gene family across different species. In addition, analysis of motifs, domains, and 3D structures further confirmed the conservatism of the AFF gene family. In particular, the gene structure of the AFF3 gene was not conserved, possibly because of intron insertion. It was also found that the neighboring genes of the Lr-AFF3 gene had a higher diversity than that in jawed vertebrates through synteny analysis. The results of the MTT and EdU experiments showed that the C-terminal homology domain (CHD) and N-terminal homology domain (NHD) of Lr-AFF3 promoted cell proliferation. In summary, our research will not only provide new insights into the origin and evolution of the AFF gene family in different species, but also provide new clues for the functions of Lr_AFF3.

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Smoking-associated gene expression alterations in nasal epithelium reveal immune impairment linked to lung cancer risk

BACKGROUND

Lung cancer is the leading cause of cancer-related death in the world. In contrast to many other cancers, a direct connection to modifiable lifestyle risk in the form of tobacco smoke has long been established. More than 50% of all smoking-related lung cancers occur in former smokers, 40% of which occur more than 15 years after smoking cessation. Despite extensive research, the molecular processes for persistent lung cancer risk remain unclear. We thus set out to examine whether risk stratification in the clinic and in the general population can be improved upon by the addition of genetic data and to explore the mechanisms of the persisting risk in former smokers.

METHODS

We analysed transcriptomic data from accessible airway tissues of 487 subjects, including healthy volunteers and clinic patients of different smoking statuses. We developed a computational model to assess smoking-associated gene expression changes and their reversibility after smoking is stopped, comparing healthy subjects to clinic patients with and without lung cancer.

RESULTS

We find persistent smoking-associated immune alterations to be a hallmark of the clinic patients. Integrating previous GWAS data using a transcriptional network approach, we demonstrate that the same immune- and interferon-related pathways are strongly enriched for genes linked to known genetic risk factors, demonstrating a causal relationship between immune alteration and lung cancer risk. Finally, we used accessible airway transcriptomic data to derive a non-invasive lung cancer risk classifier.

CONCLUSIONS

Our results provide initial evidence for germline-mediated personalized smoke injury response and risk in the general population, with potential implications for managing long-term lung cancer incidence and mortality.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia

BACKGROUND

Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs.

METHODS

We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG).

RESULTS

In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16).

CONCLUSION

The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Carrying SNP rs17506395 (T > G) in TP63 gene and CCR5Δ32 mutation associated with the occurrence of breast cancer in Burkina Faso

Genetic alterations in the TP63 (GenBank: NC_000003.12, ID: 8626) and CCR5 (receptor 5 chemokine co-receptor) (GenBank: NC_000003.12, ID: 1234) genes may increase the risk of developing breast cancer. The aim of this study was to investigate the probable involvement of polymorphisms rs17506395 in the TP63 (tumour protein 63) gene and the CCR5Δ32 mutation in the occurrence of breast cancer in Burkina Faso. This case-control study included 72 patients and 72 controls. Genotyping of SNP rs17506395 (TP63) was performed by polymerase chain reaction-restriction fragment length polymorphism, and genotyping of the CCR5Δ32 mutation was performed by allele-specific oligonucleotide polymerase chain reaction. For SNP rs17506395 (TP63), the genotypic frequencies of wild-type homozygotes (TT) and heterozygotes (TG) were, respectively, 27.72 and 72.22% in cases and 36.11 and 63.89% in controls. No mutated homozygotes (GG) were observed. For the CCR5Δ32 mutation, the genotypic frequencies of wild-type homozygotes (WT/WT) and heterozygotes (WT/Δ32) were 87.5 and 13.5%, respectively, in the cases and 89.29 and 10.71%, respectively, in the controls. No mutated homozygotes (Δ32/Δ32) were observed. None of the polymorphisms rs17506395 of the TP63 gene (OR = 1.47, 95% CI = 0.69-3.17, P = 0.284) and the CCR5Δ32 mutation (OR = 1.32, 95% CI = 0.46-3.77; P = 0.79) were associated with the occurrence of breast cancer in this study.

Progress in the relationship between GDF11 and depression

Annually, the frequency of morbidity in depression has increased progressively in response to life stressors, and there is an increasing trend toward younger morbidity. The pathogenesis of depression is complicated and includes factors such as genetic inheritance and variations in physiological functions induced by various environmental factors. Currently, drug therapy has wide adaptability in clinical practice and plays an important role in the treatment of patients with mild depression. However, the therapeutic effects of most antidepressants are typically not significant and are associated with considerable adverse effects and addiction. Therefore, it is imperative to identify the deeper mechanisms of depression and search for alternative drug targets. Growth differentiation factor 11 (GDF11) is described as an anti-ageing molecule that belongs to a member of the transforming growth factor β family. Additionally, the latest research findings suggested that GDF11 positively regulates neurogenesis and enhances neuronal activity, thereby attenuating depression-like behaviours. Although an increasing number of studies have focused on the multiple functions of GDF11 in skeletal dysplasia and carcinogenesis, its precise mechanism of action in depression remains unknown. Thus, in this review, we discuss the role of GDF11 and its mechanistic pathways in the pathogenesis of depression to develop novel therapies for depression.

Analysis of Mutational Status of IGHV, and Cytokine Polymorphisms as Prognostic Factors in Chronic Lymphocytic Leukemia: The Romanian Experience

The aim of the current study was to assess the associations between genetic risk factors (such as the mutational status of the IGHV gene and polymorphisms of the IL-10 and TNF-α genes) and CLL risk, prognosis, and overall survival. Another goal of this study was to evaluate the multivariate effect of the combination of multiple genetic risk factors (mutational status of the IGHV gene, somatic mutations, DNA CNVs, and cytokine SNPs) on the clinical characteristics and survival of patients. A total of 125 CLL patients and 239 healthy controls were included for comparative SNP analysis. IL-10 (rs1800896 and rs1800872) and TNF-α (rs361525 and rs1800750) SNPs and haplotypes were not associated with CLL risk. The absence of hypermutation in the IGHV gene was shown to be of important prognostic value, being associated with short OS. Further individual risk factors for short OS were an age above 65 years at diagnosis and the presence of somatic mutations and/or CNVs. In our multivariable analysis, the presence of somatic mutations and the IL-10 rs1800872 variant allele, and the association of CNVs with the IL-10 rs1800896 variant allele, were identified as risk factors for short OS. Moreover, the OS in unmutated IGHV patients was additionally affected (decreased) by the presence of CNVs and/or somatic mutations. Similarly, IL-10 rs1800896 modulated the OS in unmutated IGHV patients with CNVs.

Pathogen recognition pathway gene variants and inflammasome sensors gene expression in tuberculosis patients under treatment

BACKGROUND

Several epidemiological studies have suggested that genetic variations in encoding pattern recognition receptors (PRRs) genes such as Toll Like Receptors (TLRs) and their signaling products, may influence the susceptibility, severity and outcome of tuberculosis (TB). After sensing a pathogen, the cell responds producing an inflammatory response, to restrain the pathogen's successful course of infection. Herein we assessed single nucleotide polymorphisms (SNP) and gene expression from pathogen recognition and inflammasome pathways in Brazilian TB patients.

METHODS AND RESULTS

For genetic association analysis we included MYD88 and TLR4, PRRs sensing proteins. Allele distribution for MYD88 rs6853 (A > G) and TLR4 rs7873784 (C > G) presented conserved among the tested samples with statistically differential distribution in TB patients versus controls. However, when testing according to sample ethnicity (African or Caucasian-derived individuals) we identified that the rs6853 G/G genotype was associated with a lower susceptibility to TB in Caucasian population. Meanwhile, the rs7873784 G/G genotype was associated with a higher TB susceptibility in Afro-descendant ethnicity individuals. We also aimed to verify MYD88 and the inflammasome genes NLRP1 and NLRC4 expression in order to connect to active TB and/or clinical aspects.

CONCLUSIONS

We identified that inflammasome gene expression in TB patients under treatment display a similar pattern as in healthy controls, indicating that TB treatment impairs NLRP1 inflammasome activation.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

Molecular identification of missense variants in SLC3A1 gene; an approach leading to computer-aided drug design for cystinuria

Cystinuria is a rare congenital disorder characterized by the formation of cystine stones in urinary system, mainly kidneys and urinary tract. It follows the autosomal recessive inheritance pattern, where both of the parents contain the mutant allele. Cystine is an oxidized dimeric form of amino acid cysteine, shining crystal of greenish-yellow color sized greater than 5 mm. A minor genetic defect in SLC3A1 gene, downregulate the cystine transporter, rBAT protein, to absorb cystine and other dibasic amino acids in proximal tubule of nephron, causing Cystinuria. Computational and molecular analysis of SLC3A1 gene was performed to identify the deleterious missense single nucleotide variations (mSNVs) linked with Cystinuria in Pakistani population. In silico analysis of whole SLC3A1 gene nsSNPs has revealed that the exon 1, 6 and 10 are the hotspot areas, which potentially alter the protein structure and function. Three SNVs including one synonymous SNV A186C in exon 1, and two mSNVs including G314T in exon 1 and G44972A in exons 10 were identified. Both mSNVs were confirmed by ARMS PCR in all the 68 samples. The results have shown that 10% of the patients have G314T, 16% have G44972A and 74% of the patients have both of these mSNVs. Both of these mSNVs were involved in the structural and functional deterioration of rBAT protein. Additionally, computer aided drug designing tools were used to design diaminobenzylpyrimidine drug around the mutant residues which exhibit the lowest binding affinity with the target as compared to the previously reported cystine binding thiol drugs. In future, the present study could be extended to a large scale for mass screening of reported SNVs and mSNVs which, ultimately, lead to the development of knockouts for the functional studies and treatments.

Genetic variations in G-quadruplex forming sequences affect the transcription of human disease-related genes

Guanine-rich DNA strands can fold into non-canonical four-stranded secondary structures named G-quadruplexes (G4s). G4s folded in proximal promoter regions (PPR) are associated either with positive or negative transcriptional regulation. Given that single nucleotide variants (SNVs) affecting G4 folding (G4-Vars) may alter gene transcription, and that SNVs are associated with the human diseases' onset, we undertook a novel comprehensive study of the G4-Vars genome-wide (G4-variome) to find disease-associated G4-Vars located into PPRs. We developed a bioinformatics strategy to find disease-related SNVs located into PPRs simultaneously overlapping with putative G4-forming sequences (PQSs). We studied five G4-Vars disturbing in vitro the folding and stability of the G4s located into PPRs, which had been formerly associated with sporadic Alzheimer's disease (GRIN2B), a severe familiar coagulopathy (F7), atopic dermatitis (CSF2), myocardial infarction (SIRT1) and deafness (LHFPL5). Results obtained in cultured cells for these five G4-Vars suggest that the changes in the G4s affect the transcription, potentially contributing to the development of the mentioned diseases. Collectively, data reinforce the general idea that G4-Vars may impact on the different susceptibilities to human genetic diseases' onset, and could be novel targets for diagnosis and drug design in precision medicine.

Genome-Wide Association Study of Chronic Dizziness in the Elderly Identifies Loci Implicating MLLT10, BPTF, LINC01224, and ROS1

PURPOSE

Chronic age-related imbalance is a common cause of falls and subsequent death in the elderly and can arise from dysfunction of the vestibular system, an elegant neuroanatomical group of pathways that mediates human perception of acceleration, gravity, and angular head motion. Studies indicate that 27-46% of the risk of age-related chronic imbalance is genetic; nevertheless, the underlying genes remain unknown.

METHODS

The cohort consisted of 50,339 cases and 366,900 controls in the Million Veteran Program. The phenotype comprised cases with two ICD diagnoses of vertigo or dizziness at least 6 months apart, excluding acute or recurrent vertiginous syndromes and other non-vestibular disorders. Genome-wide association studies were performed as individual logistic regressions on European, African American, and Hispanic ancestries followed by trans-ancestry meta-analysis. Downstream analysis included case-case-GWAS, fine mapping, probabilistic colocalization of significant variants and genes with eQTLs, and functional analysis of significant hits.

RESULTS

Two significant loci were identified in Europeans, another in the Hispanic population, and two additional in trans-ancestry meta-analysis, including three novel loci. Fine mapping revealed credible sets of intronic single nucleotide polymorphisms (SNPs) in MLLT10 - a histone methyl transferase cofactor, BPTF - a subunit of a nucleosome remodeling complex implicated in neurodevelopment, and LINC01224 - a proto-oncogene receptor tyrosine kinase.

CONCLUSION

Despite the difficulties of phenotyping the nature of chronic imbalance, we replicated two loci from previous vertigo GWAS studies and identified three novel loci. Findings suggest candidates for further study and ultimate treatment of this common elderly disorder.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

In Silico Evaluation of Coding and Non-Coding nsSNPs in the Thrombopoietin Receptor (MPL) Proto-Oncogene: Assessing Their Influence on Protein Stability, Structure, and Function

The thrombopoietin receptor (MPL) gene is a critical regulator of hematopoiesis, and any alterations in its structure or function can result in a range of hematological disorders. Non-synonymous single nucleotide polymorphisms (nsSNPs) in MPL have the potential to disrupt normal protein function, prompting our investigation into the most deleterious MPL SNPs and the associated structural changes affecting protein-protein interactions. We employed a comprehensive suite of bioinformatics tools, including PredictSNP, InterPro, ConSurf, I-Mutant2.0, MUpro, Musitedeep, Project HOPE, STRING, RegulomeDB, Mutpred2, CScape, and CScape Somatic, to analyze 635 nsSNPs within the MPL gene. Among the analyzed nsSNPs, PredictSNP identified 28 as significantly pathogenic, revealing three critical functional domains within MPL. Ten of these nsSNPs exhibited high conservation scores, indicating potential effects on protein structure and function, while 14 were found to compromise MPL protein stability. Although the most harmful nsSNPs did not directly impact post-translational modification sites, 13 had the capacity to substantially alter the protein's physicochemical properties. Some mutations posed a risk to vital protein-protein interactions crucial for hematological functions, and three non-coding region nsSNPs displayed significant regulatory potential with potential implications for hematopoiesis. Furthermore, 13 out of 21 nsSNPs evaluated were classified as high-risk pathogenic variants by Mutpred2. Notably, amino acid alterations such as C291S, T293N, D295G, and W435C, while impactful on protein stability and function, were deemed non-oncogenic "passenger" mutations. Our study underscores the substantial impact of missense nsSNPs on MPL protein structure and function. Given MPL's central role in hematopoiesis, these mutations can significantly disrupt hematological processes, potentially leading to a variety of disorders. The identified high-risk pathogenic nsSNPs may hold promise as potential biomarkers or therapeutic targets for hematological diseases. This research lays the foundation for future investigations into the MPL gene's role in the realm of hematological health and diseases.

Activation of the transcription factor NFAT5 in the tumor microenvironment enforces CD8+ T cell exhaustion

Persistent exposure to antigen during chronic infection or cancer renders T cells dysfunctional. The molecular mechanisms regulating this state of exhaustion are thought to be common in infection and cancer, despite obvious differences in their microenvironments. Here we found that NFAT5, an NFAT family transcription factor that lacks an AP-1 docking site, was highly expressed in exhausted CD8+ T cells in the context of chronic infections and tumors but was selectively required in tumor-induced CD8+ T cell exhaustion. Overexpression of NFAT5 in CD8+ T cells reduced tumor control, while deletion of NFAT5 improved tumor control by promoting the accumulation of tumor-specific CD8+ T cells that had reduced expression of the exhaustion-associated proteins TOX and PD-1 and produced more cytokines, such as IFNɣ and TNF, than cells with wild-type levels of NFAT5, specifically in the precursor exhausted PD-1+TCF1+TIM-3-CD8+ T cell population. NFAT5 did not promote T cell exhaustion during chronic infection with clone 13 of lymphocytic choriomeningitis virus. Expression of NFAT5 was induced by TCR triggering, but its transcriptional activity was specific to the tumor microenvironment and required hyperosmolarity. Thus, NFAT5 promoted the exhaustion of CD8+ T cells in a tumor-selective fashion.

Linc2function: A Comprehensive Pipeline and Webserver for Long Non-Coding RNA (lncRNA) Identification and Functional Predictions Using Deep Learning Approaches

Long non-coding RNAs (lncRNAs), comprising a significant portion of the human transcriptome, serve as vital regulators of cellular processes and potential disease biomarkers. However, the function of most lncRNAs remains unknown, and furthermore, existing approaches have focused on gene-level investigation. Our work emphasizes the importance of transcript-level annotation to uncover the roles of specific transcript isoforms. We propose that understanding the mechanisms of lncRNA in pathological processes requires solving their structural motifs and interactomes. A complete lncRNA annotation first involves discriminating them from their coding counterparts and then predicting their functional motifs and target bio-molecules. Current in silico methods mainly perform primary-sequence-based discrimination using a reference model, limiting their comprehensiveness and generalizability. We demonstrate that integrating secondary structure and interactome information, in addition to using transcript sequence, enables a comprehensive functional annotation. Annotating lncRNA for newly sequenced species is challenging due to inconsistencies in functional annotations, specialized computational techniques, limited accessibility to source code, and the shortcomings of reference-based methods for cross-species predictions. To address these challenges, we developed a pipeline for identifying and annotating transcript sequences at the isoform level. We demonstrate the effectiveness of the pipeline by comprehensively annotating the lncRNA associated with two specific disease groups. The source code of our pipeline is available under the MIT licensefor local use by researchers to make new predictions using the pre-trained models or to re-train models on new sequence datasets. Non-technical users can access the pipeline through a web server setup.

Unveiling the immunomodulatory shift: Epithelial-mesenchymal transition Alters immune mechanisms of amniotic epithelial cells

Epithelial-mesenchymal transition (EMT) changes cell phenotype by affecting immune properties of amniotic epithelial cells (AECs). The present study shows how the response to lipopolysaccharide of cells collected pre- (eAECs) and post-EMT (mAECs) induces changes in their transcriptomics profile. In fact, eAECs mainly upregulate genes involved in antigen-presenting response, whereas mAECs over-express soluble inflammatory mediator transcripts. Consistently, network analysis identifies CIITA and Nrf2 as main drivers of eAECs and mAECs immune response, respectively. As a consequence, the depletion of CIITA and Nrf2 impairs the ability of eAECs and mAECs to inhibit lymphocyte proliferation or macrophage-dependent IL-6 release, thus confirming their involvement in regulating immune response. Deciphering the mechanisms controlling the immune function of AECs pre- and post-EMT represents a step forward in understanding key physiological events wherein these cells are involved (pregnancy and labor). Moreover, controlling the immunomodulatory properties of eAECs and mAECs may be essential in developing potential strategies for regenerative medicine applications.

MicroRNAs as powerful tool against COVID-19: Computational perspective

Severe acute respiratory syndrome coronavirus 2 is the virus that is responsible for the current pandemic, COVID-19 (SARS-CoV-2). MiRNAs, a component of RNAi technology, belong to the family of short, noncoding ssRNAs, and may be crucial in the battle against this global threat since they are involved in regulating complex biochemical pathways and may prevent viral proliferation, translation, and host expression. The complicated metabolic pathways are modulated by the activity of many proteins, mRNAs, and miRNAs working together in miRNA-mediated genetic control. The amount of omics data has increased dramatically in recent years. This massive, linked, yet complex metabolic regulatory network data offers a wealth of opportunity for iterative analysis; hence, extensive, in-depth, but time-efficient screening is necessary to acquire fresh discoveries; this is readily performed with the use of bioinformatics. We have reviewed the literature on microRNAs, bioinformatics, and COVID-19 infection to summarize (1) the function of miRNAs in combating COVID-19, and (2) the use of computational methods in combating COVID-19 in certain noteworthy studies, and (3) computational tools used by these studies against COVID-19 in several purposes. This article is categorized under: Infectious Diseases > Computational Models.

Extremely sparse models of linkage disequilibrium in ancestrally diverse association studies

Linkage disequilibrium (LD) is the correlation among nearby genetic variants. In genetic association studies, LD is often modeled using large correlation matrices, but this approach is inefficient, especially in ancestrally diverse studies. In the present study, we introduce LD graphical models (LDGMs), which are an extremely sparse and efficient representation of LD. LDGMs are derived from genome-wide genealogies; statistical relationships among alleles in the LDGM correspond to genealogical relationships among haplotypes. We published LDGMs and ancestry-specific LDGM precision matrices for 18 million common variants (minor allele frequency >1%) in five ancestry groups, validated their accuracy and demonstrated order-of-magnitude improvements in runtime for commonly used LD matrix computations. We implemented an extremely fast multiancestry polygenic prediction method, BLUPx-ldgm, which performs better than a similar method based on the reference LD correlation matrix. LDGMs will enable sophisticated methods that scale to ancestrally diverse genetic association data across millions of variants and individuals.

Multivariate genome-wide associations for immune traits in two maternal pig lines

BACKGROUND

Immune traits are considered to serve as potential biomarkers for pig's health. Medium to high heritabilities have been observed for some of the immune traits suggesting genetic variability of these phenotypes. Consideration of previously established genetic correlations between immune traits can be used to identify pleiotropic genetic markers. Therefore, genome-wide association study (GWAS) approaches are required to explore the joint genetic foundation for health biomarkers. Usually, GWAS explores phenotypes in a univariate (uv), trait-by-trait manner. Besides two uv GWAS methods, four multivariate (mv) GWAS approaches were applied on combinations out of 22 immune traits for Landrace (LR) and Large White (LW) pig lines.

RESULTS

In total 433 (LR: 351, LW: 82) associations were identified with the uv approach implemented in PLINK and a Bayesian linear regression uv approach (BIMBAM) software. Single Nucleotide Polymorphisms (SNPs) that were identified with both uv approaches (n = 32) were mostly associated with immune traits such as haptoglobin, red blood cell characteristics and cytokines, and were located in protein-coding genes. Mv GWAS approaches detected 647 associations for different mv immune trait combinations which were summarized to 133 Quantitative Trait Loci (QTL). SNPs for different trait combinations (n = 66) were detected with more than one mv method. Most of these SNPs are associated with red blood cell related immune trait combinations. Functional annotation of these QTL revealed 453 immune-relevant protein-coding genes. With uv methods shared markers were not observed between the breeds, whereas mv approaches were able to detect two conjoint SNPs for LR and LW. Due to unmapped positions for these markers, their functional annotation was not clarified.

CONCLUSIONS

This study evaluated the joint genetic background of immune traits in LR and LW piglets through the application of various uv and mv GWAS approaches. In comparison to uv methods, mv methodologies identified more significant associations, which might reflect the pleiotropic background of the immune system more accurately. In genetic research of complex traits, the SNP effects are generally small. Furthermore, one genetic variant can affect several correlated immune traits at the same time, termed pleiotropy. As mv GWAS methods consider strong dependencies among traits, the power to detect SNPs can be boosted. Both methods revealed immune-relevant potential candidate genes. Our results indicate that one single test is not able to detect all the different types of genetic effects in the most powerful manner and therefore, the methods should be applied complementary.

Genetics of hair graying with age

Hair graying is an early and obvious phenotypic and physiological trait with age in humans. Several recent advances in molecular biology and genetics have increased our understanding of the mechanisms of hair graying, which elucidate genes related to the synthesis, transport, and distribution of melanin in hair follicles, as well as genes regulating these processes above. Therefore, we review these advances and examine the trends in the genetic aspects of hair graying from enrichment theory, Genome-Wide association studies, whole exome sequencing, gene expression studies, and animal models for hair graying with age, aiming to overview the changes in hair graying at the genetic level and establish the foundation for future research. Meanwhile, by summarizing the genetics, it's of great value to explore the possible mechanism, treatment, or even prevention of hair graying with age.

Shared hotspot mutations in oncogenes position dogs as an unparalleled comparative model for precision therapeutics

Naturally occurring canine cancers have remarkable similarities to their human counterparts. To better understand these similarities, we investigated 671 client-owned dogs from 96 breeds with 23 common tumor types, including those whose mutation profile are unknown (anal sac carcinoma and neuroendocrine carcinoma) or understudied (thyroid carcinoma, soft tissue sarcoma and hepatocellular carcinoma). We discovered mutations in 50 well-established oncogenes and tumor suppressors, and compared them to those reported in human cancers. As in human cancer, TP53 is the most commonly mutated gene, detected in 22.5% of canine tumors overall. Canine tumors share mutational hotspots with human tumors in oncogenes including PIK3CA, KRAS, NRAS, BRAF, KIT and EGFR. Hotspot mutations with significant association to tumor type include NRAS G61R and PIK3CA H1047R in hemangiosarcoma, ERBB2 V659E in pulmonary carcinoma, and BRAF V588E (equivalent of V600E in humans) in urothelial carcinoma. Our findings better position canines as a translational model of human cancer to investigate a wide spectrum of targeted therapies.

Extensive Bioinformatics Analyses Reveal a Phylogenetically Conserved Winged Helix (WH) Domain (Zτ) of Topoisomerase IIα, Elucidating Its Very High Affinity for Left-Handed Z-DNA and Suggesting Novel Putative Functions

The dynamic processes operating on genomic DNA, such as gene expression and cellular division, lead inexorably to topological challenges in the form of entanglements, catenanes, knots, "bubbles", R-loops, and other outcomes of supercoiling and helical disruption. The resolution of toxic topological stress is the function attributed to DNA topoisomerases. A prominent example is the negative supercoiling (nsc) trailing processive enzymes such as DNA and RNA polymerases. The multiple equilibrium states that nscDNA can adopt by redistribution of helical twist and writhe include the left-handed double-helical conformation known as Z-DNA. Thirty years ago, one of our labs isolated a protein from Drosophila cells and embryos with a 100-fold greater affinity for Z-DNA than for B-DNA, and identified it as topoisomerase II (gene Top2, orthologous to the human UniProt proteins TOP2A and TOP2B). GTP increased the affinity and selectivity for Z-DNA even further and also led to inhibition of the isomerase enzymatic activity. An allosteric mechanism was proposed, in which topoII acts as a Z-DNA-binding protein (ZBP) to stabilize given states of topological (sub)domains and associated multiprotein complexes. We have now explored this possibility by comprehensive bioinformatic analyses of the available protein sequences of topoII representing organisms covering the whole tree of life. Multiple alignment of these sequences revealed an extremely high level of evolutionary conservation, including a winged-helix protein segment, here denoted as Zτ, constituting the putative structural homolog of Zα, the canonical Z-DNA/Z-RNA binding domain previously identified in the interferon-inducible RNA Adenosine-to-Inosine-editing deaminase, ADAR1p150. In contrast to Zα, which is separate from the protein segment responsible for catalysis, Zτ encompasses the active site tyrosine of topoII; a GTP-binding site and a GxxG sequence motif are in close proximity. Quantitative Zτ-Zα similarity comparisons and molecular docking with interaction scoring further supported the "B-Z-topoII hypothesis" and has led to an expanded mechanism for topoII function incorporating the recognition of Z-DNA segments ("Z-flipons") as an inherent and essential element. We further propose that the two Zτ domains of the topoII homodimer exhibit a single-turnover "conformase" activity on given G(ate) B-DNA segments ("Z-flipins"), inducing their transition to the left-handed Z-conformation. Inasmuch as the topoII-Z-DNA complexes are isomerase inactive, we infer that they fulfill important structural roles in key processes such as mitosis. Topoisomerases are preeminent targets of anti-cancer drug discovery, and we anticipate that detailed elucidation of their structural-functional interactions with Z-DNA and GTP will facilitate the design of novel, more potent and selective anti-cancer chemotherapeutic agents.

RNA-Seq-based discovery of genetic variants and allele-specific expression of two layer lines and broiler chicken

Recent advances in the selective breeding of broilers and layers have made poultry production one of the fastest-growing industries. In this study, a transcriptome variant calling approach from RNA-seq data was used to determine population diversity between broilers and layers. In total, 200 individuals were analyzed from three different chicken populations (Lohmann Brown (LB), n = 90), Lohmann Selected Leghorn (LSL, n = 89), and Broiler (BR, n = 21). The raw RNA-sequencing reads were pre-processed, quality control checked, mapped to the reference genome, and made compatible with Genome Analysis ToolKit for variant detection. Subsequently, pairwise fixation index (F ST) analysis was performed between broilers and layers. Numerous candidate genes were identified, that were associated with growth, development, metabolism, immunity, and other economically significant traits. Finally, allele-specific expression (ASE) analysis was performed in the gut mucosa of LB and LSL strains at 10, 16, 24, 30, and 60 weeks of age. At different ages, the two-layer strains showed significantly different allele-specific expressions in the gut mucosa, and changes in allelic imbalance were observed across the entire lifespan. Most ASE genes are involved in energy metabolism, including sirtuin signaling pathways, oxidative phosphorylation, and mitochondrial dysfunction. A high number of ASE genes were found during the peak of laying, which were particularly enriched in cholesterol biosynthesis. These findings indicate that genetic architecture as well as biological processes driving particular demands relate to metabolic and nutritional requirements during the laying period shape allelic heterogeneity. These processes are considerably affected by breeding and management, whereby elucidating allele-specific gene regulation is an essential step towards deciphering the genotype to phenotype map or functional diversity between the chicken populations. Additionally, we observed that several genes showing significant allelic imbalance also colocalized with the top 1% of genes identified by the FST approach, suggesting a fixation of genes in cis-regulatory elements.

Comprehensive behavioral analyses of mice with a glycine receptor alpha 4 deficiency

Glycine receptors (GlyRs) are ligand-gated chloride channels comprising alpha (α1-4) and β subunits. The GlyR subunits play major roles in the mammalian central nervous system, ranging from regulating simple sensory information to modulating higher-order brain function. Unlike the other GlyR subunits, GlyR α4 receives relatively little attention because the human ortholog lacks a transmembrane domain and is thus considered a pseudogene. A recent genetic study reported that the GLRA4 pseudogene locus on the X chromosome is potentially involved in cognitive impairment, motor delay and craniofacial anomalies in humans. The physiologic roles of GlyR α4 in mammal behavior and its involvement in disease, however, are not known. Here we examined the temporal and spatial expression profile of GlyR α4 in the mouse brain and subjected Glra4 mutant mice to a comprehensive behavioral analysis to elucidate the role of GlyR α4 in behavior. The GlyR α4 subunit was mainly enriched in the hindbrain and midbrain, and had relatively lower expression in the thalamus, cerebellum, hypothalamus, and olfactory bulb. In addition, expression of the GlyR α4 subunit gradually increased during brain development. Glra4 mutant mice exhibited a decreased amplitude and delayed onset of the startle response compared with wild-type littermates, and increased social interaction in the home cage during the dark period. Glra4 mutants also had a low percentage of entries into open arms in the elevated plus-maze test. Although mice with GlyR α4 deficiency did not show motor and learning abnormalities reported to be associated in human genomics studies, they exhibited behavioral changes in startle response and social and anxiety-like behavior. Our data clarify the spatiotemporal expression pattern of the GlyR α4 subunit and suggest that glycinergic signaling modulates social, startle, and anxiety-like behaviors in mice.

Comprehensive In Silico Characterization of the Coding and Non-Coding SNPs in Human Dectin-1 Gene with the Potential of High-Risk Pathogenicity Associated with Fungal Infections

The human C-type lectin domain family 7 member A (CLEC7A) gene encodes a Dectin-1 protein that recognizes beta-1,3-linked and beta-1,6-linked glucans, which form the cell walls of pathogenic bacteria and fungi. It plays a role in immunity against fungal infections through pathogen recognition and immune signaling. This study aimed to explore the impact of nsSNPs in the human CLEC7A gene through computational tools (MAPP, PhD-SNP, PolyPhen-1, PolyPhen-2, SIFT, SNAP, and PredictSNP) to identify the most deleterious and damaging nsSNPs. Further, their effect on protein stability was checked along with conservation and solvent accessibility analysis by I-Mutant 2.0, ConSurf, and Project HOPE and post-translational modification analysis using MusiteDEEP. Out of the 28 nsSNPs that were found to be deleterious, 25 nsSNPs affected protein stability. Some SNPs were finalized for structural analysis with Missense 3D. Seven nsSNPs affected protein stability. Results from this study predicted that C54R, L64P, C120G, C120S, S135C, W141R, W141S, C148G, L155P, L155V, I158M, I158T, D159G, D159R, I167T, W180R, L183F, W192R, G197E, G197V, C220S, C233Y, I240T, E242G, and Y3D were the most structurally and functionally significant nsSNPs in the human CLEC7A gene. No nsSNPs were found in the predicted sites for post-translational modifications. In the 5' untranslated region, two SNPs, rs536465890 and rs527258220, showed possible miRNA target sites and DNA binding sites. The present study identified structurally and functionally significant nsSNPs in the CLEC7A gene. These nsSNPs may potentially be used for further evaluation as diagnostic and prognostic biomarkers.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Whole exome sequencing analysis of canine urothelial carcinomas without BRAF V595E mutation: Short in-frame deletions in BRAF and MAP2K1 suggest alternative mechanisms for MAPK pathway disruption

Molecular profiling studies have shown that 85% of canine urothelial carcinomas (UC) harbor an activating BRAF V595E mutation, which is orthologous to the V600E variant found in several human cancer subtypes. In dogs, this mutation provides both a powerful diagnostic marker and a potential therapeutic target; however, due to their relative infrequency, the remaining 15% of cases remain understudied at the molecular level. We performed whole exome sequencing analysis of 28 canine urine sediments exhibiting the characteristic DNA copy number signatures of canine UC, in which the BRAF V595E mutation was undetected (UDV595E specimens). Among these we identified 13 specimens (46%) harboring short in-frame deletions within either BRAF exon 12 (7/28 cases) or MAP2K1 exons 2 or 3 (6/28 cases). Orthologous variants occur in several human cancer subtypes and confer structural changes to the protein product that are predictive of response to different classes of small molecule MAPK pathway inhibitors. DNA damage response and repair genes, and chromatin modifiers were also recurrently mutated in UDV595E specimens, as were genes that are positive predictors of immunotherapy response in human cancers. Our findings suggest that short in-frame deletions within BRAF exon 12 and MAP2K1 exons 2 and 3 in UDV595E cases are alternative MAPK-pathway activating events that may have significant therapeutic implications for selecting first-line treatment for canine UC. We developed a simple, cost-effective capillary electrophoresis genotyping assay for detection of these deletions in parallel with the BRAF V595E mutation. The identification of these deletion events in dogs offers a compelling cross-species platform in which to study the relationship between somatic alteration, protein conformation, and therapeutic sensitivity.

RNA-sequencing revealed apple pomace ameliorates expression of genes in the hypothalamus associated with neurodegeneration in female rats fed a Western diet during adolescence to adulthood

OBJECTIVES

Apple pomace, a waste byproduct of apple processing, is rich in nutrients (e.g. polyphenols and soluble fiber) with the potential to be neuroprotective. The aim of this study was to employ RNA-sequencing (RNASeq) technology to investigate diet-gene interactions in the hypothalamus of rats after feeding a Western diet calorically substituted with apple pomace.

METHODS

Adolescent (age 21-29 days) female Sprague-Dawley rats were randomly assigned (n = 8 rats/group) to consume either a purified standard diet, Western (WE) diet, or Western diet calorically substituted with 10% apple pomace (WE/AP) for 8 weeks. RNA-seq was performed (n = 5 rats/group) to determine global differentially expressed genes in the hypothalamus.

RESULTS

RNA-seq results comparing rats fed WE to WE/AP revealed 15 differentially expressed genes in the hypothalamus. Caloric substitution of WE diet with 10% apple pomace downregulated (q < 0.06) five genes implicated in brain aging and neurodegenerative disorders: synuclein alpha, phospholipase D family member 5, NADH dehydrogenase Fe-S protein 6, choline O-acetyltransferase, and frizzled class receptor 6.

DISCUSSION

Altered gene expression of these five genes suggests that apple pomace ameliorated synthesis of the neurotransmitter, acetylcholine, in rats fed a WE diet. Apple pomace, a rich source of antioxidant polyphenols and soluble fiber, has been shown to reverse nonalcoholic fatty liver disease (NAFLD). Diet-induced NAFLD decreases hepatic de novo synthesis of choline, a precursor to acetylcholine. Based on preclinical evidence, apple pomace has the potential to be a sustainable functional food for maintaining brain function and for reducing the risk of neurodegeneration.

In silico exploration and molecular dynamics of deleterious SNPs on the human TERF1 protein triggering male infertility

By limiting chromosome erosion and end-to-end fusions, telomere integrity is critical for chromosome stability and cell survival. During mitotic cycles or due to environmental stresses, telomeres become progressively shorter and dysfunctional, thus triggering cellular senescence, genomic instability and cell death. To avoid such consequences, the telomerase action, as well as the Shelterin and CST complexes, assure the telomere's protection. Telomeric repeat binding factor 1 (TERF1), which is one of the primary components of the Shelterin complex, binds directly to the telomere and controls its length and function by regulating the telomerase activity. Several reports about TERF1 gene variations have been associated with different diseases, and some of them have linked these variations to male infertility. Hence, this paper can be advantageous to investigate the association between the missense variants of the TERF1 gene and the susceptibility to male infertility. The stepwise prediction of SNPs pathogenicity followed in this study was based on stability and conservation analysis, post-translational modification, secondary structure, functional interaction prediction, binding energy evaluation and finally molecular dynamic simulation. Prediction matching among the tools revealed that out of 18 SNPs, only four (rs1486407144, rs1259659354, rs1257022048 and rs1320180267) were predicted as the most damaging and highly deleterious SNPs affecting the TERF1 protein and its molecular dynamics when interacting with the TERB1 protein by influencing the function, structural stability, flexibility and compaction of the overall complex. Interestingly, these polymorphisms should be considered during genetic screening so they can be used effectively as genetic biomarkers for male infertility diagnosis.Communicated by Ramaswamy H. Sarma.

Genetic risk of depression is different in subgroups of dietary ratio of tryptophan to large neutral amino acids

Manipulation of intake of serotonin precursor tryptophan has been exploited to rapidly induce and alleviate depression symptoms. While studies show that this latter effect is dependent on genetic vulnerability to depression, the effect of habitual tryptophan intake in the context of predisposing genetic factors has not been explored. Our aim was to investigate the effect of habitual tryptophan intake on mood symptoms and to determine the effect of risk variants on depression in those with high and low tryptophan intake in the whole genome and specifically in serotonin and kynurenine pathways. 63,277 individuals in the UK Biobank with data on depressive symptoms and tryptophan intake were included. We compared two subpopulations defined by their habitual diet of a low versus a high ratio of tryptophan to other large amino acids (TLR). A modest protective effect of high dietary TLR against depression was found. NPBWR1 among serotonin genes and POLI in kynurenine pathway genes were significantly associated with depression in the low but not in the high TLR group. Pathway-level analyses identified significant associations for both serotonin and kynurenine pathways only in the low TLR group. In addition, significant association was found in the low TLR group between depressive symptoms and biological process related to adult neurogenesis. Our findings demonstrate a markedly distinct genetic risk profile for depression in groups with low and high dietary TLR, with association with serotonin and kynurenine pathway variants only in case of habitual food intake leading to low TLR. Our results confirm the relevance of the serotonin hypothesis in understanding the neurobiological background of depression and highlight the importance of understanding its differential role in the context of environmental variables such as complexity of diet in influencing mental health, pointing towards emerging possibilities of personalised prevention and intervention in mood disorders in those who are genetically vulnerable.

Integrated-gut-liver-on-a-chip platform as an in vitro human model of non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) afflicts a significant percentage of the population; however, no effective treatments have yet been established because of the unsuitability of in vitro assays and animal experimental models. Here, we present an integrated-gut-liver-on-a-chip (iGLC) platform as an in vitro human model of the gut-liver axis (GLA) by co-culturing human gut and liver cell lines interconnected via microfluidics in a closed circulation loop, for the initiation and progression of NAFLD by treatment with free fatty acids (FFAs) for 1 and 7 days, respectively. Co-cultured Caco-2 gut-mimicking cells and HepG2 hepatocyte-like cells demonstrate the protective effects from apoptosis against FFAs treatment, whereas mono-cultured cells exhibit induced apoptosis. Phenotype and gene expression analyses reveal that the FFAs-treated gut and liver cells accumulated intracellular lipid droplets and show an increase in gene expression associated with a cellular response to copper ions and endoplasmic reticulum stress. As an in vitro human GLA model, the iGLC platform may serve as an alternative to animal experiments for investigating the mechanisms of NAFLD.

Identification of Three Human POLH Germline Variants Defective in Complementing the UV- and Cisplatin-Sensitivity of POLH-Deficient Cells

DNA polymerase (pol) η is responsible for error-free translesion DNA synthesis (TLS) opposite ultraviolet light (UV)-induced cis-syn cyclobutane thymine dimers (CTDs) and cisplatin-induced intrastrand guanine crosslinks. POLH deficiency causes one form of the skin cancer-prone disease xeroderma pigmentosum variant (XPV) and cisplatin sensitivity, but the functional impacts of its germline variants remain unclear. We evaluated the functional properties of eight human POLH germline in silico-predicted deleterious missense variants, using biochemical and cell-based assays. In enzymatic assays, utilizing recombinant pol η (residues 1-432) proteins, the C34W, I147N, and R167Q variants showed 4- to 14-fold and 3- to 5-fold decreases in specificity constants (kcat/Km) for dATP insertion opposite the 3'-T and 5'-T of a CTD, respectively, compared to the wild-type, while the other variants displayed 2- to 4-fold increases. A CRISPR/Cas9-mediated POLH knockout increased the sensitivity of human embryonic kidney 293 cells to UV and cisplatin, which was fully reversed by ectopic expression of wild-type pol η, but not by that of an inactive (D115A/E116A) or either of two XPV-pathogenic (R93P and G263V) mutants. Ectopic expression of the C34W, I147N, and R167Q variants, unlike the other variants, did not rescue the UV- and cisplatin-sensitivity in POLH-knockout cells. Our results indicate that the C34W, I147N, and R167Q variants-substantially reduced in TLS activity-failed to rescue the UV- and cisplatin-sensitive phenotype of POLH-deficient cells, which also raises the possibility that such hypoactive germline POLH variants may increase the individual susceptibility to UV irradiation and cisplatin chemotherapy.

Evolutionary History of the Poecilia picta Sex Chromosomes

The degree of divergence between the sex chromosomes is not always proportional to their age. In poeciliids, four closely related species all exhibit a male heterogametic sex chromosome system on the same linkage group, yet show a remarkable diversity in X and Y divergence. In Poecilia reticulata and P. wingei, the sex chromosomes remain homomorphic, yet P. picta and P. parae have a highly degraded Y chromosome. To test alternative theories about the origin of their sex chromosomes, we used a combination of pedigrees and RNA-seq data from P. picta families in conjunction with DNA-seq data collected from P. reticulata, P. wingei, P. parae, and P. picta. Phylogenetic clustering analysis of X and Y orthologs, identified through segregation patterns, and their orthologous sequences in closely related species demonstrates a similar time of origin for both the P. picta and P. reticulata sex chromosomes. We next used k-mer analysis to identify shared ancestral Y sequence across all four species, suggesting a single origin to the sex chromosome system in this group. Together, our results provide key insights into the origin and evolution of the poeciliid Y chromosome and illustrate that the rate of sex chromosome divergence is often highly heterogenous, even over relatively short evolutionary time frames.

The first genetic landscape of inherited retinal dystrophies in Portuguese patients identifies recurrent homozygous mutations as a frequent cause of pathogenesis

Inherited retinal diseases (IRDs) are a group of ocular conditions characterized by an elevated genetic and clinical heterogeneity. They are transmitted almost invariantly as monogenic traits. However, with more than 280 disease genes identified so far, association of clinical phenotypes with genotypes can be very challenging, and molecular diagnosis is essential for genetic counseling and correct management of the disease. In addition, the prevalence and the assortment of IRD mutations are often population-specific. In this work, we examined 230 families from Portugal, with individuals suffering from a variety of IRD diagnostic classes (270 subjects in total). Overall, we identified 157 unique mutations (34 previously unreported) in 57 distinct genes, with a diagnostic rate of 76%. The IRD mutational landscape was, to some extent, different from those reported in other European populations, including Spanish cohorts. For instance, the EYS gene appeared to be the most frequently mutated, with a prevalence of 10% among all IRD cases. This was, in part, due to the presence of a recurrent and seemingly founder mutation involving the deletion of exons 13 and 14 of this gene. Moreover, our analysis highlighted that as many as 51% of our cases had mutations in a homozygous state. To our knowledge, this is the first study assessing a cross-sectional genotype-phenotype landscape of IRDs in Portugal. Our data reveal a rather unique distribution of mutations, possibly shaped by a small number of rare ancestral events that have now become prevalent alleles in patients.

Genome-wide genotype-serum proteome mapping provides insights into the cross-ancestry differences in cardiometabolic disease susceptibility

Identification of protein quantitative trait loci (pQTL) helps understand the underlying mechanisms of diseases and discover promising targets for pharmacological intervention. For most important class of drug targets, genetic evidence needs to be generalizable to diverse populations. Given that the majority of the previous studies were conducted in European ancestry populations, little is known about the protein-associated genetic variants in East Asians. Based on data-independent acquisition mass spectrometry technique, we conduct genome-wide association analyses for 304 unique proteins in 2,958 Han Chinese participants. We identify 195 genetic variant-protein associations. Colocalization and Mendelian randomization analyses highlight 60 gene-protein-phenotype associations, 45 of which (75%) have not been prioritized in Europeans previously. Further cross-ancestry analyses uncover key proteins that contributed to the differences in the obesity-induced diabetes and coronary artery disease susceptibility. These findings provide novel druggable proteins as well as a unique resource for the trans-ancestry evaluation of protein-targeted drug discovery.

Regulatory dissection of the severe COVID-19 risk locus introgressed by Neanderthals

Individuals infected with the SARS-CoV-2 virus present with a wide variety of symptoms ranging from asymptomatic to severe and even lethal outcomes. Past research has revealed a genetic haplotype on chromosome 3 that entered the human population via introgression from Neanderthals as the strongest genetic risk factor for the severe response to COVID-19. However, the specific variants along this introgressed haplotype that contribute to this risk and the biological mechanisms that are involved remain unclear. Here, we assess the variants present on the risk haplotype for their likelihood of driving the genetic predisposition to severe COVID-19 outcomes. We do this by first exploring their impact on the regulation of genes involved in COVID-19 infection using a variety of population genetics and functional genomics tools. We then perform a locus-specific massively parallel reporter assay to individually assess the regulatory potential of each allele on the haplotype in a multipotent immune-related cell line. We ultimately reduce the set of over 600 linked genetic variants to identify four introgressed alleles that are strong functional candidates for driving the association between this locus and severe COVID-19. Using reporter assays in the presence/absence of SARS-CoV-2, we find evidence that these variants respond to viral infection. These variants likely drive the locus' impact on severity by modulating the regulation of two critical chemokine receptor genes: CCR1 and CCR5. These alleles are ideal targets for future functional investigations into the interaction between host genomics and COVID-19 outcomes.

SCAF4-related syndromic intellectual disability

The causal link between variants in the SCAF4 gene and a syndromic form of intellectual disability (ID) was established in 2020 by Fliedner et al. Since then, no additional cases have been reported. We performed exome sequencing in a 16-year-old Brazilian male presenting with ID, epilepsy, behavioral problems, speech impairment, facial dysmorphisms, heart malformations, and obesity. A de novo pathogenic variant [SCAF4(NM_020706.2):c.374_375dup(p.Glu126LeufsTer20)] was identified. This is the second study reporting the involvement of SCAF4 in syndromic ID, and the description of the patient's clinical features contributes to defining the phenotypic spectrum of this recently described Mendelian disorder.

High-throughput sequencing analysis of nuclear-encoded mitochondrial genes reveals a genetic signature of human longevity

Mitochondrial dysfunction is a well-known contributor to aging and age-related diseases. The precise mechanisms through which mitochondria impact human lifespan, however, remain unclear. We hypothesize that humans with exceptional longevity harbor rare variants in nuclear-encoded mitochondrial genes (mitonuclear genes) that confer resistance against age-related mitochondrial dysfunction. Here we report an integrated functional genomics study to identify rare functional variants in ~ 660 mitonuclear candidate genes discovered by target capture sequencing analysis of 496 centenarians and 572 controls of Ashkenazi Jewish descent. We identify and prioritize longevity-associated variants, genes, and mitochondrial pathways that are enriched with rare variants. We provide functional gene variants such as those in MTOR (Y2396Lfs*29), CPS1 (T1406N), and MFN2 (G548*) as well as LRPPRC (S1378G) that is predicted to affect mitochondrial translation. Taken together, our results suggest a functional role for specific mitonuclear genes and pathways in human longevity.

Genomic health is dependent on long-term population demographic history

Current genetic methods of population assessment in conservation biology have been challenged by genome-scale analyses due to their quantitatively novel insights. These analyses include assessments of runs-of-homozygosity (ROH), genomic evolutionary rate profiling (GERP), and mutational load. Here, we aim to elucidate the relationships between these measures using three divergent ungulates: white-tailed deer, caribou, and mountain goat. The white-tailed deer is currently expanding, while caribou are in the midst of a significant decline. Mountain goats remain stable, having suffered a large historical bottleneck. We assessed genome-wide signatures of inbreeding using the inbreeding coefficient F and %ROH (F_ROH ) and identified evolutionarily constrained regions with GERP. Mutational load was estimated by identifying mutations in highly constrained elements (CEs) and sorting intolerant from tolerant (SIFT) mutations. Our results showed that F and F_ROH are higher in mountain goats than in caribou and white-tailed deer. Given the extended bottleneck and low N_e of the mountain goat, this supports the idea that the genome-wide effects of demographic change take time to accrue. Similarly, we found that mountain goats possess more highly constrained CEs and the lowest dN/dS values, both of which are indicative of greater purifying selection; this is also reflected by fewer mutations in CEs and deleterious mutations identified by SIFT. In contrast, white-tailed deer presented the highest mutational load with both metrics, in addition to dN/dS, while caribou were intermediate. Our results demonstrate that extended bottlenecks may lead to reduced diversity and increased F_ROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations.

Multiplexed target enrichment of coding and non-coding transcriptomes enables studying Candida spp. infections from human derived samples

The study of transcriptomic interactions between host and pathogens in in vivo conditions is challenged by the low relative amounts of the pathogen RNA. Yeast opportunistic pathogens of the genus Candida can cause life-threatening systemic infections in immunocompromised patients, and are of growing medical concern. Four phylogenetically diverse species account for over 90% of Candida infections, and their specific interactions with various human tissues are still poorly understood. To enable in vivo transcriptomic analysis in these species, we designed and validated pan-Candida target capture probes to enrich protein-coding and non-coding transcriptomes. The probe-based enrichment approach outperformed enrichment based on differential lysis of host cells, and showed similar enrichment performance as an existing capture design, yet achieving better fidelity of expression levels, enabling species multiplexing and capturing of lncRNAs. In addition, we show that our probe-based enrichment strategy allows robust genotype-based identification of the infecting strain present in the sample.

Molecular characterization and expression analysis of a novel cold-inducible RNA-binding protein (CIRBP) gene in lamprey (Lethenteron reissneri)

Cold-inducible RNA-binding protein (CIRBP) responds to a wide array of cellular stresses such as cold shock, hypoxia, and inflammatory responses. However, functional studies of CIRBP in jawless vertebrates are limited. In this study, a CIRBP homolog from the jawless vertebrate lamprey (Lethenteron reissneri) was cloned and characterized (named Lr-CIRBP). The cDNA fragment of Lr-CIRBP has a 516 bp open reading frame (ORF) that encodes 171 amino acids, comprising a glycine-rich region at the C-terminal, similar to higher vertebrates but slightly shorter, and an RNA recognition motif (RRM) domain at the N-terminus. The predicted Lr-CIRBP sequence had 51.4 ~ 70.6% similarity with CIRBPs from other vertebrates. Further phylogenetic analysis revealed that Lr-CIRBP is located in the outgroup of vertebrates and is the ancestor of vertebrates. Based on real-time quantitative PCR experimental analysis, Lr-CIRBP expression was highest in leukocytes and increased significantly after multi-stimulation, peaking at 12 h. RNA interference showed that Lr-CIRBP knockdown can down-regulate the expression of inflammatory factors in Lethenteron reissneri. In conclusion, our study successfully clarifies the ancestral features and functions of CIRBP, while revealing valuable insight into how the protein is involved in the immune responses of a jawless vertebrate.

Evolution of TOP1 and TOP1MT Topoisomerases in Chordata

Type IB topoisomerases relax the torsional stress associated with DNA metabolism in the nucleus and mitochondria and constitute important molecular targets of anticancer drugs. Vertebrates stand out among eukaryotes by having two Type IB topoisomerases acting specifically in the nucleus (TOP1) and mitochondria (TOP1MT). Despite their major importance, the origin and evolution of these paralogues remain unknown. Here, we examine the molecular evolutionary processes acting on both TOP1 and TOP1MT in Chordata, taking advantage of the increasing number of available genome sequences. We found that both TOP1 and TOP1MT evolved under strong purifying selection, as expected considering their essential biological functions. Critical active sites, including those associated with resistance to anticancer agents, were found particularly conserved. However, TOP1MT presented a higher rate of molecular evolution than TOP1, possibly related with its specialized activity on the mitochondrial genome and a less critical role in cells. We could place the duplication event that originated the TOP1 and TOP1MT paralogues early in the radiation of vertebrates, most likely associated with the first round of vertebrate tetraploidization (1R). Moreover, our data suggest that cyclostomes present a specialized mitochondrial Type IB topoisomerase. Interestingly, we identified two missense mutations replacing amino acids in the Linker region of TOP1MT in Neanderthals, which appears as a rare event when comparing the genome of both species. In conclusion, TOP1 and TOP1MT differ in their rates of evolution, and their evolutionary histories allowed us to better understand the evolution of chordates.

LncRNA RP11-551L14.4 suppresses breast cancer development by inhibiting the expression of miR-4472

Background

Previous studies have been reported that long non-coding RNA (lncRNA) can regulate the expression of genes which are involved in many important cellular processes The potential role of lncRNA RP11-551L14.4 in the development of breast cancer and the possible regulatory mechanisms was investigated.

Methods

Quantitative real-time polymerase chain reaction (qRT-PCR) was conducted to analyze RP11-551L14.4 expression in 36 paired breast cancer tissues and adjacent tissues. The expression of RP11-551L14.4 in multiple breast cancer cell lines was detected by qRT-PCR. Meanwhile, overexpression of RP11-551L14.4 models was established using lentivirus in BT474 and T47D breast cancer cells. Cell counting kit-8 (CCK-8), cell colony formation and cell cycle assays were performed to detect the effects of RP11-551L14.4 on the biological function of breast cancer cells. Besides, bioinformatics techniques, dual luciferase reporter gene assay and rescue experiments were used to investigate the potential mechanisms.

Results

RP11-551L14.4 expression was negatively associated with the advanced tumor stage. Breast cancer patients with low RP11-551L14.4 expression manifested a poorer prognosis. The results of qRT-PCR showed that RP11-551L14.4 expression in breast cancer tissues was significantly lower than in adjacent tissues. Meanwhile, overexpression of RP11-551L14.4 significantly decreased the cell proliferation and cell cycle. Bioinformatics technology showed that RP11-551L14.4 could complementarily bind to miR-4472. qRT-PCR results indicated that the expression levels of miR-4472 and RP11-551L14.4 in breast cancer were negatively correlated. Luciferase reporter gene assay showed that miR-4472 remarkably decreased the relative luciferase activity of the wild-type RP11-551L14.4 vector. miR-4472 is a direct target gene of RP11-551L14.4. miR-4472 levels were reduced, and repulsive guidance molecule A (RGMA) mRNA or protein levels were increased after overexpression of RP11-551L14.4 in the breast cancer cells. miR-4472 reversed the effects caused by RP11-551L14.4 in breast cancer cells.

Conclusion

RP11-551L14.4 expression was remarkably decreased in breast cancer tissues and cells. RP11-551L14.4 may inhibit the malignant progression of breast cancer by regulating miR-4472 expression.

Structural Evaluation and Conformational Dynamics of ZNF141T474I Mutation Provoking Postaxial Polydactyly Type A

Postaxial Polydactyly (PAP) is a congenital disorder of limb abnormalities characterized by posterior extra digits. Mutations in the N-terminal region of the Zinc finger protein 141 (ZNF141) gene were recently linked with PAP type A. Zinc finger proteins exhibit similarity at their N-terminal regions due to C2-H2 type Zinc finger domains, but their functional preferences vary significantly by the binding patterns of DNA. Methods: This study delineates the pathogenic association, miss-fold aggregation, and conformational paradigm of a missense variant (c.1420C > T; p.T474I) in ZNF141 gene segregating PAP through a molecular dynamics simulations approach. Results: In ZNF141 protein, helices play a crucial role by attaching three specific target DNA base pairs. In ZNF141T474I protein, H1, H3, and H6 helices attain more flexibility by acquiring loop conformation. The outward disposition of the proximal portion of H9-helix in mutant protein occurs due to the loss of prior beta-hairpins at the C terminal region of the C2-H2 domain. The loss of hydrogen bonds and exposure of hydrophobic residues to solvent and helices turning to loops cause dysfunction of ZNF141 protein. These significant changes in the stability and conformation of the mutant protein were validated using essential dynamics and cross-correlation maps, which revealed that upon point mutation, the overall motion of the proteins and the correlation between them were completely different, resulting in Postaxial polydactyly type A. Conclusions: This study provides molecular insights into the structural association of ZNF141 protein with PAP type A. Identification of active site residues and legends offers new therapeutic targets for ZNF141 protein. Further, it reiterates the functional importance of the last residue of a protein.