Search

Abstract Background In humans, copies of the Long Interspersed Nuclear Element 1 (LINE-1) retrotransposon comprise 21% of the reference genome, and have been shown to modulate expression and produce novel splice isoforms of transcripts from genes that span or neighbor the LINE-1 insertion site. Results In this work, newly released pilot data from the 1000 Genomes Project is analyzed to detect previously unreported full length insertions of the retrotransposon LINE-1. By direct analysis of the sequence data, we have identified 22 previously unreported LINE-1 insertion sites within the sequence data reported for a mother/father/daughter trio. Conclusions It is demonstrated here that next generation sequencing data, as well as emerging high quality datasets from individual genome projects allow us to assess the amount of heterogeneity with respect to the LINE-1 retrotransposon amongst humans, and provide us with a wealth of testable hypotheses as to the impact that this diversity may have on the health of individuals and populations.

Abstract Background High-throughput molecular biology techniques yield vast amounts of data, often by detecting small portions of ribonucleotides corresponding to specific identifiers. Existing bioinformatic methodologies categorize and compare these elements using inferred descriptive annotation given this sequence information irrespective of the fact that it may not be representative of the identifier as a whole. Results All annotations, no matter the granularity, can be aligned to genomic sequences and therefore annotated by genomic intervals. We have developed AbsIDconvert, a methodology for converting between genomic identifiers by first mapping them onto a common universal coordinate system using an interval tree which is subsequently queried for overlapping identifiers. AbsIDconvert has many potential uses, including gene identifier conversion, identification of features within a genomic region, and cross-species comparisons. The utility is demonstrated in three case studies: 1) comparative genomic study mapping plasmodium gene sequences to corresponding human and mosquito transcriptional regions; 2) cross-species study of Incyte clone sequences; and 3) analysis of human Ensembl transcripts mapped by Affymetrix®; and Agilent microarray probes. AbsIDconvert currently supports ID conversion of 53 species for a given list of input identifiers, genomic sequence, or genome intervals. Conclusion AbsIDconvert provides an efficient and reliable mechanism for conversion between identifier domains of interest. The flexibility of this tool allows for custom definition identifier domains contingent upon the availability and determination of a genomic mapping interval. As the genomes and the sequences for genetic elements are further refined, this tool will become increasingly useful and accurate. AbsIDconvert is freely available as a web application or downloadable as a virtual machine at: http://bioinformatics.louisville.edu/abid/.

Abstract Background Sequencing of the approximately 1.7 billion bases of the zebrafish genome is currently underway. To date, few high resolution genetic maps exist for the zebrafish genome, based mainly on single nucleotide polymorphisms (SNPs) and short microsatellite repeats. The desire to construct a higher resolution genetic map led to the construction of a database of tandemly repeating elements within the zebrafish Zv8 assembly. Description Exact tandem repeats with a repeat length of at least three bases and a copy number of at least 10 were reported. Repeats with a total length of 250 or fewer bases and their flanking regions were masked for known vertebrate repeats. Optimal primer pairs were computationally designed in the regions flanking the detected repeats. This database of exact tandem repeats can then be used as a resource by molecular biologists with interests in experimentally testing VNTRs within a zebrafish population. Conclusions A total of 116,915 repeats with a base length of at least three nucleotides were detected. The longest of these was a 54-base repeat with fourteen tandem copies. A significant number of repeats with a base length of 18, 24, 27 and 30 were detected, many with potentially novel proline-rich coding regions. Detection of exact tandem repeats in the zebrafish genome leads to a wealth of information regarding potential polymorphic sites for VNTRs. The association of many of these repeats with potentially novel yet similar coding regions yields an exciting potential for disease associated genes. A web interface for querying repeats is available at http://bioinformatics.louisville.edu/zebrafish/. This portal allows for users to search for a repeats of a selected base size from any valid specified region within the 25 linkage groups.

Abstract Background Detection of short, subtle conserved motif regions within a set of related DNA or amino acid sequences can lead to discoveries about important regulatory domains such as transcription factor and DNA binding sites as well as conserved protein domains. In order to help assess motif detection algorithms on motifs with varying properties and levels of conservation, we have developed a computational tool, rMotifGen, with the sole purpose of generating a number of random DNA or protein sequences containing short sequence motifs. Each motif consensus can be user-defined, randomly generated, or created from a position-specific scoring matrix (PSSM). Insertions and mutations within these motifs are created according to user-defined parameters and substitution matrices. The resulting sequences can be helpful in mutational simulations and in testing the limits of motif detection algorithms. Results Two implementations of rMotifGen have been created, one providing a graphical user interface (GUI) for random motif construction, and the other serving as a command line interface. The second implementation has the added advantages of platform independence and being able to be called in a batch mode. rMotifGen was used to construct sample sets of sequences containing DNA motifs and amino acid motifs that were then tested against the Gibbs sampler and MEME packages. Conclusion rMotifGen provides an efficient and convenient method for creating random DNA or amino acid sequences with a variable number of motifs, where the instance of each motif can be incorporated using a position-specific scoring matrix (PSSM) or by creating an instance mutated from its corresponding consensus using an evolutionary model based on substitution matrices. rMotifGen is freely available at: http://bioinformatics.louisville.edu/brg/rMotifGen/.

Abstract Background Enhancements in sequencing technology have recently yielded assemblies of large genomes including rat, mouse, human, fruit fly, and zebrafish. The availability of large-scale genomic and genic sequence data coupled with advances in microarray technology have made it possible to study the expression of large numbers of sequence products under several different conditions in days where traditional molecular biology techniques might have taken months, or even years. Therefore, to efficiently study a number of gene products associated with a disease, pathway, or other biological process, it is necessary to be able to design primer pairs or oligonucleotides en masse rather than using a time consuming and laborious gene-by-gene method. Results We have developed an integrated system, MPrime, in order to efficiently calculate primer pairs or specific oligonucleotides for multiple genic regions based on a keyword, gene name, accession number, or sequence fasta format within the rat, mouse, human, fruit fly, and zebrafish genomes. A set of products created for mouse housekeeping genes from MPrime-designed primer pairs has been validated using both PCR-amplification and DNA sequencing. Conclusion These results indicate MPrime accurately incorporates standard PCR primer design characteristics to produce high scoring primer pairs for genes of interest. In addition, sequence similarity for a set of oligonucleotides constructed for the same set of genes indicates high specificity in oligo design.

The majority of malignant tumors exhibit an altered metabolic phenotype that ultimately provides the required energy and molecular precursors necessary for unregulated cell division. Within this, phosphoserine aminotransferase 1 (PSAT1) is involved in de novo serine biosynthesis and its activity promotes various biochemical processes, including one-carbon metabolism. It also directly generates α-ketoglutarate (α-KG), a Kreb cycle intermediate and epigenetic-regulating metabolite. Prior studies examining PSAT1 depletion have identified individual affected downstream pathways, such as GSK3β and E2F, in several cancer types, including non-small-cell lung cancer (NSCLC). However, global gene expression examination in response to PSAT1 loss, particularly in EGFR mutant NSCLC, has not been unexplored. Transcriptional profiling of EGFR mutant NSCLC cells with or without stable knock-down of PSAT1 identified differentially expressed genes (DEGs) enriched in several metabolic pathways required for cell division, including amino acid and nucleotide biosynthesis. Supplementation studies involving non-essential amino acids, nucleosides and α-KG partially restored defects in anchorage-independent growth due to the knockdown of PSAT1. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment analysis identified potential impacts on actin cytoskeleton arrangement and β-catenin activity, which were rescued by PSAT1 re-expression. Finally, a comparative analysis of PSAT1 DEGs against transcripts enriched in patient EGFR mutant lung tumors identified a gene signature that is associated with overall and relapse-free survival (RFS) and was able to distinguish low or high-risk populations for RFS in early-stage EGFR mutant NSCLC. Overall, investigating genes altered by PSAT1 loss confirmed known PSAT1-regulated cellular pathways, identified a previously unknown role in the mediation of cytoskeleton arrangement in EGFR mutant NSCLC cells and allowed for the characterization of a gene signature with putative predictive potential for RFS in early-stage disease., Competing Interests: The authors declare that they have no competing interests., (Copyright: © 2024 Biyik-Sit et al.)

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis., (© 2024. The Author(s).)

The Long-read RNA-Seq Genome Annotation Assessment Project (LRGASP) Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. The consortium generated over 427 million long-read sequences from cDNA and direct RNA datasets, encompassing human, mouse, and manatee species, using different protocols and sequencing platforms. These data were utilized by developers to address challenges in transcript isoform detection and quantification, as well as de novo transcript isoform identification. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. When aiming to detect rare and novel transcripts or when using reference-free approaches, incorporating additional orthogonal data and replicate samples are advised. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis., Competing Interests: Competing Interests Design of the project was discussed with Oxford Nanopore Technologies (ONT), Pacific Biosciences, and Lexogen. ONT provided partial support for flow cells and reagents. S.C-S and A.N.B. have received reimbursement for travel, accommodation, and conference fees to speak at events organized by ONT. A.N.B. is a consultant for Remix Therapeutics, Inc.

Background: Diabetic complications have been studied extensively in recent years. There are very few biomarkers in body fluids that can pinpoint a distinct diabetic complication due to insufficient known specific biomarkers for ischemia., Objective: Identifying microRNA in animal models for each complication could enable early diagnosis of a given complication if verified in humans. MicroRNA (miRNA) profiling has been done in rodent models for a number of diabetic complications, like diabetic glomerular injury, atherosclerosis, cognitive impairment, diabetic wound healing, angiopathy and other complications. Due to multiple differences between rodents and humans, the changes in rabbit skin, considered closer to humans than even pigs, may better simulate human diabetic complications of ischemia., Methods: To study the miRNA profile of rabbits in which diabetes was induced or ischemia was surgically generated, we studied whether diabetes or ischemia-induced specific miRNA could be detected. MicroRNA from the blood of diabetic rabbits and rabbits with local ischemia was collected in PAXgene Blood RNA tubes specifically designed for miRNA isolation and extracted using the PAX gene miRNA extraction kit. The isolated RNA was quality controlled using an RNA analyzer, and further, using RNA seq technology, it was analyzed for distinct miRNAs that were detected in diabetic and non-diabetic rabbits induced with ischemia., Results: A miRNA that was found to be expressed in diabetic rabbits and ischemic rabbits but not in untreated rabbits was miRNA-183. Several miRNAs were differentially expressed across comparison groups, and several upregulated miRNAs were identified being unique to each comparison. In rabbits with a potential diabetic complication of a long-term ischemic model, there was one distinct microRNA, which was highly significantly upregulated in ischemia rabbit (miRNA-133-3p). One miRNA that was highly significantly upregulated in diabetic rabbit but not in ischemic rabbits was miRNA-3074-5p. Only statistically significant results have been considered and analyzed., Conclusion: These findings could lead to a precise and timely diagnosis of a potential single diabetic complication without invasive tissue biopsies and could be a novel tool in the management of diabetic patients developing complications due to the progression of diabetes., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

Epithelial-mesenchymal transition (EMT) facilitates cancer invasion and is initiated by mesenchyme-driving transcription factors and actin cytoskeletal assembly. We show a cytoplasmic-to-nuclear transport gradient of the EMT transcription factor Zeb1 toward sites of invasion in lung adenocarcinoma (LUAD), driven by the EMT inducer Tgfb, which is expressed in M2 polarized macrophages. We show that Zeb1 binds free actin monomers and RhoA in the cytoplasm to inhibit actin polymerization, blocking cell migration and Yap1 nuclear transport. Tgfb causes turnover of the scaffold protein Rassf1a, which targets RhoA. Release of this RhoA inhibition in response to Tgfb overcomes Zeb1's block of cytoskeleton assembly and frees it for nuclear transport. A ZEB1 nuclear transport signature highlights EMT progression, identifies dedifferentiated invasive/metastatic human LUADs, and predicts survival. Blocking Zeb1 nuclear transport with a small molecule identified in this study inhibits cytoskeleton assembly, cell migration, Yap1 nuclear transport, EMT, and precancerous-to-malignant transition., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

A recent surge in the use and abuse of diverse prescribed psychotic and illicit drugs necessitates the surveillance of drug residues in source water and the associated ecological impacts of chronic exposure to the aquatic organism. Thirty-six psychotic and illicit drug residues were determined in discharged wastewater from two centralized municipal wastewater treatment facilities and two wastewater receiving creeks for seven consecutive days in Kentucky. Zebrafish (Danio rerio) larvae were exposed to the environmental relevant mixtures of all drug residues, all illicit drugs, and all prescribed psychotic drugs. The extracted RNA from fish homogenates was sequenced, and differentially expressed sequences were analyzed for known or predicted nervous system expression, and screened annotated protein-coding genes to the true environmental cocktail mixture. Illicit stimulant (cocaine and one metabolite), opioids (methadone, methadone metabolite, and oxycodone), hallucinogen (MDA), benzodiazepine (oxazepam and temazepam), carbamazepine, and all target selective serotonin reuptake inhibitors including sertraline, fluoxetine, venlafaxine, and citalopram were quantified in 100% of collected samples from both creeks. The high dose cocktail mixture exposure group revealed the largest group of differentially expressed genes: 100 upregulated and 77 downregulated (p ≤ 0.05; q ≤ 0.05). The top 20 differentially expressed sequences in each exposure group comprise 82 unique transcripts corresponding to 74% annotated genes, 7% non-coding sequences, and 19% uncharacterized sequences. Among 61 differentially expressed sequences that corresponded to annotated protein-coding genes, 23 (38%) genes or their homologs are known to be expressed in the nervous system of fish or other organisms. Several of the differentially expressed sequences are associated primarily with the immune system, including several major histocompatibility complex class I and interferon-induced proteins. Interleukin-1 beta (downregulated in this study) abnormalities are considered a risk factor for psychosis. This is the first study to assess the contributions of multiple classes of psychotic and illicit drugs in combination with developmental gene expression., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Purpose: The management of triple-negative breast cancer (TNBC) remains a significant clinical challenge due to the lack of effective targeted therapies. Inhibitors of the cyclin-dependent kinases 4 and 6 (CDK4/6) are emerging as promising therapeutic agents against TNBC; however, cells can rapidly acquire resistance through multiple mechanisms that are yet to be identified. Therefore, determining the mechanisms underlying resistance to CDK4/6 inhibition is crucial to develop combination therapies that can extend the efficacy of the CDK4/6 inhibitors or delay resistance. This study aims to identify differentially expressed genes (DEG) associated with acquired resistance to palbociclib in ER- breast cancer cells., Methods: We performed next-generation transcriptomic sequencing (RNA-seq) and pathway analysis in ER- MDA-MB-231 palbociclib-sensitive (231/pS) and palbociclib-resistant (231/pR) cells., Results: We identified 2247 up-regulated and 1427 down-regulated transcripts in 231/pR compared to 231/pS cells. DEGs were subjected to functional analysis using Gene Ontology (GO) and the KEGG database which identified many transduction pathways associated with breast cancer, including the PI3K/AKT, PTEN and mTOR pathways. Additionally, Ingenuity Pathway Analysis (IPA) revealed that resistance to palbociclib is closely associated with altered cholesterol and fatty acid biosynthesis suggesting that resistance to palbociclib may be dependent on lipid metabolic reprograming., Conclusion: This study provides evidence that lipid metabolism is altered in TNBC with acquired resistance to palbociclib. Further studies are needed to determine if the observed lipid metabolic rewiring can be exploited to overcome therapy resistance in TNBC.

In this communication, we report on the genomic surveillance of SARS-CoV-2 using wastewater samples in Jefferson County, KY. In February 2021, we analyzed seven wastewater samples for SARS-CoV-2 genomic surveillance. Variants observed in smaller catchment areas, such as neighborhood manhole locations, were not necessarily consistent when compared to associated variant results in downstream treatment plants, suggesting catchment size or population could impact the ability to detect diversity.

Acquired resistance to cyclin-dependent kinases 4 and 6 (CDK4/6) inhibition in estrogen receptor-positive (ER+) breast cancer remains a significant clinical challenge. Efforts to uncover the mechanisms underlying resistance are needed to establish clinically actionable targets effective against resistant tumors. In this study, we sought to identify differentially expressed genes (DEGs) associated with acquired resistance to palbociclib in ER+ breast cancer. We performed next-generation transcriptomic RNA sequencing (RNA-seq) and pathway analysis in ER+ MCF7 palbociclib-sensitive (MCF7/pS) and MCF7 palbociclib-resistant (MCF7/pR) cells. We identified 2183 up-regulated and 1548 down-regulated transcripts in MCF7/pR compared to MCF7/pS cells. Functional analysis of the DEGs using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database identified several pathways associated with breast cancer, including 'cell cycle', 'DNA replication', 'DNA repair' and 'autophagy'. Additionally, Ingenuity Pathway Analysis (IPA) revealed that resistance to palbociclib is closely associated with deregulation of several key canonical and metabolic pathways. Further studies are needed to determine the utility of these DEGs and pathways as therapeutics targets against ER+ palbociclib-resistant breast cancer.

Glyceollins (Glys) are produced by soy plants in response to stress and are known for their anti-estrogenic activity both in vivo and in vitro in cancer cell lines as well as peripheral tissues. Glys can also exhibit non-estrogen receptor (ER) mediated effects. The effects of Glys on gene expression in the brain are still unclear. For this study, 17-β estradiol (E2) or placebo slow-release pellets were implanted into ovariectomized CFW mice followed by 11 days of exposure to either Glys or vehicle i.p. injections. We then examined the female mouse brain transcriptome using paired-end RNA sequencing (RNA-Seq) on the Illumina GAIIx platform. The goal of this study was to compare and contrast the results obtained from RNA-Seq with the results from our previous whole brain microarray experiment, which indicated that Glys potentially act through both ER-mediated and non-ER-mediated mechanisms, exhibiting a gene expression profile distinct from E2-treated groups. Our results suggest that the transcripts regulated by both E2 and Glys alone or in combination annotated to similar pathway maps and networks in both microarray and RNA-Seq experiments. Additionally, unlike our microarray data analysis, RNA-Seq enabled the detection of treatment effects on low expression transcripts of interest (e.g., prolactin and growth hormone). Collectively, our results suggest that depending on the gene, Glys can regulate expression independently of E2 action, similarly to E2, or oppose E2's effects in the female mouse brain., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Biological sequence usually contains yet to find knowledge, and mining biological sequences usually involves a huge dataset and long computation time. Common tasks for biological sequence mining are pattern discovery, classification and clustering. The newly developed model, Plausible Neural Network (PNN), provides an intuitive and unified architecture for such a large dataset analysis. This paper introduces the basic concepts of the PNN, and explains how it is applied to biological sequence mining. The specific task of biological sequence mining, exon/intron prediction, is implemented by using PNN. The experimental results show the capability of solving biological sequence mining tasks using PNN.

Untranslated regions (UTR) play important roles in the posttranscriptional regulation of mRNA processing. There is a wealth of UTR-related information to be mined from the rapidly accumulating EST collections. A computational tool, UTR-extender, has been developed to infer UTR sequences from genomically aligned ESTs. It can completely and accurately reconstruct 72% of the 3' UTRs and 15% of the 5' UTRs when tested using 908 functionally cloned transcripts. In addition, it predicts extensions for 11% of the 5' UTRs and 28% of the 3' UTRs. These extension regions are validated by examining splicing frequencies and conservation levels. We also developed a method called polyadenylation site scan (PASS) to precisely map polyadenylation sites in human genomic sequences. A PASS analysis of 908 genic regions estimates that 40-50% of human genes undergo alternative polyadenylation. Using EST redundancy to assess expression levels, we also find that genes with short 3' UTRs tend to be highly expressed.

In the course of our efforts to build extended regions of human genomic sequence by assembling individual BAC sequences, we have encountered several instances where a region of the genome has been sequenced independently using reagents derived from two different individuals. Comparing these sequences allows us to analyze the frequency and distribution of single nucleotide polymorphisms (SNPs) in the human genome. The observed transition/transversion frequencies are consistent with a biological origin for the sequence discrepancies, and this suggests that the data produced by large sequencing centers are accurate enough to be used as the basis for SNP analysis. The observed distribution of single nucleotide polymorphisms in the human genome is not uniform. An apparent duplication in the human genome extending over more than 130 kb between chromosomes 1p34 and 16p13 is reported. Independently derived sequences covering these regions are more than 99.9% identical, indicating that this duplication event must have occurred quite recently. FISH mapping results reported by the relevant laboratories indicate that the human population may be polymorphic for this duplication. We present a population genetic theory for the expected distribution of SNPs and derive an algorithm for probabilistically segmenting genomic sequence into regions that are identical by descent (IBD) between two individuals based on this theory and the observed locations of polymorphisms. Based on these methods and a random mating model for the human population, estimates are made for the mutation rate in the human genome.

DNA sequence analysis depends on the accurate assembly of fragment reads for the determination of a consensus sequence. This report examines the possibility of analyzing multiple, independent restriction digests as a method for testing the fidelity of sequence assembly. A dynamic programming algorithm to determine the maximum likelihood alignment of error prone electrophoretic mobility data to the expected fragment mobilities given the consensus sequence and restriction enzymes is derived and used to assess the likelihood of detecting rearrangements in genomic sequencing projects. The method is shown to reliably detect errors in sequence fragment assembly without the necessity of making reference to an overlying physical map. An html form-based interface is available at http:/(/)www.ibc.wustl.edu/services/validate. html.