Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Host diversity and behavior determine patterns of interspecies transmission and geographic diffusion of avian influenza A subtypes among North American wild reservoir species.

Wild birds can carry avian influenza viruses (AIV), including those with pandemic or panzootic potential, long distances. Even though AIV has a broad host range, few studies account for host diversity when estimating AIV spread. We analyzed AIV genomic sequences from North American wild birds, including 303 newly sequenced isolates, to estimate interspecies and geographic viral transition patterns among multiple co-circulating subtypes. Our results show high transition rates within Anseriformes and Charadriiformes, but limited transitions between these orders. Patterns of transition between species were positively associated with breeding habitat range overlap, and negatively associated with host genetic distance. Distance between regions (negative correlation) and summer temperature at origin (positive correlation) were strong predictors of transition between locations. Taken together, this study demonstrates that host diversity and ecology can determine evolutionary processes that underlie AIV natural history and spread. Understanding these processes can provide important insights for effective control of AIV.

Proteomic Architecture of Human Coronary and Aortic Atherosclerosis.

BACKGOUND: The inability to detect premature atherosclerosis significantly hinders implementation of personalized therapy to prevent coronary heart disease. A comprehensive understanding of arterial protein networks and how they change in early atherosclerosis could identify new biomarkers for disease detection and improved therapeutic targets. METHODS: Here we describe the human arterial proteome and proteomic features strongly associated with early atherosclerosis based on mass spectrometry analysis of coronary artery and aortic specimens from 100 autopsied young adults (200 arterial specimens). Convex analysis of mixtures, differential dependent network modeling, and bioinformatic analyses defined the composition, network rewiring, and likely regulatory features of the protein networks associated with early atherosclerosis and how they vary across 2 anatomic distributions. RESULTS: The data document significant differences in mitochondrial protein abundance between coronary and aortic samples (coronary>>aortic), and between atherosclerotic and normal tissues (atherosclerotic<

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Genome-wide association study in Chinese identifies novel loci for blood pressure and hypertension.

Hypertension is a common disorder and the leading risk factor for cardiovascular disease and premature deaths worldwide. Genome-wide association studies (GWASs) in the European population have identified multiple chromosomal regions associated with blood pressure, and the identified loci altogether explain only a small fraction of the variance for blood pressure. The differences in environmental exposures and genetic background between Chinese and European populations might suggest potential different pathways of blood pressure regulation. To identify novel genetic variants affecting blood pressure variation, we conducted a meta-analysis of GWASs of blood pressure and hypertension in 11 816 subjects followed by replication studies including 69 146 additional individuals. We identified genome-wide significant (P < 5.0 × 10(-8)) associations with blood pressure, which included variants at three new loci (CACNA1D, CYP21A2, and MED13L) and a newly discovered variant near SLC4A7. We also replicated 14 previously reported loci, 8 (CASZ1, MOV10, FGF5, CYP17A1, SOX6, ATP2B1, ALDH2, and JAG1) at genome-wide significance, and 6 (FIGN, ULK4, GUCY1A3, HFE, TBX3-TBX5, and TBX3) at a suggestive level of P = 1.81 × 10(-3) to 5.16 × 10(-8). These findings provide new mechanistic insights into the regulation of blood pressure and potential targets for treatments.

Mutations in folate transporter genes and risk for human myelomeningocele.

The molecular mechanisms linking folate deficiency and neural tube defect (NTD) risk in offspring remain unclear. Folate transporters (SLC19A1, SLC46A1, SLC25A32, and FOLH1) and folate receptors (FOLR1, FOLR2, and FOLR3) are suggested to play essential roles in transporting folate from maternal intestinal lumen to the developing embryo. Loss of function variants in these genes may affect folate availability and contribute to NTD risk. This study examines whether variants within the folate transporter and receptor genes are associated with an increased risk for myelomeningocele (MM). Exons and their flanking intron sequences of 348 MM subjects were sequenced using the Sanger sequencing method and/or next generation sequencing to identify variants. Frequencies of alleles of single nucleotide polymorphisms (SNPs) in MM subjects were compared to those from ethnically matched reference populations to evaluate alleles' associated risk for MM. We identified eight novel variants in SLC19A1 and twelve novel variants in FOLR1, FOLR2, and FOLR3. Pathogenic variants include c.1265delG in SLC19A1 resulting in an early stop codon, four large insertion deletion variants in FOLR3, and a stop_gain variant in FOLR3. No new variants were identified in SLC46A1, SLC25A32, or FOLH1. In SLC19A1, c.80A>G (rs1051266) was not associated with our MM cohort; we did observe a variant allele G frequency of 61.7%, higher than previously reported in other NTD populations. In conclusion, we discovered novel loss of function variants in genes involved in folate transport in MM subjects. Our results support the growing evidence of associations between genes involved in folate transport and susceptibility to NTDs.

Fecal Indole as a Biomarker of Susceptibility to Cryptosporidium Infection.

Cryptosporidium causes significant diarrhea worldwide, especially among children and immunocompromised individuals, and no effective drug treatment is currently available for those who need it most. In this report, previous volunteer infectivity studies have been extended to examine the association between fecal indole and indole-producing (IP) gut microbiota on the outcome of a Cryptosporidium infection. Fecal indole concentrations (FICs) of 50 subjects and 19 taxa of common gut microbiota, including six IP taxa (11 subjects) were determined in stool samples collected before and after a challenge with Cryptosporidium oocysts. At the baseline, the mean FIC (± the standard deviation) was 1.66 ± 0.80 mM in those who became infected after a challenge versus 3.20 ± 1.23 mM in those who remained uninfected (P = 0.0001). Only 11.1% of the subjects with a FIC of >2.5 mM became infected after a challenge versus 65.2% of the subjects with a FIC of <2.5 mM. In contrast, the FICs of infected subjects at the baseline or during diarrhea were not correlated with infection intensity or disease severity. The relative abundances (percent) of Escherichia coli, Bacillus spp., and Clostridium spp. were greater ≥2.5-fold in volunteers with a baseline FIC of >2.5 mM, while those of Bacteroides pyogenes, B. fragilis, and Akkermansia muciniphila were greater in those with a baseline FIC of <2.5 mM. These data indicate that some IP bacteria, or perhaps indole alone, can influence the ability of Cryptosporidium to establish an infection. Thus, preexisting indole levels in the gut join the oocyst dose and immune status as important factors that determine the outcome of Cryptosporidium exposure.

Human epithelial Na+ channel missense variants identified in the GenSalt study alter channel activity.

Mutations in genes encoding subunits of the epithelial Na+ channel (ENaC) can cause early onset familial hypertension, demonstrating the importance of this channel in modulating blood pressure. It remains unclear whether other genetic variants resulting in subtler alterations of channel function result in hypertension or altered sensitivity of blood pressure to dietary salt. This study sought to identify functional human ENaC variants to examine how these variants alter channel activity and to explore whether these variants are associated with altered sensitivity of blood pressure to dietary salt. Six-hundred participants of the Genetic Epidemiology Network of Salt Sensitivity (GenSalt) study with salt-sensitive or salt-resistant blood pressure underwent sequencing of the genes encoding ENaC subunits. Functional effects of identified variants were examined in a Xenopus oocyte expression system. Variants that increased channel activity included three in the gene encoding the α-subunit (αS115N, αR476W, and αV481M), one in the ß-subunit (ßS635N), and one in the γ-subunit (γL438Q). One α-subunit variant (αA334T) and one γ-subunit variant (ßD31N) decreased channel activity. Several α-subunit extracellular domain variants altered channel inhibition by extracellular Na+ (Na+ self-inhibition). One variant (αA334T) decreased and one (αV481M) increased cell surface expression. Association between these variants and salt sensitivity did not reach statistical significance. This study identifies novel functional human ENaC variants and demonstrates that some variants alter channel cell surface expression and/or Na+ self-inhibition.

Genetic association of the glycine cleavage system genes and myelomeningocele.

BACKGROUND: Neural tube defects (NTDs) are one of the most common congenital birth defects, with myelomeningocele (MM) being the most severe form compatible with life. Recent studies show a link between mitochondrial folate one carbon metabolism and NTDs by means of the glycine cleavage system (GCS). We hypothesize that single nucleotide polymorphisms and novel variants in the coding regions of the GCS genes increase the risk for MM. METHODS: DNA was obtained from 96 subjects with MM born before the United States mandated folic acid fortification of grains in 1998. Primers were designed for polymerase chain reaction amplification and sequencing of all exons in the AMT gene, one of four genes in the GCS, followed by identification of single nucleotide polymorphisms and novel variants. An additional 252 MM subjects underwent whole exome sequencing to examine all four GCS genes (aminomethyltransferase, glycine dehydrogenase, glycine cleavage system protein-H, and dihydrolipoamide dehydrogenase). RESULTS: We identified six novel, heterozygous variants in the AMT gene with three predicted to be deleterious to AMT function (p.Val7Leu, p.Pro251Arg, and p.Val380Met). Five extremely rare, known heterozygous variants were found in the AMT gene and one in the GLDC gene. No novel variants in the exons of the other two GCS genes (DLD and GCSH) were identified. CONCLUSION: We identified novel and rare, known variants in two of the four GCS genes that may contribute to the development of MM. Consistent with previous findings, the current study provides additional support that genetic variations in GCS genes contribute to the risk of NTDs. Birth Defects Research (Part A) 106:847-853, 2016. © 2016 Wiley Periodicals, Inc.

Aggregate blood pressure responses to serial dietary sodium and potassium intervention: defining responses using independent component analysis.

BACKGROUND: Hypertension is a complex trait that often co-occurs with other conditions such as obesity and is affected by genetic and environmental factors. Aggregate indices such as principal components among these variables and their responses to environmental interventions may represent novel information that is potentially useful for genetic studies. RESULTS: In this study of families participating in the Genetic Epidemiology Network of Salt Sensitivity (GenSalt) Study, blood pressure (BP) responses to dietary sodium interventions are explored. Independent component analysis (ICA) was applied to 20 variables indexing obesity and BP measured at baseline and during low sodium, high sodium and high sodium plus potassium dietary intervention periods. A "heat map" protocol that classifies subjects based on risk for hypertension is used to interpret the extracted components. ICA and heat map suggest four components best describe the data: (1) systolic hypertension, (2) general hypertension, (3) response to sodium intervention and (4) obesity. The largest heritabilities are for the systolic (64%) and general hypertension (56%) components. There is a pattern of higher heritability for the component response to intervention (40-42%) as compared to those for the traditional intervention responses computed as delta scores (24%-40%). CONCLUSIONS: In summary, the present study provides intermediate phenotypes that are heritable. Using these derived components may prove useful in gene discovery applications.

Genome-wide association analysis of incident coronary heart disease (CHD) in African Americans: a short report.

African Americans have the highest rate of mortality due to coronary heart disease (CHD). Although multiple loci have been identified influencing CHD risk in European-Americans using a genome-wide association (GWAS) approach, no GWAS of incident CHD has been reported for African Americans. We performed a GWAS for incident CHD events collected during 19 years of follow-up in 2,905 African Americans from the Atherosclerosis Risk in Communities (ARIC) study. We identified a genome-wide significant SNP (rs1859023, MAF = 31%) located at 7q21 near the PFTK1 gene (HR = 0.57, 95% CI 0.46 to 0.69, p = 1.86×10(-08)), which replicated in an independent sample of over 8,000 African American women from the Women's Health Initiative (WHI) (HR = 0.81, 95% CI 0.70 to 0.93, p = 0.005). PFTK1 encodes a serine/threonine-protein kinase, PFTAIRE-1, that acts as a cyclin-dependent kinase regulating cell cycle progression and cell proliferation. This is the first finding of incident CHD locus identified by GWAS in African Americans.

Exon sequencing of PAX3 and T (brachyury) in cases with spina bifida.

BACKGROUND: Based on studies in animals and humans, PAX3 and T (brachyury) are candidate genes for spina bifida. However, neither gene has been definitively identified as a risk factor for this condition. METHODS: Sanger sequencing was used to identify variants in all PAX3 and T exons and promoter regions in 114 spina bifida cases. For known variants, allele frequencies in cases were compared with those from public databases using unadjusted odds ratios. Novel variants were genotyped in parents and assessed for predicted functional impact. RESULTS: We identified common variants in PAX3 (n = 2) and T (n = 3) for which the allele frequencies in cases were significantly different from those reported in at least one public database. We also identified novel variants in both PAX3 (n = 11) and T (n = 1) in spina bifida cases. Several of the novel PAX3 variants are predicted to be highly conserved and/or impact gene function or expression. CONCLUSION: These studies provide some evidence that common variants of PAX3 and T are associated with spina bifida. Rare and novel variants in these genes were also identified in affected individuals. However, additional studies will be required to determine whether these variants influence the risk of spina bifida.

Linkage of high-density lipoprotein-cholesterol concentrations to a locus on chromosome 9p in Mexican Americans.

High-density lipoproteins (HDLs) are anti-atherogenic lipoproteins that have a major role in transporting cholesterol from peripheral tissues to the liver, where it is removed. Epidemiologic studies have shown that low levels of high-density lipoprotein-cholesterol (HDL-C) are associated with an increased incidence of coronary heart disease and an increased mortality rate, indicating a protective role of high concentrations of HDL-C against atherogenesis and the development of coronary heart disease. HDL-C level is influenced by several genetic and nongenetic factors. Nongenetic factors include smoking, which has been shown to decrease the HDL-C level. Exercise and alcohol have been shown to increase HDL-C levels. Decreased HDL-C is often associated with other coronary heart disease risk factors such as obesity, hyperinsulinemia and insulin resistance, hypertriglyceridemia and hypertension. Although several genes have been identified for rare forms of dyslipidemia, the genes accounting for major variation in HDL-C levels have yet to be identified. Using a multipoint variance components linkage approach, we found strong evidence of linkage (lod score=3.4; P=0.00004) of a quantitative trait locus (QTL) for HDL-C level to a genetic location between markers D9S925 and D9S741 on chromosome 9p in Mexican Americans. A replication study in an independent set of Mexican American families confirmed the existence of a QTL on chromosome 9p.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE.

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

The role of the kallikrein-kinin system genes in the salt sensitivity of blood pressure: the GenSalt Study.

The current study comprehensively examined the association between common genetic variants of the kallikrein-kinin system (KKS) and blood pressure salt sensitivity. A 7-day low-sodium followed by a 7-day high-sodium dietary intervention was conducted among 1,906 Han Chinese participants recruited from 2003 to 2005. Blood pressure was measured by using a random-zero sphygmomanometer through the study. A total of 205 single nucleotide polymorphisms (SNPs) covering 11 genes of the KKS were selected for the analyses. Genetic variants of the bradykinin receptor B2 gene (BDKRB2) and the endothelin converting enzyme 1 gene (ECE1) showed significant associations with the salt-sensitivity phenotypes even after adjustment for multiple testing. Compared with the major G allele, the BDKRB2 rs11847625 minor C allele was significantly associated with increased systolic blood pressure responses to low-sodium intervention (P = 0.0001). Furthermore, a haplotype containing allele C was associated with an increased systolic blood pressure response to high-sodium intervention (P = 0.0009). Seven highly correlated ECE1 SNPs were shown to increase the diastolic blood pressure response to low-sodium intervention (P values ranged from 0.0003 to 0.002), with 2 haplotypes containing these 7 SNPs also associated with this same phenotype (P values ranged from 0.0004 to 0.002). In summary, genetic variants of the genes involved in the regulation of KKS may contribute to the salt sensitivity of blood pressure.

Genome-wide linkage and positional association study of blood pressure response to dietary sodium intervention: the GenSalt Study.

The authors conducted a genome-wide linkage scan and positional association analysis to identify the genetic determinants of salt sensitivity of blood pressure (BP) in a large family-based, dietary-feeding study. The dietary intervention was conducted among 1,906 participants in rural China (2003-2005). A 7-day low-sodium intervention was followed by a 7-day high-sodium intervention. Salt sensitivity was defined as BP responses to low- and high-sodium interventions. Signals of the logarithm of the odds to the base 10 (LOD ≥ 3) were detected at 33-42 centimorgans of chromosome 2 (2p24.3-2p24.1), with a maximum LOD score of 3.33 for diastolic blood pressure responses to high-sodium intervention. LOD scores were 2.35-2.91 for mean arterial pressure (MAP) and 0.80-1.49 for systolic blood pressure responses in this region, respectively. Correcting for multiple tests, single nucleotide polymorphism (SNP) rs11674786 (2.7 kilobases upstream of the family with sequence similarity 84, member A, gene (FAM84A)) in the linkage region was significantly associated with diastolic blood pressure (P = 0.0007) and MAP responses (P = 0.0007), and SNP rs16983422 (2.8 kilobases upstream of the visinin-like 1 gene (VSNL1)) was marginally associated with diastolic blood pressure (P = 0.005) and MAP responses (P = 0.005). An additive interaction between SNPs rs11674786 and rs16983422 was observed, with P = 7.00 × 10(-5) and P = 7.23 × 10(-5) for diastolic blood pressure and MAP responses, respectively. The authors concluded that genetic region 2p24.3-2p24.1 might harbor functional variants for the salt sensitivity of BP.

Rare PPARA variants and extreme response to fenofibrate in the Genetics of Lipid-Lowering Drugs and Diet Network Study.

OBJECTIVE: Fenofibrate, a peroxisome proliferator-activated receptor-α (PPARα) agonist, reduces triglyceride (TG) concentrations by 25-60%. Given significant interindividual variations in the TG response, we investigated the association of PPARA rare variants with treatment response in the Genetics of Lipid-Lowering Drugs and Diet Network study. METHODS: We calculated the change in the TG concentration (ΔTG) among 861 GOLDN participants treated with fenofibrate (160 mg/day) for 3 weeks. From the distribution of ΔTG adjusted for age and sex, the 150 highest and 150 lowest fenofibrate responders were selected from the tails of the distribution for PPARA resequencing. The resequencing strategy was based on VariantSEQr technology for the amplification of exons and regulatory regions. RESULTS: We identified 73 variants with an average minor allele frequency of 4.8% (range: 0.2-16%). We tested the association of rare variants located in a coding or a regulatory region (minor allele frequency<1%, 13 variants) with treatment response group by an indicator variable (presence/absence of ≥1 rare variant) using general linear mixed models to allow for adjustment for family relationship. After adjusting for baseline, fasting TG concentration carrying at least one rare variant was associated with a low fenofibrate response (odds ratio=6.46; 95% confidence interval: 1.4-30.8). Carrier status was also associated with a relative change in the total cholesterol concentration (P=0.02), but not high-density lipoprotein or low-density lipoprotein concentration. CONCLUSION: Rare, potentially functional variants in PPARA may play a role in the TG response to fenofibrate, but future experimental studies will be necessary to replicate the findings and confirm functional effects.

Common variants in the periostin gene influence development of atherosclerosis in young persons.

OBJECTIVE: We investigated the influence of genetic variants (rare and common) in the gene encoding periostin (POSTN) on atherosclerosis as measured in arterial specimens from the Pathobiological Determinants of Atherosclerosis in Youth (PDAY) study. METHODS AND RESULTS: A comprehensive survey of common POSTN variants (87 single-nucleotide polymorphisms [SNPs]) in PDAY subjects (n = 2527) identified numerous SNPs associated with raised lesions in abdominal aorta and with fatty streaks in thoracic aorta. These SNPs belonged to a small number of correlation bins that spanned the entire locus. To examine effects of rare variants, we resequenced POSTN functional regions in PDAY cases with raised lesions (n = 291) and controls with no raised lesions (n = 294). However, we found no significant associations with case-control status for carriers of POSTN rare variants using the weighted-sum method for rare variant analysis. CONCLUSIONS: We identified common variants in POSTN that are associated with arterial lesions in young persons from the PDAY study. This finding strongly supports a role for periostin in atherogenesis, as suggested by recent proteomics analysis that found abundant expression of periostin in atherosclerotic lesions. Genetic variation may influence atherosclerosis via periostin's known involvement in multiple relevant pathways, including angiogenesis, vascular remodeling, and stimulation of migration and differentiation of vascular smooth muscle cells.

Novel and extendable genotyping system for human respiratory syncytial virus based on whole-genome sequence analysis.

BACKGROUND: Human respiratory syncytial virus (RSV) is one of the leading causes of respiratory infections, especially in infants and young children. Previous RSV sequencing studies have primarily focused on partial sequencing of G gene (200-300 nucleotides) for genotype characterization or diagnostics. However, the genotype assignment with G gene has not recapitulated the phylogenetic signal of other genes, and there is no consensus on RSV genotype definition. METHODS: We conducted maximum likelihood phylogenetic analysis with 10 RSV individual genes and whole-genome sequence (WGS) that are published in GenBank. RSV genotypes were determined by using phylogenetic analysis and pair-wise node distances. RESULTS: In this study, we first statistically examined the phylogenetic incongruence, rate variation for each RSV gene sequence and WGS. We then proposed a new RSV genotyping system based on a comparative analysis of WGS and the temporal distribution of strains. We also provide an RSV classification tool to perform RSV genotype assignment and a publicly accessible up-to-date instance of Nextstrain where the phylogenetic relationship of all genotypes can be explored. CONCLUSIONS: This revised RSV genotyping system will provide important information for disease surveillance, epidemiology, and vaccine development.

Mendelian randomization supports bidirectional causality between telomere length and clonal hematopoiesis of indeterminate potential.

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.