Integrated analyses of gene expression and genetic association studies in a founder population.

Genome-wide association studies (GWASs) have become a standard tool for dissecting genetic contributions to disease risk. However, these studies typically require extraordinarily large sample sizes to be adequately powered. Strategies that incorporate functional information alongside genetic associations have proved successful in increasing GWAS power. Following this paradigm, we present the results of 20 different genetic association studies for quantitative traits related to complex diseases, conducted in the Hutterites of South Dakota. To boost the power of these association studies, we collected RNA-sequencing data from lymphoblastoid cell lines for 431 Hutterite individuals. We then used Sherlock, a tool that integrates GWAS and expression quantitative trait locus (eQTL) data, to identify weak GWAS signals that are also supported by eQTL data. Using this approach, we found novel associations with quantitative phenotypes related to cardiovascular disease, including carotid intima-media thickness, left atrial volume index, monocyte count and serum YKL-40 levels.

The combination of a genome-wide association study of lymphocyte count and analysis of gene expression data reveals novel asthma candidate genes.

Recent genome-wide association studies (GWAS) have identified a number of novel genetic associations with complex human diseases. In spite of these successes, results from GWAS generally explain only a small proportion of disease heritability, an observation termed the 'missing heritability problem'. Several sources for the missing heritability have been proposed, including the contribution of many common variants with small individual effect sizes, which cannot be reliably found using the standard GWAS approach. The goal of our study was to explore a complimentary approach, which combines GWAS results with functional data in order to identify novel genetic associations with small effect sizes. To do so, we conducted a GWAS for lymphocyte count, a physiologic quantitative trait associated with asthma, in 462 Hutterites. In parallel, we performed a genome-wide gene expression study in lymphoblastoid cell lines from 96 Hutterites. We found significant support for genetic associations using the GWAS data when we considered variants near the 193 genes whose expression levels across individuals were most correlated with lymphocyte counts. Interestingly, these variants are also enriched with signatures of an association with asthma susceptibility, an observation we were able to replicate. The associated loci include genes previously implicated in asthma susceptibility as well as novel candidate genes enriched for functions related to T cell receptor signaling and adenosine triphosphate synthesis. Our results, therefore, establish a new set of asthma susceptibility candidate genes. More generally, our observations support the notion that many loci of small effects influence variation in lymphocyte count and asthma susceptibility.

The maternal HLA-G 1597ΔC null mutation is associated with increased risk of pre-eclampsia and reduced HLA-G expression during pregnancy in African-American women.

The non-classical major histocompatibility complex molecule, human leukocyte antigen (HLA)-G, is thought to contribute to maternal immune tolerance and successful placentation during pregnancy. Genetic polymorphisms in HLA-G are known to influence expression levels as well as the relative expression of individual protein isoforms. As diminished or aberrant HLA-G expression patterns may contribute to the development of certain pregnancy complications, we sought to investigate the association between functional HLA-G polymorphisms and the risk of pre-eclampsia (PE) in African-American women. The association between maternal and fetal genotype at six HLA-G polymorphisms and risk of PE was assessed in 372 pregnancies (314 normotensive; 58 pre-eclamptic). We observed an elevated risk of PE (P = 0.00027) in pregnancies where the mother carried the 1597ΔC allele, a null allele that abolishes expression of full-length HLA-G isoforms. Furthermore, the frequency of the maternal 1597ΔC allele was highest in the subset of pre-eclamptic pregnancies that were delivered preterm, suggesting an association between the null allele and the severity of PE. We then replicated the association between higher maternal 1597ΔC allele frequency and increased severity of PE (P = 0.038) in an independent sample of 533 African-American women. Finally, to investigate the mechanistic basis of this association, we measured circulating soluble HLA-G (sHLA-G) concentrations in maternal serum collected during pregnancy in 51 healthy, normotensive African-American control women and found significantly lower levels in women carrying the 1597ΔC allele (P = 0.012). These results demonstrate that maternal HLA-G genotype is significantly associated with risk of PE in African-American women and is predictive of circulating sHLA-G levels during pregnancy.

Single cell profiling at the maternal-fetal interface reveals a deficiency of PD-L1+ non-immune cells in human spontaneous preterm labor.

The mechanisms that underlie the timing of labor in humans are largely unknown. In most pregnancies, labor is initiated at term (≥ 37 weeks gestation), but in a signifiicant number of women spontaneous labor occurs preterm and is associated with increased perinatal mortality and morbidity. The objective of this study was to characterize the cells at the maternal-fetal interface (MFI) in term and preterm pregnancies in both the laboring and non-laboring state in Black women, who have among the highest preterm birth rates in the U.S. Using mass cytometry to obtain high-dimensional single-cell resolution, we identified 31 cell populations at the MFI, including 25 immune cell types and six non-immune cell types. Among the immune cells, maternal PD1+ CD8 T cell subsets were less abundant in term laboring compared to term non-laboring women. Among the non-immune cells, PD-L1+ maternal (stromal) and fetal (extravillous trophoblast) cells were less abundant in preterm laboring compared to term laboring women. Consistent with these observations, the expression of CD274, the gene encoding PD-L1, was significantly depressed and less responsive to fetal signaling molecules in cultured mesenchymal stromal cells from the decidua of preterm compared to term women. Overall, these results suggest that the PD1/PD-L1 pathway at the MFI may perturb the delicate balance between immune tolerance and rejection and contribute to the onset of spontaneous preterm labor.

Sex-specific differences in peripheral blood leukocyte transcriptional response to LPS are enriched for HLA region and X chromosome genes.

Sex-specific differences in prevalence are well documented for many common, complex diseases, especially for immune-mediated diseases, yet the precise mechanisms through which factors associated with biological sex exert their effects throughout life are not well understood. We interrogated sex-specific transcriptional responses of peripheral blood leukocytes (PBLs) to innate immune stimulation by lipopolysaccharide (LPS) in 46 male and 66 female members of the Hutterite community, who practice a communal lifestyle. We identified 1217 autosomal and 54 X-linked genes with sex-specific responses to LPS, as well as 71 autosomal and one X-linked sex-specific expression quantitative trait loci (eQTLs). Despite a similar proportion of the 15 HLA genes responding to LPS compared to all expressed autosomal genes, there was a significant over-representation of genes with sex by treatment interactions among HLA genes. We also observed an enrichment of sex-specific differentially expressed genes in response to LPS for X-linked genes compared to autosomal genes, suggesting that HLA and X-linked genes may disproportionately contribute to sex disparities in risk for immune-mediated diseases.

Establishment of human induced trophoblast stem-like cells from term villous cytotrophoblasts.

Human trophoblast stem cells (hTSC) can be isolated from first trimester placenta but not from term placenta. Here we demonstrate that villous cytotrophoblasts (vCTB) from term placenta can be reprogrammed into induced trophoblastic stem-like cells (iTSC) by introducing sets of transcription factors. The iTSCs express TSC markers such as GATA3, TEAD4 and ELF5, and are multipotent, validated by their differentiation into both extravillous trophoblasts (EVT) and syncytiotrophoblasts (STB) in vitro and in vivo. The iTSC can be passaged indefinitely in vitro without slowing of growth. The transcriptome profile of these cells closely resembles the profile of hTSC isolated from first trimester placentae but different from the term placental vCTB from which they originated. The ability to reprogram cells from term placenta into iTSC will allow study of early gestation events which impact placental function later in gestation, including preeclampsia and spontaneous preterm birth.

Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma.

There is a life-long relationship between rhinovirus (RV) infection and the development and clinical manifestations of asthma. In this study we demonstrate that cultured primary bronchial epithelial cells from adults with asthma (n = 9) show different transcriptional and chromatin responses to RV infection compared to those without asthma (n = 9). Both the number and magnitude of transcriptional and chromatin responses to RV were muted in cells from asthma cases compared to controls. Pathway analysis of the transcriptionally responsive genes revealed enrichments of apoptotic pathways in controls but inflammatory pathways in asthma cases. Using promoter capture Hi-C we tethered regions of RV-responsive chromatin to RV-responsive genes and showed enrichment of these regions and genes at asthma GWAS loci. Taken together, our studies indicate a delayed or prolonged inflammatory state in cells from asthma cases and highlight genes that may contribute to genetic risk for asthma.

HLA-G polymorphisms and soluble HLA-G protein levels in women with recurrent pregnancy loss from Basrah province in Iraq.

HLA-G is a nonclassical, class I major histocompatibility complex (MHC) gene that exhibits immunomodulatory properties and likely plays a role in the maintenance of successful pregnancy. In this study, we investigated the role of HLA-G polymorphisms on risk for recurrent pregnancy loss (RPL) and on circulating levels of soluble (s)HLA-G in Iraqi women. DNA and plasma were obtained from blood samples collected at 9-12 weeks gestation from 50 women with RPL and 50 healthy pregnant women in Basrah province, Iraq. As measured by ELISA, median sHLA-G levels were significantly lower in the RPL cases compared to healthy controls (21.4 vs. 38.8 U/ml, respectively; P=0.025), and decreased with increasing maternal age (P=0.0051). However, HLA-G allele and haplotype frequencies did not differ significantly between cases and controls (P values ≥0.12 for all tests). In contrast, homozygosity for the C allele (CC) at a tri-allelic promoter polymorphism, -725C/G/T, was associated with lower concentrations of sHLA-G compared to the CG or CT genotypes (median levels 21.1 vs. 40.1 vs. 42.6 U/ml, respectively; P=0.0089). These results demonstrate that HLA-G genotype influences circulating sHLA-G levels during pregnancy but is not significantly associated with risk of RPL.

Increased protein-coding mutations in the mitochondrial genome of African American women with preeclampsia.

Preeclampsia occurs more frequently in women of African ancestry. The cause of this hypertensive complication is unclear, but placental oxidative stress may play a role. Because mitochondria are the major sites of oxidative phosphorylation, we hypothesized that placentas of preeclamptic pregnancies harbor mitochondrial DNA (mtDNA) mutations. Next-generation sequencing of placental mtDNA in African American preeclamptics (N = 30) and controls (N = 38) from Chicago revealed significant excesses in preeclamptics of nonsynonymous substitutions in protein-coding genes and mitochondrially encoded nicotinamide adenine dinucleotide dehydrogenase 5 gene and an increase in the substitution rate (P = .0001). Moreover, 88% of preeclamptics and 53% of controls carried at least one nonsynonymous substitution (P = .005; odds ratio [OR] = 6.36, 95% confidence interval [CI]: 1.5-39.1). These results were not replicated in a sample of African American preeclamptics (N = 162) and controls (N = 171) from Detroit. Differences in study design and heterogeneity may account for this lack of replication. Nonsynonymous substitutions in mtDNA may be risk factors for preeclampsia in some African American women, but additional studies are required to establish this relationship.

Genome-wide association study of plasma lipoprotein(a) levels identifies multiple genes on chromosome 6q.

Plasma lipoprotein(a) (Lp[a]) level is an independent risk factor of cardiovascular disease that is under strong genetic control. We conducted a genome-wide association study of plasma Lp(a) in 386 members of a founder population that adheres to a communal lifestyle, proscribes cigarette smoking, and prepares and eats meals communally. We identified associations with 77 single nucleotide polymorphisms (SNPs) spanning 12.5 Mb on chromosome 6q26-q27 that met criteria for genome-wide significance (P

Fine mapping and positional candidate studies identify HLA-G as an asthma susceptibility gene on chromosome 6p21.

Asthma affects nearly 14 million people worldwide and has been steadily increasing in frequency for the past 50 years. Although environmental factors clearly influence the onset, progression, and severity of this disease, family and twin studies indicate that genetic variation also influences susceptibility. Linkage of asthma and related phenotypes to chromosome 6p21 has been reported in seven genome screens, making it the most replicated region of the genome. However, because many genes with individually small effects are likely to contribute to risk, identification of asthma susceptibility loci has been challenging. In this study, we present evidence from four independent samples in support of HLA-G as a novel asthma and bronchial hyperresponsiveness susceptibility gene in the human leukocyte antigen region on chromosome 6p21, and we speculate that this gene might contribute to risk for other inflammatory diseases that show linkage to this region.

Variation in the HLA-G promoter region influences miscarriage rates.

The HLA-G gene is primarily expressed in placental cells that invade the maternal decidua during pregnancy. This gene encodes multiple isoforms that fulfill a variety of functions at the maternal-fetal interface throughout gestation. Recently, a null allele for the most abundant HLA-G isoform was associated with recurrent miscarriage in two independent studies, suggesting that reduced levels of the HLA-G1 protein may compromise successful pregnancy. We initiated the present study to determine whether other polymorphisms that could affect expression levels of HLA-G were associated with fetal loss in women participating in a 15-year prospective study of pregnancy outcome. We genotyped these subjects for 18 single-nucleotide polymorphisms in the 1,300 bp upstream of exon 1, 13 of which were identified as part of this study, as well as for an insertion/deletion (in/del) polymorphism in the 3' untranslated region. The 18 SNPs defined eight unique haplotypes. One polymorphism, -725C/G, was associated with fetal loss, with an increased risk for miscarriage in couples in which both partners carried the -725G allele, compared with couples not carrying this allele (odds ratio 2.76, 95% confidence interval 1.08-7.09; P=.035). Further, the G at nucleotide -725 creates a CpG dinucleotide, and we demonstrate that this CpG site is methylated on -725G alleles. Overall, this study identified extraordinary levels of variation in the 5'-upstream regulatory region of HLA-G and provides evidence for an association between a promoter-region SNP and fetal loss rates, further attesting to the novel features and critical role of this gene in pregnancy.