Genome-wide Enrichment of De Novo Coding Mutations in Orofacial Cleft Trios.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Analysis of the genomic expression profile in trisomy 18: insight into possible genes involved in the associated phenotypes.

Trisomy 18, sometimes called Edwards syndrome, occurs in about 1 in 6000 live births and causes multiple birth defects in affected infants. The extra copy of chromosome 18 causes the altered expression of many genes and leads to severe skeletal, cardiovascular and neurological systems malformations as well as other medical problems. Due to the low rate of survival and the massive genetic imbalance, little research has been aimed at understanding the molecular consequences of trisomy 18 or considering potential therapeutic approaches. Our research is the first study to characterize whole-genome expression in fibroblast cells obtained from two patients with trisomy 18 and two matched controls, with follow-up expression confirmation studies on six independent controls. We show a detailed analysis of the most highly dysregulated genes on chromosome 18 and those genome-wide. The identified effector genes and the dysregulated downstream pathways provide hints of possible genotype-phenotype relationships to some of the most common symptoms observed in trisomy 18. We also provide a possible explanation for the sex-specific differences in survival, a unique characteristic of trisomy 18. Our analysis of genome-wide expression data moves us closer to understanding the molecular consequences of the second most common human autosomal trisomy of infants who survive to term. These insights might also translate to the understanding of the etiology of associated birth defects and medical conditions among those with trisomy 18.

Genetic contributors to risk of schizophrenia in the presence of a 22q11.2 deletion.

Schizophrenia occurs in about one in four individuals with 22q11.2 deletion syndrome (22q11.2DS). The aim of this International Brain and Behavior 22q11.2DS Consortium (IBBC) study was to identify genetic factors that contribute to schizophrenia, in addition to the ~20-fold increased risk conveyed by the 22q11.2 deletion. Using whole-genome sequencing data from 519 unrelated individuals with 22q11.2DS, we conducted genome-wide comparisons of common and rare variants between those with schizophrenia and those with no psychotic disorder at age ≥25 years. Available microarray data enabled direct comparison of polygenic risk for schizophrenia between 22q11.2DS and independent population samples with no 22q11.2 deletion, with and without schizophrenia (total n = 35,182). Polygenic risk for schizophrenia within 22q11.2DS was significantly greater for those with schizophrenia (padj = 6.73 × 10-6). Novel reciprocal case-control comparisons between the 22q11.2DS and population-based cohorts showed that polygenic risk score was significantly greater in individuals with psychotic illness, regardless of the presence of the 22q11.2 deletion. Within the 22q11.2DS cohort, results of gene-set analyses showed some support for rare variants affecting synaptic genes. No common or rare variants within the 22q11.2 deletion region were significantly associated with schizophrenia. These findings suggest that in addition to the deletion conferring a greatly increased risk to schizophrenia, the risk is higher when the 22q11.2 deletion and common polygenic risk factors that contribute to schizophrenia in the general population are both present.

Open Excision of Olecranon Melorheostosis to Relieve Elbow Extension Block.

We report a 34-year-old man who presented with hyperostosis of his right elbow associated with an inability to fully extend his elbow. The radiographic assessment revealed a classical dripping candle wax appearance of his proximal ulna suggestive of melorheostosis. Computed tomography was performed to identify the impingement point and aid in surgical planning. A targeted open excision biopsy via a Boyd incision was performed to excise the exophytic component that was causing the functional block. After surgery, he achieved full elbow extension and was able to return to his usual activity.

Study of telomere length in men who carry a fragile X premutation or full mutation allele.

The fragile X premutation is defined by the expansion of the CGG trinucleotide repeat at the 5' UTR of the FMR1 gene to between 55 and 200 repeats, while repeat tracks longer than 200 are defined as full mutations. Men carrying a premutation are at increased risk for fragile X-associated tremor/ataxia syndrome (FXTAS); those with > 200 repeats have fragile X syndrome, a common genetic form of intellectual disabilities. In our study, we tested the hypothesis that men carrying a fragile X premutation or full mutation are "biologically older", as suggested by the associated age-related disorder in the presence of the fragile X premutation or the altered cellular pathology that affects both the fragile X premutation and full mutation carriers. Thus, we predicted that both groups would have shorter telomeres than men carrying the normal size repeat allele. Using linear regression models, we found that, on average, premutation carriers had shorter telomeres compared with non-carriers (n = 69 vs n = 36; p = 0.02) and that there was no difference in telomere length between full mutation carriers and non-carriers (n = 37 vs n = 29; p > 0.10). Among premutation carriers only, we also asked whether telomere length was shorter among men with vs without symptoms of FXTAS (n = 28 vs n = 38 and n = 27 vs n = 41, depending on criteria) and found no evidence for a difference (p > 0.10). Previous studies have shown that the premutation is transcribed whereas the full mutation is not, and the expanded repeat track in FMR1 transcript is thought to lead to the risk for premutation-associated disorders. Thus, our data suggest that the observed premutation-only telomere shortening may be a consequence of the toxic effect of the premutation transcript and suggest that premutation carriers are "biologically older" than men carrying the normal size allele in the same age group.

Whole genome sequencing of orofacial cleft trios from the Gabriella Miller Kids First Pediatric Research Consortium identifies a new locus on chromosome 21.

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic, and the genetic etiology of non-syndromic OFCs is only partially determined. Here, we analyze whole genome sequence (WGS) data for association with risk of OFCs in European and Colombian families selected from a multicenter family-based OFC study. This is the first large-scale WGS study of OFC in parent-offspring trios, and a part of the Gabriella Miller Kids First Pediatric Research Program created for the study of childhood cancers and structural birth defects. WGS provides deeper and more specific genetic data than using imputation on present-day single nucleotide polymorphic (SNP) marker panels. Genotypes of case-parent trios at single nucleotide variants (SNV) and short insertions and deletions (indels) spanning the entire genome were called from their sequences using human GRCh38 genome assembly, and analyzed for association using the transmission disequilibrium test. Among genome-wide significant associations, we identified a new locus on chromosome 21 in Colombian families, not previously observed in other larger OFC samples of Latin American ancestry. This locus is situated within a region known to be expressed during craniofacial development. Based on deeper investigation of this locus, we concluded that it contributed risk for OFCs exclusively in the Colombians. This study reinforces the ancestry differences seen in the genetic etiology of OFCs, and underscores the need for larger samples when studying for OFCs and other birth defects in populations with diverse ancestry.

Identification of Polycystic Kidney Disease 1 Like 1 Gene Variants in Children With Biliary Atresia Splenic Malformation Syndrome.

Biliary atresia (BA) is the most common cause of end-stage liver disease in children and the primary indication for pediatric liver transplantation, yet underlying etiologies remain unknown. Approximately 10% of infants affected by BA exhibit various laterality defects (heterotaxy) including splenic abnormalities and complex cardiac malformations-a distinctive subgroup commonly referred to as the biliary atresia splenic malformation (BASM) syndrome. We hypothesized that genetic factors linking laterality features with the etiopathogenesis of BA in BASM patients could be identified through whole-exome sequencing (WES) of an affected cohort. DNA specimens from 67 BASM subjects, including 58 patient-parent trios, from the National Institute of Diabetes and Digestive and Kidney Diseases-supported Childhood Liver Disease Research Network (ChiLDReN) underwent WES. Candidate gene variants derived from a prespecified set of 2,016 genes associated with ciliary dysgenesis and/or dysfunction or cholestasis were prioritized according to pathogenicity, population frequency, and mode of inheritance. Five BASM subjects harbored rare and potentially deleterious biallelic variants in polycystic kidney disease 1 like 1 (PKD1L1), a gene associated with ciliary calcium signaling and embryonic laterality determination in fish, mice, and humans. Heterozygous PKD1L1 variants were found in 3 additional subjects. Immunohistochemical analysis of liver from the one BASM subject available revealed decreased PKD1L1 expression in bile duct epithelium when compared to normal livers and livers affected by other noncholestatic diseases. Conclusion: WES identified biallelic and heterozygous PKD1L1 variants of interest in 8 BASM subjects from the ChiLDReN data set; the dual roles for PKD1L1 in laterality determination and ciliary function suggest that PKD1L1 is a biologically plausible, cholangiocyte-expressed candidate gene for the BASM syndrome.

PEMapper and PECaller provide a simplified approach to whole-genome sequencing.

The analysis of human whole-genome sequencing data presents significant computational challenges. The sheer size of datasets places an enormous burden on computational, disk array, and network resources. Here, we present an integrated computational package, PEMapper/PECaller, that was designed specifically to minimize the burden on networks and disk arrays, create output files that are minimal in size, and run in a highly computationally efficient way, with the single goal of enabling whole-genome sequencing at scale. In addition to improved computational efficiency, we implement a statistical framework that allows for a base by base error model, allowing this package to perform as well or better than the widely used Genome Analysis Toolkit (GATK) in all key measures of performance on human whole-genome sequences.

Identification and Repurposing of Trisubstituted Harmine Derivatives as Novel Inhibitors of Mycobacterium tuberculosis Phosphoserine Phosphatase.

Mycobacterium tuberculosis is still the deadliest bacterial pathogen worldwide and the increasing number of multidrug-resistant tuberculosis cases further complicates this global health issue. M. tuberculosis phosphoserine phosphatase SerB2 is a promising target for drug design. Besides being a key essential metabolic enzyme of the pathogen's serine pathway, it appears to be involved in immune evasion mechanisms. In this work, a malachite green-based phosphatase assay has been used to screen 122 compounds from an internal chemolibrary. Trisubstituted harmine derivatives were found among the best hits that inhibited SerB2 activity. Synthesis of an original compound helped to discuss a brief structure activity relationship evaluation. Kinetics experiments showed that the most potent derivatives inhibit the phosphatase in a parabolic competitive fashion with apparent inhibition constants ( K i ) values in the micromolar range. Their interaction modes with the enzyme were investigated through induced fit docking experiments, leading to results consistent with the experimental data. Cellular assays showed that the selected compounds also inhibited M. tuberculosis growth in vitro. Those promising results may provide a basis for the development of new antimycobacterial agents targeting SerB2.

Genome-Wide Association Study Identifies African-Specific Susceptibility Loci in African Americans With Inflammatory Bowel Disease.

BACKGROUND & AIMS: The inflammatory bowel diseases (IBD) ulcerative colitis (UC) and Crohn's disease (CD) cause significant morbidity and are increasing in prevalence among all populations, including African Americans. More than 200 susceptibility loci have been identified in populations of predominantly European ancestry, but few loci have been associated with IBD in other ethnicities. METHODS: We performed 2 high-density, genome-wide scans comprising 2345 cases of African Americans with IBD (1646 with CD, 583 with UC, and 116 inflammatory bowel disease unclassified) and 5002 individuals without IBD (controls, identified from the Health Retirement Study and Kaiser Permanente database). Single-nucleotide polymorphisms (SNPs) associated at P < 5.0 × 10-8 in meta-analysis with a nominal evidence (P < .05) in each scan were considered to have genome-wide significance. RESULTS: We detected SNPs at HLA-DRB1, and African-specific SNPs at ZNF649 and LSAMP, with associations of genome-wide significance for UC. We detected SNPs at USP25 with associations of genome-wide significance for IBD. No associations of genome-wide significance were detected for CD. In addition, 9 genes previously associated with IBD contained SNPs with significant evidence for replication (P < 1.6 × 10-6): ADCY3, CXCR6, HLA-DRB1 to HLA-DQA1 (genome-wide significance on conditioning), IL12B,PTGER4, and TNC for IBD; IL23R, PTGER4, and SNX20 (in strong linkage disequilibrium with NOD2) for CD; and KCNQ2 (near TNFRSF6B) for UC. Several of these genes, such as TNC (near TNFSF15), CXCR6, and genes associated with IBD at the HLA locus, contained SNPs with unique association patterns with African-specific alleles. CONCLUSIONS: We performed a genome-wide association study of African Americans with IBD and identified loci associated with UC in only this population; we also replicated IBD, CD, and UC loci identified in European populations. The detection of variants associated with IBD risk in only people of African descent demonstrates the importance of studying the genetics of IBD and other complex diseases in populations beyond those of European ancestry.

Identification of consensus binding sites clarifies FMRP binding determinants.

Fragile X mental retardation protein (FMRP) is a multifunctional RNA-binding protein with crucial roles in neuronal development and function. Efforts aimed at elucidating how FMRP target mRNAs are selected have produced divergent sets of target mRNA and putative FMRP-bound motifs, and a clear understanding of FMRP's binding determinants has been lacking. To clarify FMRP's binding to its target mRNAs, we produced a shared dataset of FMRP consensus binding sequences (FCBS), which were reproducibly identified in two published FMRP CLIP sequencing datasets. This comparative dataset revealed that of the various sequence and structural motifs that have been proposed to specify FMRP binding, the short sequence motifs TGGA and GAC were corroborated, and a novel TAY motif was identified. In addition, the distribution of the FCBS set demonstrates that FMRP preferentially binds to the coding region of its targets but also revealed binding along 3' UTRs in a subset of target mRNAs. Beyond probing these putative motifs, the FCBS dataset of reproducibly identified FMRP binding sites is a valuable tool for investigating FMRP targets and function.

Sex-specific hippocampal 5-hydroxymethylcytosine is disrupted in response to acute stress.

Environmental stress is among the most important contributors to increased susceptibility to develop psychiatric disorders. While it is well known that acute environmental stress alters gene expression, the molecular mechanisms underlying these changes remain largely unknown. 5-hydroxymethylcytosine (5hmC) is a novel environmentally sensitive epigenetic modification that is highly enriched in neurons and is associated with active neuronal transcription. Recently, we reported a genome-wide disruption of hippocampal 5hmC in male mice following acute stress that was correlated to altered transcript levels of genes in known stress related pathways. Since sex-specific endocrine mechanisms respond to environmental stimulus by altering the neuronal epigenome, we examined the genome-wide profile of hippocampal 5hmC in female mice following exposure to acute stress and identified 363 differentially hydroxymethylated regions (DhMRs) linked to known (e.g., Nr3c1 and Ntrk2) and potentially novel genes associated with stress response and psychiatric disorders. Integration of hippocampal expression data from the same female mice found stress-related hydroxymethylation correlated to altered transcript levels. Finally, characterization of stress-induced sex-specific 5hmC profiles in the hippocampus revealed 778 sex-specific acute stress-induced DhMRs some of which were correlated to altered transcript levels that produce sex-specific isoforms in response to stress. Together, the alterations in 5hmC presented here provide a possible molecular mechanism for the adaptive sex-specific response to stress that may augment the design of novel therapeutic agents that will have optimal effectiveness in each sex.

Genome-wide alterations in hippocampal 5-hydroxymethylcytosine links plasticity genes to acute stress.

Environmental stress is among the most important contributors to increased susceptibility to develop psychiatric disorders, including anxiety and post-traumatic stress disorder. While even acute stress alters gene expression, the molecular mechanisms underlying these changes remain largely unknown. 5-hydroxymethylcytosine (5hmC) is a novel environmentally sensitive DNA modification that is highly enriched in post-mitotic neurons and is associated with active transcription of neuronal genes. Recently, we found a hippocampal increase of 5hmC in the glucocorticoid receptor gene (Nr3c1) following acute stress, warranting a deeper investigation of stress-related 5hmC levels. Here we used an established chemical labeling and affinity purification method coupled with high-throughput sequencing technology to generate the first genome-wide profile of hippocampal 5hmC following exposure to acute restraint stress and a one-hour recovery. This approach found a genome-wide disruption in 5hmC associated with acute stress response, primarily in genic regions, and identified known and potentially novel stress-related targets that have a significant enrichment for neuronal ontological functions. Integration of these data with hippocampal gene expression data from these same mice found stress-related hydroxymethylation correlated to altered transcript levels and sequence motif predictions indicated that 5hmC may function by mediating transcription factor binding to these transcripts. Together, these data reveal an environmental impact on this newly discovered epigenetic mark in the brain and represent a critical step toward understanding stress-related epigenetic mechanisms that alter gene expression and can lead to the development of psychiatric disorders.

Analysis of FMRP mRNA target datasets reveals highly associated mRNAs mediated by G-quadruplex structures formed via clustered WGGA sequences.

Fragile X syndrome, a common cause of intellectual disability and a well-known cause of autism spectrum disorder, is the result of loss or dysfunction of fragile X mental retardation protein (FMRP), a highly selective RNA-binding protein and translation regulator. A major research priority has been the identification of the mRNA targets of FMRP, particularly as recent studies suggest an excess of FMRP targets among genes implicated in idiopathic autism and schizophrenia. Several large-scale studies have attempted to identify mRNAs bound by FMRP through several methods, each generating a list of putative target genes, leading to distinct hypotheses by which FMRP recognizes its targets; namely, by RNA structure or sequence. However, no in depth analyses have been performed to identify the level of consensus among the studies. Here, we analyze four large FMRP target datasets to generate high-confidence consensus lists, and examine all datasets for sequence elements within the target RNAs to validate reported FMRP binding motifs (GACR, ACUK and WGGA). We found GACR to be highly enriched in FMRP datasets, while ACUK was not. The WGGA pattern was modestly enriched in several, but not all datasets. The previous association between FMRP and G-quadruplexes prompted the analysis of the distribution of WGGA in the target genes. Consistent with the requirements for G-quadruplex formation, we observed highly clustered WGGA motifs in FMRP targets compared with other genes, implicating both RNA structure and sequence in the recognition motif of FMRP. In addition, we generate a list of the top 40 FMRP targets associated with FXS-related phenotypes.

Characterization of genetic loci that affect susceptibility to inflammatory bowel diseases in African Americans.

BACKGROUND & AIMS: Inflammatory bowel disease (IBD) has familial aggregation in African Americans (AAs), but little is known about the molecular genetic susceptibility. Mapping studies using the Immunochip genotyping array expand the number of susceptibility loci for IBD in Caucasians to 163, but the contribution of the 163 loci and European admixture to IBD risk in AAs is unclear. We performed a genetic mapping study using the Immunochip to determine whether IBD susceptibility loci in Caucasians also affect risk in AAs and identify new associated loci. METHODS: We recruited AAs with IBD and without IBD (controls) from 34 IBD centers in the United States; additional controls were collected from 4 other Immunochip studies. Association and admixture loci were mapped for 1088 patients with Crohn's disease, 361 with ulcerative colitis, 62 with IBD type unknown, and 1797 controls; 130,241 autosomal single-nucleotide polymorphisms (SNPs) were analyzed. RESULTS: The strongest associations were observed between ulcerative colitis and HLA rs9271366 (P = 7.5 × 10(-6)), Crohn's disease and 5p13.1 rs4286721 (P = 3.5 × 10(-6)), and IBD and KAT2A rs730086 (P = 2.3 × 10(-6)). Additional suggestive associations (P < 4.2 × 10(-5)) were observed between Crohn's disease and IBD and African-specific SNPs in STAT5A and STAT3; between IBD and SNPs in IL23R, IL12B, and C2orf43; and between ulcerative colitis and SNPs near HDAC11 and near LINC00994. The latter 3 loci have not been previously associated with IBD, but require replication. Established Caucasian associations were replicated in AAs (P < 3.1 × 10(-4)) at NOD2, IL23R, 5p15.3, and IKZF3. Significant admixture (P < 3.9 × 10(-4)) was observed for 17q12-17q21.31 (IZKF3 through STAT3), 10q11.23-10q21.2, 15q22.2-15q23, and 16p12.2-16p12.1. Network analyses showed significant enrichment (false discovery rate <1 × 10(-5)) in genes that encode members of the JAK-STAT, cytokine, and chemokine signaling pathways, as well those involved in pathogenesis of measles. CONCLUSIONS: In a genetic analysis of 3308 AA IBD cases and controls, we found that many variants associated with IBD in Caucasians also showed association evidence with these diseases in AAs; we also found evidence for variants and loci not previously associated with IBD. The complex genetic factors that determine risk for or protection against IBD in different populations require further study.

Differentially methylated plasticity genes in the amygdala of young primates are linked to anxious temperament, an at risk phenotype for anxiety and depressive disorders.

Children with an anxious temperament (AT) are at a substantially increased risk to develop anxiety and depression. The young rhesus monkey is ideal for studying the origin of human AT because it shares with humans the genetic, neural, and phenotypic underpinnings of complex social and emotional functioning. Heritability, functional imaging, and gene expression studies of AT in young monkeys revealed that the central nucleus of the amygdala (Ce) is a key environmentally sensitive substrate of this at risk phenotype. Because epigenetic marks (e.g., DNA methylation) can be modulated by environmental stimuli, these data led us to hypothesize a role for DNA methylation in the development of AT. To test this hypothesis, we used reduced representation bisulfite sequencing to examine the cross-sectional genome-wide methylation levels in the Ce of 23 age-matched monkeys (1.3 ± 0.2 years) phenotyped for AT. Because AT reflects a continuous trait-like variable, we used an analytical approach that is consistent with this biology to identify genes in the Ce with methylation patterns that predict AT. Expression data from the Ce of these same monkeys were then used to find differentially methylated candidates linked to altered gene regulation. Two genes particularly relevant to the AT phenotype were BCL11A and JAG1. These transcripts have well-defined roles in neurodevelopmental processes, including neurite arborization and the regulation of neurogenesis. Together, these findings represent a critical step toward understanding the effects of early environment on the neuromolecular mechanisms that underlie the risk to develop anxiety and depressive disorders.

Age-associated DNA methylation in pediatric populations.

DNA methylation (DNAm) plays diverse roles in human biology, but this dynamic epigenetic mark remains far from fully characterized. Although earlier studies uncovered loci that undergo age-associated DNAm changes in adults, little is known about such changes during childhood. Despite profound DNAm plasticity during embryogenesis, monozygotic twins show indistinguishable childhood methylation, suggesting that DNAm is highly coordinated throughout early development. Here we examine the methylation of 27,578 CpG dinucleotides in peripheral blood DNA from a cross-sectional study of 398 boys, aged 3-17 yr, and find significant age-associated changes in DNAm at 2078 loci. These findings correspond well with pyrosequencing data and replicate in a second pediatric population (N = 78). Moreover, we report a deficit of age-related loci on the X chromosome, a preference for specific nucleotides immediately surrounding the interrogated CpG dinucleotide, and a primary association with developmental and immune ontological functions. Meta-analysis (N = 1158) with two adult populations reveals that despite a significant overlap of age-associated loci, most methylation changes do not follow a lifelong linear pattern due to a threefold to fourfold higher rate of change in children compared with adults; consequently, the vast majority of changes are more accurately modeled as a function of logarithmic age. We therefore conclude that age-related DNAm changes in peripheral blood occur more rapidly during childhood and are imperfectly accounted for by statistical corrections that are linear in age, further suggesting that future DNAm studies should be matched closely for age.

Array-based assay detects genome-wide 5-mC and 5-hmC in the brains of humans, non-human primates, and mice.

BACKGROUND: Methylation on the fifth position of cytosine (5-mC) is an essential epigenetic mark that is linked to both normal neurodevelopment and neurological diseases. The recent identification of another modified form of cytosine, 5-hydroxymethylcytosine (5-hmC), in both stem cells and post-mitotic neurons, raises new questions as to the role of this base in mediating epigenetic effects. Genomic studies of these marks using model systems are limited, particularly with array-based tools, because the standard method of detecting DNA methylation cannot distinguish between 5-mC and 5-hmC and most methods have been developed to only survey the human genome. RESULTS: We show that non-human data generated using the optimization of a widely used human DNA methylation array, designed only to detect 5-mC, reproducibly distinguishes tissue types within and between chimpanzee, rhesus, and mouse, with correlations near the human DNA level (R(2) > 0.99). Genome-wide methylation analysis, using this approach, reveals 6,102 differentially methylated loci between rhesus placental and fetal tissues with pathways analysis significantly overrepresented for developmental processes. Restricting the analysis to oncogenes and tumor suppressor genes finds 76 differentially methylated loci, suggesting that rhesus placental tissue carries a cancer epigenetic signature. Similarly, adapting the assay to detect 5-hmC finds highly reproducible 5-hmC levels within human, rhesus, and mouse brain tissue that is species-specific with a hierarchical abundance among the three species (human > rhesus >> mouse). Annotation of 5-hmC with respect to gene structure reveals a significant prevalence in the 3'UTR and an association with chromatin-related ontological terms, suggesting an epigenetic feedback loop mechanism for 5-hmC. CONCLUSIONS: Together, these data show that this array-based methylation assay is generalizable to all mammals for the detection of both 5-mC and 5-hmC, greatly improving the utility of mammalian model systems to study the role of epigenetics in human health, disease, and evolution.

Deciphering the role of VapBC13 and VapBC26 toxin antitoxin systems in the pathophysiology of Mycobacterium tuberculosis.

The expansion of VapBC TA systems in M. tuberculosis has been linked with its fitness and survival upon exposure to stress conditions. Here, we have functionally characterized VapBC13 and VapBC26 TA modules of M. tuberculosis. We report that overexpression of VapC13 and VapC26 toxins in M. tuberculosis results in growth inhibition and transcriptional reprogramming. We have also identified various regulatory proteins as hub nodes in the top response network of VapC13 and VapC26 overexpression strains. Further, analysis of RNA protection ratios revealed potential tRNA targets for VapC13 and VapC26. Using in vitro ribonuclease assays, we demonstrate that VapC13 and VapC26 degrade serT and leuW tRNA, respectively. However, no significant changes in rRNA cleavage profiles were observed upon overexpression of VapC13 and VapC26 in M. tuberculosis. In order to delineate the role of these TA systems in M. tuberculosis physiology, various mutant strains were constructed. We show that in comparison to the parental strain, ΔvapBC13 and ΔvapBC26 strains were mildly susceptible to oxidative stress. Surprisingly, the growth patterns of parental and mutant strains were comparable in aerosol-infected guinea pigs. These observations imply that significant functional redundancy exists for some TA systems from M. tuberculosis.