Long-read sequencing in deciphering human genetics to a greater depth.

Through four decades' development, DNA sequencing has inched into the era of single-molecule sequencing (SMS), or the third-generation sequencing (TGS), as represented by two distinct technical approaches developed independently by Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT). Historically, each generation of sequencing technologies was marked by innovative technological achievements and novel applications. Long reads (LRs) are considered as the most advantageous feature of SMS shared by both PacBio and ONT to distinguish SMS from next-generation sequencing (NGS, or the second-generation sequencing) and Sanger sequencing (the first-generation sequencing). Long reads overcome the limitations of NGS and drastically improves the quality of genome assembly. Besides, ONT also contributes several unique features including ultra-long reads (ULRs) with read length above 300 kb and some close to 1 million bp, direct RNA sequencing and superior portability as made possible by pocket-sized MinION sequencer. Here, we review the history of DNA sequencing technologies and associated applications, with a special focus on the advantages as well as the limitations of ULR sequencing in genome assembly.

Nuclear receptor corepressors non-canonically drive glucocorticoid receptor-dependent activation of hepatic gluconeogenesis.

Nuclear receptor corepressors (NCoRs) function in multiprotein complexes containing histone deacetylase 3 (HDAC3) to alter transcriptional output primarily through repressive chromatin remodelling at target loci1-5. In the liver, loss of HDAC3 causes a marked hepatosteatosis largely because of de-repression of genes involved in lipid metabolism6,7; however, the individual roles and contribution of other complex members to hepatic and systemic metabolic regulation are unclear. Here we show that adult loss of both NCoR1 and NCoR2 (double knockout (KO)) in hepatocytes phenocopied the hepatomegalic fatty liver phenotype of HDAC3 KO. In addition, double KO livers exhibited a dramatic reduction in glycogen storage and gluconeogenic gene expression that was not observed with hepatic KO of individual NCoRs or HDAC3, resulting in profound fasting hypoglycaemia. This surprising HDAC3-independent activation function of NCoR1 and NCoR2 is due to an unexpected loss of chromatin accessibility on deletion of NCoRs that prevented glucocorticoid receptor binding and stimulatory effect on gluconeogenic genes. These studies reveal an unanticipated, non-canonical activation function of NCoRs that is required for metabolic health.

Short-term cold exposure induces persistent epigenomic memory in brown fat.

Deficiency of the epigenome modulator histone deacetylase 3 (HDAC3) in brown adipose tissue (BAT) impairs the ability of mice to survive in near-freezing temperatures. Here, we report that short-term exposure to mild cold temperature (STEMCT: 15°C for 24 h) averted lethal hypothermia of mice lacking HDAC3 in BAT (HDAC3 BAT KO) exposed to 4°C. STEMCT restored the induction of the thermogenic coactivator PGC-1α along with UCP1 at 22°C, which is greatly impaired in HDAC3-deficient BAT, and deletion of either UCP1 or PGC-1α prevented the protective effect of STEMCT. Remarkably, this protection lasted for up to 7 days. Transcriptional activator C/EBPß was induced by short-term cold exposure in mouse and human BAT and, uniquely, remained high for 7 days following STEMCT. Adeno-associated virus-mediated knockdown of BAT C/EBPß in HDAC3 BAT KO mice erased the persistent memory of STEMCT, revealing the existence of a C/EBPß-dependent and HDAC3-independent cold-adaptive epigenomic memory.

Hepatic Vagal Afferents Convey Clock-Dependent Signals to Regulate Circadian Food Intake.

Circadian desynchrony induced by shiftwork or jetlag is detrimental to metabolic health, but how synchronous/desynchronous signals are transmitted among tissues is unknown. Here we report that liver molecular clock dysfunction is signaled to the brain via the hepatic vagal afferent nerve (HVAN), leading to altered food intake patterns that are corrected by ablation of the HVAN. Hepatic branch vagotomy also prevents food intake disruptions induced by high-fat diet feeding and reduces body weight gain. Our findings reveal a previously unrecognized homeostatic feedback signal that relies on synchrony between the liver and the brain to control circadian food intake patterns. This identifies the hepatic vagus nerve as a therapeutic target for obesity in the setting of chrono-disruption. One Sentence Summary: The hepatic vagal afferent nerve signals internal circadian desynchrony between the brain and liver to induce maladaptive food intake patterns.

Comprehensive Cohort Analysis of Mutational Spectrum in Early Onset Breast Cancer Patients.

Early onset breast cancer (EOBC), diagnosed at age ~40 or younger, is associated with a poorer prognosis and higher mortality rate compared to breast cancer diagnosed at age 50 or older. EOBC poses a serious threat to public health and requires in-depth investigation. We studied a cohort comprising 90 Taiwanese female patients, aiming to unravel the underlying mechanisms of EOBC etiopathogenesis. Sequence data generated by whole-exome sequencing (WES) and whole-genome sequencing (WGS) from white blood cell (WBC)-tumor pairs were analyzed to identify somatic missense mutations, copy number variations (CNVs) and germline missense mutations. Similar to regular breast cancer, the key somatic mutation-susceptibility genes of EOBC include TP53 (40% prevalence), PIK3CA (37%), GATA3 (17%) and KMT2C (17%), which are frequently reported in breast cancer; however, the structural protein-coding genes MUC17 (19%), FLG (16%) and NEBL (11%) show a significantly higher prevalence in EOBC. Furthermore, the top 2 genes harboring EOBC germline mutations, MUC16 (19%) and KRT18 (19%), encode structural proteins. Compared to conventional breast cancer, an unexpectedly higher number of EOBC susceptibility genes encode structural proteins. We suspect that mutations in structural proteins may increase physical permeability to environmental hormones and carcinogens and cause breast cancer to occur at a young age.

Single cell transcriptome analysis of MCF-7 reveals consistently and inconsistently expressed gene groups each associated with distinct cellular localization and functions.

Single cell transcriptome (SCT) analysis provides superior resolution to illustrate tumor cell heterogeneity for clinical implications. We characterized four SCTs of MCF-7 using 143 housekeeping genes (HKGs) as control, of which lactate dehydrogenase B (LDHB) expression is silenced. These SCT libraries mapped to 11,423, 11,486, 10,380, and 11,306 RefSeq genes (UCSC), respectively. High consistency in HKG expression levels across all four SCTs, along with transcriptional silencing of LDHB, was observed, suggesting a high sensitivity and reproducibility of the SCT analysis. Cross-library comparison on expression levels by scatter plotting revealed a linear correlation and an 83-94% overlap in transcript isoforms and expressed genes were also observed. To gain insight of transcriptional diversity among the SCTs, expressed genes were split into consistently expressed (CE) (expressed in all SCTs) and inconsistently expressed (IE) (expressed in some but not all SCTs) genes for further characterization, along with the 142 expressed HKGs as a reference. Distinct transcriptional strengths were found among these groups, with averages of 1,612.0, 88.0 and 1.2 FPKM for HKGs, CE and IE, respectively. Comparison between CE and IE groups further indicated that expressions of CE genes vary more significantly than that of IE genes. Gene Ontology analysis indicated that proteins encoded by CE genes are mainly involved in fundamental intracellular activities, while proteins encoded by IE genes are mainly for extracellular activities, especially acting as receptors or ion channels. The diversified gene expressions, especially for those encoded by IE genes, may contribute to cancer drug resistance.

Analysis of microbial sequences in plasma cell-free DNA for early-onset breast cancer patients and healthy females.

BACKGROUND: Cell-free circulating DNA (cfDNA) is becoming a useful biopsy for noninvasive diagnosis of diseases. Microbial sequences in plasma cfDNA may provide important information to improve prognosis and treatment. We have developed a stringent method to identify microbial species via microbial cfDNA in the blood plasma of early-onset breast cancer (EOBC) patients and healthy females. Empirically, microbe-originated sequence reads were identified by mapping non-human PE reads in cfDNA libraries to microbial databases. Those mapped concordantly to unique microbial species were assembled into contigs, which were subsequently aligned to the same databases. Microbial species uniquely aligned were identified and compared across all individuals on MCRPM (Microbial CfDNA Reads Per Million quality PE reads) basis. RESULTS: The predominant microbial cfDNAs in all plasma samples examined are originated from bacteria and these bacteria were limited to only a few genera. Among those, Acinetobacter johnsonii XBB1 and low levels of Mycobacterium spp. were commonly found in all healthy females, but also present in an EOBC patient. Compared to those in healthy counterparts, bacterial species in EOBC patients are more diverse and more likely to present at high levels. Among these three EOBC patients tested, a patient who has record high titer (2,724 MCRPM) of Pseudomonas mendocina together with 8.82 MCRPM of Pannonibacter phragmitetus has passed away; another patient infected by multiple Sphingomonas species remains alive; while the third patient who has similar microbial species (Acinetobacter johnsonii XBB1) commonly seen in normal controls is having a normal life. CONCLUSIONS: Our preliminary data on the profiles of microbial cfDNA sequences suggested that it may have some prognostic value in cancer patients. Validation in larger number of patients is warranted.

T Oligo-Primed Polymerase Chain Reaction (TOP-PCR): A Robust Method for the Amplification of Minute DNA Fragments in Body Fluids.

Body fluid DNA sequencing is a powerful noninvasive approach for the diagnosis of genetic defects, infectious agents and diseases. The success relies on the quantity and quality of the DNA samples. However, numerous clinical samples are either at low quantity or of poor quality due to various reasons. To overcome these problems, we have developed T oligo-primed polymerase chain reaction (TOP-PCR) for full-length nonselective amplification of minute quantity of DNA fragments. TOP-PCR adopts homogeneous "half adaptor" (HA), generated by annealing P oligo (carrying a phosphate group at the 5' end) and T oligo (carrying a T-tail at the 3' end), for efficient ligation to target DNA and subsequent PCR amplification primed by the T oligo alone. Using DNA samples from body fluids, we demonstrate that TOP-PCR recovers minute DNA fragments and maintains the DNA size profile, while enhancing the major molecular populations. Our results also showed that TOP-PCR is a superior method for detecting apoptosis and outperforms the method adopted by Illumina for DNA amplification.

ClosIndb: A resource for computationally derived information from clostridial genomes.

Over the past few years, several clostridial genomes have been sequenced, and since then new sequencing projects are also under way. Clostridia is one of the most sequenced genera, and presently, complete genome sequences of 49 clostridial species are available in public archives. Unraveling this wealth of genomic information opens up potential avenues in clostridial research. In the present study, we have carried out in silico analysis to decipher the genomic data. Subsequently, a web resource, ClosIndb, has been developed which collates the computationally derived information associated with all clostridial genes. It features various aspects of coding regions as well as non-coding regions, such as putative orthologs, proteins physicochemical properties, operons and cis-regulatory elements. It provides users with comparative details of all clostridial proteins across the firmicutes. ClosIndb is a comprehensive resource for all completely sequenced clostridial genomes and is under constant development. ClosIndb is freely accessible at http://bif.uohyd.ac.in/closindb/.

Complete Taiwanese Macaque (Macaca cyclopis) Mitochondrial Genome: Reference-Assisted de novo Assembly with Multiple k-mer Strategy.

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.

MycoRRdb: a database of computationally identified regulatory regions within intergenic sequences in mycobacterial genomes.

The identification of regulatory regions for a gene is an important step towards deciphering the gene regulation. Regulatory regions tend to be conserved under evolution that facilitates the application of comparative genomics to identify such regions. The present study is an attempt to make use of this attribute to identify regulatory regions in the Mycobacterium species followed by the development of a database, MycoRRdb. It consist the regulatory regions identified within the intergenic distances of 25 mycobacterial species. MycoRRdb allows to retrieve the identified intergenic regulatory elements in the mycobacterial genomes. In addition to the predicted motifs, it also allows user to retrieve the Reciprocal Best BLAST Hits across the mycobacterial genomes. It is a useful resource to understand the transcriptional regulatory mechanism of mycobacterial species. This database is first of its kind which specifically addresses cis-regulatory regions and also comprehensive to the mycobacterial species. Database URL: http://mycorrdb.uohbif.in.

OrFin: A web tool for detection of putative orthologs.

UNLABELLED: Identification of ortholog is one of the important tasks to understand a novel genome. It helps to assign functional annotations, from one organism to another organism. To identify the putative ortholog, Reciprocal Best BLAST hit (RBBH) method is known to be an efficient approach. OrFin makes use of the same approach to identify pair of orthologous proteins for a given set of sequences of two species. It is a user-friendly web tool which works with user defined parameters to search RBBHs. Results are produced in both html and text format. AVAILABILITY: This web tool is freely available at http://bifl.uohyd.ac.in/orfin.