Search

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discover 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which are highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlie the coordinated expression of genes in the GNMs. Epigenetic analyses reveal that regulatory networks underlying self-renewal and pluripotency are more complex than previously realized. Genetic analyses identify thousands of regulatory variants that overlapped predicted transcription factor binding sites and are associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network are significantly enriched for regulatory variants with large effects, suggesting that they play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work bins tens of thousands of regulatory elements in hiPSCs into discrete regulatory networks, shows that pluripotency and self-renewal processes have a surprising level of regulatory complexity, and suggests that genetic factors may contribute to cell state transitions in human iPSC lines.

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Reactivation of fetal-specific genes and isoforms occurs during heart failure. However, the underlying molecular mechanisms and the extent to which the fetal program switch occurs remains unclear. Limitations hindering transcriptome-wide analyses of alternative splicing differences (i.e. isoform switching) in cardiovascular system (CVS) tissues between fetal, healthy adult and heart failure have included both cellular heterogeneity across bulk RNA-seq samples and limited availability of fetal tissue for research. To overcome these limitations, we have deconvoluted the cellular compositions of 996 RNA-seq samples representing heart failure, healthy adult (heart and arteria), and fetal-like (iPSC-derived cardiovascular progenitor cells) CVS tissues. Comparison of the expression profiles revealed that reactivation of fetal-specific RNA-binding proteins (RBPs), and the accompanied re-expression of 1,523 fetal-specific isoforms, contribute to the transcriptome differences between heart failure and healthy adult heart. Of note, isoforms for 20 different RBPs were among those that reverted in heart failure to the fetal-like expression pattern. We determined that, compared with adult-specific isoforms, fetal-specific isoforms encode proteins that tend to have more functions, are more likely to harbor RBP binding sites, have canonical sequences at their splice sites, and contain typical upstream polypyrimidine tracts. Our study suggests that compared with healthy adult, fetal cardiac tissue requires stricter transcriptional regulation, and that during heart failure reversion to this stricter transcriptional regulation occurs. Furthermore, we provide a resource of cardiac developmental stage-specific and heart failure-associated genes and isoforms, which are largely unexplored and can be exploited to investigate novel therapeutics for heart failure.

Combining single-cell cytometry datasets increases the analytical flexibility and the statistical power of data analyses. However, in many cases the full potential of co-analyses is not reached due to technical variance between data from different experimental batches. Here, we present cyCombine, a method to robustly integrate cytometry data from different batches, experiments, or even different experimental techniques, such as CITE-seq, flow cytometry, and mass cytometry. We demonstrate that cyCombine maintains the biological variance and the structure of the data, while minimizing the technical variance between datasets. cyCombine does not require technical replicates across datasets, and computation time scales linearly with the number of cells, allowing for integration of massive datasets. Robust, accurate, and scalable integration of cytometry data enables integration of multiple datasets for primary data analyses and the validation of results using public datasets.

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types.

Rationale: Interstitial lung abnormalities (ILAs) are associated with the highest genetic risk locus for idiopathic pulmonary fibrosis (IPF); however, the extent to which there are unique associations among individuals with ILAs or additional overlap with IPF is not known.Objectives: To perform a genome-wide association study (GWAS) of ILAs.Methods: ILAs and a subpleural-predominant subtype were assessed on chest computed tomography (CT) scans in the AGES (Age Gene/Environment Susceptibility), COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease [COPD]), Framingham Heart, ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-points), MESA (Multi-Ethnic Study of Atherosclerosis), and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study) studies. We performed a GWAS of ILAs in each cohort and combined the results using a meta-analysis. We assessed for overlapping associations in independent GWASs of IPF.Measurements and Main Results: Genome-wide genotyping data were available for 1,699 individuals with ILAs and 10,274 control subjects. The MUC5B (mucin 5B) promoter variant rs35705950 was significantly associated with both ILAs (P = 2.6 × 10-27) and subpleural ILAs (P = 1.6 × 10-29). We discovered novel genome-wide associations near IPO11 (rs6886640, P = 3.8 × 10-8) and FCF1P3 (rs73199442, P = 4.8 × 10-8) with ILAs, and near HTRE1 (rs7744971, P = 4.2 × 10-8) with subpleural-predominant ILAs. These novel associations were not associated with IPF. Among 12 previously reported IPF GWAS loci, five (DPP9, DSP, FAM13A, IVD, and MUC5B) were significantly associated (P

Despite the importance of understanding how variability across induced pluripotent stem cell (iPSC) lines due to non-genetic factors (clone and passage) influences their differentiation outcome, large-scale studies capable of addressing this question have not yet been conducted. Here, we differentiated 191 iPSC lines to generate iPSC-derived cardiovascular progenitor cells (iPSC-CVPCs). We observed cellular heterogeneity across the iPSC-CVPC samples due to varying fractions of two cell types: cardiomyocytes (CMs) and epicardium-derived cells (EPDCs). Comparing the transcriptomes of CM-fated and EPDC-fated iPSCs, we discovered that 91 signature genes and X chromosome dosage differences are associated with these two distinct cardiac developmental trajectories. In an independent set of 39 iPSCs differentiated into CMs, we confirmed that sex and transcriptional differences affect cardiac-fate outcome. Our study provides novel insights into how iPSC transcriptional and X chromosome gene dosage differences influence their response to differentiation stimuli and, hence, cardiac cell fate.

Rationale: Chronic obstructive pulmonary disease (COPD) has been associated with numerous genetic variants, yet the extent to which its genetic risk is mediated by variation in lung structure remains unknown.Objectives: To characterize associations between a genetic risk score (GRS) associated with COPD susceptibility and lung structure on computed tomography (CT).Methods: We analyzed data from MESA Lung (Multi-Ethnic Study of Atherosclerosis Lung Study), a U.S. general population-based cohort, and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study). A weighted GRS was calculated from 83 SNPs that were previously associated with lung function. Lung density, spatially matched airway dimensions, and airway counts were assessed on full-lung CT. Generalized linear models were adjusted for age, age squared, sex, height, principal components of genetic ancestry, smoking status, pack-years, CT model, milliamperes, and total lung volume.Measurements and Main Results: MESA Lung and SPIROMICS contributed 2,517 and 2,339 participants, respectively. Higher GRS was associated with lower lung function and increased COPD risk, as well as lower lung density, smaller airway lumens, and fewer small airways, without effect modification by smoking. Adjustment for CT lung structure, particularly small airway measures, attenuated associations between the GRS and FEV1/FVC by 100% and 60% in MESA and SPIROMICS, respectively. Lung structure (P  0.10), improved discrimination of moderate-to-severe COPD cases relative to clinical factors alone.Conclusions: A GRS associated with COPD susceptibility was associated with CT lung structure. Lung structure may be an important mediator of heritability and determinant of personalized COPD risk.

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N = 60,552), African (N = 8429), Asian (N = 9959), and Hispanic/Latino (N = 11,775) ethnicities. We identify over 50 additional loci at genome-wide significance in ancestry-specific or multiethnic meta-analyses. Using recent fine-mapping methods incorporating functional annotation, gene expression, and differences in linkage disequilibrium between ethnicities, we further shed light on potential causal variants and genes at known and newly identified loci. Several of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12. Our study highlights the utility of multiethnic and integrative genomics approaches to extend existing knowledge of the genetics of lung function and clinical relevance of implicated loci.

BackgroundIn Myanmar, civil unrest and the establishment of internally displaced person (IDP) settlements along the Myanmar-China border have impacted malaria transmission.MethodsMicrosatellite markers were used to examine source-sink dynamics for Plasmodium vivax between IDP settlements and surrounding villages in the border region. Genotypic structure and diversity were compared across the 3 years following the establishment of IDP settlements, to infer demographic history. We investigated whether human migration and landscape heterogeneity contributed to P. vivax transmission.ResultsP. vivax from IDP settlements and local communities consistently exhibited high genetic diversity within populations but low polyclonality within individuals. No apparent genetic structure was observed among populations and years. P. vivax genotypes in China were similar to those in Myanmar, and parasite introduction was unidirectional. Landscape factors, including distance, elevation, and land cover, do not appear to impede parasite gene flow.ConclusionsThe admixture of P. vivax genotypes suggested that parasite gene flow via human movement contributes to the spread of malaria both locally in Myanmar and across the international border. Our genetic findings highlight the presence of large P. vivax gene reservoirs that can sustain transmission. Thus, it is important to reinforce and improve existing control efforts along border areas.

Ethiopia is one of the few African countries where Plasmodium vivax is co-endemic with P. falciparum. Malaria transmission is seasonal and transmission intensity varies mainly by landscape and climate. Although the recent emergence of drug resistant parasites presents a major issue to malaria control in Ethiopia, little is known about the transmission pathways of parasite species and prevalence of resistant markers. This study used microsatellites to determine population diversity and gene flow patterns of P. falciparum (N = 226) and P. vivax (N = 205), as well as prevalence of drug resistant markers to infer the impact of gene flow and existing malaria treatment regimes. Plasmodium falciparum indicated a higher rate of polyclonal infections than P. vivax. Both species revealed moderate genetic diversity and similar population structure. Populations in the northern highlands were closely related to the eastern Rift Valley, but slightly distinct from the southern basin area. Gene flow via human migrations between the northern and eastern populations were frequent and mostly bidirectional. Landscape genetic analyses indicated that environmental heterogeneity and geographical distance did not constrain parasite gene flow. This may partly explain similar patterns of resistant marker prevalence. In P. falciparum, a high prevalence of mutant alleles was detected in codons related to chloroquine (pfcrt and pfmdr1) and sulfadoxine-pyrimethamine (pfdhps and pfdhfr) resistance. Over 60% of the samples showed pfmdr1 duplications. Nevertheless, no mutation was detected in pfK13 that relates to artemisinin resistance. In P. vivax, while sequences of pvcrt-o were highly conserved and less than 5% of the samples showed pvmdr duplications, over 50% of the samples had pvmdr1 976F mutation. It remains to be tested if this mutation relates to chloroquine resistance. Monitoring the extent of malaria spread and markers of drug resistance is imperative to inform policy for evidence-based antimalarial choice and interventions. To effectively reduce malaria burden in Ethiopia, control efforts should focus on seasonal migrant populations.

BackgroundRecent emergence of artemisinin-resistant P. falciparum has posed a serious hindrance to the elimination of malaria in the Greater Mekong Subregion. Parasite clearance time, a measure of change in peripheral parasitaemia in a sequence of samples taken after treatment, can be used to reflect the susceptibility of parasites or the efficiency of antimalarials. The association of genetic polymorphisms and artemisinin resistance has been documented. This study aims to examine clearance time of P. falciparum and P. vivax parasitemia as well as putative gene mutations associated with residual or recurred parasitemia in Myanmar.MethodsA total of 63 P. falciparum and 130 P. vivax samples collected from two internally-displaced populations and one surrounding village were examined for parasitemia changes. At least four samples were taken from each patient, at the first day of diagnosis up to 3 months following the initial treatment. The amount of parasite gene copy number was estimated using quantitative real-time PCR based on a species-specific region of the 18S rRNA gene. For samples that showed residual or recurred parasitemia after treatment, microsatellites were used to identify the 'post-treatment' parasite genotype and compared such with the 'pre-treatment' genotype. Mutations in genes pfcrt, pfmdr1, pfatp6, pfmrp1 and pfK13 that are potentially associated with ACT resistance were examined to identify if mutation is a factor for residual or persistent parasitemia.ResultsOver 30% of the P. falciprium infections showed delayed clearance of parasitemia after 2-3 days of treatment and 9.5% showed recurred parasitemia. Mutations in codon 876 of the pfmrp1 corroborated significance association with slow clearance time. However, no association was observed in the variation in pfmdr1 gene copy number as well as mutations of various codonsinpfatp6, pfcrt, and pfK13 with clearance time. For P. vivax, over 95% of the infections indicated cleared parasitemia at days 2-3 of treatment. Four samples were found to be re-infected with new parasite strains based on microsatellite genotypes after initial treatment.ConclusionThe appearance of P.falciparum infected samples showing delayed clearance or recurred parasitemia after treatment raises concerns on current treatment and ACT drug resistance.

Angelman syndrome (AS) is a neurological genetic disorder caused by loss of expression of the maternal copy of UBE3A in the brain. Due to brain-specific genetic imprinting at this locus, the paternal UBE3A is silenced by a long antisense transcript. Inhibition of the antisense transcript could lead to unsilencing of paternal UBE3A, thus providing a therapeutic approach for AS. However, widespread delivery of gene regulators to the brain remains challenging. Here, we report an engineered zinc finger-based artificial transcription factor (ATF) that, when injected i.p. or s.c., crossed the blood-brain barrier and increased Ube3a expression in the brain of an adult mouse model of AS. The factor displayed widespread distribution throughout the brain. Immunohistochemistry of both the hippocampus and cerebellum revealed an increase in Ube3a upon treatment. An ATF containing an alternative DNA-binding domain did not activate Ube3a. We believe this to be the first report of an injectable engineered zinc finger protein that can cause widespread activation of an endogenous gene in the brain. These observations have important implications for the study and treatment of AS and other neurological disorders.

The northern Andes, with their steep elevational and climate gradients, are home to an exceptional diversity of flora and fauna, particularly rich in avian species that have adapted to divergent ecological conditions. With this diversity comes the opportunity for parasites to exploit a wide breadth of avian hosts. However, little research has focused on examining the patterns of prevalence and lineage diversity of avian parasites in the Andes. Here, we screened a total of 428 birds from 19 species (representing nine families) and identified 133 infections of avian haemosporidia (31%), including lineages of Plasmodium, Haemoproteus, and Leucocytozoon. We document a higher prevalence of haemosporidia at higher elevations and lower temperatures, as well as an overall high diversity of lineages in the northern Andes, including the first sequences of haemosporidians reported in hummingbirds (31 sequences found in 11 species within the family Trochilidae). Double infections were distinguished using PHASE, which enables the separation of distinct parasite lineages. Results suggest that the ecological heterogeneity of the northern Andes that has given rise to a rich diversity of avian hosts may also be particularly conducive to parasite diversification and specialization.

Background: We conducted a genome-wide association study (GWAS) of subclinical interstitial lung disease (ILD), defined as high attenuation areas (HAA) on CT, in the population-based Multi-Ethnic Study of Atherosclerosis Study.Methods: We measured the percentage of high attenuation areas (HAA) in the lung fields on cardiac CT scan defined as voxels with CT attenuation values between -600 and -250 HU. Genetic analyses were performed in MESA combined across race/ethnic groups: non-Hispanic White (n = 2,434), African American (n = 2,470), Hispanic (n = 2,065) and Chinese (n = 702), as well as stratified by race/ethnicity.Results: Among 7,671 participants, regions at genome-wide significance were identified for basilar peel-core ratio of HAA in FLJ35282 downstream of ANRIL (rs7852363, P = 2.1x10-9) and within introns of SNAI3-AS1 (rs140142658, P = 9.6x10-9) and D21S2088E (rs3079677, P = 2.3x10-8). Within race/ethnic groups, 18 additional loci were identified at genome-wide significance, including genes related to development (FOXP4), cell adhesion (ALCAM) and glycosylation (GNPDA2, GYPC, GFPT1 and FUT10). Among these loci, SNP rs6844387 near GNPDA2 demonstrated nominal evidence of replication in analysis of n = 1,959 participants from the Framingham Heart Study (P = 0.029). FOXP4 region SNP rs2894439 demonstrated evidence of validation in analysis of n = 228 White ILD cases from the Columbia ILD Study compared to race/ethnicity-matched controls from MESA (one-sided P = 0.007). In lung tissue from 15 adults with idiopathic pulmonary fibrosis compared to 15 adults without lung disease. ANRIL (P = 0.001), ALCAM (P = 0.03) and FOXP4 (P = 0.046) were differentially expressed.Conclusions: Our results suggest novel roles for protein glycosylation and cell cycle disinhibition by long non-coding RNA in the pathogenesis of ILD. [ABSTRACT FROM AUTHOR]

In Africa, control programs that target primarily Plasmodium falciparum are inadequate for eliminating malaria. To learn more about prevalence and genetic variability of P. malariae in Africa, we examined blood samples from 663 asymptomatic and 245 symptomatic persons from western Kenya during June-August of 2014 and 2015. P. malariae accounted for 5.3% (35/663) of asymptomatic infections and 3.3% (8/245) of clinical cases. Among asymptomatic persons, 71% (32/45) of P. malariae infections detected by PCR were undetected by microscopy. The low sensitivity of microscopy probably results from the significantly lower parasitemia of P. malariae. Analyses of P. malariae circumsporozoite protein gene sequences revealed high genetic diversity among P. malariae in Africa, but no clear differentiation among geographic populations was observed. Our findings suggest that P. malariae should be included in the malaria elimination strategy in Africa and highlight the need for sensitive and field-applicable methods to identify P. malariae in malaria-endemic areas. [ABSTRACT FROM AUTHOR]