Development, Validation, and Implementation of an Augmented Multiwell, Multitarget Quantitative PCR-Human Papillomavirus Genotyping Analysis through Software Automation, Data Science, and Artificial Intelligence.

The value of human papillomavirus (HPV) testing for cervical cancer screening is well established, where its use as a primary screening option or as a reflex test after atypical cytology results has recently gained wide acceptance. The importance of full genotyping and viral load determination has been demonstrated to enhance the clinical understanding of the viral infection progression during follow-up or after treatment, thereby providing clinicians with supplementary tools for optimized patient management. In this study, a new analysis method for the RIATOL quantitative PCR assay was developed, validated, and implemented in the laboratory of clinical molecular pathology at AML, under national accreditation and following the International Organization for Standardization guidelines. It presents the successful validation of a high-throughput, multitarget HPV analysis method, with enhanced accuracy on both qualitative and quantitative end results. This is achieved by software standardization and automation of PCR curve analysis and interpretation, using data science and artificial intelligence. Moreover, the user-centric functionality of the platform was demonstrated to enhance both staff training and routine analysis workflows, thereby saving time and laboratory personnel resources. Overall, the integration of the FastFinder plugin semi-automatic analysis algorithm with the RIATOL real-time quantitative PCR assay proved to be a remarkable advancement in high-throughput HPV quantification, with demonstrated capability to provide highly accurate clinical-grade results and to reduce manual variability and analysis time.

Integrated genomic DNA/RNA profiling of hematologic malignancies in the clinical setting.

The spectrum of somatic alterations in hematologic malignancies includes substitutions, insertions/deletions (indels), copy number alterations (CNAs), and a wide range of gene fusions; no current clinically available single assay captures the different types of alterations. We developed a novel next-generation sequencing-based assay to identify all classes of genomic alterations using archived formalin-fixed paraffin-embedded blood and bone marrow samples with high accuracy in a clinically relevant time frame, which is performed in our Clinical Laboratory Improvement Amendments-certified College of American Pathologists-accredited laboratory. Targeted capture of DNA/RNA and next-generation sequencing reliably identifies substitutions, indels, CNAs, and gene fusions, with similar accuracy to lower-throughput assays that focus on specific genes and types of genomic alterations. Profiling of 3696 samples identified recurrent somatic alterations that impact diagnosis, prognosis, and therapy selection. This comprehensive genomic profiling approach has proved effective in detecting all types of genomic alterations, including fusion transcripts, which increases the ability to identify clinically relevant genomic alterations with therapeutic relevance.

A genome-wide association study of variants associated with acquisition of Staphylococcus aureus bacteremia in a healthcare setting.

BACKGROUND: Humans vary in their susceptibility to acquiring Staphylococcus aureus infection, and research suggests that there is a genetic basis for this variability. Several recent genome-wide association studies (GWAS) have identified variants that may affect susceptibility to infectious diseases, demonstrating the potential value of GWAS in this arena. METHODS: We conducted a GWAS to identify common variants associated with acquisition of S. aureus bacteremia (SAB) resulting from healthcare contact. We performed a logistic regression analysis to compare patients with healthcare contact who developed SAB (361 cases) to patients with healthcare contact in the same hospital who did not develop SAB (699 controls), testing 542,410 SNPs and adjusting for age (by decade), sex, and 6 significant principal components from our EIGENSTRAT analysis. Additionally, we evaluated the joint effect of the host and pathogen genomes in association with severity of SAB infection via logistic regression, including an interaction of host SNP with bacterial genotype, and adjusting for age (by decade), sex, the 6 significant principal components, and dialysis status. Bonferroni corrections were applied in both analyses to control for multiple comparisons. RESULTS: Ours is the first study that has attempted to evaluate the entire human genome for variants potentially involved in the acquisition or severity of SAB. Although this study identified no common variant of large effect size to have genome-wide significance for association with either the risk of acquiring SAB or severity of SAB, the variant (rs2043436) most significantly associated with severity of infection is located in a biologically plausible candidate gene (CDON, a member of the immunoglobulin family) and may warrant further study. CONCLUSIONS: The genetic architecture underlying SAB is likely to be complex. Future investigations using larger samples, narrowed phenotypes, and advances in both genotyping and analytical methodologies will be important tools for identifying causative variants for this common and serious cause of healthcare-associated infection.

Deficiency of asparagine synthetase causes congenital microcephaly and a progressive form of encephalopathy.

We analyzed four families that presented with a similar condition characterized by congenital microcephaly, intellectual disability, progressive cerebral atrophy, and intractable seizures. We show that recessive mutations in the ASNS gene are responsible for this syndrome. Two of the identified missense mutations dramatically reduce ASNS protein abundance, suggesting that the mutations cause loss of function. Hypomorphic Asns mutant mice have structural brain abnormalities, including enlarged ventricles and reduced cortical thickness, and show deficits in learning and memory mimicking aspects of the patient phenotype. ASNS encodes asparagine synthetase, which catalyzes the synthesis of asparagine from glutamine and aspartate. The neurological impairment resulting from ASNS deficiency may be explained by asparagine depletion in the brain or by accumulation of aspartate/glutamate leading to enhanced excitability and neuronal damage. Our study thus indicates that asparagine synthesis is essential for the development and function of the brain but not for that of other organs.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.

Mutation in TECPR2 reveals a role for autophagy in hereditary spastic paraparesis.

We studied five individuals from three Jewish Bukharian families affected by an apparently autosomal-recessive form of hereditary spastic paraparesis accompanied by severe intellectual disability, fluctuating central hypoventilation, gastresophageal reflux disease, wake apnea, areflexia, and unique dysmorphic features. Exome sequencing identified one homozygous variant shared among all affected individuals and absent in controls: a 1 bp frameshift TECPR2 deletion leading to a premature stop codon and predicting significant degradation of the protein. TECPR2 has been reported as a positive regulator of autophagy. We thus examined the autophagy-related fate of two key autophagic proteins, SQSTM1 (p62) and MAP1LC3B (LC3), in skin fibroblasts of an affected individual, as compared to a healthy control, and found that both protein levels were decreased and that there was a more pronounced decrease in the lipidated form of LC3 (LC3II). siRNA knockdown of TECPR2 showed similar changes, consistent with aberrant autophagy. Our results are strengthened by the fact that autophagy dysfunction has been implicated in a number of other neurodegenerative diseases. The discovered TECPR2 mutation implicates autophagy, a central intracellular mechanism, in spastic paraparesis.

Using ERDS to infer copy-number variants in high-coverage genomes.

Although there are many methods available for inferring copy-number variants (CNVs) from next-generation sequence data, there remains a need for a system that is computationally efficient but that retains good sensitivity and specificity across all types of CNVs. Here, we introduce a new method, estimation by read depth with single-nucleotide variants (ERDS), and use various approaches to compare its performance to other methods. We found that for common CNVs and high-coverage genomes, ERDS performs as well as the best method currently available (Genome STRiP), whereas for rare CNVs and high-coverage genomes, ERDS performs better than any available method. Importantly, ERDS accommodates both unique and highly amplified regions of the genome and does so without requiring separate alignments for calling CNVs and other variants. These comparisons show that for genomes sequenced at high coverage, ERDS provides a computationally convenient method that calls CNVs as well as or better than any currently available method.

Fine-mapping classical HLA variation associated with durable host control of HIV-1 infection in African Americans.

A small proportion of human immunodeficiency virus-1 (HIV-1) infected individuals, termed HIV-1 controllers, suppress viral replication to very low levels in the absence of therapy. Genetic investigations of this phenotype have strongly implicated variation in the class I major histocompatibility complex (MHC) region as key to HIV-1 control. We collected sequence-based classical class I HLA genotypes at 4-digit resolution in HIV-1-infected African American controllers and progressors (n = 1107), and tested them for association with host control using genome-wide single nucleotide polymorphism data to account for population structure. Several classical alleles at HLA-B were associated with host control, including B*57:03 [odds ratio (OR) = 5.1; P= 3.4 × 10(-18)] and B*81:01 (OR = 4.8; P= 1.3 × 10(-9)). Analysis of variable amino acid positions demonstrates that HLA-B position 97 is the most significant association with host control in African Americans (omnibus P = 1.2 × 10(-21)) and explains the signal of several HLA-B alleles, including B*57:03. Within HLA-B, we also identified independent effects at position 116 (omnibus P= 2.8 × 10(-15)) in the canonical F pocket, position 63 in the B pocket (P= 1.5 × 10(-3)) and the non-pocket position 245 (P= 8.8 × 10(-10)), which is thought to influence CD8-binding kinetics. Adjusting for these HLA-B effects, there is evidence for residual association in the MHC region. These results underscore the key role of HLA-B in affecting HIV-1 replication, likely through the molecular interaction between HLA-B and viral peptides presented by infected cells, and suggest that sites outside the peptide-binding pocket also influence HIV-1 control.

Copy number variation of KIR genes influences HIV-1 control.

A genome-wide screen for large structural variants showed that a copy number variant (CNV) in the region encoding killer cell immunoglobulin-like receptors (KIR) associates with HIV-1 control as measured by plasma viral load at set point in individuals of European ancestry. This CNV encompasses the KIR3DL1-KIR3DS1 locus, encoding receptors that interact with specific HLA-Bw4 molecules to regulate the activation of lymphocyte subsets including natural killer (NK) cells. We quantified the number of copies of KIR3DS1 and KIR3DL1 in a large HIV-1 positive cohort, and showed that an increase in KIR3DS1 count associates with a lower viral set point if its putative ligand is present (p = 0.00028), as does an increase in KIR3DL1 count in the presence of KIR3DS1 and appropriate ligands for both receptors (p = 0.0015). We further provide functional data that demonstrate that NK cells from individuals with multiple copies of KIR3DL1, in the presence of KIR3DS1 and the appropriate ligands, inhibit HIV-1 replication more robustly, and associated with a significant expansion in the frequency of KIR3DS1+, but not KIR3DL1+, NK cells in their peripheral blood. Our results suggest that the relative amounts of these activating and inhibitory KIR play a role in regulating the peripheral expansion of highly antiviral KIR3DS1+ NK cells, which may determine differences in HIV-1 control following infection.

SVA: software for annotating and visualizing sequenced human genomes.

SUMMARY: Here we present Sequence Variant Analyzer (SVA), a software tool that assigns a predicted biological function to variants identified in next-generation sequencing studies and provides a browser to visualize the variants in their genomic contexts. SVA also provides for flexible interaction with software implementing variant association tests allowing users to consider both the bioinformatic annotation of identified variants and the strength of their associations with studied traits. We illustrate the annotation features of SVA using two simple examples of sequenced genomes that harbor Mendelian mutations. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://www.svaproject.org.

The characterization of twenty sequenced human genomes.

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten "case" genomes from individuals with severe hemophilia A and ten "control" genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.

Host determinants of HIV-1 control in African Americans.

We performed a whole-genome association study of human immunodeficiency virus type 1 (HIV-1) set point among a cohort of African Americans (n = 515), and an intronic single-nucleotide polymorphism (SNP) in the HLA-B gene showed one of the strongest associations. We use a subset of patients to demonstrate that this SNP reflects the effect of the HLA-B*5703 allele, which shows a genome-wide statistically significant association with viral load set point (P = 5.6 x 10(-10)). These analyses therefore confirm a member of the HLA-B*57 group of alleles as the most important common variant that influences viral load variation in African Americans, which is consistent with what has been observed for individuals of European ancestry, among whom the most important common variant is HLA-B*5701.

A genomic score prognostic of outcome in trauma patients.

Traumatic injuries frequently lead to infection, organ failure, and death. Health care providers rely on several injury scoring systems to quantify the extent of injury and to help predict clinical outcome. Physiological, anatomical, and clinical laboratory analytic scoring systems (Acute Physiology and Chronic Health Evaluation [APACHE], Injury Severity Score [ISS]) are utilized, with limited success, to predict outcome following injury. The recent development of techniques for measuring the expression level of all of a person's genes simultaneously may make it possible to develop an injury scoring system based on the degree of gene activation. We hypothesized that a peripheral blood leukocyte gene expression score could predict outcome, including multiple organ failure, following severe blunt trauma. To test such a scoring system, we measured gene expression of peripheral blood leukocytes from patients within 12 h of traumatic injury. cRNA derived from whole blood leukocytes obtained within 12 h of injury provided gene expression data for the entire genome that were used to create a composite gene expression score for each patient. Total blood leukocytes were chosen because they are active during inflammation, which is reflective of poor outcome. The gene expression score combines the activation levels of all the genes into a single number which compares the patient's gene expression to the average gene expression in uninjured volunteers. Expression profiles from healthy volunteers were averaged to create a reference gene expression profile which was used to compute a difference from reference (DFR) score for each patient. This score described the overall genomic response of patients within the first 12 h following severe blunt trauma. Regression models were used to compare the association of the DFR, APACHE, and ISS scores with outcome. We hypothesized that patients with a total gene response more different from uninjured volunteers would tend to have poorer outcome than those more similar. Our data show that for measures of poor outcome, such as infections, organ failures, and length of hospital stay, this is correct. DFR scores were associated significantly with adverse outcome, including multiple organ failure, duration of ventilation, length of hospital stay, and infection rate. The association remained significant after adjustment for injury severity as measured by APACHE or ISS. A single score representing changes in gene expression in peripheral blood leukocytes within hours of severe blunt injury is associated with adverse clinical outcomes that develop later in the hospital course. Assessment of genome-wide gene expression provides useful clinical information that is different from that provided by currently utilized anatomic or physiologic scores.

Clinical utility of coronary CT angiography: coronary stenosis detection and prognosis in ambulatory patients.

INTRODUCTION: Multislice CT coronary angiography (MSCTA) accurately detects stenosis in patients undergoing coronary arteriography, but its accuracy in clinical outpatients is less certain. This study retrospectively analyzed MSCTA performance in a large outpatient cohort and examined 6-month clinical follow-up in those without invasive CA. METHODS: Patients underwent MSCTA for clinical indications including symptoms or noninvasive results being either equivocal or suspected as incorrect by referring clinicians. Standard 16-slice CT scanner techniques were used, and results were analyzed on the basis of both patient and vessel. Patients were treated medically or sent to invasive angiography on the basis of MSCTA results and judgment of referring clinicians. All invasive angiograms were analyzed using quantitative coronary angiography. Six-month clinical follow-up was determined in patients without CA. RESULTS: One thousand fifty-three consecutive patients were referred for MSCTA, resulting in 994 interpretable scans. Mean age was 58+/-13 years, 55% were male, 50% had prior noninvasive testing, and 90% had symptoms. Invasive angiography was performed in 160 patients, with significant stenoses present in 69%. MSCTA demonstrated 87% and 89% accuracy by patient- and vessel-based analysis, respectively, and was most accurate in the left main and right coronary arteries. Only two patients not referred for angiography had significant stenosis in those undergoing 6-month follow-up. CONCLUSIONS: MSCTA accurately detects obstructive coronary stenosis in clinical patients with possible cardiac symptoms, and effectively triages them for invasive angiography. Negative results are highly accurate in ruling out obstructive disease. Six-month prognosis is excellent in patients without significant disease determined by MSCT.