Genome-wide association analysis of COVID-19 mortality risk in SARS-CoV-2 genomes identifies mutation in the SARS-CoV-2 spike protein that colocalizes with P.1 of the Brazilian strain.

SARS-CoV-2 mortality has been extensively studied in relation to host susceptibility. How sequence variations in the SARS-CoV-2 genome affect pathogenicity is poorly understood. Starting in October 2020, using the methodology of genome-wide association studies (GWAS), we looked at the association between whole-genome sequencing (WGS) data of the virus and COVID-19 mortality as a potential method of early identification of highly pathogenic strains to target for containment. Although continuously updating our analysis, in December 2020, we analyzed 7548 single-stranded SARS-CoV-2 genomes of COVID-19 patients in the GISAID database and associated variants with mortality using a logistic regression. In total, evaluating 29,891 sequenced loci of the viral genome for association with patient/host mortality, two loci, at 12,053 and 25,088 bp, achieved genome-wide significance (p values of 4.09e-09 and 4.41e-23, respectively), though only 25,088 bp remained significant in follow-up analyses. Our association findings were exclusively driven by the samples that were submitted from Brazil (p value of 4.90e-13 for 25,088 bp). The mutation frequency of 25,088 bp in the Brazilian samples on GISAID has rapidly increased from about 0.4 in October/December 2020 to 0.77 in March 2021. Although GWAS methodology is suitable for samples in which mutation frequencies varies between geographical regions, it cannot account for mutation frequencies that change rapidly overtime, rendering a GWAS follow-up analysis of the GISAID samples that have been submitted after December 2020 as invalid. The locus at 25,088 bp is located in the P.1 strain, which later (April 2021) became one of the distinguishing loci (precisely, substitution V1176F) of the Brazilian strain as defined by the Centers for Disease Control. Specifically, the mutations at 25,088 bp occur in the S2 subunit of the SARS-CoV-2 spike protein, which plays a key role in viral entry of target host cells. Since the mutations alter amino acid coding sequences, they potentially imposing structural changes that could enhance viral infectivity and symptom severity. Our analysis suggests that GWAS methodology can provide suitable analysis tools for the real-time detection of new more transmissible and pathogenic viral strains in databases such as GISAID, though new approaches are needed to accommodate rapidly changing mutation frequencies over time, in the presence of simultaneously changing case/control ratios. Improvements of the associated metadata/patient information in terms of quality and availability will also be important to fully utilize the potential of GWAS methodology in this field.

Ozone exposure, vitamin C intake, and genetic susceptibility of asthmatic children in Mexico City: a cohort study.

BACKGROUND: We previously reported that asthmatic children with GSTM1 null genotype may be more susceptible to the acute effect of ozone on the small airways and might benefit from antioxidant supplementation. This study aims to assess the acute effect of ozone on lung function (FEF(25-75)) in asthmatic children according to dietary intake of vitamin C and the number of putative risk alleles in three antioxidant genes: GSTM1, GSTP1 (rs1695), and NQO1 (rs1800566). METHODS: 257 asthmatic children from two cohort studies conducted in Mexico City were included. Stratified linear mixed models with random intercepts and random slopes on ozone were used. Potential confounding by ethnicity was assessed. Analyses were conducted under single gene and genotype score approaches. RESULTS: The change in FEF(25-75) per interquartile range (60 ppb) of ozone in persistent asthmatic children with low vitamin C intake and GSTM1 null was -91.2 ml/s (p = 0.06). Persistent asthmatic children with 4 to 6 risk alleles and low vitamin C intake showed an average decrement in FEF(25-75) of 97.2 ml/s per 60 ppb of ozone (p = 0.03). In contrast in children with 1 to 3 risk alleles, acute effects of ozone on FEF25-75 did not differ by vitamin C intake. CONCLUSIONS: Our results provide further evidence that asthmatic children predicted to have compromised antioxidant defense by virtue of genetic susceptibility combined with deficient antioxidant intake may be at increased risk of adverse effects of ozone on pulmonary function.

Incorporating missingness for estimation of marginal regression models with multiple source predictors.

Multiple informant data refers to information obtained from different individuals or sources used to measure the same construct; for example, researchers might collect information regarding child psychopathology from the child's teacher and the child's parent. Frequently, studies with multiple informants have incomplete observations; in some cases the missingness of informants is substantial. We introduce a Maximum Likelihood (ML) technique to fit models with multiple informants as predictors that permits missingness in the predictors as well as the response. We provide closed form solutions when possible and analytically compare the ML technique to the existing Generalized Estimating Equations (GEE) approach. We demonstrate that the ML approach can be used to compare the effect of the informants on response without standardizing the data. Simulations incorporating missingness show that ML is more efficient than the existing GEE method. In the presence of MCAR missing data, we find through a simulation study that the ML approach is robust to a relatively extreme departure from the normality assumption. We implement both methods in a study investigating the association between physical activity and obesity with activity measured using multiple informants (children and their mothers).