ClickSeq: Replacing Fragmentation and Enzymatic Ligation with Click-Chemistry to Prevent Sequence Chimeras

SARS-CoV-2 EndoU-ribonuclease regulates RNA recombination and impacts viral fitness

Coronaviruses (CoVs) maintain large RNA genomes that frequently undergoes mutations and recombination, contributing to their evolution and emergence. In this study, we find that SARS-CoV-2 has greater RNA recombination frequency than other human CoVs. In addition, coronavirus RNA recombination primarily occurs at uridine (U)-enriched RNA sequences. Therefore, we next evaluated the role of SARS-CoV-2 NSP15, a viral endonuclease that targets uridines (EndoU), in RNA recombination and virus infection. Using a catalytically inactivated EndoU mutant (NSP15H234A), we observed attenuated viral replication in vitro and in vivo. However, the loss of EndoU activity also dysregulated inflammation resulting in similar disease in vivo despite reduced viral loads. Next-generation sequencing (NGS) demonstrated that loss of EndoU activity disrupts SARS-CoV-2 RNA recombination by reducing viral sub-genomic message but increasing recombination events that contribute to defective viral genomes (DVGs). Overall, the study demonstrates that NSP15 plays a critical role in regulating RNA recombination and SARS-CoV-2 pathogenesis.

RpoS activates formation of Salmonella Typhi biofilms and drives persistence in the gall bladder

The development of strategies for targeting the asymptomatic carriage of Salmonella Typhi in chronic typhoid patients has suffered owing to our basic lack of understanding of the molecular mechanisms that enable the formation of S. Typhi biofilms. Traditionally, studies have relied on cholesterol-attached biofilms formed by a closely related serovar, Typhimurium, to mimic multicellular Typhi communities formed on human gallstones. In long-term infections, S. Typhi adopts the biofilm lifestyle to persist in vivo and survive in the carrier state, ultimately leading to the spread of infections via the fecal-oral route of transmission. In the present work, we studied S. Typhi biofilms directly, applied targeted as well as genome-wide genetic approaches to uncover unique biofilm components that do not conform to the CsgD-dependent pathway established in S. Typhimurium. We undertook a genome-wide Tn5 mutation screen in H58, a clinically relevant multidrug resistance strain of S. Typhi, in gallstone-mimicking conditions. We generated New Generation Sequencing libraries based on the ClickSeq technology to identify the key regulators, IraP and RpoS, and the matrix components Sth fimbriae, Vi capsule and lipopolysaccharide. We discovered that the starvation sigma factor, RpoS, was required for the transcriptional activation of matrix-encoding genes in vitro, and for S. Typhi colonization in persistent infections in vivo, using a heterologous fish larval model. An rpoS null mutant failed to colonize the gall bladder in chronic zebrafish infections. Overall, our work uncovered a novel RpoS-driven, CsgD-independent paradigm for the formation of cholesterol-attached Typhi biofilms, and emphasized the role(s) of stress signaling pathways for adaptation in chronic infections. Our identification of the biofilm regulators in S. Typhi paves the way for the development of drugs against typhoid carriage, which will ultimately control the increased incidence of gall bladder cancer in typhoid carriers.

Cytoplasmic binding partners of the Integrator endonuclease INTS11 and its paralog CPSF73 are required for their nuclear function

INTS11 and CPSF73 are metal-dependent endonucleases for Integrator and pre-mRNA 3'-end processing, respectively. Here, we show that the INTS11 binding partner BRAT1/CG7044, a factor important for neuronal fitness, stabilizes INTS11 in the cytoplasm and is required for Integrator function in the nucleus. Loss of BRAT1 in neural organoids leads to transcriptomic disruption and precocious expression of neurogenesis-driving transcription factors. The structures of the human INTS9-INTS11-BRAT1 and Drosophila dIntS11-CG7044 complexes reveal that the conserved C terminus of BRAT1/CG7044 is captured in the active site of INTS11, with a cysteine residue directly coordinating the metal ions. Inspired by these observations, we find that UBE3D is a binding partner for CPSF73, and UBE3D likely also uses a conserved cysteine residue to directly coordinate the active site metal ions. Our studies have revealed binding partners for INTS11 and CPSF73 that behave like cytoplasmic chaperones with a conserved impact on the nuclear functions of these enzymes.

Optimized SMRT-UMI protocol produces highly accurate sequence datasets from diverse populations-Application to HIV-1 quasispecies

Pathogen diversity resulting in quasispecies can enable persistence and adaptation to host defenses and therapies. However, accurate quasispecies characterization can be impeded by errors introduced during sample handling and sequencing, which can require extensive optimizations to overcome. We present complete laboratory and bioinformatics workflows to overcome many of these hurdles. The Pacific Biosciences single molecule real-time platform was used to sequence polymerase-chain reaction (PCR) amplicons derived from cDNA templates tagged with unique molecular identifiers (SMRT-UMI). Optimized laboratory protocols were developed through extensive testing of different sample preparation conditions to minimize between-template recombination during PCR. The use of UMI allowed accurate template quantitation as well as removal of point mutations introduced during PCR and sequencing to produce a highly accurate consensus sequence from each template. Production of highly accurate sequences from the large datasets produced from SMRT-UMI sequencing is facilitated by a novel bioinformatic pipeline, Probabilistic Offspring Resolver for Primer IDs (PORPIDpipeline). PORPIDpipeline automatically filters and parses circular consensus reads by sample, identifies and discards reads with UMIs likely created from PCR and sequencing errors, generates consensus sequences, checks for contamination within the dataset, and removes any sequence with evidence of PCR recombination, heteroduplex formation, or early cycle PCR errors. The optimized SMRT-UMI sequencing and PORPIDpipeline methods presented here represent a highly adaptable and established starting point for accurate sequencing of diverse pathogens. These methods are illustrated through characterization of human immunodeficiency virus quasispecies in a virus transmitter-recipient pair of individuals.

A multidisciplinary approach to the identification of the protein-RNA connectome in double-stranded RNA virus capsids

How multi-segmented double-stranded RNA (dsRNA) viruses correctly incorporate their genomes into their capsids remains unclear for many viruses, including Bluetongue virus (BTV), a Reoviridae member, with a genome of 10 segments. To address this, we used an RNA-cross-linking and peptide-fingerprinting assay (RCAP) to identify RNA binding sites of the inner capsid protein VP3, the viral polymerase VP1 and the capping enzyme VP4. Using a combination of mutagenesis, reverse genetics, recombinant proteins and in vitro assembly, we validated the importance of these regions in virus infectivity. Further, to identify which RNA segments and sequences interact with these proteins, we used viral photo-activatable ribonucleoside crosslinking (vPAR-CL) which revealed that the larger RNA segments (S1-S4) and the smallest segment (S10) have more interactions with viral proteins than the other smaller segments. Additionally, using a sequence enrichment analysis we identified an RNA motif of nine bases that is shared by the larger segments. The importance of this motif for virus replication was confirmed by mutagenesis followed by virus recovery. We further demonstrated that these approaches could be applied to a related Reoviridae member, rotavirus (RV), which has human epidemic impact, offering the possibility of novel intervention strategies for a human pathogen.

Improved Nanopore full-length cDNA sequencing by PCR-suppression

Full-length transcript sequencing remains a main goal of RNA sequencing. However, even the application of long-read sequencing technologies such as Oxford Nanopore Technologies still fail to yield full-length transcript sequencing for a significant portion of sequenced reads. Since these technologies can sequence reads that are far longer than the longest known processed transcripts, the lack of efficiency to obtain full-length transcripts from good quality RNAs stems from library preparation inefficiency rather than the presence of degraded RNA molecules. It has previously been shown that addition of inverted terminal repeats in cDNA during reverse transcription followed by single-primer PCR creates a PCR suppression effect that prevents amplification of short molecules thus enriching the library for longer transcripts. We adapted this method for Nanopore cDNA library preparation and show that not only is PCR efficiency increased but gene body coverage is dramatically improved. The results show that implementation of this simple strategy will result in better quality full-length RNA sequencing data and make full-length transcript sequencing possible for most of sequenced reads.

INTS13 variants causing a recessive developmental ciliopathy disrupt assembly of the Integrator complex

Oral-facial-digital (OFD) syndromes are a heterogeneous group of congenital disorders characterized by malformations of the face and oral cavity, and digit anomalies. Mutations within 12 cilia-related genes have been identified that cause several types of OFD, suggesting that OFDs constitute a subgroup of developmental ciliopathies. Through homozygosity mapping and exome sequencing of two families with variable OFD type 2, we identified distinct germline variants in INTS13, a subunit of the Integrator complex. This multiprotein complex associates with RNA Polymerase II and cleaves nascent RNA to modulate gene expression. We determined that INTS13 utilizes its C-terminus to bind the Integrator cleavage module, which is disrupted by the identified germline variants p.S652L and p.K668Nfs*9. Depletion of INTS13 disrupts ciliogenesis in human cultured cells and causes dysregulation of a broad collection of ciliary genes. Accordingly, its knockdown in Xenopus embryos leads to motile cilia anomalies. Altogether, we show that mutations in INTS13 cause an autosomal recessive ciliopathy, which reveals key interactions between components of the Integrator complex.

Click Chemistry Enables Rapid Amplification of Full-Length Reverse Transcripts for Long-Read Third Generation Sequencing

Here we describe the development of a novel click chemistry-based method for the generation and amplification of full-length cDNA libraries from total RNA, while avoiding the need for problematic template-switching (TS) reactions. Compared with prior efforts, our method involves neither random priming nor stochastic cDNA termination, thus enabling amplification of transcripts that were previously inaccessible via related click chemistry-based RNA sequencing techniques. A key modification involving the use of PCR primers containing two overhanging 3'-nucleotides substantially improved the read-through compatibility of the 1,4-disubstituted 1,2,3-triazole-containing cDNA, where such modifications typically hinder amplification. This allowed us to more than double the possible insert size compared with the state-of-the art click chemistry-based technique, PAC-seq. Furthermore, our method performed on par with a commercially available PCR-cDNA RNA sequencing kit, as determined by Oxford Nanopore sequencing. Given the known advantages of PAC-seq, namely, suppression of PCR artifacts, we anticipate that our contribution could enable diverse applications including improved analyses of mRNA splicing variants and fusion transcripts.

Tau Modulates mRNA Transcription, Alternative Polyadenylation Profiles of hnRNPs, Chromatin Remodeling and Spliceosome Complexes

Tau protein is a known contributor in several neurodegenerative diseases, including Alzheimer’s disease (AD) and frontotemporal dementia (FTD). It is well-established that tau forms pathological aggregates and fibrils in these diseases. Tau has been observed within the nuclei of neurons, but there is a gap in understanding regarding the mechanism by which tau modulates transcription. We are interested in the P301L mutation of tau, which has been associated with FTD and increased tau aggregation. Our study utilized tau-inducible HEK (iHEK) cells to reveal that WT and P301L tau distinctively alter the transcription and alternative polyadenylation (APA) profiles of numerous nuclear precursors mRNAs, which then translate to form proteins involved in chromatin remodeling and splicing. We isolated total mRNA before and after over-expressing tau and then performed Poly(A)-ClickSeq (PAC-Seq) to characterize mRNA expression and APA profiles. We characterized changes in Gene Ontology (GO) pathways using EnrichR and Gene Set Enrichment Analysis (GSEA). We observed that P301L tau up-regulates genes associated with reactive oxygen species responsiveness as well as genes involved in dendrite, microtubule, and nuclear body/speckle formation. The number of genes regulated by WT tau is greater than the mutant form, which indicates that the P301L mutation causes loss-of-function at the transcriptional level. WT tau up-regulates genes contributing to cytoskeleton-dependent intracellular transport, microglial activation, microtubule and nuclear chromatin organization, formation of nuclear bodies and speckles. Interestingly, both WT and P301L tau commonly down-regulate genes responsible for ubiquitin-proteosome system. In addition, WT tau significantly down-regulates several genes implicated in chromatin remodeling and nucleosome organization. Although there are limitations inherent to the model systems used, this study will improve understanding regarding the nuclear impact of tau at the transcriptional and post-transcriptional level. This study also illustrates the potential impact of P301L tau on the human brain genome during early phases of pathogenesis.

Enhancing the Antiviral Potency of Nucleobases for Potential Broad-Spectrum Antiviral Therapies

Broad-spectrum antiviral therapies hold promise as a first-line defense against emerging viruses by blunting illness severity and spread until vaccines and virus-specific antivirals are developed. The nucleobase favipiravir, often discussed as a broad-spectrum inhibitor, was not effective in recent clinical trials involving patients infected with Ebola virus or SARS-CoV-2. A drawback of favipiravir use is its rapid clearance before conversion to its active nucleoside-5′-triphosphate form. In this work, we report a synergistic reduction of flavivirus (dengue, Zika), orthomyxovirus (influenza A), and coronavirus (HCoV-OC43 and SARS-CoV-2) replication when the nucleobases favipiravir or T-1105 were combined with the antimetabolite 6-methylmercaptopurine riboside (6MMPr). The 6MMPr/T-1105 combination increased the C-U and G-A mutation frequency compared to treatment with T-1105 or 6MMPr alone. A further analysis revealed that the 6MMPr/T-1105 co-treatment reduced cellular purine nucleotide triphosphate synthesis and increased conversion of the antiviral nucleobase to its nucleoside-5′-monophosphate, -diphosphate, and -triphosphate forms. The 6MMPr co-treatment specifically increased production of the active antiviral form of the nucleobases (but not corresponding nucleosides) while also reducing levels of competing cellular NTPs to produce the synergistic effect. This in-depth work establishes a foundation for development of small molecules as possible co-treatments with nucleobases like favipiravir in response to emerging RNA virus infections.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called ‘Tiled-ClickSeq’, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Application and design considerations for 3'-end sequencing using click-chemistry

Over the past 15 years, investigations into alternative polyadenylation (APA) and its function in cellular physiology and pathology have greatly expanded due to the emergent appreciation of its key role in driving transcriptomic diversity. This growth has necessitated the development of new technologies capable of monitoring cleavage and polyadenylation events genome-wide. Advancements in approaches include both the creation of computational tools to re-analyze RNA-seq to identify APA events as well as targeted sequencing approaches customized to focus on the 3'-end of mRNA. Here we describe a streamlined protocol for polyA-Click-seq (PAC-seq), which utilizes click-chemistry to create mRNA 3'-ends sequencing libraries. Importantly, we offer additional considerations not present in our previous protocols including the use of spike-ins, unique molecular identifier primers, and guidance for appropriate depth of PAC-seq. In conjunction with the companion chapter on PolyA-miner (Yalamanchili et al., 2021) to computationally analyze PAC-seq data, we provide a complete experimental pipeline to analyze mRNA 3'-end usage in eukaryotic cells.

Integrator Recruits Protein Phosphatase 2A to Prevent Pause Release and Facilitate Transcription Termination

Efficient release of promoter-proximally paused RNA Pol II into productive elongation is essential for gene expression. Recently, we reported that the Integrator complex can bind paused RNA Pol II and drive premature transcription termination, potently attenuating the activity of target genes. Premature termination requires RNA cleavage by the endonuclease subunit of Integrator, but the roles of other Integrator subunits in gene regulation have yet to be elucidated. Here we report that Integrator subunit 8 (IntS8) is critical for transcription repression and required for association with protein phosphatase 2A (PP2A). We find that Integrator-bound PP2A dephosphorylates the RNA Pol II C-terminal domain and Spt5, preventing the transition to productive elongation. Thus, blocking PP2A association with Integrator stimulates pause release and gene activity. These results reveal a second catalytic function associated with Integrator-mediated transcription termination and indicate that control of productive elongation involves active competition between transcriptional kinases and phosphatases.

Mapping RNA-capsid interactions and RNA secondary structure within virus particles using next-generation sequencing

To characterize RNA-capsid binding sites genome-wide within mature RNA virus particles, we have developed a Next-Generation Sequencing (NGS) platform: viral Photo-Activatable Ribonucleoside CrossLinking (vPAR-CL). In vPAR-CL, 4-thiouridine is incorporated into the encapsidated genomes of virus particles and subsequently UV-crosslinked to adjacent capsid proteins. We demonstrate that vPAR-CL can readily and reliably identify capsid binding sites in genomic viral RNA by detecting crosslink-specific uridine to cytidine transitions in NGS data. Using Flock House virus (FHV) as a model system, we identified highly consistent and significant vPAR-CL signals across virus RNA genome, indicating a clear tropism of the encapsidated RNA genome. Certain interaction sites coincide with previously identified functional RNA motifs. We additionally performed dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) to generate a high-resolution profile of single-stranded genomic RNA inside viral particles. Combining vPAR-CL and DMS-MaPseq reveals that the predominant RNA-capsid interaction sites favored double-stranded RNA regions. We disrupted secondary structures associated with vPAR-CL sites using synonymous mutations, resulting in varied effects to virus replication, propagation and packaging. Certain mutations showed substantial deficiency in virus replication, suggesting these RNA-capsid sites are multifunctional. These provide further evidence to support that FHV packaging and replication are highly coordinated and inter-dependent events.

The Integrator Complex Attenuates Promoter-Proximal Transcription at Protein-Coding Genes

The transition of RNA polymerase II (Pol II) from initiation to productive elongation is a central, regulated step in metazoan gene expression. At many genes, Pol II pauses stably in early elongation, remaining engaged with the 25- to 60-nt-long nascent RNA for many minutes while awaiting signals for release into the gene body. However, 15%-20% of genes display highly unstable promoter Pol II, suggesting that paused polymerase might dissociate from template DNA at these promoters and release a short, non-productive mRNA. Here, we report that paused Pol II can be actively destabilized by the Integrator complex. Specifically, we present evidence that Integrator utilizes its RNA endonuclease activity to cleave nascent RNA and drive termination of paused Pol II. These findings uncover a previously unappreciated mechanism of metazoan gene repression, akin to bacterial transcription attenuation, wherein promoter-proximal Pol II is prevented from entering productive elongation through factor-regulated termination.

The Integrator Complex Attenuates Promoter-Proximal Transcription at Protein-Coding Genes

The transition of RNA polymerase II (Pol II) from initiation to productive elongation is a central, regulated step in metazoan gene expression. At many genes, Pol II pauses stably in early elongation, remaining engaged with the 25- to 60-nt-long nascent RNA for many minutes while awaiting signals for release into the gene body. However, 15%-20% of genes display highly unstable promoter Pol II, suggesting that paused polymerase might dissociate from template DNA at these promoters and release a short, non-productive mRNA. Here, we report that paused Pol II can be actively destabilized by the Integrator complex. Specifically, we present evidence that Integrator utilizes its RNA endonuclease activity to cleave nascent RNA and drive termination of paused Pol II. These findings uncover a previously unappreciated mechanism of metazoan gene repression, akin to bacterial transcription attenuation, wherein promoter-proximal Pol II is prevented from entering productive elongation through factor-regulated termination.

Development of Poly(A)-ClickSeq as a tool enabling simultaneous genome-wide poly(A)-site identification and differential expression analysis

The use of RNA-seq as a generalized tool to measure the differential expression of genes has essentially replaced the use of the microarray. Despite the acknowledged technical advantages to this approach, RNA-seq library preparation remains mostly conducted by core facilities rather than in the laboratory due to the infrastructure, expertise and time required per sample. We have recently described two 'click-chemistry' based library construction methods termed ClickSeq and Poly(A)-ClickSeq (PAC-seq) as alternatives to conventional RNA-seq that are both cost effective and rely on straightforward reagents readily available to most labs. ClickSeq is random-primed and can sequence any (unfragmented) RNA template, while PAC-seq is targeted to poly(A) tails of mRNAs. Here, we further develop PAC-seq as a platform that allows for simultaneous mapping of poly(A) sites and the measurement of differential expression of genes. We provide a detailed protocol, descriptions of appropriate computational pipelines, and a proof-of-principle dataset to illustrate the technique. PAC-seq offers a unique advantage over other 3' end mapping protocols in that it does not require additional purification, selection, or fragmentation steps allowing sample preparation directly from crude total cellular RNA. We have shown that PAC-seq is able to accurately and sensitively count transcripts for differential gene expression analysis, as well as identify alternative poly(A) sites and determine the precise nucleotides of the poly(A) tail boundaries.

A Survey of Virus Recombination Uncovers Canonical Features of Artificial Chimeras Generated During Deep Sequencing Library Preparation

Chimeric reads can be generated by in vitro recombination during the preparation of high-throughput sequencing libraries. Our attempt to detect biological recombination between the genomes of dengue virus (DENV; +ssRNA genome) and its mosquito host using the Illumina Nextera sequencing library preparation kit revealed that most, if not all, detected host-virus chimeras were artificial. Indeed, these chimeras were not more frequent than with control RNA from another species (a pillbug), which was never in contact with DENV RNA prior to the library preparation. The proportion of chimera types merely reflected those of the three species among sequencing reads. Chimeras were frequently characterized by the presence of 1-20 bp microhomology between recombining fragments. Within-species chimeras mostly involved fragments in opposite orientations and located less than 100 bp from each other in the parental genome. We found similar features in published datasets using two other viruses: Ebola virus (EBOV; -ssRNA genome) and a herpesvirus (dsDNA genome), both produced with the Illumina Nextera protocol. These canonical features suggest that artificial chimeras are generated by intra-molecular template switching of the DNA polymerase during the PCR step of the Nextera protocol. Finally, a published Illumina dataset using the Flock House virus (FHV; +ssRNA genome) generated with a protocol preventing artificial recombination revealed the presence of 1-10 bp microhomology motifs in FHV-FHV chimeras, but very few recombining fragments were in opposite orientations. Our analysis uncovered sequence features characterizing recombination breakpoints in short-read sequencing datasets, which can be helpful to evaluate the presence and extent of artificial recombination.