Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration

μ-Opioid receptor transcriptional variants in the murine forebrain and spinal cord

Oprm1, the gene encoding the μ-opioid receptor, has multiple reported transcripts, with a variable 3' region and many alternative sequences encoding the C-terminus of the protein. The functional implications of this variability remain mostly unexplored, though a recurring notion is that it could be exploited by developing selective ligands with improved clinical profiles. Here, we comprehensively examined Oprm1 transcriptional variants in the murine central nervous system, using long-read RNAseq as well as spatial and single-cell transcriptomics. The results were validated with RNAscope in situ hybridization. We found a mismatch between transcripts annotated in the mouse genome (GRCm38/mm10) and the RNA-seq results. Sequencing data indicated that the primary Oprm1 transcript has a 3' terminus located on chr10:6,860,027, which is ∼ 9.5 kilobases downstream of the longest annotated exon 4 end. Long-read sequencing confirmed that the final Oprm1 exon included a 10.2 kilobase long 3' untranslated region, and the presence of the long variant was unambiguously confirmed using RNAscope in situ hybridization in the thalamus, striatum, cortex and spinal cord. Conversely, expression of the Oprm1 reference transcript or alternative transcripts of the Oprm1 gene was absent or close to the detection limit. Thus, the primary transcript of the Oprm1 mouse gene is a variant with a long 3' untranslated region, which is homologous to the human OPRM1 primary transcript and encodes the same conserved C-terminal amino acid sequence.

Unraveling novel mRNA transcripts of the human DNA N-glycosylase 1 (NTHL1) gene with the implementation of an innovative targeted DNA-seq assay

The human NTHL1 gene encodes a DNA glycosylase that plays a key role in the base excision repair (BER) pathway, repairing oxidative DNA damage and maintaining genome integrity. The physiological activity of NTHL1 is crucial in preventing genetic alterations that can lead to cancer. In this study, we employed an innovative targeted DNA sequencing (DNA-seq) methodology to explore the transcriptional landscape of the NTHL1 gene, revealing previously uncharacterized alternative splicing events and novel exons. Our designed approach provided significantly improved sequencing depth and coverage, enabling the identification of novel NTHL1 mRNA transcripts. Bioinformatics analysis confirmed all annotated splice junctions of the main NTHL1 transcripts (v.1 - v.3) and revealed novel mRNA transcripts (NTHL1 v.4 - v.9) derived from splicing events between annotated exons as well as mRNAs containing previously uncharacterized exons (NTHL1 v.10 - v.14). Quantitative PCR analysis highlighted a diverse expression pattern of these novel transcripts across different human cell lines, suggesting cell-specific roles and regulatory mechanisms. Notably, NTHL1 v.5 was overexpressed in luminal A breast cancer cells (MCF-7), while v.13 was prominent in triple negative (BT-20), HER2 + breast cancer (SK-BR-3), prostate, colorectal cancer cells and HEK-293 cells. Our findings suggest that specific novel NTHL1 transcripts may encode protein isoforms with distinct structural features, as indicated by ribosome profiling datasets, while others containing premature termination codons could function as long non-coding RNAs. These insights enhance our understanding of NTHL1 regulatory role and its potential as a biomarker and therapeutic target in human malignancies. This study underscores the importance of exploring the transcriptional diversity of NTHL1 to fully elucidate its role in cancer pathobiology.

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Copine C plays a role in adhesion and streaming in Dictyostelium

Copines are a family of calcium-dependent membrane-binding proteins. To study these proteins, anull mutant for cpnC was created in Dictyostelium, which has six copines genes (cpnA-cpnF). During development, cpnC- cells were able to aggregate, but did not form streams. Once aggregated into mounds, they formed large ring structures. cpnC- cells were less adherent to plastic substrates, but more adherent to other cells. These phenotypes correlated with changes in adhesion protein expression with decreased expression of SibA and increased expression of CsaA in developing cpnC- cells. We also measured the expression of RegA, a cAMP phosphodiesterase, and found that cpnC- cells have reduced RegA expression. The reduced RegA expression in cpnC- cells is most likely responsible for the observed phenotypes.

Integrated multi-omics analysis reveals molecular changes associated with chronic lipid accumulation following contusive spinal cord injury

Functional and pathological recovery after spinal cord injury (SCI) is often incomplete due to the limited regenerative capacity of the central nervous system (CNS), which is further impaired by several mechanisms that sustain tissue damage. Among these, the chronic activation of immune cells can cause a persistent state of local CNS inflammation and damage. However, the mechanisms that sustain this persistent maladaptive immune response in SCI have not been fully clarified yet. In this study, we integrated histological analyses with proteomic, lipidomic, transcriptomic, and epitranscriptomic approaches to study the pathological and molecular alterations that develop in a mouse model of cervical spinal cord hemicontusion. We found significant pathological alterations of the lesion rim with myelin damage and axonal loss that persisted throughout the late chronic phase of SCI. This was coupled by a progressive lipid accumulation in myeloid cells, including resident microglia and infiltrating monocyte-derived macrophages. At tissue level, we found significant changes of proteins indicative of glycolytic, tricarboxylic acid cycle (TCA), and fatty acid metabolic pathways with an accumulation of triacylglycerides with C16:0 fatty acyl chains in chronic SCI. Following transcriptomic, proteomic, and epitranscriptomic studies identified an increase of cholesterol and m6A methylation in lipid-droplet-accumulating myeloid cells as a core feature of chronic SCI. By characterizing the multiple metabolic pathways altered in SCI, our work highlights a key role of lipid metabolism in the chronic response of the immune and central nervous system to damage.

Analysis of the finasteride treatment and its withdrawal in the rat hypothalamus and hippocampus at whole-transcriptome level

PURPOSE

As reported in patients treated for androgenetic alopecia with finasteride (i.e., a blocker of the enzyme 5 alpha-reductase) and in an animal model, side effects affecting sexual, psychiatric, neurological, and physical domains, may occur during the treatment and persist with drug suspension. The etiopathogenesis of these side effects has been poorly explored. Therefore, we performed a genome-wide analysis of finasteride effects in the brain of adult male rat.

METHODS

Animals were treated (i.e., for 20 days) with finasteride (1mg/rat/day). 24 h after the last treatment and 1 month after drug suspension, RNA sequencing analysis was performed in hypothalamus and hippocampus. Data were analyzed by differential expression analysis and Gene-Set Enrichment Analyses (GSEA).

RESULTS

Data obtained after finasteride treatment showed that 186 genes (i.e., 171 up- and 15 downregulated) and 19 (i.e., 17 up- and 2 downregulated) were differentially expressed in the hypothalamus and hippocampus, respectively. Differential expression analysis at the drug withdrawal failed to identify dysregulated genes. Several gene-sets were enriched in these brain areas at both time points.

CONCLUSION

Some of the genes reported to be differentially expressed (i.e., TTR, DIO2, CLDN1, CLDN2, SLC4A5, KCNE2, CROT, HCRT, MARCKSL1, VGF, IRF2BPL) and GSEA, suggest a potential link with specific side effects previously observed in patients and in the animal model, such as depression, anxiety, disturbance in memory and attention, and sleep disturbance. These data may provide an important background for future experiments aimed at confirming the pathological role of these genes.

Next-Generation Sequencing and Emerging Technologies

Genetic sequencing technologies are evolving at a rapid pace with major implications for research and clinical practice. In this review, the authors provide an updated overview of next-generation sequencing (NGS) and emerging methodologies. NGS has tremendously improved sequencing output while being more time and cost-efficient in comparison to Sanger sequencing. The authors describe short-read sequencing approaches, such as sequencing by synthesis, ion semiconductor sequencing, and nanoball sequencing. Third-generation long-read sequencing now promises to overcome many of the limitations of short-read sequencing, such as the ability to reliably resolve repeat sequences and large genomic rearrangements. By combining complementary methods with massively parallel DNA sequencing, a greater insight into the biological context of disease mechanisms is now possible. Emerging methodologies, such as advances in nanopore technology, in situ nucleic acid sequencing, and microscopy-based sequencing, will continue the rapid evolution of this area. These new technologies hold many potential applications for hematological disorders, with the promise of precision and personalized medical care in the future.

Helicases DDX5 and DDX17 promote heterogeneity in HBV transcription termination in infected human hepatocytes

BACKGROUND & AIMS

Transcription termination fine-tunes gene expression and contributes to the specification of RNA function in eukaryotic cells. Transcription termination of HBV is subject to the recognition of the canonical polyadenylation signal (cPAS) common to all viral transcripts. However, the regulation of this cPAS and its impact on viral gene expression and replication is currently unknown.

METHODS

To unravel the regulation of HBV transcript termination, we implemented a 3' RACE (rapid amplification of cDNA ends)-PCR assay coupled to single molecule sequencing both in in vitro-infected hepatocytes and in chronically infected patients.

RESULTS

The detection of a previously unidentified transcriptional readthrough indicated that the cPAS was not systematically recognized during HBV replication in vitro and in vivo. Gene expression downregulation experiments demonstrated a role for the RNA helicases DDX5 and DDX17 in promoting viral transcriptional readthrough, which was, in turn, associated with HBV RNA destabilization and decreased HBx protein expression. RNA and chromatin immunoprecipitation, together with mutation of the cPAS sequence, suggested a direct role of DDX5 and DDX17 in functionally linking cPAS recognition to transcriptional readthrough, HBV RNA stability and replication.

CONCLUSIONS

Our findings identify DDX5 and DDX17 as crucial determinants of HBV transcriptional fidelity and as host restriction factors for HBV replication.

IMPACT AND IMPLICATIONS

HBV covalently closed circular (ccc)DNA degradation or functional inactivation remains the holy grail for the achievement of HBV cure. Transcriptional fidelity is a cornerstone in the regulation of gene expression. Here, we demonstrate that two helicases, DDX5 and DDX17, inhibit recognition of the HBV polyadenylation signal and thereby transcriptional termination, thus decreasing HBV RNA stability and acting as restriction factors for efficient cccDNA transcription and viral replication. The observation that DDX5 and DDX17 are downregulated in patients chronically infected with HBV suggests a role for these helicases in HBV persistence in vivo. These results open new perspectives for researchers aiming at identifying new targets to neutralise cccDNA transcription.

Fatal infection caused by a genetically distinct elephant endotheliotropic herpesvirus type 5 in a captive Asian elephant in Germany

BACKGROUND

Elephant endotheliotropic herpesvirus (EEHV) infection is the most common cause for lethal hemorrhagic disease in captive juvenile Asian elephants (Elephas maximus). Although EEHV1 is known as the most likely cause of fatal haemorrhagic disease in Asian elephants, EEHV5 was lately involved in lethal cases of haemorrhagic disease in captive elephants.

CASE PRESENTATION

Here we report the first death of a four-year old Asian elephant diagnosed with EEHV5 in Germany. Molecular diagnosis yielded detection of EEHV5 DNA in all tested tissues. Histopathological examination revealed typical features of hemorrhagic disease in all examined organs. EEHV5 was sequenced from total DNA isolated from heart tissue by Illumina and Nanopore sequencing. Sequencing data showed 3,881 variants, distributed across the entire genome, compared to the published EEHV5 sequence.

CONCLUSIONS

We have detected EEHV5 in a fatal disease case of a male Asian elephant. Whole genome sequencing revealed substantial differences of our DNA isolate compared to available EEHV5 sequences. This report of fatal haemorrhagic disease associated with EEHV5 infection should raise awareness for EEHV5 as an important elephant pathogen. Genome sequencing and downstream SNPs analysis will further encourage future research to understand genetic diversity, pathogenesis and virulence of EEHVs with respect to developing new diagnostic methods, prophylactic strategies, and implementation of surveillance and control measures.

PPARG is dispensable for bovine embryo development up to tubular stages†

Following blastocyst hatching, ungulate embryos undergo a prolonged preimplantation period termed conceptus elongation. Conceptus elongation constitutes a highly susceptible period for embryonic loss, and the embryonic requirements during this process are largely unknown, but multiple lipid compounds have been identified in the fluid nourishing the elongating conceptuses. Peroxisome proliferator-activated receptors mediate the signaling actions of prostaglandins and other lipids, and, between them, PPARG has been pointed out to play a relevant role in conceptus elongation by a functional study that depleted PPARG in both uterus and conceptus. The objective of this study has been to determine if embryonic PPARG is required for bovine embryo development. To that aim, we have generated bovine PPARG knock-out embryos in vitro using two independent gene ablation strategies and assessed their developmental ability. In vitro development to Day 8 blastocyst was unaffected by PPARG ablation, as total, inner cell mass, and trophectoderm cell numbers were similar between wild-type and knock-out D8 embryos. In vitro post-hatching development to D12 was also comparable between different genotypes, as embryo diameter, epiblast cell number, embryonic disk formation, and hypoblast migration rates were unaffected by the ablation. The development of tubular stages equivalent to E14 was assessed in vivo, following a heterologous embryo transfer experiment, observing that the development of extra-embryonic membranes and of the embryonic disk was not altered by PPARG ablation. In conclusion, PPARG ablation did not impaired bovine embryo development up to tubular stages.

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion

BACKGROUND

Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions.

RESULTS

Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP".

CONCLUSION

The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

X-Linked Bilateral Polymicrogyria With Epilepsy and Intellectual Disability Associated With a Novel KIF4A Variant

We studied three brothers and a maternal half-brother featuring global developmental delay, mild to moderate intellectual disability, epilepsy, microcephaly, and strabismus. All had bilateral perisylvian and perirolandic polymicrogyria, while some also had malformations of the hippocampus (malrotation and dysplasia), cerebellum (heterotopias and asymmetric aplasia), corpus callosum dysgenesis, and brainstem asymmetric dysplasia. Exome sequencing showed that all four patients had a novel variant (c.1597C>T:p.Leu533Phe) on the KIF4A gene on chromosome X. We discuss how this variant is possibly pathogenic and could explain the reported phenotype.

Exploring the transcription start sites and other genomic features facilitates the accurate identification and annotation of small RNAs across multiple stress conditions in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) is a pathogen that is known for its ability to persist in harsh environments and cause chronic infections. Understanding the regulatory networks of MTB is crucial for developing effective treatments. Small regulatory RNAs (sRNAs) play important roles in gene expression regulation in all kingdoms of life, and their classification based solely on genomic location can be imprecise due to the computational-based prediction of protein-coding genes in bacteria, which often neglects segments of mRNA such as 5'UTRs, 3'UTRs, and intercistronic regions of operons. To address this issue, our study simultaneously discovered genomic features such as TSSs, UTRs, and operons together with sRNAs in the M. tuberculosis H37Rv strain (ATCC 27294) across multiple stress conditions. Our analysis identified 1,376 sRNA candidates and 8,173 TSSs in MTB, providing valuable insights into its complex regulatory landscape. TSS mapping enabled us to classify these sRNAs into more specific categories, including promoter-associated sRNAs, 5'UTR-derived sRNAs, 3'UTR-derived sRNAs, true intergenic sRNAs, and antisense sRNAs. Three of these sRNA candidates were experimentally validated using 3'-RACE-PCR: predictedRNA_0240, predictedRNA_0325, and predictedRNA_0578. Future characterization and validation are necessary to fully elucidate the functions and roles of these sRNAs in MTB. Our study is the first to simultaneously unravel TSSs and sRNAs in MTB and demonstrate that the identification of other genomic features, such as TSSs, UTRs, and operons, allows for more accurate and specific classification of sRNAs.

Resolving the metabolism of monolignols and other lignin-related aromatic compounds in Xanthomonas citri

Lignin, a major plant cell wall component, has an important role in plant-defense mechanisms against pathogens and is a promising renewable carbon source to produce bio-based chemicals. However, our understanding of microbial metabolism is incomplete regarding certain lignin-related compounds like p-coumaryl and sinapyl alcohols. Here, we reveal peripheral pathways for the catabolism of the three main lignin precursors (p-coumaryl, coniferyl, and sinapyl alcohols) in the plant pathogen Xanthomonas citri. Our study demonstrates all the necessary enzymatic steps for funneling these monolignols into the tricarboxylic acid cycle, concurrently uncovering aryl aldehyde reductases that likely protect the pathogen from aldehydes toxicity. It also shows that lignin-related aromatic compounds activate transcriptional responses related to chemotaxis and flagellar-dependent motility, which might play an important role during plant infection. Together our findings provide foundational knowledge to support biotechnological advances for both plant diseases treatments and conversion of lignin-derived compounds into bio-based chemicals.

Genomic dynamics of the Lower Yellow River Valley since the Early Neolithic

The Yellow River Delta played a vital role in the development of the Neolithic civilization of China. However, the population history of this region from the Neolithic transitions to the present remains poorly understood due to the lack of ancient human genomes. This especially holds for key Neolithic transitions and tumultuous turnovers of dynastic history. Here, we report genome-wide data from 69 individuals dating to 5,410-1,345 years before present (BP) at 0.008 to 2.49× coverages, along with 325 present-day individuals collected from 16 cities across Shandong. During the Middle to Late Dawenkou period, we observed a significant influx of ancestry from Neolithic Yellow River farmers in central China and some southern Chinese ancestry that mixed with local hunter-gatherers in Shandong. The genetic heritage of the Shandong Longshan people was found to be most closely linked to the Dawenkou culture. During the Shang to Zhou Dynasties, there was evidence of genetic admixture of local Longshan populations with migrants from the Central Plain. After the Qin to Han Dynasties, the genetic composition of the region began to resemble that of modern Shandong populations. Our genetic findings suggest that the middle Yellow River Basin farmers played a role in shaping the genetic affinity of neighboring populations in northern China during the Middle to Late Neolithic period. Additionally, our findings indicate that the genetic diversity in the Shandong region during the Zhou Dynasty may be linked with their complex ethnicities.

H3T11 phosphorylation by CKII is required for heterochromatin formation in Neurospora

Heterochromatin is a key feature of eukaryotic genomes and is crucial for maintaining genomic stability. In fission yeast, heterochromatin nucleation is mainly mediated by DNA-binding proteins or the RNA interference (RNAi) pathway. In the filamentous fungus Neurospora crassa, however, the mechanism that causes the initiation of heterochromatin at the relics of repeat-induced point mutation is unknown and independent of the classical RNAi pathway. Here, we show that casein kinase II (CKII) and its kinase activity are required for heterochromatin formation at the well-defined 5-kb heterochromatin of the 5H-cat-3 region and transcriptional repression of its adjacent cat-3 gene. Similarly, mutation of the histone H3 phosphorylation site T11 also impairs heterochromatin formation at the same locus. The catalytic subunit CKA colocalizes with H3T11 phosphorylation (H3pT11) within the 5H-cat-3 domain and the deletion of cka results in a significant decrease in H3T11 phosphorylation. Furthermore, the loss of kinase activity of CKII results in a significant reduction of H3pT11, H3K9me3 (histone H3 lysine 9 trimethylation) and DNA methylation levels, suggesting that CKII regulates heterochromatin formation by promoting H3T11 phosphorylation. Together, our results establish that histone H3 phosphorylation by CKII is a critical event required for heterochromatin formation.

Characterization of sequence variations in the extended flanking regions using massively parallel sequencing in 21 A-STRs and 21 Y-STRs

In forensic genetics, utilizing massively parallel sequencing (MPS) to analyze short tandem repeats (STRs) has demonstrated several advantages compared to conventional capillary electrophoresis (CE). Due to the current technical limitations, although flanking region polymorphisms had been mentioned in several previous studies, most studies focused on the core repeat regions of STRs or the variations in the adjacent flanking regions. In this study, we developed an MPS system consisting of two sets of multiplex PCR systems to detect not only the STR core repeat regions but also to observe variants located at relatively distant positions in the flanking regions. The system contained 42 commonly used forensic STRs, including 21 autosomal STRs (A-STRs) and 21 Y-chromosomal STRs (Y-STRs), and a total of 350 male individuals from a Chinese Han population were genotyped. The length and sequence variants per locus were tallied and categorized based on length (length-based, LB), sequence without flanking region (core repeat regions sequence-based, RSB), and sequence with flanking region (core repeat and flanking regions sequence-based, FSB), respectively. Allele frequencies, Y-haplotype frequencies, and forensic parameters were calculated based on LB, RSB, and FSB, respectively, to evaluate the improvement in discrimination power, heterozygosity, and effectiveness of forensic systems. The results suggested the sequence variations have more influence on A-STRs and could improve the identification ability of MPS-STR genotyping. Concordance between MPS and CE methods was confirmed by using commercial CE-based STR kits. The impact of flanking region variations on STR genotype analysis and potential factors contributing to discordances were discussed. A total of 58 variations in the flanking regions (53 SNPs/SNVs and 5 InDels) were observed and most variations (48/58) were distributed in A-STRs. In summary, this study delved deeper into the genetic information of forensic commonly used STR and advanced the application of massively parallel sequencing in forensic genetics.

Characterisation of the erythrocyte invasion phenotype of FCB-2: A South American P. falciparum reference strain

The extent of parasite adaptive capability involved in erythrocyte invasion represents a significant challenge for the development of a Plasmodium falciparum vaccine. The parasite's geographical and populational origin may influence such adaptive behaviour; in vitro culture-adapted parasite strains are typically used for such studies. Previous studies have reported invasion phenotypes in strains from Africa and Asia and, to a lesser extent, from Latin America. This study was aimed at expanding the pool of characterised parasite strains from Latin America by describing the invasion phenotype of the P. falciparum Colombia Bogotá 2 (FCB2) strain. The FCB2 genome was sequenced and erythrocyte invasion ligand sequences were analysed and compared to other previously reported ones. RT-PCR was used for assessing Pfeba family erythrocyte invasion ligands and reticulocyte binding homologue (Pfrh) gene transcription. A flow cytometry-based erythrocyte invasion assay (using enzymatically-treated erythrocytes) was used for determining the FCB2 strain's invasion phenotype. The P. falciparum FCB2 genome sequence was analysed, bearing in mind that prolonged in vitro parasite culture may affect its genome sequence and, in some cases, lead to the deletion of certain genes; it was demonstrated that all erythrocyte invasion ligand gene sequences studied here were preserved. Comparative analysis showed that the target genome sequences were conserved whereas transcriptional analysis highlighted Pfebas and Pfrhs gene expression. Erythrocyte invasion analysis demonstrated that the FCB2 strain has a sialic acid-resistant invasion phenotype.

Genome-Wide Association Studies and QTL Mapping Reveal a New Locus Associated with Resistance to Bacterial Pustule Caused by Xanthomonas citri pv. glycines in Soybean

Bacterial pustule (BP), caused by Xanthomonas citri pv. glycines, is an important disease that, under favorable conditions, can drastically affect soybean production. We performed a genome-wide association study (GWAS) with a panel containing Brazilian and American cultivars, which were screened qualitatively and quantitatively against two Brazilian X. citri isolates (IBS 333 and IBS 327). The panel was genotyped using a genotyping by sequencing (GBS) approach, and we identified two main new regions in soybeans associated with X. citri resistance on chromosomes 6 (IBS 333) and 18 (IBS 327), different from the traditional rxp gene located on chromosome 17. The region on chromosome 6 was also detected by QTL mapping using a biparental cross between Williams 82 (R) and PI 416937 (S), showing that Williams 82 has another recessive resistance gene besides rxp, which was also detected in nine BP-resistant ancestors of the Brazilian cultivars (including CNS, S-100), based on haplotype analysis. Furthermore, we identified additional SNPs in strong LD (0.8) with peak SNPs by exploring variation available in WGS (whole genome sequencing) data among 31 soybean accessions. In these regions in strong LD, two candidate resistance genes were identified (Glyma.06g311000 and Glyma.18g025100) for chromosomes 6 and 18, respectively. Therefore, our results allowed the identification of new chromosomal regions in soybeans associated with BP disease, which could be useful for marker-assisted selection and will enable a reduction in time and cost for the development of resistant cultivars.

Mono-methylation of lysine 27 at histone 3 confers lifelong susceptibility to stress

Histone post-translational modifications are critical for mediating persistent alterations in gene expression. By combining unbiased proteomics profiling and genome-wide approaches, we uncovered a role for mono-methylation of lysine 27 at histone H3 (H3K27me1) in the enduring effects of stress. Specifically, mice susceptible to early life stress (ELS) or chronic social defeat stress (CSDS) displayed increased H3K27me1 enrichment in the nucleus accumbens (NAc), a key brain-reward region. Stress-induced H3K27me1 accumulation occurred at genes that control neuronal excitability and was mediated by the VEFS domain of SUZ12, a core subunit of the polycomb repressive complex-2, which controls H3K27 methylation patterns. Viral VEFS expression changed the transcriptional profile of the NAc, led to social, emotional, and cognitive abnormalities, and altered excitability and synaptic transmission of NAc D1-medium spiny neurons. Together, we describe a novel function of H3K27me1 in the brain and demonstrate its role as a "chromatin scar" that mediates lifelong stress susceptibility.

Characterization of mycobacteriophage Adephagia cytotoxic proteins

Mycobacterium phage Adephagia is a cluster K phage that infects Mycobacterium smegmatis and some strains of Mycobacterium pathogens. Adephagia has a siphoviral virion morphology and is temperate. Its genome is 59,646 bp long and codes for one tRNA gene and 94 predicted protein-coding genes; most genes not associated with virion structure and assembly are functionally ill-defined. Here, we determined the Adephagia gene expression patterns in lytic and lysogenic growth and used structural predictions to assign additional gene functions. We characterized 66 nonstructural genes for their toxic phenotypes when expressed in M. smegmatis, and we show that 25 of these (38%) are either toxic or strongly inhibit growth, resulting in either reduced viability or small colony sizes. Some of these genes are predicted to be involved in DNA metabolism or regulation, but others are of unknown function. We also characterize the HicAB-like toxin-antitoxin (TA) system encoded by Adephagia (gp91 and gp90, respectively) and show that the gp90 antitoxin is lysogenically expressed, abrogates gp91 toxicity, and is required for normal lytic and lysogenic growth.

Discovery of new acetamide derivatives of 5-indole-1,3,4-oxadiazol-2-thiol as inhibitors of HIV-1 Tat-mediated viral transcription

Human immunodeficiency virus-1 (HIV-1) encodes a transcriptional factor called Tat, which is critical for viral transcription. Tat-mediated transcription is highly ordered apart from the cellular manner; therefore, it is considered a target for developing anti-HIV-1 drugs. However, drugs targeting Tat-mediated viral transcription are not yet available. Our high-throughput screen of a compound library employing a dual-reporter assay identified a 1,3,4-oxadiazole scaffold against Tat and HIV-1 infection. Furthermore, a serial structure-activity relation (SAR) study performed with biological assays found 1,3,4-oxadiazole derivatives (9 and 13) containing indole and acetamide that exhibited potent inhibitory effects on HIV-1 infectivity, with half-maximal effective concentrations (EC50) of 0.17 (9) and 0.24 µM (13), respectively. The prominent derivatives specifically interfered with the viral transcriptional step without targeting other infection step(s) and efficiently inhibited the HIV-1 replication cycle in the T cell lines and peripheral blood mononuclear cells infected with HIV-1. Additionally, compared to the wild type, the compounds exhibited similar potency against anti-retroviral drug-resistant HIV-1 strains. In a series of mode-of-action studies, the compounds inhibited the ejection of histone H3 for facilitating viral transcription on the long-terminal repeat (LTR) promoter. Furthermore, SAHA (a histone deacetylase inhibitor) treatment abolished the compound potency, revealing that the compounds can possibly target Tat-regulated epigenetic modulation of LTR to inhibit viral transcription. Overall, our screening identified novel 1,3,4-oxadiazole compounds that inhibited HIV-1 Tat, and subsequent SAR-based optimization led to the derivatives 9 and 13 development that could be a promising scaffold for developing a new class of therapeutic agents for HIV-1 infection.

Single-cell detection of copy number changes reveals dynamic mechanisms of adaptation to antifungals in Candida albicans

Genomic copy number changes are associated with antifungal drug resistance and virulence across diverse fungal pathogens, but the rate and dynamics of these genomic changes in the presence of antifungal drugs are unknown. Here we optimized a dual-fluorescent reporter system in the diploid pathogen Candida albicans to quantify haplotype-specific copy number variation (CNV) and loss of heterozygosity (LOH) at the single-cell level with flow cytometry. We followed the frequency and dynamics of CNV and LOH at two distinct genomic locations in the presence and absence of antifungal drugs in vitro and in a murine model of candidiasis. Copy number changes were rapid and dynamic during adaptation to fluconazole and frequently involved competing subpopulations with distinct genotypes. This study provides quantitative evidence for the rapid speed at which diverse genotypes arise and undergo dynamic population-level fluctuations during adaptation to antifungal drugs in vitro and in vivo.

Variant graph craft (VGC): a comprehensive tool for analyzing genetic variation and identifying disease-causing variants

BACKGROUND

The variant call format (VCF) file is a structured and comprehensive text file crucial for researchers and clinicians in interpreting and understanding genomic variation data. It contains essential information about variant positions in the genome, along with alleles, genotype calls, and quality scores. Analyzing and visualizing these files, however, poses significant challenges due to the need for diverse resources and robust features for in-depth exploration.

RESULTS

To address these challenges, we introduce variant graph craft (VGC), a VCF file visualization and analysis tool. VGC offers a wide range of features for exploring genetic variations, including extraction of variant data, intuitive visualization, and graphical representation of samples with genotype information. VGC is designed primarily for the analysis of patient cohorts, but it can also be adapted for use with individual probands or families. It integrates seamlessly with external resources, providing insights into gene function and variant frequencies in sample data. VGC includes gene function and pathway information from Molecular Signatures Database (MSigDB) for GO terms, KEGG, Biocarta, Pathway Interaction Database, and Reactome. Additionally, it dynamically links to gnomAD for variant information and incorporates ClinVar data for pathogenic variant information. VGC supports the Human Genome Assembly Hg37 and Hg38, ensuring compatibility with a wide range of data sets, and accommodates various approaches to exploring genetic variation data. It can be tailored to specific user needs with optional phenotype input data.

CONCLUSIONS

In summary, VGC provides a comprehensive set of features tailored to researchers working with genomic variation data. Its intuitive interface, rapid filtering capabilities, and the flexibility to perform queries using custom groups make it an effective tool in identifying variants potentially associated with diseases. VGC operates locally, ensuring data security and privacy by eliminating the need for cloud-based VCF uploads, making it a secure and user-friendly tool. It is freely available at https://github.com/alperuzun/VGC .

PICKLE-mediated nucleosome condensing drives H3K27me3 spreading for the inheritance of Polycomb memory during differentiation

Spreading of H3K27me3 is crucial for the maintenance of mitotically inheritable Polycomb-mediated chromatin silencing in animals and plants. However, how Polycomb repressive complex 2 (PRC2) accesses unmodified nucleosomes in spreading regions for spreading H3K27me3 remains unclear. Here, we show in Arabidopsis thaliana that the chromatin remodeler PICKLE (PKL) plays a specialized role in H3K27me3 spreading to safeguard cell identity during differentiation. PKL specifically localizes to H3K27me3 spreading regions but not to nucleation sites and physically associates with PRC2. Loss of PKL disrupts the occupancy of the PRC2 catalytic subunit CLF in spreading regions and leads to aberrant dedifferentiation. Nucleosome density increase endowed by the ATPase function of PKL ensures that unmodified nucleosomes are accessible to PRC2 catalytic activity for H3K27me3 spreading. Our findings demonstrate that PKL-dependent nucleosome compaction is critical for PRC2-mediated H3K27me3 read-and-write function in H3K27me3 spreading, thus revealing a mechanism by which repressive chromatin domains are established and propagated.

Placental somatic mutation in human stillbirth and live birth: A pilot case-control study of paired placental, fetal, and maternal whole genomes

INTRODUCTION

A high frequency of single nucleotide somatic mutations in the placenta has been recently described, but its relationship to placental dysfunction is unknown.

METHODS

We performed a pilot case-control study using paired fetal, maternal, and placental samples collected from healthy live birth controls (n = 10), live births with fetal growth restriction (FGR) due to placental insufficiency (n = 7), and stillbirths with FGR and placental insufficiency (n = 11). We quantified single nucleotide and structural somatic variants using bulk whole genome sequencing (30-60X coverage) in four biopsies from each placenta. We also assessed their association with clinical and histological evidence of placental dysfunction.

RESULTS

Seventeen pregnancies had sufficiently high-quality placental, fetal, and maternal DNA for analysis. Each placenta had a median of 473 variants (range 111-870), with 95 % arising in just one biopsy within each placenta. In controls, live births with FGR, and stillbirths, the median variant counts per placenta were 514 (IQR 381-779), 582 (450-735), and 338 (245-441), respectively. After adjusting for depth of sequencing coverage and gestational age at birth, the somatic mutation burden was similar between groups (FGR live births vs. controls, adjusted diff. 59, 95 % CI -218 to +336; stillbirths vs controls, adjusted diff. -34, -351 to +419), and with no association with placental dysfunction (p = 0.7).

DISCUSSION

We confirmed the high prevalence of somatic mutation in the human placenta and conclude that the placenta is highly clonal. We were not able to identify any relationship between somatic mutation burden and clinical or histologic placental insufficiency.

ACLY and ACSS2 link nutrient-dependent chromatin accessibility to CD8 T cell effector responses

Coordination of cellular metabolism is essential for optimal T cell responses. Here, we identify cytosolic acetyl-CoA production as an essential metabolic node for CD8 T cell function in vivo. We show that CD8 T cell responses to infection depend on acetyl-CoA derived from citrate via the enzyme ATP citrate lyase (ACLY). However, ablation of ACLY triggers an alternative, acetate-dependent pathway for acetyl-CoA production mediated by acyl-CoA synthetase short-chain family member 2 (ACSS2). Mechanistically, acetate fuels both the TCA cycle and cytosolic acetyl-CoA production, impacting T cell effector responses, acetate-dependent histone acetylation, and chromatin accessibility at effector gene loci. When ACLY is functional, ACSS2 is not required, suggesting acetate is not an obligate metabolic substrate for CD8 T cell function. However, loss of ACLY renders CD8 T cells dependent on acetate (via ACSS2) to maintain acetyl-CoA production and effector function. Together, ACLY and ACSS2 coordinate cytosolic acetyl-CoA production in CD8 T cells to maintain chromatin accessibility and T cell effector function.

Unlocking lager's flavour palette by metabolic engineering of Saccharomyces pastorianus for enhanced ethyl ester production

Despite being present in trace amounts, ethyl esters play a crucial role as flavour compounds in lager beer. In yeast, ethyl hexanoate, ethyl octanoate and ethyl decanoate, responsible for fruity and floral taste tones, are synthesized from the toxic medium chain acyl-CoA intermediates released by the fatty acid synthase complex during the fatty acid biosynthesis, as a protective mechanism. The aim of this study was to enhance the production of ethyl esters in the hybrid lager brewing yeast Saccharomyces pastorianus by improving the medium chain acyl-CoA precursor supply. Through CRISPR-Cas9-based genetic engineering, specific FAS1 and FAS2 genes harbouring mutations in domains of the fatty acid synthesis complex were overexpressed in a single and combinatorial approach. These mutations in the ScFAS genes led to specific overproduction of the respective ethyl esters: overexpression of ScFAS1I306A and ScFAS2G1250S significantly improved ethyl hexanoate production and ScFAS1R1834K boosted the ethyl octanoate production. Combinations of ScFAS1 mutant genes with ScFAS2G1250S greatly enhanced predictably the final ethyl ester concentrations in cultures grown on full malt wort, but also resulted in increased levels of free medium chain fatty acids causing alterations in flavour profiles. Finally, the elevated medium chain fatty acid pool was directed towards the ethyl esters by overexpressing the esterase ScEEB1. The genetically modified S. pastorianus strains were utilized in lager beer production, and the resulting beverage exhibited significantly altered flavour profiles, thereby greatly expanding the possibilities of the flavour palette of lager beers.

Comparison of strain specific pathogenicity of Herpes Simplex Virus Type 1 by high-throughput sequencing

Herpes Simplex Virus Type 1 (HSV-1) is a widespread human pathogen known for causing a spectrum of clinical manifestations, ranging from mild cold sores to severe complications like encephalitis. Understanding the strain-specific variations of HSV-1 is crucial for elucidating its pathogenesis and developing targeted therapeutic interventions. In this multifaceted study, we investigated the strain-specific characteristics of HSV-1 using an in vivo rat model. Firstly, a pilot study was conducted to assess the capacity of three HSV-1 strains (Fisher (F), KOS (K), and MacIntyre (M)) to induce cold sores in rats. Remarkably, the F strain exhibited pronounced pathogenicity, inducing erythema, swelling, and disrupted epidermis with ulceration, distinguishing it from the K and M strains. Subsequently, the treatment capability of intravenous acyclovir injection in HSV-1 F strain-infected rats was evaluated. Acyclovir treatment resulted in a significant reduction in HSV-1 viral copy numbers in serum and dissected neuronal tissues, particularly in the spinal cord, brain, and lower lip. Lastly, whole genome sequencing data revealed that high-impact mutations occurred in the K and M strains within the UL49, US2, and US3 genes. These mutations may play a pivotal role in influencing viral replication, dissemination, pathogenesis, and infectivity. In contrast, the moderate missense variant mutations detected in the US12, US8, UL3, UL30, UL31, and UL36 genes appeared to have no effect on viral pathogenesis and infectivity, based on RT-PCR data for spinal cord, trigeminal nerve, brain, and the lower lip. These strain-specific mutations underscore the dynamic nature of HSV-1 evolution. Collectively, our findings contribute to a deeper understanding of HSV-1 strain diversity and pave the way for the development of targeted therapeutic strategies against this medically significant virus.

Haplotype-resolved genome and mapping of freezing tolerance in the wild potato Solanum commersonii

Solanum commersonii (2n = 2x = 24, 1EBN, Endosperm Balance Number), native to the southern regions of Brazil, Uruguay, and northeastern Argentina, is the first wild potato germplasm collected by botanists and exhibits a remarkable array of traits related to disease resistance and stress tolerance. In this study, we present a high-quality haplotype-resolved genome of S. commersonii. The two identified haplotypes demonstrate chromosome sizes of 706.48 and 711.55 Mb, respectively, with corresponding chromosome anchoring rates of 94.2 and 96.9%. Additionally, the contig N50 lengths are documented at 50.87 and 45.16 Mb. The gene annotation outcomes indicate that the haplotypes encompasses a gene count of 39 799 and 40 078, respectively. The genome contiguity, completeness, and accuracy assessments collectively indicate that the current assembly has produced a high-quality genome of S. commersonii. Evolutionary analysis revealed significant positive selection acting on certain disease resistance genes, stress response genes, and environmentally adaptive genes during the evolutionary process of S. commersonii. These genes may be related to the formation of diverse and superior germplasm resources in the wild potato species S. commersonii. Furthermore, we utilized a hybrid population of S. commersonii and S. verrucosum to conduct the mapping of potato freezing tolerance genes. By combining BSA-seq analysis with traditional QTL mapping, we successfully mapped the potato freezing tolerance genes to a specific region on Chr07, spanning 1.25 Mb, with a phenotypic contribution rate of 18.81%. In short, current research provides a haplotype-resolved reference genome of the diploid wild potato species S. commersonii and establishes a foundation for further cloning and unraveling the mechanisms underlying cold tolerance in potatoes.

Complete organelle genomes of the threatened aquatic species Scheuchzeria palustris (Scheuchzeriaceae): Insights into adaptation and phylogenomic placement

Scheuchzeria palustris, the only species in the Scheuchzeriaceae family, plays a crucial role in methane production and transportation, influencing the global carbon cycle and maintaining ecosystem stability. However, it is now threatened by human activities and global warming. In this study, we generated new organelle genomes for S. palustris, with the plastome (pt) measuring 158,573 bp and the mitogenome (mt) measuring 420,724 bp. We predicted 296 RNA editing sites in mt protein-coding genes (PCGs) and 142 in pt-PCGs. Notably, abundant RNA editing sites in pt-PCGs likely originated from horizontal gene transfer between the plastome and mitogenome. Additionally, we identified positive selection signals in four mt-PCGs (atp4, ccmB, nad3, and sdh4) and one pt-PCG (rps7), which may contribute to the adaptation of S. palustris to low-temperature and high-altitude environments. Furthermore, we identified 35 mitochondrial plastid DNA (MTPT) segments totaling 58,479 bp, attributed to dispersed repeats near most MTPT. Phylogenetic trees reconstructed from mt- and pt-PCGs showed topologies consistent with the APG IV system. However, the conflicting position of S. palustris can be explained by significant differences in the substitution rates of its mt- and pt-PCGs (p < .001). In conclusion, our study provides vital genomic resources to support future conservation efforts and explores the adaptation mechanisms of S. palustris.

Analyses of combined Merkel cell carcinomas with neuroblastic components suggests that loss of T antigen expression in Merkel cell carcinoma may result in cell cycle arrest and neuroblastic transdifferentiation

Merkel cell carcinoma (MCC) is an aggressive skin cancer frequently caused by genomic integration of the Merkel cell polyomavirus (MCPyV). MCPyV-negative cases often present as combined MCCs, which represent a distinctive subset of tumors characterized by association of an MCC with a second tumor component, mostly squamous cell carcinoma. Up to now, only exceptional cases of combined MCC with neuroblastic differentiation have been reported. Herein we describe two additional combined MCCs with neuroblastic differentiation and provide comprehensive morphologic, immunohistochemical, transcriptomic, genetic and epigenetic characterization of these tumors, which both arose in elderly men and appeared as an isolated inguinal adenopathy. Microscopic examination revealed biphasic tumors combining a poorly differentiated high-grade carcinoma with a poorly differentiated neuroblastic component lacking signs of proliferation. Immunohistochemical investigation revealed keratin 20 and MCPyV T antigen (TA) in the MCC parts, while neuroblastic differentiation was confirmed in the other component in both cases. A clonal relation of the two components can be deduced from 20 and 14 shared acquired point mutations detected by whole exome analysis in both combined tumors, respectively. Spatial transcriptomics demonstrated a lower expression of stem cell marker genes such as SOX2 and MCM2 in the neuroblastic component. Interestingly, although the neuroblastic part lacked TA expression, the same genomic MCPyV integration and the same large T-truncating mutations were observed in both tumor parts. Given that neuronal transdifferentiation upon TA repression has been reported for MCC cell lines, the most likely scenario for the two combined MCC/neuroblastic tumors is that neuroblastic transdifferentiation resulted from loss of TA expression in a subset of MCC cells. Indeed, DNA methylation profiling suggests an MCC-typical cellular origin for the combined MCC/neuroblastomas. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Transcriptomic analysis of Bombyx mori in its early larval stage (2nd instar) of development upon Nosema bombycis transovarial infection

The infection caused by Nosema bombycis often known as pebrine, is a devastating sericulture disease. The infection can be transmitted to the next generation through eggs laid by infected female Bombyx mori moths (transovarial) as well as with N. bombycis contaminated food (horizontal). Most diagnoses were carried out in the advanced stages of infection until the time that infection might spread to other healthy insects. Hence, early diagnosis of pebrine is of utmost importance to quarantine infected larvae from uninfected silkworm batches and stop further spread of the infection. The findings of our study provide an insight into how the silkworm larval host defence system was activated against early N. bombycis transovarial infection. The results obtained from transcriptome analysis of infected 2nd instar larvae revealed significant (adjusted P-value < 0.05) expression of 1888 genes of which 801 genes were found to be upregulated and 1087 genes were downregulated when compared with the control. Pathway analysis indicated activation of the immune deficiency (IMD) pathway, which shows a potential immune defence response against pebrine infection as well as suppression of the melanin synthesis pathway due to lower expression of prophenoloxidase activating enzyme (PPAE). Liquid chromatography mass spectrometry (LC-MS/MS) analysis of haemolymph from infected larvae shows the secretion of serpin binding protein of N. bombycis which might be involved in the suppression of the melanization pathway. Moreover, among the differentially expressed genes, we found that LPMC-61, yellow-y, gasp and osiris 9 can be utilised as potential markers for early diagnosis of transovarial pebrine infection in B. mori. Physiological as well as biochemical roles and functions of many of the essential genes are yet to be established, and enlightened research will be required to characterize the products of these genes.

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

Guidance for estimating penetrance of monogenic disease-causing variants in population cohorts

Penetrance is the probability that an individual with a pathogenic genetic variant develops a specific disease. Knowing the penetrance of variants for monogenic disorders is important for counseling of individuals. Until recently, estimates of penetrance have largely relied on affected individuals and their at-risk family members being clinically referred for genetic testing, a 'phenotype-first' approach. This approach substantially overestimates the penetrance of variants because of ascertainment bias. The recent availability of whole-genome sequencing data in individuals from very-large-scale population-based cohorts now allows 'genotype-first' estimates of penetrance for many conditions. Although this type of population-based study can underestimate penetrance owing to recruitment biases, it provides more accurate estimates of penetrance for secondary or incidental findings. Here, we provide guidance for the conduct of penetrance studies to ensure that robust genotypes and phenotypes are used to accurately estimate penetrance of variants and groups of similarly annotated variants from population-based studies.

A novel antigen biomarker for detection of high-level of Loa loa microfilaremia

BACKGROUND

Loiasis is a disease caused by the nematode Loa loa. Serious adverse events sometimes occur in people with heavy L. loa microfilaremia after ivermectin treatment. In regions of Central Africa where loiasis is endemic, this significantly impedes global elimination programs for lymphatic filariasis and onchocerciasis that use mass distribution of ivermectin. Improved diagnostic tests to identify individuals at increased risk of serious adverse events could facilitate efforts to eliminate lymphatic filariasis and onchocerciasis in this region.

METHODS AND FINDINGS

We previously identified the L. loa protein Ll-Bhp-1 in loiasis patient sera. Here, we further characterize Ll-Bhp-1 and report development of an antigen capture ELISA to detect this antigen. This assay detected Ll-Bhp-1 in 74 of 116 (63.8%) loiasis patient sera. Ll-Bhp-1 levels were significantly correlated with L. loa microfilarial counts, and the sensitivity of the assay was highest for samples from people with high counts, (94% and 100% in people with ≥20,000 and ≥50,000 microfilaria per milliliter of blood, respectively). The antigen was not detected in 112 sera from people with other filarial infections, or in 34 control sera from the USA.

CONCLUSIONS

This Ll-Bhp-1 antigen assay is specific for loiasis, and highly sensitive for identifying people with high L. loa microfilarial counts who are at increased risk for serious adverse events after ivermectin treatment. L. loa antigen detection has the potential to facilitate loiasis mapping efforts and programs to eliminate lymphatic filariasis and onchocerciasis in Central Africa.

Genome-wide association analysis identifies candidates of three bulb traits in garlic

Garlic bulbs generally possess several swelling cloves, and the swelling degree of the bulbs determines its yield and appearance quality. However, the genetic basis underlying bulb traits remains poorly known. To address this issue, we performed a genome-wide association analysis for three bulb traits: bulb weight, diameter, and height. It resulted in the identification of 51 significant associated signals from 38 genomic regions. Twelve genes from the associated regions, whose transcript abundances in the developmental bulb showed significant correlations with the investigated traits in 81 garlic accessions, were considered the candidates of the corresponding locus. We focused on five of these candidates and their variations and revealed that the promoter variations of fructose-bisphosphate aldolase-encoding Asa8G05696.1 and beta-fructofuranosidase-encoding Asa6G01167.1 are responsible for the functional diversity of these two genes in garlic population. Interestingly, our results revealed that all candidates we focused on experienced a degree of selection during garlic evolutionary history, and different genotypes of them were retained in two China-cultivated garlic groups. Taken together, these results suggest a potential involvement of those candidates in the parallel evolution of garlic bulb organs in two China-cultivated garlic groups. This study provides important insights into the genetic basis of garlic bulb traits and their evolution.

A Pilot Study Based on the Correlation Between Whole Exome and Transcriptome Reveals Potent Variants in the Indian Population of Cervical Cancer

Cervical malignancy (CC) is the 2nd most prevalent malignancy among females, leading to cancer mortality. Primary detection of CC tumors results in an improved prognosis. CC is a malignant gynecological tumor, with few treatment options. New diagnostic and therapeutic agents are required to expand patient survival and quality of life. If CC tumors can be found at an early stage, the prognosis is much brighter. New diagnostic and therapeutic agents are needed to increase patient survival and quality of life. In this work, we performed whole-exome sequencing utilizing V5 (Illumina platform) 10 samples, 5 control and 5 CC tumour tissue, and we compared the results with transcriptome studies. KMT2C variations were shown to be among the most vicious in this analysis. From an Indian viewpoint, we found a plethora of SNVs and mutations, including those with known, unknown, and possible effects on health. Based on our findings, we know that the KMT2C gene is on chr. Seven and in exon 8, all three recognized variants are missense, synonymous, coding synonymous, non-coding variants, and GnomAD MAF (- 0.05). The variation at position (7:152265091, T > A, SNV 62478356) in KMT2C is unique, potent, and pathogenic. The missense coding transcript CIQTNF maps to chromosome 7 and displays T > C SNV. In addition, we performed single strand conformational polymorphism analysis on 64 samples and further confirmed them using Sanger sequencing to understand and verify the mutations. KMT2C shows a log FC value of - 1.16. Understanding emerging harmful mutations from an Indian viewpoint is facilitated by our bioinformatics-based, extensive correlation studies of WES analysis. Potentially harmful and new mutations were found in our preliminary analysis; among these ten top mutated genes, KMT2C and CIQTNF were altered in ten cases of CC with an Indian phenotype.

Spindle architecture constrains karyotype evolution

The eukaryotic cell division machinery must rapidly and reproducibly duplicate and partition the cell's chromosomes in a carefully coordinated process. However, chromosome numbers vary dramatically between genomes, even on short evolutionary timescales. We sought to understand how the mitotic machinery senses and responds to karyotypic changes by using a series of budding yeast strains in which the native chromosomes have been successively fused. Using a combination of cell biological profiling, genetic engineering and experimental evolution, we show that chromosome fusions are well tolerated up until a critical point. Cells with fewer than five centromeres lack the necessary number of kinetochore-microtubule attachments needed to counter outward forces in the metaphase spindle, triggering the spindle assembly checkpoint and prolonging metaphase. Our findings demonstrate that spindle architecture is a constraining factor for karyotype evolution.

Single-cell EpiChem jointly measures drug-chromatin binding and multimodal epigenome

Studies of molecular and cellular functions of small-molecule inhibitors in cancer treatment, eliciting effects by targeting genome and epigenome associated proteins, requires measurement of drug-target engagement in single-cell resolution. Here we present EpiChem for in situ single-cell joint mapping of small molecules and multimodal epigenomic landscape. We demonstrate single-cell co-assays of three small molecules together with histone modifications, chromatin accessibility or target proteins in human colorectal cancer (CRC) organoids. Integrated multimodal analysis reveals diverse drug interactions in the context of chromatin states within heterogeneous CRC organoids. We further reveal drug genomic binding dynamics and adaptive epigenome across cell types after small-molecule drug treatment in CRC organoids. This method provides a unique tool to exploit the mechanisms of cell type-specific drug actions.

Qualitative and quantitative analysis of MED12 c.887G>A causing both missense and splicing variants in X-linked Ohdo syndrome

The phenotypes associated with MED12 pathogenic variants are diverse. Male patients usually have missense variants, but the effects of base substitutions on mRNA splicing have not been investigated. Here, we report a Japanese brother with intellectual disability, characteristic facial appearance with blepharophimosis, cleft palate, Fallot tetralogy, vesicoureteral reflux, and deafness. A known missense pathogenic variant was detected in MED12, NM_005120.3:c.887G>A p.(Arg296Gln), and X-linked Ohdo syndrome was diagnosed in combination with their phenotype. mRNA splicing of MED12 was evaluated qualitatively and quantitatively using long-range PCR-based targeted RNA sequencing (reverse transcribed long amplicon sequencing), and it was shown that this missense variant simultaneously causes aberrant splicing of the 42-bp in-frame deletion in exon 7, r.847_888del, which accounts for approximately 30% of the mRNAs in both siblings. The X chromosome inactivation study showed that the X chromosome carrying the mutant allele was 100% inactivated in the carrier mothers. mRNA level analysis is essential for the accurate interpretation of the effects of variants. In this case, the MED12 protein function may be reduced by more than just an amino acid substitution, resulting in the patients with the most severe phenotype of MED12-related syndrome in males.

Comprehensive molecular characterization of collecting duct carcinoma for therapeutic vulnerability

Collecting duct carcinoma (CDC) is an aggressive rare subtype of kidney cancer with unmet clinical needs. Little is known about its underlying molecular alterations and etiology, primarily due to its rarity, and lack of preclinical models. This study aims to comprehensively characterize molecular alterations in CDC and identify its therapeutic vulnerabilities. Through whole-exome and transcriptome sequencing, we identified KRAS hotspot mutations (G12A/D/V) in 3/13 (23%) of the patients, in addition to known TP53, NF2 mutations. 3/13 (23%) patients carried a mutational signature (SBS22) caused by aristolochic acid (AA) exposures, known to be more prevalent in Asia, highlighting a geologically specific disease etiology. We further discovered that cell cycle-related pathways were the most predominantly dysregulated pathways. Our drug screening with our newly established CDC preclinical models identified a CDK9 inhibitor LDC000067 that specifically inhibited CDC tumor growth and prolonged survival. Our study not only improved our understanding of oncogenic molecular alterations of Asian CDC, but also identified cell-cycle machinery as a therapeutic vulnerability, laying the foundation for clinical trials to treat patients with such aggressive cancer.

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

Leishmania major telomerase RNA knockout: From altered cell proliferation to decreased parasite infectivity

This study focuses on the biological impacts of deleting the telomerase RNA from Leishmania major (LeishTER), a parasite responsible for causing leishmaniases, for which no effective treatment or prevention is available. TER is a critical player in the telomerase ribonucleoprotein complex, containing the template sequence copied by the reverse transcriptase component during telomere elongation. The success of knocking out both LeishTER alleles was confirmed, and no off-targets were detected. LmTER-/- cells share similar characteristics with other TER-depleted eukaryotes, such as altered growth patterns and partial G0/G1 cell cycle arrest in early passages, telomere shortening, and elevated TERRA expression. They also exhibit increased γH2A phosphorylation, suggesting that the loss of LeishTER induces DNA damage signaling. Moreover, pro-survival autophagic signals and mitochondrion alterations were shown without any detectable plasma membrane modifications. LmTER-/- retained the ability to transform into metacyclics, but their infectivity capacity was compromised. Furthermore, the overexpression of LeishTER was also deleterious, inducing a dominant negative effect that led to telomere shortening and growth impairments. These findings highlight TER's vital role in parasite homeostasis, opening discussions about its potential as a drug target candidate against Leishmania.

Endogenous nege-like viral elements in arthropod genomes reveal virus-host coevolution and ancient history of two plant virus families

Negevirus is a recently proposed taxon of arthropod-infecting virus, which is associated with plant viruses of two families (Virgaviridae and Kitaviridae). Nevertheless, the evolutionary history of negevirus-host and its relationship with plant viruses remain poorly understood. Endogenous nege-like viral elements (ENVEs) are ancient nege-like viral sequences integrated into the arthropod genomes, which can serve as the molecular fossil records of previous viral infection. In this study, 292 ENVEs were identified in 150 published arthropod genomes, revealing the evolutionary history of nege-like viruses and two related plant virus families. We discovered three novel and eight strains of nege-like viruses in 11 aphid species. Further analysis indicated that 10 ENVEs were detected in six aphid genomes, and they were divided into four types (ENVE1-ENVE4). Orthologous integration and phylogenetic analyses revealed that nege-like viruses had a history of infection of over 60 My and coexisted with aphid ancestors throughout the Cenozoic Era. Moreover, two nege-like viral proteins (CP and SP24) were highly homologous to those of plant viruses in the families Virgaviridae and Kitaviridae. CP- and SP24-derived ENVEs were widely integrated into numerous arthropod genomes. These results demonstrate that nege-like viruses have a long-term coexistence with arthropod hosts and plant viruses of the two families, Virgaviridae and Kitaviridae, which may have evolved from the nege-like virus ancestor through horizontal virus transfer events. These findings broaden our perspective on the history of viral infection in arthropods and the origins of plant viruses.

IMPORTANCE

Although negevirus is phylogenetically related to plant virus, the evolutionary history of negevirus-host and its relationship with plant virus remain largely unknown. In this study, we used endogenous nege-like viral elements (ENVEs) as the molecular fossil records to investigate the history of nege-like viral infection in arthropod hosts and the evolution of two related plant virus families (Virgaviridae and Kitaviridae). Our results showed the infection of nege-like viruses for over 60 My during the arthropod evolution. ENVEs highly homologous to viral sequences in Virgaviridae and Kitaviridae were present in a wide range of arthropod genomes but were absent in plant genomes, indicating that plant viruses in these two families possibly evolved from the nege-like virus ancestor through cross-species horizontal virus transmission. Our findings provide a new perspective on the virus-host coevolution and the origins of plant viruses.

Primary unilateral macronodular adrenal hyperplasia with concomitant glucocorticoid and androgen excess and KDM1A inactivation

BACKGROUND

Primary bilateral macronodular adrenal hyperplasia (PBMAH) is a rare cause of Cushing's syndrome. Individuals with PBMAH and glucose-dependent insulinotropic polypeptide (GIP)-dependent Cushing's syndrome due to ectopic expression of the GIP receptor (GIPR) typically harbor inactivating KDM1A sequence variants. Primary unilateral macronodular adrenal hyperplasia (PUMAH) with concomitant glucocorticoid and androgen excess has never been encountered or studied.

METHODS

We investigated a woman with a large, heterogeneous adrenal mass and severe adrenocorticotropic hormone-independent glucocorticoid and androgen excess, a biochemical presentation typically suggestive of adrenocortical carcinoma. The patient presented during pregnancy (22nd week of gestation) and reported an 18-month history of oligomenorrhea, hirsutism, and weight gain. We undertook an exploratory study with detailed histopathological and genetic analysis of the resected adrenal mass and leukocyte DNA collected from the patient and her parents.

RESULTS

Histopathology revealed benign macronodular adrenal hyperplasia. Imaging showed a persistently normal contralateral adrenal gland. Whole-exome sequencing of 4 representative nodules detected KDM1A germline variants, benign NM_001009999.3:c.136G > A:p.G46S, and likely pathogenic NM_001009999.3:exon6:c.865_866del:p.R289Dfs*7. Copy number variation analysis demonstrated an additional somatic loss of the KDM1A wild-type allele on chromosome 1p36.12 in all nodules. RNA sequencing of a representative nodule showed low/absent KDM1A expression and increased GIPR expression compared with 52 unilateral sporadic adenomas and 4 normal adrenal glands. Luteinizing hormone/chorionic gonadotropin receptor expression was normal. Sanger sequencing confirmed heterozygous KDM1A variants in both parents (father: p.R289Dfs*7 and mother: p.G46S) who showed no clinical features suggestive of glucocorticoid or androgen excess.

CONCLUSIONS

We investigated the first PUMAH associated with severe Cushing's syndrome and concomitant androgen excess, suggesting pathogenic mechanisms involving KDM1A.

Five centuries of consanguinity, isolation, health, and conflict in Las Gobas: A Northern Medieval Iberian necropolis

Between the 8th and 11th centuries CE, the Iberian Peninsula underwent profound upheaval due to the Umayyad invasion against the Visigoths, resulting in population shifts and lasting demographic impacts. Our understanding of this period is hindered by limited written sources and few archaeogenetic studies. We analyzed 33 individuals from Las Gobas, a necropolis in northern Spain, spanning the 7th to 11th centuries. By combining archaeological and osteological data with kinship, metagenomics, and ancestry analyses, we investigate conflicts, health, and demography of these individuals. We reveal intricate family relationships and genetic continuity within a consanguineous population while also identifying several zoonoses indicative of close interactions with animals. Notably, one individual was infected with a variola virus phylogenetically clustering with the northern European variola complex between ~885 and 1000 CE. Last, we did not detect a significant increase of North African or Middle East ancestries over time since the Islamic conquest of Iberia, possibly because this community remained relatively isolated.

Cell cycle-regulated transcription factor AP2XII-9 is a key activator for asexual division and apicoplast inheritance in Toxoplasma gondii tachyzoite

Toxoplasma gondii is an intracellular parasitic protozoan that poses a significant risk to the fetus carried by a pregnant woman or to immunocompromised individuals. T. gondii tachyzoites duplicate rapidly in host cells during acute infection through endodyogeny. This highly regulated division process is accompanied by complex gene regulation networks. TgAP2XII-9 is a cell cycle-regulated transcription factor, but its specific role in the parasite cell cycle is not fully understood. In this study, we demonstrate that TgAP2XII-9 is identified as a nuclear transcription factor and is dominantly expressed during the S/M phase of the tachyzoite cell cycle. Cleavage Under Targets and Tagmentation (CUT&Tag) results indicate that TgAP2XII-9 targets key genes for the moving junction machinery (RON2, 4, and 8) and daughter cell inner membrane complex (IMC). TgAP2XII-9 deficiency resulted in a significant downregulation of rhoptry proteins and rhoptry neck proteins, leading to a severe defect in the invasion and egress efficiency of tachyzoites. Additionally, the loss of TgAP2XII-9 correlated with a substantial downregulation of multiple IMC and apicoplast proteins, leading to disorders of daughter bud formation and apicoplast inheritance and further contributing to the inability of cell division and intracellular proliferation. Our study reveals that TgAP2XII-9 acts as a critical S/M-phase regulator that orchestrates the endodyogeny and apicoplast division in T. gondii tachyzoites. This study contributes to a broader understanding of the complexity of the parasite's cell cycle and its key regulators.

IMPORTANCE

The intracellular apicoplast parasite Toxoplasma gondii poses a great threat to the public health. The acute infection of T. gondii tachyzoites relies on efficient invasion by forming a moving junction structure and also fast replication by highly regulated endodyogeny. This study shows that an ApiAP2 transcription factor, TgAP2XII-9, acts as an activator for the S/M-phase gene expression, including genes related to daughter buds and moving junction formation. Loss of TgAP2XII-9 results in significant growth defects and disorders in endodyogeny and apicoplast inheritance of the parasites. Our results provide valuable insights into the transcriptional regulation of the parasite cell cycle and invading machinery in T. gondii.

Polycomb proteins RING1A/B promote H2A monoubiquitination to regulate female gametophyte development in Arabidopsis

A functional female gametophyte is the basis of successful sexual reproduction in flowering plants. During female gametophyte development, the megaspore mother cell (MMC), which differentiates from a single subepidermal somatic cell in the nucellus, undergoes meiosis to produce four megaspores; only the one at the chalazal end, referred to as the functional megaspore (FM), then undergoes three rounds of mitosis and develops into a mature embryo sac. Here, we report that RING1A and RING1B (RING1A/B), two functionally redundant Polycomb proteins in Arabidopsis, are critical for female gametophyte development. Mutations of RING1A/B resulted in defects in the specification of the MMC and the FM, and in the subsequent mitosis of the FM, thereby leading to aborted ovules. Detailed analysis revealed that several genes essential for female gametophyte development were ectopically expressed in the ring1a ring1b mutant, including Argonaute (AGO) family genes and critical transcription factors. Furthermore, RING1A/B bound to some of these genes to promote H2A monoubiquitination (H2Aub). Taken together, our study shows that RING1A/B promote H2Aub modification at key genes for female gametophyte development, suppressing their expression to ensure that the development progresses correctly.

The mitochondrial orf117Sha gene desynchronizes pollen development and causes pollen abortion in Arabidopsis Sha cytoplasmic male sterility

Cytoplasmic male sterility (CMS) is of major agronomical relevance in hybrid breeding. In gametophytic CMS, abortion of pollen is determined by the grain genotype, while in sporophytic CMS, it is determined by the mother plant genotype. While several CMS mechanisms have been dissected at the molecular level, gametophytic CMS has not been straightforwardly accessible. We used the gametophytic Sha-CMS in Arabidopsis to characterize the cause and process of pollen abortion by implementing in vivo biosensing in single pollen and mitoTALEN mutagenesis. We obtained conclusive evidence that orf117Sha is the CMS-causing gene, despite distinct characteristics from other CMS genes. We measured the in vivo cytosolic ATP content in single pollen, followed pollen development, and analyzed pollen mitochondrial volume in two genotypes that differed only by the presence of the orf117Sha locus. Our results showed that the Sha-CMS is not triggered by ATP deficiency. Instead, we observed desynchronization of a pollen developmental program. Pollen death occurred independently in pollen grains at diverse stages and was preceded by mitochondrial swelling. We conclude that pollen death is grain-autonomous in Sha-CMS and propose that mitochondrial permeability transition, which was previously described as a hallmark of developmental and environmental-triggered cell death programs, precedes pollen death in Sha-CMS.