Microscopic examination of spatial transcriptome using Seq-Scope.

Spatial barcoding technologies have the potential to reveal histological details of transcriptomic profiles; however, they are currently limited by their low resolution. Here, we report Seq-Scope, a spatial barcoding technology with a resolution comparable to an optical microscope. Seq-Scope is based on a solid-phase amplification of randomly barcoded single-molecule oligonucleotides using an Illumina sequencing platform. The resulting clusters annotated with spatial coordinates are processed to expose RNA-capture moiety. These RNA-capturing barcoded clusters define the pixels of Seq-Scope that are ∼0.5-0.8 µm apart from each other. From tissue sections, Seq-Scope visualizes spatial transcriptome heterogeneity at multiple histological scales, including tissue zonation according to the portal-central (liver), crypt-surface (colon) and inflammation-fibrosis (injured liver) axes, cellular components including single-cell types and subtypes, and subcellular architectures of nucleus and cytoplasm. Seq-Scope is quick, straightforward, precise, and easy-to-implement and makes spatial single-cell analysis accessible to a wide group of biomedical researchers.

Efficient Linear dsDNA Tagging Using Deoxyuridine Excision*.

Site-specific strategies for exchanging segments of dsDNA are important for DNA library construction and molecular tagging. Deoxyuridine (dU) excision is an approach for generating 3' ssDNA overhangs in gene assembly and molecular cloning procedures. Unlike approaches that use a multi-base pair motif to specify a DNA cut site, dU excision requires only a dT→dU substitution. Consequently, excision sites can be embedded in biologically active DNA sequences by placing dU substitutions at non-perturbative positions. In this work, I describe a molecular tagging method that uses dU excision to exchange a segment of a dsDNA strand with a long synthetic oligonucleotide. The core workflow of this method, called deoxyUridine eXcision-tagging (dUX-tagging), is an efficient one-pot reaction: strategically positioned dU nucleotides are excised from dsDNA to generate a 3' overhang so that additional sequence can be appended by annealing and ligating a tagging oligonucleotide. The tagged DNA is then processed by one of two procedures to fill the 5' overhang and remove excess tagging oligo. To facilitate its widespread use, all dUX-tagging procedures exclusively use commercially available reagents. As a result, dUX-tagging is a concise and easily implemented approach for high-efficiency linear dsDNA tagging.

Seewis hantavirus in common shrew (Sorex araneus) in Sweden.

BACKGROUND: Rodent borne hantaviruses are emerging viruses infecting humans through inhalation. They cause hemorrhagic fever with renal syndrome and hemorrhagic cardiopulmonary syndrome. Recently, hantaviruses have been detected in other small mammals such as Soricomorpha (shrews, moles) and Chiroptera (bats), suggested as reservoirs for potential pandemic viruses and to play a role in the evolution of hantaviruses. It is important to study the global virome in different reservoirs, therefore our aim was to investigate whether shrews in Sweden carried any hantaviruses. Moreover, to accurately determine the host species, we developed a molecular method for identification of shrews. METHOD: Shrews (n = 198), caught during 1998 in Sweden, were screened with a pan-hantavirus PCR using primers from a conserved region of the large genome segment. In addition to morphological typing of shrews, we developed a molecular based typing method using sequencing of the mitochondrial cytochrome C oxidase I (COI) and cytochrome B (CytB) genes. PCR amplified hantavirus and shrew fragments were sequenced and phylogenetically analysed. RESULTS: Hantavirus RNA was detected in three shrews. Sequencing identified the virus as Seewis hantavirus (SWSV), most closely related to previous isolates from Finland and Russia. All three SWSV sequences were retrieved from common shrews (Sorex araneus) sampled in Västerbotten County, Sweden. The genetic assay for shrew identification was able to identify native Swedish shrew species, and the genetic typing of the Swedish common shrews revealed that they were most similar to common shrews from Russia. CONCLUSION: We detected SWSV RNA in Swedish common shrew samples and developed a genetic assay for shrew identification based on the COI and CytB genes. This was the first report of presence of hantavirus in Swedish shrews.

Application of amplicon-based targeted sequencing with the molecular barcoding system to detect uncommon minor EGFR mutations in patients with treatment-naïve lung adenocarcinoma.

BACKGROUND: In lung cancer, epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor sensitizing mutations co-existing with rare minor EGFR mutations are known as compound mutations. These minor EGFR mutations can lead to acquired resistance after EGFR tyrosine kinase inhibitor treatment, so determining the mutation status of patients is important. However, using amplicon-based targeted deep sequencing based on next-generation sequencing to characterize mutations is prone to sequencing error. We therefore assessed the benefit of incorporating molecular barcoding with high-throughput sequencing to investigate genomic heterogeneity in treatment-naïve patients who have undergone resection of their non-small cell lung cancer (NSCLC) EGFR mutations. METHODS: We performed amplicon-based targeted sequencing with the molecular barcoding system (MBS) to detect major common EGFR mutations and uncommon minor mutations at a 0.5% allele frequency in fresh-frozen lung cancer samples. RESULTS: Profiles of the common mutations of EGFR identified by MBS corresponded with the results of clinical testing in 63 (98.4%) out of 64 cases. Uncommon mutations of EGFR were detected in seven cases (10.9%). Among the three types of major EGFR mutations, patients with the G719X mutation had a significantly higher incidence of compound mutations than those with the L858R mutation or exon 19 deletion (p = 0.0052). This was validated in an independent cohort from the Cancer Genome Atlas dataset (p = 0.018). CONCLUSIONS: Our findings demonstrate the feasibility of using the MBS to establish an accurate NSCLC patient genotype. This work will help understand the molecular basis of EGFR compound mutations in NSCLC, and could aid the development of new treatment modalities.

Recent Advances in Liquid Biopsy in Precision Oncology Research.

Liquid biopsy is a minimally invasive test for cancer genetic status based on circulating tumor DNA (ctDNA), circulating tumor cells, or other tumor-derived materials in blood plasma. Although the minimal invasiveness and time resolution are attractive features of liquid biopsy, the limited amount of ctDNA in plasma poses problems. Recent developments in digital PCR and next-generation sequencing (NGS)-based technology have improved the accuracy of liquid biopsy. In particular, molecular barcoding technology in NGS-based methods, i.e., tagging of molecular barcodes to cell-free DNA before amplification, reduces technical errors by validating the consensus of sequences originating from a single molecule, leading to marked improvement of the accuracy and detection limit. However, substitutions caused by DNA damage and somatic mutations originating from normal cells are still obstacles to the sensitive detection of mutations on ctDNA. Since there have been only a few clinical applications, a deeper understanding of ctDNA biology and more advanced analytical technology are needed for the practical application of liquid biopsy.

Nx1-Seq (Well Based Single-Cell Analysis System).

Research on the hierarchical nature of cell differentiation and heterogeneity in tissues has been performed by isolating and identifying cells by the use of monoclonal antibodies, cell sorting, microdissection, and functional assays. However, it is difficult to analyze continuous changes in cell differentiation and the identification of cells for which cell markers are unclear. Furthermore, cell populations considered identical were shown to be diverse. Recently, single cell gene expression analysis was performed to help understand the complexity of cell populations. Single-cell analysis can analyze the diversity of individual cell populations as well as the tissue microenvironment, and is extremely useful for research on intercellular interactions in diseases and identifying specific marker genes. Recent advances in technology have made it possible to analyze hundreds of single cells. In this paper, we introduce our newly developed well-based single-cell transcriptome method, which includes other methods.

Quantitation of mRNA Transcripts and Proteins Using the BD Rhapsody™ Single-Cell Analysis System.

In this review, we describe the BD Rhapsody™ Single-Cell Analysis System, a platform that allows high-throughput capture of nucleic acids from single cells using a simple cartridge workflow and a multitier barcoding system. The resulting captured information can be used to generate various types of next-generation sequencing (NGS) libraries, including whole transcriptome analysis for discovery biology and targeted RNA analysis for high sensitivity transcript detection. The BD Rhapsody system can be used with emerging applications, such as BD™ AbSeq assays, to profile gene expression in both mRNA and protein level to provide ultra-high resolution analysis of single cells.

Single molecule counting and assessment of random molecular tagging errors with transposable giga-scale error-correcting barcodes.

BACKGROUND: RNA-Seq measures gene expression by counting sequence reads belonging to unique cDNA fragments. Molecular barcodes commonly in the form of random nucleotides were recently introduced to improve gene expression measures by detecting amplification duplicates, but are susceptible to errors generated during PCR and sequencing. This results in false positive counts, leading to inaccurate transcriptome quantification especially at low input and single-cell RNA amounts where the total number of molecules present is minuscule. To address this issue, we demonstrated the systematic identification of molecular species using transposable error-correcting barcodes that are exponentially expanded to tens of billions of unique labels. RESULTS: We experimentally showed random-mer molecular barcodes suffer from substantial and persistent errors that are difficult to resolve. To assess our method's performance, we applied it to the analysis of known reference RNA standards. By including an inline random-mer molecular barcode, we systematically characterized the presence of sequence errors in random-mer molecular barcodes. We observed that such errors are extensive and become more dominant at low input amounts. CONCLUSIONS: We described the first study to use transposable molecular barcodes and its use for studying random-mer molecular barcode errors. Extensive errors found in random-mer molecular barcodes may warrant the use of error correcting barcodes for transcriptome analysis as input amounts decrease.

DNA Barcoding Coupled with High Resolution Melting Analysis Enables Rapid and Accurate Distinction of Aspergillus species.

We describe a high-resolution melting (HRM) analysis method that is rapid, reproducible, and able to identify reference strains and further 40 clinical isolates of Aspergillus fumigatus (14), A. lentulus (3), A. terreus (7), A. flavus (8), A. niger (2), A. welwitschiae (4), and A. tubingensis (2). Asp1 and Asp2 primer sets were designed to amplify partial sequences of the Aspergillus benA (beta-tubulin) genes in a closed-, single-tube system. Human placenta DNA, further Aspergillus (3), Candida (9), Fusarium (6), and Scedosporium (2) nucleic acids from type strains and clinical isolates were also included in this study to evaluate cross reactivity with other relevant pathogens causing invasive fungal infections. The barcoding capacity of this method proved to be 100% providing distinctive binomial scores; 14, 34, 36, 35, 25, 15, 26 when tested among species, while the within-species distinction capacity of the assay proved to be 0% based on the aligned thermodynamic profiles of the Asp1, Asp2 melting clusters allowing accurate species delimitation of all tested clinical isolates. The identification limit of this HRM assay was also estimated on Aspergillus reference gDNA panels where it proved to be 10-102 genomic equivalents (GE) except the A. fumigatus panel where it was 103 only. Furthermore, misidentification was not detected with human genomic DNA or with Candida, Fusarium, and Scedosporium strains. Our DNA barcoding assay introduced here provides results within a few hours, and it may possess further diagnostic utility when analyzing standard cultures supporting adequate therapeutic decisions.

Molecular indexing enables quantitative targeted RNA sequencing and reveals poor efficiencies in standard library preparations.

We present a simple molecular indexing method for quantitative targeted RNA sequencing, in which mRNAs of interest are selectively captured from complex cDNA libraries and sequenced to determine their absolute concentrations. cDNA fragments are individually labeled so that each molecule can be tracked from the original sample through the library preparation and sequencing process. Multiple copies of cDNA fragments of identical sequence become distinct through labeling, and replicate clones created during PCR amplification steps can be identified and assigned to their distinct parent molecules. Selective capture enables efficient use of sequencing for deep sampling and for the absolute quantitation of rare or transient transcripts that would otherwise escape detection by standard sequencing methods. We have also constructed a set of synthetic barcoded RNA molecules, which can be introduced as controls into the sample preparation mix and used to monitor the efficiency of library construction. The quantitative targeted sequencing revealed extremely low efficiency in standard library preparations, which were further confirmed by using synthetic barcoded RNA molecules. This finding shows that standard library preparation methods result in the loss of rare transcripts and highlights the need for monitoring library efficiency and for developing more efficient sample preparation methods.

Arctic marine phytobenthos of northern Baffin Island.

Global climate change is expected to alter the polar bioregions faster than any other marine environment. This study assesses the biodiversity of seaweeds and associated eukaryotic pathogens of an established study site in northern Baffin Island (72° N), providing a baseline inventory for future work assessing impacts of the currently ongoing changes in the Arctic marine environment. A total of 33 Phaeophyceae, 24 Rhodophyceae, 2 Chlorophyceae, 12 Ulvophyceae, 1 Trebouxiophyceae, and 1 Dinophyceae are reported, based on collections of an expedition to the area in 2009, complemented by unpublished records of Robert T. Wilce and the first-ever photographic documentation of the phytobenthos of the American Arctic. Molecular barcoding of isolates raised from incubated substratum samples revealed the presence of 20 species of brown seaweeds, including gametophytes of kelp and of a previously unsequenced Desmarestia closely related to D. viridis, two species of Pylaiella, the kelp endophyte Laminariocolax aecidioides and 11 previously unsequenced species of the Ectocarpales, highlighting the necessity to include molecular techniques for fully unraveling cryptic algal diversity. This study also includes the first records of Eurychasma dicksonii, a eukaryotic pathogen affecting seaweeds, from the American Arctic. Overall, this study provides both the most accurate inventory of seaweed diversity of the northern Baffin Island region to date and can be used as an important basis to understand diversity changes with climate change.

Pollination of Vietnamese Aspidistra xuansonensis (Asparagaceae) by female Cecidomyiidi flies: larvae of pollinator feed on fertile pollen in anthers of anthetic bisexual flowers.

UNLABELLED: • PREMISE OF THE STUDY: Aspidistra is a species-rich, herbaceous monocot genus of tropical Southeast Asia. Most species are recently discovered and apparently endangered, though virtually nothing is known about their biology. Species of the genus are primarily distinguished using flower morphology, which is enormously diverse. However, the pollination process has not been directly observed in the center of diversity of the genus (N Vietnam and S China). Indirect and partly direct data on the only widely cultivated species of the genus (A. elatior) placed it among angiosperms with the most unusual pollination biology, though these data are highly controversial, suggesting pollen transfer by mollusks, crustaceans, flies, or possibly tiny soil invertebrates such as collembolans.• METHODS: Pollination of Aspidistra xuansonensis in the center of diversity of the genus was studied using visual observations and videos and light and scanning electron microscopy investigation of flowers and their pollinators. Pollinators and their larvae were molecularly barcoded.• KEY RESULTS: Aspidistra xuansonensis is pollinated by female cecidomyiid flies (gall midges). They oviposit on anthers, and larvae develop among the pollen mass. Molecular barcoding proved taxonomic identity of the larvae and the flies. The larvae neither damage floral parts nor cause gall formation, but feed on pollen grains by sucking out their content. The larvae move out of the flowers before decomposition starts. Carebara ants steal developing larvae from flowers but do not contribute to pollination.• CONCLUSIONS: More than one kind of myiophily is present in Aspidistra. Brood site pollination was documented for the first time in Aspidistra. The pollination system of A. xuansonensis differs from other kinds of brood site pollination in the exit of the larvae prior to the decomposition of floral parts.

Multi-Locus Sequence Analysis Indicates Potential Cryptic Speciation in the Chigger Mite Neoschoengastia gallinarum (Hatori, 1920) Parasitising Birds in Asia.

Neoschoengastia gallinarum is widely distributed in Asia, preferentially parasitising birds, and heavy infestations have clinical impacts on domestic fowl. In common with other trombiculid mites, the genetic diversity and potential variation in host preferences or pathology induced by N. gallinarum are poorly understood. This study aimed to unravel the geographical variation and population structure of N. gallinarum collected from galliform birds in Peninsular Malaysia and Thailand by inference from concatenated mitochondrial-encoded cytochrome c oxidase subunit I (COI), and nuclear-encoded internal transcribed spacer 2 (ITS2) and 18S ribosomal DNA gene sequences, including a comparison with previously published data from southeastern China. Our multi-locus sequence analysis revealed three monophyletic clades comprising (A) specimens from Peninsular Malaysia, (B) the samples from Thailand together with a minority of Chinese sequences, and (C) the majority of sequences from China. Similarly, most species delimitation approaches divided the specimens into three operational taxonomic units. Analysis of molecular variance revealed 96.41% genetic divergence between Malaysian and Thai populations, further supported by the absence of gene flow (Nm = 0.01). In conclusion, despite the two countries sharing a land border, populations of N. gallinarum from Peninsular Malaysia and Thailand appear to be genetically segregated and may represent distinct cryptic species.

Additional new species and new records of the genus Sticta (lichenised Ascomycota, lobarioid Peltigeraceae) from Bolivia.

Four species of the genus Sticta are described as new from Bolivia, based on morphological examination and phylogenetic analysis of the fungal ITS barcoding marker. Additionally, two species are reported as new to Bolivia (their identification confirmed by molecular data) and one previously reported species is confirmed by molecular data for the first time. Detailed morphological and anatomical descriptions are provided for all new species. Two of the new species, S.isidiolobulata Ossowska, B. Moncada, Lücking & Kukwa and S.madidiensis Ossowska, B. Moncada, Lücking & Kukwa belong to clade I, as defined in previous studies. In contrast, S.montepunkuensis Ossowska, B. Moncada, Lücking & Kukwa and S.macrolobata Ossowska, B. Moncada, Lücking & Kukwa, also described here as new to science, belong to clade III. Stictaisidiolobulata has an irregular to suborbicular thallus of medium size, with isidia developing into spathulate lobules, cyanobacterial photobiont and apothecia with entire to weakly-crenate margins. The large irregular thallus of the cyanobacteria-associated S.macrolobata has broad lobes, apothecia with verrucous to tomentose margins and cyphellae with raised margins, whereas S.madidiensis has a medium-sized, palmate to irregular thallus with a stipe, but without vegetative propagules and apothecia. Stictamontepunkuensis has large and irregular thalli with green algae as photobiont, apothecia with crenate to verrucous margins and urceolate cyphellae with a wide pore and a scabrid basal membrane. Two species, S.beauvoisii Delise and S.riparia Merc.-Díaz are reported as new to Bolivia (the latter also as new to South America) and belong to clade III. Stictatomentosa (Sw.) Ach., species confirmed from Bolivia by molecular data, belongs to clade II. Stictabeauvoisii is characterised by a smooth yellowish-brown upper surface with darker apices and abundant, marginal isidia and a brown lower surface with golden-chocolate brown primary tomentum and sparse, golden-brown rhizines. Stictariparia has a strongly branched thallus, with undulate lobes and abundant, marginal, palmate, grey to dark brown phyllidia and greyish-brown lower surface with the primary tomentum absent towards the margins. Stictatomentosa has palmate, bluish thalli with white cilia and abundant, submarginal apothecia and creamy-white lower surface with a sparse, white primary tomentum.

Single-Cell Molecular Barcoding to Decode Multimodal Information Defining Cell States.

Single-cell research has provided a breakthrough in biology to understand heterogeneous cell groups, such as tissues and organs, in development and disease. Molecular barcoding and subsequent sequencing technology insert a singlecell barcode into isolated single cells, allowing separation cell by cell. Given that multimodal information from a cell defines precise cellular states, recent technical advances in methods focus on simultaneously extracting multimodal data recorded in different biological materials (DNA, RNA, protein, etc.). This review summarizes recently developed singlecell multiomics approaches regarding genome, epigenome, and protein profiles with the transcriptome. In particular, we focus on how to anchor or tag molecules from a cell, improve throughputs with sample multiplexing, and record lineages, and we further discuss the future developments of the technology.

Strategies for improving detection of circulating tumor DNA using next generation sequencing.

Cancer has become a global health issue and liquid biopsy has emerged as a non-invasive tool for various applications. In cancer, circulating tumor DNA (ctDNA) can be detected from cell-free DNA (cfDNA) obtained from plasma and has potential for early diagnosis, treatment, resistance, minimal residual disease detection, and tumoral heterogeneity identification. However, the low frequency of ctDNA requires techniques for accurate analysis. Multitarget assay such as Next Generation Sequencing (NGS) need improvement to achieve limits of detection that can identify the low frequency variants present in the cfDNA. In this review, we provide a general overview of the use of cfDNA and ctDNA in cancer, and discuss techniques developed to optimize NGS as a tool for ctDNA detection. We also summarize the results obtained using NGS strategies in both investigational and clinical contexts.

Molecular Barcoding: A Tool to Guarantee Correct Seafood Labelling and Quality and Preserve the Conservation of Endangered Species.

The recent increase in international fish trade leads to the need for improving the traceability of fishery products. In relation to this, consistent monitoring of the production chain focusing on technological developments, handling, processing and distribution via global networks is necessary. Molecular barcoding has therefore been suggested as the gold standard in seafood species traceability and labelling. This review describes the DNA barcoding methodology for preventing food fraud and adulteration in fish. In particular, attention has been focused on the application of molecular techniques to determine the identity and authenticity of fish products, to discriminate the presence of different species in processed seafood and to characterize raw materials undergoing food industry processes. In this regard, we herein present a large number of studies performed in different countries, showing the most reliable DNA barcodes for species identification based on both mitochondrial (COI, cytb, 16S rDNA and 12S rDNA) and nuclear genes. Results are discussed considering the advantages and disadvantages of the different techniques in relation to different scientific issues. Special regard has been dedicated to a dual approach referring to both the consumer's health and the conservation of threatened species, with a special focus on the feasibility of the different genetic and genomic approaches in relation to both scientific objectives and permissible costs to obtain reliable traceability.

Nanopore adaptive sampling for targeted mitochondrial genome sequencing and bloodmeal identification in hematophagous insects.

BACKGROUND: Blood-feeding insects are important vectors for an array of zoonotic pathogens. While previous efforts toward generating molecular resources have largely focused on major vectors of global medical and veterinary importance, molecular data across a large number of hematophagous insect taxa remain limited. Advancements in long-read sequencing technologies and associated bioinformatic pipelines provide new opportunities for targeted sequencing of insect mitochondrial (mt) genomes. For engorged hematophagous insects, such technologies can be leveraged for both insect mitogenome genome assembly and identification of vertebrate blood-meal sources. METHODS: We used nanopore adaptive sampling (NAS) to sequence genomic DNA from four species of field-collected, blood-engorged mosquitoes (Aedes and Culex spp.) and one deer fly (Chrysops sp.). NAS was used for bioinformatical enrichment of mtDNA reads of hematophagous insects and potential vertebrate blood-meal hosts using publically available mt genomes as references. We also performed an experimental control to compare results of traditional non-NAS nanopore sequencing to the mt genome enrichment by the NAS method. RESULTS: Complete mitogenomes were assembled and annotated for all five species sequenced with NAS: Aedes trivittatus, Aedes vexans, Culex restuans, Culex territans and the deer fly, Chrysops niger. In comparison to data generated during our non-NAS control experiment, NAS yielded a substantially higher proportion of reference-mapped mtDNA reads, greatly streamlining downstream mitogenome assembly and annotation. The NAS-assembled mitogenomes ranged in length from 15,582 to 16,045 bp, contained between 78.1% and 79.0% A + T content and shared the anticipated arrangement of 13 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs. Maximum likelihood phylogenies were generated to further characterize each insect species. Additionally, vertebrate blood-meal analysis was successful in three samples sequenced, with mtDNA-based phylogenetic analyses revealing that blood-meal sources for Chrysops niger, Culex restuans and Aedes trivittatus were human, house sparrow (Passer domesticus) and eastern cottontail rabbit (Sylvilagus floridanus), respectively. CONCLUSIONS: Our findings show that NAS has dual utility to simultaneously molecularly identify hematophagous insects and their blood-meal hosts. Moreover, our data indicate NAS can facilitate a wide array of mitogenomic systematic studies through novel 'phylogenetic capture' methods. We conclude that the NAS approach has great potential for broadly improving genomic resources used to identify blood-feeding insects, answer phylogenetic questions and elucidate complex pathways for the transmission of vector-borne pathogens.

Genetic diversity and organic waste degrading capacity of Hermetia illucens from the evergreen forest of the Equatorial Choco lowland.

Globally, microplastics (MP) represent a growing burden for ecosystems due to their increasing presence at different trophic levels. In Ecuador, the lack of waste segregation has increased the quantity of waste, primarily organics and plastics, overloading landfills and water sources. Over time, plastics reduce in size and silently enter the food chain of animals, such as insects. The black soldier fly (BSF) larvae, Hermetia illucens (Linnaeus, 1758), is a species with devouring behavior used for waste management because of its beneficial qualities such as fly pest control, biomass production, and rapid organic waste degradation. Studies have uncovered the insect's ability to tolerate MP, and consider the possibility that they may be able to degrade polymers. For the first time in Ecuador, the present study characterized H. illucens using the sequences of different molecular markers. Finally, H. illucens' degrading capacity was evaluated in the presence of MP and decaying food residues, resembling landfill conditions.

Response of planktonic microbial assemblages to disturbance in an urban sub-tropical estuary.

Microbes are sensitive indicators of estuarine processes because they respond rapidly to dynamic disturbance events. As most of the world's population lives in urban areas and climate change-related disturbance events are becoming more frequent, estuaries bounded by cities are experiencing increasing stressors, at the same time that their ecosystem services are required more than ever. Here, using a multidisciplinary approach, we determined the response of planktonic microbial assemblages in response to seasonality and a rainfall disturbance in an urban estuary bounded by Australia's largest city, Sydney. We used molecular barcoding (16S, 18S V4 rRNA) and microscopy-based identification to compare microbial assemblages at locations with differing characteristics and urbanisation histories. Across 142 samples, we identified 8,496 unique free-living bacterial zOTUs, 8,175 unique particle associated bacterial zOTUs, and 1,920 unique microbial eukaryotic zOTUs. Using microscopy, we identified only the top <10% abundant, larger eukaryotic taxa (>10 µm), however quantification was possible. The site with the greater history of anthropogenic impact showed a more even community of associated bacteria and eukaryotes, and a significant increase in dissolved inorganic nitrogen following rainfall, when compared to the more buffered site. This coincided with a reduced proportional abundance of Actinomarina and Synechococcus spp., a change in SAR 11 clades, and an increase in the eukaryotic microbial groups Dinophyceae, Mediophyceae and Bathyoccocaceae, including a temporary dominance of the harmful algal bloom dinoflagellate Prorocentrum cordatum (syn. P. minimum). Finally, a validated hydrodynamic model of the estuary supported these results, showing that the more highly urbanised and upstream location consistently experienced a higher magnitude of salinity reduction in response to rainfall events during the study period. The best abiotic variables to explain community dissimilarities between locations were TDP, PN, modelled temperature and salinity (r = 0.73) for the free living bacteria, TP for the associated bacteria (r = 0.43), and modelled temperature (r = 0.28) for the microbial eukaryotic communities. Overall, these results show that a minor disturbance such as a brief rainfall event can significantly shift the microbial assemblage of an anthropogenically impacted area within an urban estuary to a greater degree than a seasonal change, but may result in a lesser response to the same disturbance at a buffered, more oceanic influenced location. Fine scale research into the factors driving the response of microbial communities in urban estuaries to climate related disturbances will be necessary to understand and implement changes to maintain future estuarine ecosystem services.