Ensembl variation resources

Four germline POLH variants, including two found in skin tumors, impair DNA polymerase η function and cellular tolerance to UV radiation and cisplatin

DNA polymerase (pol) η is vital for accurately replicating DNA opposite ultraviolet light (UV)-induced cyclobutane pyrimidine dimers and cisplatin-induced intrastrand purine crosslinks. While human POLH deficiency is linked to the disease xeroderma pigmentosum variant, the functional consequences of germline and somatic POLH variants remain largely unexplored. We characterized nine nonsynonymous POLH germline variants, five of which have also been found in various tumors. Enzyme activity was first assessed using recombinant pol η (residues 1-432) proteins. Variants F17S, C227Y, and R356X displayed substantially reduced or nearly abolished polymerase activity opposite cis-syn cyclobutane thymine dimer (CTD) compared to the wild-type. Cellular effects were then evaluated in POLH-deficient human embryonic kidney (HEK) 293 cells. Unlike cells transfected with wild-type POLH, cells transfected with F17S, R81C, C227Y, or R356X variants failed to rescue UV- and cisplatin-sensitivity. Interestingly, the R81C variant protein was undetectable in transfected cells. Further steady-state kinetic analysis revealed that the F17S, C227Y, and R356X variants had 3- to 5000-fold reductions in kcat/Km values for correct dATP insertion opposite CTD, while the R81C variant exhibited kinetics comparable to the wild-type enzyme. CRISPR/Cas9-mediated knock-in of the R81C variant in HEK 293T cells was associated with significantly impaired pol η protein expression and increased cisplatin sensitivity. Notably, R81C and R356X mutations have been reported in skin cancer samples. These findings suggest that R81C, F17S, C227Y, and R356X POLH variants-underexpressed or hypoactive-may be insufficient to protect cells from UV radiation and cisplatin, highlighting their potential implications for individual susceptibility to UV and cisplatin damage.

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.

A novel recurrent ARL3 variant c.209G > A p.(Gly70Glu) causes variable non-syndromic dominant retinal dystrophy with defective lipidated protein transport in human retinal stem cell models

Inherited retinal dystrophies (IRDs) are characterized by their high clinical and genetic heterogeneity. Despite significant advances in the identification of genes associated with IRDs, many individuals and families still have not received a definite molecular diagnosis. Here, we performed clinical examinations and conducted genetic testing in five families with IRD. Whole exome sequencing in the five index cases revealed a heterozygous missense variant, c.209G > A, p.(Gly70Glu) in the ARL3 gene (NM_004311.4). A de novo occurrence was demonstrated in one affected individual and autosomal dominant inheritance in nine affected individuals from four families. Their phenotypes displayed variable expressivity, and ranged from rod-cone to cone-rod dystrophy with photophobia. Human induced pluripotent stem cells (hiPSCs) were generated from dermal fibroblasts from the individual with the de novo ARL3 variant and were differentiated to retinal pigment epithelium cells (RPE) and retinal organoids. Immunofluorescence analyses in these models showed decreased INPP5E localization within the cilia of RPE and connecting cilia of retinal organoids, as well as reduced PDE6⍺ in the organoid outer segments, suggesting that the p.(Gly70Glu) variant causes IRD by defective lipidated protein transport in photoreceptors and/or RPE. This is the first study of ARL3 dysfunction in human retinal cells, highlighting its importance for retinal homeostasis, as well as a variability in the clinical presentation of ARL3-associated IRD.

Bi-allelic variants in three genes encoding distinct subunits of the vesicular AP-5 complex cause hereditary macular dystrophy

Inherited retinal diseases (IRDs) are a genetically heterogeneous group of Mendelian disorders that often lead to progressive vision loss and involve approximately 300 distinct genes. Although variants in these loci account for the majority of molecular diagnoses, other genes associated with IRD await molecular identification. In this study, we uncover bi-allelic assortments of 23 different (22 loss-of-function) variants in AP5Z1, AP5M1, and AP5B1 as independent causes of recessive IRD in members of 19 families from nine countries. Affected individuals, regardless of their genotypes, exhibit a specific form of macular degeneration, sometimes presenting in association with extraocular features. All three genes encode different subunits of the vesicular fifth adaptor protein (AP-5) complex, a component of the intracellular trafficking system involved in maintaining cellular homeostasis and ensuring the proper functioning of lysosomal pathways. The retinal pigment epithelium (RPE), a cellular monolayer located posteriorly to the neural retina, is characterized by intense lysosomal and phagocytic activity. Immunostaining of RPE cells revealed a punctate pattern of AP5Z1, AP5M1, and AP5B1 staining and co-localization with markers of late endosomes and the Golgi, suggesting a role of AP-5 in the normal physiology of this tissue. Overall, the identification of independently acting variants in three distinct proteins within the same macromolecular complex reveals AP-5 as having an important function in the preservation and maintenance of normal macular functions.

TRAPping the effects of tobacco smoking: the regulation and function of Acp5 expression in lung macrophages

Tartrate-resistant acid phosphatase [TRAP, gene acid phosphatase 5 (Acp5; gene name for TRAP)] is highly expressed in alveolar macrophages with proposed roles in lung inflammation and lung fibrosis development. We previously showed that its expression and activity are higher in lung macrophages of smokers and patients with chronic obstructive pulmonary disease (COPD), suggesting involvement in smoke-induced lung damage. In this study, we explored the function of TRAP and regulation of its different mRNA transcripts (Acp5 201-206) in lung tissue exposed to cigarette smoke to elucidate its function in alveolar macrophages. In mice exposed to cigarette smoke or air for 4-6 wk, higher Acp5 mRNA expression in lung tissue after smoking was mainly driven by transcript Acp5-202, which originates from macrophages. The expression of Acp5-202 correlated with transcription factors previously found to drive proliferation of macrophages. Treating fetal liver progenitor-derived alveolar-like macrophages [Max Planck Institute (MPI; macrophages derived from fetal liver progenitors) macrophages] with cigarette smoke extract resulted in more proliferation compared with nontreated cells. In contrast, Acp5-deficient MPI macrophages and MPI macrophages treated with a TRAP inhibitor proliferated significantly less than control macrophages. Mechanistically, this lack of proliferation after TRAP inhibition was associated with higher presence of phosphorylated Beta-catenin (β-catenin; a signaling protein) compared with nontreated controls. Phosphorylation of β-catenin is known to mark it for ubiquitination and degradation by the proteasome, preventing its activity in promoting cell proliferation. In conclusion, our findings provide strong evidence for TRAP stimulating alveolar macrophage proliferation by dephosphorylating β-catenin. By driving proliferation, TRAP likely helps sustain alveolar macrophage populations during smoke exposure, either compensating for their loss due to smoking or increasing their numbers to better manage smoke-induced damage.NEW & NOTEWORTHY This study has uncovered that the enzyme tartrate-resistant acid phosphatase (TRAP) is crucial for alveolar macrophage proliferation through a β-catenin-dependent pathway. Importantly, TRAP influences this important ability of alveolar macrophages through the Acp5-202 mRNA transcript. The increase in TRAP expression following smoke exposure suggests that it plays a key role in promoting cell renewal, potentially helping to mitigate smoke-induced lung damage.

Accessible, realistic genome simulation with selection using stdpopsim

Selection is a fundamental evolutionary force that shapes patterns of genetic variation across species. However, simulations incorporating realistic selection along heterogeneous genomes in complex demographic histories are challenging, limiting our ability to benchmark statistical methods aimed at detecting selection and to explore theoretical predictions. stdpopsim is a community-maintained simulation library that already provides an extensive catalog of species-specific population genetic models. Here we present a major extension to the stdpopsim framework that enables simulation of various modes of selection, including background selection, selective sweeps, and arbitrary distributions of fitness effects (DFE) acting on annotated subsets of the genome (for instance, exons). This extension maintains stdpopsim's core principles of reproducibility and accessibility while adding support for species-specific genomic annotations and published DFE estimates. We demonstrate the utility of this framework by benchmarking methods for demographic inference, DFE estimation, and selective sweep detection across several species and scenarios. Our results demonstrate the robustness of demographic inference methods to selection on linked sites, reveal the sensitivity of DFE-inference methods to model assumptions, and show how genomic features, like recombination rate and functional sequence density, influence power to detect selective sweeps. This extension to stdpopsim provides a powerful new resource for the population genetics community to explore the interplay between selection and other evolutionary forces in a reproducible, low-barrier framework.

Integrative In Silico Analysis to Identify Functional and Structural Impacts of nsSNPs on Programmed Cell Death Protein 1 (PD-1) Protein and UTRs: Potential Biomarkers for Cancer Susceptibility

Background: Programmed cell death protein 1 (PD-1), encoded by the PDCD1 gene, is critical in immune checkpoint regulation and cancer immune evasion. Variants in PDCD1 may alter its function, impacting cancer susceptibility and disease progression. Objectives: This study evaluates the structural, functional, and regulatory impacts of non-synonymous single-nucleotide polymorphisms (nsSNPs) in the PDCD1 gene, focusing on their pathogenic and oncogenic roles. Methods: Computational tools, including PredictSNP1.0, I-Mutant2.0, MUpro, HOPE, MutPred2, Cscape, Cscape-Somatic, GEPIA2, cBioPortal, and STRING, were used to analyze 695 nsSNPs in the PD1 protein. The analysis covered structural impacts, stability changes, regulatory effects, and oncogenic potential, focusing on conserved domains and protein-ligand interactions. Results: The analysis identified 84 deleterious variants, with 45 mapped to conserved regions like the Ig V-set domain essential for ligand-binding interactions. Stability analyses identified 78 destabilizing variants with significant protein instability (ΔΔG values). Ten nsSNPs were identified as potential cancer drivers. Expression profiling showed differential PDCD1 expression in tumor versus normal tissues, correlating with improved survival in skin melanoma but limited value in ovarian cancer. Regulatory SNPs disrupted miRNA-binding sites and transcriptional regulation, affecting PDCD1 expression. STRING analysis revealed key PD-1 protein partners within immune pathways, including PD-L1 and PD-L2. Conclusions: This study highlights the significance of PDCD1 nsSNPs as potential biomarkers for cancer susceptibility, advancing the understanding of PD-1 regulation. Experimental validation and multi-omics integration are crucial to refine these findings and enhance theraputic strategies.

An integrated DNA interactome and transcriptome profiling reveals a PU.1/enhancer RNA-mediated Feed-forward Regulatory Loop Regulating monocyte/macrophage development and innate immune functions

High expression of the myeloid master ETS transcription factor PU.1 drives the development of monocyte/macrophage (Mono/MΦ), a crucial cellular component of the innate immune system. Disruptions in normal expression patterns of PU.1 are linked to a variety myeloid malignancy and immune diseases. It is evidenced that PU.1 binds to and modulates enhancers of several myeloid genes. While noncoding RNAs transcribed from noncoding genes at the enhancers are increasingly reported to be involved in enhancer regulation, the crosstalk between PU.1 and noncoding RNAs in enhancer-mediated myeloid gene regulation in Mono/MΦ differentiation and immune response has not been systematically investigated. In this study, we interrogated the PU.1-mediated transcriptome and cistrome with our comprehensive collection of putative and verified enhancers. Among a repertoire of noncoding genes present at PU.1-bound enhancers, we discovered that PU.1 acts as a potent transcription factor inducer of the noncoding RNA LOUP , which we previously identified as an RNA inducer of PU.1. The genomic region within the LOUP locus occupied by PU.1 is characterized by the epigenetic features of a myeloid-specific super-enhancer. Targeted disruption of the PU.1-binding motifs resulted in the downregulation of LOUP promoter activity. Depletion of LOUP reduced the expression of Mono/MΦ cell markers as well as the transcriptional program associated with Mono/MΦ differentiation Mono/MΦ innate defense mechanisms, including phagocytosis, antimicrobial activity, and chemoattractant cytokine production. LOUP induces Mono/MΦ phagocytic activities. Collectively, our findings indicate that PU.1 and enhancer RNA LOUP are biomolecular components of an unidentified feed-forward loop that promotes their mutual expression, contributing to Mono/MΦ differentiation and innate immune functions. The identification of the PU.1/ LOUP regulatory circuit provides valuable insights into the mechanisms underlying cell-type and gene-specific enhancer activity and Mono/MΦ biology, as well as significant implications for advancing our understanding of immune diseases and myeloid malignancies.

Sequence-based GWAS in 180,000 German Holstein cattle reveals new candidate variants for milk production traits

BACKGROUND

Milk production traits are complex and influenced by many genetic and environmental factors. Although extensive research has been performed for these traits, with many associations unveiled thus far, due to their crucial economic importance, complex genetic architecture, and the fact that causal variants in cattle are still scarce, there is a need for a better understanding of their genetic background. In this study, we aimed to identify new candidate loci associated with milk production traits in German Holstein cattle, the most important dairy breed in Germany and worldwide. For that purpose, 180,217 cattle were imputed to the sequence level and large-scale genome-wide association study (GWAS) followed by fine-mapping and evolutionary and functional annotation were carried out to identify and prioritize new association signals.

RESULTS

Using the imputed sequence data of a large cattle dataset, we identified 50,876 significant variants, confirming many known and identifying previously unreported candidate variants for milk (MY), fat (FY), and protein yield (PY). Genome-wide significant signals were fine-mapped with the Bayesian approach that determines the credible variant sets and generates the probability of causality for each signal. The variants with the highest probabilities of being causal were further classified using external information about the function and evolution, making the prioritization for subsequent validation experiments easier. The top potential causal variants determined with fine-mapping explained a large percentage of genetic variance compared to random ones; 178 variants explained 11.5%, 104 explained 7.7%, and 68 variants explained 3.9% of the variance for MY, FY, and PY, respectively, demonstrating the potential for causality.

CONCLUSIONS

Our findings proved the power of large samples and sequence-based GWAS in detecting new association signals. In order to fully exploit the power of GWAS, one should aim at very large samples combined with whole-genome sequence data. These can also come with both computational and time burdens, as presented in our study. Although milk production traits in cattle are comprehensively investigated, the genetic background of these traits is still not fully understood, with the potential for many new associations to be revealed, as shown. With constantly growing sample sizes, we expect more insights into the genetic architecture of milk production traits in the future.

Unveiling the Digital Evolution of Molecular Tumor Boards

Molecular tumor boards (MTB) are interdisciplinary conferences involving various experts discussing patients with advanced tumors, to derive individualized treatment suggestions based on molecular variants. These discussions involve using heterogeneous internal data, such as patient clinical data, but also external resources such as knowledge databases for annotations and search for relevant clinical studies. This imposes a certain level of complexity that requires huge effort to homogenize the data and use it in a speedy manner to reach the needed treatment. For this purpose, most institutions involving an MTB are heading toward automation and digitalization of the process, hence reducing manual work requiring human intervention and subsequently time in deriving personalized treatment suggestions. The tools are also used to better visualize the patient's data, which allows a refined overview for the board members. In this paper, we present the results of our thorough literature research about MTBs, their process, the most common knowledge bases, and tools used to support this decision-making process.

Elucidating the role of MLL1 nsSNPs: Structural and functional alterations and their contribution to leukemia development

(1) BACKGROUND

The Mixed lineage leukemia 1 (MLL1) gene, located on chromosome 11q23, plays a pivotal role in histone lysine-specific methylation and is consistently associated with various types of leukemia. Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) have been tied to numerous diseases, including cancers, and have become valuable cancer biomarkers. There's a notable gap in studies probing the influence of SNPs on MLL1 protein structure, function, and subsequent modifications.

(2) METHODS

We utilized an array of bioinformatics tools, including PredictSNP, InterPro, ConSurf, I-Mutant2.0, MUpro, Musitedeep, Project HOPE, RegulomeDB, Mutpred2, and both CScape and CScape Somatic, to meticulously analyze the consequences of nsSNPs in the MLL1 gene.

(3) RESULTS

Out of 2,097 nsSNPs analyzed, 62 were determined to be significantly pathogenic by the PredictSNP tool, with ten crucial MLL1 functional domains identified using InterPro. Additionally, 50 of these nsSNPs had high conservation scores, hinting at potential effects on protein structure and function, while 32 were found to undermine MLL1 protein stability. Notably, four nsSNPs were deemed oncogenic, with two identified as cancer drivers. The nsSNP, D2724G, between the MLL1 protein's FY-rich domains, could disrupt proteolytic cleavage, altering gene expression patterns and potentially promoting cancer.

(4) CONCLUSIONS

Our research provides a comprehensive assessment of nsSNPs' impact in the MLL1 protein structure and function and consequently on leukemia development, suggesting potential avenues for personalized treatment, early detection, improved prognosis, and a deeper understanding of hematological malignancy genesis.

Effect of sequencing platforms on the sensitivity of chemical mutation detection using Hawk-Seq™

BACKGROUND

Error-corrected next-generation sequencing (ecNGS) technologies have enabled the direct evaluation of genome-wide mutations after exposure to mutagens. Previously, we reported an ecNGS methodology, Hawk-Seq™, and demonstrated its utility in evaluating mutagenicity. The evaluation of technical transferability is essential to further evaluate the reliability of ecNGS-based assays. However, cutting-edge sequencing platforms are continually evolving, which can affect the sensitivity of ecNGS. Therefore, the effect of differences in sequencing instruments on mutation data quality should be evaluated.

RESULTS

We assessed the performance of four sequencing platforms (HiSeq2500, NovaSeq6000, NextSeq2000, and DNBSEQ-G400) with the Hawk-Seq™ protocol for mutagenicity evaluation using DNA samples from mouse bone marrow exposed to benzo[a]pyrene (BP). The overall mutation (OM) frequencies per 106 bp in vehicle-treated samples were 0.22, 0.36, 0.46, and 0.26 for HiSeq2500, NovaSeq6000, NextSeq2000, and DNBSEQ-G400, respectively. The OM frequency of NextSeq2000 was significantly higher than that of HiSeq2500, suggesting the difference to be based on the platform. The relatively higher value in NextSeq2000 was a consequence of the G:C to C:G mutations in NextSeq2000 data (0.67 per 106 G:C bp), which was higher than the mean of the four platforms by a ca. of 0.25 per 106 G:C bp. A clear dose-dependent increase in G:C to T:A mutation frequencies was observed in all four sequencing platforms after BP exposure. The cosine similarity values of the 96-dimensional trinucleotide mutation patterns between HiSeq and the three other platforms were 0.93, 0.95, and 0.92 for NovaSeq, NextSeq, and DNBSeq, respectively. These results suggest that all platforms can provide equivalent data that reflect the characteristics of the mutagens.

CONCLUSIONS

All platforms sensitively detected mutagen-induced mutations using the Hawk-Seq™ analysis. The substitution types and frequencies of the background errors differed depending on the platform. The effects of sequencing platforms on mutagenicity evaluation should be assessed before experimentation.

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.