Genomics

Pseudomonas aeruginosa, is a gram-negative bacterium causes opportunistic or nosocomial infections in immunocompromised individuals. In recent years, a steady increase in human corneal infections of P. aeruginosa has been reported with increased multi-drug resistance (MDR) or extensively drug resistance (XDR). Several non-coding sRNAs, has been identified to regulate various physiological processes in P. aeruginosa, including biofilm formation, quorum sensing. However, the regulatory mechanism of sRNAs in MDR/XDR pathways of P. aeruginosa keratitis strains is not yet studied. In this study, we identified bacterial sRNAs in publicly available P. aeruginosa keratitis genomes and investigated their regulatory role in MDR/XDR pathways using bioinformatic analysis. Totally, 46 P. aeruginosa keratitis strains from different geographical regions were included. Of 46, Eight (30%) out of Twenty-seven and Nine (52) out of Nineteen P. aeruginosa strains from India and Australia were identified as not-MDR. Whereas, 10 (38%) Indian and 9 (47%) Australian strains were identified as MDR. Eight Indian strains were identified as XDR. Out of 46 strains, 23 (50%) carried ExoU, 21(45%) carried ExoS and two (5%) strains carried both ExoU and ExoS, exotoxins for their virulence. The sRNA, SPA0021 was identified in 18 MDR/XDR and 6 not-MDR strains along with UCBPP-PA14. Interestingly, majority of the imipenem resistant P. aeruginosa keratitis strains from the present study was found to be carried SPA0023 sRNA (18 out of 30 strains). The outer membrane porin protein OprD, identified as binding target of SPA0023. Negative regulation or inactivation of OprD, reported in increased imipenem resistance in P. aeruginosa. Mutation analysis revealed that SPA0023 carrying P. aeruginosa keratitis strains contains a lesser number of amino acid changes in OprD protein than other strains. These findings indicate, imipenem resistance in SPA0023 carried strains might arose from the negative regulation or inhibition of OprD by SPA0023. However, functional studies are warranted with large number of P. aeruginosa keratitis strains to confirm the negative regulation of OprD by SPA0023 and imipenem resistance.

Isoniazid (INH) is an important first-line medication for the treatment of tuberculosis. The impact that tuberculosis drug resistance has on treatment outcomes is a topic that is receiving a lot of attention these days because of the rising incidence of INH-resistant cases. Study involves a single group of patients who have been diagnosed with Isoniazid monoresistant tuberculosis. Treatment history and demographic data of the patients were obtained after informed consent. The mutation patterns of isoniazid were observed after multiplex PCR and Line Probe Assay (LPA). A total of 101 patient (M,F) records at the IRL, Puducherry were analyzed. The predominant gene responsible for TB was KATG (67.3%). The KATG Mut1 was a prime mutation observed in the present study population (58.41%). Study showed positive association with males (74%), occupation as coolie (88%), diabetes as comorbidity (33%), pulmonary tuberculosis as the TB site (98.01%), history of previous ATT intake in 43 patients (42.6%), katG mutation (67.3%), katG Mut 1 was the prime mutation (58.4%).The cure rate was high in INH high concentration resistance patients which was statistically significant (p=0.0167). INH monoresistance mutations seen in 64.3% of the patients with katG, compared to inhA (34.65%). Similar to katG mutations, inhA mutations also have MUT1 as their most frequent gene pattern. There is a significant association between males, diabetes, smoking and alcohol addictions were associated with high risk of developing high dose INH monoresistance (katG). High prevalence of recurrent tuberculosis was seen in high dose INH monoresistance tuberculosis. Patients who are microbiologically confirmed pulmonary tuberculosis and diabetes with rifampicin sensitive status needs to be checked for LPA for isoniazid sensitivity status to prevent treatment failure and relapse. It is crucial to understand the gene pattern in each of these patients since these mutations are closely associated to high or low-degree resistance to INH

Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families are known as RNA silencing machinery genes or RNAi genes. They have important activities at post-transcriptional and chromatin modification levels. They regulate gene expression relating to different stresses, growth, and development in eukaryotes. A complete cycle of gene silencing is occurred by the collaboration of these three families. However, these gene families are not yet rigorously studied in the economically important wheat genome. Our bioinformatic analysis based genome-wide identification, characterization, diversification and regulatory components of these gene families identified 7 TaDCL, 39 TaAGO and 16 TaRDR genes from wheat genome against RNAi genes of Arabidopsis thaliana. Phylogenetic analysis of wheat genome with Arabidopsis and rice RNAi genes showed that TaDCL, TaAGO and TaRDR proteins are clustered into four, eight and four subgroups respectively. Domain, motif and exon-intron structure analyses showed that the TaDCL, TaAGO and TaRDR proteins conserve identical characteristics within groups while retain diverse differences between groups. GO annotations implied that a number of biological and molecular pathways are linked to RNAi mechanism in wheat. Gene networking between transcription factors and RNAi proteins indicates that ERF is the leading family linked to maximum RNAi genes followed by MIKC-MADS, C2H2, BBR-BPC, MYB, and Dof. Cis-regulatory elements associated to RNAi genes are predicted to act as regulatory components against various environmental conditions. Expressed sequence tag analysis showed that larger numbers of RNAi genes are expressed in different tissues and organs predicted to play roles for healthy plants and grains. Expression analysis of 7 TaDCL genes using qRT-PCR showed that only TaDCL3a and TaDCL3b had root specific significant expression (p-value<0.05) with no expression in leaf validated EST results. Besides, TaDCL3b and TaDCL4 significantly prompted in drought condition indicating their potential role in drought stress tolerance. Overall results would however help researchers for in-depth biological investigation of these RNAi genes in wheat crop improvement.

Heterochromatin protein 1 (HP1) is widespread in several organisms playing a role in control of gene expression by heterochromatin formation and maintenance of silent chromatin. Schistosoma mansoni is a human parasite that is responsible for Schistosomiasis, a tropical neglected disease in the tropical and subtropical areas in the world, where the intermediate host Biomphalaria glabrata is present. In this study we attempted to investigate if the SmHP1 is enriched in S. mansoni chromatin in cercariae larvae stage, compared with another larvae stage sporocysts and its importance for S. mansoni life cycle progression and parasite fitness. We used ChIPmentation with commercial antibody ab109028 that passed in-house quality control. Our data show that S. mansoni HP1 enrichment is non-canonical with a peak at the transcription end sites of protein coding genes. We did not find strong differences in SmHP1 chromatin landscapes between sporocysts and cercariae. Knock-down of SmHP1 in schistosomula and in vivo experiments in mice unexpectedly increased parasite fitness. Our results suggest that SmHP1 may influence chromatin structure in a non-canonical way in S. mansoni stages and may play a role in regulation of parasite fitness.

SARS-CoV-2 evolution has continued to generate variants, responsible for new pandemic waves locally and globally. Varying disease presentation and severity has been ascribed to inherent variant characteristics and vaccine immunity. This study analyzed genomic data from 305 whole genome sequences from SARS-CoV-2 patients before and through the third wave in India. Delta variant was responsible for disease in patients without comorbidity(97%), while Omicron BA.2 caused disease primarily in those with comorbidity(77%). Tissue adaptation studies brought forth higher propensity of Omicron variants to bronchial tissue than lung, contrary to observation in Delta variants from Delhi. Study of codon usage pattern distinguished the prevalent variants, clustering them separately, Omicron BA.2 isolated in February grouped away from December strains, and all BA.2 after December acquired a new mutation S959P in ORF1b (44.3% of BA.2 in the study) indicating ongoing evolution. Loss of critical spike mutations in Omicron BA.2 and gain of immune evasion mutations including G142D, reported in Delta but absent in BA.1, and S371F instead of S371L in BA.1 could possibly be due to evolutionary trade-off and explain very brief period of BA.1 in December 2021, followed by complete replacement by BA.2.

Pseudomonas aeruginosa has been considered one of the major nosocomial pathogens associated with elevated morbidity and mortality worldwide. Outbreaks have been associated with few high-risk pandemic P. aeruginosa lineages, presenting a remarkable antimicrobial resistance. However, the biological features involved with the persistence and spread of such lineages among clinical settings remain to be unravel. This study reports the emergence of the ST309 P. aeruginosa lineage in South America/Brazil, more precisely, in the Amazon region. Global genomic analyses were performed with the Brazilian strain (PA834) and more 41 complete and draft ST309 genomes publicly available, giving insights about ST309 epidemiology and its resistome and mobilome. Antimicrobial susceptibility tests revealed that the Brazilian PA834 strain presented the XDR phenotype, which was mainly due to intrinsic resistance mechanisms. Genomic analyses revealed a heterogeneous distribution of acquired antimicrobial resistance genes among ST309 genomes, which included blaVIM-2, blaIMP-15 and qnrVC1, all of them associated with class 1 integrons. The mobilome mining showed the presence of Integrative and Conjugative Elements, transposons and genomic islands harbouring a huge arsenal of hevy metal resistance genes. Moreover, these elements also carried genes involved with virulence and adaptive traits. Therefore, the presence of such genes in ST309 lineage possibly accounted for the global spread and persistence of this emerging clone, and for its establishment as a pandemic lineage of clinical importance.

The limited knowledge of genomic non-coding and regulatory regions has limited our ability to decipher the genetic mechanisms underlying complex traits in pigs. In this study, we characterize the spatiotemporal landscape of putative enhancers and promoters and their target genes by combining H3K27ac targeted ChIP-Seq and RNA-Seq in fetal (day 74-75 pc) and adult (day 132-150 pn) tissues (brain, liver, heart, muscle and small intestine) sampled from Asian aboriginal Bamaxiang and European highly selected Large White pigs of both sexes. We identify 101,290 H3K27ac peaks marking 18,521 promoters and 82,769 enhancers, including peaks that are active across all tissues and developmental stages could indicate safe harbors for exogenous gene insertion, and tissue and developmental-stage specific peaks that regulate genes pathways matching tissue and developmental stage specific physiological functions. We found H3K27ac and DNA methylation in the promoter region of the XIST gene may involve in X chromosome inactivation, and demonstrate utility of the present resource to reveal regulatory patterns of known causal genes and to prioritize candidate causal variants for complex traits in pigs. We have developed a web browser to improve the accessibility of the results (http://39.108.231.116/browser/?genome=susScr11).

Idiopathic Pulmonary Fibrosis (IPF) is a chronic and progressive lung disease that leads to gradual decline in lung function. The molecular mechanism and risk factors of this disease are still obscure. Poorly understood etiology of this disease is the major obstacle in the identification of potential biomarkers and drug targets. In this study, microarray gene expression data of normal and IPF patient has been utilized for the statistical analysis of differentially expressed genes (DEGs) with a view to identifying potential molecular signatures using network-based system. Then their functional enrichment analysis revealed their predominant involvement in transcription, protein acetylation, extracellular matrix organization, apoptic process, inflammatory response etc. Protein-Protein Interaction (PPI) network revealed (UBC, PTEN, SOS1, PTK2, FGFR1, YAP1, FOXO1, RACK1, BMP4 and CD44) as hub proteins in IPF. Subsequent regulatory network analysis suggested (E2F1, STAT3, PPARG, MEF2A, FOXC1, GATA3, YY1, GATA2, NFKB1, and FOXL1) as the best regulatory transcriptional signatures and (hsa-mir-155-5p, hsa-mir-16-5p, hsa-mir-17-5p, hsa-mir-19a-3p, hsa-mir-192-5p, hsa-mir-92a-3p, hsa-mir-26b-5p, hsa-mir-335-5p, hsa-mir-124-3p, and hsa-let-7b-5p) as the best post-transcriptional signatures. This study represents proteome and RNA signatures of IPF which might be useful to uphold the present efforts in the discovery of potential biomarkers and treatments of this disease.

Background and AimsThree out of four RNA components of ribosomes are encoded by 45S rDNA loci, whose transcripts are processed into 18S, 5.8S and 26S ribosomal RNAs. The loci are organized as long head-to-tail tandem arrays of nearly identical units spanning over several megabases of sequence. Due to this peculiar structure, the number of rRNA genes, their sequence composition and expression status remain unclear, especially in complex polyploid genomes harbouring multiple loci. Here we conducted a complex study to decipher structure and activity of both major and minor rRNA loci in hexaploid bread wheat (Triticum aestivum). MethodsWe employed an original, multi-omics approach, combining chromosome flow sorting and optical mapping with transcriptome and methylome sequencing. Key ResultsThe former two techniques enabled unbiased quantification of rDNA units in particular loci of the wheat genome. Total number of rRNA genes organized in tandem arrays was 4388, with 64.1, 31.4, 3.9 and 0.7% located in short arms of chromosomes 6B, 1B, 5D and 1A, respectively. At the expression level, only 1B and 6B loci contributed to transcription at roughly 2:1 ratio. The 1B:6B ratio varied among five analysed tissues (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG contexts in the silenced 5D locus compared to the active 1B and 6B loci. ConclusionsA fine genomic organization and tissue-specific expression of rRNA loci were deciphered, for the first time, in a complex polyploid species. We documented various mechanisms of rRNA dosage control, including gene elimination and stable inactivation related to nucleolar subdominance of A and D-genome loci, and a subtle, developmentally regulated silencing of one of the major loci. The results are discussed in the context of wheat evolution and transcription regulation.

Transposable elements function as one of the major effectors in response to biological or environmental stress. Under normal conditions, host organisms deploy epigenetic and post-transcriptional machinery (histone modifications, chromatin remodelers, long non-coding RNAs (lncRNAs)) at the TE sites to contain their mobility. But many a times, the chromatin architecture undergoes TE induced changes under the effect of stress that in turn might lead to unprecedented gene expression. LncRNAs are emerging as a crucial tool in the regulation of TEs. TEs possess remarkable abilities to respond in the face of stress, ranging from undetected mutations to changing the regulatory landscape of the host. Although the relationship between stress response and TE activation/deactivation is well acknowledged but our understanding of the mechanism of regulation remains poor. This study focuses on the gene expression of THAP9, a domesticated transposon and lncRNA THAP9-AS1 (THAP9-anti sense1), which form a sense and anti-sense gene pair with a promoter overlap of approximately 350bp. The two genes exhibit different patterns of gene expression under different types of stresses in the S-phase of the cell cycle. THAP9-AS1 is always upregulated under stress whereas THAP9 exhibits both downregulation and upregulation in different stresses. Both THAP9 and THAP9-AS1 exhibit a periodic gene expression throughout the S-phase which is a characteristic of cell cycle regulated genes.

This study reports the novel allelic diversity in HKT 1;5 gene (High Affinity Potassium Transporter) in bread wheat and its phylogenetic relationship among the paralogs/orthologs of in Triticum aestivum and its wild relatives. HKT 1;5 gene is a known and pivotal gene associated with the salinity tolerance in plants upon the discrimination of K+ over Na+ in leaves without change in Na+ concentration in root. This gene was sequenced in a diverse collection of bread wheat, durum wheat, wild relatives, and ditelosomic lines. Sequence analysis in bread wheat led to the identification of four alleles, which could be distinguished by number of SNPs. Sequence comparison between monocot bread wheat and dicot Arabidopsis thaliana revealed that the HKT1;5 gene is conserved at level of exonic regions; however, the presence of transposable elements especially in intronic regions is further intriguing towards evolutionary relatedness. Two paralogous or major alleles observed in Triticum monoccum and Aegilops tauschii were further categorized as sub-alleles based on their SNPs comparison. This gene was absent in T. urartu in accordance with existing evidence, while it was found in A. speltoides (an allelic variant) with a few base pairs insertion in the exon1 region causing a frameshift mutation with an altered amino acids and genomic database mining unveiled additional alleles in this species. Ditelosomic lines with 4DL and 4DS chromosomes revealed a higher similarity with bread and durum wheat respectively. Phylogenetic studies of HKT1;5 orthologs from different Poaceae species revealed the occurrence of five different ortholog groups with taxonomic consistency. Phenotyping based salinity stress experiment distinguished the unknown lines for salinity tolerance and sensitiveness in comparison with known reference lines and possible allelic comparison was made. The salinity stress analysis further revealed that some known drought/heat tolerance lines showed slightly better salinity tolerance with mean values and variability of traits than known saline tolerant wheat line at controlled ambient.

BackgroundSpatial organization of the genome is fundamental for ensuring accurate gene expression. This process depends on the communication between gene promoters and distal cis-regulatory elements (CREs), which together make up 8% of the human genome and are supported by the chromatin structure. It is estimated that over 90% of disease-associated variants are located in the non-coding region of the genome and may affect CRE. For the cystic fibrosis transmembrane conductance regulator (CFTR) gene, a complete understanding of tissue-specific CFTR expression and regulation is missing, in particular in the pancreas. Mechanistic insights into tissue-specific expression may provide clarity on the clinical heterogeneity observed in Cystic Fibrosis and CFTR-related disorders. MethodsTo understand the role of 3D chromatin architecture in establishing tissue-specific expression of the CFTR gene, we mapped chromatin interactions via circular chromosome conformation capture (4C) and epigenomic regulation through H3K27ac and DNase Hypersensitive site I (DHS) in Capan-1 pancreatic cells. Candidate regulatory regions are validated by luciferase reporter assay and CRISPR-knock out. ResultsWe identified active regulatory regions not only around the CFTR gene but also outside the topologically associating domain (TAD). By performing functional assays, we validated our targets and revealed a cooperative effect of the -44 kb, -35 kb, +15.6 kb and 37.7 kb regions, which share common predicted transcription factor (TF) motifs. Comparative 3D genomic analysis and functional assays using the Caco-2 intestinal cell line revealed the presence of tissue-specific CREs. ConclusionBy studying the chromatin architecture of the CFTR locus in Capan-1 cells, we demonstrated the involvement of multiple CREs upstream and downstream of the CFTR gene. We also extend our analysis to compare intestinal and pancreatic cells and provide information on the tissue-specificity of CRE. These findings highlight the importance of expanding the search for causative variants beyond the gene coding sequence but also by considering the tissue-specific 3D genome.

Entamoeba is the amoebozoan parasite commonly found in the intestines of animals. E. marina is the first exception isolated from marine sediments, possibly adapting from animal intestines to the sea. However, the evolutionary process of E. marina remains uncertain due to the lack of a genome sequence. Here, we present the de novo genome and transcriptome of E. marina using Oxford Nanopore MinION and Illumina HiSeq/MiSeq. The genome of E. marina is approximately 37.5 Mbp in length and consisted of 202 contigs, which is the second longest followed by E. invadens. E. marina showed significant reduction in the major virulence-associated gene families, including cysteine proteases, lysosomal enzyme transporters, and surface galactose/N-acetylglucosamine-specific lectins, suggesting diversification, more specifically reduction of pathogenicity-related genes. Genome and RNA-seq analyses also indicated genes either conserved throughout eukaryotes or laterally transferred from prokaryotes, and potentially responsible for salt tolerance. Our study provides insights into the mechanism underlying the lifestyle changes in the evolution of parasitic eukaryotes.

The last decade has witnessed great progresses regarding the molecular basis of breast cancer with discovery of different nuclear susceptibility genes; in addition investigations and researches regarding mitochondrial DNA (mtDNA) mutations in breast cancer have been started. Mitochondrial haplogroup determinants (single nucleotide polymorphism SNP) and somatic mitochondrial mutations have recently been studied as possible risk factors for carcinogenic processes in different tissues, hence in order to identify breast cancer related SNPs and haplogroups among the population of Sulaimaniyah city/Iraq, the entire mitochondrial genome of 20-breast cancer samples and comparable controls were sequenced. Haplogrep 2.0 was used for haplogroup identification; Chi-square and Fishers exact test were applied to assess relational significance. HV haplogroup in the cancer samples appeared to be a risk factor for breast cancer compared to the most common H haplogroup in control samples with a p-values of 0.002 and 0.006 respectively and an Odd Ratio (OR) = 28.00. Besides, SNP (A8860G) was also identified as a risk factor for breast cancer as compared to other randomly selected SNPs (A750G, A1438G and C7028T) with p values {square}0.05 and OR >1.

Hypervirulent K. pneumoniae (HvKp) is typically associated with ST23 clone; however, hvKp is also emerging from clones ST11, ST15 and ST147, which are also multi-drug resistant (MDR). Here, we aimed to characterise nine novel MDR hvKp isolates harbouring mosaic plasmids simultaneously carrying antimicrobial resistance (AMR) and virulence genes. Nine HvKp isolates obtained from hospitalised patients in southern India were characterized for antimicrobial susceptibility and hypervirulence phenotypes. All nine hvKp isolates were subjected to whole genome sequencing (WGS) using Ilumina HiSeq2500 and a subset of four were sequenced using Oxford Nanopore MinION. Among the nine isolates, seven were carbapenem-resistant, two of which carried blaNDM-5 on an IncFII plasmid and five carried blaOXA-232 on a ColKP3 plasmid. The virulence determinants were encoded in a mosaic plasmid ([~]320 Kbp) that formed as a result of its insertion in a IncFIB-IncHI1B plasmid co-integrate. The mosaic plasmid carried AMR genes (aadA2, armA, blaOXA-1, msrE, mphE, sul1 and dfrA14) in addition to rmpA2, iutA and iucABCD virulence genes. Interestingly the mosaic plasmid carried its own type IV-A3 CRISPR-cas system that is likely able to target the acquisition of IncF plasmid with the help of a traL spacer. The convergence of virulence and AMR is the biggest threat among invasive K. pneumoniae infections. However, increasing reports of the presence of mosaic plasmid carrying both AMR and virulence genes suggests MDR-hvKp isolates are no longer confined to selected clones and the containment of such isolates is very challenging. IMPORTANCEKlebsiella pneumoniae is an opportunistic pathogen that commonly associated with hospital-acquired infections in the urinary tract, respiratory tract, lung, wound sites. The organism has gained notoriety by acquiring additional genetic traits to become either hypervirulent (HV) phenotype or multidrug resistant (MDR) phenotype. Though the infections by both these phenotypes were very challenging to treat, the MDR K. pneumonia (MDR-Kp) were remained in the hospital settings while HV K. pneumonia (hvKp) strains were mostly originated from the community settings. In a recent turn of events, the evolution of MDR-Kp and hvKp has converged as both clones found to carry both MDR plasmids and virulence plasmid. These convergent strains are challenging to treat and is associated with higher mortality rate. As the recent hvKp isolates harbour mosaic plasmid encoding both AMR and virulence determinants there is a need to investigate the evolution of these pathogens. The significance of our research is in characterising the novel mosaic plasmid identified in MDR-hvKp isolates that belong sequence type (ST) 2096. Tracking the possible evolution pathway of MDR-hvKPs would greatly help in the proper surveillance and management of this superbugs. RepositoriesThe whole genome sequences of the present study isolates have been deposited in GenBank, NCBI, with accession numbers CP053765 - CP053770, CP053771 - CP053780, CP058798-CP058806, JAARNO010000001.1 - JAARNO010000005.1, JAAQSG000000000, JAARNJ000000000, JAARMH000000000 and JAAQTC000000000

Genomic imprinting, a mechanism resulting in parent-of-origin expression of genes through epigenetic regulation, intersects with a broad range of biological fields including evolution, molecular genetics and epigenetics and determinism of complex traits. Although next generation sequencing technologies enable nowadays to detect imprinted genes in a genome-wide manner, a wide spectrum of this phenomena is evaluated only in humans and mice. Here, we propose to map genes showing a parental expression bias in hypothalamus, muscle and placenta in piglets around birth using an extensive and unbiased strategy that relied on reciprocal crosses, genetics reconstruction of parental phases after imputation and statistical analyses discriminating parent-of-origin from allele-specific expression. We detected 440 unique genes with a weak to exclusive parental expression bias including 114 unique genes with an imbalance ratio above 25:75. About thirty imprinted genes are common to human and/or mice and an equivalent number is shared between tissues, suggesting an overall weak conservation landscape of genomic imprinting. Interestingly, we identified novel parent-of-origin expressed genes involved in neurodevelopmental (PITRM1, Pitrilysin Metallopeptidase 1) and fetal growth (FAM20B, Glycosaminoglycan Xylosylkinase and POU6F2, POU Class 6 Homeobox 2) functions. In addition, deeper analyses of specific loci likely highlighted lineage-specific imprinted genes such as a Zinc Finger Protein 300-like gene as well as specific imprinted isoforms of COPG2 (COPI Coat Complex Subunit Gamma 2), a gene showing conflicting data in the literature. Altogether, our results bring pig as the most comprehensively and exhaustively documented species for genomic imprinting after human and mice organisms. A weak conservation of this mechanism across species and tissues suggested a distinction between a small number of core imprinted genes and others parent-of-origin expressed genes that seemed subjected to evolutionary forces for acquiring their imprinting status either in a lineage-specific or tissue-specific manner.

Pneumonia remains one of the leading causes of death worldwide, particularly amongst the elderly and young children. We performed a genome-wide meta-analysis of lifetime pneumonia diagnosis (N=266,277), that encompassed the largest collection of cases published to date. Genome-wide significant associations with pneumonia were uncovered for the first time beyond the major histocompatibility complex region, with three novel loci, including a signal fine-mapped to a cluster of mucin genes. Moreover, we demonstrated evidence of a polygenic effect of common and low frequency pneumonia associated variation impacting several other mucin genes and O-glycosylation, further suggesting a role for these processes in pneumonia pathophysiology. The pneumonia GWAS was then leveraged to identify drug repurposing opportunities, including evidence that supports the use of lipid modifying agents in the prevention and treatment of the disorder. We also propose how polygenic risk could be utilised for precision drug repurposing through pneumonia risk scores constructed using variants mapped to pathways with known drug targets. In summary, we provide novel insights into the genetic architecture of pneumonia susceptibility, with future study warranted to functionally interrogate novel association signals and evaluate the suitability of the compounds prioritised by this study as repositioning candidates.

Aquaporins (AQPs) facilitate the transport of water or other small solutes into cells in the presence of osmotic gradients. However, the current understanding of piscine AQP gene with cellular stress responses has been still limitedly exemplified. In present study, we characterized the mud loach AQP1b gene at the nucleotide and amino acid levels. We identified three AQP 1b transcript variants (mmAQP1b_tv1, mmAQP1b_tv2, and mmAQP1b_tv3). Then, we examined the AQP1b promoter region and observed several transcription factor binding sites (TFBS) for nuclear factor of activated T-cells (NFAT), SRY-box, c-AMP responsive element binding protein (CREB), GATA binding factor, and hepatic nuclear factor-1. Interestingly, mmAQP1b transcription was differentially modulated by heavy metal or immunostimulant challenge. Further studies to deepen the knowledge of fish AQP-mediated adaptation response potentially relevant to molecular pathogenesis are warranted. Summary statementWe identified mud loach AQP1b transcript variants and consensus sequences involved in stress or innate immunity in promotor region. AQP1b transcription was differentially modulated by heavy metal or immunostimulant challenge.

Recent studies have shown that three-dimensional architecture of bacterial chromatin plays an important role in gene expression regulation. However, genome topological organization in Mycobacterium tuberculosis (M. tuberculosis), the etiologic agent of tuberculosis, remains unknown. On the other hand, exact mechanism of differential pathogenesis in the canonical strains of M. tuberculosis H37Rv and H37Ra remains poorly understood in terms of their raw sequences. In this context, a detailed contact map from a Hi-C experiment is a candidate for what bridges the gap. Here we present the first comprehensive report on genome-wide contact maps between regions of H37Rv and H37Ra genomes. We tracked differences between the genome architectures of H37Rv and H37Ra, which could possibly explain the virulence attenuation in H37Ra. We confirm the existence of a differential organization between the two strains most significantly a higher Chromosome Interacting Domain (CID) size in attenuated H37Ra strain. CID boundaries are also found enriched with highly expressed genes and with higher operon density in H37Rv. Furthermore, most of the differentially expressed PE/PPE genes were present near the CID boundaries in H37Rv and not in H37Ra. Collectively our study proposes a differential genomic topological pattern between H37Rv and H37Ra, which could explain the virulence attenuation in H37Ra.

Cattle contribute to the nutritional needs and economy of a place. Performance and fitness depends on the response and adaptation to local climatic conditions. Genomic and genetic studies are important for advancing cattle breeding and availability of relevant reference genomes are essential. In the present study, the genome of a Vechur calf was sequenced on both short-read illumina and long-read nanopore platforms. Hybrid de novo assembly approach was deployed to obtain an average contigs length of 1.97 Mbp and N50 of 4.94 Mbp. Using a short-read genome sequence of the corresponding sire and dam, a haplotype-resolved genome was also assembled. In comparison to the taurine reference genome, we found 28982 autosomal structural variants, 16926990 SNVs with 883544 SNVs homozygous in the trio samples suggesting high prevalence of these SNPs in the population. Many of these SNPs have been reported to be associated with various QTLs including growth, milk yield and milk fat content which are crucial determinants of cattle production. Further, population genotype data analysis indicated that the present sample belongs to an Indian cattle breed forming a unique cluster of Bos indicus. Subsequent FST analysis revealed differentiation of Vechur cattle genome at multiple loci, especially those regions related to whole body growth and cell division especially IGF1, HMGA2, RRM2 and CD68 loci suggesting a possible role of these genes in its small stature and better disease resistance capabilities in comparison with the local crossbreeds. This provides an opportunity to select and engineer cattle breeds optimized for local conditions.