Development and validation of a high throughput SARS-CoV-2 whole genome sequencing workflow in a clinical laboratory.

Monitoring new mutations in SARS-CoV-2 provides crucial information for identifying diagnostic and therapeutic targets and important insights to achieve a more effective COVID-19 control strategy. Next generation sequencing (NGS) technologies have been widely used for whole genome sequencing (WGS) of SARS-CoV-2. While various NGS methods have been reported, one chief limitation has been the complexity of the workflow, limiting the scalability. Here, we overcome this limitation by designing a laboratory workflow optimized for high-throughput studies. The workflow utilizes modified ARTIC network v3 primers for SARS-CoV-2 whole genome amplification. NGS libraries were prepared by a 2-step PCR method, similar to a previously reported tailed PCR method, with further optimizations to improve amplicon balance, to minimize amplicon dropout for viral genomes harboring primer-binding site mutation(s), and to integrate robotic liquid handlers. Validation studies demonstrated that the optimized workflow can process up to 2688 samples in a single sequencing run without compromising sensitivity and accuracy and with fewer amplicon dropout events compared to the standard ARTIC protocol. We additionally report results for over 65,000 SARS-CoV-2 whole genome sequences from clinical specimens collected in the United States between January and September of 2021, as part of an ongoing national genomics surveillance effort.

Analytical validation and performance characteristics of a 48-gene next-generation sequencing panel for detecting potentially actionable genomic alterations in myeloid neoplasms.

Identification of genomic mutations by molecular testing plays an important role in diagnosis, prognosis, and treatment of myeloid neoplasms. Next-generation sequencing (NGS) is an efficient method for simultaneous detection of clinically significant genomic mutations with high sensitivity. Various NGS based in-house developed and commercial myeloid neoplasm panels have been integrated into routine clinical practice. However, some genes frequently mutated in myeloid malignancies are particularly difficult to sequence with NGS panels (e.g., CEBPA, CARL, and FLT3). We report development and validation of a 48-gene NGS panel that includes genes that are technically challenging for molecular profiling of myeloid neoplasms including acute myeloid leukemia (AML), myelodysplastic syndrome (MDS), and myeloproliferative neoplasms (MPN). Target regions were captured by hybridization with complementary biotinylated DNA baits, and NGS was performed on an Illumina NextSeq500 instrument. A bioinformatics pipeline that was developed in-house was used to detect single nucleotide variations (SNVs), insertions/deletions (indels), and FLT3 internal tandem duplications (FLT3-ITD). An analytical validation study was performed on 184 unique specimens for variants with allele frequencies ≥5%. Variants identified by the 48-gene panel were compared to those identified by a 35-gene hematologic neoplasms panel using an additional 137 unique specimens. The developed assay was applied to a large cohort (n = 2,053) of patients with suspected myeloid neoplasms. Analytical validation yielded 99.6% sensitivity (95% CI: 98.9-99.9%) and 100% specificity (95% CI: 100%). Concordance of variants detected by the 2 tested panels was 100%. Among patients with suspected myeloid neoplasms (n = 2,053), 54.5% patients harbored at least one clinically significant mutation: 77% in AML patients, 48% in MDS, and 45% in MPN. Together, these findings demonstrate that the assay can identify mutations associated with diagnosis, prognosis, and treatment options of myeloid neoplasms even in technically challenging genes.

Development and Validation of a 34-Gene Inherited Cancer Predisposition Panel Using Next-Generation Sequencing.

The use of genetic testing to identify individuals with hereditary cancer syndromes has been widely adopted by clinicians for management of inherited cancer risk. The objective of this study was to develop and validate a 34-gene inherited cancer predisposition panel using targeted capture-based next-generation sequencing (NGS). The panel incorporates genes underlying well-characterized cancer syndromes, such as BRCA1 and BRCA2 (BRCA1/2), along with more recently discovered genes associated with increased cancer risk. We performed a validation study on 133 unique specimens, including 33 with known variant status; known variants included single nucleotide variants (SNVs) and small insertions and deletions (Indels), as well as copy-number variants (CNVs). The analytical validation study achieved 100% sensitivity and specificity for SNVs and small Indels, with 100% sensitivity and 98.0% specificity for CNVs using in-house developed CNV flagging algorithm. We employed a microarray comparative genomic hybridization (aCGH) method for all specimens that the algorithm flags as CNV-positive for confirmation. In combination with aCGH confirmation, CNV detection specificity improved to 100%. We additionally report results of the first 500 consecutive specimens submitted for clinical testing with the 34-gene panel, identifying 53 deleterious variants in 13 genes in 49 individuals. Half of the detected pathogenic/likely pathogenic variants were found in BRCA1 (23%), BRCA2 (23%), or the Lynch syndrome-associated genes PMS2 (4%) and MLH1 (2%). The other half were detected in 9 other genes: MUTYH (17%), CHEK2 (15%), ATM (4%), PALB2 (4%), BARD1 (2%), CDH1 (2%), CDKN2A (2%), RAD51C (2%), and RET (2%). Our validation studies and initial clinical data demonstrate that a 34-gene inherited cancer predisposition panel can provide clinically significant information for cancer risk assessment.

Transcriptional profile of tuberculosis antigen-specific T cells reveals novel multifunctional features.

In latent tuberculosis infection (LTBI) spread of the bacteria is contained by a persistent immune response, which includes CD4(+) T cells as important contributors. In this study we show that TB-specific CD4(+) T cells have a characteristic chemokine expression signature (CCR6(+)CXCR3(+)CCR4(-)), and that the overall number of these cells is significantly increased in LTBI donors compared with healthy subjects. We have comprehensively characterized the transcriptional signature of CCR6(+)CXCR3(+)CCR4(-) cells and found significant differences to conventional Th1, Th17, and Th2 cells, but no major changes between healthy and LTBI donors. CCR6(+)CXCR3(+)CCR4(-) cells display lineage-specific signatures of both Th1 and Th17 cells, but also have a unique gene expression program, including genes associated with susceptibility to TB, enhanced T cell activation, enhanced cell survival, and induction of a cytotoxic program akin to CTL cells. Overall, the gene expression signature of CCR6(+)CXCR3(+)CCR4(-) cells reveals characteristics important for controlling latent TB infections.

Epigenomic analysis of primary human T cells reveals enhancers associated with TH2 memory cell differentiation and asthma susceptibility.

A characteristic feature of asthma is the aberrant accumulation, differentiation or function of memory CD4(+) T cells that produce type 2 cytokines (TH2 cells). By mapping genome-wide histone modification profiles for subsets of T cells isolated from peripheral blood of healthy and asthmatic individuals, we identified enhancers with known and potential roles in the normal differentiation of human TH1 cells and TH2 cells. We discovered disease-specific enhancers in T cells that differ between healthy and asthmatic individuals. Enhancers that gained the histone H3 Lys4 dimethyl (H3K4me2) mark during TH2 cell development showed the highest enrichment for asthma-associated single nucleotide polymorphisms (SNPs), which supported a pathogenic role for TH2 cells in asthma. In silico analysis of cell-specific enhancers revealed transcription factors, microRNAs and genes potentially linked to human TH2 cell differentiation. Our results establish the feasibility and utility of enhancer profiling in well-defined populations of specialized cell types involved in disease pathogenesis.

Predicting cell types and genetic variations contributing to disease by combining GWAS and epigenetic data.

Genome-wide association studies (GWASs) identify single nucleotide polymorphisms (SNPs) that are enriched in individuals suffering from a given disease. Most disease-associated SNPs fall into non-coding regions, so that it is not straightforward to infer phenotype or function; moreover, many SNPs are in tight genetic linkage, so that a SNP identified as associated with a particular disease may not itself be causal, but rather signify the presence of a linked SNP that is functionally relevant to disease pathogenesis. Here, we present an analysis method that takes advantage of the recent rapid accumulation of epigenomics data to address these problems for some SNPs. Using asthma as a prototypic example; we show that non-coding disease-associated SNPs are enriched in genomic regions that function as regulators of transcription, such as enhancers and promoters. Identifying enhancers based on the presence of the histone modification marks such as H3K4me1 in different cell types, we show that the location of enhancers is highly cell-type specific. We use these findings to predict which SNPs are likely to be directly contributing to disease based on their presence in regulatory regions, and in which cell types their effect is expected to be detectable. Moreover, we can also predict which cell types contribute to a disease based on overlap of the disease-associated SNPs with the locations of enhancers present in a given cell type. Finally, we suggest that it will be possible to re-analyze GWAS studies with much higher power by limiting the SNPs considered to those in coding or regulatory regions of cell types relevant to a given disease.

Memory T cells in latent Mycobacterium tuberculosis infection are directed against three antigenic islands and largely contained in a CXCR3+CCR6+ Th1 subset.

An understanding of the immunological footprint of Mycobacterium tuberculosis (MTB) CD4 T cell recognition is still incomplete. Here we report that human Th1 cells specific for MTB are largely contained in a CXCR3(+)CCR6(+) memory subset and highly focused on three broadly immunodominant antigenic islands, all related to bacterial secretion systems. Our results refute the notion that secreted antigens act as a decoy, since both secreted proteins and proteins comprising the secretion system itself are targeted by a fully functional T cell response. In addition, several novel T cell antigens were identified which can be of potential diagnostic use, or as vaccine antigens. These results underline the power of a truly unbiased, genome-wide, analysis of CD4 MTB recognition based on the combined use of epitope predictions, high throughput ELISPOT, and T cell libraries using PBMCs from individuals latently infected with MTB.

Comparative genomic reconstruction of transcriptional networks controlling central metabolism in the Shewanella genus.

BACKGROUND: Genome-scale prediction of gene regulation and reconstruction of transcriptional regulatory networks in bacteria is one of the critical tasks of modern genomics. The Shewanella genus is comprised of metabolically versatile gamma-proteobacteria, whose lifestyles and natural environments are substantially different from Escherichia coli and other model bacterial species. The comparative genomics approaches and computational identification of regulatory sites are useful for the in silico reconstruction of transcriptional regulatory networks in bacteria. RESULTS: To explore conservation and variations in the Shewanella transcriptional networks we analyzed the repertoire of transcription factors and performed genomics-based reconstruction and comparative analysis of regulons in 16 Shewanella genomes. The inferred regulatory network includes 82 transcription factors and their DNA binding sites, 8 riboswitches and 6 translational attenuators. Forty five regulons were newly inferred from the genome context analysis, whereas others were propagated from previously characterized regulons in the Enterobacteria and Pseudomonas spp.. Multiple variations in regulatory strategies between the Shewanella spp. and E. coli include regulon contraction and expansion (as in the case of PdhR, HexR, FadR), numerous cases of recruiting non-orthologous regulators to control equivalent pathways (e.g. PsrA for fatty acid degradation) and, conversely, orthologous regulators to control distinct pathways (e.g. TyrR, ArgR, Crp). CONCLUSIONS: We tentatively defined the first reference collection of ~100 transcriptional regulons in 16 Shewanella genomes. The resulting regulatory network contains ~600 regulated genes per genome that are mostly involved in metabolism of carbohydrates, amino acids, fatty acids, vitamins, metals, and stress responses. Several reconstructed regulons including NagR for N-acetylglucosamine catabolism were experimentally validated in S. oneidensis MR-1. Analysis of correlations in gene expression patterns helps to interpret the reconstructed regulatory network. The inferred regulatory interactions will provide an additional regulatory constrains for an integrated model of metabolism and regulation in S. oneidensis MR-1.

Comparative genomics of the dormancy regulons in mycobacteria.

In response to stresses, Mycobacterium cells become dormant. This process is regulated by the DosR transcription factor. In Mycobacterium tuberculosis, the dormancy regulon is well characterized and contains the dosR gene itself and dosS and dosT genes encoding DosR kinases, nitroreductases (acg; Rv3131), diacylglycerol acyltransferase (DGAT) (Rv3130c), and many universal stress proteins (USPs). In this study, we apply comparative genomic analysis to characterize the DosR regulons in nine Mycobacterium genomes, Rhodococcus sp. RHA1, Nocardia farcinica, and Saccharopolyspora erythraea. The regulons are highly labile, containing eight core gene groups (regulators, kinases, USPs, DGATs, nitroreductases, ferredoxins, heat shock proteins, and the orthologs of the predicted kinase [Rv2004c] from M. tuberculosis) and 10 additional genes with more restricted taxonomic distribution that are mostly involved in anaerobic respiration. The largest regulon is observed in M. marinum and the smallest in M. abscessus. Analysis of large gene families encoding USPs, nitroreductases, and DGATs demonstrates a mosaic distribution of regulated and nonregulated members, suggesting frequent acquisition and loss of DosR-binding sites.

Erratum.

This corrects the article on p. e1000281 in Vol. 4, PMID: 19043566. Hypermutable Non-Synonymous Sites Are Under Stronger Negative Selection.

Hypermutable non-synonymous sites are under stronger negative selection.

Mutation rate varies greatly between nucleotide sites of the human genome and depends both on the global genomic location and the local sequence context of a site. In particular, CpG context elevates the mutation rate by an order of magnitude. Mutations also vary widely in their effect on the molecular function, phenotype, and fitness. Independence of the probability of occurrence of a new mutation's effect has been a fundamental premise in genetics. However, highly mutable contexts may be preserved by negative selection at important sites but destroyed by mutation at sites under no selection. Thus, there may be a positive correlation between the rate of mutations at a nucleotide site and the magnitude of their effect on fitness. We studied the impact of CpG context on the rate of human-chimpanzee divergence and on intrahuman nucleotide diversity at non-synonymous coding sites. We compared nucleotides that occupy identical positions within codons of identical amino acids and only differ by being within versus outside CpG context. Nucleotides within CpG context are under a stronger negative selection, as revealed by their lower, proportionally to the mutation rate, rate of evolution and nucleotide diversity. In particular, the probability of fixation of a non-synonymous transition at a CpG site is two times lower than at a CpG site. Thus, sites with different mutation rates are not necessarily selectively equivalent. This suggests that the mutation rate may complement sequence conservation as a characteristic predictive of functional importance of nucleotide sites.

Deinococcus geothermalis: the pool of extreme radiation resistance genes shrinks.

Bacteria of the genus Deinococcus are extremely resistant to ionizing radiation (IR), ultraviolet light (UV) and desiccation. The mesophile Deinococcus radiodurans was the first member of this group whose genome was completely sequenced. Analysis of the genome sequence of D. radiodurans, however, failed to identify unique DNA repair systems. To further delineate the genes underlying the resistance phenotypes, we report the whole-genome sequence of a second Deinococcus species, the thermophile Deinococcus geothermalis, which at its optimal growth temperature is as resistant to IR, UV and desiccation as D. radiodurans, and a comparative analysis of the two Deinococcus genomes. Many D. radiodurans genes previously implicated in resistance, but for which no sensitive phenotype was observed upon disruption, are absent in D. geothermalis. In contrast, most D. radiodurans genes whose mutants displayed a radiation-sensitive phenotype in D. radiodurans are conserved in D. geothermalis. Supporting the existence of a Deinococcus radiation response regulon, a common palindromic DNA motif was identified in a conserved set of genes associated with resistance, and a dedicated transcriptional regulator was predicted. We present the case that these two species evolved essentially the same diverse set of gene families, and that the extreme stress-resistance phenotypes of the Deinococcus lineage emerged progressively by amassing cell-cleaning systems from different sources, but not by acquisition of novel DNA repair systems. Our reconstruction of the genomic evolution of the Deinococcus-Thermus phylum indicates that the corresponding set of enzymes proliferated mainly in the common ancestor of Deinococcus. Results of the comparative analysis weaken the arguments for a role of higher-order chromosome alignment structures in resistance; more clearly define and substantially revise downward the number of uncharacterized genes that might participate in DNA repair and contribute to resistance; and strengthen the case for a role in survival of systems involved in manganese and iron homeostasis.

Comparative genomic analysis of regulation of anaerobic respiration in ten genomes from three families of gamma-proteobacteria (Enterobacteriaceae, Pasteurellaceae, Vibrionaceae).

BACKGROUND: Gamma-proteobacteria, such as Escherichia coli, can use a variety of respiratory substrates employing numerous aerobic and anaerobic respiratory systems controlled by multiple transcription regulators. Thus, in E. coli, global control of respiration is mediated by four transcription factors, Fnr, ArcA, NarL and NarP. However, in other Gamma-proteobacteria the composition of global respiration regulators may be different. RESULTS: In this study we applied a comparative genomic approach to the analysis of three global regulatory systems, Fnr, ArcA and NarP. These systems were studied in available genomes containing these three regulators, but lacking NarL. So, we considered several representatives of Pasteurellaceae, Vibrionaceae and Yersinia spp. As a result, we identified new regulon members, functioning in respiration, central metabolism (glycolysis, gluconeogenesis, pentose phosphate pathway, citrate cicle, metabolism of pyruvate and lactate), metabolism of carbohydrates and fatty acids, transcriptional regulation and transport, in particular: the ATP synthase operon atpIBEFHAGCD, Na+-exporting NADH dehydrogenase operon nqrABCDEF, the D-amino acids dehydrogenase operon dadAX. Using an extension of the comparative technique, we demonstrated taxon-specific changes in regulatory interactions and predicted taxon-specific regulatory cascades. CONCLUSION: A comparative genomic technique was applied to the analysis of global regulation of respiration in ten gamma-proteobacterial genomes. Three structurally different but functionally related regulatory systems were described. A correlation between the regulon size and the position of a transcription factor in regulatory cascades was observed: regulators with larger regulons tend to occupy top positions in the cascades. On the other hand, there is no obvious link to differences in the species' lifestyles and metabolic capabilities.

Evolution of the NadR regulon in Enterobacteriaceae.

The NAD biosynthetic pathway and NAD transformations in E. coli and S. typhi are well characterized. Using comparative genomics methods we describe the NadR regulon in other Enterobacteriaceae, identity new candidate regulon members and demonstrate that even a very simple regulon covering an essential methabolic pathway could be different in closely related genomes.