Cohort profile: the Johns Hopkins COVID Long Study (JHCLS), a United States Nationwide Prospective Cohort Study

Purpose

COVID-19 disease continues to affect millions of individuals worldwide, both in the short and long term. The post-acute complications of SARS-CoV-2 infection, referred to as long COVID, result in diverse symptoms affecting multiple organ systems. Little is known regarding how the symptoms associated with long COVID progress and resolve over time. The Johns Hopkins COVID Long Study aims to prospectively examine the short- and long-term consequences of COVID-19 disease in individuals both with and without a history of SARS-CoV-2 infection using self-reported data collected in an online survey.

Participants

Sixteen thousand, seven hundred sixty-four adults with a history of SARS-CoV-2 infection and 799 adults without a history of SARS-CoV-2 infection who completed an online baseline survey.

Findings to date

This cohort profile describes the baseline characteristics of the Johns Hopkins COVID Long Study. Among 16,764 participants with a history of SARS-CoV-2 infection and defined long COVID status, 75% reported a good or excellent health status prior to infection, 99% reported experiencing at least one COVID-19 symptom during the acute phase of infection, 9.9% reported a hospitalization, and 63% were defined as having long COVID using the WHO definition.

Future plans

Analysis of longitudinal data will be used to investigate the progression and resolution of long COVID symptoms over time.

A hepatitis B virus (HBV) sequence variation graph improves alignment and sample-specific consensus sequence construction

Nearly 300 million individuals live with chronic hepatitis B virus (HBV) infection (CHB), for which no curative therapy is available. As viral diversity is associated with pathogenesis and immunological control of infection, improved methods to characterize this diversity could aid drug development efforts. Conventionally, viral sequencing data are mapped/aligned to a reference genome, and only the aligned sequences are retained for analysis. Thus, reference selection is critical, yet selecting the most representative reference a priori remains difficult. We investigate an alternative pangenome approach which can combine multiple reference sequences into a graph which can be used during alignment. Using simulated short-read sequencing data generated from publicly available HBV genomes and real sequencing data from an individual living with CHB, we demonstrate alignment to a phylogenetically representative 'genome graph' can improve alignment, avoid issues of reference ambiguity, and facilitate the construction of sample-specific consensus sequences more genetically similar to the individual's infection. Graph-based methods can, therefore, improve efforts to characterize the genetics of viral pathogens, including HBV, and have broader implications in host-pathogen research.

Genetic Susceptibility to Astrovirus Diarrhea in Bangladeshi Infants

Background

Astroviral infections commonly cause acute nonbacterial gastroenteritis in children globally. However, these infections often go undiagnosed outside of research settings. There is no treatment available for astrovirus, and Astroviridae strain diversity presents a challenge to potential vaccine development.

Methods

To address our hypothesis that host genetic risk factors are associated with astrovirus disease susceptibility, we performed a genome-wide association study of astrovirus infection in the first year of life from children enrolled in 2 Bangladeshi birth cohorts.

Results

We identified a novel region on chromosome 1 near the loricrin gene (LOR) associated with astrovirus diarrheal infection (rs75437404; meta-analysis P = 8.82 × 10-9; A allele odds ratio, 2.71) and on chromosome 10 near the prolactin releasing hormone receptor gene (PRLHR) (rs75935441; meta-analysis P = 1.33 × 10-8; C allele odds ratio, 4.17). The prolactin-releasing peptide has been shown to influence feeding patterns and energy balance in mice. In addition, several single-nucleotide polymorphisms in the chromosome 1 locus have previously been associated with expression of innate immune system genes PGLYRP4, S100A9, and S100A12.

Conclusions

This study identified 2 significant host genetic regions that may influence astrovirus diarrhea susceptibility and should be considered in further studies.

Genome-Wide Association Studies of Diarrhea Frequency and Duration in the First Year of Life in Bangladeshi Infants

BACKGROUND

Diarrhea is the second leading cause of death in children under 5 years old worldwide. Known diarrhea risk factors include sanitation, water sources, and pathogens but do not fully explain the heterogeneity in frequency and duration of diarrhea in young children. We evaluated the role of host genetics in diarrhea.

METHODS

Using 3 well-characterized birth cohorts from an impoverished area of Dhaka, Bangladesh, we compared infants with no diarrhea in the first year of life to those with an abundance, measured by either frequency or duration. We performed a genome-wide association analysis for each cohort under an additive model and then meta-analyzed across the studies.

RESULTS

For diarrhea frequency, we identified 2 genome-wide significant loci associated with not having any diarrhea, on chromosome 21 within the noncoding RNA AP000959 (C allele odds ratio [OR] = 0.31, P = 4.01 × 10-8), and on chromosome 8 within SAMD12 (T allele OR = 0.35, P = 4.74 × 10-7). For duration of diarrhea, we identified 2 loci associated with no diarrhea, including the same locus on chromosome 21 (C allele OR = 0.31, P = 1.59 × 10-8) and another locus on chromosome 17 near WSCD1 (C allele OR = 0.35, P = 1.09 × 10-7).

CONCLUSIONS

These loci are in or near genes involved in enteric nervous system development and intestinal inflammation and may be potential targets for diarrhea therapeutics.

Physical Health and Mental Fatigue Disability Associated with Long COVID: Baseline Results from a US Nationwide Cohort

BACKGROUND

Persistent symptoms after severe acute respiratory disease coronavirus 2 (SARS-COV-2; long COVID) occur in 10%-55% of individuals, but the impact on daily functioning and disability remains unquantified.

METHODS

To characterize disability associated with long COVID, we analyzed baseline data from an online, US-based cohort study. Adult participants included those reporting a history of COVID-19 (n = 8874) or never having COVID-19 (n = 633) without prior disability. The main outcomes were self-reported physical mobility, instrumental activities of daily living (IADL), and mental fatigue disability, assessed by measuring 5 disability components: difficulty walking a quarter mile or climbing 10 stairs (mobility), difficulty doing light or heavy housework (IADL), and Wood Mental Fatigue Inventory score (mental fatigue).

RESULTS

Of 7926 participants with long COVID, 65% were classified with at least one disability, as compared with 6% and 14% for resolved COVID and no COVID, respectively. Additionally, 22% were classified as disabled in all 3 categories. Age, prior comorbidity, increased body mass index, female sex, COVID-19 hospitalization, non-white/multi-race were associated with higher disability burden. Dizziness and heavy limbs at infection were associated with disability regardless of hospitalization. Dyspnea and tremors were associated with disability in non-hospitalized individuals. Vaccination was protective against disability.

CONCLUSIONS

We observed a high burden of new disability associated with long COVID, which has serious implications for individual and societal health. Longitudinal evaluation of COVID-19 patients is necessary to identify patterns of recovery and treatment options.

Erratum: Accelerated aging with HIV begins at the time of initial HIV infection

[This corrects the article DOI: 10.1016/j.isci.2022.104488.].

Targeting hepatitis B vaccine escape using immunogenetics in Bangladeshi infants

Hepatitis B virus (HBV) vaccine escape mutants (VEM) are increasingly described, threatening progress in control of this virus worldwide. Here we studied the relationship between host genetic variation, vaccine immunogenicity and viral sequences implicating VEM emergence. In a cohort of 1,096 Bangladeshi children, we identified human leukocyte antigen (HLA) variants associated with response vaccine antigens. Using an HLA imputation panel with 9,448 south Asian individuals DPB1*04:01 was associated with higher HBV antibody responses (p=4.5×10-30). The underlying mechanism is a result of higher affinity binding of HBV surface antigen epitopes to DPB1*04:01 dimers. This is likely a result of evolutionary pressure at the HBV surface antigen 'a-determinant' segment incurring VEM specific to HBV. Prioritizing pre-S isoform HBV vaccines may tackle the rise of HBV vaccine evasion.

Prospective Cohort Study of Cryptosporidium Infection and Shedding in Infants and Their Households

BACKGROUND

Cryptosporidium spp. are responsible for significant diarrheal morbidity and mortality in under-5 children. There is no vaccine; thus, a focus on prevention is paramount. Prior studies suggest that person-to-person spread may be an important pathway for transmission to young children. Here we describe a longitudinal cohort study of 100 families with infants to determine rates of cryptosporidiosis within households during the coronavirus disease 2019 (COVID-19) pandemic.

METHODS

Families living in Mirpur, Bangladesh, with 1 infant aged 6-8 months were enrolled and followed with weekly illness survey and stool testing for Cryptosporidium for 8 months.

RESULTS

From December 2020 to August 2021, 100 families were enrolled. Forty-four percent of index children and 35% of siblings had at least 1 Cryptosporidium infection. Shedding of Cryptosporidium occurred for a mean (standard deviation) of 19 (8.3) days in index infants, 16.1 (11.6) days in children 1-5 years, and 16.2 (12.8) days in adults. A longer duration of Cryptosporidium shedding was associated with growth faltering in infants. There was a spike in Cryptosporidium cases in May 2021, which coincided with a spike in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases in the region.

CONCLUSIONS

In this intensive, longitudinal study of Cryptosporidium infection in families we found high rates of cryptosporidiosis in infants and children, and prolonged parasite shedding, especially among malnourished children. These data support that transmission within the household is an important route of exposure for young infants and that treatment of nondiarrheal infection to interrupt person-to-person transmission within the home may be essential for preventing cryptosporidiosis in infants.

Characterizing common and rare variations in non-traditional glycemic biomarkers using multivariate approaches on multi-ancestry ARIC study

Glycated hemoglobin, fasting glucose, glycated albumin, and fructosamine are biomarkers that reflect different aspects of the glycemic process. Genetic studies of these glycemic biomarkers can shed light on unknown aspects of type 2 diabetes genetics and biology. While there exists several GWAS of glycated hemoglobin and fasting glucose, very few GWAS have focused on glycated albumin or fructosamine. We performed a multi-phenotype GWAS of glycated albumin and fructosamine from 7,395 White and 2,016 Black participants in the Atherosclerosis Risk in Communities (ARIC) study on the common variants from genotyped/imputed data. We found 2 genome-wide significant loci, one mapping to known type 2 diabetes gene (ARAP1/STARD10, p = 2.8 × 10-8) and another mapping to a novel gene (UGT1A, p = 1.4 × 10-8) using multi-omics gene mapping strategies in diabetes-relevant tissues. We identified additional loci that were ancestry-specific (e.g., PRKCA from African ancestry individuals, p = 1.7 × 10-8) and sex-specific (TEX29 locus in males only, p = 3.0 × 10-8). Further, we implemented multi-phenotype gene-burden tests on whole-exome sequence data from 6,590 White and 2,309 Black ARIC participants. Eleven genes across different rare variant aggregation strategies were exome-wide significant only in multi-ancestry analysis. Four out of 11 genes had notable enrichment of rare predicted loss of function variants in African ancestry participants despite smaller sample size. Overall, 8 out of 15 loci/genes were implicated to influence these biomarkers via glycemic pathways. This study illustrates improved locus discovery and potential effector gene discovery by leveraging joint patterns of related biomarkers across entire allele frequency spectrum in multi-ancestry analyses. Most of the loci/genes we identified have not been previously implicated in studies of type 2 diabetes, and future investigation of the loci/genes potentially acting through glycemic pathways may help us better understand risk of developing type 2 diabetes.

Defining representativeness of study samples in medical and population health research

Medical and population health science researchers frequently make ambiguous statements about whether they believe their study sample or results are representative of some (implicit or explicit) target population. This article provides a comprehensive definition of representativeness, with the goal of capturing the different ways in which a study can be representative of a target population. It is proposed that a study is representative if the estimate obtained in the study sample is generalisable to the target population (owing to representative sampling, estimation of stratum specific effects, or quantitative methods to generalise or transport estimates) or the interpretation of the results is generalisable to the target population (based on fundamental scientific premises and substantive background knowledge). This definition is explored in the context of four covid-19 studies, ranging from laboratory science to descriptive epidemiology. All statements regarding representativeness should make clear the way in which the study results generalise, the target population the results are being generalised to, and the assumptions that must hold for that generalisation to be scientifically or statistically justifiable.

Multi-ancestry sex-stratified genomic associations with HIV viral load and controller status from the ICGH

Biological sex and host genetics influence HIV pathogenesis. Females have a higher likelihood of spontaneous viral control and lower setpoint viral load (spVL). No prior studies have assessed sex-specific genetics of HIV. To address this, we performed a sex stratified genome-wide association study using data from the International Collaboration for the Genomics of HIV. Although it is the largest collection of genomic data in HIV, this multi-ethnic sample of 9,705 people is 81.3% male. We sought to identify sex-specific genetic variants and genes associated with HIV spVL and control. We confirmed associations in the HLA and CCR5 regions in males, and HLA in females. Gene-based analyses detected associations between HIV spVL and PET100 (Pvalue=8.36x10-07), PCP2 (Pvalue=8.81x10-07), XAB2 (Pvalue=1.32x10-6) and STXBP2 (Pvalue=1.65x10-4) only in males. We detected variants with a significant sex-differential effect on spVL in SDC3 and PUM1 (rs10914268,Pvalue=1.93x10-08) and PSORS1C2 (rs1265159, Pvalue=3.26x10-08) and on HIV control in SUB1 (rs687659, Pvalue=1.02×10-08), AL158151.3, PTPA and IER5L (rs4387067, Pvalue=2.07×10-09). Those variants have epigenetic and genetic interactions with relevant genes with both cis and trans effects. In summary, we identified sex-shared associations at the single variant level, sex-specific associations at the gene-based level, and genetic variants with significant differential effects between the sexes.

No Association of IFNL4 Genotype With Opportunistic Infections and Cancers Among Men With Human Immunodeficiency Virus 1 Infection

BACKGROUND

IFNL4 genetic variants that are strongly associated with clearance of hepatitis C virus have been linked to risk of certain opportunistic infections (OIs) and cancers, including Kaposi sarcoma, cytomegalovirus infection, and herpes simplex virus infection. As the interferon (IFN) λ family plays a role in response to viral, bacterial, and fungal infections, IFNL4 genotype might affect risk for a wide range of OIs/cancers.

METHODS

We examined associations between genotype for the functional IFNL4 rs368234815 polymorphism and incidence of 16 OIs/cancers among 2310 men with human immunodeficiency virus (2038 white; 272 black) enrolled in the Multicenter AIDS Cohort Study during 1984-1990. Our primary analyses used Cox proportional hazards models adjusted for self-reported racial ancestry to estimate hazard ratios with 95% confidence intervals, comparing participants with the genotypes that generate IFN-λ4 and those with the genotype that abrogates IFN-λ4. We censored follow-up at the introduction of highly effective antiretroviral therapies.

RESULTS

We found no statistically significant association between IFNL4 genotype and the incidence of Kaposi sarcoma (hazard ratio, 0.92 [95% confidence interval, .76-1.11]), cytomegalovirus infection (0.94 [.71-1.24]), herpes simplex virus infection (1.37 [.68-2.93]), or any other OI/cancer. We observed consistent results using additive genetic models and after controlling for CD4 cell count through time-dependent adjustment or restriction to participants with a low CD4 cell count.

CONCLUSIONS

The absence of associations between IFNL4 genotype and these OIs/cancers provides evidence that this gene does not affect the risk of disease from opportunistic pathogens.

Pathogen exposure misclassification can bias association signals in GWAS of infectious diseases when using population-based common control subjects

Genome-wide association studies (GWASs) have been performed to identify host genetic factors for a range of phenotypes, including for infectious diseases. The use of population-based common control subjects from biobanks and extensive consortia is a valuable resource to increase sample sizes in the identification of associated loci with minimal additional expense. Non-differential misclassification of the outcome has been reported when the control subjects are not well characterized, which often attenuates the true effect size. However, for infectious diseases the comparison of affected subjects to population-based common control subjects regardless of pathogen exposure can also result in selection bias. Through simulated comparisons of pathogen-exposed cases and population-based common control subjects, we demonstrate that not accounting for pathogen exposure can result in biased effect estimates and spurious genome-wide significant signals. Further, the observed association can be distorted depending upon strength of the association between a locus and pathogen exposure and the prevalence of pathogen exposure. We also used a real data example from the hepatitis C virus (HCV) genetic consortium comparing HCV spontaneous clearance to persistent infection with both well-characterized control subjects and population-based common control subjects from the UK Biobank. We find biased effect estimates for known HCV clearance-associated loci and potentially spurious HCV clearance associations. These findings suggest that the choice of control subjects is especially important for infectious diseases or outcomes that are conditional upon environmental exposures.

A hepatitis B virus (HBV) sequence variation graph improves sequence alignment and sample-specific consensus sequence construction for genetic analysis of HBV

Hepatitis B virus (HBV) remains a global public health concern, with over 250 million individuals living with chronic HBV infection (CHB) and no curative therapy currently available. Viral diversity is associated with CHB pathogenesis and immunological control of infection. Improved methods to characterize the viral genome at both the population and intra-host level could aid drug development efforts. Conventionally, HBV sequencing data are aligned to a linear reference genome and only sequences capable of aligning to the reference are captured for analysis. Reference selection has additional consequences, including sample-specific 'consensus' sequence construction. It remains unclear how to select a reference from available sequences and whether a single reference is sufficient for genetic analyses. Using simulated short-read sequencing data generated from full-length publicly available HBV genome sequences and HBV sequencing data from a longitudinally sampled individual with CHB, we investigate alternative graph-based alignment approaches. We demonstrate that using a phylogenetically representative 'genome graph' for alignment, rather than linear reference sequences, avoids issues of reference ambiguity, improves alignment, and facilitates the construction of sample-specific consensus sequences genetically similar to an individual's infection. Graph-based methods can therefore improve efforts to characterize the genetics of viral pathogens, including HBV, and may have broad implications in host pathogen research.

Identity by descent mapping of HCV spontaneous clearance in populations of diverse ancestry

Background

Acute infection with hepatitis C virus (HCV) affects millions of individuals worldwide. Host genetics plays a role in spontaneous clearance of the acute infection which occurs in approximately 30% of the individuals. Common variants in GPR158, genes in the interferon lambda (IFNL) cluster, and the MHC region have been associated with HCV clearance in populations of diverse ancestry. Fine mapping of those regions has identified some key variants and amino acids as potential causal variants but the role of rare variants in those regions and in the genome, in general, has not been explored. We aimed to detect haplotypes containing rare variants related to HCV clearance using identity-by-descent (IBD) haplotype sharing between unrelated cases/case pairs and case/controls pairs in 3,608 individuals with European and African ancestry.

Results

We detected 1,711,832 and 5,678,043 and individual pairs of IBD segments in the European and African ancestry individuals, respectively. As expected, individuals of African descent had more, and shorter segments compared to Europeans. We did not detect any significant IBD signals in the known associated gene regions.

Conclusions

IBD is based on sharing of haplotypes and is most powerful in populations with a shared founder or recent common ancestor. For the complex trait of HCV clearance, we used two outbred, global populations that limited our power to detect IBD associations. Overall, in this population-based sample we failed to detect rare variations associated with HCV clearance in individuals of European and African ancestry.

Rare variant analyses across multiethnic cohorts identify novel genes for refractive error

Refractive error, measured here as mean spherical equivalent (SER), is a complex eye condition caused by both genetic and environmental factors. Individuals with strong positive or negative values of SER require spectacles or other approaches for vision correction. Common genetic risk factors have been identified by genome-wide association studies (GWAS), but a great part of the refractive error heritability is still missing. Some of this heritability may be explained by rare variants (minor allele frequency [MAF] ≤ 0.01.). We performed multiple gene-based association tests of mean Spherical Equivalent with rare variants in exome array data from the Consortium for Refractive Error and Myopia (CREAM). The dataset consisted of over 27,000 total subjects from five cohorts of Indo-European and Eastern Asian ethnicity. We identified 129 unique genes associated with refractive error, many of which were replicated in multiple cohorts. Our best novel candidates included the retina expressed PDCD6IP, the circadian rhythm gene PER3, and P4HTM, which affects eye morphology. Future work will include functional studies and validation. Identification of genes contributing to refractive error and future understanding of their function may lead to better treatment and prevention of refractive errors, which themselves are important risk factors for various blinding conditions.

Physical and mental health disability associated with long-COVID: Baseline results from a US nationwide cohort

Importance

Persistent symptoms after SARS-COV-2 infection, or long-COVID, may occur in anywhere from 10-55% of those who have had COVID-19, but the extent of impact on daily functioning and disability remains unquantified.

Objective

To characterize physical and mental disability associated with long-COVID.

Design

Cross-sectional analysis of baseline data from a cohort study.

Setting

Online US nationwide survey.

Participants

Adults 18 years of age and older who live in the US who either report a history of COVID-19 illness (n=8,874) or report never having had COVID-19 (n=633).

Main Outcome and Measures

Self-reported mobility disability (difficulty walking a quarter of a mile and/or up 10 stairs, instrumental activities of daily living [IADL] disability (difficulty doing light or heavy housework), and mental fatigue as measured by the Wood Mental Fatigue Inventory (WMFI).

Results

Of 7,926 participants with long-COVID, the median age was 45 years, 84% were female, 89% self-reported white race, and 7.4% self-reported Hispanic/Latino ethnicity. Sixty-five percent of long-COVID participants were classified as having at least one disability, compared to 6% of those with resolved-COVID (n=948) and 14% of those with no-COVID (n=633). Of long-COVID participants, about 1% and 5% were classified as critically physically disabled or mentally fatigued, respectively. Age, prior comorbidity, increased BMI, female gender, hospitalization for COVID-19, non-white race, and multi-race were all associated with significantly higher disability burden. Dizziness at the time of infection (33% non-hospitalized, 39% hospitalized) was associated with all five disability components in both hospitalized and non-hospitalized groups. Heavy limbs, dyspnea, and tremors were associated with four of the five components of disability in the non-hospitalized group, and heavy limbs was associated with four of the five components in the hospitalized group. Vaccination was protective against development of disability.

Conclusion and Relevance

We observed a high burden of physical and mental disability associated with long-COVID which has serious implications for individual and societal health that may be partially mitigated by vaccination. Longitudinal characterization and evaluation of COVID-19 patients is necessary to identify patterns of recovery and treatment options.

Deep learning and social network analysis elucidate drivers of HIV transmission in a high-incidence cohort of people who inject drugs

Globally, people who inject drugs (PWID) experience some of the fastest-growing HIV epidemics. Network-based approaches represent a powerful tool for understanding and combating these epidemics; however, detailed social network studies are limited and pose analytical challenges. We collected longitudinal social (injection partners) and spatial (injection venues) network information from 2512 PWID in New Delhi, India. We leveraged network analysis and graph neural networks (GNNs) to uncover factors associated with HIV transmission and identify optimal intervention delivery points. Longitudinal HIV incidence was 21.3 per 100 person-years. Overlapping community detection using GNNs revealed seven communities, with HIV incidence concentrated within one community. The injection venue most strongly associated with incidence was found to overlap six of the seven communities, suggesting that an intervention deployed at this one location could reach the majority of the sample. These findings highlight the utility of network analysis and deep learning in HIV program design.

Benchmarking statistical methods for analyzing parent-child dyads in genetic association studies

Genetic association studies of child health outcomes often employ family-based study designs. One of the most popular family-based designs is the case-parent trio design that considers the smallest possible nuclear family consisting of two parents and their affected child. This trio design is particularly advantageous for studying relatively rare disorders because it is less prone to type 1 error inflation due to population stratification compared to population-based study designs (e.g., case-control studies). However, obtaining genetic data from both parents is difficult, from a practical perspective, and many large studies predominantly measure genetic variants in mother-child dyads. While some statistical methods for analyzing parent-child dyad data (most commonly involving mother-child pairs) exist, it is not clear if they provide the same advantage as trio methods in protecting against population stratification, or if a specific dyad design (e.g., case-mother dyads vs. case-mother/control-mother dyads) is more advantageous. In this article, we review existing statistical methods for analyzing genome-wide marker data on dyads and perform extensive simulation experiments to benchmark their type I errors and statistical power under different scenarios. We extend our evaluation to existing methods for analyzing a combination of case-parent trios and dyads together. We apply these methods on genotyped and imputed data from multiethnic mother-child pairs only, case-parent trios only or combinations of both dyads and trios from the Gene, Environment Association Studies consortium (GENEVA), where each family was ascertained through a child affected by nonsyndromic cleft lip with or without cleft palate. Results from the GENEVA study corroborate the findings from our simulation experiments. Finally, we provide recommendations for using statistical genetic association methods for dyads.

Accelerated aging with HIV begins at the time of initial HIV infection

Living with HIV infection is associated with early onset of aging-related chronic conditions, sometimes described as accelerated aging. Epigenetic DNA methylation patterns can evaluate acceleration of biological age relative to chronological age. The impact of initial HIV infection on five epigenetic measures of aging was examined before and approximately 3 years after HIV infection in the same individuals (n=102). Significant epigenetic age acceleration (median 1.9-4.8 years) and estimated telomere length shortening (all p≤ 0.001) were observed from pre-to post-HIV infection, and remained significant in three epigenetic measures after controlling for T cell changes. No acceleration was seen in age- and time interval-matched HIV-uninfected controls. Changes in genome-wide co-methylation clusters were also significantly associated with initial HIV infection (p≤ 2.0 × 10^-4). These longitudinal observations clearly demonstrate an early and substantial impact of HIV infection on the epigenetic aging process, and suggest a role for HIV itself in the earlier onset of clinical aging.

Genome-Wide Association Study of Campylobacter-Positive Diarrhea Identifies Genes Involved in Toxin Processing and Inflammatory Response

Diarrhea is responsible for the deaths of more than 500,000 children each year, many of whom reside in low-to-middle-income countries (LMICs). Additionally, children with multiple diarrheal infections early in life have increased growth stunting and malnutrition and decreased vaccine efficacy. Two bacteria that contribute to the burden of diarrhea are Campylobacter jejuni and Campylobacter coli, both are endemic in Bangladesh. However, not all children that are exposed to these pathogens, including Campylobacter, will experience diarrhea. We hypothesized that host genetics may influence susceptibility to Campylobacter infections and performed a genome-wide association study in 534 children from two independent birth cohorts in Dhaka, Bangladesh. Infants were monitored for diarrhea for the first 2 years of life and only defined as controls if all diarrheal samples in the first year were negative for Campylobacter jejuni/C. coli. Each cohort was analyzed separately under an additive model and adjusted for length-for-age z-scores at birth and 12 months, sex, water treatment, and ancestry. In a fixed effect inverse-variance weighted meta-analysis of these two cohorts, we identified a genome-wide significant region on chromosome 8 in intron 4 of the rho guanine nucleotide exchange factor 10 gene (ARHGEF10). Individuals with the G allele (rs13281104) had a 2-fold lower risk of having a Campylobacter-associated diarrheal episode than individuals with the A allele (OR 0.41, 95% CI 0.29 to 0.58, P = 3.6 × 10-7). This SNP is associated with decreased expression of the neighboring gene, CLN8, which may be involved in the transport of the cytolethal distending toxin produced by Campylobacter. IMPORTANCE Children in low-to-middle-income countries often suffer from multiple enteric infections in their first few years of life, many of which have the potential for long-lasting effects. These children are already likely to be malnourished and underweight, and chronic intestinal disturbances exacerbate these conditions. Despite public health interventions aimed at improving water, sanitation, and hygiene, enteric infections are still a leading cause of death for children under five. Previous work has included transmission dynamics, pathogen characteristics, and evaluation of interventions. Here, we examined the role of host genetic variation in susceptibility to diarrhea-associated Campylobacter infection. In our meta-analysis of two independent birth cohorts from Dhaka, Bangladesh, we found that children carrying a specific genetic variant (rs13281104, in an intron of ARHGEF10) were half as likely to have a diarrhea-associated Campylobacter infection in their first year of life. This protective effect may be achieved by decreasing gene expression and thereby impacting host-pathogen interactions and host immune response.

Phage display of environmental protein toxins and virulence factors reveals the prevalence, persistence, and genetics of antibody responses

Microbial exposures are crucial environmental factors that impact healthspan by sculpting the immune system and microbiota. Antibody profiling via Phage ImmunoPrecipitation Sequencing (PhIP-Seq) provides a high-throughput, cost-effective approach for detecting exposure and response to microbial protein products. We designed and constructed a library of 95,601 56-amino acid peptide tiles spanning 14,430 proteins with "toxin" or "virulence factor" keyword annotations. We used PhIP-Seq to profile the antibodies of ∼1,000 individuals against this "ToxScan" library. In addition to enumerating immunodominant antibody epitopes, we studied the age-dependent stability of the ToxScan profile and used a genome-wide association study to find that the MHC-II locus modulates bacterial epitope selection. We detected previously described anti-flagellin antibody responses in a Crohn's disease cohort and identified an association between anti-flagellin antibodies and juvenile dermatomyositis. PhIP-Seq with the ToxScan library is thus an effective tool for studying the environmental determinants of health and disease at cohort scale.

Increased Rate of Epigenetic Aging in Men Living With HIV Prior to Treatment

Background: Epigenetic aging is accelerated in tissues of persons living with HIV (PLWH) and may underlie the early onset of age-related illnesses. This study examines the rate-of-change in epigenetic age in PLWH following HIV infection but before HAART, using archived longitudinal samples from the Multicenter AIDS Cohort Study. Methods: DNA was isolated from cryopreserved peripheral blood mononuclear cells from 101 men living with HIV, with baseline visit <2.5 years after HIV seroconversion (Visit 1) and follow-up visit <1.5 years before the initiation of HAART (Visit 2), and 100 HIV-uninfected men matched on age and visits with comparable time intervals. DNA methylation (DNAm) age was estimated for five clocks (Pan-tissue, Extrinsic, Phenotypic, Grim, and Skin & Blood age), and a DNAm-based estimate of telomere length (DNAmTL). Multivariate linear regression models were used to examine baseline factors associated with rate-of-aging, defined as (DNAm age visit 2-DNAm age visit 1)/(age visit 2-age visit 1). Results: Epigenetic age increased approximately twice as fast in PLWH as uninfected controls (Pan-tissue, Extrinsic, and Phenotypic clocks). Shortening of DNAmTL was nearly 3-fold faster in PLWH than controls. Faster rate-of-aging was associated with HIV status (Pan-Tissue, Extrinsic, Phenotypic, and DNAmTL), white race (Extrinsic, DNAmTL), higher cumulative HIV viral load (Grim), and lower baseline DNAm age (Phenotypic, Skin & Blood). Conclusion: Epigenetic rates-of-aging were significantly faster for untreated PLWH. Our findings expand on the important impact of HIV infection on biologic aging, both in elevating epigenetic age and increasing the rate-of-aging in the years following infection.

Trans-ancestral fine-mapping of MHC reveals key amino acids associated with spontaneous clearance of hepatitis C in HLA-DQβ1

Spontaneous clearance of acute hepatitis C virus (HCV) infection is associated with single nucleotide polymorphisms (SNPs) on the MHC class II. We fine-mapped the MHC region in European (n = 1,600; 594 HCV clearance/1,006 HCV persistence) and African (n = 1,869; 340 HCV clearance/1,529 HCV persistence) ancestry individuals and evaluated HCV peptide binding affinity of classical alleles. In both populations, HLA-DQβ1Leu26 (p value_Meta = 1.24 × 10^-14) located in pocket 4 was negatively associated with HCV spontaneous clearance and HLA-DQβ1Pro55 (p value_Meta = 8.23 × 10^-11) located in the peptide binding region was positively associated, independently of HLA-DQβ1Leu26. These two amino acids are not in linkage disequilibrium (r² < 0.1) and explain the SNPs and classical allele associations represented by rs2647011, rs9274711, HLA-DQB1^∗03:01, and HLA-DRB1^∗01:01. Additionally, HCV persistence classical alleles tagged by HLA-DQβ1Leu26 had fewer HCV binding epitopes and lower predicted binding affinities compared to clearance alleles (geometric mean of combined IC₅₀ nM of persistence versus clearance; 2,321 nM versus 761.7 nM, p value = 1.35 × 10^-38). In summary, MHC class II fine-mapping revealed key amino acids in HLA-DQβ1 explaining allelic and SNP associations with HCV outcomes. This mechanistic advance in understanding of natural recovery and immunogenetics of HCV might set the stage for much needed enhancement and design of vaccine to promote spontaneous clearance of HCV infection.

Post-sequelae symptoms and comorbidities after COVID-19

The frequency, severity, and forms of symptoms months after coronavirus 2019 (COVID‐19) are poorly understood, especially in community settings. To better understand and characterize symptoms months after community‐based COVID‐19, a retrospective cohort analysis was conducted. Three hundred and twenty‐eight consecutive persons with a positive test for SARS‐CoV‐2 in the Johns Hopkins Health System, Maryland, March−May 2020, were selected for the study. Symptom occurrence and severity were measured through questionnaires. Of 328 persons evaluated, a median of 242 days (109−478 days) from the initial positive SARS‐CoV‐2 test, 33.2% reported not being fully recovered and 4.9% reported symptoms that constrained daily activities. Compared to those who reported being fully recovered, those with post‐acute sequelae were more likely to report a prior history of heart attack (p < 0.01). Among those reporting long‐term symptoms, men and women were equally represented (men = 34.8%, women = 34.6%), but only women reported symptoms that constrained daily activities, and 56% of them were caregivers. The types of new or persistent symptoms varied, and for many, included a deviation from prior COVID‐19 health, such as being less able to exercise, walk, concentrate, or breathe. A limitation is that self‐report of symptoms might be biased and/or caused by factors other than COVID‐19. Overall, even in a community setting, symptoms may persist months after COVID‐19 reducing daily activities including caring for dependents.

•

Even months after coronavirus 2019 a substantial proportion of persons continue to have symptoms that might restrict their daily activities.

•

Further research is needed to prevent this complication especially as the pandemic spreads over the world.

Analysis of antibody binding specificities in twin and SNP-genotyped cohorts reveals that antiviral antibody epitope selection is a heritable trait

Human immune responses to viral infections are highly variable, but the genetic factors that contribute to this variability are not well characterized. We used VirScan, a high-throughput epitope scanning technology, to analyze pan-viral antibody reactivity profiles of twins and SNP-genotyped individuals. Using these data, we determined the heritability and genomic loci associated with antibody epitope selection, response breadth, and control of Epstein-Barr virus (EBV) viral load. 107 EBV peptide reactivities were heritable and at least two Epstein-Barr nuclear antigen 2 (EBNA-2) reactivities were associated with variants in the MHC class II locus. We identified an EBV serosignature that predicted viral load in peripheral blood mononuclear cells and was associated with variants in the MHC class I locus. Our study illustrates the utility of epitope profiling to investigate the genetics of pathogen immunity, reports heritable features of the antibody response to viruses, and identifies specific HLA loci important for EBV epitope selection.

Association analyses of rare variants identify two genes associated with refractive error

PURPOSE

Genetic variants identified through population-based genome-wide studies are generally of high frequency, exerting their action in the central part of the refractive error spectrum. However, the power to identify associations with variants of lower minor allele frequency is greatly reduced, requiring considerable sample sizes. Here we aim to assess the impact of rare variants on genetic variation of refractive errors in a very large general population cohort.

METHODS

Genetic association analyses of non-cyclopaedic autorefraction calculated as mean spherical equivalent (SPHE) used whole-exome sequence genotypic information from 50,893 unrelated participants in the UK Biobank of European ancestry. Gene-based analyses tested for association with SPHE using an optimised SNP-set kernel association test (SKAT-O) restricted to rare variants (minor allele frequency < 1%) within protein-coding regions of the genome. All models were adjusted for age, sex and common lead variants within the same locus reported by previous genome-wide association studies. Potentially causal markers driving association at significant loci were elucidated using sensitivity analyses by sequentially dropping the most associated variants from gene-based analyses.

RESULTS

We found strong statistical evidence for association of SPHE with the SIX6 (p-value = 2.15 x 10-10, or Bonferroni-Corrected p = 4.41x10-06) and the CRX gene (p-value = 6.65 x 10-08, or Bonferroni-Corrected p = 0.001). The SIX6 gene codes for a transcription factor believed to be critical to the eye, retina and optic disc development and morphology, while CRX regulates photoreceptor specification and expression of over 700 genes in the retina. These novel associations suggest an important role of genes involved in eye morphogenesis in refractive error.

CONCLUSION

The results of our study support previous research highlighting the importance of rare variants to the genetic risk of refractive error. We explain some of the origins of the genetic signals seen in GWAS but also report for the first time a completely novel association with the CRX gene.

Spatiotemporal Phylodynamics of Hepatitis C Among People Who Inject Drugs in India

BACKGROUND AND AIMS

Implementing effective interventions for HCV requires a detailed understanding of local transmission dynamics and geospatial spread. Little is known about HCV phylodynamics, particularly among high-burden populations, such as people who inject drugs (PWID).

APPROACH AND RESULTS

We used 483 HCV sequences and detailed individual-level data from PWID across four Indian cities. Bayesian phylogeographic analyses were used to evaluate transmission hotspots and geospatial diffusion of the virus. Phylogenetic cluster analysis was performed to infer epidemiologic links and factors associated with clustering. A total of 492 HIV sequences were used to draw comparisons within the same population and, in the case of coinfections, evaluate molecular evidence for shared transmission pathways. Overall, 139/483 (28.8%) of HCV sequences clustered with a median cluster size of 3 individuals. Genetically linked participants with HCV were significantly younger and more likely to be infected with HCV subtype 3b as well as to live and inject close to one another. Phylogenetic evidence suggests likely ongoing HCV infection/reinfection with limited support for shared HIV/HCV transmission pathways. Phylogeographic analyses trace historic HCV spread back to Northeastern India and show diffusion patterns consistent with drug trafficking routes.

CONCLUSIONS

This study characterizes HCV phylodynamics among PWID in a low and middle-income country setting. Heterogeneity and recent genetic linkage of HCV across geographically disparate Indian states suggest that targeted interventions could help prevent reimportation of virus through drug trafficking routes.

Role of direct and indirect social and spatial ties in the diffusion of HIV and HCV among people who inject drugs: a cross-sectional community-based network analysis in New Delhi, India

Background:

People who inject drugs (PWID) account for some of the most explosive human immunodeficiency virus (HIV) and hepatitis C virus (HCV) epidemics globally. While individual drivers of infection are well understood, less is known about network factors, with minimal data beyond direct ties.

Methods:

2512 PWID in New Delhi, India were recruited in 2017–19 using a sociometric network design. Sampling was initiated with 10 indexes who recruited named injection partners (people who they injected with in the prior month). Each recruit then recruited their named injection partners following the same process with cross-network linkages established by biometric data. Participants responded to a survey, including information on injection venues, and provided a blood sample. Factors associated with HIV/HCV infection were identified using logistic regression.

Results:

The median age was 26; 99% were male. Baseline HIV prevalence was 37.0% and 46.8% were actively infected with HCV (HCV RNA positive). The odds of prevalent HIV and active HCV infection decreased with each additional degree of separation from an infected alter (HIV AOR: 0.87; HCV AOR: 0.90) and increased among those who injected at a specific venue (HIV AOR: 1.50; HCV AOR: 1.69) independent of individual-level factors (p<0.001). In addition, sociometric factors, for example, network distance to an infected alter, were statistically significant predictors even when considering immediate egocentric ties.

Conclusions:

These data demonstrate an extremely high burden of HIV and HCV infection and a highly interconnected injection and spatial network structure. Incorporating network and spatial data into the design/implementation of interventions may help interrupt transmission while improving efficiency.

Funding:

National Institute on Drug Abuse and the Johns Hopkins University Center for AIDS Research.

Understanding the social and spatial relationships that connect people is a key element to stop the spread of infectious diseases. These networks are particularly relevant to combat epidemics among populations that are hard to reach with public health interventions. Network-based approaches, for example, can help to stop HIV or hepatitis C from spreading amongst populations that use injectable drugs. Yet how social and geographic connections such as acquaintances, injection partners, or preferred drug use places impact the risk of infection is still poorly mapped out.

To address this question, Clipman et al. focused on people who inject drugs in New Delhi, India, a population heavily impacted by HIV and hepatitis C. Over 2500 people were recruited, each participant inviting their injection partners to also take part. The volunteers answered survey questions, including where they used drugs, and provided a blood sample to be tested.

The results showed that, even after adjusting for individual risk factors, where people used drugs and with whom affected their risk of becoming infected with HIV and hepatitis C. In terms of social ties, the likelihood of HIV and hepatitis C infection decreased by about 13% for each person separating a given individual from an infected person. However, geographical networks also had a major impact. Injecting at a popular location respectively increased the odds of HIV and hepatitis C infection by 50% and 69%. In fact, even if the participant was not using drugs at these specific places, having an injection partner who did was enough to increase the risk for disease: for each person separating an individual from the location, the likelihood of being infected with HIV and hepatitis C decreased by respectively 14% and 10%.

The results by Clipman et al. highlight how the relationships between physical spaces and social networks contribute to the spread of dangerous diseases amongst people who inject drugs. Ultimately, this knowledge may help to shape better public health interventions that would take into account the importance of geographical locations.

APOL1 variant alleles associate with reduced risk for opportunistic infections in HIV infection

Apolipoprotein L1 (APOL1), an innate immune factor against African trypanosoma brucei, inhibits HIV-1 in vitro. The impact of APOL1 G1-G2 variants on HIV-1-associated opportunistic infections (OIs) is unknown. Here, we report findings from a metaanalysis of four HIV/AIDS prospective cohorts (ALIVE, LSOCA, MACS, and WIHS) including 2066 African American participants. Using a global test combining all four cohorts, carriage of two APOL1 variant alleles is associated with a 50% reduction in odds of OI (combined OR 0.50, 95% CI 0.33-0.76). Subgroup analysis of OI etiological categories (viral, parasitic, fungal and Mycobacterial) suggests the possibility of specific protection from fungal infections (OR 0.54. 95% CI 0.32-0.93; P_{Bonferroni corrected} = 0.08). We observe an association of APOL1 variant alleles with host protection against OI in HIV-positive individuals. The study suggests a broader role of APOL1 variant alleles in innate immunity in vivo.

An et al. determine the presence of variants of the innate immune factor APOL1 in four cohorts of HIV-positive patients. The study suggests that APOL1 might confer carriers of two variant alleles protection from HIV related opportunistic infections, especially fungal infections.

Enterovirus D68 molecular and cellular biology and pathogenesis

In recent years, enterovirus D68 (EV-D68) has advanced from a rarely detected respiratory virus to a widespread pathogen responsible for increasing rates of severe respiratory illness and acute flaccid myelitis (AFM) in children worldwide. In this review, we discuss the accumulating data on the molecular features of EV-D68 and place these into the context of enterovirus biology in general. We highlight similarities and differences with other enteroviruses and genetic divergence from own historical prototype strains of EV-D68. These include changes in capsid antigens, host cell receptor usage, and viral RNA metabolism collectively leading to increased virulence. Furthermore, we discuss the impact of EV-D68 infection on the biology of its host cells, and how these changes are hypothesized to contribute to motor neuron toxicity in AFM. We highlight areas in need of further research, including the identification of its primary receptor and an understanding of the pathogenic cascade leading to motor neuron injury in AFM. Finally, we discuss the epidemiology of the EV-D68 and potential therapeutic approaches.

Mendelian Randomization Analysis of n-6 Polyunsaturated Fatty Acid Levels and Pancreatic Cancer Risk

BACKGROUND

Whether circulating polyunsaturated fatty acid (PUFA) levels are associated with pancreatic cancer risk is uncertain. Mendelian randomization (MR) represents a study design using genetic instruments to better characterize the relationship between exposure and outcome.

METHODS

We utilized data from genome-wide association studies within the Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case-Control Consortium, involving approximately 9,269 cases and 12,530 controls of European descent, to evaluate associations between pancreatic cancer risk and genetically predicted plasma n-6 PUFA levels. Conventional MR analyses were performed using individual-level and summary-level data.

RESULTS

Using genetic instruments, we did not find evidence of associations between genetically predicted plasma n-6 PUFA levels and pancreatic cancer risk [estimates per one SD increase in each PUFA-specific weighted genetic score using summary statistics: linoleic acid odds ratio (OR) = 1.00, 95% confidence interval (CI) = 0.98-1.02; arachidonic acid OR = 1.00, 95% CI = 0.99-1.01; and dihomo-gamma-linolenic acid OR = 0.95, 95% CI = 0.87-1.02]. The OR estimates remained virtually unchanged after adjustment for covariates, using individual-level data or summary statistics, or stratification by age and sex.

CONCLUSIONS

Our results suggest that variations of genetically determined plasma n-6 PUFA levels are not associated with pancreatic cancer risk.

IMPACT

These results suggest that modifying n-6 PUFA levels through food sources or supplementation may not influence risk of pancreatic cancer.

HLA tapasin independence: broader peptide repertoire and HIV control

Human leukocyte antigen (HLA) class I allotypes vary in their ability to present peptides in the absence of tapasin, an essential component of the peptide loading complex. We quantified tapasin dependence of all allotypes that are common in European and African Americans (n = 97), which revealed a broad continuum of values. Ex vivo examination of cytotoxic T cell responses to the entire HIV-1 proteome from infected subjects indicates that tapasin-dependent allotypes present a more limited set of distinct peptides than do tapasin-independent allotypes, data supported by computational predictions. This suggests that variation in tapasin dependence may impact the strength of the immune responses by altering peptide repertoire size. In support of this model, we observed that individuals carrying HLA class I genotypes characterized by greater tapasin independence progress more slowly to AIDS and maintain lower viral loads, presumably due to increased breadth of peptide presentation. Thus, tapasin dependence level, like HLA zygosity, may serve as a means to restrict or expand breadth of the HLA-I peptide repertoire across humans, ultimately influencing immune responses to pathogens and vaccines.

Genomics in the era of COVID-19: ethical implications for clinical practice and public health

Genomic studies of patients with COVID-19, or exposed to it, are underway to delineate host factors associated with variability in susceptibility, infectivity, and disease severity. Here, we highlight the ethical implications-both potential benefits and harms-of genomics for clinical practice and public health in the era of COVID-19.

Multi-ancestry fine mapping of interferon lambda and the outcome of acute hepatitis C virus infection

Clearance of acute infection with hepatitis C virus (HCV) is associated with the chr19q13.13 region containing the rs368234815 (TT/ΔG) polymorphism. We fine-mapped this region to detect possible causal variants that may contribute to HCV clearance. First, we performed sequencing of IFNL1-IFNL4 region in 64 individuals sampled according to rs368234815 genotype: TT/clearance (N = 16) and ΔG/persistent (N = 15) (genotype-outcome concordant) or TT/persistent (N = 19) and ΔG/clearance (N = 14) (discordant). 25 SNPs had a difference in counts of alternative allele >5 between clearance and persistence individuals. Then, we evaluated those markers in an association analysis of HCV clearance conditioning on rs368234815 in two groups of European (692 clearance/1 025 persistence) and African ancestry (320 clearance/1 515 persistence) individuals. 10/25 variants were associated (P < 0.05) in the conditioned analysis leaded by rs4803221 (P value = 4.9 × 10^-04) and rs8099917 (P value = 5.5 × 10^-04). In the European ancestry group, individuals with the haplotype rs368234815ΔG/rs4803221C were 1.7× more likely to clear than those with the rs368234815ΔG/rs4803221G haplotype (P value = 3.6 × 10^-05). For another nearby SNP, the haplotype of rs368234815ΔG/rs8099917T was associated with HCV clearance compared to rs368234815ΔG/rs8099917G (OR: 1.6, P value = 1.8 × 10^-04). We identified four possible causal variants: rs368234815, rs12982533, rs10612351 and rs4803221. Our results suggest a main signal of association represented by rs368234815, with contributions from rs4803221, and/or nearby SNPs including rs8099917.

A Multiancestry Sex-Stratified Genome-Wide Association Study of Spontaneous Clearance of Hepatitis C Virus

BACKGROUND

Spontaneous clearance of acute hepatitis C virus (HCV) infection is more common in women than in men, independent of known risk factors.

METHODS

To identify sex-specific genetic loci, we studied 4423 HCV-infected individuals (2903 male, 1520 female) of European, African, and Hispanic ancestry. We performed autosomal, and X chromosome sex-stratified and combined association analyses in each ancestry group.

RESULTS

A male-specific region near the adenosine diphosphate-ribosylation factor-like 5B (ARL5B) gene was identified. Individuals with the C allele of rs76398191 were about 30% more likely to have chronic HCV infection than individuals with the T allele (OR, 0.69; P = 1.98 × 10-07), and this was not seen in females. The ARL5B gene encodes an interferon-stimulated gene that inhibits immune response to double-stranded RNA viruses. We also identified suggestive associations near septin 6 and ribosomal protein L39 genes on the X chromosome. In box sexes, allele G of rs12852885 was associated with a 40% increase in HCV clearance compared with the A allele (OR, 1.4; P = 2.46 × 10-06). Septin 6 facilitates HCV replication via interaction with the HCV NS5b protein, and ribosomal protein L39 acts as an HCV core interactor.

CONCLUSIONS

These novel gene associations support differential mechanisms of HCV clearance between the sexes and provide biological targets for treatment or vaccine development.

The Role of Host Genetic Factors in Coronavirus Susceptibility: Review of Animal and Systematic Review of Human Literature

The SARS-CoV-2 pandemic raises many scientific and clinical questions. These include how host genetic factors affect disease susceptibility and pathogenesis. New work is emerging related to SARS-CoV-2; previous work has been conducted on other coronaviruses that affect different species. We reviewed the literature on host genetic factors related to coronaviruses, systematically focusing on human studies. We identified 1,832 articles of potential relevance. Seventy-five involved human host genetic factors, 36 of which involved analysis of specific genes or loci; aside from one meta-analysis, all were candidate-driven studies, typically investigating small numbers of research subjects and loci. Three additional case reports were described. Multiple significant loci were identified, including 16 related to susceptibility (seven of which identified protective alleles) and 16 related to outcomes (three of which identified protective alleles). The types of cases and controls used varied considerably; four studies used traditional replication/validation cohorts. Among other studies, 30 involved both human and non-human host genetic factors related to coronavirus, 178 involved study of non-human (animal) host genetic factors related to coronavirus, and 984 involved study of non-genetic host factors related to coronavirus, including involving immunopathogenesis. Previous human studies have been limited by issues that may be less impactful now, including low numbers of eligible participants and limited availability of advanced genomic methods; however, these may raise additional considerations. We outline key genes and loci from animal and human host genetic studies that may bear investigation in the study of COVID-19. We also discuss how previous studies may direct current lines of inquiry.

Perspectives on the Future of Epidemiology: A Framework for Training

Over the past century, the field of epidemiology has evolved and adapted to changing public health needs. Challenges include newly emerging public health concerns across broad and diverse content areas, new methods, and vast data sources. We recognize the need to engage and educate the next generation of epidemiologists and prepare them to tackle these issues of the 21st century. In this commentary, we suggest a skeleton framework upon which departments of epidemiology should build their curriculum. We propose domains that include applied epidemiology, biological and social determinants of health, communication, creativity and ability to collaborate and lead, statistical methods, and study design. We believe all students should gain skills across these domains to tackle the challenges posed to us. The aim is to train smart thinkers, not technicians, to embrace challenges and move the expanding field of epidemiology forward.

The Ethics of Precision Rationing: Human Genetics and the Need for Debate on Stratifying Access to Medication

Rising prices for new, transformative therapies are challenging health systems around the world, leading many payers and providers to begin rationing access to treatments, even in the countries that have been most resistant to doing so. This is the case for direct-acting antivirals (DAAs) for the treatment of hepatitis C virus (HCV). However, little attention has been paid to the increasing role that human genetics might play in rationing decisions. Researchers have already proposed that genetic markers associated with spontaneous HCV clearance could be used to restrict DAA access for some patients, although treatment would be medically beneficial for those patients. Would such forms of rationing present a form of genetic discrimination? And what of the public health implications of these approaches? Here we present an ethical analysis of such proposals for "precision rationing" and raise 4 key areas of concern. We argue that ethical issues arising in this area are not substantively different from the pressing ethical issues regarding rationing and discrimination more broadly, but provide important impetus for motivating broad public debate to find ethically sound ways of managing genomics and new expensive medications.

The Role of Host Genetic Factors in Coronavirus Susceptibility: Review of Animal and Systematic Review of Human Literature

BACKGROUND

The recent SARS-CoV-2 pandemic raises many scientific and clinical questions. One set of questions involves host genetic factors that may affect disease susceptibility and pathogenesis. New work is emerging related to SARS-CoV-2; previous work has been conducted on other coronaviruses that affect different species.

OBJECTIVES

We aimed to review the literature on host genetic factors related to coronaviruses, with a systematic focus on human studies.

METHODS

We conducted a PubMed-based search and analysis for articles relevant to host genetic factors in coronavirus. We categorized articles, summarized themes related to animal studies, and extracted data from human studies for analyses.

RESULTS

We identified 1,187 articles of potential relevance. Forty-five studies were related to human host genetic factors related to coronavirus, of which 35 involved analysis of specific genes or loci; aside from one meta-analysis on respiratory infections, all were candidate-driven studies, typically investigating small number of research subjects and loci. Multiple significant loci were identified, including 16 related to susceptibility to coronavirus (of which 7 identified protective alleles), and 16 related to outcomes or clinical variables (of which 3 identified protective alleles). The types of cases and controls used varied considerably; four studies used traditional replication/validation cohorts. Of the other studies, 28 involved both human and non-human host genetic factors related to coronavirus, 174 involved study of non-human (animal) host genetic factors related to coronavirus, 584 involved study of non-genetic host factors related to coronavirus, including involving immunopathogenesis, 16 involved study of other pathogens (not coronavirus), 321 involved other studies of coronavirus, and 18 studies were assigned to the other categories and removed.

KEY FINDINGS

We have outlined key genes and loci from animal and human host genetic studies that may bear investigation in the nascent host genetic factor studies of COVID-19. Previous human studies to date have been limited by issues that may be less impactful on current endeavors, including relatively low numbers of eligible participants and limited availability of advanced genomic methods.

Prevalence and Phylogenetic Characterization of Hepatitis C Virus Among Indian Men Who Have Sex With Men: Limited Evidence for Sexual Transmission

BACKGROUND

Data from high-income countries suggest increasing hepatitis C virus (HCV) prevalence/incidence among human immunodeficiency virus (HIV)-infected men who have sex with men (MSM), but limited data derive from low-and-middle-income countries.

METHODS

We recruited 4994 MSM from 5 states across India using respondent-driven sampling. Logistic regression incorporating respondent-driven sampling weights and machine learning feature selection were used to identify correlates of prevalent HCV, and Bayesian phylogenetic analysis was used to examine genetic clustering.

RESULTS

The median age was 25 years, the HIV prevalence was 7.2%, and 49.3% of participants reported recent unprotected anal intercourse. The HCV prevalence was 1.3% (95% confidence interval, 1.0%-1.6%; site range, 0.2%-3.4%) and was 3.1% in HIV-positive versus 1.1% among HIV-negative men. HCV infection was significantly associated with injection drug use (odds ratio, 177.1; 95% confidence interval, 72.7-431.5) and HIV infection (4.34; 1.88-10.05). Machine learning did not uncover any additional epidemiologic signal. Phylogenetic analysis revealed 3 clusters suggestive of linked transmission; each contained ≥1 individual reporting injection drug use.

CONCLUSIONS

We observed a low HCV prevalence in this large sample of MSM despite a high prevalence of known risk factors, reflecting either the need for a threshold of HCV for sexual transmission and/or variability in sexual practices across settings.

Case-Control Study of Cryptosporidium Transmission in Bangladeshi Households

Background

Cryptosporidium is a leading contributor to diarrheal morbidity and mortality in under-5 children worldwide. As there is no vaccine and no effective drug therapy in young children for this infection, preventing infection is critical. We undertook a pilot case-control study to define the extent of person-to-person transmission of cryptosporidiosis within an urban and a rural community in Bangladesh.

Methods

We enrolled 48 case families with a Cryptosporidium-infected child aged 6–18 months. Controls were age- and sex-matched Cryptosporidium-negative children in 12 households. Children and household members were followed for 8 weeks with weekly illness survey and stool testing with quantitative polymerase chain reaction for Cryptosporidium.

Results

In the 24 urban case families, the secondary attack rate was 35.8% (19/53) vs 0% (0/11) in controls (P = .018, χ² test). In contrast, in the 24 rural case families, the secondary attack rate was 7.8% (5/64) vs 0% (0/21) in controls (P = .19, χ² test). Genotyping by gp60 demonstrated infection with the same subspecies in 5 families, and evidence of transmission in 2. Serologic response to Cryptosporidium infection was associated with younger age, longer duration of infection, and Cryptosporidium hominis gp60_IbA9G3R2 infection.

Conclusions

In the urban site, the high rate of secondary infection and infection with the same subspecies within families suggests that person-to-person transmission is a major source of Cryptosporidium infection for young children living in this region. Molecular genotyping can be applied to determine transmission of Cryptosporidium in endemic regions. Further work is needed to understand the differences in parasite transmissibility and immunity to different genotypes.

Genomics and Infectious Diseases: Expert Perspectives on Public Health Considerations regarding Actionability and Privacy

There is growing evidence that human genetics plays a significant role in shaping human responses to infectious diseases. For instance, individuals' genetic susceptibility or resistance to infectious disease is likely to affect disease transmission. Yet little attention has been paid to the ethical, legal, and social implications of research in genomics and infectious disease, despite the unique ethical issues that arise in this arena. This article presents results from a pilot study exploring ethics in research on human genetics and response to HIV and other infectious diseases and is focused on perspectives from expert stakeholders. Whereas chairs of institutional review boards, biobank directors, and researchers in genomics and infectious disease expressed similar views about research privacy in the context of a public health emergency, they expressed different perspectives about the role that public health considerations ought to play in the return of individual results to research participants. These perspectives highlight the need to emphasize the importance of broad dialogue for helping various parties navigate the ethically complex current and future challenges of genomics and infectious disease research.

Genome-Wide Association Study of Cryptosporidiosis in Infants Implicates PRKCA

Diarrhea is a major cause of both morbidity and mortality worldwide, especially among young children. Cryptosporidiosis is a leading cause of diarrhea in children, particularly in South Asia and sub-Saharan Africa, where it is responsible for over 200,000 deaths per year. Beyond the initial clinical presentation of diarrhea, it is associated with long-term sequelae such as malnutrition and neurocognitive developmental deficits. Risk factors include poverty and overcrowding, and yet not all children with these risk factors and exposure are infected, nor do all infected children develop symptomatic disease. One potential risk factor to explain these differences is their human genome. To identify genetic variants associated with symptomatic cryptosporidiosis, we conducted a genome-wide association study (GWAS) examining 6.5 million single nucleotide polymorphisms (SNPs) in 873 children from three independent cohorts in Dhaka, Bangladesh, namely, the Dhaka Birth Cohort (DBC), the Performance of Rotavirus and Oral Polio Vaccines in Developing Countries (PROVIDE) study, and the Cryptosporidiosis Birth Cohort (CBC). Associations were estimated separately for each cohort under an additive model, adjusting for length-for-age Z-score at 12 months of age, the first two principal components to account for population substructure, and genotyping batch. The strongest meta-analytic association was with rs58296998 (P = 3.73 × 10-8), an intronic SNP and expression quantitative trait locus (eQTL) of protein kinase C alpha (PRKCA). Each additional risk allele conferred 2.4 times the odds of Cryptosporidium-associated diarrhea in the first year of life. This genetic association suggests a role for protein kinase C alpha in pediatric cryptosporidiosis and warrants further investigation.IMPORTANCE Globally, diarrhea remains one of the major causes of pediatric morbidity and mortality. The initial symptoms of diarrhea can often lead to long-term consequences for the health of young children, such as malnutrition and neurocognitive developmental deficits. Despite many children having similar exposures to infectious causes of diarrhea, not all develop symptomatic disease, indicating a possible role for human genetic variation. Here, we conducted a genetic study of susceptibility to symptomatic disease associated with Cryptosporidium infection (a leading cause of diarrhea) in three independent cohorts of infants from Dhaka, Bangladesh. We identified a genetic variant within protein kinase C alpha (PRKCA) associated with higher risk of cryptosporidiosis in the first year of life. These results indicate a role for human genetics in susceptibility to cryptosporidiosis and warrant further research to elucidate the mechanism.

The Evolving Field of Genetic Epidemiology: From Familial Aggregation to Genomic Sequencing

The field of genetic epidemiology is relatively young and brings together genetics, epidemiology, and biostatistics to identify and implement the best study designs and statistical analyses for identifying genes controlling risk for complex and heterogeneous diseases (i.e., those where genes and environmental risk factors both contribute to etiology). The field has moved quickly over the past 40 years partly because the technology of genotyping and sequencing has forced it to adapt while adhering to the fundamental principles of genetics. In the last two decades, the available tools for genetic epidemiology have expanded from a genetic focus (considering 1 gene at a time) to a genomic focus (considering the entire genome), and now they must further expand to integrate information from other “-omics” (e.g., epigenomics, transcriptomics as measured by RNA expression) at both the individual and the population levels. Additionally, we can now also evaluate gene and environment interactions across populations to better understand exposure and the heterogeneity in disease risk. The future challenges facing genetic epidemiology are considerable both in scale and techniques, but the importance of the field will not diminish because by design it ties scientific goals with public health applications.