Polygenic prediction of preeclampsia and gestational hypertension.

Preeclampsia and gestational hypertension are common pregnancy complications associated with adverse maternal and child outcomes. Current tools for prediction, prevention and treatment are limited. Here we tested the association of maternal DNA sequence variants with preeclampsia in 20,064 cases and 703,117 control individuals and with gestational hypertension in 11,027 cases and 412,788 control individuals across discovery and follow-up cohorts using multi-ancestry meta-analysis. Altogether, we identified 18 independent loci associated with preeclampsia/eclampsia and/or gestational hypertension, 12 of which are new (for example, MTHFR-CLCN6, WNT3A, NPR3, PGR and RGL3), including two loci (PLCE1 and FURIN) identified in the multitrait analysis. Identified loci highlight the role of natriuretic peptide signaling, angiogenesis, renal glomerular function, trophoblast development and immune dysregulation. We derived genome-wide polygenic risk scores that predicted preeclampsia/eclampsia and gestational hypertension in external cohorts, independent of clinical risk factors, and reclassified eligibility for low-dose aspirin to prevent preeclampsia. Collectively, these findings provide mechanistic insights into the hypertensive disorders of pregnancy and have the potential to advance pregnancy risk stratification.

Noncanonical role for Ku70/80 in the prevention of allergic airway inflammation via maintenance of airway epithelial cell organelle homeostasis.

Airway epithelial homeostasis is under constant threat due to continuous exposure to the external environment, and abnormally robust sensitivity to external stimuli is critical to the development of airway diseases, including asthma. Ku is a key nonhomologous end-joining DNA repair protein with diverse cellular functions such as VDJ recombination and telomere length maintenance. Here, we show a novel function of Ku in alleviating features of allergic airway inflammation via the regulation of mitochondrial and endoplasmic reticulum (ER) stress. We first determined that airway epithelial cells derived from both asthmatic lungs and murine asthma models demonstrate increased expression of 8-hydroxy-deoxyguanosine (8-OHdG), a marker of oxidative DNA damage. Ku protein expression was dramatically reduced in the bronchial epithelium of patients with asthma as well as in human bronchial epithelial cells exposed to oxidative stress. Knockdown of Ku70 or Ku80 in naïve mice elicited mitochondrial collapse or ER stress, leading to bronchial epithelial cell apoptosis and spontaneous development of asthma-like features, including airway hyperresponsiveness, airway inflammation, and subepithelial fibrosis. These findings demonstrate an essential noncanonical role for Ku proteins in asthma pathogenesis, likely via maintenance of organelle homeostasis. This novel function of Ku proteins may also be important in other disease processes associated with organelle stress.

Interactomics Analyses of Wild-Type and Mutant A1CF Reveal Diverged Functions in Regulating Cellular Lipid Metabolism.

Population genetic studies highlight a missense variant (G398S) of A1CF that is strongly associated with higher levels of blood triglycerides (TGs) and total cholesterol (TC). Functional analyses suggest that the mutation accelerates the secretion of very low-density lipoprotein (VLDL) from the liver by an unknown mechanism. Here, we used multiomics approaches to interrogate the functional difference between the WT and mutant A1CF. Using metabolomics analyses, we captured the cellular lipid metabolite changes induced by transient expression of the proteins, confirming that the mutant A1CF is able to relieve the TG accumulation induced by WT A1CF. Using a proteomics approach, we obtained the interactomic data of WT and mutant A1CF. Networking analyses show that WT A1CF interacts with three functional protein groups, RNA/mRNA processing, cytosolic translation, and, surprisingly, mitochondrial translation. The mutation diminishes these interactions, especially with the group of mitochondrial translation. Differential analyses show that the WT A1CF-interacting proteins most significantly different from the mutant are those for mitochondrial translation, whereas the most significant interacting proteins with the mutant are those for cytoskeleton and vesicle-mediated transport. RNA-seq analyses validate that the mutant, but not the WT, A1CF increases the expression of the genes responsible for cellular transport processes. On the contrary, WT A1CF affected the expression of mitochondrial matrix proteins and increased cell oxygen consumption. Thus, our studies confirm the previous hypothesis that A1CF plays broader roles in regulating gene expression. The interactions of the mutant A1CF with the vesicle-mediated transport machinery provide mechanistic insight in understanding the increased VLDL secretion in the A1CF mutation carriers.

Spectrum of clinical features and genetic variants in mevalonate kinase (MVK) gene of South Indian families suffering from Hyperimmunoglobulin D Syndrome.

Hyper-IgD syndrome (HIDS, OMIM #260920) is a rare autosomal recessive autoinflammatory disorder caused by pathogenic variants in the mevalonate kinase (MVK) gene. HIDS has an incidence of 1:50,000 to 1:5,000, and is thought to be prevalent mainly in northern Europe. Here, we report a case series of HIDS from India, which includes ten patients from six families who presented with a wide spectrum of clinical features such as recurrent fever, oral ulcers, rash, arthritis, recurrent diarrhea, hepatosplenomegaly, and high immunoglobulin levels. Using whole exome sequencing (WES) and/or Sanger capillary sequencing, we identified five distinct genetic variants in the MVK gene from nine patients belonging to six families. The variants were classified as pathogenic or likely pathogenic as per the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) guidelines for annotation of sequence variants. Over 70% of patients in the present study had two recurrent mutations in MVK gene i.e. a nonsynonymous variant p.V377I, popularly known as the 'Dutch mutation', along with a splicing variant c.226+2delT in a compound heterozygous form. Identity by descent analysis in two patients with the recurrent variants identified a 6.7 MB long haplotype suggesting a founder effect in the South Indian population. Our analysis suggests that a limited number of variants account for the majority of the patients with HIDS in South India. This has implications in clinical diagnosis, as well as in the development of cost-effective approaches for genetic diagnosis and screening. To our best knowledge, this is the first and most comprehensive case series of clinically and genetically characterized patients with HIDS from India.

Saliva microbiome in primary Sjögren's syndrome reveals distinct set of disease-associated microbes.

OBJECTIVE: This study systematically aims to evaluate the salivary microbiome in patients with primary Sjögren's syndrome (pSS) using 16S rRNA sequencing approach. METHODS: DNA isolation and 16S rRNA sequencing was performed on saliva of 37 pSS and 35 control (CC) samples on HiSeq 2500 platform. 16S rRNA sequence analysis was performed independently using two popular computational pipelines, QIIME and less operational taxonomic units scripts (LoTuS). RESULTS: There were no significant changes in the alpha diversity between saliva of patients and controls. However, four genera including Bifidobacterium, Lactobacillus, Dialister and Leptotrichia were found to be differential between the two sets, and common between both QIIME and LoTuS analysis pipelines (Fold change of 2 and p < .05). Bifidobacterium, Dialister and Lactobacillus were found to be enriched, while Leptotrichia was significantly depleted in pSS compared to the controls. Exploration of microbial diversity measures (Chao1, observed species and Shannon index) revealed a significant increase in the diversity in patients with renal tubular acidosis. An opposite trend was noted, with depletion of diversity in patients with steroids. CONCLUSION: Our analysis suggests that while no significant changes in the diversity of the salivary microbiome could be observed in Sjögren's syndrome compared to the controls, a set of four genera were significantly and consistently differential in the saliva of patients with pSS. Additionally, a difference in alpha diversity in patients with renal tubular acidosis and those on steroids was observed.

Role of Tmem163 in zinc-regulated insulin storage of MIN6 cells: Functional exploration of an Indian type 2 diabetes GWAS associated gene.

Genome wide association study for type 2 diabetes discovered TMEM163 as a risk locus. Perturbations in TMEM163 expression was reported to be associated with impaired intracellular zinc homeostasis. Physiological concentration of zinc is instrumental to maintain insulin storage and functionality in pancreatic ß cells. We found abundant TMEM163 expression in human pancreas, both at transcriptional and translational levels. Knockdown of endogenous Tmem163 in MIN6 cells resulted in increased intracellular zinc and total insulin content, coupled with compromised insulin secretion at high glucose stimuli. Furthermore, Tmem163 knockdown led to enhanced cellular glucose uptake. Upon next generation sequencing, one-third of the studied T2D patients were found to have a novel missense variant in TMEM163 gene. Study participants harboring this missense variant displayed a trend of higher glycemic indices. This is the first report on exploring the biological role of TMEM163 in relation to T2D pathophysiology.

A genome-wide map of circular RNAs in adult zebrafish.

Circular RNAs (circRNAs) are transcript isoforms generated by back-splicing of exons and circularisation of the transcript. Recent genome-wide maps created for circular RNAs in humans and other model organisms have motivated us to explore the repertoire of circular RNAs in zebrafish, a popular model organism. We generated RNA-seq data for five major zebrafish tissues - Blood, Brain, Heart, Gills and Muscle. The repertoire RNA sequence reads left over after reference mapping to linear transcripts were used to identify unique back-spliced exons utilizing a split-mapping algorithm. Our analysis revealed 3,428 novel circRNAs in zebrafish. Further in-depth analysis suggested that majority of the circRNAs were derived from previously well-annotated protein-coding and long noncoding RNA gene loci. In addition, many of the circular RNAs showed extensive tissue specificity. We independently validated a subset of circRNAs using polymerase chain reaction (PCR) and divergent set of primers. Expression analysis using quantitative real time PCR recapitulate selected tissue specificity in the candidates studied. This study provides a comprehensive genome-wide map of circular RNAs in zebrafish tissues.

Organellar transcriptome sequencing reveals mitochondrial localization of nuclear encoded transcripts.

Mitochondria are organelles involved in a variety of biological functions in the cell, apart from their principal role in generation of ATP, the cellular currency of energy. The mitochondria, in spite of being compact organelles, are capable of performing complex biological functions largely because of the ability to exchange proteins, RNA, chemical metabolites and other biomolecules between cellular compartments. A close network of biomolecular interactions are known to modulate the crosstalk between the mitochondria and the nuclear genome. Apart from the small repertoire of genes encoded by the mitochondrial genome, it is now known that the functionality of the organelle is highly reliant on a number of proteins encoded by the nuclear genome, which localize to the mitochondria. With exceptions to a few anecdotal examples, the transcripts that have the potential to localize to the mitochondria have been poorly studied. We used a deep sequencing approach to identify transcripts encoded by the nuclear genome which localize to the mitoplast in a zebrafish model. We prioritized 292 candidate transcripts of nuclear origin that are potentially localized to the mitochondrial matrix. We experimentally demonstrated that the transcript encoding the nuclear encoded ribosomal protein 11 (Rpl11) localizes to the mitochondria. This study represents a comprehensive analysis of the mitochondrial localization of nuclear encoded transcripts. Our analysis has provided insights into a new layer of biomolecular pathways modulating mitochondrial-nuclear cross-talk. This provides a starting point towards understanding the role of nuclear encoded transcripts that localize to mitochondria and their influence on mitochondrial function.

Identification of novel circadian transcripts in the zebrafish retina.

High fecundity, transparent embryos for monitoring the rapid development of organs and the availability of a well-annotated genome has made zebrafish a model organism of choice for developmental biology and neurobiology. This vertebrate model, which is also a favourite in chronobiology studies, shows striking circadian rhythmicity in behaviour. Here, we identify novel genes in the zebrafish genome that are expressed in the zebrafish retina. We further resolve the expression pattern over time and tentatively assign specific novel transcripts to retinal bipolar cells of the inner nuclear layer. Using chemical ablation and free run experiments, we segregate the transcripts that are rhythmic when entrained by light from those that show sustained oscillations in the absence of external cues. The transcripts reported here with rigorous annotation and specific functions in circadian biology provide the groundwork for functional characterization of novel players in the zebrafish retinal clock.

SAGE: a comprehensive resource of genetic variants integrating South Asian whole genomes and exomes.

South Asia is home to $\sim $20% of the world population and characterized by distinct ethnic, linguistic, cultural and genetic lineages. Only limited representative samples from the region have found its place in large population-scale international genome projects. The recent availability of genome scale data from multiple populations and datasets from South Asian countries in public domain motivated us to integrate the data into a comprehensive resource. In the present study, we have integrated a total of six datasets encompassing 1213 human exomes and genomes to create a compendium of 154 814 557 genetic variants and adding a total of 69 059 255 novel variants. The variants were systematically annotated using public resources and along with the allele frequencies are available as a browsable-online resource South Asian genomes and exomes. As a proof of principle application of the data and resource for genetic epidemiology, we have analyzed the pathogenic genetic variants causing retinitis pigmentosa. Our analysis reveals the genetic landscape of the disease and suggests subset of genetic variants to be highly prevalent in South Asia.

RNA secondary structure profiling in zebrafish reveals unique regulatory features.

BACKGROUND: RNA is known to play diverse roles in gene regulation. The clues for this regulatory function of RNA are embedded in its ability to fold into intricate secondary and tertiary structure. RESULTS: We report the transcriptome-wide RNA secondary structure in zebrafish at single nucleotide resolution using Parallel Analysis of RNA Structure (PARS). This study provides the secondary structure map of zebrafish coding and non-coding RNAs. The single nucleotide pairing probabilities of 54,083 distinct transcripts in the zebrafish genome were documented. We identified RNA secondary structural features embedded in functional units of zebrafish mRNAs. Translation start and stop sites were demarcated by weak structural signals. The coding regions were characterized by the three-nucleotide periodicity of secondary structure and display a codon base specific structural constrain. The splice sites of transcripts were also delineated by distinct signature signals. Relatively higher structural signals were observed at 3' Untranslated Regions (UTRs) compared to Coding DNA Sequence (CDS) and 5' UTRs. The 3' ends of transcripts were also marked by unique structure signals. Secondary structural signals in long non-coding RNAs were also explored to better understand their molecular function. CONCLUSIONS: Our study presents the first PARS-enabled transcriptome-wide secondary structure map of zebrafish, which documents pairing probability of RNA at single nucleotide precision. Our findings open avenues for exploring structural features in zebrafish RNAs and their influence on gene expression.

Egyptian tale from India: application of whole-exome sequencing in diagnosis of atypical familial Mediterranean fever.

Clinical diagnosis of autoinflammatory diseases requires a high degree of clinical suspicion and clinching molecular evidence to substantiate the diagnosis. This is more so in populations with low prevalence of these disorders. In this report, we describe the case of a young man from India with recurrent fever and persistent arthritis. The patient's forefathers were of Egyptian ancestry who practiced consanguinity. Molecular genetic analysis using whole-exome sequencing suggested the presence of variants c.443A>T:p.E148V and c.442G>C:p.E148Q in the MEFV gene, earlier independently shown to be associated with familial Mediterranean fever (FMF) in a compound heterozygous state. The variants were further confirmed by capillary sequencing. This report also highlights the application of whole exome sequencing to delineate the allelic differences in the variants apart from serving as a quick genetic screening approach for autoinflammatory diseases. To the best of our knowledge, this is the first report of a compound heterozygosity for the two well-characterized variants associated with atypical FMF in a patient.

Aptamer-Assisted Detection of the Altered Expression of Estrogen Receptor Alpha in Human Breast Cancer.

An increase in the expression of estrogen receptors (ER) and the expanded population of ER-positive cells are two common phenotypes of breast cancer. Detection of the aberrantly expressed ERα in breast cancer is carried out using ERα-antibodies and radiolabelled ligands to make decisions about cancer treatment and targeted therapy. Capitalizing on the beneficial advantages of aptamer over the conventional antibody or radiolabelled ligand, we have identified a DNA aptamer that selectively binds and facilitates the detection of ERα in human breast cancer tissue sections. The aptamer is identified using the high throughput sequencing assisted SELEX screening. Biophysical characterization confirms the binding and formation of a thermodynamically stable complex between the identified DNA aptamer (ERaptD4) and ERα (Ka = 1.55±0.298×108 M(-1); ΔH = 4.32×104±801.1 cal/mol; ΔS = -108 cal/mol/deg). Interestingly, the specificity measurements suggest that the ERaptD4 internalizes into ERα-positive breast cancer cells in a target-selective manner and localizes specifically in the nuclear region. To harness these characteristics of ERaptD4 for detection of ERα expression in breast cancer samples, we performed the aptamer-assisted histochemical analysis of ERα in tissue samples from breast cancer patients. The results were validated by performing the immunohistochemistry on same samples with an ERα-antibody. We found that the two methods agree strongly in assay output (kappa value = 0.930, p-value <0.05 for strong ERα positive and the ERα negative samples; kappa value = 0.823, p-value <0.05 for the weak/moderate ER+ve samples, n = 20). Further, the aptamer stain the ERα-positive cells in breast tissues without cross-reacting to ERα-deficient fibroblasts, adipocytes, or the inflammatory cells. Our results demonstrate a significant consistency in the aptamer-assisted detection of ERα in strong ERα positive, moderate ERα positive and ERα negative breast cancer tissues. We anticipate that the ERaptD4 aptamer targeting ERα may potentially be used for an efficient grading of ERα expression in cancer tissues.

Comparative whole-genome analysis of clinical isolates reveals characteristic architecture of Mycobacterium tuberculosis pangenome.

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.

mit-o-matic: a comprehensive computational pipeline for clinical evaluation of mitochondrial variations from next-generation sequencing datasets.

The human mitochondrial genome has been reported to have a very high mutation rate as compared with the nuclear genome. A large number of mitochondrial mutations show significant phenotypic association and are involved in a broad spectrum of diseases. In recent years, there has been a remarkable progress in the understanding of mitochondrial genetics. The availability of next-generation sequencing (NGS) technologies have not only reduced sequencing cost by orders of magnitude but has also provided us good quality mitochondrial genome sequences with high coverage, thereby enabling decoding of a number of human mitochondrial diseases. In this study, we report a computational and experimental pipeline to decipher the human mitochondrial DNA variations and examine them for their clinical correlation. As a proof of principle, we also present a clinical study of a patient with Leigh disease and confirmed maternal inheritance of the causative allele. The pipeline is made available as a user-friendly online tool to annotate variants and find haplogroup, disease association, and heteroplasmic sites. The "mit-o-matic" computational pipeline represents a comprehensive cloud-based tool for clinical evaluation of mitochondrial genomic variations from NGS datasets. The tool is freely available at http://genome.igib.res.in/mitomatic/.

Draft Genome Sequence of the Extremely Halophilic Bacterium Halomonas salina Strain CIFRI1, Isolated from the East Coast of India.

Halomonas salina strain CIFRI1 is an extremely salt-stress-tolerant bacterium isolated from the salt crystals of the east coast of India. Here we report the annotated 3.45-Mb draft genome sequence of strain CIFRI1 having 86 contigs with 3,139 protein coding loci, including 62 RNA genes.