Background DNA methylation is thought to be an important determinant of human phenotypic variation, but its inherent cell type specificity has impeded progress on this question. At exceptional genomic regions, interindividual variation in DNA methylation occurs systemically. Like genetic variants, systemic interindividual epigenetic variants are stable, can influence phenotype, and can be assessed in any easily biopsiable DNA sample. We describe an unbiased screen for human genomic regions at which interindividual variation in DNA methylation is not tissue-specific. Results For each of 10 donors from the NIH Genotype-Tissue Expression (GTEx) program, CpG methylation is measured by deep whole-genome bisulfite sequencing of genomic DNA from tissues representing the three germ layer lineages: thyroid (endoderm), heart (mesoderm), and brain (ectoderm). We develop a computational algorithm to identify genomic regions at which interindividual variation in DNA methylation is consistent across all three lineages. This approach identifies 9926 correlated regions of systemic interindividual variation (CoRSIVs). These regions, comprising just 0.1% of the human genome, are inter-correlated over long genomic distances, associated with transposable elements and subtelomeric regions, conserved across diverse human ethnic groups, sensitive to periconceptional environment, and associated with genes implicated in a broad range of human disorders and phenotypes. CoRSIV methylation in one tissue can predict expression of associated genes in other tissues. Conclusions In addition to charting a previously unexplored molecular level of human individuality, this atlas of human CoRSIVs provides a resource for future population-based investigations into how interindividual epigenetic variation modulates risk of disease. Electronic supplementary material The online version of this article (10.1186/s13059-019-1708-1) contains supplementary material, which is available to authorized users.

The birth prevalence of sex chromosome trisomies (SCT), that is individuals with an XYY, XXY or XXX sex chromosome constitution, is traditionally based on six surveys of unselected newborns carried out in the 1960s and early 1970s. All three SCTs had a prevalence of 1 in 1000 same sex births. We re-examined these prevalences based on additional cytogenetic studies of newborn surveys, spontaneous abortions, perinatal deaths and prenatal diagnoses. The more recent newborn surveys suggest there has been an increase in the prevalence of XXYs, but not of the other two SCTs since the original newborn series. The prevalence of XXYs has risen from 1.09 to 1.72 per 1000 male births (P=0.023). We suggest that such an increase, in the absence of an increase in the prevalence of XXX, is unlikely to be due to increased maternal age. As XXY is the only chromosome abnormality known where a substantial proportion ( approximately 50%) arise as the result of non-disjunction at the first paternal meiotic division, we speculate that some factor may be interfering with pairing and/or recombination of the sex bivalent at the paternal MI division.

Bardet Biedl syndrome (BBS) is a multisystem genetically heterogeneous ciliopathy that most commonly leads to obesity, photoreceptor degeneration, digit anomalies, genito-urinary abnormalities, as well as cognitive impairment with autism, among other features. Sequencing of a DNA sample from a 17-year-old female affected with BBS did not identify any mutation in the known BBS genes. Whole-genome sequencing identified a novel loss-of-function disease-causing homozygous mutation (K102*) in C8ORF37, a gene coding for a cilia protein. The proband was overweight (body mass index 29.1) with a slowly progressive rod-cone dystrophy, a mild learning difficulty, high myopia, three limb post-axial polydactyly, horseshoe kidney, abnormally positioned uterus and elevated liver enzymes. Mutations in C8ORF37 were previously associated with severe autosomal recessive retinal dystrophies (retinitis pigmentosa RP64 and cone-rod dystrophy CORD16) but not BBS. To elucidate the functional role of C8ORF37 in a vertebrate system, we performed gene knockdown in Danio rerio and assessed the cardinal features of BBS and visual function. Knockdown of c8orf37 resulted in impaired visual behavior and BBS-related phenotypes, specifically, defects in the formation of Kupffer's vesicle and delays in retrograde transport. Specificity of these phenotypes to BBS knockdown was shown with rescue experiments. Over-expression of human missense mutations in zebrafish also resulted in impaired visual behavior and BBS-related phenotypes. This is the first functional validation and association of C8ORF37 mutations with the BBS phenotype, which identifies BBS21. The zebrafish studies hereby show that C8ORF37 variants underlie clinically diagnosed BBS-related phenotypes as well as isolated retinal degeneration.

S-palmitoylation is a reversible posttranslational modification that plays an important role in regulating protein localization, trafficking, and stability. Recent studies have shown that some proteins undergo extremely rapid palmitoylation/depalmitoylation cycles after cellular stimulation supporting a direct signaling role for this posttranslational modification. Here, we investigated whether β -adrenergic stimulation of cardiomyocytes led to stimulus-dependent palmitoylation of downstream signaling proteins. We found that β -adrenergic stimulation led to rapidly increased G α s and G α i palmitoylation. The kinetics of palmitoylation was temporally consistent with the downstream production of cAMP and contractile responses. We identified the plasma membrane-localized palmitoyl acyltransferase DHHC5 as an important mediator of the stimulus-dependent palmitoylation in cardiomyocytes. Knockdown of DHHC5 showed that this enzyme is necessary for palmitoylation of G α s, G α i, and functional responses downstream of β -adrenergic stimulation. A palmitoylation assay with purified components revealed that G α s and G α i are direct substrates of DHHC5. Finally, we provided evidence that the C-terminal tail of DHHC5 can be palmitoylated in response to stimulation and such modification is important for its dynamic localization and function in the plasma membrane. Our results reveal that DHHC5 is a central regulator of signaling downstream of β -adrenergic receptors in cardiomyocytes.

Genome sequencing of the protistan parasite Entamoeba histolytica HM-1:IMSS revealed that almost all the tRNA genes are organized into tandem arrays that make up over 10% of the genome. The 25 distinct array units contain up to 5 tRNA genes each and some also encode the 5S RNA. Between adjacent genes in array units are complex short tandem repeats (STRs) resembling microsatellites. To investigate the origins and evolution of this unique gene organization, we have undertaken a genome survey to determine the array unit organization in 4 other species of Entamoeba—Entamoeba dispar, Entamoeba moshkovskii, Entamoeba terrapinae, and Entamoeba invadens—and have explored the STR structure in other isolates of E. histolytica. The genome surveys revealed that E. dispar has the same array unit organization as E. histolytica, including the presence and numerical variation of STRs between adjacent genes. However, the individual repeat sequences are completely different to those in E. histolytica. All other species of Entamoeba studied also have tandem arrays of clustered tRNA genes, but the gene composition of the array units often differs from that in E. histolytica/E. dispar. None of the other species' arrays exhibit the complex STRs between adjacent genes although simple tandem duplications are occasionally seen. The degree of similarity in organization reflects the phylogenetic relationships among the species studied. Within individual isolates of E. histolytica most copies of the array unit are uniform in sequence with only minor variation in the number and organization of the STRs. Between isolates, however, substantial differences in STR number and organization can exist although the individual repeat sequences tend to be conserved. The origin of this unique gene organization in the genus Entamoeba clearly predates the common ancestor of the species investigated to date and their function remains unclear.