Adult Human, but Not Rodent, Spermatogonial Stem Cells Retain States with a Foetal-like Signature.

Spermatogenesis involves a complex process of cellular differentiation maintained by spermatogonial stem cells (SSCs). Being critical to male reproduction, it is generally assumed that spermatogenesis starts and ends in equivalent transcriptional states in related species. Based on single-cell gene expression profiling, it has been proposed that undifferentiated human spermatogonia can be subclassified into four heterogenous subtypes, termed states 0, 0A, 0B, and 1. To increase the resolution of the undifferentiated compartment and trace the origin of the spermatogenic trajectory, we re-analysed the single-cell (sc) RNA-sequencing libraries of 34 post-pubescent human testes to generate an integrated atlas of germ cell differentiation. We then used this atlas to perform comparative analyses of the putative SSC transcriptome both across human development (using 28 foetal and pre-pubertal scRNA-seq libraries) and across species (including data from sheep, pig, buffalo, rhesus and cynomolgus macaque, rat, and mouse). Alongside its detailed characterisation, we show that the transcriptional heterogeneity of the undifferentiated spermatogonial cell compartment varies not only between species but across development. Our findings associate 'state 0B' with a suppressive transcriptomic programme that, in adult humans, acts to functionally oppose proliferation and maintain cells in a ready-to-react state. Consistent with this conclusion, we show that human foetal germ cells-which are mitotically arrested-can be characterised solely as state 0B. While germ cells with a state 0B signature are also present in foetal mice (and are likely conserved at this stage throughout mammals), they are not maintained into adulthood. We conjecture that in rodents, the foetal-like state 0B differentiates at birth into the renewing SSC population, whereas in humans it is maintained as a reserve population, supporting testicular homeostasis over a longer reproductive lifespan while reducing mutagenic load. Together, these results suggest that SSCs adopt differing evolutionary strategies across species to ensure fertility and genome integrity over vastly differing life histories and reproductive timeframes.

Near telomere-to-telomere genome assemblies of two Chlorella species unveil the composition and evolution of centromeres in green algae.

BACKGROUND: Centromeres play a crucial and conserved role in cell division, although their composition and evolutionary history in green algae, the evolutionary ancestors of land plants, remains largely unknown. RESULTS: We constructed near telomere-to-telomere (T2T) assemblies for two Trebouxiophyceae species, Chlorella sorokiniana NS4-2 and Chlorella pyrenoidosa DBH, with chromosome numbers of 12 and 13, and genome sizes of 58.11 Mb and 53.41 Mb, respectively. We identified and validated their centromere sequences using CENH3 ChIP-seq and found that, similar to humans and higher plants, the centromeric CENH3 signals of green algae display a pattern of hypomethylation. Interestingly, the centromeres of both species largely comprised transposable elements, although they differed significantly in their composition. Species within the Chlorella genus display a more diverse centromere composition, with major constituents including members of the LTR/Copia, LINE/L1, and LINE/RTEX families. This is in contrast to green algae including Chlamydomonas reinhardtii, Coccomyxa subellipsoidea, and Chromochloris zofingiensis, in which centromere composition instead has a pronounced single-element composition. Moreover, we observed significant differences in the composition and structure of centromeres among chromosomes with strong collinearity within the Chlorella genus, suggesting that centromeric sequence evolves more rapidly than sequence in non-centromeric regions. CONCLUSIONS: This study not only provides high-quality genome data for comparative genomics of green algae but gives insight into the composition and evolutionary history of centromeres in early plants, laying an important foundation for further research on their evolution.

De novo and somatic structural variant discovery with SVision-pro.

Long-read-based de novo and somatic structural variant (SV) discovery remains challenging, necessitating genomic comparison between samples. We developed SVision-pro, a neural-network-based instance segmentation framework that represents genome-to-genome-level sequencing differences visually and discovers SV comparatively between genomes without any prerequisite for inference models. SVision-pro outperforms state-of-the-art approaches, in particular, the resolving of complex SVs is improved, with low Mendelian error rates, high sensitivity of low-frequency SVs and reduced false-positive rates compared with SV merging approaches.

Haplotype-resolved assemblies and variant benchmark of a Chinese Quartet.

BACKGROUND: Recent state-of-the-art sequencing technologies enable the investigation of challenging regions in the human genome and expand the scope of variant benchmarking datasets. Herein, we sequence a Chinese Quartet, comprising two monozygotic twin daughters and their biological parents, using four short and long sequencing platforms (Illumina, BGI, PacBio, and Oxford Nanopore Technology). RESULTS: The long reads from the monozygotic twin daughters are phased into paternal and maternal haplotypes using the parent-child genetic map and for each haplotype. We also use long reads to generate haplotype-resolved whole-genome assemblies with completeness and continuity exceeding that of GRCh38. Using this Quartet, we comprehensively catalogue the human variant landscape, generating a dataset of 3,962,453 SNVs, 886,648 indels (< 50 bp), 9726 large deletions (≥ 50 bp), 15,600 large insertions (≥ 50 bp), 40 inversions, 31 complex structural variants, and 68 de novo mutations which are shared between the monozygotic twin daughters. Variants underrepresented in previous benchmarks owing to their complexity-including those located at long repeat regions, complex structural variants, and de novo mutations-are systematically examined in this study. CONCLUSIONS: In summary, this study provides high-quality haplotype-resolved assemblies and a comprehensive set of benchmarking resources for two Chinese monozygotic twin samples which, relative to existing benchmarks, offers expanded genomic coverage and insight into complex variant categories.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy.

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Fine-tuning germline mutation rates across evolution.

The germline mutation rate (GMR) sets the pace at which mutations, the raw material of evolution, are introduced into the genome. By sequencing a dataset of unprecedently broad phylogenetic scope, Bergeron et al. estimated species-specific GMR, offering numerous insights into how this parameter shapes and is shaped by life-history traits.

Personalized recurrence risk assessment following the birth of a child with a pathogenic de novo mutation.

Following the diagnosis of a paediatric disorder caused by an apparently de novo mutation, a recurrence risk of 1-2% is frequently quoted due to the possibility of parental germline mosaicism; but for any specific couple, this figure is usually incorrect. We present a systematic approach to providing individualized recurrence risk. By combining locus-specific sequencing of multiple tissues to detect occult mosaicism with long-read sequencing to determine the parent-of-origin of the mutation, we show that we can stratify the majority of couples into one of seven discrete categories associated with substantially different risks to future offspring. Among 58 families with a single affected offspring (representing 59 de novo mutations in 49 genes), the recurrence risk for 35 (59%) was decreased below 0.1%, but increased owing to parental mixed mosaicism for 5 (9%)-that could be quantified in semen for paternal cases (recurrence risks of 5.6-12.1%). Implementation of this strategy offers the prospect of driving a major transformation in the practice of genetic counselling.

Fibrillin-1 and asprosin, novel players in metabolic syndrome.

Fibrillin-1 is a major component of the extracellular microfibrils, where it interacts with other extracellular matrix proteins to provide elasticity to connective tissues, and regulates the bioavailability of TGFß family members. A peptide consisting of the C-terminal 140 amino acids of fibrillin-1 has recently been identified as a glucogenic hormone, secreted from adipose tissue during fasting and targeting the liver to release glucose. This fragment, called asprosin, also signals in the hypothalamus to stimulate appetite. Asprosin levels are correlated with many of the pathologies indicative of metabolic syndrome, including insulin resistance and obesity. Previous studies and reviews have addressed the therapeutic potential of asprosin as a target in obesity, diabetes and related conditions without considering mechanisms underlying the relationship between generation of asprosin and expression of the much larger fibrillin-1 protein. Profibrillin-1 undergoes obligatory cleavage at the cell surface as part of its assembly into microfibrils, producing the asprosin peptide as well as mature fibrillin-1. Patterns of FBN1 mRNA expression are inconsistent with the necessity for regulated release of asprosin. The asprosin peptide may be protected from degradation in adipose tissue. We present evidence for an alternative possibility, that asprosin mRNA is generated independently from an internal promoter within the 3' end of the FBN1 gene, which would allow for regulation independent of fibrillin-synthesis and is more economical of cellular resources. The discovery of asprosin opened exciting possibilities for treatment of metabolic syndrome related conditions, but there is much to be understood before such therapies could be introduced into the clinic.

Modern Solutions for Ancient Pathogens: Direct Pathogen Sequencing for Diagnosis of Lepromatous Leprosy and Cerebral Coenurosis.

Microbes unculturable in vitro remain diagnostically challenging, dependent historically on clinical findings, histology, or targeted molecular detection. We applied whole-genome sequencing directly from tissue to diagnose infections with mycobacteria (leprosy) and parasites (coenurosis). Direct pathogen DNA sequencing provides flexible solutions to diagnosis of difficult pathogens in diverse contexts.

Generation and network analysis of an RNA-seq transcriptional atlas for the rat.

The laboratory rat is an important model for biomedical research. To generate a comprehensive rat transcriptomic atlas, we curated and downloaded 7700 rat RNA-seq datasets from public repositories, downsampled them to a common depth and quantified expression. Data from 585 rat tissues and cells, averaged from each BioProject, can be visualized and queried at http://biogps.org/ratatlas. Gene co-expression network (GCN) analysis revealed clusters of transcripts that were tissue or cell type restricted and contained transcription factors implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes. Many of these clusters overlap with previous data from analysis of other species, while some (e.g. expressed specifically in immune cells, retina/pineal gland, pituitary and germ cells) are unique to these data. GCN analysis on large subsets of the data related specifically to liver, nervous system, kidney, musculoskeletal system and cardiovascular system enabled deconvolution of cell type-specific signatures. The approach is extensible and the dataset can be used as a point of reference from which to analyse the transcriptomes of cell types and tissues that have not yet been sampled. Sets of strictly co-expressed transcripts provide a resource for critical interpretation of single-cell RNA-seq data.

Contrasting gene-level signatures of selection with reproductive fitness.

Selection leaves signatures in the DNA sequence of genes, with many test statistics devised to detect its action. While these statistics are frequently used to support hypotheses about the adaptive significance of particular genes, the effect these genes have on reproductive fitness is rarely quantified experimentally. Consequently, it is unclear how gene-level signatures of selection are associated with empirical estimates of gene effect on fitness. Eukaryotic data sets that permit this comparison are very limited. Using the model plant Arabidopsis thaliana, for which these resources are available, we calculated seven gene-level substitution and polymorphism-based statistics commonly used to infer selection (dN/dS, NI, DOS, Tajima's D, Fu and Li's D*, Fay and Wu's H, and Zeng's E) and, using knockout lines, compared these to gene-level estimates of effect on fitness. We found that consistent with expectations, essential genes were more likely to be classified as negatively selected. By contrast, using 379 Arabidopsis genes for which data was available, we found no evidence that genes predicted to be positively selected had a significantly different effect on fitness than genes evolving more neutrally. We discuss these results in the context of the analytic challenges posed by Arabidopsis, one of the only systems in which this study could be conducted, and advocate for examination in additional systems. These results are relevant to the evaluation of genome-wide studies across species where experimental fitness data is unavailable, as well as highlighting an increasing need for the latter.

Rats exhibit age-related mosaic loss of chromosome Y.

Mosaic loss of the Y chromosome (LOY) is the most frequent chromosomal aberration in aging men and is strongly correlated with mortality and disease. To date, studies of LOY have only been performed in humans, and so it is unclear whether LOY is a natural consequence of our relatively long lifespan or due to exposure to human-specific external stressors. Here, we explored whether LOY could be detected in rats. We applied a locus-specific PCR and target sequencing approach that we used as a proxy to estimate LOY in 339 samples covering eleven tissues from young and old individuals. We detected LOY in four tissues of older rats. To confirm the results from the PCR screening, we re-sequenced 60 full genomes from old rats, which revealed that the Y chromosome is the sole chromosome with low copy numbers. Finally, our results suggest that LOY is associated with other structural aberrations on the Y chromosome and possibly linked to the mosaic loss of the X chromosome. This is the first report, to our knowledge, demonstrating that the patterns of LOY observed in aging men are also present in a rodent, and conclude that LOY may be a natural process in placental mammals.

Generalizable characteristics of false-positive bacterial variant calls.

Minimizing false positives is a critical issue when variant calling as no method is without error. It is common practice to post-process a variant-call file (VCF) using hard filter criteria intended to discriminate true-positive (TP) from false-positive (FP) calls. These are applied on the simple principle that certain characteristics are disproportionately represented among the set of FP calls and that a user-chosen threshold can maximize the number detected. To provide guidance on this issue, this study empirically characterized all false SNP and indel calls made using real Illumina sequencing data from six disparate species and 166 variant-calling pipelines (the combination of 14 read aligners with up to 13 different variant callers, plus four 'all-in-one' pipelines). We did not seek to optimize filter thresholds but instead to draw attention to those filters of greatest efficacy and the pipelines to which they may most usefully be applied. In this respect, this study acts as a coda to our previous benchmarking evaluation of bacterial variant callers, and provides general recommendations for effective practice. The results suggest that, of the pipelines analysed in this study, the most straightforward way of minimizing false positives would simply be to use Snippy. We also find that a disproportionate number of false calls, irrespective of the variant-calling pipeline, are located in the vicinity of indels, and highlight this as an issue for future development.

CRISPR-Cas9 Editing of Human Histone Deubiquitinase Gene USP16 in Human Monocytic Leukemia Cell Line THP-1.

USP16 is a histone deubiquitinase which facilitates G2/M transition during the cell cycle, regulates DNA damage repair and contributes to inducible gene expression. We mutated the USP16 gene in a high differentiation clone of the acute monocytic leukemia cell line THP-1 using the CRISPR-Cas9 system and generated four homozygous knockout clones. All were able to proliferate and to differentiate in response to phorbol ester (PMA) treatment. One line was highly proliferative prior to PMA treatment and shut down proliferation upon differentiation, like wild type. Three clones showed sustained expression of the progenitor cell marker MYB, indicating that differentiation had not completely blocked proliferation in these clones. Network analysis of transcriptomic differences among wild type, heterozygotes and homozygotes showed clusters of genes that were up- or down-regulated after differentiation in all cell lines. Prior to PMA treatment, the homozygous clones had lower levels than wild type of genes relating to metabolism and mitochondria, including SRPRB, encoding an interaction partner of USP16. There was also apparent loss of interferon signaling. In contrast, a number of genes were up-regulated in the homozygous cells compared to wild type at baseline, including other deubiquitinases (USP12, BAP1, and MYSM1). However, three homozygotes failed to fully induce USP3 during differentiation. Other network clusters showed effects prior to or after differentiation in the homozygous clones. Thus the removal of USP16 affected the transcriptome of the cells, although all these lines were able to survive, which suggests that the functions attributed to USP16 may be redundant. Our analysis indicates that the leukemic line can adapt to the extreme selection pressure applied by the loss of USP16, and the harsh conditions of the gene editing and selection protocol, through different compensatory pathways. Similar selection pressures occur during the evolution of a cancer in vivo, and our results can be seen as a case study in leukemic cell adaptation. USP16 has been considered a target for cancer chemotherapy, but our results suggest that treatment would select for escape mutants that are resistant to USP16 inhibitors.

CSF1R-dependent macrophages control postnatal somatic growth and organ maturation.

Homozygous mutation of the Csf1r locus (Csf1rko) in mice, rats and humans leads to multiple postnatal developmental abnormalities. To enable analysis of the mechanisms underlying the phenotypic impacts of Csf1r mutation, we bred a rat Csf1rko allele to the inbred dark agouti (DA) genetic background and to a Csf1r-mApple reporter transgene. The Csf1rko led to almost complete loss of embryonic macrophages and ablation of most adult tissue macrophage populations. We extended previous analysis of the Csf1rko phenotype to early postnatal development to reveal impacts on musculoskeletal development and proliferation and morphogenesis in multiple organs. Expression profiling of 3-week old wild-type (WT) and Csf1rko livers identified 2760 differentially expressed genes associated with the loss of macrophages, severe hypoplasia, delayed hepatocyte maturation, disrupted lipid metabolism and the IGF1/IGF binding protein system. Older Csf1rko rats developed severe hepatic steatosis. Consistent with the developmental delay in the liver Csf1rko rats had greatly-reduced circulating IGF1. Transfer of WT bone marrow (BM) cells at weaning without conditioning repopulated resident macrophages in all organs, including microglia in the brain, and reversed the mutant phenotypes enabling long term survival and fertility. WT BM transfer restored osteoclasts, eliminated osteopetrosis, restored bone marrow cellularity and architecture and reversed granulocytosis and B cell deficiency. Csf1rko rats had an elevated circulating CSF1 concentration which was rapidly reduced to WT levels following BM transfer. However, CD43hi non-classical monocytes, absent in the Csf1rko, were not rescued and bone marrow progenitors remained unresponsive to CSF1. The results demonstrate that the Csf1rko phenotype is autonomous to BM-derived cells and indicate that BM contains a progenitor of tissue macrophages distinct from hematopoietic stem cells. The model provides a unique system in which to define the pathways of development of resident tissue macrophages and their local and systemic roles in growth and organ maturation.

The Mononuclear Phagocyte System of the Rat.

The laboratory rat continues to be the model of choice for many studies of physiology, behavior, and complex human diseases. Cells of the mononuclear phagocyte system (MPS; monocytes, macrophages, and dendritic cells) are abundant residents in every tissue in the body and regulate postnatal development, homeostasis, and innate and acquired immunity. Recruitment and proliferation of MPS cells is an essential component of both initiation and resolution of inflammation. The large majority of current knowledge of MPS biology is derived from studies of inbred mice, but advances in technology and resources have eliminated many of the advantages of the mouse as a model. In this article, we review the tools available and the current state of knowledge of development, homeostasis, regulation, and diversity within the MPS of the rat.

Analysis of homozygous and heterozygous Csf1r knockout in the rat as a model for understanding microglial function in brain development and the impacts of human CSF1R mutations.

Mutations in the human CSF1R gene have been associated with dominant and recessive forms of neurodegenerative disease. Here we describe the impacts of Csf1r mutation in the rat on development of the brain. Diffusion imaging indicated small reductions in major fiber tracts that may be associated in part with ventricular enlargement. RNA-seq profiling revealed a set of 105 microglial markers depleted in all brain regions of the Csf1rko rats. There was no evidence of region or sex-specific expression of microglia-associated transcripts. Other than the microglial signature, Csf1rko had no effect on any neuronal or region-specific transcript cluster. Expression of markers of oligodendrocytes, astrocytes, dopaminergic neurons and Purkinje cells was minimally affected. However, there were defects in dendritic arborization of doublecortin-positive neurogenic precursors and expression of poly-sialylated neural cell adhesion molecule (PS-NCAM) in the dentate gyrus of the hippocampus. Heterozygous Csf1rko rats had no detectable brain phenotype. We conclude that most brain developmental processes occur normally in the absence of microglia and that CSF1R haploinsufficiency is unlikely to cause leukoencephalopathy.

Quantitative trait loci and transcriptome signatures associated with avian heritable resistance to Campylobacter.

Campylobacter is the leading cause of bacterial foodborne gastroenteritis worldwide. Handling or consumption of contaminated poultry meat is a key risk factor for human campylobacteriosis. One potential control strategy is to select poultry with increased resistance to Campylobacter. We associated high-density genome-wide genotypes (600K single nucleotide polymorphisms) of 3000 commercial broilers with Campylobacter load in their caeca. Trait heritability was modest but significant (h2 = 0.11 ± 0.03). Results confirmed quantitative trait loci (QTL) on chromosomes 14 and 16 previously identified in inbred chicken lines, and detected two additional QTLs on chromosomes 19 and 26. RNA-Seq analysis of broilers at the extremes of colonisation phenotype identified differentially transcribed genes within the QTL on chromosome 16 and proximal to the major histocompatibility complex (MHC) locus. We identified strong cis-QTLs located within MHC suggesting the presence of cis-acting variation in MHC class I and II and BG genes. Pathway and network analyses implicated cooperative functional pathways and networks in colonisation, including those related to antigen presentation, innate and adaptive immune responses, calcium, and renin-angiotensin signalling. While co-selection for enhanced resistance and other breeding goals is feasible, the frequency of resistance-associated alleles was high in the population studied and non-genetic factors significantly influenced Campylobacter colonisation.

MeDAS: a Metazoan Developmental Alternative Splicing database.

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Read trimming has minimal effect on bacterial SNP-calling accuracy.

Read alignment is the central step of many analytic pipelines that perform variant calling. To reduce error, it is common practice to pre-process raw sequencing reads to remove low-quality bases and residual adapter contamination, a procedure collectively known as 'trimming'. Trimming is widely assumed to increase the accuracy of variant calling, although there are relatively few systematic evaluations of its effects and no clear consensus on its efficacy. As sequencing datasets increase both in number and size, it is worthwhile reappraising computational operations of ambiguous benefit, particularly when the scope of many analyses now routinely incorporates thousands of samples, increasing the time and cost required. Using a curated set of 17 Gram-negative bacterial genomes, this study initially evaluated the impact of four read-trimming utilities (Atropos, fastp, Trim Galore and Trimmomatic), each used with a range of stringencies, on the accuracy and completeness of three bacterial SNP-calling pipelines. It was found that read trimming made only small, and statistically insignificant, increases in SNP-calling accuracy even when using the highest-performing pre-processor in this study, fastp. To extend these findings, >6500 publicly archived sequencing datasets from Escherichia coli, Mycobacterium tuberculosis and Staphylococcus aureus were re-analysed using a common analytic pipeline. Of the approximately 125 million SNPs and 1.25 million indels called across all samples, the same bases were called in 98.8 and 91.9 % of cases, respectively, irrespective of whether raw reads or trimmed reads were used. Nevertheless, the proportion of mixed calls (i.e. calls where <100 % of the reads support the variant allele; considered a proxy of false positives) was significantly reduced after trimming, which suggests that while trimming rarely alters the set of variant bases, it can affect the proportion of reads supporting each call. It was concluded that read quality- and adapter-trimming add relatively little value to a SNP-calling pipeline and may only be necessary if small differences in the absolute number of SNP calls, or the false call rate, are critical. Broadly similar conclusions can be drawn about the utility of trimming to an indel-calling pipeline. Read trimming remains routinely performed prior to variant calling likely out of concern that doing otherwise would typically have negative consequences. While historically this may have been the case, the data in this study suggests that read trimming is not always a practical necessity.