Head and neck squamous cell carcinomas (HNSCC) include heterogeneous group of tumors, classified according to their anatomical site. It is the sixth most prevalent cancer globally. Among South Asian countries, India accounts for 40% of HNC malignancies with significant morbidity and mortality. In the present study, we have performed exome sequencing and analysis of 51 Head and Neck squamous cell carcinoma samples. Besides known mutations in the oncogenes and tumour suppressors, we have identified novel gene signatures differentiating buccal, alveolar, and tongue cancers. Around 50% of the patients showed mutation in tumour suppressor genes TP53 and TP63. Apart from the known mutations, we report novel mutations in the genes AKT1, SPECC1, and LRP1B, which are linked with tumour progression and patient survival. A highly curated process was developed to identify survival signatures. 36 survival-related genes were identified based on the correlation of functional impact of variants identified using exome-seq with gene expression from transcriptome data (GEPIA database) and survival. An independent LASSO regression analysis was also performed. Survival signatures common to both the methods led to identification of 4 dead and 3 alive gene signatures, the accuracy of which was confirmed by performing a ROC analysis (AUC=0.79 and 0.91, respectively). Also, machine learning-based driver gene prediction tool resulted in the identification of IRAK1 as the driver (p-value = 9.7 e-08) and also as an actionable mutation. Modelling of the IRAK1 mutation showed a decrease in its binding to known IRAK1 inhibitors.

Anopheles stephensi is the most menacing malaria vector to watch for in newly urbanising parts of the world. Its fitness is reported to be a direct consequence of the vector adapting to laying eggs in over-head water tanks with street-side water puddles polluted by oil and sewage. Large frequent inversions in the genome of malaria vectors are implicated in adaptation. We report the genome assembly of a strain of An. stephensi of the type-form, collected from a construction site from Chennai (IndCh) in 2016. The genome reported here with a L50 of 4, completes the trilogy of high-resolution genomes of strains with respect to a 16.5 Mbp 2Rb genotype in An. stephensi known to be associated with adaptation to environmental heterogeneity. Unlike the reported genomes of two other strains, STE2 (2R+b/2Rb) and UCI (2Rb/2Rb), IndCh is found to be homozygous for the standard form (2R+b/2R+b). Comparative genome analysis revealed base-level details of the breakpoints and allowed extraction of 22,650 segregating SNPs for typing this inversion in populations. Whole genome sequencing of 82 individual mosquitoes from diverse geographical locations reveal that one third of both wild and laboratory populations maintain the heterozygous genotype of 2Rb. The large number of SNPs can be tailored to 1740 exonic SNPs enabling genotyping directly from transcriptome sequencing. The genome trilogy approach accelerated the study of fine structure and typing of an important inversion in An. stephensi, putting the genome resources for this understudied species on par with the extensively studied malaria vector, Anopheles gambiae. We argue that the IndCh genome is relevant for field translation work compared to those reported earlier by showing that individuals from diverse geographical locations cluster with IndCh, pointing to significant convergence resulting from travel and commerce between cities, perhaps, contributing to the survival of the fittest strain.

Grain amaranths, edible C4 dicots, produce pseudo-cereals high in lysine. Lysine being one of the most limiting essential amino acids in cereals and C4 photosynthesis being one of the most sought-after phenotypes in protein-rich legume crops, the genome of one of the grain amaranths is likely to play a critical role in crop research. We have sequenced the genome and transcriptome of Amaranthus hypochondriacus, a diploid (2n = 32) belonging to the order Caryophyllales with an estimated genome size of 466 Mb. Of the 411 linkage single-nucleotide polymorphisms (SNPs) reported for grain amaranths, 355 SNPs (86%) are represented in the scaffolds and 74% of the 8.6 billion bases of the sequenced transcriptome map to the genomic scaffolds. The genome of A. hypochondriacus, codes for at least 24,829 proteins, shares the paleohexaploidy event with species under the superorders Rosids and Asterids, harbours 1 SNP in 1,000 bases, and contains 13.76% of repeat elements. Annotation of all the genes in the lysine biosynthetic pathway using comparative genomics and expression analysis offers insights into the high-lysine phenotype. As the first grain species under Caryophyllales and the first C4 dicot genome reported, the work presented here will be beneficial in improving crops and in expanding our understanding of angiosperm evolution.

Edible/non-toxic varieties of Jatropha curcas L. are gaining increasing attention, providing both oil as biofuel feedstock or even as edible oil and the seed kernel meal as animal feed ingredient. They are a viable alternative to the limitation posed by the presence of phorbol esters in toxic varieties. Accurate genotyping of toxic/non-toxic accessions is critical to breeding management. The aim of this study was to identify SNP markers linked to seed toxicity in J. curcas. For SNP discovery, NGS technology was used to sequence the whole genomes of a toxic and non-toxic parent along with a bulk of 51 toxic and 30 non-toxic F2 plants. To ascertain the association between SNP markers and seed toxicity trait, candidate SNPs were genotyped on 672 individuals segregating for seed toxicity and two collections of J. curcas composed of 96 individuals each. In silico SNP discovery approaches led to the identification of 64 candidate SNPs discriminating non-toxic and toxic samples. These SNPs were mapped on Chromosome 8 within the Linkage Group 8 previously identified as a genomic region important for phorbol ester biosynthesis. The association study identified two new SNPs, SNP_J22 and SNP_J24 significantly linked to low toxicity with R2 values of 0.75 and 0.54, respectively. Our study released two valuable SNP markers for high-throughput, marker-assisted breeding of seed toxicity in J. curcas.

Malaria remains a major healthcare risk to growing economies like India, and a chromosome-level reference genome of Anopheles stephensi is critical for successful vector management and understanding of vector evolution using comparative genomics. We report chromosome-level assemblies of an Indian strain, STE2, and a Pakistani strain SDA-500 by combining draft genomes of the two strains using a homology-based iterative approach. The resulting assembly IndV3/PakV3 with L50 of 9/12 and N50 6.3/6.9 Mb had scaffolds long enough for building 90% of the euchromatic regions of the three chromosomes, IndV3s/PakV3s, using low-resolution physical markers and enabled the generation of the next version of genome assemblies, IndV4/PakV4, using HiC data. We have validated these assemblies using contact maps against publicly available HiC raw data from two strains including STE2 and another lab strain of An. stephensi from UCI and compare the quality of the assemblies with other assemblies made available as preprints since the submission of the manuscript. We show that the IndV3s and IndV4 assemblies are sensitive in identifying a homozygous 2Rb inversion in the UCI strain and a 2Rb polymorphism in the STE2 strain. Multiple tandem copies of CYP6a14, 4c1, and 4c21 genes, implicated in insecticide resistance, lie within this inversion locus. Comparison of assembled genomes suggests a variation of 1 in 81 positions between the UCI and STE2 lab strains, 1 in 82 between SDA-500 and UCI strain, and 1 in 113 between SDA-500 and STE2 strains of An. stephensi, which are closer than 1 in 68 variations among individuals from two other lab strains sequenced and reported here. Based on the developmental transcriptome and orthology of all the 54 olfactory receptors (ORs) to those of other Anopheles species, we identify an OR with the potential for host recognition in the genus Anopheles. A comparative analysis of An. stephensi genomes with the completed genomes of a few other Anopheles species suggests limited inter-chromosomal gene flow and loss of synteny within chromosomal arms even among the closely related species.

Recent studies have implicated the hypocretin/orexinergic system in reward-seeking behavior. Almorexant, a dual orexin/hypocretin R(1) and R(2) receptor antagonist, has proven effective in preclinical studies in promoting sleep in animal models and was in Phase III clinical trials for sleep disorders. The present study combines behavioral assays with in vitro biochemical and electrophysiological techniques to elucidate the role of almorexant in ethanol and sucrose intake. Using an operant self-administration paradigm, we demonstrate that systemic administration of almorexant decreased operant self-administration of both 20% ethanol and 5% sucrose. We further demonstrate that intra-ventral tegmental area (VTA) infusions, but not intra-substantia nigra infusions, of almorexant reduced ethanol self-administration. Extracellular recordings performed in VTA neurons revealed that orexin-A increased firing and this enhancement of firing was blocked by almorexant. The results demonstrate that orexin/hypocretin receptors in distinct brain regions regulate ethanol and sucrose mediated behaviors.

Identification of Plasmodium-resistance genes in malaria vectors remains an elusive goal despite the recent availability of high-quality genomes of several mosquito vectors. Anopheles stephensi, with its three distinctly-identifiable forms at the egg stage, correlating with varying vector competence, offers an ideal species to discover functional mosquito genes implicated in Plasmodium resistance. Recently, the genomes of several strains of An. stephensi of the type-form, known to display high vectorial capacity, were reported. Here, we report a chromosomal-level assembly of an intermediate-form of An. stephensi strain (IndInt), shown to have reduced vectorial capacity relative to a strain of type-form (IndCh). The contig level assembly with a L50 of 4 was scaffolded into chromosomes by using the genome of IndCh as the reference. The final assembly shows a heterozygous paracentric inversion, 3Li, involving 8 Mbp, which is syntenic to the extensively-studied 2La inversion implicated in Plasmodium resistance in An. gambiae involving 21 Mbp. Deep annotation of genes within the 3Li region in the IndInt assembly using the state-of-the-art protein-fold prediction and other annotation tools reveals the presence of a tumor necrosis factor-alpha (TNF-alpha) like gene, which is the homolog of the Eiger gene in Drosophila. Subsequent chromosome-wide searches revealed homologs of Wengen (Wgn) and Grindelwald (Grnd) genes, which are known to be the receptors for Eiger in Drosophila. We have identified all the genes in IndInt required for Eiger-mediated signaling by analogy to the TNF-alpha system, suggesting the presence of a functionally-active Eiger signaling pathway in IndInt. Comparative genomics of the three type-forms with that of IndInt, reveals structurally disruptive mutations in Eiger gene in all three strains of the type-form, suggesting compromised innate immunity in the type-form as the likely cause of high vectorial capacity in these strains. This is the first report of the presence of a homolog of Eiger in malaria vectors, known to be involved in cell death in Drosophila, within an inversion region in IndInt syntenic to an inversion associated with Plasmodium resistance in An. gambiae.

The hg19 assembly of the human genome is the most heavily annotated and most commonly used reference to make variant calls for individual genomes. Based on the phase 3 report of the 1000 genomes project (1000G), it is now well known that many positions in the hg19 genome represent minor alleles. Since commonly used variant call methods are developed under the assumption that hg19 reference harbors major alleles at all the ~3 billion positions, these methods mask the calls whenever an individual is homozygous to the minor allele at the respective positions. Hence, it is important to address the extent and impact of these minor alleles in hg19 from the point of view of individual genomes.

Recently, large numbers of normal human tissues have been profiled for non-coding RNAs and more than fourteen thousand long intergenic non-coding RNAs (lincRNAs) are found expressed in normal human tissues. The functional roles of these normal lincRNAs (nlincRNAs) in the regulation of protein coding genes in normal and disease biology are yet to be established. Here, we have profiled two RNA-seq datasets including cancer and matched non-neoplastic tissues from 12 individuals from diverse demography for both coding genes and nlincRNAs. We find 130 nlincRNAs significantly regulated in cancer, with 127 regulated in the same direction in the two datasets. Interestingly, according to Illumina Body Map, significant numbers of these nlincRNAs display baseline null expression in normal prostate tissues but are specific to other tissues such as thyroid, kidney, liver and testis. A number of the regulated nlincRNAs share loci with coding genes, which are either co-regulated or oppositely regulated in all cancer samples studied here. For example, in all cancer samples i) the nlincRNA, TCONS_00029157, and a neighboring tumor suppressor factor, SIK1, are both down regulated; ii) several thyroid-specific nlincRNAs in the neighborhood of the thyroid-specific gene TPO, are both up-regulated; and iii) the TCONS_00010581, an isoform of HEIH, is down-regulated while the neighboring EZH2 gene is up-regulated in cancer. Several nlincRNAs from a prostate cancer associated chromosomal locus, 8q24, are up-regulated in cancer along with other known prostate cancer associated genes including PCAT-1, PVT1, and PCAT-92. We observe that there is significant bias towards up-regulation of nlincRNAs with as high as 118 out of 127 up-regulated in cancer, even though regulation of coding genes is skewed towards down-regulation. Considering that all reported cancer associated lincRNAs (clincRNAs) are biased towards up-regulation, we conclude that this bias may be functionally relevant.

The most widely used human genome reference assembly hg19 harbors minor alleles at 2.18 million positions as revealed by 1000 Genome Phase 3 dataset. Although this is less than 2% of the 89 million variants reported, it has been shown that the minor alleles can result in 30% false positives in individual genomes, thus misleading and burdening downstream interpretation. More alarming is the fact that, significant percentage of variants that are homozygous recessive for these minor alleles, with potential disease implications, are masked from reporting.

The role of many lncRNAs in cancer remains elusive including that for a Prostate Cancer Associated Transcript 92 (PCAT92). PCAT92 shares the locus on chromosome 13 with ABCC4 gene, known to be implicated in prostate cancer. It has been shown that PCAT92 and ABCC4 are up-regulated in prostate cancer samples from multiple transcriptome datasets. Among the prostate cancer cell-lines LNCaP showed maximum overexpression of PCAT92 compared to control cell-line RWPE-1. We have shown that knockdown of PCAT92 in LNCaP cells reduces cell viability and proliferation and down-regulates ABCC4 transcript/protein expression. The shared region between PCAT92 and ABCC4 has a binding site for an oncogenic transcription factor (ZIC2) which is also upregulated in the majority of datasets studied here. ZIC2 binding to the predicted ABCC4 promoter has been confirmed using pull-down assay. Interestingly, under PCAT92 knockdown condition, there is a reduction in the ZIC2 binding to ABCC4 promoter indicating the potential involvement of PCAT92 in the recruitment of ZIC2. We have identified distinct regions on PCAT92 with potential to bind to ZIC2 non-DNA binding Zinc-finger domain and potential for triplex formation near ABCC4 promoter region, which have been experimentally validated. Together, these observations and localization in the nucleus suggests that PCAT92 may play a role in prostate cancer by increasing the local concentration of ZIC2 by forming RNA-DNA triplex near ABCC4 promoter thus helping in recruitment of ZIC2 for ABCC4 regulation.

In the age of genomics-based crop improvement, a high-quality genome of a local landrace adapted to the local environmental conditions is critically important. Grain amaranths produce highly nutritional grains with a multitude of desirable properties including C4 photosynthesis highly sought-after in other crops. For improving the agronomic traits of grain amaranth and for the transfer of desirable traits to dicot crops, a reference genome of a local landrace is necessary. Toward this end, our lab had initiated sequencing the genome of Amaranthus (A.) hypochondriacus (A.hyp_K_white) and had reported a draft genome in 2014. We selected this landrace because it is well adapted for cultivation in India during the last century and is currently a candidate for TILLING-based crop improvement. More recently, a high-quality chromosome-level assembly of A. hypochondriacus (PI558499, Plainsman) was reported. Here, we report a chromosome-level assembly of A.hyp_K_white (AhKP) using low-coverage PacBio reads, contigs from the reported draft genome of A.hyp_K_white, raw HiC data and reference genome of Plainsman (A.hyp.V.2.1). The placement of A.hyp_K_white on the phylogenetic tree of grain amaranths of known accessions clearly suggests that A.hyp_K_white is genetically distal from Plainsman and is most closely related to the accession PI619259 from Nepal (Ramdana). Furthermore, the classification of another accession, Suvarna, adapted to the local environment and selected for yield and other desirable traits, is clearly Amaranthus cruentus. A classification based on hundreds of thousands of SNPs validated taxonomy-based classification for a majority of the accessions providing the opportunity for reclassification of a few.

The identification of efficient molecular markers related to low bolting tendency is a priority in sugar beet (Beta vulgaris L.) breeding. This study aimed to identify SNP markers associated with low bolting tendency by establishing a genome-wide association study. An elaborate 3-year field trial comprising 13 sugar beet lines identified L14 as the one exhibiting the lowest bolting tendency along with an increased survival rate after autumnal sowing. For SNP discovery following phenotyping, contrasting phenotypes of 24 non-bolting and 15 bolting plants of the L14 line were sequenced by restriction site-associated DNA sequencing (RAD-seq). An association model was established with a set of 10,924 RAD-based single nucleotide polymorphism (SNP) markers. The allelic status of the most significantly associated SNPs ranked based on their differential allelic status between contrasting phenotypes (p < 0.01) was confirmed on three different validation datasets comprising diverse sugar beet lines and varieties adopting a range of SNP detection technologies. This study has led to the identification of SNP_36780842 and SNP_48607347 linked to low bolting tendency and can be used for marker-assisted breeding and selection in sugar beet.

Genome duplication event in edible dicots under the orders Rosid and Asterid, common during the oligocene period, is missing for species under the order Caryophyllales. Despite this, grain amaranths not only survived this period but display many desirable traits missing in species under rosids and asterids. For example, grain amaranths display traits like C4 photosynthesis, high-lysine seeds, high-yield, drought resistance, tolerance to infection and resilience to stress. It is, therefore, of interest to look for minor genome rearrangements with potential functional implications that are unique to grain amaranths. Here, by deep sequencing and assembly of 16 transcriptomes (86.8 billion bases) we have interrogated differential genome rearrangement unique to Amaranthus hypochondriacus with potential links to these phenotypes. We have predicted 125,581 non-redundant transcripts including 44,529 protein coding transcripts identified based on homology to known proteins and 13,529 predicted as novel/amaranth specific coding transcripts. Of the protein coding de novo assembled transcripts, we have identified 1810 chimeric transcripts. More than 30% and 19% of the gene pairs within the chimeric transcripts are found within the same loci in the genomes of A. hypochondriacus and Beta vulgaris respectively and are considered real positives. Interestingly, one of the chimeric transcripts comprises two important genes, namely DHDPS1, a key enzyme implicated in the biosynthesis of lysine, and alpha-glucosidase, an enzyme involved in sucrose catabolism, in close proximity to each other separated by a distance of 612 bases in the genome of A. hypochondriacus in a convergent configuration. We have experimentally validated that transcripts of these two genes are also overlapping in the 3' UTR with their expression negatively correlated from bud to mature seed, suggesting a potential link between the high seed lysine trait and unique genome organization.

Potassium currents play key roles in regulating motoneuron activity, including functional specializations that are important for locomotion. The thoracic and abdominal segments in the Drosophila larval ganglion have repeated arrays of motoneurons that innervate body-wall muscles used for peristaltic movements during crawling. Although abdominal motoneurons and their muscle targets have been studied in detail, owing, in part, to their involvement in locomotion, little is known about the cellular properties of motoneurons in thoracic segments. The goal of this study was to compare firing properties among thoracic motoneurons and the potassium currents that influence them. Whole-cell, patch-clamp recordings performed from motoneurons in two thoracic and one abdominal segment revealed both transient and sustained voltage-activated K(+) currents, each with Ca(++)-sensitive and Ca(++)-insensitive [A-type, voltage-dependent transient K(+) current (I(Av))] components. Segmental differences in the expression of voltage-activated K(+) currents were observed. In addition, we demonstrate that Shal contributes to I(Av) currents in the motoneurons of the first thoracic segment.

The mosquito Anopheles stephensi is a vector of urban malaria in Asia that recently invaded Africa. Studying the genetic basis of vectorial capacity and engineering genetic interventions are both impeded by limitations of a vector's genome assembly. The existing assemblies of An. stephensi are draft-quality and contain thousands of sequence gaps, potentially missing genetic elements important for its biology and evolution.

Vector management programs rely on knowledge of the biology and genetic make-up of mosquitoes. Anopheles stephensi is a major invasive urban malaria vector, distributed throughout the Indian subcontinent and Middle East, and has recently been expanding its range in Africa. With the existence of three biological forms, distinctly identifiable based on the number of ridges on eggs and varying vectorial competence, An. stephensi is a perfect species for developing isofemale lines, which can be tested for insecticide susceptibility and vectorial competence of various biological forms.