TERRA ONTseq: a long-read-based sequencing pipeline to study the human telomeric transcriptome.

The long noncoding RNA TERRA is transcribed from telomeres in virtually all eukaryotes with linear chromosomes. In humans, TERRA transcription is driven in part by promoters comprising CpG dinucleotide-rich repeats of 29 bp repeats, believed to be present in half of the subtelomeres. Thus far, TERRA expression has been analyzed mainly using molecular biology-based approaches that only generate partial and somehow biased results. Here, we present a novel experimental pipeline to study human TERRA based on long-read sequencing (TERRA ONTseq). By applying TERRA ONTseq to different cell lines, we show that the vast majority of human telomeres produce TERRA and that the cellular levels of TERRA transcripts vary according to their chromosomes of origin. Using TERRA ONTseq, we also identified regions containing TERRA transcription start sites (TSSs) in more than half of human subtelomeres. TERRA TSS regions are generally found immediately downstream from 29 bp repeat-related sequences, which appear to be more widespread than previously estimated. Finally, we isolated a novel TERRA promoter from the highly expressed subtelomere of the long arm of Chromosome 7. With the development of TERRA ONTseq, we provide a refined picture of human TERRA biogenesis and expression and we equip the scientific community with an invaluable tool for future studies.

NanoDeep: a deep learning framework for nanopore adaptive sampling on microbial sequencing.

Nanopore sequencers can enrich or deplete the targeted DNA molecules in a library by reversing the voltage across individual nanopores. However, it requires substantial computational resources to achieve rapid operations in parallel at read-time sequencing. We present a deep learning framework, NanoDeep, to overcome these limitations by incorporating convolutional neural network and squeeze and excitation. We first showed that the raw squiggle derived from native DNA sequences determines the origin of microbial and human genomes. Then, we demonstrated that NanoDeep successfully classified bacterial reads from the pooled library with human sequence and showed enrichment for bacterial sequence compared with routine nanopore sequencing setting. Further, we showed that NanoDeep improves the sequencing efficiency and preserves the fidelity of bacterial genomes in the mock sample. In addition, NanoDeep performs well in the enrichment of metagenome sequences of gut samples, showing its potential applications in the enrichment of unknown microbiota. Our toolkit is available at https://github.com/lysovosyl/NanoDeep.

Uncovering the dynamics and consequences of RNA isoform changes during neuronal differentiation.

Static gene expression programs have been extensively characterized in stem cells and mature human cells. However, the dynamics of RNA isoform changes upon cell-state-transitions during cell differentiation, the determinants and functional consequences have largely remained unclear. Here, we established an improved model for human neurogenesis in vitro that is amenable for systems-wide analyses of gene expression. Our multi-omics analysis reveals that the pronounced alterations in cell morphology correlate strongly with widespread changes in RNA isoform expression. Our approach identifies thousands of new RNA isoforms that are expressed at distinct differentiation stages. RNA isoforms mainly arise from exon skipping and the alternative usage of transcription start and polyadenylation sites during human neurogenesis. The transcript isoform changes can remodel the identity and functions of protein isoforms. Finally, our study identifies a set of RNA binding proteins as a potential determinant of differentiation stage-specific global isoform changes. This work supports the view of regulated isoform changes that underlie state-transitions during neurogenesis.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code.

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

Genomic and transcriptomic analyses of the elite rice variety Huizhan provide insight into disease resistance and heat tolerance.

The indica rice variety Huizhan shows elite traits of disease resistance and heat tolerance. However, the underlying genetic basis of these traits is not fully understood due to limited genomic resources. Here, we used Nanopore long-read and next-generation sequencing technologies to generate a chromosome-scale genome assembly of Huizhan. Comparative genomics analysis uncovered a large chromosomal inversion and expanded gene families that are associated with plant growth, development and stress responses. Functional rice blast resistance genes, including Pi2, Pib and Ptr, and bacterial blight resistance gene Xa27, contribute to disease resistance of Huizhan. Furthermore, integrated genomics and transcriptomics analyses showed that OsHIRP1, OsbZIP60, the SOD gene family, and various transcription factors are involved in heat tolerance of Huizhan. The high-quality genome assembly and comparative genomics results presented in this study facilitate the use of Huizhan as an elite parental line in developing rice varieties adapted to disease pressure and climate challenges.

Sequencing accuracy and systematic errors of nanopore direct RNA sequencing.

BACKGROUND: Direct RNA sequencing (dRNA-seq) on the Oxford Nanopore Technologies (ONT) platforms can produce reads covering up to full-length gene transcripts, while containing decipherable information about RNA base modifications and poly-A tail lengths. Although many published studies have been expanding the potential of dRNA-seq, its sequencing accuracy and error patterns remain understudied. RESULTS: We present the first comprehensive evaluation of sequencing accuracy and characterisation of systematic errors in dRNA-seq data from diverse organisms and synthetic in vitro transcribed RNAs. We found that for sequencing kits SQK-RNA001 and SQK-RNA002, the median read accuracy ranged from 87% to 92% across species, and deletions significantly outnumbered mismatches and insertions. Due to their high abundance in the transcriptome, heteropolymers and short homopolymers were the major contributors to the overall sequencing errors. We also observed systematic biases across all species at the levels of single nucleotides and motifs. In general, cytosine/uracil-rich regions were more likely to be erroneous than guanines and adenines. By examining raw signal data, we identified the underlying signal-level features potentially associated with the error patterns and their dependency on sequence contexts. While read quality scores can be used to approximate error rates at base and read levels, failure to detect DNA adapters may be a source of errors and data loss. By comparing distinct basecallers, we reason that some sequencing errors are attributable to signal insufficiency rather than algorithmic (basecalling) artefacts. Lastly, we generated dRNA-seq data using the latest SQK-RNA004 sequencing kit released at the end of 2023 and found that although the overall read accuracy increased, the systematic errors remain largely identical compared to the previous kits. CONCLUSIONS: As the first systematic investigation of dRNA-seq errors, this study offers a comprehensive overview of reproducible error patterns across diverse datasets, identifies potential signal-level insufficiency, and lays the foundation for error correction methods.

Decoding the complexity of on-target integration: characterizing DNA insertions at the CRISPR-Cas9 targeted locus using nanopore sequencing.

BACKGROUND: CRISPR-Cas9 technology has advanced in vivo gene therapy for disorders like hemophilia A, notably through the successful targeted incorporation of the F8 gene into the Alb locus in hepatocytes, effectively curing this disorder in mice. However, thoroughly evaluating the safety and specificity of this therapy is essential. Our study introduces a novel methodology to analyze complex insertion sequences at the on-target edited locus, utilizing barcoded long-range PCR, CRISPR RNP-mediated deletion of unedited alleles, magnetic bead-based long amplicon enrichment, and nanopore sequencing. RESULTS: We identified the expected F8 insertions and various fragment combinations resulting from the in vivo linearization of the double-cut plasmid donor. Notably, our research is the first to document insertions exceeding ten kbp. We also found that a small proportion of these insertions were derived from sources other than donor plasmids, including Cas9-sgRNA plasmids, genomic DNA fragments, and LINE-1 elements. CONCLUSIONS: Our study presents a robust method for analyzing the complexity of on-target editing, particularly for in vivo long insertions, where donor template integration can be challenging. This work offers a new tool for quality control in gene editing outcomes and underscores the importance of detailed characterization of edited genomic sequences. Our findings have significant implications for enhancing the safety and effectiveness of CRISPR-Cas9 gene therapy in treating various disorders, including hemophilia A.

Benchmarking short and long read polishing tools for nanopore assemblies: achieving near-perfect genomes for outbreak isolates.

BACKGROUND: Oxford Nanopore provides high throughput sequencing platforms able to reconstruct complete bacterial genomes with 99.95% accuracy. However, even small levels of error can obscure the phylogenetic relationships between closely related isolates. Polishing tools have been developed to correct these errors, but it is uncertain if they obtain the accuracy needed for the high-resolution source tracking of foodborne illness outbreaks. RESULTS: We tested 132 combinations of assembly and short- and long-read polishing tools to assess their accuracy for reconstructing the genome sequences of 15 highly similar Salmonella enterica serovar Newport isolates from a 2020 onion outbreak. While long-read polishing alone improved accuracy, near perfect accuracy (99.9999% accuracy or ~ 5 nucleotide errors across the 4.8 Mbp genome, excluding low confidence regions) was only obtained by pipelines that combined both long- and short-read polishing tools. Notably, medaka was a more accurate and efficient long-read polisher than Racon. Among short-read polishers, NextPolish showed the highest accuracy, but Pilon, Polypolish, and POLCA performed similarly. Among the 5 best performing pipelines, polishing with medaka followed by NextPolish was the most common combination. Importantly, the order of polishing tools mattered i.e., using less accurate tools after more accurate ones introduced errors. Indels in homopolymers and repetitive regions, where the short reads could not be uniquely mapped, remained the most challenging errors to correct. CONCLUSIONS: Short reads are still needed to correct errors in nanopore sequenced assemblies to obtain the accuracy required for source tracking investigations. Our granular assessment of the performance of the polishing pipelines allowed us to suggest best practices for tool users and areas for improvement for tool developers.

Molecular Epidemiology of Mayaro Virus among Febrile Patients, Roraima State, Brazil, 2018-2021.

We detected Mayaro virus (MAYV) in 3.4% (28/822) of febrile patients tested during 2018-2021 from Roraima State, Brazil. We also isolated MAYV strains and confirmed that these cases were caused by genotype D. Improved surveillance is needed to better determine the burden of MAYV in the Amazon Region.

Applications of Nanopore sequencing in precision cancer medicine.

Oxford Nanopore Technologies sequencing, also referred to as Nanopore sequencing, stands at the forefront of a revolution in clinical genetics, offering the potential for rapid, long read, and real-time DNA and RNA sequencing. This technology is currently making sequencing more accessible and affordable. In this comprehensive review, we explore its potential regarding precision cancer diagnostics and treatment. We encompass a critical analysis of clinical cases where Nanopore sequencing was successfully applied to identify point mutations, splice variants, gene fusions, epigenetic modifications, non-coding RNAs, and other pivotal biomarkers that defined subsequent treatment strategies. Additionally, we address the challenges of clinical applications of Nanopore sequencing and discuss the current efforts to overcome them.

Exploring the Genomic Landscape of Bacillus paranthracis PUMB_17 as a Proficient Phosphatidylcholine-Specific Phospholipase C Producer.

Phospholipases find versatile applications across industries, including detergent production, food modification, pharmaceuticals (especially in drug delivery systems), and cell signaling research. In this study, we present a strain of Bacillus paranthracis for the first time, demonstrating significant potential in the production of phosphatidylcholine-specific phospholipase C (PC-PLC). The investigation thoroughly examines the B. paranthracis PUMB_17 strain, focusing on the activity of PC-PLC and its purification process. Notably, the PUMB_17 strain displays extracellular PC-PLC production with high specific activity during the late exponential growth phase. To unravel the genetic makeup of PUMB_17, we employed nanopore-based whole-genome sequencing and subsequently conducted a detailed genome annotation. The genome comprises a solitary circular chromosome spanning 5,250,970 bp, featuring a guanine-cytosine ratio of 35.49. Additionally, two plasmids of sizes 64,250 bp and 5845 bp were identified. The annotation analysis reveals the presence of 5328 genes, encompassing 5186 protein-coding sequences, and 142 RNA genes, including 39 rRNAs, 103 tRNAs, and 5 ncRNAs. The aim of this study was to make a comprehensive genomic exploration that promises to enhance our understanding of the previously understudied and recently documented capabilities of Bacillus paranthracis and to shed light on a potential use of the strain in the industrial production of PC-PLC.

Closing the gap: Oxford Nanopore Technologies R10 sequencing allows comparable results to Illumina sequencing for SNP-based outbreak investigation of bacterial pathogens.

Whole-genome sequencing has become the method of choice for bacterial outbreak investigation, with most clinical and public health laboratories currently routinely using short-read Illumina sequencing. Recently, long-read Oxford Nanopore Technologies (ONT) sequencing has gained prominence and may offer advantages over short-read sequencing, particularly with the recent introduction of the R10 chemistry, which promises much lower error rates than the R9 chemistry. However, limited information is available on its performance for bacterial single-nucleotide polymorphism (SNP)-based outbreak investigation. We present an open-source workflow, Prokaryotic Awesome variant Calling Utility (PACU) (https://github.com/BioinformaticsPlatformWIV-ISP/PACU), for constructing SNP phylogenies using Illumina and/or ONT R9/R10 sequencing data. The workflow was evaluated using outbreak data sets of Shiga toxin-producing Escherichia coli and Listeria monocytogenes by comparing ONT R9 and R10 with Illumina data. The performance of each sequencing technology was evaluated not only separately but also by integrating samples sequenced by different technologies/chemistries into the same phylogenomic analysis. Additionally, the minimum sequencing time required to obtain accurate phylogenetic results using nanopore sequencing was evaluated. PACU allowed accurate identification of outbreak clusters for both species using all technologies/chemistries, but ONT R9 results deviated slightly more from the Illumina results. ONT R10 results showed trends very similar to Illumina, and we found that integrating data sets sequenced by either Illumina or ONT R10 for different isolates into the same analysis produced stable and highly accurate phylogenomic results. The resulting phylogenies for these two outbreaks stabilized after ~20 hours of sequencing for ONT R9 and ~8 hours for ONT R10. This study provides a proof of concept for using ONT R10, either in isolation or in combination with Illumina, for rapid and accurate bacterial SNP-based outbreak investigation.

A multicenter study on accuracy and reproducibility of nanopore sequencing-based genotyping of bacterial pathogens.

Nanopore sequencing has shown the potential to democratize genomic pathogen surveillance due to its ease of use and low entry cost. However, recent genotyping studies showed discrepant results compared to gold-standard short-read sequencing. Furthermore, although essential for widespread application, the reproducibility of nanopore-only genotyping remains largely unresolved. In our multicenter performance study involving five laboratories, four public health-relevant bacterial species were sequenced with the latest R10.4.1 flow cells and V14 chemistry. Core genome MLST analysis of over 500 data sets revealed highly strain-specific typing errors in all species in each laboratory. Investigation of the methylation-related errors revealed consistent DNA motifs at error-prone sites across participants at read level. Depending on the frequency of incorrect target reads, this either leads to correct or incorrect typing, whereby only minimal frequency deviations can randomly determine the final result. PCR preamplification, recent basecalling model updates and an optimized polishing strategy notably diminished the non-reproducible typing. Our study highlights the potential for new errors to appear with each newly sequenced strain and lays the foundation for computational approaches to reduce such typing errors. In conclusion, our multicenter study shows the necessity for a new validation concept for nanopore sequencing-based, standardized bacterial typing, where single nucleotide accuracy is critical.

Comprehensive analysis of m6 A methylome and transcriptome by Nanopore sequencing in clear cell renal carcinoma.

N6 -methyladenosine (m6 A) is the most prevalent epigenetic modification on eukaryotic messenger RNAs. Recent studies have focused on elucidating the key role of m6 A modification patterns in tumor progression. However, the relationship between m6 A and transcriptional regulation remains elusive. Nanopore technology enables the quantification of m6 A levels at each genomic site. In this study, a pair of tumor tissues and adjacent normal tissues from clear cell renal cell carcinoma (ccRCC) surgical samples were collected for Nanopore direct RNA sequencing. We identified 9644 genes displaying anomalous m6 A modifications, with 5343 genes upregulated and 4301 genes downregulated. Among these, 5224 genes were regarded as dysregulated genes, encompassing abnormal regulation of both m6 A modification and RNA expression. Gene Set Enrichment Analysis revealed an enrichment of these genes in pathways related to renal system progress and fatty acid metabolic progress. Furthermore, the χ2 test demonstrated a significant association between the levels of m6 A in dysregulated genes and their transcriptional expression levels. Additionally, we identified four obesity-associated genes (FTO, LEPR, ADIPOR2, and NPY5R) among the dysregulated genes. Further analyses using public databases revealed that these four genes were all related to the prognosis and diagnosis of ccRCC. This study introduced the novel approach of employing conjoint analysis of m6 A modification and RNA expression based on Nanopore sequencing to explore potential disease-related genes. Our work demonstrates the feasibility of the application of Nanopore sequencing technology in RNA epigenetic regulation research and identifies new potential therapeutic targets for ccRCC.

Real-time application of ITS and D1-D3 nanopore amplicon metagenomic sequencing in fungal infections: Enhancing fungal infection diagnostics.

While fungal infections cause considerable morbidity and mortality, the performance of the current diagnostic tests for fungal infection is low. Even though fungal metagenomics or targeted next-generation sequencing have been investigated for various clinical samples, the real-time clinical utility of these methods still needs to be elucidated. In this study, we used internal transcribed spacer (ITS) and D1-D3 ribosomal DNA nanopore amplicon metagenomic sequencing to assess its utility in patients with fungal infections. Eighty-four samples from seventy-three patients were included and categorized into 'Fungal infection,' 'Fungal colonization,' and 'Fungal contamination' groups based on the judgement of infectious disease specialists. In the 'Fungal infection' group, forty-seven initial samples were obtained from forty-seven patients. Three fungal cases detected not by the sequencing but by conventional fungal assays were excluded from the analysis. In the remaining cases, the conventional fungal assay-negative/sequencing-positive group (n=11) and conventional fungal assay-positive/sequencing-positive group (n=33) were compared. Non-Candida and non-Aspergillus fungi infections were more frequent in the conventional-negative/sequencing-positive group (p-value = 0.031). We demonstrated the presence of rare human pathogens, such as Trichosporon asahii and Phycomyces blakesleeanus. In the 'Fungal infection' group and 'Fungal colonization' group, sequencing was faster than culturing (mean difference = 4.92 days, p-value < 0.001/ mean difference = 4.67, p-value <0.001). Compared to the conventional diagnostic methods including culture, nanopore amplicon sequencing showed a shorter turnaround time and a higher detection rate for uncommon fungal pathogens.

Optimized bacterial community characterization through full-length 16S rRNA gene sequencing utilizing MinION nanopore technology.

BACKGROUND: Accurate identification of bacterial communities is crucial for research applications, diagnostics, and clinical interventions. Although 16S ribosomal RNA (rRNA) gene sequencing is a widely employed technique for bacterial taxonomic classification, it often results in misclassified or unclassified bacterial taxa. This study sought to refine the full-length 16S rRNA gene sequencing protocol using the MinION sequencer, focusing on the V1-V9 regions. Our methodological enquiry examined several factors, including the number of PCR amplification cycles, choice of primers and Taq polymerase, and specific sequence databases and workflows employed. We used a microbial standard comprising eight bacterial strains (five gram-positive and three gram-negative) in known proportions as a validation control. RESULTS: Based on the MinION protocol, we employed the microbial standard as the DNA template for the 16S rRNA gene amplicon sequencing procedure. Our analysis showed that an elevated number of PCR amplification cycles introduced PCR bias, and the selection of Taq polymerase and primer sets significantly affected the subsequent analysis. Bacterial identification at genus level demonstrated Pearson correlation coefficients ranging from 0.73 to 0.79 when assessed using BugSeq, Kraken-Silva and EPI2ME-16S workflows. Notably, the EPI2ME-16S workflow exhibited the highest Pearson correlation with the microbial standard, minimised misclassification, and increased alignment accuracy. At the species taxonomic level, the BugSeq workflow was superior, with a Pearson correlation coefficient of 0.92. CONCLUSIONS: These findings emphasise the importance of careful selection of PCR settings and a well-structured analytical framework for 16S rRNA full-length gene sequencing. The results showed a robust correlation between the predicted and observed bacterial abundances at both the genus and species taxonomic levels, making these findings applicable across diverse research contexts and with clinical utility for reliable pathogen identification.

Identification of transcriptionally-active human papillomavirus integrants through nanopore sequencing reveals viable targets for gene therapy against cervical cancer.

Integration of the human papillomavirus (HPV) genome into the cellular genome is a key event that leads to constitutive expression of viral oncoprotein E6/E7 and drives the progression of cervical cancer. However, HPV integration patterns differ on a case-by-case basis among related malignancies. Next-generation sequencing technologies still face challenges for interrogating HPV integration sites. In this study, utilizing Nanopore long-read sequencing, we identified 452 and 108 potential integration sites from the cervical cancer cell lines (CaSki and HeLa) and five tissue samples, respectively. Based on long Nanopore chimeric reads, we were able to analyze the methylation status of the HPV long control region (LCR), which controls oncogene E6/E7 expression, and to identify transcriptionally-active integrants among the numerous integrants. As a proof of concept, we identified an active HPV integrant in between RUNX2 and CLIC5 on chromosome 6 in the CaSki cell line, which was supported by ATAC-seq, H3K27Ac ChIP-seq, and RNA-seq analysis. Knockout of the active HPV integrant, by the CRISPR/Cas9 system, dramatically crippled cell proliferation and induced cell senescence. In conclusion, identifying transcriptionally-active HPV integrants with Nanopore sequencing can provide viable targets for gene therapy against HPV-associated cancers.

Accurate and comprehensive evaluation of O6-methylguanine-DNA methyltransferase promoter methylation by nanopore sequencing.

AIMS: The methylation status of the O6-methylguanine-DNA methyltransferase (MGMT) promoter region is essential in evaluating the prognosis and predicting the drug response in patients with glioblastoma. In this study, we evaluated the utility of using nanopore long-read sequencing as a method for assessing methylation levels throughout the MGMT CpG-island, compared its performance to established techniques and demonstrated its clinical applicability. METHODS: We analysed 165 samples from CNS tumours, focusing on the MGMT CpG-island using nanopore sequencing. Oxford Nanopore Technologies (ONT) MinION and PromethION flow cells were employed for single sample or barcoded assays, guided by a CRISPR/Cas9 protocol, adaptive sampling or as part of a whole genome sequencing assay. Methylation data obtained through nanopore sequencing were compared to results obtained via pyrosequencing and methylation bead arrays. Hierarchical clustering was applied to nanopore sequencing data for patient stratification. RESULTS: Nanopore sequencing displayed a strong correlation (R2 = 0.91) with pyrosequencing results for the four CpGs of MGMT analysed by both methods. The MGMT-STP27 algorithm's classification was effectively reproduced using nanopore data. Unsupervised hierarchical clustering revealed distinct patterns in methylated and unmethylated samples, providing comparable survival prediction capabilities. Nanopore sequencing yielded high-confidence results in a rapid timeframe, typically within hours of sequencing, and extended the analysis to all 98 CpGs of the MGMT CpG-island. CONCLUSIONS: This study presents nanopore sequencing as a valid and efficient method for determining MGMT promotor methylation status. It offers a comprehensive view of the MGMT promoter methylation landscape, which enables the identification of potentially clinically relevant subgroups of patients. Further exploration of the clinical implications of patient stratification using nanopore sequencing of MGMT is warranted.

Microbiome diversity and zoonotic bacterial pathogen prevalence in Peromyscus mice from agricultural landscapes and synanthropic habitat.

Rodents are key reservoirs of zoonotic pathogens and play an important role in disease transmission to humans. Importantly, anthropogenic land-use change has been found to increase the abundance of rodents that thrive in human-built environments (synanthropic rodents), particularly rodent reservoirs of zoonotic disease. Anthropogenic environments also affect the microbiome of synanthropic wildlife, influencing wildlife health and potentially introducing novel pathogens. Our objective was to examine the effect of agricultural development and synanthropic habitat on microbiome diversity and the prevalence of zoonotic bacterial pathogens in wild Peromyscus mice to better understand the role of these rodents in pathogen maintenance and transmission. We conducted 16S amplicon sequencing on faecal samples using long-read nanopore sequencing technology to characterize the rodent microbiome. We compared microbiome diversity and composition between forest and synanthropic habitats in agricultural and undeveloped landscapes and screened for putative pathogenic bacteria. Microbiome richness, diversity, and evenness were higher in the agricultural landscape and synanthropic habitat compared to undeveloped-forest habitat. Microbiome composition also differed significantly between agricultural and undeveloped landscapes and forest and synanthropic habitats. We detected overall low diversity and abundance of putative pathogenic bacteria, though putative pathogens were more likely to be found in mice from the agricultural landscape. Our findings show that landscape- and habitat-level anthropogenic factors affect Peromyscus microbiomes and suggest that landscape-level agricultural development may be important to predict zoonotic pathogen prevalence. Ultimately, understanding how anthropogenic land-use change and synanthropy affect rodent microbiomes and pathogen prevalence is important to managing transmission of rodent-borne zoonotic diseases to humans.

Identification and expression analysis of Phosphate Transporter 1 (PHT1) genes in the highly phosphorus-use-efficient Hakea prostrata (Proteaceae).

Heavy and costly use of phosphorus (P) fertiliser is often needed to achieve high crop yields, but only a small amount of applied P fertiliser is available to most crop plants. Hakea prostrata (Proteaceae) is endemic to the P-impoverished landscape of southwest Australia and has several P-saving traits. We identified 16 members of the Phosphate Transporter 1 (PHT1) gene family (HpPHT1;1-HpPHT1;12d) in a long-read genome assembly of H. prostrata. Based on phylogenetics, sequence structure and expression patterns, we classified HpPHT1;1 as potentially involved in Pi uptake from soil and HpPHT1;8 and HpPHT1;9 as potentially involved in Pi uptake and root-to-shoot translocation. Three genes, HpPHT1;4, HpPHT1;6 and HpPHT1;8, lacked regulatory PHR1-binding sites (P1BS) in the promoter regions. Available expression data for HpPHT1;6 and HpPHT1;8 indicated they are not responsive to changes in P supply, potentially contributing to the high P sensitivity of H. prostrata. We also discovered a Proteaceae-specific clade of closely-spaced PHT1 genes that lacked conserved genetic architecture among genera, indicating an evolutionary hot spot within the genome. Overall, the genome assembly of H. prostrata provides a much-needed foundation for understanding the genetic mechanisms of novel adaptations to low P soils in southwest Australian plants.