Primer3--new capabilities and interfaces

Relationship between the expression levels of myogenic regulatory factor genes and carcass characteristics in Kivircik and Hungarian Merino lambs

This study aimed to investigate the expression profiles of the myogenic regulatory genes MYOD1, MYOG, MYF5, MYF6, and MSTN in longissimus dorsi muscle, as well as the correlation of the expression levels of these genes with carcass characteristics and growth performance in the Kivircik and Hungarian Merino sheep breeds. The expression levels of the MYF5, MYF6, and MYOG genes were found to be significantly correlated with the rib proportion, the expression level of the MYOG gene was identified as being the main determinant of variations in the rib proportion in the Kivircik lambs. The regression analysis results revealed that the expression levels of the MYF5 and MSTN genes played an essential role in determining the cold carcass dressing percentage in Hungarian Merino lambs. Further, as a result of the regression analysis, the model including the expression level of the MYF6 gene demonstrated that this gene could be responsible for 36.4% of the differences observed in cold carcass weight. In conclusion, the findings of this study suggest that the expression levels of the MYF5, MYF6, and MYOG genes were associated with various carcass traits, particularly in the Kivircik breed, and these genes hold potential as markers for enhancing breed productivity.

High levels of sinigrin trigger synthesis of fatty acids in Plutella xylostella (L.)

Diamondback moth (Lepidoptera: Plutellidae; Plutella xylostella L.) is a specialist insect of the Brassicaceae family, damaging economically important crops, such as cabbage and cauliflower. Glucosinolates, also known as 'mustard oil bombs' are present in all Brassicaceae members, of which sinigrin (allyl-glucosinolate or 2-propenyl-glucosinolate) is a major aliphatic compound. During herbivory, glucosinolates are converted to toxic isothiocyanates that deter insect pests. P. xylostella possesses glucosinolate sulfatases that desulfate them. Such a conversion renders them unfit for degradation to toxic products. Changes in the larval performance prompted us for RNA sequencing to understand probable adaptation mechanism under sinigrin stress. Differentially expressed genes were found to be related to larval cuticle proteins. Further, gene ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses depict genes belonging to the categories, integral component of membrane, cellular processes and those involved in biosynthesis of fatty acids. Upregulation of cuticular genes viz. larval cuticle protein-17 (LCP-17), cuticular protein-19 (2CP-19) and ATP binding cassette transporter C7 (ABCC7), ABCC16 was validated by qRT-PCR. Liquid chromatography quadrupole time of flight mass spectrometry analysis of whole larvae feeding on sinigrin and their separated cuticle, depicted abundance of fatty acids. Changes in the topography of the larval cuticle were evident by scanning electron microscopy. Expression of PxABCH1 was corroborated to its role in the transport of cuticular lipids. Notably, molecular docking of PxABCH1 with cuticular fatty acids showed favorable binding interactions. To summarize, integrated transcriptomic and metabolomic analyses suggest that in response to a diet containing a high dose of sinigrin, P. xylostella re-programs metabolic pathways related to fatty acid biosynthesis that directly influence insect development.

Biopsy-derived organoids in personalised early breast cancer care: Challenges of tumour purity and normal cell overgrowth cap their practical utility

The ability to establish organoids composed exclusively of tumour rather than healthy cells is essential for their implementation into clinical practice. Organoids have recently emerged as a powerful tool to expand patient material in culture and generate modifiable 3D models derived from humans or animal models. For translational research, they enable the creation of model systems for an ever-increasing number of cell types and diseases. And in personalised medicine, they potentially allow for functional drug testing with high predictive power in certain settings. We found that using biopsy material from untreated, early-stage primary breast cancer patients poses significant challenges for consistently culturing tumour cells as organoids. Specifically, we observed frequent outgrowth of genetically normal, non-cancerous epithelial cells. We analysed >100 biopsy samples from early-stage breast cancer and present our large collection of >70 organoid lines. We also show methods of assessing successful tumour cell culture in a time, and cost-efficient manner, proving a high rate (>85%) of normal cell overgrowth in early-stage breast cancer organoids. Finally, we show a number of successful attempts to culture cancer organoids from mastectomy-derived tissue of advanced, metastatic breast cancer. We conclude that the usefulness of organoids from early breast cancer for translational research and personalised medicine, especially guidance of adjuvant or post-surgical maintenance therapy, is strongly limited by the low success rate of culturing cancerous cells under organoid conditions.

Design and structure of overlapping regions in PCA via deep learning

Polymerase cycling assembly (PCA) stands out as the predominant method in the synthesis of kilobase-length DNA fragments. The design of overlapping regions is the core factor affecting the success rate of synthesis. However, there still exists DNA sequences that are challenging to design and construct in the genome synthesis. Here we proposed a deep learning model based on extensive synthesis data to discern latent sequence representations in overlapping regions with an AUPR of 0.805. Utilizing the model, we developed the SmartCut algorithm aimed at designing oligonucleotides and enhancing the success rate of PCA experiments. This algorithm was successfully applied to sequences with diverse synthesis constraints, 80.4 % of which were synthesized in a single round. We further discovered structure differences represented by major groove width, stagger, slide, and centroid distance between overlapping and non-overlapping regions, which elucidated the model's reasonableness through the lens of physical chemistry. This comprehensive approach facilitates streamlined and efficient investigations into the genome synthesis.

Analysis of gill and skin microbiota in Larimichthys crocea reveals bacteria associated with cryptocaryoniasis resistance potential

Cryptocaryoniasis, caused by the ciliate parasite Cryptocaryon irritans, poses a significant threat to the large yellow croaker (Larimichthys crocea) in intensive marine aquaculture. This study explores the interaction between skin and gill microbiota and C. irritans infection, focusing on the role of commensal microbes in disease resistance. Fish were challenged with 100 theronts per gram of body weight, leading to substantial microbial dysbiosis, characterized by decreased alpha diversity and disrupted co-occurrence networks, particularly on the skin. Post-infection, Vibrio abundance significantly increased in both gills and skin, suggesting potential for secondary infections. Conversely, lower Vibrio levels correlated with higher populations of Bdellovibrio-like organisms (BALOs), which may play a beneficial role in microbial balance. Fish showed varying susceptibility, with mildly infected individuals exhibiting less histopathological damage and a stronger immune response, indicated by elevated interleukin-1β (IL-1β) and interleukin-8 (IL-8) levels. Correlation analyses revealed significant relationships between relative infection intensity (RII) and microbial composition, with certain bacteria known for anti-eukaryotic microbial properties showing negative correlations with RII. Additionally, the abundance of nitrogen-metabolizing bacteria also correlated negatively with RII. Functional predictions indicated increased bacterial genes related to denitrification and vitamin biosynthesis post-infection. Notably, Candidatus Midichloria was identified as a potential biomarker for C. irritans infection and is thought to be an endosymbiont of C. irritans, with its presence validated through PCR analysis. These findings illuminate microbial dynamics during C. irritans infection and suggest probiotic candidates for managing cryptocaryoniasis.

Quantifying Bone Collagen Fingerprint Variation Between Species

Collagen is the most ubiquitous protein in the animal kingdom and one of the most abundant proteins on Earth. Despite having a relatively repetitive amino acid sequence motif that enables its triple helical structure, in type 1 collagen, that dominates skin and bone, there is enough variation for its increasing use for the biomolecular species identification of animal tissues processed or degraded beyond the amenability of DNA-based analyses. In recent years, this has been most commonly achieved through the technique of collagen peptide mass fingerprinting (PMF) known as ZooMS (Zooarchaeology by Mass Spectrometry), applied to the analysis of tens of thousands of samples across over one hundred studies in the past decade alone. However, a robust means to quantify variation between these fingerprints remains elusive, despite being increasingly required due to the shift towards a wider range of wild fauna and those that are more distantly related from currently known sequences. This is particularly problematic in fish due to their greater sequence variation. Here we evaluate the quantification of the relative closeness of collagen fingerprints between families using ANOSIM and a modified SIMPER analysis, incorporating relative peak intensity. Our results show a clear correlation between sequence differentiation and statistical distance of PMFs, indicating that the additional complexity of type 1 collagen in fish could directly affect the efficacy of biomolecular techniques such as ZooMS. Furthermore, this multivariate statistical analysis demonstrates that PMFs in fish are substantively more distinct than those of mammalian or amphibian taxa.

Metagenomic insights into cyanotoxin dynamics in a Mexican subtropical lake

Valle de Bravo is a vital water supply for part of the metropolitan area of the Valle de Mexico megacity, providing 30% of Mexico City's water demand. This water body has experienced an acceleration in its trophic status, going from oligotrophic to eutrophic in just a few years. This temperate lake (at a tropical latitude) is in a persistent bloom dominated by a variety of co-occurring cyanobacteria, many of which have toxigenic potential based on microscopic identification, that makes it difficult or even impractical to identify the cyanotoxin producers. To unravel this complexity and directly identify the toxigenic genera, we showed that integrating classical approaches with metagenomic is required. We first characterized, from genes to metagenomes assembled genomes, the toxigenic Cyanobacteria. We found that Microcystis was the most dominant cyanobacterial genus and the sole carrier of the mcy operon, making it the only microcystin producer. We then quantified twenty-one different cyanopeptides, including twelve microcystin congeners using a high-performance liquid chromatography-high-resolution. Nine microcystins (MCs) and the emerging cyanotoxin anabaenopeptin-A and -B were found at varying concentrations throughout the year, with MC-LA being the most common and abundant. Our findings, constrained by our sampling strategy, indicate that conventional cyanotoxin biomarkers (e.g., toxin mcy genes) were not consistently reliable indicators of cyanotoxin concentrations in this freshwater system. In this study, we followed the dynamics of the cyanobacterial community and the associated cyanopeptides with unprecedented resolution. Our results have implications for better management of toxic blooms in this freshwater system, which supplies drinking water to more than 7 million people in the megalopolis of Valle de México.

Multi-omics analysis reveals CMTR1 upregulation in cancer and roles in ribosomal protein gene expression and tumor growth

BACKGROUND

CMTR1 (cap methyltransferase 1), a key nuclear mRNA cap methyltransferase, catalyzes 2'-O-methylation of the first transcribed nucleotide, a critical step in mRNA cap formation. Previous studies have implicated CMTR1 in embryonic stem cell differentiation and immune responses during viral infection; however, its role in cancer biology remains largely unexplored. This study aims to elucidate CMTR1's function in cancer progression and evaluate its potential as a novel therapeutic target in certain cancer types.

METHODS

We conducted a comprehensive multi-omics analysis of CMTR1 across various human cancers using TCGA and CPTAC datasets. Functional studies were performed using CRISPR-mediated knockout and siRNA knockdown in human and mouse basal-like breast cancer models. Transcriptomic and pathway enrichment analyses were carried out in CMTR1 knockout/knockdown models to identify CMTR1-regulated genes. In silico screening and biochemical assays were employed to identify novel CMTR1 inhibitors.

RESULTS

Multi-omics analysis revealed that CMTR1 is significantly upregulated at the mRNA, protein, and phosphoprotein levels across multiple cancer types in the TCGA and CPTAC datasets. Functional studies demonstrated that CMTR1 depletion significantly inhibits tumor growth both in vitro and in vivo. Transcriptomic analysis of CMTR1 knockout cells revealed that CMTR1 primarily regulates ribosomal protein genes and other transcripts containing 5' Terminal Oligopyrimidine (TOP) motifs. Additionally, CMTR1 affects the expression of snoRNA host genes and snoRNAs, suggesting a broader role in RNA metabolism. Mechanistic studies indicated that CMTR1's target specificity is partly determined by mRNA structure, particularly the presence of 5'TOP motifs. Finally, through in silico screening and biochemical assays, we identified several novel CMTR1 inhibitors, including N97911, which demonstrated in vitro growth inhibition activity in breast cancer cells.

CONCLUSIONS

Our findings establish CMTR1 as an important player in cancer biology, regulating critical aspects of RNA metabolism and ribosome biogenesis. The study highlights CMTR1's potential as a therapeutic target in certain cancer types and provides a foundation for developing novel cancer treatments targeting mRNA cap methylation.

Genomics assisted mapping of earliness in pea (Pisum sativum L.)

BACKGROUND

Garden pea (Pisum sativum L.), is a temperate crop belonging to the Leguminosae family. Early maturing pea varieties complete their growth cycle in ∼80-90 days and fits very well within the crop rotation of rice, wheat, and maize, thereby providing an extra source of income to the farmers. Identification of genes associated with the earliness is very important for developing early maturing pea varieties.

METHODS AND RESULTS

In the present study we investigated the genetics of earliness and identified the putative genomic regions associated with the earliness in F2 population derived from a cross between early-maturing (Matar Ageta-10) and late-maturing (Punjab-89) pea varieties using BSA-Seq approach. Genetic analysis revealed that earliness follows a monogenic recessive inheritance pattern. Two extreme phenotypic pools were constructed by identifying ten extreme early and ten extreme late plants from the F2 population, and QTL-seq analysis was performed to obtain major genomic region of 6.5 Mb located at 418.46 Mb to 424.97 Mb on chromosome 7 and has been designated as PsE7. Further, a total of 907 SNPs were identified within this 6.5 Mb genomic region of which seven SNPs were validated through KASP markers. Among these, one marker namely PS423028253 showed association with the earliness trait at distance of 1.7 cM.

CONCLUSION

This novel genomic region along with KASP marker (PS423028253) identified in this study could be used for marker-assisted selection in pea breeding programs and will aid in the identification of the candidate genes in future studies.

High de novo mutation rate in Iranian NF2-related schwannomatosis patients with a report of a novel NF2 mutation

BACKGROUND

NF2-related schwannomatosis (NF2) is a rare genetic disease that significantly impacts patients' quality of life due to the occurrence of multiple tumors within the nervous system. The high clinical heterogeneity in tumor number, location, and size makes predicting each patient's clinical outcome impossible. Genetic investigation can be crucial in diagnosis, prognosis, and management. This study aims to explore the genetic basis of eight Iranian patients with NF2.

METHODS AND RESULTS

To investigate potential genetic causes, we conducted comprehensive medical evaluations, whole-exome sequencing (WES), and multiplex ligation-dependent probe amplification (MLPA) on the probands of each family. The identified variants in the family members were confirmed using Sanger sequencing and MLPA. The variants were classified according to the American College of Medical Genetics and Genomics guidelines. Seven distinct variants linked to the NF2 gene were identified as causes of NF2-related schwannomatosis in these patients, among which the c.862_863del frameshift was a novel variant not previously reported. Seventy-five percent of these mutations were de novo. The mean diagnostic age was lower among patients with truncating mutations compared to other patients.

CONCLUSIONS

This study identified a novel mutation in the NF2 gene and showed a high rate of de novo mutations in Iranian NF2 patients. Moreover, patients with truncating mutations experienced earlier symptoms than others. Comparing the manifestations of each patient with similar mutations to previous reports expands our understanding of the phenotype of NF2. These results can provide more comprehensive insights into prognosis and early interventions.

A TTPA deletion is associated with retinopathy with vitamin E deficiency in the English Cocker Spaniel dog

Retinopathy with vitamin E deficiency is a familial disease in the English Cocker Spaniel dog breed. Ophthalmic abnormalities observed in retinopathy with vitamin E deficiency-affected English Cocker Spaniel include lipofuscin granule deposition within the tapetal fundus and subsequent retinal degeneration resulting in visual deficits. Affected dogs may also exhibit neurological signs that include ataxia and hindlimb proprioceptive deficits. In all cases, circulating plasma concentrations of α-tocopherol are low. This study sought to investigate the genetic basis of retinopathy with vitamin E deficiency in the English Cocker Spaniel breed. We undertook a genome-wide association study comprising 30 English Cocker Spaniels with normal fundic examinations aged 6 years or older (controls) and 20 diagnosed with retinopathy with vitamin E deficiency (cases) and identified a statistically associated signal on chromosome 29 (Praw = 1.909 × 10-17). Whole genome sequencing of 2 cases identified a 102 bp deletion in exon 1 of the alpha-tocopherol transfer protein gene (TTPA), truncating the protein by 34 amino acids. The c.23_124del variant segregated with retinopathy with vitamin E deficiency in a total of 30 cases and 43 controls. Variants in TTPA are causal for ataxia with vitamin E deficiency in humans which is a phenotypically similar disease to retinopathy with vitamin E deficiency. The identification of the canine variant is extremely significant as the availability of a DNA test will allow for identification of presymptomatic dogs and early therapeutic intervention which may prevent development of retinopathy and improve neurological signs. Breeders can also use the DNA test to efficiently eradicate the disease from this breed.

Identification of DAGLB variants in Japanese early-onset Parkinson's disease

Hereditary factors play a significant role in the development of Parkinson's disease and the identification of causative genes is ongoing. Biallelic variants in Diacylglycerol lipase β (DAGLB) are related to early-onset Parkinson's disease (EOPD) in the Chinese population, and have also been identified in an Algerian case. To date, no EOPD cases with DAGLB variants have been reported among Japanese patients. This study was conducted to clarify the occurrence of DAGLB variants among Japanese EOPD patients. We screened 270 patients with sporadic EOPD (male: female ratio, 1.37:1; mean age at onset ± standard deviation, 37.32 ± 7.91 years), and 276 patients with suspected autosomal recessive Parkinson's disease (ARPD, male: female ratio, 0.75:1; mean age at onset ± standard deviation, 58.86 ± 14.67 years). Genetic screening of all coding exons and flanking splicing regions was performed by Sanger sequencing. We identified two rare biallelic variants in two patients, both from consanguineous families. One variant was a homozygous frameshift variant (c.1770_1771del, p.Tyr591ProfsTer26), which was predicted to be pathogenic. The other was a missense variant (c.1444T > C, p.Tyr482His) and was predicted to be benign, with co-segregation ruled out for this variant. We identified a pathogenic variant in the DAGLB gene. Together with previous reports, these findings provide further evidence that loss-of-function variants in DAGLB are involved in EOPD in the Japanese population.

The potential of inversions to accumulate balanced sexual antagonism is supported by simulations and Drosophila experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with In(3 R)K favoring survival and In(3 L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, In(3 L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas In(3 R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Development of ferret immune repertoire reference resources and single-cell-based high-throughput profiling assays

Domestic ferrets (Mustela putorius furo) are important for modeling human respiratory diseases. However, ferret B and T cell receptors have not been completely identified or annotated, limiting immune repertoire studies. Here, we performed long-read transcriptome sequencing of ferret splenocyte and lymph node samples to obtain over 120,000 high-quality full-length immunoglobin (Ig) and T cell receptor (TCR) transcripts. We constructed a complete reference set of the constant regions of ferret Ig and TCR isotypes and chain types. We also systematically annotated germline Ig and TCR variable (V), diversity (D), joining (J), and constant (C) genes on a recent ferret reference genome assembly. We designed new ferret-specific immune repertoire profiling assays by targeting positions in constant regions without allelic diversity across 11 ferret genome assemblies and experimentally validated them using a commercially compatible single-cell-based platform. These improved resources and assays will enable future studies to fully capture ferret immune repertoire diversity.IMPORTANCEDomestic ferrets (Mustela putorius furo) are an increasingly common model organism to study human respiratory diseases such as influenza infections. However, researchers lack ferret-specific reagents and resources to study the immune system and immune response in ferrets. In this study, we developed comprehensive ferret immune repertoire reference resources and assays, which will enable more accurate analyses of the ferret immune system in the future.

Identification of quantitative trait loci (QTLs) for key cheese making phenotypes in the blue-cheese mold Penicillium roqueforti

Elucidating the genomic architecture of quantitative traits is essential for our understanding of adaptation and for breeding in domesticated organisms. Penicillium roqueforti is the mold used worldwide for the blue cheese maturation, contributing to flavors through proteolytic and lipolytic activities. The two domesticated cheese populations display very little genetic diversity, but are differentiated and carry opposite mating types. We produced haploid F1 progenies from five crosses, using parents belonging to cheese and non-cheese populations. Analyses of high-quality genome assemblies of the parental strains revealed five large translocations, two having occurred via a circular intermediate, one with footprints of Starship giant mobile elements. Offspring genotyping with genotype-by-sequencing (GBS) revealed several genomic regions with segregation distortion, possibly linked to degeneration in cheese lineages. We found transgressions for several traits relevant for cheese making, with offspring having more extreme trait values than parental strains. We identified quantitative trait loci (QTLs) for colony color, lipolysis, proteolysis, extrolite production, including mycotoxins, but not for growth rates. Some genomic regions appeared rich in QTLs for both lipid and protein metabolism, and other regions for the production of multiple extrolites, indicating that QTLs have pleiotropic effects. Some QTLs corresponded to known biosynthetic gene clusters, e.g., for the production of melanin or extrolites. F1 hybrids constitute valuable strains for cheese producers, with new traits and new allelic combinations, and allowed identifying target genomic regions for traits important in cheese making, paving the way for strain improvement. The findings further contribute to our understanding of the genetic mechanisms underlying rapid adaptation, revealing convergent adaptation targeting major gene regulators.

A non-syndromic orofacial cleft risk locus links tRNA splicing defects to neural crest cell pathologies

Orofacial clefts are the most common form of congenital craniofacial malformation worldwide. The etiology of these birth defects is multifactorial, involving genetic and environmental factors. However, in most cases, the underlying causes remain unexplained, precluding a molecular understanding of disease mechanisms. Here, we integrated genome-wide association data, targeted resequencing of case and control cohorts, tissue- and cell-type-specific epigenomic profiling, and genome architecture analyses to molecularly dissect a genomic locus associated with an increased risk of non-syndromic orofacial cleft. We found that common and rare risk variants associated with orofacial cleft intersect with an enhancer (e2p24.2) that is active in human embryonic craniofacial tissue. We mapped e2p24.2 long-range interactions to a topologically associated domain harboring MYCN, DDX1, and CYRIA. We found that MYCN and DDX1, but not CYRIA, are required during craniofacial development in chicken embryos. We investigated the role of DDX1, a key component of the tRNA splicing complex, in cranial neural crest cells (cNCCs). The loss of DDX1 in cNCCs resulted in the accumulation of unspliced tRNA fragments, depletion of mature intron-containing tRNAs, and ribosome stalling at codons decoded by these tRNAs. This was accompanied by defects in both global protein synthesis and cNCC migration. We further showed that the induction of tRNA fragments is sufficient to disrupt craniofacial development. Together, these results uncovered a molecular mechanism in which impaired tRNA splicing affects cNCCs and craniofacial development and positioned MYCN, DDX1, and tRNA processing defects as risk factors in the pathogenesis of orofacial clefts.

WP-MOD: A multi-omics and taxonomy database for woody plants

Woody plants, including trees, shrubs, and woody vines, are vital components of terrestrial ecosystems and are critical for maintaining biodiversity, regulating climate, and supporting human livelihoods. Over the past decade, the accumulation of high-throughput sequencing data, multi-omics data, and taxonomic information on woody plants has increased significantly, highlighting the need for an integrative database. Here, we present the Woody Plant Multi-Omics Database (WP-MOD, https://www.woodyplant.com), a comprehensive and user-friendly platform designed to meet the growing need for specialized resources in woody plant research. The WP-MOD integrates extensive taxonomic information and multi-omics data from 373 species across 35 orders and provides a centralized resource for the analysis and exploration of woody plant biology. The database includes high-quality reference genomes and reanalyzed data from RNA sequencing, small RNA sequencing, chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin sequencing, and bisulfite sequencing, along with 17 tools for sequence and omics analysis. The WP-MOD supports both genetic and molecular research and contributes to the conservation and sustainable management of woody plants. We believe that the WP-MOD will be an essential tool for plant science researchers.

Bi-allelic variants in BRF2 are associated with perinatal death and craniofacial anomalies

BACKGROUND

Variants in genes encoding multiple subunits of the RNA Polymerase III complex which synthesizes rRNAs, tRNAs, and other small RNAs were previously associated with neurological disorders, such as syndromic hypomyelination leukodystrophies, pontocerebellar hypoplasia, and cerebellofaciodental syndrome. One new such candidate is BRF2, which encodes a TFIIB-like factor that recruits the RNA polymerase III complex to type 3 promoters to initiate transcription of U6, RnaseP, and 7SK RNAs.

METHODS

We combined sequencing with functional analyses to investigate the effects of BRF2 variants.

RESULTS

We observe that a previously reported significant underrepresentation of double transmission of a splice variant results in recessive lethality in three large Icelandic families with multiple perinatal losses. Using data aggregation, we identified an additional seven individuals worldwide from three unrelated families carrying biallelic variants in BRF2. Affected individuals present a variable phenotype ranging from severe craniofacial anomalies with early death to intellectual disability with motor and speech development. In silico 3D modelling and functional analyses showed functional impairment of the identified variants, e.g., differences in target loci occupancy. Zebrafish knocked down for the orthologous brf2 presented with abnormal escape response, reduced swimming velocity and head size, and craniofacial malformations. These defects were complemented by the human wild-type but not mutated BRF2 mRNA further demonstrating their deleteriousness.

CONCLUSIONS

Overall, our results support the association of biallelic BRF2 variants with a novel neurodevelopmental disease and provide an additional link between RNA polymerase III, its targets and craniofacial anomalies.

Heterogeneous and novel transcript expression in single cells of patient-derived clear cell renal cell carcinoma organoids

Splicing is often dysregulated in cancer, leading to alterations in the expression of canonical and alternatively spliced isoforms. We used the multiplexed arrays sequencing (MAS-seq) protocol of PacBio to sequence full-length transcripts in patient-derived organoid (PDO) cells of clear cell renal cell carcinoma (ccRCC). The sequencing revealed a heterogeneous dysregulation of splicing across 2599 single ccRCC cells. The majority of novel transcripts could be removed with stringent filtering criteria. The remaining 31,531 transcripts (36.6% of the 86,182 detected transcripts on average) were previously uncharacterized. In contrast to known transcripts, many of the novel transcripts have cell-specific expression. Novel transcripts common to ccRCC cells belong to genes involved in ccRCC-related pathways, such as hypoxia and oxidative phosphorylation. A novel transcript of the ccRCC-related gene nicotinamide N-methyltransferase is validated using PCR. Moreover, >50% of novel transcripts possess a predicted complete protein-coding open reading frame. An analysis of the most dominant transcript-switching events between ccRCC and non-ccRCC cells shows many switching events that are cell- and sample-specific, underscoring the heterogeneity of alternative splicing events in ccRCC. Overall, our study elucidates the intricate transcriptomic architecture of ccRCC, underlying its aggressive phenotype and providing insights into its molecular complexity.

Genome-wide comparison reveals large structural variants in cassava landraces

BACKGROUND

Structural variants (SVs) are critical for plant genomic diversity and phenotypic variation. This study investigates a large, 9.7 Mbp highly repetitive segment on chromosome 12 of TMEB117, a region not previously characterized in cassava (Manihot esculenta Crantz). We aim to explore its presence and variability across multiple cassava landraces, providing insights into its genomic significance and potential implications.

RESULTS

We validated the presence of the 9.7 Mbp segment in the TMEB117 genome, distinguishing it from other published cassava genome assemblies. By mapping short-read sequencing data from 16 cassava landraces to TMEB117 chromosome 12, we observed variability in read mapping, suggesting that while all genotypes contain the insertion region, some exhibit missing segments or sequence differences. Further analysis revealed two unique genes associated with deacetylase activity, HDA14 and SRT2, within the insertion. Additionally, the MUDR-Mutator transposable element was significantly overrepresented in this region.

CONCLUSIONS

This study uncovers a large structural variant in the TMEB117 cassava genome, highlighting its variability among different genotypes. The enrichment of HDA14 and SRT2 genes and the MUDR-Mutator elements within the insertion suggests potential functional significance, though further research is needed to explore this. These findings provide important insights into the role of structural variations in shaping cassava genomic diversity.

Cathepsin X is a conserved cell death protein involved in algal response to environmental stress

Phytoplankton are responsible for half of the global photosynthesis and form vast blooms in aquatic ecosystems. Bloom demise fuels marine microbial life and is suggested to be mediated by programmed cell death (PCD) induced by diverse environmental stressors. Despite its importance, the molecular basis for algal PCD remains elusive. Here, we reveal novel PCD genes conserved across distant algal lineages using cell-to-cell heterogeneity in the response of the diatom Phaeodactylum tricornutum to oxidative stress. Comparative transcriptomics of sorted sensitive and resilient subpopulations following oxidative stress revealed genes directly linked to their contrasting fates of cell death and survival. Comparing these genes with those found in a large-scale mutant screen in the green alga Chlamydomonas reinhardtii identified functionally relevant conserved PCD gene candidates, including the cysteine protease cathepsin X/Z (CPX). CPX mutants in P. tricornutum CPX1 and C. reinhardtii CYSTEINE ENDOPEPTIDASE 12 (CEP12) exhibited resilience to oxidative stress and infochemicals that induce PCD, supporting a conserved function of these genes in algal PCD. Phylogenetic and predictive structural analyses show that CPX is highly conserved in eukaryotes, and algae exhibit strong structural similarity to human Cathepsin X/Z (CTSZ), a protein linked to various diseases. CPX is expressed by diverse algae across the oceans and correlates with upcoming demise events during toxic Pseudo-nitzschia blooms, providing support for its ecological significance. Elucidating PCD components in algae sheds light on the evolutionary origin of PCD in unicellular organisms and on the cellular strategies employed by the population to cope with stressful conditions.

Proteome of amino acids or IGF1-stimulated pacu muscle cells offers molecular insights and suggests FN1B and EIF3C as candidate markers of fish muscle growth

Study of fish skeletal muscle is essential to understand physiological or metabolic processes, and to develop programs searching for increased muscle mass and meat production. Amino acids (AA) and IGF1 stimulate processes that lead to muscle growth, but their signaling pathways and molecular regulation need further clarification in fish. We obtained the proteome of pacu (Piaractus mesopotamicus) cultured muscle cells treated with AA or IGF1, which induced the differential abundance of 67 and 53 proteins, respectively. Enrichment analyses showed that AA modulated histone methylation, cell differentiation, and metabolism, while IGF1 modulated ATP production and protein synthesis. In addition, we identified molecular networks with candidate markers that commonly regulate fish muscle cells: FN1B and EIF3C, respectively up- and down-regulated by both treatments. FN1B was related to cell proliferation, protein synthesis, and muscle repair, while EIF3C connected with negative regulators of muscle growth. Their gene expression was evaluated in pacu and Nile tilapia (Oreochromis niloticus) after nutrient manipulation, with fn1b increased during refeeding and eif3c increased during fasting in both species. Our work helps clarify the molecular regulation by AA or IGF1 and suggests that FN1B and EIF3C could be potential stimulatory and inhibitory biomarkers of fish muscle growth.

Improved isolation and PCR detection of Phytophthora agathidicida oospores from soils

Phytophthora species are eukaryotic microorganisms responsible for severe dieback and root rot in plants worldwide, impacting crops, forests, and other important ecosystems. In New Zealand, P. agathidicida leads to fatal dieback in kauri (Agathis australis), long-lived endemic trees of significant cultural and ecological importance. A critical aspect of the P. agathidicida lifecycle is the production of oospores-thick-walled spores essential for long-term survival in soil, dispersal, and disease inoculation. However, their heterogeneous distribution in soils, robust structure, and dormant state make them challenging to detect using soil baiting or DNA-based methods. Soil baiting is the basis of most current testing for P. agathidicida, but baiting-based methods have low sensitivity, are slow, and require specialised facilities. To address these challenges, we developed and validated a PCR-based method for detecting P. agathidicida oospores directly from soil. Our approach includes a technique for separating oospores from soil, improved oospore lysis and DNA extraction, and a primer pair that targets a repeat region of the P. agathidicida genome with high sensitivity and specificity. The primers amplified the target product in all tested P. agathidicida isolates without cross-reactivity against eight non-target Phytophthora species. The detection limit was 1 femtogram of P. agathidicida DNA via endpoint PCR. Performance assessment against 65 soil samples from kauri forests revealed P. agathidicida in 69% of samples compared to only 11% detected by existing methods. By eliminating the need for baiting, our assay enhances the speed, accuracy, and accessibility of testing, thereby facilitating more comprehensive monitoring and improved disease management.

IMPORTANCE

Phytophthora species are notorious plant pathogens responsible for severe dieback and root rot diseases, significantly impacting crops, forests, and irreplaceable natural ecosystems. Rapid and accurate detection of these pathogens is essential for effective disease management. In New Zealand, P. agathidicida threatens the country's endemic kauri forests. In this study, we developed and validated a PCR-based method for detecting P. agathidicida oospores in soil. Oospores are long-lived, thick-walled spores that serve as key propagules for survival in soil and the spread of disease. Their robust structure and dormant state make them particularly challenging to detect using traditional soil baiting techniques or DNA-based methods. Our method is fast, accurate, and requires minimal equipment, enabling local testing and thereby empowering communities and enhancing surveillance efforts. Although developed for P. agathidicida, this method could be adapted for other plant pathogens, potentially improving disease management across various agricultural and ecological contexts.

Secreted retropepsin-like enzymes are essential for stress tolerance and biofilm formation in Pseudomonas aeruginosa

Proteases regulate important biological functions. Here we present the structural and functional characterization of three previously uncharacterized aspartic proteases in Pseudomonas aeruginosa. We show that these proteases have structural hallmarks of retropepsin peptidases and play redundant roles for cell survival under hypoosmotic stress conditions. Consequently, we named them retropepsin-like osmotic stress tolerance peptidases (Rlo). Our research shows that while Rlo proteases are homologous to RimB, an aspartic peptidase involved in rhizosphere colonization and plant infection, they contain N-terminal signal peptides and perform distinct biological functions. Mutants lacking all three secreted Rlo peptidases show defects in antibiotic resistance, biofilm formation, and cell morphology. These defects are rescued by mutations in the inactive transglutaminase transmembrane protein RloB and the cytoplasmic ATP-grasp protein RloC, two previously uncharacterized genes in the same operon as one of the Rlo proteases. These studies identify Rlo proteases and rlo operon products as critical factors in clinically relevant processes, making them appealing targets for therapeutic strategies against Pseudomonas infections.

Development and Application of Novel SSR Markers to Assess the Genetic Diversity and Population Structure of Phacelia secunda Along an Altitudinal Gradient in the Central Chile Andes

Phacelia secunda J.F. Gmel. (Boraginaceae) is a widely distributed insect-pollinated perennial herb. In central Chile (33° S), it occurs from the sea level up to 3600 m in the Andes, exhibiting broad morphological variation. In this study, we developed and characterized novel polymorphic microsatellites for this species, using an Illimina MiSeq sequencing platform. Nineteen polymorphic loci were obtained, with alleles numbers ranging from 3 to 13 per locus (mean = 5.84). Observed (HO) and expected heterozygosities (HE) ranged from 0.050 to 0.900 and from 0.049 to 0.825, respectively. These markers were applied to assess the genetic diversity and population structure along an altitudinal spanning from 1600 to 3600 m. The highest elevation population exhibited significantly lower within-population genetic diversity compared to lower-elevation populations. Significant population differentiation was observed along the gradient. Gene flow estimates support a stepping-stone like mode of migration, with greater exchange between adjacent elevations. These new microsatellites provide a valuable tool for elucidating the influence of altitude on genetic diversity and structure, and for evaluating the roles of local adaptation and phenotypic plasticity in shaping population variation.

CREPE (CREate Primers and Evaluate): a computational tool for large-scale primer design and specificity analysis

Polymerase chain reaction (PCR) is ubiquitous in biological research labs, as it is a fast, flexible, and cost-effective technique to amplify a DNA region of interest. Novel applications often leverage this ease of implementation and parallelize reactions for multiplexed approaches, such as next-generation sequencing. However, manual primer design can be an error-prone and time-consuming process depending on the number and composition of target sites. While Primer3 has emerged as an accessible tool to solve some of these issues and increase reproducibility, additional computational pipelines are required for appropriate scaling. Moreover, this does not replace the manual confirmation of primer specificity (i.e., the assessment of off-targets). To overcome the challenges associated with large-scale primer design, we fused the functionality of Primer3 and In-Silico PCR (ISPCR); this integrated pipeline, which we call CREPE ( CRE ate P rimers and E valuate), performs primer design and specificity analysis through a custom evaluation script for any given number of target sites at scale. Its final output summarizes the lead forward and reverse primer pair for each target site, a measure of the likelihood of binding to off-targets, and additional information to aid a user's decision-making. We provide this through a customized workflow for targeted amplicon sequencing (TAS) on a 150bp paired-end Illumina platform. Experimental testing of this application on clinically relevant loci showed successful amplification for more than 90% of primers deemed acceptable by CREPE. Together, we believe that CREPE represents a useful bioinformatic tool that supports the important scaling of PCR-based applications.

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Tracking Different States of Spiked Environmental DNA Using Multiplex Digital PCR Assays

The study of microbial communities based on the total environmental DNA (eDNA) is influenced by the presence of different eDNA states, i.e., intracellular (iDNA) and extracellular DNA (exDNA), and the choice of the DNA extraction method. Although the use of spike-and-recovery controls facilitates the diagnosis of such issues, appropriate experimental setups simultaneously accounting for the different eDNA states and their bacterial origins are missing. Here, we used two single-gene deletion mutants of both Escherichia coli and Bacillus subtilis to trace exDNA and iDNA spike-ins of each selected model organism within environmental samples. Unique primer/probe sets were developed for each strain, allowing their absolute quantification using multiplex digital PCR assays. The proposed spike-and-recovery controls were successfully applied to various environments including soil, sediment, sludge and compost. While the percent recovery of spiked iDNA differed significantly between E. coli and B. subtilis, results were similar for both model organisms in the case of spiked exDNA, emphasising that the fate of DNA molecules in the environment is similar irrespective of their bacterial origin. Hence, future studies may benefit from the proposed approach to better understand methodological ambiguities related to the eDNA extraction in general as well as the separation of the different eDNA states.

Comparative analysis of esophageal gland microbes between two body sizes of Gigantopelta aegis, a hydrothermal snail from the Southwest Indian Ridge

Microbial communities within animals provide nutritional foundation and energy supply for the hydrothermal ecosystem. The peltospirid snail Gigantopelta aegis forms large aggregation in the Longqi vent field on the Southwest Indian Ridge. This endemic species is characterized by a changeable diet and morphology, especially reflected in internal organs such as remarkably enlarged esophageal glands. Here, 16S full-length rRNA gene analysis was performed to compare the variations in esophageal gland microbiota between two body size groups (small and large) of G. aegis. Phyla Proteobacteria and Bacteroidetes were the dominant featured bacteria contributing to the microbial community. No significant differences between the small and large groups were revealed by the diversity index and principal component analysis (PCA) clustering. The differences were in the relative abundance of bacteria. Compared with small-sized snails, the larger ones housed more Thiogranum (9.94% to 34.86%) and fewer Sediminibacterium (29.38% to 4.54%). Functional prediction for all of the microbiota showed that the pathways related to metabolism appeared highly abundant in smaller G. aegis. However, for the larger ones, the most distinctive pathways were those of environmental information processing. Facultative symbiotic Sulfurovum was marked as a core node in the co-occurrence network and suggested an influence on habitat selection of G. aegis in hydrothermal fields. In summary, variations in bacteria composition and potential functions possibly reflected changes in the anatomical structure and dietary habits of G. aegis. These dominant bacteria shared capabilities in nutritional supplementation and ecological niche expansion in the host, potentially a key adaptation for hydrothermal survival.IMPORTANCEDominant in the Longqi hydrothermal vent Southwest Indian Ridge, Gigantopelta aegis was observed to undergo unique and significant morphological changes and diet shifts known as cryptometamorphosis. During this process, G. aegis developed a specialized bacteria-housing organ, the esophageal gland, in the later life stages. Our research discovered variations in esophageal gland microbes between different body size groups of snails. These bacteria were closely related to the development and health of G. aegis. Full-length 16S rRNA gene analysis revealed more Thiogranum and fewer Sediminibacterium, suggesting a potential association with environmental adaptation. In the small-sized group, the potential functions were enriched in metabolism, while in larger G. aegis individuals, predictions indicated adaptive functions such as environmental information processing. Also, symbiotic Sulfurovum could be one of the factors influencing the habitat selection of G. aegis. Understanding the complex relationship between benthic macrofauna and microbes helps us describe the mechanisms of survival in extreme environments.

A portable, nanopore-based genotyping platform for near real-time detection of Puccinia graminis f. sp. tritici lineages and fungicide sensitivity

BACKGROUND

Fungal plant disease outbreaks are increasing in both scale and frequency, posing severe threats to agroecosystem stability, native biodiversity and food security. Among these, the notorious wheat stem rust fungus, Puccinia graminis f.sp. tritici (Pgt), has threatened wheat production since the earliest days of agriculture. New Pgt strains continue to emerge and quickly spread over vast distances through the airborne dispersal of asexual urediniospores, triggering extensive disease outbreaks as these exotic Pgt strains often overcome resistance in dominant crop varieties of newly affected regions. This highlights the urgent need for a point-of-care, real-time Pgt genotyping platform to facilitate early detection of emerging Pgt strains.

RESULTS

In this study, we developed a simple amplicon-based re-sequencing platform for rapid genotyping of Pgt isolates. This system is built around a core set of 276 Pgt genes that we found are highly polymorphic between Pgt isolates and showed that the sequence of these genes alone could be used to accurately type Pgt strains to particular lineages. We also developed a simplistic DNA preparation method and an automated bioinformatic pipeline, to enable these Pgt gene markers to be sequenced and analysed rapidly using the MinION nanopore sequencing device. This approach successfully enabled the typing of Pgt strains within approximately 48 h of collecting Pgt-infected wheat samples, even in resource-limited locations in Kenya and Ethiopia. In addition, we incorporated monitoring capabilities for sequence variations in Pgt genes that encode targets of the azole and succinate dehydrogenase inhibitor fungicides, enabling real-time tracking of potential shifts in fungicide sensitivity.

CONCLUSION

The newly developed Pgt Mobile And Real-time, PLant disEase (MARPLE) diagnostics platform we established, now allows precise typing of individual Pgt strains while simultaneously tracking changes in fungicide sensitivity, providing an early warning system for potential indicators of changes in the Pgt population and emerging fungicide resistance. Further integration of this Pgt MARPLE diagnostics platform into national surveillance programmes will support more informed management decisions and timely responses to Pgt disease outbreaks, helping reduce the devastating crop losses currently caused by this 'cereal killer'.

OGU: A Toolbox for Better Utilising Organelle Genomic Data

Organelle genomes serve as crucial datasets for investigating the genetics and evolution of plants and animals, genome diversity, and species identification. To enhance the collection, analysis, and visualisation of such data, we have developed a novel open-source software tool named Organelle Genome Utilities (OGU). The software encompasses three modules designed to streamline the handling of organelle genome data. The data collection module is dedicated to retrieving, validating and organising sequence information. The evaluation module assesses sequence variance using a range of methods, including novel metrics termed stem and terminal phylogenetic diversity. The primer module designs universal primers for downstream applications. Finally, a visualisation pipeline has been developed to present comprehensive insights into organelle genomes across different lineages rather than focusing solely on individual species. The performance, compatibility and stability of OGU have been rigorously evaluated through benchmarking with four datasets, including one million mixed GenBank records, plastid genomic data from the Lamiaceae family, mitochondrial data from rodents, and 308 plastid genomes sourced from various angiosperm families. Based on software capabilities, we identified 30 plastid intergenic spacers. These spacers exhibit a moderate evolutionary rate and offer practical utility comparable to coding regions, highlighting the potential applications of intergenic spacers in organelle genomes. We anticipate that OGU will substantially enhance the efficient utilisation of organelle genomic data and broaden the prospects for related research endeavours.

Alu-Mediated Deletion of FANCA in Turkish Families With Fanconi Anemia: Evidence of a Founder Effect

Fanconi anemia (FA) is a rare inherited bone marrow failure syndrome characterized by pancytopenia, increased susceptibility to malignancies, and a spectrum of congenital anomalies. Here, we report on eight affected individuals from six unrelated families with a large Alu-mediated intragenic deletion encompassing exons 6-31 in the FANCA gene, identified as a founder mutation in the Turkish population through haplotype analysis. This deletion, mediated by Alu repeat sequences, underscores the role of repetitive elements in FA pathogenesis. Clinical data revealed variable phenotypic presentations among affected individuals, highlighting the challenge of establishing genotype-phenotype correlations even in the presence of identical FANCA pathogenic variants. We carried out an easy and effective PCR-based diagnostic test for detecting this mutation, enabling precise diagnosis and genetic counseling for affected individuals and their families. This study provides valuable insights into the molecular mechanisms underlying FA pathogenesis and offers a practical approach for genetic diagnosis in affected individuals, particularly those of Turkish descent.

Cyclin-dependent kinase 13 is indispensable for normal mouse heart development

Congenital heart disease (CHD) has an incidence of approximately 1%. Over the last decade, sequencing studies including large cohorts of individuals with CHD have begun to unravel the genetic mechanisms underpinning CHD. This includes the identification of variants in cyclin-dependent kinase 13 (CDK13), in individuals with syndromic CHD. CDK13 encodes a serine/threonine protein kinase. The cyclin partner of CDK13 is cyclin K; this complex is thought to be important in transcription and RNA processing. Pathogenic variants in CDK13 cause CDK13-related disorder in humans, characterised by intellectual disability and developmental delay, recognisable facial features, feeding difficulties and structural brain defects, with 35% of individuals having CHD. To obtain a greater understanding for the role that this essential protein kinase plays in embryonic heart development, we have analysed a presumed loss of function Cdk13 transgenic mouse model (Cdk13tm1b). The homozygous mutants were embryonically lethal in most cases by E15.5. X-gal staining showed Cdk13 expression localised to developing facial regions, heart and surrounding areas at E10.5, whereas at E12.5, it was more widely present. In the E15.5 heart, staining was seen throughout. RT-qPCR showed significant reduction in Cdk13 transcript expression in homozygous compared with WT and heterozygous hearts at E10.5 and E12.5. Detailed morphological 3D analysis of embryonic and postnatal hearts was performed using high-resolution episcopic microscopy, which affords a more detailed analysis of structures such as cardiac valve leaflets and endocardial cushions, compared with more traditional histological techniques. We show that both the homozygous and heterozygous Cdk13tm1b mutants exhibit a range of CHD, including ventricular septal defects, bicuspid aortic valve, double outlet right ventricle and atrioventricular septal defects. 100% (n = 4) of homozygous hearts displayed CHD. Differential expression was seen in Cdk13tm1b homozygous mutants for two genes known to be necessary for normal heart development. The types of defects, and the presence of CHD in heterozygous mice (17.02%, n = 8/47), are consistent with the CDK13-related disorder phenotype in humans. This study provides important insights into the effects of reduced function of CDK13 in the mouse heart and contributes to our understanding of the mechanism behind this disorder as a cause of CHD.

Circular RNAs in neurological conditions - computational identification, functional validation, and potential clinical applications

Non-coding RNAs (ncRNAs) have gained significant attention in recent years due to advancements in biotechnology, particularly high-throughput total RNA sequencing. These developments have led to new understandings of non-coding biology, revealing that approximately 80% of non-coding regions in the genome possesses biochemical functionality. Among ncRNAs, circular RNAs (circRNAs), first identified in 1976, have emerged as a prominent research field. CircRNAs are abundant in most human cell types, evolutionary conserved, highly stable, and formed by back-splicing events which generate covalently closed ends. Notably, circRNAs exhibit high expression levels in neural tissue and perform diverse biochemical functions, including acting as molecular sponges for microRNAs, interacting with RNA-binding proteins to regulate their availability and activity, modulating transcription and splicing, and even translating into functional peptides in some cases. Recent advancements in computational and experimental methods have enhanced our ability to identify and validate circRNAs, providing valuable insights into their biological roles. This review focuses on recent developments in circRNA research as they related to neuropsychiatric and neurodegenerative conditions. We also explore their potential applications in clinical diagnostics, therapeutics, and future research directions. CircRNAs remain a relatively underexplored area of non-coding biology, particularly in the context of neurological disorders. However, emerging evidence supports their role as critical players in the etiology and molecular mechanisms of conditions such as schizophrenia, bipolar disorder, major depressive disorder, Alzheimer's disease, and Parkinson's disease. These findings suggest that circRNAs may provide a novel framework contributing to the molecular dysfunctions underpinning these complex neurological conditions.

The First Report of Cod Gill Poxvirus in Gills of Atlantic Cod (Gadus morhua L.) Suffering From Cardiorespiratory Disease

Atlantic cod farming experiences renewed growth in Norway, and increased awareness is essential to address emerging diseases in this species. There are few reports on gill diseases in cod, and to date, no viral gill infections of cod have been documented. In this study, we collected samples from three sequential time points in summer 2023 from farmed cod suffering from cardiorespiratory disease. We document severe heart and gill pathology, and a novel double-strand DNA virus was discovered in the gills. Through comprehensive genetic characterisation and comparative sequence analysis, this virus was classified as a new species in the genus Poxvirus, designated cod gill poxvirus (CGPV). We demonstrate disease causality with severe gill lesions as shown by histopathology and RNA scope in situ hybridisation, and poxvirus particles were identified in gill epithelial cells by transmission electron microscopy. Further, another gill pathogen not previously described in cod, Candidatus. Branchiomonas cysticola, was identified by pcr and in situ hybridisation. Our findings provide strong evidence for poxvirus in Atlantic cod and underscore the imminent threat posed to Atlantic cod farm industry.

Identification of Novel F9 Gene Variants in 143 Vietnamese Patients with Hemophilia B

Purpose

Vietnam is estimated to have approximately 30,000 hemophilia B (HB) carriers, with hundreds of new cases registered annually. However, comprehensive molecular studies on HB remain limited. Therefore, this study aimed to characterize genetic variants and assess their clinical significance in unrelated Vietnamese patients with HB.

Patients and Methods

This study included a cohort of 143 unrelated HB patients with diagnosed FIX levels. Genetic analysis of the F9 gene was performed using DNA sequencing and other molecular techniques. Variant pathogenicity was classified following ACMG/AMP guidelines, supplemented by computational predictions and clinical data.

Results

A 100% variant detection rate was achieved, identifying 83 unique variants from 143 patients. Single nucleotide variants were predominant, with missense variants accounting for 71.08%. Of the 83 unique variants, 20 novel variants were identified, including six missenses, four nonsenses, four frameshifts, two large deletions, two in-frame deletions, and two splice-site variants. The serine protease domain contained the highest proportion of variants (49.4%). Pathogenicity analysis revealed a predominance of severe phenotypes (72.03%). Among the novel variants, twelve were classified as pathogenic, one as likely pathogenic, and seven as variants of uncertain significance. A noteworthy case was the NM_000133.4:c.-21C>T promoter variant associated with HB Leyden, which demonstrated age-dependent improvements in factor IX levels.

Conclusion

This study expands the mutational spectrum of HB in the Vietnamese population and provide critical insights into genotype-phenotype correlations. The identification of novel variants enhances diagnostic precision and underscores the importance of comprehensive genomic analyses in understanding disease mechanisms.

Replenishing co-downregulated miR-100-5p and miR-125b-5p in malignant germ cell tumors causes growth inhibition through cell cycle disruption

MicroRNAs (miRNAs) are short, nonprotein-coding RNAs, and their expression is dysregulated in malignant germ cell tumors (GCTs). Here, we investigated the causes and consequences of downregulated miR-99a-5p/miR-100-5p (functionally identical) and miR-125b-5p levels in malignant GCTs regardless of age, site, or subtype. Quantitative RT-PCR was used to assess miR-99a-5p/miR-100-5p, miR-125b-5p, and associated gene expression in malignant GCT tissues/cell lines [seminoma (Sem), yolk sac tumor (YST), embryonal carcinoma (EC)]. Cells were treated with demethylating 5-azacytidine and pyrosequencing was performed. Combination miR-100-5p/miR-125b-5p mimic replenishment was used to treat malignant GCT cells. Global messenger RNA (mRNA) targets of the replenished miRNAs were identified and Metascape used to study pathway effects. We found that expression levels of miR-99a-5p/miR-100-5p and miR-125b-5p, their respective pri-miRNAs, and associated genes from chromosomes 11 and 21 (chr11/chr21) were downregulated and highly correlated in malignant GCT cells. Treatment with 5-azacytidine caused upregulation of these miRNAs, with pyrosequencing revealing hypermethylation of their chr11/chr21 loci, likely contributing to miR-100-5p/miR-125b-5p downregulation. Combination miR-100-5p/miR-125b-5p mimic replenishment resulted in growth inhibition in Sem/YST cells, with miR-100-5p/miR-125b-5p mRNA targets enriched in downregulated genes, which were involved in cell cycle (confirmed by flow cytometry) and signaling pathways. Knockdown of the miR-100-5p/miR-125b-5p target tripartite motif containing 71 (TRIM71kd) recapitulated miR-100-5p/miR-125b-5p replenishment, with growth inhibition and cell cycle disruption of Sem/YST/EC cells. Further, replenishment led to reduced lin-28 homolog A (LIN28A) levels and concomitant increases in let-7 (MIRLET7B) tumor suppressor miRNAs, creating a sustained reversion of cell phenotype. In summary, combination miR-100-5p/miR-125b-5p mimic replenishment or TRIM71kd caused growth inhibition in malignant GCT cells via cell cycle disruption. Further studies are now warranted, including mimic treatment alongside conventional platinum-based chemotherapy.

Streptococcus thermophilus CNCM I-5570 lysate counteracts the aging process in human dermal fibroblast cells by neutralizing harmful free radicals and impacting antioxidant and anti-inflammatory pathways, thus restoring their physiological functions

Previous studies have highlighted the in vitro and in vivo anti-aging potential of Streptococcus thermophilus prompting us to investigate the biomolecular mechanisms underlying its effects. We evaluated the reparative ability of S. thermophilus lysate in a hydrogen peroxide (H2O2)-induced senescence model of human dermal fibroblasts (HDFs). Cell proliferation, cell number, and senescence level were evaluated by IncuCyte® Live Cell Imager system, trypan blue dye exclusion test and β-galactosidase activity, respectively. We analyzed p21, prolyl 4-hydroxylase A1, intracellular collagen I, nuclear factor E2-related factor 2 (Nrf2), nuclear factor kappa B (NF-κB) and heme oxygenase-1 expression through western blot. Extracellular levels of collagen I, interleukin-1β, and IL-6 were assessed by ELISA. The oxidative stress markers were assayed using standard methods. The direct antioxidant activity of probiotic was quantified using multiple techniques. The presence of antioxidant genes in probiotic was detected via PCR assay. Probiotic lysate exposure increased the proliferation rate, counteracted the aging by reducing β-galactosidase activity and p21 levels, promoted collagen I synthesis and neutralized oxidative stress by activating Nrf2. The probiotic lysate inhibited the NF-κB pathway with pro-inflammatory marker downregulation. Notably, we revealed that probiotic exhibited strong free radical scavenging ability, iron-chelating properties, and significant ferric reducing power in a concentration-dependent manner. We identified seven genes with antioxidant function in its genome. Our results show that S. thermophilus lysate is efficacious in suppressing the biomolecular events associated with H2O2-induced cellular aging, thus supporting the reparative action of S. thermophilus, helpful in treating skin aging.

Combining the RCAS/tv-a retrovirus and CRISPR/Cas9 gene editing systems to generate primary mouse models of diffuse midline glioma

Diffuse midline gliomas (DMGs) are lethal brain tumors that arise in children and young adults, resulting in a median survival of less than two years. Genetically engineered mouse models (GEMMs) are critical to studying tumorigenesis and tumor-immune interactions, which may inform new treatment approaches. However, current midline glioma GEMM approaches are limited in their ability to multiplex perturbations and/or target specific cell lineages in the brain for genetic manipulation. Here, we combined the RCAS/tv-a avian retrovirus system and CRISPR/Cas9 genetic engineering to drive midline glioma formation in mice. CRISPR/Cas9-based disruption of Trp53, a tumor suppressor that is frequently disrupted in midline gliomas, along with the oncogene PDGF-B resulted in high grade tumor formation with moderate latency (median time to tumor formation of 12 weeks). We confirmed CRISPR-mediated Trp53 disruption using next-generation sequencing (NGS) and immunohistochemistry (IHC). Next, we disrupted multiple midline glioma tumor suppressor genes (Trp53, Pten, Atm, Cdkn2a) in individual mouse brains. These mini-pooled in vivo experiments generated primary midline gliomas with decreased tumor latency (median time to tumor formation of 3.6 weeks, P < 0.0001, log-rank test compared to single-plex gRNA). Quantification of gRNA barcodes and CRISPR editing events revealed that all tumors contained cells with various disruptions of all target genes and suggested a multiclonal origin for the tumors as well as stronger selection for Trp53 disruption compared to disruption of the other genes. This mouse modeling approach will streamline midline glioma research and enable complex experiments to understand tumor evolution and therapeutics.

Whole-genome sequencing, as a powerful diagnostic tool in hearing loss, reveals novel variants in PTPRQ missed by whole-exome sequencing

BACKGROUND/OBJECTIVES

Hearing loss (HL) is one of the most common congenital disorders, affecting 1-2 in 1,000 newborns. Modern genetic diagnostics using large gene panels and/or whole exome analysis (WES) can identify disease-causing mutations in 25-50 % of patients, with higher solve rates in individuals with earlier onset.

RESULTS

Here, we used whole-genome sequencing (WGS) to reanalyze 14 index patients/families who remained without genetic diagnosis by WES. We were able to identify the genetic cause of HL in 6 families ( ∼ 43 %). Two families were diagnosed with DFNB84A caused by compound heterozygous recessive mutations in PTPRQ. Three of the four underlying variants, including a structural variant, a deep intronic variant, and a splice variant, escaped detection by WES. Minigene assays confirmed the pathogenicity of the intronic and the splice variants. In addition, we used protein 3D structure prediction and rigid ligand docking to study the pathogenicity of variants that escape nonsense-mediated decay.

CONCLUSION

In our study, we present four novel variants in PTPRQ, three of which were detected only by WGS. To our knowledge, we report here the first pathogenic deep intronic PTPRQ variant causing HL. Our results suggest that the mutational spectrum of PTPRQ is not well covered by standard WES and that PTPRQ-associated hearing loss may be more frequent than previously thought. WGS provides an additional layer of information in the diagnostics of HL.

Feeding time modulates the daily rhythms of expression of digestive and metabolic enzymes in the liver, and food intake regulation and reward systems in the hypothalamus of the European sea bass (Dicentrarchus labrax)

Fish exhibit daily rhythms at the molecular level across different tissues, synchronized by zeitgebers, such as food availability. To optimize feeding, organisms align internal timekeeping systems to environmental cues. Previous studies on intermediary metabolism and the hypothalamic control of food intake in fish have underscored the significance of feeding time and daily rhythms. This study examined how feeding times-mid-light (ML) versus mid-dark (MD)-influence the rhythmic transcription of digestive and metabolic enzymes in the liver, and regulatory factors of food intake in the hypothalamus of European sea bass (Dicentrarchus labrax). It also explored the connection between food intake control and the reward system. When fish were fed at ML, genes involved in protein digestion (tryp2, tryp3, ctrl, and cpa5) exhibited daily rhythms with peaks early in the dark phase (ZT 11:17-13:36). These peaks were delayed in MD-fed fish (ZT 16:57-18:27). Pla2, a gene related to lipid metabolism, and transamination genes (c-alt, m-alt) showed rhythms only in ML-fed fish, with acrophases in the light phase (ZT 5:01-13:58), such as pyruvate kinase (pk) that peaked at ZT 6:16. Orexigenic genes (npy, orexin) had rhythms only in the MD group, with nocturnal peaks (ZT 13:09, 16:06). Conversely, reward system genes (th, bdnf) were rhythmic in ML-fed fish (ZT 17:35, 11:46), with only th retaining its rhythm in MD-fed fish (ZT 15:30). These findings suggest feeding time significantly affects rhythms in digestive and metabolic processes. They also highlight the intricate nature of food intake regulation systems, which present diverse synchronization patterns in relation to feeding time.

High-throughput screening reveals high diversity and widespread distribution of viruses in black soldier flies (Hermetia illucens)

Virus discovery in mass-reared insects is a growing topic of interest due to outbreak risks and for insect welfare concerns. In the case of black soldier flies (Hermetia illucens, BSF), pioneering bioinformatic studies have uncovered exogenous viruses from the orders Ghabrivirales and Bunyavirales, as well as endogenous viral elements from five virus families. This prompted further virome investigation of BSF metagenomes and metatranscriptomes, including from BSF individuals displaying signs and symptoms of disease. A high-throughput pipeline allowed the simultaneous investigation of 203 next generation sequencing datasets. This revealed the presence of seven viruses belonging to the families Dicistroviridae, Iflaviridae, Rhabdoviridae, Solinviviridae, Inseviridae, Lebotiviridae, and an unclassified Bunyavirales. Here we describe five viruses, which were detected in BSF from multiple origins, outlining the diversity of naturally occurring viruses associated with BSF colonies. As this viral community may also include BSF pathogens, we developed molecular detection tools which could be used for viral surveillance, both in mass-reared and wild populations of BSF.

Developing Striga resistance in sorghum by modulating host cues through CRISPR/Cas9 gene editing

KEY MESSAGE

High transformation and gene editing efficiencies in sorghum-produced, transgene-free SDN1-edited plants exhibit precise mutations, reduced germination stimulants, and enhanced resistance to Striga infection. Sorghum (Sorghum bicolor L.) is a primary food staple grain for millions in Sub-Saharan Africa (SSA). It is mainly constrained by the parasitic weed Striga, which causes up to 100% yield losses and affects over 60% of cultivable farmlands and livelihoods. In this study, CRISPR/Cas9 technology is utilized to induce mutations in core strigolactone (SL) biosynthetic genes, i.e., CCD7, CCD8, MAX1, in addition to an uncharacterized gene (DUF) in the fine-mapped 400 kb lgs1 region in sorghum to develop durable Striga resistance. Two sorghum cultivars were delivered with the expression cassettes through immature embryo-based Agrobacterium-mediated transformation. Our study demonstrated transformation and gene editing efficiencies of ~ 70 and up to 17.5% (calculated based on the numuber of established plants), respectively, in two sorghum genotypes. Subsequent analysis of homozygous E0 lines in the E1 generation confirmed stable integration of mutations for all targeted genes. Loss-of-function mutations in the CCD7, CCD8, MAX1, and DUF genes led to a significant downregulation of the expression of associated genes in the SL biosynthetic pathway. The phenotypic analysis of edited lines revealed changes in phenotypic patterns compared to wild-type plants. Analysis of root exudates showed significant reductions in SL production in edited lines compared to wild-type plants. Striga infection experiments demonstrated delayed or reduced emergence rates of Striga in edited lines with lower SL production, highlighting the potential for genetically altering SL production to control Striga infestations. This study provides insights into the functional roles of CCD7, CCD8, MAX1, and DUF genes in sorghum towards reduced and/or altered SL production and improved resistance to Striga infestations.

Chromosome-scale genome assembly of Trigonella corniculata (L.)L. (Nagauri pan /Kasuri methi), an important spice

Trigonella corniculata (L) L. or Nagauri pan /Kasuri methi, is an important spice crop with high nutraceutical potential. We report the de novo chromosome-scale assembly of T. corniculata genome using high coverage PacBio, Illumina and Hi-C reads. The assembly spans 798 Mb (Megabases) in 282 scaffolds with a scaffold N50 of 99.6 Mb. More than 98% of the sequence length is captured in eight different pseudomolecules with an average length of 98 Mb. A BUSCO score of over 97% is suggestive of the high degree of completeness and contiguity of the genome. A total of 64,801 protein-coding genes are predicted. Genome-wide Simple Sequence Repeats (99,149) have been identified and wet lab validated at forty-eight loci. The chromosome-scale genome assembly of T. corniculata and the SSR markers identified in this study will provide a strong foundation for future structural and functional genomics studies in T. corniculata and other fenugreek species.

The RSPO2 gene is associated with bilateral anterior amelia in Chihuahuas

Bilateral anterior amelia (BAA) is the congenital absence of thoracic limbs and has been reported in the Chihuahua as an autosomal recessive disorder. In some cases, the digits of the pelvic limbs can be variably affected, but otherwise, the pelvic limbs are generally spared. A GWAS performed with nine BAA affected Chihuahuas identified a significant association on chromosome 13, and homozygosity mapping delineated a 2.1 Mb chromosomal region containing the RSPO2 gene. Loss of function variants of RSPO2 in humans and cattle has been associated with the absence of all limbs. Six affected Chihuahuas were whole genome sequenced (WGS) and aligned to the CanFam4 assembly. SNVs, small indels, and structural variants within the critical interval that fitted a recessive model were investigated. Three SNVs (NC_049234.1:g.8891861C > T; NC_049234.1:g.8974204C > T and NC_049234.1:g.9789424G > A) were homozygous in five cases and absent from 3,418 genetically diverse control genome sequences, except for one Small Poodle that was heterozygous. One SNV resided in RSPO2's second intron, while the two others were intergenic. The three candidate variants were genotyped in 7 additional cases and 100 control Chihuahuas. Twelve of 13 cases were homozygous for the mutant allele, and one case was heterozygous. Controls were either homozygous for the reference allele (97%) or heterozygous (3%). Our data should facilitate genetic testing of Chihuahuas to prevent the unintentional production of BAA affected dogs. Moreover, the identification of these variants enhances understanding of RSPO2 gene function in limb development.

Application of genomic tools to study and potentially improve the upper thermal tolerance of farmed Atlantic salmon (Salmo salar)

BACKGROUND

The Atlantic salmon (Salmo salar) aquaculture industry must mitigate the impacts of rising ocean temperatures and the increased prevalence/severity of marine heat waves. Therefore, we investigated the genetic architecture and gene expression (transcriptomics) responsible for determining a salmon's upper thermal tolerance.

RESULTS

A genome-wide association study (GWAS) was conducted using fin clips of salmon from a previous incremental thermal maximum (ITMax) challenge (n = 251) and the North American 50 K SNP chip. ITMax was a highly polygenic trait with low/moderate heritability (mean SNP-based h2 = 0.20 and pedigree-based h2 = 0.25). Using data from the same fish, a separate GWAS assessed thermal-unit growth coefficient (TGC). Five significant SNPs were detected on chromosomes three and five, and high heritability estimates were calculated for TGC measured as fish grew from 12 to 20 °C (mean SNP-based h2 = 0.62 and pedigree-based h2 = 0.64). RNA-seq analyses of liver samples (n = 5-6 family-1 temperature-1) collected from the four most and four least tolerant families at 10 and 20 °C were also used to provide insights into potential mechanisms modulating this species' thermal tolerance. Between the top and bottom families, 347 and 175 differentially expressed transcripts (FDR-adjusted p < 0.01; fold-change ≥|2.0|) were identified at 10 and 20 °C, respectively. GO term enrichment analysis revealed unique responses to elevated temperature between family rankings (e.g., 'blood coagulation', 'sterol metabolic process' and 'synaptic growth at neuromuscular junction'). qPCR analyses further confirmed differences pertaining to cholesterol metabolism (lpl), inflammation (epx, elf3, ccl20), apoptosis (htra1b, htra2, anxa5b), angiogenesis (angl4, pdgfa), nervous system processes (insyn2a, kcnj11l) and heat stress (serpinh1b-1, serpinh1b-2). Three differentially expressed transcripts (i.e., ppp1r9a, gal3st1a, f5) were located in close proximity (± 120 kbp) to near-significant SNPs from the GWAS. Interestingly, ppp1r9a and gal3st1a have putative neurological functions, while f5 regulates blood coagulation.

CONCLUSIONS

These analyses provide several putative biomarkers of upper thermal tolerance in salmon that could prove valuable in helping the industry develop more temperature-tolerant fish. Further, our study supports previous reports that ITMax has low/moderate heritability in this species, and suggests that TGC at elevated temperatures is highly heritable.

Cytokine interactions and chemokine dysregulations in mastitis immunopathogenesis: insights from transcriptomic profiling of milk somatic cells in tropical Sahiwal (Bos indicus) cows

Introduction

Bovine mastitis causes a significant loss to the dairy industry by affecting the quantity and quality of milk. Addressing this challenge, the present study will leverage advanced omics techniques for early mastitis detection in early lactating Sahiwal cows (Bos indicus). This was the first differential transcriptomic study investigating the alterations in gene expression in milk somatic cells during the progression of naturally occurring mastitis in indigenous Sahiwal cows.

Methods

Cows were grouped into healthy (H), subclinical mastitis (SCM) and clinical mastitis (CM) groups by thoroughly screening them using the California Mastitis Test (CMT) and milk somatic cell counts (SCC). This was followed by detailed milk composition analysis, differential leukocyte counts (DLC), and microbiological culture.

Results

The differential gene expression of milk SCs through transcriptome profiling identified 83 and 76, up-regulated and 157 and 192 down-regulated genes in CM vs H and SCM vs H groups (log2 fold change ≥1 and ≤-1, p < 0.05) respectively. Pathway analysis revealed that upregulated genes were enriched in pathways such as phagosome activity, IL-17 signalling, Th1 and Th2 cell differentiation, while downregulated genes were linked to RIG-I-like receptor signalling, NK cell cytotoxicity, and Toll-like receptor signalling and Cytokine-cytokine receptor interactions. Notably, the study underscores the roles of chemokines CCL8, CCL2, and CXCL10 in immune cell recruitment during mastitis, where their downregulation suggests impaired mammary immune defense that governs Chemokine signalling pathways. Further, the comparative analysis with the previously available milk SCs proteome data identified the downregulation of chemokines signalling pathways during mastitis.

Discussion

Overall, this research enhances our understanding of mastitis pathogenesis and emphasizes that these targeted chemokines may boost mammary resilience through immunomodulation, genetic selection and genome editing or by utilising adjuvants in vaccine development that restore chemokine signalling offers a potential strategy to improve mastitis resistance in dairy cattle.

Functional modulation of the human gut microbiome by bacteria vehicled by cheese

Since cheese is one of the most commonly and globally consumed fermented foods, scientific investigations in recent decades have focused on determining the impact of this dairy product on human health and well-being. However, the modulatory effect exerted by the autochthonous cheese microbial community on the taxonomic composition and associated functional potential of the gut microbiota of human is still far from being fully dissected or understood. Here, through the use of an in vitro human gut-simulating cultivation model in combination with multi-omics approaches, we have shown that minor rather than dominant bacterial players of the cheese microbiota are responsible for gut microbiota modulation of cheese consumers. These include taxa from the genera Enterococcus, Bacillus, Clostridium, and Hafnia. Indeed, they contribute to expand the functional potential of the intestinal microbial ecosystem by introducing genes responsible for the production of metabolites with relevant biological activity, including genes involved in the synthesis of vitamins, short-chain fatty acids, and amino acids. Furthermore, tracing of cheese microbiota-associated bacterial strains in fecal samples from cheese consumers provided evidence of horizontal transmission events, enabling the detection of particular bacterial strains transferred from cheese to humans. Moreover, transcriptomic and metabolomic analyses of a horizontally transmitted (cheese-to-consumer) bacterial strain, i.e., Hafnia paralvei T10, cultivated in a human gut environment-simulating medium, confirmed the concept that cheese-derived bacteria may expand the functional arsenal of the consumer's gut microbiota. This highlights the functional and biologically relevant contributions of food microbes acquired through cheese consumption on the human health.IMPORTANCEDiet is universally recognized as the primary factor influencing and modulating the human intestinal microbiota both taxonomically and functionally. In this context, cheese, being a fermented food with its own microbiota, serves not only as a source of nourishment for humans, but also as a source of nutrients for the consumer's gut microbiota. Additionally, it may act as a vehicle for autochthonous food-associated microorganisms which undergo transfer from cheese to the consumer, potentially influencing host gut health. The current study highlights not only that cheese microbiota-associated bacteria can be traced in the human gut microbiota, but also that they may expand the functional repertoire of the human gut microbiota, with potentially significant implications for human health.

Distinct mutations emerge in the genome of serotype O foot-and-mouth disease virus during persistence in cattle

Like other RNA viruses, foot-and-mouth disease virus (FMDV) has a high mutation rate. After the acute phase of infection, about half of infected cattle develop a persistent FMDV infection that can last for weeks or months. During this persistent phase, the virus continues to replicate, resulting in the emergence of genomic heterogeneity. We have documented the pattern of mutations in the persistent phase by obtaining consensus-level sequences directly from oropharyngeal fluid (OPF) without prior virus isolation in culture. OPF samples were repeatedly collected from 22 experimentally infected cattle, 20 of which were virus positive in the OPF on day 21 after infection or later. We observed that during the persistent phase, the amount of non-synonymous mutations causing an amino acid change increased over time. Two amino acid changes that showed a striking increase during the persistent phase, VP3 A75T and VP2 Y79H, were present neither in the inoculum nor during the acute phase. Another amino acid change in VP3, R56C, which was previously implicated in FMDV pathogenicity, was already present in the inoculum and dominated toward the end of the trial in most samples. Several other amino acid changes occurred, particularly on the surface of VP2 around residue VP2 79. By functional analysis, we show that the persistent isolates evolve distinctly compared with cell culture adaptation but do not show signs of antigenic escape from neutralizing antibodies. In agreement with previous observations, we conclude that these amino acid changes are indeed associated with persistent infection of cattle with FMDV serotype O.

IMPORTANCE

Our research article describes the genetic changes that occur during the acute and persistent foot-and-mouth disease (FMDV) infection. This is of particular interest to understand viral dynamics within an infected population from which new viral strains could emerge. Especially FMDV, with its high antigenic diversity and very limited cross-reactivity between strains and serotypes, has already demonstrated in the past that new variants can quickly emerge and evade vaccine responses. In our study, we have observed that this dynamic evolution continues during the persistent phase. Persistently infected animals, which are clinically indistinguishable from healthy animals, also pose a reservoir for recombination. A better understanding of viral dynamics is essential for improved vaccines to prevent the emergence of antigenic variants.

A conserved lysine/arginine-rich motif is essential for the autophagic degradation of potyviral 6K1 protein and virus infection

Potyviruses possess one positive-sense single-stranded RNA genome, mainly dependent on polyprotein processing as the expression strategy. The resulting polyproteins are proteolytically processed by three virus-encoded proteases into 11 or 12 mature proteins. One such factor, 6 kDa peptide 1 (6K1), is an understudied viral factor. Its function in viral infection remains largely mysterious. This study is to reveal part of its roles by using pepper veinal mottle virus (PVMV) as the model. Alanine substitution screening analysis revealed that 15 of 17 conserved residues across potyviral 6K1 sequences are essential for PVMV infection. However, 6K1 protein is less accumulated in virus-infected cells, although P3-6K1 and 6K1-CI junctions are efficiently processed by NIa-Pro for its release, indicating that 6K1 undergoes a self-degradation event. Mutating the cleavage site to prevent NIa-Pro processing abolishes viral infection, suggesting that the generation of 6K1 along with its degradation might be important for viral multiplication. We corroborated that cellular autophagy is engaged in 6K1's degradation. Individual engineering of the 15 6K1 variants into PVMV allows their expression along with viral infection. Five of such variants, D30A, V32A, K34A, L36A, and L39A, significantly interfere with viral infection. The five residues are enclosed in a conserved lysine/arginine-rich motif; four of them appear crucial in engaging autophagy-mediated self-degradation. Based on these data, we envisaged a scenario which potyviral 6K1s interact with an unknown anti-viral component to be co-degraded by autophagy to promote viral infection.IMPORTANCEPotyvirus is the largest genus of plant-infecting RNA viruses, which encompasses socio-economically important virus species, such as Potato virus Y, Plum pox virus, and Soybean mosaic virus. Like all picorna-like viruses, potyviruses express their factors mainly via polyprotein processing. Theoretically, viral factors P3 through CP, including 6K1, should share an equivalent number of molecules. The 6K1 is small in size (~6 kDa) and conserved across potyviruses but less accumulated in virus-infected cells. This study demonstrates that cellular autophagy is engaged in the degradation of 6K1 to promote viral infection. In particular, we found a conserved lysine/arginine-rich motif in 6K1s across potyviruses that is engaged in this degradation event. This finding reveals one facet of a small protein that helps understand the pro-viral role of cellular autophagy in viral infection.