Representativeness of whole-genome sequencing approaches in England: the importance for understanding inequalities associated with SARS-CoV-2 infection.

Whole-genome sequencing (WGS) information has played a crucial role in the SARS-CoV-2 (COVID-19) pandemic by providing evidence about variants to inform public health policy. The purpose of this study was to assess the representativeness of sequenced cases compared with all COVID-19 cases in England, between March 2020 and August 2021, by demographic and socio-economic characteristics, to evaluate the representativeness and utility of these data in epidemiological analyses. To achieve this, polymerase chain reaction (PCR)-confirmed COVID-19 cases were extracted from the national laboratory system and linked with WGS data. During the study period, over 10% of COVID-19 cases in England had WGS data available for epidemiological analysis. With sequencing capacity increasing throughout the period, sequencing representativeness compared to all reported COVID-19 cases increased over time, allowing for valuable epidemiological analyses using demographic and socio-economic characteristics, particularly during periods with emerging novel SARS-CoV-2 variants. This study demonstrates the comprehensiveness of England's sequencing throughout the COVID-19 pandemic, rapidly detecting variants of concern, and enabling representative epidemiological analyses to inform policy.

Assessment of mortality and hospital admissions associated with confirmed infection with SARS-CoV-2 Alpha variant: a matched cohort and time-to-event analysis, England, October to December 2020.

BackgroundThe emergence of the SARS-CoV-2 Alpha variant in England coincided with a rapid increase in the number of PCR-confirmed COVID-19 cases in areas where the variant was concentrated.AimOur aim was to assess whether infection with Alpha was associated with more severe clinical outcomes than the wild type.MethodsLaboratory-confirmed infections with genomically sequenced SARS-CoV-2 Alpha and wild type between October and December 2020 were linked to routine healthcare and surveillance datasets. We conducted two statistical analyses to compare the risk of hospital admission and death within 28 days of testing between Alpha and wild-type infections: a matched cohort study and an adjusted Cox proportional hazards model. We assessed differences in disease severity by comparing hospital admission and mortality, including length of hospitalisation and time to death.ResultsOf 63,609 COVID-19 cases sequenced in England between October and December 2020, 6,038 had the Alpha variant. In the matched cohort analysis, we matched 2,821 cases with Alpha to 2,821 to cases with wild type. In the time-to-event analysis, we observed a 34% increased risk in hospitalisation associated with Alpha compared with wild type, but no significant difference in the risk of mortality.ConclusionWe found evidence of increased risk of hospitalisation after adjusting for key confounders, suggesting increased infection severity associated with the Alpha variant. Rapid assessments of the relative morbidity in terms of clinical outcomes and mortality associated with emerging SARS-CoV-2 variants compared with dominant variants are required to assess overall impact of SARS-CoV-2 mutations.

Genomic assessment of quarantine measures to prevent SARS-CoV-2 importation and transmission.

Mitigation of SARS-CoV-2 transmission from international travel is a priority. We evaluated the effectiveness of travellers being required to quarantine for 14-days on return to England in Summer 2020. We identified 4,207 travel-related SARS-CoV-2 cases and their contacts, and identified 827 associated SARS-CoV-2 genomes. Overall, quarantine was associated with a lower rate of contacts, and the impact of quarantine was greatest in the 16-20 age-group. 186 SARS-CoV-2 genomes were sufficiently unique to identify travel-related clusters. Fewer genomically-linked cases were observed for index cases who returned from countries with quarantine requirement compared to countries with no quarantine requirement. This difference was explained by fewer importation events per identified genome for these cases, as opposed to fewer onward contacts per case. Overall, our study demonstrates that a 14-day quarantine period reduces, but does not completely eliminate, the onward transmission of imported cases, mainly by dissuading travel to countries with a quarantine requirement.

Comparison of Next-Generation Sequencing Technologies for Comprehensive Assessment of Full-Length Hepatitis C Viral Genomes.

Affordable next-generation sequencing (NGS) technologies for hepatitis C virus (HCV) may potentially identify both viral genotype and resistance genetic motifs in the era of directly acting antiviral (DAA) therapies. This study compared the ability of high-throughput NGS methods to generate full-length, deep, HCV sequence data sets and evaluated their utility for diagnostics and clinical assessment. NGS methods using (i) unselected HCV RNA (metagenomics), (ii) preenrichment of HCV RNA by probe capture, and (iii) HCV preamplification by PCR implemented in four United Kingdom centers were compared. Metrics of sequence coverage and depth, quasispecies diversity, and detection of DAA resistance-associated variants (RAVs), mixed HCV genotypes, and other coinfections were compared using a panel of samples with different viral loads, genotypes, and mixed HCV genotypes/subtypes [geno(sub)types]. Each NGS method generated near-complete genome sequences from more than 90% of samples. Enrichment methods and PCR preamplification generated greater sequence depth and were more effective for samples with low viral loads. All NGS methodologies accurately identified mixed HCV genotype infections. Consensus sequences generated by different NGS methods were generally concordant, and majority RAVs were consistently detected. However, methods differed in their ability to detect minor populations of RAVs. Metagenomic methods identified human pegivirus coinfections. NGS provided a rapid, inexpensive method for generating whole HCV genomes to define infecting genotypes, RAVs, comprehensive viral strain analysis, and quasispecies diversity. Enrichment methods are particularly suited for high-throughput analysis while providing the genotype and information on potential DAA resistance.

Clinical Sequencing Exploratory Research Consortium: Accelerating Evidence-Based Practice of Genomic Medicine.

Despite rapid technical progress and demonstrable effectiveness for some types of diagnosis and therapy, much remains to be learned about clinical genome and exome sequencing (CGES) and its role within the practice of medicine. The Clinical Sequencing Exploratory Research (CSER) consortium includes 18 extramural research projects, one National Human Genome Research Institute (NHGRI) intramural project, and a coordinating center funded by the NHGRI and National Cancer Institute. The consortium is exploring analytic and clinical validity and utility, as well as the ethical, legal, and social implications of sequencing via multidisciplinary approaches; it has thus far recruited 5,577 participants across a spectrum of symptomatic and healthy children and adults by utilizing both germline and cancer sequencing. The CSER consortium is analyzing data and creating publically available procedures and tools related to participant preferences and consent, variant classification, disclosure and management of primary and secondary findings, health outcomes, and integration with electronic health records. Future research directions will refine measures of clinical utility of CGES in both germline and somatic testing, evaluate the use of CGES for screening in healthy individuals, explore the penetrance of pathogenic variants through extensive phenotyping, reduce discordances in public databases of genes and variants, examine social and ethnic disparities in the provision of genomics services, explore regulatory issues, and estimate the value and downstream costs of sequencing. The CSER consortium has established a shared community of research sites by using diverse approaches to pursue the evidence-based development of best practices in genomic medicine.

Assessment of the Utility of Whole Genome Sequencing of Measles Virus in the Characterisation of Outbreaks.

BACKGROUND: Measles is a highly infectious disease caused by measles virus (MeV). Despite the availability of a safe and cost-effective vaccine, measles is one of the world-leading causes of death in young children. Within Europe, there is a target for eliminating endemic measles in 2015, with molecular epidemiology required on 80% of cases for inclusion/exclusion of outbreak transmission chains. Currently, MeV is genotyped on the basis of a 450 nucleotide region of the nucleoprotein gene (N-450) and the hemagglutinin gene (H). However, this is not sufficiently informative for distinguishing endemic from imported MeV. We have developed an amplicon-based method for obtaining whole genome sequences (WGS) using NGS or Sanger methodologies from cell culture isolates or oral fluid specimens, and have sequenced over 60 samples, including 42 from the 2012 outbreak in the UK. RESULTS: Overall, NGS coverage was over 90% for approximately 71% of the samples tested. Analysis of 32 WGS excluding 3' and 5' termini (WGS-t) obtained from the outbreak indicates that the single nucleotide difference found between the two major groups of N-450 sequences detected during the outbreak is most likely a result of stochastic viral mutation during endemic transmission rather than of multiple importation events: earlier strains appear to have evolved into two distinct strain clusters in 2013, one containing strains with both outbreak-associated N-450 sequences. Additionally, phylogenetic analysis of each genomic region of MeV for the strains in this study suggests that the most information is acquired from the non-coding region located between the matrix and fusion protein genes (M/F NCR) and the N-450 genotyping sequence, an observation supported by entropy analysis across genotypes. CONCLUSIONS: We suggest that both M/F NCR and WGS-t could be used to complement the information from classical epidemiology and N-450 sequencing to address specific questions in the context of measles elimination.

Assessment of cortical and striatal involvement in 523 Huntington disease brains.

OBJECTIVE: To evaluate the relationship of striatal involvement in Huntington disease (HD) to involvement in other brain regions, CAG repeat size, onset age, and other factors. METHODS: We examined patterns of neuropathologic involvement in 664 HD brains submitted to the Harvard Brain Tissue Resource Center. Brains with concomitant Alzheimer or Parkinson changes (n = 82), more than 20% missing data (n = 46), incomplete sample submission (n = 12), or CAG repeat less than 36 (n = 1) were excluded, leaving 523 cases. Standardized ratings from 0 (absent) to 4 (severe) of gross and microscopic involvement were performed for 50 regions. Cluster analysis reduced the data to 2 main measures of involvement: striatal and cortical. RESULTS: The clusters were correlated with each other (r = 0.42) and with disease duration (striatal: r = 0.35; cortical: r = 0.31). The striatal cluster was correlated with HD repeat size (r = 0.50). The cortical cluster showed a stronger correlation with decreased brain weight (r = -0.52) than the striatal cluster (r = -0.33). The striatal cluster was correlated with younger death age (r = -0.31) and onset age (r = -0.46) while the cortical cluster was not (r = 0.09, r = -0.04, respectively). CONCLUSIONS: The 2 brain clusters had different relationships to the HD CAG repeat size, onset age, and brain weight, suggesting that neuropathologic involvement does not proceed in a strictly coupled fashion. The pattern and extent of involvement varies substantially from one brain to the next. These results suggest that regional involvement in HD brain is modified by factors which, if identified, may lend insight into novel routes to therapeutics.

Phylodynamic and phylogeographic patterns of the HIV type 1 subtype F1 parenteral epidemic in Romania.

In the late 1980s an HIV-1 epidemic emerged in Romania that was dominated by subtype F1. The main route of infection is believed to be parenteral transmission in children. We sequenced partial pol coding regions of 70 subtype F1 samples from children and adolescents from the PENTA-EPPICC network of which 67 were from Romania. Phylogenetic reconstruction using the sequences and other publically available global subtype F sequences showed that 79% of Romanian F1 sequences formed a statistically robust monophyletic cluster. The monophyletic cluster was epidemiologically linked to parenteral transmission in children. Coalescent-based analysis dated the origins of the parenteral epidemic to 1983 [1981-1987; 95% HPD]. The analysis also shows that the epidemic's effective population size has remained fairly constant since the early 1990s suggesting limited onward spread of the virus within the population. Furthermore, phylogeographic analysis suggests that the root location of the parenteral epidemic was Bucharest.

Injuries in youth football: national emergency department visits during 2001-2005 for young and adolescent players.

OBJECTIVES: Limited research exists describing youth football injuries, and many of these are confined to specific regions or communities. The authors describe U.S. pediatric football injury patterns receiving emergency department (ED) evaluation and compare injury patterns between the younger and older youth football participants. METHODS: A retrospective analysis of ED data on football injuries was performed using the National Electronic Injury Surveillance System-All Injury Program. Injury risk estimates were calculated over a 5-year period (2001-2005) using participation data from the National Sporting Goods Association. Injury types are described for young (7-11 years) and adolescent (12-17 years) male football participants. RESULTS: There were an estimated total of 1,060,823 visits to U.S. EDs for males with football-related injuries. The most common diagnoses in the younger group (7-11 years) were fracture/dislocation (29%), sprain/strain (27%), and contusion (27%). In the older group (ages 12-17 years), diagnoses included sprain/strain (31%), fracture/dislocation (29%), and contusion (23%). Older participants had a significantly higher injury risk of injury over the 5-year study period: 11.0 (95% confidence interval [CI] = 9.2 to 12.8) versus 6.1 (95% CI = 4.8 to 7.3) per 1,000 participants/year. Older participants had a higher injury risk across all categories, with the greatest disparity being with traumatic brain injury (TBI), 0.8 (95% CI = 0.6 to 1.0) versus 0.3 (95% CI = 0.2 to 0.4) per 1,000 participants/year. CONCLUSIONS: National youth football injury patterns are similar to those previously reported in community and cohort studies. Older participants have a significantly higher injury risk, especially with TBI.

Genome-wide admixture mapping for coronary artery calcification in African Americans: the NHLBI Family Heart Study.

Coronary artery calcification (CAC) is an important measure of subclinical coronary atherosclerosis and an independent predictor of coronary heart disease. To identify the genetic loci contributing to CAC, we conducted a genome-wide scan with 374 microsatellite markers by applying admixture mapping to 618 African American participants in the US National Heart, Lung, and Blood Institute Family Heart Study, in which 868 European American participants from family heart study and 157 Africans genotyped by the Marshfield Medical Genetics Center were used as the two reference founding populations for the African Americans, and a computer program based on a Markov Chain Monte Carlo algorithm, STRUCTURE 2.1, was used to estimate European and African ancestries among African Americans. A permutation test for random repeated sampling regression of CAC score on marker specific African ancestry found 22 markers statistically significant at the 0.05 level and four markers, D10S189 at 10p14, D20S159 at 20q13, D12S1294 at 12q14, and D6S1053 at 6q12, significant at the 0.01 level. D10S189 and D6S1053 were further confirmed at the 0.05 significance level by regression of CAC on allelic copy number, in which individual ancestry was used as a genetic background covariate to control possible stratification in African Americans. On the basis of the results from this and other independent studies, the location of D6S1053 at 80cM on chromosome 6 (6q12) seems to harbor a highly promising quantitative trait loci for atherosclerosis.

Assessment of automated genotyping protocols as tools for surveillance of HIV-1 genetic diversity.

BACKGROUND: The routine use of drug resistance testing provides an abundant source of HIV-1 sequence data. However, it is not clear how reliable standard genotyping of these sequences is for describing HIV-1 genetic variation and for detecting novel genetic variants and epidemiological trends. OBJECTIVES: To compare assignment of HIV-1 resistance test sequences to reference strains across commonly used genotyping protocols. METHODS: Subtype assignments were compared across three standard genotyping protocols for 10 537 resistance test sequences, representing approximately one-fifth of all reported infections in the United Kingdom. Sequences that were inconsistently genotyped across methods, or that were unassigned by at least one method, were examined for evidence of recombination using sliding-window-based approaches. RESULTS: Although agreement across methods was high for subtypes B, C and H, it was generally much lower (< 50%) for other subtypes. Disagreement between methods typically involved closely related, but epidemiologically distinct, groups or involved a significant proportion ( approximately 12%) of divergent sequences in which analysis revealed widespread evidence of recombination and a remarkable diversity of unusual recombinant forms. CONCLUSIONS: With frequent long-distance transfer of viral strains and widespread recombination between them, genetic and epidemiological relationships within HIV-1 are becoming increasingly complex. Current methods of subtype assignment vary in their ability to identify novel genetic variants and to distinguish epidemiologically distinct strains. Capturing meaningful epidemiological information from resistance test data will require a critical understanding of the methodologies used in order to appreciate the possible sources of error and misclassification.

Quality assessment of the human genome sequence.

As the final sequencing of the human genome has now been completed, we present the results of the largest examination of the quality of the finished DNA sequence. The completed study covers the major contributing sequencing centres and is based on a rigorous combination of laboratory experiments and computational analysis.

The analysis of survival data with a non-susceptible fraction and dual censoring mechanisms.

It is known that the ages of onset of many diseases are determined by both a genetic predisposition to disease as well as environmental risk factors that are capable of either triggering or hastening the onset of disease. Difficulties in modelling onset ages arise when a large fraction fail to inherit the disease-causing gene, and multiple reasons for censoring result in unobserved onset ages. We present a parametric Bayesian model that includes subjects with missing age information, non-susceptible subjects and allows for regression on risk factor information. The model is fit using Markov chain Monte Carlo simulation from the posterior distribution, and allows the simultaneous estimation of the proportion of the population at risk of disease, the mean onset age of disease, survival after disease onset, and the association of risk factors with susceptibility, onset age and survival after onset. An example employing Huntington's disease data is presented.

Genetic variability of adult body mass index: a longitudinal assessment in framingham families.

OBJECTIVE: To explore the contribution of genetics to the mean, SD, maximum value, maximum less the mean, and change over time in body mass index (BMI) and the residual of body weight after adjustment for height. BMI is frequently used as a general indicator of obesity because of its ease and reliability in ascertainment. Cross-sectional twin and family studies have shown a moderate-to-substantial genetic component for BMI. However, the contribution of genetics to the long-term average, variability, or change over time in BMI is less clear. RESEARCH METHODS AND PROCEDURES: Longitudinal data from the Framingham heart study were used to create pedigrees of age-matched individuals. Heritability estimates were derived using variance-decomposition methods on a total of 1051 individuals from 380 extended pedigrees followed for a period of 20 years. All subjects were followed from approximately age 35 to 55 years. RESULTS: Moderate heritability estimates were found for the mean BMI (h(2) = 0.37), maximum BMI (h(2) = 0.40), and the mean residual of body weight (h(2) = 0.36). Low heritability estimates (h(2) congruent with 0.20) were found for the maximum less the mean in BMI and the SDs of BMI and residual of body weight. No additive genetic contribution was found for the average change over time in BMI or the residual of body weight. DISCUSSION: These findings suggest that there is a significant genetic component for the magnitude of BMI throughout an individual's middle-adult years; however, little evidence was found for a genetic contribution to the variability or rate of change in an individual's BMI.