DNA methylation is an essential epigenetic mark whose role in gene regulation and its dependency on genomic sequence and environment are not fully understood. In this study we provide novel insights into the mechanistic relationships between genetic variation, DNA methylation and transcriptome sequencing data in three different cell-types of the GenCord human population cohort. We find that the association between DNA methylation and gene expression variation among individuals are likely due to different mechanisms from those establishing methylation-expression patterns during differentiation. Furthermore, cell-type differential DNA methylation may delineate a platform in which local inter-individual changes may respond to or act in gene regulation. We show that unlike genetic regulatory variation, DNA methylation alone does not significantly drive allele specific expression. Finally, inferred mechanistic relationships using genetic variation as well as correlations with TF abundance reveal both a passive and active role of DNA methylation to regulatory interactions influencing gene expression. DOI:http://dx.doi.org/10.7554/eLife.00523.001.

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Multiple myeloma is an incurable malignancy of plasma cells, and its pathogenesis is poorly understood. Here we report the massively parallel sequencing of 38 tumour genomes and their comparison to matched normal DNAs. Several new and unexpected oncogenic mechanisms were suggested by the pattern of somatic mutation across the data set. These include the mutation of genes involved in protein translation (seen in nearly half of the patients), genes involved in histone methylation, and genes involved in blood coagulation. In addition, a broader than anticipated role of NF-κB signalling was indicated by mutations in 11 members of the NF-κB pathway. Of potential immediate clinical relevance, activating mutations of the kinase BRAF were observed in 4% of patients, suggesting the evaluation of BRAF inhibitors in multiple myeloma clinical trials. These results indicate that cancer genome sequencing of large collections of samples will yield new insights into cancer not anticipated by existing knowledge.

Genome-wide data provide a powerful tool for inferring patterns of genetic variation and structure of human populations.

Primary biliary cirrhosis (PBC) is an uncommon autoimmune disease with a homogeneous clinical phenotype that reflects incomplete disease concordance in monozygotic (MZ) twins. We have taken advantage of a unique collection consisting of genomic DNA and mRNA from peripheral blood cells of female MZ twins (n = 3 sets) and sisters of similar age (n = 8 pairs) discordant for disease. We performed a genome-wide study to investigate differences in (i) DNA methylation (using a custom tiled four-plex array containing tiled 50-mers 19,084 randomly chosen methylation sites), (ii) copy number variation (CNV) (with a chip including markers derived from the 1000 Genomes Project, all three HapMap phases, and recently published studies), and/or (iii) gene expression (by whole-genome expression arrays). Based on the results obtained from these three approaches we utilized quantitative PCR to compare the expression of candidate genes. Importantly, our data support consistent differences in discordant twins and siblings for the (i) methylation profiles of 60 gene regions, (ii) CNV of 10 genes, and (iii) the expression of 2 interferon-dependent genes. Quantitative PCR analysis showed that 17 of these genes are differentially expressed in discordant sibling pairs. In conclusion, we report that MZ twins and sisters discordant for PBC manifest particular epigenetic differences and highlight the value of the epigenetic study of twins.

Non-additive interaction between genetic variants, or epistasis, is a possible explanation for the gap between heritability of complex traits and the variation explained by identified genetic loci. Interactions give rise to genotype dependent variance, and therefore the identification of variance quantitative trait loci can be an intermediate step to discover both epistasis and gene by environment effects (GxE). Using RNA-sequence data from lymphoblastoid cell lines (LCLs) from the TwinsUK cohort, we identify a candidate set of 508 variance associated SNPs. Exploiting the twin design we show that GxE plays a role in ∼70% of these associations. Further investigation of these loci reveals 57 epistatic interactions that replicated in a smaller dataset, explaining on average 4.3% of phenotypic variance. In 24 cases, more variance is explained by the interaction than their additive contributions. Using molecular phenotypes in this way may provide a route to uncovering genetic interactions underlying more complex traits.DOI: http://dx.doi.org/10.7554/eLife.01381.001.

Paramyotonia congenita (PC) is a rare autosomal dominant neuromuscular disorder characterized by juvenile onset and development of cold-induced myotonia after repeated activities. The disease is mostly caused by genetic mutations of the sodium channel, voltage-gated, type IV, alpha subunit (SCN4A) gene. This study intended to systematically identify the causative genetic variations of a Chinese Han PC family. Seven members of this PC family, including four patients and three healthy controls, were selected for whole exome sequencing (WES) using the Illumina HiSeq platform. Sequence variations were identified using the SoftGenetics program. The mutation R1448C of SCN4A was found to be the only causative mutation. This study applied WES technology to sequence multiple members of a large PC family and was the first to systematically confirm that the genetic change in SCN4A is the only causative variation in this PC family and the SCN4A mutation is sufficient to lead to PC.

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

Admixture mapping requires a genomewide panel of relatively evenly spaced markers that can distinguish the ancestral origins of chromosomal segments in admixed individuals. Through use of the results of the International HapMap Project and specific selection criteria, the current study has examined the ability of selected single-nucleotide polymorphisms (SNPs) to extract continental ancestry information in African American subjects and to explore parameters for admixture mapping. Genotyping of two linguistically diverse West African populations (Bini and Kanuri Nigerians, who are Niger-Congo [Bantu] and Nilo-Saharan speakers, respectively), European Americans, and African Americans validated a genomewide set of >4,000 SNP ancestry-informative markers with mean and median F(ST) values >0.59 and mean and median Fisher's information content >2.5. This set of SNPs extracted a larger amount of ancestry information in African Americans than previously reported SNP panels and provides nearly uniform coverage of the genome. Moreover, in the current study, simulations show that this more informative panel improves power for admixture mapping in African Americans when ethnicity risk ratios are modest. This is particularly important in the application of admixture mapping in complex genetic diseases for which only modest ethnicity risk ratios of relevant susceptibility genes are expected.

Using a genome-wide single nucleotide polymorphism (SNP) panel, we observed population structure in a diverse group of Europeans and European Americans. Under a variety of conditions and tests, there is a consistent and reproducible distinction between "northern" and "southern" European population groups: most individual participants with southern European ancestry (Italian, Spanish, Portuguese, and Greek) have >85% membership in the "southern" population; and most northern, western, eastern, and central Europeans have >90% in the "northern" population group. Ashkenazi Jewish as well as Sephardic Jewish origin also showed >85% membership in the "southern" population, consistent with a later Mediterranean origin of these ethnic groups. Based on this work, we have developed a core set of informative SNP markers that can control for this partition in European population structure in a variety of clinical and genetic studies.

Anti-nuclear antibodies to speckled 100 kDa (sp100) and glycoprotein 210 (gp210) are specific serologic markers of primary biliary cholangitis (PBC) of uncertain/controversial clinical or prognostic significance. To study the genetic determinants associated with sp100 and gp210 autoantibody subphenotypes, we performed a genome-wide association analysis of 930 PBC cases based on their autoantibody status, followed by a replication study in 1,252 PBC cases. We confirmed single-nucleotide polymorphisms rs492899 (P = 3.27 × 10-22 ; odds ratio [OR], 2.90; 95% confidence interval [CI], 2.34-3.66) and rs1794280 (P = 5.78 × 10-28 ; OR, 3.89; 95% CI, 3.05-4.96) in the human major histocompatibility complex (MHC) region associated with the sp100 autoantibody. However, no genetic variant was identified as being associated with the gp210 autoantibody. To further define specific classical human leukocyte antigen (HLA) alleles or amino acids associated with the sp100 autoantibody, we imputed 922 PBC cases (211 anti-sp100-positive versus 711 negative cases) using a Han Chinese MHC reference database. Conditional analysis identified that HLA-DRβ1-Asn77/Arg74, DRβ1-Ser37, and DPβ1-Lys65 were major determinants for sp100 production. For the classical HLA alleles, the strongest association was with DRB1*03:01 (P = 1.51 × 10-9 ; OR, 2.97; 95% CI, 2.06-4.29). Regression analysis with classical HLA alleles identified DRB1*03:01, DRB1*15:01, DRB1*01, and DPB1*03:01 alleles can explain most of the HLA association with sp100 autoantibody. Conclusion: This study indicated significant genetic predisposition to the sp100 autoantibody, but not the gp210 autoantibody, subphenotype in PBC patients. Additional studies will be necessary to determine if these findings have clinical significance to PBC pathogenesis and/or therapeutics.

Feed efficiency (FE) is an important trait affecting costs in swine industry. Investigation on FE-related genes in different tissues is valuable for molecular breeding. Hypothalamus is a convergent and integrated centre for multiple nutrient-related signals. The present study identified 363 differentially expressed (DE) genes and 14 DE lincRNAs in the hypothalamus of high- and low-FE Yorkshire pigs. Furthermore, 983 significantly correlated DE gene-lincRNA pairs were identified through weighted correlation network analysis (WGCNA) and Pearson correlation analysis. These DE genes were primarily enriched in the neuronal signal transduction process containing the upregulated genes of VIPR1, CCR1, CCR5, LEPR, INSR, ADRA1A, CCKAR, and ADORA3 and the downregulated genes of GRM1, GRM4, GRM5, and VIPR2, which were located in the cell membrane. These signal receptors were mainly connected to downstream Jak-STAT signaling that involved the increased genes (JAK2, STAT3, and POMC) and mTOR signaling pathway, including the decreased genes (CAMKK2, AMPK, and MTOR). STAT3 and AMPK genes also played a role in two major hypothalamic neurons of POMC and NPY/AGRP. A total of eight DE lincRNAs also participated in the potential network. In conclusion, neuronal signaling transduction-involved genes and lincRNAs were related to FE variation in pig hypothalamus.

Autotetraploid rice is a useful germplasm for polyploid rice breeding, however, low seed setting is the major barrier in commercial utilization of autotetraploid rice. Our research group has developed neo-tetraploid rice lines, which have the characteristics of high fertility and heterosis when crossed with autotetraploid rice. In the present study, re-sequencing and RNA-seq were employed to detect global DNA variations and differentially expressed genes (DEGs) during meiosis stage in three neo-tetraploid rice lines compared to their parents, respectively. Here, a total of 4109881 SNPs and 640592 InDels were detected in neo-tetraploid lines compared to the reference genome, and 1805 specific presence/absence variations (PAVs) were detected in three lines. Approximately 12% and 0.5% of the total SNPs and InDels identified in three lines were located in genic regions, respectively. A total of 28 genes, harboring at least one of the large-effect SNP and/or InDel which affect the integrity of the encoded protein, were identified in the three lines. Together, 324 specific mutation genes, including 52 meiosis-related genes and 8 epigenetics-related genes were detected in neo-tetraploid rice compared to their parents. Of these 324 genes, five meiosis-related and three epigenetics-related genes displayed differential expressions during meiosis stage. Notably, 498 specific transcripts, 48 differentially expressed transposons and 245 differentially expressed ncRNAs were also detected in neo-tetraploid rice. Our results suggested that genomic structural reprogramming, DNA variations and differential expressions of some important meiosis and epigenetics related genes might be associated with high fertility in neo-tetraploid rice.

Indoleamine 2,3-dioxygenase-1 (IDO1) is a key regulator of immune tolerance in ovarian cancer. This study investigated efficacy and safety of the IDO1 enzyme inhibitor epacadostat versus tamoxifen in patients with biochemical-only recurrence (CA-125 elevation) following complete remission after first-line chemotherapy for advanced epithelial ovarian, primary peritoneal, or fallopian tube cancer.

Common wild rice (Oryza rufipogon Griff.) is an important germplasm for rice breeding, which contains many resistance genes. Re-sequencing provides an unprecedented opportunity to explore the abundant useful genes at whole genome level. Here, we identified the nucleotide-binding site leucine-rich repeat (NBS-LRR) encoding genes by re-sequencing of two wild rice lines (i.e. Huaye 1 and Huaye 2) that were developed from common wild rice. We obtained 128 to 147 million reads with approximately 32.5-fold coverage depth, and uniquely covered more than 89.6% (> = 1 fold) of reference genomes. Two wild rice lines showed high SNP (single-nucleotide polymorphisms) variation rate in 12 chromosomes against the reference genomes of Nipponbare (japonica cultivar) and 93-11 (indica cultivar). InDels (insertion/deletion polymorphisms) count-length distribution exhibited normal distribution in the two lines, and most of the InDels were ranged from -5 to 5 bp. With reference to the Nipponbare genome sequence, we detected a total of 1,209,308 SNPs, 161,117 InDels and 4,192 SVs (structural variations) in Huaye 1, and 1,387,959 SNPs, 180,226 InDels and 5,305 SVs in Huaye 2. A total of 44.9% and 46.9% genes exhibited sequence variations in two wild rice lines compared to the Nipponbare and 93-11 reference genomes, respectively. Analysis of NBS-LRR mutant candidate genes showed that they were mainly distributed on chromosome 11, and NBS domain was more conserved than LRR domain in both wild rice lines. NBS genes depicted higher levels of genetic diversity in Huaye 1 than that found in Huaye 2. Furthermore, protein-protein interaction analysis showed that NBS genes mostly interacted with the cytochrome C protein (Os05g0420600, Os01g0885000 and BGIOSGA038922), while some NBS genes interacted with heat shock protein, DNA-binding activity, Phosphoinositide 3-kinase and a coiled coil region. We explored abundant NBS-LRR encoding genes in two common wild rice lines through genome wide re-sequencing, which proved to be a useful tool to exploit elite NBS-LRR genes in wild rice. The data here provide a foundation for future work aimed at dissecting the genetic basis of disease resistance in rice, and the two wild rice lines will be useful germplasm for the molecular improvement of cultivated rice.

Partial pollen and embryo sac sterilities are the two main reasons for low fertility in autotetraploid rice. Our previous study revealed that small RNAs changes may associate with pollen fertility in autotetraploid rice. However, knowledge on comparative analysis between the development of pollen and embryo sac by small RNAs in autotetraploid rice is still unknown. In the present study, WE-CLSM (whole-mount eosin B-staining confocal laser scanning microscopy) and high-throughput sequencing technology was employed to examine the cytological variations and to analyze small RNAs changes during pollen and embryo sac development in autotetraploid rice compared with its diploid counterpart.

Urothelial cancer-associated 1 (UCA1), an oncogenic long non-coding RNA, was aberrantly upregulated in colorectal cancer (CRC). This study aimed to further explore the clinical value of UCA1 in CRC.

Haplotype phasing of genetic variants is important for clinical interpretation of the genome, population genetic analysis and functional genomic analysis of allelic activity. Here we present phASER, an accurate approach for phasing variants that are overlapped by sequencing reads, including those from RNA sequencing (RNA-seq), which often span multiple exons due to splicing. Using diverse RNA-seq data we demonstrate that this provides more accurate phasing of rare variants compared with population-based phasing and allows phasing of variants in the same gene up to hundreds of kilobases away that cannot be obtained from DNA sequencing (DNA-seq) reads. We show that in the context of medical genetic studies this improves the resolution of compound heterozygotes. Additionally, phASER provides measures of haplotypic expression that increase power and accuracy in studies of allelic expression. In summary, phasing using RNA-seq and phASER is accurate and improves studies where rare variant haplotypes or allelic expression is needed.

Autotetraploid rice is a useful germplasm for polyploid rice breeding; however, low seed setting is a major hindrance for its utilization. Here, we reported the development of a new tetraploid rice, Huoduo1 (H1), which has the characteristic of high fertility, from crossing generations of autotetraploid rice. Cytological observations displayed the high fertility of the pollen (95.62%) in H1, a lower percentage of pollen mother cell (PMC) abnormalities, and stable chromosome configurations during the pollen development process compared with its parents. Using RNA-seq analysis, we detected 440 differentially expressed genes (DEGs) in H1 compared with its parents. Of these DEGs, 193 were annotated as pollen fertility-related genes, and 129 (~66.8%) exhibited significant up-regulation in H1 compared with the parents, including three environmentally sensitive genic male sterility genes (TMS9-1, TMS5, and CSA), one meiosis gene (RAD51D), and three tapetal-related genes (MIL2, OsAP25, and OsAP37), which were validated by qRT-PCR in this study. Two genes, TMS9-1 and TMS5, were knocked out using CRISPR/Cas9 technology, and their mutants displayed low fertility and the abnormal development of pollen. Our findings provide evidence for the regulatory mechanisms of fertility in tetraploid rice and indicated that the up-regulation of pollen fertility-related genes may contribute to the high fertility in new tetraploid rice.

The CD28 locus is associated with susceptibility to a variety of autoimmune and immune-mediated inflammatory diseases including primary sclerosing cholangitis (PSC). Previously, we linked the CD28 pathway in PSC disease pathology and found that vitamin D could maintain CD28 expression. Here, we assessed whether the PSC-associated CD28 risk variant A (rs7426056) affects CD28 expression and T cell function in healthy individuals (n = 14 AA, n = 14 AG, n = 14 GG). Homozygotes for the PSC disease risk allele (AA) showed significantly lower CD28 mRNA expression ex-vivo than either GG or AG (p < 0.001) in total peripheral blood mononuclear cells. However, the CD28 risk variant alone was not sufficient to explain CD28 protein loss on CD4+ T cells. All genotypes responded equally to vitamin D as indicated by induction of a regulatory phenotype and an increased anti-inflammatory/pro-inflammatory cytokine ratio. A genotypic effect on response to TNFα stimuli was detected, which was inhibited by vitamin D. Together our results show: (a) an altered gene expression in carriers of the susceptible CD28 variant, (b) no differences in protein levels on CD4+ T cells, and